QuantumCore

Lesson 17 · 60 min · core

The Bloch Sphere

Loading progress…

Every pure qubit state is, up to an unobservable global phase, a point on the unit sphere in R3\mathbb R^3. We derive the parametrisation ∣ψ⟩=cos⁡θ2∣0⟩+eiϕsin⁡θ2∣1⟩\lvert\psi\rangle=\cos\tfrac\theta2\lvert0\rangle+e^{i\phi}\sin\tfrac\theta2\lvert1\rangle from two arbitrary complex amplitudes, show that its projector is 12(I+r⃗⋅σ⃗)\tfrac12(I+\vec r\cdot\vec\sigma) with r⃗=(sin⁡θcos⁡ϕ,sin⁡θsin⁡ϕ,cos⁡θ)\vec r=(\sin\theta\cos\phi,\sin\theta\sin\phi,\cos\theta), prove the overlap formula ∣⟨ψ∣χ⟩∣2=12(1+r⃗ψ⋅r⃗χ)\lvert\langle\psi\vert\chi\rangle\rvert^2=\tfrac12(1+\vec r_\psi\cdot\vec r_\chi) (so orthogonal states are antipodal), explain the half-angles through the 2-to-1 map SU(2)→SO(3)SU(2)\to SO(3), and extend the picture to mixed states filling the Bloch ball.

Prerequisites: The Qubit, Density Operators: Pure vs Mixed

Objectives and key results
  • Map any normalised pair (α,β)(\alpha,\beta) to Bloch angles (θ,ϕ)(\theta,\phi) with θ∈[0,π]\theta\in[0,\pi], ϕ∈(−π,π]\phi\in(-\pi,\pi], and back.

  • Derive ∣ψ⟩⟨ψ∣=12(I+r⃗⋅σ⃗)\lvert\psi\rangle\langle\psi\rvert=\tfrac12(I+\vec r\cdot\vec\sigma) entry by entry and read off r⃗=(sin⁡θcos⁡ϕ,sin⁡θsin⁡ϕ,cos⁡θ)\vec r=(\sin\theta\cos\phi,\sin\theta\sin\phi,\cos\theta).

  • Prove ∣⟨ψ∣χ⟩∣2=12(1+r⃗ψ⋅r⃗χ)\lvert\langle\psi\vert\chi\rangle\rvert^2=\tfrac12(1+\vec r_\psi\cdot\vec r_\chi) and conclude that orthogonal states are antipodal points.

  • Explain why the sphere uses half-angles: UU and −U-U in SU(2)SU(2) give the same rotation of r⃗\vec r.

  • Place mixed states inside the ball: ∥r⃗∥≤1\lVert\vec r\rVert\le1 with eigenvalues 12(1±∥r⃗∥)\tfrac12(1\pm\lVert\vec r\rVert) and purity 12(1+∥r⃗∥2)\tfrac12(1+\lVert\vec r\rVert^2).

formula · Bloch parametrisation
∣ψ⟩=cos⁡θ2∣0⟩+eiϕsin⁡θ2∣1⟩,θ∈[0,π], ϕ∈(−π,π]\lvert\psi\rangle=\cos\tfrac\theta2\lvert0\rangle+e^{i\phi}\sin\tfrac\theta2\lvert1\rangle,\qquad\theta\in[0,\pi],\ \phi\in(-\pi,\pi]
identity · Projector as a Bloch vector
∣ψ⟩⟨ψ∣=12(I+r⃗⋅σ⃗),r⃗=(sin⁡θcos⁡ϕ, sin⁡θsin⁡ϕ, cos⁡θ)\lvert\psi\rangle\langle\psi\rvert=\tfrac12\big(I+\vec r\cdot\vec\sigma\big),\qquad\vec r=(\sin\theta\cos\phi,\ \sin\theta\sin\phi,\ \cos\theta)
theorem · Overlap formula
∣⟨ψ∣χ⟩∣2=12(1+r⃗ψ⋅r⃗χ),⟨ψ∣χ⟩=0  ⟺  r⃗χ=−r⃗ψ\lvert\langle\psi\vert\chi\rangle\rvert^2=\tfrac12\big(1+\vec r_\psi\cdot\vec r_\chi\big),\qquad\langle\psi\vert\chi\rangle=0\iff\vec r_\chi=-\vec r_\psi
theorem · Bloch ball
ρ=12(I+r⃗⋅σ⃗)⪰0  ⟺  ∥r⃗∥≤1,Trρ2=12(1+∥r⃗∥2)\rho=\tfrac12(I+\vec r\cdot\vec\sigma)\succeq0\iff\lVert\vec r\rVert\le1,\qquad\mathrm{Tr}\rho^2=\tfrac12\big(1+\lVert\vec r\rVert^2\big)

