配置間相互作用法

電子構造
原子価結合法
コールソン=フィッシャー理論
一般化された原子価結合
現代原子価結合
分子軌道法
ハートリー=フォック法
半経験的分子軌道法
メラー=プレセット法
配置間相互作用法
結合クラスター法
多配置自己無撞着場
量子化学複合手法(英語版)
量子モンテカルロ
原子軌道による線形結合法
密度汎関数理論
時間依存密度汎関数法
トーマス=フェルミ模型
オービタルフリー密度汎関数理論
電子バンド構造
ほとんど自由な電子モデル
強結合近似
マフィンティン近似
k·p摂動論
空格子近似

配置間相互作用法(はいちかんそうごさようほう、: configuration interaction method、略称: CI 法)は、量子化学において、多電子系におけるボルン-オッペンハイマー近似のもとで非相対論的シュレーディンガー方程式を解くために用いられる線形変分的ポスト-ハートリー-フォック法である。

数学的に「配置」とは、波動関数として用いられるスレイター行列式の線形結合で記述される。軌道の占有数(たとえば(1s)2(2s)2(2p)1...)の観点では、「相互作用」は異なる電子配置(状態)の混ざり合い(相互作用)を意味する。CI計算には必要なCPU時間や巨大なハードウェアが必要なため、CI法の使用は相対的に小さい系に限られる。

ハートリーフォック法では波動関数を1つのスレイター行列式で表すが、CI法では電子相関を考慮するため、スピン軌道(上付きSOで記述される)で構成される配置状態関数(CSF)の線形結合を用いる。

Ψ = I = 0 c I Φ I S O = c 0 Φ 0 S O + c 1 Φ 1 S O + {\displaystyle \Psi =\sum _{I=0}c_{I}\Phi _{I}^{\mathrm {SO} }=c_{0}\Phi _{0}^{\mathrm {SO} }+c_{1}\Phi _{1}^{\mathrm {SO} }+\dotsb }

ここで通常はΨは系の電子基底状態である。その後、変分法によって係数 c I   {\displaystyle c_{I}\ } とその時のエネルギー固有値を求める。

この展開が、適切な対称性の可能なすべての配置状態関数(CSF)を含んでいる場合、これは1粒子基底によって張られた空間で電子のシュレーディンガー方程式を正確に解くFull CI法である。上記の展開における1次項は普通はハートリー-フォック行列式である。他のCSFは、ハートリーフォック行列式から仮想軌道交換されたスピン軌道の数によって分類される。1つのスピン軌道が異なっていたならば、これを1励起行列式で記述する。2つのスピン軌道が異なっていたならば、2励起行列式である。これはCI空間と呼ばれる展開での行列式の数を制限するのに使われる。

打ち切られた (truncated) CI空間は計算時間を省くのに重要である。たとえば、CID法では2励起だけに限られる。CISD法では1励起と2励起だけに限られる。これらのCID法、CISD法は多くの場合で用いられる。

デビッドソン補正大きさについての無矛盾性を補正するために使われる。打ち切られたCI法の問題は、無限に離れた2粒子のエネルギーが1粒子のエネルギーの2倍ではないという大きさの矛盾性 (: size-inconsistency) である。

CI法は一般化行列固有値方程式へとつながる。

H c = e S c , {\displaystyle {\boldsymbol {H}}{\boldsymbol {c}}={\boldsymbol {e}}{\boldsymbol {S}}{\boldsymbol {c}},}

ここでcは係数ベクトル、eは固有値行列であり、ハミルトニアンの行列要素、重なり行列の行列要素はそれぞれ以下のようになる。

H i j = Φ i S O | H ^ | Φ j S O S i j = Φ i S O | Φ j S O {\displaystyle {\begin{aligned}{\boldsymbol {H}}_{ij}&=\langle \Phi _{i}^{\mathrm {SO} }|{\hat {H}}|\Phi _{j}^{\mathrm {SO} }\rangle \\{\boldsymbol {S}}_{ij}&=\langle \Phi _{i}^{\mathrm {SO} }|\Phi _{j}^{\mathrm {SO} }\rangle \end{aligned}}}

スレイター行列式は直交化されたスピン軌道の組から構成されるので Φ i S O | Φ j S O = δ i j {\displaystyle \langle \Phi _{i}^{\mathrm {SO} }|\Phi _{j}^{\mathrm {SO} }\rangle =\delta _{ij}} 、つまり S {\displaystyle {\boldsymbol {S}}} は恒等行列となり上記の行列の方程式は簡単な形になる。

CI法の解は、エネルギー固有値 E j {\displaystyle {\boldsymbol {E}}^{j}} と対応するエネルギー固有ベクトル c I j {\displaystyle {\boldsymbol {c}}_{I}^{j}} である。エネルギー固有値は基底状態といくつか電子励起状態のエネルギーである。よってエネルギー差(励起エネルギー)をCI法から計算することが可能である。打ち切られたCI法の励起エネルギーは一般的に高く見積もられすぎる傾向がある。なぜなら励起状態は基底状態ほど相関していないからである。

参考文献

  • Cramer, Christopher J. (2002). Essentials of Computational Chemistry. Chichester: John Wiley & Sons, Ltd.. pp. 191–232. ISBN 0-471-48552-7 
  • Sherrill, C. David; Schaefer III, Henry F. (1999). “The Configuration Interaction Method: Advances in Highly Correlated Approaches”. In Löwdin, Per-Olov. Advances in Quantum Chemistry. 34. San Diego: Academic Press. pp. 143–269. doi:10.1016/S0065-3276(08)60532-8. ISBN 0-12-034834-9 

関連項目