Monte Carlo data of the two-dimensional Ising spin glass with bimodal interactions are presented with the aim of understanding the low-temperature physics of the model. An analysis of the specific heat, spin-glass susceptibility, finite-size correlation length, and the Binder ratio is performed to try to verify a recent proposal in which for large system sizes and finite but low temperatures the effective critical exponents are identical to the critical exponents of the two-dimensional Ising spin glass with Gaussian interactions. Our results show that with present system sizes the recently proposed scenario, in which the two-dimensional Ising spin glass with bimodally distributed interactions is in the same universality class as the model with Gaussian-distributed disorder at low but finite temperatures, cannot be reliably proven.
Scaling arguments and precise simulations are used to study the square lattice ±J Ising spin glass, a prototypical model for glassy systems. Droplet theory explains, and our numerical results show, entropically stabilized long-range spin-glass order at zero temperature, which resembles the energetic stabilization of long-range order in higher-dimensional models at finite temperature. At low temperature, a temperature-dependent crossover length scale is used to predict the power-law dependence on temperature of the heat capacity and clarify the importance of disorder distributions.
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.