Lesson 2: Back to the Stern–Gerlach experiment

Keywords : Stern–Gerlach, Spin 1/2, Pauli matrices, Successive measurements, Complex amplitudes

In Theme 1, the Stern—Gerlach experiment showed us the failure of the classical description of a magnetic moment, and led us to infer some elements of a new formalism: states represented by vectors, superpositions and complex amplitudes. We now have the postulates of Lesson 1 and the mathematical tools of Theme 2. It is time to proceed in the other direction: to start from the experimental results and the postulates, and deduce the mathematical model of the apparatus, that is, the state space and the observables representing measurements along the different axes.

We shall see that this model is almost entirely determined. The results of the simplest experiments, with one or two apparatuses, fix the dimension of the state space and the form of the eigenvectors, up to a few phase conventions. The spectral theorem then provides the matrices of the observables: these are the Pauli matrices. Once constructed, these matrices allow us to calculate the results of all the other experiments in series, in particular the zz, xx, zz sequence that had defeated classical physics, as well as measurements along an arbitrary axis. Along the way, we shall see why complex numbers are indispensable to quantum mechanics.

Throughout the lesson, we are concerned only with the internal state of the atom, responsible for its magnetic moment, leaving aside its motion in space. We also assume that this state does not evolve between two apparatuses.

1. What the experiment teaches us

Let us recall the results established in Theme 1. A Stern—Gerlach apparatus is characterised by the direction of its field gradient, which we call its axis: SGz\mathrm{SG}_z denotes an apparatus with axis zz, SGx\mathrm{SG}_x an apparatus with axis xx, and so on. At the output of an apparatus, the beam splits into two channels, denoted for example by z+z+ and z−z-. By blocking one of the channels, we prepare a beam in which all the atoms have given the same result; we shall denote by ∣z+⟩\ket{z+} the state thus prepared in the z+z+ channel, and similarly for the other channels. The experiments have established the following facts.

Important (The experimental facts)
  1. (E1). Quantisation. Whatever its axis, a Stern—Gerlach apparatus produces only two spots: the component of the magnetic moment along the axis takes only two opposite values, ±μ0\pm\mu_0.
  2. (E2). Reproducibility. Two apparatuses with the same axis placed in succession give the same result: an atom selected in the z+z+ channel of a first SGz\mathrm{SG}_z apparatus emerges with certainty in the z+z+ channel of a second SGz\mathrm{SG}_z apparatus.
  3. (E3). Perpendicular axes. An atom prepared in one of the channels of an apparatus and then analysed by an apparatus with a perpendicular axis emerges in each of the two channels with probability 1/21/2. This is the case for the pairs of axes (z,x)(z,x), (z,y)(z,y) and (x,y)(x,y).

These three facts will suffice to construct the whole model. The other results of Theme 1, such as the reappearance of the z−z- channel in the zz, xx, zz sequence, or the dependence of the probabilities on the angle between two apparatuses, will then be predictions of the model, which we shall compare with experiment.

It is convenient to encode the two outcomes of an analysis by the dimensionless numbers +1+1 and −1-1, rather than by ±μ0\pm\mu_0. The observables that we shall construct, denoted by σz\sigma_z, σx\sigma_x, σy\sigma_y, will therefore have eigenvalues ±1\pm1. We shall return to physical units at the end of the lesson.

2. The state space

2.1. The observable associated with the z-axis

According to the second postulate, analysis along zz is represented by an observable, that is, a self-adjoint operator σz\sigma_z acting on the state space H\H, whose spectrum is the set of possible outcomes. According to (E1), this spectrum is {+1,−1}\{+1,-1\}. Let E+E_+ and E−E_- denote the two eigenspaces, and P+P_+ and P−P_- the orthogonal projectors onto these subspaces. From Theme 2, two eigenvectors associated with distinct eigenvalues of a self-adjoint operator are orthogonal, and the eigenspaces of a self-adjoint operator in finite dimension span the whole space. We therefore have

P++P−=1,σz=(+1) P++(−1) P−.P_++P_-=\mathbf{1},\qquad\sigma_z=(+1)\,P_++(-1)\,P_- .

The first relation also directly expresses the experimental fact that each atom emerges in one channel or the other: the sum of the probabilities of the two outcomes, ⟨ψ∣P+∣ψ⟩+⟨ψ∣P−∣ψ⟩=⟨ψ|ψ⟩\bra\psi P_+\ket\psi+\bra\psi P_-\ket\psi=\braket\psi\psi, is equal to one for every normalised state.

Fact (E2) is automatically compatible with the postulates: according to the reduction postulate, an atom that has given the result +1+1 is, immediately after the measurement, in a state in E+E_+, and a new measurement of σz\sigma_z then gives +1+1 with probability ⟨ψ∣P+∣ψ⟩=1\bra\psi P_+\ket\psi=1.

2.2. Dimension two

What is the dimension of the eigenspaces? Nothing in the postulates determines it. If E+E_+ had dimension greater than one, it would contain several different states all giving the result +1+1 with certainty; another observable, compatible with σz\sigma_z, would then be needed to distinguish them. Stern—Gerlach experiments reveal no additional internal degree of freedom of this kind: an atom selected in the z+z+ channel has fully determined behaviour in all subsequent experiments, as we shall verify. We therefore make the simplest modelling assumption: the two eigenvalues are non-degenerate. The subspaces E±E_\pm are one-dimensional, spanned by two unit vectors ∣z+⟩\ket{z+} and ∣z−⟩\ket{z-}, and the state space has dimension two.

