メインコンテンツへスキップ
このプロジェクトはベータ版です。バグを見つけましたか?GitHub Issueを開く。質問やアイデアがありますか?Discussionsで質問する — フィードバックをお待ちしています!
このページの翻訳品質はいかがですか?

表現の等価性

これまでに、量子チャネルを数学的に表現する3つの異なる方法、すなわちStinespring表現、Kraus表現、Choi表現について説明しました。 また、チャネルの定義として、チャネルとは密度行列を常に密度行列へと変換する線形写像であり、複合システムの一部にのみチャネルを適用する場合でもこの性質が成り立つことを確認しました。 このレッスンの残りは、この3つの表現が等価であり、定義を正確に捉えていることの数学的証明に充てられます。

証明の概要​

私たちの目標は、4つの命題の等価性を示すことです。まずそれらを正確に書き下しましょう。 4つの命題はすべて、レッスン全体で用いてきた同じ規約に従います。すなわち、Φ\Phi は正方行列から正方行列への線形写像であり、入力行列の行と列はシステム X\mathsf{X}(入力システム)の古典状態に対応し、出力行列の行と列はシステム Y\mathsf{Y}(出力システム)の古典状態に対応しています。

  1. Φ\Phi は X\mathsf{X} から Y\mathsf{Y} へのチャネルである。すなわち、Φ\Phi はより大きな複合システムの一部に作用する場合でも、常に密度行列を密度行列に変換する。

  2. Choi行列 J(Φ)J(\Phi) は半正定値であり、条件 Tr⁡Y(J(Φ))=IX\operatorname{Tr}_{\mathsf{Y}}(J(\Phi)) = \mathbb{I}_{\mathsf{X}} を満たす。

  3. Φ\Phi のKraus表現が存在する。すなわち、すべての入力 ρ\rho に対して等式 Φ(ρ)=∑k=0N−1AkρAk†\Phi(\rho) = \sum_{k = 0}^{N-1} A_k \rho A_k^{\dagger} が成り立ち、かつ条件 ∑k=0N−1Ak†Ak=IX\sum_{k = 0}^{N-1} A_k^{\dagger} A_k = \mathbb{I}_{\mathsf{X}} を満たす行列 A0,…,AN−1A_0,\ldots,A_{N-1} が存在する。

  4. Φ\Phi のStinespring表現が存在する。すなわち、システム W\mathsf{W} と G\mathsf{G} が存在し、(W,X)(\mathsf{W},\mathsf{X}) と (G,Y)(\mathsf{G},\mathsf{Y}) の古典状態の数が等しく、さらに (W,X)(\mathsf{W},\mathsf{X}) から (G,Y)(\mathsf{G},\mathsf{Y}) へのユニタリ演算を表すユニタリ行列 UU が存在して、Φ(ρ)=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) が成り立つ。

証明の方針は、含意の連鎖を示すことです。 リスト中の1番目の命題が2番目を含意し、2番目が3番目を含意し、3番目が4番目を含意し、4番目が1番目を含意することを証明します。 これにより、4つの命題すべてが等価である——すなわち、与えられた Φ\Phi に対してすべてが真かすべてが偽である——ことが確立されます。なぜなら、含意を推移的にたどることで任意の1つの命題から別の任意の命題へと到達できるからです。

これは複数の命題の等価性を証明する際の一般的な戦略であり、そのような文脈で有用な工夫は、含意をできるだけ証明しやすい方向に設定することです。 実際ここでもそうなっており、4つの含意のうち2つはすでに扱っています。

チャネルからChoi行列へ​

上記の番号で参照すると、最初に証明すべき含意は 1 ⇒\Rightarrow 2 です。 この含意はチャネルのChoi状態の文脈ですでに議論されています。 ここでは数学的な詳細をまとめます。

入力システム X\mathsf{X} の古典状態集合を Σ\Sigma とし、n=∣Σ∣n = \vert\Sigma\vert とします。 Φ\Phi が X\mathsf{X} の2つのコピーからなるシステムの第2のコピーに適用される状況を考えます。2つのコピーが合わせて状態

