Zum Hauptinhalt springen

Purifikationen

Definition von Purifikationen

Beginnen wir mit einer präzisen mathematischen Definition für Purifikationen.

Definition

Sei X\mathsf{X} ein System in einem Zustand, der durch eine Dichtematrix ρ\rho dargestellt wird, und sei ψ\vert\psi\rangle ein Quantenzustandsvektor eines Paares (X,Y),(\mathsf{X},\mathsf{Y}), der ρ\rho ergibt, wenn über Y\mathsf{Y} die Spur gezogen wird:

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

Der Zustandsvektor ψ\vert\psi\rangle wird dann als eine Purifikation von ρ\rho bezeichnet.

Der reine Zustand ψψ,\vert\psi\rangle\langle\psi\vert, als Dichtematrix statt als Quantenzustandsvektor ausgedrückt, wird ebenfalls häufig als Purifikation von ρ\rho bezeichnet, wenn die Gleichung in der Definition gilt – wir verwenden den Begriff jedoch im Allgemeinen für einen Quantenzustandsvektor.

Der Begriff Purifikation wird auch allgemeiner verwendet, wenn die Reihenfolge der Systeme umgekehrt ist, wenn die Systeme und Zustände andere Namen tragen (natürlich), und wenn mehr als zwei Systeme vorhanden sind. Wenn zum Beispiel ψ\vert \psi \rangle ein Quantenzustandsvektor ist, der einen reinen Zustand eines zusammengesetzten Systems (A,B,C)(\mathsf{A},\mathsf{B},\mathsf{C}) darstellt, und die Gleichung

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

für eine Dichtematrix ρ\rho gilt, die einen Zustand des Systems (A,C)(\mathsf{A},\mathsf{C}) darstellt, dann wird ψ\vert\psi\rangle weiterhin als Purifikation von ρ\rho bezeichnet.

Im Rahmen dieser Lektion konzentrieren wir uns jedoch auf die spezifische Form, die in der Definition beschrieben wird. Eigenschaften und Fakten über Purifikationen gemäß dieser Definition lassen sich typischerweise auf mehr als zwei Systeme verallgemeinern, indem man die Systeme neu anordnet und in zwei zusammengesetzte Systeme aufteilt, von denen eines die Rolle von X\mathsf{X} und das andere die Rolle von Y\mathsf{Y} übernimmt.

Existenz von Purifikationen

Seien X\mathsf{X} und Y\mathsf{Y} zwei beliebige Systeme und ρ\rho ein gegebener Zustand von X.\mathsf{X}. Wir werden beweisen, dass ein Quantenzustandsvektor ψ\vert\psi\rangle von (X,Y)(\mathsf{X},\mathsf{Y}) existiert, der ρ\rho purifiziert — was eine andere Weise ist zu sagen, dass ψ\vert\psi\rangle eine Purifikation von ρ\rho ist — vorausgesetzt, das System Y\mathsf{Y} ist groß genug. Insbesondere gilt: Wenn Y\mathsf{Y} mindestens so viele klassische Zustände hat wie X,\mathsf{X}, existiert für jeden Zustand ρ\rho notwendigerweise eine Purifikation dieser Form. Für manche Zustände ρ\rho werden weniger klassische Zustände von Y\mathsf{Y} benötigt; im Allgemeinen sind rank(ρ)\operatorname{rank}(\rho) klassische Zustände von Y\mathsf{Y} notwendig und hinreichend für die Existenz eines Quantenzustandsvektors von (X,Y),(\mathsf{X},\mathsf{Y}), der ρ\rho purifiziert.

Betrachte zunächst eine beliebige Darstellung von ρ\rho als Konvexkombination von nn reinen Zuständen, für eine beliebige positive ganze Zahl n.n.

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

In diesem Ausdruck ist (p0,,pn1)(p_0,\ldots,p_{n-1}) ein Wahrscheinlichkeitsvektor und ϕ0,,ϕn1\vert\phi_0\rangle,\ldots,\vert\phi_{n-1}\rangle sind Quantenzustandsvektoren von X.\mathsf{X}.

Eine Möglichkeit, einen solchen Ausdruck zu erhalten, ist der Spektralsatz, bei dem nn die Anzahl der klassischen Zustände von X\mathsf{X} ist, p0,,pn1p_0,\ldots,p_{n-1} die Eigenwerte von ρ\rho sind und ϕ0,,ϕn1\vert\phi_0\rangle,\ldots,\vert\phi_{n-1}\rangle orthonormale Eigenvektoren zu diesen Eigenwerten sind.

