tecktom.com

Schmidt 分解

2026年10月17日11分で読めますAI生成教科書『虎の巻』第4部 第2章 基準との同値:部分トレース・Schmidt 分解・浄化
目次
  1. 定理と証明
  2. 行列で見ると特異値分解
  3. Schmidt 分解から分かること
  4. 木下の層

AI 生成・木下が未勉強

量子力学の様々な構成方法 第4部 第2章 §4.2.2 ★★★ 2つの系の純粋状態は、両側の正規直交系の対角な和 ∑kλk∣ak⟩∣bk⟩\sum_k\sqrt{\lambda_k}| a_k \rangle | b_k \rangle に必ず書ける。両側の縮約密度演算子は同じ非零固有値 λk\lambda_k を持つ

定理と証明

定理 4.2.1(Schmidt 分解 ★★★)

∣Ψ⟩∈HA⊗HB| \Psi \rangle \in\mathcal{H}_{\mathrm A}\otimes\mathcal{H}_{\mathrm B} を単位ベクトルとし(HA\mathcal{H}_{\mathrm A} は有限次元とする)、ρ^A=Tr⁡B∣Ψ⟩⟨Ψ∣\hat\rho_{\mathrm A}=\operatorname{Tr}_{\mathrm B}| \Psi \rangle \langle \Psi | の正の固有値を λ1,…,λr\lambda_1,\dots,\lambda_r、対応する正規直交な固有ベクトルを ∣a1⟩,…,∣ar⟩| a_1 \rangle ,\dots,| a_r \rangle とする。このとき

∣bk⟩≔1λk(⟨ak∣⊗1B)∣Ψ⟩(k=1,…,r)(4.2.3)| b_k \rangle \coloneqq\frac{1}{\sqrt{\lambda_k}}\big(\langle a_k |\otimes\mathbf 1_{\mathrm B}\big)| \Psi \rangle \quad(k=1,\dots,r)\tag{4.2.3}

は HB\mathcal{H}_{\mathrm B} の正規直交系で、

∣Ψ⟩=∑k=1rλk ∣ak⟩∣bk⟩(4.2.4)| \Psi \rangle =\sum_{k=1}^r\sqrt{\lambda_k}\,| a_k \rangle | b_k \rangle\tag{4.2.4}

が成り立つ。さらに ρ^B=Tr⁡A∣Ψ⟩⟨Ψ∣=∑kλk∣bk⟩⟨bk∣\hat\rho_{\mathrm B}=\operatorname{Tr}_{\mathrm A}| \Psi \rangle \langle \Psi |=\sum_k\lambda_k| b_k \rangle \langle b_k | で、ρ^A\hat\rho_{\mathrm A} と ρ^B\hat\rho_{\mathrm B} の非零固有値は重複も込めて一致する。

証明

HA\mathcal{H}_{\mathrm A} の正規直交基底を ρ^A\hat\rho_{\mathrm A} の固有ベクトル全体 {∣ak⟩}\{| a_k \rangle \}(固有値 00 のものも含める)に取る。(4.2.3) の右辺の分子を ∣b~k⟩=(⟨ak∣⊗1)∣Ψ⟩| \tilde b_k \rangle =(\langle a_k |\otimes\mathbf 1)| \Psi \rangle と書くと、内積は

⟨b~k∣b~l⟩=⟨Ψ∣(∣ak⟩⟨al∣⊗1)∣Ψ⟩=Tr⁡[∣Ψ⟩⟨Ψ∣(∣ak⟩⟨al∣⊗1)]=Tr⁡A(ρ^A∣ak⟩⟨al∣)=⟨al∣ρ^A∣ak⟩=λkδkl\langle \tilde b_k | \tilde b_l \rangle =\langle \Psi |\big(| a_k \rangle \langle a_l |\otimes\mathbf 1\big)| \Psi \rangle =\operatorname{Tr}\big[| \Psi \rangle \langle \Psi |\big(| a_k \rangle \langle a_l |\otimes\mathbf 1\big)\big]=\operatorname{Tr}_{\mathrm A}\big(\hat\rho_{\mathrm A}| a_k \rangle \langle a_l |\big)=\langle a_l | \hat\rho_{\mathrm A} | a_k \rangle =\lambda_k\delta_{kl}

