The Free Gaussian Wave Packet and Its Spreading
~6 min read · Tier toggle changes depth (T1 plain · T2 standard · T3 rigour)
You should be comfortable with the Fourier transform relationship between position and momentum space, and you should understand that a free particle’s wave function evolves by multiplying its momentum-space amplitude by a simple phase factor.
Why this matters
A free particle does not stay localised; it spreads, and the Gaussian packet is the cleanest way to see exactly how fast and why.
In classical mechanics, a point particle with definite momentum travels in a straight line forever. In quantum mechanics, a particle with a definite momentum is delocalised over all space. To describe a particle that is *somewhere* but not *everywhere*, we must superpose many momentum states. The Gaussian wave packet is the standard choice because it keeps its Gaussian shape under free evolution, making the mathematics tractable while capturing the essential physics: the packet spreads in time. This spreading is not a failure of the theory; it is a direct consequence of the uncertainty principle and the dispersion relation for free particles. Understanding it is the bridge between the abstract Hilbert space formalism and the physical intuition of a particle moving through space.
The Idea in Plain Language
Imagine you prepare a particle in a very narrow region of space. To confine it to that small region, you must combine many different momentum components, because a single momentum state is spread out infinitely. The narrower the initial position packet, the wider the range of momenta you need. Once released, each momentum component travels at its own group velocity. Since different momenta correspond to different velocities, the components drift apart. The packet broadens. This is not diffusion due to collisions or thermal effects; it is a purely quantum-mechanical consequence of the superposition of plane waves with different wavelengths. The Gaussian shape is special because it minimises the product of position and momentum uncertainties, so it spreads in the most “efficient” way possible, and its width at any time can be calculated exactly.
Initial State and Momentum Distribution
We begin with a normalised Gaussian wave function in position space at \( t = 0 \):
Here \( \sigma_0 \) is the initial standard deviation of the position distribution. The corresponding momentum-space wave function is obtained via the Fourier transform. Because the Fourier transform of a Gaussian is another Gaussian, the momentum distribution is also Gaussian, with standard deviation \( \Delta p = \hbar / (2\sigma_0) \). This satisfies the minimum uncertainty relation \( \Delta x \, \Delta p = \hbar/2 \). The key point is that the initial localisation \( \sigma_0 \) directly determines the spread in momenta: tighter position means broader momentum, and vice versa.
Time Evolution and the Spreading Law
For a free particle, the time evolution of the momentum-space amplitude is simply a phase rotation: \( \tilde{\psi}(p, t) = \tilde{\psi}(p, 0) \exp(-i E_p t / \hbar) \), where \( E_p = p^2 / 2m \). Transforming back to position space, the wave function at time \( t \) remains Gaussian, but its width increases. The result, derived in the cited card D-012, is:
This formula tells us that at \( t = 0 \), the width is \( \sigma_0 \), and as \( t \to \infty \), the width grows linearly with time: \( \sigma(t) \approx \frac{\hbar t}{2m\sigma_0} \). The spreading is quadratic in time for small \( t \), but becomes linear for large \( t \). The characteristic spreading time is defined as \( \tau = 2m\sigma_0^2 / \hbar \). After one such time unit, the width has increased by a factor of \( \sqrt{2} \). This timescale depends only on the mass and the initial width, not on the particle’s average momentum or energy.
Physical Interpretation and Limits
The spreading can be understood as the dephasing of the momentum components. Each momentum component \( p \) travels at the group velocity \( p/m \) (its phase velocity is only half that). The spread in velocities is \( \Delta v = \Delta p / m = \hbar / (2m\sigma_0) \). Over time \( t \), the components separate by roughly \( \Delta v \, t \), which matches the asymptotic spreading rate \( \hbar t / (2m\sigma_0) \). In the classical limit, where \( \sigma_0 \) is large compared to the de Broglie wavelength, the spreading is negligible over typical observation times. For macroscopic objects, \( m \) is so large that \( \tau \) is enormous, and the packet effectively behaves like a classical particle. For microscopic particles, such as electrons, \( \tau \) can be very short, and spreading is a dominant effect in experiments like electron diffraction and atomic interferometry.
The derivation
Free Gaussian Wave Packet and Its Spreading
A free particle initially localised to width σ0 spreads, and the spreading timescale is set by the mass and the initial width alone.
A common misconception is that the spreading is caused by the uncertainty principle “forcing” the particle to move. In fact, the uncertainty principle sets the initial momentum spread, but the spreading itself is a dynamical consequence of the free-particle dispersion relation. If the particle were in a harmonic potential, the Gaussian would spread and then re-focus periodically. A free packet that starts at its minimum width never returns to it. That is not irreversibility in any deep sense: the evolution is unitary, and a packet prepared with the opposite "chirp" (a phase that makes its slow components lead) first narrows down to \( \sigma_0 \) and only then spreads. C
Worked example
Problem. An electron is prepared in a Gaussian wave packet with initial width \( \sigma_0 = 1.0 \times 10^{-10} \ \mathrm{m} \). Calculate the time \( t \) at which the packet width has doubled, i.e. \( \sigma(t) = 2\sigma_0 \). Use \( m_e = 9.11 \times 10^{-31} \ \mathrm{kg} \) and \( \hbar = 1.055 \times 10^{-34} \ \mathrm{J \cdot s} \).
- Set up the equation: \( 2\sigma_0 = \sigma_0 \sqrt{ 1 + \left( \frac{\hbar t}{2m\sigma_0^2} \right)^2 } \).
- Divide both sides by \( \sigma_0 \): \( 2 = \sqrt{ 1 + \left( \frac{\hbar t}{2m\sigma_0^2} \right)^2 } \).
- Square both sides: \( 4 = 1 + \left( \frac{\hbar t}{2m\sigma_0^2} \right)^2 \), so \( \left( \frac{\hbar t}{2m\sigma_0^2} \right)^2 = 3 \).
- Solve for \( t \): \( \frac{\hbar t}{2m\sigma_0^2} = \sqrt{3} \), hence \( t = \frac{2m\sigma_0^2 \sqrt{3}}{\hbar} \).
- Substitute values: \( t = \frac{2 \times 9.11 \times 10^{-31} \times (1.0 \times 10^{-10})^2 \times \sqrt{3}}{1.055 \times 10^{-34}} \ \mathrm{s} \).
- Calculate numerator: \( 2 \times 9.11 \times 10^{-31} \times 1.0 \times 10^{-20} \times 1.732 = 3.155 \times 10^{-50} \ \mathrm{kg \cdot m^2} \).
- Divide by \( \hbar \): \( t = \frac{3.155 \times 10^{-50}}{1.055 \times 10^{-34}} = 2.99 \times 10^{-16} \ \mathrm{s} \).
Units check. The expression \( t = 2m\sigma_0^2 \sqrt{3} / \hbar \) has units \( \mathrm{kg \cdot m^2 / (J \cdot s)} = \mathrm{kg \cdot m^2 / (kg \cdot m^2 / s^2 \cdot s)} = \mathrm{s} \), which is correct.
Three quick tests
- What is the characteristic spreading time \( \tau \) for a Gaussian packet of width \( \sigma_0 \)?
Answer
\( \tau = 2m\sigma_0^2 / \hbar \). - Why does the Gaussian packet spread?
Answer
Because it is a superposition of momentum components with different group velocities, which dephase over time. - In the classical limit, why is spreading negligible?
Answer
Because the mass \( m \) is large, making \( \tau \) very large, so the width changes slowly compared to observation times.
Problem set for PU-202 →
Each derivation in this unit ends with five graded problems with full solutions.
Lecture notes drafted with a local model and checked line by line by the editor; every number in the worked example was recomputed.