Wave Packets

As we have seen, the wavefunction of a massive particle of momentum $p$ and energy $E$ , moving in free space along the $x$ -axis, can be written

$$\psi(x,t) = \bar{\psi} {\rm e}^{-{\rm i} (\omega t-k x)},$$ (C.41)

where $k= p/\hbar$ , $\omega = E/\hbar$ , and $\bar{\psi}$ is a complex constant. Here, $\omega$ and $k$ are linked via the matter wave dispersion relation, (C.25). Expression (C.41) represents a plane wave that propagates in the $x$ -direction with the phase velocity $v_p=\omega/k$ . However, it follows from Equation (C.26) that this phase velocity is only half of the classical velocity of a massive particle.

According to the discussion in the previous section, the most reasonable physical interpretation of the wavefunction is that $\vert\psi(x,t)\vert^{ 2} dx$ is proportional to (assuming that the wavefunction is not properly normalized) the probability of finding the particle between $x$ and $x+dx$ at time $t$ . However, the modulus squared of the wavefunction (C.41) is $\vert\bar{\psi}\vert^{ 2}$ , which is a constant that depends on neither $x$ nor $t$ . In other words, the previous wavefunction represents a particle that is equally likely to be found anywhere on the $x$ -axis at all times. Hence, the fact that this wavefunction propagates at a phase velocity that does not correspond to the classical particle velocity has no observable consequences.

How can we write the wavefunction of a particle that is localized in $x$ ? In other words, a particle that is more likely to be found at some positions on the $x$ -axis than at others. It turns out that we can achieve this goal by forming a linear combination of plane waves of different wavenumbers: that is,

$$\psi(x,t) = \int_{-\infty}^{\infty} \bar{\psi}(k) {\rm e}^{-{\rm i} (\omega t-k x)} dk.$$ (C.42)

Here, $\bar{\psi}(k)$ represents the complex amplitude of plane waves of wavenumber $k$ within this combination. In writing the previous expression, we are relying on the assumption that matter waves are superposable. In other words, it is possible to add two valid wave solutions to form a third valid wave solution. The ultimate justification for this assumption is that matter waves satisfy the linear wave equation (C.24).

There is a fundamental mathematical theorem, known as Fourier's theorem, which states that if

$$f(x) = \int_{-\infty}^{\infty} \bar{f}(k) {\rm e}^{ {\rm i} k x} dk,$$ (C.43)

then

$$\bar{f}(k) = \frac{1}{2\pi}\int_{-\infty}^\infty f(x) {\rm e}^{-{\rm i} k x} dx.$$ (C.44)

Here, $\bar{f}(k)$ is known as the Fourier transform of the function $f(x)$ . We can use Fourier's theorem to find the $k$ -space function, $\bar{\psi}(k)$ , that generates any given $x$ -space wavefunction, $\psi(x)$ , at a given time.

For instance, suppose that at $t=0$ the wavefunction of our particle takes the form

$$\psi(x,0) \propto \exp\left[{\rm i} k_0 x - \frac{(x-x_0)^{ 2}}{4 ({\mit\Delta}x)^{ 2}}\right].$$ (C.45)

Thus, the initial probability distribution for the particle's $x$ -coordinate is Gaussian (see Section 2.9):

$$\vert\psi(x,0)\vert^{ 2} \propto \exp\left[- \frac{(x-x_0)^{ 2}}{2 ({\mit\Delta}x)^{ 2}}\right].$$ (C.46)

It follows that a measurement of the particle's position is most likely to yield the value $x_0$ , and very unlikely to yield a value which differs from $x_0$ by more than $3 {\mit\Delta} x$ . Thus, Equation (C.45) is the wavefunction of a particle that is initially localized in some region of $x$ -space, centered on $x=x_0$ , whose width is of order ${\mit\Delta} x$ . This type of wavefunction is known as a wave packet.

According to Equation (C.42),

$$\psi(x,0) = \int_{-\infty}^{\infty} \bar{\psi}(k) {\rm e}^{ {\rm i} k x} dk.$$ (C.47)

Hence, we can employ Fourier's theorem to invert this expression to give

$$\bar{\psi}(k)\propto \int_{-\infty}^{\infty} \psi(x,0) {\rm e}^{-{\rm i} k x} dx.$$ (C.48)

Making use of Equation (C.45), we obtain

$$\bar{\psi}(k) \propto {\rm e}^{-{\rm i} (k-k_0) x_0}\int_{-\inf... ...i} (k-k_0) (x-x_0) - \frac{(x-x_0)^{ 2}}{4 ({\mit\Delta}x)^{ 2}}\right]dx.$$ (C.49)

Changing the variable of integration to $y=(x-x_0)/ (2 {\mit\Delta} x)$ , the previous expression reduces to

