Néel-state to valence-bond-solid transition on the honeycomb lattice: evidence for deconfined criticality

Disorder Operator and Rényi Entanglement Entropy of Symmetric Mass Generation

The "symmetric mass generation" (SMG) quantum phase transition discovered in recent years has attracted great interest from both condensed matter and high energy theory communities. Here, interacting Dirac fermions acquire a gap without condensing any fermion bilinear mass term or any concomitant spontaneous symmetry breaking. It is hence beyond the conventional Gross-Neveu-Yukawa-Higgs paradigm. One important question we address in this Letter is whether the SMG transition corresponds to a true unitary conformal field theory. We employ the sharp diagnosis including the scaling of disorder operator and Rényi entanglement entropy in large-scale lattice model quantum Monte Carlo simulations. Our results strongly suggest that the SMG transition is indeed an unconventional quantum phase transition and it should correspond to a true (2+1)d unitary conformal field theory.

Deconfinement transitions in three-dimensional compact lattice Abelian Higgs models with multiple-charge scalar fields

We investigate the nature of the deconfinement transitions in three-dimensional lattice Abelian Higgs models, in which a complex scalar field of integer charge Q≥2 is minimally coupled with a compact U(1) gauge field. Their phase diagram presents two phases separated by a transition line where static charges q, with q Collapse

Key Words Collapse

MESH Headings Collapse

Grants Collapse

Affiliation(s)

Claudio Bonati Dipartimento di Fisica dell'Università di Pisa and INFN Sezione di Pisa, I-56127 Pisa, Italy
Andrea Pelissetto Dipartimento di Fisica dell'Università di Roma Sapienza and INFN Sezione di Roma I, I-00185 Roma, Italy
Ettore Vicari Dipartimento di Fisica dell'Università di Pisa, I-56127 Pisa, Italy

Collapse

Coulomb-Higgs phase transition of three-dimensional lattice Abelian Higgs gauge models with noncompact gauge variables and gauge fixing

We study the critical behavior of three-dimensional (3D) lattice Abelian Higgs (AH) gauge models with noncompact gauge variables and multicomponent complex scalar fields, along the transition line between the Coulomb and Higgs phases. Previous works that focused on gauge-invariant correlations provided evidence that, for a sufficiently large number of scalar components, these transitions are continuous and associated with the stable charged fixed point of the renormalization-group flow of the 3D AH field theory (scalar electrodynamics), in which charged scalar matter is minimally coupled with an electromagnetic field. Here we extend these studies by considering gauge-dependent correlations of the gauge and matter fields, in the presence of two different gauge fixings, the Lorenz and the axial gauge fixing. Our results for N=25 are definitely consistent with the predictions of the AH field theory and therefore provide additional evidence for the characterization of the 3D AH transitions along the Coulomb-Higgs line as charged transitions in the AH field-theory universality class. Moreover, our results give additional insights on the role of the gauge fixing at charged transitions. In particular, we show that scalar correlations are critical only if a hard Lorenz gauge fixing is imposed.

Nonequilibrium Dynamics of Deconfined Quantum Critical Point in Imaginary Time

Deconfined quantum critical point (DQCP) characterizes a kind of exotic phase transition beyond the usual Landau-Ginzburg-Wilson paradigm. Here we study the nonequilibrium imaginary-time dynamics of the DQCP in the two-dimensional J-Q_{3} model. We explicitly show the deconfinement dynamic process and identify that it is the spinon confinement length, rather than the usual correlation length, that increases proportionally to the time. Moreover, we find that, in the relaxation process, the order parameters of the Néel and the valence-bond-solid orders can be controlled by different length scales, although they satisfy the same equilibrium scaling forms. A dual dynamic scaling theory is then proposed. Our findings not only constitute a new realm of nonequilibrium criticality in DQCP, but also offer a controllable knob by which to investigate the dynamics in strongly correlated systems. Possible realizations in foreseeable quantum computers are also discussed.

Antiferromagnetism Induced by Bond Su-Schrieffer-Heeger Electron-Phonon Coupling: A Quantum Monte Carlo Study

