Přeskočit na hlavní obsah

Vícenásobné systémy a redukované stavy

Nyni se zamirime na to, jak matice hustoty fungují pro vice systému, vcetne príkladu různých typů korelací, které mohou vyjádrit, a jak je lze použít k popisu stavů izolovaných cástí složených systémů.

Vice systémů

Matice hustoty mohou reprezentovat stavy více systémů analogickým způsobem jako stavové vektory ve zjednodušené formulaci kvantové informace, přičemž se řídí stejnou základní myšlenkou, že více systémů lze nahlížet, jako by tvořily jeden složený systém. Matematicky řečeno, řádky a sloupce matic hustoty reprezentujících stavy více systémů odpovídají kartézskému součinu množin klasických stavů jednotlivých systémů.

Připomeň si například reprezentace čtyř Bellových stavů pomocí stavových vektorů.

ϕ+=1200+1211ϕ=12001211ψ+=1201+1210ψ=12011210\begin{aligned} \vert \phi^+ \rangle & = \frac{1}{\sqrt{2}} \vert 00 \rangle + \frac{1}{\sqrt{2}} \vert 11 \rangle \\[2mm] \vert \phi^- \rangle & = \frac{1}{\sqrt{2}} \vert 00 \rangle - \frac{1}{\sqrt{2}} \vert 11 \rangle \\[2mm] \vert \psi^+ \rangle & = \frac{1}{\sqrt{2}} \vert 01 \rangle + \frac{1}{\sqrt{2}} \vert 10 \rangle \\[2mm] \vert \psi^- \rangle & = \frac{1}{\sqrt{2}} \vert 01 \rangle - \frac{1}{\sqrt{2}} \vert 10 \rangle \end{aligned}

Reprezentace těchto stavů pomocí matic hustoty jsou následující.

ϕ+ϕ+=(12001200000000120012)\vert \phi^+ \rangle \langle \phi^+ \vert = \begin{pmatrix} \frac{1}{2} & 0 & 0 & \frac{1}{2}\\[2mm] 0 & 0 & 0 & 0\\[2mm] 0 & 0 & 0 & 0\\[2mm] \frac{1}{2} & 0 & 0 & \frac{1}{2} \end{pmatrix} ϕϕ=(12001200000000120012)\vert \phi^- \rangle \langle \phi^- \vert = \begin{pmatrix} \frac{1}{2} & 0 & 0 & -\frac{1}{2}\\[2mm] 0 & 0 & 0 & 0\\[2mm] 0 & 0 & 0 & 0\\[2mm] -\frac{1}{2} & 0 & 0 & \frac{1}{2} \end{pmatrix} ψ+ψ+=(00000121200121200000)\vert \psi^+ \rangle \langle \psi^+ \vert = \begin{pmatrix} 0 & 0 & 0 & 0\\[2mm] 0 & \frac{1}{2} & \frac{1}{2} & 0\\[2mm] 0 & \frac{1}{2} & \frac{1}{2} & 0\\[2mm] 0 & 0 & 0 & 0 \end{pmatrix} ψψ=(00000121200121200000)\vert \psi^- \rangle \langle \psi^- \vert = \begin{pmatrix} 0 & 0 & 0 & 0\\[2mm] 0 & \frac{1}{2} & -\frac{1}{2} & 0\\[2mm] 0 & -\frac{1}{2} & \frac{1}{2} & 0\\[2mm] 0 & 0 & 0 & 0 \end{pmatrix}

Produktové stavy

Podobně jako u stavových vektorů, tenzorové součiny matic hustoty reprezentují nezávislost mezi stavy více systémů. Například pokud je X\mathsf{X} připraven ve stavu reprezentovaném maticí hustoty ρ\rho a Y\mathsf{Y} je nezávisle připraven ve stavu reprezentovaném σ,\sigma, pak matice hustoty popisující stav (X,Y)(\mathsf{X},\mathsf{Y}) je tenzorový součin ρσ.\rho\otimes\sigma.

Používá se zde stejná terminologie jako ve zjednodušené formulaci kvantové informace: stavy tohoto tvaru se nazývají produktové stavy.

Korelované a provázané stavy

