ボース気体のソースを表示

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

== トーマス＝フェルミ近似 ==
{{see also|トーマス＝フェルミ模型}}
理想ボース気体の熱力学は、[[グランドカノニカル分布]]によって計算される。ボース気体のグランドカノニカル分布関数は次のように与えられる。
:<math>\mathcal{Z}(z,\beta,V) = \prod_i \left(1-ze^{-\beta\epsilon_i}\right)^{-g_i}</math>
この積のそれぞれの項は、固有のエネルギー {{math|''&epsilon;<sub>i</sub>''}} に相当する。{{math|''g<sub>i</sub>''}} はエネルギー {{math|''&epsilon;<sub>i</sub>''}} を持つ状態の数、{{mvar|z}} は絶対活量(または[[フガシティー]])で、[[化学ポテンシャル]]  {{mvar|&mu;}} を用いて次のように定義される。
:<math>z(\beta,\mu)= e^{\beta \mu}</math>
{{mvar|&beta;}} は次のように定義される。
:<math>\beta = \frac{1}{kT}</math>
ここで {{mvar|k}} は[[ボルツマン定数]]、{{mvar|T}} は[[温度]]である。
全ての熱力学的な量はグランドカノニカル分布関数から導出されるため、全ての熱力学的な量は３つの変数 {{mvar|z}}、{{mvar|&beta;}}（または {{mvar|T}} ）、{{mvar|V}} のみの関数として考えることができる。
全ての偏微分係数は、３つの変数のうち１つを変数とし、残りの２つは定数とすることで求められる。
ここで次のように定義される無次元の[[グランドポテンシャル]]を扱うと便利である。
:<math>\Omega=-\ln(\mathcal{Z}) = \sum_i g_i \ln\left(1-ze^{-\beta\epsilon_i}\right).</math>
平均エネルギーは準位間のエネルギー差と比べて大きいと仮定する[[トーマス＝フェルミ近似]]を適用すると、上記の和は積分で置き換えられる。
:<math>\Omega\approx \int_0^\infty \ln\left(1-ze^{-\beta E}\right)\,dg.</math>
縮退度 {{mvar|dg}} は一般的な公式によって多くの異なる状況を表現する。
:<math>dg = \frac{1}{\Gamma(\alpha)}\,\frac{E^{\,\alpha-1}}{ E_c^{\alpha}} ~dE</math>
ここで {{mvar|&alpha;}} は定数、{{math|''E''<sub>c</sub>}} は臨界エネルギー、{{mvar|&Gamma;}} は[[ガンマ関数]]である。
たとえば箱の中の質量を持つボース気体では {{math|1=''&alpha;'' = 3/2}} で、臨界エネルギーは次のように与えられる。
:<math>\frac{1}{(\beta E_c)^\alpha}=\frac{Vf}{\Lambda^3}</math>
ここで {{math|&Lambda;}} は[[熱的波長]]である。調和トラップ中の質量を持つボース気体では {{math|1=''&alpha;'' = 3}} で、臨界エネルギーは次のように与えられる。
:<math>\frac{1}{(\beta E_c)^\alpha}=\frac{f}{(\hbar\omega\beta)^3}</math>
ここで {{math|1=''V''(''r'') = ''m&omega;''<sup>2</sup>''r''<sup>2</sup> / 2}} は調和ポテンシャルである。{{math|''E''<sub>c</sub>}} は体積だけの関数である。

このグランドポテンシャルの方程式は、項別に被積分関数のテイラー級数を積分することにより、 または {{math|Li<sub>1</sub>(''z''exp(&minus;''&beta;E''))}} の[[メリン変換]]に比例するとすることにより解くことができる。ここで {{math|Li<sub>s</sub>(''x'')}} は[[多重対数関数]]である。解は次のように与えられる。
:<math>\Omega\approx-\frac{\textrm{Li}_{\alpha+1}(z)}{\left(\beta E_c\right)^\alpha}</math>
このボース気体における連続体近似の問題点は、基底状態が実質的に無視されることで、ゼロエネルギーで縮退度がゼロになることである。
この問題点は[[ボース＝アインシュタイン凝縮]]を扱うときには重大で、次章で扱う。

