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

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

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

Published Dec 21, 2017• 26 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.96.245127

Related papers

Magnetization, d -wave superconductivity, and non-Fermi-liquid behavior in a crossover from dispersive to flat bands

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

We explore the effect of inhomogeneity on electronic properties of the two-dimensional Hubbard model on a square lattice using dynamical mean-field theory (DMFT). The inhomogeneity is introduced via modulated lattice hopping such that in the extreme inhomogeneous limit the resulting geometry is a Lieb lattice, which exhibits a flat-band dispersion. The crossover can be observed in the uniform sublattice magnetization which is zero in the homogeneous case and increases with the inhomogeneity. Studying the spatially resolved frequency-dependent local self-energy, we find a crossover from Fermi-liquid to non-Fermi-liquid behavior happening at a moderate value of the inhomogeneity. This emergence of a non-Fermi liquid is concomitant of a quasiflat band. For finite doping the system with small inhomogeneity displays $d$-wave superconductivity coexisting with incommensurate spin-density order, inferred from the presence of oscillatory DMFT solutions. The $d$-wave superconductivity gets suppressed for moderate to large inhomogeneity for any finite doping while the incommensurate spin-density order still exists.

arXiv:1903.07466•Physical review. B./Physical review. B•Sep 18, 2019•16 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

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

Martin Bercx, Johannes Stephan Hofmann, Fakher F. Assaad +1 more

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.

arXiv:1610.03314•Physical review. B./Physical review. B•Jan 3, 2017•22 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

Stripes and spin-density waves in the doped two-dimensional Hubbard model: Ground state phase diagram

Hao Xu, Hao Shi, Ettore Vitali +2 more

We determine the spin and charge orders in the ground state of the doped two-dimensional (2D) Hubbard model in its simplest form, namely with only nearest-neighbor hopping and on-site repulsion. At half-filling, the ground state is known to be an antiferromagnetic Mott insulator. Doping Mott insulators is believed to be relevant to the superconductivity observed in cuprates. A variety of candidates have been proposed for the ground state of the doped 2D Hubbard model. A recent work employing a combination of several state-of-the-art numerical many-body methods established the stripe order as the ground state near $1/8$ doping at strong interactions. In this paper, we apply one of these methods, the cutting-edge constrained-path auxiliary field quantum Monte Carlo (AFQMC) method with self-consistently optimized gauge constraints, to systematically study the model as a function of doping and interaction strength. With careful finite size scaling based on large-scale computations, we map out the ground state phase diagram in terms of its spin and charge order. We find that modulated antiferromagnetic order persists from near half-filling to about $1/5$ doping. At lower interaction strengths or larger doping, these ordered states are best described as spin-density waves, with essentially delocalized holes and modest oscillations in charge correlations. When the charge correlations are stronger (large interaction or small doping), they are best described as stripe states, with the holes more localized near the node in the antiferromagnetic spin order. In both cases, we find that the wavelength in the charge correlations is consistent with so-called filled stripes.

arXiv:2112.02187•Physical Review Research•Mar 28, 2022•74 citations

Finite temperatures and flat bands: The Hubbard model on three-dimensional Lieb lattices

L. Oliveira-Lima, Julián Faúndez, Natanael C. Costa +1 more

We investigate some thermodynamic and magnetic properties of the Hubbard model on two three-dimensional extensions of the Lieb lattice: the perovskite Lieb lattice (PLL) and the layered Lieb lattice (LLL). Using determinant quantum Monte Carlo (DQMC) simulations alongside Hartree-Fock and cluster mean-field theory approaches, we analyze how flat-band degeneracy, connectivity, and lattice anisotropy influence the emergence of magnetic order. Our results show that both geometries support finite-temperature magnetic transitions, namely, ferromagnetic (FM) on the PLL and antiferromagnetic (AFM) on the LLL. Further, we have established that the critical temperature ${T}_{c}$, as a function of the uniform onsite coupling $U$, displays a maximum, which is larger in the AFM case than in the FM one, despite the absence of flat bands in the LLL. We also provide numerical evidence to show that flat bands in the PLL rapidly generate magnetic moments, but a small interorbital coordination suppresses the increase of ${T}_{c}$ at large interaction strength $U/t$. By contrast, the LLL benefits from higher connectivity, favoring magnetic order even in the absence of flat bands. The possibilities of anisotropic interlayer hoppings and inhomogeneous onsite interactions were separately explored. We have found that magnetism in the PLL is hardly affected by hopping anisotropy since the main driving mechanism is the preserved flat band; for the LLL, by contrast, spectral weight is removed from $d$ sites, which increases ${T}_{c}$ more significantly. At mean-field level, we have obtained that setting $U=0$ on $p$ sites and $U={U}_{d}\ensuremath{\ne}0$ on $d$ sites leads to a quantum critical point at some ${U}_{d}$; this behavior was not confirmed by our DQMC simulations.

