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.

Extreme Suppression of Antiferromagnetic Order and Critical Scaling in a Two-Dimensional Random Quantum Magnet.

Sr_{2}CuTeO_{6} is a square-lattice Néel antiferromagnet with superexchange between first-neighbor S=1/2 Cu spins mediated by plaquette centered Te ions. Substituting Te by W, the affected impurity plaquettes have predominantly second-neighbor interactions, thus causing local magnetic frustration. Here we report a study of Sr_{2}CuTe_{1-x}W_{x}O_{6} using neutron diffraction and µSR techniques, showing that the Néel order vanishes already at x=0.025±0.005. We explain this extreme order suppression using a two-dimensional Heisenberg spin model, demonstrating that a W-type impurity induces a deformation of the order parameter that decays with distance as 1/r^{2} at temperature T=0. The associated logarithmic singularity leads to loss of order for any x>0. Order for small x>0 and T>0 is induced by weak interplane couplings. In the nonmagnetic phase of Sr_{2}CuTe_{1-x}W_{x}O_{6}, the µSR relaxation rate exhibits quantum critical scaling with a large dynamic exponent, z≈3, consistent with a random-singlet state.

Multicritical Deconfined Quantum Criticality and Lifshitz Point of a Helical Valence-Bond Phase.

The S=1/2 square-lattice J-Q model hosts a deconfined quantum phase transition between antiferromagnetic and dimerized (valence-bond solid) ground states. We here study two deformations of this model-a term projecting staggered singlets, as well as a modulation of the J terms forming alternating "staircases" of strong and weak couplings. The first deformation preserves all lattice symmetries. Using quantum Monte Carlo simulations, we show that it nevertheless introduces a second relevant field, likely by producing topological defects. The second deformation induces helical valence-bond order. Thus, we identify the deconfined quantum critical point as a multicritical Lifshitz point-the end point of the helical phase and also the end point of a line of first-order transitions. The helical-antiferromagnetic transitions form a line of generic deconfined quantum-critical points. These findings extend the scope of deconfined quantum criticality and resolve a previously inconsistent critical-exponent bound from the conformal-bootstrap method.

Existence of a Spectral Gap in the Affleck-Kennedy-Lieb-Tasaki Model on the Hexagonal Lattice.

The S=1 Affleck-Kennedy-Lieb-Tasaki (AKLT) quantum spin chain was the first rigorous example of an isotropic spin system in the Haldane phase. The conjecture that the S=3/2 AKLT model on the hexagonal lattice is also in a gapped phase has remained open, despite being a fundamental problem of ongoing relevance to condensed-matter physics and quantum information theory. Here we confirm this conjecture by demonstrating the size-independent lower bound Δ>0.006 on the spectral gap of the hexagonal model with periodic boundary conditions in the thermodynamic limit. Our approach consists of two steps combining mathematical physics and high-precision computational physics. We first prove a mathematical finite-size criterion which gives an analytical, size-independent bound on the spectral gap if the gap of a particular cut-out subsystem of 36 spins exceeds a certain threshold value. Then we verify the finite-size criterion numerically by performing state-of-the-art DMRG calculations on the subsystem.

Monte Carlo Renormalization Flows in the Space of Relevant and Irrelevant Operators: Application to Three-Dimensional Clock Models.

We study renormalization group flows in a space of observables computed by Monte Carlo simulations. As an example, we consider three-dimensional clock models, i.e., the XY spin model perturbed by a Z_{q} symmetric anisotropy field. For q=4, 5, 6, a scaling function with two relevant arguments describes all stages of the complex renormalization flow at the critical point and in the ordered phase, including the crossover from the U(1) Nambu-Goldstone fixed point to the ultimate Z_{q} symmetry-breaking fixed point. We expect our method to be useful in the context of quantum-critical points with inherent dangerously irrelevant operators that cannot be tuned away microscopically but whose renormalization flows can be analyzed as we do here for the clock models.

Scaling and Diabatic Effects in Quantum Annealing with a D-Wave Device.

