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

Spontaneous particle-hole symmetry breaking of correlated fermions on the Lieb lattice

Martin Bercx, Johannes Stephan Hofmann, Fakher F. Assaad, Thomas C. Lang

Published Jan 3, 2017• 22 citations
arXivHTMLPDF

Abstract

We study spinless fermions with nearest-neighbor repulsive interactions ($t\text{\ensuremath{-}}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.

DOI
10.1103/physrevb.95.035108

Related papers

Temperature and doping induced instabilities of the repulsive Hubbard model on the Lieb lattice

Pramod Kumar, Tuomas I. Vanhala, Päivi Törmä

The properties of a phase at finite interactions can be significantly influenced by the underlying dispersion of the noninteracting Hamiltonian. We demonstrate this by studying the repulsive Hubbard model on the two-dimensional Lieb lattice, which has a flat band for vanishing interaction $U$. We perform real-space dynamical mean-field theory calculations at different temperatures and dopings using a continuous-time quantum Monte Carlo impurity solver. Studying the frequency dependence of the self-energy, we find that a nonmagnetic metallic region at finite temperature displays non-Fermi-liquid behavior, which is a concomitant of the flat-band singularity. At half-filling, we also find a magnetically ordered region, where the order parameter varies linearly with the interaction strength, and a strongly correlated Mott insulating phase. The double occupancy decreases sharply for small $U$, highlighting the flat-band contribution. Away from half-filling, we observe the stripe order, i.e., an inhomogeneous spin and charge density wave of finite wavelength, which turns into a sublattice ordering at higher temperatures.

arXiv:1709.09992•Physical review. B./Physical review. B•Dec 21, 2017•26 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

Ising phase transitions and thermodynamics of correlated fermions in a two-dimensional spin-dependent lattice potential

Zhuotao Xie, Yufeng Song, Yuan-Yao He

We present a numerically exact study of the Hubbard model with spin-dependent anisotropic hopping on the square lattice using auxiliary-field quantum Monte Carlo method. At half filling, the system undergoes Ising phase transitions on cooling, leading to the formation of Ising-type antiferromagnetic order for repulsive interactions and charge-density wave order for attractive interactions at finite temperatures. By elegantly implementing the sign-problem-free condition and Hubbard-Stratonovich transformations, we achieve significant improvements in precision control of the numerical calculations and obtain highly accurate results of the transition temperatures from weak to strong interactions across representative anisotropies. We further characterize the system by examining the temperature dependence of various thermodynamic properties, including the energy, double occupancy, specific heat, and charge susceptibility. Specifically, we provide unbiased numerical results of the entropy map on temperature-interaction plane, the critical entropy, and the spin, singlon and doublon correlations, all of which are directly measurable in optical lattice experiments. Away from half filling, we explore the behavior of the sign problem and investigate the possible emergence of stripe spin-density wave order in the system with repulsive interaction.

arXiv:2412.20843•Physical review. B./Physical review. B•Mar 4, 2025•6 citations

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

Kou-Han Ma, Ning-Hua Tong

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.

arXiv:2107.12256•Physical review. B./Physical review. B•Oct 11, 2021•8 citations

Metal-insulator transition in the disordered Hubbard model of the Lieb lattice

Yueqi Li, Lingyu Tian, Tianxing Ma +1 more

Using the determinant quantum Monte Carlo method, we investigate the metal-insulator transition in the interacting disordered Hubbard model of a Lieb lattice in which the system characterizes the flat band centered at the Fermi level. By choosing suitable electron densities, we ensure the weak interaction sign problem to improve the reliability of our results. It is found that disorder and on-site Coulomb repulsive interaction produce interesting effects that induce the metal-insulator transition which is impossible in the half-filled case. The density of states at the Fermi energy is still finite in the thermodynamic limit, suggesting that the system is an Anderson insulator rather than a Mott insulator. Moreover, this doping system is paramagnetic, unlike the half-filled system which is ferrimagnetic.

arXiv:2206.03890•Physical review. B./Physical review. B•Nov 28, 2022•7 citations

Interplay of flat electronic bands with Holstein phonons

Chunhan Feng, Richard T. Scalettar

Existing quantum Monte Carlo studies have investigated the properties of fermions on a Lieb (${\mathrm{CuO}}_{2}$) lattice interacting with on-site or near-neighbor electron-electron coupling. Attention has focused on the interplay of such interactions with the macroscopic degeneracy of local zero-energy modes, from which Bloch states can be formed to produce a flat band in which energy is independent of momentum. The resulting high density of states, in combination with the Stoner criterion, suggests that there should be pronounced instabilities to ordered phases. Indeed, a theorem by Lieb rigorously establishes the existence of ferrimagnetic order. Here we study the charge density wave phases induced by electron-phonon coupling on the Lieb lattice, as opposed to previous work on electron-electron interactions. Our key result is the demonstration of charge density wave phases at one third and two thirds fillings, characterized by long-range density-density correlations between doubly occupied sites on the minority or majority sublattices, and an accompanying gap. We also compute the transition temperature to the ordered phase as a function of the electron-phonon coupling.

arXiv:2009.05595•Physical review. B./Physical review. B•Dec 23, 2020•33 citations
We study spinless fermions with nearest-neighbor repulsive interactions (t\ensuremath-Vt\text{\ensuremath{-}}Vt\ensuremath-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.
The properties of a phase at finite interactions can be significantly influenced by the underlying dispersion of the noninteracting Hamiltonian. We demonstrate this by studying the repulsive Hubbard model on the two-dimensional Lieb lattice, which has a flat band for vanishing interaction UUU. We perform real-space dynamical mean-field theory calculations at different temperatures and dopings using a continuous-time quantum Monte Carlo impurity solver. Studying the frequency dependence of the self-energy, we find that a nonmagnetic metallic region at finite temperature displays non-Fermi-liquid behavior, which is a concomitant of the flat-band singularity. At half-filling, we also find a magnetically ordered region, where the order parameter varies linearly with the interaction strength, and a strongly correlated Mott insulating phase. The double occupancy decreases sharply for small UUU, highlighting the flat-band contribution. Away from half-filling, we observe the stripe order, i.e., an inhomogeneous spin and charge density wave of finite wavelength, which turns into a sublattice ordering at higher temperatures.
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 present a numerically exact study of the Hubbard model with spin-dependent anisotropic hopping on the square lattice using auxiliary-field quantum Monte Carlo method. At half filling, the system undergoes Ising phase transitions on cooling, leading to the formation of Ising-type antiferromagnetic order for repulsive interactions and charge-density wave order for attractive interactions at finite temperatures. By elegantly implementing the sign-problem-free condition and Hubbard-Stratonovich transformations, we achieve significant improvements in precision control of the numerical calculations and obtain highly accurate results of the transition temperatures from weak to strong interactions across representative anisotropies. We further characterize the system by examining the temperature dependence of various thermodynamic properties, including the energy, double occupancy, specific heat, and charge susceptibility. Specifically, we provide unbiased numerical results of the entropy map on temperature-interaction plane, the critical entropy, and the spin, singlon and doublon correlations, all of which are directly measurable in optical lattice experiments. Away from half filling, we explore the behavior of the sign problem and investigate the possible emergence of stripe spin-density wave order in the system with repulsive interaction.
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.
Using the determinant quantum Monte Carlo method, we investigate the metal-insulator transition in the interacting disordered Hubbard model of a Lieb lattice in which the system characterizes the flat band centered at the Fermi level. By choosing suitable electron densities, we ensure the weak interaction sign problem to improve the reliability of our results. It is found that disorder and on-site Coulomb repulsive interaction produce interesting effects that induce the metal-insulator transition which is impossible in the half-filled case. The density of states at the Fermi energy is still finite in the thermodynamic limit, suggesting that the system is an Anderson insulator rather than a Mott insulator. Moreover, this doping system is paramagnetic, unlike the half-filled system which is ferrimagnetic.
Existing quantum Monte Carlo studies have investigated the properties of fermions on a Lieb (CuO2{\mathrm{CuO}}_{2}CuO2​) lattice interacting with on-site or near-neighbor electron-electron coupling. Attention has focused on the interplay of such interactions with the macroscopic degeneracy of local zero-energy modes, from which Bloch states can be formed to produce a flat band in which energy is independent of momentum. The resulting high density of states, in combination with the Stoner criterion, suggests that there should be pronounced instabilities to ordered phases. Indeed, a theorem by Lieb rigorously establishes the existence of ferrimagnetic order. Here we study the charge density wave phases induced by electron-phonon coupling on the Lieb lattice, as opposed to previous work on electron-electron interactions. Our key result is the demonstration of charge density wave phases at one third and two thirds fillings, characterized by long-range density-density correlations between doubly occupied sites on the minority or majority sublattices, and an accompanying gap. We also compute the transition temperature to the ordered phase as a function of the electron-phonon coupling.