In physics, a wave packet (or wave train) is a short "burst" or "envelope" of localized wave action that travels as a unit. A wave packet can be analyzed into, or can be synthesized from, an infinite set of component sinusoidal waves of different wavenumbers, with phases and amplitudes such that they interfere constructively only over a small region of space, and destructively elsewhere. Depending on the evolution equation, the wave packet's envelope may remain constant (no dispersion, see figure) or it may change (dispersion) while propagating.
Quantum mechanics ascribes a special significance to the wave packet: it is interpreted as a "probability wave", describing the probability that a particle or particles in a particular state will be measured to have a given position and momentum. It is in this way related to the wave function. Through application of the Schrödinger equation in quantum mechanics, it is possible to deduce the time evolution of a system, similar to the process of the Hamiltonian formalism in classical mechanics. The wave packet is thus a mathematical solution to the Schrödinger equation. The area under the absolute square of the wave packet solution is interpreted as the probability density of finding the particle in a region. The dispersive character of solutions of the Schrödinger equation has played an important role in rejecting Schrödinger's original interpretation, and accepting the Born rule.
In the coordinate representation of the wave (such as the Cartesian coordinate system), the localized position of the physical object's probability is given by the position of the packet. Moreover, the narrower the spatial wave packet, and therefore the better defined the position of the wave packet, the larger the spread in the momentum of the wave. This trade-off between spread in position and spread in momentum is one example of the Heisenberg uncertainty principle.
In the early 1900s, it became apparent that classical mechanics had some major failings. Isaac Newton originally proposed the idea that light came in discrete packets, which he called "corpuscles", but the wave-like behavior of many light phenomena quickly led scientists to favor a wave description of electromagnetism. It wasn't until the 1930s that the particle nature of light really began to be widely accepted in physics. The development of quantum mechanics — and its success at explaining confusing experimental results — was at the root of this acceptance. Thus, one of the basic concepts in the formulation of quantum mechanics is that of light coming in discrete bundles called photons. The energy of light photon is a discrete function of its frequency,
The ideas of quantum mechanics continued to be developed throughout the 20th century. The picture that was developed was of a particulate world, with all phenomena and matter made of and interacting with discrete particles; however, these particles were described by a probability wave. The interactions, locations, and all of physics would be reduced to the calculations of these probability amplitude waves. The particle-like nature of the world was significantly confirmed by experiment, while, at the very same time, the wave-like phenomena could be characterized as consequences of the wave-packet nature of particles.
Basic behaviors of wave packets
As an example of propagation without dispersion, consider wave solutions to the following wave equation:
where c is the speed of the wave's propagation in a given medium. Using the physics time convention, exp(−iωt), the wave equation has plane-wave solutions
- , and
This relation between ω and k should be valid so that the plane wave is a solution to the wave equation. It is called a dispersion relation.
To simplify, consider only waves propagating in one dimension (extension to three dimensions is straightforward). Then the general solution is
in which we may take ω= kc . The first term represents a wave propagating in the positive x-direction since it is a function of x−ct only; the second term, being a function of x+ct, represents a wave propagating in the negative x-direction.
A wave packet is a localized disturbance that results from the sum of many different wave forms. If the packet is strongly localized, more frequencies are needed to allow the constructive superposition in the region of localization and destructive superposition outside the region. From the basic solutions in one dimension, a general form of a wave packet can be expressed as
As in the plane-wave case the wave packet travels to the right for ω(k)=kc (since then u(x,t)=F(x−ct)) and to the left for ω(k)=−kc (since then u(x,t) = F(x+ct)).
The factor comes from Fourier transform conventions. The amplitude A(k) contains the coefficients of the linear superposition of the plane-wave solutions. These coefficients can in turn be expressed as a function of u(x,t) evaluated at t=0 by inverting the Fourier transform relation above:
For instance, choosing
The nondispersive propagation of the real or imaginary part of this wave packet is presented in the above animation.
As an example of propagation with dispersion, consider solutions to the Schrödinger equation (with m and ħ set equal to one)
yielding as dispersion relation
Once again, restricting ourselves to one dimension the solution to the Schrödinger equation satisfying the initial condition is found according to
An impression of the dispersive behaviour of this wave packet is obtained by looking at
(note that |u(x,t)| itself is not a solution of the Schrödinger equation).
Gaussian wavepackets in quantum mechanics
The above Gaussian wavepacket, unnormalized and just centered at the origin, instead, can now be written in 3D:
where a is a positive real number, the square of the width of the wavepacket, a = 2⟨r·r⟩/3⟨1⟩ = 2 (Δx)2.
The Fourier transform is also a Gaussian in terms of the wavenumber, the k-vector, (with inverse width, 1/a = 2⟨k·k⟩/3⟨1⟩ = 2 (Δpx/ħ)2, so that Δx Δpx=ħ/2, i.e. it saturates the uncertainty relation),
Each separate wave only phase-rotates in time, so that the time dependent Fourier-transformed solution is:
The inverse Fourier transform is still a Gaussian, but the parameter a has become complex, and there is an overall normalization factor.
The integral of Ψ over all space is invariant, because it is the inner product of Ψ with the state of zero energy, which is a wave with infinite wavelength, a constant function of space. For any energy eigenstate η(x), the inner product:
only changes in time in a simple way: its phase rotates with a frequency determined by the energy of η. When η has zero energy, like the infinite wavelength wave, it doesn't change at all. The integral ∫|Ψ|2d3r is also invariant, which is a statement of the conservation of probability. Explicitly,
The width of the Gaussian is the interesting quantity which can be read off from |Ψ|2:
The width eventually grows linearly in time, as ħt /m√a, indicating wave-packet spreading.
This linear growth is a reflection of the momentum uncertainty: the wavepacket is confined to a narrow width √a, and so has a momentum which is uncertain (according to the uncertainty principle) by the amount ħ/2√a, a spread in velocity of ħ/2m√a, and thus in the future position by ħt /m√a. (The uncertainty relation is then a strict inequality, far from saturation.)
The Airy wave train
In contrast to the above Gaussian wavepacket, it has been observed that a particular wavefunction based on Airy functions, propagates freely without dispersion, maintaining its shape. It accelerates undistorted in the absence of a force field: ψ=Ai(B(x−B ³ t ²)) exp(iB ³ t (x−2B ³ t ²/3)). (For simplicity, ħ=1, m=1/2, and B is a constant, cf. nondimensionalization.)
Nevertheless, Ehrenfest's theorem is still valid in this force-free situation, because the state is both non-normalizable and has an undefined (infinite) ⟨x⟩ for all times. (To the extent that it can be defined, ⟨p⟩ =0 for all times, despite the apparent acceleration of the front.)
In phase space, this is evident in the pure state Wigner quasiprobability distribution of this wavetrain, whose shape in x and p is invariant as time progresses, but whose features accelerate to the right, in accelerating parabolas B(x−B ³ t ²) + (p/B − tB ²)² = 0,
Note the momentum distribution obtained by integrating over all x is constant.
The narrow-width limit of the Gaussian wavepacket solution discussed is the free propagator kernel K. For other differential equations, this is sometimes called the Green's function, but in quantum mechanics it is traditional to reserve the name Green's function for the time Fourier transform of K.
Returning to one dimension for simplicity, when a is the infinitesimal quantity ε, the Gaussian initial condition, rescaled so that its integral is one:
becomes a delta function, δ(x), so that its time evolution,
gives the propagator.
Note that a very narrow initial wavepacket instantly becomes infinitely wide, but with a phase which is more rapidly oscillatory at large values of x. This might seem strange—the solution goes from being localized at one point to being "everywhere" at all later times, but it is a reflection of the enormous momentum uncertainty of a localized particle. Also note that the norm of the wavefunction is infinite, but this is also correct, since the square of a delta function is divergent in the same way.
The factor involving ε is an infinitesimal quantity which is there to make sure that integrals over K are well defined. In the limit that ε→0, K becomes purely oscillatory, and integrals of K are not absolutely convergent. In the remainder of this section, it will be set to zero, but in order for all the integrations over intermediate states to be well defined, the limit ε→0 is to be only taken after the final state is calculated.
The propagator is the amplitude for reaching point x at time t, when starting at the origin, x=0. By translation invariance, the amplitude for reaching a point x when starting at point y is the same function, only now translated,
In the limit when t is small, the propagator, of course, goes to a delta function,
but only in the sense of distributions: The integral of this quantity multiplied by an arbitrary differentiable test function gives the value of the test function at zero. To see this, note that the integral over all space of K equals 1 at all times:
since this integral is the inner-product of K with the uniform wavefunction. But the phase factor in the exponent has a nonzero spatial derivative everywhere except at the origin, and so when the time is small there are fast phase cancellations at all but one point. This is rigorously true when the limit ε→0 is taken at the very end.
So the propagation kernel is the (future) time evolution of a delta function, and it is continuous, in a sense: it goes to the initial delta function at small times. If the initial wavefunction is an infinitely narrow spike at position y,
it becomes the oscillatory wave:
Now, since every function can be written as a weighted sum of such narrow spikes,
the time evolution of every function ψ0 is determined by the propagation kernel K:
Thus, this is a formal way to express the general solution. The interpretation of this expression is that the amplitude for a particle to be found at point x at time t is the amplitude that it started at y, times the amplitude that it went from y to x, summed over all the possible starting points. In other words, it is a convolution of the kernel K with the arbitrary initial condition ψ0,
Since the amplitude to travel from x to y after a time t+t ' can be considered in two steps, the propagator obeys the identity:
which can be interpreted as follows: the amplitude to travel from x to z in time t+t ' is the sum of the amplitude to travel from x to y in time t, multiplied by the amplitude to travel from y to z in time t ', summed over all possible intermediate states y. This is a property of an arbitrary quantum system, and by subdividing the time into many segments, it allows the time evolution to be expressed as a path integral.
Analytic continuation to diffusion
The spreading of wavepackets in quantum mechanics is directly related to the spreading of probability densities in diffusion. For a particle which is randomly walking, the probability density function at any point satisfies the diffusion equation:
where the factor of 2, which can be removed by a rescaling either time or space, is only for convenience.
A solution of this equation is the spreading Gaussian,
and since the integral of ρt is constant, while the width is becoming narrow at small times, this function approaches a delta function at t=0:
again only in the sense of distributions, so that
for any smooth test function f.
The spreading Gaussian is the propagation kernel for the diffusion equation and it obeys the convolution identity:
which allows diffusion to be expressed as a path integral. The propagator is the exponential of an operator H:
which is the infinitesimal diffusion operator,
A matrix has two indices, which in continuous space makes it a function of x and x '. In this case, because of translation invariance, the matrix element K only depend on the difference of the position, and a convenient abuse of notation is to refer to the operator, the matrix elements, and the function of the difference by the same name:
Translation invariance means that continuous matrix multiplication:
is essentially convolution,
The exponential can be defined over a range of ts which include complex values, so long as integrals over the propagation kernel stay convergent,
As long as the real part of z is positive, for large values of x, K is exponentially decreasing, and integrals over K are indeed absolutely convergent.
The limit of this expression for z coming close to the pure imaginary axis is the Schrödinger propagator:
and this gives a more conceptual explanation for the time evolution of Gaussians.
From the fundamental identity of exponentiation, or path integration:
holds for all complex z values where the integrals are absolutely convergent so that the operators are well defined.
Thus, quantum evolution starting from a Gaussian, which is the diffusion kernel K,
gives the time evolved state,
This illustrates the above diffusive form of the Gaussian solutions,
- Joy Manners (2000). Quantum Physics: An Introduction. CRC Press. pp. 53–56. ISBN 978-0-7503-0720-8.
- Toda, Mikito (2005). Geometric structures of phase space in multidimensional chaos.... Hoboken, New Jersey: John Wiley & Sons inc. p. 123. ISBN 0-471-70527-6.
- Quantum Mechanics, E. Abers, Pearson Ed., Addison Wesley, Prentice Hall Inc, 2004, ISBN 978-0-13-146100-0
- Leonard I. Schiff (1968). Quantum mechanics (3rd ed.). London: McGraw-Hill.
- M. V. Berry and N. L. Balazs (1979), "Nonspreading wavepackets", Am J Phys 47: 264-267, doi:10.1119/1.11855.
- From the general pedagogy web-site: Curtright, T.L., Time-dependent Wigner Functions
- Jackson, J.D. (1975). Classical Electrodynamics (2nd Ed.). New York: John Wiley & Sons, Inc. ISBN 0-471-43132-X
- Feynman, R. P., and Hibbs, A. R., Quantum Mechanics and Path Integrals, New York: McGraw-Hill, 1965, ISBN 0-07-020650-3. (Dover, 2010, ISBN 0-486-47722-3.)