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:1008.1901 · Physical Review B

Interacting fermions on the honeycomb bilayer: From weak to strong coupling

Oskar Vafek

Published Nov 3, 2010• 129 citations
arXivHTMLPDF

Abstract

Many-body instabilities of the half-filled honeycomb bilayer are studied using weak-coupling renormalization group (RG) as well as strong-coupling expansion. For spinless fermions and assuming parabolic degeneracy, there are four independent four-fermion contact couplings. While the dominant instability depends on the microscopic values of the couplings, the broken symmetry state is typically a gapped insulator with either broken inversion symmetry or broken time-reversal symmetry, with a quantized anomalous Hall effect. Under certain conditions, the dominant instability may appear in the particle-particle (pairing) channel. For some nongeneric fine-tuned initial conditions, weak-coupling RG trajectories flow into the noninteracting fixed point, although generally we find runaway flows which we associate with ordering tendencies. Additionally, a tight-binding model with nearest-neighbor hopping and nearest-neighbor repulsion is studied in weak and strong couplings and in each regime a gapped phase with inversion symmetry breaking is found. In the strong-coupling limit, the ground-state wave function is constructed for vanishing in-plane hopping but finite interplane hopping, which explicitly displays the broken inversion symmetry and a finite difference between the number of particles on the two layers. Finally, we discuss the spin-1/2 case and use Fierz identities to show that the number of independent four-fermion contact couplings is 9. The corresponding RG equations in the spin-1/2 case are also presented, and used to show that, just as in strong coupling, the most dominant weak-coupling instability of the repulsive Hubbard model (at half filling) is an antiferromagnet.

DOI
10.1103/physrevb.82.205106

Related papers

Electronic instabilities of the extended Hubbard model on the honeycomb lattice from functional renormalization

Yanick Volpez, Daniel D. Scherer, Michael M. Scherer

Interacting fermions on the half-filled honeycomb lattice with short-range repulsions have been suggested to host a variety of interesting many-body ground states, e.g., a topological Mott insulator. A number of recent studies of the spinless case in terms of exact diagonalization, the infinite density matrix renormalization group, and the functional renormalization group, however, indicate a suppression of the topological Mott insulating phase in the whole range of interaction parameters. Here, we complement the previous studies by investigating the quantum many-body instabilities of the physically relevant case of spin-1/2 fermions with onsite, nearest-neighbor, and second-nearest-neighbor repulsion. To this end, we employ the multipatch functional renormalization group for correlated fermions with refined momentum resolution observing the emergence of an antiferromagnetic spin-density wave and a charge-density wave for dominating onsite and nearest-neighbor repulsions, respectively. For dominating second-nearest neighbor interaction our results favor an ordering tendency towards a charge-modulated ground state over the topological Mott insulating state. The latter evades a stabilization as the leading instability by the additional onsite interaction.

arXiv:1606.01125•Physical review. B./Physical review. B•Oct 5, 2016•18 citations

Bond-ordered states and f -wave pairing of spinless fermions on the honeycomb lattice

Stephan Heßelmann, Daniel D. Scherer, Michael M. Scherer +1 more

Spinless fermions on the honeycomb lattice with repulsive nearest-neighbor interactions are known to harbour a quantum critical point at half-filling, with critical behavior in the Gross-Neveu (chiral Ising) universality class. The critical interaction strength separates a weak-coupling semimetallic regime from a commensurate charge-density-wave phase. The phase diagram of this basic model of correlated fermions on the honeycomb lattice beyond half-filling is, however, less well established. Here, we perform an analysis of its many-body instabilities using the functional renormalization group method with a basic Fermi surface patching scheme, which allows us to treat instabilities in competing channels on equal footing also away from half-filling. Between half-filling and the Van Hove filling, the free Fermi surface is holelike and we again find a charge-density wave instability to be dominant at large interactions. Moreover, its characteristics are those of the half-filled case. Directly at the Van Hove filling, the nesting property of the free Fermi surface stabilizes a dimerized bond-order phase. At lower filling, the free Fermi surface becomes electronlike and a superconducting instability with $f$-wave symmetry is found to emerge from the interplay of intra-unit-cell repulsion and collective fluctuations in the proximity to the charge-density wave instability. We estimate the extent of the various phases and extract the corresponding order parameters from the effective low-energy Hamiltonians.

