6.6.5 Spin-Biased MP2 Methods (SCS-MP2, SOS-MP2, and MOS-MP2)

The accuracy of MP2 calculations can be significantly improved by semi-empirically scaling the opposite-spin (OS) and same-spin (SS) correlation components with separate scaling factors, as shown by Grimme. ⁴²¹ Grimme S.
J. Chem. Phys.
(2003), 118, pp. 9095. Link Scaling with 1.2 and 0.33 (or OS and SS components) defines the SCS-MP2 method, but other parameterizations are desirable for systems involving intermolecular interactions, as in the SCS-MI-MP2 method, which uses 0.40 and 1.29 (for OS and SS components). ²⁸⁶ DiStasio, Jr. R. A., Head-Gordon M.
Mol. Phys.
(2007), 105, pp. 1073. Link

Results of similar quality for thermochemistry can be obtained by only retaining and scaling the opposite spin correlation (by 1.3), as was recently demonstrated. ⁵⁷² Jung Y. et al.
J. Chem. Phys.
(2004), 121, pp. 9793. Link Furthermore, the SOS-MP2 energy can be evaluated using the RI approximation together with a Laplace transform technique, in effort that scales only with the 4th power of molecular size. Efficient algorithms for the energy ⁵⁷² Jung Y. et al.
J. Chem. Phys.
(2004), 121, pp. 9793. Link and the analytical gradient ⁷⁵² Lochan R. C., Shao Y., Head-Gordon M.
J. Chem. Theory Comput.
(2007), 3, pp. 988. Link of this method are available since Q-Chem v. 3.0, and offer advantages in speed over MP2 for larger molecules, as well as statistically significant improvements in accuracy.

However, we note that the SOS-MP2 method does systematically underestimate long-range dispersion (for which the appropriate scaling factor is 2 rather than 1.3) but this can be accounted for by making the scaling factor distance-dependent, which is done in the modified opposite spin variant (MOS-MP2) that has recently been proposed and tested. ⁷⁵⁰ Lochan R. C., Jung Y., Head-Gordon M.
J. Phys. Chem. A
(2005), 109, pp. 7598. Link The MOS-MP2 energy and analytical gradient are also available in Q-Chem 3.0 at a cost that is essentially identical with SOS-MP2. Timings show that the 4th-order implementation of SOS-MP2 and MOS-MP2 yields substantial speedups over RI-MP2 for molecules in the 40 heavy atom regime and larger. It is also possible to customize the scale factors for particular applications, such as weak interactions, if required.

A fourth order scaling SOS-MP2/MOS-MP2 energy calculation can be invoked by setting the CORRELATION keyword to either SOSMP2 or MOSMP2. MOS-MP2 further requires the specification of the $rem variable OMEGA, which tunes the level of attenuation of the MOS operator: ⁷⁵⁰ Lochan R. C., Jung Y., Head-Gordon M.
J. Phys. Chem. A
(2005), 109, pp. 7598. Link

The recommended OMEGA value is $\omega =0.6$ bohr${}^{-1}$. ⁷⁵⁰ Lochan R. C., Jung Y., Head-Gordon M.
J. Phys. Chem. A
(2005), 109, pp. 7598. Link The fast algorithm makes use of auxiliary basis expansions and therefore, the keyword AUX_BASIS should be set consistently with the user’s choice of BASIS. Fourth-order scaling analytical gradient for both SOS-MP2 and MOS-MP2 are also available and is automatically invoked when JOBTYPE is set to OPT or FORCE. The minimum memory requirement is 3$X^2$, where $X$ = the number of auxiliary basis functions, for both energy and analytical gradient evaluations. Disk space requirement for closed shell calculations is $\sim 2OVX$ for energy evaluation and $\sim 4OVX$ for analytical gradient evaluation.

Summary of key $rem variables to be specified: