pith. sign in

arxiv: 2604.24603 · v1 · submitted 2026-04-27 · 🪐 quant-ph · nlin.CD

Quantum vs. Classical Spin: A Comparative Study of Dipolar Spin Dynamics and the Onset of Chaos

Victor Henner , Alexander Nepomnyashchy , Tatyana Belozerova This is my paper

Pith reviewed 2026-05-08 04:11 UTC · model grok-4.3

classification 🪐 quant-ph nlin.CD
keywords quantum spin dynamicsclassical spin simulationsdipolar interactionsfree induction decayquantum-classical comparisonspin chaos
0
0 comments X

The pith

Quantum and classical spins produce significantly different free induction decay despite similar overall patterns.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper compares the evolution of dipole-coupled spins by directly solving the Schrödinger equation for quantum spins against classical spin simulations. Using free induction decay as the benchmark, it shows that the two approaches yield qualitatively similar patterns but with clear discrepancies at both short and long timescales. These differences are attributed to fundamental distinctions between the quantum and classical descriptions. A reader would care because classical models are routinely employed in microscopic spin dynamics simulations, raising questions about their accuracy for dipolar systems.

Core claim

Although classical spins have long been used in microscopic spin dynamics simulations, their results differ significantly from those of quantum spins in a dipole-coupled system. When benchmarked against free induction decay, the overall patterns are qualitatively similar yet significant discrepancies appear at both short and long timescales, which can be traced to fundamental distinctions in the two descriptions.

What carries the argument

Free induction decay as the comparative metric between direct Schrödinger-equation solutions for quantum spins and classical spin simulations in a dipole-coupled system.

If this is right

  • Classical simulations cannot be assumed to reproduce quantum spin dynamics reliably at all timescales in dipolar systems.
  • The onset of chaos in spin systems may be described differently under quantum versus classical rules.
  • Precise predictions for dipolar spin evolution require quantum treatments rather than classical approximations.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Similar discrepancies could appear in other observables beyond free induction decay, such as spin correlation functions.
  • The result suggests testing whether quantum-classical differences persist or diminish as system size increases toward the thermodynamic limit.
  • This comparison framework could be applied to other interaction types, like exchange or Zeeman terms, to map where classical approximations break down.

Load-bearing premise

The chosen system size, initial conditions, and free induction decay metric provide a fair and representative comparison between the quantum and classical models without hidden approximations that favor one over the other.

What would settle it

Repeating the comparison for substantially larger numbers of spins or different initial states and finding that the short- and long-time discrepancies in free induction decay vanish would falsify the claim of significant fundamental differences.

Figures

Figures reproduced from arXiv: 2604.24603 by Alexander Nepomnyashchy, Tatyana Belozerova, Victor Henner.

