REVIEW 3 major objections 5 minor 48 references
A 97-qubit superconducting processor measures magnon spectra, lifetimes, and scattering channels in a 2D quantum magnet with precision that classical tensor-network simulations cannot match.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-02 05:34 UTC pith:N4ITM673
load-bearing objection A convincing experimental tour de force with one untested assumption in the 97-qubit mode identification; deserves serious review with requests for an overlap check and tighter fits. the 3 major comments →
Precision quantum simulation of magnon spectra and interactions
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that a 97-qubit analog-digital processor can perform high-precision spectroscopy of individual magnon modes in a 2D XY magnet, extracting linear response (temperature-dependent spectra and lifetimes) and nonlinear response (self-scattering, pump-probe, and coherent two-mode scattering). The measurements reveal that magnon decay rates vary strongly across the Brillouin zone: enhancement near van Hove singularities, and suppression for corner-localized modes; that the decay rate follows a mode-dependent power law Γ∝(ε−ε0)^η whose exponent falls below unity at high mode indices, which the authors attribute to self-consistent broadening of the scattering partners; and that c
What carries the argument
The load-bearing technique is the interleaving of analog Hamiltonian evolution with digital single-qubit Z-gates on a processor whose effective Hamiltonian (including small multi-qubit terms) has been learned to about 0.1% of the dominant coupling by optimizing cross-entropy benchmarking in 24-qubit patches and patching the result to the full device. Magnon modes are defined through the spin-wave (Holstein-Primakoff) Bogoliubov transformation of the mean-field Hamiltonian, giving mode shapes used both to excite (via phase-imprinted Z-rotations) and to read out (via weighted measurement) the retarded Green's function of each mode. Decay rates are modeled with a self-consistent Born approximat
Load-bearing premise
The 97-qubit claims rest on the assumption that the Hamiltonian learned at infinite temperature in 24-qubit patches, then patched to the full device, remains accurate at the low temperatures and 97-qubit scale used for the magnon measurements; if small long-range or temperature-dependent terms are missed, the extracted mode shapes and decay rates would mix several excitations.
What would settle it
Measure the same mode-resolved decay rates at 97 qubits using a Hamiltonian calibrated directly at the operating temperature (e.g., by learning from bitstrings of a cooled, tilted initial state close to the magnon operating point) and compare with the patched learned-Hamiltonian predictions; a mismatch in the mode-resolved decay rates beyond the reported median error, or a difference in the corner-mode pairing pattern, would falsify the largest-system claims. Alternatively, an exact state-vector simulation at 36-40 qubits with the learned Hamiltonian that fails to reproduce the measured 24-31
If this is right
- If the learned Hamiltonian is accurate at the operating temperature, the same interleaved protocol can map the full mode-resolved response function of other 2D spin models on the processor, including anisotropic and frustrated couplings.
- The observed sublinear exponent η<1 for high-energy modes establishes self-consistent broadening as a measurable feedback effect, giving a concrete target for analytic and numerical many-body theory beyond the Born approximation.
- The finding that corner modes form a nearly closed scattering manifold implies that boundary-localized magnons can act as long-lived, weakly coupled information carriers, a feature that could matter for spintronic or magnonic device design.
- The paper's evidence that MPS and PEPS cannot reach experimental precision at 97 qubits and intermediate temperatures, if correct, supplies a concrete quantum-advantage benchmark: a classically hard observable (the magnon decay rate) that a quantum processor can deliver.
Where Pith is reading between the lines
- Because the paper validates the learned Hamiltonian only at infinite temperature and in 24-qubit patches, an immediate extension is to test the same Hamiltonian at the operating temperatures via low-temperature cross-entropy benchmarking or density-matrix reconstruction; success would strengthen confidence in the 97-qubit mode-resolved results.
- The corner-mode manifold suggests a design principle the paper leaves implicit: in finite open-boundary spin lattices, modes with low effective coordination have strongly reduced scattering phase space, so geometric engineering of boundaries (e.g., notches or holes) could selectively protect or destabilize particular magnon modes.
- The power-law exponent η, extracted from an energy-density ramp, is a new spectral observable that could be measured in other platforms, such as Rydberg arrays or trapped ions, as a test of whether the self-consistent broadening mechanism is universal or specific to the XY model.
- If the learned Hamiltonian is ported to another device with the same architecture, the mode-resolved decay maps (e.g., the cross pattern in pump-probe data) could serve as a sensitive fingerprint for diagnosing small parameter drift in the effective Hamiltonian.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports a hybrid analog-digital quantum simulation of magnon excitations in a 2D XY spin-1/2 model on a superconducting processor with up to 97 qubits. The authors combine Hamiltonian learning via XEB (with per-qubit-per-cycle errors near 5e-4) with interleaved digital Z-gates to excite and measure individual Holstein-Primakoff magnon modes. They present linear-response spectra and decay rates vs temperature and mode index, reporting van Hove singularity enhancement, suppression for corner-localized modes, and a mode-dependent power-law Γ ∝ (ε - ε0)^η with η < 1 for high-energy modes. They also perform nonlinear amplitude-dependent decay measurements and two-mode pump-probe / collider experiments, including a claim that corner modes scatter almost exclusively into each other. The main experimental results are benchmarked against exact simulations up to 31 qubits (median relative errors ~0.7% in frequency and ~7.8% in decay rate) and against MPS simulations for larger systems, where MPS is reported to become inaccurate at elevated temperatures.
Significance. If the central claims hold, this is a substantial advance in quantum simulation of magnetic response functions: it demonstrates mode-resolved control of individual quasiparticles in a 2D interacting spin system, direct measurement of nonlinear magnon scattering channels, and a concrete regime where classical tensor-network methods are claimed to be insufficient. The paper is unusually thorough in its validation: machine-checked small-system comparisons, XEB fidelity analysis, detailed error modeling, and a large supplementary study of classical complexity. The pump-probe and coherent two-mode spectroscopy protocols are novel and could become broadly useful. The main risk is that the 97-qubit mode-resolved conclusions rely on identifying the measured modes with Holstein-Primakoff modes of the pure XY model, whereas the actual device Hamiltonian contains many additional terms; this untested assumption is load-bearing for the corner-mode and mode-resolved decay claims.
major comments (3)
- [§S3 A, Eq. (S25) with Eq. (S5); Figs. 2b,d,f, 3, 5] The mode shapes φ^HP_μ used both to excite and to analyze magnon modes are obtained from the pure XY Hamiltonian (Eq. S25), but the actual device is described by the much richer learned Hamiltonian of Eq. (S5), which includes density-assisted hopping, longer-range XY couplings, and density-density terms, patched from 24-qubit XEB learning. Small-system benchmarks (Nq ≤ 31) validate the full pipeline against exact simulation, but they do not test whether these HP modes remain the true quasiparticle modes of the 97-qubit patched Hamiltonian. If the quadratic part of the full learned Hamiltonian is not diagonalized by the φ^HP modes, the measured χ_μ(t) mixes several true modes. The corner-mode manifold contains only four modes, so even a small overlap leakage could substantially alter the extracted decay rates and the conclusion that corner modes scatter only into each other. I request a d
- [§S1 F, Fig. 3b,c] The central claim of sublinear power-law behavior, η < 1 for high-energy modes, is extracted from a two-parameter fit Γ ∝ (ε - ε0)^η in which ε0 is itself estimated from the same decay-rate data (Sec. S1F, ε0 = -0.57 ± 0.01). The exponent η and ε0 are strongly correlated, so the reported mode-dependent η values and the distinction between η > 1 and η < 1 need a sensitivity analysis: fix ε0 at the MPS ground-state value (-0.554) and at the Hartree-Fock value, and show how η changes for the modes displayed in Fig. 3c. Without this, the comparison to the self-consistent broadening model (which requires γ_ijkl = Γ_i + Γ_j + Γ_k) is not fully convincing, since both the data fit and the model are adjusted. Please also report confidence intervals for η that include the uncertainty in ε0.
- [§S5 A and main-text MPS comparisons (Fig. 2)] The abstract's statement that matrix-product state simulations 'become inaccurate away from these limits' is supported partly by the experimental comparisons in Fig. 2, but the quantitative computational-complexity analysis in Sec. S5 uses a simplified nearest-neighbor XY model (Eq. S45), omits thermalization, and generally adopts assumptions that the authors themselves state are favorable to classical simulation. The resource estimates (e.g., χ ~ 2^25-2^30 for 97 qubits) therefore do not directly apply to the actual learned Hamiltonian used in the experiment. I request that the authors specify whether the MPS curves in Fig. 2 were computed with the full learned Hamiltonian, and, if so, provide convergence data with bond dimension for at least one elevated-temperature case. If only the simplified model was used, the classical-complexity claims should be tempered accordingly.
minor comments (5)
- [Eq. (1)] The text below Eq. (1) defines the energy-conservation delta as δ((ω_k+ω_l-ω_i-ω_j)/γ_ijkl), but Eq. (1) writes δ(Δω_ijkl / γ_ijkl). Please make the notation consistent with the Lorentzian δ_γ(x) defined in Sec. S3E.
- [Fig. 3c] The legend labels for the experimental data, the constant-linewidth theory, and the self-consistent theory are difficult to distinguish in grayscale; please use distinct markers/colors and include error bars on the experimental η values.
- [Fig. 3e] The caption says 'Solid lines are guides to the eye,' but several curves appear to be fits. Please specify exactly which lines are fits and which are just connecting points, and add error bars for the frequency points.
- [Abstract] The phrase 'matrix-product state simulations capture the dynamics well in either small systems or at low temperatures' would benefit from specifying which Hamiltonian was used in those MPS simulations (learned vs nearest-neighbor XY).
- [§S1 G, Eq. (S1)-(S2)] The nonlinear decay model contains a large number of free parameters (a0, f, α1, α2, c, ν, t0, y0) fitted simultaneously. The authors mention that 11 parameters are fitted to 204 points; it would be helpful to report the reduced χ² or a similar goodness-of-fit statistic for representative modes to support the claim that the model is not overfitting.
Circularity Check
No significant circularity: the magnon measurements are direct experimental observables, and the theoretical models are clearly labeled approximate guides rather than fits recycled as predictions.
full rationale
The paper's central results are measured response functions (frequencies, decay rates, nonlinear decay channels) read directly from the device, not derived from the same data by construction. The Holstein-Primakoff mode shapes in Eq. (S25) define the excitation and projection basis, but they are computed from the model Hamiltonian, not fitted to the measured spectra; the measured χμ(t) is a well-defined spin correlation function even if the HP identification is imperfect. The self-consistent broadening model in Eq. (1), with γ_ijkl = Γ_j + Γ_k + Γ_l, is a parameter-free fixed-point calculation: it does not reduce to the observed η by construction, and the paper explicitly reports quantitative disagreements with experiment (Fig. S17), which is the opposite of a fitted input being renamed a prediction. The Hamiltonian learning from Ref. [37] is a methodological self-citation, but it is independently benchmarked here against exact simulation on systems up to 31 qubits (median relative error 0.007 in frequency and 0.078 in decay rate), so the load-bearing validation is not a bare self-citation. The 97-qubit extrapolation rests on an untested assumption that the HP modes remain valid for the patched learned Hamiltonian; this is a validation gap and correctness risk, not a circular equation-level reduction. No quoted step exhibits Eq. X = Eq. Y by construction or a fitted parameter presented as a prediction.
Axiom & Free-Parameter Ledger
free parameters (4)
- ε0 (ground-state energy density in Γ∝(ε−ε0)^η fit) =
-0.57 ± 0.01
- α1, α2, c, ν, t0, y0 (nonlinear decay model) =
mode-dependent values (Figs. 4c,d and S8)
- Hartree–Fock rescaling factor for Δlog ω0 =
1.2
- Lorentzian linewidth γ_ijkl in Eq. (1) =
γ = Σ_i Γ_i (self-consistent)
axioms (5)
- standard math Holstein-Primakoff expansion and bosonic Bogoliubov diagonalization define the magnon modes and mode shapes φμ_i (Eq. S25).
- domain assumption The effective spin Hamiltonian Eq. (S5) with truncation of terms ≲0.001g captures all relevant device dynamics.
- domain assumption Hamiltonian parameters learned at infinite temperature (XEB) transfer unchanged to the lower temperatures used for magnon measurements.
- domain assumption Measured energy density ε maps to a thermal magnon temperature via the mean-field relation Eq. (S41).
- domain assumption The qubit array is an isolated spin-1/2 system; coupler states and higher qubit levels are fully projected out.
read the original abstract
Quantum simulation promises to advance materials discovery by accurately simulating complex states of matter, their microscopic excitations, and macroscopic response functions. The central challenge in resolving the underlying interacting dynamics is to combine high-fidelity evolution with the sophisticated control necessary to manipulate individual quasi-particles in quantum many-body states. Here, we report on high-precision simulation of both linear and non-linear response functions in a 2D XY spin-1/2 magnet using an analog-digital superconducting processor of up to 97 qubits. By interleaving digital gates with analog evolution precisely characterized via Hamiltonian learning, we selectively excite magnons at tunable energy densities. Measuring first the linear magnon response -- a central probe in neutron-scattering experiments -- we extract temperature-dependent spectra and lifetimes. Our results reveal stark variations in magnon decay rates across the Brillouin zone, with enhancement near van Hove singularities and suppression for edge-localized modes. Next, we perform a suite of nonlinear measurements, including the study of self-scattering mechanisms, as well as pump-probe spectroscopy to directly characterize the magnon interactions. While matrix-product state simulations capture the dynamics well in either small systems or at low temperatures, their predictions become inaccurate away from these limits. This work demonstrates precise simulation of the interacting dynamics in quantum magnets, and provides key insights into quasi-particles and their microscopic scattering mechanisms.
Figures
Reference graph
Works this paper leans on
-
[1]
the PEPS bond dimensionχmust be unfeasibly large to capture the evolved states to sufficient fi- delity,
-
[2]
Let us first consider the fidelity
expectation values from low-fidelity PEPS are in- deed inaccurate. Let us first consider the fidelity. Ref. [92] describes a method to approximate the running fidelity of a PEPS during simulation of a quantum circuit as a product of the per-gate fidelities. These fidelities are computed from the magnitude of singular values thrown away, which is precisely...
-
[47]
& Lindner, N
Bairey, E., Arad, I. & Lindner, N. H. Learning a Local Hamiltonian from Local Measurements.Phys. Rev. Lett. 122, 020504 (2019)
2019
-
[48]
& Preskill, J
Huang, H.-Y., Kueng, R. & Preskill, J. Predicting many properties of a quantum system from very few measure- ments.Nat. Phys.16, 1050–1057 (2020)
2020
-
[49]
K., Choi, J., Shaw, A
Mark, D. K., Choi, J., Shaw, A. L., Endres, M. & Choi, S. Benchmarking quantum simulators using ergodic quan- tum dynamics.Phys. Rev. Lett.131, 110601 (2023)
2023
-
[50]
& Cory, D
Wiebe, N., Granade, C., Ferrie, C. & Cory, D. Quantum hamiltonian learning using imperfect quantum resources. Phys. Rev. A89, 042314 (2014)
2014
-
[51]
E., Ferrie, C., Wiebe, N
Granade, C. E., Ferrie, C., Wiebe, N. & Cory, D. G. Robust online Hamiltonian learning.New J. Phys.14, 103013 (2012)
2012
-
[52]
Broughton, M.et al.Tensorflow quantum: A software framework for quantum machine learning.arXiv preprint arXiv:2003.02989(2021)
Pith/arXiv arXiv 2003
-
[53]
Jones, T. & Gacon, J. Efficient calculation of gradients in classical simulations of variational quantum algorithms. arXiv preprint arXiv:2009.02823(2020)
Pith/arXiv arXiv 2009
-
[54]
& Schuster, D
Leung, N., Abdelhafez, M., Koch, J. & Schuster, D. Speedup for quantum optimal control from automatic differentiation based on graphics processing units.Phys. Rev. A95, 042318 (2017)
2017
-
[55]
& Held, K
Wallerberger, M. & Held, K. Trie-based ranking of quan- tum many-body states.Phys. Rev. Res.4, 033238 (2022)
2022
-
[56]
Lin, H. Q. Exact diagonalization of quantum-spin mod- els.Phys. Rev. B42, 6561–6567 (1990)
1990
-
[57]
Dormand, J. R. & Prince, P. J. A family of embedded Runge-Kutta formulae.J. Comput. Appl. Math.6, 19–26 (1980)
1980
-
[58]
Kingma, D. P. & Ba, J. Adam: A method for stochastic optimization. InInternational Conference on Learning Representations (ICLR)(2015)
2015
-
[59]
Manole, T.et al.How much can we learn from quantum random circuit sampling?arXiv preprint arXiv:2510.09919(2025)
arXiv 2025
-
[60]
Phys.14, 595–600 (2018)
Boixo, S.et al.Characterizing quantum supremacy in near-term devices.Nat. Phys.14, 595–600 (2018)
2018
-
[61]
& Primakoff, H
Holstein, T. & Primakoff, H. Field Dependence of the Intrinsic Domain Magnetization of a Ferromagnet.Phys. Rev.58, 1098–1113 (1940)
1940
-
[62]
Anderson, P. W. An approximate quantum theory of the antiferromagnetic ground state.Phys. Rev.86, 694–701 (1952)
1952
-
[63]
& Joannopoulos, J
Gomez-Santos, G. & Joannopoulos, J. D. Application of spin-wave theory to the ground state of XY quantum Hamiltonians.Phys. Rev. B36, 8707–8711 (1987)
1987
-
[64]
Colpa, J. H. P. Diagonalization of the quadratic boson hamiltonian.Phys. A: Stat. Mech. Appl.93, 327–353 (1978)
1978
-
[65]
Jonay, C., Huse, D. A. & Nahum, A. Coarse-grained dynamics of operator and state entanglement.arXiv preprint arXiv:1803.00089(2018)
Pith/arXiv arXiv 2018
-
[66]
Rakovszky, T., von Keyserlingk, C. W. & Pollmann, F. Entanglement growth after inhomogenous quenches. Phys. Rev. B.100, 125139 (2019)
2019
-
[67]
M., Verstraete, F
Schuch, N., Wolf, M. M., Verstraete, F. & Cirac, J. I. En- tropy scaling and simulability by matrix product states. Phys. Rev. Lett.100, 030504 (2008)
2008
-
[68]
& White, S
Barthel, T., Schollw¨ ock, U. & White, S. R. Spectral func- tions in one-dimensional quantum systems at finite tem- perature using the density matrix renormalization group. Phys. Rev. B79, 245101 (2009)
2009
-
[69]
Yi-Thomas, S., Ware, B., Sau, J. D. & White, C. D. Comparing numerical methods for hydrodynamics in a one-dimensional lattice spin model.Phys. Rev. B110, 134308 (2024)
2024
-
[70]
White, C. D. Effective dissipation rate in a Liouvillian- graph picture of high-temperature quantum hydrody- namics.Phys. Rev. B.107, 094311 (2023)
2023
-
[71]
& Rakovszky, T
von Keyserlingk, C., Pollmann, F. & Rakovszky, T. Op- erator backflow and the classical simulation of quantum transport.Phys. Rev. B.105, 245101 (2022)
2022
-
[72]
Rakovszky, T., von Keyserlingk, C. W. & Pollmann, F. Dissipation-assisted operator evolution method for capturing hydrodynamic transport.Phys. Rev. B.105, 075131 (2022)
2022
-
[73]
E., Cao, X., Avdoshkin, A., Scaffidi, T
Parker, D. E., Cao, X., Avdoshkin, A., Scaffidi, T. & Alt- man, E. A universal operator growth hypothesis.Phys. Rev. X.9, 041017 (2019)
2019
-
[74]
Artiaco, C., Fleckenstein, C., Aceituno Ch´ avez, D., Kvorning, T. K. & Bardarson, J. H. Efficient large-scale many-body quantum dynamics via local-information time evolution.PRX Quantum5, 020352 (2024)
2024
-
[75]
P., Trauzettel, B., Klein Kvorning, T., Bar- darson, J
Bauer, N. P., Trauzettel, B., Klein Kvorning, T., Bar- darson, J. H. & Artiaco, C. Local information flow in quantum quench dynamics.Phys. Rev. A112, 022221 (2025)
2025
-
[76]
qsim (2025)
Quantum AI team and collaborators. qsim (2025)
2025
-
[77]
Haegeman, J.et al.Time-Dependent Variational Princi- ple for Quantum Lattices.Phys. Rev. Lett.107, 070601 (2011)
2011
-
[78]
Paeckel, S.et al.Time-evolution methods for matrix- product states.Ann. Phys. (N. Y.)411, 167998 (2019)
2019
-
[79]
Cataldi, G.et al.Hilbert curve vs Hilbert space: ex- ploiting fractal 2D covering to increase tensor network efficiency.Quantum5, 556 (2021)
2021
-
[80]
& Bar Lev, Y
Kloss, B., Reichman, D. & Bar Lev, Y. Studying dynam- ics in two-dimensional quantum lattices using tree tensor network states.SciPost Phys.9, 070 (2020)
2020
-
[81]
Vovrosh, J.et al.Simulating dynamics of the two- dimensional transverse-field ising model: A comparative study of large-scale classical numerics.Phys. Rev. Res. 8, 023311 (2026)
2026
-
[82]
& White, S
Yang, M. & White, S. R. Time-dependent variational principle with ancillary krylov subspace.Phys. Rev. B 102, 094315 (2020)
2020
-
[83]
Singh, S., Pfeifer, R. N. C. & Vidal, G. Tensor network states and algorithms in the presence of a global U(1) symmetry.Phys. Rev. B83, 115125 (2011)
2011
-
[84]
Hubig, C., McCulloch, I. P. & Schollw¨ ock, U. Generic construction of efficient matrix product operators.Phys. Rev. B95, 035129 (2017)
2017
-
[85]
Thompson, A. P., Soeteman, A., Cade, C. & Niesen, I. Non-zero noise extrapolation: accurately simulat- ing noisy quantum circuits with tensor networks.arXiv preprint arXiv:2501.13237(2025)
arXiv 2025
-
[86]
Zhou, Y., Stoudenmire, E. M. & Waintal, X. What limits the simulation of quantum computers?Phys. Rev. X.10, 041038 (2020). 42
2020
-
[87]
Ayral, T.et al.Density-matrix renormalization group algorithm for simulating quantum circuits with a finite fidelity.PRX Quantum4, 020304 (2023)
2023
-
[88]
Mandr` a, S.et al.A heuristic for matrix product state simulation of out-of-equilibrium dynamics of two- dimensional transverse-field ising models.arXiv preprint arXiv:2511.23438(2025)
arXiv 2025
-
[89]
& Fishman, M
Tindall, J. & Fishman, M. Gauging tensor networks with belief propagation.SciPost Physics15, 222 (2023)
2023
-
[90]
& Chan, G
Evenbly, G., Pancotti, N., Milsted, A., Gray, J. & Chan, G. K.-L. Loop series expansions for tensor networks. Phys. Rev. Res.8, 013245 (2026)
2026
-
[91]
Gray, J.et al.Tensor network loop cluster expan- sions for quantum many-body problems.arXiv preprint arXiv:2510.05647(2025)
Pith/arXiv arXiv 2025
-
[92]
Rudolph, M. S. & Tindall, J. Simulating and sampling from quantum circuits with 2D tensor networks.arXiv preprint arXiv:2507.11424(2025)
arXiv 2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.