A. J. Diaz

Rayleigh-Taylor instability in prominences from numerical simulations including partial ionization effects

We study the Rayleigh-Taylor instability (RTI) at a prominence-corona transition region in a non-linear regime. Our aim is to understand how the presence of neutral atoms in the prominence plasma influences the instability growth rate, and the evolution of velocity, magnetic field vector and thermodynamic parameters of turbulent drops. We perform 2.5D numerical simulations of the instability initiated by a multi-mode perturbation at the corona-prominence interface using a single-fluid MHD approach including a generalized Ohms law. The initial equilibrium configuration is purely hydrostatic and contains a homogeneous horizontal magnetic field forming an angle with the direction in which the plasma is perturbed. We analyze simulations with two different orientations of the magnetic field. For each field orientation we compare two simulations, one for the pure MHD case, and one including the ambipolar diffusion in the Ohms law (AD case). Other than that, both simulations for each field orientation are identical. The numerical results in the initial stage of the instability are compared with the analytical linear calculations. We find that the configuration is always unstable in the AD case. The growth rate of the small-scale modes in the non-linear regime is up to 50% larger in the AD case than in the purely MHD case and the average velocities of flows are a few percent larger. Significant drift momenta are found at the interface between the coronal and the prominence material at all stages of the instability, produced by the faster downward motion of the neutral component with respect to the ionized component. The differences in temperature of the bubbles between the ideal and non-ideal case are also significant, reaching 30%. There is an asymmetry between large rising bubbles and small-scale down flowing fingers, favoring the detection of upward velocities in observations.

KELVIN-HELMHOLTZ INSTABILITY IN PARTIALLY IONIZED COMPRESSIBLE PLASMAS

The Kelvin-Helmholtz instability (KHI) has been observed in the solar atmosphere. Ion-neutral collisions may play a relevant role for the growth rate and evolution of the KHI in solar partially ionized plasmas such as in, e.g., solar prominences. Here, we investigate the linear phase of the KHI at an interface between two partially ionized magnetized plasmas in the presence of a shear flow. The effects of ion-neutral collisions and compressibility are included in the analysis. We obtain the dispersion relation of the linear modes and perform parametric studies of the unstable solutions. We find that, in the incompressible case, the KHI is present for any velocity shear regardless of the value of the collision frequency. In the compressible case, the domain of instability depends strongly on the plasma parameters, especially the collision frequency and the density contrast. For high collision frequencies and low density contrasts the KHI is present for super-Alfvenic velocity shear only. For high density contrasts the threshold velocity shear can be reduced to sub-Alfvenic values. For the particular case of turbulent plumes in prominences, we conclude that sub-Alfvenic flow velocities can trigger the KHI thanks to the ion-neutral coupling.

THE EFFECTS OF MAGNETIC-FIELD GEOMETRY ON LONGITUDINAL OSCILLATIONS OF SOLAR PROMINENCES

We investigate the influence of the geometry of the solar filament magnetic structure on the large-amplitude longitudinal oscillations. A representative filament flux tube is modeled as composed of a cool thread centered in a dipped part with hot coronal regions on either side. We have found the normal modes of the system and establish that the observed longitudinal oscillations are well described with the fundamental mode. For small and intermediate curvature radii and moderate to large density contrast between the prominence and the corona, the main restoring force is the solar gravity. In this full wave description of the oscillation a simple expression for the oscillation frequencies is derived in which the pressure-driven term introduces a small correction. We have also found that the normal modes are almost independent of the geometry of the hot regions of the tube. We conclude that observed large-amplitude longitudinal oscillations are driven by the projected gravity along the flux tubes and are strongly influenced by the curvature of the dips of the magnetic field in which the threads reside.

Rayleigh-Taylor instability in partially ionized compressible plasmas: One fluid approach

Aims. We study the modification of the classical criterion for the linear onset and growth rate of the Rayleigh-Taylor instability (RTI) in a partially ionized (PI) plasma in the one-fluid description by considering a generalized induction equation. Methods. The governing linear equations and appropriate boundary conditions, including gravitational terms, are derived and applied to the case of the RTI in a single interface between two partially ionized plasmas. The boundary conditions lead to an equation for the frequencies in which some have positive complex parts, marking the appearance of the RTI. We study the ambipolar term alone first, extending the result to the full induction equation later. Results. The configuration is always unstable because of the presence of a neutral species. In the classical stability regime, the growth rate is small, since the collisions prevent the neutral fluid to fully develop the RTI. For parameters in the classical instability regime, the growth rate is lowered, but the differences with the compressible MHD case are small for the considered theoretical values of the collision frequencies and diffusion coefficients for solar prominences. Conclusions. The PI modifies some aspects of the linear RTI instability, since it takes into account that neutrals do not feel the stabilizing effect of the magnetic field. For the set of parameters representative for solar prominences, our model gives the resulting timescale comparable to observed lifetimes of RTI plumes.

SPECTRO-POLARIMETRIC IMAGING REVEALS HELICAL MAGNETIC FIELDS IN SOLAR PROMINENCE FEET