Figure 1
Figure 1. Figure 1: Results of simulations of FID with 12 quantum spins. Left panel view at source ↗
Figure 2
Figure 2. Figure 2: The envelope of () x et in view at source ↗
Figure 3
Figure 3. Figure 3: Fig.3. panel a view at source ↗
Figure 4
Figure 4. Figure 4: Graph of () x et . 5 5 5, Nspins =   0.01 d p = , (0) 0.999999999 x e = . Two panels represent two runs with identical parameters. The dimensionless time is 0 tt = . The irregular behavior of () x et at large t, which changes significantly due to tiny changes of initial conditions, is in sharp contrast to the regular and reproducible behavior of the quantum system. To confirm the presence of chaos in t… view at source ↗
Figure 5
Figure 5. Figure 5: Graph of Lt() . Panel a: 5 N N N x y z = = = , (0) 0.9999 x e = , 0.01 d p = . L(30000) 0.009889 = . Panel b: 6, 1, 1 N N N x y z = = = , (0) 0.5 x e = , 0.01 d p = . L(70000) = 0.002927 The same value L(70000) = 0.002927 and practically identical graph is obtained for 12 spins, 12, 1, 1 N N N x y z = = = view at source ↗
Figure 6
Figure 6. Figure 6: ln ( ) (Lt ) versus ln( )t . (0) 0.5 x e = , 0.01 d p = . Panel a: 2, 1, 1 N N N x y z = = = , L(70000) 0.0000587 = . Panel b: 3, 1, 1 N N N x y z = = = , L(70000) 0.00006837 = . We can conclude that for N  3 spins the classical spins dynamics (for the value 0.01 d p = ) are chaotic. We also have analyzed the dependence of the Lyapunov exponent L on the dipole interactions parameter d p . For small value… view at source ↗
Figure 7
Figure 7. Figure 7: 3 N = 7 , 0.01 d p = , (0) 0.5 x e = . Panel a: linear initial distribution, panel b: random initial distribution. Solid line on panel b represents function (12) with the parameters A = 0.5 , a = 0.0009, b = 0.045. Very close results are obtained for 3 N = 5 spins and (0) 0.5 x e = with the same, as in view at source ↗
Figure 8
Figure 8. Figure 8: N =   7 7 7, 0.01 d p = , (0) 0.95 x e = . Panel a: Graph () x et . Panel b. Fourier transform of () x et . The Fourier transform in Figure 8b has the same spectral width as that of the quantum spins in view at source ↗
Figure 9
Figure 9. Figure 9: Graph () x et , 0.01 d p = , (0) 0.95 x e = . Panel a: N =   555 . Panel b: N =   333 view at source ↗
Figure 10
Figure 10. Figure 10: In the special case EM == 1, 0 , equation (B14) has three kinds of solutions, a stationary solution 1 c = 0 and aperiodic solutions 1 2,3 0 c ( ) 2 / 3cosh ( 2( ))    − =− . The stationary solution 1 c = 0 corresponds to initial conditions 12    == /2, = 0 , i.e., 12 xx (0) (0) 1 == . Solutions 2,3 c describe a homoclinic trajectory: under an arbitrary small initial disturbance, the system leave th… view at source ↗
Figure 10
Figure 10. Figure 10: 1, 1, 2 N N N x y z = = = , 0.01 d p = , (0) 0.9999 x e = References [1] S. J. Knak Jensen, O. Platz, Free-Induction-Decay Shapes in a Dipolar-Coupled Rigid Lattice of Infinite Nuclear Spins, Phys. Rev. B 7 (1973) 31. [2] S. J. Jensen and E. K. Hansen, Dynamics of classical spins with dipolar coupling in a rigid lattice at high temperatures, Phys. Rev. B 13 (1976) 1903. [3] A. A. Lundin, V. E. Zobov, Simu… view at source ↗
read the original abstract

We investigate the spin dynamics of a dipole-coupled system by comparing a direct solution of the Schrodinger equation for quantum spins with simulations of classical spins. Although classical spins have long been used in microscopic spin dynamics simulations, we demonstrate that their results differ significantly from those of quantum spins. Using Free Induction Decay as a benchmark, we find that while the overall patterns are qualitatively similar, significant discrepancies emerge at both short and long timescales. We trace these differences to fundamental distinctions in the two descriptions.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript compares quantum spin dynamics, obtained by direct solution of the Schrödinger equation, with classical spin simulations for a dipole-coupled system. It uses Free Induction Decay (FID) as the primary benchmark and reports that, while overall patterns are qualitatively similar, significant quantitative discrepancies appear at both short and long timescales. These differences are attributed to fundamental distinctions between the quantum (operator) and classical (vector) descriptions.

Significance. If the discrepancies are shown to arise from intrinsic quantum-classical differences rather than finite-size artifacts, the work would usefully caution against routine use of classical spins in microscopic dipolar simulations, especially for small clusters or when entanglement matters. The direct side-by-side comparison is a strength, but its interpretive weight depends on controls that are not yet demonstrated.

