Critical behavior of the three-dimensional Ising spin glass

Anharmonic theory of superconductivity and its applications to emerging quantum materials

The role of anharmonicity on superconductivity has often been disregarded in the past. Recently, it has been recognized that anharmonic decoherence could play a fundamental role in determining the superconducting properties (electron-phonon coupling, critical temperature, etc) of a large class of materials, including systems close to structural soft-mode instabilities, amorphous solids and metals under extreme high-pressure conditions. Here, we review recent theoretical progress on the role of anharmonic effects, and in particular certain universal properties of anharmonic damping, on superconductivity. Our focus regards the combination of microscopic-agnostic effective theories for bosonic mediators with the well-established BCS theory and Migdal-Eliashberg theory for superconductivity. We discuss in detail the theoretical frameworks, their possible implementation within first-principles methods, and the experimental probes for anharmonic decoherence. Finally, we present several concrete applications to emerging quantum materials, including hydrides, ferroelectrics and systems with charge density wave instabilities.

Proof of Single-Replica Equivalence in Short-Range Spin Glasses

We consider short-range Ising spin glasses in equilibrium at infinite system size, and prove that, for fixed bond realization and a given Gibbs state drawn from a suitable metastate, each translation and locally invariant function (for example, self-overlaps) of a single pure state in the decomposition of the Gibbs state takes the same value for all the pure states in that Gibbs state. We describe several significant applications to spin glasses.

Theory and experiments for disordered elastic manifolds, depinning, avalanches, and sandpiles

Domain walls in magnets, vortex lattices in superconductors, contact lines at depinning, and many other systems can be modeled as an elastic system subject to quenched disorder. The ensuing field theory possesses a well-controlled perturbative expansion around its upper critical dimension. Contrary to standard field theory, the renormalization group (RG) flow involves a function, the disorder correlator Δ(w), and is therefore termed the functional RG. Δ(w) is a physical observable, the auto-correlation function of the center of mass of the elastic manifold. In this review, we give a pedagogical introduction into its phenomenology and techniques. This allows us to treat both equilibrium (statics), and depinning (dynamics). Building on these techniques, avalanche observables are accessible: distributions of size, duration, and velocity, as well as the spatial and temporal shape. Various equivalences between disordered elastic manifolds, and sandpile models exist: an elastic string driven at a point and the Oslo model; disordered elastic manifolds and Manna sandpiles; charge density waves and Abelian sandpiles or loop-erased random walks. Each of the mappings between these systems requires specific techniques, which we develop, including modeling of discrete stochastic systems via coarse-grained stochastic equations of motion, super-symmetry techniques, and cellular automata. Stronger than quadratic nearest-neighbor interactions lead to directed percolation, and non-linear surface growth with additional Kardar-Parisi-Zhang (KPZ) terms. On the other hand, KPZ without disorder can be mapped back to disordered elastic manifolds, either on the directed polymer for its steady state, or a single particle for its decay. Other topics covered are the relation between functional RG and replica symmetry breaking, and random-field magnets. Emphasis is given to numerical and experimental tests of the theory.

Complexity as information in spin-glass Gibbs states and metastates: Upper bounds at nonzero temperature and long-range models

In classical finite-range spin systems, especially those with disorder such as spin glasses, a low-temperature Gibbs state may be a mixture of a number of pure or ordered states; the complexity of the Gibbs state has been defined in the past roughly as the logarithm of this number, assuming the question is meaningful in a finite system. As nontrivial pure-state structure is lost in finite size, in a recent paper [Phys. Rev. E 101, 042114 (2020)2470-004510.1103/PhysRevE.101.042114] Höller and the author introduced a definition of the complexity of an infinite-size Gibbs state as the mutual information between the pure state and the spin configuration in a finite region, and applied this also within a metastate construction. (A metastate is a probability distribution on Gibbs states.) They found an upper bound on the complexity for models of Ising spins in which each spin interacts with only a finite number of others, in terms of the surface area of the region, for all T≥0. In the present paper, the complexity of a metastate is defined likewise in terms of the mutual information between the Gibbs state and the spin configuration. Upper bounds are found for each of these complexities for general finite-range (i.e., short- or long-range, in a sense we define) mixed p-spin interactions of discrete or continuous spins (such as m-vector models), but only for T>0. For short-range models, the bound reduces to the surface area. For long-range interactions, the definition of a Gibbs state has to be modified, and for these models we also prove that the states obtained within the metastate constructions are Gibbs states under the modified definition. All results are valid for a large class of disorder distributions.

