Explicit dynamic finite-size scaling in the time-dependent correlation functions for the two-dimensional classical XY model is explored. The dynamic scaling method, proposed by us previously, is utilized for finding the time scaling variables of the model. It shows nontrivial dynamic finite-size scaling behavior near and below TKT for the total spin correlation function. Since the model is critical for all temperatures below TKT, the time scaling exponent is markedly different from that of the three-dimensional classical XY model at short times. There also appears to be another dynamic scaling region with a different time scaling exponent at T⩽TKT for a fairly wide intermediate time interval.
As a further test of the dynamic scaling method proposed by us previously, dynamic finite-size effect in an infinite-range classical antiferromagnetic spin model, which is an anisotropic version of two-sublattice Lieb-Mattis model, is studied. Application of our dynamic scaling method to the model yields time scaling exponents similar to those of the infinite-range ferromagnetic classical spin van der Waals model. With the hereby obtained time scaling variables, we can successfully collapse the time-dependent correlation functions of the model into universal curves.
A phenomenological dynamic scaling method is proposed and utilized to interpret finite-size effects in the time-dependent correlation function for the classical spin van der Waals model. We claim that the equations of motion for the classical spin van der Waals model are well defined only within a certain characteristic time. The order of magnitude for the size-dependent maximum time interval of the time-dependent correlation function, analytically obtained in the large N limit via the method of Laplace by Dekeyser and Lee, is found to be related to our characteristic time. As a by-product of this work, time scaling variables are found with which we can completely collapse the time-dependent correlation functions of the classical spin van der Waals model into one curve. \textcopyright{} 1996 The American Physical Society.
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.
Computer calculations of the time-dependent spin and energy correlation functions of the classical one-dimensional
XY model at infinite temperature are reported. Plots of
〈SziSzi(t)〉,
〈SziSzi+1(t)〉, $〈{\ensuremath{\epsilon}}^{i}{\ensuremath{\epsilon}}^{i}(t)〉$, and
(21)[〈SxiSxi(t)〉+〈SyiSyi(t)〉] out to times
Jt=9 are given (${\ensuremath{\epsilon}}^{i}$ is the energy-density operator and $\ensuremath{\hbar}J$ is the exchange integral). Comparison is made with the exact results for the spin-1/2
XY model. After an appropriate scaling of the exchange integral and normalization to the classical value at
t=0 the values of the spin-1/2 functions are close to the values of their classical counterparts for times up to $\mathrm{Jt}\ensuremath{\approx}4$.
A scaling hypothesis is set up for the correlation function G(R,T;L) of a finite-sized system, with O(n) symmetry (n\ensuremath{\ge}2), confined to geometry ${L}^{d\mathrm{\ensuremath{-}}d\mathcal{'}}$\ifmmode\times\else\texttimes\fi{}${\ensuremath{\infty}}^{d\mathcal{'}}$ (2d4, d'\ensuremath{\le}2) and subjected to periodic boundary conditions. Finite-size effects in the correlation length \ensuremath{\xi}(T;L) are predicted in the region of first-order phase transition (T
Tc) as well as the region of second-order phase transition (T\ensuremath{\simeq}
Tc); for the correlation function, the predictions relate mostly to the regime T
Tc. To test these predictions, a detailed analytical study is carried out in the case of the spherical model and all predictions are seen to be fully borne out. In the process we find that, though a single correlation length \ensuremath{\xi} governs the scaling of the variable R in all directions, the functional dependence of the correlation function on the scaled parameter R/\ensuremath{\xi} is highly anisotropic, in that it depends very differently on the components, ${\mathrm{R}}_{\ensuremath{\perp}}$ and
R?, of R pertaining to the finite and the infinite dimensions, respectively. Most importantly, while long-range order at temperatures below
Tc prevails all the way in the direction of ${\mathrm{R}}_{\ensuremath{\perp}}$ it is severely limited in the direction of
R?. In fact, for
R?\ensuremath{\gg}L, the qualitative features of the correlation function become characteristic of a d'-dimensional bulk system, though \ensuremath{\xi} still pertains to the actual, finite-sized system; the net result is that long-range order in the direction of
R? persists only to distances small in comparison with \ensuremath{\xi} which, for d'2, is found to be O(L(L/a${)}^{(d\mathrm{\ensuremath{-}}2)/(2\mathrm{\ensuremath{-}}d\mathcal{'})}$), a being the lattice constant.