One of the key predictions of Parisi's broken replica symmetry theory of spin glasses is the existence of a phase transition in an applied field to a state with broken replica symmetry. This transition takes place at the de Almeida-Thouless (AT) line in the h-T plane. We have studied this line in the power-law diluted Heisenberg spin glass in which the probability that two spins separated by a distance r interact with each other falls as 1/r^{2σ}. In the presence of a random vector field of variance h_{r}^{2} the phase transition is in the universality class of the Ising spin glass in a field. Tuning σ is equivalent to changing the dimension d of the short-range system, with the relation being d=2/(2σ-1) for σ<2/3. We have found by numerical simulations that h_{AT}^{2}∼(2/3-σ) implying that the AT line does not exist below six dimensions and that the Parisi scheme is not appropriate for spin glasses in three dimensions.
One of the key predictions of Parisi's broken replica symmetry theory of spin glasses is the existence of a phase transition in an applied field to a state with broken replica symmetry. This transition takes place at the de Almeida-Thouless (AT) line in the h-T plane. We have studied this line in the power-law diluted Heisenberg spin glass in which the probability that two spins separated by a distance r interact with each other falls as 1/r^{2σ}. In the presence of a random vector field of variance h_{r}^{2} the phase transition is in the universality class of the Ising spin glass in a field. Tuning σ is equivalent to changing the dimension d of the short-range system, with the relation being d=2/(2σ-1) for σ<2/3. We have found by numerical simulations that h_{AT}^{2}∼(2/3-σ) implying that the AT line does not exist below six dimensions and that the Parisi scheme is not appropriate for spin glasses in three dimensions.
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.
We test for the presence or absence of the de Almeida--Thouless (AT) line using one-dimensional power-law diluted Heisenberg spin-glass model, in which the rms strength of the interactions decays with distance, r as $1/{r}^{\ensuremath{\sigma}}$. It is argued that varying the power $\ensuremath{\sigma}$ is analogous to varying the space dimension d in a short-range model. For $\ensuremath{\sigma}=0.6$, which is in the mean-field regime regime, we find clear evidence for an AT line. For $\ensuremath{\sigma}=0.85$, which is in the non-mean-field regime and corresponds to a space dimension of close to 3, we find no AT line, though we cannot rule one out for very small fields. Finally for $\ensuremath{\sigma}=0.75$, which is in the non-mean-field regime but closer to the mean-field boundary, the evidence suggests that there is an AT line, though the possibility that even larger sizes are needed to see the asymptotic behavior can not be ruled out.
The de Almeida-Thouless (AT) line in Ising spin glasses is the phase boundary in the temperature T and magnetic field h plane below which replica symmetry is broken. Using perturbative renormalization group (RG) methods, we show that, when the dimension d of space is just above six, there is a multicritical point (MCP) on the AT line, which separates a low-field regime, in which the critical exponents have mean-field values, from a high-field regime, where the RG flows run away to infinite coupling strength; as d approaches six from above, the MCP approaches the zero-field critical point exponentially in 1/(d-6). Thus, on the AT line, perturbation theory for the critical properties breaks down at a sufficiently large magnetic field even above 6 dimensions, as well as for all nonzero fields when d≤6, as was known previously. We calculate the exponents at the MCP to first order in ϵ=d-6>0. The fate of the MCP as d increases from just above six to infinity is not known.
The nature of spin-glass states in a magnetic field remains a major open problem in statistical physics. The existence of the de Almeida-Thouless (dAT) transition for three-dimensional (3D) spin glasses in a field is still debated. We introduce a computational method to define the spin-glass susceptibility, which is robust against the broad tail in the overlap distribution that undermines conventional analyses. Applying this approach to the Edwards-Anderson spin-glass model in 2D and 3D, and contrasting with the 3D Ising (without disorder) and mean-field spin-glass models, we find a stark difference: the locus of susceptibility maxima bends to the right in the field-temperature plane for the Ising and 2D spin-glass cases, indicating a supercritical crossover line, but bends to the left for the mean field and 3D spin glasses-a signature of the dAT line. Finite-size scaling further suggests that the peak susceptibility diverges with system size in 3D spin glasses under a field, while saturating in 2D. These results provide direct numerical evidence for the dAT transition in 3D, supporting the replica symmetry breaking scenario.
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.