量子力学Ⅰ/平均場近似 のバックアップ(No.7)

目次 †

平均場近似による１体問題化 †

前述のように、「粒子間の相互作用がないこと」を仮定すれば、 １粒子波動関数を複数掛け合わせて、多粒子シュレーディンガー方程式を解ける。

しかし、「粒子間の相互作用がないこと」を要件とする限り、 興味のある問題にはまったく適当できない。

そこで、「擬似的に相互作用をなくすために」粒子・粒子間の相互作用を平均化して、 ポテンシャル に含めてしまう平均場近似が行われる。

ハートレーの方法 †

上で見た水素様「分子」のポテンシャル

　&math( V(\bm r_1,\bm r_2)&=\frac{e^2}{4\pi\epsilon_0}\Biggl[ \underbrace{ \frac{-1}{|\bm r_1-\bm R_1|}+ \frac{-1}{|\bm r_1-\bm R_2|}}_{電子1と原子核}+ \underbrace{ \frac{-1}{|\bm r_2-\bm R_1|}+ \frac{-1}{|\bm r_2-\bm R_2|}}_{電子2と原子核}+ \underbrace{\frac{+1}{|\bm r_1-\bm r_2|}}_{電子間相互作用} \Biggr]\\ &=V_{１体}(\bm r_1)+V_{１体}(\bm r_2)+V_{２体}(\bm r_1,\bm r_2) );

でもそうだったように、 粒子間の相互作用が「２体相互作用の重ね合わせ」で書けるとすれば、 多体問題のポテンシャルを一般に次のように表せる。

　&math( V(\bm r_1,\bm r_2,\dots,\bm r_n)= \underbrace{\sum_{j=1}^n V_j(\bm r_j)}_{１粒子ポテンシャル}+ \underbrace{\sum_{j=1}^n\sum_{k=j+1}^n V_{j,k}(\bm r_j,\bm r_k)}_{２粒子ポテンシャル} );

（相互作用ポテンシャルを２重にカウントしないため、 ２つ目の の は から始まっている）

このようなポテンシャルを仮定した場合、 例えば粒子 の感じるポテンシャルは他の粒子の位置により変化するのであるが、

　&math( v_j(\bm r_j)&=V_j(\bm r_j)+\sum_{k\ne j} V_{j,k}(\bm r_j,\bm r_k)\\ );

これを "平均場" で置き換えよう。

　&math( \overline{v_j}(\bm r_j)&= V_j(\bm r_j)+\sum_{k\ne j} \overline{ V_{j,k}}(\bm r_j)\\ );

ただし、

　

であり、これはポテンシャル を粒子 が に見いだされる確率 で重み付けして平均したものである。 番目の粒子は様々に動き回るが、その際 に存在する確率が であることから、粒子 の感じるであろう「平均的なポテンシャル」を上記のように求めたわけである。

この に対する「１体のシュレーディンガー方程式」 を解いて、粒子 に対する一体の波動関数 を得る。

　&math( \left[-\frac{\hbar^2}{2m_j}\nabla_{r_j}^2+\overline{v_j}(\bm r_j)\right]\varphi_j(\bm r_j)=\varepsilon_j\varphi_j(\bm r_j) );

このようにして、平均場近似により多体問題が１体問題に変換されたことになる。

ただし、 を求めるのに が必要で、 個々の を求めるのに が必要なので、 この方程式はそのままでは解けない。

始めに適当な を仮定して を求め、 そこから新しい を求め、、、、などと繰り返して、 「全体としてつじつまの合う（セルフコンシステントな＝自己無頓着な）」解 を得るような手順（= 自己無頓着場の方法 あるいは セルフコンシステント法）が必要となる。

そのようにして求めた から

　

を作れば、１体のハミルトニアン

　&math( &\hat h_j=-\frac{\hbar^2}{2m_j}\nabla_{r_j}+V_j(\bm r_j)+\sum_{k\ne j}\overline{V_{jk}}(\bm r_j)\\ );

をにらみつつ、

