Quantum critical behavior of a three-dimensional Ising spin glass in a transverse magnetic field

The quantum transition of the two-dimensional Ising spin glass

Quantum annealers are commercial devices that aim to solve very hard computational problems1, typically those involving spin glasses2,3. Just as in metallurgic annealing, in which a ferrous metal is slowly cooled4, quantum annealers seek good solutions by slowly removing the transverse magnetic field at the lowest possible temperature. Removing the field diminishes the quantum fluctuations but forces the system to traverse the critical point that separates the disordered phase (at large fields) from the spin-glass phase (at small fields). A full understanding of this phase transition is still missing. A debated, crucial question regards the closing of the energy gap separating the ground state from the first excited state. All hopes of achieving an exponential speed-up, compared to classical computers, rest on the assumption that the gap will close algebraically with the number of spins5-9. However, renormalization group calculations predict instead that there is an infinite-randomness fixed point10. Here we solve this debate through extreme-scale numerical simulations, finding that both parties have grasped parts of the truth. Although the closing of the gap at the critical point is indeed super-algebraic, it remains algebraic if one restricts the symmetry of possible excitations. As this symmetry restriction is experimentally achievable (at least nominally), there is still hope for the quantum annealing paradigm11-13.

Quantum critical dynamics in a 5,000-qubit programmable spin glass

Experiments on disordered alloys1-3 suggest that spin glasses can be brought into low-energy states faster by annealing quantum fluctuations than by conventional thermal annealing. Owing to the importance of spin glasses as a paradigmatic computational testbed, reproducing this phenomenon in a programmable system has remained a central challenge in quantum optimization4-13. Here we achieve this goal by realizing quantum-critical spin-glass dynamics on thousands of qubits with a superconducting quantum annealer. We first demonstrate quantitative agreement between quantum annealing and time evolution of the Schrödinger equation in small spin glasses. We then measure dynamics in three-dimensional spin glasses on thousands of qubits, for which classical simulation of many-body quantum dynamics is intractable. We extract critical exponents that clearly distinguish quantum annealing from the slower stochastic dynamics of analogous Monte Carlo algorithms, providing both theoretical and experimental support for large-scale quantum simulation and a scaling advantage in energy optimization.

Development of research network on Quantum Annealing Computation and Information using Google Scholar data

We build and analyse the network of 100 top-cited nodes (research papers and books from Google Scholar; the strength or citation of the nodes range from about 44 000 up to 100) starting in early 1980 until last year. These searched publications (papers and books) are based on Quantum Annealing Computation and Information categorized into four different sets: (A) Quantum/Transverse Field Spin Glass Model, (B) Quantum Annealing, (C) Quantum Adiabatic Computation and (D) Quantum Computation Information in the title or abstract of the searched publications. We fitted the growth in the annual number of publication ([Formula: see text]) in each of these four categories, A-D, to the form [Formula: see text] where [Formula: see text] denotes the time in years. We found the scaling time [Formula: see text] to be of the order of about 10 years for categories A and C, whereas [Formula: see text] is of the order of about 5 years for categories B and D. This article is part of the theme issue 'Quantum annealing and computation: challenges and perspectives'.

Role of quantum fluctuation in inducing ergodicity in the spin glass phase and its effect in quantum annealing

We review the numerical studies on the critical behaviour of the quantum Sherrington-Kirkpatrick (SK) spin glass model, which indicate that a quantum critical behaviour is observed up to a low but non-zero value of temperature. We revisit the numerical investigations on the spin glass order parameter distributions, which identify a low temperature along with high transverse field spin glass phase where the order parameter distribution becomes a delta function in the thermodynamic limit indicating the restoration of replica symmetry and ergodic nature of the system. In the remaining spin glass phase associated with high temperature and low transverse field, the observed distribution is broad akin to the Parisi order parameter distribution. This essentially indicates the non-ergodic behaviour of the system. We further discuss the annealing dynamics studies on the quantum SK model. Such investigations reveal the system size independence of annealing time when the annealing paths go through the ergodic spin glass region. Interestingly, when such dynamics are performed in the non-ergodic spin glass phase the annealing time becomes an increasing function of the system size. Spin autocorrelation shows faster relaxation in the ergodic spin glass region compared with that found in the non-ergodic spin glass region. This article is part of the theme issue 'Quantum annealing and computation: challenges and perspectives'.

