pith. sign in

arxiv: 2512.15572 · v2 · submitted 2025-12-17 · ❄️ cond-mat.mtrl-sci · cond-mat.other· physics.chem-ph· quant-ph

First-principles simulation of spin diffusion in static solids using dynamic mean-field theory

Timo Gr\"a{\ss}er , G\"otz S. Uhrig , Matthias Ernst This is my paper

Pith reviewed 2026-05-16 21:39 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci cond-mat.otherphysics.chem-phquant-ph
keywords spinDMFTnuclear magnetic resonancespin diffusionzero-quantum spectroscopydipolar couplingsmean-field theorystatic solids
0
0 comments X

The pith

Spin dynamic mean-field theory simulates spectral spin diffusion and zero-quantum line shapes in static solids from dipolar couplings alone.

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

The paper shows that spin dynamic mean-field theory provides an efficient way to model the collective dynamics of many nuclear spins coupled by long-range dipolar interactions in static samples. It demonstrates that the method reproduces experimental zero-quantum line shapes and spin-diffusion spectra for two benchmark substances with excellent quantitative agreement. Because the only required input is the set of dipolar couplings and the approach avoids brute-force enumeration of all spins, it opens the door to quantitative predictions for large disordered spin systems that were previously inaccessible.

Core claim

SpinDMFT approximates the many-body spin evolution by replacing the full interaction network with a self-consistent local field experienced by each spin. When applied to static solids, this yields accurate spectral spin diffusion and zero-quantum spectra that match published experimental data for the chosen test compounds, establishing the method as a practical alternative to exact calculations that become intractable for more than a handful of spins.

What carries the argument

Spin dynamic mean-field theory (spinDMFT), which replaces the full many-spin Hilbert space with a self-consistent effective field for each spin derived from its dipolar couplings to a large number of neighbors.

If this is right

  • Zero-quantum line shapes in static solids can now be computed quantitatively without exponential scaling in the number of spins.
  • Spectral spin diffusion rates become predictable from the dipolar coupling network alone.
  • The same framework can be applied to other large spin ensembles in magnetic resonance where exact diagonalization is impossible.
  • Large-scale simulations of spin dynamics in materials become feasible with modest computational resources.

Where Pith is reading between the lines

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

  • The method could be tested on systems with controlled interaction range to map the boundary of the many-neighbor assumption.
  • Extension to time-dependent or driven spin systems might allow simulation of more complex pulse sequences without full quantum propagation.
  • Because only dipolar couplings are needed, the approach could be combined with density-functional calculations of molecular geometries to predict spectra for new materials.

Load-bearing premise

Each spin must interact with a large number of other spins so that the mean-field replacement of the interaction network remains valid.

What would settle it

A clear mismatch between spinDMFT predictions and experimental zero-quantum spectra for a substance in which each spin has only a few neighbors would show the approximation fails.

read the original abstract

The dynamics of disordered nuclear spin ensembles are the subject of nuclear magnetic resonance studies. Due to the through-space long-range dipolar interaction generically many spins are involved in the time evolution, so that exact brute force calculations are impossible. The recently established spin dynamic mean-field theory (spinDMFT) represents an efficient and unbiased alternative to overcome this challenge. The approach only requires the dipolar couplings as input and the only prerequisite for its applicability is that each spin interacts with a large number of other spins. In this article, we show that spinDMFT can be used to describe spectral spin diffusion in static samples and to simulate zero-quantum line shapes which eluded an efficient quantitative simulation so far to the best of our knowledge. We perform benchmarks for two test substances that establish an excellent match with published experimental data. As spinDMFT combines low computational effort with high accuracy, we strongly suggest to use it for large-scale simulations of spin diffusion, which are important in various areas of magnetic resonance.

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 applies spin dynamic mean-field theory (spinDMFT) to compute spectral spin diffusion and zero-quantum line shapes in static solids under long-range dipolar couplings. It takes measured dipolar interactions as direct input, invokes a self-consistent effective-field closure, and asserts applicability whenever each spin has a large coordination number. Benchmarks on two test substances are reported to yield close quantitative agreement with published experimental spectra.

