We study the impact of Coulomb interactions at half-integer filling of the moiré Chern bands of twisted double-bilayer graphene and twisted mono-bilayer graphene, using unbiased exact diagonalization complemented by unrestricted Hartree-Fock calculations. For small intra-sublattice tunneling, $w_{AA}$, a non-coplanar magnetic state is found which has the same symmetries as the tetrahedral antiferromagnet of the triangular moir\'e lattice and can be thought of as a skyrmion lattice commensurate with the moiré scale. The antiferromagnetic order competes with a set of ferromagnetic, topological charge density waves, which are favored for larger $w_{AA}$ and are associated with an approximate emergent O(3) symmetry, `rotating' the different charge density wave states into each other. Direct comparison with exact diagonalization reveals that the ordered phases are accurately described within the unrestricted Hartree-Fock approximation. Exhibiting a finite charge gap and Chern number $C=\pm1$, the formation of charge density wave order which is intimately connected to a skyrmion lattice phase is consistent with recent experiments on these systems.

Quantum simulation using synthetic systems is a promising route to solve outstanding quantum many-body problems in regimes where other approaches, including numerical ones, fail. Many platforms are being developed towards this goal, in particular based on trapped ions, superconducting circuits, neutral atoms or molecules. All of which face two key challenges: (i) scaling up the ensemble size, whilst retaining high quality control over the parameters and (ii) certifying the outputs for these large systems. Here, we use programmable arrays of individual atoms trapped in optical tweezers, with interactions controlled by laser-excitation to Rydberg states to implement an iconic many-body problem, the antiferromagnetic 2D transverse field Ising model. We push this platform to an unprecedented regime with up to 196 atoms manipulated with high fidelity. We probe the antiferromagnetic order by dynamically tuning the parameters of the Hamiltonian. We illustrate the versatility of our platform by exploring various system sizes on two qualitatively different geometries, square and triangular arrays. We obtain good agreement with numerical calculations up to a computationally feasible size (around 100 particles). This work demonstrates that our platform can be readily used to address open questions in many-body physics.

We perform an extensive exact diagonalization study of interaction driven insulators in spin- and valley-polarized Moiré flat bands of twisted bilayer graphene aligned with its hexagonal boron nitride substrate. In addition to previously reported fractional Chern insulator phases, we provide compelling evidence for competing charge-density-wave phases at multiple fractional fillings of a realistic single-band model. A thorough analysis at different interlayer hopping parameters, motivated by experimental variability, and the role of kinetic energy at various Coulomb interaction strengths highlight the competition between these phases. The interplay of the single-particle and the interaction induced hole dispersion with the inherent Berry curvature of the Chern bands is intuitively understood to be the driving mechanism for the ground-state selection. The resulting phase diagram features remarkable agreement with experimental findings in a related Moiré heterostructure and affirms the relevance of our results beyond the scope of graphene based materials.

We examine a basic lattice model of interacting fermions that exhibits quadratic band crossing points (QBCPs) in the non-interacting limit. In particular, we consider spinless fermions on the honeycomb lattice with nearest neighbor hopping $t$ and third-nearest neighbor hopping $t''$, which exhibits fine-tuned QBCPs at the corners of the Brillouin zone for ${t'' = t/2}$. In this situation, the density of states remains finite at the Fermi level of the half-filled band and repulsive nearest-neighbor interactions $V$ lead to a charge-density-wave (CDW) instability at infinitesimally small $V$ in the random-phase approximation or mean-field theory. We examine the fragility of the QBCPs against dispersion renormalizations in the ${t\mbox{-}t''\mbox{-}V}$ model using perturbation theory, and find that the $t''$-value needed for the QBCPs increases with $V$ due to the hopping renormalization. However, the instability toward CDW formation always requires a nonzero threshold interaction strength, i.e., one cannot fine-tune $t''$ to recover the QBCPs in the interacting system. These perturbative arguments are supported by quantum Monte Carlo simulations for which we carefully compare the corresponding threshold scales at and beyond the QBCP fine-tuning point. From this analysis, we thus gain a quantitative microscopic understanding of the fragility of the QBCPs in this basic interacting fermion system.

