arxiv: v1 [astro-ph.sr] 27 Sep 2013

Transcription

1 Draft version September 30, 2013 Preprint typeset using L A TEX style emulateapj v. 5/2/11 MAGNETOACOUSTIC WAVES IN A PARTIALLY IONIZED TWO-FLUID PLASMA Roberto Soler 1, Marc Carbonell 2, & Jose Luis Ballester 1 1 Departament de Física, Universitat de les Illes Balears, E Palma de Mallorca, Spain and 2 Departament de Matemàtiques i Informàtica, Universitat de les Illes Balears, E Palma de Mallorca, Spain Draft version September 30, 2013 arxiv: v1 [astro-ph.sr] 27 Sep 2013 ABSTRACT Compressible disturbances propagate in a plasma in the form of magnetoacoustic waves driven by both gas pressure and magnetic forces. In partially ionized plasmas the dynamics of ionized and neutral species are coupled due to ion-neutral collisions. As a consequence, magnetoacoustic waves propagating through a partially ionized medium are affected by the ion-neutral coupling. The degree to which the behavior of the classic waves is modified depends on the physical properties of the various species and on the relative value of the wave frequency compared to the ion-neutral collision frequency. Here, we perform a comprehensive theoretical investigation of magnetoacoustic wave propagation in a partially ionized plasma using the two-fluid formalism. We consider an extensive range of values for the collision frequency, ionization ratio, and plasma β, so that the results are applicable to a wide variety of astrophysical plasmas. We determine the modification of the wave frequencies and study the frictional damping due to ion-neutral collisions. Approximate analytic expressions to the frequencies are given in the limit case of strongly coupled ions and neutrals, while numerically obtained dispersion diagrams are provided for arbitrary collision frequencies. In addition, we discuss the presence of cutoffs in the dispersion diagrams that constrain wave propagation for certain combinations of parameters. A specific application to propagation of compressible waves in the solar chromosphere is given. Subject headings: magnetic fields magnetohydrodynamics (MHD) plasmas Sun: atmosphere Sun: oscillations waves 1. INTRODUCTION The theoretical study of magnetohydrodynamic (MHD) waves in partially ionized plasmas has received increasing attention in the recent years. Energy dissipation due to ion-neutral collisions may be important in many astrophysical plasmas, e.g., in the low solar atmosphere, because of its connection to plasma heating (e.g., Khomenko & Collados 2012; Martínez-Sykora et al. 2012; Russell & Fletcher 2013). Some recent examples of theoretical works that focused on the investigation of MHD waves in partially ionized plasmas are, e.g., Pécseli & Engvold (2000); De Pontieu et al. (2001); Kumar & Roberts (2003); Khodachenko et al. (2004, 2006); Leake et al. (2005); Zaqarashvili et al. (2011, 2013); Soler et al. (2012) among many others. In the present paper, we continue the research started in Soler et al. (2013a), hereafter Paper I, about wave propagation in a partially ionized two-fluid plasma. The two-fluid formalism used here assumes that ions and electrons form together an ion-electron fluid, while neutrals form a separate fluid. The ion-electron fluid and the neutral fluid interact by means of ion-neutral collisions. Paper I was devoted to the study of Alfvén waves. Here we tackle the investigation of magnetoacoustic waves. The dynamics of magnetoacoustic waves is more involved than that of Alfvén waves even in the case of a fully ionized plasma. The reason is that Alfvén waves are incompressible and are only driven by the magnetic tension force. On the contrary, magnetoacoustic waves produce compression of the plasma and are driven by both magnetic and gas pressure forces (see, e.g., Goossens 2003). roberto.soler@uib.es, marc.carbonell@uib.es, joseluis.ballester@uib.es As a consequence, the behavior of magnetoacoustic waves is determined by the relative values of the sound and Alfvén velocities in the plasma. When the medium is partially ionized, the situation is even more complex because, on the one hand, magnetoacoustic waves also produce compression in the neutral fluid due to collisions with ions and, on the other hand, the neutral fluid is able to support its own acoustic waves that may couple with the plasma magnetoacoustic waves. Thus, the study of propagation of magnetoacoustic disturbances in a partially ionized medium is interesting from the physical point of view and challenging from the theoretical point of view. This work can be related to the previous papers by Zaqarashvili et al. (2011) and Mouschovias et al. (2011). Zaqarashvili et al. (2011) derived the dispersion relation for magnetoacoustic waves in a partially ionized twofluid plasma. They focused their study on the derivation of the two-fluid equations from the more general threefluid equations and in comparing the two-fluid results with those in the single-fluid approximation. They concluded that the single-fluid approximation breaks down when the value of the ion-neutral collision frequency is on the same order of magnitude or lower than the wave frequency. Here we proceed as in Zaqarashvili et al. (2011) and use the more general two-fluid theory instead of the single-fluid approximation. Mouschovias et al. (2011) also considered the two-fluid theory and performed a very detailed investigation of wave propagation and instabilities in interstellar molecular clouds. To realistically represent the plasma in molecular clouds, Mouschovias et al. (2011) included in addition to ion-neutral collisions other effects as, e.g., ionization and recombination and self-gravity. The main purpose

