BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Zi Yang Meng

Zi Yang Meng contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.str-elcond-mat.stat-mechhep-latquant-phcond-mat.mes-hallhep-thcond-mat.supr-concond-mat.quant-gascond-mat.dis-nnphysics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
34works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

34 published item(s)

preprint2026arXiv

Chiral Graviton Modes in Fermionic Fractional Chern Insulators

Chiral graviton modes are hallmark collective excitations of Fractional Quantum Hall (FQH) liquids. However, their existence on the lattice, where continuum symmetries that protect them from decay are lost, is still an open and urgent question, especially considering the recent advances in the realization of Fractional Chern Insulators (FCI) in transition metal dichalcogenides and rhombohedral pentalayer graphene. Here we present a comprehensive theoretical and numerical study of graviton-modes in fermionic FCI, and thoroughly demonstrate their existence. We first derive a lattice stress tensor operator in the context of the fermionic Harper-Hofstadter(HH) model which captures the graviton in the flat band limit. Importantly, we discover that such lattice stress-tensor operators are deeply connected to lattice quadrupolar density correlators, readily generalizable to generic Chern bands. We then explicitly show the adiabatic connection between FQH and FCI chiral graviton modes by interpolating from a low flux HH model to a Checkerboard lattice model that hosts a topological flat band. In particular, using state-of-the-art matrix product state and exact diagonalization simulations, we provide strong evidence that chiral graviton modes are long-lived excitations in FCIs despite the lack of continuous symmetries and the scattering with a two-magnetoroton continuum. By means of a careful finite-size analysis, we show that the lattice generates a finite but small intrinsic decay rate for the graviton mode. We discuss the relevance of our results for the exploration of graviton modes in FCI phases realized in solid state settings, as well as cold atom experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Discrete Shift and Polarization from Response to Symmetry Defects in Interacting Topological Phases

We extend the previous study of extracting crystalline symmetry-protected topological invariants to the correlated regime. We construct the interacting Hofstadter model defined on square lattice with the rotation and translation symmetry defects: disclination and dislocation. The model realizes Chern insulator and the charge density wave state as one tunes interactions. Employing the density matrix renormalization group (DMRG) method, we calculate the excess charge around the defects and find that the topological invariants remain quantized in both phases, with the topological quantity extracted to great precision. This study paves the way for utilizing matrix product state, and potentially other quantum many-body computation methods, to efficiently study crystalline symmetry defects on 2D interacting lattice systems.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Precise computation of universal corner entanglement entropy at 2+1 dimension: From Ising to Gaussian quantum critical points

Computing the subleading logarithmic term in the entanglement entropy (EE) of (2+1)d quantum many-body systems remains a significant challenge, despite its central role in revealing universal information about quantum states and quantum critical points (QCPs). Building on recent algorithmic advances that enable the stable calculation of EE as an exponential observable~\cite{zhouIncremental2024,zhangIntegral2024,liaoExtracting2024}, we develop a {\it bubble basis} projector quantum Monte Carlo (QMC) algorithm to precisely and efficiently compute the universal corner of EE at QCPs in a (2+1)d square-lattice transverse-field Ising model augmented with a four-body interaction. Turning on this interaction allows us to trace an Ising critical line, reaching the tricritical point, and then a line of first-order phase transition. In (2+1)d, the tricritical point is described by the Gaussian theory, where a theoretical calculation of the corner logarithmic term in the 2nd Rényi entropy term is available~\cite{UniversalCasini2007}. Our QMC results are in quantitative agreement with this theoretical value, providing a highly nontrivial benchmark of the algorithm. Furthermore, we also study the Rényi EE at the Ising critical line and on the first-order transition line, obtaining results consistent with theoretical expectations. These findings establish the long-sought connection between the universal values of an exactly solvable limit and those of a strongly correlated regime at (2+1)d.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Spectra of Magnetoroton and Chiral Graviton Modes of Fractional Chern Insulator

Employing the state-of-the-art time-dependent variational principle (TDVP) algorithm, we compute the spectra of charge-neutral excitations in the $ν=1/2$ (bosonic) \updated{ and $1/3$ (fermionic) fractional Chern insulator (FCI)} on the Haldane honeycomb lattice model. The magnetoroton visualized from the dynamic density structure factor acquires a minimum gap at finite momentum that can go soft with increasing interaction and give rise to a charge density wave (CDW) at the same wavevector. As the system approaches the FCI-to-CDW transition point, we observe a pronounced sharpening of the roton mode, suggesting that the magnetoroton behaves more like a quasiparticle as it softens. Notably, this occurs while the single-particle gap remains finite. Besides the magnetoroton at finite momentum, we also construct quadrupolar chiral operators in a discrete lattice and resolve the chiral graviton mode around the $Γ$ point of the Brillouin zone. Furthermore, we show the different chiralities of the gravitons of FCIs with opposite-sign Hall conductance for the first time.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Dirac fermions with plaquette interactions. I. SU(2) phase diagram with Gross-Neveu and deconfined quantum criticalities