Antiferromagnetism (AFM) such as Néel ordering is often closely related to Coulomb interactions such as Hubbard repulsion in two-dimensional (2D) systems. Whether Néel AFM ordering in two dimensions can be dominantly induced by electron-phonon couplings (EPC) has not been completely understood. Here, by employing numerically exact sign-problem-free quantum Monte Carlo (QMC) simulations, we show that bond Su-Schrieffer-Heeger (SSH) phonons with frequency ω and EPC constant λ can induce AFM ordering for a wide range of phonon frequency ω>ω_{c}. For ω<ω_{c}, a valence-bond-solid (VBS) order appears and there is a direct quantum phase transition between VBS and AFM phases at ω_{c}. The phonon mechanism of the AFM ordering is related to the fact that SSH phonons directly couple to electron hopping whose second-order process can induce an effective AFM spin exchange. Our results shall shed new light on understanding AFM ordering in correlated quantum materials.

Quantum entanglement and criticality in a one-dimensional deconfined quantum critical point

We investigate the ground-state quantum entanglement in a one-dimensional incarnation of deconfined quantum critical point by making use of the finite-size density matrix renormalization group method. We observe two distinct behaviors of two-site entanglement calculated on odd and even bonds, and the difference O_{E} is shown to obey conventional scaling relations for order parameters. Accurate deconfined critical point and associated critical exponents are numerically extracted from finite-size scaling analyses. We further notice a close similarity between O_{E} and the valence-bond-solid order parameter and same observations are also obtained for quantum coherence and trace distance. Furthermore, the deconfined critical point is suggested to possess rich quantum entanglement other than the two-site entanglement from the residual entanglement perspective. Our work explores the critical characteristics of the one-dimensional deconfined quantum critical point from the quantum information aspect and provides insights for its ground-state entanglement structure.

Higher-charge three-dimensional compact lattice Abelian-Higgs models

We consider three-dimensional higher-charge multicomponent lattice Abelian-Higgs (AH) models, in which a compact U(1) gauge field is coupled to an N-component complex scalar field with integer charge q, so that they have local U(1) and global SU(N) symmetries. We discuss the dependence of the phase diagram, and the nature of the phase transitions, on the charge q of the scalar field and the number N≥2 of components. We argue that the phase diagram of higher-charge models presents three different phases, related to the condensation of gauge-invariant bilinear scalar fields breaking the global SU(N) symmetry, and to the confinement and deconfinement of external charge-one particles. The transition lines separating the different phases show different features, which also depend on the number N of components. Therefore, the phase diagram of higher-charge models substantially differs from that of unit-charge models, which undergo only transitions driven by the breaking of the global SU(N) symmetry, while the gauge correlations do not play any relevant role. We support the conjectured scenario with numerical results, based on finite-size scaling analyses of Monte Carlo simuations for doubly charged unit-length scalar fields with small and large number of components, i.e., N=2 and N=25.

Three-dimensional monopole-free CP^{N-1} models

We investigate the phase diagram and the nature of the phase transitions of three-dimensional monopole-free CP^{N-1} models, characterized by a global U(N) symmetry, a U(1) gauge symmetry, and the absence of monopoles. We present numerical analyses based on Monte Carlo simulations for N=2, 4, 10, 15, and 25. We observe a finite-temperature transition in all cases, related to the condensation of a local gauge-invariant order parameter. For N=2 we are unable to draw any definite conclusion on the nature of the transition. The results may be interpreted in terms of either a weak first-order transition or a continuous transition with anomalously large scaling corrections. However, the results allow us to exclude that the transition belongs to the O(3) vector universality class, as it occurs in the standard three-dimensional CP^{1} model without monopole suppression. For N=4,10,and15, the transition is of first order, and significantly weaker than that observed in the presence of monopoles. For N=25 the results are consistent with a conventional continuous transition. We compare our results with the existing literature and with the predictions of different field-theory approaches. They are consistent with the scenario in which the model undergoes continuous transitions for large values of N, including N=∞, in agreement with analytic large-N calculations for the N-component Abelian-Higgs model.

