We map out the possible ordered states in bilayer graphene at the neutrality point by extending the previous renormalization group treatment of many-body instabilities to finite temperature, trigonal warping and externally applied perpendicular electric field. We were able to analytically determine all outcomes of the RG flow equations for the nine four-fermion coupling constants. While the full phase diagram exhibits a rich structure, we confirm that when forward scattering dominates, the only ordering tendency with divergent susceptibility at finite temperature is the nematic. At finite temperature, this result is stable with respect to small back and layer imbalance scattering; further increasing their strength leads to the layer antiferromagnet. We also determine conditions for other ordered states to appear and compare our results to the special cases of attractive and repulsive Hubbard models where exact results are available.
We map out the possible ordered states in bilayer graphene at the neutrality point by extending the previous renormalization group treatment of many-body instabilities to finite temperature, trigonal warping and externally applied perpendicular electric field. We were able to analytically determine all outcomes of the RG flow equations for the nine four-fermion coupling constants. While the full phase diagram exhibits a rich structure, we confirm that when forward scattering dominates, the only ordering tendency with divergent susceptibility at finite temperature is the nematic. At finite temperature, this result is stable with respect to small back and layer imbalance scattering; further increasing their strength leads to the layer antiferromagnet. We also determine conditions for other ordered states to appear and compare our results to the special cases of attractive and repulsive Hubbard models where exact results are available.
Both insulating and conducting electronic behaviors have been experimentally seen in clean bilayer graphene samples at low temperature, and there is still no consensus on the nature of the interacting ground state at half filling and in the absence of a magnetic field. Theoretically, several possibilities for the insulating ground states have been predicted for weak interaction strength. However, a recent renormalization-group calculation on a Hubbard model for charge-neutral bilayer graphene with short-range interactions suggests the emergence of low-energy Dirac fermions that would stabilize the metallic phase for weak interactions. Using a nonperturbative projective quantum Monte Carlo, we calculate the ground state for bilayer graphene using a realistic model for the Coulomb interaction that includes both short-range and long-range contributions. We find that a finite critical on-site interaction is needed to gap bilayer graphene and the transition belongs to the Gross-Neveu universality class, thereby confirming the Hubbard model expectations even in the presence of a long-range Coulomb potential, in agreement with our theoretical analysis. In addition, we also find that the critical on-site interactions necessary to destabilize the metallic ground state decrease with increasing interlayer coupling.
We study the competition of spin- and charge-density waves and their quantum multicritical behavior for the semimetal-insulator transitions of low-dimensional Dirac fermions. Employing the effective Gross-Neveu-Yukawa theory with two order parameters as a model for graphene and a growing number of other two-dimensional Dirac materials allows us to describe the physics near the multicritical point at which the semimetallic and the spin- and charge-density-wave phases meet. With the help of a functional renormalization group approach, we are able to reveal a complex structure of fixed points, the stability properties of which decisively depend on the number of Dirac fermions Nf. We give estimates for the critical exponents and observe crucial quantitative corrections as compared to the previous first-order $\ensuremath{\epsilon}$ expansion. For small Nf, the universal behavior near the multicritical point is determined by the chiral Heisenberg universality class supplemented by a decoupled, purely bosonic, Ising sector. At large Nf, a novel fixed point with nontrivial couplings between all sectors becomes stable. At intermediate Nf, including the graphene case (Nf=2), no stable and physically admissible fixed point exists. Graphene's phase diagram in the vicinity of the intersection between the semimetal, antiferromagnetic, and staggered density phases should consequently be governed by a triple point exhibiting first-order transitions.
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.
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.
We investigate the quantum many-body instabilities for electrons on the honeycomb lattice at half filling with extended interactions, motivated by a description of graphene and related materials. We employ a recently developed fermionic functional renormalization group scheme, which allows for highly resolved calculations of wave-vector dependences in the low-energy effective interactions. We encounter the expected anti-ferromagnetic spin density wave for a dominant on-site repulsion between electrons, and charge order with different modulations for dominant pure nth nearest-neighbor repulsive interactions. Novel instabilities towards incommensurate charge density waves take place when nonlocal density interactions among several bond distances are included simultaneously. Moreover, for more realistic Coulomb potentials in graphene including enough nonlocal terms there is a suppression of charge order due to competition effects between the different charge ordering tendencies, and if the on-site term fails to dominate, the semimetallic state is rendered stable. The possibility of a topological Mott insulator being the favored tendency for dominating second-nearest-neighbor interactions is not realized in our results with high momentum resolution.
In the present thesis, the properties of the entanglement spectra of the ground state of graphene bilayers Bernal-stacked is analytically investigated. In addition, the entanglement spectra of the ground state of the honeycomb lattice in the presence of superconductivity instabilities are analytically studied. We consider not only the eigenvalues of the reduced density matrix that form the entanglement spectrum, but also, its eigenstates. From these eigenstates, we construct topological quantities, such as the Berry phase and the Chern number, in order to investigate the topological properties of the entanglement spectrum. In the first part of this thesis, we present an analytical study of the graphene bilayers, mainly focusing on the effects of the trigonal warping. When the term causing the trigonal warping is neglected, the entanglement Hamiltonian obtained by tracing out one layer shows a proportionality with the energetic Hamiltonian of the remaining monolayer graphene, in the limit of strongly coupled layers. We demonstrate that this proportionality leads to an agreement of the topological quantities of these Hamiltonians. Moreover, we demonstrate that the entanglement spectrum of graphene bilayers with a trigonal warping spectrum clearly differs in the geometric shape from the energy spectrum of the remaining monolayer graphene. However, there is an agreement of the topological quantities such as Berry phase contributions to Chern numbers. In the second part of this thesis, we give a detailed study of graphene in the presence of superconductivity instabilities, mainly considering the s-wave and the chiral dx2−y2 +idxy state.We investigate the relationship between the entanglement and energy spectrum, making use of the concepts of the Chern number constructed from eigenstates of the entanglement Hamiltonian. We demonstrate that the entanglement and remaining subsystem Hamiltonians can have different topologies. These findings are illustrated by considering the entanglement Hamiltonian of the ground state of graphene with dx2−y2 + idxy superconductivity obtained by tracing out one spin direction. Our investigations are based on closed analytical expressions for the full eigensystem in the entire Brillouin zone of bilayer graphene with a trigonally warped spectrum and graphene with superconductivity instabilities.