Structure learning of Hamiltonians from real-time evolution
Ainesh Bakshi, Allen Liu, Ankur Moitra, Ewin Tang
TL;DR
The paper tackles Hamiltonian structure learning from real-time evolution for an unknown $\mathfrak{K}$-local Hamiltonian on $n$ qubits, aiming to recover all coefficients without prior knowledge of interaction terms. It introduces a recursive, bootstrapped framework that achieves Heisenberg-limited scaling, with total evolution time $t_{total}=O(\mathfrak{r}\log n/\varepsilon)$ and constant time resolution $t_{min}=\Theta(1/\mathfrak{r})$, while handling long-range and power-law decaying interactions via bounded local norms. A key innovation is a constant-time Trotterization bound enabling discrete control to replace continuous quantum control, together with a structure-learning layer that uses Goldreich–Levin–style queries on Pauli spectra to identify large coefficients efficiently. The method is applicable to both bounded-range and non-local Hamiltonians, including power-law decays that beat the standard $1/\varepsilon^2$ scaling, and it yields a fixed-parameter tractable runtime in many realistic settings. Overall, the work provides a general, scalable framework for quantum-device characterization that bridges quantum metrology, Hamiltonian simulation, and classical structure-learning paradigms.
Abstract
We study the problem of Hamiltonian structure learning from real-time evolution: given the ability to apply $e^{-\mathrm{i} Ht}$ for an unknown local Hamiltonian $H = \sum_{a = 1}^m λ_a E_a$ on $n$ qubits, the goal is to recover $H$. This problem is already well-understood under the assumption that the interaction terms, $E_a$, are given, and only the interaction strengths, $λ_a$, are unknown. But how efficiently can we learn a local Hamiltonian without prior knowledge of its interaction structure? We present a new, general approach to Hamiltonian learning that not only solves the challenging structure learning variant, but also resolves other open questions in the area, all while achieving the gold standard of Heisenberg-limited scaling. In particular, our algorithm recovers the Hamiltonian to $\varepsilon$ error with total evolution time $O(\log (n)/\varepsilon)$, and has the following appealing properties: (1) it does not need to know the Hamiltonian terms; (2) it works beyond the short-range setting, extending to any Hamiltonian $H$ where the sum of terms interacting with a qubit has bounded norm; (3) it evolves according to $H$ in constant time $t$ increments, thus achieving constant time resolution. As an application, we can also learn Hamiltonians exhibiting power-law decay up to accuracy $\varepsilon$ with total evolution time beating the standard limit of $1/\varepsilon^2$.
