The standard short-range two-dimensional Ising spin glass is numerically well accessible, in particular, because there are polynomial-time ground-state algorithms. On the other hand, in contrast to higher dimensional spin glasses, it does not exhibit a rich behavior, i.e., no ordered phase at finite temperature. Here, we investigate whether long-range correlated bonds change this behavior. This would still keep the model numerically well accessible while exhibiting a more interesting behavior. The bonds are drawn from a Gaussian distribution with a two-point correlation for bonds at distance r that decays as (1+r^{2})^{-a/2}, a≥0. We study numerically with exact algorithms the ground-state and domain-wall excitations. Our results indicate that the inclusion of bond correlations still does not lead to a spin-glass order at any finite temperature. A further analysis reveals that bond correlations have a strong effect at local length scales, inducing ferro- and antiferromagnetic domains into the system. The length scale of ferro- and antiferromagnetic order diverges exponentially as the correlation exponent approaches a critical value, a→a_{crit}=0. Thus, our results suggest that the system becomes a ferro- or antiferromagnet only in the limit a→0.
The standard short-range two-dimensional Ising spin glass is numerically well accessible, in particular, because there are polynomial-time ground-state algorithms. On the other hand, in contrast to higher dimensional spin glasses, it does not exhibit a rich behavior, i.e., no ordered phase at finite temperature. Here, we investigate whether long-range correlated bonds change this behavior. This would still keep the model numerically well accessible while exhibiting a more interesting behavior. The bonds are drawn from a Gaussian distribution with a two-point correlation for bonds at distance r that decays as (1+r^{2})^{-a/2}, a≥0. We study numerically with exact algorithms the ground-state and domain-wall excitations. Our results indicate that the inclusion of bond correlations still does not lead to a spin-glass order at any finite temperature. A further analysis reveals that bond correlations have a strong effect at local length scales, inducing ferro- and antiferromagnetic domains into the system. The length scale of ferro- and antiferromagnetic order diverges exponentially as the correlation exponent approaches a critical value, a→a_{crit}=0. Thus, our results suggest that the system becomes a ferro- or antiferromagnet only in the limit a→0.
With the help of exact ground states obtained by a polynomial algorithm we compute the domain wall energy Δ at zero temperature for the bond-random and the site-random Ising spin glass model in two dimensions. We find that in both models the stability of the ferromagnetic and the spin glass order ceases to exist at a unique concentration p c for the ferromagnetic bonds. In the vicinity of this critical point, the size and concentration dependence of the first and second moment of the domain wall energy are, for both models, described by a common finite-size scaling form. Moreover, below this concentration the stiffness exponent turns out to be slightly negative, θ S = − 0.056(6), indicating the absence of any intermediate spin glass phase at non-zero temperature.
We study ground-state properties of the two-dimensional random-bond Ising model with couplings having a concentration p∊[0,1] of antiferromagnetic and $(1\ensuremath{-}p)$ of ferromagnetic bonds. We apply an exact matching algorithm which enables us the study of systems with linear dimension L up to 700. We study the behavior of the domain-wall energies and of the magnetization. We find that the paramagnet-ferromagnet transition occurs at ${p}_{c}\ensuremath{\sim}0.103$ compared to the concentration ${p}_{n}\ensuremath{\sim}0.109$ at the Nishimori point, which means that the phase diagram of the model exhibits a re-entrance. Furthermore, we find no indications for an (intermediate) spin-glass ordering at finite temperature.
We perform Monte Carlo simulations of large two-dimensional Gaussian Ising spin glasses down to very low temperatures $\ensuremath{\beta}=1∕T=50$. Equilibration is ensured by using a cluster algorithm including Monte Carlo moves consisting of flipping fundamental excitations. We study the thermodynamic behavior using the Binder cumulant, the spin-glass susceptibility, the distribution of overlaps, the overlap with the ground state, and the specific heat. We confirm that Tc=0. All results are compatible with an algebraic divergence of the correlation length with an exponent $\ensuremath{\nu}$. We find $\ensuremath{-}1∕\ensuremath{\nu}=\ensuremath{-}0.295(40)$, which is compatible with the value for the domain-wall and droplet exponent $\ensuremath{\theta}\ensuremath{\approx}\ensuremath{-}0.29$ found previously in ground-state studies. Hence the thermodynamic behavior of this model seems to be governed by one single exponent.
The d=1 Ising ferromagnet and spin glass with long-range power-law interactions J r^-a is studied for all interaction range exponents a by a renormalization-group transformation that simultaneously projects local ferromagnetism and antiferromagnetism. In the ferromagnetic case, J>0, a finite-temperature ferromagnetic phase occurs for interaction range 0.742, the phase transition temperature discontinuously drops to zero and for a>2 there is no ordered phase above zero temperature, also as predicted by rigorous results. At the other end, on approaching a=0.74 from above, namely increasing the range of the interaction, the phase transition temperature diverges to infinity, meaning that, at all non-infinite temperatures, the system is ferromagnetically ordered. Thus, the equivalent-neighbor interactions regime is entered before (a > 0) the neighbors become equivalent, namely before the interactions become equal for all separations. The critical exponents alpha,beta, gamma,delta,eta,nu are calculated, from a large recursion matrix, varying as function of a. For antiferromagnetic J<0, all triplets of spins at all ranges have competing interactions and this highly frustrated system does not have an ordered phase. In the spin-glass system, where all couplings for all separations are randomly ferromagnetic or antiferromagnetic (with probability p), a finite-temperatures spin-glass phase is obtained in the absence of antiferromagnetic phase. In the spin-glass phase, the signature chaotic behavior under scale change occurs in a richer version than previously: In the long-range interaction of this system, the interactions at every separation become chaotic, yielding a piecewise chaotic interaction function.
The statistics of the ground-state and domain-wall energies for the two-dimensional random-bond Ising model on square lattices with independent, identically distributed bonds having probability p of ${J}_{ij}=\ensuremath{-}1$ and $(1\ensuremath{-}p)$ of Jij=+1 are studied. We are able to consider large samples of up to 3202 spins by using sophisticated matching algorithms. We study $L\ifmmode\times\else\texttimes\fi{}L$ systems, but we also consider $L\ifmmode\times\else\texttimes\fi{}M$ samples, for different aspect ratios R=L∕M. We find that the scaling behavior of the ground-state energy and its sample-to-sample fluctuations inside the spin-glass region $({p}_{c}\ensuremath{\leqslant}p\ensuremath{\leqslant}1\ensuremath{-}{p}_{c})$ are characterized by simple scaling functions. In particular, the fluctuations exhibit a cusplike singularity at pc. Inside the spin-glass region, the average domain-wall energy converges to a finite nonzero value as the sample size becomes infinite, holding R fixed. Here, large finite-size effects are visible, which can be explained for all p by a single exponent $\ensuremath{\omega}\ensuremath{\approx}2∕3$, provided higher-order corrections to scaling are included. Finally, we confirm the validity of aspect-ratio scaling for $R\ensuremath{\rightarrow}0$: the distribution of the domain-wall energies converges to a Gaussian for $R\ensuremath{\rightarrow}0$, although the domain walls of neighboring subsystems of size $L\ifmmode\times\else\texttimes\fi{}L$ are not independent.
We study domain-wall excitations in two-dimensional random-bond Ising spin systems on a square lattice with side length L, subject to two different continuous disorder distributions. In both cases an adjustable parameter allows tuning of the disorder so as to yield a transition from a spin-glass ordered ground state to a ferromagnetic ground state. We formulate an auxiliary graph-theoretical problem in which domain walls are given by undirected shortest paths with possibly negative distances. Due to the details of the mapping, standard shortest-path algorithms cannot be applied. To solve such shortest-path problems we have to apply minimum-weight perfect-matching algorithms. We first locate the critical values of the disorder parameters, where the ferromagnet to spin-glass transition occurs for the two types of the disorder. For certain values of the disorder parameters close to the respective critical point, we investigate the system size dependence of the width of the average domain-wall energy $(\ensuremath{\sim}{L}^{\ensuremath{\theta}})$ and the average domain-wall length $(\ensuremath{\sim}{L}^{{d}_{f}})$. Performing a finite-size scaling analysis for systems with a side length up to L=512, we find that both exponents remain constant in the spin-glass phase, i.e., $\ensuremath{\theta}\ensuremath{\approx}\ensuremath{-}0.28$ and ${d}_{f}\ensuremath{\approx}1.275$. This is consistent with conformal field theory, where it seems to be possible to relate the exponents from the analysis of stochastic Loewner evolutions via ${d}_{f}\ensuremath{-}1=3/[4(3+\ensuremath{\theta})]$. Finally, we characterize the transition in terms of ferromagnetic clusters of spins that form, as one proceeds from spin-glass ordered to ferromagnetic ground states.