Edwards-Anderson Ising Spin Glass Research Articles

Weighted averages in population annealing: Analysis and general framework.

Population annealing is a powerful sequential Monte Carlo algorithm designed to study the equilibrium behavior of general systems in statistical physics through massive parallelism. In addition to the remarkable scaling capabilities of the method, it allows for measurements to be enhanced by weighted averaging [J. Machta, Phys. Rev. E 82, 026704 (2010)1539-375510.1103/PhysRevE.82.026704], admitting to reduce both systematic and statistical errors based on independently repeated simulations. We give a self-contained introduction to population annealing with weighted averaging, generalize the method to a wide range of observables such as the specific heat and magnetic susceptibility and rigorously prove that the resulting estimators for finite systems are asymptotically unbiased for essentially arbitrary target distributions. Numerical results based on more than 10^{7} independent population annealing runs of the two-dimensional Ising ferromagnet and the Edwards-Anderson Ising spin glass are presented in depth. In the latter case, we also discuss efficient ways of measuring spin overlaps in population annealing simulations.

Number of thermodynamic states in the three-dimensional Edwards-Anderson spin glass

The question of the number of thermodynamic states present in the low-temperature phase of the three-dimensional Edwards-Anderson Ising spin glass is addressed by studying spin and link overlap distributions using population annealing Monte Carlo simulations. We consider overlaps between systems with the same boundary condition-which are the usual quantities measured-and also overlaps between systems with different boundary conditions, both for the full systems and also within a smaller window within the system. Our results appear to be fully compatible with a single pair of pure states such as in the droplet/scaling picture. However, our results for whether or not domain walls induced by changing boundary conditions are space filling or not are also compatible with scenarios having many thermodynamic states, such as the chaotic pairs picture and the replica symmetry breaking picture. The differing results for spin overlaps in same and different boundary conditions suggest that finite-size effects are very large for the system sizes currently accessible in low-temperature simulations.

Long-time predictability in disordered spin systems following a deep quench.

We study the problem of predictability, or "nature vs nurture," in several disordered Ising spin systems evolving at zero temperature from a random initial state: How much does the final state depend on the information contained in the initial state, and how much depends on the detailed history of the system? Our numerical studies of the "dynamical order parameter" in Edwards-Anderson Ising spin glasses and random ferromagnets indicate that the influence of the initial state decays as dimension increases. Similarly, this same order parameter for the Sherrington-Kirkpatrick infinite-range spin glass indicates that this information decays as the number of spins increases. Based on these results, we conjecture that the influence of the initial state on the final state decays to zero in finite-dimensional random-bond spin systems as dimension goes to infinity, regardless of the presence of frustration. We also study the rate at which spins "freeze out" to a final state as a function of dimensionality and number of spins; here the results indicate that the number of "active" spins at long times increases with dimension (for short-range systems) or number of spins (for infinite-range systems). We provide theoretical arguments to support these conjectures, and also study analytically several mean-field models: the random energy model, the uniform Curie-Weiss ferromagnet, and the disordered Curie-Weiss ferromagnet. We find that for these models, the information contained in the initial state does not decay in the thermodynamic limit-in fact, it fully determines the final state. Unlike in short-range models, the presence of frustration in mean-field models dramatically alters the dynamical behavior with respect to the issue of predictability.

No percolation in low temperature spin glass

We consider the Edwards-Anderson Ising Spin Glass model for temperatures $T\geq 0.$ We define notions of Boltzmann-Gibbs measure for the Edwards-Anderson spin glass at a given temperature, and of unsatisfied (frustrated) edges, and recall the notion of ground states. We prove that for low positive temperatures, in almost every spin configuration the graph formed by the unsatisfied edges is made of finite connected components. Similarly, for zero temperature, we show that in almost every ground state the graph of unsatisfied edges is a forest all of whose components are finite. In other words, for low enough temperatures the unsatisfied edges do not percolate.

Bond chaos in spin glasses revealed through thermal boundary conditions