Ground-state stability and the nature of the spin glass phase

We propose an approach toward understanding the spin glass phase at zero and low temperature by studying the stability of a spin glass ground state against perturbations of a single coupling. After reviewing the concepts of flexibility, critical droplet, and related quantities for both finite- and infinite-volume ground states, we study some of their properties and review three models in which these quantities are partially or fully understood. We also review a recent result showing the connection between our approach and that of disorder chaos. We then view four proposed scenarios for the low-temperature spin glass phase-replica symmetry breaking, scaling-droplet, TNT, and chaotic pairs-through the lens of the predictions of each scenario for the lowest-energy large-lengthscale excitations above the ground state. Using a new concept called σ-criticality, which quantifies the sensitivity of ground states to single-bond coupling variations, we show that each of these four pictures can be identified with different critical droplet geometries and energies. We also investigate necessary and sufficient conditions for the existence of multiple incongruent ground states.

Nontrivial maturation metastate-average state in a one-dimensional long-range Ising spin glass: Above and below the upper critical range

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.

One-step replica-symmetry-breaking phase below the de Almeida-Thouless line in low-dimensional spin glasses

The de Almeida-Thouless (AT) line is the phase boundary in the temperature-magnetic field plane of an Ising spin glass at which a continuous (i.e., second-order) transition from a paramagnet to a replica-symmetry-breaking (RSB) phase occurs, according to mean-field theory. Here, using field-theoretic perturbative renormalization group methods on the Bray-Roberts reduced Landau-Ginzburg-type theory for a short-range Ising spin glass in space of dimension d, we show that at nonzero magnetic field the nature of the corresponding transition is modified as follows: (a) For d-6 small and positive, with increasing field on the AT line, first, the ordered phase just below the transition becomes the so-called one-step RSB, instead of the full RSB that occurs in mean-field theory; the transition on the AT line remains continuous with a diverging correlation length. Then at a higher field, a tricritical point separates the latter transition from a quasi-first-order one, that is one at which the correlation length does not diverge, and there is a jump in part of the order parameter, but no latent heat. The location of the tricritical point tends to zero as d→6^{+}. (b) For d≤6, we argue that the quasi-first-order transition could persist down to arbitrarily small nonzero fields, with a transition to full RSB still expected at lower temperature. Whenever the quasi-first-order transition occurs, it is at a higher temperature than the AT transition would be for the same field, preempting it as the temperature is lowered. These results may explain the reported absence of a diverging correlation length in the presence of a magnetic field in low-dimensional spin glasses in some simulations and in high-temperature series expansions. We also draw attention to the similarity of the "dynamically frozen" state, which occurs at temperatures just above the quasi-first-order transition, and the "metastate-average state" of the one-step RSB phase, and discuss the issue of the number of pure states in either.

Evidence for Supersymmetry in the Random-Field Ising Model at D=5

We provide a nontrivial test of supersymmetry in the random-field Ising model at five spatial dimensions, by means of extensive zero-temperature numerical simulations. Indeed, supersymmetry relates correlation functions in a D-dimensional disordered system with some other correlation functions in a D-2 clean system. We first show how to check these relationships in a finite-size scaling calculation and then perform a high-accuracy test. While the supersymmetric predictions are satisfied even to our high accuracy at D=5, they fail to describe our results at D=4.

Morphology of renormalization-group flow for the de Almeida-Thouless-Gardner universality class

A replica-symmetry-breaking phase transition is predicted in a host of disordered media. The criticality of the transition has, however, long been questioned below its upper critical dimension, six, due to the absence of a critical fixed point in the renormalization-group flows at one-loop order. A recent two-loop analysis revealed a possible strong-coupling fixed point, but given the uncontrolled nature of perturbative analysis in the strong-coupling regime, debate persists. Here we examine the nature of the transition as a function of spatial dimension and show that the strong-coupling fixed point can go through a Hopf bifurcation, resulting in a critical limit cycle and a concomitant discrete scale invariance. We further investigate a different renormalization scheme and argue that the basin of attraction of the strong-coupling fixed point (or limit cycle) may stay finite for all dimensions.

Binomial Spin Glass