Quantum annealing: an overview

In this review, after providing the basic physical concept behind quantum annealing (or adiabatic quantum computation), we present an overview of some recent theoretical as well as experimental developments pointing to the issues which are still debated. With a brief discussion on the fundamental ideas of continuous and discontinuous quantum phase transitions, we discuss the Kibble-Zurek scaling of defect generation following a ramping of a quantum many body system across a quantum critical point. In the process, we discuss associated models, both pure and disordered, and shed light on implementations and some recent applications of the quantum annealing protocols. Furthermore, we discuss the effect of environmental coupling on quantum annealing. Some possible ways to speed up the annealing protocol in closed systems are elaborated upon: we especially focus on the recipes to avoid discontinuous quantum phase transitions occurring in some models where energy gaps vanish exponentially with the system size. This article is part of the theme issue 'Quantum annealing and computation: challenges and perspectives'.

Quantum Criticality in the 2D Quasiperiodic Potts Model

Quantum critical points in quasiperiodic magnets can realize new universality classes, with critical properties distinct from those of clean or disordered systems. Here, we study quantum phase transitions separating ferromagnetic and paramagnetic phases in the quasiperiodic q-state Potts model in 2+1D. Using a controlled real-space renormalization group approach, we find that the critical behavior is largely independent of q, and is controlled by an infinite-quasiperiodicity fixed point. The correlation length exponent is found to be ν=1, saturating a modified version of the Harris-Luck criterion.

Collective Modes at a Disordered Quantum Phase Transition

We study the collective excitations, i.e., the Goldstone (phase) mode and the Higgs (amplitude) mode, near the superfluid-Mott glass quantum phase transition in a two-dimensional system of disordered bosons. Using Monte Carlo simulations as well as an inhomogeneous quantum mean-field theory with Gaussian fluctuations, we show that the Higgs mode is strongly localized for all energies, leading to a noncritical scalar response. In contrast, the lowest-energy Goldstone mode undergoes a striking delocalization transition as the system enters the superfluid phase. We discuss the generality of these findings and experimental consequences, and we point out potential relations to many-body localization.

Griffiths-McCoy singularities in the dilute transverse-field Ising model: A numerical linked cluster expansion study

We use numerical linked cluster expansions (NLCEs) to study the site-diluted transverse-field Ising model on the square lattice at T=0. NLCE with a self-consistent mean field on the boundary of the clusters is used to obtain the ground-state magnetization, susceptibility, and structure factor as a function of transverse field h and exchange constant J. Adding site dilution to the model turns NLCE into a series expansion in the dilution parameter p. Studying the divergence of the structure factor allows us to establish the phase diagram in the h/J and p plane. By studying the magnetization of the system in a longitudinal field, we investigate the Griffiths-McCoy singularities. We find that the magnetization develops nonlinearities in the Griffiths phase with exponents that vary continuously with h. Additionally, the probability distribution of the local susceptibility develops long tails in the Griffiths phase, which is studied in terms of its moments.

Phase transitions in a programmable quantum spin glass simulator

Understanding magnetic phases in quantum mechanical systems is one of the essential goals in condensed matter physics, and the advent of prototype quantum simulation hardware has provided new tools for experimentally probing such systems. We report on the experimental realization of a quantum simulation of interacting Ising spins on three-dimensional cubic lattices up to dimensions 8 × 8 × 8 on a D-Wave processor (D-Wave Systems, Burnaby, Canada). The ability to control and read out the state of individual spins provides direct access to several order parameters, which we used to determine the lattice's magnetic phases as well as critical disorder and one of its universal exponents. By tuning the degree of disorder and effective transverse magnetic field, we observed phase transitions between a paramagnetic, an antiferromagnetic, and a spin-glass phase.

