The low temperature dynamics of the three-dimensional Ising spin-glass in zero field with a discrete bond distribution is investigated via MC simulations. The thermoremanent magnetization is found to decay algebraically and the temperature dependent exponents agree very well with the experimentally determined values. The nonequilibrium autocorrelation function C(t, t w ) shows a crossover at the waiting (or aging) time t w from algebraic quasiequilibrium decay for times t<<t w to another, faster algebraic decay for t>>t w with an exponent similar to one for the remanent magnetization.
The crossover between quasiequilibrium and nonequilibrium spin-glass dynamics has been studied in Monte Carlo simulations of two- and three-dimensional short-range Ising spin-glass systems. The spin system was quenched in zero field to a low temperature T. After equilibrating for a time
tw, a weak magnetic field was applied and the time dependence of the magnetization M(t) and the spin autocorrelation function q(t) was studied. Our results show that the relaxation rate S(t)=\ensuremath{\partial}M(t)/\ensuremath{\partial}
log10(t) exhibits a maximum at t\ensuremath{\approxeq}
tw. If the temperature is lowered or raised, immediately prior to the field application, the apparent age (
ta) of the system, as defined by the maximum in relaxation rate, is shifted to longer or shorter time scales, respectively. No significant differences between the aging behavior in two and three dimensions are found. The relation \ensuremath{\chi}(t)=[1-q(t)]/T holds only for short time scales (t\ensuremath{\ll}
tw), clear deviations are found at longer time scales (t>
tw). Since this relation should according to the fluctuation-dissipation theorem hold for a system in equilibrium, this result is interpreted as an indication of a crossover between equilibrium and nonequilibrium dynamics. Our results are qualitatively in agreement with experiments as well as with theoretical models for the nonequilibrium spin-glass dynamics.
Aging of two-dimensional (2D) Ising spin glasses with short-range interactions
Rb2${\mathrm{Cu}}_{1\mathrm{\ensuremath{-}}\mathit{x}}$
CoxF4 (
Tc=0) is studied in the critical regime via the low-frequency ac susceptibility and via the dc magnetization for x=0.22 and 0.33. After quenching to low temperatures, slow logarithmic relaxation is observed for both the in-phase and out-of-phase ac susceptibilities. Using the droplet scaling model, the spin-glass correlation length R is concluded to increase after a temperature quench as R\ensuremath{\propto}(lnt${)}^{1/\mathrm{\ensuremath{\psi}}}$, with t in units of some microscopic time. The exponent \ensuremath{\psi}=1.0\ifmmode\pm\else\textpm\fi{}0.1 compares with recent theoretical results and numerical simulations for 2D short-range Ising spin glasses. The waiting-time dependence found for the field-cooled and zero-field-cooled dc magnetizations points to a markedly slower crossover from the quasiequilibrium to the nonequilibrium regime than in 3D. From experiments with different cooling histories, finally, the equilibrium spin-glass state is inferred to be strongly dependent on the temperature.
I present an analysis of the dynamic behavior of short-range Ising spin glasses observed in stochastic simulations. The time dependence of the order parameter q(t)=〈
Sx(0)
Sx(t)〉\ifmmode\bar\else\textasciimacron\fi{}---which is the same as that of the structure factor---and the time dependence of the related dynamic correlation functions have been recorded with good statistics and very long observation times. The spin-glass model with a symmetric distribution of discrete nearest-neighbor \ifmmode\pm\else\textpm\fi{}J interactions on a simple-cubic lattice was used. Simulations were performed with a special fast computer, allowing for the first-time investigation of the equilibrium dynamics for a wide range of temperatures (0.7\ensuremath{\le}kT/J\ensuremath{\le}5.0) and lattice sizes (
83,
163,
323, and
643). I have found that the empirical formula q(t)=${\mathrm{ct}}^{\mathrm{\ensuremath{-}}x}$exp(-\ensuremath{\omega}${t}^{\ensuremath{\beta}}$) with temperature-dependent exponents x(T) and \ensuremath{\beta}(T) describes the decay very well at all temperatures above the spin-glass transition. In the spin-glass phase, only the algebraic decay q(t)=${\mathrm{ct}}^{\mathrm{\ensuremath{-}}x}$ could be observed, with different temperature dependences of the exponent x(T). The dynamic scaling hypothesis and finite-size scaling explain well the observed temperature and size dependence of the data, and the functional form of the correlation functions is com- patible with the scaling form if corrections to scaling are taken into account. The scaling behavior and the dynamic and static critical exponents found in my simulations are in reasonable agreement with recent experiments performed on insulating spin glasses, showing that despite its simplicity the discrete model of spin glasses analyzed in this work displays behavior similar to that seen in nature.
A quench of a d-dimensional spin system from a random initial configuration, {
Si(0)}, to a critical point is considered. The decay with time t of the autocorrelation with the initial condition is
q0(t)==〈
Si(0)\ensuremath{\cdot}
Si(t)〉\ensuremath{\sim}${t}_{\mathit{c}}^{\ensuremath{-}\ensuremath{\lambda}}$/z, where z is the usual dynamic critical exponent. Naively, ${\ensuremath{\lambda}}_{c}$=d, but I find ${\ensuremath{\lambda}}_{c}$simulations of pure Ising models in d=2, 3 and the \ifmmode\pm\else\textpm\fi{}J Ising spin glass in d=3. This suggests that ${\ensuremath{\lambda}}_{c}$ is a new critical exponent for nonequilibrium dynamics. For a spin glass the decay of
q0(t) is the same as that of the remanent magnetization; the exponent ${\ensuremath{\lambda}}_{c}$/z observed in the spin-glass simulation is in good agreement with a recent experimental measurement by Granberg et al.