arxiv: v1 [astro-ph.sr] 23 May 2016

Transcription

1 Astronomy& Astrophysics manuscript no. ms c ESO 2018 February 20, 2018 The role of Alfvén wave heating in solar prominences Roberto Soler 1, 2, Jaume Terradas 1, 2, Ramon Oliver 1, 2, and Jose Luis Ballester 1, 2 1 Departament de Física, Universitat de les Illes Balears, E Palma de Mallorca, Spain. 2 Institut d Aplicacions Computacionals de Codi Comunitari (IAC 3 ), Universitat de les Illes Balears, E Palma de Mallorca, Spain. roberto.soler@uib.es arxiv: v1 [astro-ph.sr] 23 May 2016 Received XXX; accepted XXX ABSTRACT Observations have shown that magnetohydrodynamic waves over a large frequency range are ubiquitous in solar prominences. The waves are probably driven by photospheric motions and may transport energy up to prominences suspended in the corona. Dissipation of wave energy can lead to heating of the cool prominence plasma, so contributing to the local energy balance within the prominence. Here we discuss the role of Alfvén wave dissipation as a heating mechanism for the prominence plasma. We consider a slab-like quiescent prominence model with a transverse magnetic field embedded in the solar corona. The prominence medium is modelled as a partially ionized plasma composed of a charged ion-electron single fluid and two separate neutral fluids corresponding to neutral hydrogen and neutral helium. Friction between the three fluids acts as a dissipative mechanism for the waves. The heating caused by externally-driven Alfvén waves incident on the prominence slab is analytically explored. We find that the dense prominence slab acts as a resonant cavity for the waves. The fraction of incident wave energy that is channelled into the slab strongly depends upon the wave period, P. Using typical prominence conditions, we obtain that wave energy trapping and associated heating are negligible when P 100 s, so that it is unlikely that those waves have a relevant influence on prominence energetics. When 1 s P 100 s the energy absorption into the slab shows several sharp and narrow peaks, that can reach up to 100%, when the incident wave frequency matches a cavity resonance of the slab. Wave heating is enhanced at those resonant frequencies. Conversely, when P 1 s cavity resonances are absent, but the waves are heavily damped by the strong dissipation. We estimate that wave heating may compensate for about 10% of radiative losses of the prominence plasma. Key words. Magnetohydrodynamics (MHD) Sun: atmosphere Sun: corona Sun: filaments, prominences Sun: oscillations Waves 1. INTRODUCTION The energy balance in solar prominences, and the understanding of the processes involved with heating and cooling of the plasma, are difficult problems that are intimately linked to the prominence formation and structure (see, e.g., Gilbert 2015; Heinzel 2015). Using the differential emission measure, Parenti & Vial (2007) found that for a prominence the integral over the temperature range of 10 4 K 10 6 K yielded radiative losses of ergs cm 2 s 1. In quiescent conditions the prominence radiation is steady, requiring a source of still unknown heating to maintain energy balance in the structure (Parenti 2014). Identifying the source or sources of prominence heating is a challenging task for both observers and theoreticians, and several heating mechanisms have been proposed (see Labrosse et al. 2010; Heinzel et al. 2010; Gilbert 2015). Incident radiation is generally accepted as the dominant prominence heating mechanism (see Gilbert 2015). The amount of radiative heating in prominences largely depends on the illumination from the surrounding solar atmosphere. For this reason, although the prominence plasma is akin to that in the chromosphere, the studies of chromospheric radiative heating (see, e.g., Carlsson & Leenaarts 2012) are not entirely applicable to the particular conditions of prominences. Because prominences are suspended in and surrounded by the corona, coronal radiation is absorbed from all sides, while they are somewhat more isolated from the photospheric radiation than the chromospheric plasma (Gilbert 2015). The radiation field plays a crucial role and has to be properly considered when evaluating the balance between radiative heating and radiative losses (see Heinzel 2015). Several works tried to construct prominence models by taking into account the balance between incident radiation and cooling. Heasley & Mihalas (1976) obtained radiativeequilibrium temperatures of about 4,600 K in slab-like prominences. Anzer & Heinzel (1999) considered similar models and concluded that some form of additional heating is necessary to balance radiative losses inside the prominence, specially in the prominence cool core. In addition, computations by Heinzel et al. (2010); Heinzel & Anzer (2012) point out that not only the radiation field is important but the accurate description of radiative losses is also crucial to estimate the radiative-equilibrium temperature in prominences. Heinzel & Anzer (2012) considered a more accurate description of prominence cooling by adding the calcium radiative losses and computed equilibrium temperatures even lower than those obtained for a pure hydrogen gas. In general, the obtained radiative-equilibrium temperatures in the works cited above are lower than what is typically inferred from observations. The values reported in the literature, as derived applying different methods, span in a range of temperatures between 7,500 K and 9,000 K in the prominence cores (see Parenti 2014). An additional, non-negligible source of heating seems to be necessary to rise the central prominence temperature up to the expected values (Labrosse et al. 2010; Heinzel et al. 2010; Heinzel & Anzer 2012). In this regard, dissipation of magneto- Article number, page 1 of 15

