Issue 
A&A
Volume 532, August 2011



Article Number  A94  
Number of page(s)  20  
Section  The Sun  
DOI  https://doi.org/10.1051/00046361/201117014  
Published online  29 July 2011 
Toward detailed prominence seismology
II. Charting the continuous magnetohydrodynamic spectrum
^{1}
FOM Institute for Plasma Physics Rijnhuizen, Association EURATOMFOM, PO Box 1207, 3430 BE Nieuwegein, The Netherlands
email: JW.S.Blokland@Rijnhuizen.nl
^{2}
Centre for Plasma Astrophysics, Mathematics Department, KU Leuven, Celestijnenlaan 200B, 3001 Heverlee, Belgium
Received: 4 April 2011
Accepted: 20 May 2011
Context. Starting from accurate magnetohydrodynamic flux rope equilibria containing prominence condensations, we initiate a systematic survey of their linear eigenoscillations. This paves the way for more detailed prominence seismology, which thus far has made dramatic simplifications about the prevailing magnetic field topologies.
Aims. To quantify the full spectrum of linear MHD eigenmodes, we require knowledge of all fluxsurface localized modes, charting out the continuous parts of the MHD spectrum. We combine analytical and numerical findings for the continuous spectrum for realistic prominence configurations, where a cool prominence is embedded in a hotter cavity, or where the flux rope contains multiple condensations supported against gravity.
Methods. The equations governing all eigenmodes for translationally symmetric, gravitating equilibria containing an axial shear flow, are analyzed, along with their fluxsurface localized limit. The analysis is valid for general 2.5D equilibria, where either density, entropy, or temperature vary from one flux surface to another. We analyze the intricate mode couplings caused by the poloidal variation in the flux rope equilibria, by performing a small gravity parameter expansion. We contrast the analytical results with continuous spectra obtained numerically.
Results. For equilibria where the density is a flux function, we show that continuum modes can be overstable, and we present the stability criterion for these convective continuum instabilities. Furthermore, for all equilibria, a fourmode coupling scheme between an Alfvénic mode of poloidal mode number m and three neighboring (m − 1,m,m + 1) slow modes is identified, occurring in the vicinity of rational flux surfaces. For realistically structured prominence equilibria, this coupling is shown to play an important role, from weak to stronger gravity parameter g values. The analytic predictions for small g are compared with numerical spectra, and progressive deviations for larger g are identified.
Conclusions. The unstable continuum modes could be relevant for shortlived prominence configurations. The gaps created by poloidal mode coupling in the continuous spectrum need further analysis, as they form preferred frequency ranges for global eigenoscillations.
Key words: Sun: filaments, prominences / instabilities / magnetohydrodynamics (MHD) / plasmas
© ESO, 2011
1. Prominence seismology from MHD modeling
Motivated by the recent interest and progress made in diagnosing oscillations in solar prominences, we have started an effort (Blokland & Keppens 2011, hereafter Paper I) to close the gap between observed prominence properties, and the models used to quantify their oscillations. From observations, increasingly detailed knowledge of thermodynamical properties in filaments is available. This relies on inverting spectroscopic data, where departures from local thermodynamic equilibrium in the radiative transfer through one and twodimensional prominence models are nowadays routinely incorporated (Labrosse et al. 2010). At the same time, high resolution observations (Lin et al. 2005), computations of threedimensional (3D) magnetohydrostatic equilibrium configurations without gravity (Aulanier et al. 2002), or modern magnetic field extrapolations from photospheric magnetograms (DeRosa et al. 2009), together with the application of spectropolarimetric techniques (Lin et al. 1998; Trujillo Bueno et al. 2002), have brought insight into the complex threedimensional topology of the surrounding magnetic field. Simultaneously, increased spatiotemporal resolution has quantified a host of oscillation modes observed in prominences, with velocity amplitudes of a few km s^{1}. Both short and longperiod signals, from minutes to hours, can be detected, and regions that oscillate coherently with similar periods allow us to quantify damping rates for seemingly global eigenoscillations (Terradas et al. 2002).
Theoretical studies in prominence seismology have meanwhile identified the most likely candidates explaining these motions, in terms of fundamental sausage, kink, and other mode types present in the MHD spectrum of magnetized slabs or cylinders. To explain the observed damping of the eigenmodes, one can on the one hand study wave damping as a result of nonadiabatic effects, such as anisotropic, thermal conduction, and radiative losses. For cylindrically structured coronal loops and prominences, with prevailing helical fields, early work by Keppens et al. (1993) demonstrated how these nonadiabatic effects can give rise to damped and/or overstable global, discrete nonadiabatic Alfvén waves. Soler et al. (2008) revisited these effects for eigenmodes of more simplified, uniform cylinder models in homogeneous magnetic fields. By using uniform internal/external conditions, the possibility of wave resonant behavior was excluded artificially. This resonant behavior, where global waves have eigenfrequencies in the range of the MHD continuous spectrum, provides an alternative means of explaining the observed damping rates. The theory for resonant MHD waves is reviewed in many papers, e.g. Goossens et al. (2011), and references therein. In other reviews of prominence oscillations (e.g. Mackay et al. 2010; Arregui & Ballester 2011), many of the accumulated theoretical findings are discussed, and a striking observation is that all these efforts dramatically simplify the underlying magnetic topology, in order to allow analytic progress. Nonuniformity, and the associated possibility of resonant absorption, is at most introduced by varying the density profiles (e.g. Soler et al. 2010), sticking to uniform field topologies. At the same time, efforts to account for the multistranded nature of actual solar filaments have also been initiated, in similarly simple magnetic fields (Luna et al. 2010).
To pave the way for more detailed prominence seismology, true MHD spectroscopy (Goedbloed et al. 1993), of the kind pioneered in fusionrelated plasma physics, is called for. This relies on accurate computations of the underlying equilibria, with full account of the complex magnetic topology, and subsequent computation of its eigenspectrum. In our accompanying Paper I, we started this effort by demonstrating our ability to compute accurate translationally symmetric flux rope equilibria containing prominence condensations, governed by the magnetohydrostationary equations balancing Lorentz forces, pressure gradients, and gravity. We here continue this effort, and use these equilibria, where we incorporate helical fields, axial shear flows, and varying plasma beta conditions consistent with observed filament properties. We here concentrate on quantifying all eigenmodes defining the continuous parts of the MHD eigenspectrum. This is a worthwhile effort in its own right, partly because the global, discrete modes are known to be organized about this continuous spectrum (see e.g. modern textbooks such as Goedbloed & Poedts 2004; Goedbloed et al. 2010), and because the possible damping of these global modes due to resonant absorption relies on knowledge of the frequency ranges that they occupy.
Charting out the continuous MHD spectrum for realistic flux rope structures with prominences has never been attempted before, in sharp contrast to the vast body of knowledge on MHD eigenmodes in axisymmetric, tokamak configurations. Our effort considers flux ropes idealized to be invariant along the filament axis, with an outer cylindrical crosssection, for simplicity. Since the gravitational stratification will act to deform a pure cylindrical flux model to one with still nested, but downwards displaced flux surfaces, the resulting poloidal variation in the equilibrium quantities will induce coupling between the different branches of the continuous MHD spectrum. These branches are labeled by their axial, k, and azimuthal, m, mode numbers. When charting the mode frequencies from the filament axis to the flux rope edge, the spectrum can help us to determine the avoided crossings, which are akin to the poloidal mode couplings discussed for equilibria without gravity, such as the nested elliptic cylinders studied by Chance et al. (1977), or configurations with noncircular crosssections and/or density variations along field lines analyzed by Poedts & Goossens (1991). We extend these latter studies to translationally invariant, gravitating equilibria with possibly shear flow v_{z}(ψ), restricted to be a flux function. The actual equilibria are presented and discussed in Paper I. The analysis of their continuous spectra presented here, benefits from the detailed knowledge of the spectrum of the more complex, axisymmetric topologies (Goedbloed 1975). The case for selfgravitating, static axisymmetric magnetized tori was treated analytically in Poedts et al. (1985). Goedbloed et al. (2004b) analyzed of the continuous MHD spectrum for magnetized accretion tori, where external gravity, nested toroidal flux surfaces, and both poloidal and toroidal (Keplerianlike) rotation is included, and identified a potential for unstable localized modes called transslow Alfvén continuum modes. The approach of Goedbloed et al. (2004b), which combined analytic findings for the equilibria, the poloidal mode couplings that were confronted with numerically computed equilibria, and eigenmodes is analogous to our efforts in this series.
Our paper is then organized as follows. In Sect. 2, we present the equations governing all eigenmodes for these gravitating, 2.5D MHD equilibria. Since we allow for axial flow, the FriemanRotenberg formalism (Frieman & Rotenberg 1960) is revisited, and using the straight field line coordinates (see Paper I) describing our equilibria, we present them in a form where a field line projection can be made. This form can directly be compared with all other cases studied previously, such as the cylindrical limit, or the analogy with tokamak configurations that contain toroidal flow. In Sect. 3, we then derive the equations for the continuous spectrum, and find that prominence equilibria where density is a flux function can also be affected by convective continuum instabilities (Blokland et al. 2007a; van der Holst et al. 2000b,a). The coupled system of differential equations can be converted to algebraic relations between poloidal Fourier harmonics, which clearly identifies how Alfvén and slow modes can couple. This coupling is then studied analytically in Sect. 4, where the small gravity parameter expansion previously used (Paper I) to analyze equilibrium effects due to gravity, reenters to analyze the leading order coupling effects. For low frequencies, we then demonstrate that a gap can arise in the Alfvén continuum branches, from a fourmode coupling scheme. In Sect. 5, we briefly summarize the numerical approach used in the remaining sections. There, the equilibria presented in Paper I are all diagnosed spectroscopically, by charting out their continuous spectrum. The analytic findings are used to check and analyze in detail the derived eigenfrequencies. This is done for cool prominences embedded in hot cavities in Sect. 6, and for complex, double layered prominences without and with axial flows in Sect. 7. Conclusions and an outlook to future work is given in Sect. 8.
2. Spectral formulation
In what follows, we analyze the equations governing linear, ideal MHD perturbations about a stationary equilibrium. The equations for the equilibrium have been presented and solved numerically in Paper I. The (timeindependent) equilibrium quantities are denoted by the magnetic field vector B, the pressure p, the density ρ, and the external gravitational potential Φ, which for the prominence application at hand are all functions of the poloidal coordinates (x,y). The background equilibrium may contain a purely longitudinal shear flow v = v_{z}(ψ)e_{z}, where the axial flow is constant on each flux surface. Each flux label ψ identifies a single flux surface, which form nested surfaces defined by the magnetic field topology within the prescribed outer flux rope shape. We first recapitulate the governing equations, and then present a representation exploiting a projection involving straight field line coordinates.
2.1. FriemanRotenberg formalism
The formalism introduced by Frieman & Rotenberg (1960) is used here to investigate the stability properties of stationary MHD equilibria. From the full set of linearized MHD equations, these authors derived one secondorder differential equation for the Lagrangian displacement field vector ξ(1)where the generalized force operator G(ξ) is (2)Here, (3)is the force operator for static equilibria without gravity, which was derived by Bernstein et al. (1958). These expressions contain the Eulerian perturbation of the total pressure (4)where γ is the ratio of specific heats, and the Eulerian perturbation of the magnetic field (5)The timedependence of the displacement field is assumed to be exponential with normal mode frequencies ω, . Using this assumption, the FriemanRotenberg equation can be written as (6)We base our derivations below on the axial symmetry of the solar prominence equilibrium by writing the axial dependence of the displacement field as , where k is the axial mode number. We then neglect the hat on and suppress the implied exp( − iωt + ikz) dependence for all linear quantities.
2.2. Fieldline projection and representation
As in all stability studies in fusion research, we need to exploit a projection based on the nested topology of the magnetic surfaces, on which the helical field lines reside. This projection ensures that on each flux surface, one can make a distinction between two wave directions: one parallel and the other perpendicular to the magnetic field. In the short wavelength limit, these directions correspond to the polarizations of the slow and Alfvén waves, respectively. To perform this projection, we adopt straight field line, nonorthogonal coordinates (ψ,ϑ,z) that are introduced in our accompanying paper (Blokland & Keppens 2011), which are known once the equilibrium is fully determined. The triad of unit vectors based on the projection is (7)where the gradient of the magnetic flux ψ is perpendicular to the flux surfaces. Using these unit vectors, the components of the displacement field ξ can be replaced by the three related variables (8)and the projected gradient operators become The second equality in the expressions for ℱ and only holds when these operators act on a component of the displacement field ξ. These expressions contain several quantities already introduced in Paper I, such as the metric elements g_{ij} and the Jacobian J for the nonorthogonal straight field line coordinates, and quantities shown to be flux functions to realize equilibrium, such as the poloidal stream function I = B_{z} and the safety factor q(ψ) = JI.
The FriemanRotenberg equation in Eq. (1)can be written as (12)by exploiting the straight field line coordinates and applying the projection defined by Eq. (8), where is the fluxsurface localized Dopplershifted eigenfrequency. Furthermore, the matrices A and B are 3 × 3 matrix operators that are also present in the case of a static equilibrium without gravity. The matrix operator G contains the elements associated with the external gravitational potential. The nonvanishing matrix elements of A and B are (13)and (14)We used the poloidal curvature of the magnetic surfaces κ_{p} from Paper I, along with the operators (15)These are related to the normal gradient operator , originally introduced by Goedbloed (1997) for the spectral analysis of static tokamak equilibria. The square brackets in the expressions in Eq. (13)for the matrix elements of A and the gradient operators defined by (15)indicate that the differential operator only acts on the term inside the bracket. This notation is also used in the expressions below. As expected, the matrices A and B can also be derived from the expressions for wave modes about static axisymmetric MHD equilibria, found in the papers by Goedbloed (1975, 1997), by formally setting the toroidal distance R to unity. The matrix G enters if there is an external gravitational potential present. The expression for this matrix is(16)where We note that the term μ represents the pressure variation on a flux surface. This matrix G can also be derived from the expression for the matrix operator R found in the paper by Blokland et al. (2007a), where axisymmetric, gravitating, magnetized accretion tori with purely toroidal flow where spectroscopically analyzed. The equivalence involves (1) setting their radial distance R = 1; (2) producing the toroidal curvature κ_{t} = 0 for a translationally symmetric solar prominence; and (3) replacing the toroidal mode number n by the axial wave number k. In that analysis, as well as in the related one for toroidally rotating tokamak plasmas by van der Holst et al. (2000a), another matrix C related to the Coriolis effect appeared, which naturally vanishes for our translationally invariant configurations.
Another interesting, related limit is the cylindrical limit known from standard textbook treatments (Goedbloed & Poedts 2004). If we ignore gravity altogether, then matrix G will be zero and the operators ℱ and will just be algebraic expressions. In this case, the spectral equation given in Eq. (12)reduces exactly to the wellknown equation presented by Hain & Lüst (1958). They rewrote this matrix form to one secondorder differential equation for the radial displacement field, assuming a ratio of specific heats γ = 1. This was generalized by Goedbloed (1971) for arbitrary γ, and rewritten, for cases involving equilibrium flow in Hameiri (1981) and Bondeson et al. (1987). A cylindrical limit involving flow and external gravity, suitable for analyzing accretion disk stability, was presented by Keppens et al. (2002), and can also be easily shown to be directly related to the above expressions.
3. Continuous MHD spectrum
In the previous section, we derived the spectral equation in Eq. (12)governing all MHD waves and instabilities in translational symmetric solar prominences. This equation is the starting point to all following MHD spectral computations. In particular, we can derive the equations for the continuous MHD spectrum by first considering modes localized on a particular flux surface ψ = ψ_{0}. We can then obtain a nonsingular eigenvalue problem in the local Alfvén and slowlike eigenfunctions η(ϑ),ζ(ϑ). This form allows us to derive a stability criterion for continuum modes, where it is found that for an axially symmetric prominence, gravity can cause the continuous spectrum to become unstable. By expanding the periodic ϑdependence of the eigenfunctions, we get an equivalent algebraic set of coupled equations.
3.1. General formalism
To derive the equation for the continuous spectrum, we assume that the normal derivative (∂/∂ψ) of the eigenfunctions is large compared to the eigenfunctions themselves and the local Dopplershifted eigenfrequency is finite. Using these assumptions, the first row of the spectral equation in Eq. (12)can be solved approximately (19)We note that X is small compared to either Y or Z. This implies that the continuum modes are mainly tangential to a particular flux surface. Furthermore, we note that ∂X/∂ψ, Y, and Z are of the same order.
Inserting the approximated solution given by Eq. (19)into the second and third row of the spectral equation in Eq. (12)results in an eigenvalue problem that is independent of the normal derivative. Hence, the reduced problem becomes nonsingular. Exploiting this property, we can write the projected displacement field components Y and Z as (20)where δ(ψ − ψ_{0}) is the Dirac delta function. Here, Y and Z have been divided into an improper (ψdependence) and proper part (ϑdependence). This kind of splitting was introduced by Goedbloed (1975) for static axisymmetric equilibria without gravity. The reduced, nonsingular eigenvalue problem becomes (21)where(22)(23)(24)In the expressions above, κ_{g} is the geodesic curvature, (25)where κ = b·∇b is the field line curvature. Furthermore, is the magnetically modified BruntVäisälä frequency projected onto a flux surface, (26)This magnetically modified BruntVäisälä frequency has also been found in various related cases. For example, Poedts et al. (1985) and Beliën et al. (1997) found this frequency for static, axisymmetric gravitationally stratified plasmas. van der Holst et al. (2000b,a) discussed its role for tokamak plasmas with purely toroidal flow. Blokland et al. (2007a) generalized the expression by including an arbitrary external gravitational potential in its application to rotating, magnetized accretion disks. When we examine the matrix g closely, we see that the poloidal variation in the pressure is present in all its matrix elements. This has the consequence that the diagonal dominance of matrix a will be destroyed, which causes possible mode coupling between Alfvén and slow modes.
The nonsingular eigenvalue problem presented in Eq. (21)is solved for each flux surface, using poloidally periodic boundary conditions for the eigenfunctions η and ζ. Typically, one solves the problem for a given axial mode number k, and we then get a set of discrete eigenvalues for each flux surface. When we quantify these eigenfrequencies for all flux surfaces, they together map out the continuous spectra. For completeness, we point out that if one uses a formulation based on the primitive variables ρ, v, p, and B, an additional entropy continuum kv_{z} is found. Goedbloed et al. (2004a) showed that this Eulerian entropy continuum, in ideal MHD, does not couple to any of the other continua, nor to the discrete stable and/or unstable modes.
3.2. Spectral properties and stability criterion
We now derive a general stability criterion for the continua. For this derivation, a Hilbert space with appropriate inner product needs to be constructed. The parallel gradient operator ℱ is a Hermitian operator under the inner product (27)for poloidally periodic functions v and w. Inspired by this property, a Hilbert space can be defined with the inner product (28)In this Hilbert space, the matrix operators a, b, and g are Hermitian operators. Setting v and w equal to the eigenfunctions (η,ζ)^{T}, one can derive a simple quadratic polynomial for the eigenfrequency from the spectral equation in Eq. (21)for the continua, namely (29)where the coefficients are Here, is the BruntVäisälä frequency projected onto a flux surface (32)The coefficient a is real and always greater than zero. As for a, the coefficient c is real but can be zero or even negative. The solutions of this second order polynomial in Eq. (29)are simply . If c is negative, a damped stable wave and an overstable mode are obtained. In this case, its absolute value is found from , while Re(ω) = kv_{z}, so the oscillation frequency is set by the local Doppler shift. By looking at the expression for c, we realize that this coefficient is equal to or greater than zero if (33)is satisfied everywhere in the plasma. Therefore, equilibria for which the entropy or the temperature (assuming γ > 1) is a flux function, will always be stable according to this criterion. However, when density is a flux function, the equilibrium variation may violate this criterion. This implies that such equilibria have unstable MHD continua, a situation identified in earlier, related work, as being caused by the convective continuum instability (CCI) (van der Holst et al. 2000b,a; Blokland et al. 2007a). This CCI instability criterion in Eq. (33)is analogous to the Schwarzschild criterion for convective instability, with one important difference. The Schwarzschild criterion deals with normal derivatives, while this criterion deals with tangential ones. This difference has also been found by various authors (Hellsten & Spies 1979; Hameiri 1983; Poedts et al. 1985; van der Holst et al. 2000a; Blokland et al. 2007a).
The relevance of this CCI instability to solar prominences may well be significant. We note that it is common practice to study the (slow) evolution of prominences carrying flux ropes by taking an equilibrium configuration, and subjecting it to e.g. photosphericboundarydriven changes, and relax it thereby repeatedly to successive equilibrium states (see e.g., Amari et al. 1999). As a result, the internal flux rope variation proceeds through a series of equilibria. For shortlived prominences, equilibria with density as a flux function may occur, and in this evolutionary process, the equilibrium can transit from the stable to the unstable regime as governed by the CCI criterion. As we later demonstrate, the CCI growth rate also varies with increasing gravity parameter g. As explained in Paper I, this gincrease could relate to an increase in size, internal density, or the weaker overall magnetic field strength of the prominence carrying flux rope.
3.3. Expansion in poloidal harmonics
We now continue to transform the two coupled differential equations for the MHD continua in Eq. (21)to an algebraic set, by expanding the eigenfunctions η and ζ in poloidal Fourier harmonics. This expansion has the advantage that we can perform a detailed investigation of how different poloidal harmonics may couple to each other. This coupling could either be between modes of the same kind, for example AlfvénAlfvén coupling, or between modes of different kind such as Alfvénslow coupling. In the next section, we show that in an actual solar prominence, more complex mode coupling is present.
Expanding the eigenfunctions η and ζ for the MHD continua in a finite Fourier series, (34)and applying the operator (35)to the spectral equation in Eq. (21), one obtains the algebraic coupled equations (36)The operator in Eq. (35)is the result of exploiting the inner product in Eq. (28)of the previously defined Hilbert space. The poloidal mode numbers m and m′ have the same range and the 2 × 2 submatrices are given by (37)The explicit expressions for the integrand functions K(ϑ), L(ϑ), M(ϑ), and N(ϑ) are These expressions are derived by exploiting the Hermitian property in Eq. (27)of the parallel gradient operator ℱ.
Once again, we consider the cylindrical limit in combination with zero gravity, which is appropriate for many earlier seismological studies of prominences. In this limit, all equilibrium quantities are independent of the poloidal angle ϑ, therefore the geodesic curvature will be zero. This means that the Alfvén and slow continuum decouple fully from each other. Furthermore, only one single harmonic is required for the expansion of the Fourier series of the eigenfunctions η and ζ. Solving the spectral equation for this case, one finds the wellknown local Dopplershifted Alfvén continua (42)and the local Dopplershifted slow continua (43)Looking at the expressions in Eqs. (38)–(41)for the submatrix elements, we can make some general statements about the MHD continua of solar prominences. In the absence of gravity (), the Alfvén and slow continuum can be coupled by the geodesic curvature present in the offdiagonal terms M(ϑ) and N(ϑ). With gravity, these offdiagonal terms can become more pronounced. The second statement is that an Alfvén continuum at a resonant surface, associated with the element K(ϑ), may shift away from the local Dopplershifted frequency kv_{z} because of the (1) geodesic curvature, involving its square ; (2) compressibility and finite pressure (γp); and (3) the pressure variation on a flux surface. All these effects suggest that an Alfvén gap can appear around the local Dopplershifted frequency and that Alfvénslow coupling should be taken into account for these gaps. This kind of coupling is investigated in the next section.
4. Small gravity expansion
In previous sections, we have derived equations for the continuous MHD spectrum of a translationally symmetric equilibrium. These derivations were completely general, and allowed us to identify an instability criterion on the one hand, as well as to identify the terms responsible for Alfvénslow type coupling. In this section, we quantify this coupling effect, using the small gravity parameter expansion for the equilibrium, similar to those discussed in Paper I. In the remainder of this paper, we assume a gravitational potential in which the prominence is embedded given by (44)where x_{0} is the location of the center of the last closed flux surface of the prominence and the role of gravity is represented by the constant g. This gravity parameter g denotes a dimensionless quantification of the relative importance of gravity on the equilibrium structure, and we explained in Paper I that this parameter can meaningfully vary by orders of magnitude for solar prominence configurations. For small g and assuming a circular crosssection of the prominence plasma, we quantified in Paper I the change induced in the equilibrium (the downward shift of the flux surfaces), and the corresponding straight field line coordinates. This kind of expansion is similar to the small inverse aspect ratio expansion for tokamak plasmas. Using this analytical result for the equilibrium, we can now demonstrate that gaps (or even instabilities) may appear in the continuous MHD spectrum because of mode coupling, as a direct consequence of the force of gravity.
4.1. Continuous MHD spectrum for equilibria weakly affected by gravity
We apply the small gravity expansion to the expressions found in Eqs. (38)–(41)for the matrix elements K(ϑ), M(ϑ), N(ϑ), and L(ϑ), respectively. For this expansion, we rely on the metric elements that are valid for straight field line coordinates exploiting the Shafranov shift Δ(r), given by Eq. (26) in our accompanying paper (Blokland & Keppens 2011). In addition, one can derive the quantities where the total magnetic field . The parameter t_{g} = 1 has been introduced to track the linear and quadratic dependence on the geodesic curvature. Furthermore, the parameter t_{e} allows us to distinguish the three options for the flux function dependence in the equilibrium, namely (51)These flux functions lead to different variants of the extended GradShafranov equation in Eq. (10) for the equilibrium, as presented in our accompanying paper (Blokland & Keppens 2011). We again note that the magnetically modified BruntVäisälä frequency from Eq. (50)can only be negative if the density is a flux function. This once more confirms that for this particular choice of the flux function, the MHD continua may become unstable as already found to be true in general above.
Applying the small gravity expansion, the eigenvalue equation in Eq. (36)can be written as (52)where We note that the first and second rows of the eigenvalue equation in Eq. (52)have been multiplied by and q^{1}, respectively, to get rid of the dimensional factors in front of the square brackets in the definitions given in Eqs. (38)–(41). The approximated expressions of the integrands of the matrix elements used in the spectral equation in Eq. (52)are now where Once more, the magnetically modified BruntVäisälä frequency (61) is only negative if the density is a flux function (t_{e} = 0).
4.2. Low frequency Δ m = 0 gap
The approximate relations just given can now be analyzed in terms of mode coupling possibilities. We focus on low frequencies, , and show that a Δm = 0 gap can be created in the neighborhood of a rational surface m + kq = 0, where q quantifies the fluxsurface specific local twist in the field lines on this surface. For such a surface, we determine the leading order coupling on an Alfvén continuum mode with the same poloidal mode number m. Close to the rational surface, we can safety neglect the third term of the expression in Eq. (55)for and Eq. (58)for . We are left with the eigenvalue equation for the MHD continua in Eq. (52), which reduces in essence to a fourmode coupling between the central Alfvén mode, and the three neighboring slow modes given by (62)where In this derivation, we made use of the finding that the integral of sin(ϑ) and sin^{2}(ϑ) over ϑ gives and , respectively. We note that for this fourmode coupling scheme, we neglected any coupling to the m ± 1 Alfvén modes, an assumption that is valid as long as the associated frequencies are far away in the spectrum. The same argument also holds for the coupling to the m ± 2 Alfvén and slow modes. By careful examination of this 4 × 4 matrix, one notices that the coupling between with either or is on the order of , while all the other types of coupling are on the order of . Because the main modes involved have the same poloidal mode number m, this coupling is known as a Δm = 0 gap possibility.
The matrix of the reduced eigenvalue problem presented in Eq. (62)is symmetric and therefore its eigenvalues will be real. This is consistent with the properties of the quadratic polynomial (29)for the eigenfrequency . If gravity drives a certain MHD continuum unstable, it means that the corresponding eigenfrequency will be smaller than zero, resulting in a complex eigenvalue . Once the eigenvalue is known, it is straightforward to determine the dominant character of the MHD mode by determining the dominant component of the eigenvector of the eigenvalue equation in Eq. (62). This analytic result will be used in what follows to verify the numerically obtained continua.
5. The spectral code PHOENIX
While the equilibria presented in Paper I were computed with the FINESSE code (Beliën et al. 2002), the stability analysis of the prominence models is subsequently done using the spectral code PHOENIX (Blokland et al. 2007b). Although this code was originally developed for analyzing axisymmetric plasma tori such as tokamaks and accretion disks, the design is flexible enough to analyze axially symmetric, gravitating plasma tubes. This possibility is exploited in the present PHOENIX version 1.2. For a given twodimensional equilibrium computed by FINESSE and converted to straight field line coordinates, we then need to select an axial wave number k and a range of poloidal wave numbers m for which we determine the complete MHD spectrum numerically. A range of poloidal wave numbers is required to allow for the anticipated mode coupling caused by the combinations of noncircular crosssections of the flux surfaces and/or the presence of a gravitational potential. Instead of solving the spectral equation in Eq. (12), PHOENIX solves the linearized version of the full set of nonlinear MHD equations. This allows us to extend the eigenspectrum computations for future incorporation of nonadiabatic effects, while it in addition will contain the trivial solution for the entropy modes known as the Eulerian entropy continuum. The discretization of these linearized MHD equations exploits a finite element method in the radial direction, while in the poloidal direction a spectral method is applied. The finite elements used are a combination of quadratic and cubic Hermite elements. This combination prevents the creation of spurious eigenvalues (Rappaz 1977). The resulting generalized nonHermitian eigenvalue problem is then solved using the JacobiDavidson method (Sleijpen & van der Vorst 1996). For such an eigenspectrum computation, one needs to supply boundary conditions at both the magnetic axis (set by symmetry) and the outer flux surface. There, we assume a perfectly conducting wall for simplicity. Strictly speaking, this condition is not applicable to true flux ropes carrying embedded solar prominences. However, since we focus on the continuous spectrum alone, the boundary conditions in the radial direction have hardly any influence, since those waves or instabilities are by definition localized on flux surfaces.
Exploiting this realization, Poedts & Schwarz (1993) introduced a more efficient method than the JacobiDavidson method to compute the continuous MHD spectrum. This method replaces the radial dependence of the cubic Hermite and the quadratic elements with log (ϵ) and 1/ϵ, respectively, on each individual flux surface. In this way, the singular behavior of the continuous modes can be approximated. Goedbloed (1975) and Pao (1975) describe in detail the reasons for this type of singular behavior in the continua for static plasma tori without gravity. Goedbloed et al. (2004b) pointed out that similarly singular behavior for the MHD continua prevails in plasma tori including toroidal plus poloidal flows and gravity. The same behavior is thus also expected for our translationally symmetric case. After replacing the elements, the eigenvalue problem becomes a small algebraic problem on each flux surface, with a size determined by the number of poloidal mode numbers taken along. For each flux surface, this can be solved using a direct QR method. All spectra below are computed in this way, for an axial wavenumber k = −1 and using seven poloidal mode numbers centered around m = 1. A larger number of poloidal Fourier harmonics can be used, which results in a more accurate solution but also in a longer computation time. The typical radial resolution used is either 101 or 1001 points, depending on how much detail is required in identifying each branch of the MHD continua.
6. Cool prominence surrounded by hot medium
The theory presented will be applied to the two equilibrium classes presented in Paper I, cool prominences embedded in hot medium and double layered prominences. We now discuss the first class. The equilibria and their associated profiles are discussed in our accompanying paper (Blokland & Keppens 2011). As described in our accompanying paper, we analyze the continuous MHD spectrum for increasing values of the gravity parameter g.
6.1. Gap creation when increasing g
Figure 1 shows part of the continuous MHD spectrum as a function of the radial flux coordinate computed by PHOENIX for a prominence equilibrium with gravity g = 0.001 and assuming that the temperature is a flux function. The frequencies are normalized to the Alfvén time , where μ_{0}, L, ρ_{mag}, and B_{mag} are the permeability, the radius of the prominence, the density on the magnetic axis, and the magnetic field on the magnetic axis, respectively. Near the magnetic axis where s = 0, each branch can be labeled immediately using the analytical expressions in Eqs. (42)and (43)for the MHD Alfvén and slow continua. At the location s ≈ 0.9, the safety factor q of the equilibrium becomes the value of one (see Paper I). At this location, the m = 1 Alfvén and the m = 1 slow continuum all become zero. At the low gravity value for Fig. 1, no true gap can be detected, and for progressively stronger gravity, we expect mode coupling to occur, ultimately inducing a possible gap or even instabilities in the continua. At the same location, the m = 0 slow and m = 2 slow continuum have the same frequency, namely approximately 0.797.
Fig. 1 Part of the continuous MHD spectrum for a cool solar prominence with gravity parameter g = 0.001 computed by PHOENIX. Shown are eigenfrequencies normalized to the Alfvén time τ_{A} as a function of the radial flux coordinate , all corresponding to flux surface localized eigenoscillations. 
We next take a closer look at the continuous MHD spectrum about the q = 1 surface, to compare the PHOENIX results with the ones predicted by the fourmode coupling scheme in Eq. (62). This comparison is presented in Fig. 2. Results from PHOENIX are represented by crosses, while the colored lines result from solving the eigenvalue equation in Eq. (62). The colors indicate the dominant characters, which are easily determined from the associated eigenvector. The colors red, green, blue, and yellow correspond to dominant (Alfvén), (slow), (slow), and (slow) characters, respectively. The plot shows that there is excellent agreement between both methods, as expected for small gravity. Furthermore, the m = 1 Alfvén continuum (red) intersects the m = 0 (blue) and m = 2 slow continuum (yellow). This intersection may lead to mode coupling, resulting in a Δm = 1 gap. The analysis of this kind of gap is left to future work. The crosses at zero frequency correspond to the Eulerian entropy continuum, as also found numerically by PHOENIX.
Fig. 2 The real part of the subspectrum of the MHD continua of a cool solar prominence equilibrium with small gravity (g = 0.001) as a function of the radial flux coordinate for axial wavenumber k = −1 and poloidal mode numbers m = {−2,...,4} . The temperature is assumed to be a flux function. The crosses are PHOENIX results, while the colored lines compare them with the fourmode coupling scheme. The colors red, green, blue, and yellow correspond to dominant , , , and components of the eigenvector, respectively. 
Fig. 3 The real part of the subspectrum of the MHD continua of a cool solar prominence equilibrium with intermediate gravity g = 0.500 as a function of the radial flux coordinate for axial wavenumber k = −1 and poloidal mode numbers m = {−2,...,4} . Temperature is assumed to be a flux function. Crosses are the PHOENIX results, while the colored lines give the fourmode coupling scheme. The colors red, green, blue, and yellow correspond to the dominant , , , and components of the eigenvector. 
Fig. 4 The real part of the subspectrum of the MHD continua of a cool solar prominence equilibrium with strong gravity (g = 1.000) as a function of the radial flux coordinate are shown for axial wavenumber k = −1 and poloidal mode numbers m = { − 2,...,4 } . Here, the temperature is assumed to be a flux function. The crosses are the PHOENIX results, while the colored lines are the results of the fourmode coupling scheme. The color red, green, blue, and yellow correspond to the dominant , , , and components of the eigenvector. 
We now repeat the same analysis, for a case with intermediate gravity g = 0.500 and still assuming that the temperature is a flux function. The results of PHOENIX results are in excellent agreement with the solutions of the fourmode coupling scheme in Eq. (62), as shown in Fig. 3. For this larger g value, the results now clearly illustrate that a gap appears in the m = 1 Alfvén continuum, while no gap appears in the m = 1 slow continuum. We note that the radial flux coordinate s of the q = 1 surface is shifted inwards relative to the previously discussed equilibrium. However, the radial location r of this surface remains almost the same (see Paper I).
If we continue to increase the gravity parameter g, for a cool prominence equilibrium where temperature is a flux function, we find progressive deviations between numerical and (approximate) analytical predictions. The plot presented in Fig. 4 shows the results from PHOENIX (crosses) and the fourmode coupling scheme (colored lines). Near q = 1, the results of the fourmode coupling scheme corresponding to the Alfvén continuum (red) show good agreement with PHOENIX, while for a slow continuum (green) there is a discrepancy. The fourmode coupling scheme indicates that a gap in the slow continuum may result, while PHOENIX does not find such a gap. The gap in the Alfvén continuum clearly becomes larger if gravity is increased. The radial location r of the q = 1 surface remains almost the same as in the two previously discussed cases, while its radial flux coordinate s is shifted outward.
Fig. 5 The frequency of the Alfvén continuum mode at the resonant q = 1 surface as a function of the gravity. The temperature is assumed to be a flux function. 
In the previous three equilibria, which only differ in terms of the increased role of gravity, we observed that the gap appearing in the Alfvén continuum becomes larger if the gravity is increased. This dependence was investigated by determining the Alfvén frequency at the q = 1 surface, for various values of the gravity parameter. The result of this parametric investigation (involving the construction of numerical equilibria, subsequently diagnosed for their MHD continua by PHOENIX) has been plotted in Fig. 5. The plot shows excellent agreement between PHOENIX and the growth of the “gap” predicted by the fourmode coupling scheme. It is seen that the Alfvén gap frequency depends nonlinearly on the gravity parameter g. Despite this excellent agreement for the Alfvén continuum gap up to g = 1, successive differences occur, e.g. as noted for the slow continua shown in Fig. 4. In Fig. 5, the horizontal axis does not use the gravity parameter g itself, but rather a scaled value. The equilibrium with gravity g = 0.001, g = 0.500, and g = 1.000 correspond to rρ_{0}g/p_{0} = 0.0003, rρ_{0}g/p_{0} = 0.154, and rρ_{0}g/p_{0} = 0.300, respectively.
6.2. Unstable continua
Fig. 6 The real part of the subspectrum of the MHD continua of a cool solar prominence equilibrium with intermediate gravity (g = 0.500) as a function of the radial flux coordinate are shown for axial wavenumber k = −1 and poloidal mode numbers m = { − 2,...,4 } . Density is assumed to be a flux function. The crosses are the PHOENIX results, while the colored lines are the results of the fourmode coupling scheme. The color red, green, blue, and yellow correspond to the dominant , , , and component of the eigenvector. 
Fig. 7 The imaginary part of the subspectrum of the MHD continua of a cool solar prominence equilibrium with intermediate gravity (g = 0.500) as a function of the radial flux coordinate are shown for axial wavenumber k = −1 and poloidal mode numbers m = {−2,...,4}. These quantify the growth rates of the CCI instabilities. Density is assumed to be a flux function. The crosses are the PHOENIX results, while the colored lines are the results of the fourmode coupling scheme. The color red, green, blue, and yellow correspond to the dominant , , , and component of the eigenvector. 
Finally, to illustrate the effect of the chosen flux function in the equilibrium on the MHD continua, we now discuss in the intermediate gravity case g = 0.500, the continuous spectra when the density or the entropy is a flux function. We recall that when the density is a flux function the continuous MHD spectrum may become unstable. The continuous MHD spectrum of the density case is presented in Fig. 6. The figure shows there is again an excellent agreement between the fourmode coupling scheme and PHOENIX. The plot shows two distinct features, that are not present when the temperature is a flux function. Firstly, no gap appears in the m = 1 Alfvén continuum. Near the q = 1 surface (s ≈ 0.90), an interesting observation can be made, where Alfvén and slow mode characters mingle. A mode at for example s = 0.89 with dominant m = 1 Alfvén character slowly changes to a dominant m = 1 slow character when it is traced closer to the q = 1 surface. The transition points are indicated by the arrows in the plot. The second distinct feature is the merger of the two m = 1 slow branches to one at s ≈ 0.852 and s ≈ 0.94. Exactly at these locations, the m = 1 slow continuum turns unstable as shown in Fig. 7. This figure quantifies the imaginary parts of the continuum eigenfrequencies, corresponding to (normalized) growth rates. This instability is the convective continuum instability (CCI) as discussed by Blokland et al. (2007a). There is again excellent agreement between PHOENIX and the fourmode coupling scheme. Furthermore, the transition points of the dominant character of the mode indicated by arrows in Fig. 7 are exactly at the same locations as those indicated in Fig. 6. From this analysis, we conclude that the CCI at the q = 1 surface has a dominant Alfvén character.
Fig. 8 The growth rate of the CCI mode at the resonant q = 1 surface as a function of the gravity parameter. The density is assumed to be a flux function. 
As in the case where temperature is a flux function, we can parametrically investigate the dependence of the growth rate of the CCI on the gravity parameter g. This analysis again involves repeated FINESSE computations for increasing values of g, all diagnosed with PHOENIX, and the extraction of the growth rate for the most unstable mode in each case. These results are presented in Fig. 8. Once more, there is excellent agreement between the two methods. Furthermore, our analysis clearly identities a linear dependence between the growth rate and the gravity parameter g. Owing to our nondimensionalization, this implies that there is a larger growth of the CCI mode for flux ropes with embedded prominences that (1) gradually grow in crosssectional extent; and/or (2) become more heavy (density increase); and/or (3) decrease in magnetic field strength.
In the context of prominence dynamics, the existence of equilibria viable to unstable continua may have relevance to the complex dynamics observed in high resolution prominence movies (Berger et al. 2008), as well as to the sudden disappearance of prominence plasma, the socalled disparitionbrusque phenomenon (see e.g., TandbergHanssen 2011, and references therein). It should be emphasized that the onset of the CCI instability will be localized on those flux surfaces where the CCI criterion is violated. As shown in Fig. 7, growth rates can vary from one flux surface to another in a given unstable case, while in addition, the maximal growth rate also varies because of gradual equilibrium changes, as demonstrated in Fig. 8. Whether this linear prediction for the onset of internal, continuumdriven dynamics leads to the rather turbulent prominence dynamics revealed by contemporary observations, or even to full disruption and the disappearance of prominence structures, must be answered by followup high resolution, nonlinear numerical simulations. This requires initializing multidimensional nonlinear simulations with the spectrally diagnosed equilibria.
6.3. Spectra for equilibria in which entropy is a flux function
We now consider the case where entropy is a flux function, for an equilibrium with intermediate gravity (g = 0.500). The resulting continuous MHD spectrum is presented in Fig. 9. There is again excellent agreement between both methods. It can be clearly seen that no gap in the Alfvén continuum appears. The continua for the strong gravity g = 1.000 case have also been computed, which also show no appearance of a gap. This can be understood by examining Eq. (29)for the eigenfrequency. When the entropy is a flux function, the BruntVäisälä frequency projected on a flux surface N_{BV,pol} is zero, thereby transforming the eigenfrequency given by Eq. (29), to one similar to a case without gravity, for which no gap can appear.
Fig. 9 The real part of the subspectrum of the MHD continua of a cool solar prominence equilibrium with intermediate gravity (g = 0.500) as a function of the radial flux coordinate are shown for axial wavenumber k = −1 and poloidal mode numbers m = { − 2,...,4 } . Here, entropy is assumed to be a flux function. The crosses are the PHOENIX results, while the colored lines are the results of the fourmode coupling scheme. The color red, green, blue, and yellow correspond to the dominant , , , and component of the eigenvector. 
7. Double layered prominences
We discuss the continuous MHD spectrum of the second equilibrium class, the double layered prominences. As for cool prominences, the equilibria and the associated profiles are discussed in our accompanying paper (Blokland & Keppens 2011). For this equilibrium class, we also investigate the influence of the gravity on the continuous MHD spectrum by slowly increasing the gravity.
7.1. Continua for increasing gravity
The first equilibrium that was analyzed with PHOENIX is the one where the temperature is assumed to be a flux function and the gravity is low at g = 0.001. Part of the MHD spectrum is shown in Fig. 10. The individual branches near the magnetic axis can be labeled using the analytical expressions for the Alfvén continuum in Eq. (42)and slow continuum in Eq. (43). The associated radial flux coordinate s of the q = 1 surface is approximately 0.927. At exactly this location, the m = 1 Alfvén and m = 1 slow continuum are zero, as expected, and the m = 0 and m = 2 slow continuum have the same frequency.
Fig. 10 A part of the continuous MHD spectrum of double layered solar prominence with gravity g = 0.001 computed by PHOENIX as a function of the radial flux coordinate . 
Fig. 11 The real part of the subspectrum of the MHD continua of a double layered solar prominence equilibrium with small gravity (g = 0.001) as a function of the radial flux coordinate are shown for axial wavenumber k = −1 and poloidal mode numbers m = { − 2,...,4 } . Here, the temperature is assumed to be a flux function. The crosses are the PHOENIX results, while the colored lines are the results of the fourmode coupling scheme. The color red, green, blue, and yellow correspond to the dominant , , , and component of the eigenvector. 
The plot shown in Fig. 11 contains the PHOENIX results (crosses) together with the solutions of the fourmode coupling scheme represented by Eq. (62)(colored lines), and illustrates the excellent agreement between the two methods. As for the cool prominences, the m = 1 Alfvén continuum (red) intersects the m = 0 (blue) and m = 2 slow continuum (yellow). This may in turn lead to mode coupling, creating a high frequency Δm = 1 gap. The analysis of these Δm = 1 gaps was not performed here, but left to future work. We note that the crosses with zero frequency correspond to the Eulerian entropy continuum.
Fig. 12 The real part of the subspectrum of the MHD continua of a double layered solar prominence equilibrium with strong gravity (g = 0.500) as a function of the radial flux coordinate are shown for axial wavenumber k = −1 and poloidal mode numbers m = { − 2,...,4 } . Here, the temperature is assumed to be a flux function. The crosses are the PHOENIX results, while the colored lines are the results of the fourmode coupling scheme. The color red, green, blue, and yellow correspond to the dominant , , , and component of the eigenvector. 
Figure 12 contains the plot for an equilibrium with intermediate gravity g = 0.500. The temperature is again assumed to be a flux function. Near the q = 1 surface (s ≈ 0.927), the results of PHOENIX and the fourmode coupling scheme show excellent agreement for the m = 1 Alfvén and slow continuum. Owing to the nonzero gravity, a gap has been created in the m = 1 Alfvén continuum. However, m = 0 and m = 2 slow continua of the coupling scheme slightly mismatch the PHOENIX results. This mismatch is because by our neglect of the coupling with m = −1 and m = 3 slow continua. In the plot, these continua are represented by the crosses with nonzero frequency in combination without a solid line. At four different radial flux coordinates, s ≈ 0.82, s ≈ 0.837, s ≈ 0.973, and s ≈ 0.978 a Δm = 1 gap was produced (note the color change of the solid lines about these radii). As expected, the fourmode coupling scheme only partly captures this kind of coupling. Therefore, about these Δm = 1 gaps there is a mismatch between the results of both methods, and a separate detailed analysis of the Δm = 1 coupling schemes appears to be required.
Fig. 13 The real part of the subspectrum of the MHD continua of a double layered solar prominence equilibrium with strong gravity (g = 1.000) as a function of the radial flux coordinate are shown for axial wavenumber k = −1 and poloidal mode numbers m = {−2,...,4}. Here, the temperature is assumed to be a flux function. The crosses are the PHOENIX results, while the colored lines are the results of the fourmode coupling scheme. The color red, green, blue, and yellow correspond to the dominant , , , and component of the eigenvector. 
The last equilibrium of this subsection that we discuss is the one with strong gravity g = 1.000. The temperature is assumed to be a flux function. For this equilibrium, one observes the onset of the double layered structure in the pressure and density as demonstrated in our accompanying paper (Blokland & Keppens 2011). The results for the m = 1 Alfvén and slow continuum show excellent agreement near the q = 1 surface (s ≈ 0.927). Owing to the stronger gravity compared to the previous case, the gap in m = 1 Alfvén continuum is larger. Near this surface, the previously mentioned mismatch of both the m = 0 and m = 2 slow continuum becomes larger because of the neglect of the coupling with both the m = −1 and m = 3 slow continuum. This coupling becomes stronger if the gravity is increased. The fourmode coupling scheme is also unable to capture the coupling at the radial flux coordinates s ∈ [0.828,0.833] and s ∈ [0.975,0.978]. This is again because we neglected the coupling with both the m = −1 and m = 3 slow continuum. Therefore, close to these radii there is a large discrepancy between PHOENIX and the coupling scheme in Eq. (62). Furthermore, we note the coupling between these two slow continua at the q = 1 surface. This kind of coupling creates a Δm = 2 gap that was neglected in our theoretical analysis.
Fig. 14 The frequency of the Alfvén continuum mode at the resonant q = 1 surface as a function of the gravity. The temperature is assumed to be a flux function. 
As for the equilibrium class of cool prominences, we observe that when the gravity is increased, the gap in the m = 1 Alfvén continuum increases. This gravity dependence of the gap was investigated in detail by computing the Alfvén frequency at the q = 1 surface for different values of the gravity. The results of this investigation are presented in Fig. 14. The results of PHOENIX and the fourmode coupling scheme show excellent agreement with each other. Furthermore, we note that the frequency varies almost linearly with the gravity parameter g. This linear dependence is very different from the clear nonlinear dependence found for cool prominences. This excellent agreement is again only valid for the Alfvén continuum. As demonstrated in Figs. 12 and 13, this does not automatically mean that it also holds for the slow continua. The dimensionless gravity parameter g = 0.001, g = 0.500, and g = 1.000 correspond to rρ_{0}g/p_{0} = 0.0009, rρ_{0}g/p_{0} = 0.4597, and rρ_{0}g/p_{0} = 0.9192, respectively.
7.2. Continua for strong gravity, with axial flow
We now discuss the last equilibrium of the second class. For this equilibrium, we assume that the temperature is a flux function, we adopt an extra strong gravity parameter value (g = 5.000), and we include an axial, shear flow. The derivation in our accompanying paper (Blokland & Keppens 2011) shows that the axial flow must be a flux function. The flow itself is described by the profile (69)where the constant A_{4} = 0.11. This corresponds to a maximum Alfvén Mach number of 0.105. Using the typical Alfvén speed of 218 km s^{1}, as mentioned in our accompanying paper (Blokland & Keppens 2011), this translates into a filamentaligned flow of at most 22.9 km s^{1}. Part of the resulting continuous MHD spectrum computed by PHOENIX is shown in Fig. 15. Instead of using 7 poloidal harmonics, we used 15 harmonics to obtain an accurate result (i.e., one that no longer changed when adding even more poloidal harmonics). This is a clear indication that the strong gravity causes strong coupling between different continua. Owing the extra strong gravity, there is no meaningful comparison between the PHOENIX results and those predicted from our theoretical analysis, which is only appropriate if the gravity is weak. In the figure, we indicate the nonzero Eulerian entropy continuum (where ω = kv_{z} labeled E) that breaks the updown symmetry about the zero frequency. There is instead now symmetry about this Eulerian entropy continuum. Furthermore, PHOENIX finds modes with zero frequency, which correspond to the solution of the generalized eigenvalue equation in combination with the perturbed magnetic field being divergence free. If one carefully examines the plot, one notices that one continuum branch intersects the Eulerian entropy continuum twice, namely at s ≈ 0.808 and at s ≈ 0.917. At these locations exactly, the safety factor q is one. Figure 14 of our accompanying paper (Blokland & Keppens 2011) clearly shows the q = 1 surface locations, as seen in a plot of safety factor q across the vertical direction. Using the radial flux coordinate s as a flux surface label where we go from the magnetic axis to the flux rope edge, there will be two q = 1 surfaces. We note that the clear coupling effects at other radial locations as well, creating gaps that may define preferred frequency ranges for global eigenoscillations.
Fig. 15 A part of the continuous MHD spectrum of a double layered solar prominence with strong gravity g = 5.000 computed by PHOENIX as a function of the radial flux coordinate . 
8. Conclusions
We presented the equations describing the complete MHD spectrum. These equations have been specialized for the continuous MHD spectrum using techniques frequently used in research on fusion plasma oscillations. From the equations for the continuous spectrum, a stability criterion was derived, showing that prominence equilibria for which the temperature or the entropy is a flux function have stable continuum modes. However, when the density is a flux function, the continua may turn unstable resulting in the convective continuum instability. Whether this instability is related to the intricate dynamics witnessed in actual prominence plasmas, or explains why certain filaments disappear abruptly, remains to be investigated by followup nonlinear computations. We have demonstrated how the onset of CCI modes in flux ropes, which carry embedded prominences, will manifest itself as truly fluxsurface localized fluctuations. Their growth rate increases when a flux rope expands (the gravity parameter g increases because of both rope widening and the associated magnetic field weakening) or when it becomes more massloaded. We note that this route to potentially violent prominence dynamics is here related to internal prominence properties. This contrasts with studies (see e.g., MacTaggert & Hood 2009) that emphasize the role of external influences, studies that emphasize the response to changes in the overlying field, or the associated pressure buildup in the surroundings by the shearing of its footpoints. We note that the nonlinear simulations that may ultimately decide between these competing or possibly cooperating aspects, will need to go beyond the zero beta, and/or no gravity assumption typically employed in these studies (Aulanier et al. 2002; Amari et al. 1999). Moreover, to capture the details of fluxsurface localized continuum modes with varying MHD character (slow versus Alfvén) from one flux surface to another, the requirements of the 3D resolution will likely only be achievable by using effectively high resolution, gridadaptive simulations. In addition, guidance on the actual magnetic field structure in prominence surroundings will be needed (such as that obtained in Schmieder et al. 2007), to distinguish between the theoretical freedom in flux profiles.
The equations for the MHD continuous spectrum have been specialized for small gravity, using an expansion in poloidal Fourier harmonics. The resulting equations have been applied to a low frequency Δm = 0 gap near a resonant q surface. Near this surface, a fourmode coupling scheme has been identified. This scheme couples one Alfvén mode and three slow modes, and can cause gaps (avoided crossings) between continuum branches as these are traced across all nested flux surfaces. These gaps appear at low frequency, and form the ideal locations for global discrete eigenoscillations. It remains to be shown which kind of coupling occurs at finite frequency, Δm = 1 gaps, and to determine whether discrete global modes exist within these frequency ranges, under realistic prominence conditions. These global modes, as well as their coupling to flux localized surface modes (with resonant damping), need further study, to address whether previously established prominence seismology results are consistent with these more realistic prominences.
The MHD spectral code PHOENIX has been used to accurately compute the continuous MHD spectrum of prominence equilibria. Two classes of equilibria have been considered, cool prominences embedded in hot medium and prominences with a double layered pressure and density structure for different values of the dimensionless gravity parameter as discussed in our accompanying paper (Blokland & Keppens 2011). The results of PHOENIX are compared with analytic results of the fourmode coupling scheme, and we have found excellent agreement for small gravity (as expected), while for strong gravity, we have identified progressive discrepancies, emphasizing the need for full numerical diagnosis. For cool prominences, this discrepancy manifested itself in the slow continuum with the same poloidal mode number as the Alfvén continuum, while for double layered prominences, differences occurred in the slow continuum with a poloidal mode number ± 1 with respect to the Alfvénic one. For both classes, when temperature is a flux function, the results demonstrate that a gap in the Alfvén continuum can appear due to the gravity. For cool prominences, the Alfvén frequency associated with the gap depends nonlinearly on the gravity parameter, while for double layered prominences it scales linearly with it.
For cool prominences for which the density or the entropy is a flux function, the continuous MHD spectrum has also been computed. In the case of the density, the continuous spectrum becomes unstable leading to convective continuum instability. Our analysis reveals that close to the q = 1 surface, the CCI mode has a dominant Alfvén character, which changes away from this surface to a dominant slow type. Furthermore, no gap in the Alfvén continuum appears, in contrast to the case when the temperature is a flux function. For the entropy case, the continuous spectrum is stable and again no gap appears in the Alfvén continuum. In the context of prominence configurations, longlived equilibria almost necessitate the temperature flux function dependence, but for active region prominences, the findings about CCI possibilities may provide important information about the onset of turbulent dynamics.
Acknowledgments
This work was carried out within the framework of the European Fusion Programme, supported by the European Communities under contract of the Association EURATOM/FOM. Views and opinions expressed herein do not necessarily reflect those of the European Commission. R.K. acknowledges financial support by project GOA/2009/009 (K.U. Leuven). The research leading to these results has received funding from the European Commission’s Seventh Framework Programme (FP7/20072013) under the grant agreement SWIFF (project No. 263340, http://www.swiff.eu).
References
 Amari, T., Luciani, J., Mikic, Z., & Linker, J. 1999, ApJ, 518, L57 [NASA ADS] [CrossRef] [Google Scholar]
 Arregui, I., & Ballester, J. 2011, Space Sci. Rev., in press [Google Scholar]
 Aulanier, G., De Vore, C., & Antiochos, S. 2002, ApJ, 567, L97 [NASA ADS] [CrossRef] [Google Scholar]
 Beliën, A., Poedts, S., & Goedbloed, J. 1997, A&A, 332, 995 [NASA ADS] [Google Scholar]
 Beliën, A., Botchev, M., Goedbloed, J., van der Holst, B., & Keppens, R. 2002, J. Comp. Phys., 182, 91 [NASA ADS] [CrossRef] [Google Scholar]
 Berger, T., Shine, R., Slater, G., et al. 2008, ApJ, 676, L89 [Google Scholar]
 Bernstein, I., Frieman, E., Kruskal, M., & Kulsrud, R. 1958, Proc. Roy. Soc. of London, A, 17 [Google Scholar]
 Blokland, J., & Keppens, R. 2011, A&A, 532, A93 (Paper I) [Google Scholar]
 Blokland, J., Keppens, R., & Goedbloed, J. 2007a, A&A, 467, 21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Blokland, J., van der Holst, B., Keppens, R., & Goedbloed, J. 2007b, J. Comp. Phys., 226, 509 [NASA ADS] [CrossRef] [Google Scholar]
 Bondeson, A., Iacono, R., & Bhattacharjee, A. 1987, Phys. Fluids, 30, 2167 [NASA ADS] [CrossRef] [Google Scholar]
 Chance, M., Green, J., Grimm, R., & Johnson, J. 1977, Nucl. Fusion, 71, 65 [NASA ADS] [CrossRef] [Google Scholar]
 DeRosa, M., Schrijver, C., Barnes, G., et al. 2009, ApJ, 696, 1780 [Google Scholar]
 Frieman, E., & Rotenberg, M. 1960, Rev. Mod. Phys., 32, 898 [NASA ADS] [CrossRef] [Google Scholar]
 Goedbloed, J. 1971, Physica, 53, 501 [NASA ADS] [CrossRef] [Google Scholar]
 Goedbloed, J. 1975, Phys. Fluids, 18, 1258 [NASA ADS] [CrossRef] [Google Scholar]
 Goedbloed, J. 1997, Transactions of Fusion Technology, 33, 115 [Google Scholar]
 Goedbloed, J., & Poedts, S. 2004, Principles of magnetohydrodynamics (Cambridge University Press) [Google Scholar]
 Goedbloed, J., Huysmans, G., Holties, H., Kerner, W., & Poedts, S. 1993, Plasma Phys. Control. Fusion, 35, B227 [NASA ADS] [CrossRef] [Google Scholar]
 Goedbloed, J., Beliën, A., van der Holst, B., & Keppens, R. 2004a, Phys. Plasmas, 11, 4332 [NASA ADS] [CrossRef] [Google Scholar]
 Goedbloed, J., Beliën, A., van der Holst, B., & Keppens, R. 2004b, Phys. Plasmas, 11, 28 [NASA ADS] [CrossRef] [Google Scholar]
 Goedbloed, J., Keppens, R., & Poedts, S. 2010, Advanced magnetohydrodynamics (Cambridge University Press) [Google Scholar]
 Goossens, M., Erdélyi, R., & Ruderman, M. 2011, Space Sci. Rev., in press [Google Scholar]
 Hain, K., & Lüst, R. 1958, Z. Naturforsch., 13a, 936 [Google Scholar]
 Hameiri, E. 1981, J. Math. Phys., 22, 2080 [NASA ADS] [CrossRef] [Google Scholar]
 Hameiri, E. 1983, Phys. Fluids, 26, 230 [NASA ADS] [CrossRef] [Google Scholar]
 Hellsten, T., & Spies, G. 1979, Phys. Fluids, 22, 743 [NASA ADS] [CrossRef] [Google Scholar]
 Keppens, R., van der Linden, R., & Goossens, M. 1993, Sol. Phys., 144, 267 [NASA ADS] [CrossRef] [Google Scholar]
 Keppens, R., Casse, F., & Goedbloed, J. 2002, ApJ, 569, L121 [NASA ADS] [CrossRef] [Google Scholar]
 Lin, H., Penn, M., & Kuhn, J. 1998, ApJ, 493, 978 [NASA ADS] [CrossRef] [Google Scholar]
 Labrosse, N., Heinzel, P., Vial, J.C., et al. 2010, Space Sci. Rev., 151, 243 [Google Scholar]
 Lin, Y., Engvold, O., Rouppe van der Voort, L., Wiik, J., & Berger, T. 2005, Sol. Phys., 226, 239 [NASA ADS] [CrossRef] [Google Scholar]
 Luna, M., Terradas, J., Oliver, R., & Ballester, J. 2010, ApJ, 716, 1371 [NASA ADS] [CrossRef] [Google Scholar]
 Mackay, D., Karpen, J., Ballester, J., Schmieder, B., & Aulanier, G. 2010, Space Sci. Rev., 151, 333 [NASA ADS] [CrossRef] [Google Scholar]
 MacTaggert, D., & Hood, A. 2009, A&A, 508, 445 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pao, Y. 1975, Nucl. Fusion, 15, 631 [CrossRef] [Google Scholar]
 Poedts, S., & Goossens, M. 1991, Sol. Phys., 133, 281 [NASA ADS] [CrossRef] [Google Scholar]
 Poedts, S., & Schwarz, E. 1993, J. Comp. Phys., 105, 165 [NASA ADS] [CrossRef] [Google Scholar]
 Poedts, S., Hermans, D., & Goossens, M. 1985, A&A, 151, 16 [NASA ADS] [Google Scholar]
 Rappaz, J. 1977, Numer. Math., 28, 15 [NASA ADS] [CrossRef] [Google Scholar]
 Schmieder, B., Gunár, S., Heinzel, P., & Anzer, U. 2007, Sol. Phys., 241, 53 [NASA ADS] [CrossRef] [Google Scholar]
 Sleijpen, G., & van der Vorst, H. 1996, SIAM J. Matrix Anal. Appl., 17, 401 [CrossRef] [Google Scholar]
 Soler, R., Oliver, R., & Ballester, J. 2008, ApJ, 684, 725 [NASA ADS] [CrossRef] [Google Scholar]
 Soler, R., Arregui, I., Oliver, R., & Ballester, J. 2010, ApJ, 722, 1778 [NASA ADS] [CrossRef] [Google Scholar]
 TandbergHanssen, E. 2011, Sol. Phys., 269, 237 [NASA ADS] [CrossRef] [Google Scholar]
 Terradas, J., MolownyHoras, R., Wiehr, E., et al. 2002, A&A, 393, 637 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Trujillo Bueno, J., Landi Degl’Innocenti, E., Collados, M., Merenda, L., & Manso Sainz, R. 2002, Nature, 415, 403 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 van der Holst, B., Beliën, A., & Goedbloed, J. 2000a, Phys. Plasmas, 7, 4208 [NASA ADS] [CrossRef] [Google Scholar]
 van der Holst, B., Beliën, A., & Goedbloed, J. 2000b, Phys. Rev. Lett., 84, 2865 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
All Figures
Fig. 1 Part of the continuous MHD spectrum for a cool solar prominence with gravity parameter g = 0.001 computed by PHOENIX. Shown are eigenfrequencies normalized to the Alfvén time τ_{A} as a function of the radial flux coordinate , all corresponding to flux surface localized eigenoscillations. 

In the text 
Fig. 2 The real part of the subspectrum of the MHD continua of a cool solar prominence equilibrium with small gravity (g = 0.001) as a function of the radial flux coordinate for axial wavenumber k = −1 and poloidal mode numbers m = {−2,...,4} . The temperature is assumed to be a flux function. The crosses are PHOENIX results, while the colored lines compare them with the fourmode coupling scheme. The colors red, green, blue, and yellow correspond to dominant , , , and components of the eigenvector, respectively. 

In the text 
Fig. 3 The real part of the subspectrum of the MHD continua of a cool solar prominence equilibrium with intermediate gravity g = 0.500 as a function of the radial flux coordinate for axial wavenumber k = −1 and poloidal mode numbers m = {−2,...,4} . Temperature is assumed to be a flux function. Crosses are the PHOENIX results, while the colored lines give the fourmode coupling scheme. The colors red, green, blue, and yellow correspond to the dominant , , , and components of the eigenvector. 

In the text 
Fig. 4 The real part of the subspectrum of the MHD continua of a cool solar prominence equilibrium with strong gravity (g = 1.000) as a function of the radial flux coordinate are shown for axial wavenumber k = −1 and poloidal mode numbers m = { − 2,...,4 } . Here, the temperature is assumed to be a flux function. The crosses are the PHOENIX results, while the colored lines are the results of the fourmode coupling scheme. The color red, green, blue, and yellow correspond to the dominant , , , and components of the eigenvector. 

In the text 
Fig. 5 The frequency of the Alfvén continuum mode at the resonant q = 1 surface as a function of the gravity. The temperature is assumed to be a flux function. 

In the text 
Fig. 6 The real part of the subspectrum of the MHD continua of a cool solar prominence equilibrium with intermediate gravity (g = 0.500) as a function of the radial flux coordinate are shown for axial wavenumber k = −1 and poloidal mode numbers m = { − 2,...,4 } . Density is assumed to be a flux function. The crosses are the PHOENIX results, while the colored lines are the results of the fourmode coupling scheme. The color red, green, blue, and yellow correspond to the dominant , , , and component of the eigenvector. 

In the text 
Fig. 7 The imaginary part of the subspectrum of the MHD continua of a cool solar prominence equilibrium with intermediate gravity (g = 0.500) as a function of the radial flux coordinate are shown for axial wavenumber k = −1 and poloidal mode numbers m = {−2,...,4}. These quantify the growth rates of the CCI instabilities. Density is assumed to be a flux function. The crosses are the PHOENIX results, while the colored lines are the results of the fourmode coupling scheme. The color red, green, blue, and yellow correspond to the dominant , , , and component of the eigenvector. 

In the text 
Fig. 8 The growth rate of the CCI mode at the resonant q = 1 surface as a function of the gravity parameter. The density is assumed to be a flux function. 

In the text 
Fig. 9 The real part of the subspectrum of the MHD continua of a cool solar prominence equilibrium with intermediate gravity (g = 0.500) as a function of the radial flux coordinate are shown for axial wavenumber k = −1 and poloidal mode numbers m = { − 2,...,4 } . Here, entropy is assumed to be a flux function. The crosses are the PHOENIX results, while the colored lines are the results of the fourmode coupling scheme. The color red, green, blue, and yellow correspond to the dominant , , , and component of the eigenvector. 

In the text 
Fig. 10 A part of the continuous MHD spectrum of double layered solar prominence with gravity g = 0.001 computed by PHOENIX as a function of the radial flux coordinate . 

In the text 
Fig. 11 The real part of the subspectrum of the MHD continua of a double layered solar prominence equilibrium with small gravity (g = 0.001) as a function of the radial flux coordinate are shown for axial wavenumber k = −1 and poloidal mode numbers m = { − 2,...,4 } . Here, the temperature is assumed to be a flux function. The crosses are the PHOENIX results, while the colored lines are the results of the fourmode coupling scheme. The color red, green, blue, and yellow correspond to the dominant , , , and component of the eigenvector. 

In the text 
Fig. 12 The real part of the subspectrum of the MHD continua of a double layered solar prominence equilibrium with strong gravity (g = 0.500) as a function of the radial flux coordinate are shown for axial wavenumber k = −1 and poloidal mode numbers m = { − 2,...,4 } . Here, the temperature is assumed to be a flux function. The crosses are the PHOENIX results, while the colored lines are the results of the fourmode coupling scheme. The color red, green, blue, and yellow correspond to the dominant , , , and component of the eigenvector. 

In the text 
Fig. 13 The real part of the subspectrum of the MHD continua of a double layered solar prominence equilibrium with strong gravity (g = 1.000) as a function of the radial flux coordinate are shown for axial wavenumber k = −1 and poloidal mode numbers m = {−2,...,4}. Here, the temperature is assumed to be a flux function. The crosses are the PHOENIX results, while the colored lines are the results of the fourmode coupling scheme. The color red, green, blue, and yellow correspond to the dominant , , , and component of the eigenvector. 

In the text 
Fig. 14 The frequency of the Alfvén continuum mode at the resonant q = 1 surface as a function of the gravity. The temperature is assumed to be a flux function. 

In the text 
Fig. 15 A part of the continuous MHD spectrum of a double layered solar prominence with strong gravity g = 5.000 computed by PHOENIX as a function of the radial flux coordinate . 

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.