Efficient calculation of the Coulomb matrix and its expansion around k=0 within the FLAPW method
From MaRDI portal
Publication:603368
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.
Recommendations
- Efficient evaluation of the matrix elements of the Coulomb potential between plane waves and Gaussians
- An algorithm for the calculation of the partial wave expansion of the Coulomb-distorted plane wave
- Fast Coulomb matrix construction via compressing the interactions between continuous charge distributions
- Efficient computation of the coupling matrix in time-dependent density functional theory
- The Jacobi-matrix method in parabolic coordinates: Expansion of Coulomb functions in parabolic Sturmians
- High-performance generation of the Hamiltonian and overlap matrices in FLAPW methods
- Efficient evaluation of Coulomb integrals: the nonlinearD- and -transformations
- Efficient calculation of Coulomb matrix elements for bilayers of confined charge carriers with arbitrary spatial separation
- Improving the convergence and estimating the accuracy of summation approximants of 1/D expansions for Coulombic systems
- Analytic spin-orbit coupling matrix element formulae in FLAPW calculations
Cites work
Cited in
(5)- High-performance generation of the Hamiltonian and overlap matrices in FLAPW methods
- Efficient evaluation of the matrix elements of the Coulomb potential between plane waves and Gaussians
- Elimination of the linearization error and improved basis-set convergence within the FLAPW method
- \texttt{OptaDOS}: a tool for obtaining density of states, core-level and optical spectra from electronic structure codes
- Analytic spin-orbit coupling matrix element formulae in FLAPW calculations
This page was built for publication: Efficient calculation of the Coulomb matrix and its expansion around k=0 within the FLAPW method
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q603368)
