トーマス＝フェルミ模型

トーマス＝フェルミ模型（トーマス＝フェルミもけい、テンプレート:Lang-en-short）とは、シュレーディンガー方程式が導入されて間もなく、それを半古典的に扱った多体系の電子構造についての量子力学的な理論のことである。テンプレート:Ill2とエンリコ・フェルミに因んで名づけられた。波動関数から離れて電子密度を用いて定式化したもので、密度汎関数理論の原型ともなった。トーマス＝フェルミ模型は、核電荷が無限大の極限においてのみ正確な結果を与える。現実的な系を考えるために近似を用いると、定量性に乏しい予言しかできず、原子の殻構造や固体のフリーデル振動のような密度についてのいくつかの一般的性質を再現することもできなくなる。しかし定性的な傾向を解析的に抽出でき、またモデルを解くことが簡単であることから、多くの分野で応用されている。トーマス＝フェルミ理論により表現された運動エネルギーは、オービタルフリー密度汎関数理論のようなより洗練された密度近似運動エネルギーの一つとしても使われている。

1927年にトーマスとフェルミは独立に、この統計的モデルを用いて原子中の電子分布を近似した。実際の電子は原子中で不均一に分布しているが、近似的に電子は微小体積要素 テンプレート:Math に（局所的に）それぞれ均一に分布しており、電子密度 テンプレート:Math は各 テンプレート:Math で異なっているとする。

運動エネルギー

基底状態にある原子中の微小体積要素 テンプレート:Math において、フェルミ球の体積 テンプレート:Math はフェルミ運動量を テンプレート:Math とすると以下のように書けるテンプレート:Sfn。

$V_{\mathrm{F}}=\frac{4}{3}\pi p_{\mathrm{F}}^3(𝒓)$

ここで テンプレート:Mvar は テンプレート:Math 中の点を表す。

したがって、この空間領域に対応する相空間上の領域の体積は次のように書ける。

$\mathrm{\Delta }{V}_{\mathrm{p}\mathrm{h}}=V_{\mathrm{F}}\mathrm{\Delta }V=\frac{4}{3}\pi p_{\mathrm{F}}^3(𝒓)\mathrm{\Delta }V$

テンプレート:Math 中の電子は均一に分布しており、この相空間での体積 テンプレート:Math（テンプレート:Mvar はプランク定数）あたり2つの電子を持つテンプレート:Sfn。テンプレート:Math 中の電子数は、

$\mathrm{\Delta }{N}_{\mathrm{p}\mathrm{h}}=\frac{2}{h^3}\mathrm{\Delta }{V}_{\mathrm{p}\mathrm{h}}=\frac{8\pi }{3h^3}{p}_{\mathrm{F}}^3(𝒓)\mathrm{\Delta }V$

一方で、テンプレート:Math 中の電子数は次のように表わされる。

$\mathrm{\Delta }N=n(𝒓)\mathrm{\Delta }V$

ここで テンプレート:Math は電子密度である。

テンプレート:Math 内の電子数と テンプレート:Math 内の電子数は等しいから、

$n(𝒓)=\frac{8\pi }{3h^3}{p}_{\mathrm{F}}^3(𝒓)$

位置 テンプレート:Mvar における テンプレート:Mvar から テンプレート:Mathの運動量をもつ電子の存在確率は、絶対零度におけるフェルミ分布から以下のように書ける。

