Solving lattice density functionals close to the Mott regime
arXiv:1403.5080 · doi:10.1103/PhysRevB.89.205130
Abstract
We study a lattice version of the local density approximation (LDA) based on Behte ansatz (BALDA). Contrary to what happens in Density Functional Theory (DFT) in the continuum and despite its name, BALDA displays some very non-local features and it has a discontinuous functional derivative. The same features prevent the convergence of the self-consistent Kohn-Sham cycle thus hindering the study of BALDA solutions close to a Mott phase or in the Coulomb blockade regime. Here we propose a numerical approach which, differently from previous works, does not introduce ad hoc parameters to smear out the singularity. Our results are relevant for all lattice models where BALDA is applied ranging from Kondo systems to harmonically trapped Hubbard fermions. As an example we apply the method to the study of a one-dimensional lattice model with Hubbard interaction and a staggered potential which can be driven from an ionic to a Mott insulating state. In the Mott regime the presence of a "vacuum" allows us to calculate the different contribution to the gap and to highlight an ultranonlocality of BALDA.
References in corpus (10)
- Many-Body Physics with Ultracold Gases
- Ultracold atomic gases in optical lattices: mimicking condensed matter physics and beyond
- Localization and delocalization errors in density functional theory and implications for band-gap prediction
- Thermodynamics of the 3D Hubbard model on approach to the Neel transition
- Discerning Incompressible and Compressible Phases of Cold Atoms in Optical Lattices
- Luther-Emery Phase and Atomic-Density Waves in a Trapped Fermion Gas
- Theory of antibound states in partially filled narrow band systems
- Successes and failures of Bethe Ansatz Density Functional Theory
- Interacting fermions in 1D disordered lattices: Exploring localization and transport properties with lattice density-functional theories
- Quantitative determination of the Hubbard model phase diagram from optical lattice experiments by two-parameter scaling