Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2503.19975 · Physical review. B./Physical review. B

Emulation of quantum correlations by classical dynamics in a spin- 1 2 Heisenberg chain

Chaebin Kim, Martin Mourigal

Published Nov 24, 2025• 3 citations
arXivHTMLPDF

Abstract

We simulate the dynamical spin structure factor (DSSF) $\mathcal{S}(q,\ensuremath{\omega})$ of the spin-1/2 Heisenberg antiferromagnetic chain using classical simulations. By employing Landau-Lifshitz Dynamics, we emulate quantum correlations through temperature-dependent corrections, including rescaling of magnetic dipoles and renormalization of exchange interactions. Our results closely match Quantum Monte-Carlo calculations for ${k}_{\mathrm{B}}T/J\ensuremath{\gtrsim}1$, extending the applicability of classical dynamics to the challenging case of gapless excitations. At higher temperatures, our simulations comply with general predictions for uncorrelated paramagnetic fluctuations in the infinite temperature limit. Entanglement witnesses derived from the quantum-equivalent DSSF act as sensitive diagnostics for the quantum-to-classical crossover. Their reliability stems from their dependence on spectral features alone, enabling classical dynamics to emulate quantum thresholds without genuine entanglement. This framework also reproduces transverse spin correlations in finite magnetic fields, in agreement with quantum simulations. Together, our results establish quantum-corrected classical dynamics as a scalable and predictive tool for interpreting scattering experiments and exploring quantum correlations in strongly correlated spin systems.

DOI
10.1103/yvvg-cg2f

Related papers

Spin dynamics of the antiferromagnetic spin- 12 chain at finite magnetic fields and intermediate temperatures

Simon Grossjohann, W. Brenig

We present a study of the dynamic structure factor of the antiferromagnetic spin-$\frac{1}{2}$ Heisenberg chain at finite temperatures and finite magnetic fields. Using quantum Monte Carlo based on the stochastic series expansion and maximum entropy methods, we evaluate the longitudinal and the transverse dynamic structure factors from vanishing magnetic fields up to and above the threshold ${B}_{c}$ for ferromagnetic saturation, as well as for high and for intermediate temperatures. We study the field-induced redistribution of spectral weight contrasting longitudinal versus transverse excitations. At finite fields below saturation incommensurate low-energy modes are found consistent with zero-temperature Bethe ansatz. The crossover between the field-induced ferromagnet above ${B}_{c}$ and the Luttinger liquid below ${B}_{c}$ is analyzed in terms of the transverse spin dynamics. Evaluating sum rules we assess the quality of the analytic continuation and demonstrate excellent consistency of the maximum entropy results.

arXiv:0811.1956•Physical Review B•Mar 10, 2009•25 citations

Evolution of spin excitations in a gapped antiferromagnet from the quantum to the high-temperature limit

M. Kenzelmann, Roger Arthur Cowley, W. J. L. Buyers +3 more

We have mapped from the quantum to the classical limit the spin excitation spectrum of the antiferromagnetic spin-1 Heisenberg chain system ${\mathrm{CsNiCl}}_{3}$ in its paramagnetic phase from $T=5$ to $200\mathrm{K}.$ Neutron scattering shows that the excitations are resonant and dispersive up to at least $T=70\mathrm{K}\ensuremath{\simeq}2.7J,$ but broaden considerably with increasing temperature. The dispersion flattens out with increasing temperature as the resonance energy $\ensuremath{\Delta}$ at the antiferromagnetic wave vector increases and the maximum in the dispersion decreases. The correlation length $\ensuremath{\xi}$ between $T=12$ and $50\mathrm{K}$ is in agreement with quantum Monte Carlo calculations for the spin-1 chain. $\ensuremath{\xi}$ is also consistent with the single mode approximation, suggesting that the excitations are short-lived single particle excitations. Below $T=12\mathrm{K}$ where three-dimensional spin correlations are important, $\ensuremath{\xi}$ is shorter than predicted and the experiment is not consistent with the random phase approximation for coupled quantum chains. At $T=200\mathrm{K},$ the structure factor and second energy moment of the excitation spectrum are in excellent agreement with the high-temperature series expansion.

arXiv:cond-mat/0112188•Physical review. B, Condensed matter•Nov 8, 2002•18 citations

Dynamical properties of the S=12 random Heisenberg chain

Yu-Rong Shu, Maxime Dupont, Dao‐Xin Yao +2 more

