the Creative Commons Attribution 4.0 License.

the Creative Commons Attribution 4.0 License.

# An excitation mechanism for discrete chorus elements in the magnetosphere

### Peter Bespalov

### Olga Savina

A beam pulsed amplifier mechanism responsible for effective amplification of short very low frequency (VLF) electromagnetic pulses is proposed. Effective amplification near the magnetic equator outside the plasmasphere is considered. A conditional growth rate of short electromagnetic pulses is calculated. Obtained results can explain some important features of the oblique electromagnetic chorus emissions without hiss-like radiation background.

Very low frequency (VLF) chorus emissions are very intense electromagnetic plasma waves that are naturally excited near the magnetic equatorial plane outside the plasmasphere (Burtis and Helliwell, 1969; Burton and Holzer, 1974; Tsurutani and Smith, 1974). Impressive experimental results of the chorus emission study were obtained within the framework of the CLUSTER project. These results have been presented in detail in many papers (e.g., Santolik, 2009). It is very important that a chorus is a succession of discrete emissions.

For an electromagnetic chorus with the wave vectors predominantly along the magnetic field, significant theoretical results were obtained. The usual kinetic cyclotron instability (Bespalov and Trakhtengerts, 1986) is sometimes possible for plasma parameters in the excitation region, but this instability typically does not have a sufficient growth rate to explain the rate of the chorus emission modification. The theory of chorus excitation based on the so-called backward wave oscillator (BWO) mechanism is now well known (see, e.g., Trakhtengerts, 1995; Trakhtengerts et al., 2007; Nunn et al., 2009). The implementation of this mechanism is closely related to the hiss-like radiation background and its dynamics. However, the chorus emissions are often recorded without the hiss-like radiation background. In Omura et al. (2008) and Fu et al. (2014) the generation process of the chorus emissions is analyzed by both theory and simulation assuming that the initial cyclotron wave growth is driven by the strong temperature anisotropy of energetic electrons. To explain the chorus spectrogram, the authors take into account inhomogeneity of the magnetic field and nonlinear wave decay, or the non-monotonic energy spectrum of particles. At present, there are significant achievements in the theoretical study of generation electrostatic chorus with the wave vectors near the resonance cone (see, e.g., Li et al., 2016a). There are extensive data on the wave normal angle measurements onboard THEMIS (Li et al., 2013) and a Van Allen probe (Li et al., 2016b).

Some problems connected with the theoretical analysis of chorus formation remain unsolved; for example, the excitation mechanism of an oblique electromagnetic chorus has not been studied. In this paper, we introduce a possible mechanism of oblique electromagnetic chorus excitation without hiss-like radiation background. This mechanism is related to the effective amplification of short electromagnetic pulses from the noise level even in a stable plasma. The amplification takes place in a suitable frequency band near the magnetic equator.

We consider small-scale wave processes near the local minimum of the magnetic
field, which typically is close to the magnetic equator where the plasma is
almost homogeneous. We use a linearized Vlasov equation for the disturbed
distribution function of energetic non-relativistic electrons *f*_{∼} in a
quasineutral cool background plasma

where *e* is the absolute value of the electron charge, *m* is the electron
mass, *c* is the speed of light in free space,
${v}_{z}=\mathrm{d}z/\mathrm{d}t$ and ${\mathbf{v}}_{\perp}=\mathrm{d}{\mathbf{r}}_{\perp}/\mathrm{d}t$ are the components of electron
velocity **v** along and across the ambient magnetic field
** B**,

*f*

_{∘}(

**v**) is the undisturbed distribution function, and

*E*_{∼}and

*B*_{∼}are the electric and magnetic field disturbances.

Let a short VLF electromagnetic pulse propagate in a homogeneous plasma along
the *z* axis directed along the magnetic field.
We assume that the short electromagnetic pulse has an envelope determined by
the function *A*(*ξ*) with a nearly step-like form, finite duration and a
unit value in the main body of the pulse

where $\mathit{\xi}=z-{v}_{\mathrm{g}z}t$, and *v*_{gz} is the
component of the group velocity along the magnetic field. The chosen pulse
shape is modeled and corresponds to the shot noise well known in electronics
(see, e.g., Rytov et al., 1989). It is convenient to replace *f*_{∼} in
Eq. (1) with an expression similar to Eq. (2):

Upon substitution of Eqs. (2) and (3) into Eq. (1) we obtain

We restrict ourselves to considering the energetic electrons for which

