Tighter multipole-based integral estimates and parallel implementation of linear-scaling AO–MP2 theory

Citations Over TimeTop 10% of 2008 papers

Abstract

Within an atomic-orbital-based (AO-based) formulation of second-order Møller-Plesset perturbation theory (MP2), we present a novel screening procedure which allows us to preselect numerically significant two-electron integrals more efficiently, especially for large basis sets. The screening is based on our recently introduced multipole-based integral estimates (MBIE) method [J. Chem. Phys., 2005, 123, 184102], that allows to exploit the 1/R(4) or 1/R(6) coupling between electronic charge distributions in transformed integral products within AO-MP2. In this way, linear scaling is attained with fully-controlled numerical accuracy. Furthermore, a parallel implementation of our linear-scaling AO-MP2 method is described, which also allows us to perform calculations with larger basis sets. First calculations reveal that for e.g. linear alkanes the scaling of the number of required transformed integral products is almost equal for 6-31G* and cc-pVTZ basis sets. Using the improved MBIE screening, the largest parallel calculation was performed for a ribozyme fragment consisting of 497 atoms and 5697 basis functions, while our largest AO-MP2 calculation was performed for a stacked DNA system (16 base pairs) comprising 1052 atoms and 10 674 basis functions on a single processor.

Related Papers

• → Distributed Multipole Analysis: Stability for Large Basis Sets(2005)618 cited
• → Rotating around the quartic angular momentum barrier in fast multipole method calculations(1996)133 cited
• → On the scaling of multipole methods for particle–particle interactions(1998)8 cited
• Multilevel fast multipole method for higher order basis functions implemented in WIPL-D Pro(2009)
• → Domain decomposition method for scattering from multiple bodies of revolution(2015)2 cited