2 2 of the present work is to obtain physical insight on the effect of ion-neutral collisions on waves without taking into account other effects that could hinder the specific role of collisions. For this reason, our approach is closer to Zaqarashvili et al. (2011) than to Mouschovias et al. (2011). As stated before, Zaqarashvili et al. (2011) put their emphasis in the derivation of the equations and the dispersion relation. As examples, Zaqarashvili et al. (2011) computed some solutions, but they restricted themselves to a specific set of parameters and did not perform an in-depth parameter study. We follow the spirit of Paper I and do not focus on a specific astrophysical plasma. Instead, we take the relevant physical quantities as, e.g., the collision frequency, the ionization fraction, the Alfvén and sound velocities, etc., as free parameters whose values can be conveniently chosen. An inconvenience of the present approach is that the space of parameters to be explored is very big. The advantage is that this approach makes possible to apply the results of this work to a wide variety of situations. Compared to Zaqarashvili et al. (2011), we cover a wider parameter range. In particular, we consider lower ionization levels and weaker magnetic fields than those used by Zaqarashvili et al. (2011). The output of this investigation will be also useful to understand the recent results by Soler et al. (2013b) concerning the behavior of slow waves propagating along a partially ionized magnetic flux tube. This paper is organized as follows. Section 2 contains the description of the equilibrium state and the basic equations. In Section 3 we explore the limit cases of uncoupled fluids and strongly coupled fluids and derive approximate expressions to the wave frequencies in those cases. Later, we investigate the general case in Section 4 by plotting dispersion diagrams when the collision frequency is varied between the uncoupled and strongly coupled regimes. We perform a specific application to magnetoacoustic waves propagating in the low solar atmosphere in Section 5. Finally, Section 6 contains some concluding remarks. 2. EQUILIBRIUM AND BASIC EQUATIONS The equilibrium configuration considered here is a uniform and unbounded partially ionized hydrogen plasma composed of ions, electrons, and neutrals. As in Paper I, the plasma dynamics is studied using the two-fluid formalism (see also, e.g., Zaqarashvili et al. 2011). Ions and electrons are considered together as forming a single fluid, while neutrals are treated as a separate fluid. The ion-electron fluid and the neutral fluid interact by means of particle collisions, which transfer momentum between the species. It is implicitly assumed in the present version of the two-fluid theory that ion-electron and ion-ion collisions occur always more frequently than ion-neutral collisions, so that ions and electrons display fluid behavior regardless of the presence of neutrals. The medium is permeated by a straight and constant magnetic field, B. We use Cartesian coordinates and conveniently choose the reference frame so that the equilibrium magnetic field is orientated along the z-direction, namely B = Bẑ, with B constant. We ignore the effect of gravity. The omission of gravity is approximately valid as long as we consider wavelengths much shorter than the gravitational stratification scale height. We also assume that there are no flows in the equilibrium. We study linear adiabatic perturbations superimposed on the equilibrium state. The equations governing the behavior of the perturbations are the same as in Paper I, namely ρ i v i t = p i + 1 µ ( b) B α in(v i v n ), (1) ρ n v n t = p n α in (v n v i ), (2) b t = (v i B), (3) p i t = γp i v i, (4) p n t = γp n v n, (5) where v i, p i, P i, and ρ i are the velocity perturbation, gas pressure perturbation, equilibrium gas pressure, and equilibrium density of the ion-electron fluid, v n, p n, P n, and ρ n are the respective quantities for the neutral fluid, b is the magnetic field perturbation, µ is the magnetic permeability, γ is the adiabatic index, and α in is the ionneutral friction coefficient. The friction coefficient, α in, plays a very important role. It determines the strength of the ion-neutral friction force and, therefore, the importance of ion-neutral coupling. Braginskii (1965) gives the expression of α in in a hydrogen plasma. As in Paper I, we take α in as a free parameter and do not use a specific expression of α in. We do so to conveniently control the strength of the ion-neutral friction force. Our aim is to keep the investigation as general as possible for the results of the present article to be applicable in a wide range of astrophysical situations. Instead of using α in, in the remaining of this article we write the equations using the ion-neutral, ν in, and neutral-ion, ν ni collision frequencies, which have a more obvious physical meaning. They are defined as ν in = α in ρ i, ν ni = α in ρ n. (6) Hence, both frequencies are related by ρ i ν in = ρ n ν ni Normal mode analysis Here we perform a normal mode analysis. The temporal dependence of the perturbations is put proportional to exp( iωt), where ω is the angular frequency. In addition, we perform a Fourier analysis of the perturbations in space. The spatial dependence of perturbations is put proportional to exp(ik x x+ik y y +ik z z), where k x, k y, and k z are the components of the wavevector in the x-, y-, and z-directions, respectively. Alfvén waves and magnetoacoustic waves are two distinctclassesofmhdwavesinauniformplasma(see, e.g., Cramer 2001; Goossens 2003). Alfvén waves are incompressible and propagate vorticity perturbations. Magnetoacoustic waves are compressible and do not produce vorticity perturbations. Alfvén waves were investigated in Paper I. Here we focus on magnetoacoustic waves. Compressibility is an appropriate quantity to describe magnetoacoustic waves (see, e.g., Lighthill 1960). By working with compressibility perturbations we are able to decouple magnetoacoustic waves from Alfvén waves without imposing any restriction on the direction of the

