# Modulation of waves due to charge-exchange collisions in magnetized partially ionized space plasma

## Abstract

A nonlinear time dependent fluid simulation model is developed that describes the evolution of magnetohydrodynamic waves in the presence of collisional and charge exchange interactions of a partially ionized plasma. The partially ionized plasma consists of electrons, ions and a significant number of neutral atoms. In our model, the electrons and ions are described by a single fluid compressible magnetohydrodynamic (MHD) model and are coupled self-consistently to the neutral gas, described by the compressible hydrodynamic equations. Both the plasma and neutral fluids are treated with different energy equations that describe thermal energy exchange processes between them. Based on our self-consistent model, we find that propagating Alfvénic and fast/slow modes grow and damp alternately through a nonlinear modulation process. The modulation appears to be robust and survives strong damping by the neutral component.

###### pacs:

96.25.Qr, 96.50.Dj,96.50.Bh,96.50.Tf,96.50.Xy,94.20.wf,94.30.cq## I Introduction

Alfvén, fast and slow mode waves are characteristics of magnetohydrodynamic (MHD) plasma. While the electromagnetic fluctuations in Alfvén waves propagate predominantly along an ambient or guide constant magnetic field, the fast/slow modes propagates isotropically. Although the linear properties of these waves are fairly well understood Shukla (1999); Shukla & Stenflo (1995); Shukla & Feix (1989), their nonlinear evolution in a complex partially ionized environment, as occurs often in space/astrophysical plasmas, remains relatively uninvestigated. In partially ionized (consisting of ions, electrons and a significant neutral gas presence) space and astrophysical plasmas, these waves can interact with the neutral gas through processes such as charge exchange, collisions, ionization and recombination etc. Hence the wave properties are expected to be modified. It is however unclear how Alfvén, fast and slow waves modify the nonlinear processes that govern plasma fluctuations in the presence of neutral dynamics. In the context of cosmic ray propagation, Kulsrud & Pearce (1969) noted that the interaction of a neutral gas and plasma can damp Alfvén waves. Neutrals interacting with plasma via a frictional drag process results in ambipolar diffusion (Oishi & Mac Low 2006). Ambipolar diffusion plays a crucial role in the dynamical evolution of the near solar atmosphere, interstellar medium, and molecular clouds and star formation. Oishi & Mac Low (2006) found that ambipolar diffusion can set a characteristic mass scale in molecular clouds. They found that structures less than the ambipolar diffusion scale are present because of the propagation of compressive slow mode MHD waves at smaller scales. Leake et al (2005) showed that the lower chromosphere contains neutral atoms, the existence of which greatly increases the efficiency of wave damping due to collisional friction momentum transfer. They noted that Alfvén waves with frequencies above 0.6Hz are completely damped and frequencies below 0.01 Hz are unaffected. They undertook a quantitative comparative study of the efficiency of the role of (ion-neutral) collisional friction, viscous and thermal conductivity mechanisms in damping MHD waves in different parts of the solar atmosphere. They found that a correct description of MHD wave damping requires the consideration of all energy dissipation mechanisms through the inclusion of the appropriate terms in the generalized Ohmâs law, the momentum, energy and induction equations. Padoan et al (2000) calculated frictional heating by ion-neutral (or ambipolar) drift in turbulent magnetized molecular clouds and showed that the ambipolar heating rate per unit volume depends on field strength for constant rms Mach number of the flow, and on the Alfvénic Mach number.

The role of ion-neutral collisions has been investigated extensively by Balsara (1996) in the context of molecular clouds. Here the momentum equation for the plasma is governed predominantly by the slow comoving massive neutrals that tend to dissipate the Cloud’s magnetic field. Because of the larger neutral to ion mass ratio, ion momentum is entirely dominated by the neutral drag. Hence the inertial terms were ignored in plasma momentum equation. Balsara (1996) found that slow waves propagate without significant damping on short wavelengths, while the fast and Alfvén waves undergo rapid damping in super-Alfvénic regimes.

In the heliospheric boundary regions, particularly the outer heliosheath Zank (1999); Shaikh & Zank (2008), the plasma is mediated by interstellar neutral gas by virtue of charge exchange. For example in the local interstellar medium (LISM), the low density plasma and neutral hydrogen (H) gas are coupled primarily through the process of charge exchange. In this process, both the total (proton+atom) momentum and energy are conserved. On sufficiently large temporal and spatial scales, a partially ionized plasma is typically regarded as equilibrated; this is the case for the LISM, but not for the outer heliosheath, for example.

