Search

Twisted double- and mono-bilayer graphene are graphene-based moir\'e materials hosting strongly correlated fermions in a gate-tunable conduction band with a topologically non-trivial character. Using unbiased exact diagonalization complemented by unrestricted Hartree-Fock calculations, we find that the strong electron-electron interactions lead to a non-coplanar magnetic state, which has the same symmetries as the tetrahedral antiferromagnet on the triangular lattice and can be thought of as a skyrmion lattice commensurate with the moir\'e scale, competing with a set of ferromagnetic, topological charge density waves featuring 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|=1$, the formation of charge density wave order which is intimately connected to a skyrmion lattice phase is consistent with recent experiments on these systems., Comment: 35 pages, 12 figures

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., Comment: Main text: 6 pages, 4 figures. Supplementary information: 10 pages, 16 figures

We perform an extensive exact diagonalization study of interaction driven insulators in spin- and valley-polarized moir\'{e} 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\'{e} heterostructure and affirms the relevance of our results beyond the scope of graphene based materials., Comment: 16 pages, 18 figures

Tang et al. [Science 361, 570 (2018)] report on the properties of Dirac fermions with both on-site and Coulomb interactions. The substantial decrease up to ~40% of the Fermi velocity of Dirac fermions with on-site interaction is inconsistent with the numerical data near the Gross-Neveu quantum critical point. This results from an inappropriate finite-size extrapolation., Comment: Submitted: October 24, 2018, Accepted: November 4, 2019

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., Comment: 9 pages, 10 figures

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_\text{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 energy 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 torus energy spectrum on finite-size clusters with periodic boundaries and extrapolate them to the thermodynamic limit. Additionally, we compute the torus energy 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 energy spectrum and we observe prominent multiplicity features related to an emergent symmetry predicted from the quantum field theory. 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., Comment: 26 pages, 14 figures

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., Comment: 6+6 pages, 12 figures

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., Comment: 9 pages, 6 figures, added data for strong Coulomb repulsion and classical Ising-limit

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+1 d 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., Comment: 7+ pages, 5 figures

Topological insulators present a bulk gap, but allow for dissipationless spin transport along the edges. These exotic states are characterized by the $Z_2$ topological invariant and are protected by time-reversal symmetry. The Kane-Mele model is one model to realize this topological class in two dimensions, also called the quantum spin Hall state. In this review, we provide a pedagogical introduction to the influence of correlation effects in the quantum spin Hall states, with special focus on the half-filled Kane-Mele-Hubbard model, solved by means of unbiased determinant quantum Monte Carlo (QMC) simulations. We explain the idea of identifying the topological insulator via $\pi$-flux insertion, the $Z_2$ invariant and the associated behavior of the zero-frequency Green's function, as well as the spin Chern number in parameter-driven topological phase transitions. The examples considered are two descendants of the Kane-Mele-Hubbard model, the generalized and dimerized Kane-Mele-Hubbard model. From the $Z_2$ index, spin Chern numbers and the Green's function behavior, one can observe that correlation effects induce shifts of the topological phase boundaries. Although the implementation of these topological quantities has been successfully employed in QMC simulations to describe the topological phase transition, we also point out their limitations as well as suggest possible future directions in using numerical methods to characterize topological properties of strongly correlated condensed matter systems., Comment: 35 pages, 16 figures, brief review for Mod. Phys. Lett B

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., Comment: 5 pages, 6 figures

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., Comment: 13 pages, 9 figures

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., Comment: 11 pages, 15 figures

We employ quantum Monte Carlo techniques to calculate the $Z_2$ 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., Comment: 7 pages, 6 figures

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., Comment: 4+ pages, 4 figures; final version

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., Comment: 4+ pages including 4 figures