2 A&A proofs: manuscript no. ms hydrodyamic (MHD) wave energy is proposed as another possible source of prominence heating to be considered in addition to the dominant radiative heating. The mechanism of wave heating in the solar atmosphere has been recently reviewed by Arregui (2015). The actual importance of prominence wave heating, however, remains to be explored in detail. High-resolution observations have unquestionably demonstrated the presence of MHD waves in solar prominences (e.g., Lin et al. 2007, 2009; Hillier et al. 2013). Also, non-thermal velocities derived from line widths have been interpreted as a signature of unresolved waves (Parenti & Vial 2007). The nonthermal velocities observed by Parenti & Vial (2007) were about 10 km s 1 for prominence core temperatures. Statistical studies performed by Hillier et al. (2013) on the properties of transverse waves in thin threads of a quiescent prominence revealed wave periods ranging from 50 s to 6000 s, and velocity amplitudes between 0.2 km s 1 and 23 km s 1, with a significant amount of wave energy found at short periods. We note that, owing to observational constraints, periods shorter than 50 s could not be detected by Hillier et al. (2013), but it is likely that such short periods are also present in prominences (see a discussion on this issue in Pécseli & Engvold 2000). Indeed, periods of about 30 s have been reported in prominences from the analysis of Doppler time series (Balthasar et al. 1993; Kolobov et al. 2008). We also note that the typical lifetime of individual threads is 20 min, so that the detection of longer periods in prominence fine structures should be interpreted with caution. The driver of the waves has been related to photospheric granular motions, so that the driver would be an agent external to the prominence itself. In this direction, Hillier et al. (2013) compared the power spectra of the prominence waves with that obtained from photospheric horizontal motions and obtained a good correlation. Thus, MHD waves might transport energy from the photosphere up to prominences suspended high in the corona. The question then arises, can this wave energy be deposited in the prominence plasma in the form of heat? The observational evidence of MHD waves in prominences and their theoretical understanding are well established (see Arregui et al. 2012; Ballester 2015). However, detailed studies of the actual role of the waves in prominence heating are scarce and their conclusions are dissimilar. For instance, Pécseli & Engvold (2000) studied the possible role of Alfvén wave dissipation in providing a supporting force for the prominence material against gravity (see also Jensen 1983, 1986) and estimated that the associated wave heating would only contribute with a minute fraction to the bulk radiative energy. Conversely, Parenti & Vial (2007) assumed that the non-thermal velocities measured in a quiescent prominence were caused by the presence of Alfvén waves and performed a simple estimation of the Alfvén wave energy flux. Then, Parenti & Vial (2007) suggested that radiative losses for prominence-corona transition region (PCTR) temperatures could be entirely compensated by Alfvén wave heating alone (see also Parenti 2014). We must note, however, that Parenti & Vial (2007) probably overestimated the wave heating rate by assuming that all the wave energy flux is converted into heat. As pointed out by Gilbert (2015), detailed investigations of the impact of wave heating in prominences are necessary. The ability of MHD waves to heat the plasma relies on the presence of a dissipative process capable of efficiently damping the waves and depositing their energy into the medium. The efficiency of a given dissipative mechanism usually depends on the temporal and spatial scales involved (see, e.g., Soler et al. 2015b). A comparative study of the role of various damping mechanisms of MHD waves was made by Khodachenko et al. (2004, 2006). In the case of plasmas with prominence conditions, Khodachenko et al. (2004, 2006) found that damping by ion-neutral collisions is dominant over other processes as, e.g., viscosity or magnetic resistivity. Thus, ion-neutral collisions can be the required mechanism for efficient dissipation of MHD waves in prominences. The damping of MHD waves in partially ionized prominence plasmas has been theoretically investigated in a number of works (e.g., Forteza et al. 2007; Soler et al. 2009b; Carbonell et al. 2010, among others), although no elaborated connection between wave damping and plasma heating has been provided so far. In the solar chromosphere, Alfvén wave damping due to ionneutral collisions is a subject of intense research (e.g., Haerendel 1992; De Pontieu et al. 2001; Leake et al. 2005; Soler et al. 2015a). It has been suggested that the heating associated to the Alfvén wave damping may be strong enough to partly compensate chromospheric radiative losses (see, e.g., Goodman 2011; Song & Vasyliūnas 2011; Russell & Fletcher 2013; Tu & Song 2013; Reep & Russell 2016; Arber et al. 2016). In view of the similarities between chromospheric and prominence plasmas, the purpose of this paper is to perform an exploratory study of the role of Alfvén waves in prominence heating. Another relevant question is whether the energy of waves driven at the phostosphere can be efficiently channelled into prominences. Basic results on the reflection and transmission of Alfvén waves incident on a contact discontinuity (see Ferraro 1954) indicate that most of the wave energy incident on prominences should be reflected because of the large prominence-tocorona density ratio. However, Hollweg (1984a) showed that the depression in the ambient Alfvén velocity caused by a dense plasma slab, such as a prominence, acts as a resonant cavity for incident Alfvén waves (see also Hollweg 1981; Ionson 1982; Zhugzhda & Locans 1982; Hollweg 1984b; Sterling & Hollweg 1984). Hollweg (1984a) applied his results to coronal loops, but the essential physics would be the same in the case of prominence slabs. Due to the existence of certain resonant frequencies associated to nearly standing waves in the cavity, large wave energy fluxes can be transmitted into the slab from the external medium, so enhancing the efficiency of wave energy transmission at those resonant frequencies. The physical origin of the cavity resonances studied by Hollweg (1984a) resides in the interference of the waves. Complex interference patterns can be formed because of the wave reflection and transmission. A resonance occurs when the waves transmitted into the cavity form a nearly standing wave and, in addition, there is a destructive interference between the reflected waves and the waves that leak out of the cavity. Thus, we note that cavity resonances have a different physical origin than the so-called Alfvén continuum resonances caused by the variation of the plasma properties across the magnetic field direction (see, e.g., Goossens et al. 2011). Here we shall use a simple slab model with a transverse magnetic field to represent a prominence embedded in the solar corona and shall explore the combined impact of cavity resonances and ionneutral dissipation on the wave heating of prominences. To the best of our knowledge, the role of cavity resonances has not been investigated before in the case of prominences. This paper is organized as follows. Section 2 contains the basic equations that govern the propagation of linear Alfvén waves in a partially ionized plasma. In Section 3 we embark on the investigation of the heating by Alfvén waves in a prominence slab model. Finally, Section 4 contains the conclusions and some ideas for future works. Article number, page 2 of 15

