相対エントロピーから情報と物理を組み立てる 第5部 第1章 §5.1.1 ★★★ 純粋状態では全体と部分系の積の相対エントロピーが部分系のエントロピーの2倍になり、それが体積に比例するか境界に比例するかで状態の性格が分かれる
格子の上に N N N 個のスピンが並んだ系の状態は、2 N 2^N 2 N 個の複素数で書ける。N = 100 N=100 N = 100 なら 2 100 ≈ 1.27 × 10 30 2^{100}\approx1.27\times10^{30} 2 100 ≈ 1.27 × 1 0 30 個で、どんな計算機の記憶にも収まらない。それでも1次元の量子スピン鎖の基底状態は、長い鎖でも計算機で高い精度で求められている(5.1.5 節の密度行列くりこみ群)。この章で扱う問いは、なぜそれができるのか、である。答えを先に書くと、物理で出てくる状態の大半は、乱数で選んだ状態に比べてもつれが極端に少なく、そのもつれは部分系の「境界」の近くにしかない。この性質を面積則 と呼ぶ。
この巻の言葉で言い直すと、この章で比べる2つの状態は、多体系の状態 ρ ^ A B \hat\rho_{AB} ρ ^ A B と、それを領域 A A A とその外 B B B に分けて周辺をかけ合わせた積状態 ρ ^ A ⊗ ρ ^ B \hat\rho_A\otimes\hat\rho_B ρ ^ A ⊗ ρ ^ B である。その間の相対エントロピーが量子相互情報量 I ( A ; B ) = D ( ρ ^ A B ∥ ρ ^ A ⊗ ρ ^ B ) I(A;B)=D(\hat\rho_{AB}\,\|\,\hat\rho_A\otimes\hat\rho_B) I ( A ; B ) = D ( ρ ^ A B ∥ ρ ^ A ⊗ ρ ^ B ) で、純粋状態ではそれが部分系の von Neumann エントロピー(エンタングルメント・エントロピー )の2倍になる。当てる性質は、第2部第2章と第3部第2章で示した Gibbs の変分原理(k B T D ( σ ^ ∥ ρ ^ β ) = F ( σ ^ ) − F e q ≥ 0 k_{\mathrm B}T\,D(\hat\sigma\,\|\,\hat\rho_\beta)=F(\hat\sigma)-F_{\mathrm{eq}}\ge0 k B T D ( σ ^ ∥ ρ ^ β ) = F ( σ ^ ) − F eq ≥ 0 )である。相手の状態 σ ^ \hat\sigma σ ^ に積状態 ρ ^ A ⊗ ρ ^ B \hat\rho_A\otimes\hat\rho_B ρ ^ A ⊗ ρ ^ B を入れると、熱状態の相互情報量が境界の大きさで押さえられる、という面積則が数行で出る。
出てくるものは3つある。1つ目は、1次元で励起にエネルギーの gap がある系の基底状態では、部分系がどれだけ長くてもエンタングルメント・エントロピーが一定の値を超えないこと(Hastings の定理)と、それを自由フェルミオンの鎖で実際に計算する方法(相関行列の方法)。2つ目は、gap が閉じた臨界点ではエントロピーが部分系の長さの対数で増え、その係数が中心電荷 c c c という普遍的な数で決まること。3つ目は、面積則を満たす状態が行列の積(行列積状態 )で少ない数のパラメータで書けることで、ここが冒頭の問いへの答えになる。最後に、学部の統計力学で習う「エントロピーは示量変数(体積に比例する)」という形と、この章の面積則がどう両立するかを、熱状態の計算で確かめる。
この章の対数は自然対数で、情報量の単位は nat とする。Pauli 行列は X X X 、Y Y Y 、Z Z Z 、恒等演算子は 1 \mathbf 1 1 と書く(第3部と同じ)。フェルミオンの生成消滅演算子は c ^ i † \hat c_i^\dagger c ^ i † 、c ^ i \hat c_i c ^ i 、2準位の二値エントロピーは h b ( p ) = − p ln p − ( 1 − p ) ln ( 1 − p ) h_{\mathrm b}(p)=-p\ln p-(1-p)\ln(1-p) h b ( p ) = − p ln p − ( 1 − p ) ln ( 1 − p ) である。
格子・領域・境界
格子 Λ \Lambda Λ の各点(サイト)に d d d 次元の量子系が1つずつ載っているとする。スピン 1 2 \frac12 2 1 なら d = 2 d=2 d = 2 である。全体の Hilbert 空間は各サイトの空間のテンソル積で、次元は d ∣ Λ ∣ d^{\lvert \Lambda\rvert } d ∣ Λ ∣ になる。∣ Λ ∣ \lvert \Lambda\rvert ∣ Λ ∣ はサイトの数である。格子のサイトを2つに分け、一方を領域 A A A 、残りを B B B とする。A A A のサイトの数 ∣ A ∣ \lvert A\rvert ∣ A ∣ を A A A の体積 、A A A のサイトと B B B のサイトを結ぶ最近接の結合の本数 ∣ ∂ A ∣ \lvert \partial A\rvert ∣ ∂ A ∣ を A A A の境界の大きさ (面積)と呼ぶ。1次元の鎖で A A A を長さ L L L の区間に取れば ∣ A ∣ = L \lvert A\rvert =L ∣ A ∣ = L 、∣ ∂ A ∣ = 2 \lvert \partial A\rvert =2 ∣ ∂ A ∣ = 2 で、境界の大きさは L L L によらない。2次元の正方格子でA A A を L × L L\times L L × L の正方形に取れば ∣ A ∣ = L 2 \lvert A\rvert =L^2 ∣ A ∣ = L 2 、∣ ∂ A ∣ = 4 L \lvert \partial A\rvert =4L ∣ ∂ A ∣ = 4 L である(図 5.1.1 )。一般に n d i m n_{\mathrm{dim}} n dim 次元の格子で一辺 L L L の領域を取ると、体積は L n d i m L^{n_{\mathrm{dim}}} L n dim 、境界は L n d i m − 1 L^{n_{\mathrm{dim}}-1} L n dim − 1 に比例する。
図 5.1.1 2次元の正方格子に置いた 4 × 4 4\times4 4 × 4 の領域 A A A (青)と、A A A と B B B を結ぶ境界の結合(赤い太線)。一辺を L L L にすると体積は L 2 L^2 L 2 、境界は 4 L 4L 4 L で増える。面積則は、部分系のエンタングルメント・エントロピーが体積でなく境界の大きさに比例することをいう
エンタングルメント・エントロピー
全体が純粋状態 ∣ Ψ ⟩ | \Psi \rangle ∣Ψ ⟩ にあるとき、比べる2つの状態は ρ ^ A B = ∣ Ψ ⟩ ⟨ Ψ ∣ \hat\rho_{AB}=| \Psi \rangle \langle \Psi | ρ ^ A B = ∣Ψ ⟩ ⟨ Ψ∣ と ρ ^ A ⊗ ρ ^ B \hat\rho_A\otimes\hat\rho_B ρ ^ A ⊗ ρ ^ B である(ρ ^ A = Tr B ρ ^ A B \hat\rho_A=\operatorname{Tr}_B\hat\rho_{AB} ρ ^ A = Tr B ρ ^ A B )。第3部第2章(第3部第2章 )で示したとおり、その相対エントロピーは量子相互情報量
I ( A ; B ) = D ( ρ ^ A B ∥ ρ ^ A ⊗ ρ ^ B ) = S ( ρ ^ A ) + S ( ρ ^ B ) − S ( ρ ^ A B ) (5.1.1) I(A;B)=D(\hat\rho_{AB}\,\|\,\hat\rho_A\otimes\hat\rho_B)=S(\hat\rho_A)+S(\hat\rho_B)-S(\hat\rho_{AB})\tag{5.1.1} I ( A ; B ) = D ( ρ ^ A B ∥ ρ ^ A ⊗ ρ ^ B ) = S ( ρ ^ A ) + S ( ρ ^ B ) − S ( ρ ^ A B ) ( 5.1.1 )
で、純粋状態では S ( ρ ^ A B ) = 0 S(\hat\rho_{AB})=0 S ( ρ ^ A B ) = 0 、Schmidt 分解から S ( ρ ^ A ) = S ( ρ ^ B ) S(\hat\rho_A)=S(\hat\rho_B) S ( ρ ^ A ) = S ( ρ ^ B ) なので
I ( A ; B ) = 2 S ( ρ ^ A ) ( ρ ^ A B が純粋のとき ) (5.1.2) I(A;B)=2S(\hat\rho_A)\qquad(\hat\rho_{AB}\ \text{が純粋のとき})\tag{5.1.2} I ( A ; B ) = 2 S ( ρ ^ A ) ( ρ ^ A B が純粋のとき ) ( 5.1.2 )
となる。第3部第4章(第3部第4章 )では、純粋状態のもつれの量は相対エントロピー・オブ・エンタングルメントでも E R = S ( ρ ^ A ) E_{\mathrm R}=S(\hat\rho_A) E R = S ( ρ ^ A ) になり、多数のコピーから取り出せる最大にもつれた対の数もこれで決まることを示した。そこで次のように定める。
定義 5.1.1(エンタングルメント・エントロピー・面積則・体積則)
純粋状態 ∣ Ψ ⟩ | \Psi \rangle ∣Ψ ⟩ の、領域 A A A についてのエンタングルメント・エントロピー を S ( ρ ^ A ) = − Tr ρ ^ A ln ρ ^ A S(\hat\rho_A)=-\operatorname{Tr}\hat\rho_A\ln\hat\rho_A S ( ρ ^ A ) = − Tr ρ ^ A ln ρ ^ A と定める。状態の族(例えば系の大きさを変えた基底状態の列)が、ある定数 c 0 c_0 c 0 について、すべての領域 A A A で
S ( ρ ^ A ) ≤ c 0 ∣ ∂ A ∣ (5.1.3) S(\hat\rho_A)\le c_0\,\lvert \partial A\rvert\tag{5.1.3} S ( ρ ^ A ) ≤ c 0 ∣ ∂ A ∣ ( 5.1.3 )
を満たすとき、その族は面積則 に従うという。S ( ρ ^ A ) S(\hat\rho_A) S ( ρ ^ A ) が ∣ A ∣ \lvert A\rvert ∣ A ∣ に比例して増えるとき、体積則 に従うという。
上限は必ず体積で押さえられる。ρ ^ A \hat\rho_A ρ ^ A は d ∣ A ∣ d^{\lvert A\rvert } d ∣ A ∣ 次元の密度行列なので、S ( ρ ^ A ) ≤ ∣ A ∣ ln d S(\hat\rho_A)\le\lvert A\rvert \ln d S ( ρ ^ A ) ≤ ∣ A ∣ ln d であり、B B B の側から見れば S ( ρ ^ A ) = S ( ρ ^ B ) ≤ ∣ B ∣ ln d S(\hat\rho_A)=S(\hat\rho_B)\le\lvert B\rvert \ln d S ( ρ ^ A ) = S ( ρ ^ B ) ≤ ∣ B ∣ ln d でもある。面積則は、この上限よりはるかに小さい値に留まる、という強い主張である。混合状態では S ( ρ ^ A ) S(\hat\rho_A) S ( ρ ^ A ) に熱的な混ざり具合も入るので、もつれの量として読めるのは純粋状態のときだけである。混合状態では I ( A ; B ) I(A;B) I ( A ; B ) を使う。これは 5.1.2 節と 5.1.6 節で扱う。
乱数で選んだ状態は体積則に従う
面積則が「強い主張」であることは、乱数で選んだ状態と比べると分かる。第2部第3章(第2部第3章 )では、d A d B d_Ad_B d A d B 次元の Hilbert 空間から一様(Haar 測度)に選んだ純粋状態について、d A ≤ d B d_A\le d_B d A ≤ d B のとき部分系のエントロピーの平均が
S ( ρ ^ A ) ‾ = ∑ k = d B + 1 d A d B 1 k − d A − 1 2 d B ≈ ln d A − d A 2 d B (5.1.4) \overline{S(\hat\rho_A)}=\sum_{k=d_B+1}^{d_Ad_B}\frac1k-\frac{d_A-1}{2d_B}\approx\ln d_A-\frac{d_A}{2d_B}\tag{5.1.4} S ( ρ ^ A ) = k = d B + 1 ∑ d A d B k 1 − 2 d B d A − 1 ≈ ln d A − 2 d B d A ( 5.1.4 )
になることを見た(D. N. Page による1 1 D. N. Page (1993), Physical Review Letters 71 , 1291–1294。 )。d A = d ∣ A ∣ d_A=d^{\lvert A\rvert } d A = d ∣ A ∣ を代入すると S ‾ ≈ ∣ A ∣ ln d \overline{S}\approx\lvert A\rvert \ln d S ≈ ∣ A ∣ ln d で、上限の体積の値にほぼ等しい。つまり「ほとんどすべての状態」は体積則に従う。基底状態が面積則に従うなら、基底状態は Hilbert 空間の中で非常に特別な、ごく一部の状態だということになる。第2部第3章はこの式を「部分系はほとんど熱的に見える」(典型性)の根拠に使った。ここでは同じ式を、基底状態が典型的な状態からどれだけ離れているかを測る基準に使う。
例題 5.1.1 一重項を並べた鎖と乱数の状態 ★★
オリジナル スピン 1 2 \frac12 2 1 を8個並べた鎖の状態 ∣ Ψ ⟩ = ∣ s ⟩ 12 ⊗ ∣ s ⟩ 34 ⊗ ∣ s ⟩ 56 ⊗ ∣ s ⟩ 78 | \Psi \rangle =| s \rangle _{12}\otimes| s \rangle _{34}\otimes| s \rangle _{56}\otimes| s \rangle _{78} ∣Ψ ⟩ = ∣ s ⟩ 12 ⊗ ∣ s ⟩ 34 ⊗ ∣ s ⟩ 56 ⊗ ∣ s ⟩ 78 を考える。∣ s ⟩ i j = 1 2 ( ∣ ↑ ↓ ⟩ − ∣ ↓ ↑ ⟩ ) | s \rangle _{ij}=\frac{1}{\sqrt2}(| \uparrow\downarrow \rangle -| \downarrow\uparrow \rangle ) ∣ s ⟩ ij = 2 1 ( ∣ ↑↓ ⟩ − ∣ ↓↑ ⟩) はサイト i i i 、j j j の一重項である。
A A A をサイト 2〜5 に取ったときと、3〜6 に取ったときの S ( ρ ^ A ) S(\hat\rho_A) S ( ρ ^ A ) を求めよ。
A A A をサイト 2〜7 に広げると S ( ρ ^ A ) S(\hat\rho_A) S ( ρ ^ A ) はどうなるか。この状態が面積則に従うことを説明せよ。
10 個のスピンの Haar 乱数の状態で A A A を左の5個に取ったときの平均 (5.1.4) を求め、上限 5 ln 2 5\ln2 5 ln 2 と比べよ。
解答 (1) 一重項は1組ごとに独立なので、ρ ^ A \hat\rho_A ρ ^ A は一重項ごとの部分系の状態の積になる。サイト 2〜5 は一重項 ( 1 , 2 ) (1,2) ( 1 , 2 ) の片割れ 2、一重項 ( 3 , 4 ) (3,4) ( 3 , 4 ) の全体、一重項 ( 5 , 6 ) (5,6) ( 5 , 6 ) の片割れ 5 を含む。丸ごと含まれる一重項は純粋状態でエントロピーに寄与しない。切られた一重項の片割れは 1 2 1 \frac12\mathbf 1 2 1 1 で、1つにつき ln 2 \ln2 ln 2 だけ寄与する。したがって
S ( ρ ^ 2 – 5 ) = 2 ln 2 = 1.386 , S ( ρ ^ 3 – 6 ) = 0 S(\hat\rho_{2\text{--}5})=2\ln2=1.386,\qquad S(\hat\rho_{3\text{--}6})=0 S ( ρ ^ 2 – 5 ) = 2 ln 2 = 1.386 , S ( ρ ^ 3 – 6 ) = 0
である。3〜6 は一重項 ( 3 , 4 ) (3,4) ( 3 , 4 ) 、( 5 , 6 ) (5,6) ( 5 , 6 ) を丸ごと含み、どれも切らない。
(2) サイト 2〜7 が切る一重項は ( 1 , 2 ) (1,2) ( 1 , 2 ) と ( 7 , 8 ) (7,8) ( 7 , 8 ) の2つだけなので、やはり 2 ln 2 2\ln2 2 ln 2 である。区間をどれだけ長くしても、切れる一重項は区間の両端の2つ以下なのでS ( ρ ^ A ) ≤ 2 ln 2 S(\hat\rho_A)\le2\ln2 S ( ρ ^ A ) ≤ 2 ln 2 が成り立つ。1次元の区間の境界は ∣ ∂ A ∣ = 2 \lvert \partial A\rvert =2 ∣ ∂ A ∣ = 2 なので、c 0 = ln 2 c_0=\ln2 c 0 = ln 2 として式 (5.1.3) を満たす。もつれが最近接の対の中にしかないので、境界をまたぐ対の数だけがエントロピーに入る。
(3) d A = d B = 2 5 = 32 d_A=d_B=2^5=32 d A = d B = 2 5 = 32 で
S ‾ = ∑ k = 33 1024 1 k − 31 64 = 2.966 , \overline{S}=\sum_{k=33}^{1024}\frac1k-\frac{31}{64}=2.966, S = k = 33 ∑ 1024 k 1 − 64 31 = 2.966 ,
上限 5 ln 2 = 3.466 5\ln2=3.466 5 ln 2 = 3.466 との差は 0.499 0.499 0.499 で、ほぼ d A 2 d B = 1 2 \frac{d_A}{2d_B}=\frac12 2 d B d A = 2 1 に等しい。乱数の状態は上限から 1 2 \frac12 2 1 nat しか離れない。一重項の鎖は同じ大きさの A A A で 0 0 0 か 2 ln 2 2\ln2 2 ln 2 なので、乱数の状態とは性格がまったく違う。20 個のスピンを半分に分ければ、乱数の状態は ln 1024 − 1 2 = 6.43 \ln1024-\frac12=6.43 ln 1024 − 2 1 = 6.43 まで増えるが、一重項の鎖は2 ln 2 2\ln2 2 ln 2 のままである。
1 D. N. Page (1993), Physical Review Letters 71 , 1291–1294。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com