tecktom.com

アンサンブルの分類:Schrödinger–HJW の定理

2026年10月17日15分で読めますAI生成教科書『虎の巻』第4部 第3章 同じ ρ̂ を作るアンサンブルは1つでない
目次
  1. 混ぜられる状態と、その重みの上限
  2. 重みのばらつきは固有値より大きい
  3. 木下の層

AI 生成・木下が未勉強

量子力学の様々な構成方法 第4部 第3章 §4.3.2 ★★★ 同じ ρ^\hat\rho を与える2つのアンサンブルは、規格化しないベクトルの組がユニタリ行列で移り合うものに限る。どの状態をどれだけの重みで混ぜられるかも、ここから決まる

前節の弦の作り方は量子ビットに限られていた。一般の次元で、同じ ρ^\hat\rho を与えるアンサンブル全部を一度に書くのが次の定理である。2部系の一方の縮約状態について Schrödinger が 1936 年に与え(Schrödinger (1936), Proc. Camb. Phil. Soc. 32, 446–452)、Hughston, Jozsa, Wootters が有限個の純粋状態のアンサンブルの完全な分類として整えた(Hughston, Jozsa, Wootters (1993), Phys. Lett. A 183, 14–18)。3人の頭文字から HJW の定理と呼ばれる。

定理 4.3.1(アンサンブルのユニタリ自由度)

{∣ψ~i⟩}i=1m\{| \tilde\psi_i \rangle \}_{i=1}^m と {∣ϕ~j⟩}j=1m\{| \tilde\phi_j \rangle \}_{j=1}^m を規格化しないベクトルの組とする(個数が違うときは短い方に零ベクトルを足して mm にそろえる)。次の2つは同値である。

  1. (a) ∑i∣ψ~i⟩⟨ψ~i∣=∑j∣ϕ~j⟩⟨ϕ~j∣\sum_i| \tilde\psi_i \rangle \langle \tilde\psi_i |=\sum_j| \tilde\phi_j \rangle \langle \tilde\phi_j |。
  2. (b) m×mm\times m のユニタリ行列 uu があって、∣ψ~i⟩=∑juij∣ϕ~j⟩| \tilde\psi_i \rangle =\sum_ju_{ij}| \tilde\phi_j \rangle(i=1,…,mi=1,\dots,m)。
証明

