ボース気体

テンプレート:物性物理学 理想ボース気体（テンプレート:Lang-en-short）とは、古典的な理想気体に類似した量子力学的な相のこと。整数値のスピンをもつボース粒子から構成され、ボース–アインシュタイン統計に従う。ボース粒子の統計力学は、サティエンドラ・ボースが光子において開拓した。アルベルト・アインシュタインは質量を持つ粒子に対してボース統計を拡張するとともに、ボース粒子の理想気体が十分に低温で凝縮し、古典的な理想気体とは挙動が異なることを示した。この凝縮はボース＝アインシュタイン凝縮と呼ばれる。

トーマス＝フェルミ近似

テンプレート:See also 理想ボース気体の熱力学は、グランドカノニカル分布によって計算される。ボース気体のグランドカノニカル分布関数は次のように与えられる。

$𝒵(z,\beta ,V)=\prod _i{(1-z{e}^{-\beta {ϵ}_i})}^{-g_i}$

この積のそれぞれの項は、固有のエネルギー テンプレート:Math に相当する。テンプレート:Math はエネルギー テンプレート:Math を持つ状態の数、テンプレート:Mvar は絶対活量(またはフガシティー)で、化学ポテンシャル テンプレート:Mvar を用いて次のように定義される。

$z(\beta ,\mu )=e^{\beta \mu }$

テンプレート:Mvar は次のように定義される。

$\beta =\frac{1}{kT}$

ここで テンプレート:Mvar はボルツマン定数、テンプレート:Mvar は温度である。 全ての熱力学的な量はグランドカノニカル分布関数から導出されるため、全ての熱力学的な量は３つの変数 テンプレート:Mvar、テンプレート:Mvar（または テンプレート:Mvar ）、テンプレート:Mvar のみの関数として考えることができる。 全ての偏微分係数は、３つの変数のうち１つを変数とし、残りの２つは定数とすることで求められる。 ここで次のように定義される無次元のグランドポテンシャルを扱うと便利である。

$\mathrm{\Omega }=-\ln (𝒵)=\sum _i{g}_i\ln (1-z{e}^{-\beta {ϵ}_i}).$

平均エネルギーは準位間のエネルギー差と比べて大きいと仮定するトーマス＝フェルミ近似を適用すると、上記の和は積分で置き換えられる。

$\mathrm{\Omega }\approx \underset{0}{\overset{\mathrm{\infty }}{\int }}\ln (1-z{e}^{-\beta E})dg.$

縮退度 テンプレート:Mvar は一般的な公式によって多くの異なる状況を表現する。

$dg=\frac{1}{\mathrm{\Gamma }(\alpha )}\frac{E^{\alpha -1}}{E_c^{\alpha }}dE$

ここで テンプレート:Mvar は定数、テンプレート:Math は臨界エネルギー、テンプレート:Mvar はガンマ関数である。 たとえば箱の中の質量を持つボース気体では テンプレート:Math で、臨界エネルギーは次のように与えられる。

$\frac{1}{(\beta E_c{)}^{\alpha }}=\frac{Vf}{{\mathrm{\Lambda }}^3}$

ここで テンプレート:Math は熱的波長である。調和トラップ中の質量を持つボース気体では テンプレート:Math で、臨界エネルギーは次のように与えられる。

$\frac{1}{(\beta E_c{)}^{\alpha }}=\frac{f}{(\hslash \omega \beta {)}^3}$

ここで テンプレート:Math は調和ポテンシャルである。テンプレート:Math は体積だけの関数である。

このグランドポテンシャルの方程式は、項別に被積分関数のテイラー級数を積分することにより、 または テンプレート:Math のメリン変換に比例するとすることにより解くことができる。ここで テンプレート:Math は多重対数関数である。解は次のように与えられる。

$\mathrm{\Omega }\approx -\frac{{\text{Li}}_{\alpha +1}(z)}{{\left(\beta E_c\right)}^{\alpha }}$

