Phases of (2+1)D SO(5) Nonlinear Sigma Model with a Topological Term on a Sphere: Multicritical Point and Disorder Phase.

Novel critical phenomena beyond the Landau-Ginzburg-Wilson paradigm have been long sought after. Among many candidate scenarios, the deconfined quantum critical point (DQCP) constitutes the most fascinating one, and its lattice model realization has been debated over the past two decades. Here we apply the spherical Landau level regularization upon the exact (2+1)D SO(5) nonlinear sigma model with a topological term to study the potential DQCP therein. We perform a density matrix renormalization group (DMRG) simulation with SU(2)_{spin}×U(1)_{charge}×U(1)_{angular-momentum} symmetries explicitly implemented. Using crossing point analysis for the critical properties of the DMRG data, accompanied by quantum Monte Carlo simulations, we accurately obtain the comprehensive phase diagram of the model and find various novel quantum phases, including Néel, ferromagnet (FM), valence bond solid (VBS), valley polarized (VP) states and a gapless quantum disordered phase occupying an extended area of the phase diagram. The VBS-disorder and Néel-disorder transitions are continuous with non-Wilson-Fisher exponents. Our results show the VBS and Néel states are separated by either a weakly first-order transition or the disordered region with a multicritical point in between, thus opening up more interesting questions on the two-decade long debate on the nature of the DQCP.

Thermodynamic Response and Neutral Excitations in Integer and Fractional Quantum Anomalous Hall States Emerging from Correlated Flat Bands.

Integer and fractional Chern insulators have been extensively explored in correlated flat band models. Recently, the prediction and experimental observation of fractional quantum anomalous Hall (FQAH) states with spontaneous time-reversal symmetry breaking have garnered attention. While the thermodynamics of integer quantum anomalous Hall (IQAH) states have been systematically studied, our theoretical knowledge on thermodynamic properties of FQAH states has been severely limited. Here, we delve into the general thermodynamic response and collective excitations of both IQAH and FQAH states within the paradigmatic flat Chern-band model with remote band considered. Our key findings include (i) in both ν=1 IQAH and ν=1/3 FQAH states, even without spin fluctuations, the charge-neutral collective excitations would lower the onset temperature of these topological states, to a value significantly smaller than the charge gap, due to band mixing and multiparticle scattering; (ii) by employing large-scale thermodynamic simulations in FQAH states in the presence of strong interband mixing between C=±1 bands, we find that the lowest collective excitations manifest as the zero-momentum excitons in the IQAH state, whereas in the FQAH state, they take the form of magnetorotons with finite momentum; (iii) the unique charge oscillations in FQAH states are exhibited with distinct experimental signatures, which we propose to detect in future experiments.

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.

Emergent Glassy Behavior in a Kagome Rydberg Atom Array.

We present large-scale quantum Monte Carlo simulation results on a realistic Hamiltonian of kagome-lattice Rydberg atom arrays. Although the system has no intrinsic disorder, intriguingly, our analyses of static and dynamic properties on large system sizes reveal emergent glassy behavior in a region of parameter space located between two valence bond solid phases. The extent of this glassy region is demarcated using the Edwards-Anderson order parameter, and its phase transitions to the two proximate valence bond solids-as well as the crossover towards a trivial paramagnetic phase-are identified. We demonstrate the intrinsically slow (imaginary) time dynamics deep inside the glassy phase and discuss experimental considerations for detecting such a quantum disordered phase with numerous nearly degenerate local minima. Our proposal paves a new route to the study of real-time glassy phenomena and highlights the potential for quantum simulation of a distinct phase of quantum matter beyond solids and liquids in current-generation Rydberg platforms.

Thermodynamic Characteristic for a Correlated Flat-Band System with a Quantum Anomalous Hall Ground State.

While the ground-state phase diagram of the correlated flat-band systems has been intensively investigated, the dynamic and thermodynamic properties of such lattice models are less explored, but it is the latter which is most relevant to the experimental probes (transport, quantum capacitance, and spectroscopy) of the quantum moiré materials such as twisted bilayer graphene. Here we show, by means of momentum-space quantum Monte Carlo and exact diagonalization, in chiral limit there exists a unique thermodynamic characteristic for the correlated flat-band model with interaction-driven quantum anomalous Hall (QAH) ground state, namely, the transition from the QAH insulator to the metallic state takes place at a much lower temperature compared with the zero-temperature single-particle gap generated by the long-range Coulomb interaction. Such low transition temperature comes from the proliferation of excitonic particle-hole excitations, which transfers the electrons across the gap between different topological bands to restore the broken time-reversal symmetry and gives rise to a pronounced enhancement in the charge compressibility. Future experiments, to verify such generic thermodynamic characteristics, are proposed.

