Skip to main content
Physics LibreTexts

6.2: Positive and Completely Positive Maps

  • Page ID
    56574
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

    ( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\id}{\mathrm{id}}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\kernel}{\mathrm{null}\,}\)

    \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\)

    \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\)

    \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    \( \newcommand{\vectorA}[1]{\vec{#1}}      % arrow\)

    \( \newcommand{\vectorAt}[1]{\vec{\text{#1}}}      % arrow\)

    \( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vectorC}[1]{\textbf{#1}} \)

    \( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

    \( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

    \( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

    \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)

    We considered the evolution of the density operator under a family of Kraus operators in Eq. (4.30):

    \[\rho \rightarrow \rho^{\prime}=\mathscr{E}(\rho)=\sum_{k} A_{k} \rho A_{k}^{\dagger},\tag{6.10}\]

    where \(\sum_{k} A_{k}^{\dagger} A_{k}=\mathbb{I}\) (that is, \(\mathscr{E}\) is trace-preserving). When \(\mathscr{E}\) transforms any positive operator into another positive operator, we call it a positive map. We may be tempted to conclude that all positive maps correspond to physically allowed transformations. After all, it maps density operators to density operators. Unfortunately, Nature (or Mathematics?) is not that tidy.

    Consider the transpose of the density operator \(\rho \rightarrow \rho^{T}\), which acts according to

    \[\rho=\sum_{i j} \rho_{i j}|i\rangle\langle j|\rightarrow \rho^{T}=\sum_{i j} \rho_{j i}| i\rangle\langle j|.\tag{6.11}\]

    You can verify immediately that the trace is preserved in this operation (check this!), and \(\rho^{T}\) is again a positive operator since the eigenvalues are identical to those of \(\rho\). For example, consider the qubit state \((|0\rangle+i|1\rangle) / \sqrt{2}\). The density operator and its transpose are

    \[\rho=\frac{1}{2}\left(\begin{array}{cc} 1 & -i \\ i & 1
    \end{array}\right) \quad \text { and } \quad \rho^{T}=\frac{1}{2}\left(\begin{array}{cc} 1 & i \\ -i & 1 \end{array}\right).\tag{6.12}\]

    The transpose therefore corresponds to the state \((|0\rangle-i|1\rangle) / \sqrt{2}\). Now consider that the qubit is part of an entangled state \((|00\rangle+|11\rangle) / \sqrt{2}\). The density operator is given by

    \[\rho=\frac{1}{2}(|00\rangle\langle 00|+| 00\rangle\langle 11|+| 11\rangle\langle 00|+| 11\rangle\langle 11|),\tag{6.13}\]

    and the partial transpose on the first qubit is

    \[\rho^{T}=\frac{1}{2}(|00\rangle\langle 00|+| 10\rangle\langle 01|+| 01\rangle\langle 10|+| 11\rangle\langle 11|).\tag{6.14}\]

    The eigenvalues of \(\rho\) are all positive, but \(\rho^{T}\) has a negative eigenvalue! So \(\rho^{T}\) cannot be a density operator. Consequently, it is not correct to say that positive maps correspond to physical processes. We need to put another restriction on maps.

    From the example of the partial transpose, we can deduce that maps must not only be positive for the system \(S\) that they act on, but also positive on larger systems that include \(S\) as a subsystem. When this is the case, we call the map completely positive. There is a very important theorem in mathematics, called Kraus’ Representation Theorem, which states that maps of the form in Eq. (6.10) with the restriction that \(\sum_{k} A_{k}^{\dagger} A_{k}=\mathbb{I}\) is a completely positive map, and moreover, that any completely positive map can be expressed in this form.

    Exercises

      1. Show that \(\sum_{k} A_{k}^{\dagger} A_{k}=\mathbb{I}\),
      2. prove that any non-Hermitian square matrix can be written as \(A+i B\), with \(A\) and \(B\) Hermitian,
      3. prove that \(L_{0}=-\frac{1}{2} \sum_{k \neq 0} L_{k}^{\dagger} L_{k}\).
    1. Consider a two-level system \((|0\rangle,|1\rangle)\) that has a dephasing process, modelled by the Lindblad operators \(L_{1}=\gamma|0\rangle\langle 1|\) and \(L_{2}=\gamma|1\rangle\langle 0|\).
      1. write down the Lindblad equation (choose \(H = 0\) for simplicity).
      2. Calculate the evolution of the pure states \(|0\rangle\) and \(|+\rangle\) at \(t=0\). Hint: write the density matrix in the Pauli matrix basis \(\{\mathbb{I}, X, Y, Z\}\). What can you say about the equilibrium state of the system?
      3. Calculate and plot the entropy \(S(\rho)\) of the state \(\rho(t)\) as a function of \(\gamma\) and \(t\).
    2. Calculate the eigenvalues of \(\rho^{T}\) in Eq. (6.14).

    This page titled 6.2: Positive and Completely Positive Maps is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by Pieter Kok via source content that was edited to the style and standards of the LibreTexts platform; a detailed edit history is available upon request.