Important (State space)
The internal state of the atom is described by a vector in H=C2\H=\C^2. The two states ∣z±⟩\ket{z\pm}, prepared in the channels of a SGz\mathrm{SG}_z apparatus, form an orthonormal basis, which we represent by the column vectors ∣z+⟩=(10),∣z−⟩=(01).\ket{z+}=\begin{pmatrix}1\\0\end{pmatrix},\qquad\ket{z-}=\begin{pmatrix}0\\1\end{pmatrix}.

Every state can be written as

∣ψ⟩=α∣z+⟩+β∣z−⟩=(αβ),∣α∣2+∣β∣2=1.\ket\psi=\alpha\ket{z+}+\beta\ket{z-}=\begin{pmatrix}\alpha\\ \beta\end{pmatrix},\qquad|\alpha|^2+|\beta|^2=1 .
(1)

The two components α\alpha and β\beta are not spatial coordinates, nor are they the components of a classical magnetic moment: they are probability amplitudes. According to the Born rule, an analysis along zz gives +1+1 with probability ∣⟨z+|ψ⟩∣2=∣α∣2|\braket{z+}{\psi}|^2=|\alpha|^2 and −1-1 with probability ∣β∣2|\beta|^2. Recall also that two vectors differing only by a global phase represent the same physical state.

Let us note a subtlety in passing: each of the vectors ∣z+⟩\ket{z+} and ∣z−⟩\ket{z-} is defined only up to a phase. If ∣z−⟩\ket{z-} is replaced by eiγ∣z−⟩e^{i\gamma}\ket{z-}, it remains a unit vector in the one-dimensional subspace E−E_-, and the projector ∣z−⟩⟨z−∣\ket{z-}\bra{z-} does not change. We shall use this freedom shortly.

2.3. The matrix of the observable along z

The spectral decomposition immediately gives the matrix of σz\sigma_z in the basis (∣z+⟩,∣z−⟩)(\ket{z+},\ket{z-}):

σz=∣z+⟩⟨z+∣−∣z−⟩⟨z−∣=(1000)−(0001)=(100−1).\begin{aligned} \sigma_z=\ket{z+}\bra{z+}-\ket{z-}\bra{z-}=\begin{pmatrix}1&0\\0&0\end{pmatrix}-\begin{pmatrix}0&0\\0&1\end{pmatrix}=\begin{pmatrix}1&0\\0&-1\end{pmatrix}. \end{aligned}
(2)

This is a real diagonal matrix: it is equal to its conjugate transpose and is therefore self-adjoint, as it should be. Its eigenvalues ±1\pm1 can be read from the diagonal.

3. The observable along x

3.1. The eigenstates

Let us now turn to the xx-axis. Analysis along xx is represented by an observable σx\sigma_x acting on the same space C2\C^2, with eigenvalues ±1\pm1 according to (E1), and whose eigenvectors ∣x+⟩\ket{x+} and ∣x−⟩\ket{x-} form another orthonormal basis. To construct it, we need only find these eigenvectors.

Write ∣x+⟩=a∣z+⟩+b∣z−⟩\ket{x+}=a\ket{z+}+b\ket{z-}. According to (E3), an atom prepared in the state ∣x+⟩\ket{x+} and analysed along zz emerges in each channel with probability 1/21/2. The Born rule therefore gives

∣⟨z+|x+⟩∣2=∣a∣2=12,∣⟨z−|x+⟩∣2=∣b∣2=12.|\braket{z+}{x+}|^2=|a|^2=\frac12,\qquad|\braket{z-}{x+}|^2=|b|^2=\frac12 .

The two coefficients have the same modulus and differ only in their phase. Since the global phase of ∣x+⟩\ket{x+} is arbitrary, we may choose aa to be real and positive, obtaining

∣x+⟩=12(∣z+⟩+eiγ∣z−⟩)\ket{x+}=\frac{1}{\sqrt2}\Bigl(\ket{z+}+e^{i\gamma}\ket{z-}\Bigr)

for some phase γ\gamma. This is where the freedom noted above comes into play: by replacing ∣z−⟩\ket{z-} with eiγ∣z−⟩e^{i\gamma}\ket{z-}, which changes nothing about the physics of the zz-axis, we can always reduce to γ=0\gamma=0. We make this choice, which completely fixes the relative phase of the vectors ∣z+⟩\ket{z+} and ∣z−⟩\ket{z-}:

∣x+⟩=12(∣z+⟩+∣z−⟩).\ket{x+}=\frac{1}{\sqrt2}\Bigl(\ket{z+}+\ket{z-}\Bigr).

The second eigenvector ∣x−⟩=c∣z+⟩+d∣z−⟩\ket{x-}=c\ket{z+}+d\ket{z-} must likewise satisfy ∣c∣2=∣d∣2=12|c|^2=|d|^2=\frac12 according to (E3), and we may choose c=12c=\frac{1}{\sqrt2}. It must also be orthogonal to ∣x+⟩\ket{x+}, since the two outcomes x+x+ and x−x- are mutually exclusive, as expressed by fact (E2) applied to the xx-axis. The condition

⟨x+|x−⟩=12(c+d)=0\braket{x+}{x-}=\frac{1}{\sqrt2}\Bigl(c+d\Bigr)=0

imposes d=−cd=-c, whence

∣x−⟩=12(∣z+⟩−∣z−⟩).\ket{x-}=\frac{1}{\sqrt2}\Bigl(\ket{z+}-\ket{z-}\Bigr).

