tecktom.com

この章のまとめ

2026年10月15日11分で読めますAI生成教科書『虎の巻』第2部 第3章 典型性と熱平衡化

AI 生成・木下が未勉強

相対エントロピーから情報と物理を組み立てる 第2部 第3章 §2.3.7 ★★★ 全体では遠い2つの状態が部分系では見分けられないことを、典型性と平衡化の2つの形でまとめ、総合問題で通して使う

この章で比べたのは、全体の状態を部分系 AA に制限したものと、ミクロカノニカル状態を同じ AA に制限したものである。古典の場合、エネルギー殻の上の一様分布と、平均エネルギーをそろえたカノニカル分布の相対エントロピーは、熱力学エントロピーの差 (Scan−Smc)/kB(S^{\mathrm{can}}-S^{\mathrm{mc}})/ k_{\mathrm B} に等しく(式 (2.3.3))、全体では 12ln⁡N\frac12\ln N の程度で増えるが、kk 個の粒子に制限すると k(k−1)/(4N2)k(k-1)/(4N^2) まで小さくなった(命題 2.3.1)。ずれの源は、励起の総数が決まっていることによる粒子の組の負の相関である。

量子の場合、全体を1つの純粋状態 ∣ψ⟩| \psi \rangle にしても部分系の縮約 ρ^A\hat\rho_A は混合状態で、全体どうしの D=ln⁡dRD=\ln d_R とは別に、部分系どうしの DD を問える。DD を2次の量 χ2\chi^2 で上から押さえる不等式(命題 2.3.2)と Haar 測度の2次のモーメント(補題 2.3.3)を組み合わせると、ランダムな ∣ψ⟩| \psi \rangle では EψD(ρ^A ∥ ρ^Amc)\mathbb E_\psi D(\hat\rho_A\,\|\,\hat\rho_A^{\mathrm{mc}}) が dA/dBeffd_A/d_B^{\mathrm{eff}} の程度に小さい(定理 2.3.4、正準典型性)。制約の無い場合は Page の式と比べて、上界が厳密な平均のちょうど2倍になることを確かめた。Pinsker の不等式から、AA で測れるすべての物理量がミクロカノニカル状態と同じ値をとる。

決まった初期状態からの時間発展では、全体の D(ρ^(t) ∥ ρ^d)=S(ρ^d)D(\hat\rho(t)\,\|\,\hat\rho_{\mathrm d})=S(\hat\rho_{\mathrm d}) は一定で、全体が平衡に近づくことはない(式 (2.3.25)・(2.3.26))。それでも部分系は、エネルギー差に縮退が無く deffd_{\mathrm{eff}} が大きければ、ほとんどすべての時刻で時間平均の状態の近くにいる(定理 2.3.5)。時間平均の状態がミクロカノニカル状態に一致するかは固有状態熱化仮説(ETH)に依存し、これは証明されていない。

この章では、部分系の DD が小さいことを「全体の DD は減らないのに、見る範囲を狭めると減って見える」という形で使った。同じ構図で第二法則を読むのが次の章(第2部第4章)で、そこでは部分系への制限の代わりに粗視化を使い、Liouville の定理と DD の単調減少を両立させる。

総合問題 2準位の粒子 20 個の典型性 ★★★

