相対エントロピーから情報と物理を組み立てる 第2部 第3章 §2.3.2 ★★★ 全体では D=lndR と最大に離れた2つの状態を、部分トレースで部分系に落として比べる。この章の主役の設定
密度行列と部分トレース(最小限)
量子系の状態は一般に密度行列 ρ^(エルミートで固有値が非負、Trρ^=1)で表す。1つの状態ベクトル ∣ψ⟩ で書ける純粋状態は ρ^=∣ψ⟩⟨ψ∣ で、それ以外を混合状態と呼ぶ。物理量 O^ の期待値は ⟨O^⟩=Tr(ρ^O^) である。全体の Hilbert 空間が HA⊗HB のとき、A だけに作用する物理量 O^A⊗1^B の期待値を正しく与える A の密度行列を縮約密度行列 ρ^A=TrBρ^ と呼ぶ。HA の正規直交基底 {∣a⟩}、HB の正規直交基底 {∣b⟩} で ∣ψ⟩=∑a,bψab∣a⟩∣b⟩ と書けば、
(ρ^A)aa′=b∑ψabψa′b∗,つまりρ^A=ΨΨ†(2.3.7)
である(Ψ は ψab を並べた dA×dB 行列)。全体が純粋状態でも、Ψ の階数が 2 以上なら ρ^A は混合状態になる。これが A と B のエンタングルメントで、詳しくは 第3部第4章 で扱う。
量子の相対エントロピーは D(ρ^∥σ^)=Trρ^(lnρ^−lnσ^) で定義し、von Neumann エントロピーは S(ρ^)=−Trρ^lnρ^ である。非負性 D≥0(等号は ρ^=σ^ だけ)は 第3部第2章 で、部分トレースを含む任意の量子チャネルで D が増えないこと(量子データ処理不等式)は第3部第3章 で示す。この章ではこの2つを結果として使う。
比べる2つの状態
N 個の2準位の粒子を、今度は量子の系として扱う。各粒子の状態空間は ∣0⟩(基底)と ∣1⟩(励起)で張られる2次元で、前の節の配置 s は全体の Hilbert 空間の正規直交基底 ∣s⟩=∣n1⟩⋯∣nN⟩ になる。励起がちょうど M 個の基底ベクトルが張る部分空間を HR(エネルギー殻、次元 dR=(MN))とし、そこへの射影を Π^R と書く。一般の系でも、HR は全体のハミルトニアンの固有値が [E,E+ΔE] に入る固有状態が張る部分空間とする。
比べる状態の一方は、HR の中の純粋状態 ∣ψ⟩ である。もう一方は、HR の上の一様な混合
ρ^mc=dRΠ^R(2.3.8)
で、これをミクロカノニカル状態と呼ぶ。全体どうしを比べると、∣ψ⟩∈HR なら Π^R∣ψ⟩=∣ψ⟩ なので
D(∣ψ⟩⟨ψ∣∥ρ^mc)=0−⟨ψ∣lndRΠ^R∣ψ⟩=lndR(2.3.9)
となる(S(∣ψ⟩⟨ψ∣)=0 で、lnρ^mc は HR の上で定数 −lndR)。lndR=Smc/kB は熱力学のエントロピーそのもので、マクロな系では 1023 の程度の数になる。1回の測定で ∣ψ⟩ への射影 ∣ψ⟩⟨ψ∣ を測れば、∣ψ⟩ なら確率 1、ρ^mc なら確率 1/dR で「はい」が出るので、全体を測れる観測者には2つの状態は完全に近い精度で区別できる。
部分系 A に制限すると、比べるのは
ρ^A=TrB∣ψ⟩⟨ψ∣,ρ^Amc=TrBρ^mc(2.3.10)
になる。2準位の模型では、ρ^mc は配置の基底で対角で、その部分トレースは前の節の周辺 pAmc を対角に並べたものになる:ρ^Amc=∑sApAmc(sA)∣sA⟩⟨sA∣。前の節の結果から、これは k≪N でカノニカル状態にほとんど等しい。したがって問題は、「1つの純粋状態から作った ρ^A が、ρ^Amc にどれだけ近いか」に帰着する。
部分トレースは量子チャネルなので、データ処理不等式から
D(ρ^A∥ρ^Amc)≤D(∣ψ⟩⟨ψ∣∥ρ^mc)=lndR(2.3.11)
が成り立つ。これは何も言っていないに等しい上界である。ρ^A は dA 次元の状態なので、そもそも D は ln(1/λmin) の程度(λmin は ρ^Amc の最小固有値)を超えにくい。問題は、典型的な ∣ψ⟩ で D(ρ^A∥ρ^Amc) が 0 に近いかどうかである。次の例題で示すように、個々の ∣ψ⟩ では近くないこともある。
例題 2.3.2 W 状態の2粒子 ★★
オリジナルN=4 個の2準位の粒子で励起を M=1 個に固定する(dR=4)。HR の中の状態
∣W⟩=21(∣1000⟩+∣0100⟩+∣0010⟩+∣0001⟩)
を考え、A を粒子 1,2、B を粒子 3,4 とする。
- (1) ρ^Amc と ρ^A を、基底 ∣00⟩,∣01⟩,∣10⟩,∣11⟩ の行列で求めよ。
- (2) D(ρ^A∥ρ^Amc) を求め、全体の D=lndR と比べよ。
解答
(1) HR の4つの基底のうち、A が ∣10⟩ のものと ∣01⟩ のものが1つずつ、A が ∣00⟩ のものが2つあるのでρ^Amc=diag(1/2,1/4,1/4,0) である(∣11⟩ は励起が2個になるので現れない)。∣W⟩ の係数 ψab は、(a,b)=(10,00),(01,00),(00,10),(00,01) で 1/2、ほかは 0 である。式 (2.3.7) から
ρ^A=1/200001/41/4001/41/400000.
(01,10) 成分が 1/4 になるのは、B の同じ状態 ∣00⟩ に A の ∣01⟩ と ∣10⟩ の両方が付いているからで、A の中に重ね合わせ(コヒーレンス)が残る。
(2) ρ^A の固有値は 1/2(固有ベクトル ∣00⟩)、1/2((∣01⟩+∣10⟩)/2)、0、0 なので S(ρ^A)=ln2。lnρ^Amc は ρ^A の台の上で対角なので −Trρ^Alnρ^Amc=21ln2+(41+41)ln4=23ln2。よって
D(ρ^A∥ρ^Amc)=−S(ρ^A)−Trρ^Alnρ^Amc=21ln2=0.347 nat
である。全体の D=ln4=1.386 nat の4分の1で、データ処理不等式 (2.3.11) を満たすが、0 に近いとは言えない。dR=4 の小さな系では、A の外に A の位相をかき消すだけの自由度が足りない。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com