メインコンテンツへスキップ

純粋化

純粋化の定義

まず、純粋化の厳密な数学的定義から始めましょう。

定義

X\mathsf{X} を密度行列 ρ\rho で表される状態にある系とし、ψ\vert\psi\rangleY\mathsf{Y} をトレースアウトすることで ρ\rho が得られるような対 (X,Y)(\mathsf{X},\mathsf{Y}) の量子状態ベクトルとします。

ρ=TrY(ψψ).\rho = \operatorname{Tr}_{\mathsf{Y}} \bigl( \vert \psi\rangle\langle\psi\vert\bigr).

このとき、状態ベクトル ψ\vert\psi\rangleρ\rho純粋化と呼ばれます。

定義の等式が成り立つとき、量子状態ベクトルではなく密度行列で表された純粋状態 ψψ\vert\psi\rangle\langle\psi\vert もまた ρ\rho の純粋化と呼ばれることが一般的ですが、本レッスンでは主に量子状態ベクトルを指す用語として使用します。

純粋化という用語は、系の順序が逆になる場合、系や状態の名前が異なる場合(もちろん)、そして3つ以上の系がある場合にも、より広い意味で使用されます。 例えば、ψ\vert \psi \rangle が複合系 (A,B,C)(\mathsf{A},\mathsf{B},\mathsf{C}) の純粋状態を表す量子状態ベクトルであり、等式

ρ=TrB(ψψ)\rho = \operatorname{Tr}_{\mathsf{B}} \bigl(\vert\psi\rangle\langle\psi\vert\bigr)

が系 (A,C)(\mathsf{A},\mathsf{C}) の状態を表す密度行列 ρ\rho に対して成り立つ場合、ψ\vert\psi\rangle はやはり ρ\rho の純粋化と呼ばれます。

ただし、このレッスンでは定義で説明した特定の形式に焦点を当てます。 この定義による純粋化に関する性質や事実は、通常、系を2つの複合系に並べ直して分割することで3つ以上の系に一般化できます。片方が X\mathsf{X} の役割を担い、もう一方が Y\mathsf{Y} の役割を担います。

純粋化の存在

X\mathsf{X}Y\mathsf{Y} を任意の2つの系とし、ρ\rhoX\mathsf{X} の与えられた状態とします。 系 Y\mathsf{Y} が十分に大きければ、ρ\rho純粋化する — つまり ψ\vert\psi\rangleρ\rho の純粋化であるということ — (X,Y)(\mathsf{X},\mathsf{Y}) の量子状態ベクトル ψ\vert\psi\rangle が存在することを示します。 特に、Y\mathsf{Y}X\mathsf{X} と同じかそれ以上の古典状態を持つならば、すべての状態 ρ\rho に対してこのような形の純粋化が必ず存在します。 一部の状態 ρ\rho では Y\mathsf{Y} の古典状態がより少なくて済みます。一般に、ρ\rho を純粋化する (X,Y)(\mathsf{X},\mathsf{Y}) の量子状態ベクトルが存在するために必要かつ十分な Y\mathsf{Y} の古典状態数は rank(ρ)\operatorname{rank}(\rho) です。

まず、任意の正の整数 nn に対して、ρ\rhonn 個の純粋状態の凸結合として表す式を考えます。

ρ=a=0n1paϕaϕa\rho = \sum_{a = 0}^{n-1} p_a \vert\phi_a\rangle\langle\phi_a\vert

この式において、(p0,,pn1)(p_0,\ldots,p_{n-1}) は確率ベクトルで、ϕ0,,ϕn1\vert\phi_0\rangle,\ldots,\vert\phi_{n-1}\rangleX\mathsf{X} の量子状態ベクトルです。

このような表現を得る一つの方法はスペクトル定理を用いることで、その場合 nnX\mathsf{X} の古典状態数、p0,,pn1p_0,\ldots,p_{n-1}ρ\rho の固有値、ϕ0,,ϕn1\vert\phi_0\rangle,\ldots,\vert\phi_{n-1}\rangle はこれらの固有値に対応する正規直交固有ベクトルになります。

実際には、和に ρ\rho のゼロ固有値に対応する項を含める必要はなく、代わりに n=rank(ρ)n = \operatorname{rank}(\rho)p0,,pn1p_0,\ldots,p_{n-1}ρ\rho の非ゼロ固有値として選ぶことができます。 これは上記の形で ρ\rho を表現するために必要な nn の最小値です。

明確にしておくと、ρ\rho を純粋状態の凸結合として表す際に、スペクトル定理から得られる表現を使う必要はありません — これはそのような表現を得る一つの方法に過ぎません。 特に、nn は任意の正の整数でよく、単位ベクトル ϕ0,,ϕn1\vert\phi_0\rangle,\ldots,\vert\phi_{n-1}\rangle が直交している必要はなく、確率 p0,,pn1p_0,\ldots,p_{n-1}ρ\rho の固有値である必要もありません。

次のようにして ρ\rho の純粋化を特定できます。

ψ=a=0n1paϕaa\vert\psi\rangle = \sum_{a = 0}^{n-1} \sqrt{p_a} \, \vert\phi_a\rangle \otimes \vert a \rangle

