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.

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.

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.

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.

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.

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.

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.

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.

Dual time scales in simulated annealing of a two-dimensional Ising spin glass.

We apply a generalized Kibble-Zurek out-of-equilibrium scaling ansatz to simulated annealing when approaching the spin-glass transition at temperature T=0 of the two-dimensional Ising model with random J=±1 couplings. Analyzing the spin-glass order parameter and the excess energy as functions of the system size and the annealing velocity in Monte Carlo simulations with Metropolis dynamics, we find scaling where the energy relaxes slower than the spin-glass order parameter, i.e., there are two different dynamic exponents. The values of the exponents relating the relaxation time scales to the system length, τ∼L^{z}, are z=8.28±0.03 for the relaxation of the order parameter and z=10.31±0.04 for the energy relaxation. We argue that the behavior with dual time scales arises as a consequence of the entropy-driven ordering mechanism within droplet theory. We point out that the dynamic exponents found here for T→0 simulated annealing are different from the temperature-dependent equilibrium dynamic exponent z_{eq}(T), for which previous studies have found a divergent behavior: z_{eq}(T→0)→∞. Thus, our study shows that, within Metropolis dynamics, it is easier to relax the system to one of its degenerate ground states than to migrate at low temperatures between regions of the configuration space surrounding different ground states. In a more general context of optimization, our study provides an example of robust dense-region solutions for which the excess energy (the conventional cost function) may not be the best measure of success.

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.

Modulated phases in a three-dimensional Maier-Saupe model with competing interactions.

This work is dedicated to the study of the discrete version of the Maier-Saupe model in the presence of competing interactions. The competition between interactions favoring different orientational ordering produces a rich phase diagram including modulated phases. Using a mean-field approach and Monte Carlo simulations, we show that the proposed model exhibits isotropic and nematic phases and also a series of modulated phases that meet at a multicritical point, a Lifshitz point. Though the Monte Carlo and mean-field phase diagrams show some quantitative disagreements, the Monte Carlo simulations corroborate the general behavior found within the mean-field approximation.

Dynamic scaling in the two-dimensional Ising spin glass with normal-distributed couplings.

We carry out simulated annealing and employ a generalized Kibble-Zurek scaling hypothesis to study the two-dimensional Ising spin glass with normal-distributed couplings. The system has an equilibrium glass transition at temperature T=0. From a scaling analysis when T→0 at different annealing velocities v, we find power-law scaling in the system size for the velocity required in order to relax toward the ground state, v∼L^{-(z+1/ν)}, the Kibble-Zurek ansatz where z is the dynamic critical exponent and ν the previously known correlation-length exponent, ν≈3.6. We find z≈13.6 for both the Edwards-Anderson spin-glass order parameter and the excess energy. This is different from a previous study of the system with bimodal couplings [Rubin et al., Phys. Rev. E 95, 052133 (2017)2470-004510.1103/PhysRevE.95.052133] where the dynamics is faster (z is smaller) and the above two quantities relax with different dynamic exponents (with that of the energy being larger). We argue that the different behaviors arise as a consequence of the different low-energy landscapes: for normal-distributed couplings the ground state is unique (up to a spin reflection), while the system with bimodal couplings is massively degenerate. Our results reinforce the conclusion of anomalous entropy-driven relaxation behavior in the bimodal Ising glass. In the case of a continuous coupling distribution, our results presented here also indicate that, although Kibble-Zurek scaling holds, the perturbative behavior normally applying in the slow limit breaks down, likely due to quasidegenerate states, and the scaling function takes a different form.

Quantum criticality with two length scales.

The theory of deconfined quantum critical (DQC) points describes phase transitions at absolute temperature T = 0 outside the standard paradigm, predicting continuous transformations between certain ordered states where conventional theory would require discontinuities. Numerous computer simulations have offered no proof of such transitions, instead finding deviations from expected scaling relations that neither were predicted by the DQC theory nor conform to standard scenarios. Here we show that this enigma can be resolved by introducing a critical scaling form with two divergent length scales. Simulations of a quantum magnet with antiferromagnetic and dimerized ground states confirm the form, proving a continuous transition with deconfined excitations and also explaining anomalous scaling at T > 0. Our findings revise prevailing paradigms for quantum criticality, with potential implications for many strongly correlated materials.

Constrained sampling method for analytic continuation.

A method for analytic continuation of imaginary-time correlation functions (here obtained in quantum Monte Carlo simulations) to real-frequency spectral functions is proposed. Stochastically sampling a spectrum parametrized by a large number of δ functions, treated as a statistical-mechanics problem, it avoids distortions caused by (as demonstrated here) configurational entropy in previous sampling methods. The key development is the suppression of entropy by constraining the spectral weight to within identifiable optimal bounds and imposing a set number of peaks. As a test case, the dynamic structure factor of the S=1/2 Heisenberg chain is computed. Very good agreement is found with Bethe ansatz results in the ground state (including a sharp edge) and with exact diagonalization of small systems at elevated temperatures.

Universal dynamic scaling in three-dimensional Ising spin glasses.

We use a nonequilibrium Monte Carlo simulation method and dynamical scaling to study the phase 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. This approach has the advantage that the equilibrium limit does not have to be strictly reached for a scaling analysis to yield critical exponents. For the dynamic exponent we obtain z=5.85(9) for bimodal couplings distribution and z=6.00(10) for the Gaussian case. Assuming universal dynamic scaling, we combine the two results and obtain z=5.93±0.07 for generic 3D Ising spin glasses.

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.