arXiv:1805.02490•Physical review. B./Physical review. B•Jul 31, 2018•10 citations

Competing electronic instabilities of extended Hubbard models on the honeycomb lattice: A functional renormalization group calculation with high-wave-vector resolution

David Sánchez de la Peña, Julian Lichtenstein, Carsten Honerkamp A

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 ${n}^{\text{th}}$ 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.

arXiv:1606.01124•Physical review. B./Physical review. B•Feb 27, 2017•37 citations

Phase diagram of the half-filled extended Hubbard model in two dimensions

Biplab Chattopadhyay, D. M. Gaitonde

We consider an extended Hubbard model of interacting fermions on a lattice. The fermion kinetic energy corresponds to a tight-binding Hamiltonian with nearest-neighbor (-t) and next-nearest-neighbor (${\mathrm{t}}^{|\mathrm{I}\mathrm{H}}$) hopping matrix elements. In addition to the on-site Hubbard interaction (U) we also consider a nearest-neighbor repulsion (V). We obtain the zero-temperature phase diagram of our model within the Hartree-Fock approximation. We consider ground states having charge- and spin-density-wave ordering as well as states with orbital antiferromagnetism or spin nematic order. The latter two states correspond to particle-hole binding with ${\mathrm{d}}_{{\mathrm{x}}^{2}\mathrm{\ensuremath{-}}{\mathrm{y}}^{2}}$ symmetry in the charge and spin channels, respectively. For ${\mathrm{t}}^{|\mathrm{I}\mathrm{H}}$=0, only the charge-density-wave and spin-density-wave states are energetically stable. For nonzero ${\mathrm{t}}^{|\mathrm{I}\mathrm{H}}$, we find that orbital antiferromagnetism (or spin nematic) order is stable over a finite portion of the phase diagram at weak coupling. This region of stability is seen to grow with increasing values of ${\mathrm{t}}^{|\mathrm{I}\mathrm{H}}$.

arXiv:cond-mat/9703213•Physical review. B, Condensed matter•Jun 15, 1997•24 citations

Ground-state phase diagram of the half-filled bilayer Hubbard model

Michael Golor, Timo Reckling, Laura Classen +2 more

Employing a combination of functional renormalization group calculations and projective determinantal quantum Monte Carlo simulations, we examine the Hubbard model on the square lattice bilayer at half filling. From this combined analysis, we obtain a comprehensive account of the ground-state phase diagram with respect to the extent of the system's metallic and (antiferromagnetically ordered) Mott-insulating as well as band-insulating regions. By means of an unbiased functional renormalization group approach, we exhibit the antiferromagnetic Mott-insulating state as the relevant instability of the free metallic state, induced by any weak finite on-site repulsion. Upon performing a careful analysis of the quantum Monte Carlo data, we resolve the difficulty of identifying this antiferromagnetic ground state for finite interlayer hopping in the weak-coupling regime, where nonmonotonic finite-size corrections are shown to relate to the two-sheeted Fermi surface structure of the metallic phase. On the other hand, quantum Monte Carlo simulations are well suited to identify the transition between the Mott-insulating phase and the band insulator in the intermediate-to-strong coupling regime. Here, we compare our numerical findings to indications for the transition region obtained from the functional renormalization group procedure.

arXiv:1409.1103•Physical Review B•Nov 18, 2014•23 citations

Topological phases of the Kitaev-Hubbard model at half filling

Jean Paul Latyr Faye, David Sénéchal, Syed R. Hassan