From two complex amplitudes to two angles#

A qubit state ∣ψ⟩=α∣0⟩+β∣1⟩\lvert\psi\rangle=\alpha\lvert0\rangle+\beta\lvert1\rangle is described by two complex numbers, that is four real parameters. Normalisation ∣α∣2+∣β∣2=1\lvert\alpha\rvert^2+\lvert\beta\rvert^2=1 removes one of them, and the measurement postulate removes another: the vectors ∣ψ⟩\lvert\psi\rangle and eiγ∣ψ⟩e^{i\gamma}\lvert\psi\rangle give identical probabilities ∣⟨m∣eiγψ⟩∣2=∣⟨m∣ψ⟩∣2\lvert\langle m\vert e^{i\gamma}\psi\rangle\rvert^2=\lvert\langle m\vert\psi\rangle\rvert^2 for every measurement, so a global phase is physically invisible. Two real parameters remain, and two real parameters are exactly what is needed to name a point on a sphere. The derivation below makes this precise.

Derivation

Derivation of ∣ψ⟩=cos⁡θ2∣0⟩+eiϕsin⁡θ2∣1⟩\lvert\psi\rangle=\cos\tfrac\theta2\lvert0\rangle+e^{i\phi}\sin\tfrac\theta2\lvert1\rangle

  1. Step 1:
    α=r0eiγ0,β=r1eiγ1,r0,r1≥0, γ0,γ1∈R\alpha=r_0e^{i\gamma_0},\quad\beta=r_1e^{i\gamma_1},\qquad r_0,r_1\ge0,\ \gamma_0,\gamma_1\in\mathbb R
    Polar form of each complex amplitude.
  2. Step 2:
    ∣ψ⟩=eiγ0(r0∣0⟩+r1ei(γ1−γ0)∣1⟩)\lvert\psi\rangle=e^{i\gamma_0}\big(r_0\lvert0\rangle+r_1e^{i(\gamma_1-\gamma_0)}\lvert1\rangle\big)
    Factor out eiγ0e^{i\gamma_0} from both terms.
  3. Step 3:
    ∣ψ⟩≃r0∣0⟩+r1eiϕ∣1⟩,ϕ:=γ1−γ0\lvert\psi\rangle\simeq r_0\lvert0\rangle+r_1e^{i\phi}\lvert1\rangle,\qquad\phi:=\gamma_1-\gamma_0
    Discard the global phase eiγ0e^{i\gamma_0}, which changes no measurement probability; only the relative phase ϕ\phi survives. Reduce ϕ\phi modulo 2π2\pi into (−π,π](-\pi,\pi].
  4. Step 4:
    r02+r12=∣α∣2+∣β∣2=1r_0^2+r_1^2=\lvert\alpha\rvert^2+\lvert\beta\rvert^2=1
    Normalisation, using ∣eiγ∣=1\lvert e^{i\gamma}\rvert=1.
  5. Step 5:
    r0=cos⁡θ2,r1=sin⁡θ2,θ2∈[0,π2]r_0=\cos\tfrac\theta2,\quad r_1=\sin\tfrac\theta2,\qquad\tfrac\theta2\in\big[0,\tfrac\pi2\big]
    (r0,r1)(r_0,r_1) is a point of the unit circle in the closed first quadrant, so it equals (cos⁡a,sin⁡a)(\cos a,\sin a) for a unique a∈[0,π/2]a\in[0,\pi/2]; write a=θ/2a=\theta/2, hence θ∈[0,π]\theta\in[0,\pi].
  6. Step 6:
    ∣ψ⟩≃cos⁡θ2∣0⟩+eiϕsin⁡θ2∣1⟩\lvert\psi\rangle\simeq\cos\tfrac\theta2\lvert0\rangle+e^{i\phi}\sin\tfrac\theta2\lvert1\rangle
    Substitute. At the poles θ=0\theta=0 or θ=π\theta=\pi the factor multiplying eiϕe^{i\phi} or the whole ∣0⟩\lvert0\rangle term vanishes, so ϕ\phi is irrelevant there (like longitude at the geographic poles); the convention is ϕ=0\phi=0.