Es ist tatsächlich nicht nötig, die Terme, die den Null-Eigenwerten von ρ\rho entsprechen, in die Summe aufzunehmen, was es uns erlaubt, alternativ n=rank(ρ)n = \operatorname{rank}(\rho) zu wählen und p0,,pn1p_0,\ldots,p_{n-1} als die von null verschiedenen Eigenwerte von ρ\rho zu nehmen. Dies ist der minimale Wert von n,n, für den ein Ausdruck der obigen Form existiert.

Um es klar zu sagen: Es ist nicht notwendig, dass die gewählte Darstellung von ρ\rho als Konvexkombination reiner Zustände aus dem Spektralsatz stammt — dies ist nur eine Möglichkeit, einen solchen Ausdruck zu erhalten. Insbesondere kann nn eine beliebige positive ganze Zahl sein, die Einheitsvektoren ϕ0,,ϕn1\vert\phi_0\rangle,\ldots,\vert\phi_{n-1}\rangle müssen nicht orthogonal sein, und die Wahrscheinlichkeiten p0,,pn1p_0,\ldots,p_{n-1} müssen keine Eigenwerte von ρ\rho sein.

Wir können nun eine Purifikation von ρ\rho wie folgt identifizieren.

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

Dabei nehmen wir an, dass die klassischen Zustände von Y\mathsf{Y} die Werte 0,,n10,\ldots,n-1 umfassen. Falls nicht, kann eine beliebige Auswahl von nn verschiedenen klassischen Zuständen von Y\mathsf{Y} für 0,,n10,\ldots,n-1 eingesetzt werden. Zu überprüfen, dass es sich tatsächlich um eine Purifikation von ρ\rho handelt, ist eine einfache Angelegenheit der Berechnung der partiellen Spur, was auf die folgenden zwei äquivalenten Weisen geschehen kann.

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

Allgemeiner gilt: Für jede orthonormale Menge von Vektoren {γ0,,γn1}\{\vert\gamma_0\rangle,\ldots,\vert\gamma_{n-1}\rangle\} ist der Quantenzustandsvektor

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

eine Purifikation von ρ.\rho.

Beispiel

Seien X\mathsf{X} und Y\mathsf{Y} beide Qubits und

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

eine Dichtematrix, die einen Zustand von X\mathsf{X} darstellt.

Mit dem Spektralsatz können wir ρ\rho schreiben als

ρ=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,

wobei ψθ=cos(θ)0+sin(θ)1.\vert \psi_{\theta} \rangle = \cos(\theta) \vert 0\rangle + \sin(\theta)\vert 1\rangle. Der Quantenzustandsvektor

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

der einen reinen Zustand des Paares (X,Y)(\mathsf{X},\mathsf{Y}) beschreibt, ist daher eine Purifikation von ρ.\rho.

Alternativ können wir schreiben

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

Dies ist eine Konvexkombination reiner Zustände, aber keine Spektralzerlegung, da 0\vert 0\rangle und +\vert +\rangle nicht orthogonal sind und 1/21/2 kein Eigenwert von ρ\rho ist. Dennoch ist der Quantenzustandsvektor

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

eine Purifikation von ρ.\rho.

Schmidt-Zerlegungen

Als nächstes besprechen wir Schmidt-Zerlegungen, d. h. Ausdrücke von Quantenzustandsvektoren von Paaren von Systemen, die eine bestimmte Form annehmen. Schmidt-Zerlegungen stehen in engem Zusammenhang mit Purifikationen und sind auch für sich genommen sehr nützlich. Wenn man über einen gegebenen Quantenzustandsvektor ψ\vert\psi\rangle eines Paares von Systemen nachdenkt, besteht der erste Schritt häufig darin, eine Schmidt-Zerlegung dieses Zustands zu identifizieren oder zu betrachten.

Definition

Sei ψ\vert \psi\rangle ein gegebener Quantenzustandsvektor eines Paares von Systemen (X,Y).(\mathsf{X},\mathsf{Y}). Eine Schmidt-Zerlegung von ψ\vert\psi\rangle ist ein Ausdruck der Form

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

wobei p0,,pr1p_0,\ldots,p_{r-1} positive reelle Zahlen mit Summe 11 sind und beide Mengen {x0,,xr1}\{\vert x_0\rangle,\ldots,\vert x_{r-1}\rangle\} und {y0,,yr1}\{\vert y_0\rangle,\ldots,\vert y_{r-1}\rangle\} orthonormal sind.

