We resolve the nature of the quantum phase transition between antiferromagnetic and valence-bond solid ground states in square-lattice spin-1/2 magnets - the putative deconfined quantum-critical point. Using quantum Monte Carlo simulations of J-Q models, in which Heisenberg exchange J competes with multi-spin interactions Qn formed by n singlet projectors, we demonstrate first-order discontinuities increasing with n. For n=2 and n=3, the first-order signatures are very weak and we extract universal scaling dimensions of a proximate multicritical point. Their values coincide closely with known bounds for an SO(5) symmetric conformal field theory (CFT) containing a relevant SO(5) singlet operator that destabilizes the generic deconfined critical point. We characterize the emergent SO(5) symmetry by its scaling dimensions of the leading irrelevant perturbations. Using the J-Q2-Q6 model, we tune the transition from weak to stronger first-order and find scaling of the coexisting order parameters with a large exponent beta ~ 0.85, in agreement with the relevant scaling dimensions of the multicritical point. The exponent that had previously been identified with the correlation length does not have a corresponding CFT level, and we explain it as pseudo-critical length scale of bubbles of one phase inside the other that are compressed in size by an SO(5)-violating irrelevant repulsive interaction between the two order parameters. The bubble scale reflects a hierarchy of excitation gaps generated by the SO(5) perturbation in the neighborhood of the first-order transition, akin to a Higgs mechanism. We argue that the multicritical point is the tip of a gapless spin liquid phase in frustrated Heisenberg models into which the J-Q models can be deformed, but the spin liquid, and perhaps the multicritical point itself, exist outside the regime amenable to sign-free simulations.
We resolve the nature of the quantum phase transition between antiferromagnetic and valence-bond solid ground states in square-lattice spin-1/2 magnets - the putative deconfined quantum-critical point. Using quantum Monte Carlo simulations of J-Q models, in which Heisenberg exchange J competes with multi-spin interactions Qn formed by n singlet projectors, we demonstrate first-order discontinuities increasing with n. For n=2 and n=3, the first-order signatures are very weak and we extract universal scaling dimensions of a proximate multicritical point. Their values coincide closely with known bounds for an SO(5) symmetric conformal field theory (CFT) containing a relevant SO(5) singlet operator that destabilizes the generic deconfined critical point. We characterize the emergent SO(5) symmetry by its scaling dimensions of the leading irrelevant perturbations. Using the J-Q2-Q6 model, we tune the transition from weak to stronger first-order and find scaling of the coexisting order parameters with a large exponent beta ~ 0.85, in agreement with the relevant scaling dimensions of the multicritical point. The exponent that had previously been identified with the correlation length does not have a corresponding CFT level, and we explain it as pseudo-critical length scale of bubbles of one phase inside the other that are compressed in size by an SO(5)-violating irrelevant repulsive interaction between the two order parameters. The bubble scale reflects a hierarchy of excitation gaps generated by the SO(5) perturbation in the neighborhood of the first-order transition, akin to a Higgs mechanism. We argue that the multicritical point is the tip of a gapless spin liquid phase in frustrated Heisenberg models into which the J-Q models can be deformed, but the spin liquid, and perhaps the multicritical point itself, exist outside the regime amenable to sign-free simulations.
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.
We introduce a quantum spin-1/2 model with many-body correlated Heisenberg-type interactions on the two-dimensional square lattice, designed so that the system can host a fourfold degenerate plaquette valencebond solid (PVBS) ground state that spontaneously breaks Z 4 symmetry.The system is sign-problem free and amenable to large-scale quantum Monte Carlo simulations, thus allowing us to carry out a detailed study of the quantum phase transition between the standard Néel antiferromagnetic (AFM) and PVBS states.We find a first-order transition, in contrast to previously studied continuous transitions from the AFM phase into a columnar valence-bond solid (CVBS) phase.The theory of deconfined quantum criticality predicts generic continuous AFM-CVBS and AFM-PVBS transitions, and, in one version of the theory, the two critical order parameters transform under SO(5) symmetry.Emergent SO(5) symmetry has indeed been observed in studies of the AFM-CVBS transition, and here we show that the first-order AFM-PVBS transition also exhibits SO(5) symmetry at the transition point.Such unexpected symmetry of the coexistence state, which implies a lack of energy barriers between the coexisting phases, has recently been observed at other first-order transitions, but the case presented here is the first example with SO(5) symmetry.The extended symmetry may indicate that the transition is connected to a deconfined critical point.We also discuss the first-order transition in the context of a recent proposal of spinons with fracton properties in the PVBS state, concluding that the fracton scenario is unlikely.Furthermore, we discover a novel type of eightfold degenerate VBS phase, arising when the PVBS state breaks a remaining Z 2 symmetry.This second phase transition, which is continuous, implies that the PVBS phase can be regarded as an intermediate "vestigial" phase, a concept recently introduced to describe multistage phase transitions involving a continuous symmetry.Here we construct a six-dimensional order parameter and also introduce a general graph-theoretic approach to describe the two-stage discrete symmetry breaking.We discuss different ways of breaking the symmetries in one or two stages at zero and finite temperatures.In the latter case, we observe fluctuation-induced first-order transitions, which are hallmarks of vestigial phase transitions.We also mention possible connections of the AFM-PVBS transition to the SO(5) theory of high-T c superconductivity.
The theory for the vanishing of N\'eel order in the spin S=1/2 square lattice antiferromagnet has been the focus of attention for many decades. A consensus appears to have emerged in recent numerical studies on the antiferromagnet with first and second neighbor exchange interactions (the ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ model): A gapless spin liquid is present for a narrow window of parameters between the vanishing of the N\'eel order and the onset of a gapped valence bond solid state. We propose a deconfined critical SU(2) gauge theory for a transition into a stable Z2 spin liquid with massless Dirac spinon excitations; on the other side of the critical point, the SU(2) spin liquid (the `$\ensuremath{\pi}$-flux' phase) is presumed to be unstable to confinement to the N\'eel phase. We identify a dangerously irrelevant coupling in the critical SU(2) gauge theory, which contributes a logarithm-squared renormalization. This critical theory is also not Lorentz invariant and weakly breaks the SO(5) symmetry which rotates between the N\'eel and valence bond solid order parameters. We also propose a distinct deconfined critical U(1) gauge theory for a transition into the same gapless Z2 spin liquid; on the other side of the critical point, the U(1) spin liquid (the `staggered flux' phase) is presumed to be unstable to confinement to the valence bond solid. This critical theory has no dangerously irrelevant coupling, dynamic critical exponent $z\ensuremath{\ne}1$, and no SO(5) symmetry. All of these phases and critical points are unified in a SU(2) gauge theory with Higgs fields and fermionic spinons which can naturally realize the observed sequence of phases with increasing J2/J1: N\'eel, gapless Z2 spin liquid, and valence bond solid.
It has been proposed that the deconfined criticality in (2+1)d---the quantum phase transition between a N\'eel antiferromagnet and a valence-bond solid (VBS)---may actually be pseudocritical, in the sense that it is a weakly first-order transition with a generically long correlation length. The underlying field theory of the transition would be a slightly complex (nonunitary) fixed point as a result of fixed points annihilation. This proposal was motivated by existing numerical results from large scale Monte Carlo simulations as well as a conformal bootstrap. However, an actual theory of such a complex fixed point, incorporating key features of the transition such as the emergent SO(5) symmetry, is so far absent. Here we propose a Wess-Zumino-Witten (WZW) nonlinear sigma model with level k=1, defined in $2+\ensuremath{\epsilon}$ dimensions, with target space ${S}^{3+\ensuremath{\epsilon}}$ and global symmetry $\text{SO}(4+\ensuremath{\epsilon})$. This gives a formal interpolation between the deconfined criticality at d=3 and the SU(2)1 WZW theory at d=2 describing the spin-1/2 Heisenberg chain. The theory can be formally controlled, at least to leading order, in terms of the inverse of the WZW level 1/k. We show that at leading order there is a fixed point annihilation at ${d}^{*}\ensuremath{\approx}2.77$, with complex fixed points above this dimension including the physical d=3 case. The pseudocritical properties such as correlation length, scaling dimensions, and the drifts of scaling dimensions as the system size increases, calculated crudely to leading order, are qualitatively consistent with existing numerics.
We develop an optimized continuous-field quantum Monte Carlo (QMC) algorithm to investigate the projected SO(5) nonlinear sigma model with a Wess-Zumino-Witten term, which describes half-filled Dirac fermions in 2+1 space-time dimensions akin to graphene and Yukawa coupled to a quintuplet of compatible mass terms. Our algorithm reduces the computational complexity to O(βNqNϕ2), yielding a speedup of a factor of Nϕ (the number of magnetic fluxes, i.e., system size) relative to prior works [1-4]. This advance enables us to simulate system sizes up to Nϕ=140 on the torus and Nϕ=59 on the sphere, far exceeding the maximum sizes previously accessed, and to map out the universal phase diagram of the model on both geometries. Most notably, we identify and characterize a critical point that separates an SO(5)-broken ordered phase at small coupling from an SO(5)-symmetric disordered phase at large coupling. The critical point becomes multicritical upon the inclusion of terms that break the SO(5) symmetry down to U(1)×SU(2), relevant for the deconfined phase transition between Néel antiferromagnetic and valence-bond-solid orders in quantum magnets. Our finding of a multicritical point in the phase diagram of the SO(5) nonlinear sigma model with Wess-Zumino-Witten term resolves the long-standing open question of its global structure, and our QMC algorithm opens a new avenue for systematic studies of projected Hamiltonians, ranging from correlated flat bands to fractional quantum (anomalous) Hall systems.
Using quantum Monte Carlo simulations, we study effects of disorder on the S=1/2 Heisenberg model with exchange constant J on the square lattice supplemented by multispin interactions Q. It was found recently [L. Liu et al., Phys. Rev. X 8, 041040 (2018)] that the ground state of this $J\ensuremath{-}Q$ model with random couplings undergoes a quantum phase transition from the N\'eel antiferromagnetic state into a randomness-induced spin-liquid-like state that is a close analog to the well known random-singlet (RS) state of the Heisenberg chain with random couplings. This 2D RS state arises from a spontaneously symmetry-broken fourfold degenerate columnar valence-bond solid that is broken up by the disorder into finite domains, with spinons localized at topological defects. The interacting spinons form a critical collective many-body state without magnetic long range order but with the mean spin-spin correlations decaying with distance r as ${r}^{\ensuremath{-}2}$, as in the one-dimensional RS state. The dynamic exponent $z\ensuremath{\ge}2$, varying continuously with the model parameters. In this work, we further investigate the properties of the RS state in the $J\ensuremath{-}Q$ model with random Q couplings. We study the temperature dependence of the specific heat and various susceptibilities for large enough systems to reach the thermodynamic limit. We also analyze the size dependence of the critical magnetic order parameter and its susceptibility in the ground state. For all these quantities, we find consistency with the conventional quantum-critical scaling laws when the condition implied by the ${r}^{\ensuremath{-}2}$ form of the spin correlations is imposed. In particular, all the different quantities can be explained by the same value of the dynamic exponent z at fixed model parameters. We argue that the RS state identified in the $J\ensuremath{-}Q$ model corresponds to a generic renormalization group fixed point that can be reached in many quantum magnets with random couplings and that it has already been observed experimentally.