
Quantum core · institution-dependent · Interactive lesson
Continuous spectra and the momentum-space wavefunction
Momentum has no normalisable eigenstates, so the tidy expansion over a discrete basis breaks down here. Trade the sum for an integral and the Kronecker delta for a Dirac one, and the expansion coefficient becomes a Fourier transform: the same state, written in momentum instead of position.
UPIV-06.07 · Lesson 7 of 9 in this unit. Position is not a completion score.Quantum core · institution-dependent. Departments differ in how much quantum mechanics they include in Physics IV or teach in a dedicated course. Follow your university’s required mathematical depth and assessment; these numbered units are not automatically optional.
Build the model
The one basis you cannot normalise, and the transform that repairs it.
The key idea
The momentum-space wavefunction φ(p) is the amplitude for momentum p — the coefficient of the δ-normalised momentum eigenfunction in the expansion of ψ, obtained as the Fourier transform of ψ(x), whose modulus squared is a probability per unit momentum.
Example. For an electron packet of width σₓ = 1.0 nm, φ(p) is a Gaussian of spread σₚ = ħ/2σₓ = 5.3 × 10⁻²⁶ kg m s⁻¹, so a momentum window of width 1.0 × 10⁻²⁷ kg m s⁻¹ at the peak holds probability |φ(0)|²Δp = 7.6 × 10⁻³.
The complete physical model and assumptions
A bound-state problem hands you a countable eigenbasis, orthonormal to a Kronecker δ, with |cₙ|² a probability. Momentum breaks that machinery. Solving −iħ du/dx = pu returns uₚ(x) ∝ e(ipx/ħ) for every real p, so the spectrum is a continuum, and |uₚ|² is the same constant everywhere — the normalisation integral diverges and no prefactor rescues it. The repair is not to abandon the basis but to weaken what normalisation means: pick the prefactor 1/√(2πħ) so that overlaps come out as δ(p − p′) instead of 1, and the eigenfunction becomes a distribution that is only ever used inside an integral. Everything else follows by replacing Σₙ with ∫dp. The coefficient cₙ becomes a function φ(p), the projection integral becomes the Fourier transform, |φ(p)|² becomes a probability density per unit momentum rather than a probability, and Parseval's theorem guarantees that the single normalisation ∫|ψ|²dx = 1 has already fixed ∫|φ|²dp = 1. The cost is that no laboratory prepares a state of sharp momentum: the price of a definite p is a state spread over all space, so what physics actually handles is a packet, and the two ways of writing that packet — ψ(x) and φ(p) — are locked together by the transform. That lock is what the uncertainty bound reports.
1Every real p is an eigenvalue, and none of them is a state
Momentum in one dimension is the operator −iħ d/dx, so the eigenvalue problem −iħ du/dx = pu integrates at once to uₚ(x) = A e(ipx/ħ). On the whole line there is no wall to impose a boundary condition. The only demand is that u stay bounded at both infinities, which forces p to be real, since a complex p makes |u| grow exponentially in one direction, and leaves every real value allowed. So the spectrum is the entire real axis, a continuum rather than a list. The trouble arrives at the next step. |uₚ(x)|² = |A|² is the same number at every point, so ∫|uₚ|²dx over the line diverges for every A ≠ 0. Read physically, the state claims the particle is exactly as likely to be found in a 1 nm window here as in a 1 nm window a light-year away, and no rescaling turns that into a total of one. This is not a defect peculiar to momentum. The same happens to the energy eigenfunctions of any unbound problem — a free particle, or a scattering state above a barrier. A continuous spectrum simply does not come with normalisable eigenfunctions.
2Delta normalisation, and the box that does the same job
Keep the eigenfunctions and weaken the condition they must satisfy. Take the overlap of two of them: ∫uₚ* uₚ′ dx = |A|²∫e(i(p′−p)x/ħ)dx = |A|² 2πħ δ(p′ − p), using the Fourier representation of the Dirac δ. Choosing |A|² = 1/2πħ makes that overlap exactly δ(p′ − p), the continuum's replacement for δₘₙ: zero for different labels, and for the same label a spike whose integral is one. The delta is a distribution, meaningful only inside an integral, and it carries the unit of 1/p, which is what forces uₚ to carry (J s)(−1/2). If distributions grate, put the system in a box of length L with periodic ends instead: momenta are then discrete, pₙ = 2πnħ/L, spaced by 2πħ/L, every eigenfunction is normalisable, and sums pass to integrals through Σₙ → (L/2πħ)∫dp. Every physical answer must come out independent of L, and letting L → ∞ returns the delta convention.
3The sum over states becomes an integral
The box shows exactly how Σₙ turns into ∫dp. Its normalised states are ψₙ = e(ipₙ x/ħ)/√L, which is √(2πħ/L) times the δ-normalised uₚ taken at p = pₙ. For a state that fits well inside the box, the projection cₙ = ∫ψₙ*ψ dx is therefore √(2πħ/L) φ(pₙ), where φ(p) = ∫uₚ*ψ dx is the same projection taken on the δ-normalised basis. The discrete probabilities are |cₙ|² = (2πħ/L)|φ(pₙ)|²: the density times the spacing. Sum them over any window and you have a Riemann sum for ∫|φ|²dp, one bin per allowed momentum. Let L grow and every |cₙ|² falls as 1/L while the number of states in the window rises as L, so the sum converges on the integral and L drops out. The expansion follows the same path, ψ = Σₙ cₙ ψₙ → ∫φ(p) uₚ(x) dp, and it works for every square-integrable ψ because the plane waves are complete: ∫uₚ(x) uₚ*(x′) dp = δ(x − x′), the continuum's version of Σₙ ψₙ(x) ψₙ*(x′) = δ(x − x′). Written out, φ(p) = (2πħ)(−1/2)∫ψ(x)e(−ipx/ħ)dx is the forward half of the transform pair derived in the position and momentum representations lesson (UPIV-05.04), and its Parseval theorem, ∫|φ|²dp = ∫|ψ|²dx = 1, is the statement that the |cₙ|² summed to one all along.
4Read |φ(p)|² as a density, and let the units enforce it
Because ∫|φ|²dp = 1 and p carries a unit, |φ(p)|² must carry the unit of 1/p, inverse kg m s⁻¹. It is a probability per unit momentum, and it becomes a probability only after an integral: P(p₁ ≤ p ≤ p₂) = ∫|φ|²dp. That the density can exceed one is no paradox but a consequence of that unit, exactly as |ψ(0)|² = κ reaches 1.0 × 10¹⁰ m⁻¹ for the exponential state of the third example, whose decay length is 0.10 nm. Take the 1 nm electron packet: σₚ = ħ/2σₓ = 5.3 × 10⁻²⁶ kg m s⁻¹, so the peak density is 1/(σₚ√(2π)) = 7.6 × 10²⁴ (kg m s⁻¹)⁻¹. Ask for the probability of a momentum inside a window of width 1.0 × 10⁻²⁷ kg m s⁻¹ at that peak and the answer is 7.6 × 10²⁴ × 1.0 × 10⁻²⁷ = 7.6 × 10⁻³. Ask for the probability of exactly p = 0 and the answer is zero, for every state — which is why a continuum is described by a density and never by a list of |cₙ|².
5What a momentum measurement can and cannot return
The measurement postulate says a reading is an eigenvalue of the observable. For momentum every real p is one, so any value where φ(p) ≠ 0 can come out; what changes is how probability attaches to it. Since no single value carries any, a real reading is always a bin: a spectrometer of resolution δp reports that p lies within δp of some p₀, with probability ∫|φ|²dp over that bin, close to |φ(p₀)|²δp when the bin is narrow. Repeated on many identical preparations, the fraction of readings in each bin, divided by δp, reproduces |φ(p)|² averaged over that bin. The spread is a property of the prepared state, not damage done by the apparatus. What a reading cannot do is leave the system in a momentum eigenstate. In the textbook projection rule — an idealised apparatus, and one interpretation of what the update means — the state afterwards is φ cut down to the reported bin and renormalised: a packet of momentum spread about δp, and so, through the transform, spread over roughly h/δp in space, 0.66 μm for δp = 1.0 × 10⁻²⁷ kg m s⁻¹. An immediate repeat returns a value inside the same bin, not the same number to every digit as a discrete level would. Sharpening the resolution towards zero sends the probability of the reading to zero and the extent of the state left behind to infinity, which is why the plane wave is a limit no measurement reaches.
6Discrete and continuous parts of one spectrum
Momentum on the line is purely continuous, but most Hamiltonians are mixed. A finite square well, or hydrogen, has a discrete set of bound levels below the threshold for escape and a continuum of unbound states above it, and neither part is complete alone. The closure relation needs both, Σₙ ψₙ(x) ψₙ*(x′) + ∫ψE(x) ψE*(x′) dE = δ(x − x′), with the continuum states δ-normalised in E and the integral also summing the two directions of travel at each E in one dimension. A normalised state then splits its probability as Σₙ |cₙ|² + ∫|c(E)|²dE = 1: an energy measurement returns a bound level with a finite probability, or finds the particle unbound with a probability that only a window of energies can carry. The same bookkeeping shows that a density belongs to its label as well as to the state. A free particle's energy E = p²/2m is reached from both p and −p, so the density per unit energy is (m/|p|)(|φ(p)|² + |φ(−p)|²) at |p| = √(2mE), which diverges at E = 0 whenever φ(0) ≠ 0 yet stays integrable. Quote momenta as k = p/ħ and the density per unit k is ħ|φ(p)|², a different number for the same state and the same outcomes. Only the probability of a window is label-free.
Key equations
Read the units, assumptions and limits before substituting.
Momentum eigenfunction
p̂ uₚ = p uₚ ⇒ uₚ(x) = e(ipx/ħ) / √(2πħ)
Every real p is an eigenvalue, so the spectrum is continuous — yet |uₚ|² = 1/2πħ everywhere, so no uₚ is a state.
p in kg m s⁻¹, x in m; the prefactor carries (J s)(−1/2), and 2πħ = h = 6.626 × 10⁻³⁴ J sDelta orthonormality
∫ uₚ*(x) uₚ′(x) dx = δ(p − p′)
The continuum's stand-in for δₘₙ. It fixes the prefactor and turns every later step into a Fourier identity.
from ∫ e(i(p′−p)x/ħ) dx = 2πħ δ(p − p′); δ(p) carries the unit of 1/p, (kg m s⁻¹)⁻¹The box that discretises the continuum
pₙ = 2πnħ/L, |cₙ|² = (2πħ/L) |φ(pₙ)|², Σₙ → (L/2πħ) ∫ dp
Each allowed momentum's probability falls as 1/L while the count in a window rises as L; only the density survives L → ∞.
periodic box of length L in m, much wider than the state; n an integer; 2πħ/L is the spacing in kg m s⁻¹Completeness of the plane waves
∫ uₚ(x) uₚ*(x′) dp = δ(x − x′)
Every square-integrable ψ can be built from plane waves, although not one plane wave is itself a state.
the continuum form of Σₙ ψₙ(x) ψₙ*(x′) = δ(x − x′); this δ carries m⁻¹The two representations
φ(p) = ∫ uₚ*(x) ψ(x) dx, ψ(x) = ∫ φ(p) uₚ(x) dp
Σₙ cₙ ψₙ has become ∫dp φ(p) uₚ — the same expansion with an integral in place of the sum.
[ψ] = m(−1/2), [φ] = (kg m s⁻¹)(−1/2); the 2πħ splits evenly between the pairBorn's rule for a continuous label
P(p₁ ≤ p ≤ p₂) = ∫ |φ(p)|² dp over [p₁, p₂], P(p = p₁) = 0
No single momentum carries a probability; windows do. Widen the window and the number can only grow.
|φ|² is a density per unit momentum, unit (kg m s⁻¹)⁻¹; over all p it integrates to 1 (Parseval, UPIV-05.04)Interactive lab
Predict first. Explore the model. Compare the result.
Explore the model
Try this. Drag L from 2 to 10 nm: the steps close onto the curve, and the shaded bin, a single allowed momentum, narrows fivefold and so carries a fifth of the probability, while the window fills with more states and its total settles on the integral. Then shrink the window towards zero width and watch its probability go with it: only an interval of momenta carries probability, never a single value.
Live results
- SPACING h/L
- 0.314 ħκ
- ONE STATE |c₀|²
- 0.200
- STATES IN WINDOW
- 7
- WINDOW ∫|φ|²dp
- 0.818
Results follow this lesson’s model and assumptions. Displayed values are rounded.
Predict → change one variable → check
Plan a change before running it. Every other input stays at the saved baseline.
Saved baseline. Box length L: 2.0 nm; Window half-width: 1.0 ħκ
This experiment is kept only while this page is open. It does not change saved course progress.
Catch the common trap
The commonest error is to read |φ(p)|² as the probability of the momentum p. It is not, and its units say so: for the 1 nm electron packet |φ|² peaks at 7.6 × 10²⁴ (kg m s⁻¹)⁻¹, a number vastly above 1 that would change again if momenta were quoted in eV/c. Only ∫|φ|²dp across a window is a probability, and for any single value of p that integral is zero. The companion error is to rescue the plane wave by normalising it — dividing by √L, or simply setting A = 1 — and then treating e(ipx/ħ) as a state some apparatus could prepare. It is a distribution, admissible only inside an integral: any box length L introduced along the way must cancel from every physical answer, and a beam of well-defined momentum is a narrow packet, never an eigenfunction.
Practice & worked examples
3 authored examples. Make an attempt, then reveal the reasoning one step at a time.
EasyExample 1
A momentum eigenfunction in one dimension satisfies −iħ du/dx = p u. Solve it, show that no constant normalises it to one, and find the constant that normalises it to a Dirac delta instead.
Sketch your answer or reasoning, then compare it with the authored solution. Your writing is not automatically graded or saved.
0 of 4 steps revealed.
- Separate and integrate: du/u = (ip/ħ)dx, so uₚ(x) = A e(ipx/ħ), and this solves the equation for every real p. The momentum spectrum is the whole real line.
- Try the usual normalisation. |uₚ(x)|² = |A|² at every x, so integrating from −L/2 to L/2 gives |A|²L, which grows without bound as L → ∞. No nonzero A can make it 1.
- Ask for delta normalisation instead: ∫uₚ*(x)uₚ′(x)dx = |A|²∫e(i(p′−p)x/ħ)dx = |A|² 2πħ δ(p′ − p).
- Setting that equal to δ(p′ − p) gives |A|² = 1/2πħ = 1/h, so A = 1/√h = 1/√(6.626 × 10⁻³⁴ J s) = 3.9 × 10¹⁶ (J s)(−1/2).
Answer. uₚ(x) = e(ipx/ħ)/√(2πħ), prefactor 3.9 × 10¹⁶ (J s)(−1/2). No choice of A normalises it to one: momentum eigenfunctions are distributions, not states.
MediumExample 2
The 1 nm electron packet of the definition example has a Gaussian momentum density centred on p = 0 with σₚ = 5.27 × 10⁻²⁶ kg m s⁻¹ (derived in the position and momentum representations lesson, UPIV-05.04), so its peak density is |φ(0)|² = 1/(σₚ√(2π)) = 7.57 × 10²⁴ (kg m s⁻¹)⁻¹. Enclose it in a periodic box of length L = 6.63 μm, chosen so that the allowed momenta are exactly 1.0 × 10⁻²⁸ kg m s⁻¹ apart. Find the probability |c₀|² of the single allowed momentum p = 0, then the probability that a measurement returns |p| ≤ 5.5 × 10⁻²⁸ kg m s⁻¹, once as a sum over box states and once as an integral. What happens to both when L is made ten times larger?
Sketch your answer or reasoning, then compare it with the authored solution. Your writing is not automatically graded or saved.
0 of 5 steps revealed.
- Spacing: pₙ = 2πnħ/L = nh/L, and h/L = 6.626 × 10⁻³⁴ J s ÷ 6.626 × 10⁻⁶ m = 1.0 × 10⁻²⁸ kg m s⁻¹. The packet, about 1 nm across, sits deep inside a box thousands of times wider, so cutting its tails off at the walls changes nothing measurable.
- Project: the normalised box states are ψₙ = e(ipₙ x/ħ)/√L, which is √(h/L) times uₚ at p = pₙ, so cₙ = ∫ψₙ*ψ dx = √(h/L) φ(pₙ) and |cₙ|² = (h/L)|φ(pₙ)|². At p = 0 that is 1.0 × 10⁻²⁸ × 7.57 × 10²⁴ = 7.6 × 10⁻⁴.
- Sum: the window |p| ≤ 5.5 × 10⁻²⁸ kg m s⁻¹ holds the eleven states n = −5 to +5. Across it the Gaussian falls by a factor exp(−(5.5 × 10⁻²⁸ ÷ 5.27 × 10⁻²⁶)²/2) = 0.99995, so every one of the eleven carries 7.57 × 10⁻⁴ and the sum is 11 × 7.57 × 10⁻⁴ = 8.3 × 10⁻³.
- Integral: the density is flat across the window, so ∫|φ|²dp = 7.57 × 10²⁴ × 1.1 × 10⁻²⁷ = 8.3 × 10⁻³. The two agree because each box state stands for a bin one spacing wide, and eleven such bins tile the window exactly: the sum is a Riemann sum for the integral.
- Now make L ten times larger. The spacing falls to 1.0 × 10⁻²⁹ kg m s⁻¹, so |c₀|² falls tenfold to 7.6 × 10⁻⁵, while the window now spans 1.1 × 10⁻²⁷ ÷ 1.0 × 10⁻²⁹ = 110 spacings. The window probability, 110 × 7.57 × 10⁻⁵ = 8.3 × 10⁻³, has not moved. As L → ∞ every single |cₙ|² goes to zero and only the density per unit momentum survives.
Answer. |c₀|² = 7.6 × 10⁻⁴ for L = 6.63 μm; P(|p| ≤ 5.5 × 10⁻²⁸ kg m s⁻¹) = 8.3 × 10⁻³ both as a sum over eleven box states and as ∫|φ|²dp. At 10L each |cₙ|² is ten times smaller and ten times as many share the window, so the probability is unchanged: L drops out, as every physical answer must.
HardExample 3
An electron is prepared in ψ(x) = √κ e(−κ|x|), with 1/κ = 0.10 nm. Project it on the momentum eigenfunctions to find φ(p) and confirm that |φ|² is normalised; find the probability that a momentum measurement returns |p| ≤ ħκ; find how often a spectrometer of resolution 1.0 × 10⁻²⁷ kg m s⁻¹ reports the bin centred on p = 0; and decide whether ⟨p²⟩ is finite despite the kink at x = 0. Take ħ = 1.055 × 10⁻³⁴ J s and mₑ = 9.11 × 10⁻³¹ kg.
Sketch your answer or reasoning, then compare it with the authored solution. Your writing is not automatically graded or saved.
0 of 5 steps revealed.
- Project: φ(p) = ∫uₚ*ψ dx = √κ (2πħ)(−1/2) ∫e(−κ|x|) e(−ipx/ħ) dx. The odd sine part integrates to zero, leaving 2∫e(−κx) cos(px/ħ) dx from 0 to ∞ = 2κ/(κ² + p²/ħ²), so φ(p) = √(2/π) (ħκ)³⁄²/(p² + ħ²κ²).
- Normalise: |φ|² = (2/π)(ħκ)³/(p² + ħ²κ²)². Substitute p = ħκ tan θ, so dp = ħκ sec²θ dθ and |φ|²dp becomes (2/π)cos²θ dθ; over −π/2 to π/2 that integrates to (2/π)(π/2) = 1. The projection carried the normalisation across, as Parseval promises.
- Window: |p| ≤ ħκ is |θ| ≤ π/4, so P = (2/π)∫cos²θ dθ from −π/4 to π/4 = (2/π)(π/4 + 1/2) = 1/2 + 1/π = 0.818, where ħκ = 1.055 × 10⁻³⁴ J s × 1.0 × 10¹⁰ m⁻¹ = 1.055 × 10⁻²⁴ kg m s⁻¹.
- Bin: the peak density is |φ(0)|² = 2/(πħκ) = 2 ÷ (π × 1.055 × 10⁻²⁴) = 6.03 × 10²³ (kg m s⁻¹)⁻¹. The bin is a thousandth of ħκ wide, so the density is flat across it and P ≈ 6.03 × 10²³ × 1.0 × 10⁻²⁷ = 6.0 × 10⁻⁴. A reading of exactly p = 0, with no bin at all, has probability zero.
- Tail: far out |φ|² falls as 1/p⁴, so p²|φ|² falls as 1/p² and ⟨p²⟩ = (2/π)(ħκ)³∫p² dp/(p² + ħ²κ²)² = (2/π)(ħκ)³ × π/(2ħκ) = ħ²κ², which is finite. The kinetic energy is ħ²κ²/2mₑ = (1.055 × 10⁻²⁴)² ÷ (2 × 9.11 × 10⁻³¹) = 6.1 × 10⁻¹⁹ J = 3.8 eV. A kink costs finite energy; only a step in ψ itself, whose |φ|² falls as slowly as 1/p², would make ⟨p²⟩ diverge.
Answer. φ(p) = √(2/π)(ħκ)³⁄²/(p² + ħ²κ²), normalised; P(|p| ≤ ħκ) = 1/2 + 1/π = 0.818 with ħκ = 1.055 × 10⁻²⁴ kg m s⁻¹; the central 1.0 × 10⁻²⁷ kg m s⁻¹ bin is reported with probability 6.0 × 10⁻⁴, and exactly p = 0 never; ⟨p²⟩ = ħ²κ² is finite, giving a kinetic energy of 6.1 × 10⁻¹⁹ J = 3.8 eV.
Quick check
Select an answer and explain it before checking.
Review & continue
The momentum-space wavefunction φ(p) is the amplitude for momentum p — the coefficient of the δ-normalised momentum eigenfunction in the expansion of ψ, obtained as the Fourier transform of ψ(x), whose modulus squared is a probability per unit momentum.
- Explain the relationship without looking at the equation.
- Check a sign, a unit and a model assumption in your working.
- Return to the lab or worked example if a step is still unclear.