オリジナルN=20N=20 個の2準位の粒子(励起エネルギー ε\varepsilon)で励起を M=5M=5 個に固定し(x=1/4x=1/4)、最初の2個を部分系 AA、残りの 18 個を BB とする。古典の分布は 2.3.1 節、量子の状態は 2.3.2 節の約束に従う。

  1. (1) (使う道具:2.3.1 節/第1部第3章の DD の定義)pAmcp^{\mathrm{mc}}_A を求め、D(pAmc ∥ pAcan)D(p^{\mathrm{mc}}_A\,\|\,p^{\mathrm{can}}_A) を計算して主要項 k(k−1)/(4N2)k(k-1)/(4N^2) と比べよ。
  2. (2) (使う道具:2.3.6 節/第1部第2章の2準位系の温度)カノニカル分布の βε\beta\varepsilon を求めよ。また AA の2個がともに励起している確率とともに基底状態にある確率の比を、ミクロカノニカル分布とカノニカル分布で比べよ。
  3. (3) (使う道具:2.3.4 節)dRd_R、Tr⁡(ρ^Bmc)2\operatorname{Tr}(\hat\rho_B^{\mathrm{mc}})^2、dBeffd_B^{\mathrm{eff}}、λmin⁡(ρ^Amc)\lambda_{\min}(\hat\rho_A^{\mathrm{mc}}) を求め、Eψ∥ρ^A−ρ^Amc∥1\mathbb E_\psi\lVert \hat\rho_A-\hat\rho_A^{\mathrm{mc}}\rVert _1 と EψD(ρ^A ∥ ρ^Amc)\mathbb E_\psi D(\hat\rho_A\,\|\,\hat\rho_A^{\mathrm{mc}}) の上界を計算せよ。DD の上界を Pinsker の不等式で ∥⋅∥1\lVert \cdot\rVert _1 に直した値とも比べよ。
  4. (4) (使う道具:2.3.3 節の式 (2.3.15))Z^i\hat Z_i を粒子 ii が基底状態なら +1+1、励起していれば −1-1 を返す物理量とする。O^=Z^1Z^2\hat O=\hat Z_1\hat Z_2 の期待値をミクロカノニカル状態とカノニカル状態で求め、ランダムな ∣ψ⟩| \psi \rangle での期待値とカノニカル状態の値との差の平均に上界を与えよ。
  5. (5) (使う道具:2.3.5 節)HR\mathcal{H}_R の中で縮退の無いハミルトニアンを考え、初期状態が HR\mathcal{H}_R のすべての固有状態に均等に広がっているとする。部分系の平衡化の上界 (2.3.28) と、全体の D(ρ^(t) ∥ ρ^d)D(\hat\rho(t)\,\|\,\hat\rho_{\mathrm d}) を求めよ。
解答

(1) (205)=15504\binom{20}{5}=15504、(185)=8568\binom{18}{5}=8568、(184)=3060\binom{18}{4}=3060、(183)=816\binom{18}{3}=816 から、2個とも基底状態の配置が 8568/15504=21/38=0.55268568/15504=21/38=0.5526、片方だけ励起した2つの配置がそれぞれ 3060/15504=15/76=0.19743060/15504=15/76=0.1974、2個とも励起した配置が 816/15504=1/19=0.0526816/15504=1/19=0.0526 である。カノニカル分布の周辺は 9/169/16、3/163/16、1/161/16 なので

D=2138ln⁡21/389/16+2⋅1576ln⁡15/763/16+119ln⁡1/191/16=1.42×10−3 nat.D=\frac{21}{38}\ln\frac{21/38}{9/16}+2\cdot\frac{15}{76}\ln\frac{15/76}{3/16}+\frac1{19}\ln\frac{1/19}{1/16}=1.42\times10^{-3}\ \text{nat}.

主要項は 2/(4⋅400)=1.25×10−32/(4\cdot400)=1.25\times10^{-3} で、値の 88% を占める(N=16N=16 の例題 2.3.1 の 85% より主要項に近い)。

(2) x=1/(eβε+1)=1/4x=1/(\mathrm{e}^{\beta\varepsilon }+1)=1/4 から eβε=3\mathrm{e}^{\beta\varepsilon }=3、βε=ln⁡3=1.099\beta\varepsilon =\ln3=1.099。第1部第2章の例題で N1/N=1/(eβε+1)N_1/N=1/(\mathrm{e}^{\beta\varepsilon }+1) を状態数から導いたのと同じ関係である。比はミクロカノニカル分布で M(M−1)(N−M)(N−M−1)=20210=221=0.0952\frac{M(M-1)}{(N-M)(N-M-1)}=\frac{20}{210}=\frac2{21}=0.0952、カノニカル分布で (1/3)2=1/9=0.111(1/3)^2=1/9=0.111 である。ミクロカノニカル分布の比はカノニカル分布の 18/21=6/7=0.85718/21=6/7=0.857 倍で、例題 2.3.6 の近似 1−(1−2x)/(Nx(1−x))=1−0.5/3.75=0.8671-(1-2x)/(Nx(1-x))=1-0.5/3.75=0.867 とおよそ合う。

(3) dR=(205)=15504d_R=\binom{20}{5}=15504。BB の配置で励起が jj 個のものは、AA に励起を 5−j5-j 個置く場合の数 (25−j)\binom{2}{5-j} 通りと組になるので、ρ^Bmc\hat\rho_B^{\mathrm{mc}} の対角成分は j=5j=5 の 85688568 個が 1/155041/15504、j=4j=4 の 30603060 個が 2/155042/15504、j=3j=3 の 816816 個が 1/155041/15504 である。よって

