ハートレー方程式の導出 のバックアップ(No.1)

更新


ハートレー方程式の導出

ハートレー法の基本方程式となるハートレー方程式を導出する。

変分法を使うので、まだ学んでいなければ次のことだけ理解しておくこと。

  • 変分原理
    • 近似パラメータを含む波動方程式を作ったとする
    • 近似パラメータを調節して真の波動方程式に最も近づけたい
    • それには波動方程式に対するエネルギー期待値を最小化するようにパラメータを調節すれば良い

目次

時間に依らないシュレーディンガー方程式

n 個の同種粒子からなる系を考える。

 &math( \hat H\Phi =\varepsilon \Phi );

ハミルトニアン \hat H は運動エネルギー T 、1体ポテンシャル V_\mathrm{1体} 、 2体ポテンシャル V_\mathrm{2体} の和として表わせる。

 &math( H&=T+V_\mathrm{1体}+V_\mathrm{2体}\\ &=\sum_i T_i + \sum_i V_1(\bm r_i,s_i) + \sum_i \sum_{j>i} V_2(\bm r_i,s_i,\bm r_j,s_j)\\ &=\sum_i T_i + \sum_i V_1(x_i) + \sum_i \sum_{j>i} V_2(x_i,x_j)\\ );*1 j>i となっているのは、同じ2つの電子に対してポテンシャルを重複して計算することのないようにするため

ここで、 \bm r_i 等は空間座標、 s_i などはスピン座標、 x_i などは空間座標とスピン座標を合わせた座標の意味で使っている。

多粒子波動関数モデル

正規化された n 個の1粒子関数 \{\phi_i\} を想定する。 繰り返しになるが x は空間座標 \bm r とスピン座標 s とを合わせた座標である。

 &math( \int d\bm r\,\sum_s\,\phi_i(\bm r,s)\phi_i(\bm r,s)= \int dx\,\phi_i(x)\phi_i(x)=1 );

多粒子波動関数を \{\phi_i\} から作られる単一のハートレー積で表すものとする。

 &math( \Phi=\phi_{p_1}(x_1) \phi_{p_2}(x_2) \cdots\phi_{p_n}(x_n) );

以下で見るとおり、この形に置くこと自体が平均場近似(電子相関および交換相互作用の無視) を仮定していることと同義となる。

エネルギーの期待値

変分法で1粒子波動関数を最適化するため、まずはエネルギーの表式を求めておく。 (変分原理によればエネルギーを最小化する関数が最良の関数である)

 &math( E=\langle H\rangle=\int d^nx\ \Phi^* H \Phi );

以下、各項毎に見ていく。

運動エネルギー

 &math( \langle T_i\rangle &=\int d^nx\ \Phi^* K_i \Phi\\ );

 &math( \phantom{\langle T_i\rangle} &=-\frac{\hbar^2}{2m} \int d^nx\, \phi_1^*(x_1) \phi_2^*(x_2) \cdots\phi_n^*(x_n) \nabla_i^2 \phi_1(x_1) \phi_2(x_2) \cdots\phi_n(x_n)\\ );

 &math( \phantom{\langle T_i\rangle} &=-\frac{\hbar^2}{2m}\int dx_i\, \phi_i^*(x_i) \nabla_i^2 \phi_i(x_i)\\ );

\nabla_i が作用するのは x_i のみなので、 j\ne i については積分が実行できて 1 が現れる。

1体エネルギー

計算は上とほぼ同様の変形により、

 &math( \langle V_i\rangle &=\int d^nx\ \Phi^* V_i \Phi\\ &=\int d^nx\, \phi_1^*(x_1) \phi_2^*(x_2) \cdots\phi_n^*(x_n)\,V_1(x_i)\,\phi_1(x_1) \phi_2(x_2) \cdots\phi_n(x_n)\\ &=\int dx_i\ \phi_i^*(x_i)\,V_1(x_i)\,\phi_i(x_i)\\ );

2体エネルギー

 &math( \langle V_{ij}\rangle &=\int d^nx\ \Phi^* V_{ij} \Phi\\ &=\int d^nx\, \phi_1^*(x_1) \phi_2^*(x_2) \cdots\phi_n^*(x_n)\,V_2(x_i,x_j)\,\phi_1(x_1) \phi_2(x_2) \cdots\phi_n(x_n)\\ &=\int dx\int dx'\, \phi_i^*(x) \phi_j^*(x') \,V_2(x,x')\,\phi_i(x) \phi_j(x') \\ );

エネルギーの最小化

