Перейти до основного вмісту

Очищення

Означення очищень

Почнімо з точного математичного означення очищень.

Означення

Нехай X\mathsf{X} — система в стані, описаному матрицею густини ρ,\rho, а ψ\vert\psi\rangle — вектор квантового стану пари (X,Y),(\mathsf{X},\mathsf{Y}), що дає ρ\rho після часткового сліду по Y\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, коли рівняння з означення виконується, — але ми здебільшого вживатимемо цей термін саме стосовно вектора квантового стану.

Термін очищення вживають і ширше: коли порядок систем змінено на протилежний, коли назви систем і станів інші (звісно), а також коли систем більше двох. Наприклад, якщо ψ\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)

виконується для матриці густини ρ,\rho, що описує стан системи (A,C),(\mathsf{A},\mathsf{C}), то ψ\vert\psi\rangle все одно називають очищенням ρ.\rho.

Утім, у цьому уроці ми зосередимося на конкретній формі з наведеного означення. Властивості й факти про очищення, відповідно до цього означення, як правило, можна узагальнити на більше двох систем, перевпорядкувавши системи й розбивши їх на дві складені системи — одна відіграє роль X,\mathsf{X}, інша — роль Y.\mathsf{Y}.

Існування очищень

Нехай X\mathsf{X} і Y\mathsf{Y} — довільні дві системи, а ρ\rho — заданий стан X.\mathsf{X}. Доведемо, що існує вектор квантового стану ψ\vert\psi\rangle пари (X,Y),(\mathsf{X},\mathsf{Y}), який очищує ρ\rho — тобто є очищенням ρ\rho — за умови, що система Y\mathsf{Y} достатньо велика. Зокрема, якщо Y\mathsf{Y} має принаймні стільки ж класичних станів, скільки X,\mathsf{X}, то очищення такого вигляду обов'язково існує для кожного стану ρ.\rho. Для деяких станів ρ\rho достатньо меншої кількості класичних станів Y\mathsf{Y}; загалом, rank(ρ)\operatorname{rank}(\rho) класичних станів Y\mathsf{Y} є необхідними й достатніми для існування вектора квантового стану пари (X,Y),(\mathsf{X},\mathsf{Y}), що очищує ρ.\rho.

Розглянемо спершу довільний запис ρ\rho у вигляді опуклої комбінації nn чистих станів, для будь-якого натурального числа n.n.

ρ=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}\rangle — вектори квантових станів системи X.\mathsf{X}.

Один із способів отримати такий запис — скористатися спектральною теоремою; тоді nn — кількість класичних станів X,\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. Це мінімальне значення n,n, при якому запис ρ\rho у наведеній формі існує.

Щоб уникнути непорозумінь: обраний запис ρ\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,,n1.0,\ldots,n-1. Якщо це не так, замість 0,,n10,\ldots,n-1 можна підставити будь-який вибір nn різних класичних станів Y.\mathsf{Y}. Перевірити, що це справді є очищенням ρ,\rho, нескладно — треба просто обчислити частковий слід, що можна зробити двома еквівалентними способами.

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} — обидва кубіти, і

ρ=(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} — додатні дійсні числа, що в сумі дають 1,1, а обидві множини {x0,,xr1}\{\vert x_0\rangle,\ldots,\vert x_{r-1}\rangle\} і {y0,,yr1}\{\vert y_0\rangle,\ldots,\vert y_{r-1}\rangle\} є ортонормованими.

Значення

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

у розкладанні Шмідта ψ\vert\psi\rangle називають коефіцієнтами Шмідта; вони визначені однозначно (з точністю до порядку) — це єдині додатні дійсні числа, що можуть входити до такого виразу для ψ.\vert\psi\rangle. Множини

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

натомість, визначені неоднозначно, і свобода вибору цих множин векторів буде пояснена нижче.

Тепер перевіримо, що заданий вектор квантового стану ψ\vert\psi\rangle справді має розкладання Шмідта, і заодно навчимося його знаходити.

Розглянемо спершу довільний (не обов'язково ортогональний) базис {x0,,xn1}\{\vert x_0\rangle, \ldots, \vert x_{n-1}\rangle\} векторного простору, що відповідає системі X.\mathsf{X}. Оскільки це базис, завжди існує єдиний набір векторів 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\}, то xa=a\vert x_a\rangle = \vert a\rangle для кожного a{0,,n1},a\in\{0,\ldots,n-1\}, і тоді

ψ=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} у стандартному базисі.