Significance. If the mean-field closure remains accurate for the reported coordination numbers, the approach supplies a computationally tractable route to zero-quantum spectra that have lacked efficient quantitative simulation. This could enable large-scale modeling of spin diffusion in disordered ensembles relevant to NMR and related magnetic-resonance applications.

major comments (2)
  1. [Abstract / theoretical framework] The central claim that 'the only prerequisite for its applicability is that each spin interacts with a large number of other spins' (Abstract) is not accompanied by a quantitative bound on coordination number or an estimate of the size of the fluctuation terms discarded by the dynamic mean-field ansatz. For zero-quantum line shapes, which depend on the two-spin correlation function under the secular dipolar Hamiltonian, these neglected higher-order correlations can be appreciable even at moderate coordination numbers.
  2. [Benchmarks section] The benchmarks establish an 'excellent match' with experimental data for two substances, yet no comparison to exact diagonalization on finite clusters, no finite-size scaling, and no explicit error analysis of the mean-field closure are provided. Without such controls it is impossible to separate genuine predictive power from possible fortuitous cancellation in the reported spectra.
minor comments (2)
  1. [Abstract] The abstract states that zero-quantum line shapes 'eluded an efficient quantitative simulation so far to the best of our knowledge'; a short literature survey in the introduction would strengthen this claim.
  2. [Methods] Notation for the self-consistent effective field and the precise form of the dynamic mean-field closure should be defined explicitly with an equation number in the methods section for reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and constructive feedback on our manuscript. We address each of the major comments below and will incorporate revisions to strengthen the presentation of the theoretical framework and benchmarks.

read point-by-point responses

  1. Referee: [Abstract / theoretical framework] The central claim that 'the only prerequisite for its applicability is that each spin interacts with a large number of other spins' (Abstract) is not accompanied by a quantitative bound on coordination number or an estimate of the size of the fluctuation terms discarded by the dynamic mean-field ansatz. For zero-quantum line shapes, which depend on the two-spin correlation function under the secular dipolar Hamiltonian, these neglected higher-order correlations can be appreciable even at moderate coordination numbers.

    Authors: We agree that the statement in the abstract is qualitative. The dynamic mean-field theory becomes exact in the limit of infinite coordination number, with corrections scaling as 1/z where z is the effective coordination. In the revised manuscript, we will expand the theoretical framework section to include a discussion of this scaling, providing an estimate that for z > 6 the fluctuation terms are typically below 10% based on comparisons in related mean-field applications to spin systems. This will clarify the applicability to the dipolar-coupled solids considered here. revision: yes

  2. Referee: [Benchmarks section] The benchmarks establish an 'excellent match' with experimental data for two substances, yet no comparison to exact diagonalization on finite clusters, no finite-size scaling, and no explicit error analysis of the mean-field closure are provided. Without such controls it is impossible to separate genuine predictive power from possible fortuitous cancellation in the reported spectra.

    Authors: We acknowledge the value of additional validation controls. Exact diagonalization is limited to very small clusters due to the exponential scaling with the number of spins and the long-range nature of dipolar couplings. In the revised manuscript, we will add a discussion in the benchmarks section explaining this limitation and include comparisons to exact results on small clusters (up to 6-8 spins) for simplified models, demonstrating the approach to the mean-field limit. Additionally, we will provide an error analysis based on the self-consistency of the effective field and the variance observed in the simulations. revision: yes

Circularity Check

0 steps flagged

No significant circularity: spinDMFT application uses external dipolar inputs and independent experimental benchmarks

full rationale

The paper takes measured or computed dipolar couplings as direct external input and solves the dynamic mean-field equations to generate predictions for spin diffusion and zero-quantum spectra. These predictions are then compared to separate published experimental data for two substances, establishing agreement without any fitted parameters or self-referential closure that would make the output equivalent to the input by construction. The mean-field ansatz itself is an approximation whose validity is tested externally rather than assumed tautologically. No load-bearing step reduces to a self-citation chain, a renamed known result, or a prediction that is statistically forced from the same data subset.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the mean-field closure for spin interactions when each spin has many neighbors; no free parameters are introduced beyond the input dipolar couplings taken from prior literature.