エネルギーの期待値は次のようになった。

 &math( E=&-\frac{\hbar^2}{2m}\sum_i \int dx\ \phi_i^*(x) \nabla^2 \phi_i(x)

  1. \sum_i \int dx\ \phi_i^*(x) V_1(x) \phi_i(x)\\ &+\sum_i\sum_{j\ne i} \iint dxdx'\phi_i^*(x)\phi_j^*(x') V_2(x,x') \phi_i(x)\phi_j(x') );

このエネルギーを最小化するような \set{\phi_i} を求めることにより、 1つのハートレー積で表現可能な波動関数の最良解を探そう。

ただし、 {\phi_i} が規格化されていることを前提としているので、

 &math( \int dx_i \phi_i^*(x)\phi_i(x)=1 );

の条件下で \phi_i を変化させて E を最小化することになる。

そこで ラグランジュの未定係数法 を使う。

正規性を表す条件式は n 個あるので、 n 個の未定係数を 2\varepsilon_{i} として、

 &math( L=E-\sum_i 2\varepsilon_i \Big[\int dx\,\phi_i^*(x)\phi_i(x)-1\Big] );

を定義し、この L \phi_i,\varepsilon_i で微分しゼロと置く。

\varepsilon_i で微分した結果をゼロと置けば正規直交条件が出てくるので、これは \{\phi_i\} として正規化された関数を用いることのみで成立する。

一方、 \phi_i(x) を変化させ、 \phi_i(x)+\delta\phi_i(x) とした時の変化を E\to E+\delta E として、 \delta E/\delta \phi_i=0 となる条件が、求める1体方程式となる。

ここで、ある演算子 \hat H がエルミートであるとき、

 &math( \langle x|\hat H|y\rangle=\langle x|\hat Hy\rangle=\langle \hat Hy|x\rangle^*=\langle y|\hat H^\dagger|x\rangle^*=\langle y|\hat H|x\rangle^* );

すなわち、

 &math( \int dx f^*(x)\hat Hg(x)+\!\int dx g^*(x)\hat Hf(x)=\int dx f^*(x)\hat Hg(x)+\Big(\int dx f^*(x)\hat Hg(x)\Big)^*=2\,\mathrm{Real}\!\int dx f^*(x)\hat Hg(x) );

のようにまとめられる。

そこで各演算子がエルミートであることを使うと、

 &math( \delta L= &-\frac{\hbar^2}{2m}\sum_i \int dx\ \delta\phi_i^*(x) \nabla^2 \phi_i(x)

  • \frac{\hbar^2}{2m}\sum_i \int dx\ \phi_i^*(x) \nabla^2 \delta\phi_i(x)\\ &+\sum_i \int dx\ \delta\phi_i^*(x) V_1(x) \phi_i(x)
  1. \sum_i \int dx\ \phi_i^*(x) V_1(x) \delta\phi_i(x)\\ &+\sum_i\sum_{j\ne i} \int dx \delta\phi_i^*(x)\Big[\int dx'V_2(x,x') |\phi_j(x')|^2\Big]\phi_i(x')
  2. \sum_i\sum_{j\ne i} \int dx \phi_i^*(x)\Big[\int dx'V_2(x,x') |\phi_j(x')|^2\Big]\delta\phi_i(x')\\ &-\sum_i 2\varepsilon_i \int dx\,\delta\phi_i^*(x)\phi_i(x)
  • \sum_i 2\varepsilon_i \int dx\,\phi_i^*(x)\delta\phi_i(x) );

 &math( \phantom{\delta L} =&\,2\,\mathrm{Real}\,\Bigg\{\sum_i \int dx\ \delta\phi_i^*(x) \bigg[-\frac{\hbar^2}{2m}\nabla^2 +V_1(x)+\sum_{j\ne i}\int dx'V_2(x,x') |\phi_j(x')|^2-\varepsilon_i\bigg]\phi_i(x)\Bigg\}\\ =&\,0 );

これが任意の \delta\phi_i(x) に対して成立するためには、 \Big[\ \ \Big]\phi_i(x) がゼロでなければならない。すなわち、

 &math( \underbrace{\Big[-\frac{\hbar^2}{2m}\nabla^2 +V_1(x)+\sum_{j\ne i}\int dx'V_2(x,x') |\phi_j(x')|^2\Big]}_{\hat H_h}\phi_i(x)=\varepsilon_i\phi_i(x) );

これが、 \phi_i(x) を求めるための1粒子方程式となる。

多粒子波動関数を単一のハートレー積の形に限定することにより、 非常に自然な流れで「平均場近似」が得られることに注目せよ。

多粒子エネルギーと1粒子エネルギー


*1 j>i となっているのは、同じ2つの電子に対してポテンシャルを重複して計算することのないようにするため

Counter: 3343 (from 2010/06/03), today: 1, yesterday: 2