We establish the universal torus low-energy spectra at the free Dirac fixed point and at the strongly coupled chiral Ising fixed point and their subtle crossover behaviour in the Gross-Neuveu-Yukawa field theory with n_D = 4 component Dirac spinors in D = (2+1) dimensions. These fixed points and the field theories are directly relevant for the long-wavelength physics of certain interacting Dirac systems, such as repulsive spinless fermions on the honeycomb lattice or $\pi$-flux square lattice. The torus spectrum has been shown previously to serve as a characteristic fingerprint of relativistic fixed points and is a powerful tool to discriminate quantum critical behaviour in numerical simulations. Here we use a combination of exact diagonalization and quantum Monte Carlo simulations of strongly interacting fermionic lattice models, to compute the critical energy spectrum on finite-size clusters with periodic boundaries and extrapolate them to the thermodynamic limit. Additionally, we compute the torus spectrum analytically using the perturbative expansion in $\epsilon$ = 4 - D, which is in good agreement with the numerical results, thereby validating the presence of the chiral Ising fixed point in the lattice models at hand. We show that the strong interaction between the spinor field and the scalar order-parameter field strongly influences the critical torus spectrum. Building on these results we are able to address the subtle crossover physics of the low-energy spectrum flowing from the chiral Ising fixed point to the Dirac fixed point, and analyze earlier flawed attempts to extract Fermi velocity renormalizations from the low-energy spectrum.

Tang et al. report a substantial decrease up to ~40% of the Fermi velocity of interacting Dirac fermions near the Gross-Neveu quantum critical point, based on analyzing one-particle excitation gaps on the honeycomb lattice from quantum Monte Carlo simulations. This conclusion is inconsistent with the numerical data and results from an inappropriate finite-size extrapolation.

We present quantum Monte Carlo simulations for the chiral Heisenberg Gross-Neveu-Yukawa quantum phase transition of relativistic fermions with $N$ = 4 Dirac spinor components subject to a repulsive, local four fermion interaction in 2+1$D$. Here we employ a two dimensional lattice Hamiltonian with a single, spin-degenerate Dirac cone, which exactly reproduces a linear energy-momentum relation for all finite size lattice momenta in the absence of interactions. This allows us to significantly reduce finite size corrections compared to the widely studied honeycomb and $\pi$-flux lattices. A Hubbard term dynamically generates a mass beyond a critical coupling of $U_c$ = 6.76(1) as the system acquires antiferromagnetic order and SU(2) spin rotational symmetry is spontaneously broken. At the quantum phase transition we extract a self-consistent set of critical exponents $\nu$ = 0.98(1), $\eta_{\phi}$ = 0.53(1), $\eta_{\psi}$ = 0.18(1), $\beta$ = 0.75(1). We provide evidence for the continuous degradation of the quasi-particle weight of the fermionic excitations as the critical point is approached from the semimetallic phase. Finally we study the effective "speed of light" of the low-energy relativistic description, which depends on the interaction U, but is expected to be regular across the quantum phase transition. We illustrate that the strongly coupled bosonic and fermionic excitations share a common velocity at the critical point.

We study spinless fermions with nearest-neighbor repulsive interactions (t-V model) on the two-dimensional three-band Lieb lattice. At half-filling, the free electronic band structure consists of a flat band at zero energy and a single cone with linear dispersion. The flat band is expected to be unstable upon inclusion of electronic correlations, and a natural channel is charge order. However, due to the three-orbital unit cell, commensurate charge order implies an imbalance of electron and hole densities and therefore doping away from half-filling. Our numerical results show that below a finite-temperature Ising transition a charge density wave with one electron and two holes per unit cell and its partner under particle-hole transformation are spontaneously generated. Our calculations are based on recent advances in auxiliary-field and continuous-time quantum Monte Carlo simulations that allow sign-free simulations of spinless fermions at half-filling. It is argued that particle-hole symmetry breaking provides a route to access levels of finite doping, without introducing a sign problem.

We revisit the effect of local interactions on the quadratic band touching (QBT) of Bernal stacked bilayer graphene models using renormalization group (RG) arguments and quantum Monte Carlo simulations of the Hubbard model. We present an RG argument which predicts, contrary to previous studies, that weak interactions do not flow to strong coupling even if the free dispersion has a QBT. Instead they generate a linear term in the dispersion, which causes the interactions to flow back to weak coupling. Consistent with this RG scenario, in unbiased quantum Monte Carlo simulations of the Hubbard model we find compelling evidence that antiferromagnetism turns on at a finite U/t, despite the U=0 hopping problem having a QBT. The onset of antiferromagnetism takes place at a continuous transition which is consistent with a dynamical critical exponent z=1 as expected for 2+1d Gross-Neveu criticality. We conclude that generically in models of bilayer graphene, even if the free dispersion has a QBT, small local interactions generate a Dirac phase with no symmetry breaking and that there is a finite-coupling transition out of this phase to a symmetry-broken state.

