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:2107.12256 · Physical review. B./Physical review. B

Interacting spinless fermions on the square lattice: Charge order, phase separation, and superconductivity

Kou-Han Ma, Ning-Hua Tong

Published Oct 11, 2021• 8 citations
arXivHTMLPDF

Abstract

We investigate the phase diagram of spinless fermions on a square lattice with nearest-neighbor interaction, using the recently developed projective truncation approximation in Green's function equation of motion. For attractive interaction, the ground state is in a homogeneous $p+ip$ superconducting (SC) phase at high or low fermion densities. Near half filling is a phase separation (PS) between the SC phases. Allowing inhomogeneous solution, we obtain $p$-wave SC domains with positive interface energy. As temperature increases, the SC phases transit into normal phases above ${T}_{\mathrm{sc}}$, generating a homogeneous normal phase (far away from $n=1/2$), or a PS between normal phases with different densities (close to $n=1/2$). Further increasing temperature to ${T}_{\mathrm{ps}}$, the PS disappears and the particle-hole symmetry of the Hamiltonian is recovered. For repulsive interaction, depending on the filling, the ground state is in charge-ordered phase (half filling), charge-disordered phase (large hole/particle doping), or PS between them (weak doping). At finite temperature, the regime of charge order phase moves to finite $V$ and extends to finite doping regime.

DOI
10.1103/physrevb.104.155116

Related papers

Superconducting phases of the square-lattice extended Hubbard model

Wei-Chih Chen, Yao Wang, Cheng-Chien Chen

We study the square-lattice extended Hubbard model with on-site $U$ and nearest-neighbor $V$ interactions by exact diagonalization. We show that nonequilibrium quench dynamics can help determine the equilibrium phase transition boundaries, which agree with the calculations of the fidelity metric, dynamical structure factor, and correlation function. At half filling, the phase diagrams in the strong-coupling regime include spin density wave and ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$-wave superconductivity at large positive $U$, charge density wave (extended ${s}^{*}$-wave superconductivity) at large positive (negative) $V$, and $s$-wave superconductivity at large negative $U$ with vanishing $V$. The energies of different particle sectors also help determine the phase separation region. With carrier doping, charge fluctuation results in a strong competition between different orders, making it more difficult to identify the leading instability on finite-size clusters. Nevertheless, the more exotic $p$-wave superconducting pairing is found to be enhanced when the system is heavily overdoped by $37.5%--50%$ holes, especially in the interaction parameter range relevant to the cuprate superconductors.

arXiv:2206.01119•Physical review. B./Physical review. B•Aug 18, 2023•23 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

Hole on a stripe in a spinless fermion model

Uzi Hizi, Christopher L. Henley

In the spinless fermion model on a square lattice with infinite nearest-neighbor repulsion, holes doped into the half-filled ordered state form stripes which, at low doping, are stable against phase separation into an ordered state and a hole-rich metal. Here we consider transport of additional holes along these stripes. The motion of a single hole on a stripe is mapped to a one-dimensional problem, a variational wave function is constructed and the energy spectrum is calculated and compared to energies obtained by exact diagonalization.

arXiv:cond-mat/0305236•Europhysics Letters (EPL)•Jan 1, 2004•2 citations

Phase diagram of interacting spinless fermions on the honeycomb lattice: A comprehensive exact diagonalization study

Sylvain Capponi, Andreas Martin Läuchli

We investigate the phase diagram of spinless fermions with nearest- and next-nearest-neighbor density-density interactions on the honeycomb lattice at half-filling. Using exact diagonalization techniques of the full Hamiltonian and constrained subspaces, combined with a careful choice of finite-size clusters, we determine the different charge orderings that occur for large interactions. In this regime, we find a two-sublattice N\'eel-like state, a charge modulated state with a tripling of the unit cell, a zigzag phase, and a charge ordered state with a 12-site unit cell we call N\'eel domain wall crystal, as well as a region of phase separation for attractive interactions. A sizable region of the phase diagram is classically degenerate, but it remains unclear whether an order-by-disorder mechanism will lift the degeneracy. For intermediate repulsion, we find evidence for a Kekul\'e or plaquette bond-order wave phase. We also investigate the possibility of a spontaneous Chern insulator phase (dubbed topological Mott insulator), as previously put forward by several mean-field studies. Although we are unable to detect convincing evidence for this phase based on energy spectra and order parameters, we find an enhancement of current-current correlations with the expected spatial structure compared to the noninteracting situation. While for the studied $t\ensuremath{-}{V}_{1}\ensuremath{-}{V}_{2}$ model, the phase transition to the putative topological Mott insulator is preempted by the phase transitions to the various ordered states, our findings might hint at the possibility for a topological Mott insulator in an enlarged Hamiltonian parameter space, where the competing phases are suppressed.

arXiv:1505.01676•Physical Review B•Aug 26, 2015•74 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

Interacting spinless fermions with disorder: the Mott transition for infinite coordination number

Götz Silvester Uhrig, Ruud Vlaming