For these particles, the first term in Eq. (4) is absent because the disturbed distribution function is nonzero only in a small vicinity of the velocity which satisfies condition (5). So, under condition (5) the interaction of a short electromagnetic pulse with energetic electrons is described by Vlasov Eq. (4) without the first term. It is very important that the reduced Eq. (4) (without the first term) is identical to the Vlasov equation for wave-particle interaction in a homogeneous plasma. Note that the particles which do not satisfy Eq. (5) exist in the plasma, but their interaction with a short electromagnetic pulse is often not important, and we do not take this interaction into account.

Now we point out the condition favorable to the existence of a short electromagnetic pulse that propagates along the magnetic field at a constant velocity without the additional phase modulation and smallest dispersion distortion of the pulse front. The solution of this problem is known. According to electrodynamics (see, e.g., Sommerfeld, 1914; Jacson, 1962) this takes place under condition ${{v}_{\mathrm{ph}}}_{z}={v}_{\mathrm{g}z}$. This relationship agrees with Eq. (5) if the following two equalities are valid simultaneously:

where *v*_{ph}_{z} is the component of the phase velocity along the
magnetic field. The particles interact with a short electromagnetic pulse
under conditions (6) as in a homogeneous plasma.

If the ion motion in a relatively dense background plasma is not important,
then the dispersion equation of electromagnetic waves for the frequency band
${\mathit{\omega}}_{\mathrm{LHF}}<\mathit{\omega}<{\mathit{\omega}}_{\mathrm{B}}$ (where *ω*_{LHF} is the
lower-hybrid frequency and *ω*_{B} is the absolute value of the
electron cyclotron frequency) in the quasi-longitudinal approximation
(Helliwell, 1965; Ginzburg, 1970) takes the well-known form

Here, *k*_{z} and *k*_{⊥} are the whistler wave vector components along
and across the magnetic field, respectively, and *ω*_{p} is the
electron plasma frequency at the magnetic equator.

Note the first key point of our analysis. It is known (Helliwell, 1995) that
according to the dispersion Eq. (7) the conditions (6) are fulfilled independently of *k*_{⊥} for two selected velocities
along the magnetic field

where *u*_{G} is the Gendrin velocity. We will take into account only one
positive velocity in the intermediate calculations.

On the plane *z*, *t*, the specific features of particles and waves, which
satisfy Eq. (8), are explained in Fig. 1. Here *z* is the axis along the
magnetic field, the start and finish lines correspond to the wave-particle
interaction region boundaries, and *t* is the time counted from the pulse
crossing the start line. The resonance particles (Eq. 3) and the short electromagnetic pulse (Eq. 2)
are situated in the domain which moves along a narrow corridor between two
parallel dashed lines. The corridor width is determined by the pulse
duration. These particles and waves in the domain form a separated plasma
subsystem.

Effective wave-particle interaction in a homogeneous plasma in the magnetic
field takes place at the resonance conditions $\mathit{\omega}-{k}_{z}{v}_{z}=s{\mathit{\omega}}_{\mathrm{B}}$, where *s* is the integer. Such a resonance condition
is part of the mentioned equalities (Eq. 8) for the
Čerenkov resonance *s*=0. We limit ourselves to this possibility only:

We expect the following two simplifications to be fulfilled: the radiation power of an individual energetic electron corresponds to the so-called dipole approximation; the specific kinetic effects like the Landau damping are not important. Both of the mentioned simplifications take place for energetic electrons with a sufficiently small dispersion over the transverse and longitudinal velocities (Ginzburg, 1970):

where the overbar means averaging and *γ* is the instability growth
rate. The first inequality makes it possible to expand the Bessel functions
for small arguments (this so-called dipole approximation is used for
obtaining the expressions of the permittivity tensor). The second inequality
makes it possible not to take into account the kinetic effects (like the
Landau damping) for instability with a large hydrodynamic growth rate.

Note the second key point for a qualitative calculation of short electromagnetic pulse amplification. According to the previous comments, the electromagnetic field with suitable frequency in the short pulse changes its value as the electromagnetic wave field in a homogeneous plasma with an electron beam described by the effective distribution function:

where *n*_{b} is the density of the resonance electron beam. It can
be shown that for the wave-particle interaction at the Čerenkov
resonance, the results of calculations do not depend on the actual form of
the distribution function over perpendicular velocities if the first
inequality (10) takes place. We select the effective distribution
function (11) with a delta function over perpendicular velocities to
simplify the further algebra. On the other hand, the delta function over
velocities along the magnetic field exactly follows from the condition
(Eq. 8) for the selected velocity.