Intrinsic Nonlinear Hall Effect and Gate-Switchable Berry Curvature Sliding in Twisted Bilayer Graphene.

Though the observation of the quantum anomalous Hall effect and nonlocal transport response reveals nontrivial band topology governed by the Berry curvature in twisted bilayer graphene, some recent works reported nonlinear Hall signals in graphene superlattices that are caused by the extrinsic disorder scattering rather than the intrinsic Berry curvature dipole moment. In this Letter, we report a Berry curvature dipole induced intrinsic nonlinear Hall effect in high-quality twisted bilayer graphene devices. We also find that the application of the displacement field substantially changes the direction and amplitude of the nonlinear Hall voltages, as a result of a field-induced sliding of the Berry curvature hotspots. Our Letter not only proves that the Berry curvature dipole could play a dominant role in generating the intrinsic nonlinear Hall signal in graphene superlattices with low disorder densities, but also demonstrates twisted bilayer graphene to be a sensitive and fine-tunable platform for second harmonic generation and rectification.

Fermion Disorder Operator at Gross-Neveu and Deconfined Quantum Criticalities.

The fermion disorder operator has been shown to reveal the entanglement information in 1D Luttinger liquids and 2D free and interacting Fermi and non-Fermi liquids emerging at quantum critical points (QCPs) [W. Jiang et al., arXiv:2209.07103]. Here we study, by means of large-scale quantum Monte Carlo simulation, the scaling behavior of the disorder operator in correlated Dirac systems. We first demonstrate the logarithmic scaling behavior of the disorder operator at the Gross-Neveu (GN) chiral Ising and Heisenberg QCPs, where consistent conformal field theory (CFT) content of the GN-QCP in its coefficient is found. Then we study a 2D monopole-free deconfined quantum critical point (DQCP) realized between a quantum-spin Hall insulator and a superconductor. Our data point to negative values of the logarithmic coefficients such that the DQCP does not correspond to a unitary CFT. Density matrix renormalization group calculations of the disorder operator on a 1D DQCP model also detect emergent continuous symmetries.

Exciton Proliferation and Fate of the Topological Mott Insulator in a Twisted Bilayer Graphene Lattice Model.

A topological Mott insulator (TMI) with spontaneous time-reversal symmetry breaking and nonzero Chern number has been discovered in a real-space effective model for twisted bilayer graphene (TBG) at 3/4 filling in the strong coupling limit [1]. However, the finite temperature properties of such a TMI state remain illusive. In this work, employing the state-of-the-art thermal tensor network and the perturbative field-theoretical approaches, we obtain the finite-T phase diagram and the dynamical properties of the TBG model. The phase diagram includes the quantum anomalous Hall and charge density wave phases at low T, and an Ising transition separating them from the high-T symmetric phases. Because of the proliferation of excitons-particle-hole bound states-the transitions take place at a significantly reduced temperature than the mean-field estimation. The exciton phase is accompanied with distinctive experimental signatures in such as in charge compressibilities and optical conductivities close to the transition. Our work explains the smearing of the many-electron state topology by proliferating excitons and opens an avenue for controlled many-body investigations on finite-temperature states in the TBG and other quantum moiré systems.

Scaling of Entanglement Entropy at Deconfined Quantum Criticality.

We develop a nonequilibrium increment method to compute the Rényi entanglement entropy and investigate its scaling behavior at the deconfined critical (DQC) point via large-scale quantum Monte Carlo simulations. To benchmark the method, we first show that, at a conformally invariant critical point of O(3) transition, the entanglement entropy exhibits universal scaling behavior of area law with logarithmic corner corrections, and the obtained correction exponent represents the current central charge of the critical theory. Then we move on to the deconfined quantum critical point, where we still observe similar scaling behavior, but with a very different exponent. Namely, the corner correction exponent is found to be negative. Such a negative exponent is in sharp contrast with the positivity condition of the Rényi entanglement entropy, which holds for unitary conformal field theories (CFTs). Our results unambiguously reveal fundamental differences between DQC and quantum critical points described by unitary CFTs.

Nematic Quantum Criticality in Dirac Systems.

