tecktom.com

相互作用の無い系ではモードごとの積に分かれる

2026年10月17日12分で読めますAI生成教科書『虎の巻』第6部 第5章 熱平衡と占有数:Bose 分布と Fermi 分布
目次
  1. Fock 空間の指数法則
  2. 分布の定理
  3. 積の対数を展開する:循環の和
  4. 木下の層

AI 生成・木下が未勉強

量子力学の様々な構成方法 第6部 第5章 §6.5.2 ★★★ Fock 空間が1粒子空間の直和に対してテンソル積に分かれるので、自由な大正準の状態はモードの積になり、Bose 分布と Fermi 分布が出る

Fock 空間の指数法則

1粒子の正規直交基底 {φα}α=1M\{\varphi_\alpha\}_{\alpha=1}^M を取る(MM は有限とし、無限の場合は極限で考える)。第2章の占有数の基底∣n1,…,nM⟩=∏α(a^α†)nα∣0⟩/∏αnα!| n_1,\dots,n_M \rangle =\prod_\alpha(\hat a_\alpha^\dagger)^{n_\alpha}| 0 \rangle /\sqrt{\prod_\alpha n_\alpha!}(フェルミ粒子では nα∈{0,1}n_\alpha\in\{0,1\} で、積の順序は第3章の約束どおり番号の小さいモードを左に置く)はFock 空間の正規直交基底である。

補題 6.5.2(Fock 空間の指数法則)

写像 ∣n1,…,nM⟩↦∣n1⟩⊗⋯⊗∣nM⟩| n_1,\dots,n_M \rangle \mapsto| n_1 \rangle \otimes\dots\otimes| n_M \rangle は、F±(CM)\mathcal F_\pm(\mathbb{C}^M) から1モードの Fock 空間の MM 個のテンソル積へのユニタリ写像で、各モードの個数演算子 n^α=a^α†a^α\hat n_\alpha=\hat a_\alpha^\dagger\hat a_\alpha を、α\alpha 番目の因子だけに働く個数演算子 1⊗⋯⊗n^⊗⋯⊗11\otimes\dots\otimes\hat n\otimes\dots\otimes1 に移す。

証明

正規直交基底を正規直交基底に1対1に移すのでユニタリである。n^α\hat n_\alpha は占有数の基底で対角で、固有値は nαn_\alpha そのものである。フェルミ粒子の場合、c^α\hat c_\alpha を単独で移すと第3章の Jordan–Wigner の文字列(符号)が付くが、n^α=c^α†c^α\hat n_\alpha=\hat c_\alpha^\dagger\hat c_\alpha では符号が2回現れて打ち消し合う。したがって両方の統計で、n^α\hat n_\alpha の像は α\alpha 番目の因子の n^\hat n になる。

これは「1粒子空間を直和 CM=C⊕⋯⊕C\mathbb{C}^M=\mathbb{C}\oplus\dots\oplus\mathbb{C} に分けると、Fock 空間はテンソル積に分かれる」という性質で、指数関数の ea+b=eaeb\mathrm{e}^{a+b}=\mathrm{e}^a\mathrm{e}^b になぞらえて指数法則と呼ばれる。

分布の定理

定理 6.5.3(自由な大正準の状態は積の状態 ★★★)

H^=∑αεαn^α\hat H=\sum_\alpha\varepsilon _\alpha\hat n_\alpha(第4章の命題で2次形式を対角化した後の形)とし、ボース粒子では μ<min⁡αεα\mu<\min_\alpha\varepsilon _\alpha とする。このとき

ρ^G=⨂αe−β(εα−μ)n^Ξα,Ξ=∏αΞα=∏α(1∓e−β(εα−μ))∓1(6.5.7)\hat\rho_{\mathrm G}=\bigotimes_\alpha\frac{\mathrm{e}^{-\beta(\varepsilon _\alpha-\mu)\hat n}}{\Xi_\alpha},\qquad \Xi=\prod_\alpha\Xi_\alpha=\prod_\alpha\bigl(1\mp\mathrm{e}^{-\beta(\varepsilon _\alpha-\mu)}\bigr)^{\mp1}\tag{6.5.7}

