Efficient calculation of the Coulomb matrix and its expansion around k=0 within the FLAPW method

Abstract

We derive formulas for the Coulomb matrix within the full-potential linearized augmented-plane-wave (FLAPW) method. The Coulomb matrix is a central ingredient in implementations of many-body perturbation theory, such as the Hartree-Fock and GW approximations for the electronic self-energy or the random-phase approximation for the dielectric function. It is represented in the mixed product basis, which combines numerical muffin-tin functions and interstitial plane waves that are here expanded with the Rayleigh formula. The resulting algorithm is very efficient in terms of both computational cost and accuracy and is superior to an implementation with the Fourier transform of the step function. In order to allow an analytic treatment of the divergence at k=0 in reciprocal space, we expand the Coulomb matrix analytically around this point without resorting to a projection onto plane waves. We then apply a basis transformation that diagonalizes the Coulomb matrix and confines the divergence to a single eigenvalue. At the same time, response matrices like the dielectric function separate into head, wings, and body with the same mathematical properties as in a plane-wave basis. As an illustration we apply the formulas to electron-energy-loss spectra for nickel at different k vectors including k=0. The convergence of the spectra towards the result at k=0 is clearly seen. Our treatment also allows to include transitions from core states that give rise to a shallow peak at high energies and lead to good agreement with experiment.

18 pages, 3 figures

References in corpus (1)

• Elimination of the linearization error in GW calculations based on the linearized augmented-plane-wave method

Cited by in corpus (22)

• Resolution-of-identity approach to Hartree-Fock, hybrid density functionals, RPA, MP2, and \textit{GW} with numeric atom-centered orbital basis functions
• Strength of effective Coulomb interactions in graphene and graphite
• Efficient implementation of the GW approximation within the all-electron FLAPW method
• Implementation of screened hybrid functionals based on the Yukawa potential within the LAPW basis set
• Self-consistent GW: All-electron implementation with localized basis functions
• Electronic structure of Pu and Am metals by self consistent relativistic GW
• Wannier-function approach to spin excitations in solids
• Wannier Function Approach to Realistic Coulomb Interactions in Layered Materials and Heterostructures
• Hybrid functionals within the all-electron FLAPW method: implementation and applications of PBE0
• Competing Coulomb and electron-phonon interactions in NbS${}_2$
• All-electron periodic $G_0{W}_0$ implementation with numerical atomic orbital basis functions: algorithm and benchmarks
• Local exact exchange potentials within the all-electron FLAPW method and a comparison with pseudopotential results
• Precise response functions in all-electron methods: Application to the optimized-effective-potential approach
• The HSE hybrid functional within the FLAPW method and its application to GdN
• Theoretical investigation into the possibility of very large moments in Fe16N2
• Lattice polarization effects on the screened Coulomb interaction $W$ of the GW approximation
• Exciton-Driven Renormalization of Quasiparticle Band Structure in Monolayer MoS2
• Analytic evaluation of the electronic self-energy in the GW approximation for two electrons on a sphere
• First-principles calculation of electronic excitations in solids with SPEX
• Position Representation of Effective Electron-Electron Interactions in Solids
• Low-scaling \textit{GW} calculations of quasi-particle energies for extended systems within the numerical atomic orbital framework
• Effect of Local Electron-Electron Correlation in Hydrogen-like Impurities in Ge