M. Collados

HEATING OF THE MAGNETIZED SOLAR CHROMOSPHERE BY PARTIAL IONIZATION EFFECTS

In this paper, we study the heating of the magnetized solar chromosphere induced by the large fraction of neutral atoms present in this layer. The presence of neutrals, together with the decrease with height of the collisional coupling, leads to deviations from the classical magnetohydrodynamic behavior of the chromospheric plasma. A relative net motion appears between the neutral and ionized components, usually referred to as ambipolar diffusion. The dissipation of currents in the chromosphere is enhanced by orders of magnitude due to the action of ambipolar diffusion, as compared with the standard ohmic diffusion. We propose that a significant amount of magnetic energy can be released to the chromosphere just by existing force-free 10-40xa0G magnetic fields there. As a consequence, we conclude that ambipolar diffusion is an important process that should be included in chromospheric heating models, as it has the potential to rapidly heat the chromosphere. We perform analytical estimations and numerical simulations to prove this idea.

Multi-layer Study of Wave Propagation in Sunspots

We analyze the propagation of waves in sunspots from the photosphere to the chromosphere using time series of co-spatial Ca II H intensity spectra (including its line blends) and polarimetric spectra of Si I 10827 and the He I 10830 multiplet. From the Doppler shifts of these lines we retrieve the variation of the velocity along the line-of-sight at several heights. Phase spectra are used to obtain the relation between the oscillatory signals. Our analysis reveals standing waves at frequencies lower than 4 mHz and a continuous propagation of waves at higher frequencies, which steepen into shocks in the chromosphere when approaching the formation height of the Ca II H core. The observed non-linearities are weaker in Ca II H than in He I lines. Our analysis suggests that the Ca II H core forms at a lower height than the He I 10830 line: a time delay of about 20 s is measured between the Doppler signal detected at both wavelengths. We fit a model of linear slow magnetoacoustic wave propagation in a stratified atmosphere with radiative losses according to Newtons cooling law to the phase spectra and derive the difference in the formation height of the spectral lines. We show that the linear model describes well the wave propagation up to the formation height of Ca II H, where non-linearities start to become very important.

MAGNETO-ACOUSTIC WAVES IN SUNSPOTS: FIRST RESULTS FROM A NEW THREE-DIMENSIONAL NONLINEAR MAGNETOHYDRODYNAMIC CODE

Waves observed in the photosphere and chromosphere of sunspots show complex dynamics and spatial patterns. The interpretation of high-resolution sunspot wave observations requires modeling of three-dimensional (3D) nonlinear wave propagation and mode transformation in the sunspot upper layers in realistic spot model atmospheres. Here, we present the first results of such modeling. We have developed a 3D nonlinear numerical code specially designed to calculate the response of magnetic structures in equilibrium to an arbitrary perturbation. The code solves the 3D nonlinear MHD equations for perturbations; it is stabilized by hyper-diffusivity terms and is fully parallelized. The robustness of the code is demonstrated by a number of standard tests. We analyze several simulations of a sunspot perturbed by pulses of different periods at a subphotospheric level, from short periods, introduced for academic purposes, to longer and realistic periods of 3 and 5 minutes. We present a detailed description of the 3D mode transformation in a non-trivial sunspot-like magnetic field configuration, including the conversion between fast and slow magneto-acoustic waves and the Alfven wave, by calculation of the wave energy fluxes. Our main findings are as follows: (1) the conversion from acoustic to the Alfven mode is only observed if the driving pulse is located out of the sunspot axis, but this conversion is energetically inefficient; (2) as a consequence of the cutoff effects and refraction of the fast magneto-acoustic mode, the energy of the evanescent waves with periods around 5 minutes remains almost completely below the level β = 1; (3) waves with frequencies above the cutoff propagate field aligned to the chromosphere and their power becomes dominating over that of evanescent 5 minute oscillations, in agreement with observations.

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.

The energy of waves in the photosphere and lower chromosphere - I. Velocity statistics