である。3つ目の等号で部分トレースの特徴づけ(第4部第1章 の命題)を使った。したがって λk>0\lambda_k>0 の kk では ∣bk⟩=∣b~k⟩/λk| b_k \rangle =| \tilde b_k \rangle /\sqrt{\lambda_k} が正規直交で、λk=0\lambda_k=0 の kk では ∥b~k∥2=0\lVert \tilde b_k\rVert ^2=0 から ∣b~k⟩=0| \tilde b_k \rangle =0 である。A の基底の完全性 ∑k∣ak⟩⟨ak∣=1A\sum_k| a_k \rangle \langle a_k |=\mathbf 1_{\mathrm A} を ∣Ψ⟩| \Psi \rangle に当てると

∣Ψ⟩=∑k(∣ak⟩⟨ak∣⊗1)∣Ψ⟩=∑k∣ak⟩⊗∣b~k⟩=∑k:λk>0λk∣ak⟩∣bk⟩| \Psi \rangle =\sum_k\big(| a_k \rangle \langle a_k |\otimes\mathbf 1\big)| \Psi \rangle =\sum_k| a_k \rangle \otimes| \tilde b_k \rangle =\sum_{k:\lambda_k>0}\sqrt{\lambda_k}| a_k \rangle | b_k \rangle

で (4.2.4) を得る。最後に、(4.2.4) から ∣Ψ⟩⟨Ψ∣=∑k,lλkλl∣ak⟩⟨al∣⊗∣bk⟩⟨bl∣| \Psi \rangle \langle \Psi |=\sum_{k,l}\sqrt{\lambda_k\lambda_l}| a_k \rangle \langle a_l |\otimes| b_k \rangle \langle b_l | で、A についてトレースを取ると⟨al∣ak⟩=δkl\langle a_l | a_k \rangle =\delta_{kl} から ρ^B=∑kλk∣bk⟩⟨bk∣\hat\rho_{\mathrm B}=\sum_k\lambda_k| b_k \rangle \langle b_k |。∣bk⟩| b_k \rangle は正規直交なので、これは ρ^B\hat\rho_{\mathrm B} の固有値分解で、非零固有値は λk\lambda_k である。

正の係数 λk\sqrt{\lambda_k} を Schmidt 係数、その個数 rr を Schmidt 階数と呼ぶ。名前は E. Schmidt が積分方程式の核 K(s,t)K(s,t) を2組の直交関数系の対角な和に展開した仕事(Schmidt (1907), Math. Ann. 63, 433)に由来する。2つの系の状態の言葉で書き直したのは Schrödinger で、同じ論文で「もつれ」の概念を導入した(Schrödinger (1935), Proc. Cambridge Phil. Soc. 31, 555)。

行列で見ると特異値分解

基底 ∣i⟩A| i \rangle _{\mathrm A}、∣j⟩B| j \rangle _{\mathrm B} で ∣Ψ⟩=∑i,jCij∣i⟩∣j⟩| \Psi \rangle =\sum_{i,j}C_{ij}| i \rangle | j \rangle と書くと、係数行列 CC は dA×dBd_{\mathrm A}\times d_{\mathrm B} の長方行列で、ρ^A\hat\rho_{\mathrm A} の行列は CC†CC^\dagger、ρ^B\hat\rho_{\mathrm B} の行列は CTC‾C^{\mathsf T}\overline{C} になる(⟨i∣ρ^A∣i′⟩=∑jCijCi′j‾\langle i | \hat\rho_{\mathrm A} | i' \rangle =\sum_jC_{ij}\overline{C_{i'j}})。定理 4.2.1 は、CC を C=∑kλk ukwkTC=\sum_k\sqrt{\lambda_k}\,\bm{u}_k\bm{w}_k^{\mathsf T}(uk\bm{u}_k、wk\bm{w}_k はそれぞれ正規直交)と書く特異値分解と同じ内容である。Schmidt 係数は CC の特異値で、その2乗が ρ^A\hat\rho_{\mathrm A} の固有値である。CC†CC^\dagger と CTC‾=C†C‾C^{\mathsf T}\overline C=\overline{C^\dagger C} が同じ非零固有値を持つことは、行列の言葉では CC†CC^\dagger と C†CC^\dagger C の非零固有値が一致するという線型代数の事実に当たる。

