2

Can't find any info on Stinespring dilation so I thought I could post here. If I have a qubit complete positive map $\Lambda$, that maps all inputs to the output $|0\rangle$, $\Lambda(\rho)=|0\rangle\langle0|$ how can I derive a Stinespring dilation for $\Lambda$ and the corresponding Kraus operators?

glS
  • 27,510
  • 7
  • 37
  • 125

2 Answers2

6

Stinespring dilation can be thought of as a way of representing an arbitrary completely positive trace preserving map $\Lambda$ on a system $A$ as a composition of two simpler maps: a unitary evolution $U_{AE}$ in a Hilbert space obtained by adjoining an auxiliary system $E$ followed the partial trace over $E$

$$ \Lambda(\rho_A) = \mathrm{tr}_E\left(U_{AE} (\rho_A \otimes |0\rangle_E\langle 0|) U_{AE}^\dagger\right) $$

where the environment is initially in a pure state, for example $|0\rangle$. Thus, in order for $\Lambda$ to return $|0\rangle_A\langle 0|$ we need $U_{AE}$ to swap the states of $A$ and $E$, i.e.

$$ U_{AE} = \text{SWAP}_{AE} = \sum_{ij} |i_A\rangle|j_E\rangle\langle j_A|\langle i_E|. $$

Let us try it

$$ \begin{align} \Lambda(\rho_A) &= \mathrm{tr}_E\left(\sum_{iji'j'} |i_A\rangle|j_E\rangle\langle j_A|\langle i_E| \left(\rho_A \otimes |0_E\rangle\langle 0_E|\right)|i'_A\rangle|j'_E\rangle\langle j'_A|\langle i'_E|\right) \\ &= \mathrm{tr}_E\left(\sum_{iji'j'} |i_A\rangle\langle j_A|\rho_A|i'_A\rangle \langle j'_A| \otimes|j_E\rangle\langle i_E|0_E\rangle\langle 0_E|j'_E\rangle\langle i'_E|\right) \\ &= \mathrm{tr}_E\left(\sum_{ji'} |0_A\rangle\langle j_A|\rho_A|i'_A\rangle \langle 0_A| \otimes|j_E\rangle\langle i'_E|\right) \\ &= \sum_{ji'} |0_A\rangle\langle j_A|\rho_A|i'_A\rangle \langle 0_A| \, \langle i'_E|j_E\rangle \\ &= \sum_j |0_A\rangle\langle j_A|\rho_A|j_A\rangle \langle 0_A| \\ &= |0_A\rangle\langle 0_A| \, \mathrm{tr}\rho_A \\ &= |0_A\rangle\langle 0_A|. \end{align} $$

Having found $U_{AE}$, we can compute Kraus operators from the relation $E_k = \langle k_E|U_{AE}|0_E\rangle$

$$ \begin{align} E_k &= \langle k_E|U_{AE}|0_E\rangle \\ &= \langle k_E| \left(\sum_{ij} |i_A\rangle|j_E\rangle\langle j_A|\langle i_E|\right)|0_E\rangle \\ &= \sum_{ij} |i_A\rangle|\langle j_A|\langle k_E|j_E\rangle\langle i_E|0\rangle \\ &= |0_A\rangle|\langle k_A| \end{align} $$

and since we are working with qubits

$$ E_0 = \begin{pmatrix} 1 & 0 \\ 0 & 0 \end{pmatrix} \quad E_1 = \begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix}. $$

Incidentally, the Kraus operators above look like those for the amplitude damping channel with decay probability $1$ which with hindsight was to be expected.

Adam Zalcman
  • 25,260
  • 3
  • 40
  • 95
3

(Stinespring) Given a basis $\{|u_k\rangle\}$ for the input space, you want an isometry $V$ such that $V |u_k\rangle=|0\rangle\otimes |u_k\rangle$ for some orthonormal set $\{|u_k\rangle\}$. Then you have $$\mathrm{Tr}_2\Big(V|u_k\rangle\!\langle u_k|V^\dagger\Big) = |0\rangle\!\langle0|,$$ and thus the same hold for any pure state $|\psi\rangle$ (because you can decompose in terms of $\{|u_k\rangle\}_k$), and thus for any input state $\rho$.

More generally, you don't need input and output bases to be the same, so any isometry of the form $V|u_k\rangle=|0\rangle\otimes|v_k\rangle$ for arbitrary bases $\{|u_k\rangle\}_k,\{|v_k\rangle\}_k$ for the respective bases will work the same. This is the most general form of such an isometry though. To see this you can notice that $V$ must be such that $V|\psi\rangle$ is a product state for all $|\psi\rangle$.

(Kraus) Denote with $A_a$ the Kraus operators. These are related to the Stinespring isometry by $A_a=(I\otimes \langle a|)V$, thus $A_a=|0\rangle\!\langle a|$.

glS
  • 27,510
  • 7
  • 37
  • 125