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

Finite-density corrections to the unitary Fermi gas: A lattice perspective from dynamical mean-field theory

A. Privitera, Massimo Capone, Claudio Castellani

Published Jan 25, 2010• 18 citations
arXivHTMLPDF

Abstract

We investigate the approach to the universal regime of the dilute unitary Fermi gas as the density is reduced to zero in a lattice model. To this end we study the chemical potential, superfluid order parameter and internal energy of the attractive Hubbard model in three different lattices with densities of states (DOSs) which share the same low-energy behavior of fermions in three-dimensional free space: a cubic lattice, a ``Bethe lattice'' with a semicircular DOS, and a ``lattice gas'' with parabolic dispersion and a sharp energy cutoff that ensures the normalization of the DOS. The model is solved using dynamical mean-field theory, that treats directly the thermodynamic limit and arbitrarily low densities, eliminating finite-size effects. At densities on the order of one fermion per site the lattice and its specific form dominate the results. The evolution to the low-density limit is smooth and it does not allow to define an unambiguous low-density regime. Such finite-density effects are significantly reduced using the lattice gas, and they are maximal for the three-dimensional cubic lattice. Even though dynamical mean-field theory is bound to reduce to the more standard static mean field in the limit of zero density due to the local nature of the self-energy and of the vertex functions, it compares well with accurate Monte Carlo simulations down to the lowest densities accessible to the latter.

DOI
10.1103/physrevb.81.014523

Related papers

Lattice approaches to dilute Fermi gases: Legacy of broken Galilean invariance

A. Privitera, Massimo Capone

In the dilute limit, the properties of fermionic lattice models with short-range attractive interactions converge to those of a dilute Fermi gas in continuum space. We investigate this connection using mean-field theory, and we show that the existence of finite lattice spacing has consequences down to very small densities. In particular, we show that the reduced translational invariance associated with the lattice periodicity has a pivotal role in the finite-density corrections to the universal zero-density limit. For a parabolic dispersion with a sharp cutoff, we provide an analytical expression for the corrections, and we find that the unavoidable cutoff contributes at leading-order to the corrections to the relevant observables. In a generic lattice we find a universal power-law behavior ${n}^{1/3}$ which leads to significant corrections already for small densities. Our results place strong constraints on lattice extrapolations of dilute Fermi gas properties.

arXiv:1112.0587•Physical Review A•Jan 30, 2012•13 citations

Zero-density-limit extrapolation of the superfluid transition temperature in a unitary atomic Fermi gas on a lattice

Qijin Chen

The superfluid transition temperature ${T}_{c}$ of a unitary Fermi gas on a three-dimensional isotropic lattice with an attractive on-site interaction is investigated as a function of density $n$, from half filling down to $5\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}7}$ per unit cell, using a pairing fluctuation theory. We show that except at very low densities (${n}^{1/3}<0.2$), where ${T}_{c}/{E}_{F}$ is linear in ${n}^{1/3}$, ${T}_{c}/{E}_{F}$ exhibits a significant higher-order nonlinear dependence on ${n}^{1/3}$. This calls for extra caution against possible significant error in the zero-density-limit value of ${T}_{c}/{E}_{F}$ from typical quantum Monte Carlo (QMC) simulations, obtained by linearly extrapolating data points at intermediate and high densities. Our result, ${T}_{c}/{E}_{F}=0.256$, at the $n=0$ limit, is close to, and should be compared with, the maximum ${T}_{c}/{E}_{F}$ near the Bose-Einstein condensation regime obtained from QMC calculations.

arXiv:1109.5327•Physical Review A•Aug 9, 2012•8 citations

Superfluid-to-insulator phase transition in a unitary Fermi gas

Nir Barnea

We study the properties of a dilute low-energy two-component Fermi gas. To this end, we approximate the Fermi gas by the Hubbard lattice Hamiltonian and solve it using the dynamical mean-field approximation. At unitarity, this approximation yields the value $\ensuremath{\xi}\ensuremath{\approx}0.44$ for the ratio between the ground-state energy per particle of the interacting and free systems $(E∕N=\ensuremath{\xi}{E}_{FG})$, in remarkable agreement with full quantum Monte Carlo simulations. Investigating the evolution of the system with temperature, we have found that below the critical temperature ${T}_{c}$, the system is a superfluid and the energy gap is decreasing monotonously. For temperatures above ${T}_{c}$, the system is an insulator and the corresponding energy gap is monotonously increasing.