2量子ビットでは CC は 2×22\times2 で、∣det⁡C∣=λ1λ2\lvert \det C\rvert =\sqrt{\lambda_1\lambda_2} である(行列式は特異値の積に位相を掛けたもの)。第1部第2章で示した「積状態 ⇔det⁡C=0\Leftrightarrow\det C=0」は、Schmidt 階数が1のときに限って積状態になる、という一般の判定の2量子ビット版である。

Schmidt 分解から分かること

定理 4.2.1 から、もつれについて3つのことがすぐに分かる。第一に、全体が純粋なら、A と B の混合の度合いは同じである。純粋度 Tr⁡ρ^A2=Tr⁡ρ^B2=∑kλk2\operatorname{Tr}\hat\rho_{\mathrm A}^2=\operatorname{Tr}\hat\rho_{\mathrm B}^2=\sum_k\lambda_k^2 も、von Neumann エントロピー S(ρ^A)=S(ρ^B)=−∑kλkln⁡λkS(\hat\rho_{\mathrm A})=S(\hat\rho_{\mathrm B})=-\sum_k\lambda_k\ln\lambda_k も一致する。この共通の値 −∑kλkln⁡λk-\sum_k\lambda_k\ln\lambda_k を純粋状態のもつれのエントロピーと呼ぶ(もつれの定量化の一般論は 相対エントロピーから情報と物理を組み立てる 第3部第4章)。次元の違う系でも同じなので、たとえば1量子ビットと巨大な環境が全体として純粋なら、環境の縮約密度演算子の非零固有値は高々2個しかない。第二に、∣Ψ⟩| \Psi \rangle が積状態であることと、Schmidt 階数が1であることと、ρ^A\hat\rho_{\mathrm A} が純粋であることは同値である。第三に、Schmidt 階数は min⁡(dA,dB)\min(d_{\mathrm A},d_{\mathrm B}) を超えない。

3つの系には Schmidt 分解は無い

3つの系 A, B, C の純粋状態を ∑kck∣ak⟩∣bk⟩∣ck⟩\sum_k c_k| a_k \rangle | b_k \rangle | c_k \rangle(各系で正規直交)の形に書けるとは限らない。定理 4.2.1 の証明は、A とそれ以外の2つに分けて使うことはできるが、3つを同時に対角にすることはできない。たとえば (∣↑↑↓⟩+∣↑↓↑⟩+∣↓↑↑⟩)/3(| \uparrow\uparrow\downarrow \rangle +| \uparrow\downarrow\uparrow \rangle +| \downarrow\uparrow\uparrow \rangle )/\sqrt3 はこの形に書けない。浄化の議論で使うのは常に「注目する系」と「それ以外全部」の2つへの分け方なので、この章では問題にならない。

例題 4.2.2 2量子ビットの状態の Schmidt 分解 ★★

オリジナル∣Ψ⟩=110(2∣↑↑⟩+∣↑↓⟩+∣↓↑⟩+2∣↓↓⟩)| \Psi \rangle =\tfrac{1}{\sqrt{10}}\big(2| \uparrow\uparrow \rangle +| \uparrow\downarrow \rangle +| \downarrow\uparrow \rangle +2| \downarrow\downarrow \rangle \big)(左が A、右が B)とする。

  1. (1) ρ^A\hat\rho_{\mathrm A} を求め、その固有値と固有ベクトルを求めよ。
  2. (2) 定理 4.2.1 の手順で ∣bk⟩| b_k \rangle を作り、Schmidt 分解を書け。
  3. (3) ∣det⁡C∣\lvert \det C\rvert と λ1λ2\sqrt{\lambda_1\lambda_2} を比べ、もつれのエントロピーを求めよ。