An analytic expression for the dynamical conductivity sigma ( omega ) of spinless fermions with local disorder and nearest-neighbour repulsion is derived on lattices with infinite coordination number Z. The model is exactly solvable in the whole parameter range assuming two possible phases: a homogeneous phase and a checkerboard charge-density wave (CDW). Away from half filling the system displays anomalous behaviour: weak particle-density fluctuations favour spontaneous symmetry breaking. First, the authors investigate the effects of this anomaly on the conductivity in the AC and DC regimes. Second, they focus on the Mott transition occurring at zero temperature at half filling. The critical exponents for the conductivity are computed for the dependence on the interaction U, the disorder gamma , the filling n, the temperature T and the frequency omega .

Journal of Physics Condensed Matter•Apr 19, 1993•10 citations
We investigate the phase diagram of spinless fermions on a square lattice with nearest-neighbor interaction, using the recently developed projective truncation approximation in Green's function equation of motion. For attractive interaction, the ground state is in a homogeneous p+ipp+ipp+ip superconducting (SC) phase at high or low fermion densities. Near half filling is a phase separation (PS) between the SC phases. Allowing inhomogeneous solution, we obtain ppp-wave SC domains with positive interface energy. As temperature increases, the SC phases transit into normal phases above Tsc{T}_{\mathrm{sc}}Tsc​, generating a homogeneous normal phase (far away from n=1/2n=1/2n=1/2), or a PS between normal phases with different densities (close to n=1/2n=1/2n=1/2). Further increasing temperature to Tps{T}_{\mathrm{ps}}Tps​, the PS disappears and the particle-hole symmetry of the Hamiltonian is recovered. For repulsive interaction, depending on the filling, the ground state is in charge-ordered phase (half filling), charge-disordered phase (large hole/particle doping), or PS between them (weak doping). At finite temperature, the regime of charge order phase moves to finite VVV and extends to finite doping regime.
We study the square-lattice extended Hubbard model with on-site UUU and nearest-neighbor VVV interactions by exact diagonalization. We show that nonequilibrium quench dynamics can help determine the equilibrium phase transition boundaries, which agree with the calculations of the fidelity metric, dynamical structure factor, and correlation function. At half filling, the phase diagrams in the strong-coupling regime include spin density wave and ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$-wave superconductivity at large positive UUU, charge density wave (extended s∗{s}^{*}s∗-wave superconductivity) at large positive (negative) VVV, and sss-wave superconductivity at large negative UUU with vanishing VVV. The energies of different particle sectors also help determine the phase separation region. With carrier doping, charge fluctuation results in a strong competition between different orders, making it more difficult to identify the leading instability on finite-size clusters. Nevertheless, the more exotic ppp-wave superconducting pairing is found to be enhanced when the system is heavily overdoped by 37.537.5%--50%37.5 holes, especially in the interaction parameter range relevant to the cuprate superconductors.
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.
In the spinless fermion model on a square lattice with infinite nearest-neighbor repulsion, holes doped into the half-filled ordered state form stripes which, at low doping, are stable against phase separation into an ordered state and a hole-rich metal. Here we consider transport of additional holes along these stripes. The motion of a single hole on a stripe is mapped to a one-dimensional problem, a variational wave function is constructed and the energy spectrum is calculated and compared to energies obtained by exact diagonalization.
We investigate the phase diagram of spinless fermions with nearest- and next-nearest-neighbor density-density interactions on the honeycomb lattice at half-filling. Using exact diagonalization techniques of the full Hamiltonian and constrained subspaces, combined with a careful choice of finite-size clusters, we determine the different charge orderings that occur for large interactions. In this regime, we find a two-sublattice N\'eel-like state, a charge modulated state with a tripling of the unit cell, a zigzag phase, and a charge ordered state with a 12-site unit cell we call N\'eel domain wall crystal, as well as a region of phase separation for attractive interactions. A sizable region of the phase diagram is classically degenerate, but it remains unclear whether an order-by-disorder mechanism will lift the degeneracy. For intermediate repulsion, we find evidence for a Kekul\'e or plaquette bond-order wave phase. We also investigate the possibility of a spontaneous Chern insulator phase (dubbed topological Mott insulator), as previously put forward by several mean-field studies. Although we are unable to detect convincing evidence for this phase based on energy spectra and order parameters, we find an enhancement of current-current correlations with the expected spatial structure compared to the noninteracting situation. While for the studied $t\ensuremath{-}{V}_{1}\ensuremath{-}{V}_{2}$ model, the phase transition to the putative topological Mott insulator is preempted by the phase transitions to the various ordered states, our findings might hint at the possibility for a topological Mott insulator in an enlarged Hamiltonian parameter space, where the competing phases are suppressed.
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.
An analytic expression for the dynamical conductivity sigma ( omega ) of spinless fermions with local disorder and nearest-neighbour repulsion is derived on lattices with infinite coordination number Z. The model is exactly solvable in the whole parameter range assuming two possible phases: a homogeneous phase and a checkerboard charge-density wave (CDW). Away from half filling the system displays anomalous behaviour: weak particle-density fluctuations favour spontaneous symmetry breaking. First, the authors investigate the effects of this anomaly on the conductivity in the AC and DC regimes. Second, they focus on the Mott transition occurring at zero temperature at half filling. The critical exponents for the conductivity are computed for the dependence on the interaction U, the disorder gamma , the filling n, the temperature T and the frequency omega .