We discuss quantum annealing of the two-dimensional transverse-field Ising model on a D-Wave device, encoded on L×L lattices with L≤32. Analyzing the residual energy and deviation from maximal magnetization in the final classical state, we find an optimal L dependent annealing rate v for which the two quantities are minimized. The results are well described by a phenomenological model with two powers of v and L-dependent prefactors to describe the competing effects of reduced quantum fluctuations (for which we see evidence of the Kibble-Zurek mechanism) and increasing noise impact when v is lowered. The same scaling form also describes results of numerical solutions of a transverse-field Ising model with the spins coupled to noise sources. We explain why the optimal annealing time is much longer than the coherence time of the individual qubits.

Quantum Phases of SrCu_{2}(BO_{3})_{2} from High-Pressure Thermodynamics.

We report heat capacity measurements of SrCu_{2}(BO_{3})_{2} under high pressure along with simulations of relevant quantum spin models and map out the (P,T) phase diagram of the material. We find a first-order quantum phase transition between the low-pressure quantum dimer paramagnet and a phase with signatures of a plaquette-singlet state below T=2 K. At higher pressures, we observe a transition into a previously unknown antiferromagnetic state below 4 K. Our findings can be explained within the two-dimensional Shastry-Sutherland quantum spin model supplemented by weak interlayer couplings. The possibility to tune SrCu_{2}(BO_{3})_{2} between the plaquette-singlet and antiferromagnetic states opens opportunities for experimental tests of quantum field theories and lattice models involving fractionalized excitations, emergent symmetries, and gauge fluctuations.

Critical Level Crossings and Gapless Spin Liquid in the Square-Lattice Spin-1/2 J_{1}-J_{2} Heisenberg Antiferromagnet.

We use the density matrix renormalization group method to calculate several energy eigenvalues of the frustrated S=1/2 square-lattice J_{1}-J_{2} Heisenberg model on 2L×L cylinders with L≤10. We identify excited-level crossings versus the coupling ratio g=J_{2}/J_{1} and study their drifts with the system size L. The lowest singlet-triplet and singlet-quintuplet crossings converge rapidly (with corrections ∝L^{-2}) to different g values, and we argue that these correspond to ground-state transitions between the Néel antiferromagnet and a gapless spin liquid, at g_{c1}≈0.46, and between the spin liquid and a valence-bond solid at g_{c2}≈0.52. Previous studies of order parameters were not able to positively discriminate between an extended spin liquid phase and a critical point. We expect level-crossing analysis to be a generically powerful tool in density matrix renormalization group studies of quantum phase transitions.

Anomalous Quantum-Critical Scaling Corrections in Two-Dimensional Antiferromagnets.

We study the Néel-paramagnetic quantum phase transition in two-dimensional dimerized S=1/2 Heisenberg antiferromagnets using finite-size scaling of quantum Monte Carlo data. We resolve the long-standing issue of the role of cubic interactions arising in the bond-operator representation when the dimer pattern lacks a certain symmetry. We find nonmonotonic (monotonic) size dependence in the staggered (columnar) dimerized model, where cubic interactions are (are not) present. We conclude that there is a new irrelevant field in the staggered model, but, at variance with previous claims, it is not the leading irrelevant field. The new exponent is ω_{2}≈1.25 and the prefactor of the correction L^{-ω_{2}} is large and comes with a different sign from that of the conventional correction with ω_{1}≈0.78. Our study highlights competing scaling corrections at quantum critical points.

Confinement in the Bulk, Deconfinement on the Wall: Infrared Equivalence between Compactified QCD and Quantum Magnets.

In a spontaneously dimerized quantum antiferromagnet, spin-1/2 excitations (spinons) are confined in pairs by strings akin to those confining quarks in non-Abelian gauge theories. The system has multiple degenerate ground states (vacua) and domain walls between regions of different vacua. For two vacua, we demonstrate that spinons on a domain wall are liberated, in a mechanism strikingly similar to domain-wall deconfinement of quarks in variants of quantum chromodynamics. This observation not only establishes a novel phenomenon in quantum magnetism, but also provides a new direct link between particle physics and condensed-matter physics. The analogy opens doors to improving our understanding of particle confinement and deconfinement by computational and experimental studies in quantum magnetism.