ここでは Y\mathsf{Y} の古典状態が 0,,n10,\ldots,n-1 を含むと仮定しています。 含まない場合は、0,,n10,\ldots,n-1 の代わりに Y\mathsf{Y}nn 個の異なる古典状態を任意に選んで代入できます。 これが確かに ρ\rho の純粋化であることの確認は、部分トレースを計算するだけの簡単な作業で、以下の2つの等価な方法で行えます。

TrY(ψψ)=a=0n1(IXa)ψψ(IXa)=a=0n1paϕaϕa=ρ\operatorname{Tr}_{\mathsf{Y}} \bigl(\vert\psi\rangle\langle\psi\vert\bigr) = \sum_{a = 0}^{n-1} (\mathbb{I}_{\mathsf{X}} \otimes \langle a\vert) \vert\psi\rangle\langle\psi\vert (\mathbb{I}_{\mathsf{X}} \otimes \vert a\rangle) = \sum_{a = 0}^{n-1} p_a \vert\phi_a\rangle\langle\phi_a\vert = \rho TrY(ψψ)=a,b=0n1papbϕaϕbTr(ab)=a=0n1paϕaϕa=ρ\operatorname{Tr}_{\mathsf{Y}} \bigl(\vert\psi\rangle\langle\psi\vert\bigr) = \sum_{a,b = 0}^{n-1} \sqrt{p_a} \sqrt{p_b} \, \vert\phi_a\rangle\langle \phi_b\vert \, \operatorname{Tr}(\vert a \rangle \langle b \vert) = \sum_{a = 0}^{n-1} p_a \, \vert\phi_a\rangle\langle \phi_a\vert = \rho

より一般的には、任意の正規直交ベクトルの組 {γ0,,γn1}\{\vert\gamma_0\rangle,\ldots,\vert\gamma_{n-1}\rangle\} に対して、量子状態ベクトル

ψ=a=0n1paϕaγa\vert\psi\rangle = \sum_{a = 0}^{n-1} \sqrt{p_a} \, \vert\phi_a\rangle \otimes \vert \gamma_a \rangle

ρ\rho の純粋化です。

X\mathsf{X}Y\mathsf{Y} がどちらも Qubit で、

ρ=(34141414)\rho = \begin{pmatrix} \frac{3}{4} & \frac{1}{4}\\[2mm] \frac{1}{4} & \frac{1}{4} \end{pmatrix}

X\mathsf{X} の状態を表す密度行列であるとします。

スペクトル定理を使って ρ\rho を次のように表せます。

ρ=cos2(π/8)ψπ/8ψπ/8+sin2(π/8)ψ5π/8ψ5π/8,\rho = \cos^2(\pi/8) \vert \psi_{\pi/8}\rangle\langle\psi_{\pi/8}\vert + \sin^2(\pi/8) \vert \psi_{5\pi/8}\rangle\langle\psi_{5\pi/8}\vert,

ここで ψθ=cos(θ)0+sin(θ)1\vert \psi_{\theta} \rangle = \cos(\theta) \vert 0\rangle + \sin(\theta)\vert 1\rangle です。 量子状態ベクトル

cos(π/8)ψπ/80+sin(π/8)ψ5π/81\cos(\pi/8) \vert \psi_{\pi/8}\rangle \otimes \vert 0\rangle + \sin(\pi/8) \vert \psi_{5\pi/8}\rangle \otimes \vert 1\rangle

は対 (X,Y)(\mathsf{X},\mathsf{Y}) の純粋状態を記述するため、ρ\rho の純粋化となります。

別の方法として、次のように書くこともできます。

ρ=1200+12++.\rho = \frac{1}{2} \vert 0\rangle\langle 0\vert + \frac{1}{2} \vert +\rangle\langle +\vert.

これは純粋状態の凸結合ですが、0\vert 0\rangle+\vert +\rangle が直交しておらず 1/21/2ρ\rho の固有値でもないため、スペクトル分解ではありません。 それでも、量子状態ベクトル

1200+12+1\frac{1}{\sqrt{2}} \vert 0 \rangle \otimes \vert 0\rangle + \frac{1}{\sqrt{2}} \vert + \rangle \otimes \vert 1\rangle

ρ\rho の純粋化です。

シュミット分解

次に、シュミット分解について説明します。これは、ある特定の形を持つの系の量子状態ベクトルの表現です。 シュミット分解は純粋化と密接に関係しており、それ自体として非常に有用です。 実際、対の系の量子状態ベクトル ψ\vert\psi\rangle について考えるとき、最初のステップはしばしばこの状態のシュミット分解を特定または考察することです。

定義

ψ\vert \psi\rangle を対の系 (X,Y)(\mathsf{X},\mathsf{Y}) の与えられた量子状態ベクトルとします。ψ\vert\psi\rangleシュミット分解とは、次の形の表現です。

ψ=a=0r1paxaya,\vert \psi\rangle = \sum_{a = 0}^{r-1} \sqrt{p_a}\, \vert x_a\rangle \otimes \vert y_a \rangle,

