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:2405.17808 · Physical Review Research

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

Zhi-Yu Xiao, Shiwei Zhang

Published Jan 27, 2025• 3 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevresearch.7.013103

Related papers

Many-Body Ground State Manifold of Flat Band Interacting Hamiltonian for Magic Angle Twisted Bilayer Graphene

Kevin D. Stubbs, Michael Ragone, Allan H. MacDonald +1 more

Abstract At a magic relative twist angle, magic angle twisted bilayer graphene (MATBG) has an octet of flat bands that can host strong correlation physics when partially filled. A key theoretical discovery in MATBG is the existence of ferromagnetic Slater determinants as exact ground states of the corresponding flat band interacting (FBI) Hamiltonian. The FBI Hamiltonian describes the behavior of electrons that interact with each other in a high-dimensional space, and is constructed from the band structure of the non-interacting Bistritzer–MacDonald model at the chiral limit. A key property of the FBI Hamiltonian for MATBG is that it is frustration-free and can be written as a sum of non-commuting terms. In this work, we provide a complete characterization of the ground state manifold of the FBI Hamiltonian, proving that it is precisely the linear span of such ferromagnetic Slater determinants.

arXiv:2503.20060•Communications in Mathematical Physics•Sep 27, 2026•1 citation

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

Exact Ground State of Interacting Electrons in Magic Angle Graphene

Simon Becker, Lin Lin, Kevin D. Stubbs

Abstract One of the most remarkable theoretical findings in magic angle twisted bilayer graphene (TBG) is the emergence of ferromagnetic Slater determinants as exact ground states for the interacting Hamiltonian at the chiral limit. This discovery provides an explanation for the correlated insulating phase which has been experimentally observed at half filling. This work is the first mathematical study of interacting models in magic angle graphene systems. These include not only TBG but also TBG-like systems featuring four flat bands per valley, and twisted trilayer graphene systems with equal twist angles. We identify symmetries of the chiral limit of the Bistritzer-MacDonald Hamiltonian that are responsible for characterizing the Hartree-Fock ground states as zero energy many-body ground states. Furthermore, for a general class of Hamiltonian, we establish criteria that the ferromagnetic Slater determinants are the unique ground states within the class of uniformly half-filled, translation invariant Slater determinants. We then demonstrate that these criteria can be explicitly verified for TBG and TBG-like systems at the chiral limit, using properties of Jacobi- $${\theta }$$ θ and Weierstrass- $${\wp }$$ ℘ functions.

arXiv:2312.15314•Communications in Mathematical Physics•May 30, 2025•6 citations

On the Hartree–Fock Ground State Manifold in Magic Angle Twisted Graphene Systems

Kevin D. Stubbs, Simon Becker, Lin Lin

Abstract. Recent experiments have shown that magic angle twisted bilayer graphene (MATBG) can exhibit correlated insulator behavior at half-filling. Seminal theoretical results towards understanding this phase in MATBG have shown that Hartree–Fock ground states (with a positive single-particle excitation gap) can be exact many-body ground states of an idealized flat band interacting (FBI) Hamiltonian. We prove that in the absence of spin and valley degrees of freedom, the only Hartree–Fock ground states of the FBI Hamiltonian for MATBG are two ferromagnetic Slater determinants. Incorporating spin and valley degrees of freedom, we provide a complete characterization of the Hartree–Fock ground state manifold, which is generated by a [Formula: see text] hidden symmetry group acting on five elements. We also introduce new tools for ruling out translation symmetry breaking in the Hartree–Fock ground state manifold, which may be of independent interest.

arXiv:2403.19890•Multiscale Modeling and Simulation•Aug 28, 2025•3 citations

Electron Correlations Study toward Understanding Twisted Bilayer Graphene: Two-Orbital Hubbard Model and Exact Numerical Calculations

Florentino Lopéz‐Urías, Francisco Sánchez-Ochoa

We investigated the honeycomb lattice using a two-orbital (α and β) Hubbard model and many-body exact calculations. The spin gap (Δ s ), charge gap (Δ c ), intrasite ( γ i i α β = ⟨ S → i α · S → i β ⟩ ), intersite ( γ i j α β = ⟨ S → i α · S → j β ⟩ ), spin–spin correlations, and orbital magnetic moments μ iα are analyzed for different band-filling ( n ), which denotes the number of electrons per site. Special attention was set in one-quarter-filled ( n = 1) and half-filled ( n = 2) bands. The results are presented by varying Coulomb repulsion ( U), Hund’s coupling ( J H ), external magnetic field ( B ), and orbital energy splitting (Δϵ). For the one-quarter-filled band, an intrasite magnetic disorder ( γ i i α β ≈ 0 ) begins to form for large values of U / t, which is interpreted as a possible signature of a resonating valence bond state. Full orbital polarization ( γ i i α β = 0.25 ) with intersite antiferromagnetic (AFM) order ( γ i j α β = − 0.25 ) was observed in the half-filled band for U / t ≫ 1, indicating that the two-orbital Hubbard model can be mapped to an AFM Heisenberg model with S = 1. The increase of J H is conductive to insulating and metallic phases for half- and one-quarter-filled bands, respectively. For a fractional n, the honeycomb lattice can exhibit spin-glass-like behavior, predominating as a disordered magnetic system characterized by random values of γ i j α β . The magnetic field induced a reduction in Δ c, which became zero at the critical external magnetic field. Twisted bilayer graphene (TBG) exhibits approximately flat electronic bands at magic angles and interesting correlated phenomena, such as unconventional superconductivity and correlated insulating states around n = 1 and n = 2, respectively. The dominance of electron–electron interactions in TBG makes the Hubbard model, in various versions (extended or multi-orbital), offer a practical, structure-scale account for investigating many-body effects in large moiré unit cells, where first-principles density-functional theory calculations are computationally challenging.

