相対エントロピーから情報と物理を組み立てる 第3部 第7章 §3.7.1 ★★★ 文字 x x x と量子状態 ρ ^ x \hat\rho_x ρ ^ x の組を1つの密度行列にまとめると、Holevo 量はその量子相互情報量、すなわち積状態との D D D になる
第1部第6章(第1部第6章 )では、古典の通信路を確率行列 W ( y ∣ x ) W(y\,|\,x) W ( y ∣ x ) で表し、その容量が相互情報量の最大値C = max p I ( X ; Y ) C=\max_pI(X;Y) C = max p I ( X ; Y ) であることを示した。この章では、送り手が文字 x x x を古典の信号でなく量子状態 ρ ^ x \hat\rho_x ρ ^ x に載せて送る場合を考える。光ファイバーの中の光パルスも、量子ビットの列も、物理としてはこの形をしている。
受け手は届いた状態を測定して文字を推定する。測定の結果 y y y が出る確率は ρ ^ x \hat\rho_x ρ ^ x と測定の選び方で決まるので、測定を1つ決めると古典の通信路W ( y ∣ x ) W(y\,|\,x) W ( y ∣ x ) が1つできる。問題は、測定をどう選んでも相互情報量 I ( X ; Y ) I(X;Y) I ( X ; Y ) がある量を超えないことで、その量が Holevo 量
χ ( { p x , ρ ^ x } ) = S ( ρ ˉ ) − ∑ x p x S ( ρ ^ x ) , ρ ˉ = ∑ x p x ρ ^ x \chi(\{p_x,\hat\rho_x\})=S(\bar\rho)-\sum_xp_xS(\hat\rho_x),\qquad\bar\rho=\sum_xp_x\hat\rho_x χ ({ p x , ρ ^ x }) = S ( ρ ˉ ) − x ∑ p x S ( ρ ^ x ) , ρ ˉ = x ∑ p x ρ ^ x
である。この章で比べる2つの状態は、文字と量子状態が相関した古典量子状態 ρ ^ X B \hat\rho_{XB} ρ ^ X B と、両者を独立とみなした積状態 ρ ^ X ⊗ ρ ^ B \hat\rho_X\otimes\hat\rho_B ρ ^ X ⊗ ρ ^ B である。Holevo 量はこの2つの量子相対エントロピーになっている。当てる性質は量子データ処理不等式(第3部第3章 )で、測定は量子チャネルの1つなので、測定の後の古典の相互情報量は測定の前の χ \chi χ を超えない。これが Holevo 限界 である。
出てくるものは3つある。第一に、d d d 次元の量子系1つで送れる古典の情報は ln d \ln d ln d を超えない(1量子ビットで1ビット)。超高密度符号(第3部第6章 )が1量子ビットで2ビットを送れることと矛盾しない理由も、ここで分かる。第二に、雑音のある量子チャネル N \mathcal N N の古典の容量が、出力の Holevo 量の最大値で与えられる(Holevo–Schumacher–Westmoreland の定理)。古典の第1部第6章の「容量は D D D の球の半径」という形も、そのまま量子に移る。第三に、光の通信路の容量が熱放射のエントロピーの式 g ( n ˉ ) g(\bar n) g ( n ˉ ) で書け、第1部第6章で古典の熱雑音の近似が崩れると予告した広帯域の極限が、ここで量子の式に置き換わる。
この章の対数は自然対数で、情報量の単位は nat とする。ビットで述べる所だけ「ビット」と書く。量子ビットの Pauli 行列は X , Y , Z X,Y,Z X , Y , Z 、恒等演算子は 1 \mathbf 1 1 、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 ) と書く。von Neumann エントロピー S ( ρ ^ ) S(\hat\rho) S ( ρ ^ ) ・量子相対エントロピー D ( ρ ^ ∥ σ ^ ) D(\hat\rho\,\|\,\hat\sigma) D ( ρ ^ ∥ σ ^ ) ・Holevo 量の定義と非負性は第3部第2章(第3部第2章 )のものをそのまま使う。
量子状態に文字を載せる
送り手は文字の集合 X \mathcal X X から文字 x x x を確率 p x p_x p x で選び、量子系 B B B を状態 ρ ^ x \hat\rho_x ρ ^ x に用意して送る。組 { p x , ρ ^ x } \{p_x,\hat\rho_x\} { p x , ρ ^ x } をアンサンブル と呼ぶ。受け手が文字を知らないとき、手元の状態は混合 ρ ˉ = ∑ x p x ρ ^ x \bar\rho=\sum_xp_x\hat\rho_x ρ ˉ = ∑ x p x ρ ^ x である。
受け手が行う測定は、一般には射影測定に限らない。結果 y y y に対応する半正定値の行列 E ^ y ≥ 0 \hat E_y\ge0 E ^ y ≥ 0 で、∑ y E ^ y = 1 \sum_y\hat E_y=\mathbf 1 ∑ y E ^ y = 1 を満たすものの組{ E ^ y } \{\hat E_y\} { E ^ y } を POVM (正作用素値測度)と呼び、状態 ρ ^ \hat\rho ρ ^ を測ったとき結果 y y y が出る確率を Tr ( E ^ y ρ ^ ) \operatorname{Tr}(\hat E_y\hat\rho) Tr ( E ^ y ρ ^ ) とする。射影測定はE ^ y \hat E_y E ^ y が互いに直交する射影子の場合である。E ^ y ≥ 0 \hat E_y\ge0 E ^ y ≥ 0 は確率が負にならないこと、∑ y E ^ y = 1 \sum_y\hat E_y=\mathbf 1 ∑ y E ^ y = 1 はどの状態でも確率の和が 1 1 1 になることを保証する。POVM は、系に補助の系をつないで全体を射影測定したものとして必ず実現できる(Naimark の拡張)。
定義 3.7.1(測定が作る古典の通信路)
アンサンブル { p x , ρ ^ x } \{p_x,\hat\rho_x\} { p x , ρ ^ x } と POVM { E ^ y } \{\hat E_y\} { E ^ y } に対して
W ( y ∣ x ) ≔ Tr ( E ^ y ρ ^ x ) (3.7.1) W(y\,|\,x)\coloneqq\operatorname{Tr}(\hat E_y\hat\rho_x)\tag{3.7.1} W ( y ∣ x ) : = Tr ( E ^ y ρ ^ x ) ( 3.7.1 )
は各 x x x について y y y の確率分布で、第1部第6章の意味の通信路になる。この通信路の入力 X X X と出力 Y Y Y の相互情報量を I ( X ; Y ) I(X;Y) I ( X ; Y ) と書く。すべての POVM についての上限
I a c c ( { p x , ρ ^ x } ) ≔ sup { E ^ y } I ( X ; Y ) (3.7.2) I_{\mathrm{acc}}(\{p_x,\hat\rho_x\})\coloneqq\sup_{\{\hat E_y\}}I(X;Y)\tag{3.7.2} I acc ({ p x , ρ ^ x }) : = { E ^ y } sup I ( X ; Y ) ( 3.7.2 )
を取り出せる情報 (accessible information)と呼ぶ。
ρ ^ x \hat\rho_x ρ ^ x が互いに直交する純粋状態なら、それらを基底とする射影測定で x x x を誤りなく読めて I ( X ; Y ) = H ( p ) I(X;Y)=H(p) I ( X ; Y ) = H ( p ) になる。直交しなければ、どの測定でも x x x を確実には読めない。I a c c I_{\mathrm{acc}} I acc は測定全体についての最大化で、一般に閉じた式を持たない。そこで、測定を選ぶ前に状態だけから計算できる上限が欲しい。
古典量子状態
文字 x x x を正規直交基底 { ∣ x ⟩ } \{| x \rangle \} { ∣ x ⟩} を持つ補助の系 X X X に書き込んでおくと、送り手の記録と受け手の系の全体を1つの密度行列
ρ ^ X B ≔ ∑ x p x ∣ x ⟩ ⟨ x ∣ ⊗ ρ ^ x (3.7.3) \hat\rho_{XB}\coloneqq\sum_xp_x\,| x \rangle \langle x |\otimes\hat\rho_x\tag{3.7.3} ρ ^ X B : = x ∑ p x ∣ x ⟩ ⟨ x ∣ ⊗ ρ ^ x ( 3.7.3 )
で表せる。X X X の部分が対角なので、X X X の部分は古典の確率変数と同じ働きをする。この形の状態を古典量子状態 (cq 状態)と呼ぶ。部分トレースは ρ ^ X = ∑ x p x ∣ x ⟩ ⟨ x ∣ \hat\rho_X=\sum_xp_x| x \rangle \langle x | ρ ^ X = ∑ x p x ∣ x ⟩ ⟨ x ∣ 、ρ ^ B = ρ ˉ \hat\rho_B=\bar\rho ρ ^ B = ρ ˉ である。
比べる相手は、文字と状態を独立とみなした積状態 ρ ^ X ⊗ ρ ^ B = ∑ x p x ∣ x ⟩ ⟨ x ∣ ⊗ ρ ˉ \hat\rho_X\otimes\hat\rho_B=\sum_xp_x| x \rangle \langle x |\otimes\bar\rho ρ ^ X ⊗ ρ ^ B = ∑ x p x ∣ x ⟩ ⟨ x ∣ ⊗ ρ ˉ である。
命題 3.7.2(Holevo 量は cq 状態の相互情報量)
cq 状態 (3.7.3) について
I ( X ; B ) ρ = D ( ρ ^ X B ∥ ρ ^ X ⊗ ρ ^ B ) = ∑ x p x D ( ρ ^ x ∥ ρ ˉ ) = S ( ρ ˉ ) − ∑ x p x S ( ρ ^ x ) = χ ( { p x , ρ ^ x } ) . (3.7.4) I(X;B)_{\rho}=D(\hat\rho_{XB}\,\|\,\hat\rho_X\otimes\hat\rho_B)=\sum_xp_x\,D(\hat\rho_x\,\|\,\bar\rho)=S(\bar\rho)-\sum_xp_xS(\hat\rho_x)=\chi(\{p_x,\hat\rho_x\}).\tag{3.7.4} I ( X ; B ) ρ = D ( ρ ^ X B ∥ ρ ^ X ⊗ ρ ^ B ) = x ∑ p x D ( ρ ^ x ∥ ρ ˉ ) = S ( ρ ˉ ) − x ∑ p x S ( ρ ^ x ) = χ ({ p x , ρ ^ x }) . ( 3.7.4 )
証明 2つの状態はどちらも X X X の基底で区画に分かれた形 ∑ x ∣ x ⟩ ⟨ x ∣ ⊗ A ^ x \sum_x| x \rangle \langle x |\otimes\hat A_x ∑ x ∣ x ⟩ ⟨ x ∣ ⊗ A ^ x をしているので、行列の関数も区画ごとに取れる。ln ( ρ ^ X B ) = ∑ x ∣ x ⟩ ⟨ x ∣ ⊗ ( ln p x 1 + ln ρ ^ x ) \ln(\hat\rho_{XB})=\sum_x| x \rangle \langle x |\otimes(\ln p_x\,\mathbf 1+\ln\hat\rho_x) ln ( ρ ^ X B ) = ∑ x ∣ x ⟩ ⟨ x ∣ ⊗ ( ln p x 1 + ln ρ ^ x ) 、ln ( ρ ^ X ⊗ ρ ^ B ) = ∑ x ∣ x ⟩ ⟨ x ∣ ⊗ ( ln p x 1 + ln ρ ˉ ) \ln(\hat\rho_X\otimes\hat\rho_B)=\sum_x| x \rangle \langle x |\otimes(\ln p_x\,\mathbf 1+\ln\bar\rho) ln ( ρ ^ X ⊗ ρ ^ B ) = ∑ x ∣ x ⟩ ⟨ x ∣ ⊗ ( ln p x 1 + ln ρ ˉ ) (p x > 0 p_x>0 p x > 0 の区画だけを考える)。差を取って ρ ^ X B \hat\rho_{XB} ρ ^ X B とのトレースを取ると、ln p x \ln p_x ln p x の項が消えて∑ x p x Tr ρ ^ x ( ln ρ ^ x − ln ρ ˉ ) = ∑ x p x D ( ρ ^ x ∥ ρ ˉ ) \sum_xp_x\operatorname{Tr}\hat\rho_x(\ln\hat\rho_x-\ln\bar\rho)=\sum_xp_xD(\hat\rho_x\,\|\,\bar\rho) ∑ x p x Tr ρ ^ x ( ln ρ ^ x − ln ρ ˉ ) = ∑ x p x D ( ρ ^ x ∥ ρ ˉ ) を得る。最後の等号は第3部第2章の Holevo 量の式そのもので、− Tr ρ ^ x ln ρ ˉ -\operatorname{Tr}\hat\rho_x\ln\bar\rho − Tr ρ ^ x ln ρ ˉ を x x x について平均すると − Tr ρ ˉ ln ρ ˉ = S ( ρ ˉ ) -\operatorname{Tr}\bar\rho\ln\bar\rho=S(\bar\rho) − Tr ρ ˉ ln ρ ˉ = S ( ρ ˉ ) にまとまることから従う。
古典の第1部第3章では、相互情報量を結合分布と周辺の積の D D D として定義した。式 (3.7.4) はその量子版で、古典の確率変数 X X X と量子系 B B B の間の相関の量が Holevo 量である。第1部第6章の I ( p , W ) = ∑ x p ( x ) D ( W x ∥ p W ) I(p,W)=\sum_xp(x)D(W_x\,\|\,pW) I ( p , W ) = ∑ x p ( x ) D ( W x ∥ p W ) (行ごとの D D D の平均)とも同じ形をしている。古典の行 W x W_x W x が量子状態 ρ ^ x \hat\rho_x ρ ^ x に、出力の平均の分布 p W pW p W が ρ ˉ \bar\rho ρ ˉ に置き換わった。
例題 3.7.1 ∣ 0 ⟩ | 0 \rangle ∣0 ⟩ と ∣ + ⟩ | + \rangle ∣ + ⟩ を載せた cq 状態 ★★
オリジナル 文字 x ∈ { 0 , 1 } x\in\{0,1\} x ∈ { 0 , 1 } を確率 1 / 2 1/2 1/2 ずつで選び、x = 0 x=0 x = 0 なら ∣ 0 ⟩ | 0 \rangle ∣0 ⟩ 、x = 1 x=1 x = 1 なら ∣ + ⟩ = ( ∣ 0 ⟩ + ∣ 1 ⟩ ) / 2 | + \rangle =(| 0 \rangle +| 1 \rangle )/\sqrt2 ∣ + ⟩ = ( ∣0 ⟩ + ∣1 ⟩) / 2 を送る。
cq 状態 ρ ^ X B \hat\rho_{XB} ρ ^ X B の固有値を求め、S ( X B ) S(XB) S ( X B ) を計算せよ。
S ( X ) S(X) S ( X ) 、S ( B ) S(B) S ( B ) を求め、I ( X ; B ) = S ( X ) + S ( B ) − S ( X B ) I(X;B)=S(X)+S(B)-S(XB) I ( X ; B ) = S ( X ) + S ( B ) − S ( X B ) が Holevo 量 χ \chi χ に一致することを確かめよ。
D ( ρ ^ x ∥ ρ ˉ ) D(\hat\rho_x\,\|\,\bar\rho) D ( ρ ^ x ∥ ρ ˉ ) を x = 0 , 1 x=0,1 x = 0 , 1 について求め、平均が χ \chi χ になることを確かめよ。
解答 (1) ρ ^ X B = 1 2 ∣ 0 ⟩ ⟨ 0 ∣ ⊗ ∣ 0 ⟩ ⟨ 0 ∣ + 1 2 ∣ 1 ⟩ ⟨ 1 ∣ ⊗ ∣ + ⟩ ⟨ + ∣ \hat\rho_{XB}=\frac12| 0 \rangle \langle 0 |\otimes| 0 \rangle \langle 0 |+\frac12| 1 \rangle \langle 1 |\otimes| + \rangle \langle + | ρ ^ X B = 2 1 ∣0 ⟩ ⟨ 0∣ ⊗ ∣0 ⟩ ⟨ 0∣ + 2 1 ∣1 ⟩ ⟨ 1∣ ⊗ ∣ + ⟩ ⟨ + ∣ は区画対角で、各区画は固有値 1 / 2 1/2 1/2 の射影子1つなので、固有値は ( 1 / 2 , 1 / 2 , 0 , 0 ) (1/2,1/2,0,0) ( 1/2 , 1/2 , 0 , 0 ) 。S ( X B ) = ln 2 S(XB)=\ln2 S ( X B ) = ln 2 。
(2) S ( X ) = ln 2 S(X)=\ln2 S ( X ) = ln 2 。ρ ^ B = ρ ˉ = 1 2 ( 1 + 1 2 X + 1 2 Z ) \hat\rho_B=\bar\rho=\frac12(\mathbf 1+\frac12X+\frac12Z) ρ ^ B = ρ ˉ = 2 1 ( 1 + 2 1 X + 2 1 Z ) の固有値は 1 2 ( 1 ± 1 / 2 ) = 0.8536 , 0.1464 \frac12(1\pm1/\sqrt2)=0.8536,\ 0.1464 2 1 ( 1 ± 1/ 2 ) = 0.8536 , 0.1464 で(第3部第2章の例題と同じ)、S ( B ) = 0.4165 S(B)=0.4165 S ( B ) = 0.4165 nat。I ( X ; B ) = ln 2 + 0.4165 − ln 2 = 0.4165 I(X;B)=\ln2+0.4165-\ln2=0.4165 I ( X ; B ) = ln 2 + 0.4165 − ln 2 = 0.4165 nat。純粋状態は S ( ρ ^ x ) = 0 S(\hat\rho_x)=0 S ( ρ ^ x ) = 0 なので χ = S ( ρ ˉ ) = 0.4165 \chi=S(\bar\rho)=0.4165 χ = S ( ρ ˉ ) = 0.4165 nat で一致する。
(3) 純粋状態 ∣ ψ ⟩ | \psi \rangle ∣ ψ ⟩ について D ( ∣ ψ ⟩ ⟨ ψ ∣ ∥ ρ ˉ ) = − ⟨ ψ ∣ ln ρ ˉ ∣ ψ ⟩ D(| \psi \rangle \langle \psi |\,\|\,\bar\rho)=-\langle \psi | \ln\bar\rho | \psi \rangle D ( ∣ ψ ⟩ ⟨ ψ ∣ ∥ ρ ˉ ) = − ⟨ ψ ∣ ln ρ ˉ ∣ ψ ⟩ 。ρ ˉ \bar\rho ρ ˉ の固有ベクトルは Bloch 球で ( X + Z ) / 2 (X+Z)/\sqrt2 ( X + Z ) / 2 の向き、すなわち ∣ 0 ⟩ | 0 \rangle ∣0 ⟩ と ∣ + ⟩ | + \rangle ∣ + ⟩ の真ん中を向く軸の上の2点で、∣ 0 ⟩ | 0 \rangle ∣0 ⟩ も ∣ + ⟩ | + \rangle ∣ + ⟩ もこの軸と角 π / 4 \pi/4 π /4 をなす。上の固有ベクトルとの重なりは cos 2 ( π / 8 ) = 0.8536 \cos^2(\pi/8)=0.8536 cos 2 ( π /8 ) = 0.8536 、下とは 0.1464 0.1464 0.1464 なのでD = − ( 0.8536 ln 0.8536 + 0.1464 ln 0.1464 ) = 0.4165 D=-(0.8536\ln0.8536+0.1464\ln0.1464)=0.4165 D = − ( 0.8536 ln 0.8536 + 0.1464 ln 0.1464 ) = 0.4165 nat が x = 0 , 1 x=0,1 x = 0 , 1 で同じ値になる。平均も 0.4165 0.4165 0.4165 nat で χ \chi χ に一致する。2つの状態が対称に置かれているので、ρ ˉ \bar\rho ρ ˉ から見た「遠さ」が等しい。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com