これまでに、量子チャネルを数学的に表現する3つの異なる方法、すなわちStinespring表現、Kraus表現、Choi表現について説明しました。
また、チャネルの定義として、チャネルとは密度行列を常に密度行列へと変換する線形写像であり、複合システムの一部にのみチャネルを適用する場合でもこの性質が成り立つことを確認しました。
このレッスンの残りは、この3つの表現が等価であり、定義を正確に捉えていることの数学的証明に充てられます。
証明の概要
私たちの目標は、4つの命題の等価性を示すことです。まずそれらを正確に書き下しましょう。
4つの命題はすべて、レッスン全体で用いてきた同じ規約に従います。すなわち、Φ \Phi Φ は正方行列から正方行列への線形写像であり、入力行列の行と列はシステム X \mathsf{X} X (入力システム)の古典状態に対応し、出力行列の行と列はシステム Y \mathsf{Y} Y (出力システム)の古典状態に対応しています。
Φ \Phi Φ は X \mathsf{X} X から Y \mathsf{Y} Y へのチャネルである。すなわち、Φ \Phi Φ はより大きな複合システムの一部に作用する場合でも、常に密度行列を密度行列に変換する。
Choi行列 J ( Φ ) J(\Phi) J ( Φ ) は半正定値であり、条件 Tr Y ( J ( Φ ) ) = I X \operatorname{Tr}_{\mathsf{Y}}(J(\Phi)) = \mathbb{I}_{\mathsf{X}} Tr Y ( J ( Φ )) = I X を満たす。
Φ \Phi Φ のKraus表現が存在する。すなわち、すべての入力 ρ \rho ρ に対して等式 Φ ( ρ ) = ∑ k = 0 N − 1 A k ρ A k † \Phi(\rho) = \sum_{k = 0}^{N-1} A_k \rho A_k^{\dagger} Φ ( ρ ) = ∑ k = 0 N − 1 A k ρ A k † が成り立ち、かつ条件 ∑ k = 0 N − 1 A k † A k = I X \sum_{k = 0}^{N-1} A_k^{\dagger} A_k = \mathbb{I}_{\mathsf{X}} ∑ k = 0 N − 1 A k † A k = I X を満たす行列 A 0 , … , A N − 1 A_0,\ldots,A_{N-1} A 0 , … , A N − 1 が存在する。
Φ \Phi Φ のStinespring表現が存在する。すなわち、システム W \mathsf{W} W と G \mathsf{G} G が存在し、( W , X ) (\mathsf{W},\mathsf{X}) ( W , X ) と ( G , Y ) (\mathsf{G},\mathsf{Y}) ( G , Y ) の古典状態の数が等しく、さらに ( W , X ) (\mathsf{W},\mathsf{X}) ( W , X ) から ( G , Y ) (\mathsf{G},\mathsf{Y}) ( G , Y ) へのユニタリ演算を表すユニタリ行列 U U U が存在して、Φ ( ρ ) = Tr G ( U ( ∣ 0 ⟩ ⟨ 0 ∣ ⊗ ρ ) U † ) \Phi(\rho) = \operatorname{Tr}_{\mathsf{G}}\bigl( U (\vert 0\rangle\langle 0 \vert \otimes \rho) U^{\dagger} \bigr) Φ ( ρ ) = Tr G ( U ( ∣0 ⟩ ⟨ 0∣ ⊗ ρ ) U † ) が成り立つ。
証明の方針は、含意の連鎖を示すことです。
リスト中の1番目の命題が2番目を含意し、2番目が3番目を含意し、3番目が4番目を含意し、4番目が1番目を含意することを証明します。
これにより、4つの命題すべてが等価である——すなわち、与えられた Φ \Phi Φ に対してすべてが真かすべてが偽である——ことが確立されます。なぜなら、含意を推移的にたどることで任意の1つの命題から別の任意の命題へと到達できるからです。
これは複数の命題の等価性を証明する際の一般的な戦略であり、そのような文脈で有用な工夫は、含意をできるだけ証明しやすい方向に設定することです。
実際ここでもそうなっており、4つの含意のうち2つはすでに扱っています。
チャネルからChoi行列へ
上記の番号で参照すると、最初に証明すべき含意は 1 ⇒ \Rightarrow ⇒ 2 です。
この含意はチャネルのChoi状態の文脈ですでに議論されています。
ここでは数学的な詳細をまとめます。
入力システム X \mathsf{X} X の古典状態集合を Σ \Sigma Σ とし、n = ∣ Σ ∣ n = \vert\Sigma\vert n = ∣Σ∣ とします。
Φ \Phi Φ が X \mathsf{X} X の2つのコピーからなるシステムの第2のコピーに適用される状況を考えます。2つのコピーが合わせて状態
∣ ψ ⟩ = 1 n ∑ a ∈ Σ ∣ a ⟩ ⊗ ∣ a ⟩ , \vert \psi \rangle = \frac{1}{\sqrt{n}} \sum_{a \in \Sigma} \vert a \rangle \otimes \vert a \rangle, ∣ ψ ⟩ = n 1 a ∈ Σ ∑ ∣ a ⟩ ⊗ ∣ a ⟩ ,
にある場合、密度行列として書くと
∣ ψ ⟩ ⟨ ψ ∣ = 1 n ∑ a , b ∈ Σ ∣ a ⟩ ⟨ b ∣ ⊗ ∣ a ⟩ ⟨ b ∣ \vert \psi \rangle \langle \psi \vert = \frac{1}{n} \sum_{a,b \in \Sigma}
\vert a\rangle\langle b \vert \otimes \vert a\rangle\langle b \vert ∣ ψ ⟩ ⟨ ψ ∣ = n 1 a , b ∈ Σ ∑ ∣ a ⟩ ⟨ b ∣ ⊗ ∣ a ⟩ ⟨ b ∣
となります。
結果は次のように書けます。
( Id ⊗ Φ ) ( ∣ ψ ⟩ ⟨ ψ ∣ ) = 1 n ∑ a , b = 0 n − 1 ∣ a ⟩ ⟨ b ∣ ⊗ Φ ( ∣ a ⟩ ⟨ b ∣ ) = J ( Φ ) n , (\operatorname{Id}\otimes \,\Phi) \bigl(\vert \psi \rangle \langle \psi \vert\bigr)
= \frac{1}{n} \sum_{a,b = 0}^{n-1} \vert a\rangle\langle b \vert \otimes \Phi\bigl(\vert a\rangle\langle b \vert\bigr)
= \frac{J(\Phi)}{n}, ( Id ⊗ Φ ) ( ∣ ψ ⟩ ⟨ ψ ∣ ) = n 1 a , b = 0 ∑ n − 1 ∣ a ⟩ ⟨ b ∣ ⊗ Φ ( ∣ a ⟩ ⟨ b ∣ ) = n J ( Φ ) ,
Φ \Phi Φ がチャネルであるという仮定から、これは密度行列でなければなりません。
すべての密度行列と同様に半正定値でなければならず、半正定値行列に正の実数を掛けても半正定値行列となります。したがって J ( Φ ) ≥ 0 J(\Phi) \geq 0 J ( Φ ) ≥ 0 です。
さらに、Φ \Phi Φ がチャネルであるという仮定のもとでトレースを保存しなければならないため、
Tr Y ( J ( Φ ) ) = ∑ a , b ∈ Σ Tr ( Φ ( ∣ a ⟩ ⟨ b ∣ ) ) ∣ a ⟩ ⟨ b ∣ = ∑ a , b ∈ Σ Tr ( ∣ a ⟩ ⟨ b ∣ ) ∣ a ⟩ ⟨ b ∣ = ∑ a ∈ Σ ∣ a ⟩ ⟨ a ∣ = I X \begin{aligned}
\operatorname{Tr}_{\mathsf{Y}} (J(\Phi))
& = \sum_{a,b\in\Sigma} \operatorname{Tr}\bigl(\Phi( \vert a\rangle\langle b \vert)\bigr) \, \vert a\rangle\langle b \vert\\
& = \sum_{a,b\in\Sigma} \operatorname{Tr}\bigl(\vert a\rangle\langle b \vert\bigr) \, \vert a\rangle\langle b \vert\\
& = \sum_{a\in\Sigma} \vert a\rangle\langle a \vert\\
& = \mathbb{I}_{\mathsf{X}}
\end{aligned} Tr Y ( J ( Φ )) = a , b ∈ Σ ∑ Tr ( Φ ( ∣ a ⟩ ⟨ b ∣ ) ) ∣ a ⟩ ⟨ b ∣ = a , b ∈ Σ ∑ Tr ( ∣ a ⟩ ⟨ b ∣ ) ∣ a ⟩ ⟨ b ∣ = a ∈ Σ ∑ ∣ a ⟩ ⟨ a ∣ = I X
が成り立ちます。
ChoiからKraus表現へ
2番目の含意は、リストの番号で参照すると 2 ⇒ \Rightarrow ⇒ 3 です。
明確にしておくと、他の命題は無視します。特に Φ \Phi Φ がチャネルであるという仮定はできません。
使えるのは、Φ \Phi Φ がChoi表現について J ( Φ ) ≥ 0 J(\Phi) \geq 0 J ( Φ ) ≥ 0 かつ
Tr Y ( J ( Φ ) ) = I X \operatorname{Tr}_{\mathsf{Y}} (J(\Phi)) = \mathbb{I}_{\mathsf{X}} Tr Y ( J ( Φ )) = I X を満たす線形写像であるという事実のみです。
しかし、これだけで Φ \Phi Φ のKraus表現
Φ ( ρ ) = ∑ k = 0 N − 1 A k ρ A k † \Phi(\rho) = \sum_{k = 0}^{N-1} A_k \rho A_k^{\dagger} Φ ( ρ ) = k = 0 ∑ N − 1 A k ρ A k †
であって、条件
∑ k = 0 N − 1 A k † A k = I X \sum_{k = 0}^{N-1} A_k^{\dagger} A_k = \mathbb{I}_{\mathsf{X}} k = 0 ∑ N − 1 A k † A k = I X
を満たすものが存在することを結論づけるのに十分です。
J ( Φ ) J(\Phi) J ( Φ ) が半正定値であるという決定的に重要な仮定から始めます。これは、
J ( Φ ) = ∑ k = 0 N − 1 ∣ ψ k ⟩ ⟨ ψ k ∣ (1) J(\Phi) = \sum_{k = 0}^{N-1} \vert \psi_k \rangle \langle \psi_k \vert
\tag{1} J ( Φ ) = k = 0 ∑ N − 1 ∣ ψ k ⟩ ⟨ ψ k ∣ ( 1 )
という形で表現できることを意味します。ベクトル ∣ ψ 0 ⟩ , … , ∣ ψ N − 1 ⟩ \vert\psi_0\rangle,\ldots,\vert\psi_{N-1}\rangle ∣ ψ 0 ⟩ , … , ∣ ψ N − 1 ⟩ の選び方は一般に複数あります。これは Φ \Phi Φ のKraus表現を選ぶ自由度と直接対応しています。
このような表現を得る1つの方法は、まずスペクトル定理を用いて
J ( Φ ) = ∑ k = 0 N − 1 λ k ∣ γ k ⟩ ⟨ γ k ∣ , J(\Phi) = \sum_{k = 0}^{N-1} \lambda_k \vert \gamma_k \rangle \langle \gamma_k \vert, J ( Φ ) = k = 0 ∑ N − 1 λ k ∣ γ k ⟩ ⟨ γ k ∣ ,
と書くことです。ここで λ 0 , … , λ N − 1 \lambda_0,\ldots,\lambda_{N-1} λ 0 , … , λ N − 1 は J ( Φ ) J(\Phi) J ( Φ ) の固有値(J ( Φ ) J(\Phi) J ( Φ ) が半正定値なので必ず非負の実数)であり、∣ γ 0 ⟩ , … , ∣ γ N − 1 ⟩ \vert\gamma_0\rangle,\ldots,\vert\gamma_{N-1}\rangle ∣ γ 0 ⟩ , … , ∣ γ N − 1 ⟩ は対応する単位固有ベクトルです。
固有値の選び方には(順番を除いて)自由度はありませんが、固有ベクトルの選び方には自由度があります。特に多重度が1より大きい固有値がある場合は顕著です。
したがって、この表現は一意ではありません——ただ1つの表現があると仮定しているだけです。
いずれにしても、固有値は非負の実数なので非負の平方根を持ち、各 k = 0 , … , N − 1 k = 0,\ldots,N-1 k = 0 , … , N − 1 に対して
∣ ψ k ⟩ = λ k ∣ γ k ⟩ \vert\psi_k\rangle = \sqrt{\lambda_k} \vert \gamma_k\rangle ∣ ψ k ⟩ = λ k ∣ γ k ⟩
を選ぶことで ( 1 ) (1) ( 1 ) の形の表現が得られます。
ただし、表現 ( 1 ) (1) ( 1 ) がこのようなスペクトル分解から得られる必要はありません。特に、ベクトル ∣ ψ 0 ⟩ , … , ∣ ψ N − 1 ⟩ \vert\psi_0\rangle,\ldots,\vert\psi_{N-1}\rangle ∣ ψ 0 ⟩ , … , ∣ ψ N − 1 ⟩ は一般には直交している必要はありません。
しかし、必要であればこれらのベクトルを直交するように選ぶこともできます。また、N N N を n m nm nm より大きくする必要はありません
(ここで n n n と m m m はそれぞれ X \mathsf{X} X と Y \mathsf{Y} Y の古典状態の数)。
次に、各ベクトル ∣ ψ 0 ⟩ , … , ∣ ψ N − 1 ⟩ \vert\psi_0\rangle,\ldots,\vert\psi_{N-1}\rangle ∣ ψ 0 ⟩ , … , ∣ ψ N − 1 ⟩ をさらに分解します。
∣ ψ k ⟩ = ∑ a ∈ Σ ∣ a ⟩ ⊗ ∣ ϕ k , a ⟩ , \vert\psi_k\rangle = \sum_{a\in\Sigma} \vert a\rangle \otimes \vert \phi_{k,a}\rangle, ∣ ψ k ⟩ = a ∈ Σ ∑ ∣ a ⟩ ⊗ ∣ ϕ k , a ⟩ ,
ここでベクトル { ∣ ϕ k , a ⟩ } \{ \vert \phi_{k,a}\rangle \} { ∣ ϕ k , a ⟩} は Y \mathsf{Y} Y の古典状態に対応する成分を持ち、各 a ∈ Σ a\in\Sigma a ∈ Σ と k = 0 , … , N − 1 k=0,\ldots,N-1 k = 0 , … , N − 1 に対して明示的に
∣ ϕ k , a ⟩ = ( ⟨ a ∣ ⊗ I Y ) ∣ ψ k ⟩ \vert \phi_{k,a}\rangle = \bigl( \langle a \vert \otimes \mathbb{I}_{\mathsf{Y}}\bigr) \vert \psi_k\rangle ∣ ϕ k , a ⟩ = ( ⟨ a ∣ ⊗ I Y ) ∣ ψ k ⟩
によって決まります。
∣ ψ 0 ⟩ , … , ∣ ψ N − 1 ⟩ \vert\psi_0\rangle,\ldots,\vert\psi_{N-1}\rangle ∣ ψ 0 ⟩ , … , ∣ ψ N − 1 ⟩ は必ずしも単位ベクトルではありませんが、これは対 ( X , Y ) (\mathsf{X},\mathsf{Y}) ( X , Y ) の量子状態ベクトルに対してシステム X \mathsf{X} X の標準基底測定を行った場合を解析するときに使うプロセスと同じです。
ここが証明のこの部分を機能させる工夫です。
次の式に従ってKraus行列 A 0 , … , A N − 1 A_0,\ldots,A_{N-1} A 0 , … , A N − 1 を定義します。
A k = ∑ a ∈ Σ ∣ ϕ k , a ⟩ ⟨ a ∣ A_k = \sum_{a\in\Sigma} \vert \phi_{k,a}\rangle\langle a \vert A k = a ∈ Σ ∑ ∣ ϕ k , a ⟩ ⟨ a ∣
この式は純粋に記号的に考えることができます。∣ a ⟩ \vert a\rangle ∣ a ⟩ が反転して ⟨ a ∣ \langle a\vert ⟨ a ∣ となり右辺に移動し、行列を形成します。
証明の検証には、この式だけで十分です。
ただし、ベクトル ∣ ψ k ⟩ \vert\psi_k\rangle ∣ ψ k ⟩ と行列 A k A_k A k の間には単純で直観的な関係があります。それは A k A_k A k をベクトル化 すると ∣ ψ k ⟩ \vert\psi_k\rangle ∣ ψ k ⟩ が得られるということです。
A k A_k A k のベクトル化とは、列を上から(一番左の列を先頭として一番右の列を末尾として)積み重ねてベクトルを形成することです。
例えば、X \mathsf{X} X と Y \mathsf{Y} Y がともにQubitであり、ある k k k の選択に対して
∣ ψ k ⟩ = α 00 ∣ 0 ⟩ ⊗ ∣ 0 ⟩ + α 01 ∣ 0 ⟩ ⊗ ∣ 1 ⟩ + α 10 ∣ 1 ⟩ ⊗ ∣ 0 ⟩ + α 11 ∣ 1 ⟩ ⊗ ∣ 1 ⟩ = ( α 00 α 01 α 10 α 11 ) , \begin{aligned}
\vert\psi_k\rangle & = \alpha_{00} \vert 0\rangle \otimes \vert 0\rangle +
\alpha_{01} \vert 0\rangle \otimes \vert 1\rangle +
\alpha_{10} \vert 1\rangle \otimes \vert 0\rangle +
\alpha_{11} \vert 1\rangle \otimes \vert 1\rangle\\[2mm]
& = \begin{pmatrix}
\alpha_{00} \\[1mm]
\alpha_{01} \\[1mm]
\alpha_{10} \\[1mm]
\alpha_{11}
\end{pmatrix},
\end{aligned} ∣ ψ k ⟩ = α 00 ∣0 ⟩ ⊗ ∣0 ⟩ + α 01 ∣0 ⟩ ⊗ ∣1 ⟩ + α 10 ∣1 ⟩ ⊗ ∣0 ⟩ + α 11 ∣1 ⟩ ⊗ ∣1 ⟩ = α 00 α 01 α 10 α 11 ,
であれば、
A k = α 00 ∣ 0 ⟩ ⟨ 0 ∣ + α 01 ∣ 1 ⟩ ⟨ 0 ∣ + α 10 ∣ 0 ⟩ ⟨ 1 ∣ + α 11 ∣ 1 ⟩ ⟨ 1 ∣ = ( α 00 α 10 α 01 α 11 ) \begin{aligned}
A_k & = \alpha_{00} \vert 0\rangle\langle 0\vert +
\alpha_{01} \vert 1\rangle\langle 0\vert +
\alpha_{10} \vert 0\rangle\langle 1\vert +
\alpha_{11} \vert 1\rangle\langle 1\vert\\[2mm]
& = \begin{pmatrix}
\alpha_{00} & \alpha_{10}\\[1mm]
\alpha_{01} & \alpha_{11}
\end{pmatrix}
\end{aligned} A k = α 00 ∣0 ⟩ ⟨ 0∣ + α 01 ∣1 ⟩ ⟨ 0∣ + α 10 ∣0 ⟩ ⟨ 1∣ + α 11 ∣1 ⟩ ⟨ 1∣ = ( α 00 α 01 α 10 α 11 )
となります。
(注意:行列のベクトル化は、行 を転置して積み重ねて列ベクトルを形成するという、やや異なる定義で使われることもあります。)
まず、このKraus行列の選択が写像 Φ \Phi Φ を正しく記述していることを確認し、その後で必要な別の条件を検証します。
整理のために、新しい写像 Ψ \Psi Ψ を次のように定義します。
Ψ ( ρ ) = ∑ k = 0 N − 1 A k ρ A k † \Psi(\rho) = \sum_{k = 0}^{N-1} A_k \rho A_k^{\dagger} Ψ ( ρ ) = k = 0 ∑ N − 1 A k ρ A k †
したがって、目標は Ψ = Φ \Psi = \Phi Ψ = Φ を検証することです。
これを行うには、両写像のChoi表現を比較します。
Choi表現は忠実なので、J ( Φ ) = J ( Ψ ) J(\Phi) = J(\Psi) J ( Φ ) = J ( Ψ ) であれば Ψ = Φ \Psi = \Phi Ψ = Φ となります。
ここで、式
∣ ψ k ⟩ = ∑ a ∈ Σ ∣ a ⟩ ⊗ ∣ ϕ k , a ⟩ と A k = ∑ a ∈ Σ ∣ ϕ k , a ⟩ ⟨ a ∣ \vert\psi_k\rangle = \sum_{a\in\Sigma} \vert a\rangle \otimes \vert \phi_{k,a}\rangle
\quad\text{と}\quad
A_k = \sum_{a\in\Sigma} \vert \phi_{k,a}\rangle\langle a \vert ∣ ψ k ⟩ = a ∈ Σ ∑ ∣ a ⟩ ⊗ ∣ ϕ k , a ⟩ と A k = a ∈ Σ ∑ ∣ ϕ k , a ⟩ ⟨ a ∣
を用い、テンソル積の双線形性によって整理しながら J ( Ψ ) J(\Psi) J ( Ψ ) を計算できます。
J ( Ψ ) = ∑ a , b ∈ Σ ∣ a ⟩ ⟨ b ∣ ⊗ ∑ k = 0 N − 1 A k ∣ a ⟩ ⟨ b ∣ A k † = ∑ a , b ∈ Σ ∣ a ⟩ ⟨ b ∣ ⊗ ∑ k = 0 N − 1 ∣ ϕ k , a ⟩ ⟨ ϕ k , b ∣ = ∑ k = 0 N − 1 ( ∑ a ∈ Σ ∣ a ⟩ ⊗ ∣ ϕ k , a ⟩ ) ( ∑ b ∈ Σ ⟨ b ∣ ⊗ ⟨ ϕ k , b ∣ ) = ∑ k = 0 N − 1 ∣ ψ k ⟩ ⟨ ψ k ∣ = J ( Φ ) \begin{aligned}
J(\Psi) & = \sum_{a,b\in\Sigma} \vert a\rangle \langle b \vert \otimes \sum_{k = 0}^{N-1} A_k \vert a\rangle \langle b \vert A_k^{\dagger}\\[2mm]
& = \sum_{a,b\in\Sigma} \vert a\rangle \langle b \vert \otimes \sum_{k = 0}^{N-1} \vert \phi_{k,a} \rangle \langle \phi_{k,b} \vert \\[2mm]
& = \sum_{k = 0}^{N-1} \biggl(\sum_{a\in\Sigma} \vert a\rangle \otimes \vert \phi_{k,a} \rangle\biggr)
\biggl(\sum_{b\in\Sigma} \langle b\vert \otimes \langle \phi_{k,b} \vert\biggr)\\[2mm]
& = \sum_{k = 0}^{N-1} \vert \psi_k \rangle \langle \psi_k \vert \\[2mm]
& = J(\Phi)
\end{aligned} J ( Ψ ) = a , b ∈ Σ ∑ ∣ a ⟩ ⟨ b ∣ ⊗ k = 0 ∑ N − 1 A k ∣ a ⟩ ⟨ b ∣ A k † = a , b ∈ Σ ∑ ∣ a ⟩ ⟨ b ∣ ⊗ k = 0 ∑ N − 1 ∣ ϕ k , a ⟩ ⟨ ϕ k , b ∣ = k = 0 ∑ N − 1 ( a ∈ Σ ∑ ∣ a ⟩ ⊗ ∣ ϕ k , a ⟩ ) ( b ∈ Σ ∑ ⟨ b ∣ ⊗ ⟨ ϕ k , b ∣ ) = k = 0 ∑ N − 1 ∣ ψ k ⟩ ⟨ ψ k ∣ = J ( Φ )
したがって、選んだKraus行列は Φ \Phi Φ を正しく記述しています。
次に A 0 , … , A N − 1 A_0,\ldots,A_{N-1} A 0 , … , A N − 1 に課せられた条件を確認します。これは仮定 Tr Y ( J ( Φ ) ) = I X \operatorname{Tr}_{\mathsf{Y}}(J(\Phi)) = \mathbb{I}_{\mathsf{X}} Tr Y ( J ( Φ )) = I X (まだ使っていません)と等価であることが分かります。
示すべき関係は次の通りです。
( ∑ k = 0 N − 1 A k † A k ) T = Tr Y ( J ( Φ ) ) (2) \Biggl( \sum_{k = 0}^{N-1} A_k^{\dagger} A_k \Biggr)^{T} = \operatorname{Tr}_{\mathsf{Y}}(J(\Phi))
\tag{2} ( k = 0 ∑ N − 1 A k † A k ) T = Tr Y ( J ( Φ )) ( 2 )
(左辺は行列の転置 を意味します。)
左辺から始めると、まず
( ∑ k = 0 N − 1 A k † A k ) T = ( ∑ k = 0 N − 1 ∑ a , b ∈ Σ ∣ b ⟩ ⟨ ϕ k , b ∣ ϕ k , a ⟩ ⟨ a ∣ ) T = ∑ k = 0 N − 1 ∑ a , b ∈ Σ ⟨ ϕ k , b ∣ ϕ k , a ⟩ ∣ a ⟩ ⟨ b ∣ \begin{aligned}
\Biggl(\sum_{k = 0}^{N-1} A_k^{\dagger} A_k\Biggr)^T
& = \Biggl(\sum_{k = 0}^{N-1} \sum_{a,b\in\Sigma} \vert b \rangle \langle \phi_{k,b} \vert \phi_{k,a} \rangle \langle a \vert\Biggr)^T\\
& = \sum_{k = 0}^{N-1} \sum_{a,b\in\Sigma} \langle \phi_{k,b} \vert \phi_{k,a} \rangle \vert a \rangle \langle b \vert
\end{aligned} ( k = 0 ∑ N − 1 A k † A k ) T = ( k = 0 ∑ N − 1 a , b ∈ Σ ∑ ∣ b ⟩ ⟨ ϕ k , b ∣ ϕ k , a ⟩ ⟨ a ∣ ) T = k = 0 ∑ N − 1 a , b ∈ Σ ∑ ⟨ ϕ k , b ∣ ϕ k , a ⟩ ∣ a ⟩ ⟨ b ∣
と観察できます。最後の等号は、転置が線形で ∣ b ⟩ ⟨ a ∣ \vert b\rangle\langle a \vert ∣ b ⟩ ⟨ a ∣ を ∣ a ⟩ ⟨ b ∣ \vert a\rangle\langle b \vert ∣ a ⟩ ⟨ b ∣ に写すことから得られます。
等式の右辺に移ると、
J ( Φ ) = ∑ k = 0 N − 1 ∣ ψ k ⟩ ⟨ ψ k ∣ = ∑ k = 0 N − 1 ∑ a , b ∈ Σ ∣ a ⟩ ⟨ b ∣ ⊗ ∣ ϕ k , a ⟩ ⟨ ϕ k , b ∣ J(\Phi) = \sum_{k = 0}^{N-1} \vert \psi_k\rangle\langle\psi_k \vert
= \sum_{k = 0}^{N-1} \sum_{a,b\in\Sigma} \vert a\rangle \langle b \vert \otimes \vert\phi_{k,a}\rangle\langle \phi_{k,b} \vert J ( Φ ) = k = 0 ∑ N − 1 ∣ ψ k ⟩ ⟨ ψ k ∣ = k = 0 ∑ N − 1 a , b ∈ Σ ∑ ∣ a ⟩ ⟨ b ∣ ⊗ ∣ ϕ k , a ⟩ ⟨ ϕ k , b ∣
なので、
Tr Y ( J ( Φ ) ) = ∑ k = 0 N − 1 ∑ a , b ∈ Σ Tr ( ∣ ϕ k , a ⟩ ⟨ ϕ k , b ∣ ) ∣ a ⟩ ⟨ b ∣ = ∑ k = 0 N − 1 ∑ a , b ∈ Σ ⟨ ϕ k , b ∣ ϕ k , a ⟩ ∣ a ⟩ ⟨ b ∣ \begin{aligned}
\operatorname{Tr}_{\mathsf{Y}}(J(\Phi))
& = \sum_{k = 0}^{N-1} \sum_{a,b\in\Sigma} \operatorname{Tr}\bigl(\vert\phi_{k,a}\rangle\langle \phi_{k,b} \vert \bigr)\,
\vert a\rangle \langle b \vert\\
& = \sum_{k = 0}^{N-1} \sum_{a,b\in\Sigma} \langle \phi_{k,b} \vert \phi_{k,a} \rangle \vert a \rangle \langle b \vert
\end{aligned} Tr Y ( J ( Φ )) = k = 0 ∑ N − 1 a , b ∈ Σ ∑ Tr ( ∣ ϕ k , a ⟩ ⟨ ϕ k , b ∣ ) ∣ a ⟩ ⟨ b ∣ = k = 0 ∑ N − 1 a , b ∈ Σ ∑ ⟨ ϕ k , b ∣ ϕ k , a ⟩ ∣ a ⟩ ⟨ b ∣
となります。
同じ結果が得られたので、等式 ( 2 ) (2) ( 2 ) が確認されました。
仮定 Tr Y ( J ( Φ ) ) = I X \operatorname{Tr}_{\mathsf{Y}} (J(\Phi)) = \mathbb{I}_{\mathsf{X}} Tr Y ( J ( Φ )) = I X から、
( ∑ k = 0 N − 1 A k † A k ) T = I X \Biggl(\sum_{k = 0}^{N-1} A_k^{\dagger} A_k\Biggr)^T = \mathbb{I}_{\mathsf{X}} ( k = 0 ∑ N − 1 A k † A k ) T = I X
が成り立ち、単位行列は自身の転置であることから、必要な条件
∑ k = 0 N − 1 A k † A k = I X \sum_{k = 0}^{N-1} A_k^{\dagger} A_k = \mathbb{I}_{\mathsf{X}} k = 0 ∑ N − 1 A k † A k = I X
が真であることが示されました。
KrausからStinespring表現へ
次に、写像のKraus表現
Φ ( ρ ) = ∑ k = 0 N − 1 A k ρ A k † \Phi(\rho) = \sum_{k = 0}^{N-1} A_k \rho A_k^{\dagger} Φ ( ρ ) = k = 0 ∑ N − 1 A k ρ A k †
であって、
∑ k = 0 N − 1 A k † A k = I X \sum_{k = 0}^{N-1} A_k^{\dagger} A_k = \mathbb{I}_{\mathsf{X}} k = 0 ∑ N − 1 A k † A k = I X
を満たすものが与えられているとします。
目標は Φ \Phi Φ のStinespring表現を見つけることです。
まず、ゴミシステム G \mathsf{G} G の古典状態集合が { 0 , … , N − 1 } \{0,\ldots,N-1\} { 0 , … , N − 1 } となるように選びます。
ただし ( W , X ) (\mathsf{W},\mathsf{X}) ( W , X ) と ( G , Y ) (\mathsf{G},\mathsf{Y}) ( G , Y ) が同じサイズになるためには、n n n が m N m N m N を割り切る必要があります。その場合、W \mathsf{W} W の古典状態集合を { 0 , … , d − 1 } \{0,\ldots,d-1\} { 0 , … , d − 1 } (d = m N / n d = mN/n d = m N / n )とすることができます。
n n n 、m m m 、N N N の任意の選択に対して m N / n mN/n m N / n が整数になるとは限りませんが、Kraus表現においていくつかの追加的な値 k k k に対して A k = 0 A_k = 0 A k = 0 を設定することで N N N を任意に増やすことができます。
したがって、m N / n mN/n m N / n が整数であること(N N N が n / gcd ( n , m ) n/\operatorname{gcd}(n,m) n / gcd ( n , m ) の倍数であることと同値)を暗黙に仮定すれば、古典状態集合が { 0 , … , N − 1 } \{0,\ldots,N-1\} { 0 , … , N − 1 } である G \mathsf{G} G を選ぶことができます。
特に N = n m N = nm N = nm の場合、W \mathsf{W} W は m 2 m^2 m 2 個の古典状態を持てます。
あとは U U U を選ぶことが残っており、次のパターンに合わせて行います。
U = ( A 0 ? ⋯ ? A 1 ? ⋯ ? ⋮ ⋮ ⋱ ⋮ A N − 1 ? ⋯ ? ) U =
\begin{pmatrix}
A_{0} & \fbox{?} & \cdots & \fbox{?} \\[1mm]
A_{1} & \fbox{?} & \cdots & \fbox{?} \\[1mm]
\vdots & \vdots & \ddots & \vdots\\[1mm]
A_{N-1} & \fbox{?} & \cdots & \fbox{?}
\end{pmatrix} U = A 0 A 1 ⋮ A N − 1 ? ? ⋮ ? ⋯ ⋯ ⋱ ⋯ ? ? ⋮ ?
このパターンはブロック行列を表しています。各ブロック(A 0 , … , A N − 1 A_{0},\ldots,A_{N-1} A 0 , … , A N − 1 および疑問符で示されたブロックを含む)は m m m 行 n n n 列です。
ブロックの行数は N N N であり、ブロックの列数は d = m N / n d = mN/n d = m N / n となります。
より数式的に表現すると、U U U を次のように定義します。
U = ∑ k = 0 N − 1 ∑ j = 0 d − 1 ∣ k ⟩ ⟨ j ∣ ⊗ M k , j = ( M 0 , 0 M 0 , 1 ⋯ M 0 , d − 1 M 1 , 0 M 1 , 1 ⋯ M 1 , d − 1 ⋮ ⋮ ⋱ ⋮ M N − 1 , 0 M N − 1 , 1 ⋯ M N − 1 , d − 1 ) \begin{aligned}
U & = \sum_{k=0}^{N-1} \sum_{j=0}^{d-1} \vert k \rangle \langle j \vert \otimes M_{k,j} \\[4mm]
& = \begin{pmatrix}
M_{0,0} & M_{0,1} & \cdots & M_{0,d-1} \\[1mm]
M_{1,0} & M_{1,1} & \cdots & M_{1,d-1} \\[1mm]
\vdots & \vdots & \ddots & \vdots\\[1mm]
M_{N-1,0} & M_{N-1,1} & \cdots & M_{N-1,d-1}
\end{pmatrix}
\end{aligned} U = k = 0 ∑ N − 1 j = 0 ∑ d − 1 ∣ k ⟩ ⟨ j ∣ ⊗ M k , j = M 0 , 0 M 1 , 0 ⋮ M N − 1 , 0 M 0 , 1 M 1 , 1 ⋮ M N − 1 , 1 ⋯ ⋯ ⋱ ⋯ M 0 , d − 1 M 1 , d − 1 ⋮ M N − 1 , d − 1
各行列 M k , j M_{k,j} M k , j は m m m 行 n n n 列を持ち、特に k = 0 , … , N − 1 k = 0,\ldots,N-1 k = 0 , … , N − 1 に対して M k , 0 = A k M_{k,0} = A_k M k , 0 = A k とします。
これはユニタリ行列でなければならず、疑問符でラベルされたブロック、すなわち j > 0 j>0 j > 0 に対する M k , j M_{k,j} M k , j はこれを念頭に置いて選ばなければなりません。ただし、U U U をユニタリにすることを除けば、疑問符のブロックは証明には関係ありません。
U U U がユニタリであるという懸念は一時的に脇に置き、Stinespring表現において入力状態 ρ \rho ρ (システム X \mathsf{X} X )から出力状態(システム Y \mathsf{Y} Y )を与える式
Tr G ( U ( ∣ 0 ⟩ ⟨ 0 ∣ W ⊗ ρ ) U † ) \operatorname{Tr}_{\mathsf{G}} \bigl( U (\vert 0\rangle \langle 0 \vert_{\mathsf{W}} \otimes \rho)U^{\dagger}\bigr) Tr G ( U ( ∣0 ⟩ ⟨ 0 ∣ W ⊗ ρ ) U † )
に注目します。
次のように書き換えることができます。
U ( ∣ 0 ⟩ ⟨ 0 ∣ ⊗ ρ ) U † = U ( ∣ 0 ⟩ ⊗ I X ) ρ ( ⟨ 0 ∣ ⊗ I X ) U † , U(\vert 0\rangle\langle 0 \vert \otimes \rho)U^{\dagger}
= U(\vert 0\rangle\otimes\mathbb{I}_{\mathsf{X}}) \rho (\langle 0\vert \otimes \mathbb{I}_{\mathsf{X}}) U^{\dagger}, U ( ∣0 ⟩ ⟨ 0∣ ⊗ ρ ) U † = U ( ∣0 ⟩ ⊗ I X ) ρ (⟨ 0∣ ⊗ I X ) U † ,
U U U の選択から、
U ( ∣ 0 ⟩ ⊗ I X ) = ∑ k = 0 N − 1 ∣ k ⟩ ⊗ A k U(\vert 0\rangle\otimes\mathbb{I}_{\mathsf{X}}) =
\sum_{k = 0}^{N-1} \vert k\rangle \otimes A_k U ( ∣0 ⟩ ⊗ I X ) = k = 0 ∑ N − 1 ∣ k ⟩ ⊗ A k
が分かります。
したがって、
U ( ∣ 0 ⟩ ⟨ 0 ∣ ⊗ ρ ) U † = ∑ j , k = 0 N − 1 ∣ k ⟩ ⟨ j ∣ ⊗ A k ρ A j † U(\vert 0\rangle\langle 0 \vert \otimes \rho)U^{\dagger}
= \sum_{j,k = 0}^{N-1} \vert k\rangle\langle j\vert \otimes A_k \rho A_j^{\dagger} U ( ∣0 ⟩ ⟨ 0∣ ⊗ ρ ) U † = j , k = 0 ∑ N − 1 ∣ k ⟩ ⟨ j ∣ ⊗ A k ρ A j †
となり、
Tr G ( U ( ∣ 0 ⟩ ⟨ 0 ∣ W ⊗ ρ ) U † ) = ∑ j , k = 0 N − 1 Tr ( ∣ k ⟩ ⟨ j ∣ ) A k ρ A j † = ∑ k = 0 N − 1 A k ρ A k † = Φ ( ρ ) \begin{aligned}
\operatorname{Tr}_{\mathsf{G}} \bigl( U (\vert 0\rangle \langle 0 \vert_{\mathsf{W}} \otimes \rho) U^{\dagger}\bigr)
& = \sum_{j,k = 0}^{N-1} \operatorname{Tr}\bigl(\vert k\rangle\langle j\vert\bigr) \, A_k \rho A_j^{\dagger} \\
& = \sum_{k = 0}^{N-1} A_k \rho A_k^{\dagger} \\
& = \Phi(\rho)
\end{aligned} Tr G ( U ( ∣0 ⟩ ⟨ 0 ∣ W ⊗ ρ ) U † ) = j , k = 0 ∑ N − 1 Tr ( ∣ k ⟩ ⟨ j ∣ ) A k ρ A j † = k = 0 ∑ N − 1 A k ρ A k † = Φ ( ρ )
が得られます。
したがって、写像 Φ \Phi Φ の正しい表現が得られました。残るは U U U をユニタリに選べることの確認です。
上記のパターンで選んだ U U U の最初の n n n 列を考えます。
これらの列だけを取り出すと、ブロック行列
( A 0 A 1 ⋮ A N − 1 ) \begin{pmatrix}
A_0\\[1mm]
A_1\\[1mm]
\vdots\\[1mm]
A_{N-1}
\end{pmatrix} A 0 A 1 ⋮ A N − 1
が得られます。
列は n n n 本あり、それぞれ X \mathsf{X} X の古典状態に対応します。各 a ∈ Σ a\in\Sigma a ∈ Σ に対してこれらの列ベクトルを ∣ γ a ⟩ \vert \gamma_a \rangle ∣ γ a ⟩ と名付けます。
ブロック行列表現と対応するこれらのベクトルの式は次の通りです。
∣ γ a ⟩ = ∑ k = 0 N − 1 ∣ k ⟩ ⊗ A k ∣ a ⟩ \vert \gamma_a\rangle = \sum_{k = 0}^{N-1} \vert k\rangle \otimes A_k \vert a \rangle ∣ γ a ⟩ = k = 0 ∑ N − 1 ∣ k ⟩ ⊗ A k ∣ a ⟩
任意の a , b ∈ Σ a,b\in\Sigma a , b ∈ Σ に対するこれらのベクトル間の内積を計算します。
⟨ γ a ∣ γ b ⟩ = ∑ j , k = 0 N − 1 ⟨ k ∣ j ⟩ ⟨ a ∣ A k † A j ∣ b ⟩ = ⟨ a ∣ ( ∑ k = 0 N − 1 A k † A k ) ∣ b ⟩ \langle \gamma_a \vert \gamma_b \rangle =
\sum_{j,k = 0}^{N-1} \langle k \vert j \rangle \, \langle a \vert A_k^{\dagger} A_j \vert b\rangle =
\langle a \vert \Biggl( \sum_{k = 0}^{N-1} A_k^{\dagger} A_k \Biggr) \vert b\rangle ⟨ γ a ∣ γ b ⟩ = j , k = 0 ∑ N − 1 ⟨ k ∣ j ⟩ ⟨ a ∣ A k † A j ∣ b ⟩ = ⟨ a ∣ ( k = 0 ∑ N − 1 A k † A k ) ∣ b ⟩
仮定
∑ k = 0 m − 1 A k † A k = I X \sum_{k = 0}^{m-1} A_k^{\dagger} A_k = \mathbb{I}_{\mathsf{X}} k = 0 ∑ m − 1 A k † A k = I X
から、n n n 個の列ベクトル { ∣ γ a ⟩ : a ∈ Σ } \{\vert\gamma_a\rangle\,:\,a\in\Sigma\} { ∣ γ a ⟩ : a ∈ Σ } が正規直交系をなすことが分かります。
⟨ γ a ∣ γ b ⟩ = { 1 a = b 0 a ≠ b \langle \gamma_a \vert \gamma_b \rangle = \begin{cases}
1 & a = b\\
0 & a\neq b
\end{cases} ⟨ γ a ∣ γ b ⟩ = { 1 0 a = b a = b
これはすべての a , b ∈ Σ a,b\in\Sigma a , b ∈ Σ に対して成り立ちます。
これにより、U U U の残りの列を埋めてユニタリ行列にすることが可能であることが示されます。
特に、Gram-Schmidtの直交化法によって残りの列を選ぶことができます。
「量子情報の基礎」の量子Circuit レッスンで、状態判別問題の文脈で同様のことが行われています。
Stinespring表現から定義へ
最後の含意は 4 ⇒ \Rightarrow ⇒ 1 です。
すなわち、システムの対 ( W , X ) (\mathsf{W},\mathsf{X}) ( W , X ) を別の対 ( G , Y ) (\mathsf{G},\mathsf{Y}) ( G , Y ) に変換するユニタリ演算が存在することを仮定し、写像
Φ ( ρ ) = Tr G ( U ( ∣ 0 ⟩ ⟨ 0 ∣ W ⊗ ρ ) U † ) \Phi(\rho) = \operatorname{Tr}_{\mathsf{G}} \bigl( U (\vert 0\rangle \langle 0 \vert_{\mathsf{W}} \otimes \rho)U^{\dagger}\bigr) Φ ( ρ ) = Tr G ( U ( ∣0 ⟩ ⟨ 0 ∣ W ⊗ ρ ) U † )
が有効なチャネルであることを示すことが目標です。
その形から Φ \Phi Φ が線形であることは明らかで、残るのは常に密度行列を密度行列に変換することの確認です。
これは比較的単純であり、重要な点はすでに議論しています。
具体的には、複合システム ( Z , X ) (\mathsf{Z},\mathsf{X}) ( Z , X ) の密度行列 σ \sigma σ から始め、追加のワークスペースシステム W \mathsf{W} W を加えると、確かに密度行列が得られます。
システムを ( W , Z , X ) (\mathsf{W},\mathsf{Z},\mathsf{X}) ( W , Z , X ) の順に並べ替えると、この状態は
∣ 0 ⟩ ⟨ 0 ∣ W ⊗ σ \vert 0\rangle\langle 0\vert_{\mathsf{W}} \otimes \sigma ∣0 ⟩ ⟨ 0 ∣ W ⊗ σ
と書けます。
次にユニタリ演算 U U U を適用します。これはすでに議論したように有効なチャネルであり、密度行列を密度行列に写します。
最後に、密度行列の部分トレースは別の密度行列です。
別の言い方をすれば、次の各操作がそれぞれ有効なチャネルであることを観察します。
初期化されたワークスペースシステムを導入すること。
ユニタリ演算を行うこと。
システムをトレースアウトすること。
最後に、チャネルの合成はまた別のチャネルになります。これは定義から直ちに分かりますが、それ自体として注目に値する事実でもあります。
これで最後の含意の証明が完了し、節の冒頭に列挙した4つの命題の等価性が確立されました。