Skip to content

9.3

Quantifying entanglement

Fitting entanglement with a measuring stick: the more mixed the reduced density matrix, the deeper the entanglement — von Neumann entropy and concurrence upgrade "is it entangled?" to "how entangled is it?".

Recommended first

After this section you should be able to

  • Decide whether a two-qubit pure state is separable or entangled, and justify it via the Schmidt decomposition
  • Compute the reduced density matrix of a given pure state and evaluate its entanglement entropy
  • Use the concurrence formula C = 2|αδ − βγ| to estimate the entanglement of a two-qubit pure state quickly
  • Explain why "globally pure, locally mixed" is the essential signature of entanglement

Last section’s CNOT built the Bell state 00+112\frac{\ket{00}+\ket{11}}{\sqrt2}, and we proved it cannot be written as a product. But “cannot be written as a product” is only a qualitative verdict. Twist the parameters a little:

ψ1=00+112,ψ2=0.9900+10.99211(9.3.1)\ket{\psi_1}=\frac{\ket{00}+\ket{11}}{\sqrt2},\qquad \ket{\psi_2}=0.99\ket{00}+\sqrt{1-0.99^2}\,\ket{11}\tag{9.3.1}

Neither can be written as a product; both are entangled. Yet ψ2\ket{\psi_2} is almost 00\ket{00} — intuition says its entanglement is “much fainter”. Intuition needs to become a number: the Bell experiments of section 9.4 and the teleportation of section 9.5 both consume entanglement, and how much is consumed, and whether there is enough, must be computable.

First, lay the foundation: separable versus entangled

The general two-qubit pure state is

ψ=α00+β01+γ10+δ11(9.3.2)\ket{\psi}=\alpha\ket{00}+\beta\ket{01}+\gamma\ket{10}+\delta\ket{11}\tag{9.3.2}

If it can be written as ab\ket{a}\otimes\ket{b} it is called separable (each qubit minding its own business); otherwise it is entangled. The criterion fits on one line: a product state’s coefficient matrix (αβγδ)\begin{pmatrix}\alpha&\beta\\\gamma&\delta\end{pmatrix} has rank 1, i.e.

αδβγ=0separable(9.3.3)\alpha\delta-\beta\gamma=0 \quad\Longleftrightarrow\quad \text{separable}\tag{9.3.3}

A non-zero determinant means entanglement — and it even hints that the “amount of entanglement” should have something to do with the size of αδβγ|\alpha\delta-\beta\gamma|. We will come back to that at the end.

Where the old tools fall short: the wavefunction has no “local readout”

Chapter 3’s state vector describes the whole system. But in an entangled state, a single qubit simply has no state vector of its own — there is no a\ket{a} that correctly predicts every experiment performed on particle A alone. This is exactly why section 3.10 introduced the density operator: everything local lives in the reduced density matrix,

ρA=TrBψψ(9.3.4)\rho_A=\mathrm{Tr}_B\,\ket{\psi}\bra{\psi}\tag{9.3.4}

(Tracing out B = averaging away all of B’s possible outcomes, with their weights.) The key observation:

Quantifying “mixedness”: Schmidt decomposition and entanglement entropy

Any bipartite pure state can be brought, via the singular value decomposition, into Schmidt form:

ψ=kλkukAvkB,λk0, kλk=1(9.3.7)\ket{\psi}=\sum_k \sqrt{\lambda_k}\,\ket{u_k}_A\ket{v_k}_B, \qquad \lambda_k\ge0,\ \sum_k\lambda_k=1\tag{9.3.7}

where uk\ket{u_k} and vk\ket{v_k} are orthonormal bases on the A and B sides respectively. For two qubits, kk runs to at most 2. The payoff is immediate:

ρA=kλkukuk,ρB=kλkvkvk(9.3.8)\rho_A=\sum_k\lambda_k\ket{u_k}\bra{u_k},\qquad \rho_B=\sum_k\lambda_k\ket{v_k}\bra{v_k}\tag{9.3.8}

The two reduced density matrices have the same eigenvalues — the Schmidt coefficients λk\lambda_k. A single non-zero λ\lambda ⇔ product state; two non-zero λ\lambda‘s ⇔ entangled. Weigh the mixedness with chapter 3’s von Neumann entropy:

E(ψ)S(ρA)=TrρAlog2ρA=kλklog2λk(9.3.9)E(\psi)\equiv S(\rho_A)=-\mathrm{Tr}\,\rho_A\log_2\rho_A =-\sum_k\lambda_k\log_2\lambda_k\tag{9.3.9}

This is the entanglement entropy (base 2, measured in bits). Three benchmark points:

StateSchmidt coefficientsEE
Product state ab\ket{a}\otimes\ket{b}1, 01,\ 000
cosχ00+sinχ11\cos\chi\ket{00}+\sin\chi\ket{11}cos2χ, sin2χ\cos^2\chi,\ \sin^2\chi0<E<10<E<1
Bell state12, 12\tfrac12,\ \tfrac1211 (maximal)

A state with entanglement entropy 1 is called maximally entangled; “spending 1 ebit” means using up one such pair — the unit in which the protocols of sections 9.5 and 9.6 keep all their accounts.

The shortcut formula: concurrence

Entropy requires a partial trace and then a diagonalisation; for two qubits there is a shortcut. The criterion determinant αδβγ\alpha\delta-\beta\gamma from the opening delivers the concurrence:

C=2αδβγ,0C1(9.3.12)C=2\,|\alpha\delta-\beta\gamma|,\qquad 0\le C\le 1\tag{9.3.12}

C=0C=0 is separable, C=1C=1 maximally entangled. Its relation to the Schmidt coefficients is C=2λ1λ2C=2\sqrt{\lambda_1\lambda_2}, and it corresponds monotonically to the entropy:

E=h ⁣(1+1C22),h(p)plog2p(1p)log2(1p)(9.3.13)E=h\!\left(\frac{1+\sqrt{1-C^2}}{2}\right),\qquad h(p)\equiv-p\log_2 p-(1-p)\log_2(1-p)\tag{9.3.13}

Check it: Φ+\ket{\Phi^+} has α=δ=12\alpha=\delta=\tfrac{1}{\sqrt2}, β=γ=0\beta=\gamma=0, so C=212=1C=2\cdot\tfrac12=1; the state cosχ00+sinχ11\cos\chi\ket{00}+\sin\chi\ket{11} gives C=sin2χC=\sin 2\chi — as χ\chi turns from 0 to 45°, the entanglement climbs smoothly from 0 to 1.

To “see” these correlations, visit Lab module 08: entanglement and Bell non-locality — the spectacle of two particles whose individual readouts are perfectly random yet strictly correlated when combined is waiting there.

What comes next

We can now weigh entanglement, but a more fundamental account is still owed. Einstein’s old objection is not answered by “the entropy equals 1”: perhaps the two particles secretly agreed on all their answers before separating — the so-called correlations being nothing but two copies of the same crib sheet? This “local hidden variable” hypothesis sounds unassailable, yet it can be hauled before the tribunal of experiment. Next section: the Bell inequality, and how quantum mechanics received its loophole-free final verdict in 2015.

Section 67 of 106 · use to turn the page