Spin glasses have competing interactions that lead to a rough energy landscape which is highly susceptible to small perturbations. These chaotic effects strongly affect numerical simulations and, as such, gaining a deeper understanding of chaos in spin glasses is of much importance. The use of thermal boundary conditions is an effective approach to study chaotic phenomena. Here, we generalize population annealing Monte Carlo, combined with thermal boundary conditions, to study bond chaos due to small perturbations in the spin-spin couplings of the three-dimensional Edwards-Anderson Ising spin glass. We show that bond and temperature-induced chaos share the same scaling exponents and that bond chaos is stronger than temperature chaos.

Zero-Temperature Fluctuations in Short-Range Spin Glasses

We consider the energy difference restricted to a finite volume for certain pairs of incongruent ground states (if they exist) in the d-dimensional Edwards-Anderson (EA) Ising spin glass at zero temperature. We prove that the variance of this quantity with respect to the couplings grows at least proportionally to the volume in any dimension greater than or equal to two. An essential aspect of our result is the use of the excitation metastate. As an illustration of potential applications, we use this result to restrict the possible structure of spin glass ground states in two dimensions.

Chaos in spin glasses revealed through thermal boundary conditions

We study the fragility of spin glasses to small temperature perturbations numerically using population annealing Monte Carlo. We apply thermal boundary conditions to a three-dimensional Edwards-Anderson Ising spin glass. In thermal boundary conditions all eight combinations of periodic versus antiperiodic boundary conditions in the three spatial directions are present, each appearing in the ensemble with its respective statistical weight determined by its free energy. We show that temperature chaos is revealed in the statistics of crossings in the free energy for different boundary conditions. By studying the energy difference between boundary conditions at free-energy crossings, we determine the domain-wall fractal dimension. Similarly, by studying the number of crossings, we determine the chaos exponent. Our results also show that computational hardness in spin glasses and the presence of chaos are closely related.

Measuring free energy in spin-lattice models using parallel tempering Monte Carlo.

An efficient and simple approach of measuring the absolute free energy as a function of temperature for spin lattice models using a two-stage parallel tempering Monte Carlo and the free energy perturbation method is discussed and the results are compared with those of population annealing Monte Carlo using the three-dimensional Edwards-Anderson Ising spin glass model as benchmark tests. This approach requires little modification of regular parallel tempering Monte Carlo codes with also little overhead. Numerical results show that parallel tempering, even though using a much less number of temperatures than population annealing, can nevertheless equally efficiently measure the absolute free energy by simulating each temperature for longer times.

Evidence against a mean-field description of short-range spin glasses revealed through thermal boundary conditions

A theoretical description of the low-temperature phase of short-range spin glasses has remained elusive for decades. In particular, it is unclear if theories that assert a single pair of pure states, or theories that are based on infinitely many pure states---such as replica symmetry breaking---best describe realistic short-range systems. To resolve this controversy, the three-dimensional Edwards-Anderson Ising spin glass in thermal boundary conditions is studied numerically using population annealing Monte Carlo. In thermal boundary conditions all eight combinations of periodic vs antiperiodic boundary conditions in the three spatial directions appear in the ensemble with their respective Boltzmann weights, thus minimizing finite-size corrections due to domain walls. From the relative weighting of the eight boundary conditions for each disorder instance a sample stiffness is defined, and its typical value is shown to grow with system size according to a stiffness exponent. An extrapolation to the large-system-size limit is in agreement with a description that supports the droplet picture and other theories that assert a single pair of pure states. The results are, however, incompatible with the mean-field replica symmetry breaking picture, thus highlighting the need to go beyond mean-field descriptions to accurately describe short-range spin-glass systems.

Self-Organized Criticality in Glassy Spin Systems Requires a Diverging Number of Neighbors