このボース気体における連続体近似の問題点は、基底状態が実質的に無視されることで、ゼロエネルギーで縮退度がゼロになることである。 この問題点はボース＝アインシュタイン凝縮を扱うときには重大で、次章で扱う。

基底状態の組み入れ

全粒子数はグランドポテンシャルから次のように与えられる。

$N=-z\frac{\partial \mathrm{\Omega }}{\partial z}\approx \frac{{\text{Li}}_{\alpha }(z)}{(\beta E_c{)}^{\alpha }}$

多重対数関数項は実で正でなければならず、最大値は テンプレート:Math のときで テンプレート:Math に等しい。ここで テンプレート:Mvar はリーマンゼータ関数である。 テンプレート:Mvar を固定すると、テンプレート:Mvar の最大値は臨界値 テンプレート:Math で、このとき以下のようになる。

$N=\frac{\zeta (\alpha )}{({\beta }_c{E}_c{)}^{\alpha }}$

これは臨界温度 テンプレート:Math に相当し、これ以下ではトーマス＝フェルミ近似は破綻する。 上記の方程式は臨界温度について解くことができ、次のようになる。

$T_c={\left(\frac{N}{\zeta (\alpha )}\right)}^{1/\alpha }\frac{E_c}{k}$

たとえば テンプレート:Math で、上述の値 テンプレート:Math を用いると、次のようになる。

$T_c={\left(\frac{N}{Vf\zeta (3/2)}\right)}^{2/3}\frac{h^2}{2\pi mk}$

さらに、ここでは臨界温度以下の結果を計算することはできない。なぜなら上記の方程式を用いた粒子数は負になるからである。 ここでの問題点は、トーマス＝フェルミ近似は基底状態の縮退度を０としていることで、これは間違っている。 凝縮を受け入れる基底状態が無いため方程式が破綻する。 しかし上記の方程式は励起状態では粒子数を比較的正確に評価しており、そこへ基底状態を単純に付け加えることは悪い近似ではないことがわかる。

$N=N_0+\frac{{\text{Li}}_{\alpha }(z)}{(\beta E_c{)}^{\alpha }}$

ここで テンプレート:Math は基底状態凝縮の粒子数で、次のように与えられる。

$N_0=\frac{g_{0}z}{1-z}$

図1：正規化された温度 テンプレート:Mvar と様々なボース気体パラメータとの関係。テンプレート:Mvar の値は3/2。実線は テンプレート:Math の場合で、点線は テンプレート:Math の場合。黒線は励起粒子の割合。青線は凝集粒子の割合。赤線は化学ポテンシャルの負数。緑線は テンプレート:Mvar の値。テンプレート:Math と仮定した。

この方程式は絶対零度まで解くことができる。 図１に テンプレート:Math におけるこの方程式の解の結果を示す。これは箱の中のボース気体に相当し、テンプレート:Math とする。 実線は テンプレート:Math の場合、点線は テンプレート:Math の場合を示す。黒線は励起粒子の割合 テンプレート:Math 、青線は凝縮粒子の割合 テンプレート:Math で、赤線は化学ポテンシャル テンプレート:Mvar にマイナス符号をつけたもの、緑線は テンプレート:Mvar の値である。 横軸は次のように定義される正規化された温度 テンプレート:Mvar である。

$\tau =\frac{T}{T_c}$

テンプレート:Mvar や テンプレート:Mvar は低温の極限で テンプレート:Mvar と線形になり、テンプレート:Math は高温の極限で テンプレート:Math と線形になることがうかがえる。 粒子数が増加すると、凝縮の割合と励起の割合は臨界温度で不連続になる。

粒子数の方程式は正規化された温度で表される。

$N=\frac{g_{0}z}{1-z}+N\frac{{\text{Li}}_{\alpha }(z)}{\zeta (\alpha )}{\tau }^{\alpha }$