== 基底状態の組み入れ ==
全[[粒子数]]はグランドポテンシャルから次のように与えられる。
:<math>N = -z\frac{\partial\Omega}{\partial z} \approx\frac{\textrm{Li}_\alpha(z)}{(\beta E_c)^\alpha}</math>
多重対数関数項は実で正でなければならず、最大値は {{math|1=''z'' = 1}} のときで {{math|''&zeta;''(''&alpha;'')}} に等しい。ここで {{mvar|&zeta;}} は[[リーマンゼータ関数]]である。
{{mvar|N}} を固定すると、{{mvar|&beta;}} の最大値は臨界値 {{math|''&beta;''<sub>c</sub>}} で、このとき以下のようになる。
:<math>N = \frac{\zeta(\alpha)}{(\beta_c E_c)^\alpha}</math>
これは臨界温度 {{math|1=''T''<sub>c</sub> = 1/''k&beta;''<sub>c</sub>}} に相当し、これ以下ではトーマス＝フェルミ近似は破綻する。
上記の方程式は臨界温度について解くことができ、次のようになる。
:<math>T_c=\left(\frac{N}{\zeta(\alpha)}\right)^{1/\alpha}\frac{E_c}{k}</math>
たとえば {{math|1=''&alpha;'' = 3/2 }} で、上述の値 {{math|''E''<sub>c</sub>}} を用いると、次のようになる。
:<math>T_c=\left(\frac{N}{Vf\zeta(3/2)}\right)^{2/3}\frac{h^2}{2\pi m k}</math>
さらに、ここでは臨界温度以下の結果を計算することはできない。なぜなら上記の方程式を用いた粒子数は負になるからである。
ここでの問題点は、トーマス＝フェルミ近似は基底状態の縮退度を０としていることで、これは間違っている。
凝縮を受け入れる基底状態が無いため方程式が破綻する。
しかし上記の方程式は励起状態では粒子数を比較的正確に評価しており、そこへ基底状態を単純に付け加えることは悪い近似ではないことがわかる。
:<math>N = N_0+\frac{\textrm{Li}_\alpha(z)}{(\beta E_c)^\alpha}</math>
ここで {{math|''N''<sub>0</sub>}} は基底状態凝縮の粒子数で、次のように与えられる。
:<math>N_0 = \frac{g_0\,z}{1-z}</math>

[[Image:Bose gas quantities.png|thumb|400px|right|図1：正規化された温度 {{mvar|&tau;}} と様々なボース気体パラメータとの関係。{{mvar|&alpha;}} の値は3/2。実線は {{math|1=''N'' = 10,000}} の場合で、点線は {{math|1=''N'' = 1,000}} の場合。黒線は励起粒子の割合。青線は凝集粒子の割合。赤線は化学ポテンシャルの負数。緑線は {{mvar|z}} の値。{{math|1=''k&epsilon;''<sub>c</sub> = 1}} と仮定した。]]

この方程式は絶対零度まで解くことができる。
図１に {{math|1=''&alpha;'' = 3/2}} におけるこの方程式の解の結果を示す。これは箱の中のボース気体に相当し、{{math|1=''k'' = ''&epsilon;''<sub>c</sub> = 1}} とする。
実線は {{math|1=''N'' = 10,000}} の場合、点線は {{math|1=''N'' = 1,000}} の場合を示す。黒線は励起粒子の割合 {{math|1 &minus; ''N''<sub>0</sub>/''N''}} 、青線は凝縮粒子の割合 {{math|''N''<sub>0</sub>/''N''}} で、赤線は化学ポテンシャル {{mvar|&mu;}} にマイナス符号をつけたもの、緑線は {{mvar|z}} の値である。 
横軸は次のように定義される正規化された温度 {{mvar|&tau;}} である。
:<math>\tau=\frac{T}{T_c}</math>
{{mvar|&mu;}} や {{mvar|z}} は低温の極限で {{mvar|&tau;<sup>&alpha;</sup>}} と線形になり、{{math|''N''<sub>0</sub>/''N''}} は高温の極限で {{math|1/''&tau;<sup>&alpha;</sup>''}} と線形になることがうかがえる。
粒子数が増加すると、凝縮の割合と励起の割合は臨界温度で不連続になる。

粒子数の方程式は正規化された温度で表される。
:<math>N = \frac{g_0\,z}{1-z}+N~\frac{\textrm{Li}_\alpha(z)}{\zeta(\alpha)}~\tau^\alpha</math>

{{mvar|N}} と {{mvar|&tau;}} が与えられると、この方程式は {{mvar|&tau;<sup>&alpha;</sup>}} について解くことができ、{{mvar|z}} についての級数解は、{{mvar|&tau;<sup>&alpha;</sup>}} のべき乗または {{mvar|&tau;<sup>&alpha;</sup>}} の逆べき乗での漸近展開としての[[級数の反転]]の方法によって得ることができる。
これらの展開から、{{math|1=''T'' = 0}} 近くでのガスの振る舞いを知ることができる。{{mvar|T}} が無限大でマクスウェル-ボルツマン分布となる。
特に我々が興味があるのは {{mvar|T}} が無限大のときで、これは上記の展開から容易に決定する。

== 熱力学 ==
粒子数に対する方程式に基底状態を加えることは、等価な基底状態の項をグランドポテンシャルに加えることに相当する。
:<math>\Omega = g_0\ln(1-z)-\frac{\textrm{Li}_{\alpha+1}(z)}{\left(\beta E_c\right)^\alpha}</math>
全ての熱力学的な性質は、グランドポテンシャルから計算される。
以下の表では低温と高温の極限、粒子数が無限大の極限での様々な熱力学的な量を示す。
厳密な結果は等号 (=) で表し、{{mvar|&tau;<sup>&alpha;</sup>}} の級数の最初の数項のみの結果は近似記号で表している。

