Search

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 \textit{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., 7+4 pages, 4+5 figures

While the ground state phase diagram of the correlated flat-band systems have 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\'e materials such as twisted bilayer graphene and transition metal dichalcogenides. Here we show, by means of momentum-space quantum Monte Carlo and exact diagonalization, there exists a unique thermodynamic characteristic for the correlated flat-band models 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 "quantum teleport" the electrons across the gap between different topological bands to restore the broken time-reversal symmetry and give rise to a pronounced enhancement in the charge compressibility. Future experiments, to verify such generic thermodynamic characteristics, are proposed., Comment: 5 pages, 3 figures with supplemental material

Although evidence of inter-valley attraction-mediated by phonon or topological fluctuations is accumulating, the origin of superconductivity in the flat-band quantum moir\'e materials remains an open question. Here, instead of attempting to pinpoint the origin of the superconductivity, we aim at identifying universal properties of moir\'e flat bands that shall emerge in the presence of inter-valley attractions. We show that by matching the interaction strength of inter-valley attraction with intra-valley repulsion, the flat-band limit becomes exactly solvable. Away from the flat-band limit, the system can be simulated via quantum Monte Carlo (QMC) methods without sign problem for any fillings. Combining analytic solutions with large-scale numerical simulations, we show that upon increasing temperature, the superconducting phase melts into a bosonic fluid of Cooper pairs with large/diverging compressibility. In contrast to flat-band attractive Hubbard models, where similar effects arise only for on-site interactions, our study indicates this physics is a universal property of moir\'e flat bands, regardless of microscopic details such as the range of interactions and/or spin-oribt couplings. At higher temperature, the boson fluid phase gives its way to a pseudo gap phase, where some Cooper pairs are torn apart by thermal fluctuations, resulting in fermion density of states inside the gap. Unlike the superconducting transition temperature, which is very sensitive to doping and twisting angles, the gap and the temperature scale of the boson fluid phase and the pseudo gap phase are found to be nearly independent of doping level and/or flat-band bandwidth. The relevance of these phases with experimental discoveries in the flat band quantum moir\'e materials is discussed.

At sufficiently low temperatures, interacting electron systems tend to develop orders. Exceptions are quantum critical point (QCP) and quantum spin liquid (QSL), where fluctuations prevent the highly entangled quantum matter to an ordered state down to the lowest temperature. While the ramification of these states may have appeared in high-temperature superconductors, ultra-cold atoms, frustrated magnets and quantum Moir\'e materials, their unbiased presence remain elusive in microscopic two-dimensional lattice models. Here, we show by means of large-scale quantum Monte Carlo simulations of correlated electrons on the $\pi$-flux square lattice subjected to extended Hubbard interaction, that a Gross-Neveu QCP separating massless Dirac fermions and an columnar valence bond solid at finite interaction, and a possible Dirac QSL at the infinite yet tractable interaction limit emerge in a coherent sequence. These unexpected novel quantum states resides in this simple-looking model, unifying ingredients including emergent symmetry, deconfined fractionalization and the dynamic coupling between emergent matter and gauge fields, will have profound implications both in quantum many-body theory and understanding of the aforementioned experimental systems., Comment: 12 pages, 7 figures

We employ strange correlators to detect 2D subsystem symmetry-protected topological (SSPT) phases which are nontrivial topological phases protected by subsystem symmetries. Specifically, we analytically construct efficient strange correlators in the 2D cluster model in the presence of a uniform magnetic field and then perform the projector Quantum Monte Carlo simulation within the quantum annealing scheme. We find that strange correlators show the long-range correlation in the SSPT phase, from which we define strange order parameters to characterize the topological phase transition between the SSPT phase at low fields and the trivial paramagnetic phase at high fields. Thus, the detection of the fully localized zero modes on the 1D physical boundary of SSPT phase has been transformed into the bulk correlation measurement about the local operators with the periodic boundary condition. We also find interesting spatial anisotropy of a strange correlator, which can be intrinsically traced back to the nature of spatial anisotropy of subsystem symmetries that protect SSPT order in the 2D cluster model. By simulating strange correlators, we, therefore, provide the first unbiased large-scale quantum Monte Carlo simulation on the easy and efficient detection in the SSPT phase and open the avenue of the investigation of the subtle yet fundamental nature of the novel interacting topological phases., 16 pages, 14 figures

The quantum dimer and loop models attract great attentions, partially because the fundamental importance in the phases and phase transitions emerging in these prototypical constrained systems, and partially due to their intimate relevance toward the on-going experiments on Rydberg atom arrays in which the blockade mechanism naturally enforces the local constraint. Here we show, by means of the sweeping cluster quantum Monte Carlo method, the complete ground state phase diagram of the fully packed quantum loop model on the square lattice. We find between the lattice nematic (LN) phase with strong dimer attraction and the staggered phase (SP) with strong dimer repulsion, there emerges a resonating plaquette (RP) phase with off-diagonal translational symmetry breaking. Such a quantum phase is separated from the LN via a first order transition and from the SP by the famous Rokhsar-Kivelson point. Our renormalization group analysis reveals the different flow directions, fully consistent with the order parameter histogram in Monte Carlo simulations. The realization and implication of our phase diagram in Rydberg experiments are proposed., 8 pages, 7 figures

As an unconventional realization of topological orders with an exotic interplay of topology and geometry, fracton (topological) orders feature subextensive topological ground state degeneracy and subdimensional excitations that are movable only within a certain subspace. It has been known in the exactly solvable three-dimensional X-cube model that universally represents the type-I fracton orders, that mobility constraints on subdimensional excitations originate from the absence of spatially deformable string-like operators. To unveil the interplay of topology and geometry, in this paper, we study the dynamical signature in the X-cube model in the presence of external Zeeman fields via large-scale quantum Monte Carlo simulation and stochastic analytic continuation. We compute both real-space correlation functions and dynamic structure factors of subdimensional excitations (i.e., fractons, lineons, and planons) in the fracton phase and their evolution into the trivial paramagnetic phase by increasing external fields. We find in the fracton phase, that the correlation functions and the spectral functions show clear anisotropy exactly caused by the underlying mobility constraints. On the other hand, the external fields successfully induce quantum fluctuations and offer mobility to excitations along the subspace allowed by mobility constraints. These numerical results provide the evolution of a dynamical signature of subdimensional particles in fracton orders, indicating that the mobility constraints on local dynamical properties of subdimensional excitations are deeply related to the existence of fracton topological order. The results will also be helpful in potential experimental identifications in spectroscopy measurements such as neutron scattering and nuclear magnetic resonance., 12 pages, 11 figures

The origin of the pseudogap behavior, found in many high-$T_c$ superconductors, remains one of the greatest puzzles in condensed matter physics. One possible mechanism is fermionic incoherence, which near a quantum critical point allows pair formation but suppresses superconductivity. Employing quantum Monte Carlo simulations of a model of itinerant fermions coupled to ferromagnetic spin fluctuations, represented by a quantum rotor, we report numerical evidence of pseudogap behavior, emerging from pairing fluctuations in a quantum-critical non-Fermi liquid. Specifically, we observe enhanced pairing fluctuations and a partial gap opening in the fermionic spectrum. However, the system remains non-superconducting until reaching a much lower temperature. In the pseudogap regime the system displays a "gap-filling" rather than "gap-closing" behavior, consistent with experimental observations. Our results provide the first unambiguous lattice model realization of a pseudogap state in a strongly correlated system, driven by superconducting fluctuations., 9+14 pages, 4+13 figures

Frustrated magnets host the promises of material realizations of new paradigm of quantum matter, while direct comparison of unbiased model calculations with experimental measurements is still very challenging. Here, we design and implement a protocol of employing many-body computation methodologies for accurate model calculation -- both equilibrium and dynamical properties -- of a frustrated rare-earth magnet TmMgGaO$_4$ (TMGO), which perfectly explains the corresponding experimental findings. Our results confirm TMGO is an ideal realization of triangular-lattice Ising model with an intrinsic transverse field. The magnetic order of TMGO is predicted to melt through two successive Kosterlitz-Thouless (KT) phase transitions, with a floating KT phase in between. The dynamical spectra calculated suggest remnant images of a vanishing magnetic stripe order that represent vortex-antivortex pairs, resembling rotons in a superfluid helium film. TMGO therefore constitutes a rare quantum magnet for realizing KT physics and we further propose experimental detections of its intriguing properties., Comment: 18 pages, 5 figures, supplementary information

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. 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 a Ising transition separating them from the high-$T$ symmetric phases. Due to 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 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 the avenue for controlled many-body investigations on finite-temperature states in the TBG and other quantum moir\'e systems., Comment: 5+15 pages, 4+12 figures

We investigate the ground state phase diagram of an extended Hubbard model with $\pi$-flux hopping term at half-filling on a square lattice, with unbiased large-scale auxiliary-field quantum Monte Carlo simulations. As a function of interaction strength, there emerges an intermediate phase which realizes two interaction-driven quantum critical points, with the first between the Dirac semimetal and an insulating phase of weak valence bond solid (VBS) order, and the second separating the VBS order and an antiferromagnetic insulating phase. These intriguing quantum critical points are respectively bestowed with Gross-Neveu and deconfined quantum criticalities, and the critical exponents $\eta_\text{VBS}=0.6(1)$ and $\eta_\text{AF}=0.58(3)$ at deconfined quantum critical point satisfy the CFT Bootstrap bound. We also investigate the dynamical properties of the spin excitation and find the spin gap open near the first transition and close at the second. The relevance of our findings in realizing deconfined quantum criticality in fermion systems and the implication to lattice models with further extended interactions such as those in quantum Moir\'e systems, are discussed., Comment: 6+2 pages, 5+2 figures

Quantum dimer models are known to host topological quantum spin liquid phases, and it has recently become possible to simulate such models with Rydberg atoms trapped in arrays of optical tweezers. Here, we present large-scale quantum Monte Carlo simulation results on an extension of the triangular lattice quantum dimer model with terms in the Hamiltonian annihilating and creating single dimers. We find distinct odd and even $${{\mathbb{Z}}}_{2}$$ Z 2 spin liquids, along with several phases with no topological order: a staggered crystal, a nematic phase, and a trivial symmetric phase with no obvious broken symmetry. We also present dynamic spectra of the phases, and note implications for experiments on Rydberg atoms.

We summarize the recent developments in the model design and computation for a few representative quantum many-body systems, encompassing quantum critical metals beyond the Hertz-Millis-Moriya framework with pseudogap and superconductivity, SYK non-Fermi-liquid with self-tuned quantum criticality and fluctuation induced superconductivity, and the flat-band quantum Moir\'e lattice models in continuum where the interplay of quantum geometry of flat-band wave function and the long-range Coulomb interactions gives rise to novel insulating phases at integer fillings and superconductivity away from them. Although the narrative choreography seems simple, we show how important the appropriate model design and their tailor-made algorithmic developments -- in other words, the scientific imagination inspired by the corresponding fast experimental developments in the aforementioned systems -- compel us to invent and discover new knowledge and insights in the sport and pastime of quantum many-body research., Comment: Invited review article on quantum many-body model design and computation. 29 pages + 27 figures. Comments and missing references are welcome

We propose a height-conserving quantum dimer model (hQDM) such that the lattice sum of its associated height field is conserved, and that it admits a Rokhsar-Kivelson (RK) point. The hQDM with minimal interaction range on the square lattice exhibits Hilbert space fragmentation and maps exactly to the XXZ spin model on the square lattice in certain Krylov subspaces. We obtain the ground-state phase diagram of hQDM via quantum Monte Carlo simulations, and demonstrate that a large portion of it is within the Krylov subspaces which admit the exact mapping to the XXZ model, with dimer ordered phases corresponding to easy-axis and easy-plane spin orders. At the RK point, the apparent dynamical exponents obtained from the single mode approximation and the height correlation function show drastically different behavior across the Krylov subspaces, exemplifying Hilbert space fragmentation and emergent glassy phenomena., Comment: 5+14 pages, 4+10 figures

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 $\epsilon$-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. Hence both experimental and numerical investigations will not be representative of the infrared fixed point, but of a crossover regime characterized by drifting exponents., Comment: 35 pages, including a 29 page supplemental information section; added notion of quasiuniversality; updated acknowledgements

The experimental discovery of the fractional Hall conductivity in two-dimensional electron gases revealed new types of quantum particles, called anyons, which are beyond bosons and fermions as they possess fractionalized exchange statistics. These anyons are usually studied deep inside an insulating topological phase. It is natural to ask whether such fractionalization can be detected more broadly, say near a phase transition from a conventional to a topological phase. To answer this question, we study a strongly correlated quantum phase transition between a topological state, called a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\mathbb{Z}}}_{2}$$\end{document}Z2 quantum spin liquid, and a conventional superfluid using large-scale quantum Monte Carlo simulations. Our results show that the universal conductivity at the quantum critical point becomes a simple fraction of its value at the conventional insulator-to-superfluid transition. Moreover, a dynamically self-dual optical conductivity emerges at low temperatures above the transition point, indicating the presence of the elusive vison particles. Our study opens the door for the experimental detection of anyons in a broader regime, and has ramifications in the study of quantum materials, programmable quantum simulators, and ultra-cold atomic gases. In the latter case, we discuss the feasibility of measurements in optical lattices using current techniques., Conventional quantum particles can break up into fractionalized excitations under the right conditions; however, their direct experimental observation is challenging. Here, the authors predict strong optical conductivity signatures of such excitations in the vicinity of a topological phase transition.