Result

∣ψ⟩≃cos⁡θ2∣0⟩+eiϕsin⁡θ2∣1⟩,θ=2arccos⁡∣α∣,ϕ=arg⁡β−arg⁡α\boxed{\lvert\psi\rangle\simeq\cos\tfrac\theta2\lvert0\rangle+e^{i\phi}\sin\tfrac\theta2\lvert1\rangle,\qquad\theta=2\arccos\lvert\alpha\rvert,\quad\phi=\arg\beta-\arg\alpha}

The projector and the Bloch vector#

Global phase disappears automatically if we work with the projector ∣ψ⟩⟨ψ∣\lvert\psi\rangle\langle\psi\rvert instead of the vector: eiγ∣ψ⟩⟨ψ∣e−iγ=∣ψ⟩⟨ψ∣e^{i\gamma}\lvert\psi\rangle\langle\psi\rvert e^{-i\gamma}=\lvert\psi\rangle\langle\psi\rvert. The four matrices I,X,Y,ZI,X,Y,Z form a basis of the real vector space of 2×22\times2 Hermitian matrices, so the projector can be expanded in them. Recall I=(1001),X=(0110),Y=(0−ii0),Z=(100−1),I=\begin{pmatrix}1&0\\0&1\end{pmatrix},\quad X=\begin{pmatrix}0&1\\1&0\end{pmatrix},\quad Y=\begin{pmatrix}0&-i\\i&0\end{pmatrix},\quad Z=\begin{pmatrix}1&0\\0&-1\end{pmatrix}, and r⃗⋅σ⃗=rxX+ryY+rzZ=(rzrx−iryrx+iry−rz)\vec r\cdot\vec\sigma=r_xX+r_yY+r_zZ=\begin{pmatrix}r_z&r_x-ir_y\\r_x+ir_y&-r_z\end{pmatrix}.

Derivation

