文摘
A procedure is presented to fit gridded molecular properties to auxiliary basis sets (ABSs) of Hermite Gaussians,analogous to the density fitting (DF) method (Dunlap; et al. J. Chem. Phys. 1979, 71, 4993). In this procedure,the ab initio calculated properties (density, electrostatic potential, and/or electric field) are fitted via a linear-or nonlinear-least-squares procedure to auxiliary basis sets (ABS). The calculated fitting coefficients fromthe numerical grids are shown to be more robust than analytic density fitting due to the neglect of the corecontributions. The fitting coefficients are tested by calculating intermolecular Coulomb and exchangeinteractions for a set of dimers. It is shown that the numerical instabilities observed in DF are caused by theattempt of the ABS to fit the core contributions. In addition, this new approach allows us to reduce thenumber of functions required to obtain an accurate fit. This results in decreased computational cost, which isshown by calculating the Coulomb energy of a 4096 water box in periodic boundary conditions. Using atomcentered Hermite Gaussians, this calculation is only 1 order of magnitude slower than conventional atom-centered point charges.