
Inphysics, awave packet (also known as awave train orwave group) is a short burst of localized wave action that travels as a unit, outlined by anenvelope. A wave packet can be analyzed into, or can be synthesized from, a potentially-infinite set of componentsinusoidal waves of differentwavenumbers, with phases and amplitudes such that they interfere constructively only over a small region of space, and destructively elsewhere.[1] Any signal of a limited width in time or space requires many frequency components around a center frequency within abandwidth inversely proportional to that width; even agaussian function is considered a wave packet because itsFourier transform is a "packet" of waves of frequencies clustered around a central frequency.[2] Each componentwave function, and hence the wave packet, are solutions of awave equation. Depending on the wave equation, the wave packet's profile may remain constant (nodispersion) or it may change (dispersion) while propagating.
Ideas related to wave packets –modulation,carrier waves,phase velocity, andgroup velocity – date from the mid-1800s. The idea of a group velocity distinct from a wave's phase velocity was first proposed byW.R. Hamilton in 1839, and the first full treatment was byRayleigh in his "Theory of Sound" in 1877.[3]
Erwin Schrödinger introduced the idea of wave packets just after publishing his famouswave equation.[4] He solved his wave equation for aquantum harmonic oscillator, introduced thesuperposition principle, and used it to show that a compact state could persist. While this work did result in the important concept ofcoherent states, the wave packet concept did not endure. The year after Schrödinger's paper,Werner Heisenberg published his paper on theuncertainty principle, showing in the process, that Schrödinger's results only applied toquantum harmonic oscillators, not for example toCoulomb potential characteristic of atoms.[4]: 829
The following year, 1927,Charles Galton Darwin exploredSchrödinger's equation for an unbound electron in free space, assuming an initialGaussian wave packet.[5] Darwin showed that at time later the position of the packet traveling at velocity would be
where is the uncertainty in the initial position.
Later in 1927Paul Ehrenfest showed that the time, for amatter wave packet of width and mass to spread by a factor of 2 was. Since is so small, wave packets on the scale of macroscopic objects, with large width and mass, double only atcosmic time scales.[6]: 49
Quantum mechanics describes the nature of atomic and subatomic systems usingSchrödinger's wave equation. The classical limit of quantum mechanics and many formulations of quantum scattering use wave packets formed from various solutions to this equation.Quantum wave packet profiles change while propagating; they show dispersion. Physicists have concluded that "wave packets would not do as representations of subatomic particles".[4]: 829
Schrodinger developed wave packets in hopes of interpreting quantum wave solutions as locally compact wave groups.[4] Such packets tradeoff position localization for spreading momentum. In the coordinate representation of the wave (such as theCartesian coordinate system), the position of the particle's localized probability is specified by the position of the packet solution. The narrower the spatial wave packet, and therefore the better localized the position of the wave packet, the larger the spread in themomentum of the wave. This trade-off between spread in position and spread in momentum is a characteristic feature of the Heisenberg uncertainty principle.
One kind of optimal tradeoff minimizes the product of position uncertainty and momentum uncertainty.[7]: 60 If we place such a packet at rest it stays at rest: the average value of the position and momentum match a classical particle. However it spreads out in all directions with a velocity given by the optimal momentum uncertainty. The spread is so fast that in the distance of once around an atom the wave packet is unrecognizable.
Particle interactions are calledscattering in physics; the wave packet concept plays an important role inquantum scattering approaches. A monochromatic (single momentum) source produces convergence difficulties in the scattering models.[8]: 150 Scattering problems also have classical limits. Whenever the scattering target (for example an atom) has a size much smaller than wave packet, the center of the wave packet follows scattering classical trajectories. In other cases, the wave packet distorts and scatters as it interacts with the target.[9]: 295
This sectionneeds additional citations forverification. Please helpimprove this article byadding citations to reliable sources in this section. Unsourced material may be challenged and removed.(October 2024) (Learn how and when to remove this message) |