major comments (2)
  1. [Abstract and Results (FID comparison)] The central claim (abstract and concluding discussion) that observed FID discrepancies trace to 'fundamental distinctions' is load-bearing yet untested against finite-N effects. Quantum exact diagonalization is feasible only for N ≲ 8–10; classical vector dynamics at the same small N need not reproduce the large-N limit in which classical approximations are normally justified. No N-scaling study, alternative initial states, or comparison to larger classical ensembles is described, leaving open the possibility that the reported differences are dominated by finite-size artifacts rather than intrinsic quantum features.
  2. [Methods and Results sections] No quantitative error analysis, specific Hamiltonian parameters (e.g., dipolar coupling strengths, lattice geometry), or explicit initial-state definitions are provided to allow independent assessment of whether the short- and long-time discrepancies are robust or sensitive to choices. The absence of these details makes it impossible to judge whether the 'significant discrepancies' exceed numerical or statistical uncertainties.
minor comments (2)
  1. [Abstract] The abstract contains no equations, system sizes, or numerical values, forcing the reader to consult the full text for even basic context; adding a concise statement of N, initial polarization, and the FID observable would improve accessibility.
  2. [Methods] Notation for the classical equations of motion and the precise definition of the FID metric should be stated explicitly (e.g., as an equation) rather than left implicit.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the careful review and for identifying key issues regarding finite-size effects and methodological details. We respond to each major comment below and indicate the changes we will incorporate in the revised manuscript.

read point-by-point responses

  1. Referee: [Abstract and Results (FID comparison)] The central claim (abstract and concluding discussion) that observed FID discrepancies trace to 'fundamental distinctions' is load-bearing yet untested against finite-N effects. Quantum exact diagonalization is feasible only for N ≲ 8–10; classical vector dynamics at the same small N need not reproduce the large-N limit in which classical approximations are normally justified. No N-scaling study, alternative initial states, or comparison to larger classical ensembles is described, leaving open the possibility that the reported differences are dominated by finite-size artifacts rather than intrinsic quantum features.

    Authors: We agree that finite-N effects must be addressed explicitly. Our direct comparison is performed at identical small N for both quantum and classical descriptions, which isolates the impact of quantum operator non-commutativity and correlations versus the classical vector model. While classical dynamics at small N does not automatically match its own large-N limit, the discrepancies we report are already present in this regime and align with the absence of quantum entanglement in the classical case. In the revision we will add classical simulations at larger N (up to N=50) to show that the short- and long-time discrepancies persist qualitatively, and we will revise the abstract and discussion to state that the differences are demonstrated within the N range accessible to exact quantum calculations rather than claiming they are proven in the thermodynamic limit. revision: partial

  2. Referee: [Methods and Results sections] No quantitative error analysis, specific Hamiltonian parameters (e.g., dipolar coupling strengths, lattice geometry), or explicit initial-state definitions are provided to allow independent assessment of whether the short- and long-time discrepancies are robust or sensitive to choices. The absence of these details makes it impossible to judge whether the 'significant discrepancies' exceed numerical or statistical uncertainties.

    Authors: We acknowledge that these details were insufficiently specified. In the revised manuscript we will add: (i) explicit initial-state definitions (e.g., fully polarized product states along a chosen axis), (ii) concrete Hamiltonian parameters (dipolar coupling strength normalized to 1, lattice geometry and nearest-neighbor distances), and (iii) quantitative error analysis including time-step convergence for classical integration and numerical tolerance for the quantum diagonalization. These additions will allow readers to verify that the reported discrepancies exceed numerical uncertainties. revision: yes

standing simulated objections not resolved
  • Exact quantum diagonalization remains computationally infeasible for N ≫ 10, so a direct quantum-classical comparison at large N cannot be provided.

Circularity Check

0 steps flagged

No circularity: direct numerical comparison without self-referential derivations

full rationale

The paper conducts a side-by-side numerical study of quantum (exact Schrödinger evolution) versus classical (vector) dipolar spin dynamics, benchmarking both against free induction decay for small N. No derivation chain, fitted parameters, or predictions are described that reduce by construction to the inputs; the central claim is an observed discrepancy attributed to operator versus vector mechanics. No self-citations are invoked as load-bearing uniqueness theorems, no ansatz is smuggled, and no known result is renamed. The work is self-contained as an empirical comparison.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract provides no information on free parameters, axioms, or new entities; ledger is therefore empty.