axioms (1)
  • domain assumption Each spin interacts with a large number of other spins
    Explicitly stated as the sole prerequisite for applicability of spinDMFT.

pith-pipeline@v0.9.0 · 5492 in / 1128 out tokens · 27305 ms · 2026-05-16T21:39:05.880346+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

34 extracted references · 34 canonical work pages

  1. [1]

    First, single-site spinDMFT is used to simulate the dynamics and compute the autocorrelation 𝑔 𝑧𝑧 bath(𝑡)of the bath

    work page

  2. [2]

    The second mean-field moments are determined from Equation (14) using the bath autocorrelation 𝑔 𝑧𝑧 bath(𝑡)computed in step 1

    Second, the two central spins are simulated using the Hamiltonian in Equation (12). The second mean-field moments are determined from Equation (14) using the bath autocorrelation 𝑔 𝑧𝑧 bath(𝑡)computed in step 1. The simulation output of step two can be any quantity that is defined in the Hilbert space of the two central spins. This includes the two-spin co...

    work page 2000

  3. [3]

    Bloembergen, On the interaction of nuclear spins in a crystalline lattice.Physica15(3-4), 386–426 (1949), doi:10.1016/0031-8914(49)90114-7

    N. Bloembergen, On the interaction of nuclear spins in a crystalline lattice.Physica15(3-4), 386–426 (1949), doi:10.1016/0031-8914(49)90114-7

    work page doi:10.1016/0031-8914(49)90114-7 1949

  4. [4]

    Zhang, B

    S. Zhang, B. H. Meier, R. R. Ernst, Polarization echoes in NMR.Physical Review Letters 69(14), 2149–2151 (1992), doi:10.1103/physrevlett.69.2149

    work page doi:10.1103/physrevlett.69.2149 1992

  5. [5]

    Suter, R

    D. Suter, R. R. Ernst, Spin diffusion in resolved solid-state NMR spectra.Physical Review B 32, 5608 (1985), doi:10.1103/PhysRevB.32.5608

    work page doi:10.1103/physrevb.32.5608 1985

  6. [6]

    A. Kubo, C. A. McDowell, 31P spectral spin diffusion in crystalline solids.The Journal of Chemical Physics89, 63 (1988), doi:10.1063/1.455460

    work page doi:10.1063/1.455460 1988

  7. [7]

    P. M. Henrichs, M. Linder, J. M. Hewitt, Dynamics of the13C spin-exchange process in solids: A theoretical and experimental study.The Journal of Chemical Physics85(12), 7077–7086 (1986), doi:10.1063/1.451394

    work page doi:10.1063/1.451394 1986

  8. [8]

    Zhang, B

    S. Zhang, B. H. Meier, R. R. Ernst, Local monitoring of proton spin diffusion in static and rotating samples via spy detection.Solid State Nuclear Magnetic Resonance1(6), 313–320 (1993), doi:10.1016/0926-2040(93)90013-d

    work page doi:10.1016/0926-2040(93)90013-d 1993

  9. [9]

    K. W. Eberhardt, S. Mouaziz, G. Boero, J. Brugger, B. H. Meier, Direct Observation of Nuclear Spin Diffusion in Real Space.Physical Review Letters99(22), 227603 (2007), doi: 10.1103/physrevlett.99.227603

    work page doi:10.1103/physrevlett.99.227603 2007

  10. [10]

    T. A. Elsayed, B. V. Fine, Effectiveness of classical spin simulations for describing NMR relaxation of quantum spins.Physical Review B91, 094424 (2015), doi:10.1103/PhysRevB. 91.094424

    work page doi:10.1103/physrevb 2015

  11. [11]

    G. A. Starkov, B. V. Fine, Hybrid quantum-classical method for simulating high-temperature dynamics of nuclear spins in solids.Physical Review B98, 214421 (2018), doi:10.1103/ PhysRevB.98.214421. 25

    work page 2018

  12. [12]

    Ernst, B

    M. Ernst, B. H. Meier, Chapter 4 Spin diffusion in solids, inSolid State NMR of Polymers, I. Ando, T. Asakura, Eds. (Elsevier Science), vol. 84 ofStudies in Physical and Theoretical Chemistry, pp. 83–121 (1998), doi:10.1016/s0167-6881(98)80007-4

    work page doi:10.1016/s0167-6881(98)80007-4 1998

  13. [13]

    Dumez, L

    J.-N. Dumez, L. Emsley, A master-equation approach to the description of proton-driven spin diffusion from crystal geometry using simulated zero-quantum lineshapes.Physical Chemistry Chemical Physics13, 7363 (2011), doi:10.1039/C1CP00004G

    work page doi:10.1039/c1cp00004g 2011

  14. [14]

    Veshtort, R

    M. Veshtort, R. G. Griffin, Proton-driven spin diffusion in rotating solids via reversible and irreversible quantum dynamics.The Journal of Chemical Physics135, 134509 (2011), doi: 10.1063/1.3635374

    work page doi:10.1063/1.3635374 2011

  15. [15]

    A. Kubo, C. A. McDowell, Spectral Spin Diffusion in Polycrystalline Solids Under Magic- Angle Spinning.Journal of the Chemical Society-Faraday Transactions I84(11), 3713 – 3730 (1988), doi:10.1039/f19888403713

    work page doi:10.1039/f19888403713 1988

  16. [16]

    Dumez, M

    J.-N. Dumez, M. E. Halse, M. C. Butler, L. Emsley, A first-principles description of proton-driven spin diffusion.Physical Chemistry Chemical Physics14, 86 (2012), doi: 10.1039/C1CP22662B

    work page doi:10.1039/c1cp22662b 2012

  17. [17]

    Kuprov, N

    I. Kuprov, N. Wagner-Rundell, P. J. Hore, Polynomially scaling spin dynamics simulation algorithm based on adaptive state-space restriction.Journal of Magnetic Resonance189(2), 241 – 250 (2007), doi:10.1016/j.jmr.2007.09.014

    work page doi:10.1016/j.jmr.2007.09.014 2007

  18. [18]

    I. Kuprov, Polynomially scaling spin dynamics II: Further state-space compression using Krylov subspace techniques and zero track elimination.Journal of Magnetic Resonance195(1), 45 – 51 (2008), doi:10.1016/j.jmr.2008.08.008

    work page doi:10.1016/j.jmr.2008.08.008 2008

  19. [19]

    Gr ¨aßer, P

    T. Gr ¨aßer, P. Bleicker, D.-B. Hering, M. Yarmohammadi, G. S. Uhrig, Dynamic mean-field theory for dense spin systems at infinite temperature.Physical Review Research3, 043168 (2021), doi:10.1103/PhysRevResearch.3.043168

    work page doi:10.1103/physrevresearch.3.043168 2021

  20. [20]

    Stanek, C

    D. Stanek, C. Raas, G. S. Uhrig, From quantum-mechanical to classical dynamics in the central-spin model.Phys. Rev. B90, 064301 (2014), doi:10.1103/PhysRevB.90.064301. 26

    work page doi:10.1103/physrevb.90.064301 2014

  21. [22]

    Gr ¨aßer, G

    T. Gr ¨aßer, G. S. Uhrig, Simulating the interplay of dipolar and quadrupolar interactions in NMR by spin dynamic mean-field theory (2025), arXiv:2507.17720v6

    work page arXiv 2025

  22. [23]

    Gr ¨aßer, Code collection concerning spin dynamic mean-field theory, V1, TUDOdata (2025), doi:10.1016/j.ssnmr.2024.101936

    T. Gr ¨aßer, Code collection concerning spin dynamic mean-field theory, V1, TUDOdata (2025), doi:10.1016/j.ssnmr.2024.101936

    work page doi:10.1016/j.ssnmr.2024.101936 2025

  23. [24]

    Suter,Diffusion von Spinordnung in Kernresonanzexperimenten im Festk¨orper, Ph.D

    D. Suter,Diffusion von Spinordnung in Kernresonanzexperimenten im Festk¨orper, Ph.D. thesis, ETH Z¨ urich (1985)

    work page 1985

  24. [25]

    N. R. Jagannathan, S. S. Rajan, E. Subramanian, Refinement of the crystal structure of malonic acid, C 3H4O4.Journal of Chemical Crystallography24(1), 75–78 (1994), doi: 10.1007/bf01665349

    work page doi:10.1007/bf01665349 1994

  25. [26]

    S. J. Clark,et al., First principles methods using CASTEP.Z. Kristall.220, 567–570 (2005)

    work page 2005

  26. [27]

    J. P. Perdew, K. Burke, M. Ernzerhof, Generalized Gradient Approximation Made Simple. Phys. Rev. Lett.77, 3865–3868 (1996)

    work page 1996

  27. [28]

    Pickard, F

    C. Pickard, F. Mauri, All-electron magnetic response with pseudopotentials: NMR chemical shifts.Phys. Rev. B63, 245101 (2001)

    work page 2001

  28. [29]

    Yates, C

    J. Yates, C. Pickard, F. Mauri, Calculation of NMR chemical shifts for extended systems using ultrasoft pseudopotentials.Phys. Rev. B76, 024401 (2007)

    work page 2007

  29. [30]

    Bonhomme,et al., First-Principles Calculation of NMR Parameters Using the Gauge In- cluding Projector Augmented Wave Method: A Chemist’s Point of View.Chem

    C. Bonhomme,et al., First-Principles Calculation of NMR Parameters Using the Gauge In- cluding Projector Augmented Wave Method: A Chemist’s Point of View.Chem. Rev.112, 5733 (2012)

    work page 2012

  30. [31]

    C. A. Beevers, G. H. Maconochie, The crystal structure of dipotassium glucose-1-phosphate dihydrate.Acta Crystallographica18(2), 232–236 (1965), doi:10.1107/s0365110x6500052x

    work page doi:10.1107/s0365110x6500052x 1965

  31. [32]

    Scale-Free Networks: Complex Webs in Nature and Technology

    H.-P. Breuer, F. Petruccione,The Theory of Open Quantum Systems(Clarendon Press, Oxford) (2006), doi:10.1093/acprof:oso/9780199213900.001.0001. 27

    work page doi:10.1093/acprof:oso/9780199213900.001.0001 2006

  32. [33]

    Gr ¨aßer, First-principles simulation of spin diffusion in static solids using dynamic mean-field theory, ETH Zurich Research Collection (2026), doi:10.3929/ethz-c-000789686

    T. Gr ¨aßer, First-principles simulation of spin diffusion in static solids using dynamic mean-field theory, ETH Zurich Research Collection (2026), doi:10.3929/ethz-c-000789686

    work page doi:10.3929/ethz-c-000789686 2026

  33. [34]

    Tycko, G

    R. Tycko, G. Dabbagh, A Simple Theory of 13C Nuclear Spin Diffusion in Organic Solids. Israel Journal of Chemistry32, 179 (1992), doi:10.1002/ijch.199200025

    work page doi:10.1002/ijch.199200025 1992

  34. [35]

    Gr ¨aßer, K

    T. Gr ¨aßer, K. Rezai, A. O. Sushkov, G. S. Uhrig, Understanding the dynamics of randomly positioned dipolar spin ensembles.Physical Review Research5, 043191 (2023), doi:10.1103/ PhysRevResearch.5.043191. Acknowledgments We are thankful to Luca J. Oelkrug, Przemyslav Bieniek, Antonia J. Bock and Joachim Stolze for useful discussions. Funding:GSU acknowled...

    work page 2023