Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Position and momentum representations

Learning objectives

After this chapter, you should be able to:

Sections 6.1–6.7 form the conceptual core. Sections 6.8–6.9 are applications to packet spreading and expectation-value dynamics; they may be taken as a second pass before the chapter’s final representation check.

6.1 From discrete labels to a continuum

Every chapter so far has worked with a spin’s discrete state space: two basis kets, a handful of amplitudes, sums instead of integrals. A particle’s position, by contrast, can take a continuum of values, and this chapter extends the same logic that worked for spin to that continuous setting. Spin states already taught us the key idea: a ket is abstract, and a column of numbers is only its representation in some chosen basis. Position simply uses a basis with a continuous label instead of a discrete one. The generalized ket ∣x⟩|x\rangle represents an ideal position outcome, and the position-space wavefunction is defined exactly as an amplitude was for spin:

ψ(x)=⟨x∣ψ⟩.\psi(x)=\langle x|\psi\rangle.

The completeness relation, which for spin was a finite sum over two basis kets, becomes an integral over the continuum of possible positions:

∫−∞∞∣x⟩⟨x∣ dx=I^,\int_{-\infty}^{\infty}|x\rangle\langle x|\,dx=\hat I,

so that

∣ψ⟩=∫−∞∞ψ(x)∣x⟩ dx.|\psi\rangle=\int_{-\infty}^{\infty}\psi(x)|x\rangle\,dx.

The basis is delta-normalized:

⟨x∣x′⟩=δ(x−x′).\langle x|x'\rangle=\delta(x-x').

The Dirac delta is not an ordinary finite function; it is defined entirely by its effect inside an integral:

∫f(x′)δ(x′−x) dx′=f(x).\int f(x')\delta(x'-x)\,dx'=f(x).

One consequence follows immediately: ideal position kets are not normalizable physical states at all. Real preparations are wave packets of finite width, as Section 6.5 develops in detail.

Optional mathematical detail: why the delta is singular

It helps to see directly why an ordinary finite function cannot play the role of ⟨x∣x′⟩\langle x|x'\rangle, rather than simply asserting it. Suppose instead that ⟨x∣x′⟩=g(x−x′)\langle x|x'\rangle=g(x-x') for some bounded function gg. Completeness would then require

⟨x′∣ψ⟩=∫g(x′−x)ψ(x) dx\langle x'|\psi\rangle=\int g(x'-x)\psi(x)\,dx

