The time-dependent behavior of the classical, ferromagnetic xy chain in a symmetry-breaking magnetic field was studied using a very fast, vectorized spin-dynamics method. The equations of motion were integrated using starting configurations determined by Monte Carlo simulations. By calculating spin-spin correlation functions and taking space-time Fourier transforms, we determined S(q,\ensuremath{\omega}) for a wide range of fields and temperatures. We identify and measure the characteristics of single spin-wave peaks and two spin-wave sum peaks, as well as central peaks due to both soliton and two spin-wave difference processes. These results are discussed in light of existing theory and experiments.
We have studied the time-dependent behavior of the classical XY chain in a symmetry breaking magnetic field using an ultrafast spin-dynamics method in which we integrate the equations of motion using starting configurations determined by Monte Carlo simulations. We calculated the three polarizations of S(q,ω) separately for temperatures from 0.1 to 1.0 J/kB and fields up to 1.5 J. The soliton contribution, the single-magnon peaks, and both the sum and difference two-magnon contributions are separated; we discuss them in detail and compare them with theoretical predictions.
We have investigated the time-dependent behavior of the classical XY chain in a symmetry-breaking anisotropy field using an ultrafast spin-dynamics method in which we integrate the equation of motion with equilibrium spin configurations determined by Monte Carlo simulations. All three polarizations of the scattering function S(q,\ensuremath{\omega}) were determined for temperature
kBT=0.2J and fields varying from 0.1J to 0.6J. Substantial differences are found as compared with the isotropic ferromagnetic chain in a transverse magnetic field. The soliton contribution, single-spin-wave peaks, and both the sum and difference two-spin-wave contributions are separated. Half-widths and integrated intensities are extracted and compared with theoretical predictions based on mapping onto the sine-Gordon model.
The classical one-dimensional XY model in a symmetry-breaking Gaussian-distributed random magnetic field was studied using an ultrafast, vectorized spin-dynamics program on a Cray Y-MP. We calculated the time- and space-displaced spin-spin correlation functions, which was then Fourier-transformed to get the neutron-scattering law S(q,\ensuremath{\omega}). We see a clear change of the different contributions when the random field is switched on. At low temperatures the random field induces significant increases in the soliton density.
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.
We present results for the static and dynamic behavior of a one-dimensional classical
XY chain. Equilibrium configurations are produced by an importance-sampling Monte Carlo method. Static properties are evaluated from these data and a spin-dynamics method is then used to evaluate time-displaced correlation functions and to calculate the scattering function $S(q,\ensuremath{\omega})$. The transverse spin-spin correlation functions decay exponentially for large times. The longitudinal spin-spin correlation functions show oscillatory behavior, especially at low temperatures. The longitudinal-scattering function displays spin-wave peaks, and the transverse-scattering function shows central peaks (relaxational behavior) which change to fluctuation-broadened spin-wave peaks for large values of $\frac{q}{\ensuremath{\kappa}}$. The energy-energy correlation function shows diffusive behavior for high temperatures and non-diffusive behavior for low temperatures.
Using Monte Carlo and spin-dynamics methods, we have investigated the dynamic behavior of the classical, antiferromagnetic XY model on a triangular lattice with linear sizes
L< 300. The temporal evolutions of spin configurations were obtained by solving numerically the coupled equations of motion for each spin using fourth-order Suzuki-Trotter decompositions of exponential operators. From space- and time-displaced spin-spin correlation functions and their space-time Fourier transforms we obtained the dynamic structure factor
S(q,w) for momentum
q and frequency $\ensuremath{\omega}.$ Below
TKT (Kosterlitz-Thouless transition), both the in-plane
(Sxx) and out-of-plane
(Szz) components of $S(\mathbf{q},\ensuremath{\omega})$ exhibit very strong and sharp spin-wave peaks. Well above
TKT, Sxx and
Szz apparently display a central peak, and spin-wave signatures are still seen in
Szz. In addition, we also observed an almost dispersionless domain-wall peak at high $\ensuremath{\omega}$ below
Tc (Ising transition), where long-range order appears in the staggered chirality. Above
Tc, the domain-wall peak disappears for all q. The line shape of these peaks is captured reasonably well by a Lorentzian form. Using a dynamic finite-size scaling theory, we determined the dynamic critical exponent
z=1.002(3). We found that our results demonstrate the consistency of the dynamic finite-size scaling theory for the characteristic frequency ${\ensuremath{\omega}}_{m}$ and the dynamic structure factor $S(\mathbf{q},\ensuremath{\omega})$ itself.