Post Hartree-Fock 法

    Loading...

## Hartree-Fock 法とスピン

### RHF 法

Restricted Hartree-Fock (RHF) 法閉殻系に適用する HF 法である。

この方法では、全波動関数 ΦHF\varPhi_\mathrm{HF} を、HF 軌道 ϕi()\phi_i(\cdot) とスピン軌道 α(), β()\alpha(\cdot), ~ \beta(\cdot) の積からなる 2N2N 次元の Slater 行列式で表す:

ΦRHF(x1,,x2N)=1(2N)!ϕ1(r1)α(σ1)ϕ1(r1)β(σ1)ϕ1(r2N)α(σ2N)ϕ1(r2N)β(σ2N)ϕN(r1)α(σ1)ϕN(r1)β(σ1)ϕN(r2N)α(σ2N)ϕN(r2N)β(σ2N)\varPhi_\mathrm{RHF}(\bm{x}_1,\ldots,\bm{x}_{2N}) = \frac{1}{\sqrt{(2N)!}} \begin{vmatrix} \phi_1(\bm{r}_1)\alpha(\sigma_1) & \phi_1(\bm{r}_1)\beta(\sigma_1) & \cdots & \phi_1(\bm{r}_{2N})\alpha(\sigma_{2N}) & \phi_1(\bm{r}_{2N})\beta(\sigma_{2N}) \\ \vdots & \vdots & \ddots & \vdots & \vdots \\ \phi_N(\bm{r}_1)\alpha(\sigma_1) & \phi_N(\bm{r}_1)\beta(\sigma_1) & \cdots & \phi_N(\bm{r}_{2N})\alpha(\sigma_{2N}) & \phi_N(\bm{r}_{2N})\beta(\sigma_{2N}) \end{vmatrix}

ここで xi:=(ri,σi)\bm{x}_i := (\bm{r}_i, \sigma_i) は電子 ii の空間座標とスピン座標の組である。

閉殻系ではすべての電子が (α,β)(\alpha, \beta) の対を作るため、各スピン軌道は同数を用意する。

### UHF 法

Unrestricted Hartree-Fock (RHF) 法開殻系に適用する HF 法である。

この方法では、NN 組の α()\alpha(\cdot)N1N-1 組の β()\beta(\cdot) を用意し、それぞれ異なる軌道関数 ϕi()\phi_i(\cdot)ψi()\psi_i(\cdot) と掛け合わせる。
同じ ii に対する ϕi\phi_iψi\psi_i は縮退しておらず、軌道エネルギーはわずかに異なる。

## 配置間相互作用(CI)法

HF 波動関数 ΦHF\varPhi_\mathrm{HF} に対して、励起配置の Slater 行列式を加えた近似を配置間相互作用(CI)法という。

電子 iiaa へ励起した配置を Φia\varPhi_i^a、電子 i,ji,ja,ba,b へ励起した配置を Φi,ja,b\varPhi_{i,j}^{a,b} \cdots とするとき、その配置関数 ΦCI\varPhi_\mathrm{CI} は、

ΦCI=C0ΦHF+iaC1Φia+i>ja,bC2Φi,ja,b+\varPhi_\mathrm{CI} = C_0\varPhi_\mathrm{HF} + \sum_i \sum_a C_1 \varPhi_i^a + \sum_{i > j} \sum_{a, b} C_2 \varPhi_{i,j}^{a,b} + \cdots

となる。ここで CkC_k は配置の重み(CI 係数)であり、これを変分法によって定める:

minimizeE[C0,C1,]=ΦCIH^ΦCI\mathrm{minimize} \quad E[C_0, C_1, \ldots] = \bra{\varPhi_\mathrm{CI}} \hat{H} \ket{\varPhi_\mathrm{CI}}

配置間相互作用法は多くの Slater 行列式を必要とするため、一般に計算コストが高い。

### CISD 法

CISD 法は、HF 配置のほかに2電子励起配置までを考慮した配置間相互作用法である。

電子 iiaa へ励起した配置を Φia\varPhi_i^a、電子 i,ji,ja,ba,b へ励起した配置を Φi,ja,b\varPhi_{i,j}^{a,b} とするとき、その配置関数 ΦCISD\varPhi_\mathrm{CISD} は、

ΦCISD=CHFΦHF+iaCiaΦia+i>ja,bCi,ja,bΦi,ja,b\varPhi_\mathrm{CISD} = C_\mathrm{HF}\varPhi_\mathrm{HF} + \sum_i \sum_a C_i^a \varPhi_i^a + \sum_{i > j} \sum_{a, b} C_{i,j}^{a,b} \varPhi_{i,j}^{a,b}

となる。

### CID 法

CID 法は、HF 配置のほかに2電子励起配置のみを考慮した配置間相互作用法である。

ΦCID=CHFΦHF+i>ja,bCi,ja,bΦi,ja,b\varPhi_\mathrm{CID} = C_\mathrm{HF}\varPhi_\mathrm{HF} + \sum_{i > j} \sum_{a, b} C_{i,j}^{a,b} \varPhi_{i,j}^{a,b}

### Full-CI 法

Full-CI 法はすべての励起配置を考慮した配置間相互作用法である。非常に計算コストが高く、十分小さな系でないと適用できない。

## Møller-Plesset 摂動(MP)法

## クラスター展開(CC)法

Discussions

記事がありません