ここで p0,,pr1p_0,\ldots,p_{r-1} は和が 11 になる正の実数で、{x0,,xr1}\{\vert x_0\rangle,\ldots,\vert x_{r-1}\rangle\}{y0,,yr1}\{\vert y_0\rangle,\ldots,\vert y_{r-1}\rangle\}両方が正規直交集合です。

ψ\vert\psi\rangle のシュミット分解における値

p0,,pr1\sqrt{p_0},\ldots,\sqrt{p_{r-1}}

シュミット係数と呼ばれ、一意に決まります(順序を除いて)— これらだけが ψ\vert\psi\rangle のそのような表現に現れることのできる正の実数です。 一方、集合

{x0,,xr1}{y0,,yr1}\{\vert x_0\rangle,\ldots,\vert x_{r-1}\rangle\} \quad\text{と}\quad \{\vert y_0\rangle,\ldots,\vert y_{r-1}\rangle\}

は一意には決まりません。これらのベクトル集合の選択における自由度は、以下の説明で明らかになります。

与えられた量子状態ベクトル ψ\vert\psi\rangle が確かにシュミット分解を持つことを検証し、その過程でシュミット分解の求め方を学びましょう。

まず、系 X\mathsf{X} に対応するベクトル空間の任意の(必ずしも直交しない)基底 {x0,,xn1}\{\vert x_0\rangle, \ldots, \vert x_{n-1}\rangle\} を考えます。 これが基底であるため、以下の等式が成り立つような z0,,zn1\vert z_0\rangle,\ldots,\vert z_{n-1}\rangle の選択が常に一意に存在します。

ψ=a=0n1xaza(1)\vert \psi\rangle = \sum_{a = 0}^{n-1} \vert x_a\rangle \otimes \vert z_a \rangle \tag{1}

例えば、{x0,,xn1}\{\vert x_0\rangle,\ldots,\vert x_{n-1}\rangle\}X\mathsf{X} に関連する標準基底とします。 X\mathsf{X} の古典状態集合が {0,,n1}\{0,\ldots,n-1\} であると仮定すると、各 a{0,,n1}a\in\{0,\ldots,n-1\} に対して xa=a\vert x_a\rangle = \vert a\rangle となり、

ψ=a=0n1aza\vert\psi\rangle = \sum_{a = 0}^{n-1} \vert a\rangle \otimes \vert z_a\rangle

が成り立ちます。ここで

za=(aIY)ψ\vert z_a \rangle = ( \langle a \vert \otimes \mathbb{I}_{\mathsf{Y}}) \vert \psi\rangle

です(各 a{0,,n1}a\in\{0,\ldots,n-1\} に対して)。 このような表現は、X\mathsf{X} の標準基底測定を考えるときに頻繁に用いられます。

重要なのは、この例での z0,,zn1\vert z_0\rangle,\ldots,\vert z_{n-1}\rangle の公式

za=(aIY)ψ\vert z_a \rangle = ( \langle a \vert \otimes \mathbb{I}_{\mathsf{Y}}) \vert \psi\rangle

が機能するのは {0,,n1}\{\vert 0\rangle,\ldots,\vert n-1\rangle\}正規直交基底であるからです。 一般に、{x0,,xn1}\{\vert x_0\rangle,\ldots,\vert x_{n-1}\rangle\} が必ずしも正規直交でない基底であっても、z0,,zn1\vert z_0\rangle,\ldots,\vert z_{n-1}\rangle は等式 (1)(1) によって一意に決まりますが、異なる公式が必要になります。 それらを求める一つの方法は、まず