Важливо зазначити, що формула

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

для векторів z0,,zn1\vert z_0\rangle,\ldots,\vert z_{n-1}\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), але потрібна інша формула. Один із способів їх знайти — спочатку визначити вектори w0,,wn1\vert w_0\rangle,\ldots,\vert w_{n-1}\rangle так, щоб

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\}, після чого

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

Для заданого базису {x0,,xn1}\{\vert x_0\rangle,\ldots,\vert x_{n-1}\rangle\} векторного простору системи X\mathsf{X} однозначно визначені вектори z0,,zn1,\vert z_0\rangle,\ldots,\vert z_{n-1}\rangle, для яких виконується рівняння (1),(1), не обов'язково задовольняють якісь особливі властивості — навіть якщо {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),

відбудеться дещо цікаве. А саме: для єдиного набору {z0,,zn1,}\{\vert z_0\rangle,\ldots,\vert z_{n-1}\rangle,\} при якому рівняння (1)(1) виконується, цей набір обов'язково буде ортогональним.

Розглянемо докладніше спектральний розклад ρ.\rho.

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

Тут власні значення ρ\rho позначено p0,,pn1p_0,\ldots,p_{n-1} з огляду на те, що ρ\rho — матриця густини, тобто вектор власних значень (p0,,pn1)(p_0,\ldots,p_{n-1}) утворює вектор ймовірностей, а {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) виконується, є ортогональним, почнемо з обчислення часткового сліду.

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=0,p_a = 0, а потім записати za=paya\vert z_a\rangle = \sqrt{p_a}\vert y_a\rangle для одиничного вектора ya\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=0.p_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.

Оскільки

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

для a=0,,r1,a=0,\ldots,r-1, визначимо одиничні вектори 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}},

так що za=paya\vert z_a\rangle = \sqrt{p_a}\vert y_a\rangle для кожного a{0,,r1}.a\in\{0,\ldots,r-1\}. Оскільки вектори {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\} визначаються однозначно.

Ситуація симетрична відносно двох систем, тому можна альтернативно обрати {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\} визначаться однозначно.

Зауваж, однак, що щойно одну з множин вибрано як множину власних векторів відповідного зведеного стану, як описано вище, інша однозначно визначена — тобто їх не можна обирати незалежно.

Хоча до цього в серії ми більше не повернемося, варто зазначити, що ненульові власні значення p0,,pr1p_0,\ldots,p_{r-1} зведеного стану TrX(ψψ)\operatorname{Tr}_{\mathsf{X}}(\vert\psi\rangle\langle\psi\vert) завжди збігаються з ненульовими власними значеннями зведеного стану TrY(ψψ)\operatorname{Tr}_{\mathsf{Y}}(\vert\psi\rangle\langle\psi\vert) для будь-якого чистого стану ψ\vert\psi\rangle пари систем (X,Y).(\mathsf{X},\mathsf{Y}).

Інтуїтивно кажучи, зведені стани X\mathsf{X} і Y\mathsf{Y} мають рівно однакову кількість випадковості, коли пара (X,Y)(\mathsf{X},\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)

для деякої матриці густини ρ,\rho, що описує стан X.\mathsf{X}. Тоді обов'язково існує унітарна операція UU на одній лише системі Y,\mathsf{Y}, що переводить перше очищення в друге:

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

для деякої матриці густини ρ,\rho, що описує стан X.\mathsf{X}.

Розглянемо спектральний розклад ρ.\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}.

Для будь-яких двох ортонормованих множин одного простору з однаковою кількістю елементів завжди існує унітарна матриця, що переводить першу множину в другу, тому можна обрати унітарну матрицю UU так, що Uua=vaU \vert u_a\rangle = \vert v_a\rangle для a=0,,r1.a = 0,\ldots,r-1. Зокрема, щоб знайти таку матрицю U,U, можна спочатку за допомогою процесу ортогоналізації Грама–Шмідта розширити ортонормовані множини до ортонормованих базисів {u0,,um1}\{\vert u_0\rangle,\ldots,\vert u_{m-1}\rangle\} і {v0,,vm1},\{\vert v_0\rangle,\ldots,\vert v_{m-1}\rangle\}, де mm — вимірність простору системи Y,\mathsf{Y}, а потім покласти

U=a=0m1vaua.U = \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}