Stavy, které nelze vyjádřit jako produktové stavy, reprezentují korelace mezi systémy. Ve skutečnosti existují různé typy korelací, které lze maticemi hustoty reprezentovat. Zde je několik příkladů.

  1. Korelované klasické stavy. Například situaci, kdy Alice a Bob sdílejí náhodný bit, můžeme vyjádřit takto:

    120000+121111=(120000000000000012)\frac{1}{2} \vert 0 \rangle \langle 0 \vert \otimes \vert 0 \rangle \langle 0 \vert + \frac{1}{2} \vert 1 \rangle \langle 1 \vert \otimes \vert 1 \rangle \langle 1 \vert = \begin{pmatrix} \frac{1}{2} & 0 & 0 & 0\\[2mm] 0 & 0 & 0 & 0\\[2mm] 0 & 0 & 0 & 0\\[2mm] 0 & 0 & 0 & \frac{1}{2} \end{pmatrix}
  2. Soubory kvantových stavů. Předpokládejme, že máme mm matic hustoty ρ0,,ρm1,\rho_0,\ldots,\rho_{m-1}, všechny reprezentující stavy systému X,\mathsf{X}, a náhodně vybereme jeden z těchto stavů podle vektoru pravděpodobnosti (p0,,pm1).(p_0,\ldots,p_{m-1}). Takový proces je reprezentován souborem (ensemble) stavů, který zahrnuje specifikaci matic hustoty ρ0,,ρm1,\rho_0,\ldots,\rho_{m-1}, stejně jako pravděpodobnosti (p0,,pm1).(p_0,\ldots,p_{m-1}). Souboru stavů můžeme přiřadit jedinou matici hustoty, popisující jak náhodný výběr k,k, tak odpovídající matici hustoty ρk,\rho_k, a to takto:

    k=0m1pkkkρk.\sum_{k = 0}^{m-1} p_k \vert k\rangle \langle k \vert \otimes \rho_k.

    Pro upřesnění, toto je stav dvojice (Y,X),(\mathsf{Y},\mathsf{X}), kde Y\mathsf{Y} reprezentuje klasický výběr kk — takže předpokládáme, že jeho množina klasických stavů je {0,,m1}.\{0,\ldots,m-1\}. Stavy tohoto tvaru se někdy nazývají klasicko-kvantové stavy.

  3. Separabilní stavy. Můžeme si představit situace, ve kterých máme klasickou korelaci mezi kvantovými stavy dvou systémů, a to takto:

    k=0m1pkρkσk.\sum_{k = 0}^{m-1} p_k \rho_k \otimes \sigma_k.

    Slovy řečeno, pro každé kk od 00 do m1m-1 platí, že s pravděpodobností pkp_k je systém vlevo ve stavu ρk\rho_k a systém vpravo ve stavu σk.\sigma_k. Stavy tohoto druhu se nazývají separabilní stavy. Tento koncept lze rozšířit i na více než dva systémy.

  4. Provázané stavy. Ne všechny stavy dvojic systémů jsou separabilní. V obecné formulaci kvantové informace se takto definuje provázanost (entanglement): stavy, které nejsou separabilní, se nazývají provázané.

    Všimni si, že tato terminologie je konzistentní s terminologií, kterou jsme používali v kurzu „Základy kvantové informace". Tam jsme řekli, že kvantové stavové vektory, které nejsou produktovými stavy, reprezentují provázané stavy — a skutečně, pro jakýkoli kvantový stavový vektor ψ,\vert\psi\rangle, který není produktovým stavem, zjistíme, že stav reprezentovaný maticí hustoty ψψ\vert\psi\rangle\langle\psi\vert není separabilní. Provázanost je pro stavy, které nejsou čisté, mnohem komplikovanější.

Redukované stavy a parciální stopa

S maticemi hustoty v kontextu více systémů můžeme udělat jednu jednoduchou, ale důležitou věc, a tou je popsat stavy, které získáme ignorováním některých systémů. Když je více systémů v kvantovém stavu a jeden nebo více z nich zahodíme nebo se rozhodneme ignorovat, stav zbývajících systémů se nazývá redukovaný stav těchto systémů. Popis redukovaných stavů pomocí matic hustoty lze snadno získat prostřednictvím zobrazení, známého jako parciální stopa, z matice hustoty popisující stav celku.

Příklad: redukované stavy pro e-bit

Předpokládejme, že máme dvojici Qubitů (A,B),(\mathsf{A},\mathsf{B}), které jsou společně ve stavu