テンプレート:Mvar と テンプレート:Mvar が与えられると、この方程式は テンプレート:Mvar について解くことができ、テンプレート:Mvar についての級数解は、テンプレート:Mvar のべき乗または テンプレート:Mvar の逆べき乗での漸近展開としての級数の反転の方法によって得ることができる。 これらの展開から、テンプレート:Math 近くでのガスの振る舞いを知ることができる。テンプレート:Mvar が無限大でマクスウェル-ボルツマン分布となる。 特に我々が興味があるのは テンプレート:Mvar が無限大のときで、これは上記の展開から容易に決定する。

熱力学

粒子数に対する方程式に基底状態を加えることは、等価な基底状態の項をグランドポテンシャルに加えることに相当する。

$\mathrm{\Omega }=g_0\ln (1-z)-\frac{{\text{Li}}_{\alpha +1}(z)}{{\left(\beta E_c\right)}^{\alpha }}$

全ての熱力学的な性質は、グランドポテンシャルから計算される。 以下の表では低温と高温の極限、粒子数が無限大の極限での様々な熱力学的な量を示す。 厳密な結果は等号 (=) で表し、テンプレート:Mvar の級数の最初の数項のみの結果は近似記号で表している。

量	一般の場合	$T\ll T_c$	$T\gg T_c$
z		$\approx \frac{\zeta (\alpha )}{{\tau }^{\alpha }}-\frac{{\zeta }^2(\alpha )}{2^{\alpha }{\tau }^{2\alpha }}$	$=1$
凝縮していない粒子の割合 $1-\frac{N_0}{N}$	$=\frac{{\text{Li}}_{\alpha }(z)}{\zeta (\alpha )}{\tau }^{\alpha }$	$={\tau }^{\alpha }$	$=1$
状態方程式 $\frac{PV\beta }{N}=-\frac{\mathrm{\Omega }}{N}$	$=\frac{{\text{Li}}_{\alpha +1}(z)}{\zeta (\alpha )}{\tau }^{\alpha }$	$=\frac{\zeta (\alpha +1)}{\zeta (\alpha )}{\tau }^{\alpha }$	$\approx 1-\frac{\zeta (\alpha )}{2^{\alpha +1}{\tau }^{\alpha }}$
ギブス自由エネルギー $G=\ln (z)$	$=\ln (z)$	$=0$	$\approx \ln \left(\frac{\zeta (\alpha )}{{\tau }^{\alpha }}\right)-\frac{\zeta (\alpha )}{2^{\alpha }{\tau }^{\alpha }}$

全ての量は、高温の極限をとると古典的な理想気体の値に近づいていく。 上記の値を用いて、他の熱力学的な量を計算することができる。 たとえば、内部エネルギーと圧力×体積との関係は、すべての温度にわたって古典的な理想気体と同じである。

$U=\frac{\partial \mathrm{\Omega }}{\partial \beta }=\alpha PV$

定積比熱においても同様である。

$C_v=\frac{\partial U}{\partial T}=k(\alpha +1)U\beta$

エントロピーは次式で与えられる。

$TS=U+PV-G$

高温の極限をとると、次式が得られる。

$TS=(\alpha +1)+\ln \left(\frac{{\tau }^{\alpha }}{\zeta (\alpha )}\right)$

これは テンプレート:Math では単なるザックール・テトローデ方程式を書き換えたものである。 デルタ相互作用をもつ1次元ボース粒子はフェルミ粒子として振る舞い、パウリの排他原理に従う。 デルタ相互作用をもつ1次元ボース粒子はベーテ仮設により厳密に解くことができる。 バルク自由エネルギーと熱力学的ポテンシャルは楊振寧によって計算された。 １次元の場合における相関関数[1]も評価された。 １次元ボース気体は量子的な非線形シュレーディンガー方程式と等価である。

関連項目

• フェルミ気体
• 箱の中の気体
• デバイ模型
• ボース＝アインシュタイン凝縮
• 黒体放射

参考文献

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