(b) ⇒\Rightarrow (a):代入して、uu の列が正規直交であること ∑iuijuij′∗=(u†u)j′j=δjj′\sum_iu_{ij}u_{ij'}^*=(u^\dagger u)_{j'j}=\delta_{jj'} を使う。

∑i∣ψ~i⟩⟨ψ~i∣=∑j,j′(∑iuijuij′∗)∣ϕ~j⟩⟨ϕ~j′∣=∑j∣ϕ~j⟩⟨ϕ~j∣.\sum_i| \tilde\psi_i \rangle \langle \tilde\psi_i |=\sum_{j,j'}\Big(\sum_iu_{ij}u_{ij'}^*\Big)| \tilde\phi_j \rangle \langle \tilde\phi_{j'} |=\sum_j| \tilde\phi_j \rangle \langle \tilde\phi_j |.

(a) ⇒\Rightarrow (b):共通の ρ^\hat\rho を ρ^=∑k=1rλk∣k⟩⟨k∣\hat\rho=\sum_{k=1}^r\lambda_k| k \rangle \langle k |(λk>0\lambda_k>0、rr は階数)と対角化し、∣k~⟩=λk∣k⟩| \tilde k \rangle =\sqrt{\lambda_k}| k \rangle と置く。

まず各 ∣ψ~i⟩| \tilde\psi_i \rangle は ρ^\hat\rho の値域(∣1⟩,…,∣r⟩| 1 \rangle ,\dots,| r \rangle の張る空間)にある。実際、値域に直交する ∣χ⟩| \chi \rangle では0=⟨χ∣ρ^∣χ⟩=∑i∣⟨χ∣ψ~i⟩∣20=\langle \chi | \hat\rho | \chi \rangle =\sum_i\lvert \langle \chi | \tilde\psi_i \rangle \rvert ^2 なので、どの ii でも ⟨χ∣ψ~i⟩=0\langle \chi | \tilde\psi_i \rangle =0 である。そこで

∣ψ~i⟩=∑k=1rcik∣k~⟩,cik=⟨k∣ψ~i⟩λk| \tilde\psi_i \rangle =\sum_{k=1}^rc_{ik}| \tilde k \rangle ,\qquad c_{ik}=\frac{\langle k | \tilde\psi_i \rangle }{\sqrt{\lambda_k}}

と展開できる。これを (a) の左辺に入れると ∑k,k′(∑icikcik′∗)∣k~⟩⟨k~′∣=∑k∣k~⟩⟨k~∣\sum_{k,k'}\big(\sum_ic_{ik}c_{ik'}^*\big)| \tilde k \rangle \langle \tilde k' |=\sum_k| \tilde k \rangle \langle \tilde k | で、∣k~⟩⟨k~′∣| \tilde k \rangle \langle \tilde k' | は線形独立だから ∑icik′∗cik=δkk′\sum_ic_{ik'}^*c_{ik}=\delta_{kk'}、つまり m×rm\times r 行列 cc の列は正規直交である(c†c=1rc^\dagger c=\mathbf 1_r。cc は等長)。正規直交な rr 本の列に m−rm-r 本の列を足して m×mm\times m のユニタリ行列 CC にでき、k>rk>r では ∣k~⟩≔0| \tilde k \rangle \coloneqq0 と置けば∣ψ~i⟩=∑k=1mCik∣k~⟩| \tilde\psi_i \rangle =\sum_{k=1}^mC_{ik}| \tilde k \rangle が成り立つ。同じ手順で ∣ϕ~j⟩=∑kDjk∣k~⟩| \tilde\phi_j \rangle =\sum_kD_{jk}| \tilde k \rangle となるユニタリ DD がとれる。DD はユニタリなので ∣k~⟩=∑jDjk∗∣ϕ~j⟩| \tilde k \rangle =\sum_jD^*_{jk}| \tilde\phi_j \rangle と逆に解け、

∣ψ~i⟩=∑kCik∑jDjk∗∣ϕ~j⟩=∑j(CD†)ij∣ϕ~j⟩.| \tilde\psi_i \rangle =\sum_kC_{ik}\sum_jD_{jk}^*| \tilde\phi_j \rangle =\sum_j(CD^\dagger)_{ij}| \tilde\phi_j \rangle .

u=CD†u=CD^\dagger はユニタリ行列の積なのでユニタリである。

証明の後半から、どのアンサンブルも固有アンサンブルから等長行列 cc で作れ、逆にどの等長行列 cc(c†c=1rc^\dagger c=\mathbf 1_r)からもアンサンブル ∣ψ~i⟩=∑kcik∣k~⟩| \tilde\psi_i \rangle =\sum_kc_{ik}| \tilde k \rangle が1つ作れることが分かる。mm 個の状態からなる ρ^\hat\rho のアンサンブル全体は、m×rm\times r の等長行列全体で媒介変数表示されている。

混ぜられる状態と、その重みの上限

定理から、特定の状態 ∣ψ⟩| \psi \rangle が ρ^\hat\rho のアンサンブルに現れうるか、現れるならどれだけの重みまで許されるかが決まる。

命題 4.3.2(混ぜられる重みの上限)

単位ベクトル ∣ψ⟩| \psi \rangle が重み p>0p>0 で ρ^\hat\rho のアンサンブルに現れるための必要十分条件は、∣ψ⟩| \psi \rangle が ρ^\hat\rho の値域にあり、かつ

p≤pmax⁡(ψ)≔1⟨ψ∣ρ^−1∣ψ⟩(4.3.7)p\le p_{\max}(\psi)\coloneqq\frac{1}{\langle \psi | \hat\rho^{-1} | \psi \rangle }\tag{4.3.7}

を満たすことである。ρ^−1\hat\rho^{-1} は値域の上での逆 ∑kλk−1∣k⟩⟨k∣\sum_k\lambda_k^{-1}| k \rangle \langle k | を表す。

証明

