tecktom.com

この章のまとめ

2026年10月15日9分で読めますAI生成教科書『虎の巻』第5部 第3章 光子とフォノンを同じ形で数える

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第5部 第3章 §5.3.7 ★★★ 光子とフォノンは、2次式の場の LL を基準振動に分けた振動子の段の数。熱の量は g(ω)g(\omega) と振動子1つの量の積分で、違いは g(ω)g(\omega) だけ

この章では、空洞の電磁場と結晶の格子振動に同じ形の Hamiltonian H^=∑αℏωα(n^α+12)\hat H=\sum_\alpha\hbar\omega_\alpha(\hat n_\alpha+\frac12) を選んだ。どちらも Lagrangian が場の変数の2次式なので、基準振動に分けると独立な調和振動子の和になり(電磁場は厳密に、格子振動は調和近似で)、各振動子を第4部第4章の方法で量子化すると、段の数 nαn_\alpha が光子やフォノンの数になる。振動子が独立なので分配関数は積に分かれ、占有数の平均は μ=0\mu=0 の Bose 分布 nB(ω)n_{\mathrm B}(\omega) になる。μ=0\mu=0 は、段の数を固定する条件が無いことの表れである。

熱力学の量はすべて、モードの密度 g(ω)g(\omega) と振動子1つの量の積の積分で書くことができる。光子では gph=Vω2/(π2c3)g_{\mathrm{ph}}=V\omega^2/(\pi^2c^3) を入れて、Planck の式・U=arVT4U=a_{\mathrm r}VT^4・F=−U/3F=-U/3・P=U/(3V)P=U/(3V)・S=4U/(3T)S=4U/(3T) が出た。P=U/(3V)P=U/(3V) は、各モードの段の数が断熱不変量で、ω∝V−1/3\omega\propto V^{-1/3} であることからも出る。フォノンでは、低い ω\omega の g(ω)g(\omega) が光子と同じ ω2\omega^2 の形なので、低温の比熱は光子気体と同じ T3T^3 則になる。モードが 3N3N 個で止まるので高温では Dulong–Petit に戻り、Debye のモデルはこの両端を正しく持つ内挿の式である。零点エネルギーは温度によらず熱の量に入らないが、壁の位置による差が Casimir 力として観測され、占有数の揺らぎは粒子の項 nBn_{\mathrm B} と波の項 nB2n_{\mathrm B}^2 の和になる。

この章では、モードの k\bm{k} に上限がある格子と、上限が無い電磁場を並べた。電磁場の零点エネルギーの発散は、高い ω\omega のモードを弱める因子を入れて計算し、結果のうち因子によらない部分だけを物理として取り出した。この操作は、第8章から第10章のループの発散と正則化で一般の形に広げる(第5部第9章)。光子とフォノンを相互作用させる(電子が光子やフォノンを吸収・放出する)ことは、第5章のゲージ原理と、光学と量子エレクトロニクスの巻 第6部で扱う(第5部第5章・光学と量子エレクトロニクス 第30章)。

総合問題 振動子の気体として、原子の列と空洞の光を並べる ★★★

オリジナル質量 mm の原子 NN 個が間隔 aa でばね定数 KK のばねにつながった1次元の列(周期境界)と、体積 VV の空洞の中の電磁場を、同じ温度 TT で比べる。ωm=2K/m\omega_{\mathrm m}=2\sqrt{K/m} とする。

  1. (1) (使う道具:5.3.1節/第2部第4章・第4部第4章)第2章の総合問題 (1)(2) で求めた、列の生成消滅演算子による Hamiltonian と零点エネルギーを、この章の ωm\omega_{\mathrm m} で書き直せ。
  2. (2) (使う道具:5.3.2節・5.3.4節/第3部第6章)列の自由エネルギーを g(ω)g(\omega) で書き、低温の比熱を求めよ。
  3. (3) (使う道具:5.3.3節/第3部第4章)空洞の光子気体の F(T,V)F(T,V) から PP, SS を求め、U=3PVU=3PV を示せ。さらに U(S,V)U(S,V) を Legendre 変換して F(T,V)F(T,V) に戻ることを確かめよ。
  4. (4) (使う道具:5.3.4節)3次元の結晶で、平均の音速 vˉ=3000 m s−1\bar v=3000\,\mathrm{m\,s^{-1}} のとき、低温で同じ体積・同じ温度のフォノンの気体と光子気体のエネルギーの比を求めよ。
  5. (5) (使う道具:5.3.3節/第1部第14章・第3部第5章)空洞の1辺をゆっくり ss 倍にすると、光子気体の温度とエントロピーはどう変わるか。温度 3000 K3000\,\mathrm K の放射が 2.72548 K2.72548\,\mathrm K まで冷えるには、ss はいくらか。
解答

(1) 〔第2部第4章・第4部第4章/第5部第2章〕第5部第2章の5.2.7節の総合問題 (1)(2) の結果で、その章の MM を mm、2ω02\omega_0 を ωm\omega_{\mathrm m} と読み替えるだけである。ωk=ωm∣sin⁡(ka/2)∣\omega_k=\omega_{\mathrm m}\lvert \sin(ka/2)\rvert で、k=0k=0 は列全体の並進なので振動子にならず、

H^=P^022m+∑k≠0ℏωk(a^k†a^k+12),E0N  →N→∞  ℏωmπ.\hat H=\frac{\hat P_0^2}{2m}+\sum_{k\ne0}\hbar\omega_k\Bigl(\hat a_k^\dagger\hat a_k+\frac12\Bigr),\qquad \frac{E_0}{N}\;\xrightarrow{N\to\infty}\;\frac{\hbar\omega_{\mathrm m}}{\pi}.