waxb={1a=b0ab\langle w_a \vert x_b \rangle = \begin{cases} 1 & a=b\\ 0 & a\neq b \end{cases}

が全ての a,b{0,,n1}a,b\in\{0,\ldots,n-1\} に対して満たされるようなベクトル w0,,wn1\vert w_0\rangle,\ldots,\vert w_{n-1}\rangle を特定することで、その後

za=(waIY)ψ\vert z_a \rangle = (\langle w_a \vert \otimes \mathbb{I}_{\mathsf{Y}}) \vert \psi\rangle

が得られます。

X\mathsf{X} に対応するベクトル空間の与えられた基底 {x0,,xn1}\{\vert x_0\rangle,\ldots,\vert x_{n-1}\rangle\} に対して、等式 (1)(1) を満たす一意に決まるベクトル z0,,zn1\vert z_0\rangle,\ldots,\vert z_{n-1}\rangle は、{x0,,xn1}\{\vert x_0\rangle,\ldots,\vert x_{n-1}\rangle\} が正規直交基底であっても特別な性質を持つとは限りません。 しかし、{x0,,xn1}\{\vert x_0\rangle, \ldots, \vert x_{n-1}\rangle\} を縮約状態

ρ=TrY(ψψ)\rho = \operatorname{Tr}_{\mathsf{Y}} \bigl( \vert \psi\rangle \langle \psi \vert \bigr)

の固有ベクトルの正規直交基底として選ぶと、興味深いことが起こります。 具体的には、等式 (1)(1) が成り立つように一意に決まる集合 {z0,,zn1}\{\vert z_0\rangle,\ldots,\vert z_{n-1}\rangle\} は、必ず直交しなければなりません。

詳しく見るために、ρ\rho のスペクトル分解を考えます。

ρ=a=0n1paxaxa\rho = \sum_{a = 0}^{n-1} p_a \vert x_a \rangle \langle x_a \vert

ここで ρ\rho は密度行列なので固有値ベクトル (p0,,pn1)(p_0,\ldots,p_{n-1}) が確率ベクトルを形成することを認識して、ρ\rho の固有値を p0,,pn1p_0,\ldots,p_{n-1} と表記しており、{x0,,xn1}\{\vert x_0\rangle,\ldots,\vert x_{n-1}\rangle\} はこれらの固有値に対応する固有ベクトルの正規直交基底です。 等式 (1)(1) が成り立つ一意な集合 {z0,,zn1}\{\vert z_0\rangle,\ldots,\vert z_{n-1}\rangle\} が必ず直交することを示すために、部分トレースの計算から始めます。

TrY(ψψ)=a,b=0n1xaxbTr(zazb)=a,b=0n1zbzaxaxb.\begin{aligned} \operatorname{Tr}_{\mathsf{Y}} (\vert\psi\rangle\langle\psi\vert) & = \sum_{a,b = 0}^{n-1} \vert x_a\rangle\langle x_b\vert \operatorname{Tr}(\vert z_a\rangle\langle z_b\vert)\\ & = \sum_{a,b = 0}^{n-1} \langle z_b\vert z_a\rangle \, \vert x_a\rangle\langle x_b\vert. \end{aligned}

この表現は ρ\rho のスペクトル分解と一致しなければなりません。 {x0,,xn1}\{\vert x_0\rangle,\ldots,\vert x_{n-1}\rangle\} が基底であるため、行列の集合

{xaxb:a,b{0,,n1}}\bigl\{ \vert x_a\rangle\langle x_b\vert \,:\, a,b\in\{0,\ldots,n-1\} \bigr\}

は線形独立であると結論付けられ、したがって

zbza={paa=b0ab,\langle z_b \vert z_a\rangle = \begin{cases} p_a & a=b\\[1mm] 0 & a\neq b, \end{cases}

が成り立ち、{z0,,zn1}\{\vert z_0\rangle,\ldots,\vert z_{n-1}\rangle\} が直交していることが確立されます。

ψ\vert\psi\rangle のシュミット分解の取得はほぼ完了しています。 (1)(1) において pa=0p_a = 0 となる項を捨て、残りの各項について単位ベクトル ya\vert y_a\rangle を用いて za=paya\vert z_a\rangle = \sqrt{p_a}\vert y_a\rangle と書くだけです。

これを行うための便利な方法は、縮約状態 ρ\rho のスペクトル分解において固有値・固有ベクトル対を自由に番号付けできるという観察から始まります — したがって固有値が降順に並んでいると仮定できます。

p0p1pn1.p_0 \geq p_1 \geq \cdots \geq p_{n-1}.

r=rank(ρ)r = \operatorname{rank}(\rho) とすると、p0,,pr1>0p_0,\ldots,p_{r-1} > 0 かつ pr==pn1=0p_r = \cdots = p_{n-1} = 0 となります。 よって、

ρ=a=0r1paxaxa\rho = \sum_{a = 0}^{r-1} p_a \vert x_a \rangle \langle x_a \vert

が得られ、量子状態ベクトル ψ\vert \psi \rangle を次のように書けます。

ψ=a=0r1xaza.\vert\psi\rangle = \sum_{a = 0}^{r-1} \vert x_a\rangle \otimes \vert z_a\rangle.

a=0,,r1a=0,\ldots,r-1 に対して

za2=zaza=pa>0\| \vert z_a \rangle \|^2 = \langle z_a \vert z_a \rangle = p_a > 0

であることから、単位ベクトル y0,,yr1\vert y_0 \rangle,\ldots,\vert y_{r-1}\rangle を次のように定義できます。

ya=zaza=zapa,\vert y_a\rangle = \frac{\vert z_a\rangle}{\|\vert z_a\rangle\|} = \frac{\vert z_a\rangle}{\sqrt{p_a}},

これにより各 a{0,,r1}a\in\{0,\ldots,r-1\} に対して za=paya\vert z_a\rangle = \sqrt{p_a}\vert y_a\rangle となります。 {z0,,zr1}\{\vert z_0\rangle, \ldots, \vert z_{r-1}\rangle\} が直交かつ非ゼロであるため、{y0,,yr1}\{\vert y_0\rangle, \ldots, \vert y_{r-1}\rangle\}正規直交集合となり、ψ\vert\psi\rangle のシュミット分解が得られました。

ψ=a=0r1paxaya\vert \psi\rangle = \sum_{a = 0}^{r-1} \sqrt{p_a}\, \vert x_a\rangle \otimes \vert y_a \rangle

ベクトル {x0,,xr1}\{\vert x_0\rangle,\ldots,\vert x_{r-1}\rangle\}{y0,,yr1}\{\vert y_0\rangle,\ldots,\vert y_{r-1}\rangle\} の選択について言えば、{x0,,xr1}\{\vert x_0\rangle,\ldots,\vert x_{r-1}\rangle\} を縮約状態 TrY(ψψ)\operatorname{Tr}_{\mathsf{Y}}(\vert\psi\rangle\langle\psi\vert) の非ゼロ固有値に対応する固有ベクトルの任意の正規直交集合として選べます(上記のように)。その場合、{y0,,yr1}\{\vert y_0\rangle,\ldots,\vert y_{r-1}\rangle\} は一意に決まります。

この状況は2つの系の間で対称なので、代わりに {y0,,yr1}\{\vert y_0\rangle,\ldots,\vert y_{r-1}\rangle\} を縮約状態 TrX(ψψ)\operatorname{Tr}_{\mathsf{X}}(\vert\psi\rangle\langle\psi\vert) の非ゼロ固有値に対応する固有ベクトルの任意の正規直交集合として選ぶこともでき、その場合 {x0,,xr1}\{\vert x_0\rangle,\ldots,\vert x_{r-1}\rangle\} が一意に決まります。

ただし、一方の集合を対応する縮約状態の固有ベクトルとして選ぶと、もう一方が決まることに注意してください — つまり、両者を独立に選ぶことはできません。

本シリーズではこれ以上登場しませんが、注目すべき点として、縮約状態 TrX(ψψ)\operatorname{Tr}_{\mathsf{X}}(\vert\psi\rangle\langle\psi\vert) の非ゼロ固有値 p0,,pr1p_0,\ldots,p_{r-1} は、(X,Y)(\mathsf{X},\mathsf{Y}) の任意の純粋状態 ψ\vert\psi\rangle に対して、縮約状態 TrY(ψψ)\operatorname{Tr}_{\mathsf{Y}}(\vert\psi\rangle\langle\psi\vert) の非ゼロ固有値と必ず一致しなければなりません。

直観的には、対 (X,Y)(\mathsf{X},\mathsf{Y}) が純粋状態にあるとき、X\mathsf{X}Y\mathsf{Y} の縮約状態にはまったく同じ量のランダム性が含まれています。 この事実はシュミット分解によって明らかになります。どちらの場合も縮約状態の固有値は純粋状態のシュミット係数の二乗と一致しなければなりません。

純粋化のユニタリ同値性

シュミット分解を使って、純粋化のユニタリ同値性として知られる純粋化に関する根本的に重要な事実を確立できます。

定理

純粋化のユニタリ同値性:X\mathsf{X}Y\mathsf{Y} を系とし、ψ\vert\psi\rangleϕ\vert\phi\rangle(X,Y)(\mathsf{X},\mathsf{Y}) の量子状態ベクトルで、ともに X\mathsf{X} の同じ状態を純粋化するものとします。記号で表すと、

TrY(ψψ)=ρ=TrY(ϕϕ)\operatorname{Tr}_{\mathsf{Y}} (\vert\psi\rangle\langle\psi\vert) = \rho = \operatorname{Tr}_{\mathsf{Y}} (\vert\phi\rangle\langle\phi\vert)

がある密度行列 ρ\rhoX\mathsf{X} の状態を表す)に対して成り立ちます。 このとき、Y\mathsf{Y} のみに作用するユニタリ操作 UU が存在し、最初の純粋化を2番目に変換します。

(IXU)ψ=ϕ.(\mathbb{I}_{\mathsf{X}} \otimes U) \vert\psi\rangle = \vert\phi\rangle.

この定理のいくつかの含意についてはレッスンを通じて議論しますが、まずこれがシュミット分解に関する先の議論からどのように導かれるかを見てみましょう。

仮定として、ψ\vert\psi\rangleϕ\vert\phi\rangle は対の系 (X,Y)(\mathsf{X},\mathsf{Y}) の量子状態ベクトルで、等式

TrY(ψψ)=ρ=TrY(ϕϕ)\operatorname{Tr}_{\mathsf{Y}} (\vert\psi\rangle\langle\psi\vert) = \rho = \operatorname{Tr}_{\mathsf{Y}} (\vert\phi\rangle\langle\phi\vert)

X\mathsf{X} の状態を表すある密度行列 ρ\rho に対して満たすとします。

ρ\rho のスペクトル分解を考えます。

ρ=a=0n1paxaxa\rho = \sum_{a = 0}^{n-1} p_a \vert x_a\rangle\langle x_a\vert

ここで {x0,,xn1}\{\vert x_0\rangle,\ldots,\vert x_{n-1}\rangle\}ρ\rho の固有ベクトルの正規直交基底です。 先述の手順に従って、ψ\vert\psi\rangleϕ\vert\phi\rangle の両方について以下の形のシュミット分解を得ることができます。

ψ=a=0r1paxauaϕ=a=0r1paxava\begin{aligned} \vert\psi\rangle & = \sum_{a = 0}^{r-1} \sqrt{p_a} \, \vert x_a\rangle \otimes \vert u_a\rangle\\[1mm] \vert\phi\rangle & = \sum_{a = 0}^{r-1} \sqrt{p_a} \, \vert x_a\rangle \otimes \vert v_a\rangle \end{aligned}

これらの式において rrρ\rho のランクで、{u0,,ur1}\{\vert u_0\rangle,\ldots,\vert u_{r-1}\rangle\}{v0,,vr1}\{\vert v_0\rangle,\ldots,\vert v_{r-1}\rangle\}Y\mathsf{Y} に対応する空間における正規直交ベクトルの集合です。

同じ空間内で同じ数の要素を持つ2つの正規直交集合に対して、常に一方を他方に変換するユニタリ行列が存在するため、a=0,,r1a = 0,\ldots,r-1 に対して Uua=vaU \vert u_a\rangle = \vert v_a\rangle となるようなユニタリ行列 UU を選べます。 特に、そのような行列 UU を見つけるには、まずグラム-シュミット直交化法を使って正規直交集合を正規直交基底 {u0,,um1}\{\vert u_0\rangle,\ldots,\vert u_{m-1}\rangle\}{v0,,vm1}\{\vert v_0\rangle,\ldots,\vert v_{m-1}\rangle\}mmY\mathsf{Y} に対応する空間の次元)に拡張し、その後