Die Werte

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

in einer Schmidt-Zerlegung von ψ\vert\psi\rangle heißen Schmidt-Koeffizienten, die eindeutig bestimmt sind (bis auf ihre Reihenfolge) — es sind die einzigen positiven reellen Zahlen, die in einem solchen Ausdruck von ψ\vert\psi\rangle auftreten können. Die Mengen

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

hingegen sind nicht eindeutig bestimmt, und die Freiheit bei der Wahl dieser Vektormengen wird in der folgenden Erklärung verdeutlicht.

Wir werden nun zeigen, dass ein gegebener Quantenzustandsvektor ψ\vert\psi\rangle tatsächlich eine Schmidt-Zerlegung besitzt, und dabei lernen, wie man eine solche findet.

Betrachte zunächst eine beliebige (nicht notwendigerweise orthogonale) Basis {x0,,xn1}\{\vert x_0\rangle, \ldots, \vert x_{n-1}\rangle\} des Vektorraums, der dem System X\mathsf{X} entspricht. Da dies eine Basis ist, gibt es immer eine eindeutig bestimmte Auswahl von Vektoren z0,,zn1,\vert z_0\rangle,\ldots,\vert z_{n-1}\rangle, für die die folgende Gleichung gilt.

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

Angenommen zum Beispiel, {x0,,xn1}\{\vert x_0\rangle,\ldots,\vert x_{n-1}\rangle\} ist die Standardbasis von X.\mathsf{X}. Unter der Annahme, dass die klassische Zustandsmenge von X\mathsf{X} gleich {0,,n1}\{0,\ldots,n-1\} ist, bedeutet dies, dass xa=a\vert x_a\rangle = \vert a\rangle für jedes a{0,,n1},a\in\{0,\ldots,n-1\}, und wir erhalten

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

wenn

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

für jedes a{0,,n1}.a\in\{0,\ldots,n-1\}. Solche Ausdrücke betrachten wir häufig, wenn wir über eine Standardbasismessung von X\mathsf{X} nachdenken.

Es ist wichtig zu beachten, dass die Formel

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

für die Vektoren z0,,zn1\vert z_0\rangle,\ldots,\vert z_{n-1}\rangle in diesem Beispiel nur funktioniert, weil {0,,n1}\{\vert 0\rangle,\ldots,\vert n-1\rangle\} eine orthonormale Basis ist. Im Allgemeinen gilt: Wenn {x0,,xn1}\{\vert x_0\rangle,\ldots,\vert x_{n-1}\rangle\} eine nicht notwendigerweise orthonormale Basis ist, sind die Vektoren z0,,zn1\vert z_0\rangle,\ldots,\vert z_{n-1}\rangle durch die Gleichung (1)(1) zwar eindeutig bestimmt, aber es wird eine andere Formel benötigt. Eine Möglichkeit, sie zu finden, besteht darin, zunächst Vektoren w0,,wn1\vert w_0\rangle,\ldots,\vert w_{n-1}\rangle zu identifizieren, so dass die Gleichung

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

für alle a,b{0,,n1}a,b\in\{0,\ldots,n-1\} erfüllt ist, woraufhin wir haben

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

Für eine gegebene Basis {x0,,xn1}\{\vert x_0\rangle,\ldots,\vert x_{n-1}\rangle\} des Vektorraums von X\mathsf{X} erfüllen die eindeutig bestimmten Vektoren z0,,zn1,\vert z_0\rangle,\ldots,\vert z_{n-1}\rangle, für die die Gleichung (1)(1) gilt, nicht notwendigerweise besondere Eigenschaften, selbst wenn {x0,,xn1}\{\vert x_0\rangle,\ldots,\vert x_{n-1}\rangle\} eine orthonormale Basis ist. Wenn wir jedoch {x0,,xn1}\{\vert x_0\rangle, \ldots, \vert x_{n-1}\rangle\} als orthonormale Basis aus Eigenvektoren des reduzierten Zustands

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

wählen, passiert etwas Interessantes. Genauer gesagt: Für die eindeutig bestimmte Kollektion {z0,,zn1},\{\vert z_0\rangle,\ldots,\vert z_{n-1}\rangle\}, für die die Gleichung (1)(1) gilt, stellt sich heraus, dass diese Kollektion orthogonal sein muss.

