Search

Using a recently implemented integration method [Krech et. al.] based on an iterative second-order Suzuki-Trotter decomposition scheme, we have performed spin dynamics simulations to study the critical dynamics of the BCC Heisenberg antiferromagnet with uniaxial anisotropy. This technique allowed us to probe the narrow asymptotic critical region of the model and estimate the dynamic critical exponent $z=2.25 \pm 0.08$. Comparisons with competing theories and experimental results are presented., Comment: Latex, 3 pages, 5 figures

We have developed a technique to accelerate the acquisition of effectively uncorrelated configurations for off-lattice models of dense polymer melts which makes use of both parallel tempering and large scale Monte Carlo moves. The method is based upon simulating a set of systems in parallel, each of which has a slightly different repulsive core potential, such that a thermodynamic path from full excluded volume to an ideal gas of random walks is generated. While each system is run with standard stochastic dynamics, resulting in an NVT ensemble, we implement the parallel tempering through stochastic swaps between the configurations of adjacent potentials, and the large scale Monte Carlo moves through attempted pivot and translation moves which reach a realistic acceptance probability as the limit of the ideal gas of random walks is approached. Compared to pure stochastic dynamics, this results in an increased efficiency even for a system of chains as short as $N = 60$ monomers, however at this chain length the large scale Monte Carlo moves were ineffective. For even longer chains the speedup becomes substantial, as observed from preliminary data for $N = 200$.

Spin-dynamics techniques can now be used to study the deterministic time-dependent behavior of magnetic systems containing over 10^5 spins with quite good accuracy. This approach will be introduced, including the theoretical foundations of the methods of analysis. Then newly developed, improved techniques based upon Suzuki-Trotter decomposition methods will be described. The current ``state-of-the-art'' will be evaluated with specific examples drawn from data on simple magnetic models. The examination of dynamic critical behavior will be highlighted but the extraction of information about excitations at low temperatures will be included., Comment: Prog. Theor. Phys. (in press)

Using Suzuki-Trotter decompositions of exponential operators we describe new algorithms for the numerical integration of the equations of motion for classical spin systems. These techniques conserve spin length exactly and, in special cases, also conserve the energy and maintain time reversibility. We investigate integration schemes of up to eighth order and show that these new algorithms can be used with much larger time steps than a well established predictor-corrector method. These methods may lead to a substantial speedup of spin dynamics simulations, however, the choice of which order method to use is not always straightforward., Comment: J. Mod. Phys. C (in press)

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., Comment: J. Phys. Soc. Japan (in press)

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\leq 60$. This system is widely recognized as an appropriate model for the magnetic properties of RbMnF$_3$. Time-evolutions of spin configurations were determined numerically from coupled equations of motion for individual spins using a new algorithm implemented by Krech {\it etal}, 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\pm 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 lineshapes obtained from our simulations with very recent experimental results for RbMnF$_3$ are presented., Comment: 30 pages, RevTex, 9 figures, to appear in PRB

Using spin dynamics simulations we predict the splitting of the longitudinal spin wave peak in all antiferromagnets with single site anisotropy into two peaks separated by twice the energy gap at the Brillouin zone center. This phenomenon has yet to be observed experimentally but can be easily investigated through neutron scattering experiments on MnF$_2$ and FeF$_2$. We have also determined that for all classical Heisenberg models the longitudinal propagative excitations are entirely multiple spin-wave in nature., Comment: four pages three figures, the last two postscript files are two parts of the third figure

We have proposed new algorithms for the numerical integration of the equations of motion for classical spin systems. In close analogy to symplectic integrators for Hamiltonian equations of motion used in Molecular Dynamics these algorithms are based on the Suzuki-Trotter decomposition of exponential operators and unlike more commonly used algorithms exactly conserve spin length and, in special cases, energy. Using higher order decompositions we investigate integration schemes of up to fourth order and compare them to a well established fourth order predictor-corrector method. We demonstrate that these methods can be used with much larger time steps than the predictor-corrector method and thus may lead to a substantial speedup of computer simulations of the dynamical behavior of magnetic materials., Comment: 9 pages RevTeX with 8 figures