SCS-MP2

Spin-component-scaled second-order Møller–Plesset perturbation theory (SCS-MP2) とは、MP2 計算をもとに、電子相関の補正をより高精度な QCISD(T) の結果と合うように調節して算出した方法。2003 年に Stephen Grimme が発表した。

Grimme, S. Improved second-order Møller–Plesset perturbation theory by separate scaling of parallel- and antiparallel-spin pair correlation energies. J. Chem. Phys. 2003, 118, 9095. DOI: 10.1063/1.1569242

MP2 計算ではより高次のMøller–Plesset法と比べて電子相関の補正にずれが大きい。この補正を、逆スピン（opposite spin, OS, alpha spin–beta spin) と 同スピン (same spin, SS, alpha spin–alpha & beta spin–beta spin) の電子相関による寄与に係数（それぞれ C_OS と C_SS ）をかけることで調節している。

計算とエネルギー算出の方法

通常のMP2 計算の出力を使って SCS-MP2 エネルギーを求めることができる。Gaussian での入力と出力の例^[1]は以下の通り。

入力（MP2 で一点計算する）：

#p MP2/6-31G** geom=check guess=read scf=(maxcyc=250,tight) SCRF=(PCM, Solvent=Methanol,read) iop(5/33=2) nosymm gfinput 5d

出力（エネルギー該当部分を抜粋）：

… SCF Done: E(RHF) = -383.641996926 A.U. after 12 cycles …中略… Spin components of T(2) and E(2): alpha-alpha T2 = 0.5736104110D-01 E2= -0.1594807297D+00 alpha-beta T2 = 0.3117674447D+00 E2= -0.9286360360D+00 beta-beta T2 = 0.5736104110D-01 E2= -0.1594807297D+00 The integrals were generated 1 times. ANorm= 0.1194357370D+01 E2 = -0.1247597495D+01 EUMP2 = -0.38488959442123D+03 (D はE に直して excel などで計算する。例えば -0.3848D+03 = -0.3848 x 103)

実際、普通のMP2計算では、

EUMP2 (in hatree) = E_RHF + E2

E2 = E_MP2^OS + E_MP2^SS

= E2_ab + E2_aa + E2_bb

というように計算がなされている。SCS-MP2 では

EUMP2 (in hatree) = E_RHF + E2_SCS-MP2

E2_SCS-MP2 = C_OSE_MP2^OS + C_SSE_MP2^SS

= 6/5 E2_ab + 1/3 (E2_aa + E2_bb)

のように計算する。C_OS = 6/5,  C_SS = 1/3, の値は QCISD(T) と一致するように決められた。（C_OS = C_SS = 1 のとき、普通のMP2と同じ値になる。 ）このパラメーターを使うとエネルギーは MP2 や B3LYP より QCISD(T) からの誤差が小さくなることが論文中で報告されている。

解釈

H2 や He などの典型的な2電子相関系でMP2 は、80%から85% 程度までしか相関エネルギーを再現できていないことを踏まえて、一重項2電子系では相関エネルギーは逆スピン二重励起によるものであるため、SCS-MP2 では C_OS を 6/5 として補っている。全体としてスピン相関のエネルギーに対する寄与を2電子相関分におさめるために、同スピン励起の係数は 1/3 に減らすことによってフィッティングをしている。

その後

分子間相互作用の評価ではCCSD(T)と比較して弱点も指摘されている^[2]。系に応じて C_OS と C_SS の改良値が提案されている（C_OS = 1.15 & C_SS = 0.75 など）。

Fink, R. F. Spin-component-scaled Møller–Plesset (SCSMP) perturbation theory: A generalization of the MP approach with improved properties. J. Chem. Phys. 2010, 133, 174113. DOI: 10.1063/1.3503041

イオン性液体（Ionic Liquid）用のパラメーターも開発されている。

Rigby, J.; Izgorodina, E. I. New SCS- and SOS-MP2 Coefficients Fitted to Semi-Coulombic Systems. J. Chem. Theory Comput. 2014, 10, 3111–3122. DOI: 10.1021/ct500309x

参考文献

1. Matsuoka, A.; Sandoval, C. A.; Uchiyama, M.; Noyori, R. Naka, H. Why p-Cymene? Conformational Effect in Asymmetric Hydrogenation of Aromatic Ketones with a η⁶-Arene/Ruthenium(II) Catalyst. Chem. Asian J. 2015, 10, 112–115. DOI: 10.1002/asia.201402979
2. Bachorz, R. A.; Bischoff, F. A.; Höfener, S.; Klopper, W.; Ottiger, P.; Leist, R.; Frey, J. A.; Leutwyler, S. Phys. Chem. Chem. Phys. 2008, 10, 2758–2766. DOI: 10.1039/b718494h

[amazonjs asin=”4785335238″ locale=”JP” title=”有機化学のための 量子化学計算入門: Gaussianの基本と有効活用のヒント”] [amazonjs asin=”0198755503″ locale=”JP” title=”Computational Chemistry (Oxford Chemistry Primers)”] [amazonjs asin=”B01MXYP48B” locale=”JP” title=”Introduction to Computational Chemistry (English Edition)”]