3 Magnetoacoustic waves in a two-fluid plasma 3 wavevector. We define i and n as the compressibility perturbations of the ion-electron fluid and the neutral fluid, respectively, i v i = ik x v x,i +ik y v y,i +ik z v z,i, (7) n v n = ik x v x,n +ik y v y,n +ik z v z,n (8) We combine Equations (1) (5) and after some algebraic manipulations we obtain the two following coupled equations involving i and n only, namely ( ω 4 ω 2 k 2( c 2 i +c 2 ) A +k 2 kzc 2 2 Ac 2 ) i i = iν in ω 3 ( i n ) iν in + ω +i(ν ni +ν in ) k2 kzc 2 2 ( A c 2 i i c 2 ) n n, (9) ( ω 2 k 2 c 2 ) n n = iν ni ω( n i ), (10) where k 2 = k 2 x +k2 y +k2 z and c A, c i, and c n are the Alfvén velocity, the ion-electron sound velocity, and the neutral sound velocity defined as c 2 A = B2 µρ i, c 2 i = γp i ρ i, c 2 n = γp n ρ n. (11) Note that the Alfvén velocity is here defined using the density of the ionized fluid only. Equations (9) and (10) are the governing equations for linear compressibility perturbations. Therefore, they govern the behavior of magnetoacoustic waves. The magnetoacoustic character of the waves is determined by the relative values of the Alfvén and sound velocities. In this work we use the following definition of the ion-electron β i, β i = c2 i c 2 = γp i A B 2 /µ. (12) Thus, using the present definition of β i we can quantify the relative importance of the Alfvén velocity, c A, and ion-electron sound velocity, c i. No restriction on the value of the neutral sound velocity, c n, has been imposed so far. However, to simplify matters and decrease the number of parameters involved, we proceed as in Zaqarashvili et al. (2011) and Soler et al. (2013b) and assume that there is a strong thermal coupling between the ion-electron fluid and the neutral fluid. As a consequence, the equilibrium temperature, T, is the same for both fluids. We use the ideal gas law to relate the equilibrium gas pressure, density, and temperature as P i = 2ρ i RT, P n = ρ n RT, (13) where R is the ideal gas constant. The factor 2 in the expression for ion-electrons is present because the partial pressures of both electrons and ions have to be summed up. In a hydrogen plasma, the partial pressures of electrons and ions are the same. From Equation (13), we deduce that the requirement that all the species have the same temperature results in imposing a relation between the sound velocities as, c 2 i = 2c 2 n. (14) Hence, using Equations (12) and (14) we see that fixing the value ofβ i isenough toestablish the relationbetween the Alfvén velocity and the two sound velocities. For the subsequent analysis, we decompose the wavevector, k, in its components parallel and perpendicularto the equilibrium magnetic field as k = k z ẑ+k 1, where k = kx 2 +ky 2 and 1 is the unity vector in the perpendicular direction to B. We define θ as the angle that forms the wavevector, k, with the equilibrium magnetic field, B. Then, we can write k z = kcosθ and k = ksinθ. For propagation perpendicular to the magnetic field, θ = π/2, whereas for propagation parallel to the magnetic field, θ = 0. Intermediate values of θ represent oblique propagation. 3. APPROXIMATIONS FOR LIMIT VALUES OF THE COLLISION FREQUENCIES Before tackling the general case for arbitrary values of the collision frequencies, it is instructive to consider two paradigmatic limit cases, namely the uncoupled case in which the two fluids do not interact, and the strongly coupled case in which the two fluids behave as a single fluid Uncoupled case Here, we remove the effect of ion-neutral collisions. We represent this situation by setting ν in = ν ni = 0. Equations (9) and (10) become ( ω 4 ω 2 k 2( c 2 i +c 2 ) A +k 4 c 2 Ac 2 i cos 2 θ ) i =0, (15) ( ω 2 k 2 c 2 ) n n =0. (16) Now, the equations governing compressibility perturbations in the two fluids are decoupled. The waves in the two fluids can be studied separately. We start by considering waves in the ion-electron fluid. For nonzero i, the solutions to Equation (15) must satisfy ω 4 ω 2 k 2( c 2 i +c2 A) +k 4 c 2 A c2 i cos2 θ = 0. (17) Equation (17) is the well-known dispersion relation of magnetoacoustic waves in a fully ionized plasma (see, e.g., Lighthill 1960). The solutions to Equation (17) are ω 2 = k 2c2 A +c2 i 2 ±k 2c2 A +c2 i 2 [ 1 4c2 A c2 i cos2 θ (c 2 A +c2 i )2 ] 1/2, (18) where the + sign is for the fast wave and the sign is for the slow wave. The magnetoacoustic character of these two distinct wave modes depends on the relative values of c A and c i (see, e.g., Goossens 2003). In the low-β i case, c 2 A c2 i and the solutions simplify to ω 2 k 2( c 2 A +c 2 i) k 2 c 2 A, (19) for the fast wave and ω 2 k 2 c2 A c2 i c 2 A +c2 i cos 2 θ k 2 c 2 i cos2 θ, (20) for the slow wave. When c 2 A c2 i the fast wave is a magnetic wave that propagates isotropically at the Alfvén velocity, c A, while the slow waveis essentiallyasound wave that is guided by the magnetic field and travels at the

