メラー=プレセット法(メラー=プレセットほう、: Møller–Plesset method)とは、計算化学の分野におけるハートリー=フォック法平均場近似)を超える、より高精度な第一原理計算手法の一つである。レイリー=シュレーディンガーの摂動論を用いて電子相関を考慮する方法であり、摂動のレベルによりMP2、MP3、MP4などと呼ばれる。

主要な着想は1934年にクリスチャン・メラーミルトン・S・プレセットによって発表された[1]

レイリー=シュレーディンガーの多体摂動論

編集

レイリー=シュレーディンガー (RS) 摂動論において、ハミルトニアンを非摂動の参照項   と摂動項   に分割する。

 

ここで、λは摂動の大きさを表すパラメータである。エネルギー と、波動関数 λについて連続的に変化するので、テイラー展開によって、

 

と書くことが出来る。これらを時間独立のシュレーディンガー方程式代入すると、以下の式が得られる。

 

この式を展開、整理して、λについて両辺の係数を比較することで、n次の摂動の式が得られる。

メラー=プレセット法

編集

レイリー=シュレーディンガーの摂動論で導かれた式は一般的な形であり、ハミルトニアンに具体的な形を与える必要がある。

メラー=プレセット法の摂動法としての振舞いはあまり良くなく、真の解に向けて単調に収束するわけではない。すなわち、取り込む摂動の次数を上げても真の解から遠ざかることがありうる。このため、メラー=プレセット法により物性値を正しく予想するためには注意を要する。

最初の定式化

編集

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補正によって与えられる。

EEHF + 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]

 

Resolution of the Identity 近似

編集

通常のMP2では4中心電子反発積分

 

の計算にコストがかかる。Resolution of the Identity (RI) 近似は、補助基底関数を導入し4中心積分を3中心積分と2中心積分に置き換えることで高速化する手法である[31]。 展開の仕方には以下の4種類がある[32]:

 

ここで、t, u, v, w は補助基底関数である。3中心および2中心の重なりおよび電子反発積分は以下の通り:

 

RI近似を適用した場合 RI-MP2 のように記される。

脚注

