Understanding the low-temperature pure state structure of spin glasses remains an open problem in the field of statistical mechanics of disordered systems. Here we study Monte Carlo dynamics, performing simulations of the growth of correlations following a quench from infinite temperature to a temperature well below the spin-glass transition temperature T_{c} for a one-dimensional Ising spin-glass model with diluted long-range interactions. In this model, the probability P_{ij} that an edge {i,j} has nonvanishing interaction falls as a power law with chord distance, P_{ij}∝1/R_{ij}^{2σ}, and we study a range of values of σ with 1/2<σ<1. We consider a correlation function C_{4}(r,t). A dynamic correlation length that shows power-law growth with time ξ(t)∝t^{1/z} can be identified in the data and, for large time t, C_{4}(r,t) decays as a power law r^{-α_{d}} with distance r when r≪ξ(t). The calculation can be interpreted in terms of the maturation metastate averaged Gibbs state, or MMAS, and the decay exponent α_{d} differentiates between a trivial MMAS (α_{d}=0), as expected in the droplet picture of spin glasses, and a nontrivial MMAS (α_{d}≠0), as in the replica-symmetry-breaking (RSB) or chaotic pairs pictures. We find nonzero α_{d} even in the regime σ>2/3 which corresponds to short-range systems below six dimensions. For σ<2/3, the decay exponent α_{d} follows the RSB prediction for the decay exponent α_{s}=3-4σ of the static metastate, consistent with a conjectured statics-dynamics relation, while it approaches α_{d}=1-σ in the regime 2/3<σ<1; however, it deviates from both lines in the vicinity of σ=2/3.
Abstract We present results of numerical simulations on a one-dimensional Ising spin glass with long-range interactions. Parameters of the model are chosen such that it is a proxy for a short-range spin glass above the upper critical dimension (i.e. in the mean-field regime). The system is quenched to a temperature well below the transition temperature and the growth of correlations is observed. The spatial decay of the correlations at distances less than the dynamic correlation length agrees quantitatively with the predictions of a static theory, the ‘metastate’, evaluated according to the replica symmetry breaking (RSB) theory. We also compute the dynamic exponent z ( T ) defined by and find that it is compatible with the mean-field value of the critical dynamical exponent for short range spin glasses.
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.
For the one-dimensional long-ranged Ising spin-glass with random couplings decaying with the distance r as J ( r ) ∼ r − σ and distributed with the Lévy symmetric stable distribution of index 1 < μ ≤ 2 (including the usual Gaussian case μ = 2), we consider the region σ > 1/ μ where the energy is extensive. We study two real space renormalization procedures at zero temperature, namely a simple box decimation that leads to explicit calculations, and a strong disorder decimation that can be studied numerically on large sizes. The droplet exponent governing the scaling of the renormalized couplings J L ∝ L θ μ ( σ ) is found to be whenever the long-ranged couplings are relevant θ μ ( σ ) ≥ −1. For the statistics of the ground state energy over disordered samples, we obtain that the droplet exponent θ μ ( σ ) governs the leading correction to extensivity of the averaged value . The characteristic scale of the fluctuations around this average is of order , and the rescaled variable is Gaussian distributed for μ = 2, or displays the negative power-law tail in 1/(− u ) 1+ μ for u → − ∞ in the Lévy case 1 < μ < 2. Finally we apply the zero-temperature renormalization procedure to the related Dyson hierarchical spin-glass model where the same droplet exponent appears.
Abstract We have studied numerically the states reached in a quench from various temperatures in the one-dimensional fully-connected Kotliar, Anderson and Stein Ising spin glass model. This is a model where there are long-range interactions between the spins which falls off as a power σ of their separation. We have made a detailed study in particular of the energies of the states reached in a quench from infinite temperature and their overlaps, including the spin glass susceptibility. In the regime where , where the model is similar to the Sherrington–Kirkpatrick model, we find that the spin glass susceptibility diverges logarithmically with increasing N , the number of spins in the system, whereas for it remains finite. We attribute the behavior for to self-organized critical behavior , where the system after the quench is close to the transition between states which have trivial overlaps and those with the non-trivial overlaps associated with replica symmetry breaking. We have also found by studying the distribution of local fields that the states reached in the quench have marginal stability but only when .
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.
Extensive Monte Carlo simulations are employed in order to study the dynamic critical behaviour of the one-dimensional Ising magnet, with algebraically decaying long-range interactions of the form 1/rd + σ, with σ = 0.75. The critical temperature, as well as the critical exponents, are evaluated from the power-law behaviour of suitable physical observables when the system is quenched from uncorrelated states, corresponding to infinite temperature, to the critical point. These results are compared with those obtained from the dynamic evolution of the system when it is annealed at the critical point from the ordered state. Also, the critical temperature in the infinite interaction limit is obtained by means of a finite-range scaling analysis of data measured with different truncated interaction ranges. All the estimated static critical exponents (γ/ν, β/ν, and 1/ν) are in good agreement with renormalization group (RG) results and previously reported numerical data obtained under equilibrium conditions. On the other hand, the dynamic exponent of the initial increase of the magnetization (θ) was close to RG predictions. However, the dynamic exponent z of the time correlation length is slightly different to the RG results probably due to the fact that it may depend on the specific dynamics used or because the two-loop expansion used in the RG analysis may be insufficient.