It is worth mentioning that the charge exchange process in its simpler (and leading order) form can be treated like a friction or viscous drag term in the fluid momentum equation, describing the relative difference in the ion and neutral fluid velocities. The drag imparted in this manner by a collision between ion and neutral also causes ambipolar diffusion, a mechanism used to describe Alfvén wave damping by cosmic rays (Kulsrud & Pearce 1969) and also discussed by Oishi & Mac Low (2006) in the context of molecular clouds.

In this Letter, we focus on understanding the propagation characteristics of MHD waves by considering collisional and charge exchange interactions thus including both simultaneously. Our model is more general and fully nonlinear, compared to that of Balsara (1996) which ignores nonlinear intertial terms in the plasma momentum equation. We retain this term and carry out fully time-dependent simulations in two dimensions. We treat plasma and neutral fluid in a self-consistent manner by employing the complete nonlinear time-dependent fluid equations for both the fluids. In Section 2, we discuss the equations of a coupled plasma-neutral model, their validity, the underlying assumptions and the normalizations. Section 3 describes results of our nonlinear, coupled, self-consistent fluid simulations. We find that the propagation characteristics of Alfvén waves are altered signifincantly by the combined action of charge exchange and collisional interaction processes. The latter leads to a modulated growth and damping of the Alfvén and fast/slow modes. In Section 4, we summarize our results and describe physical processes that may be responsible for the modulation process.

## Ii Partially ionized Model Equations

We assume that fluctuations in the plasma and neutral fluids are isotropic, homogeneous, thermally equilibrated and turbulent. A constant mean magnetic field is present. Local mean flows and nonlinear structures may subsequently be generated by self-consistently excited nonlinear instabilities. The boundary conditions are periodic. This allows us to assume an infinite partially ionized plasma where the length scale of characteristic fluctuations () is much smaller than the size of plasma or computational box () i.e. . We further assume that plasma particles interact with the neutral gas via collisions as well as charge exchange. The intrinsic magnetized waves supported by the plasma scatter plasma particles in a random manner thus sufficiently isotropizing the plasma and neutral distribution functions Monin & Yaglom (2007). The latter enables us to use fluid descriptions for both the plasma and neutral gas. Most of these assumptions are appropriate to realistic space and astrophysical turbulent flows. They allow us to use MHD and hydrodynamic descriptions for the plasma and the neutral components respectively. In the context of the LISM and the outer heliosheath plasmas, the plasma and neutral fluid remain close to thermal equilibirium and behave as Maxwellian fluids. Our model simulates a plasma-neutral fluid that is coupled via collisions and charge exchange in astrophysical plasmas. The fluid model describing nonlinear turbulent processes, in the presence of charge exchange and collision, can be described in terms of the plasma density (), velocity (), magnetic field (), pressure () components according to the conservative form

(1) |

where,

and

Note the presence of the source terms that couple the plasma self-consistently to the neutral gas. The above set of plasma equations is supplimented by and is coupled self-consistently to the neutral density (), velocity () and pressure () through a set of hydrodynamic fluid equations,

(2) |

where,

Equations (1) to (2) form an entirely self-consistent description of the coupled plasma-neutral turbulent fluid in a partially ionized medium.

Several points are worth noting. The charge-exchange momentum sources in the plasma and the neutral fluids, i.e. Eqs. (1) & (2), are described respectively by terms and . These expressions, described in Pauls et al (1995); Shaikh & Zank (2008), have the following form.

(3) |

(4) |

with

A swapping of the plasma and the neutral fluid velocities in this representation corresponds, for instance, to momentum changes (i.e. gain or loss) in the plasma fluid as a result of charge exchange with the neutral atoms (i.e. in Eq. (1)). Similarly, momentum change in the neutral fluid by virtue of charge exchange with the plasma ions is described by in Eq. (2).

The plasma-neutral collisional forces are modeled by employing the following forcing term in Eqs. (1) & (2).

(5) |

where is the collision frequency. The collisions modify plasma and neutral momentum and energy, but not the density. The negative sign before in Eq. (5) ensures conservation of total momentum. We derive modified equations of energy for both fluids following the treatment described in Landau & Lifshitz (1981). The collisional force is self-consistently calculated in our model. Note also the charge exchange source term in momentum and energy of both the fluids. In the absence of collision and charge exchange interactions, the plasma and the neutral fluid are de-coupled trivially and behave as ideal fluids. Like collisions, charge-exchange interactions modify the momentum and the energy of plasma and the neutral fluids and conserve density in both the fluids (since we neglect photoionization and recombination). Nonetheless, the volume integrated energy and the density of the entire coupled system will remain conserved in a statistical manner. The conservation processes can however be altered in the presence of any external forces. These can include large-scale random driving of turbulence due to external forces such as supernova explosions, stellar winds, etc and instabilities. Finally, the magnetic field evolution is governed by the usual induction equation, i.e. Eq. (1), that obeys the frozen-in-field theorem unless some nonlinear dissipative mechanism introduces small-scale damping.

