Skip to main content

Spin-1/2 Systems

Spin-1/2 is the smallest nontrivial quantum laboratory. Its Hilbert space is two-dimensional, but it already contains incompatible observables, sequential measurement effects, unitary rotations, density matrices, entanglement, and the difference between amplitudes and probabilities. This is why Sakurai starts with Stern-Gerlach analyzers instead of beginning with differential equations.

Sakurai's treatment is operational and experimental: beams split, recombine, and reveal the logic of spin states. Ballentine returns to spin for measurement and state preparation, emphasizing what is inferred from repeated preparations. The Gottfried-named notes use spin-1/2 as a clean illustration of the postulates and density matrices. Schiff's older notation usually reaches spin after wave mechanics, so the contrast is mainly pedagogical.

Definitions

The standard SzS_z basis is

+z=(10),z=(01).|+z\rangle=\begin{pmatrix}1\\0\end{pmatrix}, \qquad |-z\rangle=\begin{pmatrix}0\\1\end{pmatrix}.

The Pauli matrices are

σx=(0110),σy=(0ii0),σz=(1001).\sigma_x=\begin{pmatrix}0&1\\1&0\end{pmatrix},\quad \sigma_y=\begin{pmatrix}0&-i\\i&0\end{pmatrix},\quad \sigma_z=\begin{pmatrix}1&0\\0&-1\end{pmatrix}.

Spin operators are

Si=2σi,i=x,y,z.S_i={\hbar\over 2}\sigma_i,\qquad i=x,y,z.

They satisfy

[Si,Sj]=iϵijkSk,[S_i,S_j]=i\hbar\epsilon_{ijk}S_k,

and

S2=Sx2+Sy2+Sz2=324I.S^2=S_x^2+S_y^2+S_z^2={3\hbar^2\over 4}I.

A general normalized spinor can be written, up to global phase, as

n,+=cosθ2+z+eiϕsinθ2z,|\mathbf n,+\rangle =\cos{\theta\over 2}|+z\rangle +e^{i\phi}\sin{\theta\over 2}|-z\rangle,

where n\mathbf n is a unit vector on the Bloch sphere:

n=(sinθcosϕ,sinθsinϕ,cosθ).\mathbf n=(\sin\theta\cos\phi,\sin\theta\sin\phi,\cos\theta).

The observable for spin along n\mathbf n is

Sn=2nσ.S_{\mathbf n}={\hbar\over 2}\mathbf n\cdot\boldsymbol{\sigma}.

Key results

The Pauli matrices obey

σiσj=δijI+iϵijkσk.\sigma_i\sigma_j=\delta_{ij}I+i\epsilon_{ijk}\sigma_k.

This identity makes two-level calculations compact. For any unit vector n\mathbf n,

(nσ)2=I,(\mathbf n\cdot\boldsymbol{\sigma})^2=I,

so SnS_{\mathbf n} has eigenvalues ±/2\pm\hbar/2.

The projectors onto spin up and down along n\mathbf n are

Pn,±=12(I±nσ).P_{\mathbf n,\pm}={1\over 2}(I\pm \mathbf n\cdot\boldsymbol{\sigma}).

If a state has Bloch vector r\mathbf r, with density matrix

ρ=12(I+rσ),\rho={1\over 2}(I+\mathbf r\cdot\boldsymbol{\sigma}),

then

P(+n)=Tr(ρPn,+)=1+rn2.P(+\mathbf n)=\mathrm{Tr}(\rho P_{\mathbf n,+})={1+\mathbf r\cdot \mathbf n\over 2}.

For a pure state, r=1\vert \mathbf r\vert =1. For a mixed state, 0r<10\leq \vert \mathbf r\vert \lt 1. The completely unpolarized beam has ρ=I/2\rho=I/2 and gives P(+n)=1/2P(+\mathbf n)=1/2 for every analyzer direction.

Rotations are represented by