Amplitude Mode in Three-Dimensional Dimerized Antiferromagnets.

The amplitude ("Higgs") mode is a ubiquitous collective excitation related to spontaneous breaking of a continuous symmetry. We combine quantum Monte Carlo (QMC) simulations with stochastic analytic continuation to investigate the dynamics of the amplitude mode in a three-dimensional dimerized quantum spin system. We characterize this mode by calculating the spin and dimer spectral functions on both sides of the quantum critical point, finding that both the energies and the intrinsic widths of the excitations satisfy field-theoretical scaling predictions. While the line width of the spin response is close to that observed in neutron scattering experiments on TlCuCl_{3}, the dimer response is significantly broader. Our results demonstrate that highly nontrivial dynamical properties are accessible by modern QMC and analytic continuation methods.

Quantum versus classical annealing: insights from scaling theory and results for spin glasses on 3-regular graphs.

We discuss an Ising spin glass where each S=1/2 spin is coupled antiferromagnetically to three other spins (3-regular graphs). Inducing quantum fluctuations by a time-dependent transverse field, we use out-of-equilibrium quantum Monte Carlo simulations to study dynamic scaling at the quantum glass transition. Comparing the dynamic exponent and other critical exponents with those of the classical (temperature-driven) transition, we conclude that quantum annealing is less efficient than classical simulated annealing in bringing the system into the glass phase. Quantum computing based on the quantum annealing paradigm is therefore inferior to classical simulated annealing for this class of problems. We also comment on previous simulations where a parameter is changed with the simulation time, which is very different from the true Hamiltonian dynamics simulated here.

Anomalous quantum glass of bosons in a random potential in two dimensions.

We present a quantum Monte Carlo study of the "quantum glass" phase of the two-dimensional Bose-Hubbard model with random potentials at filling ρ=1. In the narrow region between the Mott and superfluid phases, the compressibility has the form κ∼exp(-b/T^{α})+c with α<1 and c vanishing or very small. Thus, at T=0 the system is either incompressible (a Mott glass) or nearly incompressible (a Mott-glass-like anomalous Bose glass). At stronger disorder, where a glass reappears from the superfluid, we find a conventional highly compressible Bose glass. On a path connecting these states, away from the superfluid at larger Hubbard repulsion, a change of the disorder strength by only 10% changes the low-temperature compressibility by more than 4 orders of magnitude, lending support to two types of glass states separated by a phase transition or a sharp crossover.

Confinement and deconfinement of spinons in two dimensions.

We use Monte Carlo methods to study spinons in two-dimensional quantum spin systems, characterizing their intrinsic size λ and confinement length Λ. We confirm that spinons are deconfined, Λ→∞ and λ finite, in a resonating valence-bond spin-liquid state. In a valence-bond solid, we find finite λ and Λ, with λ of a single spinon significantly larger than the bound state-the spinon is soft and shrinks as the bound state is formed. Both λ and Λ diverge upon approaching the critical point separating valence-bond solid and Néel ground states. We conclude that the spinon deconfinement is marginal in the lowest-energy state in the spin-1 sector, due to weak attractive spinon interactions. Deconfinement in the vicinity of the critical point should occur at higher energies.

Proximate deconfined quantum critical point in SrCu2(BO3)2.

The deconfined quantum critical point (DQCP) represents a paradigm shift in quantum matter studies, presenting a "beyond Landau" scenario for order-order transitions. Its experimental realization, however, has remained elusive. Using high-pressure 11B nuclear magnetic resonance measurements on the quantum magnet SrCu2(BO3)2, we here demonstrate a magnetic field-induced plaquette singlet to antiferromagnetic transition above 1.8 gigapascals at a notably low temperature, Tc ≃ 0.07 kelvin. First-order signatures of the transition weaken with increasing pressure, and we observe quantum critical scaling at the highest pressure, 2.4 gigapascals. Supported by model calculations, we suggest that these observations can be explained by a proximate DQCP inducing critical quantum fluctuations and emergent O(3) symmetry of the order parameters. Our findings offer a concrete experimental platform for investigation of the DQCP.

