Sari la conținutul principal

Formulări matematice ale măsurătorilor

Lecția începe cu două descrieri matematice echivalente ale măsurătorilor:

  1. Măsurătorile generale pot fi descrise prin colecții de matrice, câte una pentru fiecare rezultat posibil, într-un mod care generalizează descrierea măsurătorilor proiective.
  2. Măsurătorile generale pot fi descrise ca canale ale căror ieșiri sunt întotdeauna stări clasice (reprezentate prin matrice de densitate diagonale).

Ne vom limita atenția la măsurători cu un număr finit de rezultate posibile. Deși este posibil să definim măsurători cu infinit de rezultate posibile, acestea sunt mult mai rar întâlnite în contextul calculului și al prelucrării informației și necesită, de asemenea, matematică suplimentară (și anume teoria măsurii) pentru a fi formalizate corect.

Atenția noastră inițială se va îndrepta spre așa-numitele măsurători distructive, în care ieșirea măsurătorii este doar un rezultat clasic — fără nicio specificație a stării cuantice post-măsurătoare a sistemului măsurat. Intuitiv, ne putem imagina că o astfel de măsurătoare distruge sistemul cuantic în sine sau că sistemul este imediat eliminat odată ce măsurătoarea este efectuată. Mai târziu în lecție ne vom lărgi perspectiva și vom considera măsurătorile nedistructive, în care există atât un rezultat clasic, cât și o stare cuantică post-măsurătoare a sistemului măsurat.

Măsurători ca colecții de matrice

Fie X\mathsf{X} un sistem ce urmează a fi măsurat și presupunem, pentru simplitate, că mulțimea stărilor clasice ale lui X\mathsf{X} este {0,,n1}\{0,\ldots, n-1\} pentru un număr întreg pozitiv n,n, astfel încât matricele de densitate care reprezintă stările cuantice ale lui X\mathsf{X} sunt matrice de dimensiune n×nn\times n. Nu vom avea de fapt prea mare nevoie să ne referim la stările clasice ale lui X,\mathsf{X}, dar va fi convenabil să ne referim la n,n, numărul de stări clasice ale lui X.\mathsf{X}. Vom presupune, de asemenea, că rezultatele posibile ale măsurătorii sunt numerele întregi 0,,m10,\ldots,m-1 pentru un număr întreg pozitiv m.m.

Rețineți că folosim aceste denumiri doar pentru a simplifica lucrurile; este simplu să generalizăm tot ceea ce urmează la alte mulțimi finite de stări clasice și rezultate ale măsurătorilor, redenumindu-le după cum dorim.

Măsurători proiective

Reamintim că o măsurătoare proiectivă este descrisă de o colecție de matrice de proiecție care sumează la matricea identitate. În simboluri,

{Π0,,Πm1}\{\Pi_0,\ldots,\Pi_{m-1}\}

descrie o măsurătoare proiectivă a lui X\mathsf{X} dacă fiecare Πa\Pi_a este o matrice de proiecție de dimensiune n×nn\times n și este îndeplinită următoarea condiție.

Π0++Πm1=IX\Pi_0 + \cdots + \Pi_{m-1} = \mathbb{I}_{\mathsf{X}}

Când o astfel de măsurătoare este efectuată pe un sistem X\mathsf{X} aflat într-o stare descrisă de un vector de stare cuantică ψ,\vert\psi\rangle, fiecare rezultat aa este obținut cu probabilitate egală cu Πaψ2.\|\Pi_a\vert\psi\rangle\|^2. De asemenea, starea post-măsurătoare a lui X\mathsf{X} se obține prin normalizarea vectorului Πaψ,\Pi_a\vert\psi\rangle, dar deocamdată ignorăm starea post-măsurătoare.

Dacă starea lui X\mathsf{X} este descrisă de o matrice de densitate ρ\rho în loc de un vector de stare cuantică ψ,\vert\psi\rangle, atunci putem exprima alternativ probabilitatea de a obține rezultatul aa ca Tr(Πaρ).\operatorname{Tr}(\Pi_a \rho).

Dacă ρ=ψψ\rho = \vert \psi\rangle\langle\psi\vert este o stare pură, atunci cele două expresii sunt egale:

Tr(Πaρ)=Tr(Πaψψ)=ψΠaψ=ψΠaΠaψ=Πaψ2.\operatorname{Tr}(\Pi_a \rho) = \operatorname{Tr}(\Pi_a \vert \psi\rangle\langle\psi \vert) = \langle \psi \vert \Pi_a \vert \psi \rangle = \langle \psi \vert \Pi_a \Pi_a \vert \psi \rangle = \|\Pi_a\vert\psi\rangle\|^2.

Aici folosim proprietatea ciclică a urmei pentru a doua egalitate, iar pentru a treia egalitate folosim faptul că fiecare Πa\Pi_a este o matrice de proiecție și, prin urmare, satisface Πa2=Πa.\Pi_a^2 = \Pi_a.

În general, dacă ρ\rho este o combinație convexă

ρ=k=0N1pkψkψk\rho = \sum_{k = 0}^{N-1} p_k \vert \psi_k\rangle\langle \psi_k \vert

de stări pure, atunci expresia Tr(Πaρ)\operatorname{Tr}(\Pi_a \rho) coincide cu probabilitatea medie pentru rezultatul a,a, datorită faptului că această expresie este liniară în ρ.\rho.