arXiv:2505.16206•Physical review. B./Physical review. B•Aug 26, 2025•4 citations
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.
We explore the effect of inhomogeneity on electronic properties of the two-dimensional Hubbard model on a square lattice using dynamical mean-field theory (DMFT). The inhomogeneity is introduced via modulated lattice hopping such that in the extreme inhomogeneous limit the resulting geometry is a Lieb lattice, which exhibits a flat-band dispersion. The crossover can be observed in the uniform sublattice magnetization which is zero in the homogeneous case and increases with the inhomogeneity. Studying the spatially resolved frequency-dependent local self-energy, we find a crossover from Fermi-liquid to non-Fermi-liquid behavior happening at a moderate value of the inhomogeneity. This emergence of a non-Fermi liquid is concomitant of a quasiflat band. For finite doping the system with small inhomogeneity displays ddd-wave superconductivity coexisting with incommensurate spin-density order, inferred from the presence of oscillatory DMFT solutions. The ddd-wave superconductivity gets suppressed for moderate to large inhomogeneity for any finite doping while the incommensurate spin-density order still exists.
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.
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.
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 determine the spin and charge orders in the ground state of the doped two-dimensional (2D) Hubbard model in its simplest form, namely with only nearest-neighbor hopping and on-site repulsion. At half-filling, the ground state is known to be an antiferromagnetic Mott insulator. Doping Mott insulators is believed to be relevant to the superconductivity observed in cuprates. A variety of candidates have been proposed for the ground state of the doped 2D Hubbard model. A recent work employing a combination of several state-of-the-art numerical many-body methods established the stripe order as the ground state near 1/81/81/8 doping at strong interactions. In this paper, we apply one of these methods, the cutting-edge constrained-path auxiliary field quantum Monte Carlo (AFQMC) method with self-consistently optimized gauge constraints, to systematically study the model as a function of doping and interaction strength. With careful finite size scaling based on large-scale computations, we map out the ground state phase diagram in terms of its spin and charge order. We find that modulated antiferromagnetic order persists from near half-filling to about 1/51/51/5 doping. At lower interaction strengths or larger doping, these ordered states are best described as spin-density waves, with essentially delocalized holes and modest oscillations in charge correlations. When the charge correlations are stronger (large interaction or small doping), they are best described as stripe states, with the holes more localized near the node in the antiferromagnetic spin order. In both cases, we find that the wavelength in the charge correlations is consistent with so-called filled stripes.
We investigate some thermodynamic and magnetic properties of the Hubbard model on two three-dimensional extensions of the Lieb lattice: the perovskite Lieb lattice (PLL) and the layered Lieb lattice (LLL). Using determinant quantum Monte Carlo (DQMC) simulations alongside Hartree-Fock and cluster mean-field theory approaches, we analyze how flat-band degeneracy, connectivity, and lattice anisotropy influence the emergence of magnetic order. Our results show that both geometries support finite-temperature magnetic transitions, namely, ferromagnetic (FM) on the PLL and antiferromagnetic (AFM) on the LLL. Further, we have established that the critical temperature Tc{T}_{c}Tc​, as a function of the uniform onsite coupling UUU, displays a maximum, which is larger in the AFM case than in the FM one, despite the absence of flat bands in the LLL. We also provide numerical evidence to show that flat bands in the PLL rapidly generate magnetic moments, but a small interorbital coordination suppresses the increase of Tc{T}_{c}Tc​ at large interaction strength U/tU/tU/t. By contrast, the LLL benefits from higher connectivity, favoring magnetic order even in the absence of flat bands. The possibilities of anisotropic interlayer hoppings and inhomogeneous onsite interactions were separately explored. We have found that magnetism in the PLL is hardly affected by hopping anisotropy since the main driving mechanism is the preserved flat band; for the LLL, by contrast, spectral weight is removed from ddd sites, which increases Tc{T}_{c}Tc​ more significantly. At mean-field level, we have obtained that setting U=0U=0U=0 on ppp sites and $U={U}_{d}\ensuremath{\ne}0$ on ddd sites leads to a quantum critical point at some Ud{U}_{d}Ud​; this behavior was not confirmed by our DQMC simulations.