Ashkin-teller criticality and pseudo-first-order behavior in a frustrated Ising model on the square lattice.

We study the challenging thermal phase transition to stripe order in the frustrated square-lattice Ising model with couplings J(1) < 0 (nearest-neighbor, ferromagnetic) and J(2) > 0 (second-neighbor, antiferromagnetic) for g = J(2)/|J(1| > 1/2. Using Monte Carlo simulations and known analytical results, we demonstrate Ashkin-Teller criticality for g ≥ g*; i.e., the critical exponents vary continuously between those of the 4-state Potts model at g = g* and the Ising model for g → ∞. Thus, stripe transitions offer a route to realizing a related class of conformal field theories with conformal charge c = 1 and varying exponents. The transition is first order for g < g* = 0.67 ± 0.01, much lower than previously believed, and exhibits pseudo-first-order behavior for |g* ≤ g

Lattice model for the SU(N) Néel to valence-bond solid quantum phase transition at large N.

We generalize the SU(N=2) S=1/2 square-lattice quantum magnet with nearest-neighbor antiferromagnetic coupling (J(1)) and next-nearest-neighbor ferromagnetic coupling (J(2)) to arbitrary N. For all N>4, the ground state has valence-bond-solid order for J(2)=0 and Néel order for J(2)/J(1)≫1, allowing us access to the transition between these types of states for large N. Using quantum Monte Carlo simulations, we show that both order parameters vanish at a single quantum-critical point, whose universal exponents for large enough N (here up to N=12) approach the values obtained in a 1/N expansion of the noncompact CP(N-1) field theory. These results lend strong support to the deconfined quantum-criticality theory of the Néel-valence-bond-solid transition.

Method to characterize spinons as emergent elementary particles.

We develop a technique to directly study spinons (emergent spin S=1/2 particles) in quantum spin models in any number of dimensions. The size of a spinon wave packet and of a bound pair (a triplon) are defined in terms of wave-function overlaps that can be evaluated by quantum Monte Carlo simulations. We show that the same information is contained in the spin-spin correlation function as well. We illustrate the method in one dimension. We confirm that spinons are well-defined particles (have exponentially localized wave packet) in a valence-bond-solid state, are marginally defined (with power-law shaped wave packet) in the standard Heisenberg critical state, and are not well defined in an ordered Néel state (achieved in one dimension using long-range interactions).

Thermodynamics of a gas of deconfined bosonic spinons in two dimensions.

We consider the quantum phase transition between a Néel antiferromagnet and a valence-bond solid (VBS) in a two-dimensional system of S = 1/2 spins. Assuming that the excitations of the critical ground state are linearly dispersing deconfined spinons obeying Bose statistics, we derive expressions for the specific heat C and the magnetic susceptibility χ at low temperature T in terms of a correlation length ξ(T). Comparing with quantum Monte Carlo results for the J-Q model, which is a candidate for a deconfined Néel-VBS transition, we obtain an almost perfect consistency between C, χ, and ξ. The corresponding expressions for magnon (triplet) excitations are not internally consistent, however, lending strong support for spinon excitations in the J-Q model.

Ground states of a frustrated quantum spin chain with long-range interactions.

The ground state of a spin-1/2 Heisenberg chain with both frustration and long-range interactions is studied using Lanczos exact diagonalization. The evolution of the well-known dimerization transition of the system with short-range frustrated interactions (the J1-J2 chain) is investigated in the presence of additional unfrustrated interactions decaying with distance as 1/r(alpha). It is shown that the continuous (infinite-order) dimerization transition develops into a first-order transition between a long-range ordered antiferromagnetic state and a state with coexisting dimerization and critical spin correlations at wave number k=pi/2. The relevance of the model to real systems is discussed.