Tr(Πaρ)=k=0N1pkTr(Πaψkψk)=k=0N1pkΠaψk2\operatorname{Tr}(\Pi_a \rho) = \sum_{k = 0}^{N-1} p_k \operatorname{Tr}(\Pi_a \vert \psi_k\rangle\langle\psi_k\vert) = \sum_{k = 0}^{N-1} p_k \|\Pi_a\vert\psi_k\rangle\|^2

Măsurători generale

O descriere matematică a măsurătorilor generale se obține prin relaxarea definiției măsurătorilor proiective. Mai precis, permitem matricelor din colecția ce descrie măsurătoarea să fie matrice pozitiv semidefinite arbitrare, nu neapărat proiecții. (Proiecțiile sunt întotdeauna pozitiv semidefinite; ele pot fi definite alternativ ca matrice pozitiv semidefinite ale căror valori proprii sunt toate fie 0, fie 1.)

În particular, o măsurătoare generală a unui sistem X\mathsf{X} cu rezultatele 0,,m10,\ldots,m-1 este specificată de o colecție de matrice pozitiv semidefinite {P0,,Pm1}\{P_0,\ldots,P_{m-1}\} ale căror linii și coloane corespund stărilor clasice ale lui X\mathsf{X} și care satisfac condiția

P0++Pm1=IX.P_0 + \cdots + P_{m-1} = \mathbb{I}_{\mathsf{X}}.

Dacă sistemul X\mathsf{X} este măsurat în timp ce se află într-o stare descrisă de matricea de densitate ρ,\rho, atunci fiecare rezultat a{0,,m1}a\in\{0,\ldots,m-1\} apare cu probabilitate Tr(Paρ).\operatorname{Tr}(P_a \rho).

Așa cum este firesc să cerem, vectorul probabilităților rezultatelor

(Tr(P0ρ),,Tr(Pm1ρ))\bigl(\operatorname{Tr}(P_0 \rho),\ldots,\operatorname{Tr}(P_{m-1} \rho)\bigr)

al unei măsurători generale formează întotdeauna un vector de probabilitate, indiferent de alegerea matricei de densitate ρ.\rho. Următoarele două observații demonstrează că acesta este cazul.

  1. Fiecare valoare Tr(Paρ)\operatorname{Tr}(P_a \rho) trebuie să fie nenegativă, datorită faptului că urma produsului oricăror două matrice pozitiv semidefinite este întotdeauna nenegativă:

    Q,R0  Tr(QR)0.Q, R \geq 0 \; \Rightarrow \: \operatorname{Tr}(QR) \geq 0.

    Un mod de a argumenta acest fapt este să folosim descompunerile spectrale ale lui QQ și RR împreună cu proprietatea ciclică a urmei pentru a exprima urma produsului QRQR ca o sumă de numere reale nenegative, care trebuie, prin urmare, să fie nenegativă.

  2. Condiția P0++Pm1=IXP_0 + \cdots + P_{m-1} = \mathbb{I}_{\mathsf{X}} împreună cu liniaritatea urmei asigură că probabilitățile sumează la 1.1.

    a=0m1Tr(Paρ)=Tr(a=0m1Paρ)=Tr(Iρ)=Tr(ρ)=1\sum_{a = 0}^{m-1} \operatorname{Tr}(P_a \rho) = \operatorname{Tr}\Biggl(\sum_{a = 0}^{m-1} P_a \rho\Biggr) = \operatorname{Tr}(\mathbb{I}\rho) = \operatorname{Tr}(\rho) = 1

Exemplul 1: orice măsurătoare proiectivă

Proiecțiile sunt întotdeauna pozitiv semidefinite, deci orice măsurătoare proiectivă este un exemplu de măsurătoare generală.

De exemplu, o măsurătoare în baza standard a unui Qubit poate fi reprezentată prin {P0,P1}\{P_0,P_1\} unde

P0=00=(1000)șiP1=11=(0001).P_0 = \vert 0\rangle\langle 0\vert = \begin{pmatrix} 1 & 0 \\ 0 & 0 \end{pmatrix} \quad\text{și}\quad P_1 = \vert 1\rangle\langle 1\vert = \begin{pmatrix} 0 & 0 \\ 0 & 1 \end{pmatrix}.

Măsurarea unui Qubit în starea ρ\rho produce probabilitățile rezultatelor după cum urmează.

Prob(outcome=0)=Tr(P0ρ)=Tr(00ρ)=0ρ0Prob(outcome=1)=Tr(P1ρ)=Tr(11ρ)=1ρ1\begin{aligned} \operatorname{Prob}(\text{outcome} = 0) & = \operatorname{Tr}(P_0 \rho) = \operatorname{Tr}\bigl(\vert 0\rangle\langle 0\vert \rho\bigr) = \langle 0\vert \rho \vert 0 \rangle \\[1mm] \operatorname{Prob}(\text{outcome} = 1) & = \operatorname{Tr}(P_1 \rho) = \operatorname{Tr}\bigl(\vert 1\rangle\langle 1\vert\rho\bigr) = \langle 1 \vert \rho \vert 1 \rangle \end{aligned}

Exemplul 2: o măsurătoare ne-proiectivă a unui Qubit

Fie X\mathsf{X} un Qubit și definim două matrice după cum urmează.