4 4 ion-electron sound velocity, c i. Conversely, in the high-β i case c 2 i c 2 A and the solutions simplify to ω 2 k 2( c 2 ) A +c2 i k 2 c 2 i, (21) for the fast wave and ω 2 k 2 c2 A c2 i c 2 cos 2 θ k 2 c 2 A +c2 A cos2 θ, (22) i for the slow wave. Now the fast wave behaves as an isotropic sound wave, while the slow wave is a magnetic wave guided by the magnetic field. Weturntothewavesintheneutralfluid. Thesolutions to Equation (16) for nonzero n are ω 2 = k 2 c 2 n. (23) Equation (23) represent acoustic (or sound) waves in a gas. These waves propagate isotropically at the neutral sound velocity, c n, and are unaffected by the magnetic field. As expected, in the absence of ion-neutral collisions we consistently recover the classic magnetoacoustic waves in the ion-electron fluid and the classic acoustic waves in the neutral fluid. Therefore, three distinct waves are present in the uncoupled, collisionless case. These waves do not interact and are undamped in the absence of collisions. The waves would be damped if nonideal processes as, e.g., viscosity, resistivity, or thermal conduction, are taken into account. In the remaining of this article, we use the adjective classic to refer to the waves found in the uncoupled case Strongly coupled case Conversely to the uncoupled case, the strongly coupled limit represents the situation in which ion-electrons and neutrals behave as a single fluid. To study this case, we take the limits ν in and ν ni in Equations (9) and (10). We realize that, if ω 0, it is necessary that i = n forthe equationstoremainfinite. Thisisequivalent to assume that the two fluids move as a whole. Then, when ν in and ν ni, Equations (9) and (10) reduce to a single equation, namely ( ω 4 ω 2 k 2c2 A +c2 i +χc2 n 1+χ ( ) + k 4c2 A c 2 i +χcn) 2 (1+χ) 2 cos 2 θ i,n = 0, (24) where we use i,n to represent either i or n and the ionization fraction of the plasma, χ, is defined as χ = ρ n ρ i. (25) For nonzero i,n, the solutions to Equation (24) must satisfy ( ω 4 ω 2 k 2c2 A +c2 i +χc2 n +k 4c2 A c 2 i +χcn) 2 1+χ (1+χ) 2 cos 2 θ = 0. (26) This is the wave dispersion relation in the strongly coupled case. Its solutions are ω 2 =k 2c2 A +c2 i +χc2 n 2(1+χ) ±k 2c2 A +c2 i +χc2 n 2(1+χ) [ 1 4c2 A (c2 i +χc2 n)cos 2 θ (c 2 A +c2 i +χc2 n) 2 ] 1/2 (27), where the + sign is for the modified fast wave and the sign is for the modified slow wave. We use the adjective modified to stress that these waves are the counterparts of the classic fast and slow modes but modified by the presence of the neutral fluid. The first important difference between the uncoupled and strongly coupled cases is in the number of solutions. In the uncoupled case there are three distinct waves, namely the classic slow and fast magnetoacoustic waves and the neutral acoustic wave, but in the strongly coupled limit we only find the modified version of the slow and fast magnetoacoustic modes. The modified counterpart of the classic neutral acoustic mode is apparently absent. The question then arises, what happened to the neutral acoustic mode? To answer this question, we study in the following paragraphs the character of the modified fast and slow modes depending on the relative values of the Alfvén and sound velocities. A remark should be made before tackling the approximate study of the solutions in Equation (27). We recall that the neutral and ion-electron sound velocities are related by c 2 i = 2c 2 n. This condition means that c2 i and c 2 n are of the same order of magnitude. However, this is not the end of the story. We notice that the value of the ionization fraction, χ, also plays a role in determining the solutions to Equation (27). Specifically, χ appears multiplying c 2 n, which means that χ can increase or decrease the effective value of the neutral sound velocity. Thus, we need to compare χc 2 n with c2 A and c2 i in order to determine the nature of the solutions. The approximate study of the solutions given in Equation (27) is done considering eight typical cases: 1. Case c 2 A c2 i χc 2 n. This case corresponds to an almost fully ionized (χ 1) low-β i plasma. The solutions given in Equation (27) simplify to for the modified fast wave and ω 2 k 2 c 2 A, (28) ω 2 k 2 c 2 i cos 2 θ, (29) forthemodified slowwave. We recoverthe fastand slow waves of the uncoupled case (Equations (19) and (20)). No trace of the neutral acoustic mode remains in this limit because the amount of neutrals is negligible. 2. Case c 2 A c2 i χc 2 n. This case corresponds to a low-β i plasma with χ 1. The solutions are approximated as ω 2 k 2 c 2 A 1+χ, (30) for the modified fast wave and ω 2 k 2c2 i +χc2 n 1+χ cos2 θ, (31) for the modified slow wave. By comparing Equations (19) and (30), we see that the square of

5 Magnetoacoustic waves in a two-fluid plasma 5 the fast mode frequency is reduced by the factor (χ+1) 1 compared to the value in the uncoupled case. As a consequence, the effective Alfvén velocity in the strongly coupled limit is c A / 1+χ. This is equivalent to replace the ion density, ρ i, by the sum of the ion and neutral densities, ρ i +ρ n, in the definition of the Alfvén velocity. The same result was obtained by Kumar & Roberts (2003) and Soler et al. (2013b) for fast magnetoacoustic waves in the case c 2 n = 0. This is also the same result obtained in Paper I and Soler et al. (2012) for Alfvén waves. On the other hand, the modification of the slow mode frequency is more complicated(compare Equations (20) and (31)). The expression of the slow mode frequency involves both c i and c n. Consequently, an effective sound velocity, c eff, can be defined as c 2 eff = c2 i +χc2 n 1+χ = ρ ic 2 i +ρ nc 2 n ρ i +ρ n. (32) The effective sound velocity is the weighted average of the sound velocities of ion-electrons and neutrals, and was first obtained by Soler et al. (2013b) for slow magnetoacoustic waves propagating along magnetic flux tubes. The fact that the expression of the effective sound velocity (Equation (32)) involves both neutral and ion-electron sound velocities suggests that, in this case, the modified slow mode is the descendant of both the classic slow mode and the neutral acoustic mode. 3. Case c 2 A χc2 n c 2 i. This case corresponds to a weakly ionized (χ 1) low-β i plasma, but the magnetic field is strong enough for c 2 A to be much larger than χc 2 n. The solutions given in Equation (27) become for the modified fast wave and ω 2 k 2c2 A χ, (33) ω 2 k 2 c 2 n cos2 θ, (34) for the modified slow wave. The effective Alfvén velocity is now c A / χ. This is equivalent to replace the ion density, ρ i, by the neutral density, ρ n, in the definition of the Alfvén velocity. On the other hand, the slow wave is the descendant of the neutral acoustic wave that becomes guided by the magnetic field. Therefore, the frequencies of both fast and slow waves are determined by the density of the neutral fluid alone, which is indirectly affected by the magnetic field. 4. Case χc 2 n c 2 A c2 i. In this case we consider a very weakly ionized (χ ) low-β i plasma. The approximate frequencies turn to be for the modified fast wave and ω 2 k 2 c 2 n, (35) ω 2 k 2c2 A χ cos2 θ, (36) for the modified slow wave. The modified fast wave can straightforwardly be related to the isotropic acoustic mode of the neutral fluid (Equation (23)), while the modified slow wave is a guided magnetic wave with the effective Alfvén velocity depending onthe neutraldensity alone. Asin case(3), neither the ion-electron sound velocity nor the ion density play a role. The fact that the plasma is weakly ionized changes completely the physical nature of the solutions compared to the classic magnetoacoustic waves in a low-β i fully ionized plasma. 5. Case c 2 i c 2 A χc2 n. This corresponds to an almost fully ionized (χ 1) high-β i plasma. The solutions given in Equation (27) simplify to for the modified fast wave and ω 2 k 2 c 2 i, (37) ω 2 k 2 c 2 Acos 2 θ, (38) for the modified slow wave. We recover Equations (21) and (22) obtained in the high-β i collisionless case. No trace of the neutral acoustic mode remains in this limit. 6. Casec 2 i χc 2 n c2 A. This correspondsagainto an almost fully ionized (χ 1) high-β i plasma, but now the magnetic field is so weak that χc 2 n c 2 A. The solutions given in Equation (27) simplify to for the modified fast wave and ω 2 k 2 c 2 i, (39) ω 2 k 2 c 2 A cos2 θ, (40) for the modified slow wave. As in the previous case (5), we revert to Equations (21) and (22) obtained in the uncoupled case. Cases (5) and (6) point out that the relative values of c 2 A and χc2 n are not important as long as c 2 i remains much larger than both of them. The same approximations to the frequencies are found in cases (5) and (6). The properties of the neutral fluid are not important in both cases (5) and (6). 7. Case χc 2 n c 2 i c 2 A. This is a high-β i plasma with χ 1. The solutions given in Equation (27) simplify to ω 2 k 2c2 i +χc2 n 1+χ = k2 c 2 eff, (41) for the modified fast wave and ω 2 k 2 c 2 A 1+χ cos2 θ, (42) for the modified slow wave. Equations (41) and (42) are the modified version of Equations (21) and (22) obtained in the high-β i uncoupled case, where c eff replaces c i and c A / 1+χ replaces c A. Here we find that the modified fast mode is the descendant of both the classic fast mode and the neutral acoustic mode, while the modified slow mode is a guided magnetic wave whose frequency depends on both ion and neutral densities.