We investigate the conditions required for general spin systems with frustration and disorder to display self-organized criticality, a property which so far has been established only for the fully connected infinite-range Sherrington-Kirkpatrick Ising spin-glass model [Phys. Rev. Lett. 83, 1034 (1999)]. Here, we study both avalanche and magnetization jump distributions triggered by an external magnetic field, as well as internal field distributions in the short-range Edwards-Anderson Ising spin glass for various space dimensions between 2 and 8, as well as the fixed-connectivity mean-field Viana-Bray model. Our numerical results, obtained on systems of unprecedented size, demonstrate that self-organized criticality is recovered only in the strict limit of a diverging number of neighbors and is not a generic property of spin-glass models in finite space dimensions.

Correlations between the dynamics of parallel tempering and the free-energy landscape in spin glasses

We present the results of a large-scale numerical study of the equilibrium three-dimensional Edwards-Anderson Ising spin glass with Gaussian disorder. Using parallel tempering (replica exchange) Monte Carlo 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. 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 is examined. Our results show that autocorrelation times are directly correlated with the roughness of the free-energy landscape.

Domain growth and aging scaling in coarsening disordered systems

Using extensive Monte Carlo simulations we study aging properties of two disordered systems quenched below their critical point, namely the two-dimensional random-bond Ising model and the three-dimensional Edwards-Anderson Ising spin glass with a bimodal distribution of the coupling constants. We study the two-times autocorrelation and space-time correlation functions and show that in both systems a simple aging scenario prevails in terms of the scaling variable L(t)/L(s), where L is the time-dependent correlation length, whereas s is the waiting time and t is the observation time. The investigation of the space-time correlation function for the random-bond Ising model allows us to address some issues related to superuniversality.

Uniqueness of Ground States for Short-Range Spin Glasses in the Half-Plane

We consider the Edwards-Anderson Ising spin glass model on the half-plane $Z \times Z^+$ with zero external field and a wide range of choices, including mean zero Gaussian, for the common distribution of the collection J of i.i.d. nearest neighbor couplings. The infinite-volume joint distribution $K(J,\alpha)$ of couplings J and ground state pairs $\alpha$ with periodic (respectively, free) boundary conditions in the horizontal (respectively, vertical) coordinate is shown to exist without need for subsequence limits. Our main result is that for almost every J, the conditional distribution $K(\alpha|J)$ is supported on a single ground state pair.

Exact algorithm for sampling the two-dimensional Ising spin glass

A sampling algorithm is presented that generates spin-glass configurations of the two-dimensional Edwards-Anderson Ising spin glass at finite temperature with probabilities proportional to their Boltzmann weights. Such an algorithm overcomes the slow dynamics of direct simulation and can be used to study long-range correlation functions and coarse-grained dynamics. The algorithm uses a correspondence between spin configurations on a regular lattice and dimer (edge) coverings of a related graph: Wilson's algorithm [D. B. Wilson, Proceedings of the Eighth Symposium on Discrete Algorithms (SIAM, Philadelphia, 1997), p 258] for sampling dimer coverings on a planar lattice is adapted to generate samplings for the dimer problem corresponding to both planar and toroidal spin-glass samples. This algorithm is recursive: it computes probabilities for spins along a "separator" that divides the sample in half. Given the spins on the separator, sample configurations for the two separated halves are generated by further division and assignment. The algorithm is simplified by using Pfaffian elimination rather than Gaussian elimination for sampling dimer configurations. For n spins and given floating point precision, the algorithm has an asymptotic run-time of O(n(3/2)); it is found that the required precision scales as inverse temperature and grows only slowly with system size. Sample applications and benchmarking results are presented for samples of size up to n=128(2), with fixed and periodic boundary conditions.

Make Life Simple: Unleash the Full Power of the Parallel Tempering Algorithm

We introduce a new update scheme to systematically improve the efficiency of parallel tempering simulations. We show that, by adapting the number of sweeps between replica exchanges to the canonical autocorrelation time, the average round-trip time of a replica in temperature space can be significantly decreased. The temperatures are not dynamically adjusted as in previous attempts but chosen to yield a 50% exchange rate of adjacent replicas. We illustrate the new algorithm with results for the Ising model in two and the Edwards-Anderson Ising spin glass in three dimensions.