∣ψ⟩=1n∑a∈Σ∣a⟩⊗∣a⟩,\vert \psi \rangle = \frac{1}{\sqrt{n}} \sum_{a \in \Sigma} \vert a \rangle \otimes \vert a \rangle,

にある場合、密度行列として書くと

∣ψ⟩⟨ψ∣=1n∑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.

となります。

(Id⁡⊗ Φ)(∣ψ⟩⟨ψ∣)=1n∑a,b=0n−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},

そして Φ\Phi がチャネルであるという仮定により、これは密度行列でなければなりません。 すべての密度行列と同様に半正定値でなければならず、半正定値行列に正の実数を掛けても別の半正定値行列が得られるため、J(Φ)≥0J(\Phi) \geq 0 となります。

さらに、Φ\Phi がチャネルであるという仮定のもとでは、トレースを保存しなければなりません。したがって

Tr⁡Y(J(Φ))=∑a,b∈ΣTr⁡(Φ(∣a⟩⟨b∣)) ∣a⟩⟨b∣=∑a,b∈ΣTr⁡(∣a⟩⟨b∣) ∣a⟩⟨b∣=∑a∈Σ∣a⟩⟨a∣=IX.\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}

ChoiからKraus表現へ​

2番目の含意は、リストの文の番号を再び参照すると、2 ⇒\Rightarrow 3 です。 明確にしておくと、私たちは他の文を無視しています — 特に Φ\Phi がチャネルであるという仮定を置くことはできません。 私たちが使えるのは、Φ\Phi が線形写像であり、そのチョイ表現が J(Φ)≥0J(\Phi) \geq 0 かつ Tr⁡Y(J(Φ))=IX\operatorname{Tr}_{\mathsf{Y}} (J(\Phi)) = \mathbb{I}_{\mathsf{X}} を満たすということだけです。

しかし、これだけで、Φ\Phi が次のクラウス表現を持つと結論するのに十分です。

Φ(ρ)=∑k=0N−1AkρAk†\Phi(\rho) = \sum_{k = 0}^{N-1} A_k \rho A_k^{\dagger}

次の条件が満たされるためには

∑k=0N−1Ak†Ak=IX\sum_{k = 0}^{N-1} A_k^{\dagger} A_k = \mathbb{I}_{\mathsf{X}}

J(Φ)J(\Phi) が半正定値であるという決定的に重要な仮定から始めます。これは、

私たちは、J(Φ)J(\Phi) が半正定値であるという極めて重要な仮定から始めます。これは、次の形で表現できることを意味します。

J(Φ)=∑k=0N−1∣ψk⟩⟨ψk∣(1)J(\Phi) = \sum_{k = 0}^{N-1} \vert \psi_k \rangle \langle \psi_k \vert \tag{1}

ベクトル ∣ψ0⟩,…,∣ψN−1⟩\vert\psi_0\rangle,\ldots,\vert\psi_{N-1}\rangle の選び方によっては。 一般に、これを行う方法は複数あります — 実際、これは Φ\Phi のクラウス表現を選ぶ際の自由度を直接反映しています。

このような式を得る1つの方法は、まずスペクトル定理を使って次のように書くことです。

J(Φ)=∑k=0N−1λk∣γk⟩⟨γk∣,J(\Phi) = \sum_{k = 0}^{N-1} \lambda_k \vert \gamma_k \rangle \langle \gamma_k \vert,

ここで λ0,…,λN−1\lambda_0,\ldots,\lambda_{N-1} は J(Φ)J(\Phi) の固有値であり(J(Φ)J(\Phi) が半正定値であるため必然的に非負の実数です)、∣γ0⟩,…,∣γN−1⟩\vert\gamma_0\rangle,\ldots,\vert\gamma_{N-1}\rangle は固有値 λ0,…,λN−1\lambda_0,\ldots,\lambda_{N-1} に対応する単位固有ベクトルです。