P0=(23131313)P1=(13131323)P_0 = \begin{pmatrix} \frac{2}{3} & \frac{1}{3}\\[2mm] \frac{1}{3} & \frac{1}{3} \end{pmatrix} \qquad P_1 = \begin{pmatrix} \frac{1}{3} & -\frac{1}{3}\\[2mm] -\frac{1}{3} & \frac{2}{3} \end{pmatrix}

Ambele sunt matrice pozitiv semidefinite: sunt Hermitian și, în ambele cazuri, valorile proprii sunt 1/2±5/6,1/2 \pm \sqrt{5}/6, care sunt ambele pozitive. De asemenea, avem că P0+P1=I,P_0 + P_1 = \mathbb{I}, și prin urmare {P0,P1}\{P_0,P_1\} descrie o măsurătoare.

Dacă starea lui X\mathsf{X} este descrisă de o matrice de densitate ρ\rho și efectuăm această măsurătoare, atunci probabilitatea de a obține rezultatul 00 este Tr(P0ρ)\operatorname{Tr}(P_0 \rho) și probabilitatea de a obține rezultatul 11 este Tr(P1ρ).\operatorname{Tr}(P_1 \rho). De exemplu, dacă ρ=++,\rho = \vert + \rangle \langle + \vert, probabilitățile pentru cele două rezultate 00 și 11 sunt următoarele.

Tr(P0ρ)=Tr((23131313)(12121212))=(2312+1312)+(1312+1312)=12+13=56Tr(P1ρ)=Tr((13131323)(12121212))=(13121312)+(1312+2312)=0+16=16\begin{aligned} \operatorname{Tr}(P_0 \rho) & = \operatorname{Tr}\left( \begin{pmatrix} \frac{2}{3} & \frac{1}{3}\\[2mm] \frac{1}{3} & \frac{1}{3} \end{pmatrix} \begin{pmatrix} \frac{1}{2} & \frac{1}{2}\\[2mm] \frac{1}{2} & \frac{1}{2} \end{pmatrix} \right)\\[4mm] & = \biggl(\frac{2}{3} \cdot \frac{1}{2} + \frac{1}{3} \cdot \frac{1}{2}\biggr) + \biggl(\frac{1}{3}\cdot\frac{1}{2} + \frac{1}{3}\cdot\frac{1}{2}\biggr)\\ & = \frac{1}{2} + \frac{1}{3} = \frac{5}{6}\\[4mm] \operatorname{Tr}(P_1 \rho) & = \operatorname{Tr}\left( \begin{pmatrix} \frac{1}{3} & -\frac{1}{3}\\[2mm] -\frac{1}{3} & \frac{2}{3} \end{pmatrix} \begin{pmatrix} \frac{1}{2} & \frac{1}{2}\\[2mm] \frac{1}{2} & \frac{1}{2} \end{pmatrix} \right)\\[4mm] & = \biggl(\frac{1}{3} \cdot \frac{1}{2} - \frac{1}{3} \cdot \frac{1}{2}\biggr) + \biggl(-\frac{1}{3}\cdot\frac{1}{2} + \frac{2}{3}\cdot\frac{1}{2}\biggr)\\ & = 0 + \frac{1}{6} = \frac{1}{6} \end{aligned}

Exemplul 3: măsurătoarea tetraedrică

Definim patru vectori de stare cuantică monoQubit după cum urmează.

ϕ0=0ϕ1=130+231ϕ2=130+23e2πi/31ϕ3=130+23e2πi/31\begin{aligned} \vert\phi_0\rangle & = \vert 0 \rangle\\ \vert\phi_1\rangle & = \frac{1}{\sqrt{3}}\vert 0 \rangle + \sqrt{\frac{2}{3}} \vert 1\rangle \\ \vert\phi_2\rangle & = \frac{1}{\sqrt{3}}\vert 0 \rangle + \sqrt{\frac{2}{3}} e^{2\pi i/3} \vert 1\rangle \\ \vert\phi_3\rangle & = \frac{1}{\sqrt{3}}\vert 0 \rangle + \sqrt{\frac{2}{3}} e^{-2\pi i/3} \vert 1\rangle \end{aligned}

Aceste patru stări sunt uneori cunoscute sub denumirea de stări tetraedrice deoarece sunt vârfurile unui tetraedru regulat înscris în sfera Bloch.

Illustration of a tetrahedron inscribed in the Bloch sphere

Coordonatele carteziene ale acestor patru stări pe sfera Bloch sunt

(0,0,1),(223,0,13),(23,23,13),(23,23,13),(0,0,1),\\[2mm] \left( \frac{2\sqrt{2}}{3} , 0 , -\frac{1}{3} \right),\\[1mm] \left( -\frac{\sqrt{2}}{3} , \sqrt{\frac{2}{3}} , -\frac{1}{3} \right),\\[1mm] \left( -\frac{\sqrt{2}}{3} , -\sqrt{\frac{2}{3}} , -\frac{1}{3} \right),

ceea ce poate fi verificat exprimând reprezentările matricelor de densitate ale acestor stări ca combinații liniare de matrice Pauli.