6 6 8. Caseχc 2 n c2 i c 2 A. Thiscorrespondstoaweakly ionized (χ 1) high-β i plasma. The solutions given in Equation (27) simplify to for the modified fast wave and ω 2 k 2 c 2 n, (43) ω 2 k 2c2 A χ cos2 θ, (44) for the modified slow wave. The approximation found here are the same as in case (4). The properties of the neutral fluid completely determine the behavior of the waves. Neutrals feel the magnetic field indirectly through the collisions with ions. The results of the above approximate study are summarized in Table 1. Based on this study, we can now answer the previous question about the apparent absence of the neutral acoustic mode in a strongly coupled plasma. We conclude that the classic neutral acoustic mode and the classic ion-electron magnetoacoustic modes heavily interact when ion-neutral collisions are at work. The two resulting modes in the strongly coupled regime (modified fast mode and modified slow mode) have, in general, mixed properties and are affected by the physical conditions in the two fluids. The degree to which the properties of the classicwavesare present in the resulting waves depends on the relative values of the Alfvén and sound velocities and on the ionization fraction of the plasma. In addition, the wave frequencies in the strongly coupled limit are real as in the uncoupled case. This means thatthewavesareundamped inthelimitofhighcollision frequencies as well. The damping of magnetoacoustic waves due to ion-neutral collisions takes place for intermediate collision frequencies. As stated before, the waves would be damped if additional damping mechanisms are considered, but this falls beyond the aim of the present article. Both the ion-neutral damping and the coupling between modes are investigated in Section 4 by studying the modification of the wave frequencies when a progressive variation of the collision frequencies between the uncoupled limit and strongly coupled limit is considered. 4. NUMERICAL STUDY FOR ARBITRARY COLLISION FREQUENCIES From here on, we consider arbitrary values of ν in and ν ni. With no approximation, we combine the coupled Equations (9) and (10) to obtain separate equations for i and n as with D(ω) i =0, (45) iν ni ω D(ω) D n (ω) n=0. (46) D(ω) = D i (ω)d n (ω)+dc 2 (ω), (47) where D i (ω), D n (ω), and Dc 2 (ω) are defined as D i (ω)=ω 3 (ω +iν in ) ω 2 k 2( c 2 A +c 2 ) i ω +iν ni + ω +i(ν in +ν ni ) k4 c 2 Ac 2 i cos 2 θ, (48) D n (ω)=ω(ω +iν ni ) k 2 c 2 n, (49) D 2 c (ω)= ων in ν ni [ ω 3 (ω +i(ν in +ν ni )) ω +i(ν in +ν ni ) k 4 c 2 A c2 n cos2 θ ]. (50) The waves in the ion-electron fluid and in the neutral fluid are described by the same dispersion relation. The dispersion relation is D(ω) = 0, (51) and is equivalent to the dispersion relation given by Zaqarashvili et al. (2011) in their Equation (57), although here we use a different notation. Expressed in polynomial form, the dispersion relation is a 7th-order polynomial in ω, so that there are seven solutions (see Zaqarashvili et al. 2011). The frequencies are complex, namely ω = ω R +iω I, where ω R and ω I are the real and imaginary parts. The imaginary part of the frequency is negative and accounts for the exponential damping rate of the perturbations due to ion-neutral collisions. Before exploring the solutions to the dispersion relation, it is necessary to make a comment on the nature of the solutions. We distinguish between two different kind of solutions depending of whether ω R is zero or nonzero. (1) The first kind of solutions are those with ω R = 0, i.e., solutions that represent purely decaying disturbances. These solutions are not waves in the strict sense because they do not propagate. The solutions with ω R = 0 were called vortex modes by Zaqarashvili et al. (2011) but, in our view, the name vortex modes does not reflect the true nature of these solutions. The reason is that vortex modes in a partially ionized plasma are intrinsically related to Alfvén waves (see Paper I), not to magnetoacoustic waves. Magnetoacoustic modes do not produce vorticity perturbations (see, e.g., Goossens 2003), hence the purely imaginary solutions of the magnetoacoustic dispersion relation cannot represent fluid vorticity. Instead, we propose to physically interpret the solutions with ω R = 0 as entropy modes (see, e.g., Goedbloed & Poedts 2004; Murawski et al. 2011). Entropy modes represent perturbations of density and, therefore, compression of the plasma. Compression is an intrinsic property of magnetoacoustic modes. Entropy modes have zero frequency in the absence of collisions. Due to ion-neutral collisions, perturbations in density are damped and the entropy mode frequency acquires an imaginary part, while its real part remains zero. (2) The secondkindofsolutionsarethosewith ω R 0. Thesesolutions represent propagating waves and appear in pairs, namely ω 1 = ω R + iω I and ω 2 = ω R + iω I. Both ω 1 and ω 2 represent the same magnetoacoustic wave, but the real parts of their frequencies have opposite signs (see Carbonell et al. 2004). The different signs of the real parts account for forward (ω R > 0) and backward (ω R < 0) propagationwith respect to the direction of the wavevector, k. Since the equilibrium is static, the two directions of propagation are physically equivalent. As shown later, the number of entropy modes and propagating waves is not constant. Depending on the physical parameters considered, there can be conversion between these two kind of solutions. Unlike the limit cases studied in Section 3, no simple analytic solutions of the dispersion relation can be obtained when ν in and ν ni are arbitrary. Instead, we solve the dispersion relation numerically. The numerically ob-

