tecktom.com

授業で習った形との対応

2026年10月15日6分で読めますAI生成教科書『虎の巻』第3部 第2章 等重率とミクロカノニカル
目次
  1. 量子の状態を数える形
  2. 何を一様に数えるか
  3. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第3部 第2章 §3.2.6 ★★ 授業の「微視的状態を等確率で数える」は、相空間を h3NN!h^{3N}N! のます目に区切って数えることにあたる。何を一様に数えるかは Liouville の定理が決める

量子の状態を数える形

学部の統計力学の授業では、多くの場合、量子力学のエネルギー固有状態を数えるところから始める。その形では、孤立系で、エネルギーが EE と E+δEE+\delta E の間にある量子状態(H^\hat H の固有状態)の数を W(E)W(E) とし、そのすべてが等しい確率で実現すると仮定する。エントロピーを S=kBln⁡W(E)S= k_{\mathrm B}\ln W(E)、温度を 1/T=∂S/∂E1/T=\partial S/\partial E で定める。

この章の形との対応は、次の3点にまとめられる。第1に、授業の「量子状態の数」は、この章の W(E)=∫殻dΓW(E)=\int_{\text{殻}}\mathrm{d}\Gamma にあたる。dΓ\mathrm{d}\Gamma の分母 h3NN!h^{3N}N! が、相空間の体積を量子状態の数に換算する係数である(3.2.3節)。第2に、授業では「等確率」を仮定として置くが、古典の形では「相空間の体積 dΓ\mathrm{d}\Gamma で一様」と読み、それが Liouville の定理で時間発展と両立する唯一の選び方であることを確かめられる(3.2.2節)。量子力学では、固有状態の上で一様な密度行列 ρ^∝∑E<En≤E+δE∣n⟩⟨n∣\hat\rho\propto\sum_{E<E_n\le E+\delta E}| n \rangle \langle n | が H^\hat H と可換なので、Liouville 方程式の量子版 ∂ρ^∂t=[H^, ρ^]/(iℏ)\frac{\partial \hat\rho}{\partial t}=\left[\hat H,\,\hat\rho\right]/(\mathrm{i}\hbar) で時間に依らない。第3に、授業で δE\delta E の大きさに触れずに済ませるのは、3.2.5節で見たとおり δE\delta E が結果に関係しないからである。

量子の数え方と古典の数え方を、独立な NN 個の調和振動子(Einstein の固体の模型)で並べる。量子力学では、全部で MM 個のエネルギー量子 ℏω\hbar\omega を NN 個の振動子に分ける方法の数

Wq=(M+N−1)!M! (N−1)!(3.2.23)W_{\mathrm q}=\frac{(M+N-1)!}{M!\,(N-1)!}\tag{3.2.23}

が、エネルギー E=(M+N2)ℏωE=(M+\frac N2)\hbar\omega の状態の数である。古典では、振動子は格子点ごとに置かれていて互いに区別できるので N!N! で割らず、Φ(E)=(E/ℏω)N/N!\Phi(E)=(E/\hbar\omega)^N/N!(総合問題の (1) で導く。ここの N!N! は単体の体積から出る数で、粒子の入れ替えとは関係がない)。殻の厚みを準位の間隔 ΔE=ℏω\Delta E=\hbar\omega に取ると Wcl=g ℏω=(E/ℏω)N−1/(N−1)!W_{\mathrm{cl}}=g\,\hbar\omega=(E/\hbar\omega)^{N-1}/(N-1)! である。N=2N=2 では Wq=M+1W_{\mathrm q}=M+1、Wcl=M+1W_{\mathrm{cl}}=M+1 で厳密に一致する。N=3N=3 では Wq=(M+2)(M+1)2W_{\mathrm q}=\frac{(M+2)(M+1)}{2}、Wcl=(M+3/2)22W_{\mathrm{cl}}=\frac{(M+3/2)^2}{2} で差は 18\frac18 だけである。N=100N=100, M=106M=10^6 では ln⁡Wq\ln W_{\mathrm q} と ln⁡Wcl\ln W_{\mathrm{cl}} はどちらも 1008.6061008.606 で、小数第3位まで一致する。量子の数え上げと相空間の体積は、量子の数 MM が振動子の数 NN より十分多い(高温の)ときに一致する。

何を一様に数えるか

