Issue 
A&A
Volume 534, October 2011



Article Number  A93  
Number of page(s)  7  
Section  The Sun  
DOI  https://doi.org/10.1051/00046361/201117380  
Published online  11 October 2011 
Damping of Alfvén waves in solar partially ionized plasmas: effect of neutral helium in multifluid approach
^{1}
Space Research Institute, Austrian Academy of Sciences, Schmiedlstrasse 6, 8042 Graz, Austria
email: teimuraz.zaqarashvili@oeaw.ac.at; maxim.khodachenko@oeaw.ac.at; rucker@oeaw.ac.at
^{2} Abastumani Astrophysical Observatory at Ilia State University, University St. 2, Tbilisi, Georgia
Received: 31 May 2011
Accepted: 6 September 2011
Context. Chromospheric and prominence plasmas contain neutral atoms, which may change the plasma dynamics through collision with ions. Most of the atoms are neutral hydrogen, but a significant amount of neutral helium may also be present in the plasma with a particular temperature. Damping of MHD waves due to ion collision with neutral hydrogen has been studied well, but the effects of neutral helium are largely unknown.
Aims. We aim to study the effect of neutral helium in the damping of Alfvén waves in solar, partially ionized plasmas.
Methods. We consider a threefluid magnetohydrodynamic (MHD) approximation, where one component is electronprotonsingly ionized helium and the other two components are the neutral hydrogen and neutral helium atoms. We derive the dispersion relation of linear Alfvén waves in isothermal and homogeneous plasma. Then we solve the dispersion relation and derive the damping rates of Alfvén waves for different plasma parameters.
Results. The presence of neutral helium significantly enhances the damping of Alfvén waves compared to the damping due to neutral hydrogen at certain values of plasma temperature (10 000−40 000 K) and ionization. Damping rates have a peak near the ionneutral collision frequency, but decrease for the higher part of the wave spectrum.
Conclusions. The collision of ions with neutral helium atoms can be important for the damping of Alfvén waves in chromospheric spicules and in prominencecorona transition regions.
Key words: Sun: atmosphere / Sun: oscillations
© ESO, 2011
1. Introduction
The presence of neutral atoms may change the dynamics of astrophysical plasma due to their collisions with charged particles. Ionneutral collision may lead to damping of magnetohydrodynamic (MHD) waves in the solar photosphere/chromosphere and prominences (Khodachenko et al. 2004; Forteza et al. 2007; Zaqarashvili et al. 2011). Only part of the hydrogen atoms are ionized in the solar photosphere, but the ionization degree increases with height owing to the rise of plasma temperature. The plasma in the transition region/lower corona is almost fully ionized, but upper chromosphere, spicules, and prominences contain significant amounts of neutral atoms.
Numerous papers have studied the effects of plasma partial ionization on different aspects of the solar atmosphere (Haerendel 1992; De Pontieu & Haerendel 1998; De Pontieu et al. 2001; Khodachenko & Zaitsev 2002; James & Erdélyi 2002; James et al. 2004; Khodachenko et al. 2004; Khodachenko & Rucker 2005; Leake et al. 2005; Leake & Arber 2006; Arber et al. 2007; Forteza et al. 2007; 2008; Vranjes et al. 2008; Soler et al. 2009a,b; Gogoberidze et al. 2009; Carbonell et al. 2010; Singh & Krishan 2010; Goodman & Kazeminezhad 2010; Tsap et al. 2011). Most of these papers consider only the effects of neutral hydrogen and use a singlefluid MHD approximation. Recently, Zaqarashvili et al. (2011) have shown that the consideration of twofluid approach significantly change the dynamics of MHD waves compared to the singlefluid one when the wave frequency is comparable to or higher than the ionneutral collision frequency. They show that the damping rates of transverse waves (Alfvén, fast magnetoacoustic) due to ionneutral collision reach a maximum value at certain wavelengths, but then begin to decrease for higher frequencies. The twofluid description sheds light on a few disputed problems, such as the damping of slow magnetoacoustic waves (Forteza et al. 2007). Therefore, the multifluid approach reveals more complete dynamics for partially ionized plasmas.
On the other hand, hydrogen is not the only sort of neutral atom that may be important in the damping of MHD waves. The solar plasma may contain significant amount of neutral helium in the certain regions of the solar atmosphere. The upper chromosphere/transition region, spicules, and prominences are the regions where the presence of neutral helium may be potentially important for the damping of MHD waves. The first step towards including neutral helium has been taken by Soler et al. (2010). They used the singlefluid MHD approach and conclude that the neutral helium has no significant influence on the damping of MHD waves in the prominence plasma. However, the effect of neutral helium can be enhanced in the regions of the solar atmosphere with higher temperatures, therefore further analysis is needed.
In this paper, we study the Alfvén waves in multifluid partially ionized plasma that contain neutral hydrogen and neutral helium atoms. Particular attention is paid to the wave damping due to ionneutral collisions and the difference between damping rates in partially ionized plasma with and without neutral helium. We derive the threefluid MHD equations from initial fivefluid equations and solve the linearized equations in the simplest case of a homogeneous plasma.
2. Main equations
We study partially ionized plasma, which consists in electrons (e), protons (H^{+}), singly ionized helium (He^{ + }), neutral hydrogen (H) and neutral helium (He). We suppose that each sort of species has Maxwell velocity distribution, so they can be described as separate fluids. Below we first write the equations in fivefluid description and then perform the resulting transition to the threefluid approach.
2.1. Multifluid equations
We aim to study the dynamics of Alfvén waves, therefore we consider the incompressible plasma. We also neglect the viscosity, the heat flux, and the heat production due to collisions between particles. Then the fluid equations for each species can be written as (Braginskii 1965; Goedbloed & Poedts 2004) where m_{a}, n_{a}, p_{a}, V_{a} are the mass, the number density, the pressure, and the velocity of particles a, E is the electric field, B the magnetic field strength, R_{a} the change of impulse of particles a due to collisions with other sort of particles, e_{a} = ± 4.8 × 10^{10} statcoul is the charge of electrons, protons, and singly ionized helium (note, that e_{a} = 0 for neutral particles), and c = 2.9979 × 10^{10} cm s^{1} is the speed of light. Plasma is supposed to be quasineutral, which means n_{e} = n_{H + } + n_{He + }. The description of the system is completed by Maxwell equations, which have the forms (without displacement current) where (5)is the current density.
In the case of a Maxwell distribution in each sort of particle, R_{a} takes the form (Braginskii 1965): (6)where α_{ab} = α_{ba} are coefficients of friction between particles a and b.
For time scales longer than ionelectron and ionion collision times, the electron and ion gases can be considered as a single fluid. This significantly simplifies the equations by taking the smallness of electron mass into account with regards to the masses of ion and neutral atoms. Then the fivefluid description can be changed by a threefluid description, where one component is the charged fluid (electron+protons+singly ionized helium) and the other two components are the gases of neutral hydrogen and neutral helium.
We use the definition of total density and the total velocity of charged fluid as (7)(8)correspondingly. Then the sum of momentum equations for electron, proton, and singly ionized helium gives the equation (9)where (10)is the collision term, p = p_{e} + p_{H + } + p_{He + } is the total pressure of charged fluid, ξ_{H + } = m_{H + }n_{H + }/ρ, ξ_{He + } = m_{He + }n_{He + }/ρ are relative concentrations of ions, and w = V_{H + } − V_{He + } is the relative velocity of protons and helium ions. In deriving Eq. (9), we neglect the electron inertia. Friction of plasma with neutral atoms is mostly defined by the ionneutral collision (), so the collision of electrons with neutrals is also neglected here.
The relative velocity of heavy ions, w, is small compared to the centerofmass velocity. This can be easily seen if one uses the momentum equations for proton and singly ionized helium without neutral atoms. Subtracting the equations one may get (11)The ratio of collision and Lorentz terms is proportional to δ_{H + He + }/Ω_{H + }, where Ω_{H + } = eB/m_{H + }c is the proton gyrofrequency (~10^{5} s^{1} for the magnetic field strength of 10 G) and δ_{H + He + } is the ionion collision frequency. This ratio is much less than unity in the solar chromosphere, so the collision term can be neglected. The ratio of the second term on the lefthand side and Lorentz term is proportional to ω/Ω_{H + }, where ω is the wave characteristic frequency. This ratio is much less than unity for the time scales that are much longer than an ion gyroperiod, therefore the inertial term with w can also be neglected. The third term on the lefthand side is even smaller. Then one may get (12)which yields w ≪ V for much longer time scales than ion gyroperiod. Therefore, the terms in Eqs. (9)–(10), which contain w, are smaller than other terms. Consideration of these terms complicates the calculations, but has no significant influence on the final results as we are interested in the effects of ionneutral collisions. The terms are thus neglected in the rest of the paper.
Momentum equations for neutral hydrogen and neutral helium atoms can be written as where the collision terms are Ohm’s law is obtained from the momentum equation of electrons, and it takes the form (17)Here the terms containing the ion relative velocity, w, and electronneutral collision are again neglected.
Maxwell equation (Eq. (3)) and Ohm’s law (Eq. (17)) lead to the induction equation (18)where the battery, Hall, and Ohmic diffusion terms are neglected.
2.2. Ionneutral collision frequency
The coefficient of friction between ions and neutrals (in the case of same temperature) is calculated as (Braginskii 1965) (19)where T is the plasma temperature, m_{i} (m_{n}) the ion (neutral atom) mass, m_{in} = m_{i}m_{n}/(m_{i} + m_{n}) is reduced mass, n_{i} (n_{n}) the ion (neutral atom) number density, is ionneutral collision cross section, and k = 1.38 × 10^{16} erg K^{1} is the Boltzmann constant. Using mean atomic cross section cm^{2} and the atomic masses of hydrogen and helium, the expression can be rewritten for each collision pair as (20)where T and n are normalized by 1 K and cm^{3}, respectively.
The collision frequency between ions and neutrals needs additional discussion. It is commonly accepted that the collision frequency of ions with neutrals ν_{in} = α_{in}/m_{i}n_{i} is generally different than the collision frequency of neutrals with ions ν_{ni} = α_{in}/m_{n}n_{n} (De Pontieu et al. 2001; Vranjes et al. 2008). In the first (second) case the neutrals (ions) are supposed to have zero initial velocity, and this approximation leads to the different values for the collision frequency. However, mean collision frequency between moving ions and neutral atoms is supposed to be the same value on physical grounds. From simple equations of motions of ions and neutrals (without pressure and Lorentz terms), one can derive the equation for relative velocity between ions and neutrals (21)where V_{i} and V_{n} are velocities of ions and neutral atoms, respectively. This equation gives the collision frequency between ions and neutrals as (22)This expression significantly changes the value of collision frequency between ions and neutrals. For example, calculating the collision frequency in the parameters used in Zaqarashvili et al. (2011, Eq. (30)) gives 16 s^{1} instead of 4 s^{1}.
The chromospheric temperature of 1.6 × 10^{4} K yields the ionneutral collision frequencies (as expressed by Eq. (22)) as (23)where we used proton, neutral hydrogen, singly ionized helium, and neutral helium number densities of 9.09 × 10^{10} cm^{3}, 1.05 × 10^{10} cm^{3}, 6.9 × 10^{9} cm^{3}, and 2.5 × 10^{9} cm^{3} (Fontenla et al. 1993, model FAL933).
One may expect that the effects of ionneutral collision may be enhanced for the waves with frequencies near ionneutral collision frequency (Zaqarashvili et al. 2011). Therefore, higher frequency waves are probably more affected by ionneutral collisions.
Fig. 1 Height dependence of atmospheric parameters according to FAL933 model (Fontenla et al. 1993). Upper panel: the plasma temperature. Middle panel: blue solid line – the ratio of neutral hydrogen and electron number densities, n_{H}/n_{e}; green dashed line – the ratio of neutral helium and electron number densities, n_{He}/n_{e}. Lower panel: the ratio of neutral helium and neutral hydrogen number densities, n_{He}/n_{H}. 
2.3. Partial ionization in the solar chromosphere
Before studying the damping of linear Alfvén waves due to ionneutral collision, we briefly discuss the plasma ionization in the solar lower atmosphere. The plasma is weakly ionized in the photosphere/lower chromosphere, but becomes more and more ionized with height. The increase in temperature with height leads to the ionization of hydrogen and helium atoms. Hydrogen and helium are almost fully ionized in the solar corona, so the transition occurs near the region of sharp temperature rise.
Figure 1 shows the dependence of plasma parameters on height according to the FAL933 model (Fontenla et al. 1993). This model includes the dependence of ionization degree on heights for both hydrogen and helium. The upper panel shows the plasma temperature vs. height according to this model. The temperature minimum is located near 500 km above the basic of the photosphere, while the transition region is just above 2000 km in height. The neutral hydrogen number density is much higher than the electron number density at the lower heights, but becomes comparable near ~1900 km, which corresponds to the temperature of 9400 K (middle panel). Hydrogen atoms quickly become ionized above this height. The neutral helium number density is also higher than the electron number density at the lower heights (middle panel). They become comparable near ~1600 km, which corresponds to the temperature of ~7300 K. The ratio becomes smaller and smaller just above this height. The important parameter for this study is the ratio of neutral helium and neutral hydrogen number densities (the lower panel). This ratio stays nearly constant (~0.1) up to 1500 km height, then it quickly increases up to 0.2 at height of 2000 km. Therefore, the effect of neutral helium in the damping of Alfvén waves should become important at higher altitudes.
In the next section we study the damping of linear Alfvén waves due to ion collision with neutral hydrogen and helium atoms.
3. Damping of linear Alfvén waves
We consider the simplest case of a static and homogeneous plasma with homogeneous unperturbed magnetic field, B_{z}, directed along the z axis. We assume the wave propagation in xz plane, i.e. ∂/∂y = 0. Then the Alfvén waves are governed by the equations where v_{y} (v_{Hy}, v_{Hey}) are the perturbations of ion (neutral hydrogen, neutral helium) velocity, b_{y} is the perturbation of the magnetic field, and ρ_{0} (ρ_{H0}, ρ_{He0}) is the unperturbed ion (neutral hydrogen, neutral helium, respectively) density. Here we use the definitions (28)Collision between neutral hydrogen and neutral helium atoms, expressed by α_{HeH}, does not contribute significantly to the damping of Alfvén waves. Moreover, α_{HeH}/α_{HeH + } ~ n_{H}/n_{H + } < 1 in spicules and in prominence plasma, therefore we neglect the corresponding terms for simplicity. Then, Fourier analysis with exp [i(k_{z}z − ωt) ] gives the dispersion relation (29)where (30)The dispersion relation has four different roots: two complex solutions, which correspond to Alfvén waves damped by ionneutral collision, and two purely imaginary solutions, which correspond to damped vortex solutions of neutral hydrogen and neutral helium fluids.
To study the effects of neutral helium, the damping rates of Alfvén waves obtained from Eq. (29) should be compared with the damping rates obtained in the absence of neutral helium. The dispersion relation corresponding to Alfvén waves in twofluid plasma with only neutral hydrogen can be obtained easily from Eq. (29) (see also Zaqarashvili et al. 2011) (31)This dispersion relation has three different roots: two complex solutions, which correspond to damped Alfvén waves, and a purely imaginary solution. The purely imaginary solution is associated with the neutral fluid. This is clearly seen if one neglects the ionneutral collision, so the neutral fluid behaves independently. Then, the fluid has vortex, ω = 0, solution in the incompressible limit, which gains an imaginary part when collision is included. As a result, the purely imaginary solution of dispersion relation Eq. (31) corresponds to the damped vortex solution of neutral hydrogen fluid (Zaqarashvili et al. 2011). Then, it is straightforward that the fourth imaginary solution in Eq. (29) arises owing to the vortex solution of neutral helium fluid as mentioned above.
Figure 2 displays the damping rates of Alfvén waves in partially ionized isothermal plasma with (derived from Eq. (29)) and without (derived from Eq. (31)) neutral helium vs. normalized Alfvén frequency for three different temperatures in the chromosphere. The Alfvén frequency k_{z}v_{A} is normalized on the collision frequency between neutral hydrogen and charged heavy particles, ν_{H} = ν_{H + H} + ν_{He + H}. The values of temperature and number densities are taken from the FAL933 model (Fontenla et al. 1993). The lower panel shows the damping rates for the temperature of 8100 K with proton, neutral hydrogen, and neutral helium number densities of 1.33 × 10^{11} cm^{3}, 3.71 × 10^{11} cm^{3}, and 4.8 × 10^{10} cm^{3}, respectively. The middle panel shows the damping rates for the temperature of 10 000 K with proton, neutral hydrogen, and neutral helium number densities of 1.07 × 10^{11} cm^{3}, 8.47 × 10^{10} cm^{3} and 1.13 × 10^{10} cm^{3}, respectively. The upper panel shows the damping rates for the temperature of 16 000 K with proton, neutral hydrogen, and neutral helium number densities of 9.09 × 10^{10} cm^{3}, 1.05 × 10^{10} cm^{3}, and 2.5 × 10^{9} cm^{3}. The three different temperatures correspond to 1790 km, 1995 km, and 2015 km heights from the photosphere.
Fig. 2 Damping rate (imaginary part of frequency) of Alfvén waves, ϖ = ω/k_{z}v_{A}, vs. normalized Alfvén frequency, a = k_{z}v_{A}/ν_{H}, for three different values of chromospheric temperature. Blue crosses correspond to the damping rates due to ion collision with neutral hydrogen atoms only, while red asterisks correspond to the damping rates due to ion collision with both, neutral hydrogen and neutral helium atoms. Red (blue) solid line corresponds to the damping rate derived in the singlefluid approach (Eq. (32)) with (without) neutral helium. 
According to Fig. 2 the influence of neutral helium atoms on the damping of Alfvén waves significantly depends on the plasma temperature. The effect of neutral helium is less pronounced for a temperature of ~ 8000 K (lower panel), but it gradually increases for higher temperatures (middle and upper panels). The dependence of the damping rate on frequency is similar to the one for the neutral hydrogen: the damping rate increases at lower frequencies, reaches its maximum at a particular frequency, and then begins to decrease for higher values (Zaqarashvili et al. 2011). The damping rate of Alfvén waves has its maximum near the ionneutral collision frequency.
The results show that the neutral helium atoms significantly enhance the damping of Alfvén waves in the chromospheric plasma. For example, the upper panel of Fig. 2 corresponds to the parameters of spicules with the temperature of 16 000 K (Beckers 1968). The presence of neutral helium enhances the damping rate almost three times for lowfrequency waves and almost twice for highfrequency ones. Therefore, neutral helium atoms should be taken into account, along with neutral hydrogen, in the study of transverse waves in solar spicules (Zaqarashvili & Erdélyi 2009).
Fig. 3 Damping rate (imaginary part of frequency) of Alfvén waves, ϖ = ω/k_{z}v_{A}, vs. normalized Alfvén frequency, a = k_{z}v_{A}/ν_{H}, for the prominence plasma. Red asterisks (blue crosses) correspond to the solutions with (without) neutral helium. Red (blue) solid line corresponds to the damping rate derived in the singlefluid approach (Eq. (32)) with (without) neutral helium. 
Recently, Soler et al. (2010) have studied the damping of MHD waves due to ionneutral collision in the singlefluid description of prominence plasma including neutral helium atoms. They conclude that the neutral helium has no significant influence on the damping of MHD waves. We check the result in the threefluid MHD description by taking plasma conditions typical prominence cores with 8000 K temperature. Proton and neutral hydrogen number densities are taken as 9 × 10^{9} cm^{3} and 2.1 × 10^{10} cm^{3}, respectively (Labrosse et al. 2010). The number density of neutral helium is taken as 10 % of the neutral hydrogen, i.e. 2.1 × 10^{9} cm^{3}, and the number density of singly ionized helium is taken as 2.1 × 10^{8} cm^{3}. Figure 3 displays the damping rates of Alfvén waves in partially ionized isothermal plasma with and without neutral helium. We see that the neutral helium has an insignificant effect on the wave damping in considering parameters, especially for lowfrequency waves, as shown by Soler et al. (2010). However, the neutral helium can be important in prominencetocorona transition regions, where the plasma temperature is higher than in cores. The damping rates derived in single and multifluid approaches are similar for lowfrequency waves, but show completely different behavior when wave frequency approaches the ionneutral collision frequency.
4. Discussion
Solar chromosphere/photosphere and prominences contain significant amounts of neutral atoms, which may lead to the damping of MHD waves. In some cases, the damping can be important so should be taken into consideration. The effect of neutral hydrogen on the damping of MHD waves has been studied in both singlefluid (Khodachenko et al. 2004; Forteza et al. 2007) and multifluid approximations (Zaqarashvili et al. 2011). However, the collision of ions with neutral helium atoms is important for some cases as shown here. Hydrogen atoms quickly become ionized with increasing temperature, while the number of neutral helium may still be significant. This is clearly seen in Fig. 4, which is plotted according to the atmospheric model FAL933. The ratio of neutral helium to neutral hydrogen number densities grows significantly in the temperature interval 10 000−40 000 K.
We study the damping of Alfvén waves due to collision of ions with neutral hydrogen and neutral helium atoms. The threefluid description of plasma is used, where one component is electronprotonsingly ionized helium and the other two are neutral hydrogen and neutral helium gases. The dynamics of linear Alfvén waves in homogeneous isothermal plasma with constant unperturbed magnetic field is considered. The exact dispersion relation for the Alfvén waves is derived and solved for different plasma parameters. The damping rates due to ion collision with neutral hydrogen and neutral helium atoms are derived and compared with those obtained only for neutral hydrogen. The analysis shows that the collision between ions and neutral helium can be important at certain values of plasma temperature, when the hydrogen already begins ionization, but the neutral helium is still presented. This happens for T ~ 10 000−40 000 K, which corresponds to upper chromosphere, spicules, and prominencecorona transition regions. The presence of neutral helium enhances the wave damping a few times compared to the damping only due to ionneutral hydrogen collisions (Fig. 2, middle and upper panels).
Fig. 4 The ratio of neutral helium and neutral hydrogen number densities, n_{He}/n_{H} vs. plasma temperature according to FAL933 model (Fontenla et al. 1993). 
The collisions between ions and neutral helium atoms can be important in solar spicules, which correspond to the upper panel on Fig. 2. The plasma temperature, electron density, and ionization degree are almost constant along spicules (especially in their upper part), therefore the isothermal plasma is a good approximation there (Beckers 1968). The neutral helium atoms may significantly enhance the damping of transverse waves, which are continuously observed in spicules (Kukhianidze et al. 2006; De Pontieu et al. 2007; Zaqarashvili et al. 2007; see also recent review of Zaqarashvili & Erdélyi 2009).
The presence of neutral helium has no significant influence on Alfvén wave damping in prominence cores (Fig. 3). This makes the results of Soler et al. (2010) still valuable. They studied the damping of MHD waves in a singlefluid approach. However, the neutral helium can be important in the prominencecorona transition region, where the plasma temperature is higher than in the prominence cores.
Damping rates derived in the threefluid approach reach a peak near the ionneutral collision frequency and then decrease for higher frequencies unlike the singlefluid approach, where the damping increases linearly with increasing frequency (Figs. 2–3; see also Zaqarashvili et al. 2011). This may require the reconsideration of highfrequency Alfvén wave damping in the chromosphere. More sophisticated simulations taking the gradients of temperature and ionization degrees with height into account in the different parts of solar atmosphere are needed.
In order to understand the wave damping due to ionneutral collision, we consider the simplest case of partially ionized plasma consisting of electrons, protons, and neutral hydrogen. Then, the wellknown expression of damping rate in the singlefluid approach (Braginskii 1965; Khodachenko et al. 2004; Forteza et al. 2007) can be rewritten as (32)where ρ_{tot} is the total density (ion+neutral), ν_{in} the ionneutral collision frequency defined by Eq. (22), and ξ_{n},ξ_{i} are normalized by ρ_{tot}. This expression clearly indicates that the normalized damping rate depends on the ratio of Alfvén and ionneutral collision frequencies and the ratio of neutral and ion fluid densities. The earlier statement that the damping rate depends on the magnetic field strength can be translated as follows: the increase in the magnetic field leads the increase in the Alfvén speed, therefore the waves with the same wavenumber have a higher frequency, which is closer to ionneutral collision frequency, so the normalized damping rate increases. In fact, it is the ratio of neutral and ion fluid densities that determines the damping rate of particular wave harmonic. Highfrequency waves are damped quickly. However, this statement is valid for a lowfrequency wave spectrum below ionneutral collision frequency.
FAL933 model does not include the coronal values of plasma parameters, and the ionization degree of helium in the corona is unknown. However, it is clear that the coronal plasma contains only a small number of neutral helium, therefore the effects of neutral helium should be weaker in the corona. However, a specific influence of neutral helium on wave damping in the corona is not ruled out, so further analysis is required to make a firm conclusion.
5. Conclusions
The existence of neutral helium atoms, alongside neutral hydrogen, significantly enhances the damping of Alfvén waves at certain intervals of plasma temperature due to ionneutral collision as compared to the case of neutral hydrogen alone. Neutral helium causes the increase in damping rates few times in spicules and prominencecorona transition region at the temperature T ~ 10 000 − 40 000 K. Alfvén waves have maximal damping rates at some frequency interval peaking near the ionneutral collision frequency. At the same time, the damping rate is reduced for higher frequencies.
The study of highfrequency Alfvén wave damping in the solar chromosphere with a realistic height profile of ionization degree needs to be revised in the future, including both neutral hydrogen and neutral helium atoms.
Acknowledgments
The work was supported by the Austrian Fonds zur Förderung der wissenschaftlichen Forschung (project P21197N16) and by the European FP7PEOPLE2010IRSES269299 project SOLSPANET. T.V.Z. also acknowledges financial support from the Georgian National Science Foundation (under grant GNSF/ST09/4310). The authors thank the referee for constructive comments.
References
 Arber, T. D., Haynes, M., & Leake, J. E. 2007, ApJ, 666, 541 [NASA ADS] [CrossRef] [Google Scholar]
 Beckers, J. M. 1968, Sol. Phys., 3, 367 [NASA ADS] [Google Scholar]
 Braginskii, S. I. 1965, Rev. Plasma Phys., 1, 205 [NASA ADS] [Google Scholar]
 Carbonell, M., Forteza, P., Oliver, R., & Ballester, J. L. 2010, A&A, 515, 80 [Google Scholar]
 De Pontieu, B., & Haerendel, G. 1998, A&A, 338, 729 [NASA ADS] [Google Scholar]
 De Pontieu, B., Martens, P. C. H., & Hudson, H. S. 2001, ApJ, 558, 859 [NASA ADS] [CrossRef] [Google Scholar]
 De Pontieu, B., McIntosh, S. W., Carlsson, M., et al. 2007, Science, 318, 1574 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Fontenla, J. M., Avrett, E. H., & Loeser, R. 1993, ApJ, 355, 700 [Google Scholar]
 Forteza, P., Oliver, R., Ballester, J. L., & Khodachenko, M. L. 2007, A&A, 461, 731 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Forteza, P., Oliver, R., & Ballester, J. L. 2008, A&A, 492, 223 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Goedbloed, H., & Poedts, S. 2004, Principles of Magnetohydrodynamics (Cambridge University Press) [Google Scholar]
 Gogoberidze, G., Voitenko, Y., Poedts, S., & Goossens, M. 2009, ApJ, 706, 12 [NASA ADS] [CrossRef] [Google Scholar]
 Goodman, M. L., & Kazeminezhad, F. 2010, ApJ, 708, 268 [NASA ADS] [CrossRef] [Google Scholar]
 Haerendel G. 1992, Nature, 360, 241 [NASA ADS] [CrossRef] [Google Scholar]
 James, S. P., & Erdélyi, R. 2002, A&A, 393, L11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 James, S. P., Erdélyi, R., & De Pontieu, B. 2004, A&A, 406, 715 [Google Scholar]
 Kukhianidze, V., Zaqarashvili, T. V., & Khutsishvili, E. 2006, A&A, 449, L35 [Google Scholar]
 Khodachenko, M. L., & Rucker, H. O. 2005, Adv. Space Res., 36, 1561 [NASA ADS] [CrossRef] [Google Scholar]
 Khodachenko, M. L., & Zaitsev, V. V. 2002, Astrophys. Space Sci., 279, 389 [NASA ADS] [CrossRef] [Google Scholar]
 Khodachenko, M. L.,Arber, T. D., Rucker, H. O., & Hanslmeier, A. 2004, A&A, 422, 1073 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Labrosse, N., Heinzel, P., Vial, J.C., et al. 2010, Space Sci. Rev., 151, 243 [Google Scholar]
 Leake, J. E., & Arber, T. D. 2006, A&A, 450, 805 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Leake, J. E., Arber, T. D., & Khodachenko, M. L. 2005, A&A, 442, 1091 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Singh, K. A. P., & Krishan, V. 2010, New Astron., 15, 119 [NASA ADS] [CrossRef] [Google Scholar]
 Soler, R., Oliver, R., & Ballester, J. L. 2009a, ApJ, 699, 1553 [NASA ADS] [CrossRef] [Google Scholar]
 Soler, R., Oliver, R., & Ballester, J. L. 2009b, ApJ, 707, 662 [NASA ADS] [CrossRef] [Google Scholar]
 Soler, R., Oliver, R., & Ballester, J. L. 2010, A&A, 512, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tsap, Y. T., Stepanov, A. V., & Kopylova, Y. G. 2011, Sol. Phys., 270, 205 [NASA ADS] [CrossRef] [Google Scholar]
 Vranjes, J., Poedts, S., Pandey, B. P., & de Pontieu, B. 2008, A&A, 478, 553 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zaqarashvili, T. V., & Erdélyi, R. 2009, Space Sci. Rev., 149, 355 [NASA ADS] [CrossRef] [Google Scholar]
 Zaqarashvili, T. V., Khutsishvili, E., Kukhianidze, V., & Ramishvili, G. 2007, A&A, 474, 627 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zaqarashvili, T. V., Khodachenko, M. L., & Rucker, H. O., 2011, A&A, 529, A82 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Figures