We recover the superposition introduced in Theme 1: the state ∣x+⟩\ket{x+} is not a mixture of atoms in the state ∣z+⟩\ket{z+} and atoms in the state ∣z−⟩\ket{z-}, but a single state, an equal-weight superposition of the two.

3.2. The matrix

The spectral decomposition σx=∣x+⟩⟨x+∣−∣x−⟩⟨x−∣\sigma_x=\ket{x+}\bra{x+}-\ket{x-}\bra{x-} gives the required matrix. The bra ⟨x±∣\bra{x\pm} is the row vector of the conjugate components, which are real here, and the two projectors are

∣x+⟩⟨x+∣=12(11)(11)=12(1111),∣x−⟩⟨x−∣=12(1−1−11).\begin{aligned} \ket{x+}\bra{x+}=\frac12\begin{pmatrix}1\\1\end{pmatrix}\begin{pmatrix}1&1\end{pmatrix}=\frac12\begin{pmatrix}1&1\\1&1\end{pmatrix}, \qquad \ket{x-}\bra{x-}=\frac12\begin{pmatrix}1&-1\\-1&1\end{pmatrix}. \end{aligned}

Their sum is readily seen to be the identity. Their difference gives

σx=(0110).\begin{aligned} \sigma_x=\begin{pmatrix}0&1\\1&0\end{pmatrix}. \end{aligned}
(3)

This matrix is real and symmetric, hence self-adjoint. Its eigenvectors can be checked directly: σx(1,1)=(1,1)\sigma_x(1,1)=(1,1) and σx(1,−1)=−(1,−1)\sigma_x(1,-1)=-(1,-1). The operator σx\sigma_x exchanges the two basis vectors: σx∣z+⟩=∣z−⟩\sigma_x\ket{z+}=\ket{z-} and σx∣z−⟩=∣z+⟩\sigma_x\ket{z-}=\ket{z+}.

This result can also be read as a change of basis, in the sense of Theme 2. The change-of-basis matrix from the basis (∣z±⟩)(\ket{z\pm}) to the basis (∣x±⟩)(\ket{x\pm}) has as its columns the components of ∣x±⟩\ket{x\pm}:

S=12(111−1),\begin{aligned} S=\frac{1}{\sqrt2}\begin{pmatrix}1&1\\1&-1\end{pmatrix}, \end{aligned}

which is unitary, like any change-of-basis matrix between orthonormal bases. In the basis (∣x±⟩)(\ket{x\pm}), the matrix of the observable σx\sigma_x is diag(1,−1)\mathrm{diag}(1,-1), and we recover σx=S diag(1,−1) S†\sigma_x=S\,\mathrm{diag}(1,-1)\,S^\dagger in the basis (∣z±⟩)(\ket{z\pm}).

4. The observable along y and complex numbers

4.1. A phase imposed by experiment

Let us proceed in the same way for the yy-axis. According to (E3) applied to the pair (z,y)(z,y), the two components of ∣y+⟩\ket{y+} in the basis (∣z±⟩)(\ket{z\pm}) have modulus 12\frac{1}{\sqrt2}, and, upon fixing the global phase, we may write

∣y+⟩=12(∣z+⟩+eiδ∣z−⟩).\ket{y+}=\frac{1}{\sqrt2}\Bigl(\ket{z+}+e^{i\delta}\ket{z-}\Bigr).

This time, we can no longer eliminate the phase δ\delta by redefining ∣z−⟩\ket{z-}: that freedom has already been used for the xx-axis. The phase δ\delta therefore has physical content, and it is the pair (x,y)(x,y) that will determine it. According to (E3), an atom prepared in the state ∣y+⟩\ket{y+} and analysed along xx emerges in the x+x+ channel with probability 1/21/2. But

⟨x+|y+⟩=12(1+eiδ),∣⟨x+|y+⟩∣2=14∣1+eiδ∣2=14(2+2cos⁡δ)=1+cos⁡δ2.\braket{x+}{y+}=\frac12\Bigl(1+e^{i\delta}\Bigr),\qquad |\braket{x+}{y+}|^2=\frac14\bigl|1+e^{i\delta}\bigr|^2=\frac14\bigl(2+2\cos\delta\bigr)=\frac{1+\cos\delta}{2}.

The condition 1+cos⁡δ2=12\frac{1+\cos\delta}{2}=\frac12 imposes cos⁡δ=0\cos\delta=0, that is, δ=±π2\delta=\pm\frac\pi2 and eiδ=±ie^{i\delta}=\pm i.

Important (Why complex amplitudes are needed)
With real amplitudes, the only freedom would be the sign: eiδ=+1e^{i\delta}=+1 or −1-1. But these two choices give ∣x+⟩\ket{x+} or ∣x−⟩\ket{x-}, which do not satisfy (E3) for the pair (x,y)(x,y). It is therefore impossible to account for all three axes xx, yy, zz with real amplitudes: the Stern—Gerlach results require the complex Hilbert space C2\C^2, not R2\R^2.

It remains to choose between δ=+π2\delta=+\frac\pi2 and δ=−π2\delta=-\frac\pi2. These two choices are related by a mirror symmetry that reverses the direction of the yy-axis: they describe a right-handed or left-handed frame (x,y,z)(x,y,z). Facts (E1) to (E3) cannot distinguish them, since they involve only pairs of axes. We choose δ=+π2\delta=+\frac\pi2, which corresponds, as we shall see, to a right-handed frame; this convention will be justified in Lesson 10 by the algebra of angular momentum. Thus