ϕ0ϕ0=(1000)=I+σz2\vert \phi_0 \rangle\langle \phi_0 \vert = \begin{pmatrix} 1 & 0\\[1mm] 0 & 0 \end{pmatrix} = \frac{\mathbb{I} + \sigma_z}{2} ϕ1ϕ1=(13232323)=I+223σx13σz2\vert \phi_1 \rangle\langle \phi_1 \vert = \begin{pmatrix} \frac{1}{3} & \frac{\sqrt{2}}{3} \\[2mm] \frac{\sqrt{2}}{3} & \frac{2}{3} \end{pmatrix} = \frac{\mathbb{I} + \frac{2\sqrt{2}}{3} \sigma_x - \frac{1}{3}\sigma_z}{2} ϕ2ϕ2=(13132i6132+i623)=I23σx+23σy13σz2\vert \phi_2 \rangle\langle \phi_2 \vert = \begin{pmatrix} \frac{1}{3} & -\frac{1}{3\sqrt{2}} - \frac{i}{\sqrt{6}} \\[2mm] -\frac{1}{3\sqrt{2}} + \frac{i}{\sqrt{6}} & \frac{2}{3} \end{pmatrix} = \frac{\mathbb{I} - \frac{\sqrt{2}}{3} \sigma_x + \sqrt{\frac{2}{3}} \sigma_y - \frac{1}{3}\sigma_z}{2} ϕ3ϕ3=(13132+i6132i623)=I23σx23σy13σz2\vert \phi_3 \rangle\langle \phi_3 \vert = \begin{pmatrix} \frac{1}{3} & -\frac{1}{3\sqrt{2}} + \frac{i}{\sqrt{6}} \\[2mm] -\frac{1}{3\sqrt{2}} - \frac{i}{\sqrt{6}} & \frac{2}{3} \end{pmatrix} = \frac{\mathbb{I} - \frac{\sqrt{2}}{3} \sigma_x - \sqrt{\frac{2}{3}} \sigma_y - \frac{1}{3}\sigma_z}{2}

Aceste patru stări sunt perfect distribuite pe sfera Bloch, fiecare aflată la aceeași distanță față de celelalte trei, unghiurile dintre oricare două dintre ele fiind întotdeauna egale.

Acum să definim o măsurătoare {P0,P1,P2,P3}\{P_0,P_1,P_2,P_3\} a unui Qubit setând PaP_a după cum urmează pentru fiecare a=0,,3.a=0,\ldots,3.

Pa=ϕaϕa2P_a = \frac{\vert\phi_a\rangle\langle\phi_a\vert}{2}

Putem verifica că aceasta este o măsurătoare validă după cum urmează.

  1. Fiecare PaP_a este evident pozitiv semidefinit, fiind o stare pură împărțită la o jumătate. Adică, fiecare este o matrice Hermitian și are o valoare proprie egală cu 1/21/2 și toate celelalte valori proprii zero.
  2. Suma acestor matrice este matricea identitate: P0+P1+P2+P3=I.P_0 + P_1 + P_2 + P_3 = \mathbb{I}. Exprimarea acestor matrice ca combinații liniare de matrice Pauli face această verificare simplă.

Măsurători ca și canale

O a doua modalitate de a descrie măsurătorile în termeni matematici este ca și canale.

Informația clasică poate fi privită ca un caz special al informației cuantice, în măsura în care putem identifica stările probabilistice cu matricele de densitate diagonale. Deci, în termeni operaționali, putem privi măsurătorile ca și canale ale căror intrări sunt matrice ce descriu stările oricărui sistem măsurat și ale căror ieșiri sunt matrice de densitate diagonale ce descriu distribuția rezultatelor măsurătorii.

Vom vedea în curând că orice canal cu această proprietate poate fi întotdeauna scris într-o formă simplă, canonică, care se leagă direct de descrierea măsurătorilor ca și colecții de matrice pozitiv semidefinite. Invers, dată o măsurătoare arbitrară ca o colecție de matrice, există întotdeauna un canal valid cu proprietatea ieșirii diagonale care descrie măsurătoarea dată, conform sugestiei din paragraful anterior. Punând aceste observații împreună, concluzionăm că cele două descrieri ale măsurătorilor generale sunt echivalente.

Înainte de a merge mai departe, să fim mai preciși cu privire la măsurătoare, la modul în care o privim ca pe un canal și la ipotezele pe care le facem despre ea.

Ca și înainte, vom presupune că X\mathsf{X} este sistemul ce urmează a fi măsurat și că rezultatele posibile ale măsurătorii sunt numerele întregi 0,,m10,\ldots,m-1 pentru un număr întreg pozitiv m.m. Fie Y\mathsf{Y} sistemul care stochează rezultatele măsurătorii, deci mulțimea sa de stări clasice este {0,,m1},\{0,\ldots,m-1\}, și reprezentăm măsurătoarea ca pe un canal denumit Φ\Phi de la X\mathsf{X} la Y.\mathsf{Y}. Presupunerea noastră este că Y\mathsf{Y} este clasic — adică, indiferent de starea cu care pornim pentru X,\mathsf{X}, starea lui Y\mathsf{Y} pe care o obținem este reprezentată de o matrice de densitate diagonală.

Putem exprima în termeni matematici că ieșirea lui Φ\Phi este întotdeauna diagonală în felul următor. Mai întâi definim canalul complet de dephasing Δm\Delta_m pe Y.\mathsf{Y}.

Δm(σ)=a=0m1aσaaa\Delta_m(\sigma) = \sum_{a = 0}^{m-1} \langle a \vert \sigma \vert a\rangle \,\vert a\rangle\langle a\vert