ϕ+=1200+1211.\vert\phi^+\rangle = \frac{1}{\sqrt{2}} \vert 00 \rangle + \frac{1}{\sqrt{2}} \vert 11 \rangle.

Můžeme si představit, že Alice drží Qubit A\mathsf{A} a Bob drží B,\mathsf{B}, tedy dohromady sdílejí e-bit. Rádi bychom měli popis Alicina Qubitu A\mathsf{A} v izolaci pomocí matice hustoty, jako by se Bob rozhodl vzít svůj Qubit a vydat se ke hvězdám, aby už nikdy nebyl spatřen.

Nejprve se zamysleme nad tím, co by se stalo, kdyby se Bob někde na své cestě rozhodl změřit svůj Qubit standardním bázovým měřením. Pokud by to udělal, získal by výsledek 00 s pravděpodobností

(IA0)ϕ+2=1202=12,\bigl\| \bigl( \mathbb{I}_{\mathsf{A}} \otimes \langle 0\vert \bigr) \vert \phi^+ \rangle \bigr\|^2 = \Bigl\| \frac{1}{\sqrt{2}} \vert 0 \rangle \Bigr\|^2 = \frac{1}{2},

v takovém případě se stav Alicina Qubitu stane 0;\vert 0\rangle; a výsledek 11 by získal s pravděpodobností

(IA1)ϕ+2=1212=12,\bigl\| \bigl( \mathbb{I}_{\mathsf{A}} \otimes \langle 1\vert \bigr) \vert \phi^+ \rangle \bigr\|^2 = \Bigl\| \frac{1}{\sqrt{2}} \vert 1 \rangle \Bigr\|^2 = \frac{1}{2},

v takovém případě se stav Alicina Qubitu stane 1.\vert 1\rangle.

Takže pokud ignorujeme Bobův výsledek měření a zaměříme se na Alicin Qubit, dojdeme k závěru, že získá stav 0\vert 0\rangle s pravděpodobností 1/21/2 a stav 1\vert 1\rangle s pravděpodobností 1/2.1/2. To nás vede k popisu stavu Alicina Qubitu v izolaci pomocí matice hustoty

1200+1211=12IA.\frac{1}{2} \vert 0\rangle\langle 0\vert + \frac{1}{2} \vert 1\rangle\langle 1\vert = \frac{1}{2} \mathbb{I}_{\mathsf{A}}.

To znamená, že Alicin Qubit je v úplně smíšeném stavu. Pro upřesnění, tento popis stavu Alicina Qubitu v izolaci nezahrnuje Bobův výsledek měření; Boba zcela ignorujeme.

Nyní se může zdát, že popis matice hustoty Alicina Qubitu v izolaci, který jsme právě získali, závisí na předpokladu, že Bob změřil svůj Qubit, ale ve skutečnosti tomu tak není. To, co jsme udělali, je, že jsme využili možnost, že Bob změří svůj Qubit, abychom na základě toho, co jsme se již naučili, argumentovali, že úplně smíšený stav vznikne jako stav Alicina Qubitu. Samozřejmě nic neříká, že Bob musí svůj Qubit změřit — ale nic neříká, že to neudělá. A pokud je vzdálen na světelné roky, pak nic, co udělá nebo neudělá, nemůže nijak ovlivnit stav Alicina Qubitu při pohledu v izolaci. Jinými slovy, popis, který jsme získali pro stav Alicina Qubitu, je jediný popis konzistentní s nemožností komunikace rychlejší než světlo.

Můžeme také uvažovat stav Bobova Qubitu B,\mathsf{B}, který je rovněž úplně smíšeným stavem. Skutečně, pro všechny čtyři Bellovy stavy zjistíme, že redukovaný stav jak Alicina Qubitu, tak Bobova Qubitu je úplně smíšený stav.

Redukované stavy pro obecný kvantový stavový vektor

Nyní zobecníme právě probíraný příklad na dva libovolné systémy A\mathsf{A} a B,\mathsf{B}, nemusí to být nutně Qubity ve stavu ϕ+.\vert \phi^+\rangle. Budeme předpokládat, že množiny klasických stavů systémů A\mathsf{A} a B\mathsf{B} jsou Σ\Sigma a Γ,\Gamma, v daném pořadí. Matice hustoty ρ\rho reprezentující stav kombinovaného systému (A,B)(\mathsf{A},\mathsf{B}) má tedy řádkové a sloupcové indexy odpovídající kartézskému součinu Σ×Γ.\Sigma\times\Gamma.

