Envelope (waves)

In physics and engineering, the envelope of an oscillating signal is a smooth curve outlining its extremes. The envelope thus generalizes the concept of a constant amplitude into an instantaneous amplitude. The figure illustrates a modulated sine wave varying between an upper and a lower envelope. The envelope function may be a function of time, space, angle, or indeed of any variable.

Example: beating waves

A common situation resulting in an envelope function in both space x and time t is the superposition of two waves of almost the same wavelength and frequency: : egin{align} F(x, t) & = sin left[ 2 pi left( frac {x}{lambda - Delta lambda } - ( f + Delta f )t ight) ight] + sin left[ 2 pi left( frac {x}{lambda + Delta lambda } - ( f - Delta f )t ight) ight] \[6pt] & approx 2cos left[ 2 pi left( frac {x} {lambda_{ m mod}} - Delta f t ight) ight] sin left[ 2 pi left( frac {x}{lambda} - f t ight) ight] end{align}

which uses the trigonometric formula for the addition of two sine waves, and the approximation Δλ ≪ λ:

:frac{1}{lambda pm Delta lambda}=frac {1}{lambda} frac{1}{1pmDelta lambda / lambda }approx frac{1}{lambda}mp frac {Delta lambda}{lambda^2} .

Here the modulation wavelength λmod is given by: : lambda_{ m mod} = frac {lambda^2}{Delta lambda} . The modulation wavelength is double that of the envelope itself because each half-wavelength of the modulating cosine wave governs both positive and negative values of the modulated sine wave. Likewise the beat frequency is that of the envelope, twice that of the modulating wave, or 2Δf.

If this wave is a sound wave, the ear hears the frequency associated with f and the amplitude of this sound varies with the beat frequency.

Phase and group velocity

The argument of the sinusoids above apart from a factor 2 are: :xi_C =left( frac {x}{lambda} - f t ight) , :xi_E=left( frac {x} {lambda_{ m mod}} - Delta f t ight) , with subscripts C and E referring to the carrier and the envelope. The same amplitude F of the wave results from the same values of ξC and ξE, each of which may itself return to the same value over different but properly related choices of x and t. This invariance means that one can trace these waveforms in space to find the speed of a position of fixed amplitude as it propagates in time; for the argument of the carrier wave to stay the same, the condition is: :left( frac {x}{lambda} - f t ight) = left( frac {x+Delta x}{lambda} - f (t + Delta t) ight) , which shows to keep a constant amplitude the distance Δx is related to the time interval Δt by the so-called phase velocity vp

:v_{ m p} = frac{Delta x}{Delta t} = lambda f . On the other hand, the same considerations show the envelope propagates at the so-called group velocity vg: :v_{ m g} = frac{Delta x}{Delta t} = lambda_{ m mod}Delta f =lambda^2 frac{Delta f}{Delta lambda} .

A more common expression for the group velocity is obtained by introducing the wavevector k: :k=frac{2pi}{lambda} . We notice that for small changes Δλ, the magnitude of the corresponding small change in wavevector, say Δk, is: : Delta k = left|frac{dk}{dlambda} ight|Delta lambda = 2pi frac{Delta lambda}{lambda^2} , so the group velocity can be rewritten as: : v_{ m g}= frac {2piDelta f}{Delta k} =frac {Delta omega}{Delta k} ,

where ω is the frequency in radians/s: ω = 2f. In all media, frequency and wavevector are related by a dispersion relation, ω = ω(k), and the group velocity can be written: :v_{ m g} =frac{domega (k)}{dk} .

In a medium such as classical vacuum the dispersion relation for electromagnetic waves is: :omega = c_0 k where c0 is the speed of light in classical vacuum. For this case, the phase and group velocities both are c0.

In so-called dispersive media the dispersion relation can be a complicated function of wavevector, and the phase and group velocities are not the same. For example, for several types of waves exhibited by atomic vibrations (phonons) in GaAs, the dispersion relations are shown in the figure for various directions of wavevector k. In the general case, the phase and group velocities may have different directions.

Example: envelope function approximation

In condensed matter physics an energy eigenfunction for a mobile charge carrier in a crystal can be expressed as a Bloch wave: :psi_{nmathbf{k}}(mathbf{r})=e^{imathbf{k}cdotmathbf{r}}u_{nmathbf{k}}(mathbf{r}) ,

where n is the index for the band (for example, conduction or valence band) r is a spatial location, and k is a wavevector. The exponential is a sinusoidally varying function corresponding to a slowly varying envelope modulating the rapidly varying part of the wavefunction un,k describing the behavior of the wavefunction close to the cores of the atoms of the lattice. The envelope is restricted to k-values within a range limited by the Brillouin zone of the crystal, and that limits how rapidly it can vary with location r.

In determining the behavior of the carriers using quantum mechanics, the envelope approximation usually is used in which the Schrödinger equation is simplified to refer only to the behavior of the envelope, and boundary conditions are applied to the envelope function directly, rather than to the complete wavefunction. For example, the wavefunction of a carrier trapped near an impurity is governed by an envelope function F that governs a superposition of Bloch functions: :psi( mathbf r )= sum_{mathbf k } F( mathbf k ) e^{imathbf{kcdot r}}u_{mathbf {k}}(mathbf r ) , where the Fourier components of the envelope F(k) are found from the approximate Schrödinger equation. In some applications, the periodic part uk is replaced by its value near the band edge, say k=k0, and then: :psi( mathbf r )approx left( sum_{mathbf k } F( mathbf k ) e^{imathbf{kcdot r}} ight)u_{mathbf{k}=mathbf{k}_0}(mathbf r ) = F( mathbf r )u_{mathbf{k}=mathbf{k}_0}(mathbf r ) .

Example: diffraction patterns

Diffraction patterns from multiple slits have envelopes determined by the single slit diffraction pattern. For a single slit the pattern is given by: :I_1=I_0 sin^2left(frac {pi d sin alpha}{lambda} ight) / left(frac {pi d sin alpha }{lambda} ight)^2 , where α is the diffraction angle, d is the slit width, and λ is the wavelength. For multiple slits, the pattern is :I_q = I_1 sin^2 left( frac {qpi g sin alpha} {lambda} ight) / sin^2 left( frac{ pi g sin alpha}{lambda} ight) , where q is the number of slits, and g is the grating constant. The first factor, the single-slit result I1, modulates the more rapidly varying second factor that depends upon the number of slits and their spacing.