The underlying coupled fluid model can be non-dimensionalized straightforwardly using a typical scale-length (), density () and velocity (). The normalized plasma density, velocity, energy and the magnetic field are respectively; . The corresponding neutral fluid quantities are . The momentum and the energy charge-exchange terms, in the normalized form, are respectively . The non-dimensional temporal and spatial length-scales are . Note that we have removed bars from the set of normalized coupled model equations (1) & (2). The charge-exchange cross-section parameter (), which does not appear directly in the above set of equations, is normalized as , where the factor has dimension of (area). By defining through , we see that there exists a charge exchange mode () associated with the coupled plasma-neutral turbulent system. For a characteristic density, this corresponds physically to an area defined by the charge exchange mode being equal to (mpf). Thus the larger the area, the higher is the probability of charge exchange between plasma ions and neutral atoms. Therefore, the probability that charge exchange can directly modify those modes satisfying is high compared to modes satisfying . Since the charge exchange length-scales are much smaller than the turbulent correlation scales, this further allows many turbulent interactions amongst the nonlinear turbulent modes before they undergo at least one charge exchange. An exact quantitive form of sources due to charge exchange in our model is taken from Shaikh & Zank (2008).

## Iii Evolution of MHD waves

We now investigate the nonlinear evolution of the coupled plasma-neutral system described by Eqs. (1) & (2). Our goal is to understand the evolution of MHD waves in a partially ionized turbulent plasma, mediated by the complex coupling between the two distinguishable fluids. For this purpose, we initialize the field variables (i.e. velocity, magnetic, density, energy and pressure fields) in the plasma and neutral fluids with a random, uncorrelated, out-of-phase distribution in space. The initial ratio of the plasma and neutral density is unity. The plasma and neutral fluids initialized in this manner characterize spatially local turbulence in the partially ionized astrophysical plasma. These random fields develop spatially and temporally according to the set of Eqs. (1) & (2) and forces due to charge exchange and collision dictate nonlinear interaction processes. In ideal MHD and hydrodynamic fluids (i.e. without the coupling forces in Eqs. (1) & (2)), known quadratic conserved quantities like energy, helicity, vorticity etc predominantly govern the spectral transfer of energy across disparate characteristic scales in the inertial range. Their role is well studied within the paradigm of statistical theories of turbulence Kolmogorov (1941); Iroshnikov (1963); Kraichnan (1965); Monin & Yaglom (2007); Biskamp (2003). By contrast, the complex nature of the coupling forces in Eqs. (1) & (2) poses a formidable hurdle to deriving quadratic conserved quantities. It is not clear to us as how these quantities will influence spectral transfer in the inertial range. Furthermore, the driving forces in Eqs. (1) & (2) may presumably be operative at any length and time scales. It is unclear how this will alter the conventional dissipative and diffusive processes in a partially ionized plasma. Unlike the constant drag force (i.e. not changing in space and time) used before for a neutral fluid; and which acts to damp Alfvén waves in the partially ionized plasma (Kulsrud & Pearce 1969, Balsara 1996, Shaikh 2010), it is not clear how the complex, nonlinear, time and space dependent driving forces in Eqs. (1) & (2) will modify the propagation and interactions characteristic of Alfvén, fast and slow MHD waves or even compressive sound waves in the coupled neutral gas.

To address the above issues, we have developed a two-dimensional (2D) nonlinear fluid code to numerically integrate Eqs. (1) & (2). The 2D simulations are not only computationally less expensive (compared to a fully 3D calculation), but they offer significantly higher resolution (to compute inertial range turbulence spectra) even on moderately-sized clusters like our local Beowulf system. The spatial discretization in our code uses a discrete Fourier representation of turbulent fluctuations based on a spectral method, while we use a 4th-order Runge-Kutta method for the temporal integration. All the fluctuations are initialized isotropically with random phases and amplitudes in Fourier space. A mean constant magnetic field is assumed along the -direction. Our algorithm ensures conservation of total energy and mean fluid density per unit time in the absence of charge exchange and external random forcing. Additionally, is satisfied at each time step. Our code is massively parallelized using Message Passing Interface (MPI) libraries to facilitate higher resolution. The initial isotropic turbulent spectrum of fluctuations is chosen to be close to with random phases in all three directions. The choice of such a (or even a flatter than -2) spectrum does not influence the dynamical evolution as the final state in our simulations progresses towards fully developed turbulence.