Předpokládejme, že stav (A,B)(\mathsf{A},\mathsf{B}) je popsán kvantovým stavovým vektorem ψ,\vert\psi\rangle, takže matice hustoty popisující tento stav je ρ=ψψ.\rho = \vert\psi\rangle\langle\psi\vert. Získáme popis matice hustoty stavu systému A\mathsf{A} v izolaci, který se konvenčně značí ρA.\rho_{\mathsf{A}}. (Místo dolního indexu se někdy používá i horní index.)

Stavový vektor ψ\vert\psi\rangle lze vyjádřit ve tvaru

ψ=bΓϕbb\vert\psi\rangle = \sum_{b\in\Gamma} \vert\phi_b\rangle \otimes \vert b\rangle

pro jednoznačně určenou kolekci vektorů {ϕb:bΓ}.\{\vert\phi_b\rangle : b\in\Gamma\}. Tyto vektory lze zejména určit pomocí jednoduchého vzorce.

ϕb=(IAb)ψ\vert\phi_b\rangle = \bigl(\mathbb{I}_{\mathsf{A}} \otimes \langle b\vert\bigr)\vert\psi\rangle

Podobným uvažováním jako v předchozím příkladu s e-bitem, kdybychom měřili systém B\mathsf{B} měřením ve standardní bázi, získali bychom každý výsledek bΓb\in\Gamma s pravděpodobností ϕb2,\|\vert\phi_b\rangle\|^2, přičemž stav systému A\mathsf{A} by se stal

ϕbϕb.\frac{\vert \phi_b \rangle}{\|\vert\phi_b\rangle\|}.

Jako matici hustoty lze tento stav zapsat následovně.

(ϕbϕb)(ϕbϕb)=ϕbϕbϕb2\biggl(\frac{\vert \phi_b \rangle}{\|\vert\phi_b\rangle\|}\biggr) \biggl(\frac{\vert \phi_b \rangle}{\|\vert\phi_b\rangle\|}\biggr)^{\dagger} = \frac{\vert \phi_b \rangle\langle\phi_b\vert}{\|\vert\phi_b\rangle\|^2}

Zprůměrováním různých stavů podle pravděpodobností příslušných výsledků dospějeme k matici hustoty

ρA=bΓϕb2ϕbϕbϕb2=bΓϕbϕb=bΓ(IAb)ψψ(IAb)\rho_{\mathsf{A}} = \sum_{b\in\Gamma} \|\vert\phi_b\rangle\|^2 \frac{\vert \phi_b \rangle\langle\phi_b\vert}{\|\vert\phi_b\rangle\|^2} = \sum_{b\in\Gamma} \vert \phi_b \rangle\langle\phi_b\vert = \sum_{b\in\Gamma} \bigl(\mathbb{I}_{\mathsf{A}} \otimes \langle b\vert\bigr) \vert\psi\rangle\langle\psi\vert \bigl(\mathbb{I}_{\mathsf{A}} \otimes \vert b\rangle\bigr)

Parciální stopa

Vzorec

ρA=bΓ(IAb)ψψ(IAb)\rho_{\mathsf{A}} = \sum_{b\in\Gamma} \bigl(\mathbb{I}_{\mathsf{A}} \otimes \langle b\vert\bigr) \vert\psi\rangle\langle\psi\vert \bigl(\mathbb{I}_{\mathsf{A}} \otimes \vert b\rangle\bigr)

nás vede k popisu redukovaného stavu systému A\mathsf{A} pro libovolnou matici hustoty ρ\rho páru (A,B),(\mathsf{A},\mathsf{B}), nejen pro čistý stav.

ρA=bΓ(IAb)ρ(IAb)\rho_{\mathsf{A}} = \sum_{b\in\Gamma} \bigl( \mathbb{I}_{\mathsf{A}} \otimes \langle b \vert\bigr) \rho \bigl( \mathbb{I}_{\mathsf{A}} \otimes \vert b \rangle\bigr)

Tento vzorec musí fungovat, jednoduše díky linearitě spolu s faktem, že každou matici hustoty lze zapsat jako konvexní kombinaci čistých stavů.