Employing the recently developed momentum-space quantum Monte Carlo scheme, we study the dynamic response of single-particle and collective excitations in realistic continuum models of twisted bilayer graphene. At charge neutrality, this unbiased numerical method reveals strong competition between different symmetry breaking channels with a leading instability towards the intervalley coherent state. Single-particle spectra indicate that repulsive interactions push the fermion spectral weight away from the Fermi energy and open up an insulating gap. The spectra of collective excitations suggest an approximate valley $SU(2)$ symmetry. At low-energy, long-lived valley waves are observed, which resemble spin waves of Heisenberg ferromagnetism. At high-energy, these sharp modes quickly become over-damped, when their energy reaches the fermion particle-hole continuum., 4 pages, 2 figures with supplemental material

We develop a nonequilibrium increment method to compute the R\'enyi 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 an 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 although 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 positivity condition of the R\'enyi entanglement entropy, which holds for unitary conformal field theories(CFT). Our results unambiguously reveal fundamental differences between DQC and quantum critical points(QCPs) described by unitary CFTs., Comment: published version

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., Comment: 16 pages and 8 figures

We construct a lattice model of topological order (kagome quantum spin liquids) and solve it with unbiased quantum Monte Carlo simulations. A three-stage anyon condensation with two transitions from a $\mathbb Z_2\boxtimes\mathbb Z_2$ topological order to a $\mathbb Z_2$ topological order and eventually to a trivial symmetric phase is revealed. These results provide concrete examples of phase transitions between topological orders in quantum magnets. The designed quantum spin liquid model and its numerical solution offer a playground for further investigations on vestigial anyon condensation., Comment: 7 pages, 9 figures