Continuous Easy-Plane Deconfined Phase Transition on the Kagome Lattice

We use large scale quantum Monte Carlo simulations to study an extended Hubbard model of hard core bosons on the kagome lattice. In the limit of strong nearest-neighbor interactions at 1/3 filling, the interplay between frustration and quantum fluctuations leads to a valence bond solid ground state. The system undergoes a quantum phase transition to a superfluid phase as the interaction strength is decreased. It is still under debate whether the transition is weakly first order or represents an unconventional continuous phase transition. We present a theory in terms of an easy plane noncompact CP^{1} gauge theory describing the phase transition at 1/3 filling. Utilizing large scale quantum Monte Carlo simulations with parallel tempering in the canonical ensemble up to 15552 spins, we provide evidence that the phase transition is continuous at exactly 1/3 filling. A careful finite size scaling analysis reveals an unconventional scaling behavior hinting at deconfined quantum criticality.

Dirac Fermions with Competing Orders: Non-Landau Transition with Emergent Symmetry

We consider a model of Dirac fermions in 2+1 dimensions with dynamically generated, anticommuting SO(3) Néel and Z_{2} Kekulé mass terms that permits sign-free quantum Monte Carlo simulations. The phase diagram is obtained from finite-size scaling and includes a direct and continuous transition between the Néel and Kekulé phases. The fermions remain gapped across the transition, and our data support an emergent SO(4) symmetry unifying the two order parameters. While the bare symmetries of our model do not allow for spinon-carrying Z_{3} vortices in the Kekulé mass, the emergent SO(4) invariance permits an interpretation of the transition in terms of deconfined quantum criticality. The phase diagram also features a tricritical point at which the Néel, Kekulé, and semimetallic phases meet. The present sign-free approach can be generalized to a variety of other mass terms and thereby provides a new framework to study exotic critical phenomena.

Fermion-induced quantum critical points

A unified theory of quantum critical points beyond the conventional Landau-Ginzburg-Wilson paradigm remains unknown. According to Landau cubic criterion, phase transitions should be first-order when cubic terms of order parameters are allowed by symmetry in the Landau-Ginzburg free energy. Here, from renormalization group analysis, we show that second-order quantum phase transitions can occur at such putatively first-order transitions in interacting two-dimensional Dirac semimetals. As such type of Landau-forbidden quantum critical points are induced by gapless fermions, we call them fermion-induced quantum critical points. We further introduce a microscopic model of SU(N) fermions on the honeycomb lattice featuring a transition between Dirac semimetals and Kekule valence bond solids. Remarkably, our large-scale sign-problem-free Majorana quantum Monte Carlo simulations show convincing evidences of a fermion-induced quantum critical points for N = 2, 3, 4, 5 and 6, consistent with the renormalization group analysis. We finally discuss possible experimental realizations of the fermion-induced quantum critical points in graphene and graphene-like materials.Quantum phase transitions are governed by Landau-Ginzburg theory and the exceptions are rare. Here, Li et al. propose a type of Landau-forbidden quantum critical points induced by gapless fermions in two-dimensional Dirac semimetals.

Pyrocurrent anomalies and intrinsic magnetodielectric behavior near room temperature in Li2Ni2Mo3O12, a compound with distorted honeycomb and spin-chains

Keeping current interests to identify materials with intrinsic magnetodielectric behaviour near room temperature and with novel pyroelectric current anomalies, we report temperature and magnetic-field dependent behavior of complex dielectric permittivity and pyroelectric current for an oxide, Li₂Ni₂Mo₃O₁₂, containing magnetic ions with (distorted) honey-comb and chain arrangement and ordering magnetically below 8 K. The dielectric data reveal the existence of relaxor ferroelectricity behaviour in the range 160-240 K and there are corresponding Raman mode anomalies as well in this temperature range. Pyrocurrent behavior is also consistent with this interpretation, with the pyrocurrent peak-temperature interestingly correlating with the poling temperature. ⁷Li NMR offer an evidence for crystallographic disorder intrinsic to this compound and we therefore conclude that such a disorder is apparently responsible for the randomness of local electric field leading to relaxor ferroelectric property. Another observation of emphasis is that there is a notable decrease in the dielectric constant with the application of magnetic field to the tune of about -2.4% at 300 K, with the magnitude varying marginally with temperature. Small loss factor values validate the intrinsic behaviour of the magnetodielectric effect at room temperature.

