Collapse of the electron gas from three to two dimensions in Kohn-Sham density functional theory
Aaron D. Kaplan, Kamal Wagle, and John P. Perdew
Phys. Rev. B 98, 085147 (2018) - Published 28 August, 2018
Under pressure, a quasi-two-dimensional electron gas can collapse toward the true two-dimensional (2D) limit. In this limit, the exact exchange-correlation energy per electron has a known finite limit, but general-purpose semilocal approximate density functionals, such as the local density approximation (LDA) and the Perdew-Burke-Ernzerhof generalized gradient approximation (PBE GGA), are known to diverge to minus infinity. Here we consider a model density for a noninteracting electron gas confined to a thickness by infinite-barrier walls, with a fixed 2D density and Bohr. We estimate that LDA, PBE, and the strongly constrained and appropriately normed (SCAN) meta-GGA are accurate for the exchange-correlation energy over a wide quasi-2D range, , but not for smaller . Of these functionals, only SCAN tends to a finite limit when tends to 0. Since the noninteracting kinetic energy, treated exactly in Kohn-Sham theory, dominates in this limit within a deformable jellium model, all of the general-purpose functionals can estimate the pressure required to achieve any thickness (with SCAN and LDA better than PBE). This pressure vanishes around , where the 3D electron density is roughly that of the valence electrons in metallic potassium, and it reaches about 20 GPa at and 400 GPa at .