Im Einzelnen betrachte eine Spektralzerlegung von ρ.\rho.

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

Hier bezeichnen wir die Eigenwerte von ρ\rho mit p0,,pn1p_0,\ldots,p_{n-1} in Anerkennung der Tatsache, dass ρ\rho eine Dichtematrix ist — der Vektor der Eigenwerte (p0,,pn1)(p_0,\ldots,p_{n-1}) bildet also einen Wahrscheinlichkeitsvektor — während {x0,,xn1}\{\vert x_0\rangle,\ldots,\vert x_{n-1}\rangle\} eine orthonormale Basis von Eigenvektoren zu diesen Eigenwerten ist. Um zu sehen, dass die eindeutige Kollektion {z0,,zn1},\{\vert z_0\rangle,\ldots,\vert z_{n-1}\rangle\}, für die die Gleichung (1)(1) gilt, notwendigerweise orthogonal ist, können wir mit der Berechnung der partiellen Spur beginnen.

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}

Dieser Ausdruck muss mit der Spektralzerlegung von ρ\rho übereinstimmen. Da {x0,,xn1}\{\vert x_0\rangle,\ldots,\vert x_{n-1}\rangle\} eine Basis ist, schließen wir, dass die Menge von Matrizen

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

linear unabhängig ist, und somit folgt

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

womit gezeigt ist, dass {z0,,zn1}\{\vert z_0\rangle,\ldots,\vert z_{n-1}\rangle\} orthogonal ist.

Wir haben fast eine Schmidt-Zerlegung von ψ\vert\psi\rangle erhalten. Es bleibt noch, die Terme in (1),(1), für die pa=0p_a = 0 gilt, zu verwerfen und dann za=paya\vert z_a\rangle = \sqrt{p_a}\vert y_a\rangle für einen Einheitsvektor ya\vert y_a\rangle für jeden der verbleibenden Terme zu schreiben.

Ein bequemer Weg, dies zu tun, beginnt mit der Beobachtung, dass wir die Eigenwert/Eigenvektor-Paare in einer Spektralzerlegung des reduzierten Zustands ρ\rho beliebig nummerieren können — wir dürfen also annehmen, dass die Eigenwerte in absteigender Reihenfolge sortiert sind:

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

Mit r=rank(ρ)r = \operatorname{rank}(\rho) finden wir p0,,pr1>0p_0,\ldots,p_{r-1} > 0 und pr==pn1=0.p_r = \cdots = p_{n-1} = 0. Also haben wir

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

und wir können den Quantenzustandsvektor ψ\vert \psi \rangle schreiben als

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

Da

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

für a=0,,r1,a=0,\ldots,r-1, können wir Einheitsvektoren y0,,yr1\vert y_0 \rangle,\ldots,\vert y_{r-1}\rangle definieren als

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

so dass za=paya\vert z_a\rangle = \sqrt{p_a}\vert y_a\rangle für jedes a{0,,r1}.a\in\{0,\ldots,r-1\}. Da die Vektoren {z0,,zr1}\{\vert z_0\rangle, \ldots, \vert z_{r-1}\rangle\} orthogonal und von null verschieden sind, folgt, dass {y0,,yr1}\{\vert y_0\rangle, \ldots, \vert y_{r-1}\rangle\} eine orthonormale Menge ist, und wir haben somit eine Schmidt-Zerlegung von ψ\vert\psi\rangle erhalten.

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

Bezüglich der Wahl der Vektoren {x0,,xr1}\{\vert x_0\rangle,\ldots,\vert x_{r-1}\rangle\} und {y0,,yr1}\{\vert y_0\rangle,\ldots,\vert y_{r-1}\rangle\} können wir {x0,,xr1}\{\vert x_0\rangle,\ldots,\vert x_{r-1}\rangle\} als eine beliebige orthonormale Menge von Eigenvektoren zu den von null verschiedenen Eigenwerten des reduzierten Zustands TrY(ψψ)\operatorname{Tr}_{\mathsf{Y}}(\vert\psi\rangle\langle\psi\vert) wählen (wie wir es oben getan haben), wobei die Vektoren {y0,,yr1}\{\vert y_0\rangle,\ldots,\vert y_{r-1}\rangle\} dann eindeutig bestimmt sind.