N=1000N=1000 で和を直接計算すると E0/(Nℏωm)=0.3183E_0/(N\hbar\omega_{\mathrm m})=0.3183 で、1/π=0.31831/\pi=0.3183 と一致する。

(2) 〔2節・4節・第3部第6章〕式 (5.3.10) を g(ω)g(\omega) で書いてF=E0+kBT∫0ωmg(ω)ln⁡(1−e−βℏω) dωF=E_0+ k_{\mathrm B} T\int_0^{\omega_{\mathrm m}}g(\omega)\ln(1-\mathrm{e}^{-\beta\hbar\omega})\,\mathrm{d}\omega、g(ω)=2Nπ(ωm2−ω2)−1/2g(\omega)=\frac{2N}{\pi}(\omega_{\mathrm m}^2-\omega^2)^{-1/2}。低温では g(ω)≈g(0)=2N/(πωm)g(\omega)\approx g(0)=2N/(\pi\omega_{\mathrm m}) として、例題 5.3.4 (1) から C=2π3NkB kBT/(ℏωm)C=\frac{2\pi}{3}N k_{\mathrm B}\, k_{\mathrm B} T/(\hbar\omega_{\mathrm m})。自由エネルギーの温度による部分は ∫0∞ln⁡(1−e−x) dx=−π2/6\int_0^\infty\ln(1-\mathrm{e}^{-x})\,\mathrm{d} x=-\pi^2/6 から F−E0=−π3N(kBT)2ℏωmF-E_0=-\frac{\pi}{3}N\frac{( k_{\mathrm B} T)^2}{\hbar\omega_{\mathrm m}} で、C=−T∂2F/∂T2C=-T\partial^2F/\partial T^2 が上の CC に一致する。

(3) 〔3節・第3部第4章〕F=−ar3VT4F=-\frac{a_{\mathrm r}}3VT^4 から P=−∂F/∂V=ar3T4P=-\partial F/\partial V=\frac{a_{\mathrm r}}{3}T^4、S=−∂F/∂T=43arVT3S=-\partial F/\partial T=\frac43a_{\mathrm r}VT^3、U=F+TS=arVT4=3PVU=F+TS=a_{\mathrm r}VT^4=3PV。U(S,V)=KS4/3V−1/3U(S,V)=KS^{4/3}V^{-1/3}(K=ar−1/3(3/4)4/3K=a_{\mathrm r}^{-1/3}(3/4)^{4/3})から T=∂U/∂S=43KS1/3V−1/3T=\partial U/\partial S=\frac43KS^{1/3}V^{-1/3}。Legendre 変換 F=inf⁡S[U−TS]F=\inf_S[U-TS] の停留点はこの式を SS について解いた S=(3T/4K)3VS=(3T/4K)^3V で、代入するとU−TS=KS4/3V−1/3−43KS4/3V−1/3=−13KS4/3V−1/3=−13K(3T4K)4V=−3344T4K3VU-TS=KS^{4/3}V^{-1/3}-\frac43KS^{4/3}V^{-1/3}=-\frac13KS^{4/3}V^{-1/3}=-\frac13K\bigl(\frac{3T}{4K}\bigr)^4V=-\frac{3^3}{4^4}\frac{T^4}{K^3}V。K3=ar−1(3/4)4K^3=a_{\mathrm r}^{-1}(3/4)^4 を入れると −3344⋅4434arVT4=−ar3VT4-\frac{3^3}{4^4}\cdot\frac{4^4}{3^4}a_{\mathrm r}VT^4=-\frac{a_{\mathrm r}}3VT^4 で、元の FF に戻る。UU が SS について下に凸なので、停留点は最小である。

(4) 〔4節〕式 (5.3.26) と式 (5.3.21) の比はπ2/10⋅(ℏvˉ)−3π2/15⋅(ℏc)−3=32(cvˉ)3=32×(9.993×104)3=1.50×1015\dfrac{\pi^2/10\cdot(\hbar\bar v)^{-3}}{\pi^2/15\cdot(\hbar c)^{-3}}=\dfrac32\Bigl(\dfrac{c}{\bar v}\Bigr)^3=\dfrac32\times(9.993\times10^4)^3=1.50\times10^{15}。同じ体積と温度で、フォノンの気体は光子気体の 101510^{15} 倍のエネルギーを持つ。固体を温めるとき、空洞の光のエネルギーは比べものにならないほど小さい。

(5) 〔3節・第1部第14章・第3部第5章〕ゆっくり広げると各モードの段の数 nαn_\alpha は断熱不変量として保たれ(第1部第14章の作用変数 J=E/ωJ=E/\omega が保たれることの量子版)、ωα∝1/Lb\omega_\alpha\propto1/L_{\mathrm b} は 1/s1/s 倍になる。nα=nB(ωα)n_\alpha=n_{\mathrm B}(\omega_\alpha) が保たれるためには βℏωα\beta\hbar\omega_\alpha が一定でなければならないので T→T/sT\to T/s。エントロピーは S=43arVT3∝s3⋅s−3S=\frac43a_{\mathrm r}VT^3\propto s^3\cdot s^{-3} で変わらない。これは第3部第5章の「断熱変化=エントロピー一定」の例になっている。各モードの占有の確率分布 (5.3.12) が変わらないので、分布の言葉でもエントロピーは変わらない。3000/2.72548=1.10×1033000/2.72548=1.10\times10^3 なので、s=1.10×103s=1.10\times10^3。

木下の勉強ステータス

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

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

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