必要性:∣ψ~1⟩=p∣ψ⟩| \tilde\psi_1 \rangle =\sqrt p| \psi \rangle とすると、定理の証明の係数は c1k=p⟨k∣ψ⟩/λkc_{1k}=\sqrt p\langle k | \psi \rangle /\sqrt{\lambda_k} である。cc の列が正規直交なら、m×rm\times r 行列の1行の長さの2乗は1を超えない(cc を m×mm\times m のユニタリ CC に延ばせば、その行の長さは1)。したがって∑k∣c1k∣2=p∑k∣⟨k∣ψ⟩∣2/λk=p⟨ψ∣ρ^−1∣ψ⟩≤1\sum_k\lvert c_{1k}\rvert ^2=p\sum_k\lvert \langle k | \psi \rangle \rvert ^2/\lambda_k=p\langle \psi | \hat\rho^{-1} | \psi \rangle \le1。

十分性:p⟨ψ∣ρ^−1∣ψ⟩≤1p\langle \psi | \hat\rho^{-1} | \psi \rangle \le1 なら ρ^−p∣ψ⟩⟨ψ∣≥0\hat\rho-p| \psi \rangle \langle \psi |\ge0 を示せばよい。そうすればこの差を固有アンサンブルに分けて p∣ψ⟩⟨ψ∣p| \psi \rangle \langle \psi | と並べることで、∣ψ⟩| \psi \rangle を重み pp で含むアンサンブルができる。任意の ∣χ⟩| \chi \rangle について、∣ψ⟩| \psi \rangle が値域にあるので ⟨χ∣ψ⟩=⟨χ∣ρ^1/2ρ^−1/2∣ψ⟩\langle \chi | \psi \rangle =\langle \chi | \hat\rho^{1/2}\hat\rho^{-1/2} | \psi \rangle と書け、Cauchy–Schwarz の不等式から

∣⟨χ∣ψ⟩∣2≤⟨χ∣ρ^∣χ⟩ ⟨ψ∣ρ^−1∣ψ⟩.\lvert \langle \chi | \psi \rangle \rvert ^2\le\langle \chi | \hat\rho | \chi \rangle \,\langle \psi | \hat\rho^{-1} | \psi \rangle .

よって ⟨χ∣(ρ^−p∣ψ⟩⟨ψ∣)∣χ⟩≥⟨χ∣ρ^∣χ⟩(1−p⟨ψ∣ρ^−1∣ψ⟩)≥0\langle \chi | (\hat\rho-p| \psi \rangle \langle \psi |) | \chi \rangle \ge\langle \chi | \hat\rho | \chi \rangle \big(1-p\langle \psi | \hat\rho^{-1} | \psi \rangle \big)\ge0。

量子ビットでは逆行列が ρ^−1=2(1−r⋅σ)/(1−∣r∣2)\hat\rho^{-1}=2(\mathbf 1-\bm{r}\cdot\bm{\sigma})/(1-\lvert \bm{r}\rvert ^2)(∣r∣<1\lvert \bm{r}\rvert <1)と書けるので、Bloch ベクトル n\bm{n} の純粋状態について

pmax⁡(n)=1−∣r∣22(1−r⋅n)(4.3.8)p_{\max}(\bm{n})=\frac{1-\lvert \bm{r}\rvert ^2}{2(1-\bm{r}\cdot\bm{n})}\tag{4.3.8}

となる(⟨n∣r⋅σ∣n⟩=r⋅n\langle \bm{n} | \bm{r}\cdot\bm{\sigma} | \bm{n} \rangle =\bm{r}\cdot\bm{n} を使った)。n\bm{n} が r\bm{r} の向きに近いほど多く混ぜられ、n=r/∣r∣\bm{n}=\bm{r}/\lvert \bm{r}\rvert で固有値 (1+∣r∣)/2(1+\lvert \bm{r}\rvert )/2 に達する。これは弦の幾何でも読める。n\bm{n} を一方の端にとって r\bm{r} を通る弦を引くと、n\bm{n} の重みは r\bm{r} から反対側の端までの長さを弦の長さで割ったもの(てこの規則)で、ほかの状態をさらに n\bm{n} の近くに置いて足し込むことはできないので、これが上限になる。

重みのばらつきは固有値より大きい

固有アンサンブルは、すべてのアンサンブルの中で重みの偏りが最も大きい。固有値を大きい順に λ1≥λ2≥⋯\lambda_1\ge\lambda_2\ge\cdots(mm 個まで 00 を足す)、アンサンブルの重みを大きい順に p[1]≥p[2]≥⋯p_{[1]}\ge p_{[2]}\ge\cdots と並べると、

