Atom-partitioned multipole expansions for electrostatic potential boundary conditions
- Simulation Sciences Branch, U.S. Army Research Laboratory, Aberdeen Proving Ground, MD 21005 (United States)
- Secure Mission Solutions, a Parsons Company (United States)
Applications such as grid-based real-space density functional theory (DFT) use the Poisson equation to compute electrostatics. However, the expected long tail of the electrostatic potential requires either the use of a large and costly outer domain or Dirichlet boundary conditions estimated via multipole expansion. We find that the oft-used single-center spherical multipole expansion is only appropriate for isotropic mesh domains such as spheres and cubes. In this work, we introduce a method suitable for high aspect ratio meshes whereby the charge density is partitioned into atomic domains and multipoles are computed for each domain. While this approach is moderately more expensive than a single-center expansion, it is numerically stable and still a small fraction of the overall cost of a DFT calculation. The net result is that when high aspect ratio systems are being studied, form-fitted meshes can now be used in lieu of cubic meshes to gain computational speedup.
- OSTI ID:
- 22622232
- Journal Information:
- Journal of Computational Physics, Journal Name: Journal of Computational Physics Vol. 328; ISSN JCTPAH; ISSN 0021-9991
- Country of Publication:
- United States
- Language:
- English
Similar Records
Multipole Expansion of Coulomb Energy
Polarizable atomic multipole X-ray refinement: application to peptide crystals