アインシュタイン模型

テンプレート:統計力学 アインシュタイン模型（アインシュタインもけい、テンプレート:Lang-en-short）とは、固体の比熱 テンプレート:Mvar の温度依存性を説明するために、アインシュタインが提唱した固体の格子振動に関する模型のことである。N個の同種原子からなる結晶の格子振動を、N個の独立な３次元調和振動子とみなし、且つ全てが同じ角振動数 テンプレート:Mvar を持つとした。

アインシュタインは、1906年に執筆した論文『輻射に関するプランクの理論と比熱の理論』および1910年に執筆した論文『一原子分子からなる固体における弾性的性質と比熱の関係』でこの理論を発表した。

歴史的経緯

1907年にアインシュタインが提唱したこの理論は、歴史的に見ても大きな意味を持つ。経験則であるデュロン-プティの法則で予測される固体の比熱は、古典力学では温度に依存しないはずであった。しかし、低温での実験では、絶対零度で熱容量がゼロに変化することがわかった。温度が上がると比熱は上昇し、高温になるとデュロン-プティの法則に近づく。アインシュタインの理論は、プランクの量子仮説を用いることで、実験的に観測された比熱の温度依存性を初めて説明することに成功した。これは光電効果とともに、量子化の必要性を示す最も重要な証拠のひとつとなった。

格子比熱

格子振動のエネルギーが関与する比熱を格子比熱（テンプレート:Lang-en-short）と呼ぶ。格子振動を量子化したフォノンのエネルギーは、フォノンの波数ベクトルを テンプレート:Mvar、フォノンの振動モードを テンプレート:Mvar として、

$$E=\sum _{𝒌,\alpha }\hslash {\omega }_{𝒌,\alpha }(n_{𝒌,\alpha }+\frac{1}{2})$$

で表される。整数 テンプレート:Mvar で指定される一粒子状態に占有する粒子数の期待値がわかれば、エネルギーの平均値として内部エネルギーが与えられる。結果的に、フォノンの個数の期待値は

$$⟨n_{𝒌,\alpha }⟩=\frac{1}{e^{\hslash {\omega }_{𝒌,\alpha }/k_{B}T}-1}$$

となるので、内部エネルギー テンプレート:Mvar と格子比熱 テンプレート:Mvar はそれぞれ

$$U=\sum _{𝒌,\alpha }\hslash {\omega }_{𝒌,\alpha }(⟨n_{𝒌,\alpha }⟩+\frac{1}{2})=\sum _{𝒌,\alpha }\hslash {\omega }_{𝒌,\alpha }(\frac{1}{e^{\hslash {\omega }_{𝒌,\alpha }/k_{B}T}-1}+\frac{1}{2})$$

$$C_V=\frac{\partial U}{\partial T}=\frac{\partial }{\partial T}\sum _{𝒌,\alpha }\hslash {\omega }_{𝒌,\alpha }(\frac{1}{e^{\hslash {\omega }_{𝒌,\alpha }/k_{B}T}-1}+\frac{1}{2})$$

となる。テンプレート:Math はディラック定数、テンプレート:Mvar はボルツマン定数、テンプレート:Mvar は温度である。括弧内の1/2は温度微分で０となるので、最終的な格子比熱の表式

$$C_V=\frac{\partial }{\partial T}\sum _{𝒌,\alpha }\frac{\hslash {\omega }_{𝒌,\alpha }}{e^{\hslash {\omega }_{𝒌,\alpha }/k_{B}T}-1}$$

が得られる。

アインシュタイン模型

全ての格子振動の角振動数 テンプレート:Mvar が波数 テンプレート:Mvar に依存せず、一定値 テンプレート:Mvar とするモデルがアインシュタイン模型である。このとき、波数ベクトル テンプレート:Mvar の総和は第一ブリュアンゾーン内の逆格子点の数 テンプレート:Mvar になる。また、振動モードは１つの縦モードと２つの横モードの合計３つのモード（３次元に由来）が存在するので、総和は

$$\sum _{𝒌,\alpha }1=3N$$

を与える。つまり、このモデルは角振動数が テンプレート:Mvar の テンプレート:Mvar 個の独立した調和振動子からなる系と見なせる。これを格子比熱の表式に代入すると

$$C_V=\frac{\partial }{\partial T}\left(3N\frac{\hslash {\omega }_E}{e^{\hslash {\omega }_E/k_{B}T}-1}\right)=3N{k}_B{\left(\frac{\hslash {\omega }_E}{k_{B}T}\right)}^2\frac{e^{\hslash {\omega }_E/k_{B}T}}{(e^{\hslash {\omega }_E/k_{B}T}-1{)}^2}$$

となり、角振動数 テンプレート:Mvar を温度に換算した

$${\theta }_E=\frac{\hslash {\omega }_E}{k_B}$$

を用いると

$$C_V=3N{k}_B{\left(\frac{{\theta }_E}{T}\right)}^2\frac{e^{{\theta }_E/T}}{(e^{{\theta }_E/T}-1{)}^2}$$

と表される。この式は一般的に、「アインシュタイン模型」あるいは「アインシュタインの比熱式」と呼ばれ、テンプレート:Mvar はアインシュタイン温度と呼ばれる。

実験事実の再現と乖離

高温領域（$T\gg {\theta }_E$）では、比熱が $C_V\approx 3N{k}_B$ となり、デュロン-プティの法則に漸近することがわかる。また、低温領域（$T\ll {\theta }_E$）では指数関数的にゼロに近づく。この低温での指数関数的な温度依存性は、フォノンが確率 $e^{-\hslash {\omega }_E/k_{B}T}$ で励起することに起因している。この２点の実験事実に理論的な説明を与えたのが、アインシュタイン模型のとても大きな成果である。

しかし、アインシュタイン模型は低温領域では実験事実を正しく再現しない。上述の通り、低温領域でアインシュタイン模型は指数関数的な温度依存性を示すが、実際の実験では冪関数的な温度依存性を見せ、$T^3$ でゼロに近づく。この原因はアインシュタイン模型が低い角振動数の存在を大胆に無視したからである。フォノンの分散関係は角振動数の低い領域で波数に比例する。この点を修正したモデルはデバイ模型と呼ばれ、ピーター・デバイによって1912年に理論的根拠が与えられた。

参考文献

原論文

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

訳書

• テンプレート:Cite book

外部リンク

• テンプレート:Cite web
• テンプレート:Kotobank