アインシュタイン模型

統計力学


熱力学 · 気体分子運動論
粒子統計
マクスウェル＝ボルツマン

ボース＝アインシュタイン
フェルミ＝ディラック
パラ · エニオン · 組み紐（英語版）

アンサンブル
ミクロカノニカルアンサンブル

カノニカルアンサンブル
グランドカノニカルアンサンブル
等温定圧アンサンブル
等エンタルピー-定圧

熱力学
気体の法則（英語版） · カルノーサイクル

デュロン＝プティの法則

模型
デバイ · アインシュタイン · イジング
熱力学ポテンシャル
内部エネルギー
エンタルピー
ヘルムホルツの自由エネルギー
ギブズの自由エネルギー
グランドポテンシャル
科学者
マクスウェル · ギブズ · ボルツマン · アインシュタイン · オンサーガー · ウィルソン · 久保亮五 · カダノフ · フィッシャー · 川崎恭治 · パリージ · エドワーズ · ローレンツ · 蔵本由紀 · ジャルジンスキー
  • 表
  • 話
  • 編
  • 歴

アインシュタイン模型（アインシュタインもけい、英: Einstein model）とは、固体の比熱[注 1] CV の温度依存性を説明するために、アインシュタインが提唱した固体の格子振動に関する模型のことである。N個の同種原子からなる結晶の格子振動を、N個の独立な３次元調和振動子とみなし、且つ全てが同じ角振動数 ωE を持つとした。

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

歴史的経緯

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

格子比熱

格子振動のエネルギーが関与する比熱を格子比熱（英: lattice specific heat）と呼ぶ。格子振動を量子化したフォノンのエネルギーは、フォノンの波数ベクトルを k、フォノンの振動モードを α として、

E = ∑ k , α ℏ ω k , α ( n k , α + 1 2 ) {\displaystyle E=\sum _{{\boldsymbol {k}},\alpha }\hbar \omega _{{\boldsymbol {k}},\alpha }\left(n_{{\boldsymbol {k}},\alpha }+{\frac {1}{2}}\right)}

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

⟨ n k , α ⟩ = 1 e ℏ ω k , α / k B T − 1 {\displaystyle \left\langle n_{{\boldsymbol {k}},\alpha }\right\rangle ={\frac {1}{e^{\hbar \omega _{{\boldsymbol {k}},\alpha }/k_{B}T}-1}}}

となるので、内部エネルギー U と格子比熱 CV はそれぞれ

U = ∑ k , α ℏ ω k , α ( ⟨ n k , α ⟩ + 1 2 ) = ∑ k , α ℏ ω k , α ( 1 e ℏ ω k , α / k B T − 1 + 1 2 ) {\displaystyle U=\sum _{{\boldsymbol {k}},\alpha }\hbar \omega _{{\boldsymbol {k}},\alpha }\left(\left\langle n_{{\boldsymbol {k}},\alpha }\right\rangle +{\frac {1}{2}}\right)=\sum _{{\boldsymbol {k}},\alpha }\hbar \omega _{{\boldsymbol {k}},\alpha }\left({\frac {1}{e^{\hbar \omega _{{\boldsymbol {k}},\alpha }/k_{B}T}-1}}+{\frac {1}{2}}\right)}

C V = ∂ U ∂ T = ∂ ∂ T ∑ k , α ℏ ω k , α ( 1 e ℏ ω k , α / k B T − 1 + 1 2 ) {\displaystyle C_{V}={\partial U \over \partial T}={\partial \over \partial T}{\sum _{{\boldsymbol {k}},\alpha }\hbar \omega _{{\boldsymbol {k}},\alpha }\left({\frac {1}{e^{\hbar \omega _{{\boldsymbol {k}},\alpha }/k_{B}T}-1}}+{\frac {1}{2}}\right)}}

となる。ħ はディラック定数、kB はボルツマン定数、T は温度である。括弧内の1/2は温度微分で０となるので、最終的な格子比熱の表式

C V = ∂ ∂ T ∑ k , α ℏ ω k , α e ℏ ω k , α / k B T − 1 {\displaystyle C_{V}={\partial \over \partial T}{\sum _{{\boldsymbol {k}},\alpha }{\frac {\hbar \omega _{{\boldsymbol {k}},\alpha }}{e^{\hbar \omega _{{\boldsymbol {k}},\alpha }/k_{B}T}-1}}}}

