Home
1589 words
8 minutes
電子相関の二分法:静的相関と動的相関の定義、弊害、および解決手法

last_modified: 2026-01-21

1. 概要 (Overview)#

電子相関エネルギー EcorrE_{corr} は、厳密な非相対論的エネルギー EexactE_{exact} と、Hartree-Fock (HF) 極限エネルギー EHFE_{HF} との差として定義される。 Ecorr=EexactEHFE_{corr} = E_{exact} - E_{HF} 計算化学の実践において、この相関は物理的起源の異なる二つの成分、「静的電子相関 (Static Correlation)」と「動的電子相関 (Dynamic Correlation)」に分類して扱われることが多い。本稿では、それぞれの定義、考慮漏れによる物理的・化学的弊害、および適切な計算手法について記述する。

2. 静的電子相関 (Static Correlation)#

別名:非動的相関 (Nondynamic Correlation)、多配置性 (Multireference character)

2.1 定義 (Definition)#

基底状態の波動関数を記述するために、単一のSlater行列式では不十分であり、複数の電子配置(行列式)がほぼ等しい重みで寄与する状況に由来する相関。 主に、最高被占軌道 (HOMO) と最低空軌道 (LUMO) のエネルギー差が小さい、あるいは縮退している「準縮退 (Near-degeneracy)」系において顕著となる。これは電子間の「長距離的な回避」や「軌道の再編成」に関連する。

2.2 考慮しない場合の弊害 (Consequences of Neglect)#

静的相関を無視する(= RHFなどの単一参照法を強行する)と、定性的 (Qualitative) な誤りが生じる。

  • 結合解離の記述破綻: 閉殻分子(例: H2H_2)を解離させる際、RHF波動関数は解離極限において「イオン性項 (H+H+H^- + H^+)」と「共有結合性項 (H+HH^\cdot + H^\cdot)」を等しい重みで含んでしまう。結果として、解離エネルギーが過大評価されるだけでなく、正しい解離生成物(中性ラジカル)が得られない。
  • 遷移状態・励起状態の不適切な記述: 化学反応の遷移状態や、ビラジカル、遷移金属錯体など、フロンティア軌道が近接する系において、波動関数のトポロジー自体が誤って記述される。
  • スピン対称性の破れ: UHF(非制限HF)を用いることでエネルギー的な破綻は回避できる場合があるが、波動関数がスピン固有状態ではなくなり(スピン汚染)、物理的に純粋な状態として扱えなくなる。

2.3 解決手法 (Methodology)#

波動関数の定義自体を「多配置 (Multi-configuration)」にする必要がある。

  • CASSCF (Complete Active Space SCF): 特定の軌道空間(Active Space)内で、電子配置の全組み合わせを考慮する。静的相関の回復に最も標準的に用いられる。
  • MCSCF (Multi-Configurational SCF): CASSCFの一般形。必要な配置のみを選別して最適化する。
  • GVB (Generalized Valence Bond): 原子価結合法の概念を取り入れ、結合電子対の相関を局所的に記述する。

3. 動的電子相関 (Dynamic Correlation)#

3.1 定義 (Definition)#

電子同士が互いにクーロン反発を感じ、瞬間的に互いを避け合う運動に由来する相関。 Hartree-Fock法では、電子は「他の電子が作る平均場 (Mean Field)」の中を運動すると仮定されるため、r120r_{12} \to 0 における特異点(クーロンホール)や、瞬間的な位置相関が記述できない。これは短距離相互作用に基づく。

3.2 考慮しない場合の弊害 (Consequences of Neglect)#

動的相関を無視しても波動関数の形(定性的な性質)は正しい場合が多いが、定量的 (Quantitative) な精度が欠如する。

  • 化学精度の欠如: 結合エネルギー、活性化エネルギー、反応熱などの算出において、数 kcal/mol ~ 数十 kcal/mol の誤差が生じる。「化学精度 (Chemical Accuracy, 1\sim 1 kcal/mol)」を達成するためには、動的相関の補正が不可欠である。
  • 構造パラメータの誤差: 一般にHF法は結合長を過小評価し、振動数を過大評価する傾向がある。弱い相互作用(分散力など)は動的相関そのものであるため、これを無視するとファンデルワールス錯体などは結合しなくなる。

3.3 解決手法 (Methodology)#

単一の参照配置(通常はHF解)を出発点とし、励起配置を摂動的あるいは反復的に取り込むことで記述する。

  • Coupled Cluster (CCSD(T)): 指数関数展開により高次励起を取り込む。動的相関記述の「ゴールドスタンダード」とされる。
  • MPn (Møller-Plesset Perturbation Theory): 摂動論に基づき、MP2, MP4などで相関エネルギーを推定する。低コストだが収束性に課題がある場合も。
  • CI (Configuration Interaction): 配置間相互作用法。CISDなどが用いられるが、サイズ無矛盾性(Size Consistency)の問題があるため、現在はCC法が主流。

