Chapter 1 — Discrete outcomes in the Stern–Gerlach experiment
In the Stern–Gerlach experiment, a beam of silver atoms travels through a magnetic field that varies with position. The field deflects the atoms according to their magnetic properties. If those properties behaved like randomly oriented classical magnetic needles, the deflections would span a continuous range. The observed beam instead separates into two branches.
Otto Stern and Walther Gerlach observed this splitting in 1922. We will idealize it as a measurement with two possible outcomes, each identified by the branch an atom follows. Repeating the experiment with the same preparation gives a fraction of atoms in each branch. The quantum model predicts the probabilities governing those fractions.
Assumes: complex numbers, vectors, and matrices; no prior quantum mechanics. Introduces: the state ket, the Born rule, superposition, overall versus relative phase, projective measurement, observables, and unitary evolution. Used later in: every chapter — this is the working vocabulary for a single qubit, developed further in Chapter 3.
Complex amplitudes and interference
Choose the magnet's measurement axis as the \(z\) axis and label its two outputs \(z+\) and \(z-\). We model the internal degree of freedom responsible for this choice of output in \(\mathbb C^2\). This leaves the atom's position and its other internal degrees of freedom outside the model.
We begin with preparations represented by a single normalized vector, called pure states. Fix an orthonormal basis associated with the two \(z\) outputs. The vector's complex coordinates in this basis are called probability amplitudes. Their relation to observed frequencies is the physical rule we need to specify.
For complex conjugation and squared magnitude, we use the notation
\[ z = a + ib, \]
with real \(a\) and \(b\), \(i^2=-1\), and conjugate \(z^*=a-ib\). The squared magnitude of this number \(z\) is
\[ |z|^2 = z^* z = a^2 + b^2. \]
The squared magnitude is nonnegative; normalization will make the squared magnitudes of the two coordinates sum to one. These are both requirements for assigning probabilities to two exhaustive outcomes.
Write the amplitudes for the two outputs as \(\alpha\) and \(\beta\). In Dirac notation the state vector is written as a ket:
\[ |\psi\rangle = \begin{pmatrix} \alpha \\ \beta \end{pmatrix}. \]
Here the ket is displayed as a column in the fixed orthonormal basis. Its entries are the amplitudes for \(z+\) and \(z-\), in that order. Normalization requires
\[ |\alpha|^2+|\beta|^2=1. \]
The Born rule assigns probability \(|\alpha|^2\) to the first output and \(|\beta|^2\) to the second. The state and the chosen measurement supply the amplitudes; the Born rule converts them to probabilities. Normalization makes these probabilities sum to one, as required for the ideal measurement in which every atom gives one of the two outcomes.
A normalized column is a representative of a pure physical state. Two columns differing only by a common complex factor of magnitude one describe the same state, because this factor cancels from every Born probability.
The complex coordinates contain information beyond the two probabilities for this particular magnet. In a measurement along a different axis, contributions from the coordinates can cancel when forming the amplitude for a single output. This cancellation is an example of interference. We will see it explicitly after introducing the second measurement basis.
Basis states and superposition
The internal angular momentum tracked by this two-state model is called spin. It is an internal property, distinct from the angular momentum of the atom's motion through the apparatus. A spin component is the angular momentum measured along a chosen axis. The magnet with measurement axis \(z\) distinguishes the two possible values of that component.
The basis vectors represent preparations that give the corresponding output with certainty:
\[ |z+\rangle = \begin{pmatrix} 1 \\ 0 \end{pmatrix}, \qquad |z-\rangle = \begin{pmatrix} 0 \\ 1 \end{pmatrix}. \]
In this basis the general pure-state vector is
\[ |\psi\rangle = \alpha |z+\rangle + \beta |z-\rangle. \]
This expansion is a superposition in the \(z\) basis. Its coefficients give the amplitudes for a \(z\)-basis measurement, whose probabilities we have just specified. Spin is being described by a vector in a complex state space; these coefficients are not two simultaneous classical directions of the atom.
To compute amplitudes for other measurements, we need inner products. The conjugate transpose of a ket is written as a bra:
\[ \langle\psi| = \begin{pmatrix} \alpha^* & \beta^* \end{pmatrix}. \]
The product of a bra with a ket is the inner product \(\langle\phi|\psi\rangle\). If the first vector is a normalized basis state for a measurement outcome and the second is the prepared state, this inner product is the amplitude for that outcome. Its squared magnitude is the Born probability. In particular, normalization becomes
\[ \langle\psi|\psi\rangle = |\alpha|^2 + |\beta|^2 = 1. \]
With its standard inner product, this finite-dimensional complex vector space is a Hilbert space; completeness is automatic in finite dimension. The state vectors use its linear structure, and the measurement probabilities use its inner product.
Sequential measurements along perpendicular axes
An ideal Stern–Gerlach filter separates the two outputs of a spin-component measurement and blocks one of them. Keeping only atoms from the other output is called postselection.
We assume that each retained atom has the basis state associated with its output. This is the state-update rule for an ideal projective measurement with one-dimensional output subspaces. A repeated measurement along the same axis therefore gives the same outcome with certainty. Throughout the sequence below, we also assume that the spin state does not evolve between magnets.
Now rotate the measurement axis from \(z\) to the perpendicular \(x\) axis. The new magnet measures a different spin component. Label its outputs \(x+\) and \(x-\). In the two-state spin model, we choose their basis vectors as
\[ |x+\rangle = \frac{|z+\rangle + |z-\rangle}{\sqrt{2}}, \qquad |x-\rangle = \frac{|z+\rangle - |z-\rangle}{\sqrt{2}}. \]
This relation between the two bases is part of the physical model for perpendicular spin measurements. Both vectors have unit norm, and their inner product vanishes:
\[ \langle x+|x-\rangle = \frac12 \bigl(\langle z+| + \langle z-|\bigr) \bigl(|z+\rangle - |z-\rangle\bigr) = \frac12(1-1) = 0. \]
The vanishing inner product has a direct measurement meaning. For an atom prepared in \(|x-\rangle\), the two \(z\)-basis terms contribute equal and opposite amplitudes to the same \(x+\) outcome. They cancel, so that outcome has probability zero. Such cancellation is destructive interference.
The two terms retain the relative sign fixed by the preparation, with no record distinguishing the alternatives. They form a coherent superposition. We add their contributions to a single outcome's amplitude before taking its squared magnitude. By contrast, distinct recorded outcomes are mutually exclusive, so their probabilities add. Adding the two amplitudes for the separate outputs of a \(z\) measurement would not give the probability of detecting an atom in either output.
In step 1, pass a beam of atoms through a \(z\)-oriented magnet. Block the \(z-\) output and retain the \(z+\) output. Under the ideal measurement assumption, each survivor is prepared in
\[ |\psi_1\rangle = |z+\rangle. \]
A second \(z\)-oriented magnet would send every survivor to its \(z+\) output. The inner product of the prepared state with that output state is \(1\), so the probability is \(1\).
In step 2, instead send the survivors through a perpendicular magnet and retain only its \(x+\) output. The amplitude for passing this filter, often called a transition amplitude, is
\[ \langle x+|z+\rangle = \frac{1}{\sqrt{2}} \bigl(\langle z+| + \langle z-|\bigr)|z+\rangle = \frac{1}{\sqrt{2}}. \]
The Born rule gives the passing probability conditional on the preparation:
\[ P(x+\mid z+) = \Bigl|\frac{1}{\sqrt{2}}\Bigr|^2 = \frac12. \]
The expected surviving fraction is one half. Each retained atom, initially prepared in \(z+\), now has the normalized state
\[ |\psi_2\rangle = |x+\rangle = \frac{|z+\rangle + |z-\rangle}{\sqrt{2}}. \]
The probability describes the expected fraction passing the filter. The normalized state describes an atom conditional on having passed. Its norm is one even though only half the incoming atoms are expected to survive.
In step 3, measure the \(z\) component of these survivors and count both outputs. The amplitudes are
\[ \langle z+|x+\rangle = \frac{1}{\sqrt{2}}, \qquad \langle z-|x+\rangle = \frac{1}{\sqrt{2}}, \]
so
\[ P(z+\mid x+) = P(z-\mid x+) = \frac12. \]
After selection at the second magnet, the retained state gives a certain \(x\) result but two equally probable \(z\) results. The intervening filter has changed the preparation. This rules out a description in which the middle apparatus merely reads a pre-existing value while leaving the initially certain \(z\) value intact. The sequence alone does not rule out models in which measurement disturbs pre-existing values.
Starting from 100 atoms already prepared in \(z+\), the expected count at the retained \(x+\) output is 50. The expected counts at the two final \(z\) outputs are 25 each. These final counts are relative to the initial 100 atoms; the final probabilities of one half are conditional on passing the middle filter. Repeated runs fluctuate around the expected counts even when the ideal probabilities remain fixed.
Overall and relative phase
The two probabilities for a fixed magnet orientation generally do not specify the whole pure state. To see what information matters, first multiply both amplitudes by the same factor \(e^{i\gamma}\), with \(\gamma\) real. This changes their common phase. For any normalized measurement ket, the probability remains
\[ \bigl|\langle a|e^{i\gamma}\psi\rangle\bigr|^2 = |e^{i\gamma}|^2 \,\bigl|\langle a|\psi\rangle\bigr|^2 = \bigl|\langle a|\psi\rangle\bigr|^2. \]
Thus \(|\psi\rangle\) and \(e^{i\gamma}|\psi\rangle\) represent the same pure state. The common factor is an overall phase. Mathematically, a ray consists of all nonzero complex scalar multiples of a nonzero vector. Its normalized representatives differ only by an overall phase.
Changing one amplitude's phase while holding the other fixed can change a measurement probability. With the basis vectors fixed, the phase difference between two nonzero amplitudes is their relative phase. The previously defined states provide an example:
\[ \frac{|z+\rangle + |z-\rangle}{\sqrt{2}} = |x+\rangle, \qquad \frac{|z+\rangle - |z-\rangle}{\sqrt{2}} = |x-\rangle \]
Both states give equal probabilities for the two \(z\)-basis outcomes. Their second amplitudes differ by a minus sign, so the relative phases differ by half a turn. An \(x\)-basis measurement distinguishes them with certainty: the first state gives the plus outcome, and the second gives the minus outcome. In computing each overlap, the two contributions add or cancel according to this relative phase.
The Born rule and projection
The filters can now be described for any finite-dimensional pure-state system. Consider an ideal projective measurement with a normalized ket for each outcome. These kets \(\{|a_j\rangle\}\) form an orthonormal basis: each has norm one, distinct kets are orthogonal, and together they span the state space. The outcomes therefore form a complete set for this measurement. For a normalized input ket, the Born rule gives the probability of outcome \(j\) as
\[ P(j) = \bigl|\langle a_j|\psi\rangle\bigr|^2. \]
Born proposed the probabilistic interpretation of the wavefunction in 1926. To describe what a filter retains, we also need the component along the selected outcome ket. The orthogonal projector onto that one-dimensional subspace is
\[ \Pi_j = |a_j\rangle\langle a_j|. \]
Acting on the input ket, this projector retains its component along the selected ket and removes the orthogonal component. The squared norm of the retained component is the probability of passing the filter. Equivalently,
\[ P(j) = \langle\psi|\Pi_j|\psi\rangle. \]
These orthogonal projectors are self-adjoint and idempotent. Completeness of the outcome basis gives
\[ \sum_j \Pi_j = I. \]
Multiplying this identity on the left by the input bra and on the right by its normalized ket shows that the probabilities sum to one.
Now condition on the observed outcome \(j\), assuming it has nonzero probability. The ideal filter prepares the state represented by \(|a_j\rangle\). Applying the projector and normalizing gives the same state, up to overall phase:
\[ |\psi_j\rangle = \frac{\Pi_j|\psi\rangle}{\sqrt{\langle\psi|\Pi_j|\psi\rangle}}. \]
The denominator is the norm of the projected ket. A zero-probability outcome has no conditional state defined by this formula.
An outcome may instead correspond to a subspace of dimension greater than one. For a complete family of mutually orthogonal subspaces, use their orthogonal projectors in the same probability formula. The normalized-projection update then specifies an ideal measurement that preserves the component within the selected subspace; this assumption is part of the measurement model.
If a measurement occurs but its outcome is not retained, subsequent predictions must average over the possible conditional states with their outcome probabilities. In general, no single ket describes that ensemble. A density operator, introduced later, will represent this statistical mixture.
Operators and observables
A filter labels an output channel. To describe a measured quantity, we must also assign a numerical reading to each channel. Define an operator that multiplies each outcome component by its assigned reading. Its action on superpositions satisfies
\[ A\bigl(c_1|u\rangle + c_2|v\rangle\bigr) = c_1 A|u\rangle + c_2 A|v\rangle. \]
In the orthonormal basis used here, the adjoint \(A^\dagger\) is represented by the conjugate transpose. Assigning real readings to mutually orthogonal outcome subspaces gives an operator satisfying
\[ A = A^\dagger, \]
This is the Hermitian condition. In the ideal projective model, such an operator represents an observable: a quantity measured with those numerical outcomes.
Reading the construction backwards recovers the measurement. The spectral theorem guarantees that a Hermitian operator has real eigenvalues and a complete orthonormal eigenbasis — exactly the readings and outcome subspaces we started from. For a ket in one outcome subspace,
\[ A|a_j\rangle = a_j |a_j\rangle. \]
The normalized ket \(|a_j\rangle\) is called an eigenstate of \(A\), and preparing it makes the reading \(a_j\) certain. If several basis kets share one eigenvalue, that reading corresponds to their whole eigenspace. So an observable and its ideal measurement are two views of the same data: the operator packages the readings and the eigenbasis packages the outcomes. More general measurement procedures require a broader description than this projective one.
For the two-output magnets, first assign the dimensionless readings plus one and minus one to their respective outputs. In the fixed z basis, the resulting operators are the Pauli matrices
\[ \sigma_z = \begin{pmatrix} 1 & 0 \\ 0 & -1 \end{pmatrix}, \qquad \sigma_x = \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}. \]
Both are Hermitian. Their eigenstates are exactly the states prepared by the corresponding filters:
\[ \sigma_z |z\pm\rangle = \pm |z\pm\rangle, \qquad \sigma_x |x\pm\rangle = \pm |x\pm\rangle. \]
The Pauli matrices carry no units. To express the angular momentum measured in an ideal spin-\(1/2\) model, multiply by half the reduced Planck constant, denoted by h-bar. This universal constant sets the quantum scale of angular momentum. The physical component operators are
\[ S_z = \frac{\hbar}{2}\sigma_z, \qquad S_x = \frac{\hbar}{2}\sigma_x, \]
An \(S_z\) measurement therefore returns \(+\hbar/2\) or \(-\hbar/2\), in joule-seconds. The labels \(z+\) and \(z-\) denote the eigenstates for these respective readings.
The two magnet orientations do not admit a common eigenbasis. Their Pauli operators do not commute:
\[ [\sigma_z,\sigma_x] \equiv \sigma_z\sigma_x - \sigma_x\sigma_z \ne 0. \]
The commutator \([A,B]\) compares the compositions of \(A\) and \(B\), with \(B\) applied before \(A\) in the first product. For Hermitian operators, a common orthonormal eigenbasis exists exactly when they commute. For the particular pair \(\sigma_z\) and \(\sigma_x\), the displayed eigenstates show more: there is no common eigenvector.
Consequently, no ket predicts both a certain \(z\) outcome and a certain \(x\) outcome. The intermediate filter in the three-magnet sequence replaces one such certainty with the other. A measurement includes this conditional state change; composing observable matrices alone does not describe the filtering procedure. Neither noncommutation nor that sequence by itself rules out every model with pre-existing values. The sequence rules out a readout that leaves the original certainty undisturbed.
Repeated measurements on identically prepared atoms yield a distribution of readings. Its mean, called the expectation value of \(A\) in the normalized state \(|\psi\rangle\), is
\[ \langle A\rangle_\psi = \langle\psi|A|\psi\rangle = \sum_j a_j P(j). \]
A finite sample average estimates this mean. The mean need not be an allowed outcome of a single measurement.
For \(|x+\rangle\), a measurement of \(\sigma_z\) has readings \(+1\) and \(-1\), each with probability one half. Its expectation is \(0\), although an individual run never reports zero.
Hermiticity also guarantees reality directly: if \(q = \langle\psi|A|\psi\rangle\) and \(A = A^\dagger\), then \(q^* = q\).
Unitary time evolution
A prepared atom can evolve before reaching the next filter. Suppose the tracked system is closed over that interval: interactions with untracked degrees of freedom can be neglected, while any prescribed applied field is included in its dynamics. Quantum mechanics models this evolution by a unitary operator \(U\):
\[ U^\dagger U = U U^\dagger = I. \]
If \(|\psi'\rangle = U|\psi\rangle\), then
\[ \langle\psi'|\psi'\rangle = \langle\psi|U^\dagger U|\psi\rangle = \langle\psi|\psi\rangle. \]
The evolved ket remains normalized, so it still gives a complete probability distribution for any subsequent ideal measurement. Unitarity preserves inner products between any two evolving states as well.
To determine the evolution, specify the system's energy operator \(H\), called its Hamiltonian. It is Hermitian, with energy eigenvalues measured in joules. The law connecting this operator to the changing state is Schrödinger’s equation:
\[ i\hbar \frac{d}{dt}|\psi(t)\rangle = H|\psi(t)\rangle. \]
With time \(t\) in seconds and \(\hbar\) in joule-seconds, the differential operator on the left has the same energy units as the Hamiltonian:
\[ [\hbar\, d/dt] = (\mathrm{J\,s})(1/\mathrm{s}) = \mathrm{J} = [H]. \]
When the Hamiltonian \(H\) is constant in time, the initial-value solution is
\[ |\psi(t)\rangle = U(t)\,|\psi(0)\rangle, \qquad U(t) = e^{-iHt/\hbar}. \]
The exponent contains the dimensionless combination \(Ht/\hbar\). Since \(H\) is Hermitian, \(U(t)\) is unitary.
Consider an energy splitting between the two z-basis states, described by
\[ H = \frac{\hbar\omega}{2}\sigma_z, \]
Here \(\omega\) is an angular frequency, measured in radians per second. The two basis states acquire phases at opposite rates because their energies have opposite signs in this choice of energy reference. Starting with the state \(|x+\rangle\) gives
\[ |\psi(t)\rangle = \frac{e^{-i\omega t/2}|z+\rangle + e^{+i\omega t/2}|z-\rangle}{\sqrt{2}}. \]
Both amplitudes change with time through their phase factors. Their squared magnitudes remain one half, so the \(z\)-basis probabilities stay fixed. To detect the changing relative phase, measure in the \(x\) basis. The overlap adds the two phase-dependent contributions, giving
\[ P(x+;t) = \bigl|\langle x+|\psi(t)\rangle\bigr|^2 = \cos^2\bigl(\omega t/2\bigr). \]
Thus the probability of the plus output oscillates as the delay before measurement changes. The unchanged probabilities in the energy basis conceal a changing pure state. A measurement sensitive to the relative phase reveals that evolution.
Physical implementation of the 1922 experiment
The original apparatus sent neutral silver atoms from a heated source through a narrow beam into an inhomogeneous magnetic field. The outgoing atoms left a deposit on a collecting plate. Their positions revealed two separated branches. The sequential filters used above are an idealized extension of that experiment.
The separation arises because an atom's internal magnetic state affects the force exerted by the field. In an approximation that treats each resolved branch as a classical trajectory, the dipole interaction energy gives
\[ \mathbf F \approx \boldsymbol\nabla(\boldsymbol\mu\cdot\mathbf B), \]
The magnetic moment \(\boldsymbol\mu\) describes the atom's coupling to the magnetic field and has units of joules per tesla. The field \(\mathbf B\) is measured in tesla. Their scalar product is an energy; its spatial gradient has units of newtons. The two internal components experience different forces and can therefore reach distinguishable positions.
Stern and Gerlach published their result in 1922, three years before electron spin was proposed. The spin account is a later explanation of the observed splitting.
In the modern account, ground-state silver has electronic angular-momentum quantum number one half. The two allowed projections along the magnet axis supply the basis states of our two-dimensional model. This reduction neglects nuclear-spin structure and other atomic degrees of freedom when describing the resolved branches. The atom's position must be included as well if we want to explain how those branches form.
The magnetic field directs the two internal basis components into different spatial paths. When both components are initially present coherently, the resulting joint state generally cannot be written as an internal ket times a spatial ket: each internal component is associated with its own path. This failure to factor is called entanglement and will be developed in the next chapter. An input containing just one of the internal basis components need not become entangled in this ideal model.
The separation of paths can be described by unitary evolution of the internal and spatial degrees of freedom together. A blocker then transmits one path, or a detector registers an arrival. Conditional on a selected output, the ideal spin-state update is the projection described earlier. This is an operational model of preparation and readout; it does not derive the occurrence of an individual detector outcome from unitary evolution.
Real beams have finite width, and real devices can lose atoms or misidentify an output. Field misalignment changes which component is measured. Such effects require changes to the preparation, evolution, or detector model before its predictions are compared with counts. The Born rule remains the probability rule, but the ideal two-channel formulas alone need not describe the imperfect apparatus.
The two-component ket therefore captures a selected part of the atom's behavior. The filters already illustrate ways to prepare and read that part. To use it for quantum computation, one must also implement controlled operations and keep unwanted interactions sufficiently small over the required sequence. Those physical capabilities are additional information beyond the dimension of the state space.
Common conceptual errors
For an outcome amplitude computed from a normalized state and a normalized measurement ket, use the squared magnitude \(|c|^2 = c^*c\). Squaring the complex amplitude itself can give a negative or nonreal number.
The column \((\alpha,\beta)^T\) gives coordinates in the fixed z basis. Choosing another coordinate basis generally changes the entries without changing the physical state. Rotating a measuring magnet changes the measurement; the state can still be expressed in the original basis.
Equal probabilities in one basis do not establish equality of pure states. The two x-basis states have equal z-basis probabilities and are nevertheless distinguished by an x-oriented magnet. Overall phase is the freedom that leaves every measurement probability unchanged.
An expectation value is a distribution's mean, estimated by repeated measurements. It need not be one of the values an individual measurement can report.
In the ideal filter model, postselection conditions on an output of nonzero probability and normalizes the corresponding projection. A unitary operator on the internal state alone cannot perform this selective removal. Unitary separation of the spatial paths is an earlier, distinct part of the procedure.
A Hermitian observable specifies the readings and outcome subspaces of an ideal projective measurement. A unitary operator specifies closed-system evolution. The Pauli matrices satisfy both algebraic conditions, so an operator's physical role must be stated along with its matrix.
A numerical simulation can propagate the specified ket and evaluate Born probabilities for the chosen measurement. Computing these predictions does not reproduce the physical atom or its interaction with a detector.
A two-dimensional state space does not specify how a device is operated. Preparation, control, readout, and suppression of unwanted interactions each require a physical implementation.
Self-assessment
What specifies a pure-state experiment? Give a normalized initial ket, with overall phase irrelevant. Specify the Hamiltonian for each closed-system evolution interval. For each ideal projective measurement, give the outcome projectors and their numerical readings. State which outcomes, if any, are retained for subsequent steps.
How is an outcome probability computed? For outcome \(j\) in an orthonormal measurement basis, the amplitude is \(\langle a_j|\psi\rangle\). For the normalized input state, its squared magnitude gives \(P(j) = |\langle a_j|\psi\rangle|^2\).
Why normalize in the \(z\) basis? The condition \(|\alpha|^2 + |\beta|^2 = 1\) makes the probabilities of the two complete \(z\)-basis outcomes sum to one. It describes the ideal two-outcome measurement, without an additional loss channel.
Why require a Hermitian observable? An ideal projective observable assigns real readings to mutually orthogonal, complete outcome subspaces. An operator \(A\) constructed from those spectral data is Hermitian. A non-Hermitian operator cannot have both a real spectrum and a complete orthonormal eigenbasis, so it cannot represent an observable in this model.
How does evolution differ from postselection? Closed-system unitary evolution preserves all inner products. An ideal filter instead selects an output. Conditional on that output having nonzero probability, the normalized projection describes the atoms retained for the next operation.
What follows the sequence \(z+\), then \(x+\), then \(z\)? Prepare the first state and retain the plus output of the perpendicular filter. Among those survivors, the final outcomes \(z+\) and \(z-\) each have probability one half. If the intermediate \(x\) measurement left the prepared \(z+\) state unchanged, the final result would instead be certainly \(z+\). The sequence demonstrates the change in preparation caused by this filter.
We can now follow a pure state through closed-system evolution and predict the outcomes of ideal filters, including the conditional state passed to the next step. The next chapter applies this framework to two atoms. Some joint pure states factor into one ket for each atom; others do not. Understanding the latter requires keeping the joint amplitudes together.
Sources
[R037] M. Born, “Zur Quantenmechanik der Stoßvorgänge,” Zeitschrift für Physik 37, 863–867 (1926). DOI: 10.1007/BF01397477.
[R038] W. Gerlach and O. Stern, “Der experimentelle Nachweis der Richtungsquantelung im Magnetfeld,” Zeitschrift für Physik 9, 349–352 (1922). DOI: 10.1007/BF01326983.
[R039] B. Friedrich and D. Herschbach, “Stern and Gerlach: How a Bad Cigar Helped Reorient Atomic Physics,” Physics Today 56(12), 53–59 (2003). DOI: 10.1063/1.1650229.
[R040] P. A. M. Dirac, The Principles of Quantum Mechanics, 4th ed., Oxford University Press (1958; reissued 1981). ISBN: 978-0-19-852011-5.
[R041] E. Schrödinger, “Quantisierung als Eigenwertproblem (Vierte Mitteilung),” Annalen der Physik 386, 109–139 (1926). DOI: 10.1002/andp.19263861802.
[R004] M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information: 10th Anniversary Edition, Cambridge University Press (2010). DOI: 10.1017/CBO9780511976667.