The direct approach to gravitation and electrostatics method for periodic systems.

The direct approach to gravitation and electrostatics (DAGE) algorithm is an accurate, efficient, and flexible method for calculating electrostatic potentials. In this paper, we show that the algorithm can be easily extended to consider systems with many different kinds of periodicities, such as crystal lattices, surfaces, or wires. The accuracy and performance are nearly the same for periodic and aperiodic systems. The electrostatic potential for semiperiodic systems, namely defects in crystal lattices, can be obtained by combining periodic and aperiodic calculations. The method has been applied to an ionic model system mimicking NaCl, and to a corresponding covalent model system.

[1]  Dage Sundholm,et al.  A non-iterative numerical solver of Poisson and Helmholtz equations using high-order finite-element functions , 2005 .

[2]  Lucas Visscher,et al.  Calculation of local excitations in large systems by embedding wave-function theory in density-functional theory. , 2008, Physical chemistry chemical physics : PCCP.

[3]  Benny G. Johnson,et al.  THE CONTINUOUS FAST MULTIPOLE METHOD , 1994 .

[4]  Robert J. Harrison,et al.  Singular operators in multiwavelet bases , 2004, IBM J. Res. Dev..

[5]  Joachim Sauer,et al.  Point defects in CaF2 and CeO2 investigated by the periodic electrostatic embedded cluster method. , 2009, The Journal of chemical physics.

[6]  Hideo Sekino,et al.  Basis set limit Hartree-Fock and density functional theory response property evaluation by multiresolution multiwavelet basis. , 2008, The Journal of chemical physics.

[7]  Stefan Goedecker,et al.  Efficient solution of Poisson's equation with free boundary conditions. , 2006, The Journal of chemical physics.

[8]  Dage Sundholm,et al.  Parallel implementation of a direct method for calculating electrostatic potentials. , 2007, The Journal of chemical physics.

[9]  Leslie Greengard,et al.  A fast algorithm for particle simulations , 1987 .

[10]  Leif Laaksonen,et al.  A Numerical Hartree-Fock Program for Diatomic Molecules , 1996 .

[11]  A. Warshel,et al.  Frozen density functional approach for ab initio calculations of solvated molecules , 1993 .

[12]  Gregory Beylkin,et al.  Multiresolution quantum chemistry: basic theory and initial applications. , 2004, The Journal of chemical physics.

[13]  R. Orlando,et al.  Convergence properties of the cluster model in the study of local perturbations in ionic systems. The case of bulk defects in MgO , 1994 .

[14]  P. P. Ewald.  Die Berechnung optischer und elektrostatischer Gitterpotentiale , 1921 .

[15]  Reinhold Schneider,et al.  Daubechies wavelets as a basis set for density functional pseudopotential calculations. , 2008, The Journal of chemical physics.

[16]  Stefan Goedecker,et al.  Efficient and accurate three-dimensional Poisson solver for surface problems. , 2007, The Journal of chemical physics.

[17]  Robert J. Harrison,et al.  Multiresolution Quantum Chemistry in Multiwavelet Bases , 2003, International Conference on Computational Science.

[18]  James R. Chelikowsky,et al.  Real-space pseudopotential method for first principles calculations of general periodic and partially periodic systems , 2008 .

[19]  D. Sundholm,et al.  Universal method for computation of electrostatic potentials. , 2005, The Journal of chemical physics.

[20]  Kimihiko Hirao,et al.  A linear-scaling spectral-element method for computing electrostatic potentials. , 2008, The Journal of chemical physics.

[21]  L Greengard,et al.  Fast Algorithms for Classical Physics , 1994, Science.

[22]  P. Pyykkö,et al.  Two‐dimensional fully numerical solutions of molecular Schrödinger equations. I. One‐electron molecules , 1983 .