Die Situation ist symmetrisch zwischen den beiden Systemen, daher können wir alternativ {y0,,yr1}\{\vert y_0\rangle,\ldots,\vert y_{r-1}\rangle\} als eine beliebige orthonormale Menge von Eigenvektoren zu den von null verschiedenen Eigenwerten des reduzierten Zustands TrX(ψψ)\operatorname{Tr}_{\mathsf{X}}(\vert\psi\rangle\langle\psi\vert) wählen, wobei dann die Vektoren {x0,,xr1}\{\vert x_0\rangle,\ldots,\vert x_{r-1}\rangle\} eindeutig bestimmt sind.

Beachte jedoch: Sobald eine der Mengen als Eigenvektormenge des entsprechenden reduzierten Zustands gewählt wurde, ist die andere bestimmt — sie können also nicht unabhängig voneinander gewählt werden.

Obwohl es in dieser Reihe nicht wieder vorkommen wird, ist es bemerkenswert, dass die von null verschiedenen Eigenwerte p0,,pr1p_0,\ldots,p_{r-1} des reduzierten Zustands TrX(ψψ)\operatorname{Tr}_{\mathsf{X}}(\vert\psi\rangle\langle\psi\vert) stets mit den von null verschiedenen Eigenwerten des reduzierten Zustands TrY(ψψ)\operatorname{Tr}_{\mathsf{Y}}(\vert\psi\rangle\langle\psi\vert) für jeden reinen Zustand ψ\vert\psi\rangle eines Paares von Systemen (X,Y)(\mathsf{X},\mathsf{Y}) übereinstimmen müssen.

Anschaulich gesagt enthalten die reduzierten Zustände von X\mathsf{X} und Y\mathsf{Y} genau gleich viel Zufälligkeit, wenn sich das Paar (X,Y)(\mathsf{X},\mathsf{Y}) in einem reinen Zustand befindet. Diese Tatsache wird durch die Schmidt-Zerlegung offenbart: In beiden Fällen müssen die Eigenwerte der reduzierten Zustände mit den Quadraten der Schmidt-Koeffizienten des reinen Zustands übereinstimmen.

Unitäre Äquivalenz von Purifikationen

Wir können Schmidt-Zerlegungen nutzen, um eine grundlegend wichtige Tatsache über Purifikationen zu zeigen, die als unitäre Äquivalenz von Purifikationen bekannt ist.

Theorem

Unitäre Äquivalenz von Purifikationen: Seien X\mathsf{X} und Y\mathsf{Y} Systeme, und seien ψ\vert\psi\rangle und ϕ\vert\phi\rangle Quantenzustandsvektoren von (X,Y),(\mathsf{X},\mathsf{Y}), die beide denselben Zustand von X\mathsf{X} purifizieren. In Formeln:

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

für eine Dichtematrix ρ,\rho, die einen Zustand von X\mathsf{X} darstellt. Dann muss eine unitäre Operation UU auf Y\mathsf{Y} allein existieren, die die erste Purifikation in die zweite überführt:

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

Wir werden einige Implikationen dieses Theorems im weiteren Verlauf der Lektion besprechen, aber zunächst wollen wir sehen, wie es aus unserer vorherigen Diskussion über Schmidt-Zerlegungen folgt.

Unsere Annahme ist, dass ψ\vert\psi\rangle und ϕ\vert\phi\rangle Quantenzustandsvektoren eines Paares von Systemen (X,Y)(\mathsf{X},\mathsf{Y}) sind, die die Gleichung

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

für eine Dichtematrix ρ,\rho, die einen Zustand von X\mathsf{X} darstellt, erfüllen.

Betrachte eine Spektralzerlegung von ρ.\rho.

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

Hier ist {x0,,xn1}\{\vert x_0\rangle,\ldots,\vert x_{n-1}\rangle\} eine orthonormale Basis von Eigenvektoren von ρ.\rho. Indem wir der zuvor beschriebenen Vorgehensweise folgen, erhalten wir Schmidt-Zerlegungen sowohl für ψ\vert\psi\rangle als auch für ϕ\vert\phi\rangle in der folgenden Form.

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

In diesen Ausdrücken ist rr der Rang von ρ\rho und {u0,,ur1}\{\vert u_0\rangle,\ldots,\vert u_{r-1}\rangle\} und {v0,,vr1}\{\vert v_0\rangle,\ldots,\vert v_{r-1}\rangle\} sind orthonormale Mengen von Vektoren im Raum, der Y\mathsf{Y} entspricht.