U(Rn(θ))=exp(iθ2nσ)=Icosθ2i(nσ)sinθ2.U(R_{\mathbf n}(\theta))=\exp\left(-{i\theta\over 2}\mathbf n\cdot\boldsymbol{\sigma}\right) =I\cos{\theta\over 2}-i(\mathbf n\cdot\boldsymbol{\sigma})\sin{\theta\over 2}.

The half-angle is essential. A 2π2\pi rotation changes a spinor sign, although all single-state probabilities remain unchanged. Sakurai uses this feature to emphasize that spin is not a tiny classical spinning ball.

In a magnetic field with Hamiltonian

H=μB=γSB,H=-\boldsymbol{\mu}\cdot\mathbf B=-\gamma\mathbf S\cdot\mathbf B,

the Bloch vector precesses about B\mathbf B at the Larmor frequency ω=γB\omega=\vert \gamma\vert B up to sign convention. This connects the algebra of Pauli matrices to magnetic resonance and two-level dynamics.

Visual

DirectionSpin-up ket in SzS_z basisSpin-down ket in SzS_z basis
zz(10)\begin{pmatrix}1\\0\end{pmatrix}(01)\begin{pmatrix}0\\1\end{pmatrix}
xx12(11){1\over\sqrt2}\begin{pmatrix}1\\1\end{pmatrix}12(11){1\over\sqrt2}\begin{pmatrix}1\\-1\end{pmatrix}
yy12(1i){1\over\sqrt2}\begin{pmatrix}1\\i\end{pmatrix}12(1i){1\over\sqrt2}\begin{pmatrix}1\\-i\end{pmatrix}

Worked example 1: Sequential Stern-Gerlach analyzers

Problem. A beam is prepared in +z\vert +z\rangle. It passes through an SxS_x analyzer, and only the +x+x output is kept. That beam then passes through an SzS_z analyzer. What are the final SzS_z probabilities?

Method.

  1. Express +z\vert +z\rangle in the xx basis:
+z=12+x+12x.|+z\rangle={1\over \sqrt2}|+x\rangle+{1\over \sqrt2}|-x\rangle.
  1. The probability to pass the +x+x selection is
P(+x+z)=+x+z2=12.P(+x|+z)=|\langle +x|+z\rangle|^2={1\over 2}.
  1. Conditional on being selected, the state is now +x\vert +x\rangle.

  2. Expand +x\vert +x\rangle in the zz basis:

+x=12+z+12z.|+x\rangle={1\over \sqrt2}|+z\rangle+{1\over \sqrt2}|-z\rangle.
  1. Therefore
P(+z+x)=P(z+x)=12.P(+z|+x)=P(-z|+x)={1\over 2}.

Checked answer. The final SzS_z result is not certainly +/2+\hbar/2. The intermediate incompatible SxS_x selection changed the state. This is the central lesson of Sakurai's sequential Stern-Gerlach discussion.

Worked example 2: Spin precession in a z-directed magnetic field

Problem. A spin starts in +x\vert +x\rangle and evolves under

H=ω2σz.H={\hbar\omega\over 2}\sigma_z.

Find the probability of measuring +x+x at time tt.

Method.

  1. Write the initial state:
+x=12(11).|+x\rangle={1\over \sqrt2}\begin{pmatrix}1\\1\end{pmatrix}.
  1. The evolution operator is
U(t)=eiHt/=(eiωt/200eiωt/2).U(t)=e^{-iHt/\hbar} =\begin{pmatrix}e^{-i\omega t/2}&0\\0&e^{i\omega t/2}\end{pmatrix}.
  1. Evolve the state:
ψ(t)=12(eiωt/2eiωt/2).|\psi(t)\rangle={1\over \sqrt2}\begin{pmatrix}e^{-i\omega t/2}\\e^{i\omega t/2}\end{pmatrix}.
  1. Project onto +x\vert +x\rangle:
+xψ(t)=12(eiωt/2+eiωt/2)=cosωt2.\begin{aligned} \langle +x|\psi(t)\rangle &={1\over 2}\left(e^{-i\omega t/2}+e^{i\omega t/2}\right)\\ &=\cos{\omega t\over 2}. \end{aligned}
  1. Square the modulus:
P(+x,t)=cos2ωt2.P(+x,t)=\cos^2{\omega t\over 2}.

Checked answer. At t=0t=0, the probability is 11. At ωt=π\omega t=\pi, the spin has rotated to x\vert -x\rangle and the probability is 00.

Code

import numpy as np

sx = np.array([[0, 1], [1, 0]], dtype=complex)
sz = np.array([[1, 0], [0, -1]], dtype=complex)
plus_x = np.array([1, 1], dtype=complex) / np.sqrt(2)

omega = 1.0
for t in [0, np.pi / 2, np.pi]:
u = np.diag([np.exp(-1j * omega * t / 2), np.exp(1j * omega * t / 2)])
psi_t = u @ plus_x
prob_plus_x = abs(np.vdot(plus_x, psi_t)) ** 2
print(t, prob_plus_x, np.vdot(psi_t, sx @ psi_t).real, np.vdot(psi_t, sz @ psi_t).real)

Common pitfalls

  • Calling spin-up along one axis the same state as spin-up along another axis. +z\vert +z\rangle and +x\vert +x\rangle are different states.
  • Forgetting that Stern-Gerlach selection is state preparation, not passive observation.
  • Losing the factor of 1/21/2 in spin rotations. Spinors rotate with eiθσ/2e^{-i\theta\sigma/2}.
  • Treating the sign change under a 2π2\pi spinor rotation as directly observable for a single isolated spin state. It becomes observable through interference.
  • Using Pauli matrices without the /2\hbar/2 factor when computing physical angular momentum.
  • Confusing a mixed unpolarized beam with a coherent superposition. ρ=I/2\rho=I/2 is not the same preparation as (+z+z)/2(\vert +z\rangle+\vert -z\rangle)/\sqrt2.
  • Expecting repeated incompatible measurements to preserve earlier information. An SxS_x selection erases definite SzS_z preparation in the projective model.

Spin-1/2 calculations are often short enough that students skip the physical story, but the story is what prevents sign and basis errors. Always name the preparation axis, the measurement axis, and whether a beam is selected or recombined. A selected +x+x beam is a new preparation. A recombined pair of beams can restore interference if relative phase information has been preserved. Sakurai's sequential Stern-Gerlach examples are built around exactly this distinction.

The Bloch-sphere picture is powerful but has limits. Every pure spin-1/2 state corresponds to a point on the sphere, and every mixed spin-1/2 state corresponds to a point inside it. That picture does not mean the spin is a tiny arrow already pointing in that direction before measurement. It means the state assigns probabilities for all possible spin-axis measurements by the rule P(+n)=(1+rn)/2P(+\mathbf n)=(1+\mathbf r\cdot\mathbf n)/2. Ballentine's ensemble language is a useful guardrail: the vector summarizes preparation statistics, not necessarily hidden classical orientation.

When using Pauli matrices, keep track of whether you are computing dimensionless matrix algebra or physical angular momentum. The matrices σi\sigma_i have eigenvalues ±1\pm1, while SiS_i has eigenvalues ±/2\pm\hbar/2. Rotation operators use σi/2\sigma_i/2 in the exponent for spin-1/2, and Hamiltonians use the physical magnetic moment, which may introduce signs through charge and gyromagnetic ratio. Many wrong Larmor-precession answers are only sign-convention mistakes, but those signs matter when comparing to an experiment.

Spin also previews later themes. Two spin-1/2 systems decompose into singlet and triplet sectors, giving the simplest nontrivial example of angular-momentum addition. The singlet gives the cleanest Bell-correlation example. Mixed spin beams are the easiest density-matrix examples. For that reason, mastering this page pays off repeatedly.

Connections