∣y+⟩=12(∣z+⟩+i∣z−⟩).\ket{y+}=\frac{1}{\sqrt2}\Bigl(\ket{z+}+i\ket{z-}\Bigr).

The second eigenvector must have components of modulus 12\frac{1}{\sqrt2} and be orthogonal to ∣y+⟩\ket{y+}. Writing ∣y−⟩=12(∣z+⟩+d∣z−⟩)\ket{y-}=\frac{1}{\sqrt2}(\ket{z+}+d\ket{z-}), the inner product is 12(1+(−i)d)\frac12(1+(-i)d), since the bra ⟨y+∣\bra{y+} has conjugate components 12(1,−i)\frac{1}{\sqrt2}(1,-i). It vanishes for d=−id=-i:

∣y−⟩=12(∣z+⟩−i∣z−⟩).\ket{y-}=\frac{1}{\sqrt2}\Bigl(\ket{z+}-i\ket{z-}\Bigr).

We can also check that ∣⟨x±|y−⟩∣2=12|\braket{x\pm}{y-}|^2=\frac12, as required by (E3).

4.2. The matrix

The projectors are calculated as before, taking care with conjugation in the bras:

∣y+⟩⟨y+∣=12(1i)(1−i)=12(1−ii1),∣y−⟩⟨y−∣=12(1i−i1).\begin{aligned} \ket{y+}\bra{y+}=\frac12\begin{pmatrix}1\\ i\end{pmatrix}\begin{pmatrix}1&-i\end{pmatrix}=\frac12\begin{pmatrix}1&-i\\ i&1\end{pmatrix}, \qquad \ket{y-}\bra{y-}=\frac12\begin{pmatrix}1&i\\-i&1\end{pmatrix}. \end{aligned}

Their sum is the identity, and their difference gives

σy=∣y+⟩⟨y+∣−∣y−⟩⟨y−∣=(0−ii0).\begin{aligned} \sigma_y=\ket{y+}\bra{y+}-\ket{y-}\bra{y-}=\begin{pmatrix}0&-i\\ i&0\end{pmatrix}. \end{aligned}
(4)

This matrix is not real. Let us check that it is self-adjoint by taking its transpose and then conjugating the coefficients:

σyT=(0i−i0),σy†=(σyT)∗=(0−ii0)=σy.\begin{aligned} \sigma_y^{\mathsf T}=\begin{pmatrix}0&i\\-i&0\end{pmatrix},\qquad \sigma_y^\dagger=\bigl(\sigma_y^{\mathsf T}\bigr)^*=\begin{pmatrix}0&-i\\ i&0\end{pmatrix}=\sigma_y . \end{aligned}

This is indeed the case, as guaranteed by its construction as a real linear combination of orthogonal projectors. Finally, its eigenvectors can be checked: σy(1,i)=(−i⋅i, i⋅1)=(1,i)\sigma_y(1,i)=(-i\cdot i,\,i\cdot1)=(1,i) and σy(1,−i)=(−1,i)=−(1,−i)\sigma_y(1,-i)=(-1,i)=-(1,-i).

5. The Pauli matrices

5.1. The result

We have constructed the three observables associated with the three axes. They are the matrices introduced by Pauli in 1927.

Important (Pauli matrices)
In the basis (∣z+⟩,∣z−⟩)(\ket{z+},\ket{z-}), analyses along the xx, yy and zz axes are represented by

σx=(0110),σy=(0−ii0),σz=(100−1),\begin{aligned} \sigma_x=\begin{pmatrix}0&1\\1&0\end{pmatrix},\qquad \sigma_y=\begin{pmatrix}0&-i\\ i&0\end{pmatrix},\qquad \sigma_z=\begin{pmatrix}1&0\\0&-1\end{pmatrix}, \end{aligned}
(5)

with eigenvectors

∣x±⟩=∣z+⟩±∣z−⟩2,∣y±⟩=∣z+⟩±i∣z−⟩2,\ket{x\pm}=\frac{\ket{z+}\pm\ket{z-}}{\sqrt2},\qquad \ket{y\pm}=\frac{\ket{z+}\pm i\ket{z-}}{\sqrt2},
(6)

for the eigenvalues ±1\pm1. This result is unique up to phase conventions.

Let us summarise what comes from experiment and what comes from the postulates. Facts (E1) to (E3) fixed the eigenvalues and the moduli of the components of the eigenvectors, as well as the relative phase ±π/2\pm\pi/2 between the components of ∣y+⟩\ket{y+}. The postulates provided the framework: a Hilbert space, self-adjoint observables and the Born rule for translating measured probabilities into inner products. Theme 2 provided the construction tool: the spectral decomposition A^=∑aa Pa\hat A=\sum_aa\,P_a. Three conventions remain: the phase of ∣z−⟩\ket{z-}, fixed by the choice ∣x+⟩=(∣z+⟩+∣z−⟩)/2\ket{x+}=(\ket{z+}+\ket{z-})/\sqrt2, the global phases of the eigenvectors, and the sign of δ\delta, which is related to the orientation of the frame. None of these conventions changes a measurable probability.

5.2. Algebraic properties