We study dynamical properties at finite temperature $(T)$ of Heisenberg spin chains with random antiferromagnetic exchange couplings, which realize the random singlet phase in the low-energy limit, using three complementary numerical methods: exact diagonalization, matrix-product-state algorithms, and stochastic analytic continuation of quantum Monte Carlo results in imaginary time. Specifically, we investigate the dynamic spin structure factor $S(q,\ensuremath{\omega})$ and its $\ensuremath{\omega}\ensuremath{\rightarrow}0$ limit, which are closely related to inelastic neutron scattering and nuclear magnetic resonance (NMR) experiments (through the spin-lattice relaxation rate $1/{T}_{1})$. Our study reveals a continuous narrow band of low-energy excitations in $S(q,\ensuremath{\omega})$, extending throughout the $q$ space, instead of being restricted to $q\ensuremath{\approx}0$ and $q\ensuremath{\approx}\ensuremath{\pi}$ as found in the uniform system. Close to $q=\ensuremath{\pi}$, the scaling properties of these excitations are well captured by the random-singlet theory, but disagreements also exist with some aspects of the predicted $q$ dependence further away from $q=\ensuremath{\pi}$. Furthermore we also find spin diffusion effects close to $q=0$ that are not contained within the random-singlet theory but give non-negligible contributions to the mean $1/{T}_{1}$. To compare with NMR experiments, we consider the distribution of the local relaxation rates $1/{T}_{1}$. We show that the local $1/{T}_{1}$ values are broadly distributed, approximately according to a stretched exponential. The mean $1/{T}_{1}$ first decreases with $T$, but below a crossover temperature it starts to increase and likely diverges in the limit of a small nuclear resonance frequency ${\ensuremath{\omega}}_{0}$. Although a similar divergent behavior has been predicted and experimentally observed for the static uniform susceptibility, this divergent behavior of the mean $1/{T}_{1}$ has never been experimentally observed. Indeed, we show that the divergence of the mean $1/{T}_{1}$ is due to rare events in the disordered chains and is concealed in experiments, where the typical $1/{T}_{1}$ value is accessed.

arXiv:1712.01701•Physical review. B./Physical review. B•Mar 28, 2018•31 citations

Reconstructing the spatial structure of quantum correlations in materials

Allen O. Scheie, Pontus Laurell, Elbio R. Dagotto +2 more

Quantum correlations are a fundamental property of quantum many-body states. Yet they remain experimentally elusive, hindering certification of genuine quantum behavior, especially in quantum materials. Here we show that the momentum-dependent dynamical susceptibility measured via inelastic neutron scattering enables the systematic reconstruction of a general family of quantum correlation functions, which express the degree of quantum coherence in the fluctuations of two spins at arbitrary mutual distance. Using neutron scattering data on the compound KCuF$_3$ $\unicode{x2014}$ a system of weakly coupled $S=1/2$ Heisenberg chains $\unicode{x2014}$ and of numerically exact quantum Monte Carlo data, we show that quantum correlations possess a radically different spatial structure with respect to conventional correlations. Indeed, they exhibit a new emergent length scale $\unicode{x2014}$ the quantum coherence length $\unicode{x2014}$ which is finite at any finite temperature (including when long-range magnetic order develops). Moreover, we show theoretically that coupled Heisenberg spin chains exhibit a form of quantum monogamy, with a trade-off between quantum correlations along and transverse to the spin chains. These results highlight real-space quantum correlators as an informative, model-independent means of probing the underlying quantum state of real quantum materials.

arXiv:2306.11723•Jun 20, 2023•0 citations

Finite temperature properties of alternating ferro-antiferromagnetic Heisenberg chains

Niklas Casper, Wolfram Brenig

We present a study of the spin-1/2 Heisenberg chain with alternating ferro- and antiferromagnetic exchange, focusing on the role of the exchange couplings to cover both, dimer and Haldane limit. Employing a complementary combination of perturbation theory and quantum Monte Carlo simulation, we report results for the magnetic susceptibility as well as the dynamic structure factor over a wide range of coupling constants and for different temperatures to extract the spin gap. For a small ferromagnetic coupling, we find good agreement between perturbation theory and quantum Monte Carlo. For arbitrary exchange couplings we show that the dynamic structure factor, obtained from quantum Monte Carlo, scales between triplons and a Haldane chain spectrum. Finally, we contrast our findings for the spin gap versus the exchange couplings against existing literature.

arXiv:2001.04696•Jan 14, 2020•0 citations

Ground State of the S = 1 / 2 Heisenberg Spin Chain with Random Ferromagnetic and Antiferromagnetic Couplings

Sibei Li, Hui Shao, Anders W. Sandvik

