Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userPlane wave basis set correction methods for RPA correlation energies
78
0
0.0
(
0
)
Added by
Stefan Riemelmoser MSc
Publication date
2020
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Electronic correlation energies from the random-phase approximation converge slowly with respect to the plane wave basis set size. We study the conditions, under which a short-range local density functional can be used to account for the basis set incompleteness error. Furthermore, we propose a one-shot extrapolation scheme based on the Lindhard response function of the homogeneous electron gas. The different basis set correction methods are used to calculate equilibrium lattice constants for prototypical solids of different bonding types.
rate research
Read More
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 Hedins 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.
In order to increase the accuracy of the linearized augmented plane wave method (LAPW) we present a new approach where the plane wave basis function is augmented by two different atomic radial components constructed at two different linearization energies corresponding to two different electron bands (or energy windows). We demonstrate that this case can be reduced to the standard treatment within the LAPW paradigm where the usual basis set is enriched by the basis functions of the tight binding type, which go to zero with zero derivative at the sphere boundary. We show that the task is closely related with the problem of extended core states which is currently solved by applying the LAPW method with local orbitals (LAPW+LO). In comparison with LAPW+LO, the number of supplemented basis functions in our approach is doubled, which opens up a new channel for the extension of the LAPW and LAPW+LO basis sets. The appearance of new supplemented basis functions absent in the LAPW+LO treatment is closely related with the existence of the $dot{u}_l-$component in the canonical LAPW method. We discuss properties of additional tight binding basis functions and apply the extended basis set for computation of electron energy bands of lanthanum (face and body centered structures) and hexagonal close packed lattice of cadmium. We demonstrate that the new treatment gives lower total energies in comparison with both canonical LAPW and LAPW+LO, with the energy difference more pronounced for intermediate and poor LAPW basis sets.
The constrained electron density method of embedding a Kohn-Sham system in a substrate system (first described by P. Cortona, Phys. Rev. B {bf 44}, 8454 (1991) and T.A. Wesolowski and A. Warshel, J. Phys. Chem {bf 97}, 8050 (1993)) is applied with a plane-wave basis and both local and non-local pseudopotentials. This method divides the electron density of the system into substrate and embedded electron densities, the sum of which is the electron density of the system of interest. Coupling between the substrate and embedded systems is achieved via approximate kinetic energy functionals. Bulk aluminium is examined as a test case for which there is a strong interaction between the substrate and embedded systems. A number of approximations to the kinetic-energy functional, both semi-local and non-local, are investigated. It is found that Kohn-Sham results can be well reproduced using a non-local kinetic energy functional, with the total energy accurate to better than 0.1 eV per atom and good agreement between the electron densities.
Coupled-cluster theory with single and double excitations (CCSD) is a promising ab initio method for the electronic structure of three-dimensional metals, for which second-order perturbation theory (MP2) diverges in the thermodynamic limit. However, due to the high cost and poor convergence of CCSD with respect to basis size, applying CCSD to periodic systems often leads to large basis set errors. In a common composite method, MP2 is used to recover the missing dynamical correlation energy through a focal-point correction, but the inadequacy of MP2 for metals raises questions about this approach. Here we describe how high-energy excitations treated by MP2 can be downfolded into a low-energy active space to be treated by CCSD. Comparing how the composite and downfolding approaches perform for the uniform electron gas, we find that the latter converges more quickly with respect to the basis set size. Nonetheless, the composite approach is surprisingly accurate because it removes the problematic MP2 treatment of double excitations near the Fermi surface. Using the method to estimate the CCSD correlation energy in the combined complete basis set and thermodynamic limits, we find CCSD recovers over 90% of the exact correlation energy at $r_s=4$. We also test the composite and downfolding approaches with the random-phase approximation used in place of MP2, yielding a method that is more effective but more expensive.
Accurately describing excited states within Kohn-Sham (KS) density functional theory (DFT), particularly those which induce ionization and charge transfer, remains a great challenge. Common exchange-correlation (xc) approximations are unreliable for excited states owing, in part, to the absence of a derivative discontinuity in the xc energy ($Delta$), which relates a many-electron energy difference to the corresponding KS energy difference. We demonstrate, analytically and numerically, how the relationship between KS and many-electron energies leads to the step structures observed in the exact xc potential, in four scenarios: electron addition, molecular dissociation, excitation of a finite system, and charge transfer. We further show that steps in the potential can be obtained also with common xc approximations, as simple as the LDA, when addressed from the ensemble perspective. The article therefore highlights how capturing the relationship between KS and many-electron energies with advanced xc approximations is crucial for accurately calculating excitations, as well as the ground-state density and energy of systems which consist of distinct subsystems.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