固有値の選択には自由度がありません(順序を除いて)が、固有ベクトルの選択には自由度があることに注意してください。特に多重度が1より大きい固有値がある場合はそうです。 したがって、これは J(Φ)J(\Phi) の一意な表現ではありません — 私たちは単にそのような表現の1つを持っていると仮定しているだけです。 いずれにせよ、固有値は非負の実数であるため、非負の平方根を持ちます。したがって、次を選ぶことができます。

∣ψk⟩=λk∣γk⟩\vert\psi_k\rangle = \sqrt{\lambda_k} \vert \gamma_k\rangle

各 k=0,…,N−1k = 0,\ldots,N-1 に対して、これにより式 (1)(1) の形の式が得られます。

しかし、式 (1)(1) がこのようにスペクトル分解から来ることは必須ではなく、特にベクトル ∣ψ0⟩,…,∣ψN−1⟩\vert\psi_0\rangle,\ldots,\vert\psi_{N-1}\rangle は一般に直交している必要はありません。 ただし、望むならこれらのベクトルを直交するように選ぶことができ、さらに NN を nmnm より大きくする必要は決してないことは注目に値します (nn と mm はそれぞれ X\mathsf{X} と Y\mathsf{Y} の古典状態の数を表すことを思い出してください)。

次に、ベクトル ∣ψ0⟩,…,∣ψN−1⟩\vert\psi_0\rangle,\ldots,\vert\psi_{N-1}\rangle のそれぞれをさらに次のように分解できます。

∣ψk⟩=∑a∈Σ∣a⟩⊗∣ϕk,a⟩,\vert\psi_k\rangle = \sum_{a\in\Sigma} \vert a\rangle \otimes \vert \phi_{k,a}\rangle,

ここで、ベクトル {∣ϕk,a⟩}\{ \vert \phi_{k,a}\rangle \} は Y\mathsf{Y} の古典状態に対応するエントリを持ち、次の方程式によって明示的に決定できます。

∣ϕk,a⟩=(⟨a∣⊗IY)∣ψk⟩\vert \phi_{k,a}\rangle = \bigl( \langle a \vert \otimes \mathbb{I}_{\mathsf{Y}}\bigr) \vert \psi_k\rangle

各 a∈Σa\in\Sigma および k=0,…,N−1k=0,\ldots,N-1 に対して。 ∣ψ0⟩,…,∣ψN−1⟩\vert\psi_0\rangle,\ldots,\vert\psi_{N-1}\rangle は必ずしも単位ベクトルではありませんが、これはペア (X,Y)(\mathsf{X},\mathsf{Y}) の量子状態ベクトルが与えられたときに、システム X\mathsf{X} 上で標準基底測定が行われた場合に何が起こるかを解析するために使うのと同じプロセスです。

さて、この証明のこの部分を機能させるトリックにたどり着きました。 次の方程式に従って、クラウス行列 A0,…,AN−1A_0,\ldots,A_{N-1} を定義します。

Ak=∑a∈Σ∣ϕk,a⟩⟨a∣A_k = \sum_{a\in\Sigma} \vert \phi_{k,a}\rangle\langle a \vert

この公式は純粋に記号的に考えることができます。∣a⟩\vert a\rangle は実質的にひっくり返って ⟨a∣\langle a\vert を形成し、右辺に移動して行列を形成します。 証明を検証する目的では、この公式だけで十分です。

しかし、ベクトル ∣ψk⟩\vert\psi_k\rangle と行列 AkA_k の間には単純で直感的な関係があります。それは、AkA_k をベクトル化すると ∣ψk⟩\vert\psi_k\rangle が得られるということです。 AkA_k をベクトル化するとは、列を(一番左の列を一番上にして一番右の列を一番下にして)積み重ねてベクトルを形成することを意味します。 例えば、X\mathsf{X} と Y\mathsf{Y} が両方ともキュビットで、ある kk の選択に対して次のようになるとします。

∣ψ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}

then

Ak=α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}

(注意:行列のベクトル化は、時々わずかに異なる方法で定義されることがあります。それは、行列の行を転置して積み重ね、列ベクトルを形成するという方法です。)

