No Arabic abstract
We discuss generation of series expansions for Ising spin-glasses with a symmetric $pm J$ (i.e. bimodal) distribution on d-dimensional hypercubic lattices using linked-cluster methods. Simplifications for the bimodal distribution allow us to go to higher order than for a general distribution. We discuss two types of problem, one classical and one quantum. The classical problem is that of the Ising spin glass in a longitudinal magnetic field, $h$, for which we obtain high temperature series expansions in variables $tanh(J/T)$ and $tanh(h/T)$. The quantum problem is a $T=0$ study of the Ising spin glass in a transverse magnetic field $h_T$ for which we obtain a perturbation theory in powers of $J/h_T$. These methods require (i) enumeration and counting of textit{all} connected clusters that can be embedded in the lattice up to some order $n$, and (ii) an evaluation of the contribution of each cluster for the quantity being calculated, known as the weight. We discuss a general method that takes the much smaller list (and count) of all no free-end (NFE) clusters on a lattice up to some order $n$, and automatically generates all other clusters and their counts up to the same order. The weights for finite clusters in both cases have a simple graphical interpretation that allows us to proceed efficiently for a general configuration of the $pm J$ bonds, and at the end perform suitable disorder averaging. The order of our computations is limited by the weight calculations for the high-temperature expansions of the classical model, while they are limited by graph counting for the $T=0$ quantum system. Details of the calculational methods are presented.
The strongest evidence for superiority of quantum annealing on spin glass problems has come from comparing simulated quantum annealing using quantum Monte Carlo (QMC) methods to simulated classical annealing [G. Santoro et al., Science 295, 2427(2002)]. Motivated by experiments on programmable quantum annealing devices we revisit the question of when quantum speedup may be expected for Ising spin glass problems. We find that even though a better scaling compared to simulated classical annealing can be achieved for QMC simulations, this advantage is due to time discretization and measurements which are not possible on a physical quantum annealing device. QMC simulations in the physically relevant continuous time limit, on the other hand, do not show superiority. Our results imply that care has to be taken when using QMC simulations to assess quantum speedup potential and are consistent with recent arguments that no quantum speedup should be expected for two-dimensional spin glass problems.
The stability of spin-glass (SG) phase is analyzed in detail for a fermionic Ising SG (FISG) model in the presence of a magnetic transverse field $Gamma$. The fermionic path integral formalism, replica method and static approach have been used to obtain the thermodynamic potential within one step replica symmetry breaking ansatz. The replica symmetry (RS) results show that the SG phase is always unstable against the replicon. Moreover, the two other eigenvalues $lambda_{pm}$ of the Hessian matrix (related to the diagonal elements of the replica matrix) can indicate an additional instability to the SG phase, which enhances when $Gamma$ is increased. Therefore, this result suggests that the study of the replicon can not be enough to guarantee the RS stability in the present quantum FISG model, especially near the quantum critical point. In particular, the FISG model allows changing the occupation number of sites, so one can get a first order transition when the chemical potential exceeds a certain value. In this region, the replicon and the $lambda_{pm}$ indicate instability problems for the SG solution close to all range of first order boundary.
We investigate the performance of the recently proposed stationary Fokker-Planck sampling method considering a combinatorial optimization problem from statistical physics. The algorithmic procedure relies upon the numerical solution of a linear second order differential equation that depends on a diffusion-like parameter D. We apply it to the problem of finding ground states of 2d Ising spin glasses for the +-J-Model. We consider square lattices with side length up to L=24 with two different types of boundary conditions and compare the results to those obtained by exact methods. A particular value of D is found that yields an optimal performance of the algorithm. We compare this optimal value of D to a percolation transition, which occurs when studying the connected clusters of spins flipped by the algorithm. Nevertheless, even for moderate lattice sizes, the algorithm has more and more problems to find the exact ground states. This means that the approach, at least in its standard form, seems to be inferior to other approaches like parallel tempering.
We consider the paramagnetic phase of the random transverse-field Ising spin chain and study the dynamical properties by numerical methods and scaling considerations. We extend our previous work [Phys. Rev. B 57, 11404 (1998)] to new quantities, such as the non-linear susceptibility, higher excitations and the energy-density autocorrelation function. We show that in the Griffiths phase all the above quantities exhibit power-law singularities and the corresponding critical exponents, which vary with the distance from the critical point, can be related to the dynamical exponent z, the latter being the positive root of [(J/h)^{1/z}]_av=1. Particularly, whereas the average spin autocorrelation function in imaginary time decays as [G]_av(t)~t^{-1/z}, the average energy-density autocorrelations decay with another exponent as [G^e]_av(t)~t^{-2-1/z}.
We use a non-equilibrium simulation method to study the spin glass transition in three-dimensional Ising spin glasses. The transition point is repeatedly approached at finite velocity $v$ (temperature change versus time) in Monte Carlo simulations starting at a high temperature. The normally problematic critical slowing-down is not hampering this kind of approach, since the system equilibrates quickly at the initial temperature and the slowing-down is merely reflected in the dynamic scaling of the non-equilibrium order parameter with $v$ and the system size. The equilibrium limit does not have to be reached. For the dynamic exponent we obtain $z = 5.85(9)$ for bimodal couplings distribution and $z=6.00(10)$ for the Gaussian case, thus supporting universal dynamic scaling (in contrast to recent claims of non-universal behavior).