arxiv: v1 [physics.plasm-ph] 3 Nov 2015

Transcription

1 Mon. Not. R. Astron. So. 000, Printed 4 November 2015 MN LATEX style file v2.2 Charge exhange in fluid desription of partially ionized plasmas arxiv: v1 [physis.plasm-ph] 3 Nov 2015 J. Vranjes 1, M. Kono 2 and M. Luna 1 1 Instituto de Astrofísia de Canarias, La Laguna, Tenerife, Spain and Departamento de Astrofísia, Universidad de la Laguna, La Laguna, Tenerife, Spain 2 Faulty of Poliy Studies, Chuo University, Tokyo, Japan Aepted 1988 Deember 15. Reeived 1988 Deember 14; in original form 1988 Otober 11 1 INTRODUCTION Charge exhange also known as harge transfer, or eletron apture is a proess in whih an ion and an atom exhange eletrons in the proess of ollision. In noble gasses the harge exhange is more dominant than the usual elasti sattering Raizer In hydrogen plasmas, like in the ase of the lower solar atmosphere, the harge exhange takes plae at a very high rate as well Vranjes & Krsti 2013; Vranjes & Kono 2014; see more about it further in the text. Nevertheless, in the past the harge exhange has been ompletely ignored in the solar plasma literature. This in spite of the fat that the theory has been well known for more than 60 years Dalgarno 1958; Raizer 1991; Brandsen & MDowell 1991; Eihler 2005, and it has been thoroughly studied in the atomi physis and in the general plasma theory Krsti & Shultz 1998, 1999a,b. In fat, to the best of our knowledge, our two referenes given above are jvranjes@yahoo.om kono@fps.huo-u.a.jp mluna@ia.es ABSTRACT The effets of harge exhange on waves propagating in weakly ionized plasmas are disussed. It is shown that for low-frequeny proesses, ions and neutrals should be treated as a single fluid with some effetive harge on all of them. We have derived a new momentum equation whih should be used in suh an environment. As a result, the low-frequeny magneti waves an propagate even if partiles are not magnetized, whih is entirely due to the harge exhange and the fat that it is not possible to separate partiles into two different populations as harged and neutral speies. So there an be no frition fore between ions and neutrals in the usual sense. The mean fore per partile is proportional to the ionization ratio n i /. Regarding the appliation of the theory to the Alfvén wave propagation in the lower solar atmosphere, the results predit that the plane of displaement of the fluid must hange by 90 degrees when an Alfvén wave propagates from the area where partiles are un-magnetized photosphere to the area where they are magnetized hromosphere. Beause of the most aurate ross setions whih we have here, it is possible to very aurately determine altitudes at whih suh rotation of the Alfvén wave takes plae. Key words: Plasmas; waves; Sun: photosphere; Sun: hromosphere the only published works where the harge exhange has been studied in detail in the ontext of the solar plasma. In what follows we are going to show that the harge exhange may be a ruial fator in explaining the propagation of waves in weakly ionized plasmas, in partiular in the appliation to the lower solar atmosphere. The reason for fousing onto the harge exhange is the following. Using standard desriptions, in the lower solar atmosphere the ion frition with neutrals turns out to be so strong that any magneti perturbation is almost instantly destroyed Vranjes & Kono 2014; Vranjes So an enormously strong magneti field would be needed to have a magneti wave propagating. The term magneti waves is used here to desribe the Alfvén wave, and any other wave whih requires magnetized partiles. Related to this, there is a problem of magnetization. As shown in Vranjes & Krsti 2013, the proton ollision frequeny an exeed 10 9 Hz. As suh, it is several orders of magnitude above the proton gyrofrequeny even if we assume strong magneti fields of the order of 0.1 T. Under suh onditions the motion of an individual partile is similar to brownian motion. Though eletrons an be magnetized and there may be a relative drift between the protons and eletrons in the diretion transverse to the magneti field RAS