Based on the path integral formulation of the reduced density matrix, we develop a scheme to overcome the exponential growth of computational complexity in reliably extracting low-lying entanglement spectrum from quantum Monte Carlo simulations. We test the method on the Heisenberg spin ladder with long entangled boundary between two chains and the results support the Li and Haldane's conjecture on entanglement spectrum of topological phase. We then explain the conjecture via the wormhole effect in the path integral and show that it can be further generalized for systems beyond gapped topological phases. Our further simulation results on the bilayer antiferromagnetic Heisenberg model with 2D entangled boundary across the $(2+1)$D O(3) quantum phase transition clearly demonstrate the correctness of the wormhole picture. Finally, we state that since the wormhole effect amplifies the bulk energy gap by a factor of $\beta$, the relative strength of that with respect to the edge energy gap will determine the behavior of low-lying entanglement spectrum of the system., Comment: 9 pages, 5 figures

We consider the effect of the coupling between 2D quantum rotors near an XY ferromagnetic quantum critical point and spins of itinerant fermions. We analyze how this coupling affects the dynamics of rotors and the self-energy of fermions.A common belief is that near a $q=0$ ferromagnetic transition, fermions induce an $\Omega/q$ Landau damping of rotors (i.e., the dynamical critical exponent is $z=3$) and Landau overdamped rotors give rise to non-Fermi liquid fermionic self-energy $\Sigma\propto \omega^{2/3}$. This behavior has been confirmed in previous quantum Monte Carlo (QMC) studies.Here we show that for the XY case the behavior is different.We report the results of large scale quantum Monte Carlo simulations,which show that at small frequencies $z=2$ and $\Sigma\propto \omega^{1/2}$. We argue that the new behavior is associated with the fact that a fermionic spin is by itself not a conserved quantity due to spin-spin coupling to rotors, and a combination of self-energy and vertex corrections replaces $1/q$ in the Landau damping by a constant. We discuss the implication of these results to experiments., Comment: 12 pages, 5 figures