Comment on “Ultrametricity in the Edwards-Anderson Model”

In a recent interesting Letter Contucci {\it et al.} have investigated several properties of the three-dimensional (3d) Edwards-Anderson (EA) Ising spin glass. They claim to have found strong numerical evidence for the presence of a complex ultrametric structure similar to the one described by the replica symmetry breaking solution of the mean field model. We illustrate by numerical simulations that the relations used by Contucci {\it et al.} as evidence for an ultrametric structure in the 3d EA model are fulfilled to similar accuracy in the two-dimensional EA model, which is well-described by the droplet picture and has no spin glass phase at finite temperature. We conclude that the data presented in the Contucci {\it et al.} Letter is not sufficient to dismiss the possibility that, e.g., the droplet model might describe the behavior of the 3d EA model.

Finite versus zero-temperature hysteretic behavior of spin glasses: Experiment and theory

We present experimental results attempting to fingerprint nonanalyticities in the magnetization curves of spin glasses found by Katzgraber et al. [Phys. Rev. Lett. 89, 257202 (2002)] via zero-temperature Monte Carlo simulations of the Edwards-Anderson Ising spin glass. Our results show that the singularities at zero temperature due to the reversal-field memory effect are washed out by the finite temperatures of the experiments. The data are analyzed via the first order reversal curve (FORC) magnetic fingerprinting method. The experimental results are supported by Monte Carlo simulations of the Edwards-Anderson Ising spin glass at finite temperatures which agree qualitatively very well with the experimental results. This suggests that the hysteretic behavior of real Ising spin-glass materials is well described by the Edwards-Anderson Ising spin glass. Furthermore, reversal-field memory is a purely zero-temperature effect.

Scaling Analysis of Domain-Wall Free Energy in the Edwards-Anderson Ising Spin Glass in a Magnetic Field

The stability of the spin-glass phase against a magnetic field is studied in the three- and four-dimensional Edwards-Anderson Ising spin glasses. Effective couplings J(eff) and effective fields H(eff) associated with length scale L are measured by a numerical domain-wall renormalization-group method. The results obtained by scaling analysis of the data strongly indicate the existence of a crossover length beyond which the spin-glass order is destroyed by field H. The crossover length well obeys a power law of H which diverges as H --> 0 but remains finite for any nonzero H, implying that the spin-glass phase is absent even in an infinitesimal field. These results are well consistent with the droplet theory for short-range spin glasses.

Universality in three-dimensional Ising spin glasses: A Monte Carlo study

We study universality in three-dimensional Ising spin glasses by large-scale Monte Carlo simulations of the Edwards-Anderson Ising spin glass for several choices of bond distributions, with particular emphasis on Gaussian and bimodal interactions. A finite-size scaling analysis suggests that three-dimensional spin glasses obey universality.

Temperature Chaos and Bond Chaos in Edwards-Anderson Ising Spin Glasses: Domain-Wall Free-Energy Measurements

Domain-wall free energy sigmaF, entropy sigmaS, and the correlation function C(temp) of sigmaF are measured independently in the four-dimensional +/-J Edwards-Anderson (EA) Ising spin glass. The stiffness exponent theta, the fractal dimension of domain walls d(s), and the chaos exponent zeta are extracted from the finite-size scaling analysis of sigmaF, sigmaS, and C(temp), respectively, well inside the spin-glass phase. The three exponents are confirmed to satisfy the scaling relation zeta = d(s)/2 - theta derived by the droplet theory within our numerical accuracy. We also study bond chaos induced by random variation of bonds, and find that the bond and temperature perturbations yield the universal chaos effects described by a common scaling function and the chaos exponent. These results strongly support the appropriateness of the droplet theory for the description of chaos effect in the EA Ising spin glasses.