To establish a unified framework for studying both discrete and continuous coupling distributions, we introduce the binomial spin glass, a class of models where the couplings are sums of m identically distributed Bernoulli random variables. In the continuum limit m→∞, the class reduces to one with Gaussian couplings, while m=1 corresponds to the ±J spin glass. We demonstrate that for short-range Ising models on d-dimensional hypercubic lattices the ground-state entropy density for N spins is bounded from above by (sqrt[d/2m]+1/N)ln2, and further show that the actual entropies follow the scaling behavior implied by this bound. We thus uncover a fundamental noncommutativity of the thermodynamic and continuous coupling limits that leads to the presence or absence of degeneracies depending on the precise way the limits are taken. Exact calculations of defect energies reveal a crossover length scale L^{*}(m)∼L^{κ} below which the binomial spin glass is indistinguishable from the Gaussian system. Since κ=-1/(2θ), where θ is the spin-stiffness exponent, discrete couplings become irrelevant at large scales for systems with a finite-temperature spin-glass phase.

Analysis and optimization of population annealing

Population annealing is an easily parallelizable sequential Monte Carlo algorithm that is well suited for simulating the equilibrium properties of systems with rough free-energy landscapes. In this work we seek to understand and improve the performance of population annealing. We derive several useful relations between quantities that describe the performance of population annealing and use these relations to suggest methods to optimize the algorithm. These optimization methods were tested by performing large-scale simulations of the three-dimensional (3D) Edwards-Anderson (Ising) spin glass and measuring several observables. The optimization methods were found to substantially decrease the amount of computational work necessary as compared to previously used, unoptimized versions of population annealing. We also obtain more accurate values of several important observables for the 3D Edwards-Anderson model.

Triviality of the ground-state metastate in long-range Ising spin glasses in one dimension

We consider the one-dimensional model of a spin glass with independent Gaussian-distributed random interactions, which have mean zero and variance 1/|i-j|^{2σ}, between the spins at sites i and j for all i≠j. It is known that, for σ>1, there is no phase transition at any nonzero temperature in this model. We prove rigorously that, for σ>3/2, any translation-covariant Newman-Stein metastate for the ground states (i.e., the frequencies with which distinct ground states are observed in finite-size samples in the limit of infinite size, for given disorder) is trivial and unique. In other words, for given disorder and asymptotically at large sizes, the same ground state, or its global spin flip, is obtained (almost) always. The proof consists of two parts: One is a theorem (based on one by Newman and Stein for short-range two-dimensional models), valid for all σ>1, that establishes triviality under a convergence hypothesis on something similar to the energies of domain walls and the other (based on older results for the one-dimensional model) establishes that the hypothesis is true for σ>3/2. In addition, we derive heuristic scaling arguments and rigorous exponent inequalities which tend to support the validity of the hypothesis under broader conditions. The constructions of various metastates are extended to all values σ>1/2. Triviality of the metastate in bond-diluted power-law models for σ>1 is proved directly.

Numerical Construction of the Aizenman-Wehr Metastate

Chaotic size dependence makes it extremely difficult to take the thermodynamic limit in disordered systems. Instead, the metastate, which is a distribution over thermodynamic states, might have a smooth limit. So far, studies of the metastate have been mostly mathematical. We present a numerical construction of the metastate for the d=3 Ising spin glass. We work in equilibrium, below the critical temperature. Leveraging recent rigorous results, our numerical analysis gives evidence for a dispersed metastate, supported on many thermodynamic states.

Restoration of dimensional reduction in the random-field Ising model at five dimensions

The random-field Ising model is one of the few disordered systems where the perturbative renormalization group can be carried out to all orders of perturbation theory. This analysis predicts dimensional reduction, i.e., that the critical properties of the random-field Ising model in D dimensions are identical to those of the pure Ising ferromagnet in D-2 dimensions. It is well known that dimensional reduction is not true in three dimensions, thus invalidating the perturbative renormalization group prediction. Here, we report high-precision numerical simulations of the 5D random-field Ising model at zero temperature. We illustrate universality by comparing different probability distributions for the random fields. We compute all the relevant critical exponents (including the critical slowing down exponent for the ground-state finding algorithm), as well as several other renormalization-group invariants. The estimated values of the critical exponents of the 5D random-field Ising model are statistically compatible to those of the pure 3D Ising ferromagnet. These results support the restoration of dimensional reduction at D=5. We thus conclude that the failure of the perturbative renormalization group is a low-dimensional phenomenon. We close our contribution by comparing universal quantities for the random-field problem at dimensions 3≤D<6 to their values in the pure Ising model at D-2 dimensions, and we provide a clear verification of the Rushbrooke equality at all studied dimensions.