Operace prováděná na ρ\rho k získání ρA\rho_{\mathsf{A}} v této rovnici je známá jako parciální stopa, a přesněji řečeno, parciální stopa se provádí přes B,\mathsf{B}, neboli B\mathsf{B} se vystopuje. Tato operace se značí TrB,\operatorname{Tr}_{\mathsf{B}}, takže můžeme psát

TrB(ρ)=bΓ(IAb)ρ(IAb).\operatorname{Tr}_{\mathsf{B}} (\rho) = \sum_{b\in\Gamma} \bigl( \mathbb{I}_{\mathsf{A}} \otimes \langle b \vert\bigr) \rho \bigl( \mathbb{I}_{\mathsf{A}} \otimes \vert b \rangle\bigr).

Můžeme také definovat parciální stopu přes A,\mathsf{A}, takže je to systém A,\mathsf{A}, který se vystopuje místo B,\mathsf{B}, takto.

TrA(ρ)=aΣ(aIB)ρ(aIB)\operatorname{Tr}_{\mathsf{A}} (\rho) = \sum_{a\in\Sigma} \bigl(\langle a \vert\otimes\mathbb{I}_{\mathsf{B}}\bigr) \rho \bigl(\vert a \rangle\otimes\mathbb{I}_{\mathsf{B}}\bigr)

To nám dává popis matice hustoty ρB\rho_{\mathsf{B}} stavu systému B\mathsf{B} v izolaci místo systému A.\mathsf{A}.

Abychom shrnuli, pokud (A,B)(\mathsf{A},\mathsf{B}) je libovolný pár systémů a máme matici hustoty ρ\rho popisující stav (A,B),(\mathsf{A},\mathsf{B}), redukované stavy systémů A\mathsf{A} a B\mathsf{B} jsou následující.

ρA=TrB(ρ)=bΓ(IAb)ρ(IAb)ρB=TrA(ρ)=aΣ(aIB)ρ(aIB)\begin{aligned} \rho_{\mathsf{A}} & = \operatorname{Tr}_{\mathsf{B}}(\rho) = \sum_{b\in\Gamma} \bigl( \mathbb{I}_{\mathsf{A}} \otimes \langle b \vert\bigr) \rho \bigl( \mathbb{I}_{\mathsf{A}} \otimes \vert b \rangle\bigr)\\[2mm] \rho_{\mathsf{B}} & = \operatorname{Tr}_{\mathsf{A}}(\rho) = \sum_{a\in\Sigma} \bigl( \langle a \vert \otimes \mathbb{I}_{\mathsf{B}}\bigr) \rho \bigl( \vert a \rangle\otimes \mathbb{I}_{\mathsf{B}} \bigr) \end{aligned}

Pokud ρ\rho je matice hustoty, pak ρA\rho_{\mathsf{A}} a ρB\rho_{\mathsf{B}} budou také nutně matice hustoty.

Tyto pojmy lze přirozeným způsobem zobecnit na libovolný počet systémů místo dvou. Obecně můžeme do dolního indexu matice hustoty ρ\rho vložit názvy libovolných systémů, které si zvolíme, a popsat tak redukovaný stav právě těchto systémů. Například pokud A,\mathsf{A}, B\mathsf{B} a C\mathsf{C} jsou systémy a ρ\rho je matice hustoty popisující stav (A,B,C),(\mathsf{A},\mathsf{B},\mathsf{C}), pak můžeme definovat

ρAC=TrB(ρ)=bΓ(IAbIC)ρ(IAbIC)ρC=TrAB(ρ)=aΣbΓ(abIC)ρ(abIC)\begin{aligned} \rho_{\mathsf{AC}} & = \operatorname{Tr}_{\mathsf{B}}(\rho) = \sum_{b\in\Gamma} \bigl( \mathbb{I}_{\mathsf{A}} \otimes \langle b \vert \otimes \mathbb{I}_{\mathsf{C}} \bigr) \rho \bigl( \mathbb{I}_{\mathsf{A}} \otimes \vert b \rangle \otimes \mathbb{I}_{\mathsf{C}} \bigr) \\[2mm] \rho_{\mathsf{C}} & = \operatorname{Tr}_{\mathsf{AB}}(\rho) = \sum_{a\in\Sigma} \sum_{b\in\Gamma} \bigl( \langle a \vert \otimes \langle b \vert \otimes \mathbb{I}_{\mathsf{C}} \bigr) \rho \bigl( \vert a \rangle \otimes \vert b \rangle \otimes \mathbb{I}_{\mathsf{C}} \bigr) \end{aligned}