Context. The energy source powering the solar chromosphere is still undetermined, but leaves its traces in observed intensities. Aims. We investigate the statistics of the intensity distributions as a function of the wavelength for Caii Ha nd the Caii IR line at 854.2 nm to estimate the energy content in the observed intensity fluctuations. Methods. We derived the intensity variations at different heights of the solar atmosphere, as traced by the line wings and line cores of the two spectral lines. We converted the observed intensities to absolute energy units employing reference profiles calculated in non-local thermal equilibrium (NLTE). We also converted the intensity fluctuations to corresponding brightness temperatures assuming LTE. Results. The root-mean-square (rms) fluctuations of the emitted intensity are about 0.6 (1.2) W m −2 ster −1 pm −1 near the core of the Caii IR line at 854.2 nm (Caii H), corresponding to relative intensity fluctuations of about 20% (30%). For the line wing, we find rms values of about 0.3 W m −2 ster −1 pm −1 for both lines, corresponding to relative fluctuations below 5%. The relative rms values show a local minimum for wavelengths forming at a height of about 130 km, but otherwise increase smoothly from the wing to the core, i.e., from photosphere to chromosphere. The corresponding rms brightness temperature fluctuations are below 100 K for the photosphere and up to 500 K in the chromosphere. The skewness of the intensity distributions is close to zero in the outer line wing and positive throughout the rest of the line spectrum, owing to the frequent occurrence of high-intensity events. The skewness shows a pronounced local maximum at locations with photospheric magnetic fields for wavelengths in-between those of the line wing and the line core (z ≈ 150−300 km), and a global maximum at the very core (z ≈ 1000 km) for both magnetic and field-free locations. Conclusions. The energy content of the intensity fluctuations is insufficient to create a chromospheric temperature rise that would be similar to the one in most reference models of the solar atmosphere. The increase in the rms fluctuations with height indicates the presence of upwardly propagating acoustic waves of increasing oscillation amplitude. The intensity and temperature variations indicate that there is a clear increase in dynamical activity from photosphere towards the chromosphere, but the variations fall short of the magnitude predicted by fully dynamical chromospheric models by a factor of about five. The enhanced skewness between the photosphere and lower solar chromosphere at magnetic locations is indicative of a mechanism that acts solely on magnetized plasma.

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.

Magnetohydrodynamic wave propagation from the subphotosphere to the corona in an arcade-shaped magnetic field with a null point

The aim of this work is to study the energy transport by means of MHD waves propagating in quiet Sun magnetic topology from layers below the surface to the corona. Upward propagating waves find obstacles, such as the equipartition layer with plasma b=1 and the transition region, and get converted, reflected and refracted. Understanding the mechanisms by which MHD waves can reach the corona can give us information about the solar atmosphere and the magnetic structures. We carry out two-dimensional numerical simulations of wave propagation in a magnetic field structure that consists of two vertical flux tubes separated by an arcade shaped magnetic field. This configuration contains a null point in the corona, that significantly modifies the behaviour of the waves. We describe in detail the wave propagation through the atmosphere under different driving conditions. We also present the spatial distribution of the mean acoustic and magnetic energy fluxes and the spatial distribution of the dominant frequencies in the whole domain. We conclude that the energy reaches the corona preferably along vertical magnetic fields, inside the flux tubes, and it has an acoustic nature. Most of the magnetic energy keeps concentrated below the transition region due to the refraction of the magnetic waves and the continuous conversion of acoustic-like waves into fast magnetic waves in the equipartition layer located in the photosphere. However, part of the magnetic energy reaches the low corona when propagating in the region where the arcades are located, but waves are sent back downwards to the lower atmosphere at the null point surroundings. This phenomenon, together with the reflection and refraction of waves in the TR and the lower turning point, act as a re-feeding of the atmosphere. In the frequency distribution, we find that high frequency waves can reach the corona outside the vertical flux tubes.