We study scaling behavior of the disorder parameter, defined as the expectation value of a symmetry transformation applied to a finite region, at the deconfined quantum critical point in (2+1)$d$ in the $J$-$Q_3$ model via large-scale quantum Monte Carlo simulations. We show that the disorder parameter for U(1) spin rotation symmetry exhibits perimeter scaling with a logarithmic correction associated with sharp corners of the region, as generally expected for a conformally-invariant critical point. However, for large rotation angle the universal coefficient of the logarithmic corner correction becomes negative, which is not allowed in any unitary conformal field theory. We also extract the current central charge from the small rotation angle scaling, whose value is much smaller than that of the free theory., Comment: 20 pages, 11 figures; v2 improved measurement on disorder operator; v3 additional analysis on data fitting

Simulating lattice gauge theory (LGT) Hamiltonian and its nontrivial states by programmable quantum devices has attracted numerous attention in recent years. Rydberg atom arrays constitute one of the most rapidly developing arenas for quantum simulation and quantum computing. The Z2 LGT and topological order has been realized in experiments while the U(1) LGT is being worked hard on the way. States of LGT have local constraint and are fragmented into several winding sectors with topological protection. It is therefore difficult to prepare a state in certain sector for experiments, and it is also an important task for quantum topological memory. Here, we propose a protocol of sweeping quantum annealing (SQA) for searching the state within target topological sector. With the Monte Carlo method, we show that this SQA has linear time complexity of size with applications to the antiferromagnetic transverse field Ising model, which has emergent U(1) gauge fields. This SQA protocol can be realized easily on quantum simulation platforms such as Rydberg atoms and superconducting circuits. We expect this approach would provide a generic recipe for resolving the topological hindrances in quantum optimization and the preparation of quantum topological state.

We develop a nonequilibrium increment method in quantum Monte Carlo simulations to obtain the Rényi entanglement entropy of various quantum many-body systems with high efficiency and precision. To demonstrate its power, we show the results on a few important yet difficult (2 + 1)d quantum lattice models, ranging from the Heisenberg quantum antiferromagnet with spontaneous symmetry breaking, the quantum critical point with O(3) conformal field theory (CFT) to the toric code $${{\mathbb{Z}}}_{2}$$ Z 2 topological ordered state and the Kagome $${{\mathbb{Z}}}_{2}$$ Z 2 quantum spin liquid model with frustration and multi-spin interactions. In all these cases, our method either reveals the precise CFT data from the logarithmic correction or extracts the quantum dimension in topological order, from the dominant area law in finite-size scaling, with very large system sizes, controlled errorbars, and minimal computational costs. Our method, therefore, establishes a controlled and practical computation paradigm to obtain the difficult yet important universal properties in highly entangled quantum matter.

The fluctuations of massless Dirac fermion can not only turn a first-order bosonic phase transition (in the Landau sense) to a quantum critical point, but also work reversely to enhance the first-order transition itself, depending on the implementation of finite size effects in the coupling corrections. Here, we report a case study of the latter by employing quantum Monte Carlo simulation upon a lattice model in which the bosonic part featuring the Landau-Devonshire first-order phase transition and Yukawa coupled to the Dirac fermions. We find that the parameter range for the first-order phase transition becomes larger as the Yukawa coupling increases and the microscopic mechanism of this phenomena is revealed, at a quantitative level, as the interplay between the critical fluctuations and the finite-size effects. Moreover, the scaling behavior at the separation point between the first-order and the continuous phase transitions is found to belong to the chiral tricritical Gross-Neveu universality. Our result demonstrates that the interplay of massless Dirac fermions, critical fluctuations and the finite size effects could trigger a plethora of interesting phenomena and therefore great care is called for when making generalizations., 9 pages, 6 figures

The Berezinskii-Kosterlitz-Thouless (BKT) mechanism, building upon proliferation of topological defects in 2D systems, is the first example of phase transition beyond the Landau-Ginzburg paradigm of symmetry breaking. Such a topological phase transition has long been sought yet undiscovered directly in magnetic materials. Here, we pin down two transitions that bound a BKT phase in an ideal 2D frustrated magnet TmMgGaO$_4$, via nuclear magnetic resonance under in-plane magnetic fields, which do not disturb the low-energy electronic states and allow BKT fluctuations to be detected sensitively. Moreover, by applying out-of-plane fields, we find a critical scaling behaviour of the magnetic susceptibility expected for the BKT transition. The experimental findings can be explained by quantum Monte Carlo simulations applied on an accurate triangular-lattice Ising model of the compound which hosts a BKT phase. These results provide a concrete example for the BKT phase and offer an ideal platform for future investigations on the BKT physics in magnetic materials., Comment: Version as accepted by Nature Communications, main text: 7 pages, 3 figures; supplementary information: 5 pages, 7 figures

Over the past several years, reliable quantum Monte Carlo results for the charge density wave transition temperature ${T}_{\mathrm{cdw}}$ of the half-filled two-dimensional Holstein model in square and honeycomb lattices have become available for the first time. Exploiting the further development of numerical methodology, here we present results in three dimensions, which are made possible through the use of Langevin evolution of the quantum phonon degrees of freedom. In addition to determining ${T}_{\mathrm{cdw}}$ from the scaling of the charge correlations, we also examine the nature of charge order at general wave vectors for different temperatures, couplings, and phonon frequencies, and the behavior of the spectral function and specific heat.