arXiv:0803.2293•Physical Review A•Nov 24, 2008•8 citations

Universal quantum behavior of interacting fermions in one-dimensional traps: From few particles to the trap thermodynamic limit

Adriano Angelone, Massimo Campostrini B, Ettore Vicari

We investigate the ground-state properties of trapped fermion systems described by the Hubbard model with an external confining potential. We discuss the universal behaviors of systems in different regimes: from few particles, i.e., in dilute regimes, to the trap thermodynamic limit. The asymptotic trap-size (TS) dependence in the dilute regime (increasing the trap size $\ensuremath{\ell}$ keeping the particle number $N$ fixed) is described by a universal TS scaling controlled by the dilute fixed point associated with the metal-to-vacuum quantum transition. This scaling behavior is numerically checked by DMRG simulations of the one-dimensional (1D) Hubbard model. In particular, the particle density and its correlations show crossovers among different regimes: for strongly repulsive interactions they approach those of a spinless Fermi gas, for weak interactions those of a free Fermi gas, and for strongly attractive interactions they match those of a gas of hard-core bosonic molecules. The large-$N$ limit keeping the ratio $N/\ensuremath{\ell}$ fixed corresponds to a 1D trap thermodynamic limit. We address issues related to the accuracy of the local density approximation (LDA). We show that the particle density approaches its LDA in the large-$\ensuremath{\ell}$ limit. When the trapped system is in the metallic phase, corrections at finite $\ensuremath{\ell}$ are $O({\ensuremath{\ell}}^{\ensuremath{-}1})$ and oscillating around the center of the trap. They become significantly larger at the boundary of the fermion cloud, where they get suppressed as $O({\ensuremath{\ell}}^{\ensuremath{-}1/3})$ only. This anomalous behavior arises from the nontrivial scaling at the metal-to-vacuum transition occurring at the boundaries of the fermion cloud.

arXiv:1401.3514•Physical Review A•Feb 28, 2014•25 citations

Low-density phase diagram of the three-dimensional electron gas

Sam Azadi, N. D. Drummond

Variational and diffusion quantum Monte Carlo methods are employed to investigate the zero-temperature phase diagram of the three-dimensional homogeneous electron gas at very low density. Fermi fluid and body-centered cubic Wigner crystal ground-state energies are determined using Slater-Jastrow-backflow and Slater-Jastrow many-body wave functions at different densities and spin polarizations in finite simulation cells. Finite-size errors are removed using twist-averaged boundary conditions and extrapolation of the energy per particle to the thermodynamic limit of infinite system size. Unlike previous studies, our results show that the electron gas undergoes a first-order quantum phase transition directly from a paramagnetic fluid to a body-centered cubic crystal at density parameter ${r}_{\text{s}}=86.6(7)$, with no region of stability for an itinerant ferromagnetic fluid.

arXiv:2201.08743•Physical review. B./Physical review. B•Jun 23, 2022•31 citations

Exact solution of the Bose-Hubbard model on the Bethe lattice

Guilhem Semerjian, Marco Tarzia, Francesco Zamponi

The exact solution of a quantum Bethe lattice model in the thermodynamic limit amounts to solve a functional self-consistent equation. In this paper we obtain this equation for the Bose-Hubbard model on the Bethe lattice, under two equivalent forms. The first one, based on a coherent-state path integral, leads in the large connectivity limit to the mean-field treatment of Fisher et al. [Phys. Rev. B 40, 546 (1989)] at the leading order, and to the bosonic dynamical mean field theory as a first correction, as recently derived by Byczuk and Vollhardt [Phys. Rev. B 77, 235106 (2008)]. We obtain an alternative form of the equation using the occupation number representation, which can be easily solved with an arbitrary numerical precision, for any finite connectivity. We thus compute the transition line between the superfluid and Mott insulator phases of the model, along with thermodynamic observables and the space and imaginary-time dependence of correlation functions. The finite connectivity of the Bethe lattice induces a richer physical content with respect to its infinitely connected counterpart: a notion of distance between sites of the lattice is preserved, and the bosons are still weakly mobile in the Mott insulator phase. The Bethe lattice construction can be viewed as an approximation to the finite-dimensional version of the model. We show indeed a quantitatively reasonable agreement between our predictions and the results of Quantum Monte Carlo simulations in two and three dimensions.