pith-pipeline@v0.9.0 · 5384 in / 1007 out tokens · 24230 ms · 2026-05-08T04:11:02.684300+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

16 extracted references · 16 canonical work pages

  1. [1]

    S. J. Knak Jensen, O. Platz, Free-Induction-Decay Shapes in a Dipolar-Coupled Rigid Lattice of Infinite Nuclear Spins, Phys. Rev. B 7 (1973) 31

    work page 1973

  2. [2]

    S. J. Jensen and E. K. Hansen, Dynamics of classical spins with dipolar coupling in a rigid lattice at high temperatures, Phys. Rev. B 13 (1976) 1903

    work page 1976

  3. [3]

    A. A. Lundin, V. E. Zobov, Simulation of the nuclear magnetic system of a crystal by the system of classical magnetic moments, J. Magn. Reson. 26 (1977) 229

    work page 1977

  4. [4]

    C. Tang, J. S. Waugh, Dynamics of classical spins on a lattice: Spin diffusion, Phys. Rev. B 45 (1992) 748

    work page 1992

  5. [5]

    T. A. Elsayed, B. V. Fine, Effectiveness of classical spin simulations for describing NMR relaxation of quantum spins, Phys. Rev. B 91 (2015) 094424

    work page 2015

  6. [6]

    Henner, K

    V. Henner, K. Klots, T. Belozerova, Simulation of Pake doublet with classical spins and correspondence between the quantum and classical approaches, Eur. Phys. J. B 89 (2016) 264

    work page 2016

  7. [7]

    Henner, A

    V. Henner, A. Klots, A. Nepomnyashchy, T. Belozerova, Correspondence Principle for Spin Systems: Simulations of Free Induction Decay with Classical and Quantum Spins , Applied Magnetic Resonance, 52 (2021) 859

    work page 2021

  8. [8]

    Tsiberkin, T

    K. Tsiberkin, T. Belozerova, V. Henner, Simulation of free induction decay at low-temperature with spin waves and classical spins, Eur. Phys. J. B 92 (2019) 140

    work page 2019

  9. [9]

    Davis, I

    C. Davis, I. Kaganov, V. Henner, Superradiation in magnetic resonance , Phys. Rev . B 62 (2000) 12328-12337

    work page 2000

  10. [10]

    Henner, H

    V. Henner, H. Desvaux, T. Belozerova, D. Marion, P. Kharebov, A. Klots, Collective effects due to dipolar fields as the origin of the extremely random behavior in hyperpolarized NMR maser: A theoretical and numerical study, J. Chem. Phys. 139 (2013) 144111

    work page 2013

  11. [11]

    Kharebov, V

    P. Kharebov, V. Henner, V. Yukalov, Optimal conditions for magnetization reversal of nanocluster assemlies with random properties, J. Appl. Phys. 113 (2013) 043902. 27

    work page 2013

  12. [12]

    Davis, V

    C. Davis, V. Henner, A. Tchernatinsky, I. Kaganov, Spin-system radio -frequency superradiation: A phenomenological study and comparison with numeric simulations , Phys. Rev. B 72 (2005) 054406-054415

    work page 2005

  13. [13]

    de Wijn, B

    A.S. de Wijn, B. Hess and B.V. Fine, Largest Lyapunov Exponents for Lattices of Interacting Classical Spins, Phys. Rev. Lett. 109 (2012) 034101

    work page 2012

  14. [14]

    G. E. Pake, Nuclear Resonance Absorption in Hydrated Crystals: Fine Structure of the Proton Line, J. Chem. Phys. 16 (1948) 327-336

    work page 1948

  15. [15]

    Abragam, Principles of Nuclear Magnetism (Clarendon Press, Oxford, 1961)

    A. Abragam, Principles of Nuclear Magnetism (Clarendon Press, Oxford, 1961)

    work page 1961

  16. [16]

    Pikovsky and A

    A. Pikovsky and A. Politi, Lyapunov Exponents: A Tool to Explore Complex Dynamics , (Cambridge University Press, 2016)

    work page 2016