2 2 J. Vranjes, M. Kono, and M. Luna This auses the Hall effet whih onsequently may affet protons as well although they remain unmagnetized, ν i > Ω i. This is in fat manifested as the Hall term in Eq. 31 later in the text. It is therefore neessary to address the issue of magneti waves in suh a weakly ionized and highly ollisional environment from a ompletely different perspetive, by applying the harge exhange effets in a self-onsistent manner. In the present ase this implies that partiles hange their identity so frequently that it beomes impossible to distinguish what is a harged partile and what is a neutral. In pratial terms, they have a dual nature all the time and, in average, they all take part in eletromagneti perturbations. This turns out to be ruial, the frition effetively vanishes and eletromagneti perturbations an propagate. But there is a prie for this: the effetive eletromagneti fore of a wave ats on all partiles ions and neutrals at any time, so the total mean fore felt by an individual partile an be drastially redued. All these effets are desribed in the text below. 2 ESTABLISHING FACTS The ondutivity in plasmas is usually dominated by eletrons [σ = e 2 n e/m eν e, where the notation is obvious, and ν e is the eletron total ollision frequeny] beause of their muh higher mobility as ompared to ions. But the ollisions in the solar photosphere are enormously frequent, and the eletron ollision frequeny may be over Hz Vranjes & Krsti So the ondutivity is very low, and it is in the range 1 30 S/m. Compare this with the ondutivity of the terrestrial sea water, whih is around 5 S/m, i.e., similar to the photosphere. We know that ordinary fluid dynamis is perfetly able to desribe phenomena in terrestrial oeans. So the question is why suh a fluid theory is not good enough in the ase of the solar photosphere? We shall ome to this question again in the text below. In the ase of ollisions of ions with their parental atoms, the quantum-mehanial effet of indistinguishability applies. A sketh of this effet is given in Fig. 1. When an ion and an atom ollide, and we detet the partiles after ollisions, there is no way to know whih is whih unless they are somehow additionally labeled. In other words, we may have the usual elasti sattering as well as the harge exhange. This is the meaning of the indistinguishability effet aused by the harge exhange, whih therefore must be onsistently inluded in the theory Krsti & Shultz 1998, 1999a,b. This implies that the usual ollision ross setions for elasti sattering, momentum transfer, and visosity must inlude the harge exhange effet. Suh ross setions are provided for the solar hydrogen plasma in our reent papers Vranjes & Krsti 2013; Vranjes & Kono Consequene 1: The important onlusion whih follows from this is that there will be only one frition fore term in the ion momentum equation, instead of two or even more as an frequently be seen in the plasma physis literature. The same holds of ourse for the momentum equation of neutrals. This all holds for relatively low energies, roughly speaking below 1 ev. Note that this energy range in fat overs most ases of interest, beause at higher energies the H + H + Figure 1. Whih is whih? Collision of a hydrogen atom and a proton in the presene of harge exhange, demonstrating the impossibility to distinguish partiles exiting the ollision zone. amount of neutrals may be negligible and frition plays no important role. Let us now see more quantitatively how important the hange exhange is. We shall first disuss the ross setion for this proess. A lot of details are available in papers by Krsti and his group Krsti & Shultz 1999a,b, and those have been used in our papers as well Vranjes & Krsti 2013; Vranjes & Kono In Fig. 2 we give the momentum transfer ross setion affeted by the harge transfer line 1, and the harge transfer alone line 2. The total orret value for the momentum transfer is inreased due to the harge exhange roughly by a fator 2. Note that the harge transfer affets the other ross setions as well, i.e., those for elasti sattering and visosity Vranjes & Krsti The variations of the lines in Fig. 2 at low energy is the usual quantum mehanial effet. In the ase of the solar photosphere and hromosphere, in many situations it may be appropriate to use the approximate value for the momentum transfer σ mt 200 a.u., where a.u. is the square of the Bohr radius m 2. Using a value whih varies with the altitude i.e., with the temperature an be done but for some problems [like the Alfvén wave propagation studied by Vranjes 2014] this will not have a profound effet on the results. Consequene 2: Using the values for the momentum transfer presented by line 1 in Fig. 2, will most onsistently inlude the effets of harge transfer in the frition fore term in the momentum equation. Having so aurate ross setions we an now alulate the ollision frequeny for harge exhange alone. The reason for this will beome lear soon. So we use the usual formula ν i,ex = σ exn nv T i. Here, n n is the number density of neutral atoms, v 2 T i = κt i/m i, and for σ ex we use data from the line 2 in Fig. 2. For plasma parameters in the lower solar atmosphere we use the well known data for the quiet solar atmosphere Fontenla et al The result is presented in Fig. 3. Note that we are in a very narrow low-energy range with those quantum flutuations, and this is why the full line has those small bumps on it. The graph shows that, for an average proton, the harge exhange ollision frequeny hanges from Hz at the altitude h = 0 km, to around 300 Hz at h = 2000 km. For omparison we give the proton gyro frequeny Ω i = eb/m i for B = B 0 exp[ h/250], and for the two starting values of H H 2002 RAS, MNRAS 000, 1 10

