Spin-glass (SG) and chiral-glass (CG) orderings in three-dimensional (3D) Heisenberg spin glass with and without magnetic anisotropy are studied by using large-scale off-equilibrium Monte Carlo simulations. A characteristic time of relaxation, which diverges at a transition temperature in the thermodynamic limit, is obtained as a function of the temperature and the system size. Based on the finite-size-scaling analysis for the relaxation time, it is found that in the isotropic Heisenberg spin glass, the CG phase transition occurs at a finite temperature, while the SG transition occurs at a lower temperature, which is compatible with zero. Our results of the anisotropic case support the chirality scenario for the phase transitions in the 3D Heisenberg spin glasses.
An extensive equilibrium Monte Carlo simulation is performed on the three-dimensional isotropic Heisenberg spin-glass (SG) model with the random nearest-neighbor Gaussian coupling, with particular interest in its chiral-glass (CG) and SG orderings. For this model, the possibility of the spin-chirality decoupling, i.e., the CG order setting in at a higher temperature than that of the SG order, was suggested earlier but still remains controversial. We simulate the model up to the maximum size (linear dimension)
L=48 under both periodic and open boundary conditions (BCs). In locating the CG and SG transition temperatures
TCG and
TSG with the $L\ensuremath{\rightarrow}\ensuremath{\infty}$ extrapolation, a variety of independent physical quantities under both BCs are computed and utilized to get larger number of degrees of freedom (NDF). Because of the large NDF up to
NDF=43, we succeed in obtaining stable and accurate estimates of the CG and SG transition temperatures, ${T}_{\mathrm{CG}}=0.142\ifmmode\pm\else\textpm\fi{}0.001$ and ${T}_{\mathrm{SG}}=0.{131}_{\ensuremath{-}0.006}^{+0.001}$. No sign of the size crossover is observed. For larger
L, the CG correlation length progressively outgrows the SG correlation length at low temperatures. These results provide strong numerical support for the spin-chirality decoupling. The critical exponents associated with the CG and SG transitions are evaluated using the finite-size scaling with the scaling correction. For the CG transition, we get the CG exponents, ${\ensuremath{\nu}}_{\mathrm{CG}}=1.36\ifmmode\pm\else\textpm\fi{}0.10$ and ${\ensuremath{\eta}}_{\mathrm{CG}}=0.49\ifmmode\pm\else\textpm\fi{}0.10$, consistent with the corresponding experimental exponents of canonical SG. Implications for the chirality scenario of experimental SG ordering are discussed.
The statics-dynamics correspondence in spin glasses relate nonequilibrium results on large samples (the experimental realm) with equilibrium quantities computed on small systems (the typical arena for theoretical computations). Here we employ statics-dynamics equivalence to study the Ising spin-glass critical behavior in three dimensions. By means of Monte Carlo simulation, we follow the growth of the coherence length (the size of the glassy domains), on lattices too large to be thermalized. Thanks to the large coherence lengths we reach, we are able to obtain accurate results in excellent agreement with the best available equilibrium computations. To do so, we need to clarify the several physical meanings of the dynamic exponent close to the critical temperature.
The statics-dynamics correspondence in spin glasses relate non-equilibrium results on large samples (the experimental realm) with equilibrium quantities computed on small systems (the typical arena for theoretical computations). Here we employ statics-dynamics equivalence to study the Ising spin-glass critical behavior in three dimensions. By means of Monte Carlo simulation, we follow the growth of the coherence length (the size of the glassy domains), on lattices too large to be thermalized. Thanks to the large coherence lengths we reach, we are able to obtain accurate results in excellent agreement with the best available equilibrium computations. To do so, we need to clarify the several physical meanings of the dynamic exponent close to the critical temperature.
It is shown, by means of Monte Carlo simulation and finite size scaling analysis, that the Heisenberg spin glass undergoes a finite-temperature phase transition in three dimensions. There is a single critical temperature, at which both a spin glass and a chiral glass ordering develop. The Monte Carlo algorithm, adapted from lattice gauge theory simulations, makes it possible to thermalize lattices of size L = 32, larger than in any previous spin-glass simulation in three dimensions. High accuracy is reached thanks to the use of the Marenostrum supercomputer. The large range of system sizes studied allows us to consider scaling corrections.
As an analysis of equilibrium phase transitions, the nonequilibrium relaxation method is extended to the spin glass (SG) transition. The $\ifmmode\pm\else\textpm\fi{}J$ Ising SG model is analyzed for three-dimensional (cubic) lattices up to the linear size of
L=127 and for four-dimensional (hypercubic) lattice up to
L=41. These sizes of systems are quite large as compared with those calculated, so far, by equilibrium simulations. As a dynamical order parameter, we calculate the clone correlation function (CCF) ${Q(t,t}_{\mathrm{w}})\ensuremath{\equiv}[〈{S}_{i}^{(1)}{(t+t}_{\mathrm{w}}){S}_{i}^{(2)}{(t+t}_{\mathrm{w}}){〉}^{\mathrm{F}}],$ which is a spin correlation of two replicas produced after the waiting time
tw from a simple starting state. It is found that the CCF shows an exponential decay in the paramagnetic phase, and a power-law decay after aginglike development $(t\ensuremath{\gg}{t}_{\mathrm{w}})$ in the SG phase. This provides a reliable upper bound of the transition temperature
Tg. It is also found that a scaling relation, ${Q(t,t}_{\mathrm{w}}{)=t}_{\mathrm{w}}^{\ensuremath{-}{\ensuremath{\lambda}}_{q}}\overline{q}{(t/t}_{\mathrm{w}}),$ holds just around the transition point providing the lower bound of
Tg. Together with these two bounds, we propose a new dynamical way for the estimation of
Tg from much larger systems. In the SG phase, the power-law behavior of the CCF for $t\ensuremath{\gg}{t}_{\mathrm{w}}$ suggests that the SG phase in short-range Ising models has a rugged phase space.