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法の数学的定式化、歴史的背景、およびベンチマーク結果について詳述する。
幾何構造最適化における初期ヘシアン推定の数理的基礎と拡張:Schlegel Hessianについて
2026-01-25
幾何構造最適化の収束効率を決定づける初期ヘシアン推定法(Schlegel Hessian)について、H. Bernhard Schlegelによる1984年の提唱から1997年の重元素への拡張に至るまでの理論的背景、数理的アルゴリズム、および経験的パラメータ決定のプロセスを詳述する。原子価座標系を用いた力場構築と座標変換の数学的定式化に焦点を当てる。
冗長内部座標系における自動化された鞍点探索アルゴリズムの数理と実装:Sellaについて
2026-01-25
ポテンシャルエネルギー曲面(PES)上の一次鞍点(遷移状態)を探索するための新規アルゴリズム『Sella』について、その数学的基礎、座標変換の幾何学的処理、およびNull-Space SQPを用いた制約付き最適化手法を詳述する。特に、自動化を阻害する直線結合角問題へのダミー原子を用いた対処法と、反復的ヘシアン対角化による計算コスト削減効果に焦点を当てる。
振動解析における幾何学的基礎: 質量荷重ヘシアンの導出と配置空間の平坦性に関する考察
2026-01-24
分子振動解析において中心的役割を果たす質量荷重座標系への変換について、その数理的導出、計量テンソルの定義、および空間の平坦性(曲率の不在)に関する幾何学的正当性を詳細に解説する。