arXiv:0904.3075•Physical Review B•Jul 31, 2009•42 citations
We investigate the approach to the universal regime of the dilute unitary Fermi gas as the density is reduced to zero in a lattice model. To this end we study the chemical potential, superfluid order parameter and internal energy of the attractive Hubbard model in three different lattices with densities of states (DOSs) which share the same low-energy behavior of fermions in three-dimensional free space: a cubic lattice, a Bethe lattice'' with a semicircular DOS, and a lattice gas'' with parabolic dispersion and a sharp energy cutoff that ensures the normalization of the DOS. The model is solved using dynamical mean-field theory, that treats directly the thermodynamic limit and arbitrarily low densities, eliminating finite-size effects. At densities on the order of one fermion per site the lattice and its specific form dominate the results. The evolution to the low-density limit is smooth and it does not allow to define an unambiguous low-density regime. Such finite-density effects are significantly reduced using the lattice gas, and they are maximal for the three-dimensional cubic lattice. Even though dynamical mean-field theory is bound to reduce to the more standard static mean field in the limit of zero density due to the local nature of the self-energy and of the vertex functions, it compares well with accurate Monte Carlo simulations down to the lowest densities accessible to the latter.
In the dilute limit, the properties of fermionic lattice models with short-range attractive interactions converge to those of a dilute Fermi gas in continuum space. We investigate this connection using mean-field theory, and we show that the existence of finite lattice spacing has consequences down to very small densities. In particular, we show that the reduced translational invariance associated with the lattice periodicity has a pivotal role in the finite-density corrections to the universal zero-density limit. For a parabolic dispersion with a sharp cutoff, we provide an analytical expression for the corrections, and we find that the unavoidable cutoff contributes at leading-order to the corrections to the relevant observables. In a generic lattice we find a universal power-law behavior n1/3{n}^{1/3}n1/3 which leads to significant corrections already for small densities. Our results place strong constraints on lattice extrapolations of dilute Fermi gas properties.
The superfluid transition temperature Tc{T}_{c}Tc​ of a unitary Fermi gas on a three-dimensional isotropic lattice with an attractive on-site interaction is investigated as a function of density nnn, from half filling down to $5\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}7}$ per unit cell, using a pairing fluctuation theory. We show that except at very low densities (n1/3<0.2{n}^{1/3}<0.2n1/3<0.2), where Tc/EF{T}_{c}/{E}_{F}Tc​/EF​ is linear in n1/3{n}^{1/3}n1/3, Tc/EF{T}_{c}/{E}_{F}Tc​/EF​ exhibits a significant higher-order nonlinear dependence on n1/3{n}^{1/3}n1/3. This calls for extra caution against possible significant error in the zero-density-limit value of Tc/EF{T}_{c}/{E}_{F}Tc​/EF​ from typical quantum Monte Carlo (QMC) simulations, obtained by linearly extrapolating data points at intermediate and high densities. Our result, Tc/EF=0.256{T}_{c}/{E}_{F}=0.256Tc​/EF​=0.256, at the n=0n=0n=0 limit, is close to, and should be compared with, the maximum Tc/EF{T}_{c}/{E}_{F}Tc​/EF​ near the Bose-Einstein condensation regime obtained from QMC calculations.
We study the properties of a dilute low-energy two-component Fermi gas. To this end, we approximate the Fermi gas by the Hubbard lattice Hamiltonian and solve it using the dynamical mean-field approximation. At unitarity, this approximation yields the value $\ensuremath{\xi}\ensuremath{\approx}0.44$ for the ratio between the ground-state energy per particle of the interacting and free systems $(E∕N=\ensuremath{\xi}{E}_{FG})$, in remarkable agreement with full quantum Monte Carlo simulations. Investigating the evolution of the system with temperature, we have found that below the critical temperature Tc{T}_{c}Tc​, the system is a superfluid and the energy gap is decreasing monotonously. For temperatures above Tc{T}_{c}Tc​, the system is an insulator and the corresponding energy gap is monotonously increasing.
We investigate the ground-state properties of trapped fermion systems described by the Hubbard model with an external confining potential. We discuss the universal behaviors of systems in different regimes: from few particles, i.e., in dilute regimes, to the trap thermodynamic limit. The asymptotic trap-size (TS) dependence in the dilute regime (increasing the trap size $\ensuremath{\ell}$ keeping the particle number NNN fixed) is described by a universal TS scaling controlled by the dilute fixed point associated with the metal-to-vacuum quantum transition. This scaling behavior is numerically checked by DMRG simulations of the one-dimensional (1D) Hubbard model. In particular, the particle density and its correlations show crossovers among different regimes: for strongly repulsive interactions they approach those of a spinless Fermi gas, for weak interactions those of a free Fermi gas, and for strongly attractive interactions they match those of a gas of hard-core bosonic molecules. The large-NNN limit keeping the ratio $N/\ensuremath{\ell}$ fixed corresponds to a 1D trap thermodynamic limit. We address issues related to the accuracy of the local density approximation (LDA). We show that the particle density approaches its LDA in the large-$\ensuremath{\ell}$ limit. When the trapped system is in the metallic phase, corrections at finite $\ensuremath{\ell}$ are $O({\ensuremath{\ell}}^{\ensuremath{-}1})$ and oscillating around the center of the trap. They become significantly larger at the boundary of the fermion cloud, where they get suppressed as $O({\ensuremath{\ell}}^{\ensuremath{-}1/3})$ only. This anomalous behavior arises from the nontrivial scaling at the metal-to-vacuum transition occurring at the boundaries of the fermion cloud.
Variational and diffusion quantum Monte Carlo methods are employed to investigate the zero-temperature phase diagram of the three-dimensional homogeneous electron gas at very low density. Fermi fluid and body-centered cubic Wigner crystal ground-state energies are determined using Slater-Jastrow-backflow and Slater-Jastrow many-body wave functions at different densities and spin polarizations in finite simulation cells. Finite-size errors are removed using twist-averaged boundary conditions and extrapolation of the energy per particle to the thermodynamic limit of infinite system size. Unlike previous studies, our results show that the electron gas undergoes a first-order quantum phase transition directly from a paramagnetic fluid to a body-centered cubic crystal at density parameter rs=86.6(7){r}_{\text{s}}=86.6(7)rs​=86.6(7), with no region of stability for an itinerant ferromagnetic fluid.
The exact solution of a quantum Bethe lattice model in the thermodynamic limit amounts to solve a functional self-consistent equation. In this paper we obtain this equation for the Bose-Hubbard model on the Bethe lattice, under two equivalent forms. The first one, based on a coherent-state path integral, leads in the large connectivity limit to the mean-field treatment of Fisher et al. [Phys. Rev. B 40, 546 (1989)] at the leading order, and to the bosonic dynamical mean field theory as a first correction, as recently derived by Byczuk and Vollhardt [Phys. Rev. B 77, 235106 (2008)]. We obtain an alternative form of the equation using the occupation number representation, which can be easily solved with an arbitrary numerical precision, for any finite connectivity. We thus compute the transition line between the superfluid and Mott insulator phases of the model, along with thermodynamic observables and the space and imaginary-time dependence of correlation functions. The finite connectivity of the Bethe lattice induces a richer physical content with respect to its infinitely connected counterpart: a notion of distance between sites of the lattice is preserved, and the bosons are still weakly mobile in the Mott insulator phase. The Bethe lattice construction can be viewed as an approximation to the finite-dimensional version of the model. We show indeed a quantitatively reasonable agreement between our predictions and the results of Quantum Monte Carlo simulations in two and three dimensions.