7 Magnetoacoustic waves in a two-fluid plasma 7 tained dispersion diagrams are plotted as functions of the averaged collision frequency, ν, defined as ν = ρ iν in +ρ n ν ni 2 = α in. (52) ρ i +ρ n ρ i +ρ n Because of the complex behavior displayed by the solutions and the large number of parameters involved, we organize the presentation of the results and their discussion as follows. First, we investigate the specific situation in which the wave propagation is strictly perpendicular to the magnetic field direction. This paradigmatic case is helpful to get physical insight of the behavior of the various solutions and to illustrate some general results that repeatedly appear in this investigation. Slow modes are absent from the discussion for perpendicular propagation because they are unable to propagate across the magnetic field. Subsequently, we study the other limit case, i.e., wave propagation parallel to the magnetic field, so that slow modes are added to the discussion. Finally, the general case of oblique propagation is taken into account Perpendicular propagation We start by studying the case of perpendicular propagation to the magnetic field. Hence we set θ = π/2. The dispersion relation(equation(51)) reduces to a 5th-order polynomial because the slow magnetoacoustic wave is absent. Conversely, the fast magnetoacoustic mode and the neutral acoustic mode can propagate across the field and so they remain when θ = π/2. We compute the real and imaginary parts of the frequency of the various solutions of the dispersion relation as functions of ν. Frequencies are expressed in units of kc i so that the results can be applied to the astrophysical plasma of interest by providing appropriate numeric values of k and c i. The averaged collision frequency is varied four orders of magnitude between ν/kc i = 0.01 and ν/kc i = 100. These two values are chosen to represent the uncoupled and strongly coupled limits, respectively. Beyond these two values of ν/kc i the behavior of the solutions remains essentially unaltered. Figures 1 and 2 show the results for β i = 0.04 (low-β i case) and β i = 25 (high-β i case), respectively. Three values of χ are considered in each case, namely χ = 0.2, 2, and 20, which cover the cases of largely ionized, moderately ionized, and weakly ionized plasmas, respectively. First, we analyze the real part of the frequency. The comparisonof the numericallyobtained ω R with the analytic approximations of Section 3 helps us identifying the various propagating waves in the limit values of ν. The modes are labeled and plotted with different line styles in Figures 1 and 2. When ν/kc i = 0.01 we find the classic fast magnetoacoustic (solid black line) and neutral acoustic (dashed blue line) modes, while for ν/kc i = 100 only the modified fast wave is present as a propagating solution, in agreement with Section 3. As expected, the modes do not display the properties of the classic waves when the averaged collision frequency takes intermediate values. To understand how the solutions change due to the effect of collisions we continuously follow the various solutions when ν/kc i increases from ν/kc i = 0.01 to ν/kc i = 100. We find that one of the two solutions present at the low collision frequency limit has a cutoff when ν/kc i reaches a certain value. At the cutoff, the forward (ω R > 0) and backward (ω R < 0) waves of this specific solution merge and their ω R becomes zero. On the contrary, the other solution does not have a cutoff and eventually becomes the modified fast wave at the high collision frequency limit. The location of the cutoff is important because it determines the number of distinct waves that are able to propagate for a certain combination of parameters. The effect of the various parameters on the location of the cutoff can be studied as in Paper I by taking advantage of the fact that the dispersion relation is a polynomial in ω. So, we can compute the polynomial discriminant of the dispersion relation. By definition, the discriminant is zero when the dispersion relation has a multiple solution. Since the cutoff takes place when a forward and a backward propagating wave merge, the dispersion relation has necessarily a double root at the cutoff. The discriminant is a function of χ, ν/kc i, and β i. Hence, the zeros of the discriminant inform us about the relation between these parameters at the cutoff. We omit here the expression of the discriminant because it is rather cumbersome and can be straightforwardly obtained from the dispersion relation using standard algebraic methods. The number of propagating solutions for a given set of parameters is represented in Figure 3 in the χ ν/kc i plane for three different values of β i. The zero of the discriminant is plotted with a line. Below this line two propagating waves are possible, namely the classic fast and acoustic modes, whileabovethelineonlythemodifiedfastwaveis present as a propagating wave. Importantly, we find that a third region, where propagating waves are not possible, appears when β i 1 (see the shaded area in Figure 3 when β i = 10 3 ). This forbidden region takes place for large ionization ratio (weakly ionized plasmas), very low β i (strong magnetic fields) and relatively high collision frequency (strongly coupled plasmas). This could be the case of the plasma in intense magnetic tubes of the low solar atmosphere. The physical reason for the existence of wave cutoffs is qualitativelythe sameasdiscussedinpaperiforthe case of Alfvén waves (see also, e.g., Kulsrud & Pearce 1969; Pudritz 1990; Kamaya & Nishi 1998; Mouschovias et al. 2011). There is competition between the ion-neutral friction force and the wave restoring forces. Friction becomes the dominant force and overcomes the remaining forces at the cutoff. Hence, waves are unable to propagate when the strength of friction greatly exceeds those of the restoring forces. The cutoff found here corresponds to the situation in which ions and neutrals are sufficiently coupled for the neutral fluid to be significantly affected by the magnetic field. As a consequence, the neutral fluid isnot able to support its ownpure acousticmodes, which start to behave as slow-like magnetoacoustic waves. For θ = π/2 slow modes cannot propagate and so these solutions are cut off. We turn to the imaginary part of the frequency. For presentationpurposes, we plot the absolute value ofω I in logarithmic scale. The line styles used in the graphics of ω I are the same as in the plots of ω R so that the various solutions can be easily identified. The inspection of ω I in Figures 1 and 2 informs us of the existence of a third purely imaginary mode (red dotted line). This mode is the descendant of the undamped entropy mode, which