Entry-wise computation of ∣ψ⟩⟨ψ∣\lvert\psi\rangle\langle\psi\rvert

  1. Step 1:
    ∣ψ⟩⟨ψ∣=(cos⁡θ2eiϕsin⁡θ2)(cos⁡θ2e−iϕsin⁡θ2)\lvert\psi\rangle\langle\psi\rvert=\begin{pmatrix}\cos\frac\theta2\\e^{i\phi}\sin\frac\theta2\end{pmatrix}\begin{pmatrix}\cos\frac\theta2&e^{-i\phi}\sin\frac\theta2\end{pmatrix}
    The bra is the conjugate transpose of the ket, so eiϕe^{i\phi} becomes e−iϕe^{-i\phi}.
  2. Step 2:
    =(cos⁡2θ2e−iϕsin⁡θ2cos⁡θ2eiϕsin⁡θ2cos⁡θ2sin⁡2θ2)=\begin{pmatrix}\cos^2\frac\theta2&e^{-i\phi}\sin\frac\theta2\cos\frac\theta2\\e^{i\phi}\sin\frac\theta2\cos\frac\theta2&\sin^2\frac\theta2\end{pmatrix}
    Column times row: entry (r,c)(r,c) is ψrψc‾\psi_r\overline{\psi_c}.
  3. Step 3:
    cos⁡2θ2=12(1+cos⁡θ),sin⁡2θ2=12(1−cos⁡θ)\cos^2\tfrac\theta2=\tfrac12(1+\cos\theta),\qquad\sin^2\tfrac\theta2=\tfrac12(1-\cos\theta)
    Add and subtract cos⁡2θ2+sin⁡2θ2=1\cos^2\tfrac\theta2+\sin^2\tfrac\theta2=1 and cos⁡2θ2−sin⁡2θ2=cos⁡θ\cos^2\tfrac\theta2-\sin^2\tfrac\theta2=\cos\theta (double-angle formula), then halve.
  4. Step 4:
    sin⁡θ2cos⁡θ2=12sin⁡θ\sin\tfrac\theta2\cos\tfrac\theta2=\tfrac12\sin\theta
    Double-angle formula 2sin⁡θ2cos⁡θ2=sin⁡θ2\sin\tfrac\theta2\cos\tfrac\theta2=\sin\theta.
  5. Step 5:
    ∣ψ⟩⟨ψ∣=12(1+cos⁡θsin⁡θ (cos⁡ϕ−isin⁡ϕ)sin⁡θ (cos⁡ϕ+isin⁡ϕ)1−cos⁡θ)\lvert\psi\rangle\langle\psi\rvert=\tfrac12\begin{pmatrix}1+\cos\theta&\sin\theta\,(\cos\phi-i\sin\phi)\\\sin\theta\,(\cos\phi+i\sin\phi)&1-\cos\theta\end{pmatrix}
    Substitute the three identities and Euler's formula e±iϕ=cos⁡ϕ±isin⁡ϕe^{\pm i\phi}=\cos\phi\pm i\sin\phi.
  6. Step 6:
    =12(1+rzrx−iryrx+iry1−rz)=12(I+rxX+ryY+rzZ)=\tfrac12\begin{pmatrix}1+r_z&r_x-ir_y\\r_x+ir_y&1-r_z\end{pmatrix}=\tfrac12\big(I+r_xX+r_yY+r_zZ\big)
    Compare with 12(I+r⃗⋅σ⃗)\tfrac12(I+\vec r\cdot\vec\sigma) entry by entry: rz=cos⁡θr_z=\cos\theta, rx=sin⁡θcos⁡ϕr_x=\sin\theta\cos\phi, ry=sin⁡θsin⁡ϕr_y=\sin\theta\sin\phi.
  7. Step 7:
    ∥r⃗∥2=sin⁡2θ(cos⁡2ϕ+sin⁡2ϕ)+cos⁡2θ=1\lVert\vec r\rVert^2=\sin^2\theta(\cos^2\phi+\sin^2\phi)+\cos^2\theta=1
    Pythagoras twice: every pure state lies on the unit sphere.

Result

∣ψ⟩⟨ψ∣=12(I+r⃗⋅σ⃗),r⃗=(sin⁡θcos⁡ϕ, sin⁡θsin⁡ϕ, cos⁡θ),rj=⟨ψ∣σj∣ψ⟩\boxed{\lvert\psi\rangle\langle\psi\rvert=\tfrac12\big(I+\vec r\cdot\vec\sigma\big),\qquad\vec r=(\sin\theta\cos\phi,\ \sin\theta\sin\phi,\ \cos\theta),\quad r_j=\langle\psi\vert\sigma_j\vert\psi\rangle}

The last form in the result follows from Tr(σjσk)=2δjk\mathrm{Tr}(\sigma_j\sigma_k)=2\delta_{jk} and Tr σj=0\mathrm{Tr}\,\sigma_j=0: ⟨ψ∣σj∣ψ⟩=Tr(σj∣ψ⟩⟨ψ∣)=12Trσj+12∑krkTr(σjσk)=rj\langle\psi\vert\sigma_j\vert\psi\rangle=\mathrm{Tr}(\sigma_j\lvert\psi\rangle\langle\psi\rvert)=\tfrac12\mathrm{Tr}\sigma_j+\tfrac12\sum_kr_k\mathrm{Tr}(\sigma_j\sigma_k)=r_j. So the Bloch vector is the vector of Pauli expectation values — which is how the simulator computes it.

