Local correlation methods rely on the assumption that electron correlation is nearsighted. In this work, we develop a method to alleviate this assumption. This new method is demonstrated by calculating the random phase approximation (RPA) correlation energies in several one‐dimensional model systems. In this new method, the first step is to approximately decompose the RPA correlation energy to the nearsighted and farsighted components based on the wavelength decomposition of electron correlation developed by Langreth and Perdew. The short‐wavelength (SW) component of the RPA correlation energy is then considered to be nearsighted, and the long‐wavelength (LW) component of the RPA correlation energy is considered to be farsighted. The SW RPA correlation energy is calculated using a recently developed local correlation method: the embedded cluster density approximation (ECDA). The LW RPA correlation energy is calculated globally based on the system's Kohn‐Sham orbitals. This new method is termed
We review the theory and application of adiabatic exchange–correlation (xc)-kernels for ab initio calculations of ground state energies and quasiparticle excitations within the frameworks of the adiabatic connection fluctuation dissipation theorem and Hedin’s equations, respectively. Various different xc-kernels, which are all rooted in the homogeneous electron gas, are introduced but hereafter we focus on the specific class of renormalized adiabatic kernels, in particular the rALDA and rAPBE. The kernels drastically improve the description of short-range correlations as compared to the random phase approximation (RPA), resulting in significantly better correlation energies. This effect greatly reduces the reliance on error cancellations, which is essential in RPA, and systematically improves covalent bond energies while preserving the good performance of the RPA for dispersive interactions. For quasiparticle energies, the xc-kernels account for vertex corrections that are missing in the GW self-energy. In this context, we show that the short-range correlations mainly correct the absolute band positions while the band gap is less affected in agreement with the known good performance of GW for the latter. The renormalized xc-kernels offer a rigorous extension of the RPA and GW methods with clear improvements in terms of accuracy at little extra computational cost.
more » « less- NSF-PAR ID:
- 10153968
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- npj Computational Materials
- Volume:
- 5
- Issue:
- 1
- ISSN:
- 2057-3960
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract λ ‐ECDA, whereλ indicates the wavelength decomposition. The performance ofλ ‐ECDA is examined on a one‐dimensional model system: aH24 chain, in which the RPA correlation energy is highly nonlocal. In this model system, a softened Coulomb interaction is used to describe the electron‐electron and electron‐ion interactions, and slightly stronger nuclear charges (1.2 ) are assigned to the pseudo‐H atoms. Bond stretching energies, RPA correlation potentials, and Kohn‐Sham eigenvalues predicted bye λ ‐ECDA are in good agreement with the benchmarks when the clusters are made reasonably large. We find that the LW RPA correlation energy is critical for obtaining accurate prediction of the RPA correlation potential, even though the LW RPA correlation energy contributes to only a few percent of the total RPA correlation energy. -
more » « less
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.
-
de Jong, Bert ; Nieminen, Risto (Ed.)A kernel polynomial method is developed to calculate the random phase approximation (RPA) correlation energy. In the method, the RPA correlation energy is formulated in terms of the matrix that is the product of the Coulomb potential and the density linear response functions. The integration over the matrix's eigenvalues is calculated by expanding the density of states of the matrix in terms of the Chebyshev polynomials. The coefficients in the expansion are obtained through stochastic sampling. Since it is often the energy difference between two systems that is of much interest in practice, another focus of this work is to develop a correlated sampling scheme to accelerate the convergence of the stochastic calculations of the RPA correlation energy difference between two similar systems. The scheme is termed the atom-based correlated sampling (ACS). The performance of ACS is examined by calculating the isomerization energy of acetone to 2-propenol and the energy of the water–gas shift reaction. Using ACS, the convergences of these two examples are accelerated by 3.6 and 4.5 times, respectively. The methods developed in this work are expected to be useful for calculating RPA-level reaction energies for the reactions that take place in local regions, such as calculating the adsorption energies of molecules on transition metal surfaces for modeling surface catalysis.more » « less
-
Blaschke, David (Ed.)Two-point fermionic propagators in strongly-correlated media are considered with an emphasis on the dynamical interaction kernels of their equations of motion (EOM). With the many-body Hamiltonian confined by a two-body interaction, the EOMs for the two-point fermionic propagators acquire the Dyson form and, before taking any approximation, the interaction kernels decompose into the static and dynamical (time-dependent) contributions. The latter translate to the energy-dependent and the former map to the energy-independent terms in the energy domain. We dwell particularly on the energy-dependent terms, which generate long-range correlations while making feedback on their short-range static counterparts. The origin, forms, and various approximations for the dynamical kernels of one-fermion and two-fermion propagators, most relevant in the intermediate-coupling regime, are discussed. Applications to the electromagnetic dipole response of 68,70-Ni and low-energy quadrupole response of 114,116,124-Sn are presented.more » « less
-
more » « less
Understanding, predicting, and ultimately controlling exciton band structure and exciton dynamics are central to diverse chemical and materials problems. Here, we have developed a first-principles method to determine exciton dispersion and exciton–phonon interaction in semiconducting and insulating solids based on time-dependent density functional theory. The first-principles method is formulated in planewave bases and pseudopotentials and can be used to compute exciton band structures, exciton charge density, ionic forces, the non-adiabatic coupling matrix between excitonic states, and the exciton–phonon coupling matrix. Based on the spinor formulation, the method enables self-consistent noncollinear calculations to capture spin-orbital coupling. Hybrid exchange-correlation functionals are incorporated to deal with long-range electron–hole interactions in solids. A sub-Hilbert space approximation is introduced to reduce the computational cost without loss of accuracy. For validations, we have applied the method to compute the exciton band structure and exciton–phonon coupling strength in transition metal dichalcogenide monolayers; both agree very well with the previous GW-Bethe–Salpeter equation and experimental results. This development paves the way for accurate determinations of exciton dynamics in a wide range of solid-state materials.