We review analytical and numerical studies of correlated insulating states in twisted bilayer graphene, focusing on real-space lattice models constructions and their unbiased quantum many-body solutions. We show that by constructing localized Wannier states for the narrow bands, the projected Coulomb interactions can be approximated by interactions of cluster charges with assisted nearest neighbor hopping terms. With the interaction part only, the Hamiltonian is $SU(4)$ symmetric considering both spin and valley degrees of freedom. In the strong coupling limit where the kinetic terms are neglected, the ground states are found to be in the $SU(4)$ manifold with degeneracy. The kinetic terms, treated as perturbation, break this large $SU(4)$ symmetry and propel the appearance of intervalley coherent state, quantum topological insulators and other symmetry-breaking insulating states. We first present the theoretical analysis of moir\'e lattice model construction and then show how to solve the model with large-scale quantum Monte Carlo simulations in an unbiased manner. We further provide potential directions such that from the real-space model construction and its quantum many-body solutions how the perplexing yet exciting experimental discoveries in the correlation physics of twisted bilayer graphene can be gradually understood. This review will be helpful for the readers to grasp the fast growing field of the model study of twisted bilayer graphene., Comment: 15 pages, 12 figures, invited review article

Quantum Monte Carlo (QMC) simulations of correlated electron systems provide unbiased information about system behavior at a quantum critical point (QCP) and can verify or disprove the existing theories of non-Fermi liquid (NFL) behavior at a QCP. However, simulations are carried out at a finite temperature, where quantum-critical features are masked by finite temperature effects. Here we present a theoretical framework within which it is possible to separate thermal and quantum effects and extract the information about NFL physics at $T=0$. We demonstrate our method for a specific example of 2D fermions near a Ising-ferromagnetic QCP. We show that one can extract from QMC data the zero-temperature form of fermionic self-energy $\Sigma (\omega)$ even though the leading contribution to the self-energy comes from thermal effects. We find that the frequency dependence of $\Sigma (\omega)$ agrees well with the analytic form obtained within the Eliashberg theory of dynamical quantum criticality, and obeys $\omega^{2/3}$ scaling at low frequencies. Our results open up an avenue for QMC studies of quantum-critical metals., Comment: 10 pages, 8 figures

Among the quantum many-body models that host anyon excitation and topological orders, quantum dimer models (QDM) provide a unique playground for studying the relation between single-anyon and multi-anyon continuum spectra. However, as the prototypical correlated system with local constraints, the generic solution of QDM at different lattice geometry and parameter regimes is still missing due to the lack of controlled methodologies. Here we obtain, via the newly developed sweeping cluster quantum Monte Carlo algorithm, the excitation spectra in different phases of the triangular lattice QDM. Our results reveal the single vison excitations inside the $Z_2$ quantum spin liquid (QSL) and its condensation towards the $\sqrt{12}\times\sqrt{12}$ valence bond solid (VBS), and demonstrate the translational symmetry fractionalization and emergent O(4) symmetry at the QSL-VBS transition. We find the single vison excitations, whose convolution qualitatively reproduces the dimer spectra, are not free but subject to interaction effects throughout the transition. The nature of the VBS with its O(4) order parameters are unearthed in full scope. Our approach opens the avenue for generic solution of the static and dynamic properties of QDMs and has relevance towards the realization and detection of fractional excitations in programmable quantum simulators., 4+6 pages,3+2 figures

Since the discovery of the pseudogap and Fermi arc states in underdoped cuprates, the understanding of such non-Fermi-liquid states and the associated violation of Luttinger's theorem have been the central theme in correlated electron systems. However, still lacking is a well-accepted theoretical framework to unambiguously explain these metallic states that are clearly beyond Landau's Fermi liquid and Luttinger's theorem of a Fermi surface and electron filling. Here, we design a lattice model of orthogonal metals with fermion and Ising matter fields coupled to topological order and, by solving the model via unbiased quantum Monte Carlo simulation at generic electron fillings, find that the system gives birth to phenomena of the Fermi arc and pseudogap in the single-particle spectrum that go beyond the Luttinger sum rule with broken Fermi surface but no symmetry breaking. The pseudogap and Fermi arcs coexist with a background of a deconfined Z2 gauge field, and we further find that the confinement transition of the gauge field triggers a superconductivity instability and that the hopping of the gauge-neutral fermions brings the "large" Fermi surface back from the Fermi arc state. Our unbiased numerical results provide a concrete model realization and theoretical framework for the coupling between gauge field and fermions and, in the process, generate the rich phenomena of the pseudogap, the Fermi arc, and superconductivity in generic correlated electron systems., 14 pages, 9 figures

Twisted bilayer graphene (TBG) provides a unique framework to elucidate the interplay between strong correlations and topological phenomena in two-dimensional systems. The existence of multiple electronic degrees of freedom -- charge, spin, and valley -- gives rise to a plethora of possible ordered states and instabilities. Identifying which of them are realized in the regime of strong correlations is fundamental to shed light on the nature of the superconducting and correlated insulating states observed in the TBG experiments. Here, we use unbiased, sign-problem-free quantum Monte Carlo simulations to solve an effective interacting lattice model for TBG at charge neutrality. Besides the usual cluster Hubbard-like repulsion, this model also contains an assisted hopping interaction that emerges due to the non-trivial topological properties of TBG. Such a non-local interaction fundamentally alters the phase diagram at charge neutrality, gapping the Dirac cones even for infinitesimally small interaction. As the interaction strength increases, a sequence of different correlated insulating phases emerge, including a quantum valley Hall state with topological edge states, an intervalley-coherent insulator, and a valence bond solid. The charge-neutrality correlated insulating phases discovered here provide the sought-after reference states needed for a comprehensive understanding of the insulating states at integer fillings and the proximate superconducting states of TBG., 15 pages, 9 figures, 2 tables

