- 1Department of Astronomy, University of Wisconsin, Madison, WI, United States
- 2Centro de Investigación en Astronomía, Universidad Bernardo O’Higgins, Santiago, Chile
- 3Institute for Advanced Study, Hubble Fellow, Princeton, NJ, United States
Alfvénic component of MHD turbulence damps Alfvénic waves. The consequences of this effect are important for many processes, from cosmic ray (CR) propagation to launching outflows and winds in galaxies and other magnetized systems. We discuss the differences in the damping of the streaming instability by turbulence and the damping of a plane parallel wave. The former takes place in the system of reference aligned with the local direction of magnetic field along which CRs stream. The latter is in the reference frame of the mean magnetic field and traditionally considered in plasma studies. We also compare the turbulent damping of streaming instability with ion-neutral collisional damping, which becomes the dominant damping effect at a sufficiently low ionization fraction. Numerical testing and astrophysical implications are also discussed.
1 Propagation of Alfvén Waves in MHD Turbulence
Astrophysical media are turbulent and magnetized (see a collection of relevant reviews in [1]). The propagation of Alfvén waves in turbulent magnetized media is an important astrophysical problem that influences fundamental astrophysical processes (see e.g., [2–5]). This review focuses on the damping of Alfvén waves in MHD turbulence. The Alfvén waves can arise from instabilities induced by cosmic rays (CRs), e.g., from the streaming of CRs [6–9], and the gyroresonance instability related to the compression of magnetic field and CRs (see [10]). They can also be generated by large scale perturbations of magnetic field (see [11] and ref. therein [12]).
Turbulent damping suppresses its growth and affects the streaming speed of CRs. As a result, turbulent damping of streaming instability is important for studies on the diffusion and acceleration of CRs in shocks, galaxies, and galaxy clusters ([13–15], [3,16–18]), stellar wind launching (e.g., [19,20]), and galaxy evolution (e.g., [21]).
It should be noted that the well-known study of Alfvén wave damping by turbulent plasmas performed by Silimon and Sudan [22] employed an unrealistic model of isotropic MHD turbulence. Later, turbulent damping of Alfvén waves was mentioned as a process for suppressing CR streaming instability in Yan and Lazarian [23], henceforth YL02. This process was quantified by Farmer and Goldreich [24], henceforth FG04, where the Goldreich and Sridhar [25]; henceforth, GS95 model of Alfvénic turbulence with scale-dependent anisotropy was adopted. The limitation of the aforementioned study was that for the calculations it was assumed that turbulence is injected isotropically with the turbulent velocity uL exactly equal to the Alfvén velocity VA, i.e., Alfvén Mach number MA equal to unity. In addition, only turbulent damping of streaming instability was considered.
Following the study in Lazarian [26], we will seperately discuss the turbulent damping of Alfvén waves that are generated by streaming instability and by large-scale magnetic perturbations. We will demonstrate the strong dependence of turbulent damping on MA in various turbulence regimes and astrophysical media with different levels of medium magnetization. In §2 we provide the derivation of the Alfvénic turbulent scaling. In §3 we describe the turbulent damping of Alfvén waves generated by streaming instability in the reference system aligned with the local direction of turbulent magnetic field. In §4 we discuss the turbuelnt damping of Alfvén waves induced by large-scale magnetic perturbations in a global system of reference. We compare the turbulent damping with ion-neutral collisional damping of streaming instability in a partially ionized medium in §5. The numerical testing of the theoretical predictions is provided in §6. The discussion of the astrophysical implications on propagation of CRs in galaxies and launching of winds follows in §7. The summary is given in §8.
2 Derivation of Alfvénic Turbulent Scaling
In Alfvénic turbulence the relative perturbations of velocities and magnetic fields are related as follows:
where Bl is the fluctuation of the magnetic field B at scale l, BL is the fluctuation of the magnetic field at the driving scale L of turbulence. Correspondingly, ul is the turbulent velocity fluctuation at the scale l and uL is the turbulent velocity at L. MA = uL/VA is the Alfvén Mach number.
One way to understand the non-linear interactions of Alfvén waves within the MHD turbulent cascade is to consider colliding Alfvén wave packets with parallel scales l‖ and perpendicular scales l⊥. The collision of a wave packet induces an energy change.
where the term in brackets manifests the change of the energy of a wave packet induced by its interaction with the oppositely moving Alfvén wave packet. The time of this interaction is equal to the time of the passage of these wave packets through each other. As the size of the packet is l‖, the interaction time is simply Δt ∼ l‖/VA.
The rate of turbulent energy cascade is related to the rate of structure change of the oppositely moving wave packet. The latter is ul/l⊥. As a result, Eq. 2 provides
The fractional change of packet energy taking place per collision is ΔE/E. This characterises the strength of the nonlinear turbulent interaction:
In Eq. 4, f is the ratio of the shearing rate of the wave packet, i.e., ul/l⊥, to its propagation rate, i.e., VA/l‖.
One can identify two distinct cases. If f ≪ 1, the shearing rate is significantly smaller than the propagation rate, and the cascade presents a random walk process. Therefore
steps are required for the energy cascade, and therefore the cascading time is
ℵ > 1 corresponds to the weak turbulent cascade. Naturally, ℵ cannot become less than unity. Therefore, the limiting case is ℵ ≈ 1. This is the case of strong MHD turbulence.
Traditionally, the wavevectors are defined in the system of reference related to the mean field. However, the system of reference related to a wave packet with given parallel and perpendicular dimensions is more relevant when dealing with strong MHD turbulence. We take this into account by considering Alfvén wave packets having the dispersion relation ω = VA|k‖|, where we use
In weak turbulence, the decrease of l⊥ while l‖ does not change signifies the increase of the energy change per collision. This forces ℵ to be of the order of unity. In this case one gets
in strong turbulence, which signifies the cascading time being equal to the wave period
The cascade of turbulent energy satisfies the relation [27]:
which for the hydrodynamic cascade provides
where the relation for the cascading time tcas ≈ l/ul is employed.
For the weak turbulent cascade with ℵ≫ 1, we have (LV99)
where Eqs 6, 8 are used. The isotropic turbulence injection at scale L results in the second relation in Eq. 10. Taking into account that for the weak turbulence l‖ is constant, it is easy to see that Eq. 10 provides
which is different from the hydrodynamic
It was shown in LV99 that for turbulence with isotropic injection at scale L with VL < VA the transition to the strong regime corresponding to ℵ ≈ 1 happens at the scale
As a result, the inertial range of weak turbulence is limited, i.e.,
which follows from ℵ ≈ 1 condition given by Eqs 4, 5.
The scaling relations for the strong turbulence with VL < VA can be easily obtained. The turbulence is strong and cascades over one wave period, which according to Eq. 7 is equal to l⊥/ul. Substituting the latter in Eq. 8 one gets
The latter energy cascading rate is analogous to that in an ordinary hydrodynamic Kolmogorov cascade. However, this cascading takes place in the direction perpendicular to the local direction of the magnetic field2.
This strong MHD turbulence cascade starts at ltrans and its injection velocity is given by Eq. 13. This provides another way to obtain the Alfvénic turbulent scaling in strong turbulence regime (LV99).
which can be rewritten in terms of the injection velocity uL (see Eq. 15).
Substituting this in Eq. 7 we get the relation between the parallel and perpendicular scales of the eddies (LV99):
The relations Eqs 15, 17 reduce to the GS95 scaling for transAlfvénic turbulence if MA ≡ 1.
In the opposite case we deal with superAlfvénic turbulence, i.e., with uL > VA. As a result, at scales close to the injection scale the turbulence is essentially hydrodynamic as the influence of magnetic forces is marginal. Therefore, the velocity is Kolmogorov.
The magnetic field becomes more important at smaller scales and the cascade changes its nature at the scale.
at which the turbulent velocity becomes equal to the Alfvén velocity [28]. The rate of cascade for l < lA is:
Unlike the case of subAlfvénic turbulence, the case of superAlfvénic turbulence can be reduced to the case of transAlfvénic turbulence, but with lA acting as the injection scale. At scales l < lA.
The relations for subAlfvénic and superAlfvénic tubulence that we obtain above coincide with the expressions first obtained in Lazarian and Vishniac [29] using a different approach. These expressions will be used below in our discussion on turbulent damping of Alfvén waves.
3 Turbulent Damping of Streaming Instability
Linear Alfvén waves undergo non-linear cascading when they propagate through Alfvénic turbulence. This process is of MHD nature and the non-linear damping of Alfvén waves does not depend on plasma microphysics. The interaction between CR-driven Aflven waves and turbulence is similar to that of oppositely moving wave packets of turbulent cascade.
The Alfvén waves emitted parallel to the local magnetic field experiences the least distortions from the oppositely moving eddies. Thus the least distorted Alfvén waves are those with the largest value of l⊥. Indeed, the larger l⊥, the longer time it takes for the evolution of the oppositely moving wave packets. For instance, for strong GS95 turbulence the time corresponds to
The case of Alfvén waves parallel to the local direction of magnetic field corresponds to streaming and gyroresonance instabilities. In what follows, we will focus on the streaming instability. The dispersion of magnetic field directions with respect to the mean magnetic field determines the corresponding l⊥. Naturally, the turbulent damping of Alfvén waves is different for weak turbulence and strong turbulence. Thus we will separately discuss turbulent damping of streaming instability in different turbulence regimes.
3.1 Streaming Instability and Local System of Reference
The streaming instability of CRs happens as CR particles moving in one direction scatter back from a magnetic field perturbation and thus increase the amplitude of the perturbation. The induced perturbations are Alfvén waves. If the Alfvén waves are severely damped, the CR particles can stream freely along the magnetic field.
Physically, the generation of Alfvén waves takes place as CRs stream along the local magnetic field. During the process the sampling scale for the magnetic field is the CR Larmor radius rL. In this setting one should consider the process in the system of reference related to the local direction3 of the wondering magnetic field (LV99, [31–33]).
In the direction parallel to the local magnetic field, the growth rate of the streaming instability is (see [7]):
where ΩB = eB/mc is the nonrelativistic gyrofrequency, ncr is the number density of CRs with gyroradius rL > λ = γmc2/eB, and γ is the Lorentz factor. If the growth rate given by Eq. 23 is less than the rate of turbulent damping, the streaming instability is suppressed.
3.2 Damping by SubAlfvénic Strong Turbulence
Our first approach is based on calculating the distortion of Alfvén waves by MHD turbulence as the waves propagate along magnetic field. The cause of the wave distortion is the field line wandering over angle θx. This angle is determined by the amplitude of magnetic field fluctuations δBx that are induced by turbulent eddies with perpendicular scale x. One can see that the distortion induced during the time t is
where the fluctuation induced by turbulence evolves as
In the above expression ux denotes the velocity corresponding to the magnetic field fluctuation δBx. The time t in Eq. 25 is chosen to be less than the eddy turnover time x/ux. As a result, the ratio reflects the partial sampling of the magnetic perturbation by the wave. By using the velocity scaling of strong subAlfvénic turbulence for ux in Eq. 25, it is easy to rewrite Eq. 24 as
The wave damping corresponds to the “resonance condition” δx = λ, where λ is the wavelength. Inserting this in Eq. 26 we obtain the perpendicular scale of the “resonance” magnetic fluctuations that distort the Alfvén waves:
The time required to damp the Alfvén waves is equal to the turnover time of the “resonant” eddy:
This provides the rate of non-linear damping of the Alfvén waves,
or
where the subsscript “s” denotes “strong turbulence.” For transAlfvénic turbulence, i.e., MA = 1, this result was obtained in FG04. The square of the Alfvén Mach number dependence presented in Eq. 30 means a significant change of the damping rate compared to the transAlfvénic case4.
If the injection of turbulence is isotropic, the maximal perpendicular scale of strong subAlfvénic motions is
The streaming CRs generate Alfvén waves at a scale comparable to the gyroradius rL. Thus it requires that
which is a notable limitation on CR energy if MA is small. The CRs with larger energies interact with weak turbulence as we will discuss in Section 3.3.
Due to the importance of turbulent damping of streaming instability, it is advantageous to provide another derivation of Eq. 30. This alternative derivation is based on the picture of propagating wave packets that we used while obtaining Eq. 3. Consider two oppositely moving Alfvén wave packets with the perpendicular scale
Using the scaling in Eqs 15, 25, we derive the “resonant” perpendicular scale x:
This can be used to determine the rate of damping defined as ΓsubA,s ≈ ux/x. This coincides with the earlier result given by Eq. 30. Then the maximal wavelength of the non-linearly damped Alfvén waves can be obtained from Eq. 33 if the scale ltrans is used instead of x, i.e.
Naturally, the latter coincides with the result given by Eq. 31. The minimal scale of non-linearly damped waves depends on the perpendicular scale of the smallest Alfvénic eddies lmin. The full range of rL for which turbulent damping is essential can be obtained by using Eq. 33 and the scaling of strong turbulence given by Eq. 15:
The value of lmin depends on the particular damping process of MHD turbulent cascade, which can be relatively large in a weakly ionized gas (see [34]). Due to the differences of rL for protons and electrons, Eq. 36 presents a possible situation when the streaming instability of CR electrons is not damped by turbulence, while it is damped for CR protons.
We note that the turbulent damping of streaming instability for
The latter clearly illustrates the inefficiency of damping when turbulence has the perpendicular scale larger than the “resonant” scale.
3.3 Damping by SubAlfvénic Weak Turbulence
In many instances the weak turbulence is not important. It has a limited inertial range and transfers to strong turbulence at smaller scales. However, as we show below, this may not be true for wave damping by turbulence. For wavelengths longer than λmax ,s the wave is non-linearly damped through interactions with the wave packets of the weak turbulence, having perpendicular scales given by Eq. 33. Naturally, the scaling of weak turbulence given by Eq. 11 should be used. This provides the relation between the Alfvén wave wavelength and the perpendicular scale of the “resonant” weak turbulence perturbation.
This delivers the perpendicular scale
According to Eqs 5, 6, the weak turbulence packets cascade ℵ times slower compared to the case of strong turbulence. Taking into account that the parallel scale of weak turbulence wave packets is equal to the injection scale L, we have
The rate of turbulent damping of the wave is therefore
which gives
where the subsscript “w” denotes “weak turbulence.” Note that compared to the case of damping given by Eq. 30 we now have a stronger dependence on MA, as well as a different scaling with the wavelength λ.
The maximal wavelength of the Alfvén waves that is cascaded by the weak cascade we derive by substituting l⊥ = L, i.e., using the energy injection scale in Eq. 38. This gives:
Thus for CRs that generate Alfvén waves, their rL should satisfy
to interact with weak Alfvénic turbulence. The underlying assumption here is that
Waves with λ > λmax ,w can interact with the turbulent motions at the injection scale L. The cascade of such waves is induced by the largest wave packets at a rate
which does not depend on wavelength. Physically, this means that all waves in the range LMA < λ < L decay at the same rate that is determined by the restructuring of the magnetic field at the injection scale.
The above expression is valid for λ < L. In the case of λ ≫ L the rate is reduced due to the random walk, which results in a factor (L/λ)2, i.e.,
The latter result is relevant for the damping induced by turbulence injected at scales smaller than the wavelength.
In terms of the dependence of damping rate on λ for subAlfvénic turbulence, we observe that the dependence becomes stronger with the increase of λ up to λ = LMA. For λ less than
In terms of the dependence of turbulent damping rate on MA, it changes from
3.4 Damping by SuperAlfvénic Turbulence
As we discussed earlier, if turbulence is superAlfvénic, at large scales the effects of magnetic field are marginal and turbulence is hydrodynamic-like. However, the turbulent velocity decreases with the decreasing scale and at a scale lA becomes equal to the Alfvén velocity. This scale can be considered as the injection scale of transAlfvénic turbulence. Therefore, the case of Alfvén wave damping by superAlfvénic turbulence at scales less than lA can be related to the case of damping by transAlfvénic turbulence considered in FG04. Indeed, a simple substitution of L by lA provides the required rate of magnetic structure evolution on scales less than lA. This gives:
Treating lA as the effective injection scale and using Eq. 31, it is easy to obtain the maximal wavelength up to which our treatment of the non-linear damping is applicable:
Associating λ with rL, we define the corresponding range of rL
assuming that the minimal/damping scale of turbulent motions lmin is less than
For Alfvén waves with λ larger than that given by Eq. 48 and therefore for
where we use Eq. 18.
Similar to the case of sub-Alfvénic turbulence, in superAlfvénic case, we observe the change of the rate of Alfvén wave damping changing from λ−1/2 for short wavelengths to λ−2/3 for λ longer than
3.5 Other Forms of Presenting Our Results
The scaling of weak turbulence is different from that of strong turbulence that starts at the transition scale
The cascading rate of the weak turbulence is given by Eq. 10 and we can write it as
This reflects the decrease of energy dissipation by
The peculiar feature of Eq. 52 is that if one formally substitutes instead of ϵw the cascading rate of strong turbulence, one will get the expression in FG04. This is exactly the universality of expressions that we sought. Nevertheless, this analogy is only formal as the cascading rate for weak turbulence is
For wavelengths in the range
The expression given by Eq. 53 demonstrates a slower damping rate in comparison to Eq. 52. The decrease of damping rate by the factor
For damping of Alfvén waves generated by CRs with larger rL, i.e. LMA < rL < L (see Eq. 38) we obtain:
For superAlfvénic turbulence at scales less than lA, by expressing MA from Eq. 20 and substituting it in Eq. 47, we obtain
Formally, the above expression coincides with the expression for the damping by subAlfvénic strong turbulence given by Eq. 52. Nevertheless, the subAlfvénic turbulence demonstrates the significant reduction of the cascading rate compared to the transAlfvénic turbulence. On the contrary, the superAlfvénic strong MHD turbulence corresponds to a significant increase of dissipation rate in comparison with the transAlfvénic case. Thus, for the same injection scale L and the same injection velocity VL, the damping of Alfvén waves depends on the magnetization of media. At a lower magnetization, e.g., for superAlfvénic turbulence, the damping of Alfvén waves is more efficient than that at a higher medium magnetization, i.e., for the subAlfvénic case.
As we discussed earlier, in superAlfvénic turbulence, the long Alfvén waves with λ larger than
where the hydrodynamic dissipation rate is
Below we present a few more forms of presenting the damping rates that we obtained above. For instance, it could be sometimes useful to rewrite the expressions given by Eqs 30, 42 in terms of λmax ,s given by Eq. 35. We remind the reader that the physical meaning of λmax ,s is the longest wavelength that still interacts with strong turbulent cascade. Then,
and
It is easy to see that Eq. 57 demonstrates that the damping by strong MHD turbulence ΓsubA,s happens faster than the Alfvén crossing rate of the injection scale eddies. In the case of weak turbulence, Eq. 58 demonstrates that ΓsubA,w is slower than the above rate.
4 Turbulent Damping of Alfvén Waves Generated in the Global System of Reference
The turbulent damping of Alfvén waves generated by streaming CRs is an important special case of turbulent damping as the streaming instability induces Alfvén waves that are aligned with the local direction of magnetic field. Another case arises if we consider the damping of a flux of Alfvén waves generated by an extended source. The difference between the two cases is that in the latter setting the waves are generated irrespectively to the local direction of magnetic field. Therefore, such Alfvén waves should be viewed in the global system of reference related to the mean magnetic field. As a result, our earlier treatment of the Alfvén wave damping by MHD turbulence should be modified.
4.1 Case of Strong SubAlfvénic Turbulence
Consider an Alfvén wave generated at an angle θ ≫ δB/B with respect to the global mean magnetic field. In this situation it is natural to disregard the dispersion of angles that arises from magnetic wandering induced by turbulence6. To distinguish these two cases we use sin θ instead of sin θx in Eq. 33. In this case the perpendicular scale of eddies that the waves interact with is given by:
For strong turbulence the rate of the wave damping is equal to the turnover rate of subAlfvénic eddies. Therefore using Eq. 59, we find
This provides the non-linear damping rate of an Alfvén wave moving at the angle θ with respect to the mean field.
Using the expression of weak turbulent cascading rate ϵw (see Eq. 10), one can write:
The turbulent damping given by Eq. 61 is applicable to
where lmin is the perpendicular damping scale, and
Naturally, the adopted approximation θ ≫ δB/B fails if the wave is launched parallel to the mean magnetic field. The directions of the local magnetic field deviates from the mean field and this makes the actual θ0 different from zero. In the global system of reference the dispersion is dominated by the magnetic field variations presented at the injection scale (see [33]). Therefore
Substituting this into Eq. 60 we get
The above expression is different from Eqs 30, 52. The difference stems from the different properties of Alfvén waves generated in the local system versus the global system of reference. The damping rate in Eq. 64 is applicable to the range of wavelength.
the latter result trivially follows from Eqs 62, 63.
4.2 The Case of Weak SubAlfvénic Turbulence
For weak subAlfvénic turbulence, in the case θ ≫ δB/B, one should use Eq. 59 to relate λ to the scale of perpendicular motions that the wave strongly non-linearly interacts with. To obtain the damping rate, Eq. 41 should be used:
where we use the weak cascading rate ϵw. The range of wavelength for this type of damping is
The last inequality is obtained by substituting the maximal perpendicular eddy scale LMA for x in Eq. 59.
In the case of Alfvén wave propagation along the mean magnetic field, one should use Eq. 63 to get
Using Eqs 63, 67, we find the range of wavelength that is subject to the turbulent damping:
4.3 Other Cases
After illustrating the difference of non-linear damping for waves generated in the local reference system of magnetic field and in the global reference system of the mean field, we can provide results for other cases. More detailed discussion was presented in Lazarian [26]. For instance, for superAlfvénic turbulence, there is
where in superAlfvénic turbulence angle θ varies from one turbulent eddy of size lA to another. As a result, the corresponding averaging over such changing directions should be performed. For the random distribution of the relevant directions, the corresponding geometric factor is ⟨ sin2/3θ⟩ = 3/5.
On scales larger than lA, MHD turbulence is marginally affected by magnetic fields. As a result, no difference between local and global frames is present in terms of Alfvén wave damping. This difference also disappears for the damping by turbulent fluctuations at the injection scale.
4.4 Summary of Main Results in Sections 3 and 4 on Turbulent Damping
Some of our results for non-linear turbulent damping of Alfvén waves in different turbulence regimes are summarized in Table 1.
We show results relevant both to the damping of waves in the local system of reference, e.g., corresponding to the waves generated by streaming instability (fifth column in Table 1 with the name “Instability damping rate”), and the damping of waves generated by external sources parallel to the mean magnetic field (sixth column in Table 1 with the name “Wave damping rate”). The table illustrates that the rate of damping and the ranges of wavelengths for which damping is applicable are very different for the two situations. At the first glance, this seems strange. However, the difference stems from the fact that in the case of streaming instability the waves are aligned with the local magnetic field, while the waves generated by an extended source are sent parallel to the mean magnetic field.
We did not cover in Table 1 the general case of Alfvén waves generated at an arbitrary angle relative to the mean magnetic field, as well as damping of Alfvén waves by outer-scale turbulence. It is also necessary to stress again the important role of weak turbulence for the suppression of streaming instability at low MA. While the weak turbulence is frequently disregarded due to its limited inertial range
5 Ion-Neutral Collisional Damping of Streaming Instability
In the presence of partial ionization, an additional effect of damping by ion-neutral collisions becomes important. This effect was discussed originally by Kulsrud and Pearce [7] for Alfvén waves. The damping of turbulent motions in partially ionized gas was recently summarized in Xu and Lazarian [34].
In the presence of neutrals, a slippage between them and ions induces the dissipation. In a mostly neutral medium, at wave frequencies ω = VAk‖ less than the neutral-ion collisional frequency νni, both species move together and the dissipation is minimal. As the wave frequency increases, not all neutrals get the chance to collide with ions and the relative motions of ions and neutrals induce significant dissipation. For strongly coupled ions and neutrals, the ion-neutral collisional (IN) damping rate is [7].
where ξn = ρn/ρ, and ρn and ρ are the neutral and total mass densities. For weakly coupled ions and neutrals with ω = VAik‖ > νin, where VAi is the Alfvén speed in ions and νin is the ion-neutral collisional frequency, there is
We note that both turbulent and wave motions are subject to the IN damping. Strong Alfvénic turbulence injected in the strong coupling regime cannot cascade into the weak coupling regime due to the severe damping effect [35,36].
IN damping is sensitive to the ionization fraction and becomes weak at a high ionization fraction. For strongly coupled ions and neutrals with VAk‖ < νin, ΓIN is still given by Eq. 71. For decoupled ions with VAik‖ > νin, there is [36].
where χ = ρn/ρi and ρi is the ion mass density. Furthermore, when neutrals are also decoupled from ions with VAik‖ > νni, the above expression is reduced to Eq. 72. Because of the weak damping effect, Alfvénic cascade in a highly ionized medium is not dissipated by IN damping [18].
Naturally, to understand whether turbulent damping or IN damping is more important for damping the streaming instability, ΓIN should be compared with the turbulent damping rate Γ that we provided earlier. This comparison has been recently carried out in detail by Xu and Lazarian [18]. Here we selectively review some of their results.
In a weakly ionized interstellar medium, e.g., molecular clouds, CR-driven Alfvén waves are likely in the weak coupling regime with
where B0 is the mean magnetic field strength, ne and nH are number densities of electrons and atomic hydrogen, and ECR is the CR energy. As already mentioned above, strong Alfvénic turbulence injected at a large scale in the strong coupling regime is severely damped and its cascade cannot persist in the weak coupling regime. Therefore, there is
So the damping of streaming instability in a weakly ionized medium is dominated by IN damping.
In a highly ionized interstellar medium, e.g., the warm ionized medium, CR-generated Alfvén waves are still in the weak coupling regime and have
To have the turbulent damping dominate over IN damping, there should be
which can be rewritten as
for superAlfvénic turbulence, where ni and nn are the number densities of ions and neutrals, and
for subAlfvénic turbulence. We see that the condition in Eq. 78 is naturally satisfied for superAlfvénic turbulence. In a highly ionized medium, as the IN damping is weak, streaming instability is predominantly damped by the turbulent damping.
6 Numerical Testing of Turbulent Damping of Alfvén Waves
Numerical testing of Lazarian [26] is essential in a variety of regimes. By using 3D MHD turbulence simulations [33], the results of numerical testing on turbulent damping of externally driven Alfvén waves are presented in Figure 1. The observed scaling is consistent with Lazarian [26] predictions, but inconsistent with FG04 prediction.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/702799/fphy-10-702799-HTML-r1/image_m/fphy-10-702799-g001.jpg)
FIGURE 1. The damping time-scale Γ−1 of Alfvén waves that are injected at k‖ = 10 in 3D MHD turbulence, where the parallel direction is chosen with respect to the mean magnetic field. In one approach the Alfvén wave energy Ew decays in the turbulent medium over the time scale τ1 = ln(E(t1)/E(t2))/(t2 − t1). The values of τ1 are given by triangular symbols. In the other approach the wave energy is continuously injected at k‖ = 10 until it reaches a saturation level Ew. The corresponding damping time scale is given by τ2 = Ew/ϵdriving, where ϵdriving is the wave energy injection rate. τ2 is denoted by diamond symbols. The two measurements are both consistent with k−2/3 scaling. From Cho and Lazarian (in prep)’.
The reason for this difference arises from the global reference frame adopted in the numerical experiment. Launching of Alfvén waves with respect to the local direction of magnetic field is complicated in turbulent fluid. Therefore, the testing presented in Figure 1 was carried out with Alfvén waves launched with respect to the mean magnetic field. This is the setting corresponding to turbulent damping of Alfvén waves generated in the global system of reference that we considered in §4. As a result, the numerical simulations confirmed the scaling of inverse of damping rate Γ−1, i.e., damping time scale, which is measured at different λ as
7 Astrophysical Implications
7.1 Propagation of CRs
For decades the study on CR propagation was performed within a simple model, the so-called “leaky box model” (see [38]). In this model Galactic CRs propagate freely within the partially ionized disk of the Galaxy. The Alfvén waves experience damping in the partially ionized gas [7,34,36,39] and thus the streaming instability is suppressed. On the contrary, in fully ionized plasmas of the Galactic halo, the damping of Alfvén waves is significantly reduced and the streaming instability is present. Therefore, in this classical simplistic picture that ignores turbulence, Galactic CRs stream freely through the Galactic disk and are scattered backwards in the Galactic halo.
This classical “leaky box model” is problematic, as it is well known now that the Galactic disk is not fully filled with partially ionized gas. In fact, a significant fraction of the Galactic disk material is warm ionized gas [40,41]. Therefore, CRs cannot zoom through the Galactic disk due to the streaming instability.
FG04 quantified the idea of turbulent damping of streaming instability mentioned in Yan and Lazarian [23] and came to a paradoxical conclusion by applying their theory to the propagation of CRs in the Galaxy. By assuming homogeneous transAlfvénic turbulence in the Galaxy, they found significant turbulent damping of streaming instability and thus poor confinement of CRs. This would entail problems with explaining, e.g., the observed isotropy of CRs and their residence time in the Galaxy.
In Lazarian [26] the gist of the “leaky box model” was preserved, but instead of damping by ion-neutral collisional friction, the study appealed to the turbulent damping of streaming instability in the Galactic disk and proposed a “turbulent leaky box model.” Different from FG04, by considering inhomogeneous turbulence properties in the Galaxy and the strong MA dependence of turbulent damping, they found that the damping by weak subAlfvénic turbulence is marginal in the Galactic halo and thus CRs, even at high energies, can still be confined by streaming instability.
In a recent study by Xu and Lazarian [18], they identified the important role of turbulent damping of streaming instability in the warm ionized medium (WIM). Figure 2 shows the diffusion coefficient D of streaming CRs. The MA dependence comes from both turbulent damping of streaming instability and wandering of turbulent magnetic field lines. In particular, the smaller D in superAlfvénic turbulence is caused by the tangling of turbulent magnetic fields, which results in an effective mean free path lA of the CRs streaming along turbulent magnetic fields [42].
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/702799/fphy-10-702799-HTML-r1/image_m/fphy-10-702799-g002.jpg)
FIGURE 2. Diffusion coefficient vs. ECR of streaming CRs in super and subAlfvénic turbulence in the WIM. From Xu and Lazarian [18].
The MA-dependent diffusion of CRs is important for a realistic modeling of inhomogeneous CR diffusion in the Galaxy [43]. The actual values of MA in the Galaxy can be measured from observations using a newly developed gradient technique ([44], see also [45]) or with more traditional magnetic field and turbulent velocity measurements.
7.2 Launching of Winds and Heating
While the damping of Alfvén waves by turbulence is an accepted process in the field of CR research, we would like to point out that the turbulent damping of Alfvén waves can be responsible for many fundamental astrophysical processes. For instance, different processes of damping were discussed for heating of stellar corona by Alfvén waves, as well as for launching of stellar winds (see [19,46–50]). It is clear that the turbulent damping of Alfvén waves can be very important in these settings. More recently, launching galactic winds by turbulent damping of the Alfvén waves generated by galactic activity was considered in Suzuki and Lazarian [51]. Accounting for the dependence of turbulent damping on MA is important for the quantitative modeling of the process. A similar process is important for launching winds from other types of active disk systems, e.g., circumstellar disks.
Apart from launching galactic winds by turbulent damping of Alfvén waves generated by the galaxy, the turbulent damping of streaming instability also plays a very important role in coupling CRs and magnetized galactic matter. The pressure of CRs in galactic settings is significant and it can modify interstellar dynamics. Galactic winds driven by CRs present an important example of this modification.
In general, the importance of galactic winds is easy to understand. For galaxies of the Milky Way luminosity, about 20 percent of baryons are accounted for when matching the observed luminosity to the halo mass function. Observing absorption lines in spectra of background quasars testifies for the efficient expulsion of galactic baryons from the galaxies. In fact there is evidence that galaxies with significant star formation can drive mass outflows up to 10 times the rate of star formation [52].
Numerical simulations have demonstrated that CRs indeed influence the generation of global outflows and the local structure of the interstellar medium (ISM) (see [53]). The exact properties of the simulated outflows depend sensitively on how CR transport is modeled. Recent simulations by Holguin et al. [54] employed Lazarian [26] model of turbulent damping and obtained the results that differ significantly from the earlier modeling in, e.g., Ruszkowski et al. [53]. The difference stemmed from the fact that the earlier calculations employed the model by FG04, which is only applicable to transAlfvénic turbulence, i.e., MA = 1. However, the actual MA of gas can vary significantly in simulations.
The results of the numerical simulations in Holguin et al. [54] are presented in Figure 3. Some of the implications include, fist of all, when turbulent damping of CR streaming instability is included, there is an increase of star formation rate, and the increase is more significant at a higher level of turbulence. The reason is that the turbulent damping increases the average CR streaming speed. This allows CRs to leave the dense mid-plane, reducing the pressure support from CRs to the gas. As a result, the gas in the disk collapses and stars form more efficiently. Furthermore, the higher efficiency of star formation results in more CRs produced in the mid-plane. The increased streaming speed of CRs leads to a more extended CR distribution away from the mid-plane. It is also important that the escape of CRs from the dense regions allows them to interact with lower-density gas. This widens the gas distribution in height and accelerates the gas to form CR-driven galactic winds.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/702799/fphy-10-702799-HTML-r1/image_m/fphy-10-702799-g003.jpg)
FIGURE 3. Simulations of the galactic ISM evolution in the presence of star formation and CR driven outflows. The figure shows the gas (ni) and CR (ncr) density slices ± 5 kpc along z direction perpendicular to the midplane obtained in two simulations over time 200 Myr. The CR streaming is affected by turbulent damping of streaming instability with the turbulent velocity σ = 10 km/s. The results obtained in the absence of turbulent damping on the left side of each pair of plots are clearly different from those with turbulent damping on the right side. The distribution of both gas and CRs is more extended in the presence of turbulent damping. From Holguin et al. [54].
In addition, the theory of Alfvén wave damping by turbulence suggests that Alfvén waves can propagate across longer distances in highly magnetized regions of solar atmosphere (small MA) compared to the regions with higher MA. This prediction can be observationally tested. This effect should be accounted in both modelling of solar wind launching and modelling of plasma heating. For instance, it is likely that the turbulent damping can be important in order to explain the observed “unexpected” damping of Alfvén waves in the regions above the Sun’s polar coronal holes [55].
8 Summary
Alfvén waves are damped in turbulent media and the damping depends on the Alfvén Mach number MA of the turbulence. At the same wavelength, the wave damping depends on whether the waves are generated in the local reference system of magnetic eddies by the CR streaming or they are injected at an angle relative to the large-scale mean magnetic field from an extended astrophysical source. The latter is, e.g., the case of the Alfvén waves arising from magnetic reconnection, or oscillations in accretion disks and stellar atmospheres. The difference in their damping rates arises from the difference between the local and global systems of reference where the Alfvén waves are generated.
The dependence of damping rate on the wavelength λ of the Alfvén waves in the local system of reference is λ−1/2, as opposed to a stronger dependence λ−2/3 for the waves in the global reference system.
The turbulent damping also depends on whether Alfvén waves interact with weak or strong Alfvénic turbulence. For MA < 1, the turbulence from the injection scale L to the scale
In a partially ionized gas, the turbulent damping still dominates the damping of streaming instability when the ionization fraction is sufficiently high, e.g., in the warm ionized medium [18]. In star burst galaxies, the ionization fraction is low and the ion-neutral collisional damping can be more important (e.g., [56]).
The turbulent damping of streaming instability has important implications on propagation of CRs in the Galaxy, star formation, coupling between CRs and magnetized gas and thus driving galactic winds. In addition, the turbulent damping of Alfvén waves results in heating of the medium and transfer of the momentum from Alfvénic flux to the medium. The latter is also important for launching winds.
Author Contributions
The article was devised and written by AL and SX.
Funding
The research is supported by NASA TCAN 144AAG1967. SX acknowledges the support for this work provided by NASA through the NASA Hubble Fellowship grant # HST-HF2-51 473.001-A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Incorporated, under NASA contract NAS5-26 555.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
AL thank Jungyeon Cho for discussing the paper.
Footnotes
1Using the relation
2There is an intuitive way of presenting the Alfvénic cascade in terms of eddies mixing the magnetic field in the direction perpendicular to the magnetic field surrounding the eddies. The existence of such magnetic eddies is possible due to the fact that, as shown in LV99, the turbulent magnetic reconnection happens within one eddy turnover. As a result, the existence of magnetic field does not constrain magnetic eddies, if they are aligned with the magnetic field in their vicinity, i.e., with the local magnetic field. This eddy representation of MHD turbulence vividly demonstrates the importance of the local system of reference, where l⊥ and l‖ are defined.
3The fact that MHD turbulence is formulated in terms of the local quantities is required for describing the interaction of MHD turbulence with CRs. Indeed, perturbations in the local system of reference are exactly what CRs interact with.
4We note that in FG04 the injection scale for turbulence was defined not as the actual injection scale, but the scale at which the turbulent velocity becomes equal to the Alfvén one. Such scale does not exist for subAlfvénic turbulence.
5Simple estimates demonstrate that the interactions with smaller and larger turbulent scales are subdominant compared with the interaction with the “resonant” scale.
6In the case θ ∼ δB/B, one should average the final expressions over the θ dispersion that arises from magnetic field wandering.
References
1. Lazarian A, de Gouveia Dal Pino EM, Melioli C Magnetic Fields in Diffuse Media. In: Astrophysics and Space Science Library, Vol. 407. 978-3-662-44624-9. Berlin: Springer-Verlag Berlin Heidelberg (2015). doi:10.1007/978-3-662-44625-6
2. Uhlig M, Pfrommer C, Sharma M, Nath BB, Enßlin TA, Springel V Galactic Winds Driven by Cosmic ray Streaming. MNRAS (2012) 423:2374–96. doi:10.1111/j.1365-2966.2012.21045.x
3. Wiener J, Oh SP, Guo F Cosmic ray Streaming in Clusters of Galaxies. MNRAS (2013) 434:2209–28. doi:10.1093/mnras/stt1163
4. van der Holst B, Sokolov IV, Meng X, Jin M, Manchester WB, Tóth G, et al. Alfven Wave Solar Model (AWSoM): Coronal Heating. ApJ (2014) 782:81. doi:10.1088/0004-637x/782/2/81
5. Lynch BJ, Edmondson JK, Li Y Interchange Reconnection Alfvén Wave Generation. Sol Phys (2014) 289:3043–58. doi:10.1007/s11207-014-0506-x
6. Lerche I Unstable Magnetosonic Waves in a Relativistic Plasma. ApJ (1967) 147:689. doi:10.1086/149045
7. Kulsrud R, Pearce WP The Effect of Wave-Particle Interactions on the Propagation of Cosmic Rays. ApJ (1969) 156:445. doi:10.1086/149981
9. Skilling J Cosmic Rays in the Galaxy: Convection or Diffusion? ApJ (1971) 170:265. doi:10.1086/151210
10. Lazarian A, Beresnyak A Cosmic ray Scattering in Compressible Turbulence. MNRAS (2006) 373:1195–202. doi:10.1111/j.1365-2966.2006.11093.x
12. Suzuki TK Evolution of Solar-type Stellar Winds. Astron nachr (2013) 334:81–4. doi:10.1002/asna.201211751
13. Bell AR The Acceleration of Cosmic Rays in Shock Fronts - I. Monthly Notices R Astronomical Soc (1978) 182:147–56. doi:10.1093/mnras/182.2.147
14. Kulsrud RM Plasma Physics for Astrophysics / Russell M. Kulsrud. Princeton, N.J., USA: Princeton University Press (2005). (Princeton series in astrophysics).
17. Badruddin and Kumar A. Study of the Cosmic-Ray Modulation during the Passage of ICMEs and CIRs. Sol Phys (2016) 291:559–80. doi:10.1007/s11207-015-0843-4
19. Suzuki TK, Inutsuka S-i Making the Corona and the Fast Solar Wind: A Self-Consistent Simulation for the Low-Frequency Alfvén Waves from the Photosphere to 0.3 AU. ApJ (2005) 632:L49–L52. doi:10.1086/497536
21. Hopkins PF, Squire J, Chan TK, Quataert E, Ji S, Kereš S, et al. MNRAS (2021) 501:4184. doi:10.1093/mnras/staa3691
23. Yan H, Lazarian A Scattering of Cosmic Rays by Magnetohydrodynamic Interstellar Turbulence. Phys Rev Lett (2002) 89:281102. doi:10.1103/physrevlett.89.281102
24. Farmer AJ, Goldreich P Wave Damping by Magnetohydrodynamic Turbulence and its Effect on Cosmic‐Ray Propagation in the Interstellar Medium. ApJ (2004) 604:671–4. doi:10.1086/382040
25. Goldreich P, Sridhar S Toward a Theory of Interstellar Turbulence. 2: Strong Alfvenic Turbulence. ApJ (1995) 438:763. doi:10.1086/175121
26. Lazarian A Damping of Alfvén Waves by Turbulence and its Consequences: from Cosmic-ray Streaming to Launching Winds. ApJ (2016) 833:131. doi:10.3847/1538-4357/833/2/131
28. Lazarian A Enhancement and Suppression of Heat Transfer by MHD Turbulence. ApJ (2006) 645:L25–L28. doi:10.1086/505796
29. Lazarian A, Vishniac ET Reconnection in a Weakly Stochastic Field. ApJ (1999) 517:700–18. doi:10.1086/307233
31. Cho J, Vishniac ET The Anisotropy of Magnetohydrodynamic Alfvenic Turbulence. ApJ (2000) 539:273–82. doi:10.1086/309213
33. Cho J, Lazarian A, Vishniac ET Simulations of Magnetohydrodynamic Turbulence in a Strongly Magnetized Medium. Vishniac ETApJ (2002) 564:291–301. doi:10.1086/324186
34. Xu S, Lazarian A Magnetohydrodynamic Turbulence and Turbulent Dynamo in Partially Ionized Plasma. New J Phys (2017) 19:065005. doi:10.1088/1367-2630/aa6ec9
35. Xu S, Lazarian A, Yan H The Line Width Difference of Neutrals and Ions Induced by Mhd Turbulence. ApJ (2015) 810:44. doi:10.1088/0004-637x/810/1/44
36. Xu S, Yan H, Lazarian A Damping of Magnetohydrodynamic Turbulence in Partially Ionized Plasma: Implications for Cosmic Ray Propagation. ApJ (2016) 826:166. doi:10.3847/0004-637x/826/2/166
37. Cho J, Lazarian A Compressible Magnetohydrodynamic Turbulence: Mode Coupling, Scaling Relations, Anisotropy, Viscosity-Damped Regime and Astrophysical Implications. MNRAS (2003) 345:325–39. doi:10.1046/j.1365-8711.2003.06941.x
38. Longair MS In: MS Longair, editor. High Energy Astrophysics. Cambridge, UK: Cambridge University Press (2011).
39. Lithwick Y, Goldreich P Compressible Magnetohydrodynamic Turbulence in Interstellar Plasmas. ApJ (2001) 562:279–96. doi:10.1086/323470
40. McKee CF, Ostriker JP A Theory of the Interstellar Medium - Three Components Regulated by Supernova Explosions in an Inhomogeneous Substrate. ApJ (1977) 218:148. doi:10.1086/155667
41. Draine BT Physics of the Interstellar and Intergalactic Medium by Bruce T. Draine. 978-0-691-12214-4. Princeton, N.J., USA: Princeton University Press (2011).
42. Brunetti G, Lazarian A Compressible Turbulence in Galaxy Clusters: Physics and Stochastic Particle Re-acceleration. Monthly Notices R Astronomical Soc (2007) 378:245–75. doi:10.1111/j.1365-2966.2007.11771.x
45. Xu S, Hu Y Measuring Magnetization with Rotation Measures and Velocity Centroids in Supersonic MHD Turbulence. ApJ (2021) 910:88. doi:10.3847/1538-4357/abe403
46. Verdini A, Velli M, Oughton S Propagation and Dissipation of Alfvén Waves in Stellar Atmospheres Permeated by Isothermal Winds. A&A (2005) 444:233–44. doi:10.1051/0004-6361:20052748
47. Evans RM, Opher M, Jatenco-Pereira V, Gombosi TI Surface Alfvén Wave Damping in a Three-Dimensional Simulation of the Solar Wind. IApJ (2009) 703:179–86. doi:10.1088/0004-637x/703/1/179
48. Vidotto AA, Jatenco-Pereira V Alfvén Waves as a Driving Mechanism in Stellar Winds. Adv Space Res (2010) 46:509–13. doi:10.1016/j.asr.2008.12.025
49. Verdini A, Velli M, Matthaeus WH, Oughton S, Dmitruk P A Turbulence-Driven Model for Heating and Acceleration of the Fast Wind in Coronal Holes. ApJ (2010) 708:L116–L120. doi:10.1088/2041-8205/708/2/l116
53. Ruszkowski M, Yang H-YK, Zweibel E Global Simulations of Galactic Winds Including Cosmic-ray Streaming. ApJ (2017) 834:208. doi:10.3847/1538-4357/834/2/208
54. Holguin F, Ruszkowski M, Lazarian A, Farber R, Yang H-YK Role of Cosmic-ray Streaming and Turbulent Damping in Driving Galactic Winds. MNRAS (2019) 490:1271–82. doi:10.1093/mnras/stz2568
55. Hahn M, Landi E, Savin DW Evidence of Wave Damping at Low Heights in a Polar Coronal Hole. ApJ (2012) 753:36. doi:10.1088/0004-637x/753/1/36
Keywords: cosmic rays, magnetic fields, turbulence, interstellar medium, galactic winds
Citation: Lazarian A and Xu S (2022) Damping of Alfvén Waves in MHD Turbulence and Implications for Cosmic Ray Streaming Instability and Galactic Winds. Front. Phys. 10:702799. doi: 10.3389/fphy.2022.702799
Received: 30 April 2021; Accepted: 11 January 2022;
Published: 05 May 2022.
Edited by:
Luca Comisso, Columbia University, United StatesReviewed by:
Takeru Suzuki, The University of Tokyo, JapanBogdan Hnat, University of Warwick, United Kingdom
Copyright © 2022 Lazarian and Xu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alex Lazarian, alazarian@facstaff.wisc.edu