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:cond-mat/0301526 · Physical review. B, Condensed matter

Critical dynamics of the simple-cubic Heisenberg antiferromagnet RbMnF3: Extrapolation to q=0

Shan-Ho Tsai, David Paul Landau

Published Mar 17, 2003• 23 citations
arXivHTMLPDF

Abstract

Monte Carlo and spin dynamics simulations have been used to study the dynamic critical behavior of ${\mathrm{RbMnF}}_{3},$ treated as a classical Heisenberg antiferromagnet on a simple cubic lattice. In an attempt to understand the difference in the value of the dynamic critical exponent z between experiment and theory, we have used larger lattice sizes than in our previous simulations to better probe the asymptotic critical region in momentum. We estimate $z=1.49\ifmmode\pm\else\textpm\fi{}0.03,$ in good agreement with the renormalization-group theory and dynamic scaling predictions. In addition, the central peak in the dynamic structure factor at ${T}_{c},$ seen in experiments and previous simulations, but absent in the renormalization-group and mode-coupling theories, is shown to be solely in the longitudinal component.

DOI
10.1103/physrevb.67.104411

Related papers

Spin dynamics simulations of excitations and critical dynamics in RbMnF_3

D. P. Landau, Shan-Ho Tsai, Alex Bunker

Spin-dynamics simulations have been used to investigate the dynamic behavior of RbMnF_3, treating it as a classical Heisenberg antiferromagnet on a simple cubic lattice. Time-evolutions of spin configurations were determined numerically from coupled equations of motion for individual spins using a new algorithm which is based on Suzuki-Trotter decompositions of exponential operators. The dynamic structure factor was calculated from the space- and time-displaced spin-spin correlation function. The crossover from hydrodynamic to critical behavior of the dispersion curve and spin-wave half-width was studied as the temperature was increased towards the critical value. The dynamic critical exponent was estimated to be z=(1.43\pm 0.03), which is slightly lower than the dynamic scaling prediction, but in good agreement with a recent experimental value. Comparisons are made of both the dispersion curve and the lineshapes obtained from our simulations with very recent experimental results for RbMnF_3 are presented.

arXiv:cond-mat/9912374•Dec 20, 1999•0 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

Spin-dynamics simulations of the magnetic dynamics of RbMnF3 and direct comparison with experiment

Shan-Ho Tsai, Alex Bunker, David Paul Landau

Spin-dynamics techniques have been used to perform large-scale simulations of the dynamic behavior of the classical Heisenberg antiferromagnet in simple cubic lattices with linear sizes $L<~60.$ This system is widely recognized as an appropriate model for the magnetic properties of ${\mathrm{RbMnF}}_{3}.$ Time evolutions of spin configurations were determined numerically from coupled equations of motion for individual spins using an algorithm implemented by Krech et al., which is based on fourth-order Suzuki-Trotter decompositions of exponential operators. The dynamic structure factor was calculated from the space- and time-displaced spin-spin correlation function. The crossover from hydrodynamic to critical behavior of the dispersion curve and spin-wave half-width was studied as the temperature was increased towards the critical temperature. The dynamic critical exponent was estimated to be $z=(1.43\ifmmode\pm\else\textpm\fi{}0.03),$ which is slightly lower than the dynamic scaling prediction, but in good agreement with a recent experimental value. Direct, quantitative comparisons of both the dispersion curve and the line shapes obtained from our simulations with very recent experimental results for ${\mathrm{RbMnF}}_{3}$ are presented.

arXiv:cond-mat/9910429•Physical review. B, Condensed matter•Jan 1, 2000•28 citations

CRITICAL DYNAMICS IN THE 3D SPIN-1/2 HEISENBERG MODEL: A DECOUPLED CELL MONTE CARLO STUDY

Cynthia J. Sisson

The three-dimensional spin-1/2 Heisenberg model on a simple cubic lattice is studied for ferromagnetic and antiferromagnetic interactions using the Decoupled Cell Method for quantum Monte Carlo. Results for the relaxation time τL are determined for both ferromagnetic and antiferromagnetic systems and found to be similar to those found for the classical (s → ∞) Heisenberg model. The scaling of τL with system size is used to extract the dynamical critical exponent z for the two systems. The values of z = 1.98 ± 0.12 for the ferromagnet and z = 1.94 ± 0.09 for the antiferromagnet are in good agreement with theoretical predictions and previous Monte Carlo studies of the classical Heisenberg model.

International Journal of Modern Physics C•Jun 1, 1996•1 citation

High-accuracy Monte Carlo study of the three-dimensional classical Heisenberg ferromagnet

Pawel Peczak, Alan M. Ferrenberg, D. P. Landau

Using extensive Monte Carlo simulations, we study the equilibrium properties of the simple-cubic, classical Heisenberg ferromagnet. We employ very long runs for L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L lattices to obtain high-precision data for the magnetization probability distribution. Using finite-size scaling for L\ensuremath{\le}24 and an optimized multiple-histogram data analysis, we obtain an accurate value of the inverse critical temperature J/${\mathit{k}}_{\mathit{B}}$${\mathit{T}}_{\mathit{c}}$=0.6929\ifmmode\pm\else\textpm\fi{}0.0001, which is higher than previously accepted estimates. Calculated values of various static exponents are in excellent agreement with renormalization-group and \ensuremath{\epsilon}-expansion predictions.

Physical review. B, Condensed matter•Mar 1, 1991•286 citations

Spin-dynamics study of the dynamic critical behavior of the three-dimensional classical Heisenberg ferromagnet

Kun Chen, David Paul Landau