We explain additionally the expression for the effective distribution
function (Eq. 11). Assume that a short non-spreading
pulse propagates at a constant velocity (${v}_{\mathrm{ph}{}_{z}}={v}_{\mathrm{g}z}={u}_{\mathrm{G}}$) along the
magnetic field in a plasma with an arbitrary undisturbed distribution
function *f*_{∘}(**v**). The wave-particle interaction within a
short pulse is determined by the characteristic time of the electron velocity
variation *τ*. For correct accounting of the wave-particle interaction,
it is necessary to know the distribution function averaged over the timescale
*τ*. Let us consider the effective distribution function inside the
pulse. Inside a current pulse localization there is a beam of electrons which
have entered the interaction region together with the pulse and move jointly
with it. The contribution of these electrons to the effective distribution
function is proportional to *n*_{b}*δ*(*v*_{z}−*u*_{G}).
The suprathermal electrons with other velocities or flight moments do not
make an appreciable contribution to the effective distribution function
within the pulse since they cross the pulse in a narrow corridor in Fig. 1
too quickly on a timescale *τ* or do not have contact with the pulse.
Thus, we take the effective distribution function (Eq. 11) to analyze the evolution of a short electromagnetic pulse.

We expect the quasineutral plasma to consist of three fractions: unmoved
protons; cool electrons; and a weak electron beam along the magnetic field
without thermal dispersion, described by the effective distribution function
(Eq. 11). So, the complicated geophysical problem is
reduced to a typical problem of plasma physics (e.g., Akhiezer et al., 1975).
We now consider the permittivity tensor *ε*_{αβ} of
such a medium taking into account that the resonance beam density is
relatively small. The beam-related terms in the permittivity tensor produce
any effect only if they have a resonance pole at $\mathit{\omega}-{k}_{z}{u}_{\mathrm{G}}=\mathrm{0}$ (see Eq. 9). There is one such resonance term only
in the component *ε*_{33}. The permittivity tensor with the
mentioned resonance term has the following non-zero components:

where *n*_{p} is the background plasma density. The plasma
dispersion equation (Ginzburg, 1970) for a plane electromagnetic wave in
which all the variables are proportional to $\mathrm{exp}\mathit{\left\{}i\right(-\mathit{\omega}t+{k}_{\perp}x+{k}_{z}z\left)\mathit{\right\}}$ takes the well-known form

This determinant in the quasi-longitudinal approximation is reduced to the following dispersion equation, $(\mathit{\omega}-{\mathit{\omega}}_{w}({k}_{z},\mathit{\theta}\left)\right)(\mathit{\omega}-{k}_{z}{u}_{\mathrm{G}}{)}^{\mathrm{2}}=({n}_{\mathrm{b}}/\mathrm{4}{n}_{\mathrm{p}}){\mathit{\omega}}^{\mathrm{3}}{\mathrm{sin}}^{\mathrm{2}}\mathit{\theta}$, and finally, in accordance with Eqs. (7) and (8), we have

where the absolute value symbols take into account two selected velocities (Eq. 8) along the magnetic field. Recall that the dispersion Eq. (7) was found in the same quasi-longitudinal approximation. Exactly in accordance with Eq. (7) we have

where $\mathrm{cos}\mathit{\theta}={k}_{z}/({k}_{z}^{\mathrm{2}}+{k}_{\perp}^{\mathrm{2}}{)}^{\mathrm{1}/\mathrm{2}}$, *θ* being the
angle between the wave vector and the magnetic field.

The dispersion Eq. (14) is a cubic equation for the frequency *ω*.
This equation has three roots, one of which (see Fig. 2) for conditions close
to optimal (Eq. 8) corresponds to the short electromagnetic pulse instability
with the conditional growth rate

According to Eq. (16), the conditional growth rate of the short
electromagnetic pulse instability is maximum for cos^{2}*θ*_{BPA}=0.6. So, according to Eqs. (8) and (16), it is easy to find the following
optimal wave characteristics favorable to the maximum conditional growth rate
of a short electromagnetic pulse

Note that the resonance beam density *n*_{b} is significantly less
than the density of energetic electrons, but according to the additional
calculations, it is enough for explanation of the experimentally
recorded values of *γ*∼10^{2} s^{−1} in the chorus excitation region. Now it is possible to
estimate an electromagnetic pulse duration *t*_{p}. The beam pulsed
amplifier mechanism is effective for electromagnetic pulses with the duration
${t}_{\mathrm{p}}\le \mathrm{10}/{\mathit{\gamma}}_{\mathrm{BPA}}\sim \mathrm{0.1}$ s at the linear
stage of the pulse evolution.