U=a=0m1vauaU = \sum_{a = 0}^{m-1} \vert v_a\rangle\langle u_a\vert

とすることができます。

すると次のことがわかります。

(IXU)ψ=a=0r1paxaUua=a=0r1paxava=ϕ,\begin{aligned} (\mathbb{I}_{\mathsf{X}} \otimes U) \vert\psi\rangle & = \sum_{a = 0}^{r-1} \sqrt{p_a} \, \vert x_a\rangle \otimes U \vert u_a\rangle\\ & = \sum_{a = 0}^{r-1} \sqrt{p_a} \, \vert x_a\rangle \otimes \vert v_a\rangle\\ & = \vert\phi\rangle, \end{aligned}

これで証明が完了しました。

以下では、純粋化のユニタリ同値性に関連する多くの興味深い例と含意のうちのいくつかを示します。 後ほど、フィデリティの文脈でウールマンの定理として知られる、さらに重要な例が登場します。

超高密度符号化

超高密度符号化プロトコルでは、アリスとボブはeビットを共有します。すなわち、アリスはQubit A\mathsf{A} を、ボブはQubit B\mathsf{B} を保持しており、ペア (A,B)(\mathsf{A},\mathsf{B}) は全体として ϕ+\vert\phi^{+}\rangle ベル状態にあります。 このプロトコルは、アリスが自分のQubit A\mathsf{A} にユニタリ操作を施すことで、この共有状態を4つのベル状態 ϕ+,\vert\phi^+\rangle, ϕ,\vert\phi^-\rangle, ψ+,\vert\psi^+\rangle, ψ\vert\psi^-\rangle のいずれかに変換できることを説明するものです。 アリスがこれを行った後、A\mathsf{A} をボブに送り、ボブはペア (A,B)(\mathsf{A},\mathsf{B}) に対する測定を行い、どのベル状態にあるかを確認します。

