解析力学から様々な理論分野を組み立てる 第3部 第5章 §3.5.5 ★★ 量子力学の断熱定理は量子数を保つので、エネルギー以下の固有状態の数と各準位の占有確率が保たれ、エントロピーが変わらない。断熱消磁はその応用
Ehrenfest の断熱仮説
量子力学では、Φ \Phi Φ にあたるのはエネルギー E E E 以下の固有状態の数である。自由度1の Bohr–Sommerfeld の量子化条件 J = ( n + 1 2 ) ℏ J=(n+\frac12)\hbar J = ( n + 2 1 ) ℏ (第1部第14章 )を式 (3.5.6) に入れると Φ = n + 1 2 \Phi=n+\frac12 Φ = n + 2 1 で、n n n 番目の準位のエネルギーまでに n + 1 n+1 n + 1 個の状態があることとほぼ一致する。P. Ehrenfest は、パラメータをゆっくり変えると量子数 n n n が保たれる、すなわち古典の断熱不変量が量子化される量になる、という原理を前期量子論の柱に据えた1 1 P. Ehrenfest(1916)、Annalen der Physik 51、327頁。 。量子力学ではこれが断熱定理の形で示される。固有値 E n ( λ ) E_n(\lambda) E n ( λ ) がほかの固有値から離れていて、λ \lambda λ の変化がゆっくりなら、初めに n n n 番目の固有状態にあった系は、λ \lambda λ が変わった後も n n n 番目の固有状態にとどまる(近似の数理 第7章 ・第4部第9章 )。
近似の裏付け
型 B(極限) | 小さいもの ℏ ∣ λ ˙ ∣ ∥ ∂ H ^ / ∂ λ ∥ / g 2 → 0 \hbar\lvert \dot\lambda\rvert \,\lVert \partial\hat H/\partial\lambda\rVert /g^2\to0 ℏ ∣ λ ˙ ∣ ∥ ∂ H ^ / ∂ λ ∥ / g 2 → 0 (g g g は準位の間隔) | 証明 条件付き(固有値が孤立し、変化のあいだ間隔 g g g が 0 0 0 に近づかない。Kato 1950) | 行き先 近似の数理 第7章
占有確率が保たれる
系が混合状態 ρ ^ = ∑ n p n ∣ n ( λ ) ⟩ ⟨ n ( λ ) ∣ \hat\rho=\sum_np_n| n(\lambda) \rangle \langle n(\lambda) | ρ ^ = ∑ n p n ∣ n ( λ )⟩ ⟨ n ( λ ) ∣ にあるとき、断熱定理により各固有状態は自分の番号を保って動くので、占有確率 p n p_n p n は変わらない。したがって ρ ^ \hat\rho ρ ^ の固有値の組が変わらず、エントロピー
S = − k B ∑ n p n ln p n (3.5.20) S=- k_{\mathrm B}\sum_np_n\ln p_n\tag{3.5.20} S = − k B n ∑ p n ln p n ( 3.5.20 )
は一定である。ここで保たれているのは、固有状態の番号と、それぞれの番号に乗っている確率である。古典の Hertz の定理が「エネルギー以下の体積」を保つのと同じく、量子の断熱変化は「何番目の状態か」を保つ。
準位がすべて同じ倍率で動く場合、すなわち E n ( λ ) = s ( λ ) ε n E_n(\lambda)=s(\lambda)\varepsilon _n E n ( λ ) = s ( λ ) ε n の形のときは、初めにカノニカル分布 p n ∝ e − E n / ( k B T ) p_n\propto\mathrm{e}^{-E_n/( k_{\mathrm B} T)} p n ∝ e − E n / ( k B T ) にあった系は、λ \lambda λ が変わった後もカノニカル分布で、温度だけが
T ∝ s ( λ ) (3.5.21) T\propto s(\lambda)\tag{3.5.21} T ∝ s ( λ ) ( 3.5.21 )
と変わる。p n ∝ e − s ε n / ( k B T ) p_n\propto\mathrm{e}^{-s\varepsilon _n/( k_{\mathrm B} T)} p n ∝ e − s ε n / ( k B T ) が変わらないためには s / T s/T s / T が一定であればよいからである。1辺 a a a の箱の中の粒子では E n ∝ a − 2 E_n\propto a^{-2} E n ∝ a − 2 なので T ∝ a − 2 = V − 2 / 3 T\propto a^{-2}=V^{-2/3} T ∝ a − 2 = V − 2/3 となり、古典の式 (3.5.15) の T V 2 / 3 = TV^{2/3}= T V 2/3 = 一定と一致する。
断熱消磁
磁場 B B B の中のスピン 1 2 \frac12 2 1 の常磁性体を、互いに相互作用しない N N N 個の磁気モーメント(大きさ μ B \mu_{\mathrm B} μ B )の集まりとして扱う。各スピンのエネルギーは ∓ μ B B \mp\mu_{\mathrm B}B ∓ μ B B で、準位は B B B に比例して動く(s ( λ ) = B s(\lambda)=B s ( λ ) = B )。温度 T T T のカノニカル分布では、x = μ B B / ( k B T ) x=\mu_{\mathrm B}B/( k_{\mathrm B} T) x = μ B B / ( k B T ) と置いて
S N k B = ln ( 2 cosh x ) − x tanh x (3.5.22) \frac{S}{N k_{\mathrm B}}=\ln(2\cosh x)-x\tanh x\tag{3.5.22} N k B S = ln ( 2 cosh x ) − x tanh x ( 3.5.22 )
で、エントロピーは B / T B/T B / T だけの関数である。ゆっくり磁場を下げるとエントロピーが一定なので、式 (3.5.21) のとおり B / T B/T B / T が一定に保たれ、温度が磁場に比例して下がる。これが断熱消磁で、低温を作る方法として使われる。図 3.5.3 に2つの磁場でのエントロピーを示す。数値の例は第3部第4章 の3.4.7 節の総合問題 (5)(1 K 1\,\mathrm K 1 K ・1 T 1\,\mathrm T 1 T から磁場を下げる)にある。
図 3.5.3 スピン 1 2 \frac12 2 1 の常磁性体のエントロピー(式 (3.5.22) 、黒の点線は B = 1 T B=1\,\mathrm T B = 1 T 、青の点線は B = 0.1 T B=0.1\,\mathrm T B = 0.1 T )。赤の矢印は 1 K 1\,\mathrm K 1 K で B = 1 T B=1\,\mathrm T B = 1 T から 0.1 T 0.1\,\mathrm T 0.1 T へゆっくり下げる変化で、エントロピーが一定のまま温度が 0.1 K 0.1\,\mathrm K 0.1 K へ下がる。
例題 3.5.5 箱の中の粒子をゆっくり広げる ★★
オリジナル 1辺 L L L の立方体の箱に、互いに相互作用しない質量 m m m の粒子が入っている。1粒子の準位は E n = h 2 8 m L 2 ( n x 2 + n y 2 + n z 2 ) E_{\bm{n}}=\dfrac{h^2}{8mL^2}(n_x^2+n_y^2+n_z^2) E n = 8 m L 2 h 2 ( n x 2 + n y 2 + n z 2 ) (n x , n y , n z = 1 , 2 , … n_x,n_y,n_z=1,2,\dots n x , n y , n z = 1 , 2 , … )である。初め温度 T 0 = 300 K T_0=300\,\mathrm K T 0 = 300 K のカノニカル分布にあり、熱の出入りを断って L L L をゆっくり 2 2 2 倍にする。
(1) 量子の断熱定理で占有確率が保たれるとして、変化の後もカノニカル分布で、温度が T 1 = T 0 / 4 T_1=T_0/4 T 1 = T 0 /4 になることを示せ。式 (3.5.21) の s s s はいくらか。
(2) 粒子が 1 m o l 1\,\mathrm{mol} 1 mol の単原子理想気体として、圧力が何倍になるかを求め、古典の式 (3.5.15) の P V 5 / 3 = PV^{5/3}= P V 5/3 = 一定と比べよ。
(3) 気体がされた仕事を求めよ。R = 8.314 J m o l − 1 K − 1 R=8.314\,\mathrm{J\,mol^{-1}\,K^{-1}} R = 8.314 J mo l − 1 K − 1 とする。
解答 (1) すべての準位が同じ倍率 s = ( L 0 / L ) 2 s=(L_0/L)^2 s = ( L 0 / L ) 2 で動き、L = 2 L 0 L=2L_0 L = 2 L 0 で s = 1 / 4 s=1/4 s = 1/4 である。粒子が相互作用しないので、多粒子の準位も1粒子の準位の和としてすべて s s s 倍になる。準位の順番は入れ替わらず、断熱定理で各状態の占有確率 p ∝ e − E / ( k B T 0 ) p\propto\mathrm{e}^{-E/( k_{\mathrm B} T_0)} p ∝ e − E / ( k B T 0 ) はそのまま運ばれる。変化の後の準位 s E sE s E で同じ確率を書くと e − s E / ( k B T 1 ) \mathrm{e}^{-sE/( k_{\mathrm B} T_1)} e − s E / ( k B T 1 ) で、T 1 = s T 0 = 75 K T_1=sT_0=75\,\mathrm K T 1 = s T 0 = 75 K のカノニカル分布になる。この結論は準位の倍率だけから出たので、準位の間隔が k B T k_{\mathrm B} T k B T と同じ程度の低温(古典の数え方が使えない所)でも、粒子の統計によらず成り立つ。
(2) V = L 3 V=L^3 V = L 3 は 8 8 8 倍、T T T は 1 / 4 1/4 1/4 倍なので、P = N k B T / V P=N k_{\mathrm B} T/V P = N k B T / V は 1 4 ⋅ 1 8 = 1 32 \frac14\cdot\frac18=\frac1{32} 4 1 ⋅ 8 1 = 32 1 倍である。P V 5 / 3 PV^{5/3} P V 5/3 は 1 32 ⋅ 8 5 / 3 = 32 32 = 1 \frac1{32}\cdot8^{5/3}=\frac{32}{32}=1 32 1 ⋅ 8 5/3 = 32 32 = 1 倍で変わらず、T V 2 / 3 TV^{2/3} T V 2/3 も 1 4 ⋅ 8 2 / 3 = 1 \frac14\cdot8^{2/3}=1 4 1 ⋅ 8 2/3 = 1 倍である。量子の「準位の倍率」と古典の「Φ \Phi Φ が一定」が、高温で同じ式を与えている。
(3) 熱の出入りが無いので、仕事はエネルギーの変化 3 2 R ( T 1 − T 0 ) = 3 2 × 8.314 × ( − 225 ) = − 2.81 × 10 3 J \frac32R(T_1-T_0)=\frac32\times8.314\times(-225)=-2.81\times10^3\,\mathrm J 2 3 R ( T 1 − T 0 ) = 2 3 × 8.314 × ( − 225 ) = − 2.81 × 1 0 3 J である。負なので、気体は外に 2.81 × 10 3 J 2.81\times10^3\,\mathrm J 2.81 × 1 0 3 J の仕事をした。
1 P. Ehrenfest(1916)、Annalen der Physik 51、327頁。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com