The physics of strongly correlated quantum particles within a flat band was originally explored as a route to itinerant ferromagnetism and, indeed, a celebrated theorem by Lieb rigorously establishes that the ground state of the repulsive Hubbard model on a bipartite lattice with an unequal number of sites in each sublattice must have nonzero spin $S$ at half filling. Recently, there has been interest in Lieb geometries due to the possibility of topological insulator, nematic, and Bose-Einstein condensed (BEC) phases. In this paper, we extend the understanding of the attractive Hubbard model on the Lieb lattice by using determinant quantum Monte Carlo to study real space charge and pair correlation functions not addressed by the Lieb theorems. Specifically, our results show unusual charge and charge transfer signatures within the flat band, and a reduction in pairing order at $\ensuremath{\rho}=2/3$ and $\ensuremath{\rho}=4/3$, the points at which the flat band is first occupied and then completely filled. We compare our results to the case of flat bands in the Kagome lattice and demonstrate that the behavior observed in the two cases is rather different.The physics of strongly correlated quantum particles within a flat band was originally explored as a route to itinerant ferromagnetism and, indeed, a celebrated theorem by Lieb rigorously establishes that the ground state of the repulsive Hubbard model on a bipartite lattice with an unequal number of sites in each sublattice must have nonzero spin $S$ at half filling. Recently, there has been interest in Lieb geometries due to the possibility of topological insulator, nematic, and Bose-Einstein condensed (BEC) phases. In this paper, we extend the understanding of the attractive Hubbard model on the Lieb lattice by using determinant quantum Monte Carlo to study real space charge and pair correlation functions not addressed by the Lieb theorems. Specifically, our results show unusual charge and charge transfer signatures within the flat band, and a reduction in pairing order at $\ensuremath{\rho}=2/3$ and $\ensuremath{\rho}=4/3$, the points at which the flat band is first occupied and then completely filled. We compare our results to the case of flat bands in the Kagome lattice and demonstrate that the behavior observed in the two cases is rather different.

arXiv:1404.5482Physical Review BSep 11, 2014111 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.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.16206Physical review. B./Physical review. BAug 26, 20254 citations

Quantum monte carlo simulations for high-Tc superconductors

A. Muramatsu, G. Dopf, J. Wagner +2 more

Advances in Solid State PhysicsSep 11, 20071 citation

S

at half filling. Recently, there has been interest in Lieb geometries due to the possibility of topological insulator, nematic, and Bose-Einstein condensed (BEC) phases. In this paper, we extend the understanding of the attractive Hubbard model on the Lieb lattice by using determinant quantum Monte Carlo to study real space charge and pair correlation functions not addressed by the Lieb theorems. Specifically, our results show unusual charge and charge transfer signatures within the flat band, and a reduction in pairing order at $\ensuremath{\rho}=2/3$ and $\ensuremath{\rho}=4/3$, the points at which the flat band is first occupied and then completely filled. We compare our results to the case of flat bands in the Kagome lattice and demonstrate that the behavior observed in the two cases is rather different.

{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

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.