$$\bar{\psi}(k) \propto {\rm e}^{-{\rm i} k x_0 - \beta^{ 2}/4}\int_{-\infty}^{\infty} {\rm e}^{-(y-y_0)^{ 2}} dy,$$ (C.50)

where $\beta = 2 (k-k_0) {\mit\Delta}x$ and $y_0 = - {\rm i} \beta/2$ . The integral in the previous equation is now just a number, as can easily be seen by making the second change of variable $z=y-y_0$ . Hence, we deduce that

$$\bar{\psi}(k) \propto \exp\left[-{\rm i} k x_0 - \frac{(k-k_0)^{ 2}}{4 ({\mit\Delta}k)^{ 2}}\right],$$ (C.51)

where

$${\mit\Delta} k = \frac{1}{2 {\mit\Delta} x}.$$ (C.52)

If $\vert\psi(x,0)\vert^{ 2} dx$ is proportional to the probability of a measurement of the particle's position yielding a value in the range $x$ to $x+dx$ at time $t=0$ then it stands to reason that $\vert\bar{\psi}(k)\vert^{ 2} dk$ is proportional to the probability of a measurement of the particle's wavenumber yielding a value in the range $k$ to $k+dk$ . (Recall that $p = \hbar k$ , so a measurement of the particle's wavenumber, $k$ , is equivalent to a measurement of the particle's momentum, $p$ ). According to Equation (C.51),

$$\vert\bar{\psi}(k)\vert^{ 2} \propto \exp\left[- \frac{(k-k_0)^{ 2}}{2 ({\mit\Delta}k)^{ 2}}\right].$$ (C.53)

This probability distribution is a Gaussian in $k$ -space. [See Equation (C.46).] Hence, a measurement of $k$ is most likely to yield the value $k_0$ , and very unlikely to yield a value that differs from $k_0$ by more than $3 {\mit\Delta}k$ .

We have just seen that a wave packet with a Gaussian probability distribution of characteristic width ${\mit\Delta} x$ in $x$ -space [see Equation (C.46)] is equivalent to a wave packet with a Gaussian probability distribution of characteristic width ${\mit\Delta} k$ in $k$ -space [see Equation (C.53)], where

$${\mit\Delta}x {\mit\Delta} k = \frac{1}{2}.$$ (C.54)

This illustrates an important property of wave packets. Namely, in order to construct a packet that is highly localized in $x$ -space (i.e., with small ${\mit\Delta} x$ ) we need to combine plane waves with a very wide range of different $k$ -values (i.e., with large ${\mit\Delta} k$ ). Conversely, if we only combine plane waves whose wavenumbers differ by a small amount (i.e., if ${\mit\Delta} k$ is small) then the resulting wave packet is highly extended in $x$ -space (i.e., ${\mit\Delta} x$ is large).

According to standard wave theory, a wave packet made up of a superposition of plane waves that is strongly peaked around some central wavenumber $k_0$ propagates at the group velocity,

$$v_g = \frac{d\omega(k_0)}{dk},$$ (C.55)

rather than the phase velocity, $v_p = (\omega/k)_{k_0}$ , assuming that all of the constituent plane waves satisfy a dispersion relation of the form $\omega=\omega(k)$ . For the case of matter waves, the dispersion relation is (C.25). Thus, the associated group velocity is

$$v_g = \frac{\hbar k_0}{m} = \frac{p}{m},$$ (C.56)

where $p=\hbar k_0$ . This velocity is identical to the classical velocity of a (non-relativistic) massive particle. We conclude that the matter wave dispersion relation (C.25) is perfectly consistent with classical physics, as long as we recognize that particles must be identified with wave packets (which propagate at the group velocity) rather than plane waves (which propagate at the phase velocity).

Standard wave theory also implies that the spatial extent of a wave packet of initial extent $({\mit\Delta} x)_0$ grows, as the packet evolves in time, like

$${\mit\Delta} x \simeq ({\mit\Delta} x)_0 + \frac{d^{ 2}\omega(k_0)}{dk^{ 2}} \frac{t}{({\mit\Delta} x)_0},$$ (C.57)

where $k_0$ is the packet's central wavenumber. Thus, it follows from the matter wave dispersion relation, (C.25), that the width of a particle wave packet grows in time as

$${\mit\Delta} x \simeq ({\mit\Delta} x)_0 + \frac{\hbar}{m} \frac{t}{({\mit\Delta} x)_0}.$$ (C.58)

For example, if an electron wave packet is initially localized in a region of atomic dimensions (i.e., $\Delta x\sim 10^{-10} {\rm m}$ ) then the width of the packet doubles in about $10^{-16} {\rm s}$ .