Für zwei beliebige orthonormale Mengen im selben Raum mit gleicher Anzahl von Elementen gibt es stets eine unitäre Matrix, die die erste Menge in die zweite überführt, also können wir eine unitäre Matrix UU so wählen, dass Uua=vaU \vert u_a\rangle = \vert v_a\rangle für a=0,,r1.a = 0,\ldots,r-1. Um eine solche Matrix UU zu finden, können wir zunächst das Gram-Schmidt-Orthogonalisierungsverfahren verwenden, um unsere orthonormalen Mengen zu orthonormalen Basen {u0,,um1}\{\vert u_0\rangle,\ldots,\vert u_{m-1}\rangle\} und {v0,,vm1}\{\vert v_0\rangle,\ldots,\vert v_{m-1}\rangle\} zu erweitern, wobei mm die Dimension des Raums entsprechend Y\mathsf{Y} ist, und dann

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

zu nehmen.

Wir finden nun, dass

(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}

womit der Beweis abgeschlossen ist.

Hier sind einige von vielen interessanten Beispielen und Implikationen, die mit der unitären Äquivalenz von Purifikationen zusammenhängen. Ein weiteres entscheidend wichtiges Beispiel werden wir später in der Lektion im Kontext der Fidelität sehen, bekannt als Uhlmanns Theorem.

Superdense Coding

Im Superdense-Coding-Protokoll teilen sich Alice und Bob ein E-Bit, das bedeutet: Alice hält ein Qubit A,\mathsf{A}, Bob hält ein Qubit B,\mathsf{B}, und zusammen befindet sich das Paar (A,B)(\mathsf{A},\mathsf{B}) im Bell-Zustand ϕ+.\vert\phi^{+}\rangle. Das Protokoll beschreibt, wie Alice diesen gemeinsamen Zustand in einen der vier Bell-Zustände, ϕ+,\vert\phi^+\rangle, ϕ,\vert\phi^-\rangle, ψ+\vert\psi^+\rangle und ψ,\vert\psi^-\rangle, überführen kann, indem sie eine unitäre Operation auf ihr Qubit A\mathsf{A} anwendet. Sobald sie das getan hat, schickt sie A\mathsf{A} an Bob, der dann eine Messung auf dem Paar (A,B)(\mathsf{A},\mathsf{B}) durchführt, um herauszufinden, in welchem Bell-Zustand er sich befindet.

Für alle vier Bell-Zustände ist der reduzierte Zustand von Bobs Qubit B\mathsf{B} der vollständig gemischte Zustand.

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}

Durch die unitäre Äquivalenz von Purifikationen schließen wir unmittelbar, dass für jeden Bell-Zustand eine unitäre Operation auf Alices Qubit A\mathsf{A} allein existieren muss, die ϕ+\vert\phi^+\rangle in den gewählten Bell-Zustand überführt. Obwohl dies die genauen Details des Protokolls nicht offenbart, impliziert die unitäre Äquivalenz von Purifikationen unmittelbar, dass Superdense Coding möglich ist.

Wir können auch schlussfolgern, dass Verallgemeinerungen des Superdense Codings auf größere Systeme stets möglich sind, vorausgesetzt, wir ersetzen die Bell-Zustände durch eine beliebige orthonormale Basis von Purifikationen des vollständig gemischten Zustands.

Kryptografische Implikationen

Die unitäre Äquivalenz von Purifikationen hat Auswirkungen auf die Implementierung kryptografischer Primitive mittels Quanteninformation. Zum Beispiel zeigt die unitäre Äquivalenz von Purifikationen, dass es unmöglich ist, eine ideale Form des Bit-Commitments mit Quanteninformation zu implementieren.

Das Bit-Commitment-Primitiv umfasst zwei Teilnehmer, Alice und Bob (die einander nicht vertrauen), und verläuft in zwei Phasen.

  • Die erste Phase ist die Commit-Phase, in der sich Alice auf einen Binärwert b{0,1}b\in\{0,1\} festlegt. Dieses Commitment muss bindend sein, d. h. Alice kann ihre Meinung nicht ändern, sowie verbergend, d. h. Bob kann nicht erkennen, auf welchen Wert sich Alice festgelegt hat.
  • Die zweite Phase ist die Enthüllungsphase, in der der von Alice festgelegte Bit Bob bekannt wird und Bob überzeugt sein soll, dass es tatsächlich der zugesagte Wert war, der enthüllt wurde.