to equal ψ(x′)\psi(x') for every state ∣ψ⟩|\psi\rangle. But a bounded kernel convolves and smooths a function as it integrates against it; it cannot reproduce the function’s exact value at a single point. The only way it could would be for the kernel to be concentrated entirely at x=x′x=x' with unit weight there, and no ordinary function is simultaneously concentrated at one point and normalized to integrate to one. Only the singular limit of an infinitely narrow, infinitely tall spike—the Dirac delta—has the sifting property we actually need. This is the same logic that made a discrete basis’s Kronecker delta δmn\delta_{mn} work for spin: orthonormality is exactly what lets us read off a coefficient by taking an inner product. The continuum version of that same requirement forces delta-function normalization in place of unit inner products, and it is exactly what pushes the generalized kets ∣x⟩|x\rangle just outside the space of normalizable states.

6.2 The continuous Born rule

The probability of finding the particle between aa and bb is

P(a≤x≤b)=∫ab∣ψ(x)∣2 dx.P(a\le x\le b)=\int_a^b|\psi(x)|^2\,dx.

Normalization requires

∫−∞∞∣ψ(x)∣2 dx=1.\int_{-\infty}^{\infty}|\psi(x)|^2\,dx=1.

Unlike a discrete probability, ∣ψ(x)∣2|\psi(x)|^2 is a probability density, with units of inverse length. This has a striking consequence: the probability of finding the particle at any one exact point is zero, even at a point where the density itself is nonzero.

This is a genuine conceptual shift from the spin chapters, not merely a change of notation. There, ∣⟨+z∣ψ⟩∣2|\langle+z|\psi\rangle|^2 was already, by itself, a probability. Here one more step—multiplying by an interval dxdx, or integrating over a finite range—is required before a number between 0 and 1 appears. Forgetting this step is the most common bookkeeping error when moving from discrete to continuous quantum mechanics: a “large” value of ∣ψ(x)∣2|\psi(x)|^2 at some point says only that a narrow interval there carries appreciable probability, not that the point itself is likely.

A probability-density curve has the area between positions a and b shaded. A second graph shows that interval probability approaches zero as the detector interval shrinks to an exact point.

Figure 1:In a continuous measurement, height is a density and area is probability. Shrinking the detector interval to zero removes its area even when the density at the selected point remains nonzero.

Example 6.1: normalize an exponential packet

Let ψ(x)=Ae−κ∣x∣\psi(x)=Ae^{-\kappa|x|} with κ>0\kappa>0. Because ∣x∣|x| makes the integrand symmetric about the origin, the integral over the whole line splits into twice the integral over x>0x>0, where ∣x∣=x|x|=x and the exponential is easy to integrate directly:

1=∫−∞∞∣A∣2e−2κ∣x∣ dx=2∣A∣2∫0∞e−2κx dx=2∣A∣2[−e−2κx2κ]0∞=∣A∣2κ,1=\int_{-\infty}^\infty|A|^2e^{-2\kappa|x|}\,dx =2|A|^2\int_0^\infty e^{-2\kappa x}\,dx =2|A|^2\left[\frac{-e^{-2\kappa x}}{2\kappa}\right]_0^\infty =\frac{|A|^2}{\kappa},

so one convenient choice is A=κA=\sqrt\kappa. For ⟨x⟩\langle x\rangle, the integrand x∣ψ(x)∣2=∣A∣2xe−2κ∣x∣x|\psi(x)|^2=|A|^2xe^{-2\kappa|x|} is an odd function times an even one, hence odd overall, and an odd function integrated over a range symmetric about the origin always integrates to zero; so ⟨x⟩=0\langle x\rangle=0 without any further calculation. For ⟨x2⟩\langle x^2\rangle, the same even/odd split as before applies, and the standard integral ∫0∞x2e−2κx dx=2/(2κ)3=1/(4κ3)\int_0^\infty x^2e^{-2\kappa x}\,dx=2/(2\kappa)^3=1/(4\kappa^3) (obtained by integrating by parts twice, or by differentiating ∫0∞e−sxdx=1/s\int_0^\infty e^{-sx}dx=1/s twice with respect to ss) gives

⟨x2⟩=2κ∫0∞x2e−2κx dx=2κ⋅14κ3=12κ2.\langle x^2\rangle =2\kappa\int_0^\infty x^2e^{-2\kappa x}\,dx =2\kappa\cdot\frac{1}{4\kappa^3} =\frac{1}{2\kappa^2}.

Thus Δx=⟨x2⟩−⟨x⟩2=1/(2κ)\Delta x=\sqrt{\langle x^2\rangle-\langle x\rangle^2}=1/(\sqrt2\kappa).

Concept check 6.1

Is ∣ψ(x0)∣2|\psi(x_0)|^2 the probability that a position measurement gives exactly x0x_0?

Answer

No. It is a density. A finite probability requires integration over a finite interval, and the density carries units of inverse length.

Screenshot of the Quantum Wave Interference simulation

Figure 2:Send electrons, photons, or atoms one at a time through a double slit and watch each arrive as a single dot, yet build up the interference pattern predicted by ∣ψ(x)∣2|\psi(x)|^2 over many trials. Cover one slit with a detector and the pattern collapses to two bands, the continuous-variable version of Chapter 1’s Rule 3.

Interactive simulation: Quantum Wave Interference

6.3 Operators in the position representation

Position acts by multiplication:

⟨x∣x^∣ψ⟩=xψ(x).\langle x|\hat x|\psi\rangle=x\psi(x).

Momentum acts by differentiation:

⟨x∣p^∣ψ⟩=−iℏdψdx.\langle x|\hat p|\psi\rangle =-i\hbar\frac{d\psi}{dx}.

Apply both operator orders explicitly to an arbitrary function ψ(x)\psi(x), rather than taking the commutator below on faith. The term x^p^ ψ\hat x\hat p\,\psi first differentiates, then multiplies by xx:

x^p^ ψ(x)=x(−iℏdψdx)=−iℏ xdψdx.\hat x\hat p\,\psi(x)=x\left(-i\hbar\frac{d\psi}{dx}\right) =-i\hbar\,x\frac{d\psi}{dx}.

The term p^x^ ψ\hat p\hat x\,\psi first multiplies by xx, then differentiates the product xψ(x)x\psi(x) using the ordinary product rule:

p^x^ ψ(x)=−iℏddx(xψ(x))=−iℏ ψ(x)−iℏ xdψdx.\hat p\hat x\,\psi(x)=-i\hbar\frac{d}{dx}\big(x\psi(x)\big) =-i\hbar\,\psi(x)-i\hbar\,x\frac{d\psi}{dx}.

Subtracting, the term proportional to dψ/dxd\psi/dx—the piece generated by the derivative landing on ψ\psi itself—is identical in both orderings and cancels exactly, leaving only the extra term the product rule generates when the derivative instead lands on the factor of xx:

[x^,p^] ψ(x)=x^p^ ψ−p^x^ ψ=−iℏ xdψdx−(−iℏ ψ−iℏ xdψdx)=iℏ ψ(x).[\hat x,\hat p]\,\psi(x)=\hat x\hat p\,\psi-\hat p\hat x\,\psi =-i\hbar\,x\frac{d\psi}{dx} -\left(-i\hbar\,\psi-i\hbar\,x\frac{d\psi}{dx}\right) =i\hbar\,\psi(x).

Since ψ(x)\psi(x) was an arbitrary function, this is an operator identity, true independent of any particular state:

[x^,p^]ψ(x)=iℏψ(x),[x^,p^]=iℏI^.[\hat x,\hat p]\psi(x)=i\hbar\psi(x), \qquad [\hat x,\hat p]=i\hbar\hat I.

The expectation values are

⟨x⟩=∫ψ∗(x)xψ(x) dx,⟨p⟩=∫ψ∗(x)(−iℏddx)ψ(x) dx.\langle x\rangle=\int\psi^*(x)x\psi(x)\,dx, \qquad \langle p\rangle=\int\psi^*(x) \left(-i\hbar\frac{d}{dx}\right)\psi(x)\,dx.

Boundary behavior matters here in a way it did not for spin.

Quick calculation: integration by parts and the boundary term

Hermiticity of p^\hat p would require ⟨ϕ∣p^ψ⟩=⟨p^ϕ∣ψ⟩\langle\phi|\hat p\psi\rangle =\langle\hat p\phi|\psi\rangle for every pair of states ϕ,ψ\phi,\psi. Writing the left side as an integral and integrating by parts once,

∫−∞∞ϕ∗(x)(−iℏdψdx)dx=−iℏ[ϕ∗(x)ψ(x)]−∞∞+iℏ∫−∞∞dϕ∗dxψ(x) dx.\int_{-\infty}^{\infty}\phi^*(x)\left(-i\hbar\frac{d\psi}{dx}\right)dx =-i\hbar\Big[\phi^*(x)\psi(x)\Big]_{-\infty}^{\infty} +i\hbar\int_{-\infty}^{\infty}\frac{d\phi^*}{dx}\psi(x)\,dx.

The remaining integral on the right is exactly ⟨p^ϕ∣ψ⟩\langle\hat p\phi|\psi \rangle, the quantity symmetry of p^\hat p requires the left side to equal. Thus the two sides agree when the boundary form [ϕ∗ψ]−∞∞\big[\phi^*\psi\big]_{-\infty}^{\infty} vanishes for every pair in the operator’s domain. Square integrability alone does not guarantee a pointwise limit at infinity; one also needs suitable regularity and boundary behavior. The standard momentum domain on the whole line is chosen so these conditions hold.

An unbounded operator is therefore not fully defined by a differential expression alone: its domain is part of its definition. A further subtlety appears on a finite interval. Requiring wavefunctions to vanish at two hard walls makes the boundary form vanish, so −iℏd/dx-i\hbar d/dx is symmetric on that domain, but it is not self-adjoint there and does not define the well’s momentum observable. The operator −ℏ2d2/dx2/(2m)-\hbar^2d^2/dx^2/(2m) with the hard-wall boundary conditions does define a self-adjoint Hamiltonian. Chapter 7 uses that Hamiltonian directly.

Concept check 6.2

If ψ(x)\psi(x) is multiplied by the constant phase eiγe^{i\gamma}, do ⟨x⟩\langle x\rangle or ⟨p⟩\langle p\rangle change?

Answer

No. The phase cancels from all expectation values. A position-dependent phase, such as eik0xe^{ik_0x}, can change momentum statistics.

Example 6.2: phase carries mean momentum

That last remark is worth demonstrating explicitly, since it is easy to mistake “the phase doesn’t matter” for a blanket statement rather than one about constant phase specifically. Let ψ(x)=f(x)eik0x\psi(x)=f(x)e^{ik_0x}, where ff is real, normalized, and vanishes at infinity. Then

⟨p⟩=∫fe−ik0x[−iℏddx(feik0x)]dx=ℏk0,\langle p\rangle =\int f e^{-ik_0x}\left[-i\hbar \frac{d}{dx}(fe^{ik_0x})\right]dx =\hbar k_0,

because ∫ff′ dx=0\int f f'\,dx=0. The density ∣ψ∣2=f2|\psi|^2=f^2 does not reveal k0k_0 at all; only the spatial phase carries that information. As in the spin interferometer of Chapter 4, phase information becomes observable only through a suitable measurement.

6.4 Momentum states and Fourier transforms

Momentum eigenstates satisfy p^∣p⟩=p∣p⟩\hat p|p\rangle=p|p\rangle. In position space,

⟨x∣p⟩=12πℏeipx/ℏ.\langle x|p\rangle=\frac{1}{\sqrt{2\pi\hbar}}e^{ipx/\hbar}.

The momentum-space wavefunction is

ϕ(p)=⟨p∣ψ⟩=12πℏ∫−∞∞e−ipx/ℏψ(x) dx,\phi(p)=\langle p|\psi\rangle =\frac{1}{\sqrt{2\pi\hbar}} \int_{-\infty}^{\infty}e^{-ipx/\hbar}\psi(x)\,dx,

with inverse

ψ(x)=12πℏ∫−∞∞eipx/ℏϕ(p) dp.\psi(x)=\frac{1}{\sqrt{2\pi\hbar}} \int_{-\infty}^{\infty}e^{ipx/\hbar}\phi(p)\,dp.

These are basis changes, not dynamical evolution: exactly the same distinction Chapter 2 drew between rewriting a spin ket in a new basis and performing an actual measurement. Normalization is preserved either way: ∫∣ψ(x)∣2dx=∫∣ϕ(p)∣2dp=1\int|\psi(x)|^2dx=\int|\phi(p)|^2dp=1.

A broad packet in position generally requires a narrow range of wave numbers; a narrow position packet requires many. This Fourier fact underlies the position–momentum uncertainty relation, and separating what is and is not quantum here matters. As a mathematical statement about pairs of Fourier partners, it has no quantum content of its own: it is true of classical wave pulses too. What quantum mechanics adds is the physical claim that ϕ(p)\phi(p), not merely ψ(x)\psi(x), is a genuine probability amplitude for an equally real, independently measurable observable. A classical wave packet’s Fourier width is a statement about its mathematical shape; a quantum wavefunction’s Fourier width is simultaneously a statement about the distribution of a different physical quantity’s measured values.

Four plots compare Fourier-paired Gaussian densities. A narrow position distribution corresponds to a broad momentum distribution, while a broad position distribution corresponds to a narrow momentum distribution.

Figure 3:Position and momentum descriptions are Fourier partners. Localizing the same state more tightly in one representation necessarily spreads it across more basis states in the other.

Screenshot of the Fourier: Making Waves simulation

Figure 4:Add harmonics one at a time and watch a broad range of wave numbers synthesize a narrow, localized packet in position, or remove harmonics and watch the reverse. The spread displayed for each domain is exactly the trade-off behind equation (27) below.

Interactive simulation: Fourier: Making Waves

Concept check 6.3

Does a momentum eigenfunction describe a localized physical particle?

Answer

No. Its magnitude is constant over all space and it cannot be normalized in the ordinary sense. It is a generalized basis state; physical beams are modeled by normalizable superpositions with a range of momenta.

6.5 Gaussian wave packets

A normalized Gaussian centered at x0x_0 with mean wave number k0k_0 can be written

ψ(x)=1(2πσx2)1/4exp⁡[−(x−x0)24σx2+ik0x].\psi(x)=\frac{1}{(2\pi\sigma_x^2)^{1/4}} \exp\left[-\frac{(x-x_0)^2}{4\sigma_x^2}+ik_0x\right].

It has

⟨x⟩=x0,Δx=σx,⟨p⟩=ℏk0,Δp=ℏ2σx.\langle x\rangle=x_0,\quad \Delta x=\sigma_x,\quad \langle p\rangle=\hbar k_0,\quad \Delta p=\frac{\hbar}{2\sigma_x}.

Thus

Δx Δp=ℏ2.\Delta x\,\Delta p=\frac{\hbar}{2}.

The Gaussian saturates the Robertson bound associated with equation (17). This uncertainty is not, however, caused by a badly calibrated position detector kicking the particle around. It describes the spreads found in two genuinely different ensembles, each prepared identically in the same state, one measured for position and the other for momentum.

Optional mathematical extension: why the Gaussian saturates the bound

For the canonical pair, equality in the Cauchy–Schwarz argument behind the uncertainty relation can occur only when the two centered state vectors are related by

(A^−⟨A⟩)∣ψ⟩=iλ(B^−⟨B⟩)∣ψ⟩(\hat A-\langle A\rangle)|\psi\rangle =i\lambda(\hat B-\langle B\rangle)|\psi\rangle

for some real number λ\lambda. It is worth solving this condition explicitly for A^=x^\hat A=\hat x, B^=p^\hat B=\hat p, to see directly why the minimizing states turn out to be Gaussian and not some other shape. In the position representation, with ⟨p⟩\langle p\rangle momentarily set to zero for simplicity,

(x−⟨x⟩)ψ(x)=iλ(−iℏdψdx)=λℏdψdx,(x-\langle x\rangle)\psi(x) =i\lambda\left(-i\hbar\frac{d\psi}{dx}\right) =\lambda\hbar\frac{d\psi}{dx},

a first-order differential equation with solution

ψ(x)∝exp⁡[(x−⟨x⟩)22λℏ].\psi(x)\propto \exp\left[\frac{(x-\langle x\rangle)^2}{2\lambda\hbar}\right].

Normalizability requires the exponent to decay, so λ\lambda must be negative; writing λ=−2σx2/ℏ\lambda=-2\sigma_x^2/\hbar reproduces exactly the Gaussian profile of equation (25) (restoring ⟨p⟩=ℏk0\langle p\rangle=\hbar k_0 simply reinserts the phase eik0xe^{ik_0x}, which affects neither Δx\Delta x nor Δp\Delta p). No other line shape solves equation (28): any state that is not Gaussian must have Δx Δp>ℏ/2\Delta x\,\Delta p>\hbar/2 strictly. So Gaussian wave packets, far from being merely a convenient textbook example, are the states of least joint position–momentum uncertainty: there is no better one to be found.

Example 6.3: overlap of two wave packets

How distinguishable are two Gaussian packets prepared with the same width σx\sigma_x and mean wave number k0k_0, but different centers x1x_1 and x2x_2? Their overlap is

⟨ψ2∣ψ1⟩=12πσx2∫−∞∞exp⁡[−(x−x1)24σx2−(x−x2)24σx2]dx,\langle\psi_2|\psi_1\rangle =\frac{1}{\sqrt{2\pi\sigma_x^2}} \int_{-\infty}^{\infty} \exp\left[-\frac{(x-x_1)^2}{4\sigma_x^2} -\frac{(x-x_2)^2}{4\sigma_x^2}\right]dx,

where the common phase eik0xe^{ik_0x} has canceled between ψ2∗\psi_2^* and ψ1\psi_1. Writing xˉ=(x1+x2)/2\bar x=(x_1+x_2)/2, the identity (x−x1)2+(x−x2)2=2(x−xˉ)2+12(x1−x2)2(x-x_1)^2+(x-x_2)^2=2(x-\bar x)^2+\tfrac12(x_1-x_2)^2 separates the exponent into a piece that integrates to the normalization and a piece that depends only on the separation:

⟨ψ2∣ψ1⟩=exp⁡[−(x1−x2)28σx2].\langle\psi_2|\psi_1\rangle =\exp\left[-\frac{(x_1-x_2)^2}{8\sigma_x^2}\right].

Packets separated by many widths (∣x1−x2∣≫σx|x_1-x_2|\gg\sigma_x) are nearly orthogonal and behave as reliably distinguishable preparations; packets separated by less than a width overlap substantially, and no measurement can sort individual particles between them with certainty. This is the continuous counterpart of the Stern–Gerlach statement that only exactly orthogonal spin states are perfectly distinguishable: here “exactly orthogonal” is replaced by “infinitely far apart,” and any finite separation always leaves some residual overlap.

Concept check 6.4

Two Gaussian packets of width σx\sigma_x are centered 6σx6\sigma_x apart. Are they reliably distinguishable by an ideal position measurement?

Answer

Yes, to excellent approximation. Equation (32) gives an overlap of exp⁡(−36/8)=exp⁡(−4.5)≈0.011\exp(-36/8)=\exp(-4.5)\approx0.011, so the two states are nearly orthogonal and a position measurement almost always lands cleanly in the region associated with the correct packet.

6.6 Momentum generates translations

Define the translation operator

T^(a)=e−iap^/ℏ.\hat T(a)=e^{-ia\hat p/\hbar}.

Using p^=−iℏ d/dx\hat p=-i\hbar\,d/dx, the exponent becomes −iap^/ℏ=−a d/dx-ia\hat p/\hbar=-a\,d/dx, so in the position representation the exponential’s defining power series—the same construction reviewed in Chapter 4 for a matrix exponential, now applied to the differentiation operator d/dxd/dx instead—gives

⟨x∣T^(a)∣ψ⟩=e−a d/dxψ(x)=∑n=0∞(−a)nn!dnψdxn=ψ(x)−a ψ′(x)+a22 ψ′′(x)−⋯ .\langle x|\hat T(a)|\psi\rangle =e^{-a\,d/dx}\psi(x) =\sum_{n=0}^\infty\frac{(-a)^n}{n!}\frac{d^n\psi}{dx^n} =\psi(x)-a\,\psi'(x)+\frac{a^2}{2}\,\psi''(x)-\cdots.

This is exactly the Taylor series of ψ(x−a)\psi(x-a) expanded about the point xx, so

⟨x∣T^(a)∣ψ⟩=ψ(x−a).\langle x|\hat T(a)|\psi\rangle=\psi(x-a).

The translated packet is centered a distance aa to the right: its value at xx equals the old packet’s value at x−ax-a. Infinitesimally,

T^(a)≈I^−iaℏp^.\hat T(a)\approx\hat I-\frac{ia}{\hbar}\hat p.

Momentum is therefore the generator of spatial translations, just as the Hamiltonian is the generator of time translations and spin components generate rotations. The pattern is always the same: a Hermitian observable exponentiates into a unitary operator that implements a continuous symmetry, and the observable’s expectation value is conserved exactly when that symmetry leaves the Hamiltonian unchanged. Chapter 9 will show that orbital and spin angular momentum generate rotations in precisely this sense, so that the rotation operators already used for spin in Chapter 4 are a special case of a single unifying structure that also produces (33) here.

Conjugating position by the translation gives

T^†(a)x^T^(a)=x^+aI^.\hat T^\dagger(a)\hat x\hat T(a)=\hat x+a\hat I.

This operator statement makes the geometry independent of any particular wavefunction: it is true for every state at once, not just the ones we happen to write down explicitly.

Concept check 6.5

Does translating a state by aa change the magnitude of its momentum-space wavefunction?

Answer

No. Since T^(a)=e−iap^/ℏ\hat T(a)=e^{-ia\hat p/\hbar} commutes with p^\hat p, translation multiplies ϕ(p)\phi(p) by the pure phase e−iap/ℏe^{-iap/\hbar} and leaves ∣ϕ(p)∣2|\phi(p)|^2 unchanged. Sliding a packet in space redistributes nothing among its momentum components; only where it is found changes, not how fast it is found to be moving.

6.7 Measurement resolution and wave packets

An ideal sharp position measurement uses projectors ∣x⟩⟨x∣|x\rangle\langle x| only as a mathematical limit. A detector pixel covering a region RR corresponds instead to

P^R=∫R∣x⟩⟨x∣ dx,P(R)=⟨ψ∣P^R∣ψ⟩.\hat P_R=\int_R|x\rangle\langle x|\,dx, \qquad P(R)=\langle\psi|\hat P_R|\psi\rangle.

Selecting that detector outcome produces a state localized to the region, after normalization. A narrower region generally broadens the momentum distribution. This is the continuous version of an analyzer preparing an outcome state, with finite detector resolution replacing an ideal one-dimensional projector.

6.8 Application: free-particle evolution and packet spreading

For V=0V=0, each momentum component is an energy eigenstate with E(p)=p2/(2m)E(p)=p^2/(2m). The momentum wavefunction evolves as

ϕ(p,t)=e−ip2t/(2mℏ)ϕ(p,0).\phi(p,t)=e^{-ip^2t/(2m\hbar)}\phi(p,0).

Its magnitude is unchanged, so the momentum distribution stays fixed for all time. The phase, however, is nonlinear in pp, and transforming back to position shows that this nonlinear phase causes a packet to spread. For the minimum-uncertainty Gaussian of equation (25),

Δx(t)=σx1+(ℏt2mσx2)2.\Delta x(t) =\sigma_x\sqrt{1+ \left(\frac{\hbar t}{2m\sigma_x^2}\right)^2}.

The center, meanwhile, moves classically:

⟨x⟩t=x0+ℏk0mt.\langle x\rangle_t=x_0+\frac{\hbar k_0}{m}t.

A narrow initial packet spreads faster than a broad one, because a narrow packet in position necessarily contains a broader range of momenta—and therefore a broader range of velocities—than a broad one does.

Three Gaussian probability densities at successive times move to the right while becoming wider and lower. Their centers travel uniformly but their spatial uncertainty increases.

Figure 5:Free evolution advances the packet center at the group velocity while its constituent momentum components separate. The fixed momentum spread therefore appears as a growing position spread.

Group and phase velocity

A plane-wave component has phase velocity vph=ω/k=ℏk/(2m)v_{\mathrm{ph}}=\omega/k=\hbar k/(2m), while the packet envelope itself travels at the different speed

vg=dωdk=ℏkm=pm.v_{\mathrm g}=\frac{d\omega}{dk}=\frac{\hbar k}{m}=\frac{p}{m}.

It is the group velocity that matches the particle’s classical velocity. Confusing the two gives a factor-of-two error for the nonrelativistic dispersion relation ω=ℏk2/(2m)\omega=\hbar k^2/(2m). The phase velocity is not meaningless—it governs how the fine ripples inside the envelope move—but it is not what a particle detector downstream registers as “how fast the particle arrived.” For a narrow-band packet, the carrier oscillations do not appear separately in ∣ψ∣2|\psi|^2, whereas the envelope locates the probability packet. Its group velocity is therefore the speed to compare with the packet’s classical motion.

6.9 Application: Ehrenfest’s theorem

For

H^=p^22m+V(x^),\hat H=\frac{\hat p^2}{2m}+V(\hat x),

the expectation-value equation from Chapter 4 gives

d⟨x⟩dt=⟨p⟩m,d⟨p⟩dt=−⟨dVdx⟩.\frac{d\langle x\rangle}{dt}=\frac{\langle p\rangle}{m}, \qquad \frac{d\langle p\rangle}{dt} =-\left\langle\frac{dV}{dx}\right\rangle.

These resemble Hamilton’s equations, and it is tempting to conclude that quantum expectation values always follow classical trajectories. They do not, in general: ⟨V′(x^)⟩≠V′(⟨x⟩)\langle V'(\hat x)\rangle\ne V'(\langle x\rangle) except in special cases. The center of a narrow packet in a slowly varying potential only approximately follows a classical trajectory.

There is one important exception, and the reason it arises deserves a closer look. For quadratic potentials, the equations close exactly for the means. The reason is purely algebraic: the discrepancy ⟨V′(x^)⟩−V′(⟨x⟩)\langle V'(\hat x)\rangle-V'(\langle x\rangle) is controlled by the higher moments of x^−⟨x⟩\hat x-\langle x\rangle, through a Taylor expansion of V′V' about ⟨x⟩\langle x\rangle. Every term in that expansion beyond the first involves at least ⟨(x^−⟨x⟩)2⟩\langle(\hat x-\langle x\rangle)^2\rangle, or a higher power still. If V′V' is linear in xx—equivalently, if VV is at most quadratic—that Taylor expansion has no higher terms at all, and the mean-value equations close exactly, regardless of how wide the packet is. For any other potential, this gap is generally nonzero and grows with the packet’s spread.

Ehrenfest’s theorem does not turn a quantum state into a classical point. The packet retains its uncertainties, may spread, and may even split into reflected and transmitted pieces. Classical behavior emerges only for suitable states, observables, potentials, and time scales, not as a general guarantee.

Example 6.4: constant force

For V(x)=−FxV(x)=-Fx, equation (44) gives

⟨p⟩t=⟨p⟩0+Ft,⟨x⟩t=⟨x⟩0+⟨p⟩0mt+F2mt2.\langle p\rangle_t=\langle p\rangle_0+Ft,\qquad \langle x\rangle_t=\langle x\rangle_0+ \frac{\langle p\rangle_0}{m}t+\frac{F}{2m}t^2.

The mean follows the classical constant-acceleration trajectory exactly, even though the position distribution need not remain narrow.

Example 6.5: the harmonic oscillator

For V(x)=12mω2x2V(x)=\tfrac12m\omega^2x^2, studied in full in Chapter 7, V′(x)=mω2xV'(x)=m\omega^2x is linear, so equation (44) closes exactly:

d⟨x⟩dt=⟨p⟩m,d⟨p⟩dt=−mω2⟨x⟩.\frac{d\langle x\rangle}{dt}=\frac{\langle p\rangle}{m}, \qquad \frac{d\langle p\rangle}{dt}=-m\omega^2\langle x\rangle.

Differentiating the first equation and substituting the second gives

d2⟨x⟩dt2=−ω2⟨x⟩,\frac{d^2\langle x\rangle}{dt^2}=-\omega^2\langle x\rangle,

the classical simple-harmonic equation of motion, satisfied exactly by ⟨x⟩t\langle x\rangle_t and ⟨p⟩t\langle p\rangle_t for any initial state—a narrow packet, a broad one, or a superposition of many energy eigenstates. Together with Example 6.4, this shows that exactly classical mean motion is not a generic feature of quantum dynamics; it is a special consequence of having a potential no steeper than quadratic. The width of the packet is a separate question entirely: for the free particle it grows without bound (equation (40)), while a harmonic packet prepared in the right shape returns periodically to its initial width, oscillating in step with ⟨x⟩t\langle x\rangle_t.

Concept check 6.6

For the harmonic oscillator, does Ehrenfest’s theorem guarantee that a widely spread-out state behaves like a classical particle in every respect?

Answer

No. It guarantees only that the means ⟨x⟩t\langle x\rangle_t and ⟨p⟩t\langle p\rangle_t trace the exact classical trajectory. The state can remain highly uncertain, and observables other than these first moments—such as Δx(t)\Delta x(t) itself—need not follow any classical law.

6.10 A representation is not a physical location

In momentum space the same abstract operators act differently:

⟨p∣p^∣ψ⟩=pϕ(p),⟨p∣x^∣ψ⟩=iℏdϕdp.\langle p|\hat p|\psi\rangle=p\phi(p),\qquad \langle p|\hat x|\psi\rangle=i\hbar\frac{d\phi}{dp}.

Nothing has happened to the particle when one changes from ψ(x)\psi(x) to ϕ(p)\phi(p). One has simply changed the basis used to describe the same ket. A physical position measurement, by contrast, changes conditional predictions through state update. This distinction is the continuous analog of rewriting a spin ket in the xx basis versus actually sending it through an xx analyzer.

Summary

Exercises

Exercises 1–12 build the core representation, normalization, Fourier, and translation skills. Exercises 13–17 accompany the dynamics applications. Exercises 18–34 are deeper derivations and synthesis; a first course should select from them rather than assign the entire bank at once.

  1. Normalize ψ(x)=A\psi(x)=A on 0<x<L0<x<L and zero elsewhere. Find ⟨x⟩\langle x\rangle and Δx\Delta x.

  2. For the exponential packet in Example 6.1, find the probability that ∣x∣<1/κ|x|<1/\kappa.

  3. Verify equation (17) by applying both operator orders to an arbitrary differentiable function.

  4. Show by integration by parts that ⟨p⟩\langle p\rangle is real when ψ(x)\psi(x) vanishes at infinity.

  5. Normalize xe−κxxe^{-\kappa x} on x>0x>0 and zero elsewhere.

  6. Find ⟨p⟩\langle p\rangle for a real, normalized wavefunction that vanishes at the boundaries of its domain.

  7. Verify the four moments stated for equation (25).

  8. Fourier transform a Gaussian by completing the square and confirm Δp=ℏ/(2σx)\Delta p=\hbar/(2\sigma_x).

  9. Show that T^(a)T^(b)=T^(a+b)\hat T(a)\hat T(b)=\hat T(a+b) and that T^†(a)=T^(−a)\hat T^\dagger(a) =\hat T(-a).

  10. Derive T^†x^T^=x^+aI^\hat T^\dagger\hat x\hat T=\hat x+a\hat I to first order in aa using the commutator.

  11. A detector reports only whether x<0x<0 or x≥0x\ge0. Write its two projectors and the conditional wavefunction after the x≥0x\ge0 result.

  12. Explain how two wavefunctions can have identical position densities but different momentum distributions. Give an explicit pair.

Dynamics applications

  1. Derive equation (40) by evolving the Gaussian in momentum space, or verify its short- and long-time limits.

  2. Find the time at which the width of a free Gaussian has doubled.

  3. Derive both Ehrenfest equations using commutators.

  4. Give a potential for which ⟨V′(x^)⟩=V′(⟨x⟩)\langle V'(\hat x)\rangle=V'(\langle x\rangle) for every state and one for which it is generally false.

  5. Verify the momentum-space representation of x^\hat x by differentiating the Fourier kernel.

Further extensions and synthesis

  1. Suppose ⟨x∣x′⟩=g(x−x′)\langle x|x'\rangle=g(x-x') for a bounded function gg with ∫g(x) dx=1\int g(x)\,dx=1. Show that ∫g(x′−x)ψ(x) dx\int g(x'-x)\psi(x)\,dx generally differs from ψ(x′)\psi(x') for a wavefunction that varies on the scale over which gg is spread, and explain why only the Dirac-delta limit removes this discrepancy for every ψ\psi.

  2. Using equation (32), find the separation ∣x1−x2∣|x_1-x_2|, in units of σx\sigma_x, at which ∣⟨ψ2∣ψ1⟩∣2=12|\langle\psi_2|\psi_1\rangle|^2=\tfrac12.

  3. Repeat the overlap calculation of Example 6.3 for two Gaussian packets with the same center but different mean wave numbers k1≠k2k_1\ne k_2, and show that the result depends on k1−k2k_1-k_2 and σx\sigma_x in a way that parallels equation (32).

  4. Solve equation (28) for A^=x^\hat A=\hat x, B^=p^\hat B=\hat p keeping ⟨p⟩\langle p\rangle general from the start, and confirm that the phase eik0xe^{ik_0x} appears with ⟨p⟩=ℏk0\langle p\rangle=\hbar k_0.

  5. Verify by direct substitution that the harmonic-oscillator Ehrenfest equations of Example 6.5 combine to d2⟨x⟩/dt2=−ω2⟨x⟩d^2\langle x\rangle/dt^2=-\omega^2\langle x\rangle.

  6. Identify the mathematical feature shared by the free particle and the harmonic oscillator that allows Ehrenfest’s theorem to close exactly for both, and explain why a cubic potential V(x)=cx3V(x)=cx^3 does not share it.

  7. A particle is confined with certainty to a finite interval [0,L][0,L]. Explain using the Fourier relation why its momentum-space wavefunction cannot be concentrated at a single value of pp.

  8. Verify ⟨x2⟩\langle x^2\rangle and ⟨p2⟩\langle p^2\rangle for the Gaussian packet of equation (25) by direct integration, and confirm the stated Δx\Delta x and Δp\Delta p.

  9. Two Gaussian packets share a center x0x_0 but have different widths σ1\sigma_1 and σ2\sigma_2. Compute their overlap and show that it depends only on the ratio σ1/σ2\sigma_1/\sigma_2.

  10. Show that a momentum eigenstate is unchanged in magnitude, up to an overall phase, under any translation T^(a)\hat T(a), and interpret this as the statement that momentum eigenstates carry no distinguishing spatial location.

  11. Show directly from its integral definition that P^R\hat P_R in Section 6.7 is Hermitian and idempotent (P^R2=P^R\hat P_R^2=\hat P_R).

  12. Model a detector pixel of width ε\varepsilon centered at x0x_0 by P^ε=∫x0−ε/2x0+ε/2∣x⟩⟨x∣ dx\hat P_\varepsilon=\int_{x_0-\varepsilon/2}^{x_0+\varepsilon/2} |x\rangle\langle x|\,dx. Show that P(R)/ε→∣ψ(x0)∣2P(R)/\varepsilon\to|\psi(x_0)|^2 as ε→0\varepsilon\to0, recovering the probability density from a finite-resolution measurement.

  13. Find the time at which Δx(t)=2 σx\Delta x(t)=\sqrt2\,\sigma_x for a free Gaussian packet, in terms of mm, σx\sigma_x, and ℏ\hbar.

  14. Compute the group velocity for the relativistic dispersion relation ω=c2k2+(mc2/ℏ)2\omega=\sqrt{c^2k^2+(mc^2/\hbar)^2}, and show that it reduces to p/mp/m when ℏk≪mc\hbar k\ll mc.

  15. A student measures momentum, obtains a result, and then asks for the “position-space picture” of the resulting state by Fourier transforming. Explain what has and has not been physically established about position by this point, contrasting it with an actual position measurement.

  16. A real, even wavefunction satisfies ψ(x)=ψ(−x)\psi(x)=\psi(-x). Using the Fourier relation between ψ\psi and ϕ\phi, show that ϕ(p)\phi(p) is real and even, and use that fact to show ⟨p⟩=0\langle p\rangle=0 without evaluating the momentum integral directly.

  17. Prove that T^(a)\hat T(a) is unitary directly from its exponential definition and the Hermiticity of p^\hat p.