3 Soler et al.: Alfvén wave heating in prominences 2. BASIC EQUATIONS Here we give the general governing equations of linear Alfvén wave propagation in a homogeneous partially ionized prominence plasma with a straight and constant magnetic field, B. We assume that the plasma constituents are hydrogen and helium. Because of the relatively low prominence core temperature (T 10 4 K), we assume that hydrogen is only partially ionized while, for simplicity, helium is assumed to be fully neutral. The abundance of ionized helium at prominence core temperatures is very small, so that its influence on the behavior of Alfvén waves is expected to be of minor importance Multi-fluid plasma For the study of wave propagation, we consider the so-called multi-fluid formalism (see, e.g., Zaqarashvili et al. 2011b). It is assumed that ions and electrons form a ion-electron single fluid, while neutral hydrogen and neutral helium form two separate neutral fluids. We denote byρ i,ρ H, andρ He the density of ions, neutral hydrogen, and neutral helium, respectively. Hence the total density is ρ ρ i +ρ H +ρ He, (1) where the mass density of electrons is neglected owing to the very small electron mass. The hydrogen ionization ratio is ρ i ξ i =. (2) ρ i +ρ H We considerξ i to be an arbitrary parameter whose value is between 0<ξ i < 1. The helium to hydrogen ratio of abundances by number is A He A H = n He n i + n H, (3) where n i, n H, and n He are the number densities of ions, neutral hydrogen, and neutral helium, respectively. We also consider A He /A H to be an arbitrary parameter. Hence, the densities of the individual fluids as functions of the total density are ξ i ρ i = ρ, 1+4A He /A H (4) ρ H = 1 ξ i ρ, 1+4A He /A H (5) ρ He = 4A He /A H ρ. 1+4A He /A H (6) The various fluids exchange momentum by means of particle collisions. Electron-neutral collisions are of minor importance for the present study and are not taken into account here. The friction coefficient determines the strength of the coupling between fluids. The expression of the friction coefficient for collisions involving neutrals in the approximation of small velocity drift is (e.g., Braginskii 1965; Draine 1986) α ββ = ρ βρ β 4 m β + m β 3 σ 8k B T ββ πm ββ (7) whereβandβ stand for i, H, or He, T is the temperature (we assume the same temperature for all the components of the plasma), m β is the particle mass of species β, m ββ = m β m β / ( m β + m β ) is the reduced mass, kb is Boltzmann s constant, andσ ββ is the momentum-transfer cross section. We take σ ih m 2 and σ ihe m 2 from Vranjes & Krstic (2013), andσ HHe m 2 from Lewkow et al. (2012) Governing equations of linear waves To study linear waves, we assume small perturbations superimposed on the background plasma and linearize the general, non-linear multi-fluid equations (see, e.g., Zaqarashvili et al. 2011b,a; Khomenko et al. 2014a). For simplicity, we ignore the effect of gravity and also assume that there are no equilibrium flows. The linearized equations considered here for the study of Alfvén waves are ρ i v i t ρ H v H t v He ρ He t b t = 1 µ ( b) B α ih (v i v H ) α ihe (v i v He ), (8) = α ih (v H v i ) α HHe (v H v He ), (9) = α ihe (v He v i ) α HHe (v He v H ), (10) = (v i B), (11) where v i, v H and v He are the velocity perturbations of ions, neutral hydrogen, and neutral helium, respectively, b is the magnetic field perturbation, µ is the magnetic permeability, and the rest of parameters have been defined before. In these equations we have not considered gas pressure perturbations because they are irrelevant for Alfvén waves. Other effects neglected in these equations are viscosity and magnetic resistivity, whose role on the wave damping has been shown to be less important than that of ion-neutral collisions in a partially ionized prominence plasma (Khodachenko et al. 2004, 2006). The induction equation (Equation (11)) is derived from the generalized Ohm s Law obtained after neglecting electron inertia. The full expression of the generalized Ohm s Law in a partially ionized plasma can be found in, e.g., Zaqarashvili et al. (2011b); Khomenko et al. (2014a); Soler et al. (2015a). Here, we use the ideal version of Ohm s Law, in which the term corresponding to electron drift (i.e., Hall s effect) is neglected. This term has been shown to be unimportant concerning the behavior of MHD waves in a partially ionized prominence plasma (see Soler et al. 2009c). Specifically, Hall s effect is negligible for wave frequencies lower than the ion cyclotron frequency, namelyω i = eb/m i, with B the magnetic field strength. For a typical prominence magnetic field strength of 10 G, we haveω i 96 khz, which is much higher than the frequency of the observed waves in prominences (see, e.g., Hillier et al. 2013) Wave energy and heating rate An expression for the evolution of the energy of the linear waves can be obtained by following a similar procedure to that described in, e.g., Goedbloed & Poedts (2004) and Walker (2005) for the case of ideal MHD waves in a fully ionized plasma. We take the dot product of Equations (8), (9), and (10) with v i, v H, and v He, respectively, and also take the dot product of Equation (11) with b/µ. Next, we add the equations and the resulting expression governing the wave energy evolution can be cast as U t + Π= Q, (12) Article number, page 3 of 15

