next up previous
Next: The Brillouin precursor Up: Electromagnetic wave propagation in Previous: The method of stationary

The group velocity

The point of stationary phase, defined by $\partial\phi/\partial\omega
=0$, satisfies the condition
\begin{displaymath}
\frac{c}{v_g} = \frac{ct}{x},
\end{displaymath} (768)

where
\begin{displaymath}
v_g = \frac{d \omega}{d k}
\end{displaymath} (769)

is conventionally termed the group velocity. Thus, the signal seen at position $x$ and time $t$ is dominated by the frequency range whose group velocity $v_g$ is equal to $x/t$. In this respect, the signal incident at the surface of the medium ($x=0$) at time $t=0$ can be said to propagate through the medium at the group velocity $v_g(\omega)$.

The simple one-resonance dispersion relation (4.86) yields

\begin{displaymath}
\frac{c}{v_g} \simeq n(\omega) \left[1
+\frac{\omega^2}{\ome...
...+\frac{\omega^2}
{\omega^2-\omega_0^{~2}-\omega_p^{~2}}\right]
\end{displaymath} (770)

in the limit $g\rightarrow 0$, where
\begin{displaymath}
n(\omega) =\frac{c k}{\omega} =
\sqrt{\frac{\omega_0^{~2} + \omega_p^{~2}-\omega^2}{\omega_0^{~2}-\omega^2}}.
\end{displaymath} (771)

The variation of $c/v_g$ and the refractive index $n$ with frequency is sketched in Fig. 12. With $g=0$ the group velocity is less than $c$ for all $\omega$, except for $\omega_0<\omega<
\omega_1\equiv\sqrt{\omega_0^{~2} +\omega_p^{~2}}$, where it is purely imaginary. Note that the refractive index is also complex in this frequency range. The phase velocity $v_p=c/n$ is subluminal for $\omega<\omega_0$, imaginary for $\omega_0
\leq \omega\leq \omega_1$, and superluminal for $\omega>\omega_1$.

Figure 12: The typical variation of the functions $c/v_g(\omega )$ and $n(\omega )$. Here, $\omega_1=(\omega_0^{~2} + \omega_p^{~2})^{1/2}$.
\begin{figure}\epsfysize =2.5in
\centerline{\epsffile{group.eps}}\end{figure}

The frequency range which contributes to the amplitude at time $t$ is determined graphically by finding the intersection of a horizontal line with ordinate $ct/x$ with the solid curve in Fig. 12. There is no crossing of the two curves for $t<t_0 \equiv x/c$, thus no signal arrives before this time. For times immediately following $t_0$ the point of stationary phase is seen to be at $\omega\rightarrow\infty$. In this large $\omega$ limit the point of stationary phase is given by

\begin{displaymath}
\omega_s \simeq \omega_p \sqrt{\frac{t_0}{2(t-t_0)}}.
\end{displaymath} (772)

Note that $\omega=-\omega_s$ is also a point of stationary phase. It is easily demonstrated that
\begin{displaymath}
\phi_s \simeq -2\sqrt{\xi(t-t_0)},
\end{displaymath} (773)

and
\begin{displaymath}
\phi_s'' \simeq -2\,\frac{(t-t_0)^{3/2}}{\xi^{1/2}},
\end{displaymath} (774)

with
\begin{displaymath}
F(\omega_s) \simeq \frac{t-t_0}{\tau\,\xi}.
\end{displaymath} (775)

Here, $\xi$ is given by Eq. (4.109). The stationary phase approximation (4.133) gives
\begin{displaymath}
f_s \simeq \sqrt{\frac{\pi \,\xi^{1/2} }
{(t-t_0)^{3/2}}} \,...
... e}^{-2{\rm i} \sqrt{\xi(t-t_0)}
+3\pi {\rm i}/4} +{\rm c.c.},
\end{displaymath} (776)

where c.c. denotes the complex conjugate of the preceding term (this contribution comes from the second point of stationary phase located at $\omega=-\omega_s$). The above expression reduces to
\begin{displaymath}
f_s \simeq \frac{ 2\sqrt{\pi}}{\tau} \,\frac{(t-t_0)^{1/4}}{\xi^{3/4}}
\cos\!\left[2\sqrt{\xi(t-t_0)} - 3\pi/4\right].
\end{displaymath} (777)

It is easily demonstrated that the above formula is the same as the expression (4.118) for the Sommerfeld precursor in the large argument limit $t-t_0 \gg 1/\xi$. Thus, the method of stationary phase yields an expression for the Sommerfeld precursor which is accurate at all times except those immediately following the first arrival of the signal.


next up previous
Next: The Brillouin precursor Up: Electromagnetic wave propagation in Previous: The method of stationary
Richard Fitzpatrick 2002-05-18