まず、このクラウス行列の選択が写像 Φ\Phi を正しく記述していることを確認し、その後もう一方の必要条件を確認します。 話を整理するために、新しい写像 Ψ\Psi を次のように定義しましょう。

Ψ(ρ)=∑k=0N−1AkρAk†\Psi(\rho) = \sum_{k = 0}^{N-1} A_k \rho A_k^{\dagger}

したがって、私たちの目標は Ψ=Φ\Psi = \Phi を検証することです。

これを行う方法は、これらの写像のチョイ表現を比較することです。 チョイ表現は忠実であるため、J(Φ)=J(Ψ)J(\Phi) = J(\Psi) のとき、かつそのときに限り Ψ=Φ\Psi = \Phi となります。 この時点で、次の式を使って J(Ψ)J(\Psi) を単純に計算できます。

∣ψk⟩=∑a∈Σ∣a⟩⊗∣ϕk,a⟩andAk=∑a∈Σ∣ϕk,a⟩⟨a∣\vert\psi_k\rangle = \sum_{a\in\Sigma} \vert a\rangle \otimes \vert \phi_{k,a}\rangle \quad\text{and}\quad A_k = \sum_{a\in\Sigma} \vert \phi_{k,a}\rangle\langle a \vert

テンソル積の双線形性と合わせて簡略化します。

J(Ψ)=∑a,b∈Σ∣a⟩⟨b∣⊗∑k=0N−1Ak∣a⟩⟨b∣Ak†=∑a,b∈Σ∣a⟩⟨b∣⊗∑k=0N−1∣ϕk,a⟩⟨ϕk,b∣=∑k=0N−1(∑a∈Σ∣a⟩⊗∣ϕk,a⟩)(∑b∈Σ⟨b∣⊗⟨ϕk,b∣)=∑k=0N−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}

したがって、私たちのクラウス行列は Φ\Phi を正しく記述しています。

残るは A0,…,AN−1A_0,\ldots,A_{N-1} に対する必要条件をチェックすることです。これは、まだ使っていない仮定 Tr⁡Y(J(Φ))=IX\operatorname{Tr}_{\mathsf{Y}}(J(\Phi)) = \mathbb{I}_{\mathsf{X}} と同値であることが分かります。 私たちが示すのは、次の関係です。

(∑k=0N−1Ak†Ak)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=0N−1Ak†Ak)T=(∑k=0N−1∑a,b∈Σ∣b⟩⟨ϕk,b∣ϕk,a⟩⟨a∣)T=∑k=0N−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}

最後の等式は、転置が線形であり、∣b⟩⟨a∣\vert b\rangle\langle a \vert を ∣a⟩⟨b∣\vert a\rangle\langle b \vert に写すという事実から導かれます。

方程式の右辺に移ると、次のようになります。

J(Φ)=∑k=0N−1∣ψk⟩⟨ψk∣=∑k=0N−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

したがって、

Tr⁡Y(J(Φ))=∑k=0N−1∑a,b∈ΣTr⁡(∣ϕk,a⟩⟨ϕk,b∣) ∣a⟩⟨b∣=∑k=0N−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}

同じ結果が得られたので、方程式 (2)(2) が検証されました。 仮定 Tr⁡Y(J(Φ))=IX\operatorname{Tr}_{\mathsf{Y}} (J(\Phi)) = \mathbb{I}_{\mathsf{X}} により、次のことが分かります。

(∑k=0N−1Ak†Ak)T=IX\Biggl(\sum_{k = 0}^{N-1} A_k^{\dagger} A_k\Biggr)^T = \mathbb{I}_{\mathsf{X}}

が真であることが示されました。

∑k=0N−1Ak†Ak=IX\sum_{k = 0}^{N-1} A_k^{\dagger} A_k = \mathbb{I}_{\mathsf{X}}

KrausからStinespring表現へ​

次に、写像のKraus表現

Φ(ρ)=∑k=0N−1AkρAk†\Phi(\rho) = \sum_{k = 0}^{N-1} A_k \rho A_k^{\dagger}

であって、