　&math( &\hat H\Phi(\bm r_1,\bm r_2,\dots,\bm r_n)\\ &=\left[-\sum_j\frac{\hbar^2}{2m_j}\nabla_{r_j}+\sum_jV_j(\bm r_j)+\sum_j\sum_{k>j}V_{jk}(\bm r_j,\bm r_k)\right]\varphi_1(\bm r_1)\varphi_2(\bm r_2)\dots\varphi_n(\bm r_n)\\ &\sim\sum_j\left[-\frac{\hbar^2}{2m_j}\nabla_{r_j}+V_j(\bm r_j)+\sum_{k>j}\overline{V_{jk}}(\bm r_j)\right]\varphi_1(\bm r_1)\varphi_2(\bm r_2)\dots\varphi_n(\bm r_n)\\ &=\sum_j\left[\hat h_j-\sum_{k<j}\overline{V_{jk}}(\bm r_j)\right]\varphi_1(\bm r_1)\varphi_2(\bm r_2)\dots\varphi_n(\bm r_n)\\ &\sim\underbrace{\left[\sum_j\varepsilon_j-\sum_j\sum_{k<j}\overline{V_{jk}}\right]}_{E}\varphi_1(\bm r_1)\varphi_2(\bm r_2)\dots\varphi_n(\bm r_n)\\ &=E\Psi(\bm r_1,\bm r_2,\dots,\bm r_n) );

のように、多粒子ハミルトニアンの(近似的な)固有関数になっていることが分かる。

ここで、

　&math( \overline{V_{jk}}=\iint V_{jk}(\bm r_j,\bm r_k)\,|\varphi_j(\bm r_j)|^2\,|\varphi_k(\bm r_k)|^2\,d\bm r_j\,d\bm r_k );

とした。

このときのエネルギー固有は、

　

であるから、１体問題のエネルギー固有値の単純な足し算にはなっていない。 これは、２粒子ポテンシャルのエネルギーを と とで２回取り込んでしまっているためであり、上式の２つ目のシグマはこの分を差し引いているものである。

上記のようにして求まる解が「多体波動関数を１つのハートレー積で表示する場合」の最適解となることを「変分法」（量子力学II で学ぶ）を用いて確かめることができる。

ハートレー・フォックの方法 †

多体波動関数を「１つのハートレー積で表示する場合」の最適解を求める方法がハートレー法だったのに対して、多体波動関数を「１つのスレイター行列で表示する場合」の最適解を求める方法がハートレー・フォック法である。*1「最適解」の意味は変分法を学ぶところで理解せよ

ハートレー・フォック法ではハートレー法とは異なり、１粒子方程式（「フォック方程式」と呼ばれる）はすべての粒子に対して共通となる。

フォック方程式には、ハートレー法でも現れた

• １粒子ポテンシャル
• 平均化された２粒子ポテンシャル

の他に、

• 波動関数が反対称化されたことを反映した「交換ポテンシャル」

と呼ばれる項が現れる。

この項は、波動関数の反対称性から生まれるパウリの排他律により、 同じスピンを持つ電子が互いに避け合う効果のために、 反対称化しない場合に比べて電子・電子間の斥力エネルギーが低下する影響（交換相互作用）を表している。

パウリの排他律により複数の粒子が同じ固有関数状態を取ることはできないため、 フォック方程式の固有関数をエネルギーの低い方から 個集めて作ったスレイター行列が、求める最適解となる。

この１粒子波動関数は正規直交系に選べるため、理論的にはハートレー法よりも見通しが良くなる。 （実際に数値計算等で求めるのはより難しくなっている）

また、この形では（仮想的な）励起状態やイオン化状態を容易に表せることもメリットとなる。

固有関数をエネルギーの低い方から順に 個集める代りに、 いくつか飛ばしながら 個集めれば、 それは低エネルギー状態の電子を高エネルギー状態へ励起した、 (仮想的な)「励起状態」を表す。

当然、固有関数の取り方を変えればフォック方程式に現れるポテンシャル形状も変化するため、 取り方を変えるたびに自己無撞着計算をやり直さなければならず、 それにともなって固有関数の具体的な形も変化するが、 それでも最終的に「励起状態」に対応する解が求まる。

電子数を増やしたり減らしたりする場合も、対応するポテンシャル項を付けたり、 減らしたりすれば任意の に対して電子 個の解を得られる。

&uml( start repeat

 :平均場ポテンシャル;
 :   フォック方程式   ;
 :正規直交固有関数系;
 :任意のn個を対称化 ;

repeat while (収束まだ？) stop );

ハートレー・フォック法での解法をまとめると上記のようになる。

適当なポテンシャルから始めて、

1. １体方程式の解として、フォック方程式の正規直交固有関数系を得る
2. 任意の数の１体波動方程式を選び、それらを対称化して多体波動関数を作る
   • ここでの選び方により、任意の数の粒子の任意の基底／励起状態を作れる
3. 多体波動関数から平均場ポテンシャルを求めて、
   • 1. で仮定したものと等しければ終了
   • 等しくなければポテンシャルを適当に更新して 1. から再スタート

全体としてセルフコンシステントな解が得られるまで繰り返す

2. において、任意の数の粒子を含む任意の基底／励起状態を作れることが、 さらに進んだ議論において重要となってくる。

質問・コメント †