The six cardinal states. Plugging the special angles into the parametrisation gives the eigenstates of the three Pauli matrices:

stateamplitudesθ\thetaϕ\phir⃗\vec r
∣0⟩\lvert0\rangle(1,0)(1,0)0000(0,0,1)(0,0,1)
∣1⟩\lvert1\rangle(0,1)(0,1)π\pi00(0,0,−1)(0,0,-1)
∣+⟩\lvert+\rangle12(1,1)\tfrac1{\sqrt2}(1,1)π2\tfrac\pi200(1,0,0)(1,0,0)
∣−⟩\lvert-\rangle12(1,−1)\tfrac1{\sqrt2}(1,-1)π2\tfrac\pi2π\pi(−1,0,0)(-1,0,0)
∣+i⟩\lvert{+i}\rangle12(1,i)\tfrac1{\sqrt2}(1,i)π2\tfrac\pi2π2\tfrac\pi2(0,1,0)(0,1,0)
∣−i⟩\lvert{-i}\rangle12(1,−i)\tfrac1{\sqrt2}(1,-i)π2\tfrac\pi2−π2-\tfrac\pi2(0,−1,0)(0,-1,0)

Circuit · 1 qubit

Quantum circuit on 1 qubitH on q0; S on q0q0HS

Preparing ∣+i⟩\lvert{+i}\rangle with HH then SS: H∣0⟩=12(∣0⟩+∣1⟩)H\lvert0\rangle=\tfrac1{\sqrt2}(\lvert0\rangle+\lvert1\rangle) moves r⃗\vec r from +z+z to +x+x, and S=diag(1,i)S=\mathrm{diag}(1,i) multiplies the ∣1⟩\lvert1\rangle amplitude by i=eiπ/2i=e^{i\pi/2}, rotating the azimuth by π/2\pi/2 about zz to +y+y: SH∣0⟩=12(∣0⟩+i∣1⟩)=∣+i⟩SH\lvert0\rangle=\tfrac1{\sqrt2}(\lvert0\rangle+i\lvert1\rangle)=\lvert{+i}\rangle.

∣ψ⟩=12∣0⟩+12 i∣1⟩\lvert\psi\rangle = \tfrac{1}{\sqrt{2}}\lvert 0\rangle + \tfrac{1}{\sqrt{2}}\,i\lvert 1\rangle

Quick check

Which Bloch vector belongs to ∣ψ⟩=32∣0⟩+i2∣1⟩\lvert\psi\rangle=\tfrac{\sqrt3}{2}\lvert0\rangle+\tfrac{i}{2}\lvert1\rangle?

Select one answer

Overlaps: orthogonal means antipodal#

Theorem

Overlap formula

For pure qubit states with Bloch vectors r⃗ψ\vec r_\psi and r⃗χ\vec r_\chi, ∣⟨ψ∣χ⟩∣2=12(1+r⃗ψ⋅r⃗χ)=cos⁡2Θ2\lvert\langle\psi\vert\chi\rangle\rvert^2=\tfrac12(1+\vec r_\psi\cdot\vec r_\chi)=\cos^2\tfrac\Theta2, where Θ\Theta is the angle between the two vectors. In particular ⟨ψ∣χ⟩=0\langle\psi\vert\chi\rangle=0 if and only if r⃗χ=−r⃗ψ\vec r_\chi=-\vec r_\psi.