∑k=0N−1Ak†Ak=IX.\sum_{k = 0}^{N-1} A_k^{\dagger} A_k = \mathbb{I}_{\mathsf{X}}.

目標は Φ\Phi のStinespring表現を見つけることです。

まず、ゴミシステム G\mathsf{G} の古典状態集合が {0,…,N−1}\{0,\ldots,N-1\} となるように選びます。 ただし (W,X)(\mathsf{W},\mathsf{X}) と (G,Y)(\mathsf{G},\mathsf{Y}) が同じサイズになるためには、nn が mNm N を割り切る必要があります。その場合、W\mathsf{W} の古典状態集合を {0,…,d−1}\{0,\ldots,d-1\}(d=mN/nd = mN/n)とすることができます。

nn、mm、NN の任意の選択に対して、mN/nmN/n が整数であるとは限らないため、実際には G\mathsf{G} をその古典状態集合が {0,…,N−1}\{0,\ldots,N-1\} になるように自由に選べるわけではありません。 しかし、望むだけの追加の kk の値に対して Ak=0A_k = 0 を選ぶことで、クラウス表現において NN を常に任意に増やすことができます。

したがって、mN/nmN/n が整数であること(NN が n/gcd⁡(n,m)n/\operatorname{gcd}(n,m) の倍数であることと同値)を暗黙に仮定すれば、古典状態集合が {0,…,N−1}\{0,\ldots,N-1\} である G\mathsf{G} を選ぶことができます。 特に N=nmN = nm の場合、W\mathsf{W} は m2m^2 個の古典状態を持てます。

あとは UU を選ぶことが残っており、次のパターンに合わせて行います。

U=(A0?⋯?A1?⋯?⋮⋮⋱⋮AN−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}

このパターンはブロック行列を表しています。各ブロック(A0,…,AN−1A_{0},\ldots,A_{N-1} および疑問符で示されたブロックを含む)は mm 行 nn 列です。 ブロックの行数は NN であり、ブロックの列数は d=mN/nd = mN/n となります。

より数式的に表現すると、UU を次のように定義します。

U=∑k=0N−1∑j=0d−1∣k⟩⟨j∣⊗Mk,j=(M0,0M0,1⋯M0,d−1M1,0M1,1⋯M1,d−1⋮⋮⋱⋮MN−1,0MN−1,1⋯MN−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}

各行列 Mk,jM_{k,j} は mm 行 nn 列を持ち、特に k=0,…,N−1k = 0,\ldots,N-1 に対して Mk,0=AkM_{k,0} = A_k とします。

これはユニタリ行列でなければならず、疑問符でラベルされたブロック、すなわち j>0j>0 に対する Mk,jM_{k,j} はこれを念頭に置いて選ばなければなりません。ただし、UU をユニタリにすることを除けば、疑問符のブロックは証明には関係ありません。

UU がユニタリであるという懸念を一時的に無視して、次の式に注目しましょう。

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)

これは、私たちのスタインスプリング表現に対して、X\mathsf{X} の入力状態 ρ\rho が与えられたときの Y\mathsf{Y} の出力状態を記述します。 代わりに次のように書くことができます。

U(∣0⟩⟨0∣⊗ρ)U†=U(∣0⟩⊗IX)ρ(⟨0∣⊗IX)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},

UU の選択から、

U(∣0⟩⊗IX)=∑k=0N−1∣k⟩⊗Ak.U(\vert 0\rangle\otimes\mathbb{I}_{\mathsf{X}}) = \sum_{k = 0}^{N-1} \vert k\rangle \otimes A_k.

が分かります。

U(∣0⟩⟨0∣⊗ρ)U†=∑j,k=0N−1∣k⟩⟨j∣⊗AkρAj†,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},

and so

Tr⁡G(U(∣0⟩⟨0∣W⊗ρ)U†)=∑j,k=0N−1Tr⁡(∣k⟩⟨j∣) AkρAj†=∑k=0N−1AkρAk†=Φ(ρ).\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}

したがって、写像 Φ\Phi の正しい表現が得られました。残るは UU をユニタリに選べることを検証することです。

