多配置自己無撞着場のソースを表示

{{電子構造論}}
'''多配置自己無撞着場'''（たはいちじこむどうちゃくば、{{lang-en|Multi-configurational Self-consistent Field}}、略称: '''MCSCF'''）は、[[量子化学]]において、全電子[[波動関数]]の取り扱いが[[ハートリー＝フォック法]]や[[密度汎関数理論]]では不十分な場合（低励起状態で擬縮退している[[分子]]の[[基底状態]]や、[[化学結合|結合]]が乖離している場合など）に、分子の質的に正しい参照状態を生成するための[[量子化学]]の手法である。MCSCFは、[[配置状態関数]]（CSF）または[[スレーター行列式]]の[[線形結合]]を使用して、原子または分子の正確な電子波動関数を近似する。さらにMCSCFでは、CSFまたはスレーター行列式と分子軌道の基底関数の両方の係数を変化させて、可能な限り低いエネルギーの全電子波動関数を得る。この方法は、[[配置間相互作用法]]（分子軌道係数は変化させず、CSFまたはスレーター行列式の展開係数を最適化する）とハートリー＝フォック法（スレーター行列式は1つだけだが、分子軌道係数を変化させる）を組み合わせたものと考えることができる。

MCSCF波動関数は、[[多参照配置間相互作用]]（MRCI）や、[[完全活性空間摂動論]]（CASPT2）などの多参照摂動論の参照状態としてよく用いられる。これらの手法は、非常に複雑な化学的状況を扱うことができ、計算機の性能が許せば、他の手法が失敗した場合においても、分子の基底状態や励起状態を計算しうる強力な手法である。

