The spin dynamics of two dimensional XY ferromagnet are reexamined at low temperatures in the framework of the Mori continued fraction formalism using the Gaussian approximation. In this formalism, the terms on denominators for the Laplace transform of the relaxation function R(t) are related to the frequency moments of the relaxation shape function R(q,w). In the Gaussian approximation scheme, we truncate the continued fraction for R(t) on the second stage. Adopting this approximation, we calculated up to the sixth moment. The moments are writing in terms the static spin correlation functions. In the estimation of the fourth and sixth moment at finite temperature, the four and six-spin correlation functions may be approximated by a sum of products of appropriate pair correlation functions(mode-mode decoupling). In this work we calculated the static spin correlation functions for the expressions to the fourth and sixth moments, needed in the study of the dynamics.
The spin dynamics of two dimensional XY ferromagnet are reexamined at low temperatures in the framework of the Mori continued fraction formalism using the Gaussian approximation. In this formalism, the terms on denominators for the Laplace transform of the relaxation function R(t) are related to the frequency moments of the relaxation shape function R(q,w). In the Gaussian approximation scheme, we truncate the continued fraction for R(t) on the second stage. Adopting this approximation, we calculated up to the sixth moment. The moments are writing in terms the static spin correlation functions. In the estimation of the fourth and sixth moment at finite temperature, the four and six-spin correlation functions may be approximated by a sum of products of appropriate pair correlation functions(mode-mode decoupling). In this work we calculated the static spin correlation functions for the expressions to the fourth and sixth moments, needed in the study of the dynamics.
The static two- and four-spin correlations for the one-dimensional easy-plane Heisenberg spin system with a ferromagnetic coupling, like CsNiF 3 , are investigated in the classical limit by means of the transfer matrix method. They are found to consist essentially of two components which both depend exponentially on the distance between sites. Both in the high- and low-temperature regions, only one component remains and this component still dominates at intermediate temperatures. Analytic forms for the correlation functions are given for asymptotically low temperatures. These are obtained by calculating the second-frequency moment of the dynamical structure factor in three different ways, and by requiring the T=0 limit of the moment to coincide with the square of the spin wave frequency. Generalisations of the Tomita-Mashiyama recursion relations (1972) for the isotropic case are given, valid for low temperatures, Monte Carlo calculations of the correlation functions have also been performed and some results of these are presented as well.
The exact short-time expansion coefficients of the two-time spin-pair correlation function σ(Rif,t)=≪siz(t)sfz(0)≫c at infinite temperature are used in order to estimate the spin diffusion constants of the diluted Heisenberg and XY magnets of spin 1/2. The used coefficients are those up to the term of t10 for the Heisenberg magnet on the linear chain, and those up to the term of t8 for the Heisenberg and XY magnets on the square, sc and bcc lattices. Our calculations are based on the “linear extrapolation approximation”, in which unknown coefficients in Mori's continued fraction representation for ~Is(k), the Laplace-time Fourier-space transform of σ(Rif,t), are linearly extrapolated from the last two of the known coefficients. The obtained values of the diffusion constant and their concentration dependences show a reasonable trend of convergence for successive approximations in a range of concentration not close to the percolation limit. Discussions are also given on the diffusion constants obtained by the “two-parameter Gaussian approximation in time” applied to the memory function.
The infinite-temperature space- and time-dependent spin-correlation functions grx(t) are studied for the one-dimensional XY model. Numerical calculations are performed to obtain the exact autocorrelation function g0x(t) for chains containing 5, 7, and 9 spins (S=21). This yields exact results for the first 16 moments of the frequency autocorrelation function of the infinite chain, and estimates for a few of the higher moments as well. The analysis suggests that g0x(t) for the infinite chain is identical to $\frac{\mathrm{exp}(\ensuremath{-}{J}^{2}{t}^{2})}{4}$. We show that grx(t) for $r\ensuremath{\ne}0$ vanishes identically for all values of time, implying a wave-vector-independent relaxation shape function. Our result for g0x(t) is compared with that obtained by Huber for the classical ($S=\ensuremath{\infty}$) chain.
Using spin-dynamics techniques we have performed large-scale computer simulations of the dynamic behavior of the classical three component XY model (i.e., the anisotropic limit of an easy-plane Heisenberg ferromagnet), on square lattices of size up to 1922, for several temperatures below, at, and above TKT. The temporal evolution of spin configurations was determined numerically from coupled equations of motion for individual spins by a fourth-order predictor-corrector method, with initial spin configurations generated by a hybrid Monte Carlo algorithm. The neutron-scattering function S(q,\ensuremath{\omega}) was calculated from the resultant space-time displaced spin-spin correlation function. Pronounced spin-wave peaks were found both in the in-plane and the out-of-plane scattering function over a wide range of temperatures. The in-plane scattering function Sxx also has a large number of clear but weak additional peaks, which we interpret to come from two-spin-wave scattering. In addition, we observed a small central peak in Sxx, even at temperatures well below the phase transition. We used dynamic finite-size scaling theory to extract the dynamic critical exponent z. We find z=1.00(4) for all T\ensuremath{\le}TKT, in excellent agreement with theoretical predictions, although the shape of S(q,\ensuremath{\omega}) is not well described by current theory. \textcopyright{} 1996 The American Physical Society.
The two-time spin-pair correlation function σ(R,t) is expressed as a product of a Gaussian distribution function and a power series with respect to time. The width of the Gaussian distribution function is determined by the second derivative with respect to time of the autocorrelation function σ(0,t). The coefficients of the power series are determined up to terms of order t8 for the square and s.c. isotropic Heisenberg and XY magnets of spin ½ at infinite temperature, and up to terms of t10 for the linear magnets. The expressions obtained by truncating the power series to the exactly known terms give a very good fit at short times to the exact expression for the linear XY model, and to the results of the computer simulation due to Windsor for the s.c. Heisenberg ferromagnet. The Fourier-time transforms of the expressions thus obtained are shown to be the Gram-Charlier expansions. Thus we conclude that the Gram-Charlier expansion of the Fourier-time transform of σ(R,t),σ̃(R,ω), and that of its Fourier-space transform I(k,t),S(k,ω), is especially useful in the cases when the short-time behavior of σ(R,t) or I(k,t) is considered important; generally speaking, this is the case at large k or at ω which is not very small. When the width of the Gaussian distribution function is determined by the second derivative of I(k,t), a convergent result is not obtained for small values of k. For larger k, the convergence is as good as for the case when the width is determined by the autocorrelation function.
The dynamical correlation function Gr(t)==〈Sr(t)\ensuremath{\cdot}S0(0)〉 and dynamical structure factor Sq(\ensuremath{\omega}) are calculated by the modified spin-wave theory for the low-dimensional quantum Heisenberg ferromagnets at low temperature. We use the Dyson-Maleev transformation, ideal spin-wave states, and the rotational averaging. Sq(\ensuremath{\omega}) satisfies the dynamic scaling relation. The explicit form of the scaling function is obtained. The classical limit of our results is compared with a molecular-dynamics calculation. The agreement is surprisingly good.