The proposed beam pulsed amplifier (BPA) mechanism produces effective
amplification of short electromagnetic pulses with frequency close to
*ω*_{BPA} in the VLF frequency band near the magnetic equator. The
key conditions (Eqs. 8 and 17) take place outside the plasmasphere, where *u*_{G}
is close to the typical velocity of energetic non-relativistic electrons. In
this region the ratio $({n}_{\mathrm{b}}/{n}_{\mathrm{p}}{)}^{\mathrm{1}/\mathrm{3}}$ is not so small. This leads to
the amplification of short electromagnetic pulses from the noise level and
the formation of discrete VLF emissions.

The threshold of the BPA mechanism is mainly determined by the kinetic
Čerenkov damping (*s*=0) of electromagnetic waves in the epithermal
background plasma. The magnitude of the decay rate is considerably less than
the growth rate *γ*_{BPA} in the chorus excitation region.

By the further analysis it is possible to explain the gap between the lower
and upper chorus frequency bands which are well known experimentally (e.g.,
Bell et al., 2009). Actually, the short electromagnetic pulse excitation
takes place for frequency bands both below and above *ω*_{BPA}
because of the spectral pulse distortion and the resonance beam velocity
decrease as the pulse moves through the wave-particle interaction region.

Note that electromagnetic signals with a smooth envelope are not amplified due to BPA mechanisms during their propagation through the region near the magnetic equator. Therefore, it is possible to explain the excitation of chorus emissions without a hiss-like background.

BPA mechanism is not connected directly with the energetic electron
precipitation into the ionosphere because this mechanism is responsible for
the distribution function modification in the velocity space only far enough
from the loss cone. On the other hand, exactly after the amplification, a
strong electromagnetic pulse during its propagation to the ionosphere can
interact with more energetic electrons due to the cyclotron resonance
(*s*=1). Such an interaction produces precipitation of energetic electrons
into the ionosphere that is conjugate to the electromagnetic pulse.

The paper is theoretical and no new experimental data are used. All figure data are obtained from numerical calculation in MATLAB codes. Corresponding parameters are listed in the text.

PB proposed and analyzed the BPA mechanism, and wrote the paper. OS analyzed the BPA mechanism, and wrote the paper.

The authors declare that they have no conflict of interest.

Peter Bespalov was supported by RSF, project no. 16-12-10528 (Sect. 2), the Russian Ministry of Science and
Education, project no. 14.Z50.31.0007 (Sect. 3), and Fundamental Research Program no. 28 (numerical calculations).

The topical editor, Matina Gkioulidou, thanks one anonymous referee for help in evaluating this paper.

Akhiezer, A. I., Akhiezer, I. A., Polovin, R. V., Sitenko, A. G., and Stepanov, K. N.: Plasma electrodynamics. V.1. Linear theory, Pergamon Press, Oxford, UK, New York, USA, 432 pp., 1975.

Bell, T. F., Inan, U. S., Haque, N., and Pickett, J. S.: Source regions of banded chorus, Geophys. Res. Lett., 36, L11101, https://doi.org/10.1029/2009GL037629, 2009.

Bespalov, P. A. and Trakhtengerts, V. Yu.: Cyclotron instability of the Earth radiation belts, Reviews of Plasma Physics, Consultants Bureau, New York, USA, 10, 155–292, 1986.

Burtis, W. J. and Helliwell, R. A.: Banded chorus? A new type of VLF radiation observed in the magnetosphere by OGO 1 and OGO 3, J. Geophys. Res., 74, 3002–3010, https://doi.org/10.1029/JA074i011p03002, 1969.

Burton, R. K. and Holzer, R. E.: The origin and propagation of chorus in the outer magnetosphere, J. Geophys. Res., 79, 1014–1023, https://doi.org/10.1029/JA079i007p01014, 1974.

Fu, X., Cowee, M. M., Friedel, R. H., Funsten, H. O., Gary, S. P., Hospodarsky, G. B., Kletzing, C., Kurth, W., Larsen, B. A., Liu, K., MacDonald, E. A., Min, K., Reeves, G. D., Skoug, R. M., and Winske, D.: Whistler anisotropy instabilities as the source of banded chorus: Van Allen Probes observations and Particle-in-Cell simulations, J. Geophys. Res., 119, 8288–8298, https://doi.org/10.1002/2014JA020364, 2014.