Without dispersion the wave packet maintains its shape as it propagates.As an example of propagationwithout dispersion, consider wave solutions to the followingwave equation fromclassical physics
wherec is the speed of the wave's propagation in a given medium.
Using the physics time convention,e−iωt, the wave equation hasplane-wave solutions
where the relation between theangular frequencyω andangular wave vectork is given by thedispersion relation:such that. This relation should be valid so that the plane wave is a solution to the wave equation. As the relation islinear, the wave equation is said to benon-dispersive.
To simplify, consider the one-dimensional wave equation withω(k)= ±kc. Then the general solution iswhere the first and second term represent a wave propagating in the positive respectively negativex-direction.
A wave packet is a localized disturbance that results from the sum of many differentwave 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.[10] From the basic one-dimensional plane-wave solutions, a general form of a wave packet can be expressed aswhere the amplitudeA(k), containing the coefficients of thewave superposition, follows from taking theinverse Fourier transform of a "sufficiently nice" initial waveu(x,t) evaluated att = 0:and comes fromFourier transform conventions.
For example, choosing
we obtain
and finally
The nondispersive propagation of the real or imaginary part of this wave packet is presented in the above animation.


By contrast, in the case of dispersion, a wave changes shape during propagation. For example, thefree Schrödinger equation ,has plane-wave solutions of the form:where is a constant and the dispersion relation satisfies[11][12]with the subscripts denotingunit vector notation. As the dispersion relation is non-linear, the free Schrödinger equation isdispersive.
In this case, the wave packet is given by:where once again is simply the Fourier transform of. If (and therefore) is aGaussian function, the wave packet is called aGaussian wave packet.[13]
For example, the solution to the one-dimensional free Schrödinger equation (with2Δx,m, andħ set equal to one) satisfying the initial condition representing a wave packet localized in space at the origin as a Gaussian function, is seen to be
An impression of the dispersive behavior of this wave packet is obtained by looking at the probability density:It is evident that this dispersive wave packet, while moving with constant group velocityko, is delocalizing rapidly: it has awidth increasing with time as√ 1 + 4t2 → 2t, so eventually it diffuses to an unlimited region of space.


The above dispersive Gaussian wave packet, unnormalized and just centered at the origin, instead, att=0, can now be written in 3D, now in standard units:[14][15]The Fourier transform is also a Gaussian in terms of the wavenumber, thek-vector,Witha and its inverse adhering to theuncertainty relation such thatcan be considered thesquare of the width of the wave packet, whereas its inverse can be written as
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 now the parametera 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 anyenergy 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.
For a given, the phase of the wave function varies with position as. It variesquadratically with position, which means that it is different from multiplication by a linearphase factor as is the case of imparting a constant momentum to the wave packet. In general, the phase of a gaussian wave packet has both a linear term and a quadratic term. The coefficient of the quadratic term begins by increasing from towards as the gaussian wave packet becomes sharper, then at the moment of maximum sharpness, the phase of the wave function varies linearly with position. Then the coefficient of the quadratic term increases from towards, as the gaussian wave packet spreads out again.
The integral∫ |Ψ|2d3r is also invariant, which is a statement of the conservation of probability.[16] Explicitly,wherer is the distance from the origin, the speed of the particle is zero, and width given bywhich is√a at (arbitrarily chosen) timet = 0 while eventually growing linearly in time, asħt/(m√a), indicatingwave-packet spreading.[17]
For example, if an electron wave packet is initially localized in a region of atomic dimensions (i.e.,10−10 m) then the width of the packet doubles in about10−16 s. Clearly, particle wave packets spread out very rapidly indeed (in free space):[18] For instance, after1 ms, the width will have grown to about a kilometer.
This linear growth is a reflection of the (time-invariant) momentum uncertainty: the wave packet is confined to a narrowΔx =√a/2, and so has a momentum which is uncertain (according to the uncertainty principle) by the amountħ/√2a, a spread in velocity ofħ/m√2a, and thus in the future position byħt /m√2a. The uncertainty relation is then a strict inequality, very far from saturation, indeed! The initial uncertaintyΔxΔp =ħ/2 has now increased by a factor ofħt/ma (for larget).