The three Pauli matrices have remarkable properties, which can be verified by direct calculation.

  • They are self-adjoint: σj†=σj\sigma_j^\dagger=\sigma_j for j=x,y,zj=x,y,z, as we have checked.
  • Their square is the identity: for example, σy2=(0−ii0)(0−ii0)=((−i)(i)00(i)(−i))=(1001).\begin{aligned} \sigma_y^2=\begin{pmatrix}0&-i\\ i&0\end{pmatrix}\begin{pmatrix}0&-i\\ i&0\end{pmatrix}=\begin{pmatrix}(-i)(i)&0\\0&(i)(-i)\end{pmatrix}=\begin{pmatrix}1&0\\0&1\end{pmatrix}. \end{aligned} This is the algebraic expression of the fact that their eigenvalues are ±1\pm1: a self-adjoint operator whose square is the identity has eigenvalues λ\lambda satisfying λ2=1\lambda^2=1.
  • They have zero trace and determinant −1-1, which again gives their eigenvalues as the roots of λ2−1=0\lambda^2-1=0.
  • The product of two different matrices is proportional to the third. For example, σxσy=(0110)(0−ii0)=(i00−i)=i σz,σyσx=(−i00i)=−i σz.\begin{aligned} \sigma_x\sigma_y=\begin{pmatrix}0&1\\1&0\end{pmatrix}\begin{pmatrix}0&-i\\ i&0\end{pmatrix}=\begin{pmatrix}i&0\\0&-i\end{pmatrix}=i\,\sigma_z, \qquad \sigma_y\sigma_x=\begin{pmatrix}-i&0\\0&i\end{pmatrix}=-i\,\sigma_z . \end{aligned} Similarly, σyσz=iσx\sigma_y\sigma_z=i\sigma_x and σzσx=iσy\sigma_z\sigma_x=i\sigma_y, by cyclic permutation of (x,y,z)(x,y,z). With the choice δ=−π/2\delta=-\pi/2, we would have obtained σxσy=−iσz\sigma_x\sigma_y=-i\sigma_z: it is in this sense that our convention corresponds to the right-handed frame.

It follows from the last property that two different Pauli matrices anticommute, σxσy=−σyσx\sigma_x\sigma_y=-\sigma_y\sigma_x, and that their commutators are

[σx,σy]=2iσz,[σy,σz]=2iσx,[σz,σx]=2iσy.[\sigma_x,\sigma_y]=2i \sigma_z,\qquad[\sigma_y,\sigma_z]=2i \sigma_x,\qquad[\sigma_z,\sigma_x]=2i \sigma_y .
(7)

None of these observables commutes with another. According to Theme 2, they therefore have no common eigenbasis; and one can check directly that they do not even have a single eigenvector in common, since the eigenvectors of each have components of modulus 12\frac{1}{\sqrt2} in the eigenbasis of each of the others. This is the mathematical expression of the incompatibility observed in Theme 1: no state gives a certain outcome for two different axes.

5.3. Physical units

The observables σj\sigma_j encode the outcomes as ±1\pm1. For the electron in the silver atom, we shall see in Lesson 10 that the magnetic moment is associated with spin, an angular momentum whose components are

Sj=ℏ2 σj,with eigenvalues ±ℏ2.S_j=\frac\hbar2\,\sigma_j, \qquad\text{with eigenvalues }\pm\frac\hbar2 .

Relations (7) then become [Sx,Sy]=iℏSz[S_x,S_y]=i\hbar S_z and its cyclic permutations: these are the commutation relations of angular momentum, and this is the origin of our convention for the sign of δ\delta. The magnetic moment is proportional to the spin, μ=γS\boldsymbol\mu=\gamma\mathbf S, with γ≃−e/me\gamma\simeq-e/m_e for the electron; its components therefore take the values ±∣γ∣ℏ/2=±eℏ/(2me)\pm|\gamma|\hbar/2=\pm e\hbar/(2m_e), and we recover the Bohr magneton measured by Stern and Gerlach. Since the sign of γ\gamma is negative, the channel of the apparatus in which the magnetic moment is positive corresponds to a negative spin: identifying the channels with the spin states requires a convention, which we shall assume has been fixed.

6. Experiments in series

6.1. The calculation rule

Now that the model has been constructed, we can calculate the result of any succession of apparatuses. We need only apply the postulates step by step. Consider an atom in the state ∣ψ⟩\ket\psi, which passes through a first apparatus where the outcome aa, with projector PaP_a, is selected, and then through a second where the outcome bb, with projector PbP_b, is observed.

  1. According to the Born rule, the atom emerges in channel aa with probability p(a)=⟨ψ∣Pa∣ψ⟩=∥Paψ∥2p(a)=\bra\psi P_a\ket\psi=\norm{P_a\psi}^2.
  2. According to the reduction postulate, the state of the atoms selected in this channel is ∣ψa⟩=Pa∣ψ⟩/∥Paψ∥\ket{\psi_a}=P_a\ket\psi/\norm{P_a\psi}.
  3. These atoms enter the second apparatus and emerge in channel bb with conditional probability p(b∣a)=∥Pbψa∥2=∥PbPaψ∥2/∥Paψ∥2p(b\mid a)=\norm{P_b\psi_a}^2=\norm{P_bP_a\psi}^2/\norm{P_a\psi}^2.

The probability of the complete path, relative to the beam entering the first apparatus, is the product of the two:

p(a,b)=p(a)p(b∣a)=∥PbPaψ∥2.p(a,b)=p(a) p(b\mid a)=\norm{P_bP_a\psi}^2 .
(8)

This formula generalises to any number of apparatuses. Conditional probabilities, relative to the beam entering a given apparatus, must be carefully distinguished from probabilities relative to the initial beam: the former are measured by counting the impacts at the output of an apparatus, the latter by taking into account all the atoms lost at the selection stages.

6.2. A single axis: ZZ