In anschaulichen, operativen Begriffen sollte die erste Phase des Bit-Commitments so funktionieren, als ob Alice einen Binärwert auf ein Stück Papier schreibt, das Papier in einem Safe einschließt und den Safe an Bob übergibt, während sie den Schlüssel behält. Alice hat sich auf den auf dem Papier geschriebenen Binärwert festgelegt, weil der Safe in Bobs Besitz ist (bindend), aber da Bob den Safe nicht öffnen kann, kann er nicht erkennen, auf welchen Wert sich Alice festgelegt hat (verbergend). Die zweite Phase sollte so funktionieren, als ob Alice Bob den Schlüssel zum Safe übergibt, so dass er ihn öffnen kann, um den Wert zu enthüllen, auf den sich Alice festgelegt hat.

Wie sich herausstellt, ist es unmöglich, ein perfektes Bit-Commitment-Protokoll allein durch Quanteninformation zu implementieren, da dies der unitären Äquivalenz von Purifikationen widerspricht. Hier ist eine übergeordnete Zusammenfassung eines Arguments, das dies zeigt.

Zunächst können wir annehmen, dass Alice und Bob nur unitäre Operationen durchführen oder neue initialisierte Systeme einführen, während das Protokoll ausgeführt wird. Die Tatsache, dass jeder Kanal eine Stinespring-Darstellung hat, erlaubt uns, diese Annahme zu machen.

Am Ende der Commit-Phase des Protokolls hält Bob ein zusammengesetztes System, das sich in einem von zwei Quantenzuständen befinden muss: ρ0,\rho_0, wenn Alice sich auf den Wert 00 festgelegt hat, und ρ1,\rho_1, wenn Alice sich auf den Wert 11 festgelegt hat. Damit das Protokoll perfekt verbergend ist, sollte Bob nicht in der Lage sein, zwischen diesen beiden Zuständen zu unterscheiden — es muss also ρ0=ρ1\rho_0 = \rho_1 gelten. (Andernfalls gäbe es eine Messung, die diese Zustände probabilistisch unterscheiden kann.)

Da Alice und Bob jedoch nur unitäre Operationen verwendet haben, muss der Zustand aller am Protokoll beteiligten Systeme zusammen nach der Commit-Phase in einem reinen Zustand sein. Sei insbesondere ψ0\vert\psi_0\rangle der reine Zustand aller am Protokoll beteiligten Systeme, wenn Alice sich auf 00 festlegt, und ψ1\vert\psi_1\rangle der reine Zustand aller am Protokoll beteiligten Systeme, wenn Alice sich auf 11 festlegt. Schreiben wir A\mathsf{A} und B\mathsf{B} für Alices und Bobs (möglicherweise zusammengesetzte) Systeme, dann gilt

ρ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}

Angesichts der Voraussetzung ρ0=ρ1\rho_0 = \rho_1 für ein perfekt verbergendes Protokoll stellen wir fest, dass ψ0\vert\psi_0\rangle und ψ1\vert\psi_1\rangle Purifikationen desselben Zustands sind — und somit muss nach der unitären Äquivalenz von Purifikationen eine unitäre Operation UU auf A\mathsf{A} allein existieren, so dass

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

Alice kann also ihr Commitment von 00 auf 11 ändern, indem sie UU auf A\mathsf{A} anwendet, oder von 11 auf 00 durch Anwenden von U,U^{\dagger}, und so scheitert das betrachtete hypothetische Protokoll vollständig, bindend zu sein.

Hughston-Jozsa-Wootters-Theorem

Die letzte Implikation der unitären Äquivalenz von Purifikationen, die wir in diesem Teil der Lektion besprechen, ist das folgende Theorem, bekannt als das Hughston-Jozsa-Wootters-Theorem. (Dies ist tatsächlich eine leicht vereinfachte Version des unter diesem Namen bekannten Theorems.)

Theorem

Hughston-Jozsa-Wootters: Seien X\mathsf{X} und Y\mathsf{Y} Systeme und sei ϕ\vert\phi\rangle ein Quantenzustandsvektor des Paares (X,Y).(\mathsf{X},\mathsf{Y}). Sei außerdem NN eine beliebige positive ganze Zahl, (p0,,pN1)(p_0,\ldots,p_{N-1}) ein Wahrscheinlichkeitsvektor und ψ0,,ψN1\vert\psi_0\rangle,\ldots,\vert\psi_{N-1}\rangle Quantenzustandsvektoren, die Zustände von X\mathsf{X} darstellen, so dass

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.

