箱の中の気体のソースを表示

本項では、量子力学における'''箱の中の量子的な理想気体'''について述べる。すなわち、容器に多数の分子が入っており、[[熱力学的平衡#統計力学|熱化]]のプロセスで一瞬に行われる衝突を除けば、分子どうしの相互作用を行わない系である。この系の平衡状態における性質を調べるには、無限の深さの[[井戸型ポテンシャル]]に置かれた量子的粒子についての結果を用いることができる。

この単純なモデルは、質量をもつ理想[[フェルミ気体]]や、質量を持つ理想[[ボース気体]]、質量をもたないボース気体として扱うことが可能な[[黒体放射]]などの様々な量子理想気体だけでなく、古典的な[[理想気体]]も記述することができる。黒体放射における熱化は、フォトンおよび熱平衡状態にある物体との間の相互作用により促進されると仮定される。

[[マクスウェル＝ボルツマン統計]]または[[ボース＝アインシュタイン統計]]または[[フェルミ＝ディラック統計]]の結果を用い、箱の大きさが無限大だとすると、[[トーマス＝フェルミ近似]]によりエネルギー状態の縮退度は微分として、状態の総和は積分として表現される。
これにより気体の熱力学的な性質は[[分配関数]]や[[グランドカノニカル分配関数]]を用いて計算できる。
ここではいくつかの簡単な例を示す。

==状態の縮退におけるトーマス＝フェルミ近似==
{{Main|トーマス＝フェルミ模型}}
無限の深さを持つ3次元[[井戸型ポテンシャル]]では、粒子に質量がある場合とない場合どちらにおいても、量子数の組[''n''<sub>''x''</sub>,&nbsp;''n''<sub>''y''</sub>,&nbsp;''n''<sub>''z''</sub>]によって粒子の状態の一覧表を作ることができる。
運動量の大きさは次のように与えられる。
:<math>p=\frac{h}{2L}\sqrt{n_x^2+n_y^2+n_z^2} \qquad \qquad n_x,n_y,n_z=1,2,3,\ldots </math>
ここで''h ''は[[プランク定数]]、''L ''は箱の1辺の長さである。
粒子の可能な状態それぞれは自然数の３次元格子の点として考えることができる。
原点から任意の点までの距離は、次のように書ける。
:<math>n=\sqrt{n_x^2+n_y^2+n_z^2}=\frac{2Lp}{h}</math>

それぞれの量子数の組は''f&nbsp;''個の状態を与えるとする。
ここで''f&nbsp;''は粒子の内部自由度で、衝突によって変化する。
たとえば、スピン1/2の粒子では''f=2''で、上向きと下向きそれぞれのスピン状態について1個の状態を数える。
''n ''が大きい場合、運動量の大きさが''p ''以下の状態数は、近似的に
:<math>
g=\left(\frac{f}{8}\right)  \frac{4}{3}\pi n^3 
 = \frac{4\pi f}{3} \left(\frac{Lp}{h}\right)^3
</math>
これはちょうど半径''n&nbsp;''の球の体積の''f&nbsp;''倍を８で割ったものである。なぜなら正の''n<sub>i&nbsp;</sub>''を持つ球の1/8のみを考慮したからである。
連続体近似を用いると、''p&nbsp;''から''p+dp&nbsp;''の運動量を持つ状態の数は、
:<math>
dg=\frac{\pi}{2}~f n^2\,dn =  \frac{4\pi fV}{h^3}~ p^2\,dp
</math>
ここで''V=L<sup>3</sup>&nbsp;''は箱の体積である。
このような連続体近似を用いると、''n<sub>i&nbsp;</sub>=1''の基底状態を含む低エネルギー状態の特徴を描写できなくなることに注意しなければならない。
このことは多くの場合では問題にはならないが、ボース＝アインシュタイン凝縮を考える際は、気体の大半が基底状態または基底状態付近にあり、低エネルギー状態を扱えるかどうかが重要となる。

連続体近似を用いないと、エネルギーε<sub>i </sub>の粒子の数は、次のように与えられる。
:<math>  N_i = \frac{g_i}{\Phi(\epsilon_i)}</math>
ここで
:{|
|-
|<math> g_i</math>,  &nbsp; 状態''i''の[[縮退度]]
|-
|&nbsp;
|-
|<math> \Phi(\epsilon_i) = 
\begin{cases} 
  e^{\beta(\epsilon_i-\mu)},  & \mbox{マ  ク  ス  ウ  ェ  ル  ＝  ボ  ル  ツ  マ  ン  統  計  に  従  う  粒  子  } \\
  e^{\beta(\epsilon_i-\mu)}-1, & \mbox{ボ  ー  ス  ＝  ア  イ  ン  シ  ュ  タ  イ  ン  統  計  に  従  う  粒  子  }\\ 
  e^{\beta(\epsilon_i-\mu)}+1, & \mbox{フ  ェ  ル  ミ  ＝  デ  ィ  ラ  ッ  ク  統  計  に  従  う  粒  子  }\\
\end{cases}</math>
|-
|ここでβ = ''1/kT&nbsp;'',  [[ボルツマン定数]] ''k'', [[温度]] ''T'',  [[化学ポテンシャル]] ''μ'' . 
|-
|([[マクスウェル＝ボルツマン統計]], [[ボース＝アインシュタイン統計]], [[フェルミ＝ディラック統計]]を参照)
|}

連続体近似を用いると、''E''&nbsp; から''E+dE''&nbsp;のエネルギーを持つ粒子数''dN<sub>E</sub>''は、
:<math>dN_E= \frac{dg_E}{\Phi(E)} </math>
:ここで<math>dg_E</math>&nbsp;は''E''&nbsp; から''E+dE''&nbsp;のエネルギーを持つ状態数である。

==エネルギー分布==
前項から導出された結果を用いると、箱の中の気体におけるいくつかの分布を決定できる。
粒子系において、変数<math>A</math>の分布<math>P_A</math>は、<math>A</math>から<math>A+dA</math>の値をもつ粒子の割合を表す<math>P_AdA</math>から定義される。

:<math>P_A~dA = \frac{dN_A}{N} = \frac{dg_A}{N\Phi_A}</math>
ここで
:<math>dN_A</math>, &nbsp;<math>A</math>から<math>A+dA</math>の値を持つ粒子数
:<math>dg_A</math>, &nbsp;<math>A</math>から<math>A+dA</math>の値を持つ状態数
:<math>\Phi_A^{-1}</math>,  &nbsp;<math>A</math>の値を持つ状態が粒子に占有されている確率
: <math>N</math>, &nbsp;&nbsp; &nbsp; 全粒子数
これは次を満たす。
:<math>\int_A P_A~dA = 1</math>
運動量分布<math>P_p</math>において、<math>p</math> から<math>p+dp</math>の運動量をもつ粒子の割合は、
:<math>P_p~dp = \frac{Vf}{N}~\frac{4\pi}{h^3\Phi_p}~p^2dp</math>
またエネルギー分布<math>P_E</math>において、<math>E</math> から<math>E+dE</math>のエネルギーを持つ粒子の割合は、
:<math>P_E~dE = P_p\frac{dp}{dE}~dE</math>
箱の中の粒子（と自由粒子）において、エネルギー <math>E</math>と運動量 <math>p</math>との関係は、質量がある粒子とない粒子では異なっている。
質量のある粒子では、
:<math> E=\frac{p^2}{2m}</math>
質量のない粒子では、
:<math>E=pc\,</math>
ここで<math>m</math>は粒子の質量、<math>c</math>は光速である。これらの関係を用いると、
*  質量のある粒子では
:<math>\begin{alignat}{2}
 dg_E & = \quad \ \left(\frac{Vf}{\Lambda^3}\right)
\frac{2}{\sqrt{\pi}}~\beta^{3/2}E^{1/2}~dE \\
 P_E~dE & = \frac{1}{N}\left(\frac{Vf}{\Lambda^3}\right)
\frac{2}{\sqrt{\pi}}~\frac{\beta^{3/2}E^{1/2}}{\Phi(E)}~dE \\
\end{alignat}
</math>
ここでΛは気体の[[熱的波長]]である。
:<math>\Lambda =\sqrt{\frac{h^2 \beta }{2\pi m}}</math>
これは重要な量である。なぜならΛが粒子間距離 <math>(V/N)</math><sup>''1/3''</sup>のオーダーのときは、量子的な効果が支配し始め、気体はマクスウェル＝ボルツマン気体とは見なせなくなるからである。

*  質量のない粒子では
:<math>\begin{alignat}{2}
 dg_E & = \quad \ \left(\frac{Vf}{\Lambda^3}\right)
\frac{1}{2}~\beta^3E^2~dE \\
 P_E~dE & = \frac{1}{N}\left(\frac{Vf}{\Lambda^3}\right)
\frac{1}{2}~\frac{\beta^3E^2}{\Phi(E)}~dE \\
\end{alignat}
</math>
ここでのΛは質量のない粒子の熱的波長である。
:<math>\Lambda = \frac{ch\beta}{2\,\pi^{1/3}}</math>

==具体例==
以下ではいくつかの具体例を示す。

===質量のあるマクスウェル＝ボルツマン粒子===
この場合は、
:<math>\Phi(E)=e^{\beta(E-\mu)}</math>
エネルギー分布関数を積分し''N''について解くと、
:<math>N = \left(\frac{Vf}{\Lambda^3}\right)\,\,e^{\beta\mu}</math>
これらを元々のエネルギー分布関数に代入すると、
:<math>P_E~dE = 2 \sqrt{\frac{\beta^3 E}{\pi}}~e^{-\beta E}~dE</math>
これは古典的な[[マクスウェル=ボルツマン分布]]から得られる結果と同じである。
その他の結果は[[理想気体]]を参照。

===質量のあるボース＝アインシュタイン粒子===
この場合は、
:<math>\Phi(E)=\frac{e^{\beta E}}{z}-1\,</math>
:ここで&nbsp;&nbsp;<math> z=e^{\beta\mu}.\,</math>
エネルギー分布関数を積分し''N''について解くと[[粒子数]]が得られる。
:<math>N = \left(\frac{Vf}{\Lambda^3}\right)\textrm{Li}_{3/2}(z)</math>
ここでLi<sub>s</sub>(z) は[[多重対数関数]]、Λは[[熱的波長]]である。
多重対数関数の項は正かつ実でなければならず、これは''z&nbsp;''が０から１に増加すると多重対数関数の項は０からζ(3/2)になることを意味する。
温度が０まで下がっていくとΛは増加していき、最終的には''z=1''で臨界値Λ<sub>c </sub>に行き着き、次のようになる。
:<math>N = \left(\frac{Vf}{\Lambda_c^3}\right)\zeta(3/2)</math>
ここで<math>\zeta</math> は[[リーマンゼータ関数]]を表す。
Λ=Λ<sub>c</sub>での温度は臨界温度である。
臨界温度以下では、上記の粒子数についての方程式は解を持たない。
臨界温度はボース＝アインシュタイン凝縮が起こり始める温度である。
上述の通り、連続体近似の問題点は基底状態が無視されていることである。
しかし上記の粒子数についての方程式は、励起状態のボース粒子の数はかなりうまく表現していることがわかり、よって
:<math>
N=\frac{g_0
z}{1-z}+\left(\frac{Vf}{\Lambda^3}\right)\textrm{Li}_{3/2}(z)
</math>
ここで付け加えられた項は基底状態の粒子数である（基底状態エネルギーは無視されていた）。
この方程式は温度０まで成立する。
その他の結果は[[ボース気体]]を参照。

===質量のないボース＝アインシュタイン粒子（黒体放射など）===
質量のない粒子では、前述した質量のない粒子についてのエネルギー分布関数を用いなければならない。
この関数を振動数分布関数に変換すると便利である。
:<math>
P_\nu~d\nu = \frac{h^3}{N}\left(\frac{Vf}{\Lambda^3}\right)
\frac{1}{2}~\frac{\beta^3\nu^2}{e^{(h\nu-\mu)/kT}-1}~d\nu
</math>
ここでΛは質量のない粒子の熱的波長である。
このときエネルギースペクトル密度（単位体積あたり単位振動数あたりのエネルギー）は、
:<math>U_\nu~d\nu = \left(\frac{N\,h\nu}{V}\right) P_\nu~d\nu =  \frac{4\pi f h\nu^3 }{c^3}~\frac{1}{e^{(h\nu-\mu)/kT}-1}~d\nu</math>
その他の熱力学的パラメータは、質量のある粒子の場合と同じように導出される。
たとえば振動数分布関数を積分し、''N''について解くと粒子数が得られる。
:<math>N=\frac{16\,\pi V}{c^3h^3\beta^3}\,\mathrm{Li}_3\left(e^{\mu/kT}\right)</math>
最も一般的な質量のないボース気体は[[黒体]]における[[光子気体]]である。
黒体の空洞を箱と考えると、光子は壁によって断続的に吸収・再放出される。
この場合、光子の数は保存されない。
[[ボース＝アインシュタイン統計]]の導出において粒子数の制限が取り除かれると、これは化学ポテンシャル(''μ'')が０である状況と実質的に同じである。
さらに光子は２つのスピン状態をもつので、''f''の値は２である。
このときエネルギースペクトル密度は、
:<math>U_\nu~d\nu = \frac{8\pi h\nu^3 }{c^3}~\frac{1}{e^{h\nu/kT}-1}~d\nu </math>
これはまさに黒体放射の[[プランクの法則]]におけるエネルギースペクトル密度である。
この手続きを質量のないマクスウェル＝ボルツマン粒子で実行すると、プランクの法則を高温や低密度で近似した[[ヴィーンの放射法則]]が得られる。

ある状況では、光子を含む反応により光子数が保存される（たとえば[[発光ダイオード]]や「白い」空洞）。
これらのケースでは、光子分布関数は非ゼロ化学ポテンシャルを含んでいる(Hermann 2005)。

その他の質量のないボース気体として、[[熱容量]]における[[デバイ模型]]がある。
このとき箱の中の[[フォノン]]気体を考えるが、フォノンの速度は光速より小さく、箱の各軸で波長に最大値が存在する点でフォトンの場合とは異なる。
これは相空間にわたる積分を無限の範囲まで実行することができないことを意味し、[[多重対数関数]]の代わりに[[デバイ関数]]で表されるようになる。

===質量のあるフェルミ＝ディラック粒子（金属中の電子など）===
この場合、以下が適用される。
:<math>\Phi(E)=e^{\beta(E-\mu)}+1\,</math>
エネルギー分布関数を積分することで、
:<math>N=\left(\frac{Vf}{\Lambda^3}\right)\left[-\textrm{Li}_{3/2}(-z)\right]</math>
ここでもLi<sub>s</sub>(z)は[[多重対数関数]]で、Λは[[ド・ブロイの熱的波長]]である。
そのほかの結果は[[フェルミ気体]]を参照。

== 参考文献 ==
* {{cite journal| last = Herrmann| first = F.|author2=Würfel, P.
 |date=August 2005
 | title = Light with nonzero chemical potential| journal = American Journal of Physics
 | volume = 73| issue = 8| pages = 717–723| doi = 10.1119/1.1904623
 | url = http://scitation.aip.org/journals/doc/AJPIAS-ft/vol_73/iss_8/717_1.html
 | accessdate = 2006-11-20|bibcode = 2005AmJPh..73..717H }}
* {{cite book |last=Huang |first=Kerson |authorlink=
 |title=Statistical Mechanics |year=1967 |publisher=[[John Wiley & Sons]] |location=New York}}
* {{cite book |last=Isihara |first=A.
 |title=Statistical Physics |year=1971 |publisher=Academic Press |location=New York}}
* {{cite book |last=Landau |first=L. D.|author2=E. M. Lifshitz
 |title=Statistical Physics |edition=3rd Edition Part 1 |year=1996
 |publisher=Butterworth-Heinemann |location=Oxford}}
* {{cite journal| last = Yan| first = Zijun| year = 2000
 | title = General thermal wavelength and its applications| journal = Eur. J. Phys.
 | volume = 21| pages = 625–631|doi=10.1088/0143-0807/21/6/314
 | url = http://www.iop.org/EJ/article/0143-0807/21/6/314/ej0614.pdf
 | format = PDF| accessdate = 2006-11-20|bibcode = 2000EJPh...21..625Y| issue = 6 }}
* Vu-Quoc, L., [http://clesm.mae.ufl.edu/wiki.pub/index.php/Configuration_integral_%28statistical_mechanics%29 Configuration integral (statistical mechanics)], 2008. this wiki site is down; see [https://web.archive.org/web/20120428193950/http://clesm.mae.ufl.edu/wiki.pub/index.php/Configuration_integral_%28statistical_mechanics%29 this article in the web archive on 2012 April 28].
<!---
{{cite web |url=http://clesm.mae.ufl.edu/wiki.pub/index.php/Configuration_integral_(statistical_mechanics) |title=Configuration_integral_(statistical_mechanics) |accessdate=2008-10-12 |last=Vu-Quoc |first=Loc  }}
--->
{{デフォルトソート:はこのなかのきたい}}
[[Category:統計力学]]