Let us prepare a beam in the state ∣z+⟩\ket{z+} and analyse it along zz. The amplitudes ⟨z+|z+⟩=1\braket{z+}{z+}=1 and ⟨z−|z+⟩=0\braket{z-}{z+}=0 give p(z+)=1p(z+)=1 and p(z−)=0p(z-)=0. We recover fact (E2), which was one of the inputs to the construction.

6.3. Two perpendicular axes: ZX

Let us now analyse the state ∣z+⟩\ket{z+} along xx. We can decompose ∣z+⟩\ket{z+} in the eigenbasis of σx\sigma_x: adding the expressions for ∣x+⟩\ket{x+} and ∣x−⟩\ket{x-} gives

∣z+⟩=12(∣x+⟩+∣x−⟩),\ket{z+}=\frac{1}{\sqrt2}\Bigl(\ket{x+}+\ket{x-}\Bigr),

whence p(x+∣z+)=p(x−∣z+)=12p(x+\mid z+)=p(x-\mid z+)=\frac12. Again we recover one of the inputs to the construction. Note that the expectation value of σx\sigma_x in the state ∣z+⟩\ket{z+} is zero: this can be seen from the probabilities, or directly, ⟨z+∣σx∣z+⟩=⟨z+∣(σx∣z+⟩)=⟨z+|z−⟩=0\bra{z+}\sigma_x\ket{z+}=\bra{z+}\bigl(\sigma_x\ket{z+}\bigr)=\braket{z+}{z-}=0.

6.4. Three apparatuses: ZXZ

Here is now a genuine prediction. Starting from the state ∣z+⟩\ket{z+}, let us select the x+x+ channel of a SGx\mathrm{SG}_x apparatus, and then analyse along zz again. After selection, the state is ∣x+⟩=12(∣z+⟩+∣z−⟩)\ket{x+}=\frac{1}{\sqrt2}(\ket{z+}+\ket{z-}), and the last apparatus gives

p(z+∣x+)=∣⟨z+|x+⟩∣2=12,p(z−∣x+)=∣⟨z−|x+⟩∣2=12.p(z+\mid x+)=|\braket{z+}{x+}|^2=\frac12,\qquad p(z-\mid x+)=|\braket{z-}{x+}|^2=\frac12 .

The z−z- channel, which was empty in the ZZ experiment, reappears: this is exactly what was observed in Theme 1, and what had defeated the picture of a simple classical sorting process. Relative to the initial beam, the probabilities of the two complete paths are, according to (8),

p(x+,z±)=∥Pz±Px+∣z+⟩∥2=∣⟨z±|x+⟩⟨x+|z+⟩∣2=12⋅12=14.p(x+,z\pm)=\norm{P_{z\pm}P_{x+}\ket{z+}}^2=\Bigl|\braket{z\pm}{x+}\braket{x+}{z+}\Bigr|^2=\frac12\cdot\frac12=\frac14 .

Each final output channel contains one quarter of the initial atoms, the remaining half having been stopped by the x−x- filter.

The model accounts for this result without any additional assumption: the x+x+ filter does not merely sort atoms that already have a value of μx\mu_x and a value of μz\mu_z; it prepares a new state, ∣x+⟩\ket{x+}, which retains no memory of the z+z+ outcome obtained previously. This is a consequence of the reduction postulate and the non-commutation of σx\sigma_x and σz\sigma_z.

6.5. ZYZ: same probabilities, different states

Let us replace the intermediate apparatus with one whose axis is yy, still selecting the positive channel. The last apparatus receives the state ∣y+⟩=12(∣z+⟩+i∣z−⟩)\ket{y+}=\frac{1}{\sqrt2}(\ket{z+}+i\ket{z-}), and

p(z+∣y+)=∣12∣2=12,p(z−∣y+)=∣i2∣2=12.p(z+\mid y+)=\Bigl|\frac{1}{\sqrt2}\Bigr|^2=\frac12,\qquad p(z-\mid y+)=\Bigl|\frac{i}{\sqrt2}\Bigr|^2=\frac12 .

The ZXZ and ZYZ sequences therefore give exactly the same results. Yet the states ∣x+⟩\ket{x+} and ∣y+⟩\ket{y+} are different: an apparatus with axis xx gives x+x+ with certainty for the former, and x+x+ or x−x- with probability 12\frac12 for the latter. Thus two states can have the same statistics for one axis and differ for another. What distinguishes them is the relative phase between their components along ∣z+⟩\ket{z+} and ∣z−⟩\ket{z-}, which is 11 in one case and ii in the other: a relative phase has physical consequences, whereas a global phase does not.

6.6. Selecting is not the same as doing nothing

In the ZXZ sequence, what would happen if neither of the two channels of the intermediate apparatus were blocked? If the atoms are detected, or if their trajectory makes it possible in principle to know which channel they passed through, the apparatus performs a measurement along xx, and the last apparatus receives a mixture of atoms in the states ∣x+⟩\ket{x+} and ∣x−⟩\ket{x-}, in equal proportions: each gives z±z\pm with probability 12\frac12. But if the two beams emerging from the intermediate apparatus could be perfectly recombined, without any information about the path taken, the state ∣z+⟩=12(∣x+⟩+∣x−⟩)\ket{z+}=\frac{1}{\sqrt2}(\ket{x+}+\ket{x-}) would be reconstituted, and the last apparatus would give z+z+ with certainty. A Stern—Gerlach apparatus therefore measures nothing until the path taken has been recorded in one way or another. Such recombination is extremely difficult to achieve with atoms; but this is exactly the situation in the interferometers studied in Theme 1, where the two paths are recombined coherently. We shall return to this connection between measurement and information when discussing composite systems.