3 Charge exhange in fluid desription of partially ionized plasmas 3 ross setion p-h [a.u.] proton energy [ev] lab. frame mom. transfer + h. exh. 2 - harge exhange alone proton energy [ev] CM frame Figure 2. Cross setions for harge transfer line 2, and for the momentum transfer with the effets of harge transfer inluded line 1. Atomi unit is the square of the Bohr radius, a.u= r 2 b = m 2. Data adapted from Krsti & Shultz 1999a; Vranjes & Krsti h.ex. ollision frequeny [10 6 Hz] E-3 harge exhange gyro frequeny B 0 =0.1 T gyro frequeny B 0 =0.01 T altitude [km] Figure 3. Proton ollision frequeny for harge exhange ν i,ex in the lower solar atmosphere in terms of the altitude. Straight dot-line and dash-line are gyro frequeny for B = 0.1, 0.01 T, respetively. the field B 0 = 0.1, 0.01 T. They are presented by the dashed and dotted lines. It is seen that protons are unmagnetized in a large part of the lower solar atmosphere. We stress that the ollision frequeny for total elasti sattering whih neessarily inludes the harge transfer as well is in fat higher by about a fator 4. So the magnetization is in fat even weaker than in Fig. 3, and it extends to even higher altitudes Vranjes & Krsti These numbers have a profound physial meaning. They suggest that a proton in the photosphere hanges its identity more than hundred million times per seond! The nature of the proess implies that Hydrogen atoms also hange identity and this happens with the frequeny ν n,ex = ν i,exn i/n n. 1 2 Table 1. Charge exhange ollision frequeny for neutral hydrogen at several altitudes in the solar atmosphere. h [km] ν ni,ex [Hz] For neutrals, the harge-exhange frequeny Vranjes & Krsti 2013 at several altitude is presented in Table 1. It is obvious that for both protons and Hydrogen atoms the typial harge exhange time τ ex is usually many orders of magnitude shorter than any realisti period of Alfvén waves whih an be expeted in suh an environment, τ ex τ A. 1 In view of the results presented in Fig. 3 and in Table 1, we ome to the following: Consequene 3: Both ions and atoms in the lower solar atmosphere hange identities so frequently that it is impossible to say whih is whih. For any pratial purposes, i.e., for time sales far exeeding this identity-hange time sale 1/ν ex the largest being 1/8 s for neutrals at the temperature minimum, it is ompletely obvious that, in average, all partiles will respond to eletromagneti perturbations in nearly the same manner. The nature of all partiles is intermittent; they spend one part of the time in harged state and one part in neutral state. Consequene 4: This further implies that for proesses with typial time sales far exeeding 1/ν ex, the effetive ion plasma density must be n ef =. This fat must be taken into aount, for example, in the estimates of the Alfvén wave flux from the lower solar atmosphere. These issues will be disussed in more details farther in the text. 3 EFFECTIVE CHARGE ON PARTICLES To proeed, we have to find the effetive or mean harge on partiles, whih is diretly aused by the harge exhange effet. This an be done as follows. For eah individual partile in the system there is a harateristi time interval τ = τ 1 + τ 2 within whih it passes through two yles: i when it is harged τ 1 interval, and, ii when it is in the neutral state τ 2 interval. After these two yles the proess is repeated, so we an fous on the τ interval only. Beause of the harge flip-flop, an effetive mean harge of suh a partile is q ef = 1 τ τ 0 qft, τ = 1 ν in + 1 ν ni, where we have a step-funtion ft = 1 in the starting interval τ 1 0, 1/ν in, and ft = 0 in the remaining interval of the yle τ 2 1/ν in, 1/ν ni. Note, the ollision frequenies here desribe the harge-exhange. This yields an effetive partile harge in the interval τ ex: q ef = q q = ν in/ν ni 1 + n n/n i Without neutrals i.e., in the upper layers we have q ef = q, as it should be. But for a weak ionization the effetive harge 2002 RAS, MNRAS 000, 1 10

4 4 J. Vranjes, M. Kono, and M. Luna q ef /q E-3 1E-4 1E-5 1E x [km] Figure 4. Normalized effetive harge on a partile in a proton- Hydrogen mixture, in terms of the altitude. is very small. This is due to the fat that, in average, a partile spends muh more time in the neutral state. So in the photosphere q ef /q is a very small parameter whih is spatially dependent, and it hanges rapidly with the altitude. At the temperature minimum h = 490 km its value is q ef /q 10 6, and in suh an environment slow eletromagneti perturbations are not expeted to onsiderably affet the partile dynamis. On the other hand, at the altitudes h = 1580, 2018 km, the values of the parameter are, respetively, q ef /q = 0.06, This normalized effetive harge is presented in terms of the altitude in Fig. 4. Consequene 5: The effetive harge applies to eah and every ion and atom in the system. Therefore for any proess with the harateristi time τ τ all partiles are subjet to eletromagneti fores in the same manner. This implies that the usual separation of the population into ions reating on eletromagneti perturbation and neutrals whih initially remain stati in the presene of suh perturbations and later move due to frition by ions does not hold any longer. The total eletromagneti fore ating on all partiles in a unit volume is proportional to q ef n i,ef. The mean motion of all partiles is thus the same and in pratial sense there is no frition in the system. The single fluid of ions and neutrals as a whole an now be weakly affeted by the magneti field as frition between them is absorbed into the one fluid. So, as long as the effets of the harge exhange are not inluded in a self-onsistent manner, every analytial model will fail in desribing eletromagneti perturbations in suh an environment. 4 THE ISSUE OF MAGNETIZATION The ollisions of various kinds, both elasti and inelasti still remain as a fat, and harged partiles may be unmagnetized. Following the usual wisdom they should not be diretly affeted by the magneti field, but see the omment about the Hall effet earlier in the text, Se. 1. This is the ase in the lower solar atmosphere and it may be the same in any other partially ionized plasma. In the photosphere, the total elasti sattering ollision frequeny may far exeed the gyro-frequeny and this may even happen with eletrons, as it is demonstrated in our reent referenes Vranjes & Krsti 2013; Vranjes & Kono However, it is laimed that magneti waves do exist and that they propagate in the photosphere and hromosphere Jess et al. 2009; Kostik & Khomenko If this is a fat then we are definitely still missing a omplete piture, and the understanding of magneti waves in this weakly ionized environment remains elusive and unexplained. The purpose of this work is to get some physial insight into the physis of waves in suh an environment, and to find mehanisms whih might allow for their existene. 4.1 Rotation of partile displaement plane in the presene of an Alfvén wave In general, we may have two possible ases: with magnetized, and with unmagnetized partiles. Let us assume that the bakground magneti field is in z-diretion, and it is perturbed in y-diretion. This implies a perturbed eletri field in x-diretion in aordane with Faraday law. This geometry is presented in Fig. 5, and more details are given below in the text Unmagnetized partiles. In the ase of unmagnetized partiles the ollision frequeny ν j above the gyro-frequeny Ω j the effet of the Lorentz fore is negligible. The partiles are never able to perform a gyro-rotation, therefore their motion annot have a drift harater. In other words, the partile predominantly moves in the diretion of fore, whih is the eletri field fore in the present ase. But in view of the disussion presented in the text above, due to the harge exhange, this applies to all of them. We are dealing with a mean fore ating on the whole fluid. This motion is depited in Fig. 5 by the lower part of the graph Magnetized partiles When partiles are magnetized, ν j < Ω j, their dominant motion is due to the E B-drift Chen 1988, and this drift motion is the same for both eletrons and ions. The motion in the diretion of fore is only due to polarization drift Vranjes & Kono Compared with the leading drift motion, the polarization drift is very small and the ratio of the two displaements is of the order of ω/ω j 1. So when an Alfvén wave propagates from below through the upper solar layers, there must be a 90 degrees hange of the displaement plane of the fluid as it passes through regions with unmagnetized and magnetized partiles. This is depited in Fig. 5, ompare the lower and the upper graphs desribing the partile speed. This rotation of the displaement plane happens only one, when partiles hange their feature from being un-magnetized to being magnetized. Observers should be able to detet suh a rotation. This motion of magnetized partiles an be desribed by starting from the momentum equation of the ion fluid, after keeping only terms essential for this disussion. Note that we are now in the frequeny range when the ollisions are relatively weak and we shall assume that the harge exhange plays no big role, so we have the usual frition fore in the ion equation: m i d v i dt = e E + vi B m iν ie v i v e 2002 RAS, MNRAS 000, 1 10

