Abstract Studying strong electron correlations has been an essential driving force for pushing the frontiers of condensed matter physics. In particular, in the vicinity of correlation-driven quantum phase transitions (QPTs), quantum critical fluctuations of multiple degrees of freedom facilitate exotic many-body states and quantum critical behaviours beyond Landau’s framework1. Recently, moiré heterostructures of van der Waals materials have been demonstrated as highly tunable quantum platforms for exploring fascinating, strongly correlated quantum physics2–22. Here we report the observation of tunable quantum criticalities in an experimental simulator of the extended Hubbard model with spin–valley isospins arising in chiral-stacked twisted double bilayer graphene (cTDBG). Scaling analysis shows a quantum two-stage criticality manifesting two distinct quantum critical points as the generalized Wigner crystal transits to a Fermi liquid by varying the displacement field, suggesting the emergence of a critical intermediate phase. The quantum two-stage criticality evolves into a quantum pseudo criticality as a high parallel magnetic field is applied. In such a pseudo criticality, we find that the quantum critical scaling is only valid above a critical temperature, indicating a weak first-order QPT therein. Our results demonstrate a highly tunable solid-state simulator with intricate interplay of multiple degrees of freedom for exploring exotic quantum critical states and behaviours.
Two-dimensional van der Waals materials offer a highly tunable platform for engineering electronic band structures and interactions. By employing techniques such as twisting, gating, or applying pressure, these systems enable precise control over the electronic excitation spectrum. In moir\'e bilayer graphene, the tunability facilitates the transition from a symmetric Dirac semimetal phase through a quantum critical point into an interaction-induced long-range ordered phase with a finite band gap. At charge neutrality, the ordered state proposed to emerge from twist-angle tuning is the Kramers intervalley-coherent insulator. In this case, the transition falls into the quantum universality class of the relativistic Gross-Neveu-XY model in
2+1 dimensions. Here, we refine estimates for the critical exponents characterizing this universality class using an expansion around the lower critical space-time dimension of two. We compute the order-parameter anomalous dimension ${\ensuremath{\eta}}_{\ensuremath{\varphi}}$ and the correlation-length exponent $\ensuremath{\nu}$ at one-loop order, and the fermion anomalous dimension ${\ensuremath{\eta}}_{\ensuremath{\psi}}$ at two-loop order. Combining these results with previous findings from the expansion around the upper critical dimension, we obtain improved estimates for the universal exponents in
2+1 dimensions via Pad\'e interpolation. For
Nf=4 four-component Dirac fermions, relevant to moir\'e bilayer graphene, we estimate $1/\ensuremath{\nu}=0.916(5), {\ensuremath{\eta}}_{\ensuremath{\varphi}}=0.926(13)$, and ${\ensuremath{\eta}}_{\ensuremath{\psi}}=0.0404(13)$. For
Nf=2, potentially relevant to recent tetralayer
WSe2 experiments, the Gross-Neveu-XY fixed point may be unstable owing to a fixed-point collision at
Nf,c, with ${N}_{\text{f,c}}=1+\sqrt{2}+\mathcal{O}(\ensuremath{\epsilon})$ in the expansion around the lower critical dimension.
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.
Motivated by recent experimental studies that have found signatures of a correlated insulator phase and tuning superconductivity in twisted bilayer graphene, we study the temperature-dependent conductivity, the spin correlation and the superconducting pairing correlation within a two-orbital Hubbard model on an emergent honeycomb lattice. Evaluation of the temperature dependence of the conductivity demonstrates that there is a metal-insulator transition, and the Mott phase at strong coupling is accompanied by antiferromagnetic order. The electronic correlation drives a
d+id superconducting pairing to be dominant over a wide filling region. All of the dc conductivity, the spin correlation and the superconductivity are suppressed as the interlayer coupling strength increases, and the critical
Uc for metal-insulator transition is also reduced. Our intensive numerical results reveal that the TBG should be a uniquely tunable platform for exploring correlated states.
Twisted transition metal dichalcogenides offer an exceptionally tunable moiré platform for studying correlation physics beyond conventional condensed matter systems. In particular, the intriguing interplay between the displacement field and the twist angle remains to be fully resolved. In this paper, we use large-scale density matrix renormalization group simulations to study the minimal moiré Hubbard model on a triangular lattice at half-filling, where the displacement field effect is captured by a spin-dependent staggered flux. We find that the displacement field significantly enriches the triangular Hubbard phase diagram in several qualitative ways. It rapidly destabilizes the chiral spin liquid phase beyond a narrow weak-field regime, induces pronounced chiral correlations in the strong-coupling
120∘-antiferromagnetic phase, and stabilizes incommensurate spin-density wave phases at weaker coupling. We further find signatures of a continuous transition between the chiral spin liquid and chiral antiferromagnetic phases at a finite displacement field, potentially driven by spinon condensation. Our results uncover rich displacement-field-driven many-body physics and provide useful guidance for future experiments in moiré superlattice systems.
Twisted van der Waals heterostructures have latterly received prominent attention for their many remarkable experimental properties and the promise that they hold for realizing elusive states of matter in the laboratory. We propose that these systems can, in fact, be used as a robust quantum simulation platform that enables the study of strongly correlated physics and topology in quantum materials. Among the features that make these materials a versatile toolbox are the tunability of their properties through readily accessible external parameters such as gating, straining, packing and twist angle; the feasibility to realize and control a large number of fundamental many-body quantum models relevant in the field of condensed-matter physics; and finally, the availability of experimental readout protocols that directly map their rich phase diagrams in and out of equilibrium. This general framework makes it possible to robustly realize and functionalize new phases of matter in a modular fashion, thus broadening the landscape of accessible physics and holding promise for future technological applications.
Moiré systems provide a rich platform for studies of strong correlation physics. Recent experiments on heterobilayer transition metal dichalcogenide Moiré systems are exciting in that they manifest a relatively simple model system of an extended Hubbard model on a triangular lattice. Inspired by the prospect of the hetero-transition metal dichalcogenide Moiré system's potential as a solid-state-based quantum simulator, we explore the extended Hubbard model on the triangular lattice using the density matrix renormalization group. Specifically, we explore the two-dimensional phase space spanned by the key tuning parameters in the extended Hubbard model, namely, the kinetic energy strength and the further-range Coulomb interaction strengths. We find competition between Fermi fluid, chiral spin liquid, spin density wave, and charge order. In particular, our finding of the optimal further-range interaction for the chiral correlation presents a tantalizing possibility.