Abstract.A simplified approach to treating the electron correlation energy is suggested in which only the alpha-beta component of the second order Møller-Plesset energy is evaluated, and then scaled by an empirical factor which is suggested to be 1.3. This scaled opposite spin second order energy (SOS-MP2) yields results for relative energies and derivative properties that are statistically improved over the conventional MP2 method. Furthermore, the SOS-MP2 energy can be evaluated without the 5th order computational steps associated with MP2 theory, even without exploiting any spatial locality. A 4th order algorithm is given for evaluating the opposite spin MP2 energy using auxiliary basis expansions, and a Laplace approach, and timing comparisons are given. * These authors contributed equally to this work. † To whom correspondence should be addressed. E-mail: mhg@cchem.berkeley.edu 2
1.Introduction.The most popular electronic structure method for application to systems with large numbers of electrons is density functional theory (DFT) [1,2]. However DFT methods at present completely neglect the dispersion interactions [3] that give rise to base pair stacking and other long-range correlation effects (for example the TCNE dimer dianion [4]). Novel workarounds are being explored for dispersion interactions of monomers [5][6] or ordered layers and surfaces [7,8], but do not presently apply to molecular systems. More empirical modifications of standard functionals have also been developed to improve non-bonded interactions [9,10]. Also we note that present-day DFT methods are somewhat suspect for reaction barriers. Standard functionals tend to underestimate activation energies [11], largely as a consequence of the self-interaction issue [12].The simplest electronic structure alternative to DFT that can correctly treat dispersion and hydrogen-bonding interactions is second order Møller-Plesset theory (MP2) [13]. MP2 theory is capable of quite accurately treating long-range dispersion interactions [14], as well as the dispersion, polarization and covalency effects associated with hydrogen bonding (for instance in water clusters [15]). However, MP2 has several significant drawbacks: First is relatively high computational cost, even with the best standard algorithms. Second is the need for quite large atomic orbital basis sets in order to obtain good results [16], which can further reduce the upper limit on system size. Third is the fact that poor results can be obtained for open shell systems [17], in contrast to the good behavior for closed shell molecules [18].There has been significant progress in addressing the steep cost increase of MP2calculations with molecular size in recent years. Three main types of developments can be 3 identified. First are methods that reduce the prefactor without changing the underlying scaling, such as "resolution of the identity" methods [19,20] or the pseudo-spectral approach [21], and others [22]. Second are methods that attempt to exploit "underlying locality" in the MP2problem, w...