Search for: All records

In density-functional theory, the exchange–correlation (XC) energy can be defined exactly through the coupling-constant (λ) averaged XC hole n̄xc(r,r′), representing the probability depletion of finding an electron at r′ due to an electron at r. Accurate knowledge of n̄xc(r,r′) has been crucial for developing XC energy density-functional approximations and understanding their performance for molecules and materials. However, there are very few systems for which accurate XC holes have been calculated since this requires evaluating the one- and two-particle reduced density matrices for a reference wave function over a range of λ while the electron density remains fixed at the physical (λ = 1) density. Although the coupled-cluster singles and doubles (CCSD) method can yield exact results for a two-electron system in the complete basis set limit, it cannot capture the electron–electron cusp using finite basis sets. Focusing on Hooke’s atom as a two-electron model system for which certain analytic solutions are known, we examine the effect of this cusp error on the XC hole calculated using CCSD. The Lieb functional is calculated at a range of coupling constants to determine the λ-integrated XC hole. Our results indicate that, for Hooke’s atoms, the error introduced by the description of the electron–electron cusp using Gaussian basis sets at the CCSD level is negligible compared to the basis set incompleteness error. The system-, angle-, and coupling-constant-averaged XC holes are also calculated and provide a benchmark against which the Perdew–Burke–Ernzerhof and local density approximation XC hole models are assessed.

Realistic description of competing phases in complex quantum materials has proven extremely challenging. For example, much of the existing density-functional-theory-based first-principles framework fails in the cuprate superconductors. Various many-body approaches involve generic model Hamiltonians and do not account for the interplay between the spin, charge, and lattice degrees of freedom. Here, by deploying the recently constructed strongly constrained and appropriately normed (SCAN) density functional, we show how the landscape of competing stripe and magnetic phases can be addressed on a first-principles basis both in the parent insulator YBa₂Cu₃O₆and the near-optimally doped YBa₂Cu₃O₇as archetype cuprate compounds. In YBa₂Cu₃O₇, we find many stripe phases that are nearly degenerate with the ground state and may give rise to the pseudogap state from which the high-temperature superconducting state emerges. We invoke no free parameters such as the HubbardU, which has been the basis of much of the existing cuprate literature. Lattice degrees of freedom are found to be crucially important in stabilizing the various phases.