Perpendicular Wave Propagation

Let us now consider wave propagation, through a warm plasma, perpendicular to the equilibrium magnetic field. For perpendicular propagation, $k_\parallel\rightarrow 0$, and, hence, from Equation (7.85), $\xi_n\rightarrow \infty$. Making use of the asymptotic expansions (7.86) and (7.87), the matrix $T_{\alpha\beta}$ simplifies considerably. The dispersion relation can again be written in the form (7.93), where

$\displaystyle M_{11}$ $\displaystyle = 1 - \sum_s\frac{{\mit\Pi}_s^{2}}{\omega}
\frac{{\rm e}^{-\lambd...
...!\sum_{n=-\infty,\infty} \frac{n^2\,I_n(
\lambda_s)}{\omega-n\,{\mit\Omega}_s},$ (7.102)
$\displaystyle M_{12}$ $\displaystyle = -M_{21} =- {\rm i}\sum_s \frac{{\mit\Pi}_s^{2}}{\omega}
\,\,{\r...
...rac{n\,\left[I_n'(\lambda_s)
-I_n(\lambda_s)\right]}{\omega-n\,{\mit\Omega}_s},$ (7.103)
$\displaystyle M_{22}$ $\displaystyle = 1 -\frac{k_\perp^{2}\,c^2}{\omega^2} \- \sum_s \frac{{\mit\Pi}_...
...
I_n(\lambda_s) - 2\,\lambda_s^{2}\,I_n'(\lambda_s)}{\omega-n\,{\mit\Omega}_s},$ (7.104)
$\displaystyle M_{33}$ $\displaystyle = 1 - \frac{k_\perp^{2}\,c^2}{\omega^2} - \!
\sum_s \frac{{\mit\P...
...a_s} \!\sum_{n=-\infty,\infty}
\frac{I_n(\lambda_s)}{\omega-n\,{\mit\Omega}_s},$ (7.105)

and $M_{13} = M_{31} = M_{23}=M_{32}=0$. Here,

$\displaystyle \lambda_s = \frac{(k_\perp\, \rho_s)^2}{2},$ (7.106)

where $\rho_s=v_s/\vert{\mit\Omega}_s\vert$ is the species-$s$ gyroradius.

The first root of the dispersion relation (7.93) is

$\displaystyle n_\perp^{2}=\frac{k_\perp^{2}\,c^2}{\omega^2} = 1 - \!
\sum_s \fr...
...a_s} \!\sum_{n=-\infty,\infty}
\frac{I_n(\lambda_s)}{\omega-n\,{\mit\Omega}_s},$ (7.107)

with the eigenvector $(0,\,0,\,E_z)$. This dispersion relation obviously corresponds to the electromagnetic plasma wave, or ordinary mode, discussed in Section 5.10. However, in a warm plasma, the dispersion relation for the ordinary mode is strongly modified by the introduction of resonances (where the refractive index, $n_\perp$, becomes infinite) at all the harmonics of the cyclotron frequencies:

$\displaystyle \omega_{n\,s} = n\,{\mit\Omega}_s,$ (7.108)

where $n$ is a non-zero integer. These resonances are a finite gyroradius effect. In fact, they originate from the variation of the wave phase across a gyro-orbit (Cairns 1985). Thus, in the cold plasma limit, $\lambda_s\rightarrow 0$, in which the gyroradii shrink to zero, all of the resonances disappear from the dispersion relation. In the limit in which the wavelength, $\lambda$, of the wave is much larger than a typical gyroradius, $\rho_s$, the relative amplitude of the $n$th harmonic cyclotron resonance, as it appears in the dispersion relation (7.107), is approximately $(\rho_s/\lambda)^{\vert n\vert}$ [see Equations (7.88) and (7.106)]. It is clear, therefore, that, in this limit, only low-order resonances [i.e., $n\sim {\cal O}(1)$] couple strongly into the dispersion relation, and high-order resonances (i.e., $\vert n\vert \gg 1$) can effectively be neglected. As $\lambda\rightarrow
\rho_s$, the high-order resonances become increasingly important, until, when $\lambda\lesssim \rho_s$, all of the resonances are of approximately equal strength. Because the ion gyroradius is generally much larger than the electron gyroradius, it follows that the ion cyclotron harmonic resonances are generally more important than the electron cyclotron harmonic resonances.

Observe that the cyclotron harmonic resonances appearing in the dispersion relation (7.107) are of zero width in frequency space: that is, they are just like the resonances that appear in the cold-plasma limit. Actually, this is just an artifact of the fact that the waves we are studying propagate exactly perpendicular to the equilibrium magnetic field. It is clear, from an examination of Equations (7.83) and (7.85), that the cyclotron harmonic resonances originate from the zeros of the plasma dispersion functions. Adopting the usual rule that substantial damping takes place whenever the arguments of the dispersion functions are less than or of order unity, it follows that the cyclotron harmonic resonances lead to significant damping whenever

$\displaystyle \omega - \omega_{n\,s} \lesssim k_\parallel\,v_s.$ (7.109)

Thus, the cyclotron harmonic resonances possess a finite width in frequency space provided the parallel wavenumber, $k_\parallel$, is non-zero: that is, provided the wave does not propagate exactly perpendicular to the magnetic field.

The appearance of the cyclotron harmonic resonances in a warm plasma is of great practical importance in plasma physics, because it greatly increases the number of resonant frequencies at which waves can transfer energy to the plasma. In magnetic fusion experiments, these resonances are routinely exploited to heat plasmas via externally launched electromagnetic waves (Stix 1992; Swanson 2003).

Figure: 7.8 Dispersion relation for electron Bernstein waves in a warm plasma for which $\omega _{UH}/\vert{\mit \Omega }_e\vert=2.5$.
\includegraphics[height=3.5in]{Chapter07/fig7_8.eps}

The other roots of the dispersion relation (7.93) satisfy

  $\displaystyle \left(1 - \sum_s\frac{{\mit\Pi}_s^{2}}{\omega}
\frac{{\rm e}^{-\l...
...ga-n\,{\mit\Omega}_s}
\right) \left(1 -\frac{k_\perp^{2}\,c^2}{\omega^2}\right.$    
  $\displaystyle \phantom{=}\left.-
\sum_s \frac{{\mit\Pi}_s^{2}}{\omega} \frac{{\...
...mbda_s) - 2\,\lambda_s^{2}\,I_n'(\lambda_s)}
{\omega-n\,{\mit\Omega}_s}
\right)$    
  $\displaystyle = \left(\sum_s \frac{{\mit\Pi}_s^{2}}{\omega}
\,\,{\rm e}^{-\lamb...
...eft[I_n'(\lambda_s)
-I_n(\lambda_s)\right]}{\omega-n\,{\mit\Omega}_s}\right)^2,$ (7.110)

with the eigenvector $(E_x,\, E_y,\, 0)$. In the cold plasma limit, $\lambda_s\rightarrow 0$, this dispersion relation reduces to that of the extraordinary mode discussed in Section 5.10. This mode, for which $\lambda_s\ll 1$, unless the plasma possesses a thermal velocity approaching the velocity of light, is little affected by thermal effects, except close to the cyclotron harmonic resonances, $\omega=\omega_{n\,s}$, where small thermal corrections are important because of the smallness of the denominators in the previous dispersion relation (Cairns 1985).

However, another mode also exists. In fact, if we look for a mode with a phase-velocity much less than the velocity of light (i.e., $c^2\,k_\perp^2/\omega^2\gg 1$) then it is clear from (7.102)–(7.105) that the dispersion relation is approximately

$\displaystyle 1 - \sum_s\frac{{\mit\Pi}_s^{2}}{\omega}
\frac{{\rm e}^{-\lambda_...
...m_{n=-\infty,\infty} \frac{n^2\,I_n(
\lambda_s)}{\omega-n\,{\mit\Omega}_s} = 0,$ (7.111)

and the associated eigenvector is $(E_x, \,0, \,0)$. The new waves, which are called Bernstein waves—after I.B. Bernstein, who first discovered them (Bernstein 1958)—are a type of slowly propagating, longitudinal, electrostatic wave.

Let us consider electron Bernstein waves, for the sake of definiteness. Neglecting the contribution of the ions, which is reasonable provided that the wave frequencies are sufficiently high, the dispersion relation (7.111) reduces to

$\displaystyle 1 - \frac{{\mit\Pi}_e^{2}}{\omega}
\frac{{\rm e}^{-\lambda_e}}{\l...
...m_{n=-\infty,\infty} \frac{n^2\,I_n(
\lambda_e)}{\omega-n\,{\mit\Omega}_e} = 0.$ (7.112)

In the limit $\lambda_e\rightarrow 0$ (with $\omega\neq n\,{\mit\Omega}_e$), only the $n=\pm 1$ terms survive in the previous expression. [See Equation (7.88).] In fact, because $I_{\pm 1}(\lambda_e)/\lambda_e \rightarrow
1/2$ as $\lambda_e\rightarrow 0$, the dispersion relation yields

$\displaystyle \omega^2\rightarrow {\mit\Pi}_e^{2} + {\mit\Omega}_e^{2}.$ (7.113)

It follows that there is a Bernstein wave whose frequency asymptotes to the upper hybrid frequency (see Section 5.10) in the limit $k_\perp\rightarrow 0$. For other non-zero values of $n$, we have $I_n(\lambda_e)/\lambda_e
\rightarrow 0$ as $\lambda_e\rightarrow 0$. However, a solution to Equation (7.111) can be obtained if $\omega\rightarrow n\,{\mit\Omega}_e$ at the same time. Similarly, as $\lambda_e\rightarrow\infty$, we have ${\rm e}^{-\lambda_e}\,I_n(\lambda_e)\rightarrow
0$ (Abramowitz and Stegun 1965). In this case, a solution can only be obtained if $\omega\rightarrow n\,{\mit\Omega}_e$, for some $n$, at the same time. The complete solution to Equation (7.111) is plotted in Figure 7.8 for a case where the upper hybrid frequency lies between $2\,\vert{\mit\Omega}_e\vert$ and $3\,\vert{\mit\Omega}_e\vert$. In fact, wherever the upper hybrid frequency lies, the Bernstein modes above and below it behave like those shown in the diagram.

At small values of $k_\perp$, the phase-velocity becomes large, and it is no longer legitimate to neglect the extraordinary mode (Cairns 1985). A more detailed examination of the complete dispersion relation shows that the extraordinary mode and the Bernstein mode cross over near the harmonics of the cyclotron frequency to give the pattern shown in Figure 7.9. Here, the dashed line shows the cold plasma extraordinary mode.

Figure: 7.9 Dispersion relation for extraordinary/electron Bernstein waves in a warm plasma for which $\omega _{UH}/\vert{\mit \Omega }_e\vert=2.5$ and $v_{t\,e}/c=0.2$. The dashed line indicates the cold plasma extraordinary mode.
\includegraphics[height=3.5in]{Chapter07/fig7_9.eps}

In a lower frequency range, a similar phenomena occurs at the harmonics of the ion cyclotron frequency, producing ion Bernstein waves, with somewhat similar properties to electron Bernstein waves. Note, however, that while the ion contribution to the dispersion relation can be neglected for high-frequency waves, the electron contribution cannot be neglected for low-frequency waves, so there is not a complete symmetry between the two types of Bernstein waves.