Efficient calculation of the Coulomb matrix and its expansion around k=0 within the FLAPW method
Friedrich, Christoph
Schindlmayr, Arno
Blügel, Stefan
ddc:530
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 constructed from products of FLAPW basis functions. The interstitial plane waves 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. Without additional approximations, 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 (EELS) for nickel at different k vectors including k=0. The convergence of the spectra towards the result at k=0 is clearly seen. Our all-electron treatment also allows to include transitions from 3s and 3p core states in the EELS spectrum that give rise to a shallow peak at high energies and lead to good agreement with experiment.
Elsevier
2009
info:eu-repo/semantics/article
doc-type:article
text
http://purl.org/coar/resource_type/c_6501
https://ris.uni-paderborn.de/record/18636
Friedrich C, Schindlmayr A, Blügel S. Efficient calculation of the Coulomb matrix and its expansion around k=0 within the FLAPW method. <i>Computer Physics Communications</i>. 2009;180(3):347-359. doi:<a href="https://doi.org/10.1016/j.cpc.2008.10.009">10.1016/j.cpc.2008.10.009</a>
eng
info:eu-repo/semantics/altIdentifier/doi/10.1016/j.cpc.2008.10.009
info:eu-repo/semantics/altIdentifier/issn/0010-4655
info:eu-repo/semantics/altIdentifier/wos/000264735800002
info:eu-repo/semantics/altIdentifier/arxiv/0811.2363
info:eu-repo/semantics/closedAccess