In a recent article T. Grover [Phys. Rev. Lett. 111, 130402 (2013)] introduced a simple method to compute Renyi entanglement entropies in the realm of the auxiliary field quantum Monte Carlo algorithm. Here, we further develop this approach and provide a stabilization scheme to compute higher order Renyi entropies and an extension to access the entanglement spectrum. The method is tested on systems of correlated topological insulators.

We study the quantum phases of fermions with an explicit SU(N)-symmetric, Heisenberg-like nearest-neighbor flavor exchange interaction on the honeycomb lattice at half-filling. Employing projective (zero temperature) quantum Monte Carlo simulations for even values of N, we explore the evolution from a weak-coupling semimetal into the strong-coupling, insulating regime. Furthermore, we compare our numerical results to a saddle-point approximation in the large-N limit. From the large-N regime down to the SU(6) case, the insulating state is found to be a columnar valence bond crystal, with a direct transition to the semimetal at weak, finite coupling, in agreement with the mean-field result in the large-N limit. At SU(4) however, the insulator exhibits a subtly different valence bond crystal structure, stabilized by resonating valence bond plaquettes. In the SU(2) limit, our results support a direct transition between the semimetal and an antiferromagnetic insulator.

Electronic states at the ends of a narrow armchair nanoribbon give rise to a pair of non-locally entangled spins. We propose two experiments to probe these magnetic states, based on magnetometry and tunneling spectroscopy, in which correlation effects lead to a striking, nonlinear response to external magnetic fields. On the basis of low-energy theories that we derive here, it is remarkably simple to assess these nonlinear signatures for magnetic edge states. The effective theories are especially suitable in parameter regimes where other methods such as quantum Monte-Carlo simulations are exceedingly difficult due to exponentially small energy scales. The armchair ribbon setup discussed here provides a promisingly well-controlled (both experimentally and theoretically) environment for studying the principles behind edge magnetism in graphene-based nano-structures.

We describe a method for deriving effective low-energy theories of electronic interactions at graphene edges. Our method is applicable to general edges of honeycomb lattices (zigzag, chiral, and even disordered) as long as localized low-energy states (edge states) are present. The central characteristic of the effective theories is a dramatically reduced number of degrees of freedom. As a consequence, the solution of the effective theory by exact diagonalization is feasible for reasonably large ribbon sizes. The quality of the involved approximations is critically assessed by comparing the correlation functions obtained from the effective theory with numerically exact quantum Monte-Carlo calculations. We discuss effective theories of two levels: a relatively complicated fermionic edge state theory and a further reduced Heisenberg spin model. The latter theory paves the way to an efficient description of the magnetic features in long and structurally disordered graphene edges beyond the mean-field approximation.

We investigate chiral graphene nanoribbons using projective quantum Monte Carlo simulations within the local Hubbard model description and study the effects of electron-electron interactions on the electronic and magnetic properties at the ribbon edges. Static and dynamical properties are analyzed for nanoribbons of varying width and edge chirality, and compared to a self-consistent Hartee-Fock mean-field approximation. Our results show that for chiral ribbons of sufficient width, the spin correlations exhibit exceedingly long correlation lengths, even between zigzag segments that are well separated by periodic armchair regions. Characteristic enhancements in the magnetic correlations for distinct ribbon widths and chiralities are associated with energy gaps in the tight-binding limit of such ribbons. We identify specific signatures in the local density of states and low- energy modes in the local spectral function which directly relate to enhanced electronic correlations along graphene nanoribbons and which can be accessed scanning tunneling spectroscopy.

We employ quantum Monte Carlo techniques to calculate the Z₂ topological invariant in a two-dimensional model of interacting electrons that exhibits a quantum spin Hall topological insulator phase. In particular, we consider the parity invariant for inversion-symmetric systems, which can be obtained from the bulk's imaginary-time Green's function after an appropriate continuation to zero frequency. This topological invariant is used here in order to study the trivial-band to topological-insulator transitions in an interacting system with spin-orbit coupling and an explicit bond dimerization. We discuss the accessibility and behavior of this topological invariant within quantum Monte Carlo simulations.

Using a combination of quantum Monte Carlo simulations, functional renormalization group calculations and mean-field theory, we study the Hubbard model on the Bernal-stacked honeycomb bilayer at half-filling as a model system for bilayer graphene. The free bands consisting of two Fermi points with quadratic dispersions lead to a finite density of states at the Fermi level, which triggers an antiferromagnetic instability that spontaneously breaks sublattice and spin rotational symmetry once local Coulomb repulsions are introduced. Our results reveal an inhomogeneous participation of the spin moments in the ordered ground state, with enhanced moments at the three-fold coordinated sites. Furthermore, we find the antiferromagnetic ground state to be robust with respect to enhanced interlayer couplings and extended Coulomb interactions.