4 A&A proofs: manuscript no. ms where U,Π, and Q are given by U = 1 2 ρ i v i ρ H v H ρ He v He µ b 2, (13) Π = 1 µ [(B b) v i (v i b) B], (14) Q = α ih v i v H 2 +α ihe v i v He 2 +α HHe v H v He 2. (15) Equation (12) shows that the wave energy is not conserved. We discuss the meaning of the various terms in Equation (12). The first term on the left-hand side accounts for the temporal evolution of the wave energy density, U, which is the sum of the kinetic energies of ions, neutral hydrogen, and neutral helium, plus the magnetic energy. The second term on the left-hand side is the divergence of the wave energy flux, Π, that informs us about the amount of energy that propagates with the wave. We note that only the velocity of ions appears in the expression of Π, so that Equation (14) is formally identical to that in fully ionized plasmas (see Goedbloed & Poedts 2004). The right-hand side of Equation (12) depends on the quantity Q (Equation (15)), which accounts for the loss of wave energy owing to dissipation due to collisions (see also Braginskii 1965; Draine 1986). It is implicit in Equation (12) that, as the wave propagates, its energy is absorbed into the plasma, where it is eventually transformed into heat. The energy lost by the wave must be gained by the plasma and converted into internal energy. This could be seen in the full version of the energy equation, in which the changes of the background plasma were also included (see Goedbloed & Poedts 2004, for the energy equation in a fully ionized plasma). The full energy equation remains a conservation equation. Wave energy is converted into internal energy of the background plasma, but the sum is constant. The plasma heating itself cannot be studied with the present linear analysis, but Equation (12) informs about the rate at which the wave energy is lost. We can therefore identify Q as the wave heating rate Alfvén wave dispersion relation We consider monochromatic, parallel propagating Alfvén waves and set the perturbations as proportional to exp ( iωt + iks), whereωis the wave frequency, k is the wavenumber in the direction of the background magnetic field, and s denotes the coordinate along the background magnetic field lines. The frequency,ω, is assumed to be real. Because of wave damping, the wavenumber, k, is complex. The imaginary part of the wavenumber defines the characteristic damping length of the waves. For incompressible Alfvén waves, Equations (8) (11) can be recast as follows, ωv i = kb b iν ih (v i v H ) iν ihe (v i v He ), (16) µρ i ωv H = iν Hi (v H v i ) iν HHe (v H v He ), (17) ωv He = iν Hei (v He v i ) iν HeH (v He v H ), (18) ωb = kb v i, (19) whereν ih,ν Hi,ν ihe,ν Hei,ν HHe, andν HeH are collision frequencies defined as ν ih ν ihe ν HHe = α ih ρ i, ν Hi = α ih ρ H, (20) = α ihe, ν Hei = α ihe (21) ρ i ρ He = α HHe, ν HeH = α HHe. (22) ρ H ρ He Article number, page 4 of 15 From Equations (17) and (18) we find that the velocities of neutral hydrogen and neutral helium are related to the ion velocity by v H = v He = iν Hi (ω+iν He ) ν HHe ν Hei v i, (ω+iν H ) (ω+iν He )+ν HHe ν HeH (23) iν Hei (ω+iν H ) ν HeH ν Hi v i, (ω+iν H ) (ω+iν He )+ν HHe ν HeH (24) withν H =ν Hi +ν HHe andν He =ν Hei +ν HeH the neutral hydrogen and neutral helium total collision frequencies, respectively. We combine Equation (16) with Equations (19), (23), and (24) to remove the other perturbations and arrive at an equation for v i alone (not given here). The condition that v i 0 provides us with the dispersion relation, namely k 2 = ωω col, (25) c 2 A where c 2 A is the square of the ion Alfvén velocity (which only includes the ion density) given by c 2 A = B2 µρ i, (26) andω col is a redefined, complex-valued frequency that incorporates the effects of collisions with neutrals and is given by Ω col ω+iν i + ν ihν Hi (ω+iν He )+iν ih ν HHe ν Hei (ω+iν H ) (ω+iν He )+ν HHe ν HeH + ν iheν Hei (ω+iν H )+iν ihe ν HeH ν Hi (ω+iν H ) (ω+iν He )+ν HHe ν HeH, (27) whereν i =ν ih +ν ihe is the ion total collision frequency. In the case that the plasma is fully ionized,ω col =ω and Equation (25) simplifies to the well-known dispersion relation of ideal Alfvén waves, namely k 2 = ω2. (28) c 2 A The wavenumber is real in a fully ionized ideal plasma Wave energy deposition We use Equation (14) to compute the Alfvén wave energy flux. Because of the Fourier dependence exp ( iωt + iks) the wave variables are complex and the physical meaningful part resides in the real part. In addition, to drop the oscillatory behavior of the waves and retain the net energy flux, we averageπ over one wave period, P = 2π/ω. Thus, the time-averaged energy flux, Π, for a forward-propagating wave is Π = B2 Re (k) 2µ ω V2 exp [ 2Im (k) s] ê B, (29) where ê B is the unit vector in the direction of the background magnetic field, Re and Im denote the real and imaginary parts, respectively, k is given by the principal square root of Equation (25), and V is the amplitude of the ion velocity perturbation at the arbitrary position s=0. Both Re (k) and Im (k) are positive. The units of Π are W m 2. Equation (29) agrees with Equation (8) of Hollweg (1984a).

