
1. Introduction
Long ago, Alfven first postulated (derived) the existence of Magneto-Hydrodynamic Waves―MHD waves―today called Alfven wave. For the history sake one can have a look at: reference [1].
Alfven (global or surface) wave is now very well accepted as a possible plasma Eigen mode for magnetized plasmas, it can be excited/present in fusion plasmas or space/astrophysical plasmas. As it is very well known the Alfven wave is a MHD Eigen mode where the ions provide the necessary inertia and the magnetic field lines support the effective tension that creates a restoring force for any equilibrium perturbation. Those modes can be excited by twisting the magnetic field lines (by an imposed dc electric current, for instance) or by compressing the plasma. The lines and the plasma are frozen together and the neutralizing electrons plasma is considered one fluid with ions [1] [2] [3] [4].
In past Alfven waves were not very well accepted as a possible plasma mode.
At that time, the scarce plasma physics community did not accept this finding promptly. However, after the observation of those plasma modes in solar wind, laboratory and magnetosphere/ionosphere, and behaving exactly as Alfven had postulated, the Alfven wave became then a respectable reality. After that many names were assigned to this mode, in different plasma physics communities; in space, in astrophysics, in technology and in fusion, inducing sometimes the bad idea that a new kind of “Eigen” modes were discovered. But, they are no other mode, they are with no caveat the well-known Alfven’s Alfven wave [1] [2] [3] [4] [5].
Alfven derived his new wave mode, using the very basic MHD equations. With a more elaborate calculation than Alfven’s original model, using for instance the plasma dielectric properties, any plasma mode can be obtained including the original Alfven wave.
In general, due to its intrinsic complicated dielectric properties plasmas can support infinite different Eigen modes. The scenarios get even more complicated in complex geometries such as toroidal and dipolar―as found in fusion machines, planetary magnetospheres and pulsar magnetospheres. Due to the intrinsic complications, the Alfven wave can be very diverse related to the properties of dispersivity, polarization and dissipation in different background plasmas, but it keeps being Alfven wave, with no caveat. Indeed, no new dramatic physics is found, just new geometric and thermal effects due to thermal dynamics and the curved background magnetic field; this curved magnetic field lines force the electrons and ions to follow curved trajectories and also permit particles to be trapped in magnetic traps.
In this paper we discuss the reason for so many names for the same MHD shear Alfven wave earlier derived by Alfven [6] - [31].
2. Alfven Wave Zoo
The importance to understand the Alfvenic plasma motions is of paramount importance to determine plasma stability, plasma heating, non-inductive current drive, and particle acceleration.
Therefore, before moving to any application Alfven wave-particle interaction processes, it is in order to discuss a bit more on what type of Alfven wave people are talking about and how they name it. This means that it is necessary to explore a little further the “animals” present in the Alfven wave’s zoo.
Moving now inside the Alfven wave’s zoo, we learn that they can be by nature of two types: shear and compressional, but in structure they can be global or surface. All of them can be derived from MHD equations with appropriate boundary conditions―it means bounded plasma or infinite plasma―as Alfven did. In magnetospheres, for instance, all types of these modes are observed, of course with different names used such as Alfven wave, MHD Alfven wave, kinetic Alfven wave, shear Alfven wave, inertial Alfven wave, compressional Alfven wave, Alfven surface wave and electromagnetic ion cyclotron wave (EMICW). Those names are used depending on the characteristic of the local plasma: MHD (cold) or kinetic (lukewarm/warm/hot), current carrying or not current carrying, homogeneous or nonhomogeneous. Geometrically speaking, for practical purposes, one can access the Alfven Eigen modes via Cartesian, cylindrical, toroidal or dipolar systems, since the geometrical model matters, it means that this can introduce modifications at their dispersion relation and one should not confuse those novel properties induced by real new physical effects. Complex geometry can aloud, for instance, different Alfven Eigen modes to couple M. Mahajan [11].
The global Alfven wave also called current carrying Alfven wave, TAE and GAE, should also be observed in some regions of the current carrying Space/Solar/Astrophysical plasmas.
If one is interested to study the particle acceleration processes in plasmas, they need to access the dissipation (energy and momentum) aspects of Alfven waves, therefore it is necessary to introduce kinetic corrections into the Alfven wave dispersion relation, in order to show how the non-dissipative MHD modes can be dissipated via the collisionless Cherenkov damping (Landau damping and TTMP―Transit Time Magnetic Pumping). Those dissipations arise due to the Alfven wave kinetic coupling to plasma electrostic modes giving raise to wave electric field aligned with the ambient magnetic field.
The Landau damping is related to the wave electric field and the Transit-Time Magnetic Pumping due to the wave magnetic field. At the linear regime, the plasma electromagnetic field is small, and therefore can be considered a first order expansion parameter to expand at all equations of the MHD system. If the electromagnetic field has finite amplitude the quadratic terms plays important role for plasma heating, current drive, and particle acceleration. Plasma heating occurs due to the wave energy thermalization into the bulk plasma and particle acceleration are more related to only the resonant electrons. If the thermalization time is smaller than the acceleration time, perceptible acceleration will be observed [15] - [29].
In order to move further, it is interesting to understand a bit more why people use using so many different names for the same Alfven mode? Is this just due to the tradition in plasma physics to create names and more names meaning just exactly the same thing, physically speaking? The point is, there are indeed only two Alfven modes, that is, the shear mode and the compressional mode in different geometries and derived with/without kinetic corrections. They have nothing really physically new from the original Alfven’s Alfven wave derived so many years ago. They are not a new Alfven wave so to say.
We can find in the literature all sorts of names for the Alfven shear mode: For instance, oblique Alfven wave, dispersive Alfven wave, Inertial Alfven wave, discrete Alfven wave, etc.. All of them refer to the same shear Alfven mode, of course in a certain plasma and frequency region particular conditions.
In order to move further with this conjecture, we consider a coordinate system in which z is the direction of the ambient magnetic field. In our model the x, y, and z coordinates correspond to the east west, north south, and vertical directions, respectively.
The general dispersion relation for the shear Alfven branch in homogenous uniformly magnetized plasmas is given by the elegant expression [12]:
(1)
where
,
and c are the total refraction index, the xx and zz
components of the dielectric tensor, the wave number, the speed of light and the wave frequency, respectively. This simple relation (1) allows two uncoupled independent modes, a shear mode and a compressible one.
This dispersion relation (Equation (1)) is as complicated or as complete as the dielectric tensor components, which are dependent on the plasma geometry and dimension models used. The model can be one fluid―MHD, two―fluids, n―fluids, hybrid (fluid ions + kinetic electrons or vice versa), kinetic, cold, lukewarm, hot, low beta, high beta, homogeneous, non-homogeneous, relativistic, non-relativistic, cartesian, cypindrical, toroidal or dipolar―1D, 2D or 3D.
From all those models, we can derive all the different modes that are part of the Alfven wave’s known and unknown zoo.
Note that the perpendicular wave vector component of the shear Alfven wave
is related to the kinetic correction terms
(LD),
(TTMP) and to the
ion Larmor radius correction that come from the component
. This correction is finite if we use a full kinetic plasma model. The information on the Landau wave-particle interaction (growth or dissipation) is described by the anti-Hermitian part of
and the finite ion Larmor correction is contained in the Hermitian part of
.
From this conjecture, it is clearly seen that only the electron mass correction is not enough to explain the necessary strong auroral electric potentials [thousands of volts], using the shear Inertial Alfven wave (IAW) as source mechanism. Indeed, the combination of finite electron mass and finite ion Larmor radius is the necessary ingredient to guarantee enough electric field intensity to explain the acceleration of electrons in aurora. If the ion Larmor radius is negligible the field aligned potential will be not zero of course, but too small.
The kinetic Alfven wave (kAw) has no transit-time magnetic pumping and cross term contributions described in
and
. The explanation is that this mode has finite
but very small
.
From Equation (1), we describe:
Now we can play around with Equation (1) to name the Alfven wave as found in the literature:
2.1. Alfven Waves Considering the Target Plasma Cold and as a Fluid [
It Means That the Wave Magnetic Field Aligned Electric Field Is Short Circuited in Those Plasmas
]
We get the regular well-known shear Alfven wave also called MHD Alfven wave [13].
For small wave number k and
not near to zero we get,
(1.1)
It can be seen how the ions dynamic introduces dispersivity into the shear Alfven wave dispersion relation. So, here, we just have basically the shear Alfven wave dispersion relation corrected for finite perpendicular wavelength and finite ion cyclotron frequency effects. If we assume in (1’) that
), we have.
or
(1.2)
For
, we get a more compact form
(1.3)
This mode is also called electromagnetic ion cyclotron wave (EMICW).
For parallel propagation (
) propagation, (1.3) reads,
.
and for
, we get from (1.3) the well known MHD shear Alfven wave dispersion relation.
.
This mode is non-dispersive and non-evanescent by nature, it is a natural “eigenmode” of magnetized homogeneous infinite plasmas. Of course it makes not much sense to talk about “Eigen” modes in an infinite plasma (not bounded), but we use this term here just to show that this mode is excitable in this type of plasma. The MHD shear Alfven wave has no electric field and no magnetic field aligned with the local background magnetic field and this mode is global by nature. It can be easily excited by internal or surface currents, Kelvin Helmoltz instability (sheared plasma flow) [28] and also creating a shear at the ambient magnetic field just by twisting the magnetic field lines somehow. For this mode,
.
The wave number
is undetermined (degenerated) for this mode. As mentioned by S. M. Mahajan [11], “this degeneracy with respect to
has profound consequences.” It leads to a continuous spectrum for shear Alfven waves in inhomogeneous plasma, even in a bounded system. To “regenerate”
we have to introduce an additional effect (for instance finite electron mass or finite ion Larmor radius), and in doing that we can get
explicitly at the shear Alfven wave dispersion relation. The group velocity for the MHD shear
Alfven wave mode is
, which means that all shear modes travel with
the same velocity, so a packet of shear Alfven waves do not get distorted during propagation, this feature creates a very interesting property for the shear Alfvenic turbulence in cold plasmas.
2.2. Shear Alfven with Finite Electron Inertia [
] and Cold Ions [
] Effects
In this scenario the shear Alfven wave is called by some authors inertial Alfven wave (IAW). I consider this a very misleading name to the shear Alfven wave, since this is nothing more than a special case of Equation (1), it is indeed a shear Alfven wave after all. The dispersion relation for the IAW is:
(1.4)
As said, the Inertial Alfven Wave (1.4) is just the shear Alfven wave derived from MHD equations, but to get the IAW a partial kinetic or two fluids corrections have to be considered, to bring out its propagation, polarization, and dissipation properties. What is the gain considering this correction? That is what one gets; a small magnetic field aligned wave electric field and so an explicitly determined
as we will see below.
Of course this wave could be also called kinetic Alfven wave, finite electron mass Alfven wave or inertial Alfven wave; however it is nothing more than the shear Alfven mode in the kinetic theory context. This mode has the non-vanishing small first order electric field component, aligned to the background magnetic field, since now
is finite implying that
.
This mode is evanescent for certain ranges of plasma densities [
], when
, and therefore it can be excited only by a tunneled surface wave electric
field, internal current sources, turbulence or by resonant mode conversion where
. The group velocity for this mode is:
.
This shows that the source of the shear Alfven wave dispersivity is to consider the electron inertia finite, which is neglected at the MHD one fluid equations―ions and electrons are not distinct.
2.3. Shear Alfven Wave with Quasi-Cold (Cold Ion and Warm Electrons) [
] [
] and Finite Electron Inertia [
] [
] [
]
In this case Equation (1) leads to:
(1.5)
where now the “Larmor radius” is emulated by
, since
.
This mode could be easily called Inertial Alfven wave with finite electron temperature correction and finite ion cyclotron frequency correction. If
, we have also
and so we recover the so called Inertial Alfven wave independent of the value of
.
The factor
[
] is due to the electron temperature correction and it plays a similar role as in the finite ion Larmor radius correction. The factor
is due to the electron inertia. The relative size of these two corrections can be evaluated by observing that
.
For
larger than
the finite Larmor radius correction dominates and vice versa. If
, we have:
and the propagation is allowed in the large vA region, whereas when
, the reverse is true and the dispersion relation (1.5)
becomes:
.
These cases can be already obtained in the MHD approximation using the modified Ohm’s Law, however in the MHD model it is not possible to account for properly the appearance of the mode dissipation (Landau damping).
2.4. Shear Alfven Wave with Relatively Cold Plasma [lukewarm Ions
and Warm Electrons
] [
]
In this case the electron inertia becomes important and the dispersion relation reads ―which is the Equation (1.5):
.
The hot low
[
] small Larmor radius plasma shear Alfven wave
Assuming that
, we find for (1.5):
(1.6)
This mode is already labeled as kinetic shear Alfven wave.
2.5. Shear Alfven Wave in Hot Low βplasmas [
] and Finite Ion Larmor Radius Plasma
In this case the shear Alfven wave has the dispersion relation described also by Equation (1), that is,
, (1.7)
with,
(1.8)
In (1.8),
is the modified Bessel’s function of zeroth order.
In (1.7) the Larmor radius is kept finite, no expansion is made. This means that this is the full kinetic shear Alfven wave. This dispersion relation reproduces, taking the appropriate limits, both; the Inertial Alfven Wave and the Kinetic Alfven Wave and, of course, all the intermediates nameless modes.
And so, the kinetic Alfven wave is just the same shear Alfven wave derived from MHD equations, but considering the kinetic approach. It is still the shear Alfven wave, but now both finite ion Larmor radius and finite electrons mass (inertia) corrections are considered to derive its propagation, polarization and dissipation properties.
The name kinetic Alfven wave “appears” when the dispersion relation of the shear Alfven wave contains the finite electron mass and the finite ion Larmor radius corrections, this name was first used by Akira Hasegawa [16] in the context of inhomogeneous plasmas. However, this dispersion relation can easily be obtained from the Alfven’s dielectric tensor/dispersion relation, for homogeneous plasmas, as we have shown. In high beta plasmas the kinetic Alfven wave is coupled to the ion acoustic wave.
Now, we consider the ion Larmor radius finite, no expansion is made. This means that this is the full shear kinetic Alfven wave dispersion relation. This dispersion relation contemplates, taking the appropriate limits, both; the so called Inertial Alfven Wave (IAW) and the kAw.
When I = 0 we recover the dispersion relation for the MHD shear Alfven wave, it means when
.
Since,
. (1.7’)
Depending on the signal of
the wave can be evanescent. This is directly
related to the background magnetic field and plasma density. When this happens
(
) where
. This must always be checked before we can say that we have observed a shear Alfven wave. This is more complex in non-homogeneous and layered plasma where surface wave is normally present and has propagation property not so different as the shear Alfven at the field aligned background magnetic field direction.
The group velocity is given by,
. (1.9)
Equation (1.9) shows how dispersive is the kAw and again the dispersivity is due to the electron inertia and due to the additional term related to the finite corrections of the ion Larmor radius.
Since aurora is highly populated with current systems, it is in order to ask how background currents affect the shear Alfven wave dispersion relation. The Current carrying shear Alfven wave (called in nuclear fusion by global Alfven wave, GAE, TAE et cetera) is nothing more than the well-known shear Alfven mode where a local dc current system is present (toroidal, dipolar et cetera)―for aurora dynamics the current might be the dc current system aligned with the ambient magnetic field. The current creates a poloidal dc magnetic field (emulated by the y-direction in our coordinate system), which modifies directly the known shear Alfven wave dispersion relation derived from the non-current carrying MHD plasma.
What we gain considering current carrying plasmas? First of all, the current system creates a dispersivity. The current carrying mode can be excited directly (without the need of surface wave mode conversation) by any shear or compression in the ambient magnetic field. It propagates with certain level of dissipation, due to a small but finite magnetic field aligned electron field and magnetic field, very similar situation to the compressional case. This mode gets some compressive property due to
, which is undetermined (degenerated) in the case of a shear Alfven mode derived from pure MHD equations. In this way, this mode can be dissipated by either Landau damping or TTMP or by both depending on the situation
2.6. The Alfven Surface Wave
There are two types of Alfven surface wave (Asw), the shear-type (non compressive) [22] and the compressional [21]. As I have already mentioned, the Asw are natural eigenmodes of sharp density boundaries magnetized plasmas, and evanescent eigenmodes of smooth density boundaries magnetized plasmas [21] [22] [27]. For the smooth-boundary condition, surface waves can very efficiently couple energy to the kAw via the resonant mode conversion. The Asw can be easily confused with the kinetic Alfven wave if the polarization properties of the observed turbulence are not looked at carefully, and this is due to its frequency being close to the bulk shear Alfven wave frequency, i.e.;
(2)
(1.9’)
The group velocity (1.9’) of the surface wave (2) is slightly larger than that of the shear ideal Alfven wave, however, again, the all modes at the wave packet travel with the same constant group velocity
,at this approximation.
In order to introduce the final clarification to this issue, now, if one looks at Stix’s book [14] one can see that the left-hand cut-off, the right-hand cut-off, the Alfven resonance, the surface wave, the fast magnetosonic wave and the kinetic Alfven wave are closely related. Stix just plots the perpendicular index of
refraction
versus plasma density.
He gets the following interesting results and final results to this Alfven wave ensemble name problem:
1) Starting from low density side, the square of the refraction index plotted in this case for the compressional Alfven waves according to cold-plasma (two-fluid) theory shows that for increasing density
evolves from an evanescent (
) mode to a propagating one at the left-hand cut-off point
. However, it diverges at the Alfven resonance
. Starting from higher density side
evolves from a propagating mode to an evanescent mode exactly at the Alfven resonance
, so both compressional modes (low frequency and high frequency) are evanescent at
.
2) Incorporating now, the finite-Larmor-radius corrections (but zero electron temperature), from the hot-plasma theory, to the dispersion relation. In doing so, Stix has gotten for this case (called “cool”) a short-wavelength propagating mode on the low density side of the resonance and he called it “surface wave”. This new situation does not show any singular behaviour at the Alfven resonance region. It means that the surface modes comes from low density side as evanescent up to close to the left-hand cut-off region and then propagate back to the lower density side. The fast mode propagates from the higher density side to the lower density side and just after the right-hand cut-off it becomes an evanescent mode.
3) Incorporating into the dispersion relation both; the finite-Larmor-radius and finite electron temperature corrections, using the hot-plasma theory, Stix has gotten for this case (called “warm”) a surprising short-wavelength propagating mode on the higher density side of the resonance. This mode propagation on the higher density side is labelled by other people kinetic Alfven wave. Since it is somehow related to the previous lower density propagating modes, it is evanescent for plasma densities where we have the match
.
Going further with this analysis, Stix mention that the finite electron temperature has a remarkable effect similar to ion temperature in changing the wave “cool” wave behavior to “warm” wave behavior. And going further he says that whereas raising the ion temperature increases the finite-Larmor-radius effects, raising the electron temperature so that changes sign of the zz dielectric tensor component
.
The Alfven shear mode has only
when derived from MHD equations. However, if kinetic corrections are included a perpendicular component
appears and the dispersion relation becomes
(3)
Note that the electron dynamic along the magnetic field lines gives rise to a finite field aligned electric field
, in both Alfven modes namely the shear mode and the compressional mode described by (3). For the compressional mode this electric field is related to the wave field aligned magnetic field
through
.
It is important to note that the Alfven surface wave and kinetic Alfven wave are coupled in non-homogeneous kinetic plasmas, in this kinetic plasma model the differential equation that describes the electric field in the direction of the non-homogeneity is of fourth order and the coupling is due to the finite ion Larmor radius kinetic correction. Surface waves in inhomogeneous plasmas are evanescent in the direction of the plasma non homogeneity. For sharp nonhomogeneity jump the surface wave is a propagating mode in the direction of the ambient magnetic field but it is evanescent in the direction of the plasma nonhomogeneity. So, the electric field in the direction of the plasma gradient is
localized at a distance inside the plasma of
.
Due to the nature of the surface wave and the kinetic Alfven the former can couple energy to the latter. Since, as we have seen, the surface wave is an excitable eigenmode at the lower density side of the plasma and the kinetic Alfven is evanescent at the lower density side of the plasma and cannot be excited directly at this region, the connection between surface waves and kinetic Alfven can guarantee the energy coupling from surface to kinetic Alfven wave, what is called in kinetic theory resonant mode conversion.
Here,
, (4)
where,
,
,
,
.
It is interesting what Hollweg [30] says about the kinetic Alfven wave: “the wave becomes strongly compressive when
is of the order of the ion inertial length. Thus, in low-β plasma, the kinetic Alfven wave can be compressive at values of
for which the dispersion relation departs only slightly from that of the usual MHD Alfven wave. The compression is accompanied by a magnetic field fluctuation
such that the total pressure perturbation
. Thus the wave undergoes transit-time damping as well as Landau damping; the two effects are comparable if the ion thermal speed is of the order of the Alfven speed. We find that the transverse electric field is elliptically polarized but rotating in the electron sense; this surprising behaviour of the polarization of the Alfven branch was discovered numerically by Gary [31] ”.
It is also interesting that Holweg [30] explicitly shows how the kinetic Alfven wave takes on some properties of the large-
limit of the slow mode, as can also be seen from (4).
3. Conclusion
The Alfven wave is, with no caveat, the most interesting eigenmode for magnetized plasmas. They can be externally excited to heat and drive non-inductive current in fusion plasmas and also in space and astrophysical plasmas. They can accelerate ions and electrons, induce magnetic reconnections, and other remarkable effects in space. These modes can also couple to each other and with other different plasma wave modes, and therefore they can excite other frequency range modes such as ion acoustic waves and whistler waves. Alfven waves can exhibit so interesting behaviours that even now after so many years of its discovery it is still fascinating and it is challenging plasma physicist, from the different plasma communities; as from astrophysics, stellar, space, fusion, and technological areas, to still study and enjoy them [11] [14] [15] [17] - [24].
Acknowledgements
The authors (A. S. de Assis and C. E. da Silva) would like to thank Yahweh for granting us the beauties of knowledge and a learning spirit.