From our time-dependent nonlinear fluid simulations, we determine the quantative evolution of Alfvén and fast/slow waves in a partially ionized environment. For this purpose, we distinguish Alfvénic and non-Alfvénic, i.e. corresponding to the compressional or due to slow and fast magnetosonic modes, contributions to the turbulent velocity fluctuations. Thus, we introduce diagnostics that distinguish between the modes corresponding to Alfvénic and slow/fast magnetosonic fluctuations. Since the Alfvénic fluctuations are transverse Biskamp (2003); Shaikh & Zank (2008), the propagation wave vector is orthogonal to the velocity field fluctuations i.e. , and the average spectral energy contained in these (shear Alfvénic modes) fluctuations can be computed as Shaikh & Zank (2006)

(6) |

The above relationship yields a finite spectral contribution from the characteristic turbulent Alfvénic modes. It is essentially a modal energy associated with the Alfvénic fluctuations that obey Eq. (6) in the inertial range turbulent spectrum. On the other hand, fast/slow (i.e. compressive) magnetosonic modes propagate longitudinally along the velocity field fluctuations, i.e. and thus carry a finite component of energy corresponding only to the part of the velocity field, which can be determined from the following relationship

(7) |

The expression for essentially describes the modal energy contained in the non-solenoidal component of the MHD turbulent modes.

The evolution of an average mode associated with Alfvénic and compressive waves is shown in Fig. (1). We vary charge exchange strength and collision parameter in our simulation to examine their effects on the propagation of Alfvénic and compressive modes in a partially ionized plasma. Clearly, the two processes operate on different time and length scales and are self-consistently modeled in our simulations. We find that charge exchange and collisional interactions jointly modulate Alfvén and fast/slow modes. This is shown in Fig. (1) for modes in a two dimensional box of length . The other parameters in our simulations are charge exchange wavenumber , fixed time step , and collision parameter . The background normalized constant magnetic field . Our simulations are fully nonlinear and the ratio of the mean and fluctuating magnetic fields .

The modulation process depicted in Fig. (1) can be understood on the basis of charge exchange and collisional interactions, together with convective nonlinearities. During the modulation, there exists a rapid onset of the average mode associated with Alfvénic waves. This onset is triggered essentially by linear instability process. Our linear stability analysis (not described here) indicates that the underlying coupled plasma-neutral system possesses several unstable modes. These modes account for the growing as well as damping of Alfvénic and fast/slow compressive waves. As illustrated in Fig. (1), a concurrent damping of the fast/slow compressive mode preceding the growth of Alfvénic mode occurs. When the two modes reach their extremal rise or fall, they reverse their behavior. The fast/slow compressive mode begins to rise at the expense of damping of Alfvénic mode. This process continues non periodically and repeats itself in time thereby exhibiting a nearly oscillatory behavior. It follows from Fig. (1) that the two modes, i.e. Alfvénic and fast/slow compressive, regulate each other in a predator-prey manner. Dynamically, the Alfvénic mode grows at the expense of fast/slow mode. When Alfvénic mode reaches its maximum amplitude, where the fast/slow mode remains at its lowest magnitude, the latter eats up the Alfvénic mode and vice versa. This process is repetitive. It is also noteworthy that the overall amplitude of both the modes decays in time and regrows later.

The growth and damping of wave energy can further be understood from the energy transfer rates () between the plasma and neutral fluids. This is shown in Fig. (2). It is evident from the figure that the total energy of the plasma increases at the expense of the neutral. When the waves damp, the total energy of the plasma decreases while it increases for the neutrals. The energy transfer rates also exhibit a periodic behavior following the modulation of energy in Alfvén and fast/slow modes. Nonetheless, the energy transfer rates try to approach nearly zero i.e. , but are not entirely zero due to the modulation process. This is an interesting result that differs from pure (i.e. uncoupled) plasma or neutral fluid turbulence Monin & Yaglom (2007); Biskamp (2003). The latter exhibits thereby leading to a constant transfer of energy across inertial range fluctuations. This further demonstrates that the total energy of the neutral and plasma fluid is nearly constant during the steady state (See Fig 2). What we find in our simulations is a modulated-energy-transfer process that is completely different from the constant energy transfer rate. The energy transfer rates in Fig. (2) thus provide a self-consistent description of the modulation process observed in our simulations.

## Iv Conclusion