で、平均の占有数は

⟨a^α†a^α′⟩=δαα′ 1eβ(εα−μ)∓1(ボースは nB(εα)、フェルミは f(εα))(6.5.8)\langle \hat a_\alpha^\dagger\hat a_{\alpha'} \rangle =\delta_{\alpha\alpha'}\,\frac{1}{\mathrm{e}^{\beta(\varepsilon _\alpha-\mu)}\mp1} \qquad\bigl(\text{ボースは }n_{\mathrm B}(\varepsilon _\alpha)\text{、フェルミは }f(\varepsilon _\alpha)\bigr)\tag{6.5.8}

である。異なるモードの占有数は独立で、⟨n^αn^α′⟩=⟨n^α⟩⟨n^α′⟩\langle \hat n_\alpha\hat n_{\alpha'} \rangle =\langle \hat n_\alpha \rangle \langle \hat n_{\alpha'} \rangle(α≠α′\alpha\ne\alpha')が成り立つ。

証明

K^=∑α(εα−μ)n^α\hat K=\sum_\alpha(\varepsilon _\alpha-\mu)\hat n_\alpha は互いに可換な項の和なので e−βK^=∏αe−β(εα−μ)n^α\mathrm{e}^{-\beta\hat K}=\prod_\alpha\mathrm{e}^{-\beta(\varepsilon _\alpha-\mu)\hat n_\alpha} で、補題 6.5.2 の写像で各因子は別々のテンソル因子に働く。テンソル積のトレースは各因子のトレースの積なので Ξ=∏αΞα\Xi=\prod_\alpha\Xi_\alpha になり、Ξα\Xi_\alpha は例題 6.5.1 の1モードの値 (1−e−yα)−1(1-\mathrm{e}^{-y_\alpha})^{-1}(ボース)、1+e−yα1+\mathrm{e}^{-y_\alpha}(フェルミ)である(yα=β(εα−μ)y_\alpha=\beta(\varepsilon _\alpha-\mu))。対角の期待値 ⟨n^α⟩\langle \hat n_\alpha \rangle は α\alpha 番目の因子だけの計算になり、例題 6.5.1 の nˉ\bar n を与える。非対角 α≠α′\alpha\ne\alpha' では、a^α†a^α′\hat a_\alpha^\dagger\hat a_{\alpha'} は占有数の基底の状態を、nαn_\alpha を1増やし nα′n_{\alpha'} を1減らした別の基底の状態へ移す(フェルミでは符号が付くだけ)。ρ^G\hat\rho_{\mathrm G} はこの基底で対角なので、トレースの対角成分は 0 になる。独立性は積の状態の性質そのものである。

(6.5.8) が Bose 分布と Fermi 分布である。分母の ∓1\mp1 の違いは、1つのモードに入れる個数の上限が無いか1か、すなわち交換関係か反交換関係か(第3章)の違いだけで決まっている。eβ(ε−μ)≫1\mathrm{e}^{\beta(\varepsilon -\mu)}\gg1 では両方が e−β(ε−μ)\mathrm{e}^{-\beta(\varepsilon -\mu)} に近づき、粒子が同じモードに重なることがほとんど無くなるので統計の違いが消える。これが Maxwell–Boltzmann 分布の極限である。3つを並べたのが図 6.5.1 である。

図

図 6.5.1 1つのモードの平均の占有数を y=β(ε−μ)y=\beta(\varepsilon -\mu) の関数として並べた(点線は厳密な Bose・Fermi 分布、灰色の実線は Maxwell–Boltzmann の近似)。Bose 分布は y→0+y\to0^+ で発散し、Fermi 分布は 1 を超えない。y≳3y\gtrsim3 では3つがほぼ重なる。

積の対数を展開する:循環の和

(6.5.7) の対数を逃散能 zz で展開すると、第2章の置換の和との関係が見える。

命題 6.5.4(循環展開 ★★★)

1粒子の分配関数を Z1(β)≔∑αe−βεαZ_1(\beta)\coloneqq\sum_\alpha\mathrm{e}^{-\beta\varepsilon _\alpha} とする。z e−βεmin⁡<1z\,\mathrm{e}^{-\beta\varepsilon _{\min}}<1 で

ln⁡Ξ=∑ℓ=1∞(±1)ℓ+1ℓ zℓ Z1(ℓβ)(6.5.9)\ln\Xi=\sum_{\ell=1}^\infty\frac{(\pm1)^{\ell+1}}{\ell}\,z^\ell\,Z_1(\ell\beta)\tag{6.5.9}

が成り立つ。特に Ξ\Xi の z2z^2 の係数は Z2±=12[Z1(β)2±Z1(2β)]Z_2^\pm=\frac12\bigl[Z_1(\beta)^2\pm Z_1(2\beta)\bigr] で、第2章の式(置換の和)と一致する。

証明

(6.5.7) から ln⁡Ξ=∓∑αln⁡(1∓ze−βεα)\ln\Xi=\mp\sum_\alpha\ln(1\mp z\mathrm{e}^{-\beta\varepsilon _\alpha}) である。∣w∣<1\lvert w\rvert <1 で −ln⁡(1−w)=∑ℓwℓ/ℓ-\ln(1-w)=\sum_\ell w^\ell/\ell を使うと、ボースでは ∑α∑ℓzℓe−ℓβεα/ℓ\sum_\alpha\sum_\ell z^\ell\mathrm{e}^{-\ell\beta\varepsilon _\alpha}/\ell、フェルミでは ln⁡(1+w)=∑ℓ(−1)ℓ+1wℓ/ℓ\ln(1+w)=\sum_\ell(-1)^{\ell+1}w^\ell/\ell から符号 (−1)ℓ+1(-1)^{\ell+1} が付く。α\alpha の和を先に取ると ∑αe−ℓβεα=Z1(ℓβ)\sum_\alpha\mathrm{e}^{-\ell\beta\varepsilon _\alpha}=Z_1(\ell\beta) で、(6.5.9) を得る。Ξ=exp⁡(ln⁡Ξ)\Xi=\exp(\ln\Xi) を z2z^2 まで展開すると 1+zZ1(β)+z2[12Z1(β)2±12Z1(2β)]+O(z3)1+zZ_1(\beta)+z^2\bigl[\frac12Z_1(\beta)^2\pm\frac12Z_1(2\beta)\bigr]+\mathcal{O}(z^3) で、(6.5.3) と比べれば Z2±Z_2^\pm が読める。

ℓ\ell 番目の項は、第2章で長さ ℓ\ell の巡回置換が与えた Tr⁡(PγA^⊗ℓ)=Tr⁡A^ℓ\operatorname{Tr}(P_\gamma\hat A^{\otimes\ell})=\operatorname{Tr}\hat A^\ell と同じ量である。ℓ\ell 個の粒子が順に入れ替わって輪になると、1粒子が温度 1/ℓ1/\ell 倍の熱浴にいるのと同じ Z1(ℓβ)Z_1(\ell\beta) を与え、フェルミ粒子では輪の長さが偶数のとき符号が −1-1 になる。正準の ZNZ_N では置換の分解(どの長さの輪がいくつあるか)を全部数える必要があったが、大正準では ln⁡Ξ\ln\Xi が輪1つずつの和になり、組み合わせの計算が消える。虚時間の経路積分では、この輪は「ℓ\ell 回巻いてから閉じる経路」であり、解析力学から様々な理論分野を組み立てる 第5部第1章が入れ替えた経路の輪として同じ式を出している。

例題 6.5.2 3つの準位のフェルミ粒子を数え上げる ★★

オリジナルスピンを持たないフェルミ粒子が入る準位が 0, Δε, 2Δε0,\ \Delta\varepsilon ,\ 2\Delta\varepsilon の3つだけある。μ=Δε\mu=\Delta\varepsilon、βΔε=ln⁡2\beta\Delta\varepsilon =\ln2 とする。

  1. (1) Fock 空間の8つの基底の状態を書き、Ξ\Xi を状態の和で求めよ。定理 6.5.3 の積と一致することを確かめよ。
  2. (2) 3つの準位の占有数と ⟨N^⟩\langle \hat N \rangle を求めよ。
  3. (3) 粒子数が N=0,1,2,3N=0,1,2,3 である確率を求めよ。
解答

(1) 状態は占有数 (n1,n2,n3)∈{0,1}3(n_1,n_2,n_3)\in\{0,1\}^3 の8つで、重みは e−β∑α(εα−μ)nα\mathrm{e}^{-\beta\sum_\alpha(\varepsilon _\alpha-\mu)n_\alpha} である。εα−μ=−Δε, 0, +Δε\varepsilon _\alpha-\mu=-\Delta\varepsilon ,\,0,\,+\Delta\varepsilon なので、各準位の因子は eβΔε=2\mathrm{e}^{\beta\Delta\varepsilon }=2、11、e−βΔε=12\mathrm{e}^{-\beta\Delta\varepsilon }=\frac12 で、

Ξ=1⏟(000)+2+1+12⏟N=1+2⋅1+2⋅12+1⋅12⏟N=2+2⋅1⋅12⏟N=3=1+72+72+1=9.\Xi=\underbrace{1}_{(000)}+\underbrace{2+1+\tfrac12}_{N=1}+\underbrace{2\cdot1+2\cdot\tfrac12+1\cdot\tfrac12}_{N=2}+\underbrace{2\cdot1\cdot\tfrac12}_{N=3} =1+\tfrac72+\tfrac72+1=9.

積の公式では (1+2)(1+1)(1+12)=9(1+2)(1+1)(1+\frac12)=9 で一致する。8項を展開した形が積の展開そのものになっていることが、指数法則の中身である。

(2) f=1/(ey+1)f=1/(\mathrm{e}^{y}+1) に ey=12,1,2\mathrm{e}^y=\frac12,1,2 を入れて f=23, 12, 13f=\frac23,\ \frac12,\ \frac13、⟨N^⟩=32\langle \hat N \rangle =\frac32 である。μ\mu と同じ高さの準位がちょうど半分埋まり、上下が 12\frac12 を中心に対称に並ぶ。

(3) 状態の和を NN ごとに分けて Ξ\Xi で割ると、P(0)=19P(0)=\frac19、P(1)=7/29=718P(1)=\frac{7/2}{9}=\frac7{18}、P(2)=718P(2)=\frac7{18}、P(3)=19P(3)=\frac19 である。大正準の状態は粒子数の違う状態の混合で、⟨N^⟩=32\langle \hat N \rangle =\frac32 の周りに Var(N^)=∑αfα(1−fα)=29+14+29=2536\mathrm{Var}(\hat N)=\sum_\alpha f_\alpha(1-f_\alpha)=\frac29+\frac14+\frac29=\frac{25}{36} だけ揺らぐ(∑NN2P(N)−94=718+2818+1−94=2536\sum_NN^2P(N)-\frac94=\frac7{18}+\frac{28}{18}+1-\frac94=\frac{25}{36} と一致する)。N=1N=1 の中の3つの状態の比 2:1:122:1:\frac12 は、1粒子の正準分布 e−βεα\mathrm{e}^{-\beta\varepsilon _\alpha} の比そのもので、6.5.6 節で示す「NN で切ると正準の状態になる」ことの最小の例になっている。数値は verify_ch6-05.py で確かめた。

1つのフェルミのモードは、空と占有の2つの状態を持つ2準位系である(第3章の Jordan–Wigner の対応)。μ=0\mu=0 でエネルギー ℏω0\hbar\omega_0 のモードを熱浴に置くと、占有の確率と空の確率の比は e−βℏω0\mathrm{e}^{-\beta\hbar\omega_0} で、共通例題②の2準位系の熱平衡(第4部第4章の例題の詳細つり合い Pe/Pg=e−βℏω0P_e/P_g=\mathrm{e}^{-\beta\hbar\omega_0})と同じ式になる。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