Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2607.27745

Trion Excitations in Twisted Bilayer Graphene: A Quantum Monte Carlo Study

Shibo Shan, Cheng Huang, Patrick Ledwith, Zi Yang Meng

cond-mat.str-elquant-phSubmitted Jul 30, 2026• Announced Jul 31, 2026
arXivHTMLPDF

Abstract

Determining the nature of charge carriers is a fundamental goal in the study of strongly correlated electron systems. Here, we employ the continuous-field momentum-space quantum Monte Carlo method to reveal exotic "Dirac trion" excitations in the finite-temperature normal state of twisted bilayer graphene. While the ground state is a symmetry-breaking insulator with gapped ($\sim$ 20 meV) electron-like excitations, we show that a small temperature ($\sim$ 3 meV), well below the interaction scale, drives the system into a strongly fluctuating symmetric normal state. We demonstrate that this normal state hosts gapless excitations consisting of three-particle bound states, two electrons and one hole, that are exactly orthogonal to the higher-energy electrons at the zero-momentum gapless point. These Dirac trions have the remarkable property of being arbitrarily light despite being composed of heavy constituents, and their spectra can be easily tuned by varying the twist angle and interlayer hopping strength. Our unbiased quantum many-body computation sheds light on the Dirac trions in a projected correlated flat-band setting and opens the door for further investigation of many-body excitations in strongly correlated topological bands beyond Landau levels.

Comments
8 pages, 4 figures in the main text, 5 pages, 1 figure in the Supplemental Material

Related papers

Momentum Space Quantum Monte Carlo on Twisted Bilayer Graphene

Xu Zhang, Gaopei Pan, Yi Zhang +2 more

We report an implementation of the momentum space quantum Monte Carlo (QMC) method on the interaction model for the twisted bilayer graphene (TBG). 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 when either the two valleys 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.

arXiv:2105.07010•Chinese Physics Letters•Jul 1, 2021•75 citations

Fermionic Monte Carlo Study of a Realistic Model of Twisted Bilayer Graphene

Johannes Stephan Hofmann, Eslam Khalaf, Ashvin Vishwanath +2 more

The rich phenomenology of twisted bilayer graphene (TBG) near the magic angle is believed to arise from electron correlations in topological flat bands.An unbiased approach to this problem is highly desirable, but also particularly challenging, given the multiple electron flavors, the topological obstruction to defining tight-binding models, and the long-ranged Coulomb interactions.While numerical simulations of realistic models have thus far been confined to zero temperature, typically excluding some spin or valley species, analytic progress has relied on fixed point models away from the realistic limit.Here, we present unbiased Monte Carlo simulations of realistic models of magic-angle TBG at charge neutrality.We establish the absence of a sign problem for this model in a momentumspace approach and describe a computationally tractable formulation that applies even on breaking chiral symmetry and including band dispersion.Our results include (i) the emergence of an insulating Kramers intervalley coherent ground state in competition with a correlated semimetal phase, (ii) detailed temperature evolution of order parameters and electronic spectral functions that reveal a "pseudogap" regime, in which gap features are established at a higher temperature than the onset of order, and (iii) predictions for electronic tunneling spectra and their evolution with temperature.Our results pave the way towards uncovering the physics of magic-angle graphene through exact simulations of over a hundred electrons across a wide temperature range.

arXiv:2105.12112•Physical Review X•Mar 30, 2022•73 citations

Angle-tuned Gross-Neveu quantum criticality in twisted bilayer graphene

Cheng Huang, Nikolaos Parthenios, Maksim Ulybyshev +4 more

Abstract Twisted bilayber graphene (TBG) has revealed fascinating quantum many-body states, and unstrained chargeneutral TBG at 1.08° is understood, both experimentally and theoretically , as a correlated insulator due to the interplay of long-range Coulomb interactions and the quantum metrics of flat bands. However, the fate of the ground state as one continuously tunes the twisted angle is still largely unknown. Here, theoretically, by employing a newly developed momentum-space continuous-field quantum Monte Carlo method that can fully take into account the long-range Coulomb interactions and quantum metrics of flat bands with system sizes that were not possible before , we show that charge-neutral TBG realizes an angle-tuned quantum phase transition with Gross-Neveu criticality at Θc ~ 1.20(1)°, from the analyses of single-particle spectra, free energy and the order parameter of inter-valley coherence. This shows the exciting possibility that by tuning away from 1.08°, pristine TBG can mimic a fundamental phase transition of Dirac fermions.

