Fire electrons one at a time at a pair of slits and record where each lands: a dot here, a dot there, apparently at random — and after ten thousand dots, the fringes of Young’s experiment, as if each electron had passed through both slits and interfered with itself. The last chapter of the Year 1 volume met this strangeness: light comes in photons, matter has a wavelength λ=h/p, and position and momentum cannot both be sharp. What it did not give is the equation — the rule that says how the wave associated with a particle evolves, as Newton’s law says how a trajectory evolves. That equation is Schrödinger’s (1926), and with it the quantum description becomes a calculable theory: one writes a wave functionψ(x,t), whose squared modulus is the probability of finding the particle, and one solves a linear partial differential equation. This chapter states the equation, explains what ψ means and what it does not, finds its stationary states and the free particle’s wave packets, and shows where Newton’s mechanics re-emerges — for everything heavier than a molecule.
A transmission electron microscope: electrons accelerated to a hundred kilovolts have a wavelength of a few picometres, and their wave function, focused by magnetic lenses, images the atoms of a crystal.
30.1 What the Year 1 volume left us
Proposition 30.1(Particles are waves)
A photon of frequency ν has the energy E=hν=ℏω and the momentum p=h/λ=ℏk; a material particle of momentum p behaves as a wave of wavelength λ=h/p (de Broglie), with E=ℏω and p=ℏk as well — electron diffraction by crystals, neutron and even molecule interferometry confirm it; and no state can have a position and a momentum sharper than ΔxΔp≥ℏ/2 (Heisenberg). With ℏ=h/2π=1.05×10−34Js: an electron at 100eV has λ=0.12nm, the size of an atom; a dust grain of 1×10−15kg at 1mm/s has λ=7×10−16m, smaller than a nucleus — it behaves classically.
Proof. Recalled from the Year 1 volume; the task now is the dynamics of the wave. ∎
An electron diffraction pattern recorded in an electron microscope (a Kikuchi pattern from a crystal): electrons of a few picometres’ wavelength, diffracted by the atomic planes — matter behaving as a wave (NIST).
30.2 The wave function
Definition 30.2(Wave function and Born’s rule)
The state of a particle moving along x is described by a complex wave functionψ(x,t) such that
dP=∣ψ(x,t)∣2dx
is the probability of finding the particle between x and x+dx at the time t (Born’s rule); hence the normalisation∫−∞∞∣ψ∣2dx=1, and the expectation values ⟨x⟩=∫x∣ψ∣2dx, Δx2=⟨x2⟩−⟨x⟩2. Wave functions obey the superposition principle: if ψ1 and ψ2 are possible states, so is αψ1+βψ2 (normalised), and the probabilities then contain the interference term 2Re(α∗βψ1∗ψ2). A global phase factor eiθ changes nothing; a relative phase between two superposed components changes the interference.
Remark 30.3(What ψ is not)
ψ is not a classical field spread in space like E, nor a cloud of matter: the electron is always found whole, at one point. ψ is an amplitude of probability, the object that interferes; only ∣ψ∣2 is observed, and only statistically, by repeating the experiment on identically prepared particles. The dots of the double slit are single electrons; the fringes are ∣ψ1+ψ2∣2.
30.3 The Schrödinger equation
Theorem 30.4(Schrödinger equation)
A particle of mass m in the potential energy V(x) has a wave function obeying
iℏ∂t∂ψ=−2mℏ2∂x2∂2ψ+V(x)ψ.
It is linear (superposition), first order in time (the wave function at one instant determines it at all later times), and complex (ψ cannot be real: the i is essential). It is a postulate, justified by its consequences.
Proof. A motivation rather than a proof. A free particle of sharp momentum p=ℏk and energy E=ℏω should be a plane wave ψ=Aei(kx−ωt), for which iℏ∂tψ=ℏωψ=Eψ and −(ℏ2/2m)∂x2ψ=(ℏ2k2/2m)ψ=(p2/2m)ψ. The classical relation E=p2/2m is thus reproduced by the equation iℏ∂tψ=−(ℏ2/2m)∂x2ψ, and the natural generalisation for E=p2/2m+V adds Vψ. Notice that the equation must be linear to hold for superpositions of plane waves, and first order in time with an i so that e−iωt and not cosωt is the time dependence — a real wave equation would give E2∝p4. ∎
Proposition 30.5(Conservation of probability)
The probability density ρ=∣ψ∣2 and the probability current
j=mℏIm(ψ∗∂x∂ψ)=2imℏ(ψ∗∂xψ−ψ∂xψ∗)
obey ∂tρ+∂xj=0: probability is neither created nor destroyed, and ∫∣ψ∣2dx stays equal to 1. For a plane wave Aeikx, j=∣A∣2ℏk/m=ρv: density times velocity, as for any flow.
Proof.∂t(ψ∗ψ)=ψ∗∂tψ+ψ∂tψ∗; from the equation and its conjugate, ∂tψ=(iℏ/2m)∂x2ψ−(i/ℏ)Vψ, ∂tψ∗=−(iℏ/2m)∂x2ψ∗+(i/ℏ)Vψ∗; the V terms cancel and the rest is (iℏ/2m)(ψ∗∂x2ψ−ψ∂x2ψ∗)=−∂xj. ∎
Proposition 30.6(Momentum, expectation values and the classical limit)
The momentum of a state is obtained by the operator p^=−iℏ∂x (which gives ℏk on a plane wave): ⟨p⟩=∫ψ∗(−iℏ∂xψ)dx, and the energy by H^=−(ℏ2/2m)∂x2+V. The expectation values obey (Ehrenfest)
dtd⟨x⟩=m⟨p⟩,dtd⟨p⟩=−⟨dxdV⟩:
the centre of the wave packet follows Newton’s law exactly when V is at most quadratic, and approximately whenever the packet is narrow compared with the scale over which V varies — the classical limit, which holds for any macroscopic object because ℏ is so small.
Proof. Differentiate ⟨x⟩=∫x∣ψ∣2, use ∂tρ=−∂xj and integrate by parts: d⟨x⟩/dt=∫jdx=⟨p⟩/m. The second relation is the same computation one order higher (admitted here). ∎
30.4 Stationary states
Theorem 30.7(Stationary states)
If V does not depend on time, the equation has solutions of the separated form
ψ(x,t)=φ(x)e−iEt/ℏ,−2mℏ2φ′′+Vφ=Eφ
(the time-independent Schrödinger equation): the stationary states, of well-defined energy E, whose probability density ∣φ∣2 does not move. The energies E for which φ is acceptable (bounded, normalisable for a bound state) are the energy levels; the general solution is a superposition ψ=∑ncnφn(x)e−iEnt/ℏ, and a superposition of two levels has a density that oscillates at the Bohr frequencyν21=(E2−E1)/h — the frequency of the light the atom emits or absorbs.
Proof. Insert φ(x)f(t): iℏf′/f=(−ℏ2φ′′/2m+Vφ)/φ, a function of t equal to a function of x, hence a constant E; f=e−iEt/ℏ. For two levels, ∣c1φ1e−iE1t/ℏ+c2φ2e−iE2t/ℏ∣2=∣c1φ1∣2+∣c2φ2∣2+2Re(c1∗c2φ1∗φ2e−i(E2−E1)t/ℏ). Completeness of the stationary states is admitted. ∎
A stationary state: the spatial profile φ(x) (here the ground state of a well) multiplied by a phasor turning at the rate E/ℏ; real and imaginary parts rotate into each other while ∣ψ∣2 stands still.
A superposition of the two lowest states of a well: the density sloshes from one side to the other at the Bohr frequency(E2−E1)/h — an oscillating charge, hence (Chapter 17) a radiating dipole at exactly the frequency of the spectral line.
30.5 The free particle: plane waves and wave packets
matter waves are dispersive; the phase velocity has no physical meaning, the group velocity is the particle’s velocity. A plane wave is not normalisable — a particle of perfectly sharp momentum is spread over all space — so a real particle is a wave packet, a superposition of plane waves over a band Δk, localised over Δx∼1/Δk, moving at vg and spreading: a Gaussian packet of initial width σ0 has at time t the width
σ(t)=σ01+(2mσ02ℏt)2,τ=ℏ2mσ02,
and satisfies ΔxΔp=ℏ/2 at t=0 (the minimum allowed) and more later.
Proof. The dispersion relation follows from the equation; vg as in Chapter 8. The spreading formula comes from the Fourier synthesis of the packet — each component eikx picks up the phase e−iℏk2t/2m, and the Gaussian integral with a complex coefficient gives a Gaussian of growing width; the Fourier integral is the tool of the Year 3 mathematics volume and we admit the result. Physically: the packet contains momenta spread over Δp=ℏ/2σ0, hence velocities spread over Δv=ℏ/2mσ0, and after t the faster components have outrun the slower by Δvt — the formula says exactly this. ∎
Left: the free particle’s dispersion curve, with the group velocity (tangent) and the phase velocity (chord to the origin), half of it. Right: a Gaussian packet moving at vg and spreading — its width 2 times larger after τ=2mσ02/ℏ, its peak lower, its area constant.
Example 30.9(Who spreads)
An electron localised within σ0=1nm: τ=2mσ02/ℏ=1.7×10−14s — after one nanosecond its packet is 60µm wide: an electron left alone does not stay put; bound in an atom it does not spread because the potential holds it (a stationary state). A dust grain of 1×10−15kg localised to 1µm: τ=2×107s, eight months; a bullet: longer than the age of the universe. The quantum spreading is real, and invisible for anything one can see.
Remark 30.10(Energy and time)
A state that exists only during Δt — an excited atom of lifetime τ, a wave train of length cτ — is a superposition of energies spread over ΔE∼ℏ/Δt: the Fourier relation of Chapter 18 between the duration of a train and its spectral width, now applied to E=ℏω. An atomic level of lifetime 10ns has a width ΔE=7×10−8eV, a natural linewidth of 16MHz; the Doppler width of Chapter 23 (1.5GHz) hides it in a gas, but not in a cold atom or an ion at rest.
Method 30.11(Working with wave functions)
(1) Normalise; probabilities are ∫∣ψ∣2 over the region. (2) Time-independent V: look for φ(x)e−iEt/ℏ, solve the stationary equation, impose the boundary conditions — the levels. (3) Superpose stationary states for the evolution; two levels beat at (E2−E1)/h. (4) Free motion: vg=p/m, spreading time 2mσ02/ℏ. (5) Currents: j=(ℏ/m)Im(ψ∗ψ′); for Aeikx+Be−ikx, j=(∣A∣2−∣B∣2)ℏk/m (no cross term). (6) Classical limit: Ehrenfest, and λ against the size of the apparatus.
30.6 Exercises
Exercise 30.1★
De Broglie wavelengths: an electron at 100eV, 10keV, 100keV (use p=2mE, then say why the last needs a correction); a thermal neutron (25meV, m=1.67×10−27kg); a C60 molecule (m=1.2×10−24kg) at 200m/s; a tennis ball. Which of these show interference?
Solution
Solution of Exercise 30.1.
h/2mE: 0.123nm, 12.3pm, 3.9pm (at 100keV the kinetic energy is a fifth of mc2; the relativistic momentum gives 3.7pm); neutron 0.18nm; C602.8pm; tennis ball ∼10−34m. All but the ball have shown interference.
Exercise 30.2★
ψ(x)=Ae−x2/2a2: A (use ∫e−u2du=π); ⟨x⟩, Δx; probability of finding the particle in ∣x∣<a (erf(1)=0.84) and in ∣x∣<2a (0.9995).
Solution
Solution of Exercise 30.2.
A=(πa2)−1/4; ⟨x⟩=0, Δx=a/2; 0.84; 0.9995.
Exercise 30.3★
Check that φ(x)e−iEt/ℏ solves the equation when φ solves the stationary one. Frequency of the phasor for E=1eV; for a superposition of two levels 2eV apart, frequency and wavelength of the oscillation of ∣ψ∣2 — and of the light emitted.
Solution
Solution of Exercise 30.3.
Substitution gives Eφ=−(ℏ2/2m)φ′′+Vφ. E/h=2.4×1014Hz; 4.8×1014Hz, 620nm: the red light the transition emits.
Exercise 30.4★
Spreading time τ=2mσ02/ℏ and width after 1ns and 1s: an electron with σ0=1nm and σ0=1mm; a proton with σ0=1nm; a grain of 1×10−12kg with σ0=1µm. Why do objects stay where they are put?
Solution
Solution of Exercise 30.4.
Electron, 1nm: τ=1.7×10−14s, 60µm after 1ns, 60km after 1s; 1mm: τ=17s, unchanged after 1ns, +0.2% after 1s; proton, 1nm: τ=3.2×10−11s, 30nm, 30m; grain: τ=600 years. For visible objects τ is astronomical — and their surroundings keep localising them.
Exercise 30.5★★
Current. (a) Derive the continuity equation. (b) j for Aeikx, for Asinkx, for Aeikx+Be−ikx. (c) Why is the absence of a cross term in the last case the statement that reflected and incident fluxes simply subtract? (d) Show that for a real wave functionj=0: a bound stationary state carries no current.
Solution
Solution of Exercise 30.5.
(a) Proposition 30.5. (b) ∣A∣2ℏk/m; 0; (∣A∣2−∣B∣2)ℏk/m. (c) The incident and reflected fluxes subtract with no interference term: transmission and reflection coefficients can be defined from fluxes. (d) ψ∗ψ′ real, its imaginary part zero.
Exercise 30.6★★
A Gaussian packet.ψ(x,0)=Aeik0xe−x2/4σ2. (a) A, ⟨x⟩, Δx. (b) ⟨p⟩ with p^=−iℏ∂x. (c) Admitting Δp=ℏ/2σ, check ΔxΔp=ℏ/2. (d) Prove d⟨x⟩/dt=⟨p⟩/m from the continuity equation.
Electrons in a tube. Electrons accelerated through 100V cross 30cm. (a) Speed, wavelength, time of flight. (b) Phase and group velocities; which is the electron’s? (c) A packet initially 10µm wide: spreading time, width at arrival. (d) Diffraction by a 0.3mm aperture on the way: angular spread and spot on the screen. Is the ray picture valid?
Solution
Solution of Exercise 30.7.
(a) 5.9×106m/s, 0.123nm, 51ns. (b) vφ=3.0×106m/s, vg=5.9×106m/s: the latter. (c) τ=1.7µs≫51ns: unchanged. (d) λ/a=4×10−7rad, 0.1µm: rays are fine.
Exercise 30.8★★
Phases. (a) Show that ψ and eiθψ give the same predictions. (b) For ψ=(ψ1+eiαψ2)/2, write ∣ψ∣2 and show that α matters. (c) In the double slit, what are ψ1 and ψ2, and what fixes their relative phase at a point of the screen? (d) Why does detecting which slit the electron passed destroy the fringes (qualitative)?
Solution
Solution of Exercise 30.8.
(a) ∣ψ∣2 and every expectation value are unchanged. (b) 21(∣ψ1∣2+∣ψ2∣2)+Re(eiαψ1∗ψ2). (c) The waves from the two slits; their relative phase is k(r2−r1). (d) The detector ends in different states for the two paths; the two components no longer interfere, and the pattern is the sum of two single-slit patterns.
Exercise 30.9★★
Ehrenfest. (a) For V=21mω2x2, show that ⟨x⟩ obeys exactly ⟨x⟩¨=−ω2⟨x⟩. (b) For V=mgx, ⟨x⟩ falls like a stone. (c) For a general V, when is ⟨V′(x)⟩≈V′(⟨x⟩)? (d) A 1g bead in a bowl: estimate the size of its packet needed for it to be “quantum” and compare with an atom.
Solution
Solution of Exercise 30.9.
(a) ⟨V′⟩=mω2⟨x⟩: exact. (b) ⟨V′⟩=mg: exact. (c) When the packet is narrow against the scale on which V′′ varies. (d) Its wavelength at 0.1m/s is 10−29m; the packet would have to be as wide as the bowl to matter. An atom’s thermal wavelength, 0.1nm, is its own size.
Exercise 30.10★★★
Sloshing. In a box 0<x<L with impenetrable walls the stationary states are φn=2/Lsin(nπx/L), En=n2h2/8mL2 (Chapter 31). Take ψ=(φ1e−iE1t/ℏ+φ2e−iE2t/ℏ)/2. (a) ∣ψ∣2 as a function of time. (b) Show ⟨x⟩=L/2−(16L/9π2)cosω21t (use ∫0Lxsin(πx/L)sin(2πx/L)dx=−8L2/9π2). (c) Electron, L=1nm: E2−E1, the frequency and wavelength of the oscillation. (d) What does an oscillating charge do (Chapter 17)? Conclude on spectral lines.
Solution
Solution of Exercise 30.10.
(a) 21(φ12+φ22)+φ1φ2cosω21t. (b) L/2+cosω21t∫xφ1φ2=L/2−(16L/9π2)cosω21t. (c) 3h2/8mL2=1.13eV, 2.7×1014Hz, 1.1µm. (d) It radiates at ν21: the line; a pure stationary state has no oscillating dipole and does not radiate.
Exercise 30.11★★★
Spreading, derived. The packet ψ(x,0)∝e−x2/4σ02 is the superposition ∫g(k)eikxdk with g(k)∝e−σ02k2 (admitted: the Fourier transform of a Gaussian is a Gaussian, width 1/2σ0 in k). (a) Write ψ(x,t) by giving each component its phase e−iℏk2t/2m. (b) Complete the square in the exponent and admit ∫e−αk2+βkdk=π/αeβ2/4α for complex α with positive real part: show ∣ψ∣2 is a Gaussian of width σ(t)=σ01+(ℏt/2mσ02)2. (c) Interpret τ via Δv=Δp/m. (d) Why does a packet of photons in vacuum not spread?
Solution
Solution of Exercise 30.11.
(a) ψ∝∫e−σ02k2ei(kx−ℏk2t/2m)dk. (b) α=σ02+iℏt/2m, β=ix: ∣ψ∣2∝exp[−x2/2σ02(1+(ℏt/2mσ02)2)]. (c) Δv=ℏ/2mσ0, and Δvτ=σ0. (d) ω=ck: no dispersion, all components at c.
Exercise 30.12★★★
Energy–time. (a) A state decaying as e−t/2τe−iE0t/ℏ (probability e−t/τ) is a superposition of energies: admitting that its energy distribution is a Lorentzian of full width Γ=ℏ/τ, give Γ for τ=10ns (atom), 1×10−8s, 1×10−23s (a hadronic resonance, in MeV). (b) Natural linewidth of the 10ns level in frequency; compare with the Doppler width 1.5GHz. (c) How does one see the natural width (cold atoms, trapped ions)? (d) Relate to the coherence length of the emitted wave train.
Solution
Solution of Exercise 30.12.
(a) Γ=ℏ/τ: 6.6×10−8eV; 6.6×10−8eV; 66MeV. (b) 1/2πτ=16MHz, a hundred times less than Doppler. (c) In cold atoms and trapped ions the Doppler width vanishes and the natural width is reached — the basis of frequency standards. (d) cτ=3m, the train’s length.
30.7 Problem: The electron in a tube and the spreading of a wave packet
Problem 30.1
Weekend problem — when is an electron a particle, when a wave
Part I — The cathode-ray tube. Electrons leave a hot cathode with negligible energy, are accelerated through U=20kV, pass a 0.3mm aperture, and hit a screen 40cm away.
Speed and momentum (non-relativistic); de Broglie wavelength.
The kinetic energy is 4% of mc2: is the non-relativistic treatment acceptable for this problem?
Time of flight to the screen.
Diffraction by the aperture: angular spread ∼λ/a and the spot it produces on the screen; compare with a pixel (0.3mm).
Heisenberg at the aperture: Δpy≳ℏ/2a gives the same spread — check.
A packet 0.3mm wide along the beam: spreading time 2mσ02/ℏ and width at the screen.
Deflection plates apply a transverse force: by Ehrenfest, how does ⟨y⟩ move? Conclude that the tube is a classical device.
Phase velocity of the electron wave; why is it no paradox that it differs from the electron’s speed?
20. At ν21: the spectral line; a stationary state has no oscillating dipole and is stable; the superposition radiates its way down to φ1.
21. Not definite: ⟨E⟩=(E1+E2)/2=0.94eV; a measurement gives E1 or E2, each with probability one half.
22. Random: where each electron lands; not random: the pattern ∣ψ∣2 that the landings build.
23.j=∣A∣2ℏk/m: with ∣A∣2=n electrons per unit length, nv electrons per second.
24. First order so that ψ(x,0) is a complete description that fixes the future; complex so that a definite energy has the single time dependence e−iEt/ℏ and E=p2/2m comes out.
25.ψ is the amplitude whose square is the probability; it obeys Schrödinger’s linear equation; stationary states carry a phasor and fixed densities; a free packet moves at p/m and spreads in 2mσ02/ℏ; Ehrenfest gives Newton back for the large.