5 To check Equation (29) we consider a fully ionized plasma, so that Im (k)=0 and Re (k)=ω/c A. Then, Equation (29) consistently reverts to the well-known expression of the ideal Alfvén wave energy flux (e.g., Walker 2005), namely Π = 1 2 ρ iv 2 c A ê B. (30) Equation (29) evidences that the net energy transported by the Alfvén wave exponentially decreases as the wave propagates through the plasma. This behavior is a consequence of the damping due to collisions. A characteristic energy absorption exponential length scale, L A, can be defined as Soler et al.: Alfvén wave heating in prominences L A = 1 2Im (k). (31) The wave energy is absorbed, i.e., deposited into the plasma in a length scale of the order of L A. In order to estimate L A, we consider typical physical conditions in quiescent prominences. We takeρ= kg m 3, T= 8000 K, and B=10 G. In addition, we take A He /A H = 0.1 (Gilbert et al. 2002). Figure 1 displays the energy absorption length, L A, computed from Equation (31) as a function of the wave period, P, for various values of the hydrogen ionization ratio,ξ i. As expected, the lower the hydrogen ionization ratio, the stronger the wave damping because of the larger presence of neutrals in the plasma. Waves with short periods are very efficiently damped in the prominence plasma. The corresponding energy absorption lengths for these short-period waves are rather small. For instance, L A is between 100 km and 1000 km for waves with periods between 1 s and 10 s, and L A keeps decreasing when P is further reduced. Conversely, the wave damping, and so the energy absorption, is less efficient for waves with longer periods. We obtain that L A increases as the wave period is increased. For instance, for wave periods of the order of 100 s, the energy absorption length is as large as 10 6 km. In practice, this result suggests that the energy deposition associated to waves with periods of the order of minutes is not significant over realistic length scales in solar prominences (see, e.g, Engvold 2015). We shall confirm this result later. To determine the relative importance of neutral helium on the wave damping, we have repeated the computation of Figure 1 but when collisions with neutral helium are neglected. For the sake of simplicity, these results are not shown here. We obtain that collisions with helium are of minor importance in determining the value of L A because collisions of ions with neutral hydrogen is the dominant damping mechanism of the waves. 3. WAVE HEATING IN A PROMINENCE SLAB In Section 2.5 we have learnt that short-period Alfvén waves can efficiently damp in prominence plasmas. Here, we investigate whether the energy of Alfvén waves incident on a solar prominence can be deposited in the prominence medium because of this damping. To do so, we represent the prominence by a slab with a transverse magnetic field embedded in the solar corona and assume the presence of Alfvén waves that are incident on the prominence-corona interface (see Figure 2). In this work we do not care about the driver of the waves, which has been associated to photospheric motions (Hillier et al. 2013). We assume, however, that the wave driver is external to the prominence. Here we are interested in studying the interaction of the incident waves with the dense and cool, partially ionized slab. Fig. 1. Energy absorption length, L A, as a function of the wave period, P=2π/ω, for an Alfvén wave propagating in a prominence plasma when different values of the hydrogen ionization ratio,ξ i, are assumed. The meaning of the various linestyles is indicated within the figure. The parameters used in the computations areρ= kg m 3, T = 8000 K, B=10 G, and A He /A H = 0.1. We note that logarithmic scale has been used in both axes. The incident waves will be partly transmitted into the prominence slab and partly reflected back to the corona. The waves that penetrate into the prominence will travel back and forth through the whole prominence slab, and part of their energy will be leaked again to the corona. The superposition of these transmitted and reflected waves can cause complex interference patterns. Previous investigations by Hollweg (1984a) showed that the dense slab acts as a resonant cavity for the incident Alfvén waves. Cavity resonances may appear when the waves transmitted into the slab interfere so that a nearly standing wave is formed and, at the same time, the waves reflected back to the corona and the waves that leak out of the slab interfere destructively. In such a case, most of the incident wave energy is channelled into the slab. In addition, since the waves are dissipated in the partially ionized prominence plasma because of ion-neutral collisions, the wave energy fed into the cavity resonances may efficiently contribute to the heating of the plasma. Essentially, the condition for the excitation of cavity resonances is that the frequency of the incident waves matches a natural frequency (eigenmode) of the prominence. The eigenmodes of a prominence modelled as a dense plasma slab have been investigated in the literature (e.g., Joarder & Roberts 1992; Oliver et al. 1993), where the dependence of the natural frequencies on the physical and geometrical properties of the slab was explored. Here, we shall give an approximate expression for the Alfvén cavity frequencies in Section 3.4. The excitation of an internal Alfvén eigenmode naturally produces large wave amplitudes within the prominence slab. In the stationary state, the energy of the wave driver must be distributed among the internal and external plasmas. Hence, if a cavity resonance is hit, most of the available energy goes to the internal plasma and, consequently, the amount of energy in the exterior is much smaller. Physically, the low wave amplitude in the external medium can be understood in terms of the destructive interference mentioned before. Article number, page 5 of 15