Critical and Griffiths-McCoy singularities in quantum Ising spin glasses on d-dimensional hypercubic lattices: A series expansion study

We study the ±J transverse-field Ising spin-glass model at zero temperature on d-dimensional hypercubic lattices and in the Sherrington-Kirkpatrick (SK) model, by series expansions around the strong-field limit. In the SK model and in high dimensions our calculated critical properties are in excellent agreement with the exact mean-field results, surprisingly even down to dimension d=6, which is below the upper critical dimension of d=8. In contrast, at lower dimensions we find a rich singular behavior consisting of critical and Griffiths-McCoy singularities. The divergence of the equal-time structure factor allows us to locate the critical coupling where the correlation length diverges, implying the onset of a thermodynamic phase transition. We find that the spin-glass susceptibility as well as various power moments of the local susceptibility become singular in the paramagnetic phase before the critical point. Griffiths-McCoy singularities are very strong in two dimensions but decrease rapidly as the dimension increases. We present evidence that high enough powers of the local susceptibility may become singular at the pure-system critical point.

Classical-to-quantum crossover in the critical behavior of the transverse-field Sherrington-Kirkpatrick spin glass model

We study the critical behavior of the Sherrington-Kirkpatrick model in transverse field (at finite temperature) using Monte Carlo simulation and exact diagonalization (at zero temperature). We determine the phase diagram of the model by estimating the Binder cumulant. We also determine the correlation length exponent from the collapse of the scaled data. Our numerical studies here indicate that critical Binder cumulant (indicating the universality class of the transition behavior) and the correlation length exponent cross over from their "classical" to "quantum" values at a finite temperature (unlike the cases of pure systems, where such crossovers occur at zero temperature). We propose a qualitative argument supporting such an observation, employing a simple tunneling picture.

Spin glasses: redux: an updated experimental/materials survey

This article reviews the 40+ year old spin-glass field and one of its earliest model interpretations as a spin density wave. Our description is from an experimental phenomenological point of view with emphasis on new spin glass materials and their relation to topical problems and strongly correlated materials in condensed matter physics. We first simply define a spin glass (SG), give its basic ingredients and explain how the spin glasses enter into the statistical mechanics of classical phase transitions. We then consider the four basic experimental properties to solidly characterize canonical spin glass behavior and introduce the early theories and models. Here the spin density wave (SDW) concept is used to explain the difference between a short-range SDW, i.e. a SG and, in contrast, a long-range SDW, i.e. a conventional magnetic phase transition. We continue with the present state of SG, its massive computer simulations and recent proposals of chiral glasses and quantum SG. We then collect and mention the various SG 'spin-off's'. A major section uncovers the fashionable unconventional materials that display SG-like freezing and glassy ground states, such as (high temperature) superconductors, heavy fermions, intermetallics and Heuslers, pyrochlor and spinels, oxides and chalogenides and exotics, e.g. quasicrystals. Some conclusions and future directions complete the review.

Bose-Einstein condensation in quantum glasses

The role of geometrical frustration in strongly interacting bosonic systems is studied with a combined numerical and analytical approach. We demonstrate the existence of a novel quantum phase featuring both Bose-Einstein condensation and spin-glass behavior. The differences between such a phase and the otherwise insulating "Bose glasses" are elucidated.

Quantum and classical glass transitions in LiHoxY1-xF4

When performed in the proper low-field, low-frequency limits, measurements of the dynamics and the nonlinear susceptibility in the model Ising magnet in a transverse field LiHo(x)Y(1-x)F(4) prove the existence of a spin-glass transition for x=0.167 and 0.198. The classical behavior tracks for the two concentrations, but the behavior in the quantum regime at large transverse fields differs because of the competing effects of quantum entanglement and random fields.