a analogicky pro další volby systémů.

Alternativní popis parciální stopy

Alternativní způsob, jak popsat zobrazení parciální stopy TrA\operatorname{Tr}_{\mathsf{A}} a TrB,\operatorname{Tr}_{\mathsf{B}}, je, že jsou to jedinečná lineární zobrazení splňující vzorce

TrA(MN)=Tr(M)NTrB(MN)=Tr(N)M.\begin{aligned} \operatorname{Tr}_{\mathsf{A}}(M \otimes N) & = \operatorname{Tr}(M) N \\[2mm] \operatorname{Tr}_{\mathsf{B}}(M \otimes N) & = \operatorname{Tr}(N) M. \end{aligned}

V těchto vzorcích jsou NN a MM čtvercové matice příslušných rozměrů: řádky a sloupce MM odpovídají klasickým stavům systému A\mathsf{A} a řádky a sloupce NN odpovídají klasickým stavům systému B.\mathsf{B}.

Tato charakterizace parciální stopy je nejen fundamentální z matematického hlediska, ale v některých situacích také umožňuje rychlé výpočty. Uvažuj například tento stav páru Qubitů (A,B).(\mathsf{A},\mathsf{B}).

ρ=120000+1211++\rho = \frac{1}{2} \vert 0\rangle\langle 0\vert \otimes \vert 0\rangle\langle 0\vert + \frac{1}{2} \vert 1\rangle\langle 1\vert \otimes \vert +\rangle\langle +\vert

Pro výpočet redukovaného stavu ρA\rho_{\mathsf{A}} například můžeme využít linearitu spolu s faktem, že 00\vert 0\rangle\langle 0\vert a ++\vert +\rangle\langle +\vert mají jednotkovou stopu.

ρA=TrB(ρ)=12Tr(00)00+12Tr(++)11=1200+1211\rho_{\mathsf{A}} = \operatorname{Tr}_{\mathsf{B}}(\rho) = \frac{1}{2} \operatorname{Tr}\bigl(\vert 0\rangle\langle 0\vert\bigr)\, \vert 0\rangle\langle 0\vert + \frac{1}{2} \operatorname{Tr}\bigl(\vert +\rangle\langle +\vert\bigr) \vert 1\rangle\langle 1\vert = \frac{1}{2} \vert 0\rangle\langle 0\vert + \frac{1}{2} \vert 1\rangle\langle 1\vert

Redukovaný stav ρB\rho_{\mathsf{B}} lze vypočítat podobně.

ρB=TrA(ρ)=12Tr(00)00+12Tr(11)++=1200+12++\rho_{\mathsf{B}} = \operatorname{Tr}_{\mathsf{A}}(\rho) = \frac{1}{2} \operatorname{Tr}\bigl(\vert 0\rangle\langle 0\vert\bigr)\, \vert 0\rangle\langle 0\vert + \frac{1}{2} \operatorname{Tr}\bigl(\vert 1\rangle\langle 1\vert\bigr) \vert +\rangle\langle +\vert = \frac{1}{2} \vert 0\rangle\langle 0\vert + \frac{1}{2} \vert +\rangle\langle +\vert

Parciální stopa pro dva Qubity

Parciální stopu lze také popsat explicitně pomocí matic. Zde to provedeme jen pro dva Qubity, ale lze to zobecnit i na větší systémy. Předpokládejme, že máme dva Qubity (A,B),(\mathsf{A},\mathsf{B}), takže libovolnou matici hustoty popisující stav těchto dvou Qubitů lze zapsat jako

ρ=(α00α01α02α03α10α11α12α13α20α21α22α23α30α31α32α33)\rho = \begin{pmatrix} \alpha_{00} & \alpha_{01} & \alpha_{02} & \alpha_{03}\\[2mm] \alpha_{10} & \alpha_{11} & \alpha_{12} & \alpha_{13}\\[2mm] \alpha_{20} & \alpha_{21} & \alpha_{22} & \alpha_{23}\\[2mm] \alpha_{30} & \alpha_{31} & \alpha_{32} & \alpha_{33} \end{pmatrix}

pro nějakou volbu komplexních čísel {αjk:0j,k3}.\{\alpha_{jk} : 0\leq j,k\leq 3\}.