8 8 becomes damped by collisions when ν/kc i increases. The damping rate of this nonpropagating mode monotonically increases with ν/kc i. Regarding the propagating waves, we see that the solution that has the cutoff in ω R displays a bifurcation in ω I at the same location. Two purely imaginary branches emerge at the bifurcation. We identify these new branches as entropy modes. On the one hand, the upper branch is heavily damped and follows the behavior of the original entropy mode plotted with a red dotted line. On the other hand, the lower branch is less damped and becomes an undamped entropy mode when ν/kc i. Finally, the remaining solution whose ω I does not have a bifurcation corresponds to the propagating wave that becomes the modified fast wave when ν/kc i 1. The damping rate of this wave is nonmonotonic. Its ω I is maximal, i.e., the damping is most efficient, when ν ω R. This is the same result as for Alfvén waves (see Paper I) Parallel propagation We move to the case of parallel propagation to the magnetic field. We set θ = 0 and perform the same computations as in the previous section. These results are displayed in Figures 4 and 5, where we have considered the same parameters as before. As expected, the slow magnetoacoustic wave is now present (dash-dotted green line). As in the perpendicular propagation case, we analyze Figures 4 and 5 by observing how the modes present at the low collision frequency limit evolve as we increase the collision frequency toward the high collision frequency limit. In the low-β i case (Figure 4), the slow magnetoacoustic wave and the neutral acoustic wave strongly interact. There is a change of character between these two waves depending on the value of χ. When χ 1, i.e., largely ionized plasma, the neutral acoustic mode has a cutoff and the classic slow mode becomes the modified slow mode in the limit of large ν/kc i. On the contrary, the situationisreversedwhenχ 1,sothatthesolutionthat is cut offis the classicslow mode while the neutral acoustic mode is the solution that becomes the modified slow mode. These results are consistent with the approximate study of Section 3.2. On the other hand, the behavior of the fast magnetoacoustic wave in the low-β i case is very similar to that of the Alfvén waves studied in Paper I. It is well-known that in a low-β i plasma fast waves behave as Alfvén waves for parallel propagation to the magnetic field (see, e.g., Goossens 2003). As in the case of Alfvén waves, fast magnetoacoustic waves are nonpropagating in an interval of ν/kc i when χ > 8 (see extensive details in Paper I for the case of Alfvén waves). This nonpropagating interval can be seen in Figure 4 in the results for χ = 20. In this interval the fast wave becomes a purely imaginary solution. The boundaries of the nonpropagating interval can be obtained from Equation (20) of Paper I. Using the present notation, the boundaries of the nonpropagating interval are [ χ 2 +20χ 8 8(χ+1) 3 3/2 ] 1/2, ν 2χ = ± χ1/2 (χ 8) kc i β 1/2 i (χ+1) 8(χ+1) 3 (53) where the and + signs correspond to the left and right boundaries of the nonpropagating range. In the high-β i case (Figure 5), the slow magnetoacoustic wave displays little interaction with the other two solutions. In this case, the behavior of the fast magnetoacoustic wave and the neutral acoustic wave is similar to that found when θ = π/2, while the slow wave behaves as an Alfvén wave (Paper I). Now, it is the slow wave and not the fast wave that has a nonpropagating interval of ν/kc i. Equation (53) holds for the slow wave in the high-β i case too. Since Equation (53) depends upon β 1/2 i, the nonpropagating interval is now shifted towards smaller ν/kc i than in the low-β i case (compare the results for χ = 20 in Figure 4 and 5). We study in more detail the presence of cutoffs and nonpropagating intervals by using again the discriminant of the dispersion relation. This is shown in Figure 6 for β i = 0.04 and β i = 25. As in the case of perpendicular propagation (see Figure 3) we find two zones in the χ ν/kc i planeseparatedbyacutoff(blacksolidline). These separate zones correspond to the region where the classic slow, fast, and acoustic waves live and the region where the modified slow and fast modes exist. In addition, we find the presence of the nonpropagating interval previously discussed (area between the red lines). In the lowβ i case(figure6a) the nonpropagatingintervalislocated above the boundary between the two regions. Within the nonpropagating interval, the modified slow wave is the only propagating solution. On the contrary, in the high-β i case (Figure 6b) the nonpropagating interval intersects the boundary between the two regions. Hence the number of possible scenarios increases, as schematically indicated in Figure 6b. It is appropriate to compare our results with those plotted in Figure 3 of Zaqarashvili et al. (2011) in order to check whether the present work is consistent with that previous study. In our notation, the parameters used by Zaqarashvili et al. (2011) are χ = 1, β i = 0.25, and θ = 0. These parameters are close to those used in Figure 4(c) (d). To perform a proper comparison we should use the same notation as Zaqarashvili et al.(2011). They used two dimensionless quantities, w and a, that correspond to the wave frequency and to the inverse of the collision frequency, respectively. These two quantities are expressed in our notation as βi (1+χ) ω w=, (54) cosθ kc i a= 2cosθ βi (1+χ) kc i ν. (55) We solve the dispersion relation using the above parameters and reproduce in Figure 7 the results shown in Figure 3 of Zaqarashvili et al. (2011). A perfect agreement between our results and those of Zaqarashvili et al. (2011) is obtained. With the help of the discriminant of the dispersion relation, we compute that the cutoff of the neutral acoustic wave takes place at ν/kc i 0.84, which corresponds to a 3.4 in the notation of Zaqarashvili et al. (2011). Since χ < 8, the fast mode does not have a nonpropagating interval. Thus, our results confirm and fully agree with the specific example studied by Zaqarashvili et al. (2011) Oblique propagation