4つのベル状態すべてに対して、ボブのQubit B\mathsf{B} の縮約状態は完全混合状態です。

TrA(ϕ+ϕ+)=TrA(ϕϕ)=TrA(ψ+ψ+)=TrA(ψψ)=I2\operatorname{Tr}_{\mathsf{A}}(\vert\phi^+\rangle\langle\phi^+\vert) = \operatorname{Tr}_{\mathsf{A}}(\vert\phi^-\rangle\langle\phi^-\vert) = \operatorname{Tr}_{\mathsf{A}}(\vert\psi^+\rangle\langle\psi^+\vert) = \operatorname{Tr}_{\mathsf{A}}(\vert\psi^-\rangle\langle\psi^-\vert) = \frac{\mathbb{I}}{2}

浄化のユニタリ等価性から、各ベル状態に対して、アリスのQubit A\mathsf{A} のみに作用するユニタリ操作が存在し、ϕ+\vert\phi^+\rangle を目的のベル状態に変換できることが直ちに結論づけられます。 これはプロトコルの詳細を明示するわけではありませんが、浄化のユニタリ等価性によって超高密度符号化が可能であることがすぐに示唆されます。

また、ベル状態をより大きな系における完全混合状態の浄化の正規直交基底に置き換えれば、超高密度符号化をより大きな系へ一般化することは常に可能であるとも結論づけられます。

暗号への応用

浄化のユニタリ等価性は、量子情報を用いた暗号プリミティブの実装に関して示唆を与えます。 例えば、浄化のユニタリ等価性によって、量子情報だけを用いて理想的な形式のビットコミットメントを実装することは不可能であることが明らかになります。

ビットコミットメントプリミティブには、(互いに信頼していない)アリスとボブの2人の参加者が関与し、2つのフェーズがあります。

  • 最初のフェーズはコミットフェーズであり、アリスは二値 b{0,1}b\in\{0,1\} にコミットします。 このコミットメントは拘束力を持つ必要があります。つまり、アリスは考えを変えられず、また秘匿性も必要です。つまり、ボブはアリスがどの値にコミットしたかを知ることができません。
  • 2番目のフェーズは開示フェーズであり、アリスがコミットしたビットがボブに明かされ、ボブはそれが本当にコミットされた値であることを確信できるようになります。