Parciální stopa přes první systém má následující vzorec.

TrA(α00α01α02α03α10α11α12α13α20α21α22α23α30α31α32α33)=(α00α01α10α11)+(α22α23α32α33)=(α00+α22α01+α23α10+α32α11+α33)\operatorname{Tr}_{\mathsf{A}} \begin{pmatrix} \alpha_{00} & \alpha_{01} & \alpha_{02} & \alpha_{03}\\[2mm] \alpha_{10} & \alpha_{11} & \alpha_{12} & \alpha_{13}\\[2mm] \alpha_{20} & \alpha_{21} & \alpha_{22} & \alpha_{23}\\[2mm] \alpha_{30} & \alpha_{31} & \alpha_{32} & \alpha_{33} \end{pmatrix} = \begin{pmatrix} \alpha_{00} & \alpha_{01} \\[2mm] \alpha_{10} & \alpha_{11} \end{pmatrix} + \begin{pmatrix} \alpha_{22} & \alpha_{23}\\[2mm] \alpha_{32} & \alpha_{33} \end{pmatrix} = \begin{pmatrix} \alpha_{00} + \alpha_{22} & \alpha_{01} + \alpha_{23}\\[2mm] \alpha_{10} + \alpha_{32} & \alpha_{11} + \alpha_{33} \end{pmatrix}

Jeden způsob, jak o tomto vzorci přemýšlet, začíná pohledem na 4×44\times 4 matice jako na 2×22\times 2 blokové matice, kde každý blok je 2×2.2\times 2. To znamená,

ρ=(M0,0M0,1M1,0M1,1)\rho = \begin{pmatrix} M_{0,0} & M_{0,1} \\[1mm] M_{1,0} & M_{1,1} \end{pmatrix}

pro

M0,0=(α00α01α10α11),M0,1=(α02α03α12α13),M1,0=(α20α21α30α31),M1,1=(α22α23α32α33).M_{0,0} = \begin{pmatrix} \alpha_{00} & \alpha_{01} \\[2mm] \alpha_{10} & \alpha_{11} \end{pmatrix}, \quad M_{0,1} = \begin{pmatrix} \alpha_{02} & \alpha_{03} \\[2mm] \alpha_{12} & \alpha_{13} \end{pmatrix}, \quad M_{1,0} = \begin{pmatrix} \alpha_{20} & \alpha_{21} \\[2mm] \alpha_{30} & \alpha_{31} \end{pmatrix}, \quad M_{1,1} = \begin{pmatrix} \alpha_{22} & \alpha_{23} \\[2mm] \alpha_{32} & \alpha_{33} \end{pmatrix}.

Pak máme

TrA(M0,0M0,1M1,0M1,1)=M0,0+M1,1.\operatorname{Tr}_{\mathsf{A}}\begin{pmatrix} M_{0,0} & M_{0,1} \\[1mm] M_{1,0} & M_{1,1} \end{pmatrix} = M_{0,0} + M_{1,1}.

Zde je vzorec pro případ, kdy se trasuje přes druhý systém místo prvního.

TrB(α00α01α02α03α10α11α12α13α20α21α22α23α30α31α32α33)=(Tr(α00α01α10α11)Tr(α02α03α12α13)Tr(α20α21α30α31)Tr(α22α23α32α33))=(α00+α11α02+α13α20+α31α22+α33)\operatorname{Tr}_{\mathsf{B}} \begin{pmatrix} \alpha_{00} & \alpha_{01} & \alpha_{02} & \alpha_{03}\\[2mm] \alpha_{10} & \alpha_{11} & \alpha_{12} & \alpha_{13}\\[2mm] \alpha_{20} & \alpha_{21} & \alpha_{22} & \alpha_{23}\\[2mm] \alpha_{30} & \alpha_{31} & \alpha_{32} & \alpha_{33} \end{pmatrix} = \begin{pmatrix} \operatorname{Tr} \begin{pmatrix} \alpha_{00} & \alpha_{01}\\[1mm] \alpha_{10} & \alpha_{11} \end{pmatrix} & \operatorname{Tr} \begin{pmatrix} \alpha_{02} & \alpha_{03}\\[1mm] \alpha_{12} & \alpha_{13} \end{pmatrix} \\[4mm] \operatorname{Tr} \begin{pmatrix} \alpha_{20} & \alpha_{21}\\[1mm] \alpha_{30} & \alpha_{31} \end{pmatrix} & \operatorname{Tr} \begin{pmatrix} \alpha_{22} & \alpha_{23}\\[1mm] \alpha_{32} & \alpha_{33} \end{pmatrix} \end{pmatrix} = \begin{pmatrix} \alpha_{00} + \alpha_{11} & \alpha_{02} + \alpha_{13}\\[2mm] \alpha_{20} + \alpha_{31} & \alpha_{22} + \alpha_{33} \end{pmatrix}