編集
  1. ^ Møller, Christian; Plesset, Milton S. (1934). “Note on an Approximation Treatment for Many-Electron Systems”. Phys. Rev. 46 (7): 618–622. Bibcode1934PhRv...46..618M. doi:10.1103/PhysRev.46.618. http://link.aps.org/abstract/PR/v46/p618. "この文献は数学に複数の軽微であるが、面倒な問題を含んでいる。MP摂動理論のn次への簡潔な導出については、よい量子力学の教科書を参照されたい。" 
  2. ^ Head-Gordon, Martin; Pople, John A.; Frisch, Michael J. (1988). “MP2 energy evaluation by direct methods”. Chemical Physics Letters 153 (6): 503–506. Bibcode1988CPL...153..503H. doi:10.1016/0009-2614(88)85250-3. 
  3. ^ 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. 
  4. ^ 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. 
  5. ^ 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. 
  6. ^ 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. 
  7. ^ 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. Bibcode2000JChPh.112.9213L. doi:10.1063/1.481764. 
  8. ^ 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. Bibcode1984JChPh..81.5031H. doi:10.1063/1.447489. 
  9. ^ 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. 
  10. ^ Gordon, Mark S.; Schmidt, Michael W.; Chaban, Galina M.; Glaesemann, Kurt R.; Stevens, Walter J.; Gonzalez, Carlos (1999). “A natural orbital diagnostic for multiconfigurational character in correlated wave functions”. J. Chem. Phys. 110 (9): 4199–4207. Bibcode1999JChPh.110.4199G. doi:10.1063/1.478301. https://works.bepress.com/mark_gordon/216/download/. 
  11. ^ Helgaker, Trygve; Poul Jorgensen; Jeppe Olsen (2000). Molecular Electronic Structure Theory. Wiley. ISBN 978-0-471-96755-2 
  12. ^ Glaesemann, Kurt R.; Schmidt, Michael W. (2010). “On the Ordering of Orbital Energies in High-Spin ROHF”. The Journal of Physical Chemistry A 114 (33): 8772–8777. Bibcode2010JPCA..114.8772G. doi:10.1021/jp101758y. PMID 20443582. 
  13. ^ Crawford, T. Daniel; Schaefer, Henry F.; Lee, Timothy J. (1996). “On the energy invariance of open-shell perturbation theory with respect to unitary transformations of molecular orbitals”. The Journal of Chemical Physics 105 (3): 1060. Bibcode1996JChPh.105.1060C. doi:10.1063/1.471951. https://zenodo.org/record/1232940. 
  14. ^ 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. Bibcode1980PhRvA..22.2392H. doi:10.1103/PhysRevA.22.2392. 
  15. ^ Amos, Roger D.; Andrews, Jamie S.; Handy, Nicholas C.; Knowles, Peter J. (1991). “Open-shell Møller—Plesset perturbation theory”. Chemical Physics Letters 185 (3–4): 256–264. Bibcode1991CPL...185..256A. doi:10.1016/S0009-2614(91)85057-4. 
  16. ^ Knowles, Peter J.; Andrews, Jamie S.; Amos, Roger D.; Handy, Nicholas C.; Pople, John A. (1991). “Restricted Møller—Plesset theory for open-shell molecules”. Chemical Physics Letters 186 (2–3): 130–136. Bibcode1991CPL...186..130K. doi:10.1016/S0009-2614(91)85118-G. 
  17. ^ 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. Bibcode1991CPL...187...21L. doi:10.1016/0009-2614(91)90478-R. 
  18. ^ Murray, Christopher; Davidson, Ernest R. (1991). “Perturbation theory for open shell systems”. Chemical Physics Letters 187 (5): 451–454. Bibcode1991CPL...187..451M. doi:10.1016/0009-2614(91)80281-2. 
  19. ^ Lee, Timothy J.; Jayatilaka, Dylan (1993). “An open-shell restricted Hartree—Fock perturbation theory based on symmetric spin orbitals”. Chemical Physics Letters 201 (1–4): 1–10. Bibcode1993CPL...201....1L. doi:10.1016/0009-2614(93)85024-I. https://zenodo.org/record/1253842. 
  20. ^ Kozlowski, P. M.; Davidson, Ernest R. (1994). “Construction of open shell perturbation theory invariant with respect to orbital degeneracy”. Chemical Physics Letters 226 (5–6): 440–446. Bibcode1994CPL...226..440K. doi:10.1016/0009-2614(94)00763-2. 
  21. ^ 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. Bibcode1992JChPh..97.6509M. doi:10.1063/1.463680. 
  22. ^ 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. 
  23. ^ 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. 
  24. ^ 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 
  25. ^ Nakano, Haruyuki (1993). “Quasidegenerate perturbation theory with multiconfigurational self-consistent-field reference functions”. The Journal of Chemical Physics 99 (10): 7983–7992. Bibcode1993JChPh..99.7983N. doi:10.1063/1.465674. 
  26. ^ Granovsky, A. A. (2011). “Extended multi-configuration quasi-degenerate perturbation theory: The new approach to multi-state multi-reference perturbation theory”. J. Chem. Phys. 134 (21): 214113. Bibcode2011JChPh.134u4113G. doi:10.1063/1.3596699. PMID 21663350. 
  27. ^ 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 
  28. ^ Glaesemann, Kurt R.; Gordon, Mark S.; Nakano, Haruyuki (1999). “A study of FeCO+ with correlated wavefunctions”. Physical Chemistry Chemical Physics 1 (6): 967–975. Bibcode1999PCCP....1..967G. doi:10.1039/a808518h. 
  29. ^ 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. 
  30. ^ 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. 
  31. ^ 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. 
  32. ^ 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. 

関連記事

編集