6 A&A proofs: manuscript no. ms Corona Prominence Corona Assuming a typical helium abundance of A He /A H = 0.1, Equation (34) results in µ c 0.6. On the other hand, in the prominence we consider the hydrogen ionization degree to be arbitrary and helium to be fully neutral. Hence, µ p = 1+4 A He A H 1+ξ i + A He A H. (35) For instance, if A He /A H = 0.1 andξ i = 0.5, Equation (35) gives µ p = Fig. 2. Sketch of the prominence slab model Slab model B We assume that the equilibrium state is composed of a uniform partially ionized plasma slab of width d, representing a solar prominence, embedded into a uniform and unbounded fully ionized medium representing the solar coronal environment. When necessary, we shall use the subscripts p and c to refer to the plasma in the prominence slab and in the corona, respectively. We use Cartesian coordinates and set the coordinate frame so that the equilibrium is infinite in the y- and z-directions, while the prominence slab is bounded in the x-direction. The center of the slab is located at x=0. The equilibrium is permeated by a straight magnetic field along the x-direction, namely B= Bê x, so that the magnetic field is transverse to the prominence slab. Similar slab models have been used in the past to investigate global oscillations (eigenmodes) of prominences (see, e.g., Joarder & Roberts 1992; Oliver et al. 1993; Anzer 2009; Soler et al. 2009a). A sketch of the slab model can be seen in Figure 2. Since the equilibrium magnetic field is uniform, the thermal pressure is required to be continuous across the prominencecorona interface. We use the ideal gas law to compute the coronal, P c, and prominence, P p, equilibrium pressures as P c = ρ crt c µ c, P p = ρ prt p µ p, (32) while the requirement that the thermal pressure is continuous across the prominence-corona interface gives the relation between the densities and temperatures of the two plasmas as ρ c ρ p = T p/ µ p T c / µ c. (33) In these equations, R is the ideal gas constant, and µ c and µ p are the mean atomic weight in the corona and in the prominence, respectively. On the one hand, the coronal plasma is fully ionized so that µ c = 1+4 A He A H 2+3 A He A H. (34) Article number, page 6 of Matching conditions at the prominence-corona boundary We consider Alfvén waves propagating in the corona that are incident on the prominence-corona interface. Since both the prominence slab and the corona are homogeneous, the Alfvén wave dispersion relations (Equation (25) in the prominence and Equation (28) in the corona) are locally satisfied. The incident waves are partly transmitted into the prominence slab and partly reflected back to the corona. Here we find the velocity amplitudes of the transmitted and reflected waves in terms of the amplitudes of the incident waves. The process is similar to that followed by Hollweg (1984a). In the steady state, and assuming that the waves are plane polarized along the y-direction, we can write the ion velocity perturbation as I 1 exp (ik c x)+r 1 exp ( ik c x), if x< d 2, v iy (x)= T 1 exp ( ik p x ) + T 2 exp ( ik p x ), if x < d 2, (36) I 2 exp ( ik c x)+r 2 exp (ik c x), if x> d 2, where I 1 and I 2 are the amplitudes of the incident waves, R 1 and R 2 are the amplitudes of the waves reflected back to the corona, and T 1 and T 2 are the amplitudes of the waves transmitted into the prominence slab. In addition, k c and k p are the parallel wavenumber in the corona and in the prominence, respectively. On the one hand, k c is given by taking the positive square root of the right-hand side of Equation (28) and using coronal conditions. On the other hand, k p is given by taking the principal square root of the right-hand side of Equation (25) and using prominence conditions. We assume that the amplitudes of the incident waves, I 1 and I 2, are arbitrary. Then, the amplitudes R 1, R 2, T 1, and T 2 are determined by imposing appropriate boundary conditions at x=±d/2. In the case of an Alfvén wave incident on a contact discontinuity, the boundary conditions are the continuity of the tangential components of the magnetic and electric fields (see, e.g., Goedbloed & Poedts 2004). In our case, these conditions simply reduce to the requirements [[ ]] [[ ]] viy viy = 0, = 0, (37) x where [[X]] denotes the variation of the quantity X across the prominence corona-interface. We impose the boundary conditions given in Equation (37) at x = ±d/2 and, after some algebraic manipulations, we find that the amplitudes of the reflected and transmitted waves as func-

