REVIEW 3 major objections 5 minor 47 references
This paper shows that neural quantum states, propagated by time-dependent variational methods, reproduce 1H NMR spectra of small molecules with spectral error below 10^-3, and that an exact interaction-frame transformation reduces the requi
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-04 01:03 UTC pith:LNGNMZLL
load-bearing objection A solid proof-of-principle that NQS dynamics can match exact NMR spectra for small spins, with an exact interaction-frame trick that empirically cuts integration steps 8–40×; the theoretical justification is sloppier than the numerics, and the 14-spin demo is a benign disconnected case. the 3 major comments →
Neural Quantum States for Nuclear Magnetic Resonance Spectroscopy
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 discovery is that the steep step-count growth in NQS-NMR simulation is numerical stiffness from the chemical-shift term H0 = Σ ω_i I_z,i, not a failure of the neural-network ansatz. Moving to the interaction frame of H0 — the Dirac picture, with |Ψ(t)⟩ = U0†(t)|Ψ_I(t)⟩ — removes the large one-body term from the generator exactly, leaving only coupling terms dressed by phases oscillating at chemical-shift differences. The transformation is exact and leaves the spectrum unchanged: for the RBM ansatz it amounts to a time-dependent drift of the visible biases, a_i(t) = a_i(0) − i(ω_i/2)t, and peak positions are restored analytically. In practice, this cuts the number of integration s
What carries the argument
The load-bearing object is the interaction-frame generator of Eq. (4), H_I(t) = 2π Σ_{i<j} J_{ij} [I_z,i I_z,j + (1/2)(e^{i(ω_i−ω_j)t} I_+^{(i)} I_-^{(j)} + h.c.)], obtained by moving into the rotating frame of the chemical-shift term H0 = Σ ω_i I_z,i. This transformation is exact: the spectrum is unchanged, and for the RBM the chemical-shift evolution becomes an analytic drift of the visible biases, a_i(t) = a_i(0) − i(ω_i/2)t. The generator removes the large one-body term that forces tiny integration steps, leaving a coupling-scale generator whose time dependence the state follows slowly. The RBM wavefunction of Eq. (6) supplies the variational manifold, and TDVP or p-tVMC supplies the pro
Load-bearing premise
The advantage rests on the assumption that the only numerically stiff part of the Hamiltonian is the one-body chemical-shift term, so that after applying it exactly the remaining time-dependent coupling generator can be integrated with large steps; if the oscillating phases between shifts still force tiny steps, or if strongly entangled states demand exponentially many hidden units, the 14-spin feasibility does not extend to realistic molecules.
What would settle it
Run interaction-frame TDVP on a synthetic 12-proton network with all spins coupled in one connected cluster (no isolated blocks), using α = 4, and compare the resulting spectrum to exact Krylov propagation at T = 1 s: if the step count needed for MSE < 10^-3 grows at the same steep rate as in the rotating frame, or if α must scale steeply with N to hold the error, the central claim fails.
If this is right
- At the benchmarked sizes, accuracy is limited by integration and optimization error, not by the RBM's expressive capacity, shifting attention to integrators and update schemes rather than network architecture.
- The interaction-frame transformation is exact and spectrum-preserving, so any future NQS or non-NQS propagator can adopt it as a cheap preconditioner for the chemical-shift term.
- The 14-spin sucrose run is a feasibility point: it works because the coupling network splits into 7-, 5-, and 2-spin clusters, so the hardest part is the 7-spin block.
- Spectral MSE in the 10^-4 to 10^-3 range is below typical experimental linewidths, so the accuracy is adequate for practical comparison with experiment if it persists at larger N.
- For p-tVMC, the per-step infidelity provides an adaptive convergence diagnostic that TDVP lacks, potentially making optimization budgets self-tuning at larger scale.
Where Pith is reading between the lines
- The interaction-frame idea is independent of the neural-network ansatz: tensor-network and restricted state-space NMR simulations that suffer from the same chemical-shift stiffness could absorb the one-body term analytically as well.
- The sucrose result is a best-case large-system point; a molecule with 10–20 protons in a single connected coupling network is the natural next stress test, and volume-law entanglement could make the RBM fail exactly where the interaction frame has done its job.
- Because the visible-bias drift is exact, the same trick should transfer directly to p-tVMC, combining the frame's step reduction with p-tVMC's error diagnostic.
- One can test the entanglement-growth hypothesis directly by computing bipartite von Neumann entropy along the trajectory for a synthetic 10-spin connected network; if the worst-cut entropy approaches volume-law scaling, expect the RBM hidden-unit requirement to grow quickly.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper benchmarks Restricted Boltzmann Machine (RBM) neural quantum states (NQS) for simulating nuclear magnetic resonance spectra of small coupled proton spin systems. Using time-dependent variational principle (TDVP) and projected time-dependent variational Monte Carlo (p-tVMC), the authors compute 1H NMR spectra for four molecules with 2–5 spins and compare against exact diagonalization. They report spectral MSEs below 10^-3 for TDVP at optimized parameters, and find that the dominant cost is the number of integration steps, which grows with the chemical-shift bandwidth. An interaction-frame transformation is introduced that removes the one-body chemical-shift term exactly; the authors claim this reduces the required integration steps roughly eightfold (3 spins) to an order of magnitude or more (4–5 spins), and enables a Monte Carlo simulation of sucrose (14 spins). The paper also analyzes entanglement entropy and discusses limitations, including disconnected coupling clusters in sucrose.
Significance. If the results hold, the paper provides a credible proof-of-concept that NQS can reproduce exact NMR spectra at accuracy sufficient for practical comparison, and it identifies a potentially useful numerical trick (interaction-frame propagation) for reducing integration cost. The benchmarks are non-circular (compared with exact spectra), include multi-seed statistics, and the entanglement analysis is a valuable caveat. The claim of a general speed-up and the extension to a 14-spin molecule are, however, not fully supported by the theoretical argumentation and the specific choice of a factorizable test case. The paper is careful and honest in many of its stated limitations, which is a strength.
major comments (3)
- [§II.B, Eq. (4)] The claim that the interaction frame allows larger steps because 'the rate at which the state changes is bounded by the norm of the generator, not by how rapidly its phases oscillate' is not a correct criterion for step-size selection in a second-order integrator for a time-dependent generator. Local error of Heun's method involves derivatives of H_I(t), so the oscillatory phases at Δω must be resolved. For the reported 4-spin run (N_t≈200, h=5×10^-3, stiffness 36.9), the phase increment Δω h is O(1–10), which should make the trapezoidal approximation of the flip-flop terms unreliable. The paper needs either a local error analysis for this time-dependent Hamiltonian or an empirical convergence study (e.g., comparing with a smaller h or a higher-order integrator) to support the step-reduction claim.
- [§IV.A, Table I vs. p-tVMC paragraph] The TDVP/p-tVMC comparison is not apples-to-apples: TDVP used molecule-specific optimal α (Table I) while p-tVMC used α=4 and n_iter=50 (or 75) throughout. The reported p-tVMC MSEs (2.6×10^-3 to 3.7×10^-2) are substantially worse than TDVP's. If the paper aims to demonstrate that both methods 'reproduce the same features', the quantitative comparison is incomplete. An accuracy-matched comparison (e.g., tuning α and n_iter per molecule, or reporting MSE at matching wall time) would strengthen the claim.
- [§IV.B, Fig. 7] The sucrose (N=14) feasibility point is a weak demonstration of scalability because the J-coupling network splits into three disconnected clusters (7, 5, and 2 spins). As Fig. 7 shows, the state remains a product over clusters, so the effective problem size is 7 spins. The abstract's statement that the method 'enables a 14-spin molecule (sucrose) to be accurately propagated' is misleading without this caveat. A connected 14-spin system with genuine entanglement would be a more convincing test of the NQS approach at larger N.
minor comments (5)
- [§IV.B] The sentence 'The successful simulation of the 14-proton molecule points works in favor of NQ.' is malformed; it should read 'points work in favor of NQS' or similar.
- [§III.C] The dense propagator for p-tVMC is precomputed and cached, which limits the method to N≤12; this is mentioned but could be emphasized more clearly in the abstract and conclusions.
- [§II.B] The 'stiffness' ratio max_{ij}|ω_i−ω_j|/(2πJ_ij) is used to quantify the numerical difficulty, but for a multi-spin system the relevant frequency is the spectral gap of H0, not a single pair ratio. Consider clarifying or using a more conventional definition.
- [§I/References] The paper would benefit from a brief comparison with existing restricted state-space methods (e.g., Spinach) to contextualize the claimed advantages; the discussion section mentions this as future work, but not even a qualitative comparison is given.
- [§IV.A] The p-tVMC results are reported only for α=4 and n_iter=50 (and one n_iter=75 point for the 5-spin system). Since the paper argues p-tVMC is competitive, a more thorough parameter study (e.g., n_iter sweep for each molecule) would be useful.
Circularity Check
No significant circularity: spectra are benchmarked against independently computed exact references; the interaction-frame transformation is exact and standard, and the speed-up claim is empirical.
full rationale
The paper's central claim is that NQS (RBM) with TDVP and p-tVMC reproduces 1H NMR spectra for 2-5 spins, benchmarked against exact diagonalization / sparse Krylov references. This is an external, independent benchmark: the Hamiltonian parameters come from experimental databases, the reference spectra are computed by dense exponentiation or Krylov propagation, and the reported MSEs measure the variational/integration error against those references. Choosing the timestep or hidden-unit ratio by minimizing MSE is a standard numerical benchmark procedure, not fitting a parameter to a target and then calling it a prediction. The interaction-frame transformation of Eq. (4) is presented explicitly as the standard Dirac picture, with |Psi(t)> = U0^dag(t)|Psi_I(t)> reproducing the rotating-frame state exactly; the paper does not redefine any quantity in terms of the claimed result. The reported step-count reduction is an empirical observation from the same exact-reference comparison, not a conclusion forced by definition. The only self-citations are [23] and [27], co-authored by B. Xing, but they are background references for p-tVMC variants and ground-state search; the core p-tVMC method is attributed to [22] (Sinibaldi et al.), and no load-bearing conclusion rests on a self-citation or on an imported uniqueness theorem. The skeptical concern about the Heun integrator needing to resolve oscillatory phases in H_I(t) is a legitimate numerical-analysis and correctness concern, but it is not circularity: it challenges whether the stated error analysis is valid, not whether the derivation reduces to its inputs. The paper itself includes caveats (second-order integrator, sucrose as a feasibility point, open questions about strongly entangled systems) that further indicate the claims are presented as empirical benchmarks rather than as definitionally forced results.
Axiom & Free-Parameter Ledger
free parameters (4)
- hidden unit ratio α =
α=4 for p-tVMC; α=10 (N=2,3,5), α=1 (N=4) for TDVP; α=2 for sucrose
- timestep Δt =
0.005 s (N=2), 0.002 s (N=3), 2.5×10^-4 s (N=4), 1.0×10^-4 s (N=5)
- QGT regularization ϵ =
0.01
- apodization decay λ =
3.0 s^-1
axioms (4)
- domain assumption Variational Monte Carlo with 4096 samples estimates QGT/forces and infidelity gradients with small enough bias for the central accuracy claim.
- domain assumption The effective NMR Hamiltonian (1) is the correct model for the simulated spectra.
- domain assumption Spectral MSE over a normalized region is a meaningful error metric.
- domain assumption The RBM can represent the time-evolved states with the chosen α (expressivity assumption).
read the original abstract
Predicting a nuclear magnetic resonance (NMR) spectrum from first principles requires propagating a quantum state of dimension $2^N$ for $N$ coupled spins, which becomes intractable beyond larger $N$. We benchmark Neural Quantum States (NQS), a class of variational quantum states expressed as an artificial neural network, as an alternative representation for this problem. Using two propagation methods, the Time-Dependent Variational Principle (TDVP) and projected time-dependent Variational Monte Carlo (p-tVMC), we compute the $^1$H spectra of four ($2 \to 5$ spins) experimentally parameterized molecules. TDVP reproduces all line positions and intensities with average spectral mean squared errors of $<10^{-3}$; p-tVMC reproduces the same features, with accuracy determined by its per-step optimization parameters. One dominant obstacle to larger systems is the steep growth of the number of integration steps with spectral bandwidth, which can be removed by propagating in the interaction frame of the chemical-shifted Hamiltonian. Retaining the same accuracy, this reduces the number of integration steps roughly eightfold for the 3-spin system and by at least an order of magnitude for the 4- and 5-spin systems, and it enables a 14-spin molecule (sucrose) to be accurately propagated via Monte Carlo sampling.
Figures
Reference graph
Works this paper leans on
-
[1]
R. R. Ernst, G. Bodenhausen, and A. Wokaun,Prin- ciples of Nuclear Magnetic Resonance in One and Two Dimensions(Clarendon Press, Oxford, 1987)
1987
-
[2]
M. H. Levitt,Spin Dynamics: Basics of Nuclear Mag- netic Resonance, 2nd ed. (John Wiley & Sons, 2008)
2008
-
[3]
Helgaker, M
T. Helgaker, M. Jaszu´ nski, and K. Ruud, Ab initio meth- ods for the calculation of NMR shielding and indirect spin–spin coupling constants, Chemical Reviews99, 293 (1999)
1999
-
[4]
M. W. Lodewyk, M. R. Sieber, and D. J. Tantillo, Com- putational prediction of 1H and 13C chemical shifts: A useful tool for natural product, mechanistic, and syn- thetic organic chemistry, Chemical Reviews112, 1839 (2012)
2012
-
[5]
Wolinski, J
K. Wolinski, J. F. Hinton, and P. Pulay, Efficient im- plementation of the gauge-independent atomic orbital method for NMR chemical shift calculations, Journal of the American Chemical Society112, 8251 (1990)
1990
-
[6]
M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone, G. A. Petersson, H. Nakatsuji,et al., Gaussian 16, revi- sion C.01 (2016)
2016
-
[7]
Neese, F
F. Neese, F. Wennmohs, U. Becker, and C. Riplinger, The ORCA quantum chemistry program package, The Journal of Chemical Physics152, 224108 (2020)
2020
-
[8]
D. V. Savostyanov, S. V. Dolgov, J. M. Werner, and I. Kuprov, Exact NMR simulation of protein-size spin systems using tensor train formalism, Physical Review B 90, 085139 (2014)
2014
-
[9]
Kuprov, N
I. Kuprov, N. Wagner-Rundell, and P. J. Hore, Polyno- mially scaling spin dynamics simulation algorithm based on adaptive state-space restriction, Journal of Magnetic Resonance189, 241 (2007)
2007
-
[10]
H. J. Hogben, M. Krzystyniak, G. T. P. Charnock, P. J. Hore, and I. Kuprov, Spinach – a software library for sim- ulation of spin dynamics in large spin systems, Journal 12 of Magnetic Resonance208, 179 (2011)
2011
-
[11]
J. A. Jones, Controlling NMR spin systems for quantum computation, Progress in Nuclear Magnetic Resonance Spectroscopy140–141, 49 (2024)
2024
-
[12]
Lloyd, Universal quantum simulators, Science273, 1073 (1996)
S. Lloyd, Universal quantum simulators, Science273, 1073 (1996)
1996
-
[13]
Preskill, Quantum computing in the NISQ era and beyond, Quantum2, 79 (2018)
J. Preskill, Quantum computing in the NISQ era and beyond, Quantum2, 79 (2018)
2018
-
[14]
worst cut
In this case, the coupling structure can further sim- plify molecules. More broadly, the appeal of NQS for larger systems is that, unlike tensor networks, they are not bounded by connectivity or entanglement across a cut [43]. They can remain efficient as that entanglement grows. Two caveats bound this comparison. First, both frames used the same second-o...
-
[15]
A. Khedri, P. Stadler, K. Bark, M. Lodi, R. Reiner, N. Vogt, M. Marthaler, and J. Lepp¨ akangas, The im- pact of noise on the simulation of NMR spectroscopy on NISQ devices, arXiv preprint arXiv:2404.18903 10.48550/arXiv.2404.18903 (2024)
-
[16]
K. R. Fratuset al., Can a quantum computer simulate nuclear magnetic resonance spectra better than a classical one?, arXiv preprint arXiv:2508.06448 10.48550/arXiv.2508.06448 (2025)
-
[17]
Or´ us, A practical introduction to tensor networks: Matrix product states and projected entangled pair states, Annals of Physics349, 117 (2014)
R. Or´ us, A practical introduction to tensor networks: Matrix product states and projected entangled pair states, Annals of Physics349, 117 (2014)
2014
-
[18]
S. R. White and A. E. Feiguin, Real-time evolution us- ing the density matrix renormalization group, Physical Review Letters93, 076401 (2004)
2004
-
[19]
Vidal, Efficient simulation of one-dimensional quan- tum many-body systems, Physical Review Letters93, 040502 (2004)
G. Vidal, Efficient simulation of one-dimensional quan- tum many-body systems, Physical Review Letters93, 040502 (2004)
2004
-
[20]
Carleo and M
G. Carleo and M. Troyer, Solving the quantum many- body problem with artificial neural networks, Science 355, 602 (2017), seminal work introducing Neural Quan- tum States
2017
-
[21]
Schmitt and M
M. Schmitt and M. Heyl, Quantum many-body dynamics in two dimensions with artificial neural networks, Physi- cal Review Letters125, 100503 (2020)
2020
-
[22]
X. Yuan, S. Endo, Q. Zhao, Y. Li, and S. C. Benjamin, Theory of variational quantum simulation, Quantum3, 191 (2019)
2019
-
[23]
Sinibaldi, C
A. Sinibaldi, C. Giuliani, G. Carleo, and F. Vicentini, Un- biasing time-dependent variational monte carlo by pro- jected quantum evolution, Quantum7, 1131 (2023)
2023
-
[24]
Zhang, B
W. Zhang, B. Xing, X. Xu, and D. Poletti, Paths to- wards time evolution with larger neural-network quan- tum states, Computer Physics Communications312, 109577 (2025)
2025
-
[25]
Gao and L.-M
X. Gao and L.-M. Duan, Efficient representation of quan- tum many-body states with deep neural networks, Na- ture Communications8, 662 (2017)
2017
-
[26]
Hibat-Allah, M
M. Hibat-Allah, M. Ganahl, L. E. Hayward, R. G. Melko, and J. Carrasquilla, Recurrent neural network wave func- tions, Physical Review Research2, 023358 (2020)
2020
-
[27]
D. Pfau, J. S. Spencer, A. G. Matthews, and W. M. C. Foulkes, Ab initio solution of the many-electron Schr¨ odinger equation with deep neural networks, Physi- cal Review Research2, 033429 (2020)
2020
-
[28]
Zhang, X
W. Zhang, X. Xu, Z. Wu, V. Balachandran, and D. Po- letti, Ground state search by local and sequential updates of neural network quantum states, Phys. Rev. B107, 165149 (2023)
2023
-
[29]
Dashti, W
H. Dashti, W. M. Westler, M. Tonelli, J. R. Wedell, J. L. Markley, and H. R. Eghbalnia, Spin system modeling of nuclear magnetic resonance spectra for applications in metabolomics and small molecule screening, Analytical Chemistry89, 12201 (2017), gISSMO database
2017
-
[30]
HQS Quantum Simulations GmbH, HQS quantum sim- ulations: Struqture and cheshire NMR parameter database (2023), accessed: 2024
2023
-
[31]
Bloch, Nuclear induction, Physical Review70, 460 (1946)
F. Bloch, Nuclear induction, Physical Review70, 460 (1946)
1946
-
[32]
E. L. Hahn, Spin echoes, Physical Review80, 580 (1950)
1950
-
[33]
Abragam,The Principles of Nuclear Magnetism(Ox- ford University Press, Oxford, 1961)
A. Abragam,The Principles of Nuclear Magnetism(Ox- ford University Press, Oxford, 1961)
1961
-
[34]
C. P. Slichter,Principles of Magnetic Resonance, 3rd ed. (Springer-Verlag, Berlin, 1990)
1990
-
[35]
Kubo, Statistical-mechanical theory of irreversible processes
R. Kubo, Statistical-mechanical theory of irreversible processes. I. general theory and simple applications to magnetic and conduction problems, Journal of the Phys- ical Society of Japan12, 570 (1957)
1957
-
[36]
Becca and S
F. Becca and S. Sorella,Quantum Monte Carlo Ap- proaches for Correlated Systems(Cambridge University Press, Cambridge, 2017)
2017
-
[37]
A. D. McLachlan, A variational solution of the time- dependent Schr¨ odinger equation, Molecular Physics8, 39 (1964)
1964
-
[38]
Haegeman, J
J. Haegeman, J. I. Cirac, T. J. Osborne, I. Piˇ zorn, H. Ver- schelde, and F. Verstraete, Time-dependent variational principle for quantum lattices, Physical Review Letters 107, 070601 (2011)
2011
-
[39]
Sorella, Green function Monte Carlo with stochastic reconfiguration, Physical Review Letters80, 4558 (1998)
S. Sorella, Green function Monte Carlo with stochastic reconfiguration, Physical Review Letters80, 4558 (1998)
1998
-
[40]
Stokes, J
J. Stokes, J. Izaac, N. Killoran, and G. Carleo, Quantum natural gradient, Quantum4, 269 (2020)
2020
-
[41]
Hackl, T
L. Hackl, T. Guaita, T. Shi, J. Haegeman, E. Demler, and J. I. Cirac, Geometry of variational methods: dynamics of closed quantum systems, SciPost Physics9, 48 (2020)
2020
-
[42]
D. P. Kingma and J. Ba, Adam: A method for stochas- tic optimization, arXiv preprint arXiv:1412.6980 (2015), presented at ICLR 2015
Pith/arXiv arXiv 2015
-
[43]
Vicentini, D
F. Vicentini, D. Hofmann, A. Szab´ o, D. Wu, C. Roth, C. Giuliani, G. Pescia, J. Nys, V. Vargas-Calder´ on, N. Astrakhantsev, and G. Carleo, NetKet 3: Machine learning toolbox for many-body quantum systems, Sci- Post Physics Codebases , 7 (2022)
2022
-
[44]
D.-L. Deng, X. Li, and S. Das Sarma, Quantum entan- glement in neural network states, Physical Review X7, 021021 (2017)
2017
-
[45]
Chen and M
A. Chen and M. Heyl, Empowering deep neural quantum states through efficient optimization, Nature Physics20, 1476 (2024)
2024
-
[46]
L. L. Viteritti, R. Rende, and F. Becca, Transformer vari- ational wave functions for frustrated quantum spin sys- tems, Physical Review Letters130, 236401 (2023)
2023
-
[47]
Breuer and F
H.-P. Breuer and F. Petruccione,The Theory of Open Quantum Systems(Oxford University Press, 2002)
2002
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.