Home
1167 words
6 minutes
相対論的量子化学:スカラー効果、スピン軌道相互作用、およびECP

last_modified: 2026-01-21

1. 概要 (Overview)#

非相対論的量子化学(シュレーディンガー方程式)は、光速 cc を無限大 (cc \to \infty) と仮定している。 しかし、原子番号 ZZ が大きい重原子の内殻電子は、原子核の強い引力により光速に近い速度(vZ/137cv \approx Z/137 c)で運動する。これにより、質量増大などの相対論効果が無視できなくなる。 計算化学において、相対論効果は主にスカラー相対論効果 (Scalar Relativistic Effects)スピン軌道相互作用 (Spin-Orbit Coupling, SOC) の二つに分類して扱われる。

2. 物理的現象と化学的帰結 (Phenomena and Consequences)#

2.1 スカラー相対論効果#

電子の質量 mm が速度 vv の増加に伴い増大する効果 (m=m0/1(v/c)2m = m_0 / \sqrt{1-(v/c)^2}) と、ダーウィン項(電子の位置のゆらぎ)などのスピンに依存しない補正。

  • s, p軌道の収縮: 質量増大により内殻のs軌道が原子核に引き寄せられ、収縮する。直交性の制約により、外殻のs, p軌道も連動して収縮する。
    • : 金 (AuAu) の色が黄色である理由(5d-6sギャップの縮小)、水銀 (HgHg) が液体である理由(6s軌道の安定化による結合性低下)。
  • d, f軌道の拡張: 内殻s, p軌道の収縮に伴い、原子核電荷の遮蔽効果が高まるため、d, f軌道は逆に外側へ広がり、エネルギー的に不安定化(上昇)する。

2.2 スピン軌道相互作用 (SOC)#

電子のスピン角運動量 s\vec{s} と軌道角運動量 l\vec{l} が磁気的に相互作用する効果。ベクトル的な補正であり、スピン対称性を破る。

  • 縮退の分裂 (Splitting): 非相対論では縮退している軌道(例: p軌道, d軌道)が、全角運動量量子数 j=l±sj = l \pm s に応じてエネルギー分裂を起こす。
    • 重要性: 重原子を含む系の励起状態、項間交差 (Intersystem Crossing)、燐光 (Phosphorescence) の計算において支配的な役割を果たす。

3. 実装階層:ハミルトニアンの近似 (Hierarchy of Hamiltonians)#

計算コストと精度のバランスに応じた、3つの主要なアプローチが存在する。

3.1 レベル1: 有効内殻ポテンシャル (Effective Core Potential, ECP)#

内殻電子を明示的に扱わず、パラメータ化されたポテンシャル(擬ポテンシャル)に置き換える手法。

  • 原理: 内殻電子の相対論効果をポテンシャルの中に「焼き込む」ことで、価電子のみを非相対論的シュレーディンガー方程式で計算する。
  • 利点: 計算コストが劇的に低下する。基底関数系も小さくて済む(LanL2DZ, SDDなど)。
  • 欠点: 内殻物性(X線吸収など)は計算不可。コア近傍のトポロジーは失われる。

3.2 レベル2: スカラー相対論的全電子計算 (Scalar Relativistic All-electron)#

全電子を扱いながら、相対論効果のうち「スカラー項」のみを取り込んだハミルトニアンを用いる。

  • Douglas-Kroll-Hess (DKH) 法: 外部ポテンシャルの特異点を取り除くユニタリ変換を用いる。DKH2, DKH4などの次数がある。
  • ZORA (Zeroth-Order Regular Approximation): ディラック方程式の小成分を正則化して近似する。内殻付近の記述に優れる。
  • 用途: 重原子を含む構造最適化や結合エネルギー計算の標準。

3.3 レベル3: 4成分/2成分ディラック方程式 (4-component / 2-component Dirac)#

スピン軌道相互作用を含めた完全な相対論的取り扱い。

  • 4成分 (Dirac-Coulomb): 電子(大成分スピノル)と陽電子(小成分スピノル)を露わに扱う。コストは非相対論の数十倍~百倍以上。
  • 2成分 (X2C, BSSなど): 陽電子成分を厳密に分離し、電子成分のみのハミルトニアンを構築する。4成分とほぼ同等の精度で低コスト。
  • 用途: 精密な分光予測(g値、NMR遮蔽定数、励起エネルギー)。

4. 結論 (Conclusion)#

  • 第1~3周期 (H - Ar): 相対論効果は無視可能(シュレーディンガー方程式で十分)。
  • 第4周期 (K - Kr): 遷移金属の精密計算にはスカラー効果あるいはECPが推奨される。
  • 第5周期以降 (Rb - ): 相対論効果の考慮は必須である。これを無視した計算結果(非相対論)は、定性的にも誤り(結合長の過大な誤差、誤った電子配置)を含むリスクが高い。
  • 励起状態・分光: 重原子が関与する項間交差を扱う場合は、SOCを含めた計算(2成分以上または摂動論的SOC)が必要となる。
相対論的量子化学:スカラー効果、スピン軌道相互作用、およびECP
https://ss0832.github.io/posts/20260121_compchem_relativistic_effect/
Author
ss0832
Published at
2026-01-21
License
CC BY-NC-SA 4.0

Related Posts

【計算化学】ガウス基底関数系の数理と歴史:多中心積分の克服から現代の多種多様な選択肢へ
2025-12-30
量子化学計算の根幹を支えるガウス基底関数系(Gaussian Basis Sets)について、その数理的背景と歴史的発展を包括的に解説する。STOからGTOへの転換点となったBoysの定理、Pople系・Dunning系という古典的潮流に加え、現代DFTのデファクトスタンダードであるKarlsruhe系(def2)やDFT特化型のJensen系(pc-n)、重原子計算に不可欠なECP基底(LanL2DZ, SDD)までを網羅する。計算コストと精度の狭間で進化を続けてきた基底関数の設計思想を詳述する。
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)によって提案された、冗長内部座標多様体上の測地線を利用した構造最適化手法の数理的定式化、アルゴリズムの実装、およびその有効性について詳細に解説する。従来のニュートン法との幾何学的差異と、微分幾何学的アプローチによる収束性の改善に焦点を当てる。