Efficient numerical methods for the random-field Ising model: Finite-size scaling, reweighting extrapolation, and computation of response functions

It was recently shown [Phys. Rev. Lett. 110, 227201 (2013)PRLTAO0031-900710.1103/PhysRevLett.110.227201] that the critical behavior of the random-field Ising model in three dimensions is ruled by a single universality class. This conclusion was reached only after a proper taming of the large scaling corrections of the model by applying a combined approach of various techniques, coming from the zero- and positive-temperature toolboxes of statistical physics. In the present contribution we provide a detailed description of this combined scheme, explaining in detail the zero-temperature numerical scheme and developing the generalized fluctuation-dissipation formula that allowed us to compute connected and disconnected correlation functions of the model. We discuss the error evolution of our method and we illustrate the infinite limit-size extrapolation of several observables within phenomenological renormalization. We present an extension of the quotients method that allows us to obtain estimates of the critical exponent α of the specific heat of the model via the scaling of the bond energy and we discuss the self-averaging properties of the system and the algorithmic aspects of the maximum-flow algorithm used.

Thermodynamic Identities and Symmetry Breaking in Short-Range Spin Glasses

We present a technique to generate relations connecting pure state weights, overlaps, and correlation functions in short-range spin glasses. These are obtained directly from the unperturbed Hamiltonian and hold for general coupling distributions. All are satisfied in phases with simple thermodynamic structure, such as the droplet-scaling and chaotic pairs pictures. If instead nontrivial mixed-state pictures hold, the relations suggest that replica symmetry is broken as described by a Derrida-Ruelle cascade, with pure state weights distributed as a Poisson-Dirichlet process.

Spin glasses: redux: an updated experimental/materials survey

This article reviews the 40+ year old spin-glass field and one of its earliest model interpretations as a spin density wave. Our description is from an experimental phenomenological point of view with emphasis on new spin glass materials and their relation to topical problems and strongly correlated materials in condensed matter physics. We first simply define a spin glass (SG), give its basic ingredients and explain how the spin glasses enter into the statistical mechanics of classical phase transitions. We then consider the four basic experimental properties to solidly characterize canonical spin glass behavior and introduce the early theories and models. Here the spin density wave (SDW) concept is used to explain the difference between a short-range SDW, i.e. a SG and, in contrast, a long-range SDW, i.e. a conventional magnetic phase transition. We continue with the present state of SG, its massive computer simulations and recent proposals of chiral glasses and quantum SG. We then collect and mention the various SG 'spin-off's'. A major section uncovers the fashionable unconventional materials that display SG-like freezing and glassy ground states, such as (high temperature) superconductors, heavy fermions, intermetallics and Heuslers, pyrochlor and spinels, oxides and chalogenides and exotics, e.g. quasicrystals. Some conclusions and future directions complete the review.

Short-range Ising spin glasses: the metastate interpretation of replica symmetry breaking

Parisi's formal replica-symmetry-breaking (RSB) scheme for mean-field spin glasses has long been interpreted in terms of many pure states organized ultrametrically. However, the early version of this interpretation, as applied to the short-range Edwards-Anderson model, runs into problems because as shown by Newman and Stein (NS) it does not allow for chaotic size dependence, and predicts non-self-averaging that cannot occur. NS proposed the concept of the metastate (a probability distribution over infinite-size Gibbs states in a given sample that captures the effects of chaotic size dependence) and a nonstandard interpretation of the RSB results in which the metastate is nontrivial and is responsible for what was called non-self-averaging. In this picture, each state drawn from the metastate has the ultrametric properties of the old theory, but when the state is averaged using the metastate, the resulting mixed state has little structure. This picture was constructed so as to agree both with the earlier RSB results and with rigorous results. Here we use the effective field theory of RSB, in conjunction with the rigorous definitions of pure states and the metastate in infinite-size systems, to show that the nonstandard picture follows directly from the RSB mean-field theory. In addition, the metastate-averaged state possesses power-law correlations throughout the low-temperature phase; the corresponding exponent ζ takes the value 4 according to the field theory in high dimensions d, and describes the effective fractal dimension of clusters of spins. Further, the logarithm of the number of pure states in the decomposition of the metastate-averaged state that can be distinguished if only correlations in a window of size W can be observed is of order W(d-ζ). These results extend the nonstandard picture quantitatively; we show that arguments against this scenario are inconclusive. More generally, in terms of Parisi's function q(x), if q(0)≠∫(0)(1)dxq(x), then the metastate is nontrivial. In an Appendix, we also prove rigorously that the metastate-averaged state of the Sherrington-Kirkpatrick model is a uniform distribution on all spin configurations at all temperatures.

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.