MAGNETOACOUSTIC WAVE ENERGY FROM NUMERICAL SIMULATIONS OF AN OBSERVED SUNSPOT UMBRA

We aim at reproducing the height dependence of sunspot wave signatures obtained from spectropolarimetric observations through three-dimensional MHD numerical simulations. A magnetostatic sunspot model based on the properties of the observed sunspot is constructed and perturbed at the photosphere, introducing the fluctuations measured with the Si I {lambda}10827 line. The results of the simulations are compared with the oscillations observed simultaneously at different heights from the He I {lambda}10830 line, the Ca II H core, and the Fe I blends in the wings of the Ca II H line. The simulations show a remarkable agreement with the observations. They reproduce the velocity maps and power spectra at the formation heights of the observed lines, as well as the phase and amplification spectra between several pairs of lines. We find that the stronger shocks at the chromosphere are accompanied with a delay between the observed signal and the simulated one at the corresponding height, indicating that shocks shift the formation height of the chromospheric lines to higher layers. Since the simulated wave propagation matches very well the properties of the observed one, we are able to use the numerical calculations to quantify the energy contribution of the magnetoacoustic waves to themorexa0» chromospheric heating in sunspots. Our findings indicate that the energy supplied by these waves is too low to balance the chromospheric radiative losses. The energy contained at the formation height of the lowermost Si I {lambda}10827 line in the form of slow magnetoacoustic waves is already insufficient to heat the higher layers, and the acoustic energy which reaches the chromosphere is around 3-9 times lower than the required amount of energy. The contribution of the magnetic energy is even lower.«xa0less

ON THE ROBUSTNESS OF THE PENDULUM MODEL FOR LARGE-AMPLITUDE LONGITUDINAL OSCILLATIONS IN PROMINENCES

Large-amplitude longitudinal oscillations (LALOs) in prominences are spectacular manifestations of the solar activity. In such events nearby energetic disturbances induce periodic motions on filaments with displacements comparable to the size of the filaments themselves and with velocities larger than 20 km/s. The pendulum model, in which the gravity projected along a rigid magnetic field is the restoring force, was proposed to explain these events. However, it can be objected that in a realistic situation where the magnetic field reacts to the mass motion of the heavy prominence, the simplified pendulum model could be no longer valid. We have performed non-linear time-dependent numerical simulations of LALOs considering a dipped magnetic field line structure. In this work we demonstrate that for even relatively weak magnetic fields the pendulum model works very well. We therefore validate the pendulum model and show its robustness, with important implications for prominence seismology purposes. With this model it is possible to infer the geometry of the dipped field lines that support the prominence.

Deep probing of the photospheric sunspot penumbra: no evidence of field-free gaps

Context. Some models for the topology of the magnetic field in sunspot penumbrae predict regions free of magnetic fields or with only dynamically weak fields in the deep photosphere. Aims. We aim to confirm or refute the existence of weak-field regions in the deepest photospheric layers of the penumbra. Methods. We investigated the magnetic field at logu2009 τ 5 = 0 is by inverting spectropolarimetric data of two different sunspots located very close to disk center with a spatial resolution of approximately 0.4−0.45′′. The data have been recorded using the GRIS instrument attached to the 1.5-m solar telescope GREGOR at the El Teide observatory. The data include three Fe iu2009lines around 1565 nm, whose sensitivity to the magnetic field peaks half a pressure scale height deeper than the sensitivity of the widely used Fe iu2009spectral line pair at 630 nm. Before the inversion, the data were corrected for the effects of scattered light using a deconvolution method with several point spread functions. Results. At logu2009 τ 5 = 0 we find no evidence of regions with dynamically weak ( B < 500 Gauss) magnetic fields in sunspot penumbrae. This result is much more reliable than previous investigations made on Fe iu2009lines at 630 nm. Moreover, the result is independent of the number of nodes employed in the inversion, is independent of the point spread function used to deconvolve the data, and does not depend on the amount of stray light (i.e., wide-angle scattered light) considered.