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.
An efficient Monte Carlo method is extended to evaluate directly domain-wall free energy for randomly frustrated spin systems. Using the method, critical phenomena of spin-glass phase transition are investigated in the 4d+/-J Ising model under the replica boundary condition. Our values of the critical temperature and exponent, obtained by finite-size scaling, are in good agreement with those of the standard Monte Carlo and the series expansion studies. In addition, two exponents, the stiffness exponent and the fractal dimension of the domain wall, which characterize the ordered phase, are obtained. The latter value is larger than d-1, indicating that the domain wall is really rough in the 4d Ising spin-glass phase.
The ground-state properties of quasi-one-dimensional (Q1D) Ising spin glass are investigated using an exact numerical approach and analytical arguments. A set of coupled recursive equations for the ground-state energy are introduced and solved numerically. For various types of coupling distribution, we obtain accurate results for magnetization, particularly in the presence of a weak external magnetic field. We show that in the weak magnetic field limit, similar to the 1D model, magnetization exhibits a singular power-law behavior with divergent susceptibility. Remarkably, the spectrum of magnetic exponents is markedly different from that of the 1D system even in the case of two coupled chains. The magnetic exponent makes a crossover from being dependent on a distribution function to a constant value independent of distribution. We provide an analytic theory for these observations by extending the Chen-Ma argument to the Q1D case. We derive an analytical formula for the exponent which is in perfect agreement with the numerical results.
A new method to investigate the two-point correlation function decay in a finite size random spin system and at finite temperatures is proposed. It is applicable to systems with cyclic boundary conditions and in any local or uniform magnetic field. The efficiency of this method is presented in the case of an isotropic XY spin glass model with Gaussian bond distribution in external field. Using unitary U(1) gauge transformation and operator inequalities we have obtained explicit upper bounds for generalized spin-glass susceptibility κ(T,R) that describes off-diagonal transverse magnetization in two dimensional random XY spin lattices at finite temperatures. It is found that the slowest possible decay for susceptibility in random XY spin glass is of the Kosterlitz–Thouless-type with power-law decay in 2D, and the exponential type in 1D, at low temperatures. These upper bounds rule out the possibility of the corresponding magnetic ordering in 1D and 2D isotropic XY spin glasses at finite temperatures.
Exact ground states of two-dimensional Ising spin glasses with Gaussian and bimodal $(\ifmmode\pm\else\textpm\fi{}J)$ distributions of the disorder are calculated using a ``matching''algorithm, which allows large system sizes of up to
N=4802 spins to be investigated. We study domain walls induced by two rather different types of boundary condition changes, and, in each case, analyze the system-size dependence of an appropriately defined ``defect energy,'' which we denote by $\ensuremath{\Delta}E.$ For Gaussian disorder, we find a power-law behavior $\ensuremath{\Delta}E\ensuremath{\sim}{L}^{\ensuremath{\theta}},$ with $\ensuremath{\theta}=\ensuremath{-}0.266(2)$ and $\ensuremath{\theta}=\ensuremath{-}0.282(2)$ for the two types of boundary condition changes. These results are in reasonable agreement with each other, allowing for small systematic effects. They also agree well with earlier work on smaller sizes. The negative value indicates that two dimensions is below the lower critical dimension
dc. For the $\ifmmode\pm\else\textpm\fi{}J$ model, we obtain a different result, namely that the domain-wall energy saturates at a nonzero value for $\stackrel{\ensuremath{\rightarrow}}{L}\ensuremath{\infty},$ so $\ensuremath{\theta}=0,$ indicating that the lower critical dimension for the $\ifmmode\pm\else\textpm\fi{}J$ model is exactly
dc=2.We study the magnetization of the ground states of Ising spin glasses as a function of the magnetic field. The grids are two-dimensional with Gaussian bond distribution. Finding a ground state in this case is a difficult combinatorial problem. We present a method that guarantees finding a true ground state. Although it is unlikely that one can have an algorithm with a polynomial bound for the running time, in all cases we tried, our method found a solution in a reasonable amount of time. We used this to obtain estimations of the critical exponents.