▎ 摘　　要

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.