We investigate nematic quantum phase transitions in two different Dirac fermion models. The models feature twofold and fourfold, respectively, lattice rotational symmetries that are spontaneously broken in the ordered phase. Using negative-sign-free quantum Monte Carlo simulations and an ε-expansion renormalization group analysis, we show that both models exhibit continuous phase transitions. In contrast to generic Gross-Neveu dynamical mass generation, the quantum critical regime is characterized by large velocity anisotropies, with fixed-point values being approached very slowly. Both experimental and numerical investigations will not be representative of the infrared fixed point, but of a quasiuniversal regime where the drift of the exponents tracks the velocity anisotropy.

Itinerant quantum critical point with fermion pockets and hotspots.

Metallic quantum criticality is among the central themes in the understanding of correlated electronic systems, and converging results between analytical and numerical approaches are still under review. In this work, we develop a state-of-the-art large-scale quantum Monte Carlo simulation technique and systematically investigate the itinerant quantum critical point on a 2D square lattice with antiferromagnetic spin fluctuations at wavevector [Formula: see text]-a problem that resembles the Fermi surface setup and low-energy antiferromagnetic fluctuations in high-Tc cuprates and other critical metals, which might be relevant to their non-Fermi-liquid behaviors. System sizes of [Formula: see text] ([Formula: see text]) are comfortably accessed, and the quantum critical scaling behaviors are revealed with unprecedented high precision. We found that the antiferromagnetic spin fluctuations introduce effective interactions among fermions and the fermions in return render the bare bosonic critical point into a different universality, different from both the bare Ising universality class and the Hertz-Mills-Moriya RPA prediction. At the quantum critical point, a finite anomalous dimension [Formula: see text] is observed in the bosonic propagator, and fermions at hotspots evolve into a non-Fermi liquid. In the antiferromagnetically ordered metallic phase, fermion pockets are observed as the energy gap opens up at the hotspots. These results bridge the recent theoretical and numerical developments in metallic quantum criticality and can serve as the stepping stone toward final understanding of the 2D correlated fermions interacting with gapless critical excitations.

Amplitude Mode in Quantum Magnets via Dimensional Crossover.

We investigate the amplitude (Higgs) mode associated with longitudinal fluctuations of the order parameter at the continuous spontaneous symmetry breaking phase transition. In quantum magnets, due to the fast decay of the amplitude mode into low-energy Goldstone excitations, direct observation of this mode represents a challenging task. By focusing on a quasi-one-dimensional geometry, we circumvent the difficulty and investigate the amplitude mode in a system of weakly coupled spin chains with the help of quantum Monte Carlo simulations, stochastic analytic continuation, and a chain-mean field approach combined with a mapping to the field-theoretic sine-Gordon model. The amplitude mode is observed to emerge in the longitudinal spin susceptibility in the presence of a weak symmetry-breaking staggered field. A conventional measure of the amplitude mode in higher dimensions, the singlet bond mode, is found to appear at a lower than the amplitude mode frequency. We identify these two excitations with the second (first) breather of the sine-Gordon theory, correspondingly. In contrast to higher-dimensional systems, the amplitude and bond order fluctuations are found to carry significant spectral weight in the quasi-1D limit.

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.

Valence Bond Orders at Charge Neutrality in a Possible Two-Orbital Extended Hubbard Model for Twisted Bilayer Graphene.

An extended Hubbard model on a honeycomb lattice with two orbitals per site at charge neutrality is investigated with unbiased large-scale quantum Monte Carlo simulations. The Fermi velocity of the Dirac fermions is renormalized as the cluster charge interaction increases, until a mass term emerges and a quantum phase transition from Dirac semimetal to valence bond solid (VBS) insulator is established. The quantum critical point is discovered to belong to the 3D N=4 Gross-Neveu chiral XY universality with the critical exponents obtained at high precision. Further enhancement of the interaction drives the system into two different VBS phases, the properties and transition between them are also revealed. Since the model is related to magic-angle twisted bilayer graphene, our results may have relevance towards the symmetry breaking order at the charge neutrality point of the material, and associate the wide range of universal strange metal behavior around it with quantum critical fluctuations.

Charge-Density-Wave Transitions of Dirac Fermions Coupled to Phonons.

The spontaneous generation of charge-density-wave order in a Dirac fermion system via the natural mechanism of electron-phonon coupling is studied in the framework of the Holstein model on the honeycomb lattice. Using two independent and unbiased quantum Monte Carlo methods, the phase diagram as a function of temperature and coupling strength is determined. It features a quantum critical point as well as a line of thermal critical points. Finite-size scaling appears consistent with fermionic Gross-Neveu-Ising universality for the quantum phase transition and bosonic Ising universality for the thermal phase transition. The critical temperature has a maximum at intermediate couplings. Our findings motivate experimental efforts to identify or engineer Dirac systems with sufficiently strong and tunable electron-phonon coupling.