The Kitaev-Hubbard model of interacting fermions is defined on the honeycomb lattice and, at strong coupling, interpolates between the Heisenberg model and the Kitaev model. It is basically a Hubbard model with ordinary hopping $t$ and spin-dependent hopping ${t}^{\ensuremath{'}}$. We study this model in the weak to intermediate coupling regime, at half filling, using the cellular dynamical impurity approximation (CDIA), an approach related to dynamical mean field theory but based on Potthoff's variational principle. We identify four phases in the $(U,{t}^{\ensuremath{'}})$ plane: two semimetallic phases with different numbers of Dirac points, an antiferromagnetic insulator, and an algebraic spin liquid. The last two are separated by a first-order transition. These four phases all meet at a single point and could be realized in cold atom systems.

arXiv:1506.00501•Physical Review B•Mar 27, 2014•15 citations
Many-body instabilities of the half-filled honeycomb bilayer are studied using weak-coupling renormalization group (RG) as well as strong-coupling expansion. For spinless fermions and assuming parabolic degeneracy, there are four independent four-fermion contact couplings. While the dominant instability depends on the microscopic values of the couplings, the broken symmetry state is typically a gapped insulator with either broken inversion symmetry or broken time-reversal symmetry, with a quantized anomalous Hall effect. Under certain conditions, the dominant instability may appear in the particle-particle (pairing) channel. For some nongeneric fine-tuned initial conditions, weak-coupling RG trajectories flow into the noninteracting fixed point, although generally we find runaway flows which we associate with ordering tendencies. Additionally, a tight-binding model with nearest-neighbor hopping and nearest-neighbor repulsion is studied in weak and strong couplings and in each regime a gapped phase with inversion symmetry breaking is found. In the strong-coupling limit, the ground-state wave function is constructed for vanishing in-plane hopping but finite interplane hopping, which explicitly displays the broken inversion symmetry and a finite difference between the number of particles on the two layers. Finally, we discuss the spin-1/2 case and use Fierz identities to show that the number of independent four-fermion contact couplings is 9. The corresponding RG equations in the spin-1/2 case are also presented, and used to show that, just as in strong coupling, the most dominant weak-coupling instability of the repulsive Hubbard model (at half filling) is an antiferromagnet.
Interacting fermions on the half-filled honeycomb lattice with short-range repulsions have been suggested to host a variety of interesting many-body ground states, e.g., a topological Mott insulator. A number of recent studies of the spinless case in terms of exact diagonalization, the infinite density matrix renormalization group, and the functional renormalization group, however, indicate a suppression of the topological Mott insulating phase in the whole range of interaction parameters. Here, we complement the previous studies by investigating the quantum many-body instabilities of the physically relevant case of spin-1/2 fermions with onsite, nearest-neighbor, and second-nearest-neighbor repulsion. To this end, we employ the multipatch functional renormalization group for correlated fermions with refined momentum resolution observing the emergence of an antiferromagnetic spin-density wave and a charge-density wave for dominating onsite and nearest-neighbor repulsions, respectively. For dominating second-nearest neighbor interaction our results favor an ordering tendency towards a charge-modulated ground state over the topological Mott insulating state. The latter evades a stabilization as the leading instability by the additional onsite interaction.
Spinless fermions on the honeycomb lattice with repulsive nearest-neighbor interactions are known to harbour a quantum critical point at half-filling, with critical behavior in the Gross-Neveu (chiral Ising) universality class. The critical interaction strength separates a weak-coupling semimetallic regime from a commensurate charge-density-wave phase. The phase diagram of this basic model of correlated fermions on the honeycomb lattice beyond half-filling is, however, less well established. Here, we perform an analysis of its many-body instabilities using the functional renormalization group method with a basic Fermi surface patching scheme, which allows us to treat instabilities in competing channels on equal footing also away from half-filling. Between half-filling and the Van Hove filling, the free Fermi surface is holelike and we again find a charge-density wave instability to be dominant at large interactions. Moreover, its characteristics are those of the half-filled case. Directly at the Van Hove filling, the nesting property of the free Fermi surface stabilizes a dimerized bond-order phase. At lower filling, the free Fermi surface becomes electronlike and a superconducting instability with fff-wave symmetry is found to emerge from the interplay of intra-unit-cell repulsion and collective fluctuations in the proximity to the charge-density wave instability. We estimate the extent of the various phases and extract the corresponding order parameters from the effective low-energy Hamiltonians.
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{n}^{\text{th}}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.
We consider an extended Hubbard model of interacting fermions on a lattice. The fermion kinetic energy corresponds to a tight-binding Hamiltonian with nearest-neighbor (-t) and next-nearest-neighbor (t∣IH{\mathrm{t}}^{|\mathrm{I}\mathrm{H}}t∣IH) hopping matrix elements. In addition to the on-site Hubbard interaction (U) we also consider a nearest-neighbor repulsion (V). We obtain the zero-temperature phase diagram of our model within the Hartree-Fock approximation. We consider ground states having charge- and spin-density-wave ordering as well as states with orbital antiferromagnetism or spin nematic order. The latter two states correspond to particle-hole binding with ${\mathrm{d}}_{{\mathrm{x}}^{2}\mathrm{\ensuremath{-}}{\mathrm{y}}^{2}}$ symmetry in the charge and spin channels, respectively. For t∣IH{\mathrm{t}}^{|\mathrm{I}\mathrm{H}}t∣IH=0, only the charge-density-wave and spin-density-wave states are energetically stable. For nonzero t∣IH{\mathrm{t}}^{|\mathrm{I}\mathrm{H}}t∣IH, we find that orbital antiferromagnetism (or spin nematic) order is stable over a finite portion of the phase diagram at weak coupling. This region of stability is seen to grow with increasing values of t∣IH{\mathrm{t}}^{|\mathrm{I}\mathrm{H}}t∣IH.
Employing a combination of functional renormalization group calculations and projective determinantal quantum Monte Carlo simulations, we examine the Hubbard model on the square lattice bilayer at half filling. From this combined analysis, we obtain a comprehensive account of the ground-state phase diagram with respect to the extent of the system's metallic and (antiferromagnetically ordered) Mott-insulating as well as band-insulating regions. By means of an unbiased functional renormalization group approach, we exhibit the antiferromagnetic Mott-insulating state as the relevant instability of the free metallic state, induced by any weak finite on-site repulsion. Upon performing a careful analysis of the quantum Monte Carlo data, we resolve the difficulty of identifying this antiferromagnetic ground state for finite interlayer hopping in the weak-coupling regime, where nonmonotonic finite-size corrections are shown to relate to the two-sheeted Fermi surface structure of the metallic phase. On the other hand, quantum Monte Carlo simulations are well suited to identify the transition between the Mott-insulating phase and the band insulator in the intermediate-to-strong coupling regime. Here, we compare our numerical findings to indications for the transition region obtained from the functional renormalization group procedure.
The Kitaev-Hubbard model of interacting fermions is defined on the honeycomb lattice and, at strong coupling, interpolates between the Heisenberg model and the Kitaev model. It is basically a Hubbard model with ordinary hopping ttt and spin-dependent hopping ${t}^{\ensuremath{'}}$. We study this model in the weak to intermediate coupling regime, at half filling, using the cellular dynamical impurity approximation (CDIA), an approach related to dynamical mean field theory but based on Potthoff's variational principle. We identify four phases in the $(U,{t}^{\ensuremath{'}})$ plane: two semimetallic phases with different numbers of Dirac points, an antiferromagnetic insulator, and an algebraic spin liquid. The last two are separated by a first-order transition. These four phases all meet at a single point and could be realized in cold atom systems.