Minimum spanning trees and random resistor networks in d dimensions

We consider minimum-cost spanning trees, both in lattice and Euclidean models, in d dimensions. For the cost of the optimum tree in a box of size L , we show that there is a correction of order L(theta) , where theta< or =0 is a universal d -dependent exponent. There is a similar form for the change in optimum cost under a change in boundary condition. At nonzero temperature T , there is a crossover length xi approximately T(-nu) , such that on length scales larger than xi, the behavior becomes that of uniform spanning trees. There is a scaling relation theta=-1/nu, and we provide several arguments that show that nu and -1/theta both equal nu(perc) , the correlation length exponent for ordinary percolation in the same dimension d , in all dimensions d> or =1 . The arguments all rely on the close relation of Kruskal's greedy algorithm for the minimum spanning tree, percolation, and (for some arguments) random resistor networks. The scaling of the entropy and free energy at small nonzero T , and hence of the number of near-optimal solutions, is also discussed. We suggest that the Steiner tree problem is in the same universality class as the minimum spanning tree in all dimensions, as is the traveling salesman problem in two dimensions. Hence all will have the same value of theta=-3/4 in two dimensions.

Why temperature chaos in spin glasses is hard to observe

The overlap length of a three-dimensional Ising spin glass on a cubic lattice with Gaussian interactions has been estimated numerically by transfer matrix methods and within a Migdal-Kadanoff renormalization group scheme. We find that the overlap length is large, explaining why it has been difficult to observe spin glass chaos in numerical simulations and experiment.

Hierarchical approach for computing spin glass ground states

We describe a numerical algorithm for computing spin glass ground states with a high level of reliability. The proposed method uses a population based search and applies optimization on multiple scales. Benchmarks are given leading to estimates of the performance on large lattices.

Ground-state degeneracies of Ising spin glasses on diamond hierarchical lattices

The total number of ground states for short-range Ising spin glasses, defined on diamond hierarchical lattices of fractal dimensions d=2, 3, 4, 5, and 2.58, is estimated by means of analytic calculations (three last hierarchy levels of the d=2 lattice) and numerical simulations (lower hierarchies for d=2 and all remaining cases). It is shown that in the case of continuous probability distributions for the couplings, the number of ground states is finite in the thermodynamic limit. However, for a bimodal probability distribution (+/-J with probabilities p and 1-p, respectively), the average number of ground states is maximum for a wide range of values of p around p=1 / 2 and depends on the total number of sites at hierarchy level n, Nn. In this case, for all lattices investigated, it is shown that the ground-state degeneracy behaves like exp[h(d)Nn], in the limit Nn large, where h(d) is a positive number which depends on the lattice fractal dimension. The probability of finding frustrated cells at a given hierarchy level n, Fn(p), is calculated analytically (three last hierarchy levels for d=2 and the last hierarchy of the d=3 lattice, with 0<or=p<or=1), as well as numerically (all other cases, with p=1 / 2). Except for d=2, in which case Fn(1 / 2) increases by decreasing the hierarchy level, all other dimensions investigated present an exponential decrease in Fn(1 / 2) for decreasing values of n. For d=2 our results refer to the paramagnetic phase, whereas for all other dimensions considered [which are greater than the lower critical dimension dl (dl approximately 2.5)], our results refer to the spin-glass phase at zero temperature; in the latter cases h(d) increases with the fractal dimension. For n>>1, only the last hierarchies contribute significantly to the ground-state degeneracy; such a dominant behavior becomes stronger for high fractal dimensions. The exponential increase of the number of ground states with the total number of sites is in agreement with the mean-field picture of spin glasses.