Solar prominences are clouds of cool plasma levitating above the solar surface and insulated from the million-degree corona by magnetic fields. They form in regions of complex magnetic topology, characterized by non-potential fields, which can evolve abruptly, disintegrating the prominence and ejecting magnetized material into the heliosphere. However, their physics is not yet fully understood because mapping such complex magnetic configurations and their evolution is extremely challenging, and must often be guessed by proxy from photometric observations.Using state-of-the-art spectro-polarimetric data, we reconstruct the structure of the magnetic field in a prominence. We find that prominence feet harbor helical magnetic fields connecting the prominence to the solar surface below.

Effect of partial ionization on wave propagation in solar magnetic flux tubes

Observations show that waves are ubiquitous in the solar atmosphere and may play an important role for plasma heating. The study of waves in the solar corona is usually based on linear ideal magnetohydrodynamics (MHD) for a fully ionized plasma. However, the plasma in the photosphere and the chromosphere is only partially ionized. Here we theoretically investigate the impact of partial ionization on MHD wave propagation in cylindrical flux tubes in a two-fluid model. We derive the general dispersion relation that takes into account the effects of neutral-ion collisions and the neutral gas pressure. We assumed the neutral-ion collision frequency to be an arbitrary parameter. Specific results for transverse kink modes and slow magnetoacoustic modes are shown. We find that the wave frequencies only depend on the properties of the ionized fluid when the neutral-ion collision frequency is much lower that the wave frequency. For high collision frequencies that realistically represent the solar atmosphere, ions and neutrals behave as a single fluid with an effective density corresponding to the sum of densities of fluids plus an effective sound velocity computed as the average of the sound velocities of ions and neutrals. The MHD wave frequencies are modified accordingly. The neutral gas pressure can be neglected when studying transverse kink waves but it has to be included for a consistent description of slow magnetoacoustic waves. The MHD waves are damped by neutral-ion collisions. The damping is most efficient when the wave frequency and the collision frequency are on the same order of magnitude. For high collision frequencies slow magnetoacoustic waves are more efficiently damped than transverse kink waves. In addition, we find the presence of cut-offs for certain combinations of parameters that cause the waves to become non-propagating.

BAYESIAN ANALYSIS OF MULTIPLE HARMONIC OSCILLATIONS IN THE SOLAR CORONA

The detection of multiple mode harmonic kink oscillations in coronal loops enables us to obtain information on coronal density stratification and magnetic field expansion using seismology inversion techniques. The inference is based on the measurement of the period ratio between the fundamental mode and the first overtone and theoretical results for the period ratio under the hypotheses of coronal density stratification and magnetic field expansion of the wave guide. We present a Bayesian analysis of multiple mode harmonic oscillations for the inversion of the density scale height and magnetic flux tube expansion under each of the hypotheses. The two models are then compared using a Bayesian model comparison scheme to assess how plausible each one is given our current state of knowledge.

The effects of magnetic-field geometry on longitudinal oscillations of solar prominences: Cross-sectional area variation for thin tubes

Solar prominences are subject to both field-aligned (longitudinal) and transverse oscillatory motions, as evidenced by an increasing number of observations. Large-amplitude longitudinal motions provide valuable information on the geometry of the filament-channel magnetic structure that supports the cool prominence plasma against gravity. Our pendulum model, in which the restoring force is the gravity projected along the dipped field lines of the magnetic structure, best explains these oscillations. However, several factors can influence the longitudinal oscillations, potentially invalidating the pendulum model. The aim of this work is to study the influence of large-scale variations in the magnetic field strength along the field lines, i.e., variations of the cross-sectional area along the flux tubes supporting prominence threads. We studied the normal modes of several flux tube configurations, using linear perturbation analysis, to assess the influence of different geometrical parameters on the oscillation properties. We found that the influence of the symmetric and asymmetric expansion factors on longitudinal oscillations is small.}{We conclude that the longitudinal oscillations are not significantly influenced by variations of the cross-section of the flux tubes, validating the pendulum model in this context.

Evershed flow observed in neutral and singly ionized iron lines

The amplitudes of the Evershed flow are measured using pairs of carefully selected FeI and FeII spectral lines located close in wavelength and registered simultaneously. A sunspot belonging to the NOAA 11582 group was scanned using the spectrograph of the German Vacuum Tower Telescope (Observatorio del Teide, Tenerife). Velocities were extracted from intensity profiles using the lambda-meter technique. The formation heights of the observed spectral lines were calculated using semi-empirical models of a bright and dark penumbral filament taking into account the sunspot location at the limb. Our objective is to compare azimuthally averaged amplitudes of the Evershed flow extracted from neutral and ion lines. We find measurable differences in the radial component of the flow. All five pairs of lines show the same tendency, with a few hundred m/s larger amplitude of the flow measured from FeI lines compared to FeII lines. This tendency is preserved at all photospheric heights and radial distances in the penumbra. We discuss the possible origin of this effect.