1 概要
ハートリー・フォック法は、多電子系を近似的に扱うための量子化学・量子力学の計算法である。電子同士の相互作用を、各電子が感じる平均的な場として組み込み、波動関数とエネルギーを求める。原子、分子、固体の電子構造を記述する基礎理論として広く利用され、より精密な電子相関を導入する手法の出発点にもなる。
1.1 多電子系の問題
多電子系では、電子が多数存在するため、正確な波動関数を直接求めることがきわめて難しい。各電子は原子核からの引力だけでなく、他の電子からの反発も受けるため、運動は強く結び付いている。この相互依存性が、厳密解の計算を複雑にしている。
1.2 平均場近似としての位置づけ
この方法は、複雑な電子間相互作用を個々の電子に対する平均的効果へ置き換える。つまり、各電子は他の電子の詳細な瞬間配置ではなく、全体として平均化されたポテンシャルの中を動くとみなされる。こうした近似により、現実的な計算量で多電子問題を扱える。
1.3 量子化学における役割
量子化学では、ハートリー・フォック法は基礎的な参照計算として重要である。分子軌道の概形や結合の特徴を把握する際に使われるほか、相関エネルギーを別途補う諸手法の基準にもなる。計算化学の多くの方法は、この枠組みを起点として構築されてきた。
2 基本原理
2.1 多電子波動関数
多電子系では、全電子の状態を表す波動関数が必要になる。ハートリー・フォック法では、この波動関数を単一の反対称化された形で近似する。各電子の軌道を組み合わせて全体系を表現するのが特徴である。
2.1.1 反対称性の要請
電子はフェルミ粒子であるため、二つの電子を入れ替えると波動関数は符号が反転しなければならない。この性質はパウリの排他原理と対応している。したがって、同じ量子状態に複数の電子が入ることはできない。
2.1.2 スレーター行列式
反対称性を満たす代表的な表現がスレーター行列式である。これは複数の一電子軌道を行列式にまとめた形式で、電子の交換に対して自動的に符号が変わる。ハートリー・フォック法では、この行列式が基本的な試行波動関数となる。
2.2 自己無撞着の考え方
各電子が感じる平均場は、他の電子の分布に依存する。そのため、最初に仮定した場から計算した軌道を用いて新しい場を作り、これを何度も更新する必要がある。この反復的な手続きが自己無撞着である。
2.2.1 有効ハミルトニアン
多電子ハミルトニアンを直接扱う代わりに、電子ごとの有効ハミルトニアンを用いる。そこには核からの引力項に加え、他電子によるクーロン効果と交換効果が含まれる。結果として、一電子問題の形に近い方程式へ整理される。
2.2.2 反復計算
計算は、初期軌道を与えることから始まる。そこから得られた密度やポテンシャルを使って新しい軌道を解き、差が十分小さくなるまで繰り返す。最終的に、入力と出力が整合した状態が収束解として採用される。
2.3 エネルギー最小化の原理
この方法は、試行波動関数のもとで全エネルギーを最小にするように軌道を決める。変分原理に基づくため、得られるエネルギーは厳密値の上限として振る舞う。したがって、近似の質はエネルギーと軌道の両面から評価される。
3 ハートリー・フォック方程式
3.1 方程式の導出
スレーター行列式を変分原理に代入し、軌道についてエネルギーを最小化すると、ハートリー・フォック方程式が導かれる。これは各軌道に対する一連の連立方程式であり、自己無撞着条件を満たすように定まる。実際には積分方程式あるいは行列方程式として扱われることが多い。
3.2 クーロン項
クーロン項は、ある電子が他の電子の電荷分布から受ける静電反発を表す。これは平均場の主要な部分であり、電子密度に応じて変化する。分子や原子の空間的な広がりを反映する要素でもある。
3.3 交換項
交換項は、電子の反対称性に由来する量子力学的効果である。古典的な電荷反発とは異なり、同じスピンを持つ電子の分布に特有の補正を与える。これにより、単純な静電模型では再現できない構造が現れる。
3.4 フォック演算子
フォック演算子は、各電子が従う有効な一電子演算子である。核引力、クーロン相互作用、交換相互作用をまとめて含み、軌道の方程式を簡潔に表す。自己無撞着計算では、この演算子を繰り返し更新していく。
4 主な変種
4.1 制限ハートリー・フォック法
制限ハートリー・フォック法では、対をなす電子が同じ空間軌道を共有する。主として閉殻系の記述に向いており、電子がすべて対になっている場合に適する。計算構造が比較的単純で、広く使われている。
4.1.1 閉殻系への適用
閉殻系では、各軌道に反対スピンの電子が一組ずつ入る。制限法はこの状況を自然に表現できるため、基底状態の多くで有効である。分子や原子の安定な電子配置の初期近似としても重要である。
4.1.2 対称性の扱い
この手法では、空間的対称性を保ちながら計算を進めやすい。対称性を反映した軌道を用いることで、計算の効率が向上する。加えて、得られる解の解釈もしやすい。
4.2 非制限ハートリー・フォック法
非制限ハートリー・フォック法は、スピン上向きと下向きで異なる軌道を許す。電子対が崩れた系や、磁性をもつ状態の近似に適している。対称性の制約が少ないため、より柔軟な記述が可能である。
4.2.1 開殻系への適用
開殻系では、未対電子が存在する。非制限法は、このような場合に各スピン成分を別々に扱えるので、実際の電子配置に近い描像を与える。励起状態やラジカルにも用いられる。
4.2.2 スピン汚染
この方法では、理想的なスピン量子数が混入してしまうことがある。これをスピン汚染と呼ぶ。計算結果の解釈には注意が必要であり、必要に応じて補正や別手法が検討される。
4.3 制限開殻ハートリー・フォック法
制限開殻ハートリー・フォック法は、開殻系を扱いつつ、特定の対称性やスピン構造を保つよう設計された中間的な枠組みである。制限法と非制限法の利点を折衷する形で用いられる。
4.3.1 一般開殻系の扱い
この形式では、対電子と未対電子を区別しながら、同一の電子配置に対する秩序ある記述を行う。一般的な開殻分子の基底状態や一部の励起状態を安定に扱える。
4.3.2 適用条件
適用には、対象系の電子配置や対称性の条件を確認する必要がある。すべての開殻問題に最適というわけではなく、系の性質に応じて制限法や非制限法と使い分けられる。
5 数値計算と実装
5.1 基底関数系
実装では、軌道を有限個の基底関数の線形結合で表す。基底の選び方は計算精度と負荷に大きく影響する。実用上は、目的に応じた適切な基底セットを選ぶことが重要である。
5.1.1 原子軌道基底
原子軌道基底は、原子中心に局在した関数で構成される。原子や局所的な電子分布の表現に適しており、直感的な解釈がしやすい。分子計算では、原子ごとの組み合わせとして使われることが多い。
5.1.2 分子軌道基底
分子軌道基底では、分子全体に広がる関数を用いる。これにより、結合性や反結合性の特徴が明瞭になる。電子の非局在化を扱うときに有利である。
5.2 行列要素の計算
ハートリー・フォック計算では、基底関数間の行列要素を数値的に評価する。これには一電子項と二電子項の計算が含まれる。計算コストの多くは、後者の取り扱いに費やされる。
5.2.1 一電子積分
一電子積分は、運動エネルギーや核引力に対応する。比較的標準的な形を持ち、基底関数が決まれば計算できる。系の骨格となる情報を与える部分である。
5.2.2 二電子積分
二電子積分は、電子間反発を表すための中心的な量である。組み合わせの数が急増するため、実装上の負担が大きい。効率化のために、対称性や近似技術がしばしば利用される。
5.3 収束判定と反復法
自己無撞着計算では、反復のたびにエネルギーや密度が変化する。これらの変化が十分小さくなった時点で収束とみなす。収束条件の設定は、安定な計算結果を得るうえで重要である。
5.3.1 密度行列の更新
密度行列は、占有軌道から構成される。新しい軌道が得られるたびにこれを更新し、次の反復の出発点とする。こうして、電子分布が徐々に自己整合的な形へ近づく。
5.3.2 収束加速法
収束が遅い場合には、混合手法や補助的なアルゴリズムが使われる。これらは反復の振動や発散を抑え、安定化に寄与する。大規模計算では特に重要である。
6 応用
6.1 原子の電子構造
原子計算では、殻構造や軌道エネルギーの概略を与える。多電子原子に対しても、基本的な電子配置の理解に役立つ。スペクトルやイオン化傾向を考える際の出発点にもなる。
6.2 分子軌道の解析
分子では、結合の形成や電子の分布を調べるために用いられる。得られた軌道は、化学結合の性質や反応性の傾向を理解する手がかりとなる。定性的な議論と数値解析の両方に利用される。
6.3 固体物理への拡張
固体では、周期境界条件や結晶対称性を考慮して拡張される。バンド構造の初期近似としても利用されるが、電子相関が強い系では追加の補正が必要になることが多い。材料計算の基礎的枠組みの一つである。
6.4 反応解析の基礎
化学反応の解析では、反応物と生成物の電子構造を比較する際の基準になる。遷移状態探索や反応座標の検討においても、初期近似として重要である。反応機構の理解を支える土台となる。
7 限界と改良
7.1 電子相関の不足
この方法では、電子が互いの瞬間的な運動を細かく避け合う効果を十分に表しきれない。したがって、相関エネルギーが落ちる。強い相関を示す系では、精度が不十分になることがある。
7.2 相関補正手法
不足分を補うため、さまざまな相関補正法が導入される。これらはハートリー・フォック解を基準として、より現実的なエネルギーや波動関数を与える。実際の計算では、目的に応じて選択される。
7.2.1 摂動法との併用
摂動論を組み合わせると、平均場からのずれを段階的に評価できる。比較的少ない追加計算で精度を高められるため、広く使われている。単純な近似と高度な理論の橋渡しにもなる。
7.2.2 多配置法との比較
多配置法では、複数の電子配置を同時に考慮する。これにより、単一行列式では表しにくい相関を扱いやすくなる。一方で計算負荷は増すため、対象系との兼ね合いが重要である。
7.3 より高精度な手法への発展
ハートリー・フォック法は、より高度な量子化学手法の基盤として位置付けられる。そこから、相関を精密に扱う多様な理論へ発展した。現代の計算化学では、基礎近似と高精度法を組み合わせて用いることが多い。
8 関連する理論
8.1 ハートリー法との違い
ハートリー法は、反対称性を明示的に組み込まない点で異なる。これに対し、ハートリー・フォック法は交換効果を反映し、フェルミ粒子としての電子の性質をより忠実に扱う。したがって、後者のほうが物理的に洗練されている。
8.2 密度汎関数法との関係
密度汎関数法は、電子密度を中心に記述する別系統の理論である。両者は平均場的な側面を共有するが、相関の扱い方が異なる。実務上は、ハートリー・フォック解が比較や参照に使われることが多い。
8.3 多体摂動論との関係
多体摂動論では、平均場状態を出発点として、相互作用の補正を系統的に加える。ハートリー・フォック法は、その摂動展開の基準状態として自然に現れる。理論的な一貫性を保つうえで重要な位置を占める。
9 歴史
9.1 提案の背景
この方法は、多電子原子の計算を現実的に行う必要から発展した。量子力学の成立後、原子や分子の電子構造を扱うための近似法が求められたことが背景にある。計算可能性と物理的妥当性の両立が目標であった。
9.2 フォックによる拡張
フォックは、反対称性と交換効果を理論に取り入れて手法を拡張した。これにより、電子の量子統計に即した記述が完成度を増した。名称にもその貢献が反映されている。
9.3 量子化学への定着
その後、計算機の発展とともに、分子や固体への応用が急速に広がった。現在では、電子構造計算の標準的な出発点として定着している。教育、研究、実務の各場面で基本理論として扱われる。