We study the two-dimensional Kane-Mele-Hubbard model at half filling by means of quantum Monte Carlo simulations. We present a refined phase boundary for the quantum spin liquid. The topological insulator at finite Hubbard interaction strength is adiabatically connected to the groundstate of the Kane-Mele model. In the presence of spin-orbit coupling, magnetic order at large Hubbard U is restricted to the transverse direction. The transition from the topological band insulator to the antiferromagnetic Mott insulator is in the universality class of the three-dimensional XY model. The numerical data suggest that the spin liquid to topological insulator and spin liquid to Mott insulator transitions are both continuous.

We investigate the edge-state magnetism of graphene nanoribbons using projective quantum Monte Carlo simulations and a self-consistent mean-field approximation of the Hubbard model. The static magnetic correlations are found to be short ranged. Nevertheless, the correlation length increases with the width of the ribbon such that already for ribbons of moderate widths we observe a strong trend towards mean-field-type ferromagnetic correlations at a zigzag edge. These correlations are accompanied by a dominant low-energy peak in the local spectral function and we propose that this can be used to detect edge-state magnetism by scanning tunneling microscopy. The dynamic spin structure factor at the edge of a ribbon exhibits an approximately linearly dispersing collective magnonlike mode at low energies that decays into Stoner modes beyond the energy scale where it merges into the particle-hole continuum.

We consider the Kane-Mele model with spin-orbit coupling supplemented by a Hubbard U term. On the basis of projective auxiliary field quantum Monte Carlo simulations on lattice sizes up to 15 x 15, we map out the phase diagram. The quantum spin-liquid state found in the Hubbard model is shown to be robust against weak spin-orbit interaction, and is not adiabatically connected to the spin-Hall insulating state. Beyond a critical value of U > U_c both states are unstable toward magnetic ordering. Within the quantum spin-Hall state we study the spin, charge and single-particle dynamics of the helical Luttinger liquid by retaining the Hubbard interaction only on the edge of a ribbon. The Hubbard interaction greatly suppresses charge currents along the edge, promotes edge magnetism, but leaves the single-particle signatures of the helical liquid intact.

At sufficiently low temperatures, condensed-matter systems tend to develop order. An exception are quantum spin-liquids, where fluctuations prevent a transition to an ordered state down to the lowest temperatures. While such states are possibly realized in two-dimensional organic compounds, they have remained elusive in experimentally relevant microscopic two-dimensional models. Here, we show by means of large-scale quantum Monte Carlo simulations of correlated fermions on the honeycomb lattice, a structure realized in graphene, that a quantum spin-liquid emerges between the state described by massless Dirac fermions and an antiferromagnetically ordered Mott insulator. This unexpected quantum-disordered state is found to be a short-range resonating valence bond liquid, akin to the one proposed for high temperature superconductors. Therefore, the possibility of unconventional superconductivity through doping arises. We foresee its realization with ultra-cold atoms or with honeycomb lattices made with group IV elements.

It is shown that the semimetallic state of the two-dimensional honeycomb lattice with a point-like Fermi surface is unstable towards a canted antiferromagnetic insulator upon application of an in-plane magnetic field. This instability is already present at the mean-field level; the magnetic field shifts the up- and the down-spin cones in opposite directions thereby generating a finite density of states at the Fermi surface and a perfect nesting between the up- and the down-spin Fermi sheets. This perfect nesting triggers a canted antiferromagnetic insulating state. Our conclusions, based on mean-field arguments, are confirmed by auxiliary field projective quantum Monte Carlo methods on lattices up to 12 x 12 unit cells.

We extend the weak-coupling diagrammatic determinantal algorithm to projective schemes as well as to the inclusion of phonon degrees of freedom. The projective approach provides a very efficient algorithm to access zero temperature properties. To implement phonons, we integrate them out in favor of a retarded density-density interaction and simulate the resulting purely electronic action with the weak-coupling diagrammatic determinantal algorithm. Both extensions are tested within the dynamical mean field approximation for the Hubbard and Hubbard-Holstein models.

We investigate charge ordering in a quarter-filled ladder at finite temperature by determinantal Quantum Monte Carlo. The sign problem is moderate in a wide range of model parameters relevant for NaV2O5. The charge order parameter exhibits a crossover as a function of inverse temperature on finite systems. Above a critical nearest neighbor Coulomb repulsion, the correlation length grows exponentially with inverse temperature, indicative of the ordered phase at T=0. We find a clear single-particle gap manifesting itself in a flat n(\mu) dependence at large nearest neighbor Coulomb repulsion.