We study the Heisenberg S=1/2 chain with random ferro- and antiferromagnetic couplings using quantum Monte Carlo simulations at ultra-low temperatures, converging to the ground state. Finite-size scaling of correlation functions and excitation gaps demonstrate an exotic critical state in qualitative agreement with previous strong-disorder renormalization group calculations but with scaling exponents depending on the coupling distribution. We find dual scaling regimes of the transverse correlations versus the distance, with an L independent form C(r)=r^{-μ} for r≪L and C(r,L)=L^{-η}f(r/L) for r/L>0, where μ>η and the scaling function is delivered by our analysis. These results are at variance with previous spin-wave and density-matrix renormalization group calculations, thus highlighting the power of unbiased quantum Monte Carlo simulations.

arXiv:2408.12911•Physical Review Letters•Feb 24, 2025•2 citations
We simulate the dynamical spin structure factor (DSSF) S(q,\ensuremathω)\mathcal{S}(q,\ensuremath{\omega})S(q,\ensuremathω) of the spin-1/2 Heisenberg antiferromagnetic chain using classical simulations. By employing Landau-Lifshitz Dynamics, we emulate quantum correlations through temperature-dependent corrections, including rescaling of magnetic dipoles and renormalization of exchange interactions. Our results closely match Quantum Monte-Carlo calculations for kBT/J\ensuremath≳1{k}_{\mathrm{B}}T/J\ensuremath{\gtrsim}1kB​T/J\ensuremath≳1, extending the applicability of classical dynamics to the challenging case of gapless excitations. At higher temperatures, our simulations comply with general predictions for uncorrelated paramagnetic fluctuations in the infinite temperature limit. Entanglement witnesses derived from the quantum-equivalent DSSF act as sensitive diagnostics for the quantum-to-classical crossover. Their reliability stems from their dependence on spectral features alone, enabling classical dynamics to emulate quantum thresholds without genuine entanglement. This framework also reproduces transverse spin correlations in finite magnetic fields, in agreement with quantum simulations. Together, our results establish quantum-corrected classical dynamics as a scalable and predictive tool for interpreting scattering experiments and exploring quantum correlations in strongly correlated spin systems.
We present a study of the dynamic structure factor of the antiferromagnetic spin-12\frac{1}{2}21​ Heisenberg chain at finite temperatures and finite magnetic fields. Using quantum Monte Carlo based on the stochastic series expansion and maximum entropy methods, we evaluate the longitudinal and the transverse dynamic structure factors from vanishing magnetic fields up to and above the threshold Bc{B}_{c}Bc​ for ferromagnetic saturation, as well as for high and for intermediate temperatures. We study the field-induced redistribution of spectral weight contrasting longitudinal versus transverse excitations. At finite fields below saturation incommensurate low-energy modes are found consistent with zero-temperature Bethe ansatz. The crossover between the field-induced ferromagnet above Bc{B}_{c}Bc​ and the Luttinger liquid below Bc{B}_{c}Bc​ is analyzed in terms of the transverse spin dynamics. Evaluating sum rules we assess the quality of the analytic continuation and demonstrate excellent consistency of the maximum entropy results.
We have mapped from the quantum to the classical limit the spin excitation spectrum of the antiferromagnetic spin-1 Heisenberg chain system CsNiCl3{\mathrm{CsNiCl}}_{3}CsNiCl3​ in its paramagnetic phase from T=5T=5T=5 to 200K.200\mathrm{K}.200K. Neutron scattering shows that the excitations are resonant and dispersive up to at least $T=70\mathrm{K}\ensuremath{\simeq}2.7J,$ but broaden considerably with increasing temperature. The dispersion flattens out with increasing temperature as the resonance energy $\ensuremath{\Delta}$ at the antiferromagnetic wave vector increases and the maximum in the dispersion decreases. The correlation length $\ensuremath{\xi}$ between T=12T=12T=12 and 50K50\mathrm{K}50K is in agreement with quantum Monte Carlo calculations for the spin-1 chain. $\ensuremath{\xi}$ is also consistent with the single mode approximation, suggesting that the excitations are short-lived single particle excitations. Below T=12KT=12\mathrm{K}T=12K where three-dimensional spin correlations are important, $\ensuremath{\xi}$ is shorter than predicted and the experiment is not consistent with the random phase approximation for coupled quantum chains. At T=200K,T=200\mathrm{K},T=200K, the structure factor and second energy moment of the excitation spectrum are in excellent agreement with the high-temperature series expansion.
We study dynamical properties at finite temperature (T)(T)(T) of Heisenberg spin chains with random antiferromagnetic exchange couplings, which realize the random singlet phase in the low-energy limit, using three complementary numerical methods: exact diagonalization, matrix-product-state algorithms, and stochastic analytic continuation of quantum Monte Carlo results in imaginary time. Specifically, we investigate the dynamic spin structure factor $S(q,\ensuremath{\omega})$ and its $\ensuremath{\omega}\ensuremath{\rightarrow}0$ limit, which are closely related to inelastic neutron scattering and nuclear magnetic resonance (NMR) experiments (through the spin-lattice relaxation rate 1/T1)1/{T}_{1})1/T1​). Our study reveals a continuous narrow band of low-energy excitations in $S(q,\ensuremath{\omega})$, extending throughout the qqq space, instead of being restricted to $q\ensuremath{\approx}0$ and $q\ensuremath{\approx}\ensuremath{\pi}$ as found in the uniform system. Close to $q=\ensuremath{\pi}$, the scaling properties of these excitations are well captured by the random-singlet theory, but disagreements also exist with some aspects of the predicted qqq dependence further away from $q=\ensuremath{\pi}$. Furthermore we also find spin diffusion effects close to q=0q=0q=0 that are not contained within the random-singlet theory but give non-negligible contributions to the mean 1/T11/{T}_{1}1/T1​. To compare with NMR experiments, we consider the distribution of the local relaxation rates 1/T11/{T}_{1}1/T1​. We show that the local 1/T11/{T}_{1}1/T1​ values are broadly distributed, approximately according to a stretched exponential. The mean 1/T11/{T}_{1}1/T1​ first decreases with TTT, but below a crossover temperature it starts to increase and likely diverges in the limit of a small nuclear resonance frequency ${\ensuremath{\omega}}_{0}$. Although a similar divergent behavior has been predicted and experimentally observed for the static uniform susceptibility, this divergent behavior of the mean 1/T11/{T}_{1}1/T1​ has never been experimentally observed. Indeed, we show that the divergence of the mean 1/T11/{T}_{1}1/T1​ is due to rare events in the disordered chains and is concealed in experiments, where the typical 1/T11/{T}_{1}1/T1​ value is accessed.
Quantum correlations are a fundamental property of quantum many-body states. Yet they remain experimentally elusive, hindering certification of genuine quantum behavior, especially in quantum materials. Here we show that the momentum-dependent dynamical susceptibility measured via inelastic neutron scattering enables the systematic reconstruction of a general family of quantum correlation functions, which express the degree of quantum coherence in the fluctuations of two spins at arbitrary mutual distance. Using neutron scattering data on the compound KCuF3_33​ $\unicode{x2014}$ a system of weakly coupled S=1/2S=1/2S=1/2 Heisenberg chains $\unicode{x2014}$ and of numerically exact quantum Monte Carlo data, we show that quantum correlations possess a radically different spatial structure with respect to conventional correlations. Indeed, they exhibit a new emergent length scale $\unicode{x2014}$ the quantum coherence length $\unicode{x2014}$ which is finite at any finite temperature (including when long-range magnetic order develops). Moreover, we show theoretically that coupled Heisenberg spin chains exhibit a form of quantum monogamy, with a trade-off between quantum correlations along and transverse to the spin chains. These results highlight real-space quantum correlators as an informative, model-independent means of probing the underlying quantum state of real quantum materials.
We present a study of the spin-1/2 Heisenberg chain with alternating ferro- and antiferromagnetic exchange, focusing on the role of the exchange couplings to cover both, dimer and Haldane limit. Employing a complementary combination of perturbation theory and quantum Monte Carlo simulation, we report results for the magnetic susceptibility as well as the dynamic structure factor over a wide range of coupling constants and for different temperatures to extract the spin gap. For a small ferromagnetic coupling, we find good agreement between perturbation theory and quantum Monte Carlo. For arbitrary exchange couplings we show that the dynamic structure factor, obtained from quantum Monte Carlo, scales between triplons and a Haldane chain spectrum. Finally, we contrast our findings for the spin gap versus the exchange couplings against existing literature.
We study the Heisenberg S=1/2 chain with random ferro- and antiferromagnetic couplings using quantum Monte Carlo simulations at ultra-low temperatures, converging to the ground state. Finite-size scaling of correlation functions and excitation gaps demonstrate an exotic critical state in qualitative agreement with previous strong-disorder renormalization group calculations but with scaling exponents depending on the coupling distribution. We find dual scaling regimes of the transverse correlations versus the distance, with an L independent form C(r)=r^{-μ} for r≪L and C(r,L)=L^{-η}f(r/L) for r/L>0, where μ>η and the scaling function is delivered by our analysis. These results are at variance with previous spin-wave and density-matrix renormalization group calculations, thus highlighting the power of unbiased quantum Monte Carlo simulations.