解答

(1) 係数行列は C=110(2112)C=\tfrac{1}{\sqrt{10}}\begin{pmatrix}2&1\\1&2\end{pmatrix} で、ρ^A=CC†=110(5445)=12(1+45σx)\hat\rho_{\mathrm A}=CC^\dagger=\tfrac1{10}\begin{pmatrix}5&4\\4&5\end{pmatrix}=\tfrac12\big(\mathbf 1+\tfrac45\sigma_x\big)。固有ベクトルは σx\sigma_x の固有状態 ∣±⟩| \pm \rangle で、固有値は λ+=12(1+45)=910\lambda_+=\tfrac12(1+\tfrac45)=\tfrac{9}{10}、λ−=110\lambda_-=\tfrac1{10}。

(2) ⟨±∣⊗1\langle \pm |\otimes\mathbf 1 を当てる。⟨+∣↑⟩=⟨+∣↓⟩=12\langle + | \uparrow \rangle =\langle + | \downarrow \rangle =\tfrac1{\sqrt2} なので

(⟨+∣⊗1)∣Ψ⟩=110⋅12(2∣↑⟩+∣↓⟩+∣↑⟩+2∣↓⟩)=310∣+⟩,(\langle + |\otimes\mathbf 1)| \Psi \rangle =\frac{1}{\sqrt{10}}\cdot\frac{1}{\sqrt2}\big(2| \uparrow \rangle +| \downarrow \rangle +| \uparrow \rangle +2| \downarrow \rangle \big)=\frac{3}{\sqrt{10}}| + \rangle ,

λ+=310\sqrt{\lambda_+}=\tfrac{3}{\sqrt{10}} で割って ∣b+⟩=∣+⟩| b_+ \rangle =| + \rangle。同様に (⟨−∣⊗1)∣Ψ⟩=120(2∣↑⟩+∣↓⟩−∣↑⟩−2∣↓⟩)=110∣−⟩(\langle - |\otimes\mathbf 1)| \Psi \rangle =\tfrac{1}{\sqrt{20}}(2| \uparrow \rangle +| \downarrow \rangle -| \uparrow \rangle -2| \downarrow \rangle )=\tfrac{1}{\sqrt{10}}| - \rangle で ∣b−⟩=∣−⟩| b_- \rangle =| - \rangle。したがって

∣Ψ⟩=310∣+⟩∣+⟩+110∣−⟩∣−⟩.| \Psi \rangle =\frac{3}{\sqrt{10}}| + \rangle | + \rangle +\frac{1}{\sqrt{10}}| - \rangle | - \rangle .

展開すると 1210(4∣↑↑⟩+2∣↑↓⟩+2∣↓↑⟩+4∣↓↓⟩)\tfrac{1}{2\sqrt{10}}\big(4| \uparrow\uparrow \rangle +2| \uparrow\downarrow \rangle +2| \downarrow\uparrow \rangle +4| \downarrow\downarrow \rangle \big) で元に戻る。

(3) det⁡C=110(4−1)=310\det C=\tfrac{1}{10}(4-1)=\tfrac{3}{10}、λ+λ−=9/100=310\sqrt{\lambda_+\lambda_-}=\sqrt{9/100}=\tfrac3{10} で一致する。もつれのエントロピーは−910ln⁡910−110ln⁡110=0.0948+0.2303=0.325 nat-\tfrac9{10}\ln\tfrac9{10}-\tfrac1{10}\ln\tfrac1{10}=0.0948+0.2303=0.325\ \mathrm{nat}。最大の ln⁡2=0.693\ln2=0.693 の半分に満たない。

木下の勉強ステータス

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

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

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