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

Journal of Applied Physics

Computer simulation of the approach to equilibrium in the classical Heisenberg chain

Sven P. Rudin, P. Heller

Published Mar 1, 1981• 0 citations
Publisher

Abstract

Computer studies of the approach to equilibrium in the classical 3-component Heisenberg chain, starting from a configuration with a large staggered magnetization (i. e., far from equilibrium) are reported. The equal time transverse and reduced longitudinal correlations g⊥(ℓ, t) and g∥(ℓ, t), respectively, were monitored for times long enough to show the approach to equilibrium. Results for g⊥ confirm predictions by Nelson and Fisher of wave-like behavior with a propagation speed equal to twice the zero-temperature spin wave velocity. The equilibration of g∥ was found to be slower and diffusive in nature. The fitted diffusivity appears to be comparable to the value of the spin diffusivity for the corresponding equilibrium problem.

DOI
10.1063/1.329588

Related papers

Computer studies of spin and energy transport in one-dimensional Heisenberg magnets

N.A. Lurie, David L. Huber, Martin Blume

We have studied the dynamics of the classical Heisenberg chain at finite and infinite temperatures by computer simulation of an array of 4000 spins. The evaluation of the spin-spin and energy-energy wave-vector- and time-dependent correlation functions out to time $10{J}^{\ensuremath{-}1}$ made it possible to determine the spin and energy diffusion constants. At infinite temperature we find ${D}_{S}=(1.33\ifmmode\pm\else\textpm\fi{}0.10)JS{a}^{2}$ and ${D}_{E}=(3\ifmmode\pm\else\textpm\fi{}1)JS{a}^{2}$. Values of ${D}_{S}$ have also been determined, with less precision, for ferromagnetic and antiferromagnetic interactions at temperatures $\frac{{k}_{B}T}{|J|}=1.0 \mathrm{and} 0.5$. Our results are compared with various theoretical estimates. In addition, we calculate the time-dependent spin and energy self-correlation functions.

Physical review. B, Solid state•Mar 1, 1974•99 citations

Spin transport in the XXZ model at high temperatures: Classical dynamics vs. quantum S=1/2 autocorrelations

Robin Steinigeweg

The transport of magnetization is analyzed for the classical Heisenberg chain at and especially above the isotropic point. To this end, the Hamiltonian equations of motion are solved numerically for initial states realizing harmonic-like magnetization profiles of small amplitude and with random phases. Above the isotropic point, the resulting dynamics is observed to be diffusive in a hydrodynamic regime starting at comparatively small times and wave lengths. In particular, hydrodynamic regime and diffusion constant are both found to be in quantitative agreement with close-to-equilibrium results from quantum S =1/2 autocorrelations at high temperatures. At the isotropic point, the resulting dynamics turns out to be non-diffusive at the considered times and wave lengths.

arXiv:1110.6610•Europhysics Letters (EPL)•Mar 1, 2012•26 citations

Low-Temperature Spin-Wave Dynamics in Three-Component Classical Heisenberg Chains

P. Heller, Martin Blume

A detailed quantitative study of the low-temperature spin-wave dynamics has been made for the classical Heisenberg-coupled chain using computer simulation. Results for the relaxation rates are presented and compared with dynamic scaling and other predictions.

Physical Review Letters•Oct 10, 1977•38 citations

Critical dynamics of the body-centered-cubic classical Heisenberg antiferromagnet

Alex Bunker, Kun Chen, David Paul Landau

We have used spin dynamics techniques to perform large-scale simulations of the dynamic behavior of the L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L body-centered-cubic classical Heisenberg antiferromagnet with L\ensuremath{\le}48 at a range of temperatures above and below as well as at the critical point ${\mathit{T}}_{\mathit{c}}$. The temporal evolutions of the spin configurations were determined numerically from coupled equations of motion for individual spins by a fourth-order predictor-corrector method, with initial spin configurations generated by Monte Carlo simulations. The neutron scattering function S(q,\ensuremath{\omega}) was calculated from the space- and time-displaced spin-spin correlation function. We used a previously developed dynamic finite-size scaling theory to extract the dynamic critical exponent z from S(q,\ensuremath{\omega}) at ${\mathit{T}}_{\mathit{c}}$. Our results are in agreement with the theoretical prediction of z=1.5 and with experimental results; however, we find that the asymptotic regime was only entered at L\ensuremath{\approxeq}30. In the analysis of the form of the transverse and longitudinal components of S(q,\ensuremath{\omega}) we found that a central diffusion peak appears below ${\mathit{T}}_{\mathit{c}}$ predominantly in the longitudinal component and remains present through and above ${\mathit{T}}_{\mathit{c}}$. The transverse component of the spin-wave peak is Lorentzian below ${\mathit{T}}_{\mathit{c}}$ but for T\ensuremath{\geqslant}${\mathit{T}}_{\mathit{c}}$ is described best by a more complex functional form. Below ${\mathit{T}}_{\mathit{c}}$ we see evidence of multiple spin-wave peaks in the longitudinal component. \textcopyright{} 1996 The American Physical Society.

Physical review. B, Condensed matter•Oct 1, 1996•31 citations

Time-dependent behavior of classical spin chains at infinite temperature

R. W. Gerling, David Paul Landau