Acest canal este analog canalului complet de dephasing al unui Qubit Δ\Delta din lecția anterioară. Ca aplicație liniară, anulează toate intrările ne-diagonale ale unei matrice de intrare și lasă diagonala neschimbată.

Iar acum, un mod simplu de a exprima că o matrice de densitate dată σ\sigma este diagonală este prin ecuația σ=Δm(σ).\sigma = \Delta_m(\sigma). Cu alte cuvinte, anularea tuturor intrărilor ne-diagonale ale unei matrice de densitate nu are niciun efect dacă și numai dacă intrările ne-diagonale erau deja zero. Canalul Φ\Phi satisface, prin urmare, presupunerea noastră — că Y\mathsf{Y} este clasic — dacă și numai dacă

Φ(ρ)=Δm(Φ(ρ))\Phi(\rho) = \Delta_m(\Phi(\rho))

pentru orice matrice de densitate ρ\rho ce reprezintă o stare a lui X.\mathsf{X}.

Echivalența formulărilor

De la canale la matrice

Fie un canal de la X\mathsf{X} la Y\mathsf{Y} cu proprietatea că

Φ(ρ)=Δm(Φ(ρ))\Phi(\rho) = \Delta_m(\Phi(\rho))

pentru orice matrice de densitate ρ.\rho. Aceasta poate fi exprimată alternativ după cum urmează.

Φ(ρ)=a=0m1aΦ(ρ)aaa(1)\Phi(\rho) = \sum_{a = 0}^{m-1} \langle a \vert \Phi(\rho) \vert a\rangle\, \vert a\rangle\langle a \vert \tag{1}

Ca orice canal, putem exprima Φ\Phi în formă Kraus pentru un anumit mod de alegere a matricelor Kraus A0,,AN1.A_0,\ldots,A_{N-1}.

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

Aceasta ne oferă o expresie alternativă pentru intrările diagonale ale lui Φ(ρ) ⁣:\Phi(\rho)\!:

aΦ(ρ)a=k=0N1aAkρAka=k=0N1Tr(AkaaAkρ)=Tr(Paρ)\begin{aligned} \langle a \vert \Phi(\rho) \vert a\rangle & = \sum_{k = 0}^{N-1} \langle a \vert A_k \rho A_k^{\dagger} \vert a\rangle \\ & = \sum_{k = 0}^{N-1} \operatorname{Tr}\bigl( A_k^{\dagger} \vert a\rangle\langle a \vert A_k \rho\bigr)\\ & = \operatorname{Tr}\bigl(P_a\rho\bigr) \end{aligned}

pentru

Pa=k=0N1AkaaAk.P_a = \sum_{k = 0}^{N-1} A_k^{\dagger} \vert a\rangle\langle a \vert A_k.

Astfel, pentru aceleași matrice P0,,Pm1P_0,\ldots,P_{m-1} putem exprima canalul Φ\Phi după cum urmează.

Φ(ρ)=a=0m1Tr(Paρ)aa\Phi(\rho) = \sum_{a = 0}^{m-1} \operatorname{Tr}(P_a \rho) \vert a\rangle\langle a\vert

Această expresie este consistentă cu descrierea noastră a măsurătorilor generale în termeni de matrice, deoarece vedem că fiecare rezultat al măsurătorii apare cu probabilitatea Tr(Paρ).\operatorname{Tr}(P_a \rho).

Acum să observăm că cele două proprietăți cerute colecției de matrice {P0,,Pm1}\{P_0,\ldots,P_{m-1}\} pentru a descrie o măsurătoare generală sunt într-adevăr satisfăcute. Prima proprietate este că toate sunt matrice pozitiv semidefinite. Un mod de a vedea acest lucru este să observăm că, pentru orice vector ψ\vert \psi\rangle cu intrări în corespondență cu starea clasică a lui X,\mathsf{X}, avem

ψPaψ=k=0N1ψAkaaAkψ=k=0N1aAkψ20.\langle \psi \vert P_a \vert \psi\rangle = \sum_{k = 0}^{N-1} \langle \psi \vert A_k^{\dagger} \vert a\rangle\langle a \vert A_k\vert \psi\rangle = \sum_{k = 0}^{N-1} \bigl\vert\langle a \vert A_k\vert \psi\rangle\bigr\vert^2 \geq 0.

A doua proprietate este că dacă sumăm aceste matrice obținem matricea identitate.

a=0m1Pa=a=0m1k=0N1AkaaAk=k=0N1Ak(a=0m1aa)Ak=k=0N1AkAk=IX\begin{aligned} \sum_{a = 0}^{m-1} P_a & = \sum_{a = 0}^{m-1} \sum_{k = 0}^{N-1} A_k^{\dagger} \vert a\rangle\langle a \vert A_k \\ & = \sum_{k = 0}^{N-1} A_k^{\dagger} \Biggl(\sum_{a = 0}^{m-1} \vert a\rangle\langle a \vert\Biggr) A_k \\ & = \sum_{k = 0}^{N-1} A_k^{\dagger} A_k \\ & = \mathbb{I}_{\mathsf{X}} \end{aligned}

Ultima egalitate urmează din faptul că Φ\Phi este un canal, deci matricele sale Kraus trebuie să satisfacă această condiție.

De la matrice la canale

Acum să verificăm că pentru orice colecție {P0,,Pm1}\{P_0,\ldots,P_{m-1}\} de matrice pozitiv semidefinite care satisfac P0++Pm1=IX,P_0 + \cdots + P_{m-1} = \mathbb{I}_{\mathsf{X}}, aplicația definită prin