Ginzburg, V. L.: The propagation of electromagnetic waves in plasmas, Pergamon Press, New York, USA, 2nd edn., 615 pp., 1970.

Helliwell, R. A.: Whistlers and related ionospheric phenomena, Stanford Univercity Press, Stanford, CA, USA, 349 pp., 1965.

Helliwell, R. A.: The role of the Gendrin mode of VLF propagation in the generation of magnetospheric emissions, Geophys. Res. Lett., 22, 2095–2098, https://doi.org/10.1029/95GL02003, 1995.

Jacson, J. D.: Classical electrodynamics, John Wiley & Sons, New York, USA, London, UK, 641 pp., 1962.

Li, W., Bortnik, J., Thorne, R. M., Cully, C. M., Chen, L., Angelopoulos, V., Nishimura, Y., Tao, J. B., Bonnell, J .W., and LeContel, O.: Characteristics of the Poynting flux and wave normal vectors of whistler-mode waves observed on THEMIS, J. Geophys. Res., 118, 1461–1471, https://doi.org/10.1002/jgra.50176, 2013.

Li, W., Mourenas, D., Artemyev, A. V., Bortnik, J., Thorne, R. M., Kletzing, C. A., Kurth, W. S., Hospodarsky, G. B., Reeves, G. D., Funsten, H. O., and Spence, H. E.: Unraveling the excitation mechanisms of highly oblique lower band chorus waves, Geophys. Res. Lett., 43, 8867–8875, https://doi.org/10.1002/2016GL070386, 2016a.

Li, W., Santolik, O., Bortnik, J., Thorne, R. M., Kletzing, C. A., Kurth, W. S., and Hospodarsky, G. B.: New chorus wave properties near the equator from Van Allen Probes wave observations, Geophys. Res. Lett., 43, 4725, https://doi.org/10.1002/2016GL068780, 2016b.

Nunn, D., Santolik, O., Rycroft, M., and Trakhtengerts, V.: On the numerical modelling of VLF chorus dynamical spectra, Ann. Geophys., 27, 2341–2359, https://doi.org/10.5194/angeo-27-2341-2009, 2009.

Omura, Y., Katoh, Y., and Summers, D.: Theory and simulation of the generation of whistler-mode chorus, J. Geophys. Res., 113, A04223, https://doi.org/10.1029/2007JA012622, 2008.

Rytov, S. M., Kravtsov, Yu. A., and Tatarskii, V. I.: Principles of Statistical Radiophysics 3. Elements of Random Fields, Springer-Verlag, Berlin Heidelberg, Germany, 249 pp., 1989.

Santolik, O., Gurnett, D. A., Pickett, J. S., Chum, J., and Cornilleau-Wehrlin, N.: Oblique propagation of whistler mode waves in the chorus source region, J. Geophys. Res., 114, A00F03, https://doi.org/10.1029/2009JA014586, 2009.

Sommerfeld, A.: Über die Fortpflanzung des Lichtes in dispergierenden Medien, Annalen der Physik, 349, 177–202, https://doi.org/10.1002/andp.19143491002, 1914.

Trakhtengerts, V. Y.: Magnetosphere cyclotron maser: Backward wave oscillator generation regime, J. Geophys. Res., 100, 17205–17210, https://doi.org/10.1029/95JA00843, 1995.

Trakhtengerts, V. Y., Demekhov, A. G., Titova, E. E., Kozelov, B. V., Santolik, O., Macusova, E., Garnett, D., Pickett, J. S., Rycroft, M. J., and Nunn, D.: Formation of VLF chorus frequency spectrum: Cluster data and comparision with the backward wave oscillator model, Geophys. Res. Lett., 34, L02104. https://doi.org/10.1029/2006GL027953, 2007.

Tsurutani, B. T. and Smith, E. J.: Postmidnight chorus: A substorm phenomenon, J. Geophys. Res., 79, 118–127, https://doi.org/10.1029/JA079i001p00118, 1974.

- Abstract
- Introduction
- Wave-particle interaction under special conditions
- The beam pulsed amplifier mechanism of the oblique electromagnetic chorus excitation due to the conditional instability of a short electromagnetic pulse
- Conclusions and discussion
- Data availability
- Author contributions
- Competing interests
- Acknowledgements
- References

- Abstract
- Introduction
- Wave-particle interaction under special conditions
- The beam pulsed amplifier mechanism of the oblique electromagnetic chorus excitation due to the conditional instability of a short electromagnetic pulse
- Conclusions and discussion
- Data availability
- Author contributions
- Competing interests
- Acknowledgements
- References