Efficient implementation of the fast multipole method

Robust Tensor Hypercontraction of the Particle-Particle Ladder Term in Equation-of-Motion Coupled Cluster Theory

One method of representing a high-rank tensor as a (hyper-)product of lower-rank tensors is the tensor hypercontraction (THC) method of Hohenstein et al. This strategy has been found to be useful for reducing the polynomial scaling of coupled-cluster methods by representation of a four-dimensional tensor of electron-repulsion integrals in terms of five two-dimensional matrices. Pierce et al. have already shown that the application of a robust form of THC to the particle-particle ladder (PPL) term reduces the cost of this term in couple-cluster singles and doubles (CCSD) from O ( N 6 ) to O ( N 5 ) with negligible errors in energy with respect to the density-fitted variant. In this work, we have implemented the least-squares variant of THC (LS-THC) which does not require a nonlinear tensor factorization, including the robust form (R-LS-THC), for the calculation of the excitation and electron attachment energies using equation-of-motion coupled cluster methods EOMEE-CCSD and EOMEA-CCSD, respectively. We have benchmarked the effect of the R-LS-THC-PPL approximation on excitation energies using the comprehensive QUEST database and the accuracy of electron attachment energies using the NAB22 database. We find that errors on the order of 1 meV are achievable with a reduction in total calculation time of approximately 5 ×.

Robust Approximation of Tensor Networks: Application to Grid-Free Tensor Factorization of the Coulomb Interaction

Approximation of a tensor network by approximating (e.g., factorizing) one or more of its constituent tensors can be improved by canceling the leading-order error due to the constituents' approximation. The utility of such robust approximation is demonstrated for robust canonical polyadic (CP) approximation of a (density-fitting) factorized two-particle Coulomb interaction tensor. The resulting algebraic (grid-free) approximation for the Coulomb tensor, closely related to the factorization appearing in pseudospectral and tensor hypercontraction approaches, is efficient and accurate, with significantly reduced rank compared to the naive (nonrobust) approximation. Application of the robust approximation to the particle-particle ladder term in the coupled-cluster singles and doubles reduces the size complexity from O (N⁶) to O (N⁵) with robustness ensuring negligible errors in chemically relevant energy differences using CP ranks approximately equal to the size of the density-fitting basis.

The grid-based fast multipole method – a massively parallel numerical scheme for calculating two-electron interaction energies

A grid-based fast multipole method has been developed for calculating two-electron interaction energies for non-overlapping charge densities.

Intermediate electrostatic field for the elongation method

A simple way to improve the accuracy of the fragmentation methods is proposed. The formalism was applied to the elongation (ELG) method at restricted open-shell Hartree-Fock (ROHF) level of theory. The α-helix conformer of polyglycine was taken as a model system. The modified ELG method includes a simplified electrostatic field resulting from point-charge distribution of the system's environment. In this way the long-distance polarization is approximately taken into account. The field attenuates during the ELG process to eventually disappear when the final structure is reached. The point-charge distributions for each ELG step are obtained from charge sensitivity analysis (CSA) in force-field atoms resolution. The presence of the intermediate field improves the accuracy of ELG calculations. The errors in total energy and its kinetic and potential contributions are reduced by at least one-order of magnitude. In addition the SCF convergence of ROHF scheme is improved.

Further analysis and comparative study of intermolecular interactions using dimers from the S22 database

Accurate MP2 and CCSD(T) complete basis set (CBS) interaction energy curves (14 points for each curve) have been obtained for 20 of the dimers reported in the S22 set and analytical Morse curves have been fitted that can be used in developing updated density functional theory (DFT) and force field models. The magnitude and the effect of the basis set superposition error (BSSE) were carefully investigated. We found that going up to aug-cc-pVDZ and aug-cc-pVTZ basis sets is enough to obtain accurate CBS MP2 energies when BSSE corrected values are used but aug-cc-pVTZ and aug-cc-pVQZ basis sets are needed when the BSSE uncorrected total energies are used in CBS extrapolations. MP2 interaction energies with smaller basis sets such as 6-31G* are found to have very little dispersion energy and that the true source of dispersion attributed attractive interactions is almost entirely due to BSSE. MP2 and CCSD(T) CBS interaction energies are found to be very close to one another if aromatic systems are not involved. Comparative analyses have been performed with semiempirical and ab initio methods utilizing the moderate in size but affordable 6-31G* basis set both of which can be readily applied to macromolecular systems. The new M06-2X and M06-L DFT functionals were found to be more accurate than all methods tested herein. Interaction energy curves using the SG1 grid showed discontinuities for several of the dimer systems but this problem disappeared when finer DFT numerical grids were used.

Extension of adaptive tree code and fast multipole methods to high angular momentum particle charge densities

The development and implementation of a tree code (TC) and fast multipole method (FMM) for the efficient, linear-scaling calculation of long-range electrostatic interactions of particle distributions with variable shape and multipole character are described. The target application of these methods are stochastic boundary molecular simulations with polarizable force fields and/or combined quantum mechanical/molecular mechanical potentials. Linear-scaling is accomplished through the adaptive decomposition of the system into a hierarchy of interacting particle sets. Two methods for effecting this decomposition are evaluated: fluc-splitting and box-splitting, for which the latter is demonstrated to be generally more accurate. In addition, a generalized termination criterion is developed that delivers optimal performance at fixed error tolerance that, in the case of quadrupole-represented Drude water, effects a speed-up by a factor of 2-3 relative to a multipole-independent termination criteria. The FMM is shown to be approximately 2-3 times faster than the TC, independent of the system size and multipole order of the particles. The TC and FMM are tested for a variety of static and polarizable water systems, and for the the 70S ribosome functional complex containing an assembly of transfer and messenger RNAs.

A Linear Scaling Study of Solvent−Solute Interaction Energy of Drug Molecules in Aqua Solution

Solvent-solute interaction energies for three well-known drug molecules in water solution are computed at the Hartree-Fock and B3LYP density functional theory levels using a linear scaling technique, which allows one to explicitly include in the model water molecules up to 14 A away from the solute molecule. The dependence of calculated interaction energies on the amount of included solvent has been examined. It is found that it is necessary to account for water molecules within an 8 A radius around the drug molecule to reach the saturated solvent interaction level. Effects of electron correlation and basis set on solvent-solute interaction energies are discussed.