що завершує доведення.

Нижче наведено лише кілька з багатьох цікавих прикладів і наслідків, пов'язаних з унітарною еквівалентністю очищень. Ще один критично важливий — в контексті точності — ми побачимо пізніше в уроці; він відомий як теорема Ulmann.

Надгусте кодування

У протоколі надгустого кодування Аліса і Боб спільно використовують е-біт: Аліса тримає кубіт A,\mathsf{A}, Боб тримає кубіт B,\mathsf{B}, і разом пара (A,B)(\mathsf{A},\mathsf{B}) перебуває в стані Белла ϕ+.\vert\phi^{+}\rangle. Протокол описує, як Аліса може перевести цей спільний стан у будь-який із чотирьох станів Белла — ϕ+,\vert\phi^+\rangle, ϕ,\vert\phi^-\rangle, ψ+\vert\psi^+\rangle і ψ\vert\psi^-\rangle — застосувавши унітарну операцію до свого кубіта A.\mathsf{A}. Після цього вона відправляє A\mathsf{A} Бобу, і Боб виконує вимірювання пари (A,B),(\mathsf{A},\mathsf{B}), щоб з'ясувати, у якому стані Белла вона перебуває.

Для всіх чотирьох станів Белла зведений стан кубіта Боба 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}

З унітарної еквівалентності очищень безпосередньо випливає, що для кожного стану Белла має існувати унітарна операція на одному лише кубіті Аліси A,\mathsf{A}, що переводить ϕ+\vert\phi^+\rangle у вибраний стан Белла. Хоча це не розкриває точних деталей протоколу, унітарна еквівалентність очищень безпосередньо означає, що надгусте кодування можливе.

Можна також зробити висновок, що узагальнення надгустого кодування на більші системи завжди можливі за умови, що стани Белла замінити будь-яким ортонормованим базисом із очищень цілком змішаного стану.

Криптографічні наслідки

Унітарна еквівалентність очищень має наслідки для реалізації криптографічних примітивів за допомогою квантової інформації. Зокрема, унітарна еквівалентність очищень доводить, що неможливо реалізувати ідеальну форму зобов'язання бітом (bit commitment) за допомогою квантової інформації.

Примітив зобов'язання бітом передбачає двох учасників — Аліса і Боб (які не довіряють один одному) — і складається з двох фаз.

  • Перша фаза — це фаза зобов'язання (commit), протягом якої Аліса бере зобов'язання щодо бінарного значення b{0,1}.b\in\{0,1\}. Це зобов'язання має бути прив'язуючим (binding) — тобто Аліса не може передумати — а також приховуючим (concealing) — тобто Боб не може дізнатися, яке значення вибрала Аліса.
  • Друга фаза — це фаза розкриття (reveal), у якій біт, зафіксований Алісою, стає відомим Бобу, і Боб має бути переконаний, що саме зафіксоване значення й було розкрито.