9 Magnetoacoustic waves in a two-fluid plasma 9 Finally, here we consider the case in which the direction of wave propagation is in between the limit cases studied in the previous two Subsections. Hence, we use θ = π/4. The results are displayed in Figures 8 and 9. The behavior of the various solutions is more complex than in the previous limits and can be understood as a combination of the results for perpendicular and parallel propagation. Compared to the case with θ = 0, the acoustic and fast modes display similar behavior, while the slow mode is the solution whose behavior is most altered, specially in thelow-β i case(comparefigures4and8). Bothslowand fast modes display nonpropagating intervals for oblique propagation, whereas only the fast mode have a nonpropagating interval for parallel propagation. The fast mode nonpropagating interval is also approximately described by Equation (53), which points out that the location of this interval does not depend on θ. Weproceedasbeforeandshowin Figure10thevarious possible scenarios for wave propagation in the χ ν/kc i plane of parameters. The high-β i case (Figure 10b) is verysimilartotheequivalentresultforθ = 0(Figure6b), so that no additional comments are needed. However, in the low-β i case, the presence of nonpropagating intervals for both slow and fast modes and their intersections cause the presence of six possible scenarios (Figure 10a). The result represented in Figure 10a clearly points out the high complexity of the interactions between the various waves in a partially ionized plasma, which do not occur in fully ionized plasmas (Mouschovias et al. 2011). The existence of two nonpropagating intervals in the low-β i case is consistent with the results by Soler et al. (2013b) of slow waves propagating along a partially ionized magnetic flux tube. Soler et al. (2013b) also found the presence of two nonpropagating intervals, see their Figure 5. Here we see that this double nonpropagating interval is a result also present for magnetoacoustic waves in a homogeneous, partially ionized medium. Thus, we can now conclude that this result is not caused by the geometry of the waveguide, but it is an intrinsic property of magnetoacoustic waves propagating obliquely to the magnetic field in a partially ionized plasma. The slow magnetoacoustic waves in a magnetic flux tube studied by Soler et al. (2013b) behave as the waves studied here for oblique propagation. 5. APPLICATION TO THE SOLAR CHROMOSPHERE Here we perform a specific application to magnetoacoustic waves propagating in the partially ionized solar chromosphere. This application is motivated by the recent observations of ubiquitous compressive waves in the low solar atmosphere (Morton et al. 2011, 2012). To represent the chromospheric plasma we use the same simplified model considered in Paper I. We use the quiet sun model C of Vernazza et al. (1981), hereafter VALC model, to account for the variation of physical parameters with height from the photosphere to the base of the corona. The expression of α in is taken after Braginskii (1965), namely α in = 1 2 ρ i ρ n 16kB T σ in, (56) m i πm i wherem i istheion(proton)mass,k B isboltzmann sconstant, and σ in isthe collisioncrosssection. Equation(56) considers only hydrogen and ignores the influence of heavier species. In Equation (56) it is implicitly assumed that the ion and neutral masses are approximately equal. For the collision cross section of protons with neutral hydrogen atoms we consider the typical value of σ in m 2 (used in, e.g., Khodachenko et al. 2004; Leake et al. 2005; Arber et al. 2007; Soler et al. 2009; Khomenko & Collados 2012, among others). Here we must note that, in the application done in Paper I for the case of Alfvén waves, we took σ in m 2 based on the hard sphere collision model (see, e.g., Braginskii 1965; Zaqarashvili et al. 2013). In a recent paper, Vranjes & Krstic (2013) claim that the realistic value of σ in in the solar chromosphere is about two orders of magnitude higher than that estimated in the hard sphere model. The value σ in m 2 used here is closer to the value proposed by Vranjes & Krstic (2013) than to the hard sphere value. We adopt the magnetic field strength model by Leake & Arber (2006), which aims to represent the field strength in a chromospheric vertical magnetic flux tube expanding with height, namely B = B ph ( ρ ρ ph ) 0.3, (57) where ρ = ρ i + ρ n is the total density, and B ph and ρ ph = kg m 3 are the magnetic field strength and the total density, respectively, at the photospheric level. The variation of ρ with height is taken from the VALC model. Accordingly, the magnetic field strength decreases with height. Regarding the value of B ph, we consider two possible scenarios: an active region and the quiet Sun. For the active region case, we set B ph = 1.5 kg, so that B 100 G at 1,000 km and B 20 G at 2,000 km above the photosphere. For the quiet Sun case, we set B ph = 100 G, so that B 7 G at 1,000 km and B 1 G at 2,000 km above the photosphere. As explained in Paper I, the physical parameters in this simplified model of the low solar atmosphere depend on the vertical direction, while in the previous theoretical analysis all the parameters are taken constant. Here we perform a local analysis and use the dispersion relation derived for a homogeneous plasma. We use the physical parameters at a given height to locally solve the dispersion relation at that height. A limitation of the present approach is that it ignores the possible presence of cutoff frequencies due to gravitational stratification (see, e.g., Roberts 2006). The condition for this method to be approximately valid is that the wavelength, λ = 2π/k, is much shorter than the stratification scale height. This is fulfilled by the wavelengths used in this analysis. Figure 11 displays the variation with height above the photosphere of the three physical quantities relevant for the behavior of the waves: the ionization fraction, χ, the ionized fluid β i, and the averaged ion-neutral collision frequency, ν. Only β i is affected by the magnetic field scenario considered. First of all, we see that χ ranges several orders of magnitude from very weakly ionized plasma at the low levels of the chromosphere to fully ionized plasma when the transition region to the solar corona is reached. Full ionization of hydrogen takes place