{| border="1" cellpadding="4" cellspacing="0" align="center" style="margin: 0 0 1em; border-color:#ccc"
|-
! align="center" |量
! align="center" |一般の場合
! align="center" |<math>T \ll T_c\,</math>
! align="center" |<math>T \gg T_c\,</math>
|-
! align="center" |z
! align="center" |
! align="center" |<math>\approx \frac{\zeta(\alpha)}{\tau^\alpha}
-\frac{\zeta^2(\alpha)}{2^\alpha\tau^{2\alpha}}</math>
! align="center" |<math>=1\,</math>
|-
! align="center" |凝縮していない粒子の割合<br><math>1-\frac{N_0}{N}\,</math>
! align="center" |<math>=\frac{\textrm{Li}_{\alpha}(z)}{\zeta(\alpha)}
\,\tau^\alpha</math>
! align="center" |<math>=\tau^\alpha\,</math>
! align="center" |<math>=1\,</math>
|-
! align="center" |状態方程式<br><math>\frac{PV\beta}{N}=
-\frac{\Omega}{N}\,</math>
! align="center" |<math>=
\frac{\textrm{Li}_{\alpha\!+\!1}(z)}{\zeta(\alpha)}\,\tau^\alpha</math>
! align="center" |<math>=
\frac{\zeta(\alpha\!+\!1)}{\zeta(\alpha)}\,\tau^\alpha</math>
! align="center" |<math>\approx
1-\frac{\zeta(\alpha)}{2^{\alpha\!+\!1}\tau^\alpha}</math>
|-
! align="center" |ギブス自由エネルギー<br><math>G=\ln(z)\,</math>
! align="center" |<math>=\ln(z)\,</math>
! align="center" |<math>=0\,</math>
! align="center" |<math>\approx
\ln\left(\frac{\zeta(\alpha)}{\tau^\alpha}\right)
-\frac{\zeta(\alpha)}{2^{\alpha}\tau^\alpha}</math>
|}

全ての量は、高温の極限をとると古典的な理想気体の値に近づいていく。
上記の値を用いて、他の熱力学的な量を計算することができる。
たとえば、内部エネルギーと圧力×体積との関係は、すべての温度にわたって古典的な理想気体と同じである。
:<math>U=\frac{\partial \Omega}{\partial \beta}=\alpha PV</math>
定積比熱においても同様である。
:<math>C_v=\frac{\partial U}{\partial T}=k(\alpha+1)\,U\beta</math>
エントロピーは次式で与えられる。
:<math>TS=U+PV-G\,</math>
高温の極限をとると、次式が得られる。
:<math>TS=(\alpha+1)+\ln\left(\frac{\tau^\alpha}{\zeta(\alpha)}\right)</math>
これは {{math|1=''&alpha;'' = 3/2}} では単なる[[ザックール・テトローデ方程式]]を書き換えたものである。
デルタ相互作用をもつ1次元ボース粒子はフェルミ粒子として振る舞い、[[パウリの排他原理]]に従う。
デルタ相互作用をもつ1次元ボース粒子は[[ベーテ仮設]]により厳密に解くことができる。
バルク自由エネルギーと熱力学的ポテンシャルは[[楊振寧]]によって計算された。
１次元の場合における相関関数[https://books.google.com/books?id=kaZ0pKIHhxAC&dq=quantum+inverse+scattering+method&printsec=frontcover#v=onepage&q&f=false]も評価された。
１次元ボース気体は量子的な[[非線形シュレーディンガー方程式]]と等価である。

== 関連項目 ==
* [[フェルミ気体]]
* [[箱の中の気体]]
* [[デバイ模型]]
* [[ボース＝アインシュタイン凝縮]]
* [[黒体放射]]

== 参考文献 ==
{{reflist}}
* {{cite book |last=Huang |first=Kerson |title=Statistical Mechanics 
 |year=1967 |publisher=John Wiley and 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, 3rd Edition Part 1 
 |year=1996 |publisher=Butterworth-Heinemann |location=Oxford  }}
* {{cite book |last=Pethick |first=C. J.|author2=H. Smith 
 |title=Bose–Einstein Condensation in Dilute Gases |year=2004 
 |publisher=[[Cambridge University Press]] |location=Cambridge  }}
* {{cite journal | last = Yan | first = Zijun | year = 2000
 | title = General Thermal Wavelength and its Applications | journal = Eur. J. Phys
 | volume = 21 | pages = 625–631 
 | url = http://www.iop.org/EJ/article/0143-0807/21/6/314/ej0614.pdf | format = PDF
 | doi = 10.1088/0143-0807/21/6/314|bibcode = 2000EJPh...21..625Y | issue = 6 }}

{{Normdaten}}
{{デフォルトソート:ほおすきたい}}
[[Category:物性物理学]]
[[Category:量子力学]]
[[Category:熱力学]]
[[Category:サティエンドラ・ボース]]
[[Category:物理学のエポニム]]