An ultrafast, vectorized spin-dynamics method is used to study the time-dependent properties of classical XY and Heisenberg spin chains at infinite temperature. The decay of the energy-energy and spin-spin correlation functions is oscillatory for short times and at long times is consistent with classical diffusion, although the approach to the asymptotic behavior is extremely slow. We have also calculated S(q,\ensuremath{\omega}) and find clear indication of spin-wave peaks in both models.

Physical review. B, Condensed matter•Nov 1, 1990•44 citations

Propagating Energy Modes in the Classical Heisenberg Chain: "Magnons" and "Second Magnons"

J M Loveluck, Ewald Balcar

The ferromagnetically coupled classical Heisenberg chain in an applied magnetic field has been studied by computer simulation. The results indicate the presence of a second collective mode, in addition to the damped spin-wave-like modes which have previously been observed in the absence of a magnetic field. For intermediate wavelengths, the mode manifests itself by well-defined oscillations in the energy-density correlation function, and by a second peak in the spectrum of the longitudinal spin-density correlation functions.

Physical Review Letters•Jun 4, 1979•32 citations
Computer studies of the approach to equilibrium in the classical 3-component Heisenberg chain, starting from a configuration with a large staggered magnetization (i. e., far from equilibrium) are reported. The equal time transverse and reduced longitudinal correlations g⊥(ℓ, t) and g∥(ℓ, t), respectively, were monitored for times long enough to show the approach to equilibrium. Results for g⊥ confirm predictions by Nelson and Fisher of wave-like behavior with a propagation speed equal to twice the zero-temperature spin wave velocity. The equilibration of g∥ was found to be slower and diffusive in nature. The fitted diffusivity appears to be comparable to the value of the spin diffusivity for the corresponding equilibrium problem.
We have studied the dynamics of the classical Heisenberg chain at finite and infinite temperatures by computer simulation of an array of 4000 spins. The evaluation of the spin-spin and energy-energy wave-vector- and time-dependent correlation functions out to time $10{J}^{\ensuremath{-}1}$ made it possible to determine the spin and energy diffusion constants. At infinite temperature we find ${D}_{S}=(1.33\ifmmode\pm\else\textpm\fi{}0.10)JS{a}^{2}$ and ${D}_{E}=(3\ifmmode\pm\else\textpm\fi{}1)JS{a}^{2}$. Values of DS{D}_{S}DS​ have also been determined, with less precision, for ferromagnetic and antiferromagnetic interactions at temperatures kBT∣J∣=1.0and0.5\frac{{k}_{B}T}{|J|}=1.0 \mathrm{and} 0.5∣J∣kB​T​=1.0and0.5. Our results are compared with various theoretical estimates. In addition, we calculate the time-dependent spin and energy self-correlation functions.
The transport of magnetization is analyzed for the classical Heisenberg chain at and especially above the isotropic point. To this end, the Hamiltonian equations of motion are solved numerically for initial states realizing harmonic-like magnetization profiles of small amplitude and with random phases. Above the isotropic point, the resulting dynamics is observed to be diffusive in a hydrodynamic regime starting at comparatively small times and wave lengths. In particular, hydrodynamic regime and diffusion constant are both found to be in quantitative agreement with close-to-equilibrium results from quantum S =1/2 autocorrelations at high temperatures. At the isotropic point, the resulting dynamics turns out to be non-diffusive at the considered times and wave lengths.
We have used spin dynamics techniques to perform large-scale simulations of the dynamic behavior of the L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L body-centered-cubic classical Heisenberg antiferromagnet with L\ensuremath{\le}48 at a range of temperatures above and below as well as at the critical point Tc{\mathit{T}}_{\mathit{c}}Tc​. The temporal evolutions of the spin configurations were determined numerically from coupled equations of motion for individual spins by a fourth-order predictor-corrector method, with initial spin configurations generated by Monte Carlo simulations. The neutron scattering function S(q,\ensuremath{\omega}) was calculated from the space- and time-displaced spin-spin correlation function. We used a previously developed dynamic finite-size scaling theory to extract the dynamic critical exponent z from S(q,\ensuremath{\omega}) at Tc{\mathit{T}}_{\mathit{c}}Tc​. Our results are in agreement with the theoretical prediction of z=1.5 and with experimental results; however, we find that the asymptotic regime was only entered at L\ensuremath{\approxeq}30. In the analysis of the form of the transverse and longitudinal components of S(q,\ensuremath{\omega}) we found that a central diffusion peak appears below Tc{\mathit{T}}_{\mathit{c}}Tc​ predominantly in the longitudinal component and remains present through and above Tc{\mathit{T}}_{\mathit{c}}Tc​. The transverse component of the spin-wave peak is Lorentzian below Tc{\mathit{T}}_{\mathit{c}}Tc​ but for T\ensuremath{\geqslant}Tc{\mathit{T}}_{\mathit{c}}Tc​ is described best by a more complex functional form. Below Tc{\mathit{T}}_{\mathit{c}}Tc​ we see evidence of multiple spin-wave peaks in the longitudinal component. \textcopyright{} 1996 The American Physical Society.
The ferromagnetically coupled classical Heisenberg chain in an applied magnetic field has been studied by computer simulation. The results indicate the presence of a second collective mode, in addition to the damped spin-wave-like modes which have previously been observed in the absence of a magnetic field. For intermediate wavelengths, the mode manifests itself by well-defined oscillations in the energy-density correlation function, and by a second peak in the spectrum of the longitudinal spin-density correlation functions.