4. 両者の統合的解決 (Integrated Approaches)#

現実の複雑な化学系(遷移金属触媒、光化学反応など)では、静的相関と動的相関の両方が重要となる。これらを同時に扱う手法は「多参照相関手法 (Multi-Reference Correlation Methods)」と呼ばれる。

手法カテゴリ代表的手法アプローチ
MR-PTCASPT2, NEVPT2CASSCF波動関数(静的)を参照とし、摂動論で動的相関を加える。
MR-CIMRCI多参照波動関数に対してCIを行う。極めて高コストだが高精度。
MR-CCMR-CCSD, Mk-MRCC多参照Coupled Cluster法。理論的に複雑で実装・収束が困難な場合が多い。

5. 結論 (Conclusion)#

  • 静的相関は「正しい出発点(参照関数)」を得るための要件であり、欠落すると定性的に誤った物理(誤った解離極限など)を導く。
  • 動的相関は「精密なエネルギー」を得るための要件であり、欠落すると定量的に誤った値(大きな誤差)を導く。

計算対象が「平衡構造付近の閉殻有機分子」であれば動的相関の考慮(CCSD(T)等)のみで十分であるが、「結合解離、遷移状態、遷移金属」を含む場合は静的相関の考慮(CASSCF等)を前提としなければならない。

電子相関の二分法:静的相関と動的相関の定義、弊害、および解決手法
https://ss0832.github.io/posts/20260121_electron_correlation/
Author
ss0832
Published at
2026-01-21
License
CC BY-NC-SA 4.0

Related Posts

結合クラスター法(Coupled-Cluster Theory)の理論的枠組みと『ゴールドスタンダード』への道程:CCSD(T)法の本質と限界
2026-01-13
R. J. BartlettとM. Musiałによる包括的レビュー(Rev. Mod. Phys. 79, 291, 2007)に基づき、現代量子化学における最高精度の標準的手法である結合クラスター法(CC法)を解説。指数関数的アンザッツによるサイズ示量性の保証、CCSD(T)の統計的誤差(結合長・振動数・エネルギー)、および単参照性を前提とした限界について、理論化学の視点から詳述する。
ポテンシャル曲面交差上の最小エネルギー点を特定するための直接法:射影演算子を用いた制約なし最適化アルゴリズム
2026-01-04
M. J. Bearpark, M. A. Robb, H. B. Schlegelによる論文 *Chem. Phys. Lett.* **1994**, *223*, 269-274 の包括的解説。光化学反応の機構解明に不可欠な円錐交差(Conical Intersection)および項間交差の最低エネルギー点を、Lagrange未定乗数法を用いずに探索する「直接法(Direct Method)」の理論的導出、アルゴリズムの実装、およびベンチマーク結果について詳述する。
分岐平面更新法(UBP法)による最小エネルギー円錐交差の探索:非断熱結合ベクトルを用いない効率的アルゴリズムの理論と実装
2026-01-04
Satoshi Maeda, Koichi Ohno, Keiji Morokumaによる論文 *J. Chem. Theory Comput.* **2010**, *6*, 1538-1545 の詳細な学術的解説。最小エネルギー円錐交差(MECI)探索において計算コストのボトルネックとなる非断熱結合ベクトル(CDV)の算出を回避し、分岐平面の更新によって効率的な最適化を実現するUBP法の数学的定式化、歴史的背景、およびベンチマーク結果について詳述する。
Kabsch アルゴリズムによる最小二乗重ね合わせと RMSD 算出:数学的導出・歴史的背景・実利的応用
2026-03-16
分子構造の比較において中心的な役割を果たす Kabsch アルゴリズムについて、特異値分解(SVD)に基づく数学的導出、歴史的成立過程、実装上の注意点(キラリティ処理を含む)、および計算化学・構造生物学における実利的成果を包括的に解説する。
O1NumHess: O(1) 勾配評価によるセミ数値的ヘシアン行列算出法の数理的構造と実装
2026-02-15
分子系のヘシアン行列算出において、原子数に依存しない定数オーダー O(1) の勾配計算回数で高精度な近似を実現する O1NumHess アルゴリズムについて、その数理的背景(Off-Diagonal Low-Rank 性質)、実装の詳細、および計算精度とコストに関する実証的評価を包括的に解説する。
密度汎関数理論における交換相関汎関数の系統的最適化:B97汎関数の数学的構造と有限級数展開による拡張性
2026-02-04
Axel D. Beckeによる1997年の論文 (J. Chem. Phys. 107, 8554) に基づき、B97汎関数の設計思想とその数学的定式化について詳述する。特に、勾配補正項の有限級数展開(冪級数展開)を用いた系統的な最適化手法、半無限区間の有限変数への写像、およびG2テストセットを用いたパラメータ決定プロセスに焦点を当てる。また、この手法が示唆するGGA/正確交換混合形式の精度の限界と、物理的妥当性と過剰適合(オーバーフィッティング)の境界に関する考察を行う。