が得られる。

アインシュタイン模型

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

∑ k , α 1 = 3 N {\displaystyle \sum _{{\boldsymbol {k}},\alpha }1=3N}

を与える。つまり、このモデルは角振動数が ωE の 3N 個の独立した調和振動子からなる系と見なせる。これを格子比熱の表式に代入すると

C V = ∂ ∂ T ( 3 N ℏ ω E e ℏ ω E / k B T − 1 ) = 3 N k B ( ℏ ω E k B T ) 2 e ℏ ω E / k B T ( e ℏ ω E / k B T − 1 ) 2 {\displaystyle C_{V}={\partial \over \partial T}\left(3N{\frac {\hbar \omega _{E}}{e^{\hbar \omega _{E}/k_{B}T}-1}}\right)=3Nk_{B}\left({\frac {\hbar \omega _{E}}{k_{B}T}}\right)^{2}{\frac {e^{\hbar \omega _{E}/k_{B}T}}{(e^{\hbar \omega _{E}/k_{B}T}-1)^{2}}}}

となり、角振動数 ωE を温度に換算した

θ E = ℏ ω E k B {\displaystyle \theta _{E}={\hbar \omega _{E} \over k_{B}}}

を用いると

C V = 3 N k B ( θ E T ) 2 e θ E / T ( e θ E / T − 1 ) 2 {\displaystyle C_{V}=3Nk_{B}\left({\frac {\theta _{E}}{T}}\right)^{2}{\frac {e^{\theta _{E}/T}}{(e^{\theta _{E}/T}-1)^{2}}}}

と表される。この式は一般的に、「アインシュタイン模型」あるいは「アインシュタインの比熱式」と呼ばれ、θE はアインシュタイン温度と呼ばれる。

実験事実の再現と乖離

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

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

脚注

[脚注の使い方]

注釈

  1. ^ 熱容量は示量性の物理量であり、これを質量やモル数などの示量性の量で割って示強性にしたものが比熱である。本項では、熱容量と呼ぶべきところを、慣習に従って比熱と呼ぶ。
  2. ^ 論文が雑誌に掲載され刊行されたのは1907年である。
  3. ^ 論文が雑誌に掲載され刊行されたのは1911年である。

出典

参考文献

原論文

  • Einstein, A. (November 9, 1906). “Die Plancksche Theorie der Strahlung und die Theorie der spezifischen Wärme [輻射に関するプランクの理論と比熱の理論]” (German) (PDF). Annalen der Physik. Ser. 4 (Weinheim: Wiley-VCH Verlag（ドイツ語版、英語版）) 327 (1): 180–190. Bibcode: 1906AnP...327..180E. doi:10.1002/andp.19063270110. ISSN 0003-3804. LCCN 50-13519. OCLC 5854993. http://www.physik.uni-augsburg.de/annalen/history/einstein-papers/1907_22_180-190.pdf. 
  • Einstein, A. (November 30, 1910). “Eine Beziehung zwischen dem elastischen Verhalten und der Spezifischen Wärme mit einatomigem Molekül [一原子分子からなる固体における弾性的性質と比熱の関係]” (German) (PDF). Annalen der Physik. Ser. 4 (Weinheim: Wiley-VCH Verlag（ドイツ語版、英語版）) 339 (1): 170–174. Bibcode: 1911AnP...339..170E. doi:10.1002/andp.19113390110. ISSN 0003-3804. LCCN 50-13519. OCLC 5854993. http://www.physik.uni-augsburg.de/annalen/history/einstein-papers/1911_34_170-174.pdf. 

訳書

外部リンク

  • Shinkai, Hisaaki (2005年7月9日). “Albert Einsteinの生涯と代表論文 年表”. Einstein 1905「数理科学/宇宙物理」研究の最新情報と「研究者/研究生活」に関する情報サイト. 2016年11月18日閲覧。
  • ブリタニカ国際大百科事典 小項目事典『アインシュタインの比熱式』 - コトバンク
  • 表
  • 話
  • 編
  • 歴
キャリア
著作物（英語版）
家族
アインシュタイン賞
その他
カテゴリ カテゴリ
スタブアイコン

この項目は、物理学に関連した書きかけの項目です。この項目を加筆・訂正などしてくださる協力者を求めています（プロジェクト:物理学／Portal:物理学）。

  • 表示
  • 編集