Pomocí blokových matic ve tvaru podobném předchozímu máme tento vzorec.

TrB(M0,0M0,1M1,0M1,1)=(Tr(M0,0)Tr(M0,1)Tr(M1,0)Tr(M1,1))\operatorname{Tr}_{\mathsf{B}} \begin{pmatrix} M_{0,0} & M_{0,1} \\[1mm] M_{1,0} & M_{1,1} \end{pmatrix} = \begin{pmatrix} \operatorname{Tr}(M_{0,0}) & \operatorname{Tr}(M_{0,1}) \\[1mm] \operatorname{Tr}(M_{1,0}) & \operatorname{Tr}(M_{1,1}) \end{pmatrix}

Popis těchto funkcí pomocí blokových matic lze přirozeným a přímým způsobem rozšířit na systémy větší než Qubit​y.

Na závěr lekce aplikujme tyto vzorce na tentýž stav, který jsme uvažovali výše.

ρ=120000+1211++=(120000000001414001414).\rho = \frac{1}{2} \vert 0\rangle \langle 0 \vert \otimes \vert 0\rangle \langle 0 \vert + \frac{1}{2} \vert 1\rangle \langle 1 \vert \otimes \vert +\rangle \langle + \vert = \begin{pmatrix} \frac{1}{2} & 0 & 0 & 0\\[2mm] 0 & 0 & 0 & 0 \\[2mm] 0 & 0 & \frac{1}{4} & \frac{1}{4}\\[2mm] 0 & 0 & \frac{1}{4} & \frac{1}{4} \end{pmatrix}.

Redukovaný stav prvního systému A\mathsf{A} je

TrB(120000000001414001414)=(Tr(12000)Tr(0000)Tr(0000)Tr(14141414))=(120012)\operatorname{Tr}_{\mathsf{B}} \begin{pmatrix} \frac{1}{2} & 0 & 0 & 0\\[2mm] 0 & 0 & 0 & 0\\[2mm] 0 & 0 & \frac{1}{4} & \frac{1}{4}\\[2mm] 0 & 0 & \frac{1}{4} & \frac{1}{4} \end{pmatrix} = \begin{pmatrix} \operatorname{Tr} \begin{pmatrix} \frac{1}{2} & 0\\[1mm] 0 & 0 \end{pmatrix} & \operatorname{Tr} \begin{pmatrix} 0 & 0\\[1mm] 0 & 0 \end{pmatrix} \\[4mm] \operatorname{Tr} \begin{pmatrix} 0 & 0\\[1mm] 0 & 0 \end{pmatrix} & \operatorname{Tr} \begin{pmatrix} \frac{1}{4} & \frac{1}{4}\\[2mm] \frac{1}{4} & \frac{1}{4} \end{pmatrix} \end{pmatrix} = \begin{pmatrix} \frac{1}{2} & 0\\[2mm] 0 & \frac{1}{2} \end{pmatrix}

a redukovaný stav druhého systému B\mathsf{B} je

TrA(120000000001414001414)=(12000)+(14141414)=(34141414).\operatorname{Tr}_{\mathsf{A}} \begin{pmatrix} \frac{1}{2} & 0 & 0 & 0\\[2mm] 0 & 0 & 0 & 0\\[2mm] 0 & 0 & \frac{1}{4} & \frac{1}{4}\\[2mm] 0 & 0 & \frac{1}{4} & \frac{1}{4} \end{pmatrix} = \begin{pmatrix} \frac{1}{2} & 0\\[1mm] 0 & 0 \end{pmatrix} + \begin{pmatrix} \frac{1}{4} & \frac{1}{4}\\[2mm] \frac{1}{4} & \frac{1}{4} \end{pmatrix} = \begin{pmatrix} \frac{3}{4} & \frac{1}{4}\\[2mm] \frac{1}{4} & \frac{1}{4} \end{pmatrix}.