tecktom.com

比べる2つの状態:純粋状態の部分系とミクロカノニカル状態の部分系

2026年10月15日8分で読めますAI生成教科書『虎の巻』第2部 第3章 典型性と熱平衡化
目次
  1. 密度行列と部分トレース(最小限)
  2. 比べる2つの状態
  3. 木下の層

AI 生成・木下が未勉強

相対エントロピーから情報と物理を組み立てる 第2部 第3章 §2.3.2 ★★★ 全体では D=ln⁡dRD=\ln d_R と最大に離れた2つの状態を、部分トレースで部分系に落として比べる。この章の主役の設定

密度行列と部分トレース(最小限)

量子系の状態は一般に密度行列 ρ^\hat\rho(エルミートで固有値が非負、Tr⁡ρ^=1\operatorname{Tr}\hat\rho=1)で表す。1つの状態ベクトル ∣ψ⟩| \psi \rangle で書ける純粋状態は ρ^=∣ψ⟩⟨ψ∣\hat\rho=| \psi \rangle \langle \psi | で、それ以外を混合状態と呼ぶ。物理量 O^\hat O の期待値は ⟨O^⟩=Tr⁡(ρ^O^)\langle \hat O \rangle =\operatorname{Tr}(\hat\rho\hat O) である。全体の Hilbert 空間が HA⊗HB\mathcal{H}_A\otimes\mathcal{H}_B のとき、AA だけに作用する物理量 O^A⊗1^B\hat O_A\otimes\hat 1_B の期待値を正しく与える AA の密度行列を縮約密度行列 ρ^A=Tr⁡Bρ^\hat\rho_A=\operatorname{Tr}_B\hat\rho と呼ぶ。HA\mathcal{H}_A の正規直交基底 {∣a⟩}\{| a \rangle \}、HB\mathcal{H}_B の正規直交基底 {∣b⟩}\{| b \rangle \} で ∣ψ⟩=∑a,bψab∣a⟩∣b⟩| \psi \rangle =\sum_{a,b}\psi_{ab}| a \rangle | b \rangle と書けば、

(ρ^A)aa′=∑bψabψa′b∗,つまりρ^A=ΨΨ†(2.3.7)(\hat\rho_A)_{aa'}=\sum_b\psi_{ab}\psi_{a'b}^*,\qquad\text{つまり}\quad\hat\rho_A=\Psi\Psi^\dagger\tag{2.3.7}

である(Ψ\Psi は ψab\psi_{ab} を並べた dA×dBd_A\times d_B 行列)。全体が純粋状態でも、Ψ\Psi の階数が 22 以上なら ρ^A\hat\rho_A は混合状態になる。これが AA と BB のエンタングルメントで、詳しくは 第3部第4章 で扱う。

量子の相対エントロピーは D(ρ^ ∥ σ^)=Tr⁡ρ^(ln⁡ρ^−ln⁡σ^)D(\hat\rho\,\|\,\hat\sigma)=\operatorname{Tr}\hat\rho(\ln\hat\rho-\ln\hat\sigma) で定義し、von Neumann エントロピーは S(ρ^)=−Tr⁡ρ^ln⁡ρ^S(\hat\rho)=-\operatorname{Tr}\hat\rho\ln\hat\rho である。非負性 D≥0D\ge0(等号は ρ^=σ^\hat\rho=\hat\sigma だけ)は 第3部第2章 で、部分トレースを含む任意の量子チャネルで DD が増えないこと(量子データ処理不等式)は第3部第3章 で示す。この章ではこの2つを結果として使う。

比べる2つの状態

NN 個の2準位の粒子を、今度は量子の系として扱う。各粒子の状態空間は ∣0⟩| 0 \rangle(基底)と ∣1⟩| 1 \rangle(励起)で張られる2次元で、前の節の配置 ss は全体の Hilbert 空間の正規直交基底 ∣s⟩=∣n1⟩⋯∣nN⟩| s \rangle =| n_1 \rangle \cdots| n_N \rangle になる。励起がちょうど MM 個の基底ベクトルが張る部分空間を HR\mathcal{H}_R(エネルギー殻、次元 dR=(NM)d_R=\binom NM)とし、そこへの射影を Π^R\hat\Pi_R と書く。一般の系でも、HR\mathcal{H}_R は全体のハミルトニアンの固有値が [E,E+ΔE][E,E+\Delta E] に入る固有状態が張る部分空間とする。