In summary, we developed a two dimensional self-consistent model of plasma and neutral fluids that is coupled through charge exchange and collisional interactions. Notably, we find that these interactions not only modify the convective nonlinear interactions in the coupled Eqs. (1) & (2), but they also influence the propagation characteristics of MHD waves significantly. To determine the evolution of characteristic MHD waves, we distinguish the Alfvénic and non-Alfvénic contributions of MHD modes in the partially ionized turbulent environment. It is found that charge exchange and collisional interactions lead to a modulation of Alfvénic and fast/slow compressive modes. This is unlike purely collisional processes described by Kulsrud & Pearce (1969) and Balsara (1996) that tend to damp Alfvénic modes. By contrast, inclusion of the charge exchange interactions leads to an entirely different scenario and it alters the propagation characteristics of MHD waves. Charge exchange, in the presence of collisions and convective nonlinear interactions, leads to an alternate growth and damping (i.e. modulation) of Alfvénic and fast/slow compressive modes. It is further noted that the modulation of waves is mediated by the conservation of momentum because of the characteristic wave velocity that is involved in Eqs. (6) & (7). It is the momentum transfer function () that is responsible for regulating the rise/fall of the amplitude of the two waves and maintains the conservation process. This point is further consistent with Fig (1). We finally comment on the conservation of momentum and energy in our simulations. It seems from Fig (2) that the term hovers around zero, but does not fully converge to zero. We learn that due to the small scale self-consistent dissipation (or numerical dissipation), it is difficult to observe a 100% conservation. We continue to work to improve this conservation in our model.

Our results should find application to a variety of astrophysical environments in which a partially ionized plasma is typical. Examples include the outer heliosheath formed by the interactions of the solar wind with the local interstellar medium, the magnetic collapse of molecular clouds and star formation and the general transfer of energy in partially ionized plasmas surrounding other stellar systems Wood et al (2007).

## Acknowledgments

We acknowledge the partial support of NASA grants NNX09AB40G, NNX07AH18G, NNG05EC85C, NNX09AG63G, NNX08AJ21G, NNX09AB24G, NNX09AG29G, and NNX09AG62G.

### Footnotes

- Email : dastgeer.shaikh@uah.edu

### References

- Balsara, D., ApJ., 1996, 465, 775.
- Biskamp, D. ’Magnetohydrodynamic Turbulence’ Published by Cambridge University Press, 2003
- Iroshnikov, P. S., Astron. Zh. 40, 742 (1963).
- Kolmogorov, A. N. Dokl. Acad. Sci. URSS 30, 301 (1941).
- Kraichnan, R. H., Phys. Fluids 8, 1385 (1965).
- Kulsrud, R., and Pearce, 1969, ApJ 156, 445.
- Leake, J., Arber, T. D., M. L. Khodachenko, M. L., 2005, Astronomy and Astrophysics, 442, 1091.
- Monin, A. S. & Yaglom, A. M. Statistical Fluid Mechanics: Mechanics of Turbulence, Volume 1 & 2, Dover Publication, Inc., 2007
- Oishi, J. S. and Mac Low, M., 2006, ApJ, 638, 281.
- Padoan, P., Zweibel, E., and Nordlund, A., 2000, ApJ, 540, 332.
- Pauls, H. L., Zank, G. P., and Williams, L. L., 1995, J. Geophys. Res. A11, 21595.
- Pitaevskii, L. P., and Lifshitz, E. M., in Landau and Lifshitz Course of theoretical Physics, Volume 10, Physical Kinetics, Elsevier Inc. 1981.
- Shaikh, D., and Zank, G. P., 2008, The Astrophysical Journal, Volume 688, Issue 1, pp. 683-694.
- Shaikh, D., and Zank, G. P., 2006, The Astrophysical Journal, Volume 640, Issue 2, pp. L195-L198.
- Shaikh, D. and Shukla, P. K., 2009, Physical Review Letters, vol. 102, Issue 4, id. 045004
- Shukla, P. K.; Stenflo, L.; Bingham, R. Physics of Plasmas, Volume 6, Issue 5, pp. 1677-1680 (1999).
- Shukla, P. K.; Stenflo, L. Phys. Scr. T, Vol. T60, p. 32 - 35, 1995
- Shukla, P. K.; Feix, G., Solar Physics (ISSN 0038-0938), vol. 123, no. 1, 1989, p. 117-125, 1989.
- Wood, B. E., Harper, G. M.; Muller, H. R., Heerikhuisen, J., Zank, G. P., ApJ, 655, 946, 2007.
- Zank, G. P., 1999, Space Sci. Rev., 89, 413-688, 1999.