A gaussian 2D quantum wave function:
where[19]
In contrast to the above Gaussian wave packet, which moves at constant group velocity, and always disperses, there exists a wave function based onAiry functions, that propagates freely without envelope dispersion, maintaining its shape, and accelerates in free space:[20]where, for simplicity (andnondimensionalization), choosingħ = 1,m = 1/2, andB an arbitrary constant results in
There is no dissonance withEhrenfest's theorem 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 could be defined,⟨p⟩ = 0 for all times, despite the apparent acceleration of the front.)
The Airy wave train is the only dispersionless wave in one dimensional free space.[21] In higher dimensions, other dispersionless waves are possible.[22]
Inphase space, this is evident in thepure stateWigner quasiprobability distribution of this wavetrain, whose shape inx andp is invariant as time progresses, but whose features accelerate to the right, in accelerating parabolas. The Wigner function satisfiesThe three equalities demonstrate three facts:
Note the momentum distribution obtained by integrating over allx is constant. Since this is theprobability density in momentum space, it is evident that the wave function itself is not normalizable.
The narrow-width limit of the Gaussian wave packet solution discussed is the freepropagator kernelK. For other differential equations, this is usually called theGreen's function,[23] but in quantum mechanics it is traditional to reserve the name Green's function for the time Fourier transform ofK.
Returning to one dimension for simplicity, withm and ħ set equal to one, whena is the infinitesimal quantityε, the Gaussian initial condition, rescaled so that its integral is one,becomes adelta function,δ(x), so that its time evolution,yields the propagator.
Note that a very narrow initial wave packet instantly becomes infinitely wide, but with a phase which is more rapidly oscillatory at large values ofx. This might seem strange—the solution goes from being localized at one point to being "everywhere" atall later times, but it is a reflection of the enormousmomentum uncertainty of a localized particle, as explained above.
Further note that the norm of the wave function is infinite, which is also correct, since the square of adelta function is divergent in the same way.
The factor involvingε is an infinitesimal quantity which is there to make sure that integrals overK are well defined. In the limit thatε → 0,K becomes purely oscillatory, and integrals ofK are not absolutely convergent. In the remainder of this section, itwill 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 pointx at timet, when starting at the origin,x=0. By translation invariance, the amplitude for reaching a pointx when starting at pointy is the same function, only now translated,
In the limit whent is small, the propagator goes to a delta functionbut only in the sense ofdistributions: The integral of this quantity multiplied by an arbitrary differentiabletest function gives the value of the test function at zero.
To see this, note that the integral over all space ofK equals 1 at all times,since this integral is the inner-product ofK with the uniform wave function. 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 wave function is an infinitely narrow spike at positiony,it becomes the oscillatory wave,
Now, since every function can be written as a weighted sum of such narrow spikes,the time evolution ofevery functionψ0 is determined by this propagation kernelK,
Thus, this is a formal way to express thefundamental solution orgeneral solution. The interpretation of this expression is that the amplitude for a particle to be found at pointx at timet is the amplitude that it started aty, times the amplitude that it went fromy tox,summed over all the possible starting points. In other words, it is aconvolution of the kernelK with the arbitrary initial conditionψ0,
Since the amplitude to travel fromx toy after a timet+t' can be considered in two steps, the propagator obeys the composition identity,which can be interpreted as follows: the amplitude to travel fromx toz in timet+t' is the sum of the amplitude to travel fromx toy in timet, multiplied by the amplitude to travel fromy toz in timet', summed overall 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 apath integral.[24]
The spreading of wave packets in quantum mechanics is directly related to the spreading of probability densities indiffusion. For a particle which israndomly walking, the probability density function satisfies thediffusion equation[25]where the factor of 2, which can be removed by rescaling either time or space, is only for convenience.
A solution of this equation is the time-varyingGaussian functionwhich is a form of theheat kernel. Since the integral ofρt is constant while the width is becoming narrow at small times, this function approaches a delta function att=0,again only in the sense of distributions, so thatfor anytest functionf.
The time-varying Gaussian is the propagation kernel for the diffusion equation and it obeys theconvolution identity,which allows diffusion to be expressed as a path integral. The propagator is the exponential of an operatorH,which is the infinitesimal diffusion operator,
A matrix has two indices, which in continuous space makes it a function ofx andx'. In this case, because of translation invariance, the matrix elementK 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 ofts which include complex values, so long as integrals over the propagation kernel stay convergent,As long as the real part ofz is positive, for large values ofx,K is exponentially decreasing, and integrals overK are indeed absolutely convergent.
The limit of this expression forz approaching the pure imaginary axis is the above Schrödinger propagator encountered,which illustrates the above time evolution of Gaussians.
From the fundamental identity of exponentiation, or path integration,holds for all complexz values, where the integrals are absolutely convergent so that the operators are well defined.
Thus, quantum evolution of a Gaussian, which is the complex diffusion kernelK,amounts to the time-evolved state,
This illustrates the above diffusive form of the complex Gaussian solutions,