Series expansions for the Ising spin glass in general dimension

de Almeida-Thouless instability in short-range Ising spin glasses

We use high-temperature series expansions to study the ±J Ising spin glass in a magnetic field in d-dimensional hypercubic lattices for d=5-8 and in the infinite-range Sherrington-Kirkpatrick (SK) model. The expansions are obtained in the variable w=tanh^{2}J/T for arbitrary values of u=tanh^{2}h/T complete to order w^{10}. We find that the scaling dimension Δ associated with the ordering-field h^{2} equals 2 in the SK model and for d≥6. However, in agreement with the work of Fisher and Sompolinsky [Phys. Rev. Lett. 54, 1063 (1985)PRLTAO0031-900710.1103/PhysRevLett.54.1063], there is a violation of scaling in a finite field, leading to an anomalous h-T dependence of the de Almeida-Thouless (AT) [J. Phys. A 11, 983 (1978)JPHAC50305-447010.1088/0305-4470/11/5/028] line in high dimensions, whereas scaling is restored as d→6. Within the convergence of our series analysis, we present evidence supporting an AT line in d≥6. In d=5, the exponents γ and Δ are substantially larger than mean-field values, but we do not see clear evidence for the AT line in d=5.

Ising spin glasses in dimension five

Ising spin-glass models with bimodal, Gaussian, uniform, and Laplacian interaction distributions in dimension five are studied through detailed numerical simulations. The data are analyzed in both the finite-size scaling regime and the thermodynamic limit regime. It is shown that the values of critical exponents and of dimensionless observables at criticality are model dependent. Models in a single universality class have identical values for each of these critical parameters, so Ising spin-glass models in dimension five with different interaction distributions each lie in different universality classes. This result confirms conclusions drawn from measurements in dimension four and dimension two.

Bimodal and Gaussian Ising spin glasses in dimension two

An analysis is given of numerical simulation data to size L=128 on the archetype square lattice Ising spin glasses (ISGs) with bimodal (±J) and Gaussian interaction distributions. It is well established that the ordering temperature of both models is zero. The Gaussian model has a nondegenerate ground state and thus a critical exponent η≡0, and a continuous distribution of energy levels. For the bimodal model, above a size-dependent crossover temperature T(*)(L) there is a regime of effectively continuous energy levels; below T(*)(L) there is a distinct regime dominated by the highly degenerate ground state plus an energy gap to the excited states. T(*)(L) tends to zero at very large L, leaving only the effectively continuous regime in the thermodynamic limit. The simulation data on both models are analyzed with the conventional scaling variable t=T and with a scaling variable τ(b)=T(2)/(1+T(2)) suitable for zero-temperature transition ISGs, together with appropriate scaling expressions. The data for the temperature dependence of the reduced susceptibility χ(τ(b),L) and second moment correlation length ξ(τ(b),L) in the thermodynamic limit regime are extrapolated to the τ(b)=0 critical limit. The Gaussian critical exponent estimates from the simulations, η=0 and ν=3.55(5), are in full agreement with the well-established values in the literature. The bimodal critical exponents, estimated from the thermodynamic limit regime analyses using the same extrapolation protocols as for the Gaussian model, are η=0.20(2) and ν=4.8(3), distinctly different from the Gaussian critical exponents.

Non-self-averaging in Ising spin glasses and hyperuniversality

Ising spin glasses with bimodal and Gaussian near-neighbor interaction distributions are studied through numerical simulations. The non-self-averaging (normalized intersample variance) parameter U_{22}(T,L) for the spin glass susceptibility [and for higher moments U_{nn}(T,L)] is reported for dimensions 2,3,4,5, and 7. In each dimension d the non-self-averaging parameters in the paramagnetic regime vary with the sample size L and the correlation length ξ(T,L) as U_{nn}(β,L)=[K_{d}ξ(T,L)/L]^{d} and so follow a renormalization group law due to Aharony and Harris [Phys. Rev. Lett. 77, 3700 (1996)PRLTAO0031-900710.1103/PhysRevLett.77.3700]. Empirically, it is found that the K_{d} values are independent of d to within the statistics. The maximum values [U_{nn}(T,L)]_{max} are almost independent of L in each dimension, and remarkably the estimated thermodynamic limit critical [U_{nn}(T,L)]_{max} peak values are also practically dimension-independent to within the statistics and so are "hyperuniversal." These results show that the form of the spin-spin correlation function distribution at criticality in the large L limit is independent of dimension within the ISG family. Inspection of published non-self-averaging data for three-dimensional Heisenberg and XY spin glasses the light of the Ising spin glass non-self-averaging results show behavior which appears to be compatible with that expected on a chiral-driven ordering interpretation but incompatible with a spin-driven ordering scenario.

Evidence for nonuniversal scaling in dimension-four Ising spin glasses

The critical behavior of the Binder cumulant for Ising spin glasses in dimension four is studied through simulation measurements. Data for the bimodal interaction model are compared with those for the Laplacian interaction model. Special attention is paid to scaling corrections. The limiting infinite size value at criticality for this dimensionless variable is a parameter characteristic of a universality class. This critical limit is estimated to be equal to 0.523(3) in the bimodal model and to 0.473(3) in the Laplacian model.

High-temperature series expansions for the q-state Potts model on a hypercubic lattice and critical properties of percolation

We present results for the high-temperature series expansions of the susceptibility and free energy of the q-state Potts model on a D-dimensional hypercubic lattice ZD for arbitrary values of q. The series are up to order 20 for dimension D<or=3, order 19 for D<or=5, and up to order 17 for arbitrary D. Using the q-->1 limit of these series, we estimate the percolation threshold pc and critical exponent gamma for bond percolation in different dimensions. We also extend the 1/D expansion of the critical coupling for arbitrary values of q up to order D-9.

Stiffness of the Edwards-Anderson model in all dimensions

A comprehensive description in all dimensions is provided for the scaling exponent y of low-energy excitations in the Ising spin glass introduced by Edwards and Anderson. A combination of extensive numerical as well as theoretical results suggest that its lower critical dimension is exactly dl=5/2. Such a result would be an essential feature of any complete theory of low-temperature spin glass order and imposes a constraint that may help to distinguish between theories.