1167 words
6 minutes
スピンの取り扱い:制限法 (RHF) と非制限法 (UHF) の定義と弊害
last_modified: 2026-01-21
1. 概要 (Overview)
Hartree-Fock法およびDFTにおいて、スピン(アップ)電子と スピン(ダウン)電子の空間軌道をどのように扱うかによって、手法は大きく二つに分類される。
- 制限法 (Restricted): とが同じ空間軌道を共有すると仮定する。
- 非制限法 (Unrestricted): とが異なる空間軌道を占有することを許容する。
この選択は、基底状態のエネルギー精度だけでなく、波動関数の物理的な対称性(スピン固有状態かどうか)に直結する。
2. 制限閉殻 Hartree-Fock (RHF: Restricted Closed-Shell HF)
2.1 定義
電子数が偶数の閉殻系(全ての電子が対になっている状態)に対し、電子と電子に同一の空間軌道関数 を割り当てる手法。
2.2 特徴
- 計算コスト: 求めるべき変分パラメータ(軌道係数)が1セットで済むため、計算が高速である。
- スピン対称性: 全スピン角運動量演算子 の固有関数となることが保証される。
- 一重項 (Singlet) ならば、厳密に となる。
- 適用限界: 結合解離の記述に失敗する(H2 Catastrophe)。 水素分子 () を引き剥がしていくと、RHF波動関数は無限遠でも「イオン性項 ()」と「共有結合性項 ()」を等しい重みで含んでしまうため、解離エネルギーが物理的に誤った高い値に収束する。
3. 非制限 Hartree-Fock (UHF: Unrestricted HF)
3.1 定義 (Pople-Nesbet方程式)
電子と電子に対し、互いに独立した異なる空間軌道関数 を割り当てる手法。 これにより、変分自由度が増加し、エネルギーは常に となる。
3.2 物理的効果:スピン分極 (Spin Polarization)
不対電子(ラジカル)が存在する場合、その不対電子のスピンとの交換相互作用により、対を作っている内殻電子であっても と で感じるポテンシャルが異なる。 UHFはこの違いを反映し、同じ軌道上の電子対を空間的に分離(分極)させることができる。
3.3 特徴と弊害
- 結合解離の成功: の解離において、電子を片方の原子に、電子をもう片方の原子に局在化させることができるため、正しい解離極限エネルギー () へと収束する。
- スピン汚染 (Spin Contamination):
これが最大の欠点である。UHF波動関数は、一般に の固有関数ではなくなる。
本来のスピン多重度(例:二重項 Doublet)に、より高次のスピン状態(例:四重項 Quartet)が混入してしまう。
- 判定基準: 期待値 が理論値から 10%以上 乖離している場合、その計算結果(特にエネルギーや構造)は物理的に信頼できない可能性がある。
4. 制限開殻 Hartree-Fock (ROHF: Restricted Open-Shell HF)
4.1 定義
開殻系(ラジカル等)に対して、対になっている電子は制限軌道(同じ空間軌道)に入れ、不対電子のみを単独の軌道に配置する手法。
4.2 特徴
- メリット: の固有関数であることが保証され、スピン汚染がゼロである。スペクトル計算など、スピン対称性が重要な場合に必須となる。
- デメリット: 一つの軌道決定方程式(Fock演算子)で解けるRHF/UHFと異なり、閉殻部と開殻部で異なる演算子を扱う必要があり、実装と収束計算が複雑になる。また、スピン分極効果(Core Polarization)を記述できないため、UHFよりエネルギーが高くなることがある。
結論 (Selection Guide)
- 通常の安定有機分子(一重項): 迷わず RHF (R-DFT) を使用する。
- ラジカル、遷移金属、結合解離過程:
まず UHF (U-DFT) を使用する。
- 結合解離エネルギーを求めるにはUHFが必須である(RHFでは解離しない)。
- 出力された
Spin S2の値を確認し、スピン汚染が深刻であれば ROHF の使用、あるいは事後的なスピン射影法 (Spin Projection) を検討する。
- DFTの場合: DFTは電子相関を含んでいるため、HF法に比べてスピン汚染の影響が軽微であることが多い。したがって、開殻系計算では U-DFT が事実上の標準である。
スピンの取り扱い:制限法 (RHF) と非制限法 (UHF) の定義と弊害
https://ss0832.github.io/posts/20260121_compchem_rks_uks/ Related Posts
非制限法におけるスピン汚染(Spin Contamination)の解析:S^2 期待値の評価とBroken Symmetry解の物理的意義
2026-01-02
量子化学計算(UHF/UDFT)において生じるスピン汚染の理論的背景と、全スピン角運動量 S^2 に基づく定量的な判定基準を詳説する。特に、数値的なアーチファクトと物理的実体を持つBroken Symmetry(対称性の破れ)解との峻別、およびGaussian等の計算パッケージにおける波動関数の安定性解析(stable=opt)の挙動について、SzaboやJensenらの成書を参照しつつ論じる。
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/正確交換混合形式の精度の限界と、物理的妥当性と過剰適合(オーバーフィッティング)の境界に関する考察を行う。
デカルト座標系における角度ポテンシャル勾配の特異点除去:数値的正則化手法の数理的導出と実装
2026-02-02
分子動力学および構造最適化計算において、3点間角度ポテンシャルの評価時に発生する特異点(ゼロ除算および桁落ち)の問題を数学的に解析し、テイラー展開を用いた正則化手法(Regularization)を詳細に論じる。Szabo-OstlundのBoys関数処理との類似性、誤差解析、およびPythonによる参照実装を含め、学術的・技術的観点から包括的に解説する。
冗長内部座標系における測地線アプローチを用いた幾何構造最適化の加速
2026-01-29
Hermesら(2021)によって提案された、冗長内部座標多様体上の測地線を利用した構造最適化手法の数理的定式化、アルゴリズムの実装、およびその有効性について詳細に解説する。従来のニュートン法との幾何学的差異と、微分幾何学的アプローチによる収束性の改善に焦点を当てる。