We systematically test the performance of several Monte Carlo update schemes for the $(2+1)$d XY phase transition of quantum rotor model. By comparing the local Metropolis (LM), LM plus over-relaxation (OR), Wolff-cluster (WC), hybrid Monte Carlo (HM), hybrid Monte Carlo with Fourier acceleration (FA) scheme, it is clear that among the five different update schemes, at the quantum critical point, the WC and FA schemes acquire the smallest autocorrelation time and cost the least amount of CPU hours in achieving the same level of relative error, and FA enjoys a further advantage of easily implementable for more complicated interactions such as the long-range ones. These results bestow one with the necessary knowledge of extending the quantum rotor model, which plays the role of ferromagnetic/antiferromagnetic critical bosons or Z$_2$ topological order, to more realistic and yet challenging models such as Fermi surface Yukawa-coupled to quantum rotor models., Comment: 12 pages, 8 figures

We systematically study the low-temperature specific heats for the two-dimensional kagome antiferromagnet, Cu$_{3}$Zn(OH)$_6$FBr. The specific heat exhibits a $T^{1.7}$ dependence at low temperatures and a shoulder-like feature above it. We construct a microscopic lattice model of $Z_2$ quantum spin liquid and perform large-scale quantum Monte Carlo simulations to show that the above behaviors come from the contributions from gapped anyons and magnetic impurities. Surprisingly, we find the entropy associated with the shoulder decreases quickly with grain size $d$, although the system is paramagnetic to the lowest temperature. While this can be simply explained by a core-shell picture in that the contribution from the interior state disappears near the surface, the 5.9-nm shell width precludes any trivial explanations. Such a large length scale signifies the coherence length of the nonlocality of the quantum entangled excitations in quantum spin liquid candidate, similar to Pippard's coherence length in superconductors. Our approach therefore offers a new experimental probe of the intangible quantum state of matter with topological order., Comment: 6 pages, 3 figures

The coupling between fermionic matter and gauge fields plays a fundamental role in our understanding of nature, while at the same time posing a challenging problem for theoretical modeling. In this situation, controlled information can be gained by combining different complementary approaches. Here, we study a confinement transition in a system of $N_f$ flavors of interacting Dirac fermions charged under a U(1) gauge field in 2+1 dimensions. Using Quantum Monte Carlo simulations, we investigate a lattice model that exhibits a continuous transition at zero temperature between a gapless deconfined phase, described by three-dimensional quantum electrodynamics, and a gapped confined phase, in which the system develops valence-bond-solid order. We argue that the quantum critical point is in the universality class of the QED$_3$-Gross-Neveu-XY model. We study this field theory within a $1/N_f$ expansion in fixed dimension as well as a renormalization group analysis in $4-\epsilon$ space-time dimensions. The consistency between numerical and analytical results is revealed from large to intermediate flavor number., Comment: 10 pages, 9 figures; v2: additional data and explanations, corrected erroneous interpretation of dimer scaling dimension

Magic-angle twisted bilayer graphene has recently become a thriving material platform realizing correlated electron phenomena taking place within its topological flat bands. Several numerical and analytical methods have been applied to understand the correlated phases therein, revealing some similarity with the quantum Hall physics. In this work, we provide a Mott-Hubbard perspective for the TBG system. Employing the large-scale density matrix renormalization group on the lattice model containing the projected Coulomb interactions only, we identify a first-order quantum phase transition between the insulating stripe phase and the quantum anomalous Hall state with the Chern number of $\pm 1$. Our results not only shed light on the mechanism of the quantum anomalous Hall state discovered at three-quarters filling, but also provide an example of the topological Mott insulator, i.e., the quantum anomalous Hall state in the strong coupling limit., Comment: 6+7 pages, 3+7 figures

In recent years, new phases of matter that are beyond the Landau paradigm of symmetry breaking have been accumulating, and to catch up with this fast development, new notions of global symmetry are introduced. Among them, the higher-form symmetry, whose symmetry charges are spatially extended, can be used to describe topologically ordered phases as the spontaneous breaking of the symmetry, and consequently unify the unconventional and conventional phases under the same conceptual framework. However, such conceptual tools have not been put into quantitative tests except for certain solvable models, therefore limiting their usage in the more generic quantum many-body systems. In this work, we study ${\mathbb{Z}}_{2}$ higher-form symmetry in a quantum Ising model, which is dual to the global (zero-form) Ising symmetry. We compute the expectation value of the Ising disorder operator, which is a nonlocal order parameter for the higher-form symmetry, analytically in free scalar theories and through unbiased quantum Monte Carlo simulations for the interacting fixed point in $(2+1)d$. From the scaling form of this extended object, we confirm that the higher-form symmetry is indeed spontaneously broken inside the paramagnetic, or quantum disordered phase (in the Landau sense), but remains symmetric in the ferromagnetic or ordered phase. At the Ising critical point, we find that the disorder operator also obeys a ``perimeter'' law scaling with possibly multiplicative power-law corrections. We discuss examples where both the global zero-form symmetry and the dual higher-form symmetry are preserved, in systems with a codimension-1 manifold of gapless points in momentum space. These results provide nontrivial working examples of higher-form symmetry operators, including the direct computation of one-form order parameter in an interacting conformal field theory, and open the avenue for their generic implementation in quantum many-body systems.

Non-Fermi liquids (NFL) are a class of strongly interacting gapless fermionic systems without long-lived quasiparticle excitations. An important group of NFL model features itinerant fermions coupled to soft bosonic fluctuations near a quantum-critical point (QCP), and are widely believed to capture the essential physics of many unconventional superconductors. However numerically the direct observation of a canonical NFL behavior in such systems, characterized by a power-law form in the Green's function, has been elusive. Here we consider a Sachdev-Ye-Kitaev (SYK)-like model with random Yukawa interaction between critical bosons and fermions (dubbed Yukawa-SYK model). We show it is immune from minus-sign problem and hence can be solved exactly via large-scale quantum Monte Carlo simulation beyond the large-$N$ limit accessible to analytical approaches. Our simulation demonstrates the Yukawa-SYK model features "self-tuned quantum criticality", namely the system is critical independent of the bosonic bare mass. We put these results to test at finite $N$, and our unbiased numerics reveal clear evidence of these exotic quantum-critical NFL properties -- the power-law behavior in Green's function of fermions and bosons -- which propels the theoretical understanding of critical Planckian metals and unconventional superconductors., Comment: 13 pages, 8 figures

