We report the results of Monte Carlo simulations on several spin glass models at low temperatures. By using the parallel tempering (Exchange Monte Carlo) technique we are able to equilibrate down to low temperatures, for moderate sizes, and hence the data should not be affected by critical fluctuations. Our results for short range models are consistent with a picture proposed earlier that there are large scale excitations which cost only a finite energy in the thermodynamic limit, and these excitations have a surface whose fractal dimension is less than the space dimension. For the infinite range Viana-Bray model, our results obtained for a similar number of spins are consistent with standard replica symmetry breaking.
We investigate the balanced
M=4,
p=4 spin-glass model for a one-dimensional long-range proxy for the finite dimensional short-range
p-spin glass model to examine the nature of the glass transition beyond mean-field theory. We perform large-scale Monte Carlo equilibrated simulations for both fully connected and power-law diluted versions of the model. The critical temperatures extracted from the finite-size scaling (FSS) analysis of spin-glass susceptibility are in good agreement with theoretical predictions for
σ=0,0.25, and 0.55. For these values of the long-range exponent
σ (which is the power of the decrease of the interactions between the spins with their separation), one might have expected that mean-field theory would provide a good description of the system. However, the spin-overlap distribution and the value of the
λ-parameter do not provide numerical evidence for a one-step replica symmetry breaking (1RSB) phase transition. Instead, our results indicate a direct transition from the paramagnetic state to a full replica symmetry broken phase, with a renormalized value of
λ≡ω2/ω1<1 suggesting a continuous FRSB transition, despite this ratio being equal to 2 at mean-field level. A value of
λ>1 is required for the discontinuous 1RSB transition. We argue that strong finite-size effects and closely spaced transition temperatures remove the expected 1RSB transition for the system sizes which we can study. For values of the exponent
σ=0.85, which roughly corresponds to a three dimensional system, we find that the renormalized value of
λ is again less than 1, with no signs of either the 1RSB transition or the continuous FRSB transition, suggesting that the Kauzmann temperature
TK in three dimensions might be zero and the complete absence of phase transitions in structural glasses.
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.
Spin glasses have been the subject of intense study and considerable controversy for decades, and the low-temperature phase of short-range spin glasses is still poorly understood. Our main goal is to improve our understanding in this area and find an answer to the following question: Are there only a single pair or a countable infinity of pure states in the low temperature phase of the EA spin glass? To that aim we first start by introducing spin glasses and provide a brief history of their research, then proceed to describe our method of simulation, the parallel tempering Monte Carlo algorithm. Next, we present the results of a large-scale numerical study of the equilibrium three-dimensional Edwards-Anderson Ising spin glass with Gaussian disorder. In order to understand how the parallel tempering algorithm works, we measure various static, as well as dynamical quantities, such as the autocorrelation times and round-trip times for the parallel tempering Monte Carlo method. We examine the correlation between static and dynamic observables for _ 5000 disorder realizations and up to 1000 spins down to temperatures at 20% of the critical temperature, and our results show that autocorrelation times are directly correlated with the roughness of the free-energy landscape. In the following chapters, the three- and four-dimensional Edwards-Anderson and mean-field Sherrington-Kirkpatrick Ising spin glasses are studied again via large-scale Monte Carlo simulations at low temperatures, deep within the spinglass phase. Performing a careful statistical analysis of several thousand independent disorder realizations and using an observable that detects peaks in the overlap distribution, we show that the Sherrington-Kirkpatrick and Edwards-Anderson models have a distinctly different low-temperature behavior. We arrive to the following conclusion: The structure of the spin-glass overlap distribution for the Edwards-Anderson model suggests that its low-temperature phase has only a single pair of pure states. Finally we present results for several new observables, along with a few preliminary studies and suggestions for future research.
Spin glasses are the paradigm of complex systems. These materials present really slow dynamics. However, the nature of the spin glass phase in finite dimensional systems is still controversial. Different theories describing the low temperature phase have been proposed: droplet, replica symmetry breaking and chaotic pairs. We present analytical studies of critical properties of spin glasses, in particular, critical exponents at and below the phase transition, existence of a phase transition in a magnetic field, computation of the lower critical dimension (in presence/absence of a magnetic field). We also introduce some rigorous results based on the concept of metastate. Finally, we report some numerical results regarding the construction of the Aizenman-Wehr metastate, scaling of the correlation functions in the spin glass phase and existence of a phase transition in a field, confronting these results with the predictions of different theories.