1. Solution to the Modified Helmholtz Equation for Arbitrary Periodic Charge Densities
- Author
-
Miriam Winkelmann, Edoardo Di Napoli, Daniel Wortmann, and Stefan Blügel
- Subjects
Helmholtz equation ,materials science ,Materials Science (miscellaneous) ,Green functions technique ,Biophysics ,FOS: Physical sciences ,General Physics and Astronomy ,02 engineering and technology ,01 natural sciences ,symbols.namesake ,electronic structure methods ,0103 physical sciences ,partial differential equations ,ddc:530 ,Physical and Theoretical Chemistry ,010306 general physics ,Fourier series ,Mathematical Physics ,density functional theory ,Condensed Matter - Materials Science ,Partial differential equation ,Mathematical analysis ,Yukawa potential ,Materials Science (cond-mat.mtrl-sci) ,Charge density ,Mathematical Physics (math-ph) ,Computational Physics (physics.comp-ph) ,021001 nanoscience & nanotechnology ,electrostatics ,lcsh:QC1-999 ,symbols ,Poisson's equation ,0210 nano-technology ,Multipole expansion ,Physics - Computational Physics ,Bessel function ,lcsh:Physics - Abstract
We present a general method for solving the modified Helmholtz equation without shape approximation for an arbitrary periodic charge distribution, whose solution is known as the Yukawa potential or the screened Coulomb potential. The method is an extension of Weinert's pseudo-charge method [M. Weinert, J. Math. Phys. 22, 2433 (1981)] for solving the Poisson equation for the same class of charge density distributions. The inherent differences between the Poisson and the modified Helmholtz equation are in their respective radial solutions. These are polynomial functions, for the Poisson equation, and modified spherical Bessel functions, for the modified Helmholtz equation. This leads to a definition of a modified pseudo-charge density and modified multipole moments. We have shown that Weinert's convergence analysis of an absolutely and uniformly convergent Fourier series of the pseudo-charge density is transferred to the modified pseudo-charge density. We conclude by illustrating the algorithmic changes necessary to turn an available implementation of the Poisson solver into a solver for the modified Helmholtz equation., submitted to the Journal of Mathematical Physics
- Published
- 2021
- Full Text
- View/download PDF