$F_{𝒓}(p)\mathrm{d}p=\{\begin{array}{cc}\frac{4\pi p^2}{(4/3)\pi p_{\mathrm{F}}^3(𝒓)}\mathrm{d}p& p\le p_{\mathrm{F}}(𝒓)\\ 0& \text{otherwise}\end{array}$

電子を質量 テンプレート:Math の質点とみなして古典的に運動エネルギーを計算することにすると、原子中の電子の位置 テンプレート:Mvar における単位体積あたりの運動エネルギーは、

$\begin{array}{cc}t(𝒓)& =\int \frac{p^2}{2m_{\mathrm{e}}}n(𝒓)F_{𝒓}(p)\mathrm{d}p\\ & =n(𝒓)\underset{0}{\overset{p_{\mathrm{F}}(𝒓)}{\int }}\frac{p^2}{2m_e}\frac{4\pi p^2}{(4/3)\pi p_{\mathrm{F}}^3(𝒓)}\mathrm{d}p\\ & =C_{\mathrm{F}}[n(𝒓){]}^{5/3}\end{array}$

ここで テンプレート:Math は先ほどの テンプレート:Math で表したものであり、テンプレート:Math は、

$C_{\mathrm{F}}=\frac{3h^2}{10m_{\mathrm{e}}}{\left(\frac{3}{8\pi }\right)}^{2/3}=\frac{3{\hslash }^2}{10m_{\mathrm{e}}}{\left(3{\pi }^2\right)}^{2/3}$

単位体積あたりの運動エネルギー テンプレート:Math を全空間で積分すると、電子の全運動エネルギーが得られるテンプレート:Sfn。

$T=C_{\mathrm{F}}\int [n(𝒓){]}^{5/3}{\mathrm{d}}^{3}r$

よってトーマス＝フェルミ模型によって、電子の全エネルギーは空間的に変化する電子密度 テンプレート:Math のみで表せることが示された。この電子の運動エネルギー表現と、原子核-電子相互作用と電子-電子相互作用の古典的な表現(どちらも電子密度で表せる)とを合わせることで、原子のエネルギーを計算できる。

ポテンシャルエネルギー

原子中の電子のポテンシャルエネルギーは、古典的には正電荷である原子核のと電子との間のクーロン引力によるので、以下のように書ける。

$U_{\mathrm{e}\mathrm{N}}=\int n(𝒓)V_{\mathrm{N}}(𝒓){\mathrm{d}}^{3}r$

ここで テンプレート:Math は位置 テンプレート:Mvar における電子のポテンシャルエネルギーで、原子核の電場によるものである。位置 テンプレート:Math を中心とする電荷 テンプレート:Mvar（テンプレート:Mvar は電気素量、テンプレート:Mvar は正整数）の原子核の場合、

$V_{\mathrm{N}}(𝒓)=-\frac{Z{e}^2}{r}$

同様に、電気的な相互反発による電子のポテンシャルエネルギーは以下のように書ける。

$U_{\mathrm{e}\mathrm{e}}=\frac{1}{2}{e}^2\iint \frac{n(𝒓)n({𝒓}^{\prime })}{|𝒓-{𝒓}^{\prime }|}{\mathrm{d}}^{3}r{\mathrm{d}}^3{r}^{\prime }$

全エネルギー

電子の全エネルギーは、運動エネルギーとポテンシャルエネルギーの和であるテンプレート:Sfn。

$\begin{array}{cc}E& =T+U_{\mathrm{e}\mathrm{N}}+U_{\mathrm{e}\mathrm{e}}\\ & =C_{\mathrm{F}}\int [n(𝒓){]}^{5/3}{\mathrm{d}}^{3}r+\int n(𝒓)V_{\mathrm{N}}(𝒓){\mathrm{d}}^{3}r+\frac{1}{2}{e}^2\iint \frac{n(𝒓)n({𝒓}^{\prime })}{|𝒓-{𝒓}^{\prime }|}{\mathrm{d}}^{3}r{\mathrm{d}}^3{r}^{\prime }\end{array}$

誤差と改良

トーマス＝フェルミ方程式による運動エネルギーは近似に過ぎず、またパウリの原理による原子の交換エネルギーも考慮していない。交換エネルギー項は1928年にポール・ディラックによって付け加えられた。

しかしトーマス＝フェルミ＝ディラック理論は依然として不正確である。誤差の大部分は運動エネルギー部分によるもので、その次に大きいのが電子相関を完全に無視したことによる誤差である。1962年にエドワード・テラーはトーマス＝フェルミ理論では分子結合を記述できないことを示した。トーマス＝フェルミ理論で計算した分子のエネルギーは、構成原子のエネルギーの和より高くなる。しかし一般的に結合長が均一に増加し、ばらばらの原子の集まりにより近い状態になると、分子の全エネルギーは減少するテンプレート:Sfn。これは運動エネルギー表現を改良することで克服することができるテンプレート:Sfn。その処方の一つとしてコーン・シャム法が挙げられる。

トーマス＝フェルミの運動エネルギーは、ヴァイツゼッカー相関（1935年）を付け加えることで改良でき、トーマス＝フェルミ＝ディラック＝ヴァイツゼッカー密度汎関数理論(TFDW-DFT)と呼ばれる。

$T_{\mathrm{W}}=\frac{1}{8}\frac{{\hslash }^2}{m}\int \frac{|\mathrm{\nabla }n(\overrightarrow{r}){|}^2}{n(\overrightarrow{r})}dr$

参考文献

• テンプレート:Cite book
• テンプレート:Cite book
• テンプレート:Cite book

関連項目

• 縮退した状態でのトーマス＝フェルミ近似