7 Soler et al.: Alfvén wave heating in prominences tions of I 1 and I 2 are R 1 = R 2 = T 1 = T 2 = 4k p k c p k2 c) I (k p +k c) i(k2 sin ( k (k p +k c) 2 p d ) I 1 [( ) ] ( ) 2 kp k 1 c ( k p +k c exp i2kp d ) exp i kp k c d, (38) 4k p k c p k 2 c) I (k p +k c) i(k2 sin ( k (k p +k c) 2 p d ) I 2 [( ) ] ( ) 2 kp k 1 c ( k p +k c exp i2kp d ) exp i kp k c d, (39) 2k c k p +k c I 1 + 2k c(k p k c) 1 exp ( ik (k p +k c) 2 p d ) I 2 ( ) 2 kp k c ( k p +k c exp i2kp d ) exp 2k c k p +k c I 2 + 2k c(k p k c) 1 exp ( ik (k p +k c) 2 p d ) I 1 ( ) 2 kp k c ( k p +k c exp i2kp d ) exp [ i ( ) ] d k p k c,(40) 2 [ i ( ) ] d k p k c.(41) 2 We note that Hollweg (1984a) only considered wave incidence on the left interface of the slab, so his case can be recovered by setting I 2 = 0 in our expressions. The amplitudes of the reflected and transmitted waves are complex valued and depend upon the wave frequency through the expression of the wavenumbers. It is worth noting that in the absence of dissipation the amplitudes of the reflected and transmitted Alfvén waves on a contact discontinuity are real and independent of the frequency (see, e.g., Ferraro 1954). In general, the consideration of dissipative mechanisms causes the amplitudes of the reflected and transmitted waves to depend upon the wave frequency Energy absorption coefficient The fraction of incident wave energy that is reflected back to the corona is determined by the coefficient of energy reflectivity defined as R= Π R 1 + Π R2 Π I1 + Π I2, (42) where Π I1 and Π I2 denote the moduli of the time-averaged energy flux of the incident waves, and Π R1 and Π R2 are the corresponding quantities of the reflected waves. The energy flux is given by Equation (29). In terms of the wave amplitudes, the expression of the reflectivity simply becomes R= R 1R 1 + R 2R 2 I 1 I1 + I 2I2, (43) where the superscript denotes the complex conjugate. We assume that, in general, the amplitudes of the incident waves, I 1 and I 2, are complex valued, so that we can consider an arbitrary phase relation between the two incident waves. The values ofr range betweenr=0 (no reflection) andr=1 (full reflection). In addition, the fraction of incident wave energy that is transmitted and trapped into the prominence slab can be computed by invoking conservation of energy. So, the fraction of energy absorbed into the prominence is A=1 R. (44) The values ofarange betweena=0 (no energy absorption) and A = 1 (complete absorption). To avoid confusion, we must warn the reader that the present definitions of the reflection and absorption coefficients are different from those used by Hollweg (1984a). We define the reflectivity as the fraction of incident energy that returns to the corona, and the absorption as the fraction of incident energy that gets trapped within the slab. Hollweg (1984a) only considered incidence on one boundary of the slab and defined the transmissivity as the fraction of incident wave energy that emerges from the opposite slab boundary. To perform a numerical application, we take B=10 G and A He /A H = 0.1. In the prominence plasma, we use the same physical conditions as before, namelyρ p = kg m 3 and T p = 8000 K. We assume the density in the corona to be ρ c =ρ p /200= kg m 3, and the coronal temperature is computed from the pressure balance condition (Equation (33)). For instance, ifξ i = 0.5, we get T c 1.1 MK. We note, however, that the coronal temperature plays no role in the following results. In addition, we take d= 5,000 km as the width of the prominence slab. Figure 3 displays the coefficient of energy absorption, A, for a wide range of incident wave periods (in logarithmic scale). The results for three different values ofξ i have been stacked one on top of the other for better comparison. In these computations we considered I 1 = I 2. The logarithmic scale used in Figure 3 evidences the existence of three different regimes depending upon the wave period: The large-period range, in which the energy trapping into the prominence slab is negligible and practically all the wave energy returns to the corona. This regime is independent of the value ofξ i and happens roughly for P 100 s. The intermediate-period range, roughly for 1 s P 100 s, which is characterized by the presence of sharp absorption peaks that can reach up toa 1. The location of the absorption peaks is independent ofξ i. However, the height, width, and number of peaks seem to be affected byξ i. We can relate the presence of these absorption peaks with the cavity resonances studied by Hollweg (1984a). We shall explore the physics behind cavity resonances in Section 3.4. The short-period range, approximately for P 1 s, where the absorption stabilizes at about A 0.3 and slowly grows when P decreases. No absorption peaks are present here and the value ofξ i does not have an important effect. The results of Figure 3 indicate that it is unlikely that waves with P 100 s impact on prominence energetics. The energy of these large-period waves is not trapped within the prominence. On the other hand, although the energy trapping of short-period waves with P 1 s is more important, most of the energy is reflected back to the corona. Conversely, the existence of large absorption peaks in the intermediate range of periods indicates that those resonances may efficiently trap energy within the slab. Cavity resonances are, therefore, excellent candidates to be involved in strong wave heating. Importantly, we note that the coefficient of energy absorption as a function of the frequency is an intrinsic property of the model for a given set of parameters. The energy absorption does not depend on the velocity amplitude of the incident waves Cavity resonances Hollweg (1984a) provided an instructive study about the nature of cavity resonances. We can follow Hollweg (1984a) to understand the physics behind cavity resonances. According to Article number, page 7 of 15

8 A&A proofs: manuscript no. ms Fig. 3. Coefficient of energy absorption,a, as a function of the wave period, P=2π/ω, for various prominence hydrogen ionization ratio:ξ i = 0.9 (top graphs),ξ i = 0.5 (mid graphs), andξ i = 0.1 (bottom graphs). We considered I 1 = I 2 and d= 5,000 km. We note that the horizontal axis is in logarithmic scale. The parameters used in the computations are given in the text. Hollweg (1984a), a resonance occurs when the waves transmitted into the slab interfere so that a nearly standing wave is formed. We can obtain an approximation to the resonant-cavity frequencies by neglecting the effect of dissipation (i.e., assuming a perfect coupling between the fluids). The approximation can be cast as ω n nπ d B µρp, with n=1, 2, 3... (45) To derive this approximation we have assumed that the standing Alfvén wave is perfectly trapped within the slab. Equation (45) corresponds to the frequencies of the internal Alfvén eigenmodes of the slab. The eigenmodes of a prominence slab embedded in the solar corona have been investigated by, e.g., Joarder & Roberts (1992) and Oliver et al. (1993). In addition to the internal eigenmodes, there are external and hybrid modes that depend on the line-tying conditions at the base of the corona. Since in the present model the corona is assumed to be unbounded, external and hybrid modes are absent. Anyway, neither hybrid nor external modes are expected to play an important role in prominence heating. Hybrid modes have smaller frequencies (longer periods) than all internal modes, so that their damping by ion-neutral collisions would be inefficient. In turn, the amplitude of external modes and their associated energy within the prominence slab is very small. Article number, page 8 of 15

9 Soler et al.: Alfvén wave heating in prominences Fig. 4. Coefficient of energy absorption, A, as a function of the incident wave frequency (in Hz), f = ω/2π. (a) Results when the incident waves on both sides of the slab are in phase (I 1 = I 2, black solid line) and anti-phase (I 1 = I 2, dashed blue line). (b) Results when I 2 = I 1 e iπ/2. In both panels, the vertical dotted lines denote the approximate cavity resonance frequencies according to Equation (45). In all cases we usedξ i = 0.2 and the remaining parameters are the same as in Figure 3. Figure 4(a) shows the coefficient of energy absorption as a function of the incident wave frequency (in Hz), f = ω/2π, for a particular set of parameters given in the caption of the figure. In that figure, we use the frequency and not the period in the horizontal axis to better distinguish the resonance peaks. We use vertical lines to indicate the location of the first 31 resonances approximately given by Equation (45). We see that the location of the absorption peaks agrees very well with the approximate formula. In Figure 4(a) we consider two different cases, namely I 1 = I 2 and I 1 = I 2. The first case assumes that incident waves on both sides are in phase, while in the second case the incident waves are in anti-phase. These two cases are considered to illustrate that the symmetry of the incident waves determines whether their interference within the slab is constructive or destructive, so that only the resonances with the appropriate symmetry are excited. We see in Figure 4(a) that only odd resonances (i.e., n=1, 3, 5...) are excited when I 1 = I 2, while only even resonances (i.e., n=2, 4, 6...) appear when I 1 = I 2. The two cases considered in Figure 4(a) represent paradigmatic situations. In reality, we should not expect any special relation between the amplitudes and phases of the waves incident on both sides of the slab. Therefore, odd and even resonances would be naturally excited by an arbitrary driver. This is seen in Figure 4(b), where we took I 2 = I 1 e iπ/2. We note that the height Article number, page 9 of 15

10 A&A proofs: manuscript no. ms of the absorption peaks in Figure 4(b) is essentially half that of the peaks in Figure 4(a). This is so because the total number of resonances is twice that of the paradigmatic cases of Figure 4(a). Hence, wave energy absorption is now distributed among twice as many absorption peaks. Equation (45) also informs us about the dependence of the resonant frequencies on the prominence density,ρ p, the magnetic field strength, B, and the width of the prominence slab, d. Figure 5 shows the coefficient of energy absorption for various values of the slab width (panel a) and the prominence density (panel b). Both results are in agreement with Equation (45). The decrease/increase of d causes the resonant peaks to move to shorter/longer periods. Changing the density has a similar impact, but the density also influences the overall reflectivity of the slab, as happens in the case of a single discontinuity (Ferraro 1954). The smaller the prominence density, the smaller the reflectivity, and so the larger the fraction of energy that gets trapped within the slab. This effect is specially noticeable in the short-period regime, but the density does not influence much the quality of the resonances. Hollweg (1984a) showed that, in the absence of dissipation, leakage of the nearly standing wave to the external medium due to imperfect trapping determines the width and the maxima of the resonance peaks. The more important the leakage, the broader the resonance width and the smaller the resonance peak. Then, Hollweg (1984a) further added the effect of an arbitrary dissipation mechanism and obtained that dissipation has a similar effect as that of leakage on the shape of the resonances. The results of Hollweg (1984a) can be used to qualitatively explain the shape of the absorption resonances obtained in the present work. Again, we use Figure 4 to support the discussion. Dissipation by ion-neutral collisions increases as the wave frequency increases. Accordingly, Hollweg (1984a) suggests that the width of the resonances should grow as the frequency of the resonance increases. Figure 4 clearly demonstrates this effect: the resonance peaks get wider with increasing frequency. Concerning the peaks height, we see that the first six resonance peaks (i.e., those with the lowest frequencies) grow when the frequency increase. Conversely, the height of the resonance peaks with n>6 decrease as the frequency increases. According to Hollweg (1984a), dissipation would be responsible for the decrease of the resonance peaks as the frequency increases. In turn, the increase of the height of the low-frequency resonances could be attributed to leakage becoming less efficient as the frequency increases. A side effect of the decrease of the resonance height when the frequency increases is that resonances are in practice suppressed for high enough frequencies or, equivalently, short enough periods. This explains why resonance peaks are absent in the shortperiod range of Figure 3. The very strong dissipation that occurs for high frequencies (short periods) has the negative consequence of suppressing the cavity resonances and so reducing the efficiency of wave energy absorption into the slab. Asξ i increases, dissipation becomes more important and the effective number of resonances decreases Associated heating The presence of cavity resonances in the intermediate range of periods dramatically enhances the transmission and trapping of wave energy into the prominence slab compared to the classic result of wave incidence on a simple plasma discontinuity (Ferraro 1954). The purpose here is to determine whether the ion-neutral dissipation of the wave energy can provide significant plasma heating. The heating rate caused by the Alfvén wave dissipation, Q, is given by Equation (15). Using the ion velocity amplitude of the waves transmitted into the slab given by Equations (40) and (41), we compute the heating rate within the slab. In order to retain the net heating, we average Q over one wave period to drop the oscillatory dependence. After some algebraic manipulations, we obtain the expression of the time-averaged heating, namely Q = 1 2 ρ pλ { T 1 T 1 exp[ 2Im ( k p ) x ] + T2 T 2 exp[ 2Im ( k p ) x ]}, (46) with ξ i ν ih Λ = 1 q 2 1 ξ i ν ihe 1+4A He /A H p + 1 q A He /A H p + (1 ξ i)ν HHe q 1 q A He /A H p, (47) p = (ω+iν H ) (ω+iν He )+ν HHe ν HeH, (48) q 1 = iν Hi (ω+iν He ) ν HHe ν Hei, (49) q 2 = iν Hei (ω+iν H ) ν HeH ν Hi. (50) Equation (46) evidences that the heating profile within the prominence slab is the superposition of two exponential functions that are maximum at the slab boundaries. Hence, Q is a function of position. The spatial dependence of the heating rate can be dropped by spatially averaging Q within the slab. The result is Q vol = 1 d/2 Q dx=ρ p Λ sinh[ Im ( ) ] k p d d d/2 Im ( T 1 ) T 1 + T 2 T 2. k p d 2 (51) The units of Q vol are W m 3, so that Q vol is the timeaveraged, spatially-averaged volumetric heating rate. The volumetric heating, Q vol, depends on the parameters of the model, on the wave frequency/period, and on the velocity amplitude of the incident waves. In fact, as shown in Section 3.3, for a given model parameters the intrinsic energy absorption depends only on the wave frequency. However, to compute the actual heating information about the power spectrum of the waves is also needed. The spectral weighting function relating the amplitude of the waves with their frequency needs to be known, but the power spectrum of the waves that are incident on prominences is not known. This poses an important limitation to estimate the actual heating rate, since the only alternative left is assuming ad hoc the amplitude of the waves. A first possibility is assuming a flat power spectrum and using a constant velocity amplitude. This is done in Figure 6, where we represent the volumetric heating as a function of the period when a constant velocity amplitude of 10 km s 1 is used for the incident waves. We use this velocity because it corresponds to the non-thermal velocities for T 10 4 K measured by Parenti & Vial (2007). Cavity resonances are clearly seen in Figure 6 as specific periods for which wave heating is enhanced. We consider two different values ofξ i to explore the effect of ionneutral dissipation. Reducingξ i increases the heating efficiency of the cavity resonances with long periods but reduces the efficiency of the resonances with short periods. Conversely, heating in the range of short periods (P<1 s), where resonances are not present, shows little dependence onξ i. Article number, page 10 of 15