直感的・操作的に言えば、ビットコミットメントの最初のフェーズは、アリスが二値を紙に書き、その紙を金庫に入れてボブに渡し、自分は鍵を手元に置くようなものとして機能するべきです。 アリスは紙に書いた二値にコミットしています。なぜなら、金庫はボブの手元にあるので(拘束力)、ボブは金庫を開けられないためアリスがどの値にコミットしたかわかりません(秘匿性)。 第2フェーズは、アリスがボブに鍵を渡し、ボブが金庫を開けてアリスがコミットした値を確認できるようなものとして機能するべきです。

結局のところ、量子情報のみによって完全なビットコミットメントプロトコルを実装することは不可能です。これは浄化のユニタリ等価性に矛盾するためです。 以下はそれを示す議論の概要です。

まず、プロトコルの実行中にアリスとボブがユニタリ操作のみを行うか、初期化された新しい系を導入するのみと仮定できます。 すべてのチャネルがStinespring表現を持つという事実により、この仮定が可能になります。

プロトコルのコミットフェーズの終了時に、ボブは何らかの複合系を保持しており、その状態はアリスが 00 にコミットした場合は ρ0\rho_0、アリスが 11 にコミットした場合は ρ1\rho_1 の2つの量子状態のいずれかでなければなりません。 プロトコルが完全に秘匿性を持つためには、ボブがこれら2つの状態を区別できないことが必要なので、ρ0=ρ1\rho_0 = \rho_1 でなければなりません。 (そうでなければ、これらの状態を確率的に識別できる測定が存在することになります。)

しかし、アリスとボブがユニタリ操作のみを使用したことから、コミットフェーズ後のプロトコルに関与するすべての系の状態は純粋状態でなければなりません。 具体的に、アリスが 00 にコミットするときのプロトコルに関与するすべての系の純粋状態を ψ0\vert\psi_0\rangle、アリスが 11 にコミットするときを ψ1\vert\psi_1\rangle とします。 アリスとボブの(複合系の可能性がある)系をそれぞれ A\mathsf{A}B\mathsf{B} と表すと、

ρ0=TrA(ψ0ψ0)ρ1=TrA(ψ1ψ1).\begin{aligned} \rho_0 & = \operatorname{Tr}_{\mathsf{A}}(\vert\psi_0\rangle\langle\psi_0\vert)\\[1mm] \rho_1 & = \operatorname{Tr}_{\mathsf{A}}(\vert\psi_1\rangle\langle\psi_1\vert). \end{aligned}

完全に秘匿性を持つプロトコルに必要な ρ0=ρ1\rho_0 = \rho_1 という条件から、ψ0\vert\psi_0\rangleψ1\vert\psi_1\rangle は同じ状態の浄化であることがわかります。したがって、浄化のユニタリ等価性から、A\mathsf{A} のみに作用するユニタリ操作 UU が存在して

(UIB)ψ0=ψ1.(U\otimes\mathbb{I}_{\mathsf{B}})\vert\psi_0\rangle = \vert\psi_1\rangle.

が成り立ちます。したがって、アリスは A\mathsf{A}UU を適用することで 00 から 11 へのコミットメントを自由に変更でき、または UU^{\dagger} を適用することで 11 から 00 へ変更できます。つまり、この仮想的なプロトコルは拘束力を全く持たないことになります。

Hughston-Jozsa-Woottersの定理

本レッスンのこのパートで、浄化のユニタリ等価性の最後の応用として、Hughston-Jozsa-Wootters定理として知られる以下の定理について説明します。 (これは、実際にはその名で知られる定理のわずかに簡略化された表現です。)

定理

Hughston-Jozsa-Wootters: X\mathsf{X}Y\mathsf{Y} を系とし、ϕ\vert\phi\rangle をペア (X,Y)(\mathsf{X},\mathsf{Y}) の量子状態ベクトルとします。 また、NN を任意の正の整数、(p0,,pN1)(p_0,\ldots,p_{N-1}) を確率ベクトル、ψ0,,ψN1\vert\psi_0\rangle,\ldots,\vert\psi_{N-1}\rangleX\mathsf{X} の状態を表す量子状態ベクトルとし、

TrY(ϕϕ)=a=0N1paψaψa.\operatorname{Tr}_{\mathsf{Y}}\bigl(\vert\phi\rangle\langle\phi\vert\bigr) = \sum_{a = 0}^{N-1} p_a \vert\psi_a\rangle\langle\psi_a\vert.

が成り立つとします。 Y\mathsf{Y} 上の(一般の)測定 {P0,,PN1}\{P_0,\ldots,P_{N-1}\} が存在して、(X,Y)(\mathsf{X},\mathsf{Y}) が状態 ϕ\vert\phi\rangle にあるときに Y\mathsf{Y} に対してこの測定を行うと、以下の2つの条件が満たされます:

  1. 各測定結果 a{0,,N1}a\in\{0,\ldots,N-1\} は確率 pap_a で現れます。
  2. 測定結果 aa が得られたとき、X\mathsf{X} の状態は ψa\vert\psi_a\rangle になります。