Proof (6 steps)
  1. Step 1:
    ∣⟨ψ∣χ⟩∣2=⟨ψ∣χ⟩⟨χ∣ψ⟩=Tr(∣ψ⟩⟨ψ∣ ∣χ⟩⟨χ∣)=Tr(ρσ)\lvert\langle\psi\vert\chi\rangle\rvert^2=\langle\psi\vert\chi\rangle\langle\chi\vert\psi\rangle=\mathrm{Tr}\big(\lvert\psi\rangle\langle\psi\rvert\,\lvert\chi\rangle\langle\chi\rvert\big)=\mathrm{Tr}(\rho\sigma)
    ⟨ψ∣χ⟩‾=⟨χ∣ψ⟩\overline{\langle\psi\vert\chi\rangle}=\langle\chi\vert\psi\rangle, then the trace identity ⟨a∣B∣a⟩=Tr(B∣a⟩⟨a∣)\langle a\vert B\vert a\rangle=\mathrm{Tr}(B\lvert a\rangle\langle a\rvert) with ρ=∣ψ⟩⟨ψ∣\rho=\lvert\psi\rangle\langle\psi\rvert, σ=∣χ⟩⟨χ∣\sigma=\lvert\chi\rangle\langle\chi\rvert.
  2. Step 2:
    Tr(ρσ)=14Tr[(I+∑irψ,iσi)(I+∑jrχ,jσj)]\mathrm{Tr}(\rho\sigma)=\tfrac14\mathrm{Tr}\Big[\big(I+\textstyle\sum_ir_{\psi,i}\sigma_i\big)\big(I+\sum_jr_{\chi,j}\sigma_j\big)\Big]
    Insert the Bloch forms of both projectors.
  3. Step 3:
    =14[Tr I+∑irψ,iTr σi+∑jrχ,jTr σj+∑i,jrψ,irχ,jTr(σiσj)]=\tfrac14\Big[\mathrm{Tr}\,I+\textstyle\sum_ir_{\psi,i}\mathrm{Tr}\,\sigma_i+\sum_jr_{\chi,j}\mathrm{Tr}\,\sigma_j+\sum_{i,j}r_{\psi,i}r_{\chi,j}\mathrm{Tr}(\sigma_i\sigma_j)\Big]
    Expand the product and use linearity of the trace.
  4. Step 4:
    =14[2+0+0+∑i,jrψ,irχ,j 2δij]=12(1+r⃗ψ⋅r⃗χ)=\tfrac14\Big[2+0+0+\textstyle\sum_{i,j}r_{\psi,i}r_{\chi,j}\,2\delta_{ij}\Big]=\tfrac12\big(1+\vec r_\psi\cdot\vec r_\chi\big)
    Tr I=2\mathrm{Tr}\,I=2, the Paulis are traceless, and σiσj=δijI+iϵijkσk\sigma_i\sigma_j=\delta_{ij}I+i\epsilon_{ijk}\sigma_k gives Tr(σiσj)=2δij\mathrm{Tr}(\sigma_i\sigma_j)=2\delta_{ij}.
  5. Step 5:
    12(1+cos⁡Θ)=cos⁡2Θ2\tfrac12(1+\cos\Theta)=\cos^2\tfrac\Theta2
    Both vectors have unit length, so r⃗ψ⋅r⃗χ=cos⁡Θ\vec r_\psi\cdot\vec r_\chi=\cos\Theta; then the half-angle identity.
  6. Step 6:
    ⟨ψ∣χ⟩=0  ⟺  r⃗ψ⋅r⃗χ=−1  ⟺  r⃗χ=−r⃗ψ\langle\psi\vert\chi\rangle=0\iff\vec r_\psi\cdot\vec r_\chi=-1\iff\vec r_\chi=-\vec r_\psi
    For unit vectors, Cauchy–Schwarz gives r⃗ψ⋅r⃗χ≥−1\vec r_\psi\cdot\vec r_\chi\ge-1 with equality exactly when they point in opposite directions.

∎

Quick check

What is ∣⟨0∣ψ⟩∣2\lvert\langle0\vert\psi\rangle\rvert^2 for a state with θ=2π3\theta=\tfrac{2\pi}{3}? Give a decimal.

tolerance ±0.001

Why half-angles, and the Bloch ball#

Derivation