7. An arbitrary axis

7.1. The observable

Finally, consider an apparatus whose axis has an arbitrary direction, specified by the unit vector

u=(sin⁡θcos⁡φ, sin⁡θsin⁡φ, cos⁡θ),\mathbf u=(\sin\theta\cos\varphi,\ \sin\theta\sin\varphi,\ \cos\theta),

where θ\theta is the angle with the zz-axis and φ\varphi the azimuth, measured from the xx-axis in the (x,y)(x,y) plane. How do we construct the corresponding observable σu\sigma_{\mathbf u}? Classically, the component of a magnetic moment along u\mathbf u is the linear combination u⋅μ=uxμx+uyμy+uzμz\mathbf u\cdot\boldsymbol\mu=u_x\mu_x+u_y\mu_y+u_z\mu_z. We assume that the same holds for quantum expectation values, ⟨σu⟩=ux⟨σx⟩+uy⟨σy⟩+uz⟨σz⟩\langle\sigma_{\mathbf u}\rangle=u_x\langle\sigma_x\rangle+u_y\langle\sigma_y\rangle+u_z\langle\sigma_z\rangle in every state. Since two self-adjoint operators that have the same expectation value in every state are equal, this imposes

σu=u⋅σ=sin⁡θcos⁡φσx+sin⁡θsin⁡φσy+cos⁡θσz=(cos⁡θe−iφsin⁡θeiφsin⁡θ−cos⁡θ).\begin{aligned} \sigma_{\mathbf u}=\mathbf u\cdot\boldsymbol\sigma=\sin\theta\cos\varphi \sigma_x+\sin\theta\sin\varphi \sigma_y+\cos\theta \sigma_z =\begin{pmatrix}\cos\theta&e^{-i\varphi}\sin\theta\\ e^{i\varphi}\sin\theta&-\cos\theta\end{pmatrix}. \end{aligned}
(9)

To obtain the matrix, we have combined sin⁡θ(cos⁡φ σx+sin⁡φ σy)\sin\theta(\cos\varphi\,\sigma_x+\sin\varphi\,\sigma_y), whose upper-right coefficient is sin⁡θ(cos⁡φ−isin⁡φ)=sin⁡θ e−iφ\sin\theta(\cos\varphi-i\sin\varphi)=\sin\theta\,e^{-i\varphi}. This matrix is self-adjoint, since it is a linear combination with real coefficients of self-adjoint matrices; this can also be checked directly, as its two off-diagonal coefficients are conjugates of one another.

7.2. A first check: quantisation along every axis

The preceding assumption is acceptable only if it respects fact (E1): an apparatus with an arbitrary axis must also give the two outcomes ±1\pm1. This is not obvious a priori, since σu\sigma_{\mathbf u} is a combination of the three matrices. Let us calculate its square. On expanding, the cross terms cancel in pairs owing to anticommutation, for example uxuy(σxσy+σyσx)=0u_xu_y(\sigma_x\sigma_y+\sigma_y\sigma_x)=0, leaving

σu2=(ux2+uy2+uz2) 1=1.\sigma_{\mathbf u}^2=(u_x^2+u_y^2+u_z^2)\,\mathbf{1}=\mathbf{1}.

The eigenvalues of σu\sigma_{\mathbf u} therefore satisfy λ2=1\lambda^2=1, and since its trace is zero, they are +1+1 and −1-1. The model does indeed predict quantisation along every axis, in accordance with experiment. We can also calculate the characteristic polynomial of the matrix (9) directly: (cos⁡θ−λ)(−cos⁡θ−λ)−sin⁡2θ=λ2−1(\cos\theta-\lambda)(-\cos\theta-\lambda)-\sin^2\theta=\lambda^2-1.

7.3. The eigenstates

The eigenvectors of σu\sigma_{\mathbf u} are

∣u+⟩=cos⁡θ2∣z+⟩+eiφsin⁡θ2∣z−⟩,∣u−⟩=−e−iφsin⁡θ2∣z+⟩+cos⁡θ2∣z−⟩.\ket{u+}=\cos\frac\theta2 \ket{z+}+e^{i\varphi}\sin\frac\theta2 \ket{z-},\qquad \ket{u-}=-e^{-i\varphi}\sin\frac\theta2 \ket{z+}+\cos\frac\theta2 \ket{z-}.
(10)

Let us check the first, setting c=cos⁡θ2c=\cos\frac\theta2 and s=sin⁡θ2s=\sin\frac\theta2, and using cos⁡θ=c2−s2\cos\theta=c^2-s^2 and sin⁡θ=2sc\sin\theta=2sc. The first component of σu∣u+⟩\sigma_{\mathbf u}\ket{u+} is

cos⁡θ c+e−iφsin⁡θ eiφs=(c2−s2)c+2s2c=c (c2+s2)=c,\cos\theta\,c+e^{-i\varphi}\sin\theta\,e^{i\varphi}s=(c^2-s^2)c+2s^2c=c\,(c^2+s^2)=c,

and the second is

eiφsin⁡θ c−cos⁡θ eiφs=eiφ(2sc2−(c2−s2)s)=eiφs (c2+s2)=eiφs.e^{i\varphi}\sin\theta\,c-\cos\theta\,e^{i\varphi}s=e^{i\varphi}\bigl(2sc^2-(c^2-s^2)s\bigr)=e^{i\varphi}s\,(c^2+s^2)=e^{i\varphi}s .