We investigate the ground state phase diagram of an extended Hubbard model with $π$-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 $η_\text{VBS}=0.6(1)$ and $η_\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é systems, are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Dynamical properties of collective excitations in twisted bilayer Graphene

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Evolution of Dynamical Signature in the X-cube Fracton Topological Order

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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é systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fermion sign bounds theory in quantum Monte Carlo simulation

Sign problem in fermion quantum Monte Carlo (QMC) simulation appears to be an extremely hard problem. Traditional lore passing around for years tells people that when there is a sign problem, the average sign in QMC simulation approaches zero exponentially fast with the space-time volume of the configurational space. We, however, analytically show this is not always the case and manage to find physical bounds for the average sign. Our understanding is based on a direct connection between the sign bounds and a well-defined partition function of reference system and could distinguish when the bounds have the usual exponential scaling, and when they are bestowed on an algebraic scaling at low temperature limit. We analytically explain such algebraic sign problems found in flat band moiré lattice models at low temperature limit. At finite temperature, a domain size argument based on sign bounds also explains the connection between sign behavior and finite temperature phase transition. Sign bounds, as a well-defined observable, may have ability to ease or even make use of the sign problem.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Height-conserving quantum dimer models

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Measuring Rényi entanglement entropy with high efficiency and precision in quantum Monte Carlo simulations

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 $\Z_2$ topological ordered state and the Kagome $\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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Monte Carlo study of the pseudogap and superconductivity emerging from quantum magnetic fluctuations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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. Hence both experimental and numerical investigations will not be representative of the infrared fixed point, but of a crossover regime characterized by drifting exponents.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Network-Initialized Monte Carlo Based on Generative Neural Networks

We design generative neural networks that generate Monte Carlo configurations with complete absence of autocorrelation from which only short Markov chains are needed before making measurements for physical observables, irrespective of the system locating at the classical critical point, fermionic Mott insulator, Dirac semimetal, or quantum critical point. We further propose a network-initialized Monte Carlo scheme based on such neural networks, which provides independent samplings and can accelerate the Monte Carlo simulations by significantly reducing the thermalization process. We demonstrate the performance of our approach on the two-dimensional Ising and fermion Hubbard models, and expect it can systematically speed up the Monte Carlo simulations especially for the very challenging many-electron problems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Possible Dirac quantum spin liquid in a kagome quantum antiferromagnet YCu$_3$(OH)$_6$Br$_2$[Br$_{x}$(OH)$_{1-x}$]