The 2-to-1 map SU(2)→SO(3)SU(2)\to SO(3)

  1. Step 1:
    Rn^(α)=e−iα n^⋅σ⃗/2=cos⁡α2 I−isin⁡α2 n^⋅σ⃗R_{\hat n}(\alpha)=e^{-i\alpha\,\hat n\cdot\vec\sigma/2}=\cos\tfrac\alpha2\,I-i\sin\tfrac\alpha2\,\hat n\cdot\vec\sigma
    (n^⋅σ⃗)2=I(\hat n\cdot\vec\sigma)^2=I for a unit vector n^\hat n, so the power series of the exponential splits into even (cosine) and odd (sine) parts.
  2. Step 2:
    U 12(I+r⃗⋅σ⃗) U†=12(I+(Mr⃗)⋅σ⃗),Mjk=12Tr(σjUσkU†)U\,\tfrac12(I+\vec r\cdot\vec\sigma)\,U^\dagger=\tfrac12\big(I+(M\vec r)\cdot\vec\sigma\big),\qquad M_{jk}=\tfrac12\mathrm{Tr}\big(\sigma_jU\sigma_kU^\dagger\big)
    UσkU†U\sigma_kU^\dagger is Hermitian and traceless, so it expands in the Paulis with real coefficients MjkM_{jk} (found using Tr(σjσl)=2δjl\mathrm{Tr}(\sigma_j\sigma_l)=2\delta_{jl}). MM preserves ∥r⃗∥\lVert\vec r\rVert because purity is unitarily invariant, and det⁡M=1\det M=1 by continuity from U=IU=I, so M∈SO(3)M\in SO(3).
  3. Step 3:
    Rn^(α) ⟼ rotation of r⃗ by angle α about n^R_{\hat n}(\alpha)\ \longmapsto\ \text{rotation of }\vec r\text{ by angle }\alpha\text{ about }\hat n
    E.g. for Rz(α)=diag(e−iα/2,eiα/2)R_z(\alpha)=\mathrm{diag}(e^{-i\alpha/2},e^{i\alpha/2}) the coherence ρ10\rho_{10} picks up eiαe^{i\alpha}, i.e. ϕ↦ϕ+α\phi\mapsto\phi+\alpha: the Bloch angle is the full α\alpha while the matrix uses α/2\alpha/2.
  4. Step 4:
    Rn^(2π)=cos⁡π I−isin⁡π n^⋅σ⃗=−IR_{\hat n}(2\pi)=\cos\pi\,I-i\sin\pi\,\hat n\cdot\vec\sigma=-I
    A full turn of the Bloch vector gives −I-I on the Hilbert space; only a turn by 4π4\pi gives +I+I.
  5. Step 5:
    (−U) ρ (−U)†=UρU†(-U)\,\rho\,(-U)^\dagger=U\rho U^\dagger
    The signs cancel, so UU and −U-U produce the same rotation MM: the map U↦MU\mapsto M is exactly 2-to-1. The factor −1-1 is a global phase, invisible on the sphere, which is why the amplitudes carry θ/2\theta/2 while the sphere shows θ\theta.

Result

SU(2)→SO(3),  ±U↦M(2-to-1),Rn^(2π)=−I\boxed{SU(2)\to SO(3),\ \ \pm U\mapsto M\quad\text{(2-to-1)},\qquad R_{\hat n}(2\pi)=-I}
Theorem

Mixed states fill the Bloch ball

Every qubit density operator can be written uniquely as ρ=12(I+r⃗⋅σ⃗)\rho=\tfrac12(I+\vec r\cdot\vec\sigma) with r⃗∈R3\vec r\in\mathbb R^3, and such an operator is positive semidefinite if and only if ∥r⃗∥≤1\lVert\vec r\rVert\le1. Its eigenvalues are 12(1±∥r⃗∥)\tfrac12(1\pm\lVert\vec r\rVert) and its purity is Trρ2=12(1+∥r⃗∥2)\mathrm{Tr}\rho^2=\tfrac12(1+\lVert\vec r\rVert^2); pure states are exactly the surface ∥r⃗∥=1\lVert\vec r\rVert=1 and I/2I/2 is the centre.