Thus σu∣u+⟩=∣u+⟩\sigma_{\mathbf u}\ket{u+}=\ket{u+}, and similarly one checks that σu∣u−⟩=−∣u−⟩\sigma_{\mathbf u}\ket{u-}=-\ket{u-}. We recover the special cases already constructed: θ=0\theta=0 gives ∣z+⟩\ket{z+}, θ=π2\theta=\frac\pi2 and φ=0\varphi=0 give ∣x+⟩\ket{x+}, and θ=π2\theta=\frac\pi2 and φ=π2\varphi=\frac\pi2 give ∣y+⟩\ket{y+}. For an arbitrary axis in the (x,y)(x,y) plane, we obtain ∣u+⟩=12(∣z+⟩+eiφ∣z−⟩)\ket{u+}=\frac{1}{\sqrt2}(\ket{z+}+e^{i\varphi}\ket{z-}), which is exactly the form inferred in Theme 1 to describe the continuous range of preparations ∣u+⟩\ket{u+}.

7.4. A prediction: the cosine-squared law

Let us prepare an atom in the state ∣z+⟩\ket{z+} and analyse it along u\mathbf u. The Born rule gives

p(u+∣z+)=∣⟨u+|z+⟩∣2=cos⁡2θ2,p(u−∣z+)=sin⁡2θ2.p(u+\mid z+)=|\braket{u+}{z+}|^2=\cos^2\frac\theta2,\qquad p(u-\mid z+)=\sin^2\frac\theta2 .
(11)

The probability depends only on the angle θ\theta between the two axes, and not on the azimuth, as required by rotational symmetry about the zz-axis. It is 11 for identical axes, 12\frac12 for perpendicular axes and 00 for opposite axes: an atom prepared in the z+z+ channel emerges with certainty in the −- channel of a reversed apparatus, which is evident since reversing the apparatus amounts to exchanging its two channels. Between these values, the cos⁡2(θ/2)\cos^2(\theta/2) law is a prediction of the model, confirmed experimentally. We deduce the expectation value ⟨σu⟩=cos⁡2θ2−sin⁡2θ2=cos⁡θ\langle\sigma_{\mathbf u}\rangle=\cos^2\frac\theta2-\sin^2\frac\theta2=\cos\theta, which is the projection onto u\mathbf u of a unit vector directed along zz: on average, everything behaves as for the projection of a classical magnetic moment, even though each individual measurement gives only ±1\pm1.

For a zz, u\mathbf u, zz sequence with selection of the u+u+ channel, formula (8) gives the final probabilities, relative to the initial beam,

p(u+,z+)=cos⁡2θ2⋅cos⁡2θ2=cos⁡4θ2,p(u+,z−)=cos⁡2θ2 sin⁡2θ2.p(u+,z+)=\cos^2\frac\theta2\cdot\cos^2\frac\theta2=\cos^4\frac\theta2, \qquad p(u+,z-)=\cos^2\frac\theta2\,\sin^2\frac\theta2 .

For θ=π2\theta=\frac\pi2, we recover the values 14\frac14 of the ZXZ sequence.

7.5. The beam emerging from the oven

One final consequence deserves emphasis. After a choice of global phase that makes α\alpha real and positive, every normalised state α∣z+⟩+β∣z−⟩\alpha\ket{z+}+\beta\ket{z-} can be written in the form (10): it suffices to set α=cos⁡θ2\alpha=\cos\frac\theta2 and β=eiφsin⁡θ2\beta=e^{i\varphi}\sin\frac\theta2, which is always possible since ∣α∣2+∣β∣2=1|\alpha|^2+|\beta|^2=1. In other words, every pure state is the state ∣u+⟩\ket{u+} for some direction u\mathbf u, and therefore gives a certain outcome for an apparatus oriented along that direction. The pure states of the system are thus in correspondence with directions in space, that is, with the points of a sphere: this is the Bloch sphere, which we shall study in Lesson 10.

Yet the beam emerging from the oven gives two spots of equal intensity whatever the orientation of the apparatus. From the preceding argument, no pure state can have this property. The atoms emerging from the oven are therefore not all in the same pure state: they form a statistical mixture of atoms in different, randomly oriented states. To describe such a beam, we shall have to generalise the notion of state and introduce mixed states, which we shall encounter in Theme 4.

8. Summary

Let us summarise the approach followed in this lesson, which is the one used to construct the quantum model of any system.

  • The possible outcomes of a measurement determine the spectrum of the corresponding observable: here, two values for each axis.
  • Reproducibility experiments show that the states prepared by selection are eigenstates and, together with a simplifying assumption, determine the dimension of the state space: here, two.
  • The probabilities measured in experiments in series determine, through the Born rule, the moduli of the inner products between the eigenstates of the different observables, and ultimately the eigenvectors themselves, up to phase conventions. It was at this stage that complex numbers proved indispensable.
  • The spectral decomposition then reconstructs each observable from its eigenvalues and its eigenprojectors: this gives the Pauli matrices.
  • The model can then predict the results of other experiments: the reappearance of the z−z- channel in the ZXZ sequence, quantisation along every axis, the cos⁡2(θ/2)\cos^2(\theta/2) law, and the impossibility of describing the beam emerging from the oven by a pure state.

The following lessons will apply the same approach to richer systems whose state space is infinite-dimensional: a particle in space. For further discussion of two-level calculations and measurement, see the MIT notes, course 8.321, session 3[1], as well as J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics[2], Chapter 1, whose presentation of the Stern—Gerlach experiments inspired this lesson.

9. References