Es existiert eine (allgemeine) Messung {P0,,PN1}\{P_0,\ldots,P_{N-1}\} auf Y,\mathsf{Y}, so dass die folgenden zwei Aussagen wahr sind, wenn diese Messung auf Y\mathsf{Y} durchgeführt wird, wenn (X,Y)(\mathsf{X},\mathsf{Y}) im Zustand ϕ\vert\phi\rangle ist:

  1. Jedes Messergebnis a{0,,N1}a\in\{0,\ldots,N-1\} tritt mit Wahrscheinlichkeit pap_a auf.
  2. Bedingt auf das Messergebnis aa wird der Zustand von X\mathsf{X} zu ψa.\vert\psi_a\rangle.

Anschaulich gesagt besagt dieses Theorem, dass — solange wir einen reinen Zustand zweier Systeme haben — für jede Möglichkeit, über den reduzierten Zustand des ersten Systems als Konvexkombination reiner Zustände nachzudenken, eine Messung des zweiten Systems existiert, die diese Denkweise zur Realität macht. Beachte, dass die Zahl NN nicht notwendigerweise durch die Anzahl der klassischen Zustände von X\mathsf{X} oder Y\mathsf{Y} begrenzt ist. Zum Beispiel könnte N=1.000.000N = 1.000.000 sein, während X\mathsf{X} und Y\mathsf{Y} Qubits sind.

Wir werden dieses Theorem mithilfe der unitären Äquivalenz von Purifikationen beweisen, beginnend mit der Einführung eines neuen Systems Z\mathsf{Z} mit der klassischen Zustandsmenge {0,,N1}.\{0,\ldots,N-1\}. Betrachte die folgenden zwei Quantenzustandsvektoren des Tripels (X,Y,Z).(\mathsf{X},\mathsf{Y},\mathsf{Z}).

γ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}

Der erste Vektor γ0\vert\gamma_0\rangle ist einfach der gegebene Quantenzustandsvektor ϕ\vert\phi\rangle tensoriert mit 0\vert 0\rangle für das neue System Z.\mathsf{Z}. Für den zweiten Vektor γ1\vert\gamma_1\rangle haben wir im Wesentlichen einen Quantenzustandsvektor, der das Theorem trivial machen würde — zumindest wenn Y\mathsf{Y} durch Z\mathsf{Z} ersetzt würde — da eine Standardbasismessung auf Z\mathsf{Z} offensichtlich jedes Ergebnis aa mit Wahrscheinlichkeit pap_a liefert und bedingt auf dieses Ergebnis der Zustand von X\mathsf{X} zu ψa\vert\psi_a\rangle wird.

Indem wir das Paar (Y,Z)(\mathsf{Y},\mathsf{Z}) als ein einziges zusammengesetztes System betrachten, über das gespurt werden kann, um X\mathsf{X} zu erhalten, stellen wir fest, dass wir zwei verschiedene Purifikationen des Zustands

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

identifiziert haben.

Konkret haben wir für die erste:

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

und für die zweite:

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}

Es muss daher eine unitäre Operation UU auf (Y,Z)(\mathsf{Y},\mathsf{Z}) existieren, die

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

erfüllt, nach der unitären Äquivalenz von Purifikationen.

Mithilfe dieser unitären Operation UU können wir eine Messung implementieren, die die Anforderungen des Theorems erfüllt, wie das folgende Diagramm veranschaulicht. Mit anderen Worten: Wir führen das neue System Z\mathsf{Z} im Zustand 0\vert 0\rangle ein, wenden UU auf (Y,Z)(\mathsf{Y},\mathsf{Z}) an, was den Zustand von (X,Y,Z)(\mathsf{X},\mathsf{Y},\mathsf{Z}) von γ0\vert\gamma_0\rangle in γ1\vert\gamma_1\rangle überführt, und messen dann Z\mathsf{Z} mit einer Standardbasismessung, was — wie wir bereits gesehen haben — das gewünschte Verhalten ergibt.

Ein Quantenschaltkreis zur Implementierung einer Messung für das HSW-Theorem

Das gepunktete Rechteck in der Abbildung stellt eine Implementierung dieser Messung dar, die als Kollektion positiv semidefiniter Matrizen {P0,,PN1}\{P_0,\ldots,P_{N-1}\} wie folgt beschrieben werden kann.

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)