Φ(ρ)=a=0m1Tr(Paρ)aa\Phi(\rho) = \sum_{a = 0}^{m-1} \operatorname{Tr}(P_a \rho) \vert a \rangle\langle a\vert

este într-adevăr un canal valid de la X\mathsf{X} la Y.\mathsf{Y}.

Un mod de a face acest lucru este să calculăm reprezentarea Choi a acestei aplicații.

J(Φ)=b,c=0n1bcΦ(bc)=b,c=0n1a=0m1bcTr(Pabc)aa=b,c=0n1a=0m1bbPaTccaa=a=0m1PaTaa\begin{aligned} J(\Phi) & = \sum_{b,c = 0}^{n-1} \vert b \rangle \langle c \vert \otimes \Phi(\vert b \rangle \langle c \vert)\\[1mm] & = \sum_{b,c = 0}^{n-1} \sum_{a = 0}^{m-1} \vert b \rangle \langle c \vert \otimes \operatorname{Tr}(P_a \vert b \rangle \langle c \vert) \vert a \rangle\langle a\vert\\[1mm] & = \sum_{b,c = 0}^{n-1} \sum_{a = 0}^{m-1} \vert b \rangle \langle b \vert P_a^T \vert c \rangle \langle c \vert \otimes \vert a \rangle\langle a\vert\\[1mm] & = \sum_{a = 0}^{m-1} P_a^T \otimes \vert a \rangle\langle a\vert \end{aligned}

Transpusa fiecărui PaP_a este introdusă pentru a treia egalitate deoarece

cPab=bPaTc.\langle c \vert P_a \vert b\rangle = \langle b \vert P_a^T \vert c\rangle.

Aceasta permite apariția expresiilor bb\vert b \rangle \langle b \vert și cc,\vert c \rangle \langle c \vert, care se simplifică la matricea identitate prin sumarea peste b,b, respectiv c.c.

Prin presupunerea că P0,,Pm1P_0,\ldots,P_{m-1} sunt pozitiv semidefinite, la fel sunt și P0T,,Pm1T.P_0^{T},\ldots,P_{m-1}^{T}. În particular, transpusa unei matrice Hermitian rezultă în altă matrice Hermitian, iar valorile proprii ale oricărei matrice pătrate și ale transpusei sale coincid întotdeauna. Rezultă că J(Φ)J(\Phi) este pozitiv semidefinit. Calculând urma parțială a sistemului de ieșire Y\mathsf{Y} (care este sistemul din dreapta) obținem

TrY(J(Φ))=a=0m1PaT=IXT=IX,\operatorname{Tr}_{\mathsf{Y}} (J(\Phi)) = \sum_{a = 0}^{m-1} P_a^T = \mathbb{I}_{\mathsf{X}}^T = \mathbb{I}_{\mathsf{X}},

și astfel concluzionăm că Φ\Phi este un canal.

Măsurători parțiale

Fie mai multe sisteme care sunt colectiv într-o stare cuantică și o măsurătoare generală se efectuează pe unul dintre sisteme. Aceasta produce unul dintre rezultatele măsurătorii, selectat aleatoriu în conformitate cu probabilitățile determinate de măsurătoare și de starea sistemului înainte de măsurătoare. Starea rezultantă a sistemelor rămase va depinde, în general, de rezultatul măsurătorii obținut.

Să examinăm cum funcționează acest lucru pentru o pereche de sisteme (X,Z)(\mathsf{X},\mathsf{Z}) când sistemul X\mathsf{X} este măsurat. (Numim sistemul din dreapta Z\mathsf{Z} deoarece vom lua Y\mathsf{Y} ca sistem ce reprezintă ieșirea clasică a măsurătorii când o privim ca pe un canal.) Putem apoi generaliza cu ușurință la situația în care ordinea sistemelor este inversată, precum și la trei sau mai multe sisteme.

Fie starea lui (X,Z)(\mathsf{X},\mathsf{Z}) înainte de măsurătoare descrisă de o matrice de densitate ρ,\rho, pe care o putem scrie după cum urmează.

ρ=b,c=0n1bcρb,c\rho = \sum_{b,c = 0}^{n-1} \vert b\rangle\langle c\vert \otimes \rho_{b,c}

În această expresie presupunem că stările clasice ale lui X\mathsf{X} sunt 0,,n1.0,\ldots,n-1.

Vom presupune că măsurătoarea în sine este descrisă de colecția de matrice {P0,,Pm1}.\{P_0,\ldots,P_{m-1}\}. Această măsurătoare poate fi descrisă alternativ ca un canal Φ\Phi de la X\mathsf{X} la Y,\mathsf{Y}, unde Y\mathsf{Y} este un sistem nou cu mulțimea de stări clasice {0,,m1}.\{0,\ldots,m-1\}. Mai precis, acțiunea acestui canal poate fi exprimată după cum urmează.

Φ(ξ)=a=0m1Tr(Paξ)aa\Phi(\xi) = \sum_{a = 0}^{m-1} \operatorname{Tr}(P_a \xi)\, \vert a \rangle \langle a \vert

Probabilitățile rezultatelor