∑i=1kp[i]≤∑i=1kλi(k=1,…,m)(4.3.9)\sum_{i=1}^kp_{[i]}\le\sum_{i=1}^k\lambda_i\qquad(k=1,\dots,m)\tag{4.3.9}

が成り立つ(pp は λ\lambda に多数化されるという)。理由は定理の証明から読める。pi=⟨ψ~i∣ψ~i⟩=∑k∣Cik∣2λkp_i=\langle \tilde\psi_i | \tilde\psi_i \rangle =\sum_k\lvert C_{ik}\rvert ^2\lambda_k で、行列 Bik=∣Cik∣2B_{ik}=\lvert C_{ik}\rvert ^2 は CC がユニタリなので行の和も列の和も1(二重確率行列)である。二重確率行列は置換行列の凸結合である(Birkhoff の定理)ので、p=Bλp=B\lambda は λ\lambda の成分を並べ替えたものの平均で、上位 kk 個の和は λ\lambda の上位 kk 個の和を超えない。Shannon エントロピー H(p)=−∑ipiln⁡piH(p)=-\sum_ip_i\ln p_i は凹関数で並べ替えで変わらないので、この平均の形から

H(p)≥H(λ)=S(ρ^)=−Tr⁡(ρ^ln⁡ρ^)(4.3.10)H(p)\ge H(\lambda)=S(\hat\rho)=-\operatorname{Tr}(\hat\rho\ln\hat\rho)\tag{4.3.10}

を得る。von Neumann エントロピーは、ρ^\hat\rho を純粋状態の混合として書くときの「どれを混ぜたかの不確かさ」の最小値で、最小は固有アンサンブルで達成される。第4部第4章(第4部第4章)で熱平衡のエントロピーを S(ρ^)S(\hat\rho) で書くとき、この最小性がその意味づけの1つになる。

図 4.3.2 は、図 4.3.1 の ρ^=diag(34,14)\hat\rho=\mathrm{diag}(\tfrac34,\tfrac14) について、弦の向き φ\varphi を 00 から π\pi まで回したときの2点アンサンブルの H(p)H(p) である。φ=0,π\varphi=0,\pi(固有アンサンブル)で最小値 S(ρ^)=hb(14)=0.5623S(\hat\rho)=h_{\mathrm b}(\tfrac14)=0.5623 をとり、φ=π/2\varphi=\pi/2(横の弦、重み 12\tfrac12 ずつ)で ln⁡2=0.6931\ln2=0.6931 に達する。

図

図 4.3.2 ρ^=diag(34,14)\hat\rho=\mathrm{diag}(\tfrac34,\tfrac14) を Bloch 球の xzxz 面内の弦(zz 軸から角 φ\varphi)で2点に分けたときの重みの Shannon エントロピー H(p)H(p)(実線、(4.3.4)・(4.3.5) から計算)。点線は S(ρ^)S(\hat\rho)。どの向きでも H(p)≥S(ρ^)H(p)\ge S(\hat\rho) で、等号は固有アンサンブル(φ=0,180∘\varphi=0,180^\circ)だけで成り立つ。

例題 4.3.2 移り合うユニタリと、混ぜられる重みの上限 ★★

オリジナルρ^=diag(34,14)\hat\rho=\mathrm{diag}(\tfrac34,\tfrac14)(基底 ∣↑⟩,∣↓⟩| \uparrow \rangle ,| \downarrow \rangle)について答えよ。

  1. (1) 固有アンサンブルの規格化しないベクトル ∣ψ~1⟩=32∣↑⟩| \tilde\psi_1 \rangle =\tfrac{\sqrt3}{2}| \uparrow \rangle、∣ψ~2⟩=12∣↓⟩| \tilde\psi_2 \rangle =\tfrac12| \downarrow \rangle と、(4.3.6) の ∣ϕ~±⟩=12∣a±⟩| \tilde\phi_\pm \rangle =\tfrac{1}{\sqrt2}| a_\pm \rangle を結ぶユニタリ行列を求めよ。
  2. (2) ∣↑⟩| \uparrow \rangle、∣+⟩| + \rangle、∣a+⟩| a_+ \rangle をこの ρ^\hat\rho のアンサンブルに混ぜられる重みの上限を、(4.3.7) と (4.3.8) の両方で求めよ。
  3. (3) (1) の2つのアンサンブルの重みについて (4.3.10) を確かめよ。
