Atom-partitioned multipole expansions for electrostatic potential boundary conditions

详细信息    查看全文

• 作者：M. Lee^a ; ^{michael.s.lee131.civ@mail.mil" class="auth_mail" title="E-mail the corresponding author} ; K. Leiter^a ; C. Eisner^b ; ^a ; J. Knap^a
• 关键词：Finite element method ; Density functional theory ; Spherical harmonics ; Self-consistent field ; Poisson equation ; Partial differential equation ; Quantum chemistry
• 刊名：Journal of Computational Physics
• 出版年：2017
• 出版时间：1 January 2017
• 年：2017
• 卷：328
• 期：Complete
• 页码：344-353
• 全文大小：1055 K

文摘

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.