arXiv:2412.11382•Nature Communications•Aug 4, 2025•15 citations

Correlation effects in magic-angle twisted bilayer graphene: An auxiliary-field quantum Monte Carlo study

Zhi-Yu Xiao, Shiwei Zhang

Magic-angle twisted bilayer graphene (MATBG) presents a fascinating platform for investigating the effects of electron interactions in topological flat bands. The Bistritzer-MacDonald (BM) model provides a simplified quantitative description of the flat bands. Introducing long-range Coulomb interactions leads to an interacting BM (IBM) Hamiltonian, a momentum-space continuum description which offers a very natural starting point for many-body studies of MATBG. Accurate and reliable many-body computations in the IBM model are challenging, however, and have been limited mostly to special fillings or smaller lattice sizes. We employ a state-of-the-art auxiliary-field quantum Monte Carlo (AFQMC) method to study the IBM model, which constrains the sign problem to enable accurate treatment of large system sizes. We determine ground-state properties and quantify errors compared to mean-field theory calculations. Our calculations identify correlated metal states and their competition with the insulating Kramers intervalley-coherent state at both half-filling and charge neutrality. Additionally, we investigate one- and three-quarter fillings, and examine the effect of many-body corrections beyond single Slater determinant solutions. We discuss the effect that details of the IBM Hamiltonian have on the results, including different forms of double-counting corrections, and the need to establish and precisely specify many-body Hamiltonians to allow more direct and quantitative comparisons with experiments in MATBG.

arXiv:2405.17808•Physical Review Research•Jan 27, 2025•3 citations

Fractional Chern Insulators in Twisted Bilayer Optical Lattices

Yan-Bin Yang, Jiong-Hao Wang, Shicheng Ma +1 more

Twisted bilayer materials provide a versatile platform for realizing correlated topological states. Motivated by the recent experimental realization of atomic Bose-Einstein condensates in twisted bilayer optical lattices, we investigate the emergence of fractional topological phases in such systems. For experimentally realistic parameters, the single-particle spectrum hosts a pair of quasi-degenerate nearly flat moiré bands. We show using self-consistent Hartree-Fock calculations that atomic interactions spontaneously break time-reversal symmetry, lift the quasi-degeneracy, and generate an isolated nearly flat band with a nonzero Chern number. At fractional filling, exact diagonalization of the Hamiltonian projected to this band reveals a threefold degenerate ground-state manifold and a characteristic particle entanglement spectrum, providing strong evidence for a fractional Chern insulator. The resulting many-body gap is of order \(10^{-3}E_R\), corresponding to a nanokelvin energy scale accessible to state-of-the-art cold-atom experiments. Our scheme requires neither spin-orbit coupling, as in transition metal dichalcogenides, nor a finite magnetic field, as in existing twisted-bilayer-graphene experiments, providing a highly tunable route to strongly correlated topological phases in twisted bilayer optical lattices.

arXiv:2609.09733•Sep 9, 2026•0 citations

Realization of topological Mott insulator in a twisted bilayer graphene lattice model

Bin-Bin Chen, Yuan Da Liao, Ziyu Chen +4 more

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