解答

(1) ∣ϕ~±⟩=12(32∣↑⟩±12∣↓⟩)=12(∣ψ~1⟩±∣ψ~2⟩)| \tilde\phi_\pm \rangle =\tfrac1{\sqrt2}\big(\tfrac{\sqrt3}{2}| \uparrow \rangle \pm\tfrac12| \downarrow \rangle \big)=\tfrac1{\sqrt2}\big(| \tilde\psi_1 \rangle \pm| \tilde\psi_2 \rangle \big) なので、

(∣ϕ~+⟩∣ϕ~−⟩)=12(111−1)(∣ψ~1⟩∣ψ~2⟩).\begin{pmatrix}| \tilde\phi_+ \rangle \\ | \tilde\phi_- \rangle \end{pmatrix}=\frac1{\sqrt2}\begin{pmatrix}1&1\\1&-1\end{pmatrix}\begin{pmatrix}| \tilde\psi_1 \rangle \\ | \tilde\psi_2 \rangle \end{pmatrix}.

係数行列は Hadamard 行列で、ユニタリである。

(2) ρ^−1=diag(43,4)\hat\rho^{-1}=\mathrm{diag}(\tfrac43,4)。∣↑⟩| \uparrow \rangle:⟨↑∣ρ^−1∣↑⟩=43\langle \uparrow | \hat\rho^{-1} | \uparrow \rangle =\tfrac43 で pmax⁡=34p_{\max}=\tfrac34。∣+⟩| + \rangle:12(43+4)=83\tfrac12(\tfrac43+4)=\tfrac83 で pmax⁡=38p_{\max}=\tfrac38。∣a+⟩| a_+ \rangle:34⋅43+14⋅4=2\tfrac34\cdot\tfrac43+\tfrac14\cdot4=2 で pmax⁡=12p_{\max}=\tfrac12。(4.3.8) では ∣r∣2=14\lvert \bm{r}\rvert ^2=\tfrac14、1−∣r∣2=341-\lvert \bm{r}\rvert ^2=\tfrac34 で、n=(0,0,1)\bm{n}=(0,0,1) なら r⋅n=12\bm{r}\cdot\bm{n}=\tfrac12、pmax⁡=3/42⋅1/2=34p_{\max}=\tfrac{3/4}{2\cdot1/2}=\tfrac34。n=(1,0,0)\bm{n}=(1,0,0) なら r⋅n=0\bm{r}\cdot\bm{n}=0、pmax⁡=38p_{\max}=\tfrac38。n=(32,0,12)\bm{n}=(\tfrac{\sqrt3}2,0,\tfrac12) なら r⋅n=14\bm{r}\cdot\bm{n}=\tfrac14、pmax⁡=3/43/2=12p_{\max}=\tfrac{3/4}{3/2}=\tfrac12。両方の式が一致する。∣↑⟩| \uparrow \rangle と ∣a+⟩| a_+ \rangle の上限が (1) のアンサンブルの重みにちょうど等しいのは、どちらも r\bm{r} を通る弦の端で、残りを反対側の1点に置き切っているからである。∣+⟩| + \rangle は、重み 12\tfrac12 の ∣↑⟩| \uparrow \rangle を足しても ρ^\hat\rho には届かず、38\tfrac38 が限界である。

(3) 固有アンサンブルは H=hb(14)=−14ln⁡14−34ln⁡34=0.5623H=h_{\mathrm b}(\tfrac14)=-\tfrac14\ln\tfrac14-\tfrac34\ln\tfrac34=0.5623 で、これが S(ρ^)S(\hat\rho)。横の弦は H=ln⁡2=0.6931H=\ln2=0.6931。0.6931≥0.56230.6931\ge0.5623 で (4.3.10) が成り立つ。重みで見ると (12,12)(\tfrac12,\tfrac12) は (34,14)(\tfrac34,\tfrac14) に多数化される(12≤34\tfrac12\le\tfrac34)。

木下の勉強ステータス

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

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

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