Role of Noether's Theorem at the Deconfined Quantum Critical Point.

Noether's theorem is one of the fundamental laws of physics, relating continuous symmetries and conserved currents. Here we explore the role of Noether's theorem at the deconfined quantum critical point (DQCP), which is a quantum phase transition beyond the Landau-Ginzburg-Wilson paradigm. It was expected that a larger continuous symmetry could emerge at the DQCP, which, if true, should lead to conserved current at low energy. By identifying the emergent current fluctuation in the spin excitation spectra, we can quantitatively study the current-current correlation in large-scale quantum Monte Carlo simulations. Our results reveal the conservation of the emergent current, as signified by the vanishing anomalous dimension of the current operator, and hence provide supporting evidence for the emergent symmetry at the DQCP. Our study demonstrates an elegant yet practical approach to detect emergent symmetry by probing the spin excitation, which could potentially guide the ongoing experimental search for the DQCP in quantum magnets.

Propagation and Localization of Collective Excitations on a 24-Qubit Superconducting Processor.

Superconducting circuits have emerged as a powerful platform of quantum simulation, especially for emulating the dynamics of quantum many-body systems, because of their tunable interaction, long coherence time, and high-precision control. Here in experiments, we construct a Bose-Hubbard ladder with a ladder array of 20 qubits on a 24-qubit superconducting processor. We investigate theoretically and demonstrate experimentally the dynamics of single- and double-excitation states with distinct behaviors, indicating the uniqueness of the Bose-Hubbard ladder. We observe the linear propagation of photons in the single-excitation case, satisfying the Lieb-Robinson bounds. The double-excitation state, initially placed at the edge, localizes; while placed in the bulk, it splits into two single-excitation modes spreading linearly toward two boundaries, respectively. Remarkably, these phenomena, studied both theoretically and numerically as unique properties of the Bose-Hubbard ladder, are represented coherently by pairs of controllable qubits in experiments. Our results show that collective excitations, as a single mode, are not free. This work paves the way to simulation of exotic logic particles by subtly encoding physical qubits and exploration of rich physics by superconducting circuits.

Dynamics of Topological Excitations in a Model Quantum Spin Ice.

We study the quantum spin dynamics of a frustrated XXZ model on a pyrochlore lattice by using large-scale quantum Monte Carlo simulation and stochastic analytic continuation. In the low-temperature quantum spin ice regime, we observe signatures of coherent photon and spinon excitations in the dynamic spin structure factor. As the temperature rises to the classical spin ice regime, the photon disappears from the dynamic spin structure factor, whereas the dynamics of the spinon remain coherent in a broad temperature window. Our results provide experimentally relevant, quantitative information for the ongoing pursuit of quantum spin ice materials.

Dynamical Signature of Symmetry Fractionalization in Frustrated Magnets.

The nontriviality of quantum spin liquids (QSLs) typically manifests in the nonlocal observables that signify their existence; however, this fact actually casts a shadow on detecting QSLs with experimentally accessible probes. Here, we provide a solution by unbiasedly demonstrating a dynamical signature of anyonic excitations and symmetry fractionalization in QSLs. Employing large-scale quantum Monte Carlo simulation and stochastic analytic continuation, we investigate the extended XXZ model on the kagome lattice, and find out that, across the phase transitions from Z_{2} QSLs to different symmetry breaking phases, spin spectral functions can reveal the presence and condensation of emergent anyonic spinon and vison excitations, in particular, the translational symmetry fractionalization of the latter, which can be served as the dynamical signature of the seemingly ephemeral QSLs in spectroscopic techniques such as inelastic neutron or resonance (inelastic) x-ray scatterings.

Quantum Spin Liquid with Even Ising Gauge Field Structure on Kagome Lattice.

Employing large-scale quantum Monte Carlo simulations, we study the extended XXZ model on the kagome lattice. A Z_{2} quantum spin liquid phase with effective even Ising gauge field structure emerges from the delicate balance among three symmetry-breaking phases including stripe solid, staggered solid, and ferromagnet. This Z_{2} spin liquid is stabilized by an extended interaction related to the Rokhsar-Kivelson potential in the quantum dimer model limit. The phase transitions from the staggered solid to a spin liquid or ferromagnet are found to be first order and so is the transition between the stripe solid and ferromagnet. However, the transition between a spin liquid and ferromagnet is found to be continuous and belongs to the 3D XY^{*} universality class associated with the condensation of spinons. The transition between a spin liquid and stripe solid appears to be continuous and associated with the condensation of visons.