上記のパターンに従って選ばれたとき、UU の最初の nn 列を考えます。 これらの列だけを取ると、ブロック行列が得られます。

(A0A1⋮AN−1).\begin{pmatrix} A_0\\[1mm] A_1\\[1mm] \vdots\\[1mm] A_{N-1} \end{pmatrix}.

X\mathsf{X} の各古典状態に対して1つずつ、nn 個の列があります。ベクトルとして、各 a∈Σa\in\Sigma に対してこれらの列を ∣γa⟩\vert \gamma_a \rangle と名付けましょう。 上記のブロック行列表現と一致させることができるこれらのベクトルの公式を以下に示します。

∣γa⟩=∑k=0N−1∣k⟩⊗Ak∣a⟩\vert \gamma_a\rangle = \sum_{k = 0}^{N-1} \vert k\rangle \otimes A_k \vert a \rangle

それでは、これらのベクトルのうちの任意の2つ、つまり a,b∈Σa,b\in\Sigma の任意の選択に対応するものの内積を計算してみましょう。

⟨γa∣γb⟩=∑j,k=0N−1⟨k∣j⟩ ⟨a∣Ak†Aj∣b⟩=⟨a∣(∑k=0N−1Ak†Ak)∣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∈Σa,b\in\Sigma に対するこれらのベクトル間の内積を計算します。

∑k=0m−1Ak†Ak=IX\sum_{k = 0}^{m-1} A_k^{\dagger} A_k = \mathbb{I}_{\mathsf{X}}

nn 個の列ベクトル {∣γa⟩ : a∈Σ}\{\vert\gamma_a\rangle\,:\,a\in\Sigma\} が正規直交集合を形成することが分かります。

⟨γa∣γb⟩={1a=b0a≠b\langle \gamma_a \vert \gamma_b \rangle = \begin{cases} 1 & a = b\\ 0 & a\neq b \end{cases}

から、nn 個の列ベクトル {∣γa⟩ : a∈Σ}\{\vert\gamma_a\rangle\,:\,a\in\Sigma\} が正規直交系をなすことが分かります。

これは、UU の残りの列を埋めることで、それがユニタリ行列になり得ることを意味します。 特に、グラム・シュミットの直交化過程を使って残りの列を選ぶことができます。 「Basics of quantum information」のQuantum circuitsレッスンで、状態識別問題の文脈で同様のことが行われました。

Stinespring表現から定義へ​

最後の含意は 4 ⇒\Rightarrow 1 です。 つまり、一対のシステム (W,X)(\mathsf{W},\mathsf{X}) を一対のシステム (G,Y)(\mathsf{G},\mathsf{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)

が有効なチャネルであることを結論することです。 その形から、Φ\Phi が線形であることは明らかです。残るは、それが常に密度行列を密度行列に変換することを検証することです。 これはかなり単純明快であり、重要なポイントはすでに議論済みです。

特に、複合システム (Z,X)(\mathsf{Z},\mathsf{X}) の密度行列 σ\sigma から始めて、追加の作業用システム W\mathsf{W} を付け加えると、確実に密度行列が残ります。 便宜上システム (W,Z,X)(\mathsf{W},\mathsf{Z},\mathsf{X}) を並べ替えると、この状態を次のように書くことができます。

∣0⟩⟨0∣W⊗σ.\vert 0\rangle\langle 0\vert_{\mathsf{W}} \otimes \sigma.

次にユニタリ演算 UU を適用します。これはすでに議論したように有効なチャネルであり、密度行列を密度行列に写します。 最後に、密度行列の部分トレースは別の密度行列です。

別の言い方をすれば、次の各操作がそれぞれ有効なチャネルであることを観察します。

  1. 初期化されたワークスペースシステムを導入すること。

  2. ユニタリ演算を行うこと。

  3. システムをトレースアウトすること。

最後に、チャネルの合成はまた別のチャネルになります。これは定義から直ちに分かりますが、それ自体として注目に値する事実でもあります。

これで最後の含意の証明が完了し、節の冒頭に列挙した4つの命題の等価性が確立されました。