The Journal of Physical Chemistry C•Feb 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
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.
Abstract At a magic relative twist angle, magic angle twisted bilayer graphene (MATBG) has an octet of flat bands that can host strong correlation physics when partially filled. A key theoretical discovery in MATBG is the existence of ferromagnetic Slater determinants as exact ground states of the corresponding flat band interacting (FBI) Hamiltonian. The FBI Hamiltonian describes the behavior of electrons that interact with each other in a high-dimensional space, and is constructed from the band structure of the non-interacting Bistritzer–MacDonald model at the chiral limit. A key property of the FBI Hamiltonian for MATBG is that it is frustration-free and can be written as a sum of non-commuting terms. In this work, we provide a complete characterization of the ground state manifold of the FBI Hamiltonian, proving that it is precisely the linear span of such ferromagnetic Slater determinants.
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 One of the most remarkable theoretical findings in magic angle twisted bilayer graphene (TBG) is the emergence of ferromagnetic Slater determinants as exact ground states for the interacting Hamiltonian at the chiral limit. This discovery provides an explanation for the correlated insulating phase which has been experimentally observed at half filling. This work is the first mathematical study of interacting models in magic angle graphene systems. These include not only TBG but also TBG-like systems featuring four flat bands per valley, and twisted trilayer graphene systems with equal twist angles. We identify symmetries of the chiral limit of the Bistritzer-MacDonald Hamiltonian that are responsible for characterizing the Hartree-Fock ground states as zero energy many-body ground states. Furthermore, for a general class of Hamiltonian, we establish criteria that the ferromagnetic Slater determinants are the unique ground states within the class of uniformly half-filled, translation invariant Slater determinants. We then demonstrate that these criteria can be explicitly verified for TBG and TBG-like systems at the chiral limit, using properties of Jacobi- θ{\theta }θ θ and Weierstrass- ℘{\wp }℘ ℘ functions.
Abstract. Recent experiments have shown that magic angle twisted bilayer graphene (MATBG) can exhibit correlated insulator behavior at half-filling. Seminal theoretical results towards understanding this phase in MATBG have shown that Hartree–Fock ground states (with a positive single-particle excitation gap) can be exact many-body ground states of an idealized flat band interacting (FBI) Hamiltonian. We prove that in the absence of spin and valley degrees of freedom, the only Hartree–Fock ground states of the FBI Hamiltonian for MATBG are two ferromagnetic Slater determinants. Incorporating spin and valley degrees of freedom, we provide a complete characterization of the Hartree–Fock ground state manifold, which is generated by a [Formula: see text] hidden symmetry group acting on five elements. We also introduce new tools for ruling out translation symmetry breaking in the Hartree–Fock ground state manifold, which may be of independent interest.
We investigated the honeycomb lattice using a two-orbital (α and β) Hubbard model and many-body exact calculations. The spin gap (Δ s ), charge gap (Δ c ), intrasite ( γ i i α β = ⟨ S → i α · S → i β ⟩ ), intersite ( γ i j α β = ⟨ S → i α · S → j β ⟩ ), spin–spin correlations, and orbital magnetic moments μ iα are analyzed for different band-filling ( n ), which denotes the number of electrons per site. Special attention was set in one-quarter-filled ( n = 1) and half-filled ( n = 2) bands. The results are presented by varying Coulomb repulsion ( U), Hund’s coupling ( J H ), external magnetic field ( B ), and orbital energy splitting (Δϵ). For the one-quarter-filled band, an intrasite magnetic disorder ( γ i i α β ≈ 0 ) begins to form for large values of U / t, which is interpreted as a possible signature of a resonating valence bond state. Full orbital polarization ( γ i i α β = 0.25 ) with intersite antiferromagnetic (AFM) order ( γ i j α β = − 0.25 ) was observed in the half-filled band for U / t ≫ 1, indicating that the two-orbital Hubbard model can be mapped to an AFM Heisenberg model with S = 1. The increase of J H is conductive to insulating and metallic phases for half- and one-quarter-filled bands, respectively. For a fractional n, the honeycomb lattice can exhibit spin-glass-like behavior, predominating as a disordered magnetic system characterized by random values of γ i j α β . The magnetic field induced a reduction in Δ c, which became zero at the critical external magnetic field. Twisted bilayer graphene (TBG) exhibits approximately flat electronic bands at magic angles and interesting correlated phenomena, such as unconventional superconductivity and correlated insulating states around n = 1 and n = 2, respectively. The dominance of electron–electron interactions in TBG makes the Hubbard model, in various versions (extended or multi-orbital), offer a practical, structure-scale account for investigating many-body effects in large moiré unit cells, where first-principles density-functional theory calculations are computationally challenging.
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.