arXiv:2011.07602•Nature Communications•Sep 16, 2021•73 citations
Determining the nature of charge carriers is a fundamental goal in the study of strongly correlated electron systems. Here, we employ the continuous-field momentum-space quantum Monte Carlo method to reveal exotic "Dirac trion" excitations in the finite-temperature normal state of twisted bilayer graphene. While the ground state is a symmetry-breaking insulator with gapped (∼\sim∼ 20 meV) electron-like excitations, we show that a small temperature (∼\sim∼ 3 meV), well below the interaction scale, drives the system into a strongly fluctuating symmetric normal state. We demonstrate that this normal state hosts gapless excitations consisting of three-particle bound states, two electrons and one hole, that are exactly orthogonal to the higher-energy electrons at the zero-momentum gapless point. These Dirac trions have the remarkable property of being arbitrarily light despite being composed of heavy constituents, and their spectra can be easily tuned by varying the twist angle and interlayer hopping strength. Our unbiased quantum many-body computation sheds light on the Dirac trions in a projected correlated flat-band setting and opens the door for further investigation of many-body excitations in strongly correlated topological bands beyond Landau levels.
We report an implementation of the momentum space quantum Monte Carlo (QMC) method on the interaction model for the twisted bilayer graphene (TBG). 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 when either the two valleys 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.
The rich phenomenology of twisted bilayer graphene (TBG) near the magic angle is believed to arise from electron correlations in topological flat bands.An unbiased approach to this problem is highly desirable, but also particularly challenging, given the multiple electron flavors, the topological obstruction to defining tight-binding models, and the long-ranged Coulomb interactions.While numerical simulations of realistic models have thus far been confined to zero temperature, typically excluding some spin or valley species, analytic progress has relied on fixed point models away from the realistic limit.Here, we present unbiased Monte Carlo simulations of realistic models of magic-angle TBG at charge neutrality.We establish the absence of a sign problem for this model in a momentumspace approach and describe a computationally tractable formulation that applies even on breaking chiral symmetry and including band dispersion.Our results include (i) the emergence of an insulating Kramers intervalley coherent ground state in competition with a correlated semimetal phase, (ii) detailed temperature evolution of order parameters and electronic spectral functions that reveal a "pseudogap" regime, in which gap features are established at a higher temperature than the onset of order, and (iii) predictions for electronic tunneling spectra and their evolution with temperature.Our results pave the way towards uncovering the physics of magic-angle graphene through exact simulations of over a hundred electrons across a wide temperature range.
Abstract Twisted bilayber graphene (TBG) has revealed fascinating quantum many-body states, and unstrained chargeneutral TBG at 1.08° is understood, both experimentally and theoretically , as a correlated insulator due to the interplay of long-range Coulomb interactions and the quantum metrics of flat bands. However, the fate of the ground state as one continuously tunes the twisted angle is still largely unknown. Here, theoretically, by employing a newly developed momentum-space continuous-field quantum Monte Carlo method that can fully take into account the long-range Coulomb interactions and quantum metrics of flat bands with system sizes that were not possible before , we show that charge-neutral TBG realizes an angle-tuned quantum phase transition with Gross-Neveu criticality at Θc ~ 1.20(1)°, from the analyses of single-particle spectra, free energy and the order parameter of inter-valley coherence. This shows the exciting possibility that by tuning away from 1.08°, pristine TBG can mimic a fundamental phase transition of Dirac fermions.
Magic-angle twisted bilayer graphene (MATBG) presents a fascinating platform for investigating the effects of electron interactions in topological flat bands. The Bistritzer-MacDonald (BM) model provides a simplified quantitative description of the flat bands. Introducing long-range Coulomb interactions leads to an interacting BM (IBM) Hamiltonian, a momentum-space continuum description which offers a very natural starting point for many-body studies of MATBG. Accurate and reliable many-body computations in the IBM model are challenging, however, and have been limited mostly to special fillings or smaller lattice sizes. We employ a state-of-the-art auxiliary-field quantum Monte Carlo (AFQMC) method to study the IBM model, which constrains the sign problem to enable accurate treatment of large system sizes. We determine ground-state properties and quantify errors compared to mean-field theory calculations. Our calculations identify correlated metal states and their competition with the insulating Kramers intervalley-coherent state at both half-filling and charge neutrality. Additionally, we investigate one- and three-quarter fillings, and examine the effect of many-body corrections beyond single Slater determinant solutions. We discuss the effect that details of the IBM Hamiltonian have on the results, including different forms of double-counting corrections, and the need to establish and precisely specify many-body Hamiltonians to allow more direct and quantitative comparisons with experiments in MATBG.
Twisted bilayer materials provide a versatile platform for realizing correlated topological states. Motivated by the recent experimental realization of atomic Bose-Einstein condensates in twisted bilayer optical lattices, we investigate the emergence of fractional topological phases in such systems. For experimentally realistic parameters, the single-particle spectrum hosts a pair of quasi-degenerate nearly flat moiré bands. We show using self-consistent Hartree-Fock calculations that atomic interactions spontaneously break time-reversal symmetry, lift the quasi-degeneracy, and generate an isolated nearly flat band with a nonzero Chern number. At fractional filling, exact diagonalization of the Hamiltonian projected to this band reveals a threefold degenerate ground-state manifold and a characteristic particle entanglement spectrum, providing strong evidence for a fractional Chern insulator. The resulting many-body gap is of order 10−3ER10^{-3}E_R10−3ER​, corresponding to a nanokelvin energy scale accessible to state-of-the-art cold-atom experiments. Our scheme requires neither spin-orbit coupling, as in transition metal dichalcogenides, nor a finite magnetic field, as in existing twisted-bilayer-graphene experiments, providing a highly tunable route to strongly correlated topological phases in twisted bilayer optical lattices.
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 ±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.