Avem în vedere o măsurătoare a sistemului X,\mathsf{X}, deci probabilitățile cu care se obțin diferitele rezultate ale măsurătorii pot depinde numai de ρX,\rho_{\mathsf{X}}, starea redusă a lui X.\mathsf{X}. În particular, probabilitatea pentru fiecare rezultat a{0,,m1}a\in\{0,\ldots,m-1\} de a apărea poate fi exprimată în trei moduri echivalente.

Tr(PaρX)=Tr(PaTrZ(ρ))=Tr((PaIZ)ρ)\operatorname{Tr}\bigl( P_a \rho_{\mathsf{X}}\bigr) = \operatorname{Tr}\bigl( P_a \operatorname{Tr}_{\mathsf{Z}}(\rho)\bigr) = \operatorname{Tr}\bigl( (P_a \otimes \mathbb{I}_{\mathsf{Z}}) \rho \bigr)

Prima expresie reprezintă în mod natural probabilitatea de a obține rezultatul aa pe baza a ceea ce știm deja despre măsurătorile unui singur sistem. Pentru a obține a doua expresie folosim pur și simplu definiția ρX=TrZ(ρ).\rho_{\mathsf{X}} = \operatorname{Tr}_{\mathsf{Z}}(\rho).

Pentru a obține a treia expresie este nevoie de mai multă gândire — și cei care studiază sunt încurajați să se convingă că este adevărată. Iată un indiciu: Echivalența dintre a doua și a treia expresie din ecuația anterioară nu depinde de faptul că ρ\rho este o matrice de densitate sau că fiecare PaP_a este pozitiv semidefinit. Încearcă să o demonstrezi mai întâi pentru produse tensoriale de forma ρ=MN\rho = M\otimes N și apoi concluzionează că trebuie să fie adevărat în general prin liniaritate.

Deși echivalența primei și a treia expresii din ecuația anterioară poate să nu fie imediată, ea are sens. Pornind de la o măsurătoare pe X,\mathsf{X}, definim efectiv o măsurătoare a lui (X,Z),(\mathsf{X},\mathsf{Z}), în care pur și simplu eliminăm Z\mathsf{Z} și măsurăm X.\mathsf{X}. Ca orice măsurătoare, această nouă măsurătoare poate fi descrisă de o colecție de matrice și nu este surprinzător că această măsurătoare este descrisă de colecția

{P0IZ,,Pm1IZ}.\{P_0\otimes\mathbb{I}_{\mathsf{Z}}, \ldots, P_{m-1}\otimes\mathbb{I}_{\mathsf{Z}}\}.

Stări condiționate de rezultatele măsurătorii

Dacă vrem să determinăm nu doar probabilitățile pentru diferitele rezultate, ci și starea rezultantă a lui Z\mathsf{Z} condiționată de fiecare rezultat al măsurătorii, putem apela la descrierea prin canal a măsurătorii. În particular, să examinăm starea pe care o obținem când aplicăm Φ\Phi lui X\mathsf{X} și nu facem nimic lui Z.\mathsf{Z}.

(ΦIdZ)(ρ)=b,c=0n1Φ(bc)ρb,c=a=0m1b,c=0n1Tr(Pabc)aaρb,c=a=0m1aab,c=0n1Tr(Pabc)ρb,c=a=0m1aab,c=0n1TrX((PaIZ)(bcρb,c))=a=0m1aaTrX((PaIZ)ρ)\begin{aligned} (\Phi\otimes\operatorname{Id}_{\mathsf{Z}})(\rho) & = \sum_{b,c = 0}^{n-1} \Phi(\vert b\rangle\langle c\vert) \otimes \rho_{b,c}\\ & = \sum_{a = 0}^{m-1} \sum_{b,c = 0}^{n-1} \operatorname{Tr}(P_a \vert b\rangle\langle c\vert) \,\vert a\rangle \langle a \vert \otimes \rho_{b,c}\\ & = \sum_{a = 0}^{m-1} \vert a\rangle \langle a \vert \otimes \sum_{b,c = 0}^{n-1} \operatorname{Tr}(P_a \vert b\rangle\langle c\vert) \rho_{b,c}\\ & = \sum_{a = 0}^{m-1} \vert a\rangle \langle a \vert \otimes \sum_{b,c = 0}^{n-1} \operatorname{Tr}_{\mathsf{X}}\bigl((P_a\otimes\mathbb{I}_{\mathsf{Z}}) (\vert b\rangle\langle c\vert\otimes\rho_{b,c})\bigr)\\ & = \sum_{a = 0}^{m-1} \vert a\rangle \langle a \vert \otimes \operatorname{Tr}_{\mathsf{X}}\bigl((P_a \otimes \mathbb{I}_{\mathsf{Z}}) \rho\bigr) \end{aligned}

Rețineți că aceasta este o matrice de densitate datorită faptului că Φ\Phi este un canal, deci fiecare matrice TrX((PaIZ)ρ)\operatorname{Tr}_{\mathsf{X}}\bigl((P_a \otimes \mathbb{I}_{\mathsf{Z}}) \rho) este neapărat pozitiv semidefinită.

Un ultim pas transformă această expresie în una care relevă ce căutăm.

a=0m1Tr((PaIZ)ρ)aaTrX((PaIZ)ρ)Tr((PaIZ)ρ)\sum_{a = 0}^{m-1} \operatorname{Tr}\bigl((P_a \otimes \mathbb{I}_{\mathsf{Z}}) \rho)\, \vert a\rangle \langle a \vert \otimes \frac{\operatorname{Tr}_{\mathsf{X}}\bigl((P_a \otimes \mathbb{I}_{\mathsf{Z}}) \rho)}{\operatorname{Tr}\bigl((P_a \otimes \mathbb{I}_{\mathsf{Z}}) \rho)}

