メラー=プレセット法(メラー=プレセットほう、英: Møller–Plesset method)とは、計算化学の分野におけるハートリー=フォック法(平均場近似)を超える、より高精度な第一原理計算手法の一つである。レイリー=シュレーディンガーの摂動論を用いて電子相関を考慮する方法であり、摂動のレベルによりMP2、MP3、MP4などと呼ばれる。
主要な着想は1934年にクリスチャン・メラーとミルトン・S・プレセットによって発表された[1]。
レイリー=シュレーディンガーの摂動論で導かれた式は一般的な形であり、ハミルトニアンに具体的な形を与える必要がある。
メラー=プレセット法の摂動法としての振舞いはあまり良くなく、真の解に向けて単調に収束するわけではない。すなわち、取り込む摂動の次数を上げても真の解から遠ざかることがありうる。このため、メラー=プレセット法により物性値を正しく予想するためには注意を要する。
最初の定式化
MP-エネルギー補正は「シフトした」フォック演算子として定義された非摂動ハミルトニアン
と「補正ポテンシャル」として定義された摂動ハミルトニアン
を使ってレイリー=シュレーディンガー(RS)摂動論から得ることができる。
上式において、規格化されたスレイター行列式Φ0はフォック演算子の最低固有状態である。
ここで、Nは対象としている分子中の電子の数(エネルギー中の係数2は、それぞれの軌道が逆向きスピンを持つ電子対によって占有されている事実に起因する)、は通常の電子ハミルトニアン、は1電子フォック演算子、εiは二重に占有された空間軌道φiに属する軌道エネルギーである。
スレイター行列式Φ0はの固有状態であるため、以下が示される。
すなわち、0次エネルギーはΦ0に関するの期待値、ハートリー–フォックエネルギーである。同様に、「この定式化での」MP1エネルギーは
となる。したがって、意味のある最初の補正はMP2エネルギーに現われる。
閉殻分子に対するMP2式を得るために、2次のRS-PT式は2電子励起スレイター行列式を基にして書かれる(ブリルアンの定理のため、1電子励起スレイター行列式は寄与しない)。N-電子行列要素を単純化するためにスレイター–コンドン則を適用した後、MP2エネルギー(スピンを除いてブラ-ケット記法でまとめる)は
となる。上式において、𝜑iおよび𝜑jは正準被占軌道、𝜑aおよび𝜑bは仮想(空)軌道である。量εi、εj、εa、およびεbは対応する軌道エネルギーである。全電子エネルギーはハートリー–フォックエネルギー足す2次MP補正によって与えられる。
- E ≈ EHF + EMP2
別の定式化
等価な式はハミルトニアンのわずかに異なる分割によって得られる。0次と1次の寄与でのエネルギー項の分割が元々の定式化とは異なるものの、2次おび高次エネルギー補正は2つの分割法で同一の結果を与える。この定式化は(現在MP法の大口ユーザーである)化学者によって一般的に使われている。
この違いはハートリー–フォック理論でよく知られている以下の事実(ハートリー–フォックエネルギーは被占軌道エネルギーの和ではない)が原因である。
この定式化では、非摂動ハミルトニアンとしてフォック演算子を取る。
疑いもなく、この分割では
となる。
明白に、この定式化では、メラー=プレセット理論はEMP1 ≠ 0を満たさない。0次のMP方程式の解は軌道エネルギーの和であり、0次の値に1次の補正を足すことでハートリー–フォックエネルギーが得られる。元々の定式化と同じく、最初の意味のある補正は2次のエネルギーである。繰り返し述べるが、2次および高次の補正は両方の定式化において同じである。
まとめると、この定式化で得られる0次、1次、2次のエネルギーの物理的な意味は以下の通り。
- 0次:フォック演算子の固有値の和、すなわち、各分子軌道エネルギーの和に相当する。これは電子間反発を二重に計算するため、実用性は皆無である。
- 1次:電子間反発のエネルギー補正に相当する。よって、MP1のエネルギーはハートリー=フォックエネルギーに等しい。
- 2次:2電子励起配置に由来する電子相関の補正エネルギーに相当する。ポスト-ハートリー=フォックとして意味を成す最小の摂動項であるため、実用面ではMP2以上のみが用いられる。
2次(MP2)[2]、3次(MP3)[3][4]、および4次(MP4)[5]のメラー=プレセット計算が小さな系を計算する際に使われる標準的レベルであり、多くの計算化学コードに実装されている。さらに高次のMP計算(一般的にはMP5のみ[6])は一部のコードで可能である。しかしながら、それらの計算コストのためほとんど使われることはない。
MP摂動論の系統的研究では、高次で必ずしも収束的な理論ではないことが示されている。収束は、化学系あるいは基底関数系に依って、遅かったり、速かったり、振動したり、規則正しかったり、高度に不規則だったり、単に収束しなかったりする[7]。1次および2次(MP2)波動関数についての密度行列は「応答密度」と呼ばれる種類であり、より一般的な「期待値密度」とは異なる[8][9]。応答密度行列の固有値(MP2自然軌道の占有数)はしたがって2よりも大きかったり、負にもなり得る。非物理的数は発散摂動展開のしるしである[10]。
加えて、MP3およびMP4レベルで計算された様々な重要な分子の性質は、小分子についてでさえも、MP2レベルのものよりも良くない[11]。
開殻分子については、MPn-理論は非制限ハートリー=フォック(UHF)参照関数に対してのみ直接的に適用可能である(UHF状態は一般にフォック演算子の固有ベクトルではないため)。しかしながら、得られるエネルギーはしばしば深刻なスピン汚染に悩まされ、大きな誤差を生む。可能なより良い代替法は制限開殻ハートリー=フォック法(ROHF法)に基づいたMP2様の手法の1つを使うことである。残念ながら、ROHF波動関数の任意性のため、多くのROHFに基づくMP2様の手法が存在する[12][13](例えば、HCPT[14]、ROMP[15]、RMP[16]〔ROHF-MBPT2とも呼ばれる〕[17]、OPT1およびOPT2[18]、ZAPT[19]、IOPT[20]等[21][22])。ROHFに基づくMP2様の理論の一部は2次を超えるそれらの摂動密度およびエネルギーにおいてスピン汚染に悩まされる[23]。
これらの手法、ハートリー–フォック、非制限ハートリー–フォック、制限ハートリー–フォックは単一の行列式波動関数を使用する。多配置自己無撞着場(MCSCF)法は複数の行列式を使用し、非摂動演算子のために使うことができるが、一意的ではなく、完全活性空間摂動論(CASPT2)[24]、多配置準縮退摂動論(MCQDPT)[25][26]といった非常にたくさんの手法が開発されている[27]。残念なことに、MCSCFに基づく手法は摂動級数の発散がないわけではない[28]。
通常のMP2相関エネルギーは同スピン(平行スピン、same spin、SS)と逆スピン(反平行スピン、opposite spin、OS)成分の和である:
グリメは各項の寄与をスケーリングする手法(spin-component scaling MP2, SCS-MP2)を提案した[29]:
スケーリング係数はQCISD/QCISD(T)の結果と合うように決定され、これにより精度が改善される。さらには同スピン成分を完全に無視する手法(scaled opposite spin MP2, SOS-MP2)も提案されている[30]:
通常のMP2では4中心電子反発積分
の計算にコストがかかる。Resolution of the Identity (RI) 近似は、補助基底関数を導入し4中心積分を3中心積分と2中心積分に置き換えることで高速化する手法である[31]。
展開の仕方には以下の4種類がある[32]:
ここで、t, u, v, w は補助基底関数である。3中心および2中心の重なりおよび電子反発積分は以下の通り:
RI近似を適用した場合 RI-MP2 のように記される。
Pople, J. A.; Seeger, R.; Krishnan, R. (1977). “Variational configuration interaction methods and comparison with perturbation theory”. International Journal of Quantum Chemistry 12 (S11): 149–163. doi:10.1002/qua.560120820. Pople, John A.; Binkley, J. Stephen; Seeger, Rolf (1976). “Theoretical models incorporating electron correlation”. International Journal of Quantum Chemistry 10 (S10): 1–19. doi:10.1002/qua.560100802. Krishnan, Raghavachari; Pople, John A. (1978). “Approximate fourth-order perturbation theory of the electron correlation energy”. International Journal of Quantum Chemistry 14 (1): 91–100. doi:10.1002/qua.560140109. Raghavachari, Krishnan.; Pople, John A.; Replogle, Eric S.; Head-Gordon, Martin (1990). “Fifth order Moeller-Plesset perturbation theory: comparison of existing correlation methods and implementation of new methods correct to fifth order”. The Journal of Physical Chemistry 94 (14): 5579–5586. doi:10.1021/j100377a033. Leininger, Matthew L.; Allen, Wesley D.; Schaeferd, Henry F.; Sherrill, C. David (2000). “Is Moller–Plesset perturbation theory a convergent ab initio method?”. J. Chem. Phys. 112 (21): 9213–9222. Bibcode: 2000JChPh.112.9213L. doi:10.1063/1.481764. Handy, Nicholas C.; Schaefer, Henry F. (1984). “On the evaluation of analytic energy derivatives for correlated wave functions”. The Journal of Chemical Physics 81 (11): 5031. Bibcode: 1984JChPh..81.5031H. doi:10.1063/1.447489. Wiberg, Kenneth B.; Hadad, Christopher M.; Lepage, Teresa J.; Breneman, Curt M.; Frisch, Michael J. (1992). “Analysis of the effect of electron correlation on charge density distributions”. The Journal of Physical Chemistry 96 (2): 671. doi:10.1021/j100181a030. Helgaker, Trygve; Poul Jorgensen; Jeppe Olsen (2000). Molecular Electronic Structure Theory. Wiley. ISBN 978-0-471-96755-2 Hubač, Ivan; Čársky, Petr (1980). “Correlation energy of open-shell systems. Application of the many-body Rayleigh-Schrödinger perturbation theory in the restricted Roothaan-Hartree-Fock formalism”. Physical Review A 22 (6): 2392–2399. Bibcode: 1980PhRvA..22.2392H. doi:10.1103/PhysRevA.22.2392. Lauderdale, Walter J.; Stanton, John F.; Gauss, Jürgen; Watts, John D.; Bartlett, Rodney J. (1991). “Many-body perturbation theory with a restricted open-shell Hartree—Fock reference”. Chemical Physics Letters 187 (1–2): 21–28. Bibcode: 1991CPL...187...21L. doi:10.1016/0009-2614(91)90478-R. Murray, Christopher W.; Handy, Nicholas C. (1992). “Comparison and assessment of different forms of open shell perturbation theory”. The Journal of Chemical Physics 97 (9): 6509. Bibcode: 1992JChPh..97.6509M. doi:10.1063/1.463680. Murray, Christopher; Davidson, Ernest R. (1992). “Different forms of perturbation theory for the calculation of the correlation energy”. International Journal of Quantum Chemistry 43 (6): 755. doi:10.1002/qua.560430604. Fletcher, Graham D; Gordon, Mark S; Bell, Robert S (2002). “Gradient of the ZAPT2 energy”. Theoretical Chemistry Accounts: Theory, Computation, and Modeling (Theoretica Chimica Acta) 107 (2): 57. doi:10.1007/s00214-001-0304-z. Roos, Bjrn O; Andersson, Kerstin; Flscher, Markus P; Malmqvist, Per-ke; Serrano-Andrs, Luis; Pierloot, Kristin; Merchn, Manuela (1996). “Multiconfigurational Perturbation Theory: Applications in Electronic Spectroscopy”. Advances in Chemical Physics. Advances in Chemical Physics. 93. pp. 219. doi:10.1002/9780470141526.ch5. ISBN 978-0-470-14152-6 Nakano, Haruyuki (1993). “Quasidegenerate perturbation theory with multiconfigurational self-consistent-field reference functions”. The Journal of Chemical Physics 99 (10): 7983–7992. Bibcode: 1993JChPh..99.7983N. doi:10.1063/1.465674. Davidson, Ernest R.; Jarzecki, A. A. (1999). K. Hirao. ed. Recent Advances in Multireference Methods. World Scientific. pp. 31–63. ISBN 978-981-02-3777-6 Stefan Grimme (2003). “Improved second-order Møller–Plesset perturbation theory by separate scaling of parallel- and antiparallel-spin pair correlation energies”. J. Chem. Phys. 118: 9095–9102. doi:10.1063/1.1569242. Yousung Jung; Rohini C. Lochan; Anthony D. Dutoi; Martin Head-Gordon (2004). “Scaled opposite-spin second order Møller–Plesset correlation energy: An economical electronic structure method”. J. Chem. Phys. 121: 9793–9802. doi:10.1063/1.1809602. Martin Feyereisen; George Fitzgerald; Andrew Komornicki (1993). “Use of approximate integrals in ab initio theory. An application in MP2 energy calculations”. Chemical Physics Letters 208: 359–363. doi:10.1016/S0009-2614(98)00862-8. O. Vahtras; J. Almlöf; M. W. Feyereisen (1993). “Integral approximations for LCAO-SCF calculations”. Chemical Physics Letters 5–6: 514–518. doi:10.1016/0009-2614(93)89151-7.