The deconfined quantum critical point (DQCP) was originally proposed as a continuous transition between two spontaneous symmetry breaking phases in 2D spin-1/2 systems. While great efforts have been spent on the DQCP for 2D systems, both theoretically and numerically, ambiguities among the nature of the transition are still not completely clarified. Here we shift the focus to a recently proposed 1D incarnation of DQCP in a spin-1/2 chain. By solving it with the variational matrix product state in the thermodynamic limit, a continuous transition between a valence-bond solid phase and a ferromagnetic phase is discovered. The scaling dimensions of various operators are calculated and compared with those from field theoretical description. At the critical point, two emergent O(2) symmetries are revealed, and the associated conserved current operators with exact integer scaling dimensions are determined with scrutiny. Our findings provide the low-dimensional analog of DQCP where unbiased numerical results are in perfect agreement with the controlled field theoretical predictions and have extended the realm of the unconventional phase transition as well as its identification with the advanced numerical methodology., 9 pages + references + appendix, 9 figures

Thanks to the development in quantum Monte Carlo technique, the compact U(1) lattice gauge theory coupled to fermionic matter at (2+1)D is now accessible with large-scale numerical simulations, and the ground state phase diagram as a function of fermion flavor ($N_f$) and the strength of gauge fluctuations is mapped out~\cite{Xiao2018Monte}. Here we focus on the large fermion flavor case ($N_f=8$) to investigate the dynamic properties across the deconfinement-to-confinement phase transition. In the deconfined phase, fermions coupled to the fluctuating gauge field to form U(1) spin liquid with continua in both spin and dimer spectral functions, and in the confined phase fermions are gapped out into valence bond solid phase with translational symmetry breaking and gapped spectra. The dynamical behaviors provide supporting evidence for the existence of the U(1) deconfined phase and could shine light on the nature of the U(1)-to-VBS phase transition which is of the QED$_3$-Gross-Neveu chiral O(2) universality whose properties still largely unknown., 11 pages, 6 figures

Electron-electron interactions play an important role in graphene and related systems and can induce exotic quantum states, especially in a stacked bilayer with a small twist angle(1-7). For bilayer graphene where the two layers are twisted by the 'magic angle', flat band and strong many-body effects lead to correlated insulating states and superconductivity(4-7). In contrast to monolayer graphene, the band structure of untwisted bilayer graphene can be further tuned by a displacement field(8-10), providing an extra degree of freedom to control the flat band that should appear when two bilayers are stacked on top of each other. Here, we report the discovery and characterization of displacement field-tunable electronic phases in twisted double bilayer graphene. We observe insulating states at a half-filled conduction band in an intermediate range of displacement fields. Furthermore, the resistance gap in the correlated insulator increases with respect to the in-plane magnetic fields and we find that the g factor, according to the spin Zeeman effect, is 2, indicating spin polarization at half-filling. These results establish twisted double bilayer graphene as an easily tunable platform for exploring quantum many-body states., Placing two Bernal-stacked graphene bilayers on top of each other with a small twist angle gives correlated states. As the band structure can be tuned by an electric field, this platform is a more varied setting to study correlated electrons.

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., 4+3 pages, 4+2 figures

We report an implementation of the momentum space quantum Monte Carlo (QMC) method on the interaction model for the twisted bilayer graphene (TBG) at integer fillings. The long-range Coulomb repulsion is treated exactly with the flat bands, spin and valley degrees of freedom of electrons taking into account. We prove the absence of the minus sign problem for QMC simulation at integer fillings when either the two valley or the two spin degrees of freedom are considered. By taking the realistic parameters of the twist angle and interlayer tunnelings into the simulation, we benchmark the QMC data with the exact band gap obtained at the chiral limit, to reveal the insulating ground states at the charge neutrality point (CNP). Then, with the exact Green's functions from QMC, we perform stochastic analytic continuation to obtain the first set of single-particle spectral function for the TBG model at CNP. Our momentum space QMC scheme therefore offers the controlled computation pathway for systematic investigation of the electronic states in realistic TBG model at various electron fillings., 4 pages, 2 figures with supplemental material

Recent experiments [J. Guo et al., Phys. Rev. Lett.124,206602 (2020)] on thermodynamic properties of the frustrated layered quantum magnet SrCu$_2$(BO$_3$)$_2$ -- the Shastry-Sutherland material -- have provided strong evidence for a low-temperature phase transition between plaquette-singlet and antiferromagnetic order as a function of pressure. Further motivated by the recently discovered unusual first-order quantum phase transition with an apparent emergent O(4) symmetry of the antiferromagnetic and plaquette-singlet order parameters in a two-dimensional "checkerboard J-Q" quantum spin model [B. Zhao et al., Nat. Phys. 15, 678 (2019)], we here study the same model in the presence of weak inter-layer couplings. Our focus is on the evolution of the emergent symmetry as the system crosses over from two to three dimensions and the phase transition extends from strictly zero temperature in two dimensions up to finite temperature as expected in SrCu$_2$(BO$_3$)$_2$. Using quantum Monte Carlo simulations, we map out the phase boundaries of the plaquette-singlet and antiferromagnetic phases, with particular focus on the triple point where these two order phases meet the paramagnetic phase for given strength of the inter-layer coupling. All transitions are first-order in the neighborhood of the triple points. We show that the emergent O(4) symmetry of the coexistence state breaks down clearly when the interlayer coupling becomes sufficiently large, but for a weak coupling, of the magnitude expected experimentally, the enlarged symmetry can still be observed at the triple point up to significant length scales. Thus, it is likely that the plaquette-singlet to antiferromagnetic transition in SrCu$_2$(BO$_3$)$_2$ exhibits remnants of emergent O(4) symmetry, which should be observable due to additional weakly gapped Goldstone modes., Comment: 14 pages, 8 figures

