We have studied the diluted Heisenberg spin glass model in a three-component random field for the commonly used one-dimensional long-range model where the probability that two spins separated by a distance r interact with one another falls as 1/r^{2σ}, for two values of σ, 0.75 and 0.85. No de Almeida-Thouless line is expected at these σ values. The spin glass correlation length ξ_{SG} varies with the random field as expected from the Imry-Ma argument and the droplet scaling picture of spin glasses. However, when ξ_{SG} becomes comparable to the system size L, there are departures which we attribute to the features deriving from the replica symmetry breaking picture of spin glasses. For the case σ=0.85 these features go away for system sizes with L>L^{*}, where L^{*} is large (≈4000-8000 lattice spacings). In the case of σ=0.75 we have been unable to study large enough systems to determine its value of L^{*}. We sketch a renormalization group scenario to explain how these features could arise. In this scenario finite size effects on the droplet scaling picture in low-dimensional spin glasses produce some aspects of Parisi's replica symmetry breaking theory of the Sherrington-Kirkpatrick model.
Despite decades of effort, our understanding of low-temperature phase of spin glass models with short-range interactions remains incomplete. Replica symmetry breaking (RSB) theory, based on the solution of the Sherrington-Kirkpatrick mean-field model, predicts many pure states; meanwhile, competing theories of short-range systems, such as the droplet picture, predict a single pair of pure states related by time-reversal symmetry, analogously to the ferromagnet. Since RSB certainly holds for the mean-field (infinite-range) model, it is interesting to study short-range models in high dimensions to observe whether RSB also holds here; however, computer simulations of short-range models in high dimensions are difficult because the number of spins to equilibrate grows so rapidly with the linear size of the system.A relatively recent idea which has been fruitful is to instead study one-dimensional models with long-range (power-law) interactions, which are argued to have the same critical behavior as corresponding short-range models in high dimensions, but for which simulations for a range of sizes (crucial for finite-size scaling analysis) are feasible. For these one-dimensional long-range (1DLR) models, we fill in a previously unexplored region of parameter space where the interactions become sufficiently long-range that they must be rescaled with the system size to maintain the thermodynamic limit. We find strong evidence that detailed behavior of the 1DLR models everywhere in this "nonextensive regime" is identical to that of the Sherrington-Kirkpatrick model, lending support to a recent conjecture.In an attempt to distinguish the RSB and droplet pictures, we study recently-proposed observables based on the statistics of individual disorder samples, rather than simply averaging over the disorder as is most frequently done in previous studies. We compare Monte Carlo results for 1DLR models which are proxies for short-range models in 3, 4, and 10 dimensions with previously-obtained data for the 3D and 4D short-range models and the SK model. For one statistic, which is expected to sharply distinguish between the two pictures in the thermodynamic limit, we find that larger system sizes than those currently feasible to simulate are needed to obtain an unambiguous result. We also find that two other recently-proposed statistics, the median of the cumulative overlap distribution and the "typical" overlap distribution, are not particularly helpful in distinguishing between the RSB and droplet pictures.If there are many pure states in the spin-glass phase, we need to carry out some sort of statistical average over them to obtain the thermodynamics. One such prescription for doing this is called the "metastate." Motivated by similarities between the average over pure states specified by the metastate theory and that presumably generated by the nonequilibrium dynamics, we study a 1DLR model which is a proxy for a short-range model in
d=8 dimensions and measure the evolution of dynamical correlations. We find that the spatial decay of the correlations at distances less than the dynamical correlation length
xi(t) agrees quantitatively with the predictions of the metastate theory, evaluated according to the RSB picture. We also compute the dynamic exponent defined by
xi(t)proptot1/z(T) and find that it is compatible with the mean-field value of the critical dynamical exponent for short-range spin glasses.Finally, we present a unified view of finite-size scaling (FSS) in dimensions
d above the upper critical dimension
du, for both free and periodic boundary conditions. For
d>du, a dangerous irrelevant variable is responsible for both the violation of hyperscaling and the violation of "standard" FSS. We find that the modified hyperscaling proposed to allow for this applies only to
veck=vec0 fluctuations, while standard FSS applies to
veckneqvec0 fluctuations. Hence the exponent
eta describing the power-law decay of correlations at criticality is unambiguously
eta=0. With free boundary conditions, the finite-size "shift" is greater than the rounding. Nonetheless, using
T−TL, where
TL is the finite-size pseudocritical temperature, as the scaling variable, the data do collapse onto a scaling form that includes the behavior both at
TL, where the susceptibility
chi diverges like
Ld/2, and the bulk
Tc, where it diverges like
L2. We support these claims with data from large-scale simulations of the five-dimensional Ising model.
We study in Ising spin glasses the finite-size effects near the spin-glass transition in zero field and at the de Almeida-Thouless transition in a field by Monte Carlo methods and by analytical approximations. In zero field, the finite-size scaling function associated with the spin-glass susceptibility of the Sherrington-Kirkpatrick mean-field spin-glass model is of the same form as that of one-dimensional spin-glass models with power-law long-range interactions in the regime where they can be a proxy for the Edwards-Anderson short-range spin-glass model above the upper critical dimension. We also calculate a simple analytical approximation for the spin-glass susceptibility crossover function. The behavior of the spin-glass susceptibility near the de Almeida-Thouless transition line has also been studied, but here we have only been able to obtain analytically its behavior in the asymptotic limit above and below the transition. We have also simulated the one-dimensional system in a field in the non-mean-field regime to illustrate that when the Imry-Ma droplet length scale exceeds the system size one can then be erroneously lead to conclude that there is a de Almeida-Thouless transition even though it is absent.
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.
We study the existence of a line of transitions of an Ising spin glass in a magnetic field---known as the de Almeida-Thouless line---using one-dimensional power-law diluted Ising spin-glass models. We choose the power-law exponent to have values that approximately correspond to three- and four-dimensional nearest-neighbor systems and perform a detailed finite-size scaling analysis of the data for large linear system sizes, using both a new approach proposed recently [Phys. Rev. Lett. 103, 267201 (2009)], as well as traditional approaches. Our results for the model corresponding to a three-dimensional system are consistent with there being no de Almeida-Thouless line, although the new finite-size scaling approach does not rule one out. For the model corresponding to four space dimensions, the new and traditional finite-size scaling analyses give conflicting results, indicating the need for a better understanding of finite-size scaling of spin glasses in a magnetic field.
We study the de Almeida-Thouless (AT) line in the one-dimensional power-law diluted XY spin-glass model, in which the probability that two spins separated by a distance r interact with each other, decays as 1/r^{2σ}. Tuning the exponent σ is equivalent to changing the space dimension of a short-range model. We develop a heat bath algorithm to equilibrate XY spins; using this in conjunction with the standard parallel tempering and overrelaxation sweeps, we carry out large-scale Monte Carlo simulations. For σ=0.6, which is in the mean-field regime above six dimensions-it is similar to being in 10 dimensions-we find clear evidence for an AT line. For σ=0.75 and σ=0.85, which are in the non-mean-field regime and similar to four and three dimensions, respectively, our data is like that found in previous studies of the Ising and Heisenberg spin glasses when reducing the temperature at fixed field. For σ=0.75, there is evidence from finite-size-scaling studies for an AT transition but for σ=0.85, the evidence for a transition is nonexistent. We have also studied these systems at fixed temperature varying the field and discovered that at both σ=0.75 and at σ=0.85 there is evidence of an AT transition! Confusingly, the correlation length and spin-glass susceptibility as a function of the field are both entirely consistent with the predictions of the droplet picture and hence the nonexistence of an AT line. In the usual finite-size critical point scaling studies used to provide evidence for an AT transition, there is seemingly good evidence for an AT line at σ=0.75 for small values of the system size N, which is strengthening as N is increased, but for N>2048 the trend changes and the evidence then weakens as N is further increased. We have also studied with fewer bond realizations the system at σ=0.70, which is the analog of a system with short-range interactions just below six dimensions, and found that it is similar in its behavior to the system at σ=0.75 but with larger finite-size corrections. The evidence from our simulations points to the complete absence of the AT line in dimensions outside the mean-field region and to the correctness of the droplet picture. Previous simulations which suggested there was an AT line can be attributed to the consequences of studying systems which are just too small. The collapse of our data to the droplet scaling form is poor for σ=0.75 and to some extent also for σ=0.85, when the correlation length becomes of the order of the length of the system, due to the existence of excitations which only cost a free energy of O(1), just as envisaged in the TNT picture of the ordered state of spin glasses. However, for the case of σ=0.85 we can provide evidence that for larger system sizes, droplet scaling will prevail even when the correlation length is comparable to the system size.