比べる状態の一方は、HR\mathcal{H}_R の中の純粋状態 ∣ψ⟩| \psi \rangle である。もう一方は、HR\mathcal{H}_R の上の一様な混合

ρ^mc=Π^RdR(2.3.8)\hat\rho_{\mathrm{mc}}=\frac{\hat\Pi_R}{d_R}\tag{2.3.8}

で、これをミクロカノニカル状態と呼ぶ。全体どうしを比べると、∣ψ⟩∈HR| \psi \rangle \in\mathcal{H}_R なら Π^R∣ψ⟩=∣ψ⟩\hat\Pi_R| \psi \rangle =| \psi \rangle なので

D(∣ψ⟩⟨ψ∣ ∥ ρ^mc)=0−⟨ψ∣ln⁡Π^RdR∣ψ⟩=ln⁡dR(2.3.9)D(| \psi \rangle \langle \psi |\,\|\,\hat\rho_{\mathrm{mc}})=0-\langle \psi |\ln\frac{\hat\Pi_R}{d_R}| \psi \rangle =\ln d_R\tag{2.3.9}

となる(S(∣ψ⟩⟨ψ∣)=0S(| \psi \rangle \langle \psi |)=0 で、ln⁡ρ^mc\ln\hat\rho_{\mathrm{mc}} は HR\mathcal{H}_R の上で定数 −ln⁡dR-\ln d_R)。ln⁡dR=Smc/kB\ln d_R=S^{\mathrm{mc}}/ k_{\mathrm B} は熱力学のエントロピーそのもので、マクロな系では 102310^{23} の程度の数になる。1回の測定で ∣ψ⟩| \psi \rangle への射影 ∣ψ⟩⟨ψ∣| \psi \rangle \langle \psi | を測れば、∣ψ⟩| \psi \rangle なら確率 11、ρ^mc\hat\rho_{\mathrm{mc}} なら確率 1/dR1/d_R で「はい」が出るので、全体を測れる観測者には2つの状態は完全に近い精度で区別できる。

部分系 AA に制限すると、比べるのは

ρ^A=Tr⁡B∣ψ⟩⟨ψ∣,ρ^Amc=Tr⁡Bρ^mc(2.3.10)\hat\rho_A=\operatorname{Tr}_B| \psi \rangle \langle \psi |,\qquad\hat\rho_A^{\mathrm{mc}}=\operatorname{Tr}_B\hat\rho_{\mathrm{mc}}\tag{2.3.10}

になる。2準位の模型では、ρ^mc\hat\rho_{\mathrm{mc}} は配置の基底で対角で、その部分トレースは前の節の周辺 pAmcp^{\mathrm{mc}}_A を対角に並べたものになる:ρ^Amc=∑sApAmc(sA)∣sA⟩⟨sA∣\hat\rho_A^{\mathrm{mc}}=\sum_{s_A}p^{\mathrm{mc}}_A(s_A)| s_A \rangle \langle s_A |。前の節の結果から、これは k≪Nk\ll N でカノニカル状態にほとんど等しい。したがって問題は、「1つの純粋状態から作った ρ^A\hat\rho_A が、ρ^Amc\hat\rho_A^{\mathrm{mc}} にどれだけ近いか」に帰着する。

部分トレースは量子チャネルなので、データ処理不等式から

D(ρ^A ∥ ρ^Amc)≤D(∣ψ⟩⟨ψ∣ ∥ ρ^mc)=ln⁡dR(2.3.11)D(\hat\rho_A\,\|\,\hat\rho_A^{\mathrm{mc}})\le D(| \psi \rangle \langle \psi |\,\|\,\hat\rho_{\mathrm{mc}})=\ln d_R\tag{2.3.11}

が成り立つ。これは何も言っていないに等しい上界である。ρ^A\hat\rho_A は dAd_A 次元の状態なので、そもそも DD は ln⁡(1/λmin⁡)\ln(1/\lambda_{\min}) の程度(λmin⁡\lambda_{\min} は ρ^Amc\hat\rho_A^{\mathrm{mc}} の最小固有値)を超えにくい。問題は、典型的な ∣ψ⟩| \psi \rangle で D(ρ^A ∥ ρ^Amc)D(\hat\rho_A\,\|\,\hat\rho_A^{\mathrm{mc}}) が 00 に近いかどうかである。次の例題で示すように、個々の ∣ψ⟩| \psi \rangle では近くないこともある。