Fig. 1 Height dependence of atmospheric parameters according to FAL933 model (Fontenla et al. 1993). Upper panel: the plasma temperature. Middle panel: blue solid line – the ratio of neutral hydrogen and electron number densities, n_{H}/n_{e}; green dashed line – the ratio of neutral helium and electron number densities, n_{He}/n_{e}. Lower panel: the ratio of neutral helium and neutral hydrogen number densities, n_{He}/n_{H}. 

In the text 
Fig. 2 Damping rate (imaginary part of frequency) of Alfvén waves, ϖ = ω/k_{z}v_{A}, vs. normalized Alfvén frequency, a = k_{z}v_{A}/ν_{H}, for three different values of chromospheric temperature. Blue crosses correspond to the damping rates due to ion collision with neutral hydrogen atoms only, while red asterisks correspond to the damping rates due to ion collision with both, neutral hydrogen and neutral helium atoms. Red (blue) solid line corresponds to the damping rate derived in the singlefluid approach (Eq. (32)) with (without) neutral helium. 

In the text 
Fig. 3 Damping rate (imaginary part of frequency) of Alfvén waves, ϖ = ω/k_{z}v_{A}, vs. normalized Alfvén frequency, a = k_{z}v_{A}/ν_{H}, for the prominence plasma. Red asterisks (blue crosses) correspond to the solutions with (without) neutral helium. Red (blue) solid line corresponds to the damping rate derived in the singlefluid approach (Eq. (32)) with (without) neutral helium. 

In the text 
Fig. 4 The ratio of neutral helium and neutral hydrogen number densities, n_{He}/n_{H} vs. plasma temperature according to FAL933 model (Fontenla et al. 1993). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.