Translation-Invariant Parent Hamiltonians of Valence Bond Crystals

We present a general method to construct translation-invariant and SU(2) symmetric antiferromagnetic parent Hamiltonians of valence bond crystals (VBCs). The method is based on a canonical mapping transforming S=1/2 spin operators into a bilinear form of a new set of dimer fermion operators. We construct parent Hamiltonians of the columnar and the staggered VBCs on the square lattice, for which the VBC is an eigenstate in all regimes and the exact ground state in some region of the phase diagram. We study the departure from the exact VBC regime upon tuning the anisotropy by means of the hierarchical mean field theory and exact diagonalization on finite clusters. In both Hamiltonians, the VBC phase extends over the exact regime and transits to a columnar antiferromagnet (CAFM) through a window of intermediate phases, revealing an intriguing competition of correlation lengths at the VBC-CAFM transition. The method can be readily applied to construct other VBC parent Hamiltonians in different lattices and dimensions.

Quantum phase diagram of distorted J 1 - J 2 Heisenberg S = 1/2 antiferromagnet in honeycomb lattice: a modified spin wave study

Using the modified spin wave method, we study the [Formula: see text] Heisenberg model with first and second neighbor antiferromagnetic exchange interactions. For a symmetric S  =  1/2 model, with the same couplings for all the equivalent neighbors, we find three phases in terms of the frustration parameter [Formula: see text]: (1) a commensurate collinear ordering with staggered magnetization (Néel.I state) for [Formula: see text], (2) a magnetically gapped disordered state for [Formula: see text], preserving all the symmetries of the Hamiltonian and lattice, which by definition is a quantum spin liquid (QSL) state and (3) a commensurate collinear ordering in which two out of the three nearest neighbor magnetizations are antiparallel and the remaining pair are parallel (Néel.II state), for [Formula: see text]. We also explore the phase diagram of a distorted [Formula: see text] model with S  =  1/2. Distortion is introduced as an inequality of one nearest neighbor coupling with the other two. This yields a richer phase diagram by the appearance of a new gapped QSL, a gapless QSL and also a valence bond crystal phase in addition to the previous three phases found for the undistorted model.

Necessary Condition for Emergent Symmetry from the Conformal Bootstrap

We use the conformal bootstrap program to derive the necessary conditions for emergent symmetry enhancement from discrete symmetry (e.g., Z_{n}) to continuous symmetry [e.g., U(1)] under the renormalization group flow. In three dimensions, in order for Z_{2} symmetry to be enhanced to U(1) symmetry, the conformal bootstrap program predicts that the scaling dimension of the order parameter field at the infrared conformal fixed point must satisfy Δ_{1}>1.08. We also obtain the similar necessary conditions for Z_{3} symmetry with Δ_{1}>0.580 and Z_{4} symmetry with Δ_{1}>0.504 from the simultaneous conformal bootstrap analysis of multiple four-point functions. As applications, we show that our necessary conditions impose severe constraints on the nature of the chiral phase transition in QCD, the deconfinement criticality in Néel valence bond solid transitions, and anisotropic deformations in critical O(n) models. We prove that some fixed points proposed in the literature are unstable under the perturbation that cannot be forbidden by the discrete symmetry. In these situations, the second-order phase transition with enhanced symmetry cannot happen.

Orbitally induced hierarchy of exchange interactions in the zigzag antiferromagnetic state of honeycomb silver delafossite Ag3Co2SbO6