Proof (5 steps)
  1. Step 1:
    ρ=12(I+r⃗⋅σ⃗),rj=Tr(ρσj)∈R\rho=\tfrac12\big(I+\vec r\cdot\vec\sigma\big),\qquad r_j=\mathrm{Tr}(\rho\sigma_j)\in\mathbb R
    {I,X,Y,Z}\{I,X,Y,Z\} is a real basis of 2×22\times2 Hermitian matrices; Trρ=1\mathrm{Tr}\rho=1 fixes the coefficient of II to 12\tfrac12, and Tr(σjσk)=2δjk\mathrm{Tr}(\sigma_j\sigma_k)=2\delta_{jk} extracts rjr_j.
  2. Step 2:
    (r⃗⋅σ⃗)2=∑j,krjrkσjσk=∥r⃗∥2I(\vec r\cdot\vec\sigma)^2=\textstyle\sum_{j,k}r_jr_k\sigma_j\sigma_k=\lVert\vec r\rVert^2I
    σjσk+σkσj=2δjkI\sigma_j\sigma_k+\sigma_k\sigma_j=2\delta_{jk}I, so the antisymmetric parts cancel in the symmetric sum.
  3. Step 3:
    spec(r⃗⋅σ⃗)={+∥r⃗∥,−∥r⃗∥} ⇒ λ±=12(1±∥r⃗∥)\text{spec}(\vec r\cdot\vec\sigma)=\{+\lVert\vec r\rVert,-\lVert\vec r\rVert\}\ \Rightarrow\ \lambda_\pm=\tfrac12\big(1\pm\lVert\vec r\rVert\big)
    Eigenvalues square to ∥r⃗∥2\lVert\vec r\rVert^2 and sum to Tr(r⃗⋅σ⃗)=0\mathrm{Tr}(\vec r\cdot\vec\sigma)=0 (for r⃗≠0\vec r\ne0); adding II and halving shifts them.
  4. Step 4:
    ρ⪰0  ⟺  λ−=12(1−∥r⃗∥)≥0  ⟺  ∥r⃗∥≤1\rho\succeq0\iff\lambda_-=\tfrac12\big(1-\lVert\vec r\rVert\big)\ge0\iff\lVert\vec r\rVert\le1
    A Hermitian matrix is positive semidefinite iff its smallest eigenvalue is non-negative.
  5. Step 5:
    Trρ2=λ+2+λ−2=14[(1+∥r⃗∥)2+(1−∥r⃗∥)2]=12(1+∥r⃗∥2)\mathrm{Tr}\rho^2=\lambda_+^2+\lambda_-^2=\tfrac14\big[(1+\lVert\vec r\rVert)^2+(1-\lVert\vec r\rVert)^2\big]=\tfrac12\big(1+\lVert\vec r\rVert^2\big)
    Purity is the sum of squared eigenvalues; the cross terms cancel. It equals 1 iff ∥r⃗∥=1\lVert\vec r\rVert=1 (pure) and 12\tfrac12 at r⃗=0\vec r=0 (ρ=I/2\rho=I/2).

∎

Quick check

A qubit has Bloch vector r⃗=(0.3,0,0.4)\vec r=(0.3,0,0.4). Which statement is correct?

Select one answer

References

  1. M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, 10th anniversary ed., Cambridge University Press (2010), §1.2 and §4.2.doi:10.1017/CBO9780511976667
  2. F. Bloch, Nuclear Induction, Physical Review 70, 460 (1946).doi:10.1103/PhysRev.70.460
  3. J. Preskill, Lecture Notes for Physics 219: Quantum Computation, Chapter 2.theory.caltech.edu

Hands-on labs

  • Hit the target on the Bloch sphere

    2 steps · auto-graded

    Place the state ∣+i⟩=12(∣0⟩+i∣1⟩)\lvert{+i}\rangle=\tfrac1{\sqrt2}(\lvert0\rangle+i\lvert1\rangle) on the Bloch sphere, first by choosing (θ,ϕ)(\theta,\phi) directly and then by preparing it from ∣0⟩\lvert0\rangle with at most two gates from {H,S}\{H,S\}.

    Open lab