5 Charge exhange in fluid desription of partially ionized plasmas 5 a magnetized x v B 0 unmagnetized E a z v k b B 1 b y E=E 1 osk z - t Figure 5. Twist in the motion of the fluid element in the presene of an Alfvén wave whih oupies layers of unmagnetized plasma the lower region of the piture and magnetized one the upper part of the piture. The motion within a half period of the wave is presented. m iν in v i v n. 3 Ion-eletron frition is omitted here only for simpliity, the mass differene is suh that frition with eletrons is of a seondary importane. So to explain the essential physis whih is in fous here, we keep only the leading order effets. The omplete frition fore has indeed been studied in our reent works Vranjes & Kono 2014; Vranjes Performing just the vetor transformations, from Eq. 3 the linearized perpendiular speed of ions beomes: v i = α i 1 e z E B e z v i + νin E1 0 Ω i t Ω i B 0 νin v i + iν2 in ν ni v Ω 2 i t Ω 2 i iνinνni e z v i. 4 i ω n Ω iω n Here, ω n = ω + iν ni, and α i = 1/1 + ν 2 in/ω 2 i 1. We have used the momentum equation for neutrals assuming that they are set into motion by frition with ions only Vranjes & Kono 2014, whih gave: v n = iνni ω + iν ni v i. 5 For small ratio ν 2 in/ω 2 i magnetized partiles, in average ollisions will not signifiantly affet the partiles. Their motion will have drift harater, and the dominant term in Eq. 4 is the first term the E B-drift. Partiles move predominantly perpendiular to the eletri fore, therefore their dynamis is mainly in the y- diretion. This leading term may be replaed into the remaining terms on the right-hand side of the reurrent formula 4, and they beome determined through the eletri fore. Note, the seond term is the polarization drift, whih is in the diretion of the eletri field, but this is small in view of the small parameter /t/ω i ω /Ω i. For negligible ollisions we would have only these two drifts: v i = v E e y + v p e x. The motion is essentially two-dimensional although the orbits are strongly extended ellipses. When partiles are unmagnetized, writing the perpendiular speed in the shape 4 will make no sense beause Eq. 4 is a reurrent formula only when v E e y is the dominant drift. The motion annot have drift harater in the ase of unmagnetized partiles, so in this ase we are supposed to disuss the starting equation 3. First, we note that the only ative fore in Eq. 3 is due to eletri field. For unmagnetized partiles in the presene of harge exhange, frition terms will vanish in any ase, as we laimed earlier. The remaining Lorentz fore only ats on partile or fluid element that is already in motion. In other words, the momentum equation ontaining only the Lorentz fore on the right-hand side would have a trivial solution v i = 0. Now, when partiles are not magnetized their dominant motion will be the diret one, in the diretion of the eletri field vetor. Consequently, only after a fluid element or a partile start moving in the x-diretion, it will start feeling the effets of the Lorentz fore term. For a osωt-type dependene of the eletri field E, the maximum displaement in the y-diretion is expeted when the speed is maximum, that is when E is maximum. After that the partile displaement in the y-diretion beomes smaller and eventually it moves only in the diretion of the field E after half-period of the wave. After this, the eletri field diretion is hanged and the partile or fluid element is moving bak along the same trajetory. So the motion is again two dimensional as in the usual multiomponent desription, but the dominant motion is in the diretion of fore, and not perpendiular to it as in the ase of magnetized speies. To the best of our knowledge, the twist presented in Fig. 5 has never been disussed or suggested in the literature. This is a new result whih may be appliable even in ompletely different situations, i.e., in fully ionized plasmas with frequent ollisions. The usual wisdom about Alfvén wave propagation Chen 1988 implies a dominant drift motion of partiles, and as suh it annot explain their presene in an environment with unmagnetized speies. 5 MOMENTUM EQUATION IN THE PRESENCE OF CHARGE EXCHANGE In every moment the harge quasineutrality is satisfied, so in the present ase this ondition reads en e en i = q ef n ef, where the effetive number of harged partiles is n ef =. For low frequeny proesses, with harateristi time sales τ τ, the orret momentum equation in the presene of harge-exhange and for a plasma with ions and their parental atoms reads: m in d v dt = q ef n ef E + v B qni E + v B RAS, MNRAS 000, 1 10