We report the revised crystal structure, static and dynamic magnetic properties of quasi-two dimensional honeycomb-lattice silver delafossite Ag3Co2SbO6. The magnetic susceptibility and specific heat data are consistent with the onset of antiferromagnetic long range order at low temperatures with Néel temperature TN ∼ 21.2 K. In addition, the magnetization curves revealed a field-induced (spin-flop type) transition below TN in moderate magnetic fields. The GGA+U calculations show the importance of the orbital degrees of freedom, which maintain a hierarchy of exchange interaction in the system. The strongest antiferromagnetic exchange coupling was found in the shortest Co-Co pairs and is due to direct and superexchange interaction between the half-filled xz + yz orbitals pointing directly to each other. The other four out of six nearest neighbor exchanges within the cobalt hexagon are suppressed, since for these bonds the active half-filled orbitals turned out to be parallel and do not overlap. The electron spin resonance (ESR) spectra reveal a broad absorption line attributed to the Co(2+) ion in an octahedral coordination with an average effective g-factor g = 2.40 ± 0.05 at room temperature and show strong divergence of the ESR parameters below ∼150 K, which implies an extended region of short-range correlations. Based on the results of magnetic and thermodynamic studies in applied fields, we propose a magnetic phase diagram for the new honeycomb-lattice delafossite.

Kernel method for corrections to scaling

Scaling analysis, in which one infers scaling exponents and a scaling function in a scaling law from given data, is a powerful tool for determining universal properties of critical phenomena in many fields of science. However, there are corrections to scaling in many cases, and then the inference problem becomes ill-posed by an uncontrollable irrelevant scaling variable. We propose a new kernel method based on Gaussian process regression to fix this problem generally. We test the performance of the new kernel method for some example cases. In all cases, when the precision of the example data increases, inference results of the new kernel method correctly converge. Because there is no limitation in the new kernel method for the scaling function even with corrections to scaling, unlike in the conventional method, the new kernel method can be widely applied to real data in critical phenomena.

Plaquette ordered phase and quantum phase diagram in the spin-1/2 J(1)-J(2) square Heisenberg model

We study the spin-1/2 Heisenberg model on the square lattice with first- and second-neighbor antiferromagnetic interactions J(1) and J(2), which possesses a nonmagnetic region that has been debated for many years and might realize the interesting Z(2) spin liquid. We use the density matrix renormalization group approach with explicit implementation of SU(2) spin rotation symmetry and study the model accurately on open cylinders with different boundary conditions. With increasing J(2), we find a Néel phase and a plaquette valence-bond (PVB) phase with a finite spin gap. From the finite-size scaling of the magnetic order parameter, we estimate that the Néel order vanishes at J(2)/J(1)≃0.44. For 0.5<J(2)/J(1)<0.61, we find dimer correlations and PVB textures whose decay lengths grow strongly with increasing system width, consistent with a long-range PVB order in the two-dimensional limit. The dimer-dimer correlations reveal the s-wave character of the PVB order. For 0.44<J(2)/J(1)<0.5, spin order, dimer order, and spin gap are small on finite-size systems, which is consistent with a near-critical behavior. The critical exponents obtained from the finite-size spin and dimer correlations could be compatible with the deconfined criticality in this small region. We compare and contrast our results with earlier numerical studies.

Vector boson excitations near deconfined quantum critical points

We show that the Néel states of two-dimensional antiferromagnets have low energy vector boson excitations in the vicinity of deconfined quantum critical points. We compute the universal damping of these excitations arising from spin-wave emission. Detection of such a vector boson will demonstrate the existence of emergent topological gauge excitations in a quantum spin system.

Fate of CPN-1 fixed points with q monopoles

We present an extensive quantum Monte Carlo study of the Néel to valence-bond solid (VBS) phase transition on rectangular- and honeycomb-lattice SU(N) antiferromagnets in sign-problem-free models. We find that in contrast to the honeycomb lattice and previously studied square-lattice systems, on the rectangular lattice for small N, a first-order Néel-VBS transition is realized. On increasing N≥4, we observe that the transition becomes continuous and with the same universal exponents as found on the honeycomb and square lattices (studied here for N=5, 7, 10), providing strong support for a deconfined quantum critical point. Combining our new results with previous numerical and analytical studies, we present a general phase diagram of the stability of CPN-1 fixed points with q monopoles.