A low-scaling method for second order Møller–Plesset calculations

A new low-order scaling algorithm for second-order Moller–Plesset perturbation theory (MP2) is described. The method employs localized orbitals. However, in contrast to other local MP2 schemes reported in the literature, all atomic orbitals with significant contribution to the correlation energy for a given pair or orbitals are retained. Our approach introduces no significant approximation if all significant pairs are taken into account. However, the method is intended for the description of correlation between spatially close (“strong”) pairs; for distant pairs, more economical alternative methods can be used. The algorithm, including a new implementation of the simultaneous two-index integral transformation scheme, is described in detail. Test calculations demonstrate asymptotic linear scaling for large systems. Strong-pair MP2 calculations with more that 1800 contracted basis functions performed on a single PC are reported. As a consequence of its lower scaling, the calculation of the MP2 energy is less expensive than the calculations of the Hartree–Fock energy for large systems. This has been demonstrated for strong pair correlation but is expected to hold also for full MP2 calculations.

[1]  Hans-Joachim Werner,et al.  Local perturbative triples correction (T) with linear cost scaling , 2000 .

[2]  Guntram Rauhut,et al.  Integral transformation with low‐order scaling for large local second‐order Møller–Plesset calculations , 1998 .

[3]  Jan Almlöf,et al.  Laplace transform techniques in Mo/ller–Plesset perturbation theory , 1992 .

[4]  Svein Saebo,et al.  Avoiding the integral storage bottleneck in LCAO calculations of electron correlation , 1989 .

[5]  Philippe Y. Ayala,et al.  Linear scaling second-order Moller–Plesset theory in the atomic orbital basis for large molecular systems , 1999 .

[6]  J. Pople,et al.  Self‐consistent molecular orbital methods. XX. A basis set for correlated wave functions , 1980 .

[7]  P. C. Hariharan,et al.  The influence of polarization functions on molecular orbital hydrogenation energies , 1973 .

[8]  Martin Head-Gordon,et al.  Noniterative local second order Mo/ller–Plesset theory: Convergence with local correlation space , 1998 .

[9]  Georg Hetzer,et al.  Multipole approximation of distant pair energies in local MP2 calculations , 1998 .

[10]  Georg Hetzer,et al.  Low-order scaling local electron correlation methods. I. Linear scaling local MP2 , 1999 .

[11]  Peter Pulay,et al.  Localizability of dynamic electron correlation , 1983 .

[12]  Philippe Y. Ayala,et al.  Linear scaling coupled cluster and perturbation theories in the atomic orbital basis , 1999 .

[13]  Hans-Joachim Werner,et al.  Local treatment of electron correlation in coupled cluster theory , 1996 .

[14]  Peter Pulay,et al.  Orbital-invariant formulation and second-order gradient evaluation in Møller-Plesset perturbation theory , 1986 .

[15]  Martin Schütz,et al.  Low-order scaling local electron correlation methods. III. Linear scaling local perturbative triples correction (T) , 2000 .

[16]  Peter Pulay,et al.  An efficient reformulation of the closed‐shell self‐consistent electron pair theory , 1984 .

[17]  Peter R. Taylor,et al.  Integral processing in beyond-Hartree-Fock calculations , 1987 .

[18]  Richard A. Friesner,et al.  Pseudospectral localized Mo/ller–Plesset methods: Theory and calculation of conformational energies , 1995 .

[19]  P Pulay,et al.  Local Treatment of Electron Correlation , 1993 .

[20]  H. Schaefer.  Methods of Electronic Structure Theory , 1977 .

[21]  J. E. Boggs,et al.  An electron-rich, three-coordinated sulfur atom : a correlated ab initio study of 1,6-dioxa-6a-thiapentalene , 1992 .

[22]  Peter Pulay,et al.  Comparison of the boys and Pipek–Mezey localizations in the local correlation approach and automatic virtual basis selection , 1993, J. Comput. Chem..

[23]  Wilfried Meyer,et al.  Theory of self‐consistent electron pairs. An iterative method for correlated many‐electron wavefunctions , 1976 .

[24]  Peter Pulay,et al.  The tautomers of uracil : a local correlation treatment , 1993 .

[25]  Martin Head-Gordon,et al.  Non-iterative local second order Møller–Plesset theory , 1998 .

[26]  Georg Hetzer,et al.  Low-order scaling local correlation methods II: Splitting the Coulomb operator in linear scaling local second-order Møller–Plesset perturbation theory , 2000 .

[27]  Peter Pulay,et al.  The local correlation treatment. II. Implementation and tests , 1988 .

[28]  Peter Pulay,et al.  Local configuration interaction: An efficient approach for larger molecules , 1985 .

[29]  P. Pulay.  The force constants of benzene: Local many‐body perturbation theory vs new experiment , 1986 .

[30]  Richard A. Friesner,et al.  Pseudospectral localized generalized Mo/ller–Plesset methods with a generalized valence bond reference wave function: Theory and calculation of conformational energies , 1997 .

[31]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[32]  Hans-Joachim Werner,et al.  Low-order scaling local electron correlation methods. IV. Linear scaling local coupled-cluster (LCCSD) , 2001 .

[33]  S. Saebo.  Strategies for electron correlation calculations on large molecular systems , 1992 .

[34]  Yixiang Cao,et al.  Correlated ab Initio Electronic Structure Calculations for Large Molecules , 1999 .

[35]  Peter Pulay,et al.  Fourth‐order Mo/ller–Plessett perturbation theory in the local correlation treatment. I. Method , 1987 .