Recently, quantum simulation of low-dimensional lattice gauge theories (LGTs) has attracted many interests, which may improve our understanding of strongly correlated quantum many-body systems. Here, we propose an implementation to approximate $\mathbb{Z}_2$ LGT on superconducting quantum circuits, where the effective theory is a mixture of a LGT and a gauge-broken term. Using matrix product state based methods, both the ground state properties and quench dynamics are systematically investigated. With an increase of the transverse (electric) field, the system displays a quantum phase transition from a disordered phase to a translational symmetry breaking phase. In the ordered phase, an approximate Gauss law of the $\mathbb{Z}_2$ LGT emerges in the ground state. Moreover, to shed light on the experiments, we also study the quench dynamics, where there is a dynamical signature of the spontaneous translational symmetry breaking. The spreading of the single particle of matter degree is diffusive under the weak transverse field, while it is ballistic with small velocity for the strong field. Furthermore, due to the emergent Gauss law under the strong transverse field, the matter degree can also exhibit confinement dynamics which leads to a strong suppression of the nearest-neighbor hopping. Our results pave the way for simulating the LGT on superconducting circuits, including the quantum phase transition and quench dynamics., 9 Pages, 11 Figures

This review summarizes recent developments in the study of fermionic quantum criticality, focusing on new progress in numerical methodologies, especially quantum Monte Carlo methods, and insights that emerged from recently large-scale numerical simulations. Quantum critical phenomena in fermionic systems have attracted decades of extensive research efforts, partially lured by their exotic properties and potential technology applications and partially awaked by the profound and universal fundamental principles that govern these quantum critical systems. Due to the complex and non-perturbative nature, these systems belong to the most difficult and challenging problems in the study of modern condensed matter physics, and many important fundamental problems remain open. Recently, new developments in model design and algorithm improvements enabled unbiased large-scale numerical solutions to be achieved in the close vicinity of these quantum critical points, which paves a new pathway towards achieving controlled conclusions through combined efforts of theoretical and numerical studies, as well as possible theoretical guidance for experiments in heavy-fermion compounds, Cu-based and Fe-based superconductors, ultra-cold fermionic atomic gas, twisted graphene layers, etc., where signatures of fermionic quantum criticality exist., Comment: 23 pages, 13 figures, invited Topical Review article, comments on content and references are welcome

Orthogonal metal is a new quantum metallic state that conducts electricity but acquires no Fermi surface (FS) or quasiparticles, and hence orthogonal to the established paradigm of Landau's Fermi-liquid (FL). Such a state may hold the key of understanding the perplexing experimental observations of quantum metals that are beyond FL, i.e., dubbed non-Fermi-liquid (nFL), ranging from the Cu- and Fe-based oxides, heavy fermion compounds to the recently discovered twisted graphene heterostructures. However, to fully understand such an exotic state of matter, at least theoretically, one would like to construct a lattice model and to solve it with unbiased quantum many-body machinery. Here we achieve this goal by designing a 2D lattice model comprised of fermionic and bosonic matter fields coupled with dynamic $\mathbb Z_2$ gauge fields, and obtain its exact properties with sign-free quantum Monte Carlo simulations. We find that as the bosonic matter fields become disordered, with the help of deconfinement of the $\mathbb Z_2$ gauge fields, the system reacts with changing its nature from the conventional normal metal with an FS to an orthogonal metal of nFL without FS and quasiparticles and yet still responds to magnetic probe like an FL. Such a quantum phase transition from a normal metal to an orthogonal metal, with its electronic and magnetic spectral properties revealed, is calling for the establishment of new paradigm of quantum metals and their transition with conventional ones., Comment: 7 + 2 pages, 4 + 1 figures

Deconfined quantum critical points govern continuous quantum phase transitions at which fractionalized (deconfined) degrees of freedom emerge. Here we study dynamical signatures of the fractionalized excitations in a quantum magnet (the easy-plane J-Q model) that realize a deconfined quantum critical point with emergent O(4) symmetry. By means of large-scale quantum Monte Carlo simulations and stochastic analytic continuation of imaginary-time correlation functions, we obtain the dynamic spin structure factors in the $S^{x}$ and $S^{z}$ channels. In both channels, we observe broad continua that originate from the deconfined excitations. We further identify several distinct spectral features of the deconfined quantum critical point, including the lower edge of the continuum and its form factor on moving through the Brillouin Zone. We provide field-theoretical and lattice model calculations that explain the overall shapes of the computed spectra, which highlight the importance of interactions and gauge fluctuations to explaining the spectral-weight distribution. We make further comparisons with the conventional Landau O(2) transition in a different quantum magnet, at which no signatures of fractionalization are observed. The distinctive spectral signatures of the deconfined quantum critical point suggest the feasibility of its experimental detection in neutron scattering and nuclear magnetic resonance experiments., 13 pages, 8 figures and 2 appendices

The Mott insulator is the quintessential strongly correlated electronic state. We obtain complete insight into the physics of the two-dimensional Mott insulator by extending the slave-fermion (holon-doublon) description to finite temperatures. We first benchmark its predictions against state-of-the-art quantum Monte Carlo simulations, demonstrating quantitative agreement. Qualitatively, the short-ranged spin fluctuations both induce holon-doublon bound states and renormalize the charge sector to form the Hubbard bands. The Mott gap is understood as the charge gap renormalized downwards by these spin fluctuations. As temperature increases, the Mott gap closes before the charge gap, causing a pseudogap regime to appear naturally during the melting of the Mott insulator., 12 pages, 10 figures