直感的に言えば、この定理は、2つの系の純粋状態がある場合、最初の系の縮約状態を純粋状態の凸結合として考えるいかなる方法に対しても、2番目の系の測定を行うことで、最初の系についてのその考え方を現実のものにできる、ということを示しています。 NN は必ずしも X\mathsf{X}Y\mathsf{Y} の古典状態数で制限されないことに注意してください。 例えば、X\mathsf{X}Y\mathsf{Y} がQubitであっても N=1,000,000N = 1,000,000 であることもあり得ます。

この定理を浄化のユニタリ等価性を用いて証明します。まず、古典状態集合が {0,,N1}\{0,\ldots,N-1\} である新しい系 Z\mathsf{Z} を導入します。 トリプル (X,Y,Z)(\mathsf{X},\mathsf{Y},\mathsf{Z}) の次の2つの量子状態ベクトルを考えます。

γ0=ϕXY0Zγ1=a=0N1paψaX0YaZ\begin{aligned} \vert\gamma_0\rangle & = \vert\phi\rangle_{\mathsf{XY}}\otimes\vert 0\rangle_{\mathsf{Z}}\\[1mm] \vert\gamma_1\rangle & = \sum_{a = 0}^{N-1} \sqrt{p_a}\, \vert\psi_a\rangle_{\mathsf{X}} \otimes \vert 0\rangle_{\mathsf{Y}} \otimes \vert a\rangle_{\mathsf{Z}} \end{aligned}

最初のベクトル γ0\vert\gamma_0\rangle は、与えられた量子状態ベクトル ϕ\vert\phi\rangle に新しい系 Z\mathsf{Z}0\vert 0\rangle をテンソル積したものです。 2番目のベクトル γ1\vert\gamma_1\rangle は、少なくとも Y\mathsf{Y}Z\mathsf{Z} に置き換えれば定理を自明にするような量子状態ベクトルです。なぜなら、Z\mathsf{Z} に対する標準基底測定を行うと、各結果 aa が確率 pap_a で得られ、その結果が得られたとき X\mathsf{X} の状態は ψa\vert\psi_a\rangle になることが明らかだからです。

ペア (Y,Z)(\mathsf{Y},\mathsf{Z}) をひとつの複合系として考えてトレースアウトすると X\mathsf{X} が残り、状態

ρ=a=0N1paψaψa.\rho = \sum_{a = 0}^{N-1} p_a \vert\psi_a\rangle\langle\psi_a\vert.

の2つの異なる浄化を特定したことになります。

具体的に、最初の浄化については

TrYZ(γ0γ0)=TrY(ϕϕ)=ρ\operatorname{Tr}_{\mathsf{YZ}} (\vert\gamma_0\rangle\langle\gamma_0\vert) = \operatorname{Tr}_{\mathsf{Y}} (\vert\phi\rangle\langle\phi\vert) = \rho

であり、2番目の浄化については

TrYZ(γ1γ1)=a,b=0N1papbψaψaTr(00ab)=a=0N1paψaψa=ρ.\begin{aligned} \operatorname{Tr}_{\mathsf{YZ}} (\vert\gamma_1\rangle\langle\gamma_1\vert) & = \sum_{a,b = 0}^{N-1} \sqrt{p_a}\sqrt{p_b} \, \vert\psi_a\rangle\langle\psi_a\vert \operatorname{Tr}(\vert 0\rangle\langle 0\vert \otimes \vert a\rangle\langle b\vert)\\ & = \sum_{a = 0}^{N-1} p_a \, \vert\psi_a\rangle\langle\psi_a\vert\\ & = \rho. \end{aligned}

です。

したがって、浄化のユニタリ等価性により、(Y,Z)(\mathsf{Y},\mathsf{Z}) 上のユニタリ操作 UU が存在して

(IXU)γ0=γ1(\mathbb{I}_{\mathsf{X}} \otimes U) \vert \gamma_0 \rangle = \vert\gamma_1\rangle

を満たします。

このユニタリ操作 UU を用いて、定理の要件を満たす測定を次の図のように実装できます。 言葉で説明すると、0\vert 0\rangle 状態に初期化した新しい系 Z\mathsf{Z} を導入し、(Y,Z)(\mathsf{Y},\mathsf{Z})UU を適用することで (X,Y,Z)(\mathsf{X},\mathsf{Y},\mathsf{Z}) の状態を γ0\vert\gamma_0\rangle から γ1\vert\gamma_1\rangle に変換し、次に Z\mathsf{Z} を標準基底測定で測定します。これにより、すでに確認した通り望ましい振る舞いが得られます。

HJW定理の測定の量子Circuit実装

図の点線の矩形はこの測定の実装を表しており、次のように正半定値行列の集合 {P0,,PN1}\{P_0,\ldots,P_{N-1}\} として記述できます。

Pa=(IY0)U(IYaa)U(IY0)P_a = (\mathbb{I}_{\mathsf{Y}} \otimes \langle 0\vert) U^{\dagger} (\mathbb{I}_{\mathsf{Y}} \otimes \vert a\rangle\langle a \vert)U (\mathbb{I}_{\mathsf{Y}} \otimes \vert 0\rangle)