В інтуїтивно-операційних термінах перша фаза зобов'язання бітом має працювати так, ніби Аліса записує бінарне значення на папері, замикає папір у сейфі і передає сейф Бобу, залишаючи ключ собі. Аліса взяла зобов'язання щодо бінарного значення, записаного на папері, бо сейф знаходиться у Боба (тобто це прив'язуюче), а оскільки Боб не може відкрити сейф, він не може дізнатися, яке значення вибрала Аліса (тобто це приховуюче). Друга фаза має працювати так, ніби Аліса передає Бобу ключ від сейфа, щоб він міг відкрити його й розкрити значення, до якого зобов'язалася Аліса.

Як виявляється, неможливо реалізувати досконалий протокол зобов'язання бітом лише засобами квантової інформації, оскільки це суперечить унітарній еквівалентності очищень. Нижче наведено стислий виклад аргументу, що це доводить.

Для початку можна припустити, що в ході виконання протоколу Аліса і Боб виконують лише унітарні операції або вводять нові ініціалізовані системи. Це припущення правомірне завдяки тому, що кожний канал має представлення Штайнспрінга.

Наприкінці фази зобов'язання протоколу Боб має у своєму розпорядженні деяку складену систему, яка повинна перебувати в одному з двох квантових станів: ρ0\rho_0 — якщо Аліса зобов'язалася щодо значення 00, і ρ1\rho_1 — якщо щодо значення 1.1. Щоб протокол був досконало приховуючим, Боб не повинен уміти розрізняти ці два стани — тобто має виконуватися ρ0=ρ1.\rho_0 = \rho_1. (Інакше існувало б вимірювання, яке ймовірнісно розрізняє ці стани.)

Однак оскільки Аліса і Боб використовували лише унітарні операції, стан усіх систем, задіяних у протоколі, після фази зобов'язання має бути чистим. Зокрема, припустимо, що ψ0\vert\psi_0\rangle — це чистий стан усіх систем протоколу, коли Аліса зобов'язується щодо 00, а ψ1\vert\psi_1\rangle — чистий стан усіх систем, коли Аліса зобов'язується щодо 1.1. Позначивши 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 є очищеннями одного й того самого стану — а отже, за унітарною еквівалентністю очищень, існує унітарна операція UU над A\mathsf{A} така, що

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

Таким чином, Аліса може вільно змінити своє зобов'язання з 00 на 11, застосувавши UU до A,\mathsf{A}, або з 11 на 00, застосувавши UU^{\dagger}, тому розглянутий гіпотетичний протокол абсолютно не є прив'язуючим.

Теорема Г'юстона-Йожи-Вуттерса

Останній наслідок унітарної еквівалентності очищень, який ми розглянемо в цій частині уроку, — це теорема, відома як теорема Г'юстона-Йожи-Вуттерса. (Насправді це дещо спрощена формулювання теореми, що носить цю назву.)

Теорема

Г'юстон-Йожи-Вуттерс: Нехай 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}\rangle — вектори квантових станів, що описують стани системи X\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.

Існує (загальне) вимірювання {P0,,PN1}\{P_0,\ldots,P_{N-1}\} над Y\mathsf{Y} таке, що при виконанні цього вимірювання над Y\mathsf{Y}, коли (X,Y)(\mathsf{X},\mathsf{Y}) перебуває у стані ϕ\vert\phi\rangle, виконуються такі два твердження:

  1. Кожний результат вимірювання a{0,,N1}a\in\{0,\ldots,N-1\} з'являється з імовірністю pap_a.
  2. За умови отримання результату вимірювання aa стан системи X\mathsf{X} стає ψa.\vert\psi_a\rangle.

Інтуїтивно ця теорема стверджує: якщо ми маємо чистий стан двох систем, то для будь-якого способу представлення зведеного стану першої системи як опуклої комбінації чистих станів існує вимірювання другої системи, яке фактично робить цей спосіб мислення про першу систему реальністю. Зауваж, що число NN не обов'язково обмежене кількістю класичних станів X\mathsf{X} чи Y.\mathsf{Y}. Наприклад, може бути N=1000000N = 1\,000\,000, тоді як X\mathsf{X} і Y\mathsf{Y} є кубітами.

Доведемо цю теорему, використовуючи унітарну еквівалентність очищень, починаючи з введення нової системи Z\mathsf{Z}, множина класичних станів якої — {0,,N1}.\{0,\ldots,N-1\}. Розглянемо такі два вектори квантових станів трійки (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}

Перший вектор γ0\vert\gamma_0\rangle — це просто заданий вектор квантового стану ϕ\vert\phi\rangle, доповнений тензорним добутком з 0\vert 0\rangle нової системи Z.\mathsf{Z}. Щодо другого вектора γ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.

Зокрема, для першого маємо

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

а для другого маємо

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}

Тому за унітарною еквівалентністю очищень має існувати унітарна операція UU над (Y,Z)(\mathsf{Y},\mathsf{Z}), що задовольняє

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

Використовуючи цю унітарну операцію UU, ми можемо реалізувати вимірювання, яке задовольняє вимоги теореми, як показано на діаграмі нижче. Словами: ми вводимо нову систему Z\mathsf{Z}, ініціалізовану у стан 0\vert 0\rangle, застосовуємо UU до (Y,Z)(\mathsf{Y},\mathsf{Z}), що переводить стан (X,Y,Z)(\mathsf{X},\mathsf{Y},\mathsf{Z}) з γ0\vert\gamma_0\rangle у γ1\vert\gamma_1\rangle, а потім вимірюємо Z\mathsf{Z} стандартним базисним вимірюванням, яке, як ми вже бачили, дає потрібну поведінку.

A quantum circuit implementation of a measurement for the HSW theorem

Пунктирний прямокутник на рисунку зображує реалізацію цього вимірювання, яке можна описати як сукупність позитивно напіввизначених матриць {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)