Tr⁡(ρ^Bmc)2=8568+3060×4+816155042=21624155042,dBeff=15504221624=11116.\operatorname{Tr}(\hat\rho_B^{\mathrm{mc}})^2=\frac{8568+3060\times4+816}{15504^2}=\frac{21624}{15504^2},\qquad d_B^{\mathrm{eff}}=\frac{15504^2}{21624}=11116 .

BB の配置の総数 8568+3060+816=124448568+3060+816=12444 より少し小さい。λmin⁡(ρ^Amc)=1/19\lambda_{\min}(\hat\rho_A^{\mathrm{mc}})=1/19(2個とも励起)。定理 2.3.4 から

Eψ∥ρ^A−ρ^Amc∥1≤411116=0.0190,EψD(ρ^A ∥ ρ^Amc)≤ln⁡(1+1911116)=1.71×10−3 nat.\mathbb E_\psi\lVert \hat\rho_A-\hat\rho_A^{\mathrm{mc}}\rVert _1\le\sqrt{\frac{4}{11116}}=0.0190,\qquad \mathbb E_\psi D(\hat\rho_A\,\|\,\hat\rho_A^{\mathrm{mc}})\le\ln\Bigl(1+\frac{19}{11116}\Bigr)=1.71\times10^{-3}\ \text{nat}.

Pinsker の不等式と E2D≤2ED\mathbb E\sqrt{2D}\le\sqrt{2\mathbb ED} を使うと Eψ∥⋅∥1≤2×1.71×10−3=0.0585\mathbb E_\psi\lVert \cdot\rVert _1\le\sqrt{2\times1.71\times10^{-3}}=0.0585 で、直接の評価 0.01900.0190 の約3倍ゆるい。DD の上界は λmin⁡=1/19\lambda_{\min}=1/19 で割った分だけ不利になり、Pinsker の不等式でさらに平方根を取るからである。

(4) Z^1Z^2\hat Z_1\hat Z_2 は2個が同じ状態なら +1+1、異なれば −1-1 なので、ミクロカノニカル状態で 2138−2⋅1576+119=419=0.2105\frac{21}{38}-2\cdot\frac{15}{76}+\frac1{19}=\frac4{19}=0.2105、カノニカル状態で (1−2x)2=1/4(1-2x)^2=1/4 である。差 3/76=0.03953/76=0.0395 は、pAmcp^{\mathrm{mc}}_A と pAcanp^{\mathrm{can}}_A の ∥⋅∥1\lVert \cdot\rVert _1(0.03950.0395)と一致する。∥Z^1Z^2∥∞=1\lVert \hat Z_1\hat Z_2\rVert _\infty=1 なので、式 (2.3.15) と三角不等式から

Eψ∣⟨Z^1Z^2⟩ψ−14∣≤Eψ∥ρ^A−ρ^Amc∥1+∥ρ^Amc−ρ^Acan∥1≤0.0190+0.0395=0.058.\mathbb E_\psi\lvert \langle \hat Z_1\hat Z_2 \rangle _\psi-\tfrac14\rvert \le\mathbb E_\psi\lVert \hat\rho_A-\hat\rho_A^{\mathrm{mc}}\rVert _1+\lVert \hat\rho_A^{\mathrm{mc}}-\hat\rho_A^{\mathrm{can}}\rVert _1\le0.0190+0.0395=0.058 .

この大きさの系では、カノニカル状態からのずれの大部分は、∣ψ⟩| \psi \rangle を1つに選んだことによるゆらぎよりも、ミクロカノニカル分布とカノニカル分布の違い(有限の NN の補正)による。

(5) deff=dR=15504d_{\mathrm{eff}}=d_R=15504 なので、上界は dA3/deff=64/15504=0.0642\sqrt{d_A^3/d_{\mathrm{eff}}}=\sqrt{64/15504}=0.0642。全体の D(ρ^(t) ∥ ρ^d)=ln⁡15504=9.649D(\hat\rho(t)\,\|\,\hat\rho_{\mathrm d})=\ln15504=9.649 nat で、時間によらない。部分系 AA は長い時間のほとんどで時間平均の状態から ∥⋅∥1\lVert \cdot\rVert _1 で平均 0.0640.064 以内にいるが、全体の状態は時間平均の状態から 9.69.6 nat 離れたまま動かない。

木下の勉強ステータス

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

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

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