6 6 J. Vranjes, M. Kono, and M. Luna Here, n = n ef and v denotes the mean speed of both ions and their parental atoms. We have omitted the frition with eletrons only for simpliity beause it is not relevant for the harge exhange disussion. The number of eletrons is and their distribution funtion are not hanged in the proess. But this term should be kept in a proper multi-omponent desription. The equation shows that in an ions-atoms mixture in any moment the total EM fore per unit volume applies only on partiles that are in that moment harged. Obviously, in the limit of vanishing neutrals fully ionized plasma, Eq. 6 redues to the usual ion momentum equation. Note that this equation an be derived from first priniples, starting from kineti equations for a single heavy omponent ions plus neutrals and assuming that all of them have the effetive harge q ef, and that the number density n ef = n. Two alternative ways of its formal derivation are presented in Ses. 6, 7 below. 6 FORMAL DERIVATION OF MOMENTUM EQUATION In a plasma where harge exhange ours frequently, ions and neutrals are replaing eah other in a very short time interval, so the equations of motion for ions and neutrals may be unified for the jth partile as: d x j dt = vj, 7 d v j dt = θt e m [ E xj, t + v j B x j, t ]. 8 Here j = 1, 2,, N i, N i + 1,, N i + N n, and θt is a random funtion defined by { 1 for t Si θt =, 9 0 for t S n, and S i and S n are the harged state and neutral state, respetively. The masses of ions and neutrals are assumed to be equal, whih is appropriate for a proton-hydrogen plasma. Note that the number densities of ions and neutrals are preserved and the total harge density is as well. The eletromagneti fields are determined self-onsistently by the Maxwell equations. A mirosopi density distribution funtion for the unified ion-neutral partiles is introdued by F x, v, t = 1 V N i +N n j=1 δ x x jtδ v v jt, 10 whih allows us to express the urrent and density. We an further write t F x, v, t = 1 V = 1 V = N i +N n i=1 N i +N n i=1 d xi dt [ v x θt e m δ x xitδ v vit t x d vi dt δ x x itδ v v it v ] E + v B F α x, v, t. v So this yields F α t + v Fα x + θt e m E + v B F α v = Sine the density distribution funtion F is not ontinuous and depends on the initial onditions of partiles { x j0, v j0}, we introdue an ensemble averaged distribution funtion over the most probable distribution of the initial onditions as f x, v, t =< F x, v, t >, 12 where < F > denotes an ensemble average over the initial onditions. Then the partile disreteness is expressed by δf x, v, t = F x, v, t < F x, v, t >. 13 Obviously < δf x, v, t >= 0. Sine E and B depend on eletrons as well, we an also introdue marosopi and mirosopi fields as E =< E > +δ E, B =< B > +δ B. To simplify notation, in the following we use E and B instead of < E > and < B >. Taking the ensemble average we obtain: f f + v t x + θt e m E + v B f f v =, 14 t where f = < θt e δe + v δ B δf >. 15 t m v The right hand side in Eq. 14 represents the interation among disrete partiles through the eletromagneti fields, that is, ollisional effets due to the disreteness. When the right hand side is omitted, the remaining equation is a generalized Vlasov equation for plasmas with frequent harge exhange phenomena. Introduing a mean density and a mean veloity as n = fd v =, n v = vfd v, and taking the first and seond moments of the kineti equation without the ollision terms for simpliity gives n + n v = 0, t 16 ne n v + n v : v = θt E + v B. t m 17 When the two obtained equations are ombined we obtain: v t + v v = θt e E + v B. 18 m We may take an average over the time longer than the harateristi time of harge exhange, and in this ase n i θt =. Eq. 18 is the same as Eq. 6. This is a ompletely new result appliable for all proesses satisfying the ondition 1. Yet another way of deriving the ommon momentum equation for ions and their parental atoms is presented in the next setion RAS, MNRAS 000, 1 10

7 7 AN ALTERNATIVE DERIVATION OF THE COMMON MOMENTUM EQUATION 6 We shall show that the momentum equation 6, whih desribes the ommon dynamis of all heavy partiles ions and neutrals in the presene of frequent harge exhange, an be formally derived even by starting from the two separate momentum equations for ions and neutrals. f i ei [ + v fi+ E x, t + v f i B x, t ] t m i v = fi t f n t + v fn = fn t Charge exhange in fluid desription of partially ionized plasmas 7, 19, 20 where f α x, v, t is an ensemble averaged distribition funtion of α speies, f α x, v, t = 1 N α V < δ x x jtδ v v jt >, 21 j=1 and f α/t is a ollision integral. In general eletrons are involved in the ollision integral whih may be negleted for the present ase and for harge exhange the density and momentum are onserved, that is, fα t d v = 0, v fα t d v = 0. Here we introdue a two-body distribution funtion by f in x, x, v, v, t = f i x, v, t f n x, v, t, and we get { t + v + v + ei = where fin, t fin t [ } E x, t + v B x, t ] f in m i v fi fn =f n + f i, 22 t t f in x, x, v, v, t d xd x d vd v = n in n. 23 For a two-speies system with the same number of partiles the two-body distribution funtion is expressed due to the one-to-one orrespondene by f x, x, v, v, t = nf x, v, t < δ x x δ v v > + C x, x, v, v, t, 24 where the seond term of RHS is a orrelation funtion. However, in the ase studied here, in general the numbers of ions and neutrals are different and Eq. 24 is not straightforwardly appliable to the present ase. So now we introdue a enter of density oordinate system through R = ni x + nn x, U = n i v + n n v, r = 2 ni x nn x, u = 2 ni v nn v, This transforms Eqs. 22, 24 into { t + U R + u r + ei m i [ n i ni + n n E r R + 2n i RAS, MNRAS 000, 1 10, t ] ni + nn + u U + B 2n i 2 ni + n n R r +, t 2n i 2 } U + 2 u f in = fin, t [ f ind xd vd x d v nin n = F in R, U, t ni nn ni < δ r 2 nn R δ u 2 U > +C in R, U, r, u, t ] d Rd Ud rd u, 25 where we have used x x ni + nn2 ni nn = r 2 R. 4n in n The volume element of the phase spae is preserved as f in x, x, v, v, td x x d vd v =F α R, U, r, u, td Rd Ud rd u. The oeffiient of F in R, U, t in the RHS of Eq. 25 is uniquely determined by the ondition that it is symmetri with respet to the speies and proportional to N. From Eqs. 23, 25 we have F ind Rd U =. 26 When harge exhange is so frequent that ions and neutrals are supposed to behave together, the seond term of the RHS of Eq. 25 is negleted ompared with the first term. Integrating Eq. 25 over r and u yields [ t + U R + F ] m i U F in R, U, t = 0, 27 [ ] n i F = e i E R, n i + n t + U B R, t, n where we have used fin d ud r = t Eq. 27 is a generalized Vlasov kineti equation for the speies with an effetive harge q ef whih we introdued earlier in the text, see Eq. 2. So now, following the wellknown proedure, it diretly yields all fluid equations as zero-momentum ontinuity equation, the first momentum equation of motion 6, et. 8 APPLICATION TO ALFVÉN WAVES In Fig. 5, a plain-polarized wave is assumed just to explain the basi features, whih is fine most of the time for large wavelengths kλ i 1, λ i = /ω pi, ω 2 pi = e 2 n i0/ε 0m i. But in general any diretion of the perturbed perpendiular field is possible i.e., we should allow for a irular polarization. So we may have B 1 = B x e x + B y e y and similar for the eletri field. Using momentum equation for inertialess eletrons we have E = v e B. 29

8 8 J. Vranjes, M. Kono, and M. Luna In the present ase, for the time sales determined by 1, the urrent is produed by all heavy speies and by eletrons, j = q ef n ef v en e v e, whih with the help of Ampère law yields v e = v 1 µ 0en i B. 30 This is used in Eq. 29 and the resulting expression is further used in the Faraday law whih finally yields B t = v B 1 µ 0en i [ B B ]. 31 We now linearize our new momentum equation 6, and we eliminate the eletri field with the help of the Faraday law E x = ωb y/k, E y = ωb x/k, v = ieni m iωn [ ωby k + v yb 0 e x ωbx k v xb 0 e y ]. 32 We further linearize Eq. 31, and we have a losed set for B 1, v 1. This yields the following dispersion equation: [ ] ω 4 2ω 2 k 2 va nn0 k 2 λ 2 i + k 4 va 4 = n i0 2 Here, v 2 a = B 2 0/[µ 0m in i0 + n n0] 34 is the Alfvén speed whih inludes both ions and neutrals. This an be written as: { ω 2 = k 2 va nn0 k 2 λ 2 i ± 1 + nn0 1/2 kλi n i0 2 n i0 [ ] } nn0 k 2 λ 2 1/2 i. 35 n i0 4 Setting n n0 = 0 yields a solution well-known from the literature Almaguer 1992, with orretions ontaining the small parameter kλ i whih gives so alled left and right irularly polarized waves. In the usual AW limit kλ i 1 this yields ω 2 = k 2 v 2 a. 36 In any partiular moment, all partiles both harged and unharged take part in low-frequeny perturbations, and this is refleted in the expression for the AW speed. 9 PROVING THE THEORY I. There are numerous reports about the presene of Alfvén waves in the lower solar atmosphere. If this is taken as a fat then this might be understood as a plausibility of the theory presented here beause the frition, when alulated in the usual way, destroys them instantly. II. The presene of the waves an only be explained with the help of harge exhange, as desribed with our Fig. 5, whih is important for the theory presented here. If suh a rotation of the Alfvén wave plane ould be verified by observations, this alone would mean the verifiation of the theory as well. To help observers in deteting the twist presented in Fig. 5, here we give some more details about possible loations altitudes where suh a transformation an take plae. So the wave features depend on the ratio ν i/ω i. 37 Here, ν i is the most dominant ollision frequeny for ions. As shown in Vranjes & Krsti 2013, this an be either ν in or ν ii, and this is dependent on the altitude. Now, important to stress is that in both ases the ollision frequeny ontains the ross setion for total elasti sattering. For ion-neutral ollisions this ross setion is determined self-onsistently and it ontains effets of the harge exhange. From Fig. 1 in Vranjes & Krsti 2013 it may be seen that elasti sattering ross setion is roughly by a fator 4 greater than the harge exhange ross setion whih we used in Fig. 3 in the present work. So using very aurate values from Vranjes & Krsti 2013, and for the magneti field of the same shape as used previously in the text, B = B 0 exp[ h/250], we an now hek at whih altitudes the ratio 37 is equal to unity. The wave is expeted to hange its polarization plane around this altitude. From the graphs in Vranjes & Krsti 2013 it may be seen that even assuming an exeptionally strong magneti field B 0 = 0.1 T, we have ν in > Ω i up to roughly h = 1400 km. But from this altitude and above, Coulomb ollisions beome dominant so ν ii > Ω i and protons remain unmagnetized. A big drop in the Coulomb ollision frequeny takes plae only somewhere above h = 2400 km. So only at these altitudes, or higher, the mentioned twist of the wave plane should be expeted. This all holds on ondition that the assumed variation of the magneti field is valid. But even if we keep the magneti field B 0 = 0.1 T onstant with altitude, the protons will remain unmagnetized for some 800 kilometers. Therefore the hange of the wave plane annot possibly happen below this altitude. We believe these should be lear indiators for observers; the wave twist should be taking plae in the range km in the lower atmosphere. With some more aurate values for the magneti field the altitude an be pinpointed far more preisely. III. The estimated Alfvén wave speed whih an be found in some referenes Jess et al is around 22 km/s, with the magneti field whih they assumed to be 0.1 T. Suh a speed an only be obtained if the number density involved [see Eq. 34 in the previous setion] has the value n = m 3. We note that the ion number density in the photosphere is several orders of magnitude below this value Fontenla et al. 1993, so the density whih enters the expression for the Alfvén speed is in fat the density of neutrals n n0, as expeted from our analysis. This looks like a rather onvining proof of our theory. However, it is neessary to stress that suh an expression for the Alfvén speed with the neutral density has been used in the past in the usual MHD desriptions as well. In fat, it an even be obtained analytially under ertain onditions Kulsrud & Peare In suh a modeling there are two propagation windows for Alfvén waves, determined either by the amount of ollisions or equivalently by the wave length. In other words, there is a range of parameters for whih the wave annot propagate and this has been known for more than half a entury Kulsrud & Peare We have disussed these issues in detail in our reent work Vranjes & Kono But MHD theory annot provide a self-onsistent explanation for the issue of magne 2002 RAS, MNRAS 000, 1 10

9 Charge exhange in fluid desription of partially ionized plasmas 9 tization of partiles and the presene of waves; it operates with guiding enter of mass and not with atual partiles. Even drift motion annot be desribed within suh a rough model. The magnetized partiles are an underlying assumption in MHD theory; in other words, plasma is by definition affeted by the magneti field. Contrary to this, our theory inludes the harge exhange, it is based on the single ondition 1, and the theory works regardless if partiles are magnetized or not, hene the predited twist of the wave plane in Fig. 5. Note also that there are no two propagation windows in our theory, the wave propagates for any parameters. IV. Yet another pratial way of proving our theory may be dedued by analyzing displaement speed of the plasma involved in wave motion. As disussed in Se. 4, the nature of displaement depends on magnetization. In the upper layers with magnetized partiles the leading order speed is due the E B-drift v E = E1 = ω B B 0 k B 0 In this expression it might be appropriate to set ω/k = a, where a = B 0/µ 0m in i0. This beause the harged speies are subjet to drift motion, they are magnetized, whih implies that the number of neutrals is onsiderably redued so the Alfvén speed ontains ion density only. In the lower layers our new momentum equation 6 applies. The dominant motion must be due to the eletri field fore. The speed of this diret motion from Eq. 6 is: v dir = eni Ωi n i0 ω B 1 E1 =. 39 m iωn ω n i0 + n n0 k B 0 In this ase we know that ω/k = v a, where v a is given by Eq. 34. Assume now that the wave oupies the whole region, with magnetized partiles in the upper layers and with the environment where the harge exhange plays the main role in the lower layers, and that the frequeny is the same in both regions. If we further assume that the magneti field perturbation is the same whih may be a very rough assumption, this yields: v dir v E = Ωi ω n i0 n i0 + n n0 3/2. 40 Obviously, both terms on the right-hand side are altitude dependent. The minimum value of the seond term is 10 9 but it inreases with altitude. The first term an be huge, so the ratio of the two speeds an have any value. For example, taking parameters from Jess et al. 2009, the speed ratio at the temperature minimum may beome of the order on unity. But the estimate formula 40 is rough due to various reasons; one obvious is that the ratio B 1/B 0 annot be the same in both regions. The equilibrium magneti field hanges with the altitude and the value for B 1 is hard to guess. But it is definitely possible to find differenes in the displaement speeds for the two regions by using the expressions 38, 39 as they are, without further assumptions, and then alulate the speed by using the parameters from observations. V. The Alfvén speed introdued in Eq. 34 inludes the neutral density, whih hanges for about 6 orders of magnitude. Being under the square root, this implies that the speed hanges by a fator 1000 if the magneti field is onstant. Though the magneti field hanges as well, but there exist models for its variation, so the resulting hanges in the speed ould be measured. 10 SUMMARY The analysis presented here leads to the following onlusions: a For low-frequeny proesses, partiles ions and neutrals in plasma where harge exhange is taking plae should be treated as a single fluid with an effetive harge as alulated in the text. b In a plasma like the solar photosphere, the lowfrequeny magneti waves an propagate even if partiles are not magnetized. This is due to the fat that it is not possible to separate partiles into two different populations as harged and neutral speies, therefore there an be no frition fore between ions and neutrals in the usual sense of its meaning. Beause of the harge exhange, all partiles take part in low-frequeny eletromagneti perturbations. This is the reason why eletrodynami theory is needed in an environment like the photosphere, and not a pure fluid theory as in the terrestrial seas, in spite of the fat that the ondutivity in these two environments is pratially the same. The frequent harge exhange indues one fluid of ions and neutrals and as suh, ollisions between ions and neutrals are absorbed in suh a ollisionless one fluid. So the harateristi frequeny whih needs to be ompared with the ylotron frequeny is no longer the ollision frequeny but the fluid osillation frequeny. In this sense one fluid is magnetized or magneti although the harge density is redued. Eq. 6 is no longer an equation for single partile. d The effets desribed here are prinipally different from the lassi MHD interpretation whih predits two propagation windows for the Alfvén waves in a weakly ionized environment. Suh a behavior within MHD desription is aused by the frition between ions and neutrals. The harge exhange effets disussed here have a ompletely different nature, and the analysis presented in the work suggests, in average, the absene of ion-neutral frition. We expet that the flux of Alfvén waves eventually oming from the photosphere ould be onsiderably modified, and their role in the heating of upper layers should be reonsidered. e The theory presented here explains why magneti waves an still propagate in an environment in whih ions are obviously not magnetized. In appliation to the lower solar atmosphere, the results presented here predit that the plane of displaement of the fluid must hange by 90 degrees when an Alfvén wave propagates from the area where partiles are un-magnetized photosphere to the area where they are magnetized hromosphere. Providing some aurate values for the magneti field, in view of the most aurate ross setions we are using here, and the orresponding ollision frequenies, we are able to very aurately determine altitudes at whih partile magnetization hanges. In the present study we have been foused on harge exhange as the most dominant of all inelasti ollisions, whih is in the same time a rather speifi inelasti effet beause 2002 RAS, MNRAS 000, 1 10

10 10 J. Vranjes, M. Kono, and M. Luna it preserves the number density of all partiles involved in it. In appliation to the lower solar atmosphere, other inelasti ollisions, like radiative reombination and 3-body reombination, are known to be far less frequent Vranjes & Poedts 2006, with maximum frequenies of the order of a few tens of Hz only. ACKNOWLEDGMENTS J.V. aknowledges a the finanial support from the Spanish Ministry of Eonomy and Competitiveness MINECO under the 2011 Severo Ohoa Program MINECO SEV , and b fruitful disussions with M. Collados and E. Khomenko. M. Luna aknowledges the support by the Spanish Ministry of Eonomy and Competitiveness through projets AYA , AYA and AYA P. This work ontributes to the deliverables identified in FP7 European Researh Counil grant agreement , Magneti Connetivity through the Solar Partially Ionized Atmosphere PI: E. Khomenko. REFERENCES Almaguer J. A., 1992, Phys. Fluids B, 4, 3443 Brandsen B. H., MDowell M. R. C., 1992, Charge Exhange And The Theory of Ion-Atom Collisions. Clarendon Press, Oxford Chen F. F., 1988, Introdution to Plasma Physis and Controlled Fusion. Plenum Press, New York Dalgarno A., 1958, Phil. Trans. R. So. Lond. A, 250, 426 Eihler J., 2005, Letures On Ion-Atom Collisions. Elsevier, Amsterdam Fontenla J. M., Avrett E. H., Loeser R., 1993, ApJ, 406, 319 Glassgold A. E., Krsti P. S., Shultz D. R., 2005, ApJ, 621, 808 Haerendel G., 1992, Nature, 360, 241 Jess D. B., Mathioudakis M., Erdélyi R., Crokett P. J., Keenan F. P., Christian D. J., 2009, Siene, 323, 1582 Kostik R., Khomenko E., 2013, A&A, 559, A107 Krsti P. S., Shultz D. R., 1998, Atomi and Plasma- Material Interation Data for Fusion Vol. 8. IAEA, Vienna Krsti P. S., Shultz D. R., 1999a, J. Phys. B: At. Mol. Opt. Phys., 32, 3485 Krsti P. S., Shultz D. R., 1999b, Phys. Rev. A, 60, 2118 Kulsrud R., Peare W. P., 1969, ApJ., 156, 445 Raizer Y. P., 1991, Gas Disharge Physis. Springer-Verlag, Berlin Heidelberg Vranjes J., Poedts S., 2006, Phys. Lett. A, 348, 346 Vranjes J., Krsti P. S., 2013, A&A, 554, A22 Vranjes J., Kono M., 2014, Phys. Plasmas, 21, Vranjes J., 2014, MNRAS, 445, 1614 This paper has been typeset from a TEX/ L A TEX file prepared by the author RAS, MNRAS 000, 1 10