We propose a change of style for numerical estimations of physical quantities from measurements to inferences. We estimate the most probable quantities for all the parameter region simultaneously by using the raw data cooperatively. Estimations with higher precisions are made possible. We can obtain a physical quantity as a continuous function, which is processed to obtain another quantity. We applied the method to the Heisenberg spin-glass model in three dimensions. A dynamic correlation-length scaling analysis suggests that the spin-glass and the chiral-glass transitions occur at the same temperature with a common exponent ν. The value is consistent with the experimental results. We explained a spin-chirality separation problem by a size-crossover effect.
A spin-glass transition has been investigated for a long time but we have not reached a conclusion yet due to difficulties in the simulation studies. They are slow dynamics, strong finite-size effects, and sample-to-sample dependencies. We found that a size of the spin-glass order reaches a lattice boundary within a very short Monte Carlo step. A competition between the spin-glass order and a boundary condition causes these difficulties. Once the boundary effect was removed, physical quantities exhibited quite normal behaviors. They became self-averaging in a limit of large replica numbers. These findings suggest that the nonequilibrium relaxation method is a good choice for solving the difficulties if a lattice size and a replica number are set sufficiently large. A dynamic scaling analysis on nonequilibrium relaxation functions gave a result that the spin-glass transition and the chiral-glass transition occurs at the same temperature in the Heisenberg model in three dimensions. The estimated critical exponent ν agrees with the experimental result.
The finite-size-scaling method in the equilibrium Monte Carlo (MC) simulations and the finite-time-scaling method in the nonequilibrium relaxation simulations are compromised. MC time data of various physical quantities are scaled by the MC time data of the dynamic correlation length. It corresponds to changing the system size in the finite-size-scaling method. This scaling method is tested in the three-dimensional ferromagnetic Ising model and in the three-dimensional $\ifmmode\pm\else\textpm\fi{}J$ Ising spin-glass model. The transition temperature and the critical exponents are obtained. We also comment on the definition of the dynamic correlation length in the nonequilibrium relaxation process. The Ornstein-Zernike formula is not always appropriate.
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 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.
Using single cluster flip Monte Carlo simulations we accurately determine new finite size scaling functions which are expressed only in terms of the variable x=${\ensuremath{\xi}}_{\mathit{L}}$/L, where ${\ensuremath{\xi}}_{\mathit{L}}$ is the correlation length in a finite system of size L. Data for the d=2 and d=3 Ising models, taken at different temperatures and for different size lattices, show excellent data collapse over the entire range of scaling variable for susceptibility and correlation length. From these finite size scaling functions we can estimate critical temperatures and exponents with rather high accuracy even though data are not obtained extremely close to the critical point. The bulk values of the renormalized four-point coupling constant are accurately measured and show strong evidence for hyperscaling. \textcopyright{} 1996 The American Physical Society.