Acesta este un exemplu de stare clasic-cuantică,

a=0m1p(a)aaσa,\sum_{a = 0}^{m-1} p(a)\, \vert a\rangle\langle a\vert \otimes \sigma_a,

cum am văzut în lecția despre Matrice de densitate. Pentru fiecare rezultat al măsurătorii a{0,,m1},a\in\{0,\ldots,m-1\}, avem cu probabilitatea

p(a)=Tr((PaIZ)ρ)p(a) = \operatorname{Tr}\bigl((P_a \otimes \mathbb{I}_{\mathsf{Z}}) \rho)

Y\mathsf{Y} se află în starea clasică aa\vert a \rangle \langle a \vert și Z\mathsf{Z} se află în starea

σa=TrX((PaIZ)ρ)Tr((PaIZ)ρ).(2)\sigma_a = \frac{\operatorname{Tr}_{\mathsf{X}}\bigl((P_a \otimes \mathbb{I}_{\mathsf{Z}}) \rho)}{\operatorname{Tr}\bigl((P_a \otimes \mathbb{I}_{\mathsf{Z}}) \rho)}. \tag{2}

Adică, aceasta este matricea de densitate pe care o obținem normalizând

TrX((PaIZ)ρ)\operatorname{Tr}_{\mathsf{X}}\bigl((P_a \otimes \mathbb{I}_{\mathsf{Z}}) \rho)

prin împărțire la urma sa. (Formal vorbind, starea σa\sigma_a este definită numai când probabilitatea p(a)p(a) este nenulă; când p(a)=0p(a) = 0 această stare este irelevantă, deoarece se referă la un eveniment discret care apare cu probabilitate zero.)

Firesc, probabilitățile rezultatelor sunt consistente cu observațiile noastre anterioare.

Pe scurt, acesta este ceea ce se întâmplă când măsurătoarea {P0,,Pm1}\{P_0,\ldots,P_{m-1}\} este efectuată pe X\mathsf{X} când (X,Z)(\mathsf{X},\mathsf{Z}) se află în starea ρ.\rho.

  1. Fiecare rezultat aa apare cu probabilitatea p(a)=Tr((PaIZ)ρ).p(a) = \operatorname{Tr}\bigl((P_a \otimes \mathbb{I}_{\mathsf{Z}}) \rho).
  2. Condiționat de obținerea rezultatului a,a, starea lui Z\mathsf{Z} este reprezentată de matricea de densitate σa\sigma_a arătată în ecuația (2),(2), care se obține normalizând TrX((PaIZ)ρ).\operatorname{Tr}_{\mathsf{X}}\bigl((P_a \otimes \mathbb{I}_{\mathsf{Z}}) \rho).

Generalizare

Putem adapta această descriere la alte situații, cum ar fi atunci când ordinea sistemelor este inversată sau când există trei sau mai multe sisteme. Conceptual este simplu, deși formulele pot deveni greoaie.

În general, dacă avem rr sisteme X1,,Xr,\mathsf{X}_1,\ldots,\mathsf{X}_r, starea sistemului compus (X1,,Xr)(\mathsf{X}_1,\ldots,\mathsf{X}_r) este ρ,\rho, și măsurătoarea {P0,,Pm1}\{P_0,\ldots,P_{m-1}\} este efectuată pe Xk\mathsf{X}_k, se întâmplă următoarele.

  1. Fiecare rezultat aa apare cu probabilitatea

    p(a)=Tr((IX1IXk1PaIXk+1IXr)ρ).p(a) = \operatorname{Tr}\bigl((\mathbb{I}_{\mathsf{X}_1}\otimes \cdots \otimes\mathbb{I}_{\mathsf{X}_{k-1}} \otimes P_a \otimes \mathbb{I}_{\mathsf{X}_{k+1}} \otimes \cdots \otimes\mathbb{I}_{\mathsf{X}_r}) \rho\bigr).
  2. Condiționat de obținerea rezultatului a,a, starea lui (X1,,Xk1,Xk+1,,Xr)(\mathsf{X}_1,\ldots,\mathsf{X}_{k-1},\mathsf{X}_{k+1},\ldots,\mathsf{X}_r) este reprezentată de următoarea matrice de densitate.

    TrXk((IX1IXk1PaIXk+1IXr)ρ)Tr((IX1IXk1PaIXk+1IXr)ρ)\frac{\operatorname{Tr}_{\mathsf{X}_k}\bigl((\mathbb{I}_{\mathsf{X}_1}\otimes \cdots \otimes\mathbb{I}_{\mathsf{X}_{k-1}} \otimes P_a \otimes \mathbb{I}_{\mathsf{X}_{k+1}} \otimes \cdots \otimes\mathbb{I}_{\mathsf{X}_r}) \rho\bigr)}{\operatorname{Tr}\bigl((\mathbb{I}_{\mathsf{X}_1}\otimes \cdots \otimes\mathbb{I}_{\mathsf{X}_{k-1}} \otimes P_a \otimes \mathbb{I}_{\mathsf{X}_{k+1}} \otimes \cdots \otimes\mathbb{I}_{\mathsf{X}_r}) \rho\bigr)}