G3(MP2)-RAD theory - G3(MP2)B3 for open shell systems

G3(MP2)-RAD is a variant of G3(MP2)B3 theory optimized for open shell systems. The most important changes concern the use of URCCSD(T) instead of UQCISD(T) single point calculations and the replacement of all UMP2 by restricted open shell ROMP2 single point calculations. The G3(MP2)-RAD energy at 0 degree Kelvin E₀(G3(MP2)-RAD) is defined as:

E₀(G3(MP2)-RAD) = E[URCCSD(T)/6-31G(d)//UB3LYP/6-31G(d)]
        + DE(G3MP2large)
        + DE(HLC)
        + ZPE
        + DE(SO)

The definition of the components being:

DE(G3MP2large) = E[ROMP2(FC)/G3MP2large//UB3LYP/6-31G(d)] - E[ROMP2(FC)/6-31G(d)//UB3LYP/6-31G(d)]

DE(HLC) = -An(beta) - B(n(alpha) - n(beta))
                     A = 9.413 mHartrees; B = 3.969 mHartrees (for molecules)
                     A = 9.438 mHartrees; B = 1.888 mHartrees (for atoms)
                     n(alpha) = No. of alpha valence electrons
                     n(beta) = No. of beta valence electrons

ZPE = 0.9806 * ZPE[B3LYP/6-31G(d)]

The necessary energies can be calculated most efficiently in the following sequence:

• Optimization and frequency calculation at the UB3LYP/6-31G(d) level of theory
• URCCSD(T)/6-31G(d)//UB3LYP/6-31G(d) single point (available in MOLPRO)
• ROMP2(FC)/G3MP2large//UB3LYP/6-31G(d) single point

Comments:

• The basis sets used in G3(MP2)-RAD are identical to those used in G3(MP2)B3, and use the (5d,7f) format for polarization functions in the ROMP2(FC)/G3MP2large single point calculation. A local copy of the G3MP2large basis set can be found here.
• The higher level correction (HLC) parameters A and B have been reoptimized to give the best possible performance for the G2/97 test set.
• Spin orbit correction terms E(SO) (mainly of experimental origin) are added only for atoms
• The mean absolute deviation for the full G2/97 test set for G3(MP2)-RAD is 5.17 kJ/mol. This compares to values of 5.23 for G3(MP2)B3, 4.26 for G3, 4.14 for G3B3, 4.02 for G3X, 3.96 for G3-RAD, and 3.65 for G3X-RAD.

Literature:

• D. J. Henry, C. J. Parkinson, L. Radom,
  "An Assessment of the Performance of High-Level Theoretical Procedures in the Computation of the Heats of Formation of Small Open-Shell Molecules"
  J. Phys. Chem. A 2002, 106, 7927 - 7936.
• D. J. Henry, M. B. Sullivan, L. Radom,
  "G3-RAD and G3X-RAD: Modified Gaussian-3 (G3) and Gaussian-3X (G3X) procedures for radical thermochemistry"
  J. Chem. Phys. 2003, 118, 4849 - 4860.
• L. A. Curtiss, K. Raghavachari, P. C. Redfern, V. Rassolov, J. A. Pople,
  "Gaussian-3 (G3) theory for molecules containing first and second-row atoms"
  J. Chem. Phys. 1998, 109, 7764 - 7776.
• L. A. Curtiss, K. Raghavachari,
  "G2 Theory"
  The Encyclopedia of Computational Chemistry, P. v. R. Schleyer (editor-in-chief),
  John Wiley & Sons Ltd, Athens, USA, 1998, 2, 1104 - 1114.
• A. G. Baboul, L. A. Curtiss, P. C. Redfern, K. Raghavachari,
  "Gaussian-3 theory using density functional geometries and zero-point energies"
  J. Chem. Phys. 1999, 110, 7650 - 7657.