== 概要 ==
最も単純な単結合である[[水素#水素分子|H<sub>2</sub>分子]]の場合を例に説明する。この場合、[[分子軌道]]は常に2つの原子核AとB上にそれぞれ位置する[[原子軌道]]によく似た関数 χ<sub>''iA''</sub>およびχ<sub>''iB''</sub>を用いて描くことができる。

: <math>\varphi_i = N_i(\chi_{iA} \pm \chi_{iB}),</math>

ここで、''N''<sub>''i''</sub>は規格化定数である。H<sub>2</sub>分子の平衡核間距離での基底状態の波動関数は (''φ''<sub>1</sub>)<sup>2</sup>という配置で支配されており、これは分子軌道''φ''<sub>1</sub>がほぼ二重に占有されていることを意味する。ハートリー＝フォック（HF）モデルでは、この軌道が二重に占有されていると「仮定」しているため、全体の波動関数は次のようになる。

: <math>\Phi_1 = \varphi_1(\mathbf{r}_1) \varphi_1(\mathbf{r}_2) \Theta_{2,0},</math>

ここで、 <math>\Theta_{2,0}</math> は、2つの電子の一重項 (''S''&nbsp;=&nbsp;0) のスピン関数である。この場合の分子軌道 ''φ''<sub>1</sub> は、両原子の[[s軌道|1s原子軌道]]の和、すなわち''N''<sub>1</sub>(1s<sub>A</sub>&nbsp;+&nbsp;1s<sub>B</sub>) となる。上の式を原子軌道に展開すると、次のようになる。

: <math>\Phi_1 = N_1^2 \left[ 1s_A(\mathbf{r}_1) 1s_A(\mathbf{r}_2) + 1s_A(\mathbf{r}_1) 1s_B(\mathbf{r}_2) + 1s_B(\mathbf{r}_1) 1s_A(\mathbf{r}_2) + 1s_B(\mathbf{r}_1) 1s_B(\mathbf{r}_2) \right] \Theta_{2,0}.</math>

このHFモデルは、 H<sub>2</sub>の結合長が約0.735&nbsp;Å（実験値は0.746&nbsp;Å）、結合エネルギーが350&nbsp;kJ/mol（84&nbsp;kcal/mol）（実験値は432&nbsp;kJ/mol [103&nbsp;kcal/mol]）<ref>{{cite journal |authors=[[ゲルハルト・ヘルツベルク|G. Herzberg]], A. Monfils |title=The dissociation energies of the H<sub>2</sub>, HD, and D<sub>2</sub> molecules |journal=[[Journal of Molecular Spectroscopy|J. Mol. Spec.]] |volume=5 |issue=1–6 |page=482–498 |year=1961}}</ref>と、平衡核間距離付近において H<sub>2</sub>分子を合理的に記述している。このように、HFモデルは、分子構造が平衡核間距離にある閉殻系を非常によく記述する。しかし、結合が乖離しているような状況では、片方の原子に2つの電子が残ったスレーター行列式も波動関数に含めることが重要になる。これは、H + Hよりもはるかに大きなエネルギーを持つH<sup>+</sup> + H<sup>−</sup>への解離に対応している。

そのため、HFモデルは解離過程の開殻的な状態を記述するためには使用できない。この問題を解決する最も簡単な方法は、Ψ<sub>1</sub>に含まれるそれぞれの項の前に係数を導入することである。

: <math>\Psi_1 = C_\text{ion} \Phi_\text{ion} + C_\text{cov} \Phi_\text{cov},</math>

これは[[化学結合]]を記述する[[原子価結合法]]の基本となる。係数''C''<sub>ion</sub>と''C''<sub>cov</sub>を変化させると、波動関数は正しい形になり、乖離極限では ''C''<sub>ion</sub> = 0、平衡状態では''C''<sub>ion</sub>と''C''<sub>cov</sub>が同程度になる。しかし、この記述は非直交基底関数を使用しているため、数学的構造が複雑化する。その代わり、MCSCFでは、直交する分子軌道を用いることで、多配置の波動関数を構成している。反結合軌道を導入すると
: <math>\phi_2 = N_2 (1s_A - 1s_B),</math>
となり、H<sub>2</sub>の全波動関数は、結合軌道と反結合軌道から構成される配置状態関数の線形結合として書くことができる。

: <math>\Psi_\text{MC} = C_1 \Phi_1 + C_2 \Phi_2,</math>

ここで、Φ<sub>2</sub>は[[電子配置]] (φ<sub>2</sub>)<sup>2</sup> である。このH<sub>2</sub>化学結合の多配置波動関数による記述では、平衡核間距離に近い状態では''C''<sub>1</sub>&nbsp;=&nbsp;1、 ''C''<sub>2</sub>&nbsp;=&nbsp;0となり、大きく結合が乖離している場合には''C''<sub>1</sub>と''C''<sub>2</sub>の絶対値がほぼ同じになる
<ref>
  {{cite book
  | last = McWeeny
  | first = Roy
  | authorlink = ロイ・マクウィーニー
  | title = [[チャールズ・クールソン|Coulson]]'s Valence
  | publisher = Oxford University Press
  | year = 1979
  | location = Oxford
  | pages = 124–129
  | isbn = 0-19-855145-2}} 
</ref>。

== 完全活性空間SCF法（{{vanchor|CASSCF}}） ==
特に重要なMCSCF法は、完全活性空間SCF法（CASSCF）である。CASSCFでは、波動関数を構成するCSFに[[完全活性空間]]を指定する、すなわち、活性空間に属する電子すべての配置を考慮する（完全最適化反応空間（FORS-MCSCF）とも呼ばれる）。例えば、[[一酸化窒素|NO]]分子に対してCASSCF(11,8) を定義すると、11個の価電子を8個の分子軌道（すなわち16個のスピン軌道）に詰めて構成されるすべての配置状態関数を線形結合に含めることになる<ref>
{{cite book
  | last = Jensen
  | first = Frank
  | title = Introduction to Computational Chemistry
  | url = https://archive.org/details/introductiontoco00jens_022
  | url-access = limited
  | publisher = John Wiley and Sons
  | year = 2007
  | pages = [https://archive.org/details/introductiontoco00jens_022/page/n153 133]–158
  | location = Chichester, England
  | isbn = 0-470-01187-4}}
</ref><ref>{{cite book
  | last = Cramer  | first = Christopher J.
  | title = Essentials of Computational Chemistry
  | publisher = John Wiley & Sons, Ltd.  | year = 2002  | location = Chichester
  | pages = 191–232  | isbn = 0-471-48552-7}}</ref>。

== 制限活性空間SCF（{{vanchor|RASSCF}}） ==
CASSCFでは活性軌道の数が増えるとCSFの数も増え、計算コストも増加するため、使用するCSFをより少なくすることが望ましい場合がある。このための1つの方法は、活性空間を分割し、その部分空間における電子の数の最大数と最小数を制限することである。RASSCFはこのような方法の一つであり、活性空間を3つの部分空間に分割し各部分空間の電子の最大数と最小数を指定する。例えば、強く占有された部分空間からは1電子または2電子の励起のみに制限（最小数を最大数より2だけ小さく設定）する、部分空間内の電子数を最大2個に制限する、などの指定方法がありうる。

==脚注==
{{reflist}}

==参考文献==
* {{cite book
  | last = Cramer
  | first = Christopher J.
  | title = Essentials of Computational Chemistry
  | publisher = John Wiley and Sons
  | year = 2002
  | location = Chichester
  | isbn = 0-471-48552-7
}}

== 関連項目 ==

* [[ハートリー＝フォック方程式]]
* [[ダグラス・ハートリー]]
* [[ウラジミール・フォック]]
* [[量子化学および固体物理計算ソフトの一覧]]

{{DEFAULTSORT:たはいちしこむとうちゃくは}}
[[Category:量子化学]]