例題 2.3.2 W 状態の2粒子 ★★

オリジナルN=4N=4 個の2準位の粒子で励起を M=1M=1 個に固定する(dR=4d_R=4)。HR\mathcal{H}_R の中の状態

∣W⟩=12(∣1000⟩+∣0100⟩+∣0010⟩+∣0001⟩)| \mathrm W \rangle =\frac12\bigl(| 1000 \rangle +| 0100 \rangle +| 0010 \rangle +| 0001 \rangle \bigr)

を考え、AA を粒子 1,21,2、BB を粒子 3,43,4 とする。

  1. (1) ρ^Amc\hat\rho_A^{\mathrm{mc}} と ρ^A\hat\rho_A を、基底 ∣00⟩,∣01⟩,∣10⟩,∣11⟩| 00 \rangle ,| 01 \rangle ,| 10 \rangle ,| 11 \rangle の行列で求めよ。
  2. (2) D(ρ^A ∥ ρ^Amc)D(\hat\rho_A\,\|\,\hat\rho_A^{\mathrm{mc}}) を求め、全体の D=ln⁡dRD=\ln d_R と比べよ。
解答

(1) HR\mathcal{H}_R の4つの基底のうち、AA が ∣10⟩| 10 \rangle のものと ∣01⟩| 01 \rangle のものが1つずつ、AA が ∣00⟩| 00 \rangle のものが2つあるのでρ^Amc=diag(1/2,1/4,1/4,0)\hat\rho_A^{\mathrm{mc}}=\mathrm{diag}(1/2,1/4,1/4,0) である(∣11⟩| 11 \rangle は励起が2個になるので現れない)。∣W⟩| \mathrm W \rangle の係数 ψab\psi_{ab} は、(a,b)=(10,00),(01,00),(00,10),(00,01)(a,b)=(10,00),(01,00),(00,10),(00,01) で 1/21/2、ほかは 00 である。式 (2.3.7) から

ρ^A=(1/200001/41/4001/41/400000).\hat\rho_A=\begin{pmatrix}1/2&0&0&0\\0&1/4&1/4&0\\0&1/4&1/4&0\\0&0&0&0\end{pmatrix}.

(01,10)(01,10) 成分が 1/41/4 になるのは、BB の同じ状態 ∣00⟩| 00 \rangle に AA の ∣01⟩| 01 \rangle と ∣10⟩| 10 \rangle の両方が付いているからで、AA の中に重ね合わせ(コヒーレンス)が残る。

(2) ρ^A\hat\rho_A の固有値は 1/21/2(固有ベクトル ∣00⟩| 00 \rangle)、1/21/2((∣01⟩+∣10⟩)/2(| 01 \rangle +| 10 \rangle )/\sqrt2)、00、00 なので S(ρ^A)=ln⁡2S(\hat\rho_A)=\ln2。ln⁡ρ^Amc\ln\hat\rho_A^{\mathrm{mc}} は ρ^A\hat\rho_A の台の上で対角なので −Tr⁡ρ^Aln⁡ρ^Amc=12ln⁡2+(14+14)ln⁡4=32ln⁡2-\operatorname{Tr}\hat\rho_A\ln\hat\rho_A^{\mathrm{mc}}=\frac12\ln2+\bigl(\frac14+\frac14\bigr)\ln4=\frac32\ln2。よって

D(ρ^A ∥ ρ^Amc)=−S(ρ^A)−Tr⁡ρ^Aln⁡ρ^Amc=12ln⁡2=0.347 natD(\hat\rho_A\,\|\,\hat\rho_A^{\mathrm{mc}})=-S(\hat\rho_A)-\operatorname{Tr}\hat\rho_A\ln\hat\rho_A^{\mathrm{mc}}=\tfrac12\ln2=0.347\ \text{nat}

である。全体の D=ln⁡4=1.386D=\ln4=1.386 nat の4分の1で、データ処理不等式 (2.3.11) を満たすが、00 に近いとは言えない。dR=4d_R=4 の小さな系では、AA の外に AA の位相をかき消すだけの自由度が足りない。

木下の勉強ステータス

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

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

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