Using spin-dynamics techniques we have performed large-scale computer 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 ferromagnet with L\ensuremath{\le}40 in the vicinity of the critical point ${\mathit{T}}_{\mathit{c}}$. The temporal evolutions of 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 methods. The space- and time-displaced spin-spin correlation functions and their space-time Fourier transforms were calculated to determine the neutron-scattering functions. We developed a dynamic finite-size scaling theory for the neutron-scattering function at ${\mathit{T}}_{\mathit{c}}$ and used it to extract the dynamic critical exponent z. Within our resolution limit, the value of z was estimated to be 2.478(28), in excellent agreement with the dynamic scaling prediction, and dynamic scaling was found to be valid for momentum transfer q up to 0.4\ensuremath{\pi} in the (1,0,0) direction.

Physical review. B, Condensed matter•Feb 1, 1994•88 citations
Monte Carlo and spin dynamics simulations have been used to study the dynamic critical behavior of RbMnF3,{\mathrm{RbMnF}}_{3},RbMnF3​, treated as a classical Heisenberg antiferromagnet on a simple cubic lattice. In an attempt to understand the difference in the value of the dynamic critical exponent z between experiment and theory, we have used larger lattice sizes than in our previous simulations to better probe the asymptotic critical region in momentum. We estimate z=1.49\ifmmode±\else\textpm\fi0.03,z=1.49\ifmmode\pm\else\textpm\fi{}0.03,z=1.49\ifmmode±\else\textpm\fi0.03, in good agreement with the renormalization-group theory and dynamic scaling predictions. In addition, the central peak in the dynamic structure factor at Tc,{T}_{c},Tc​, seen in experiments and previous simulations, but absent in the renormalization-group and mode-coupling theories, is shown to be solely in the longitudinal component.
Spin-dynamics simulations have been used to investigate the dynamic behavior of RbMnF_3, treating it as a classical Heisenberg antiferromagnet on a simple cubic lattice. Time-evolutions of spin configurations were determined numerically from coupled equations of motion for individual spins using a new algorithm which is based on Suzuki-Trotter decompositions of exponential operators. The dynamic structure factor was calculated from the space- and time-displaced spin-spin correlation function. The crossover from hydrodynamic to critical behavior of the dispersion curve and spin-wave half-width was studied as the temperature was increased towards the critical value. The dynamic critical exponent was estimated to be z=(1.43\pm 0.03), which is slightly lower than the dynamic scaling prediction, but in good agreement with a recent experimental value. Comparisons are made of both the dispersion curve and the lineshapes obtained from our simulations with very recent experimental results for RbMnF_3 are presented.
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.
Spin-dynamics techniques have been used to perform large-scale simulations of the dynamic behavior of the classical Heisenberg antiferromagnet in simple cubic lattices with linear sizes L< 60.L<~60.L< 60. This system is widely recognized as an appropriate model for the magnetic properties of RbMnF3.{\mathrm{RbMnF}}_{3}.RbMnF3​. Time evolutions of spin configurations were determined numerically from coupled equations of motion for individual spins using an algorithm implemented by Krech et al., which is based on fourth-order Suzuki-Trotter decompositions of exponential operators. The dynamic structure factor was calculated from the space- and time-displaced spin-spin correlation function. The crossover from hydrodynamic to critical behavior of the dispersion curve and spin-wave half-width was studied as the temperature was increased towards the critical temperature. The dynamic critical exponent was estimated to be $z=(1.43\ifmmode\pm\else\textpm\fi{}0.03),$ which is slightly lower than the dynamic scaling prediction, but in good agreement with a recent experimental value. Direct, quantitative comparisons of both the dispersion curve and the line shapes obtained from our simulations with very recent experimental results for RbMnF3{\mathrm{RbMnF}}_{3}RbMnF3​ are presented.
The three-dimensional spin-1/2 Heisenberg model on a simple cubic lattice is studied for ferromagnetic and antiferromagnetic interactions using the Decoupled Cell Method for quantum Monte Carlo. Results for the relaxation time τL are determined for both ferromagnetic and antiferromagnetic systems and found to be similar to those found for the classical (s → ∞) Heisenberg model. The scaling of τL with system size is used to extract the dynamical critical exponent z for the two systems. The values of z = 1.98 ± 0.12 for the ferromagnet and z = 1.94 ± 0.09 for the antiferromagnet are in good agreement with theoretical predictions and previous Monte Carlo studies of the classical Heisenberg model.
Using extensive Monte Carlo simulations, we study the equilibrium properties of the simple-cubic, classical Heisenberg ferromagnet. We employ very long runs for L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L lattices to obtain high-precision data for the magnetization probability distribution. Using finite-size scaling for L\ensuremath{\le}24 and an optimized multiple-histogram data analysis, we obtain an accurate value of the inverse critical temperature J/kB{\mathit{k}}_{\mathit{B}}kB​Tc{\mathit{T}}_{\mathit{c}}Tc​=0.6929\ifmmode\pm\else\textpm\fi{}0.0001, which is higher than previously accepted estimates. Calculated values of various static exponents are in excellent agreement with renormalization-group and \ensuremath{\epsilon}-expansion predictions.
Using spin-dynamics techniques we have performed large-scale computer 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 ferromagnet with L\ensuremath{\le}40 in the vicinity of the critical point Tc{\mathit{T}}_{\mathit{c}}Tc​. The temporal evolutions of 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 methods. The space- and time-displaced spin-spin correlation functions and their space-time Fourier transforms were calculated to determine the neutron-scattering functions. We developed a dynamic finite-size scaling theory for the neutron-scattering function at Tc{\mathit{T}}_{\mathit{c}}Tc​ and used it to extract the dynamic critical exponent z. Within our resolution limit, the value of z was estimated to be 2.478(28), in excellent agreement with the dynamic scaling prediction, and dynamic scaling was found to be valid for momentum transfer q up to 0.4\ensuremath{\pi} in the (1,0,0) direction.