授業の形では「状態」が初めから数えられる物として与えられるので、「何で測って一様か」という問いは表に出ない。古典の形ではこの問いが残り、答えは正準変数 (q,p)(q,p) の体積である。位置と速度 (q,q˙)(q,\dot q) で一様に数えると、質量が場所によって変わる系では違う答えになる。次の例題で確かめる。

例題 3.2.8 位置で質量が変わる粒子 ★★★

オリジナル区間 0≤q≤a0\le q\le a を、両端の壁で跳ね返りながら自由に動く粒子を考える。Lagrangian は L=12m(q)q˙2L=\frac12m(q)\dot q^2、m(q)=m0(1+q/a)2m(q)=m_0(1+q/a)^2 とする(位置で慣性が変わる。曲がった管の中のビーズや、回転する座標の自由度に現れる形)。エネルギーを EE とする。

  1. (1) 正準運動量 pp と Hamiltonian を書き、殻の上で (q,p)(q,p) について一様な分布での qq の確率密度 P(q)P(q) を求めよ。区間の右半分 q>a/2q>a/2 にいる確率を求めよ。
  2. (2) 殻の上で (q,q˙)(q,\dot q) について一様な分布(速度で測って一様)では、右半分にいる確率はいくつか。
  3. (3) 実際の運動で、右半分にいる時間の割合を求め、(1) と (2) のどちらと一致するかを述べよ。
解答

(1) p=∂L/∂q˙=m(q)q˙p=\partial L/\partial\dot q=m(q)\dot q、H=p2/(2m(q))H=p^2/(2m(q))。殻の上の重みは ∫δ(E−H) dp=∑±1/∣∂H/∂p∣=2m(q)/∣p∣=2m(q)/(2E)\int\delta(E-H)\,\mathrm{d} p=\sum_\pm1/\lvert \partial H/\partial p\rvert =2m(q)/\lvert p\rvert =2\sqrt{m(q)/(2E)} で、m(q)∝1+q/a\sqrt{m(q)}\propto1+q/a に比例する。∫0a(1+q/a) dq=32a\int_0^a(1+q/a)\,\mathrm{d} q=\frac32a で割って P(q)=23a(1+q/a)P(q)=\frac{2}{3a}(1+q/a)。右半分の確率は

23a∫a/2a(1+qa)dq=23a(a2+3a8)=712≈0.583.\frac{2}{3a}\int_{a/2}^{a}\Bigl(1+\frac qa\Bigr)\mathrm{d} q=\frac{2}{3a}\Bigl(\frac a2+\frac{3a}{8}\Bigr)=\frac{7}{12}\approx0.583.

(2) 速度で測ると、E=12m(q)q˙2E=\frac12m(q)\dot q^2 の殻の上の重みは ∫δ(E−12mq˙2) dq˙=2/(m∣q˙∣)=2/2Em(q)∝1/(1+q/a)\int\delta(E-\frac12m\dot q^2)\,\mathrm{d}\dot q=2/(m\lvert \dot q\rvert )=2/\sqrt{2Em(q)}\propto1/(1+q/a) である。∫0adq1+q/a=aln⁡2\int_0^a\frac{\mathrm{d} q}{1+q/a}=a\ln2、∫a/2adq1+q/a=aln⁡43\int_{a/2}^a\frac{\mathrm{d} q}{1+q/a}=a\ln\frac43 なので、右半分の確率は ln⁡(4/3)/ln⁡2≈0.415\ln(4/3)/\ln2\approx0.415 になる。

(3) 速さは ∣q˙∣=2E/m(q)∝1/(1+q/a)\lvert \dot q\rvert =\sqrt{2E/m(q)}\propto1/(1+q/a) で、dq\mathrm{d} q を通るのにかかる時間 dq/∣q˙∣\mathrm{d} q/\lvert \dot q\rvert は 1+q/a1+q/a に比例する。時間の割合は (1) と同じ 7/127/12 である。質量の大きい右側では粒子がゆっくり動くので長く滞在し、その重みを正しく数えるのは正準運動量 pp で測った体積の方である。速度で一様に数えると、重い所を軽く、軽い所を重く数えてしまう。(q,q˙)(q,\dot q) の体積は Hamilton の流れで保たれない(dq dp=m(q) dq dq˙\mathrm{d} q\,\mathrm{d} p=m(q)\,\mathrm{d} q\,\mathrm{d}\dot q で、係数 m(q)m(q) が軌道に沿って変わる)ので、Liouville の定理の保証が無い。

木下の勉強ステータス

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

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

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