相対エントロピーから情報と物理を組み立てる 第5部 第1章 §5.1.6 ★★ 学部の統計力学の「エントロピーは体積に比例する」は熱状態の部分系のエントロピーの話で、相互情報量は熱状態でも面積則に従う
エントロピーは示量変数である
学部の統計力学では、エントロピーは示量変数で、系を2倍にすると2倍になる、と習う。例えば単原子分子の理想気体の Sackur–Tetrode の式
S = N k B [ ln V N λ T 3 + 5 2 ] , λ T = h 2 π m k B T (5.1.19) S=Nk_{\mathrm B}\Bigl[\ln\frac{V}{N\lambda_T^3}+\frac52\Bigr],\qquad\lambda_T=\frac{h}{\sqrt{2\pi mk_{\mathrm B}T}}\tag{5.1.19} S = N k B [ ln N λ T 3 V + 2 5 ] , λ T = 2 π m k B T h ( 5.1.19 )
は N N N と V V V を同じ割合で増やすと S S S も同じ割合で増える。これと面積則は矛盾しない。式 (5.1.19) は熱状態という混合状態のエントロピーで、この章の言葉では「部分系 A A A の von Neumann エントロピー S ( ρ ^ A ) S(\hat\rho_A) S ( ρ ^ A ) に k B k_{\mathrm B} k B を掛けたもの」である。熱状態の ρ ^ A \hat\rho_A ρ ^ A はほぼ A A A だけの熱状態に近く(第2部第3章の典型性)、そのエントロピーは体積に比例する。一方、A A A と B B B の間の相関の量 I ( A ; B ) I(A;B) I ( A ; B ) は、定理 5.1.2 により境界の大きさで押さえられる。2つを合わせると、熱状態の部分系のエントロピーは
S ( ρ ^ A ) = s ∣ A ∣ + O ( ∣ ∂ A ∣ ) (5.1.20) S(\hat\rho_A)=s\lvert A\rvert +\mathcal{O}(\lvert \partial A\rvert )\tag{5.1.20} S ( ρ ^ A ) = s ∣ A ∣ + O (∣ ∂ A ∣) ( 5.1.20 )
という形になる(s s s は1サイトあたりのエントロピー密度を k B k_{\mathrm B} k B で割ったもの)。体積の項が熱力学のエントロピー、境界の項に A A A と B B B の相関が入る。純粋状態にすると体積の項が消え、境界の項だけが残る。
これを自由フェルミオンの一様な鎖で確かめる。温度 T T T の熱状態の相関行列は、式 (5.1.14) の占有の階段を Fermi 分布 f ( ε ) = 1 / ( e β ε + 1 ) f(\varepsilon )=1/(\mathrm{e}^{\beta\varepsilon }+1) f ( ε ) = 1/ ( e β ε + 1 ) に替えたC i j = ∫ − π π d k 2 π f ( ε ( k ) ) e i k ( i − j ) C_{ij}=\int_{-\pi}^{\pi}\frac{\mathrm{d} k}{2\pi}f(\varepsilon (k))\mathrm{e}^{\mathrm{i} k(i-j)} C ij = ∫ − π π 2 π d k f ( ε ( k )) e i k ( i − j ) である(化学ポテンシャル 0 0 0 、半分の占有)。k B T = 0.5 t k_{\mathrm B}T=0.5\,t k B T = 0.5 t で式 (5.1.13) から S ( L ) S(L) S ( L ) を計算すると、図 5.1.4 のように L L L に比例して増え、傾きはエントロピー密度 s = ∫ d k 2 π h b ( f ( ε ( k ) ) ) = 0.2973 s=\int\frac{\mathrm{d} k}{2\pi}h_{\mathrm b}(f(\varepsilon (k)))=0.2973 s = ∫ 2 π d k h b ( f ( ε ( k ))) = 0.2973 に一致した。S ( L ) − s L S(L)-sL S ( L ) − s L は L L L によらず 0.475 0.475 0.475 に落ち着き、これが境界の分である。同じ温度で、64 サイトの環を半分ずつに分けた熱状態の相互情報量は I ( A ; B ) = 0.950 I(A;B)=0.950 I ( A ; B ) = 0.950 で、16 サイトの環でも 32 サイトの環でも同じ値だった。境界は結合2本で、1本の項 − t ( c ^ i † c ^ i + 1 + h.c. ) -t(\hat c_i^\dagger\hat c_{i+1}+\text{h.c.}) − t ( c ^ i † c ^ i + 1 + h.c. ) のノルムは t t t なので、上界 (5.1.8) は 2 × 2 t / ( 0.5 t ) = 8 2\times2t/(0.5\,t)=8 2 × 2 t / ( 0.5 t ) = 8 である。I I I は上界より1桁小さい。
図 5.1.4 一様な自由フェルミオンの鎖で、長さ L L L の区間のエントロピー。T = 0 T=0 T = 0 (青)は対数でしか増えないが、k B T = 0.5 t k_{\mathrm B}T=0.5\,t k B T = 0.5 t (赤)は L L L に比例して増える(体積則)。ただし A A A と B B B の相互情報量はこの温度でも O ( 1 ) O(1) O ( 1 ) に留まる。値は相関行列の方法による厳密な計算
一重項・Sommerfeld 展開・ブラックホール
学部の量子力学で習う2つのスピン 1 2 \frac12 2 1 の一重項 1 2 ( ∣ ↑ ↓ ⟩ − ∣ ↓ ↑ ⟩ ) \frac1{\sqrt2}(| \uparrow\downarrow \rangle -| \downarrow\uparrow \rangle ) 2 1 ( ∣ ↑↓ ⟩ − ∣ ↓↑ ⟩) は、片方を取り出すと 1 2 1 \frac12\mathbf 1 2 1 1 で S = ln 2 S=\ln2 S = ln 2 になる。例題 5.1.1 の鎖のエントロピーは、境界をまたぐ一重項の数にこの ln 2 \ln2 ln 2 を掛けたものだった。
低温の自由電子の熱容量で習う Sommerfeld 展開も、この章の対数の係数とつながる。低温ではエントロピー密度が T T T に比例し、S / V = π 2 3 k B 2 T D ( ε F ) S/V=\frac{\pi^2}{3}k_{\mathrm B}^2T\,D(\varepsilon _{\mathrm F}) S / V = 3 π 2 k B 2 T D ( ε F ) になる(D ( ε F ) D(\varepsilon _{\mathrm F}) D ( ε F ) は Fermi 準位の状態密度)。1次元の臨界系では、これが共形場理論の s = π c k B T 3 ℏ v s=\frac{\pi c\,k_{\mathrm B}T}{3\hbar v} s = 3ℏ v π c k B T (v v v は Fermi 速度)と同じ式になり、比熱の係数にも中心電荷 c c c が現れる(次の例題)。
最後に、面積則という名前は、ブラックホールのエントロピーが地平面の面積に比例する Bekenstein–Hawking の式 S B H = k B A H 4 ℓ P 2 S_{\mathrm{BH}}=\frac{k_{\mathrm B}A_{\mathrm H}}{4\ell_{\mathrm P}^2} S BH = 4 ℓ P 2 k B A H (A H A_{\mathrm H} A H は地平面の面積、ℓ P \ell_{\mathrm P} ℓ P は Planck 長。定義は第7部第1章 の冒頭)を意識したものである。L. Bombelli、R. K. Koul、J. Lee、R. D. Sorkin1 1 L. Bombelli, R. K. Koul, J. Lee, R. D. Sorkin (1986), Physical Review D 34 , 373–383。 と Srednicki(1993)は、真空の量子場を球の内と外に分けたときのエンタングルメント・エントロピーが球の表面積に比例することを示し、S B H S_{\mathrm{BH}} S BH の起源の候補とした。格子定数を短くすると係数が発散するので、そのままでは S B H S_{\mathrm{BH}} S BH の値を与えないが、形は同じである。この話は第7部(第7部第1章 ・第7部第2章 )で扱う。
近似の裏付け
型 B(極限) | 小さいもの k B T / t → 0 k_{\mathrm B}T/t\to0 k B T / t → 0 (温度がバンド幅に比べて低い) | 証明 証明済み(Fermi 準位の近くで状態密度が滑らかなら、Sommerfeld 展開は誤差の評価つきの漸近展開になる) | 行き先 近似の数理 第7章
例題 5.1.6 低温のエントロピー密度と中心電荷 ★★
オリジナル 一様な鎖の1粒子のエネルギーは ε ( k ) = − 2 t cos k \varepsilon (k)=-2t\cos k ε ( k ) = − 2 t cos k (格子定数 1 1 1 、− π < k ≤ π -\pi<k\le\pi − π < k ≤ π )、半分の占有で Fermi 準位は ε F = 0 \varepsilon _{\mathrm F}=0 ε F = 0 である。
1サイトあたりの状態密度 D ( ε F ) D(\varepsilon _{\mathrm F}) D ( ε F ) を求め、Sommerfeld 展開で1サイトあたりのエントロピー s s s (k B k_{\mathrm B} k B で割ったもの)を k B T / t k_{\mathrm B}T/t k B T / t の1次まで求めよ。
Fermi 速度 v = 1 ℏ ∣ d ε d k ∣ k = π / 2 v=\frac1\hbar\lvert \frac{\mathrm{d} \varepsilon }{\mathrm{d} k}\rvert _{k=\pi/2} v = ℏ 1 ∣ d k d ε ∣ k = π /2 を求め、s = π c k B T 3 ℏ v s=\frac{\pi c\,k_{\mathrm B}T}{3\hbar v} s = 3ℏ v π c k B T と比べて c c c を読み取れ。
k B T = 0.05 t k_{\mathrm B}T=0.05\,t k B T = 0.05 t で (1) の値を求め、数値積分 s = ∫ d k 2 π h b ( f ( ε ( k ) ) ) = 0.02622 s=\int\frac{\mathrm{d} k}{2\pi}h_{\mathrm b}(f(\varepsilon (k)))=0.02622 s = ∫ 2 π d k h b ( f ( ε ( k ))) = 0.02622 と比べよ。k B T = 0.5 t k_{\mathrm B}T=0.5\,t k B T = 0.5 t ではどうか。
解答 (1) Fermi 準位は k = ± π 2 k=\pm\frac\pi2 k = ± 2 π の2点にあり、それぞれで ∣ d ε d k ∣ = 2 t \lvert \frac{\mathrm{d} \varepsilon }{\mathrm{d} k}\rvert =2t ∣ d k d ε ∣ = 2 t である。1サイトあたりの状態密度は k k k の密度 1 2 π \frac{1}{2\pi} 2 π 1 を使って
D ( ε F ) = 2 × 1 2 π × 1 2 t = 1 2 π t . D(\varepsilon _{\mathrm F})=2\times\frac{1}{2\pi}\times\frac{1}{2t}=\frac{1}{2\pi t}. D ( ε F ) = 2 × 2 π 1 × 2 t 1 = 2 π t 1 .
Sommerfeld 展開 s = π 2 3 k B T D ( ε F ) s=\frac{\pi^2}{3}k_{\mathrm B}T\,D(\varepsilon _{\mathrm F}) s = 3 π 2 k B T D ( ε F ) に入れて s = π k B T 6 t s=\frac{\pi k_{\mathrm B}T}{6t} s = 6 t π k B T である。
(2) v = 2 t ℏ v=\frac{2t}{\hbar} v = ℏ 2 t (格子定数1あたり)。π c k B T 3 ℏ v = π c k B T 6 t \frac{\pi c\,k_{\mathrm B}T}{3\hbar v}=\frac{\pi c\,k_{\mathrm B}T}{6t} 3ℏ v π c k B T = 6 t π c k B T が (1) と一致するのは c = 1 c=1 c = 1 のときで、5.1.4 節で T = 0 T=0 T = 0 のエンタングルメント・エントロピーの傾きから読み取った値と同じである。
(3) k B T = 0.05 t k_{\mathrm B}T=0.05\,t k B T = 0.05 t では π × 0.05 6 = 0.02618 \frac{\pi\times0.05}{6}=0.02618 6 π × 0.05 = 0.02618 で、数値積分の 0.02622 0.02622 0.02622 と 0.2 % 0.2\% 0.2% 以内で一致する。k B T = 0.5 t k_{\mathrm B}T=0.5\,t k B T = 0.5 t では1次の式は 0.262 0.262 0.262 を与えるが、数値積分は 0.297 0.297 0.297 で、1割以上ずれる。温度がバンド幅 4 t 4t 4 t に比べて小さくないので、Fermi 準位の近くだけで決まる Sommerfeld 展開の前提が成り立たない。
1 L. Bombelli, R. K. Koul, J. Lee, R. D. Sorkin (1986), Physical Review D 34 , 373–383。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com