We studied the magnetic properties of YCu$_3$(OH)$_6$Br$_2$[Br$_{1-x}$(OH)$_{x}$] ($x$ = 0.33), where Cu$^{2+}$ ions form two-dimensional kagome layers. There is no magnetic order down to 50 mK while the Curie-Weiss temperature is on the order of -100 K. At zero magnetic field, the low-temperature specific heat shows a $T^2$ dependence. Above 2 T, a linear temperature dependence term in specific heat emerges, and the value of $γ= C/T$ increases linearly with the field. Furthermore, the magnetic susceptibility tends to a constant value at $T = 0$. Our results suggest that the magnetic ground state of YCu$_3$(OH)$_6$Br$_2$[Br$_{1-x}$(OH)$_{x}$] is consistent with a Dirac quantum-spin-liquid state with a linearly dispersing spinon strongly coupled to an emergent gauge field, which has long been theoretically proposed as a candidate ground state in the two-dimensional kagome Heisenberg antiferromagnetic system.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quantum Simulation of Lattice Gauge Theories on Superconducting Circuits: Quantum Phase Transition and Quench Dynamics

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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 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ényi 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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Solving quantum rotor model with different Monte Carlo techniques

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The dynamical exponent of a quantum critical itinerant ferromagnet: a Monte Carlo study

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 $Ω/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 $Σ\propto ω^{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 $Σ\propto ω^{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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Thermodynamic characteristic for correlated flat-band system with quantum anomalous Hall ground state

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é 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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Topological Disorder Parameter

We introduce a many-body topological invariant, called the topological disorder parameter (TDP), to characterize gapped quantum phases with global internal symmetry in (2+1)d. TDP is defined as the constant correction that appears in the ground state expectation value of a partial symmetry transformation applied to a connected spatial region $M$, the absolute value of which scales generically as $\exp(-αl+γ)$ where $l$ is the perimeter of $M$ and $γ$ is the TDP. Motivated by a topological quantum field theory interpretation of the operator, we show that $e^γ$ can be related to the quantum dimension of the symmetry defect, and provide a general formula for $γ$ when the entanglement Hamiltonian of the topological phase can be described by a (1+1)d conformal field theory (CFT). A special case of TDP is equivalent to the topological Rényi entanglement entropy when the symmetry is the cyclic permutation of the replica of the gapped phase. We then investigate several examples of lattice models of topological phases, both analytically and numerically, in particular when the assumption of having a CFT edge theory is not satisfied. We also consider an example of partial translation symmetry in Wen's plaquette model and show that the result can be understood using the edge CFT. Our results establish a new tool to detect quantum topological order.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Correlated Insulating Phases in the Twisted Bilayer Graphene

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é 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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Correlation-induced insulating topological phases at charge neutrality in twisted bilayer graphene

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Fermion enhanced first-order phase transition and chiral Gross-Neveu tricritical point

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Vestigial anyon condensation in kagome quantum spin liquids

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Confinement transition in the QED$_3$-Gross-Neveu-XY universality class

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-ε$ space-time dimensions. The consistency between numerical and analytical results is revealed from large to intermediate flavor number.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Correlated states in twisted double bilayer graphene

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. For bilayer graphene where the two layers are twisted by a "magic angle", flat band and strong many-body effects lead to correlated insulating states and superconductivity. In contrast to monolayer graphene, the band structure of untwisted bilayer graphene can be further tuned by a displacement field, 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 such 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 spin Zeeman effect is ~2, indicating spin polarization at half filling. These results establish the twisted double bilayer graphene as an easily tunable platform for exploring quantum many-body states.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Designer Monte Carlo Simulation for Gross-Neveu Transition

In this manuscript, we study quantum criticality of Dirac fermions via large-scale numerical simulations, focusing on the Gross-Neveu-Yukawa(GNY) chiral-Ising quantum critical point with critical bosonic modes coupled with Dirac fermions. We show that finite-size effects at this quantum critical point can be efficiently minimized via model design, which maximizes the ultraviolet cutoff and at the same time places the bare control parameters closer to the nontrivial fixed point to better expose the critical region. Combined with the efficient self-learning quantum Monte Carlo algorithm, which enables non-local update of the bosonic field, we find that moderately-large system size (up to $16\times 16$) is already sufficient to produce robust scaling behavior and critical exponents.The conductance of the Dirac fermions is also calculated and its frequency dependence is found to be consistent with the scaling behavior predicted by the conformal field theory. The methods and model-design principles developed for this study can be generalized to other fermionic QCPs, and thus provide a promising direction for controlled studies of strongly-correlated itinerant systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Identification of non-Fermi liquid fermionic self-energy from quantum Monte Carlo data

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 $Σ(ω)$ even though the leading contribution to the self-energy comes from thermal effects. We find that the frequency dependence of $Σ(ω)$ agrees well with the analytic form obtained within the Eliashberg theory of dynamical quantum criticality, and obeys $ω^{2/3}$ scaling at low frequencies. Our results open up an avenue for QMC studies of quantum-critical metals.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Kosterlitz-Thouless Melting of Magnetic Order in the Triangular Quantum Ising Material TmMgGaO$_4$

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnetism of Finite Graphene Samples: Mean-Field Theory compared with Exact Diagonalization and Quantum Monte Carlo Simulation

The magnetic properties of graphene on finite geometries are studied using a self-consistent mean-field theory of the Hubbard model. This approach is known to predict ferromagnetic edge states close to the zig-zag edges in single-layer graphene quantum dots and nanoribbons. In order to assess the accuracy of this method, we perform complementary exact diagonalization and quantum Monte Carlo simulations. We observe good quantitative agreement for all quantities investigated provided that the Coulomb interaction is not too strong.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Metal to Orthogonal Metal Transition

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantum Many-Body Simulations of the 2D Fermi-Hubbard Model in Ultracold Optical Lattices

Understanding quantum many-body states of correlated electrons is one main theme in modern condensed matter physics. Given that the Fermi-Hubbard model, the prototype of correlated electrons, has been recently realized in ultracold optical lattices, it is highly desirable to have controlled numerical methodology to provide precise finite-temperature results upon doping, to directly compare with experiments. Here, we demonstrate the exponential tensor renormalization group (XTRG) algorithm [Phys. Rev. X 8, 031082 (2018)], complemented with independent determinant quantum Monte Carlo (DQMC) offer a powerful combination of tools for this purpose. XTRG provides full and accurate access to the density matrix and thus various spin and charge correlations, down to unprecedented low temperature of few percents of the fermion tunneling energy scale. We observe excellent agreement with ultracold fermion measurements at both half-filling and finite-doping, including the sign-reversal behavior in spin correlations due to formation of magnetic polarons, and the attractive hole-doublon and repulsive hole-hole pairs that are responsible for the peculiar bunching and antibunching behavior of the antimoments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantum phases of SrCu2(BO3)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 inter-layer 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.

0 citations0 reviewsNo code linkedOpen access