Structure learning of Hamiltonians from real-time evolution
Pith reviewed 2026-05-24 01:14 UTC · model grok-4.3
The pith
An algorithm recovers the full structure and parameters of an unknown local Hamiltonian from real-time evolution queries using total time O(log n / ε).
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We present a new general approach to Hamiltonian learning that solves the structure-learning variant by recovering the Hamiltonian to ε error with total evolution time O(log n / ε). The algorithm does not need to know the Hamiltonian terms in advance, works for any H where the sum of terms interacting with a qubit has bounded norm, and evolves according to H in constant time t increments. As an application it also learns Hamiltonians with power-law decay to accuracy ε with total evolution time that beats the usual 1/ε² scaling.
What carries the argument
A general approach to structure learning that identifies unknown interaction terms by querying the evolution operator in fixed-length constant-time increments under a per-qubit bounded-norm condition on interactions.
If this is right
- The method extends Hamiltonian learning beyond the short-range setting to any Hamiltonian obeying the bounded per-qubit interaction norm.
- It achieves Heisenberg-limited scaling for the harder structure-learning problem rather than only parameter estimation.
- It enables learning of power-law decaying interactions with total evolution time better than the standard 1/ε² limit.
- All queries use constant-time evolution steps, removing the need for arbitrarily short time resolutions.
Where Pith is reading between the lines
- Experimental platforms could apply the same constant-time queries even when the underlying interaction graph is completely unknown in advance.
- The bounded-norm condition may be checkable in advance on physical devices by measuring local energy scales.
- The approach could be combined with existing short-range learning routines to handle hybrid short- and long-range systems.
Load-bearing premise
The sum of all terms interacting with any single qubit has bounded norm.
What would settle it
Construct a family of Hamiltonians on n qubits where the per-qubit interaction norm grows with n and check whether the algorithm still recovers the Hamiltonian to error ε using only O(log n / ε) total evolution time.
Figures
read the original 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 \lambda_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, $\lambda_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$.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a new algorithm for structure learning of an unknown local Hamiltonian H = sum λ_a E_a from real-time evolution oracles e^{-i H t}. It claims to recover H to ε error using total evolution time O(log n / ε), without prior knowledge of the interaction terms E_a, for any Hamiltonian where the sum of norms of terms interacting with a qubit is bounded; the algorithm uses constant-time t increments and extends to power-law decaying interactions with scaling better than the standard 1/ε² limit.
Significance. If the central claims and derivations hold, the result would be significant for quantum Hamiltonian learning: it resolves the structure-learning variant (previously open without term knowledge) at Heisenberg-limited scaling, provides a general method beyond short-range interactions, and achieves constant time resolution. The power-law application is a concrete strength. No machine-checked proofs or reproducible code are mentioned in the provided description.
minor comments (2)
- [Abstract] Abstract: the phrase 'gold standard of Heisenberg-limited scaling' would benefit from an explicit comparison (in the introduction or §1) to the best prior bounds for the structure-learning setting to clarify the improvement.
- The bounded-norm condition is stated as the operating regime; a brief remark in the main theorem statement on whether the O(log n / ε) bound is tight or contains hidden constants dependent on the norm bound would improve clarity.
Simulated Author's Rebuttal
We thank the referee for their positive assessment of our work on Hamiltonian structure learning and for recommending minor revision. The report correctly identifies the key contributions, including the resolution of the structure-learning problem at Heisenberg-limited scaling without prior knowledge of the interaction terms, the extension beyond short-range interactions, and the application to power-law decaying Hamiltonians. Since the report lists no specific major comments or requested changes, we have no points to address point-by-point at this stage. We will incorporate any minor editorial suggestions during revision.
Circularity Check
No significant circularity identified
full rationale
The paper's central claim is an algorithm recovering H to ε error in O(log n / ε) total evolution time from the e^{-iHt} oracle, under the explicit operating assumption that the sum of norms of terms interacting with any qubit is bounded (and without prior knowledge of the E_a terms). This bounded-norm condition is stated as the regime of applicability rather than derived internally; the scaling, constant-t increments, and extension beyond short-range are presented as consequences of the oracle access and assumption. No self-definitional reductions, fitted inputs renamed as predictions, or load-bearing self-citation chains appear in the provided abstract or description. The derivation is therefore self-contained against external benchmarks given in the problem statement.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The unknown Hamiltonian satisfies that the sum of terms interacting with any qubit has bounded norm.
- domain assumption Access to the real-time evolution operator e^{-i H t} for chosen t is available as an oracle.
Forward citations
Cited by 2 Pith papers
-
Score-matching-based Structure Learning for Temporal Data on Networks
PICK adds a parent-finding subroutine for leaf nodes to speed up pruning in score-matching causal discovery, extending it from i.i.d. data to static and temporal network data.
-
Stability of digital and analog quantum simulations under noise
Rigorous worst- and average-case error bounds show comparable worst-case scaling for digital and analog quantum simulators under perturbative noise, with distinct average-case error cancellation and concentration boun...
Reference graph
Works this paper leans on
-
[1]
Nearly optimal lattice simulation by product formulas
arXiv:2201.11211 [stat.ML](page 14). [CS19] Andrew M. Childs and Yuan Su. “Nearly optimal lattice simulation by product formulas”. In: Physical Review Letters 123.5 (Aug. 2019), p. 050503. ISSN : 1079-
-
[2]
Theory of Trotter error with commutator scaling
DOI: 10.1103/physrevlett.123.050503 . arXiv: 1901.00564 [quant-ph] (page 5). [CSTWZ21] Andrew M. Childs, Yuan Su, Minh Tran, Nathan Wiebe, and Shuchen Zhu. “Theory of Trotter error with commutator scaling”. In: Physical Review X 11.1 (Feb. 2021), p. 011020. ISSN : 2160-3308. DOI: 10.1103/PhysRevX.11.011020 . arXiv:1912.08854 [quant-ph](pages 8, 9). [DDMPR...
-
[3]
Scalable reconstruction of unitary processes and Hamiltonians
ISSN : 0031-9007. DOI:10.1103/physrevlett.71.1355 (page 1). [HBCP15] M. Holzäpfel, T. Baumgratz, M. Cramer, and M. B. Plenio. “Scalable reconstruc- tion of unitary processes and hamiltonians”. In: Physical Review A 91.4 (Apr. 2015), p. 042129. ISSN : 1094-1622. DOI:10.1103/physreva.91.042129. quant-ph: 1411.6379 (arXiv) (page 13). [HCP23] Hsin-Yuan Huang,...
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physrevlett.71.1355 2015
-
[4]
Gradient Descent Learns Linear Dynamical Systems
arXiv: 1609.05191 [cs.LG]. URL: http://jmlr.org/papers/v19/16-465.html (page 14). [HTFS23] Hsin-Yuan Huang, Yu Tong, Di Fang, and Yuan Su. “Learning many-body Hamiltonians with Heisenberg-limited scaling”. In: Physical Review Letters 130.20 (May 2023), p. 200403. DOI: 10.1103/physrevlett.130.200403 . arXiv: 2210.03030 [quant-ph](pages 2–6, 13–15, 42). [KL...
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1103/physrevlett.130.200403 2023
-
[5]
Cambridge University Press, doi:10.1017/CBO9781139814782
ISBN : 9781107471542. DOI:10.1017/cbo9781139814782. arXiv:2105.10386 [cs.DM] (pages 11, 35). [OKSM24] Tatsuki Odake, Hlér Kristjánsson, Akihito Soeda, and Mio Murao. “Higher- order quantum transformations of Hamiltonian dynamics”. In: Physical Review Research 6.1 (Mar. 2024), p. l012063.ISSN : 2643-1564. DOI:10.1103/physrevresearch. 6.l012063. arXiv:2303....
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.