Significance of the hyperfine interactions in the phase diagram of LiHoxY1-xF4

We consider the quantum magnet at LiHo(x)Y(1-x)F(4) at x = 0.167. Experimentally the spin glass to paramagnet transition in this system was studied as a function of the transverse magnetic field and temperature, showing peculiar features: for example, (i) the spin glass order is destroyed much faster by thermal fluctuations than by the transverse field; and (ii) the cusp in the nonlinear susceptibility signaling the glass state decreases in size at lower temperature. Here we show that the hyperfine interactions of the Ho atom must dominate in this system, and that along with the transverse inter-Ho dipolar interactions they dictate the structure of the phase diagram. The experimental observations are shown to be natural consequences of this.

Phase diagram and critical exponents of a dissipative Ising spin chain in a transverse magnetic field

We consider a one-dimensional Ising model in a transverse magnetic field coupled to a dissipative heat bath. The phase diagram and the critical exponents are determined from extensive Monte Carlo simulations. It is shown that the character of the quantum phase transition is radically altered from the corresponding nondissipative model and the double well coupled to a dissipative heat bath with linear friction. Spatial couplings and the dissipative dynamics combine to form a new quantum criticality which is independent of dissipation strength.

Stochastic series expansion method for quantum Ising models with arbitrary interactions

A quantum Monte Carlo algorithm for the transverse Ising model with arbitrary short- or long-range interactions is presented. The algorithm is based on sampling the diagonal matrix elements of the power-series expansion of the density matrix (stochastic series expansion), and avoids the interaction summations necessary in conventional methods. In the case of long-range interactions, the scaling of the computation time with the system size N is therefore reduced from N2 to N ln(N). The method is tested on a one-dimensional ferromagnet in a transverse field, with interactions decaying as 1/r(2).

Cluster Monte Carlo algorithm for the quantum rotor model

We propose a highly efficient "worm"-like cluster Monte Carlo algorithm for the quantum rotor model in the link-current representation. We explicitly prove detailed balance for the algorithm even in the presence of disorder. For the pure quantum rotor model with mu=0, the algorithm yields high- precision estimates for the critical point K(c)=0.333 05(5) and the correlation length exponent nu=0.670(3). For the disordered case, mu=1 / 2+/-1 / 2, we find nu=1.15(10).

Disorder induced cross-over effects at quantum critical points

Critical properties of quantum spin chains with varying degrees of disorder are studied at zero temperature by analytical and extensive density matrix renormalization methods. Generally the phase diagram is found to contain three phases. The weak disorder regime, where the critical behavior is controlled by the fixed points of the pure system, and the strong disorder regime, which is attracted by an infinite randomness fixed point, are separated by an intermediate disorder regime, where dynamical scaling is anisotropic and the static and dynamical exponents are disorder dependent.

Dynamic scaling in diluted systems: deactivation through thermal dilution

Activated scaling is confirmed to hold in transverse field-induced phase transitions of randomly diluted Ising systems. Quantum Monte Carlo calculations have been made not just at the percolation threshold (pc) but also well below and above it. We follow the evolution of the activated scaling at zero temperature in the phase transition from ferromagnetic to quantum Griffiths phase (p>pc) at the phase boundary (p=pc) and for transitions inside the nonferromagnetic quantum Griffiths phase (p<pc). A novel deactivation phenomenon inside the nonferromagnetic Griffiths-McCoy phase (p<pc) is observed using a thermal (in contrast to random) dilution of the system.

From second to first order transitions in a disordered quantum magnet

We study the spin-glass transition in a disordered quantum model. In a region of the phase diagram quantum effects are small and the phase transition is second order, as in the classical case. In another region, quantum fluctuations drive the transition to first order. Across the first order line the susceptibility is discontinuous and shows hysteresis. Our findings qualitatively reproduce observations on LiHoxY1-xF4. We also discuss a marginally stable spin-glass state and derive some results previously obtained from the real-time dynamics of the model coupled to a bath.