Magnetohydrodynamic (MHD) kink waves have been observed frequently in solar coronal flux tubes, which makes them a great tool for seismology of the solar corona. Here, the effect of viscosity is studied on the evolution of kink waves. To this aim, we solve the initial value problem for the incompressible linearized viscous MHD equations in a radially inhomogeneous flux tube in the limit of long wavelengths. Using a modal expansion technique the spatio-temporal behavior of the perturbations is obtained. We confirm that for large Reynolds numbers representative of the coronal plasma the decrement in the amplitude of the kink oscillations is due to the resonant absorption mechanism that converts the global transverse oscillation to rotational motions in the inhomogeneous layer of the flux tube. We show that viscosity suppresses the rate of phase mixing of the perturbations in the inhomogeneous region of the flux tube and prevents the continuous building up of small scales in the system once a sufficiently small scale is reached. The viscous dissipation function is calculated to investigate plasma heating by viscosity in the inhomogeneous layer of the flux tube. For Reynolds numbers of the order of 10^6-10^8, the energy of the kink wave is transformed into heat in two to eight periods of the kink oscillation. For larger and more realistic Reynolds numbers, heating happens, predominantly, after the global kink oscillation is damped, and no significant heating occurs during the observable transverse motion of the flux tube.

Observations show that transverse magnetohydrodynamic (MHD) waves and flows are often simultaneously present in magnetic loops of the solar corona. The waves are resonantly damped in the Alfv\'en continuum because of plasma and/or magnetic field nonuniformity across the loop. The resonant damping is relevant in the context of coronal heating, since it provides a mechanism to cascade energy down to the dissipative scales. It has been theoretically shown that the presence of flow affects the waves propagation and damping, but most of the studies rely on the unjustified assumption that the transverse nonuniformity is confined to a boundary layer much thinner than the radius of the loop. Here we present a semi-analytic technique to explore the effect of flow on resonant MHD waves in coronal flux tubes with thick nonuniform boundaries. We extend a published method, which was originally developed for a static plasma, in order to incorporate the effect of flow. We allowed the flow velocity to continuously vary within the nonuniform boundary from the internal velocity to the external velocity. The analytic part of the method is based on expressing the wave perturbations in the thick nonuniform boundary of the loop as a Frobenius series that contains a singular term accounting for the Alfv\'en resonance, while the numerical part of the method consists of solving iteratively the transcendental dispersion relation together with the equation for the Alfv\'en resonance position. As an application of this method, we investigated the impact of flow on the phase velocity and resonant damping length of MHD kink waves. We consistently recover results in the thin boundary approximation obtained in previous studies. We have extended those results to the case of thick boundaries. We also explored the error associated with the use of the thin boundary approximation beyond its regime of applicability.

In the solar atmosphere, Alfv\'en waves are believed to play an important role in the transfer of energy from the photosphere to the corona and solar wind, and in the heating of the chromosphere. We perform numerical computations to investigate energy transport and dissipation associated with torsional Alfv\'en waves propagating in magnetic flux tubes that expand from the photosphere to the corona in quiet-Sun conditions. We place a broadband driver at the photosphere that injects a wave energy flux of $10^7$ erg cm$^{-2}$ s$^{-1}$ and consider Ohm's magnetic diffusion and ion-neutral collisions as dissipation mechanisms. We find that only a small fraction of the driven flux, $\sim 10^5$ erg cm$^{-2}$ s$^{-1}$, is able to reach coronal heights, but it may be sufficient to partly compensate the total coronal energy loss. The frequency of maximal transmittance is $\sim 5$ mHz for a photospheric field strength of 1 kG and is shifted to smaller/larger frequencies for weaker/stronger fields. Lower frequencies are reflected at the transition region, while higher frequencies are dissipated producing enough heat to balance chromospheric radiative losses. Heating in the low and middle chromosphere is due to Ohmic dissipation, while ion-neutral friction dominates in the high chromosphere. Ohmic diffusion is enhanced by phase mixing because of the expansion of the magnetic field. This effect has the important consequence of increasing the chromospheric dissipation and, therefore, reducing the energy flux that reaches the corona. We provide empirical fits of the transmission coefficient that could be used as input for coronal models.

The multi-fluid modelling of high-frequency waves in partially ionized plasmas has shown that the behavior of magnetohydrodynamics waves in the linear regime is heavily influenced by the collisional interaction between the different species that form the plasma. Here, we go beyond linear theory and study large-amplitude waves in partially ionized plasmas using a nonlinear multi-fluid code. It is known that in fully ionized plasmas, nonlinear Alfv\'en waves generate density and pressure perturbations. Those nonlinear effects are more pronounced for standing oscillations than for propagating waves. By means of numerical simulations and analytical approximations, we examine how the collisional interaction between ions and neutrals affects the nonlinear evolution. The friction due to collisions dissipates a fraction of the wave energy, which is transformed into heat and consequently rises the temperature of the plasma. As an application, we investigate frictional heating in a plasma with physical conditions akin to those in a solar quiescent prominence.

There are observational evidences for the existence of twisted magnetic field in the solar corona. This inspires us to investigate the effect of a twisted magnetic field on the evolution of magnetohydrodynamic (MHD) kink waves in coronal loops. To this aim, we solve the incompressible linearized MHD equations in a magnetically twisted nonuniform coronal flux tube in the limit of long wavelengths. Our results show that a twisted magnetic field can enhance or diminish the rate of phase-mixing of the Alfv\'{e}n continuum modes and the decay rate of the global kink oscillation depending on the twist model and the sign of the longitudinal ($k_z$) and azimuthal ($m$) wavenumbers. Also our results confirm that in the presence of a twisted magnetic field, when the sign of one of the two wavenumbers $m$ and $k_z$ is changed, the symmetry with respect to the propagation direction is broken. Even a small amount of twist can have an important impact on the process of energy cascade to small scales.

Partially ionized plasmas are found across the Universe in many different astrophysical environments. They constitute an essential ingredient of the solar atmosphere, molecular clouds, planetary ionospheres and protoplanetary disks, among other environments, and display a richness of physical effects which are not present in fully ionized plasmas. This review provides an overview of the physics of partially ionized plasmas, including recent advances in different astrophysical areas in which partial ionization plays a fundamental role. We outline outstanding observational and theoretical questions and discuss possible directions for future progress.

The presence of neutral species in a plasma has been shown to greatly affect the properties of magnetohydrodynamic waves. For instance, the interaction between ions and neutrals through momentum transfer collisions causes the damping of Alfv\'en waves and alters their oscillation frequency and phase speed. When the collision frequencies are larger than the frequency of the waves, single-fluid magnetohydrodynamic approximations can accurately describe the effects of partial ionization, since there is a strong coupling between the various species. However, at higher frequencies, the single-fluid models are not applicable and more complex approaches are required. Here, we use a five-fluid model with three ionized and two neutral components, which takes into consideration Hall's current and Ohm's diffusion in addition to the friction due to collisions between different species. We apply our model to plasmas composed of hydrogen and helium, and allow the ionization degree to be arbitrary. By means of the analysis of the corresponding dispersion relation and numerical simulations, we study the properties of small-amplitude perturbations. We discuss the effect of momentum transfer collisions on the ion-cyclotron resonances and compare the importance of magnetic resistivity, ion-neutral and ion-ion collisions on the wave damping at various frequency ranges. Application to partially ionized plasmas of the solar atmosphere are performed.

Ideal MHD provides an accurate description of low-frequency Alfv\'en waves in fully ionized plasmas. However, higher frequency waves in many plasmas of the solar atmosphere cannot be correctly described by ideal MHD and a more accurate model is required. Here, we study the properties of small-amplitude incompressible perturbations in both the low and the high frequency ranges in plasmas composed of several ionized species. We use a multi-fluid approach and take into account the effects of collisions between ions and the inclusion of Hall's term in the induction equation. Through the analysis of the corresponding dispersion relations and numerical simulations we check that at high frequencies ions of different species are not as strongly coupled as in the low frequency limit. Hence, they cannot be treated as a single fluid. In addition, elastic collisions between the distinct ionized species are not negligible for high frequency waves since an appreciable damping is obtained. Furthermore, Coulomb collisions between ions remove the cyclotron resonances and the strict cut-off regions that are present when collisions are not taken into account. The implications of these results for the modelling of high-frequency waves in solar plasmas are discussed.

Observations have shown that magnetohydrodynamic waves over a large frequency range are ubiquitous in solar prominences. The waves are probably driven by photospheric motions and may transport energy up to prominences suspended in the corona. Dissipation of wave energy can lead to heating of the cool prominence plasma, so contributing to the local energy balance within the prominence. Here we discuss the role of Alfv\'en wave dissipation as a heating mechanism for the prominence plasma. We consider a slab-like quiescent prominence model with a transverse magnetic field embedded in the solar corona. The prominence medium is modelled as a partially ionized plasma composed of a charged ion-electron single fluid and two separate neutral fluids corresponding to neutral hydrogen and neutral helium. Friction between the three fluids acts as a dissipative mechanism for the waves. The heating caused by externally-driven Alfv\'en waves incident on the prominence slab is analytically explored. We find that the dense prominence slab acts as a resonant cavity for the waves. The fraction of incident wave energy that is channelled into the slab strongly depends upon the wave period, $P$. Using typical prominence conditions, we obtain that wave energy trapping and associated heating are negligible when $P \gtrsim 100$ s, so that it is unlikely that those waves have a relevant influence on prominence energetics. When $1$ s $\lesssim P \lesssim 100$ s the energy absorption into the slab shows several sharp and narrow peaks, that can reach up to 100%, when the incident wave frequency matches a cavity resonance of the slab. Wave heating is enhanced at those resonant frequencies. Conversely, when $P \lesssim 1$ s cavity resonances are absent, but the waves are heavily damped by the strong dissipation. We estimate that wave heating may compensate for about 10% of radiative losses of the prominence plasma.

This document was created by the Solar Simulations for the Atacama Large Millimeter Observatory Network (SSALMON) in preparation of the first regular observations of the Sun with the Atacama Large Millimeter/submillimeter Array (ALMA), which are anticipated to start in ALMA Cycle 4 in October 2016. The science cases presented here demonstrate that a large number of scientifically highly interesting observations could be made already with the still limited solar observing modes foreseen for Cycle 4 and that ALMA has the potential to make important contributions to answering long-standing scientific questions in solar physics. With the proposal deadline for ALMA Cycle 4 in April 2016 and the Commissioning and Science Verification campaign in December 2015 in sight, several of the SSALMON Expert Teams composed strategic documents in which they outlined potential solar observations that could be feasible given the anticipated technical capabilities in Cycle 4. These documents have been combined and supplemented with an analysis, resulting in recommendations for solar observing with ALMA in Cycle 4. In addition, the detailed science cases also demonstrate the scientific priorities of the solar physics community and which capabilities are wanted for the next observing cycles. The work on this White Paper effort was coordinated in close cooperation with the two international solar ALMA development studies led by T. Bastian (NRAO, USA) and R. Brajsa, (ESO). This document will be further updated until the beginning of Cycle 4 in October 2016. In particular, we plan to adjust the technical capabilities of the solar observing modes once finally decided and to further demonstrate the feasibility and scientific potential of the included science cases by means of numerical simulations of the solar atmosphere and corresponding simulated ALMA observations.

The Atacama Large Millimeter/submillimeter Array (ALMA) is a new powerful tool for observing the Sun at high spatial, temporal, and spectral resolution. These capabilities can address a broad range of fundamental scientific questions in solar physics. The radiation observed by ALMA originates mostly from the chromosphere - a complex and dynamic region between the photosphere and corona, which plays a crucial role in the transport of energy and matter and, ultimately, the heating of the outer layers of the solar atmosphere. Based on first solar test observations, strategies for regular solar campaigns are currently being developed. State-of-the-art numerical simulations of the solar atmosphere and modeling of instrumental effects can help constrain and optimize future observing modes for ALMA. Here we present a short technical description of ALMA and an overview of past efforts and future possibilities for solar observations at submillimeter and millimeter wavelengths. In addition, selected numerical simulations and observations at other wavelengths demonstrate ALMA's scientific potential for studying the Sun for a large range of science cases.

The temporal evolution of a solar prominence inserted in a three-dimensional magnetic flux rope is investigated numerically. Using the model of Titov Demoulin (1999) under the regime of weak twist, the cold and dense prominence counteracts gravity by modifying the initially force-free magnetic configuration. In some cases a quasi-stationary situation is achieved after the relaxation phase, characterized by the excitation of standing vertical oscillations. These oscillations show a strong attenuation with time produced by the mechanism of continuum damping due to the inhomogeneous transition between the prominence and solar corona. The characteristic period of the vertical oscillations does not depend strongly on the twist of the flux rope. Nonlinearity is the responsible for triggering the Kelvin-Helmholtz instability associated to the vertical oscillations and that eventually produces horizontal structures. Contrary to other configurations in which the longitudinal axis of the prominence is permeated by a perpendicular magnetic field, like in unsheared arcades, the orientation of the prominence along the flux rope axis prevents the development of Rayleigh-Taylor instabilities and therefore the appearance of vertical structuring along this axis.

Coronal rain clumps and prominence knots are dense condensations with chromospheric to transition region temperatures that fall down in the much hotter corona. Their typical speeds are in the range 30--150~km~s$^{-1}$ and of the order of 10--30~km~s$^{-1}$, respectively, i.e., they are considerably smaller than free fall velocities. These cold blobs contain a mixture of ionized and neutral material that must be dynamically coupled in order to fall together, as observed. We investigate this coupling by means of hydrodynamic simulations in which the coupling arises from the friction between ions and neutrals. The numerical simulations presented here are an extension of those of \citet{oliver2014} to the partially ionized case. We find that, although the relative drift speed between the two species is smaller than 1~m~s$^{-1}$ at the blob center, it is sufficient to produce the forces required to strongly couple charged particles and neutrals. The ionization degree has no discernible effect on the main results of our previous work for a fully ionized plasma: the condensation has an initial acceleration phase followed by a period with roughly constant velocity and, in addition, the maximum descending speed is clearly correlated with the ratio of initial blob to environment density.

In this paper we show that the phase mixing of continuum Alfv\'{e}n waves and/or continuum slow waves in magnetic structures of the solar atmosphere as, e.g., coronal arcades, can create the illusion of wave propagation across the magnetic field. This phenomenon could be erroneously interpreted as fast magnetosonic waves. The cross-field propagation due to phase mixing of continuum waves is apparent because there is no real propagation of energy across the magnetic surfaces. We investigate the continuous Alfv\'{e}n and slow spectra in 2D Cartesian equilibrium models with a purely poloidal magnetic field. We show that apparent superslow propagation across the magnetic surfaces in solar coronal structures is a consequence of the existence of continuum Alfv\'{e}n waves and continuum slow waves that naturally live on those structures and phase mix as time evolves. The apparent cross-field phase velocity is related to the spatial variation of the local Alfv\'{e}n/slow frequency across the magnetic surfaces and is slower than the Alfv\'{e}n/sound velocities for typical coronal conditions. Understanding the nature of the apparent cross-field propagation is important for the correct analysis of numerical simulations and the correct interpretation of observations.

The spatial variation of physical quantities, such as the mass density, across solar atmospheric waveguides governs the timescales and spatial scales for wave damping and energy dissipation. The direct measurement of the spatial distribution of density, however, is difficult and indirect seismology inversion methods have been suggested as an alternative. We applied Bayesian inference, model comparison, and model-averaging techniques to the inference of the cross-field density structuring in solar magnetic waveguides using information on periods and damping times for resonantly damped magnetohydrodynamic (MHD) transverse kink oscillations. Three commonly employed alternative profiles were used to model the variation of the mass density across the waveguide boundary. Parameter inference enabled us to obtain information on physical quantities such as the Alfv\'en travel time, the density contrast, and the transverse inhomogeneity length scale. The inference results from alternative density models were compared and their differences quantified. Then, the relative plausibility of the considered models was assessed by performing model comparison. Our results indicate that the evidence in favor of any of the three models is minimal, unless the oscillations are strongly damped. In such a circumstance, the application of model-averaging techniques enables the computation of an evidence-weighted inference that takes into account the plausibility of each model in the calculation of a combined inversion for the unknown physical parameters.

Damped transverse oscillations of magnetic loops are routinely observed in the solar corona. This phenomenon is interpreted as standing kink magnetohydrodynamic waves, which are damped by resonant absorption owing to plasma inhomogeneity across the magnetic field. The periods and damping times of these oscillations can be used to probe the physical conditions of the coronal medium. Some observations suggest that interaction between neighboring oscillating loops in an active region may be important and can modify the properties of the oscillations compared to those of an isolated loop. Here we theoretically investigate resonantly damped transverse oscillations of interacting non-uniform coronal loops. We provide a semi-analytic method, based on the T-matrix theory of scattering, to compute the frequencies and damping rates of collective oscillations of an arbitrary configuration of parallel cylindrical loops. The effect of resonant damping is included in the T-matrix scheme in the thin boundary approximation. Analytic and numerical results in the specific case of two interacting loops are given as an application.

Dissipation of magnetohydrodynamic (MHD) wave energy has been proposed as a viable heating mechanism in the solar chromospheric plasma. Here, we use a simplified one-dimensional model of the chromosphere to theoretically investigate the physical processes and the spatial scales that are required for the efficient dissipation of Alfv\'en waves and slow magnetoacoustic waves. We consider the governing equations for a partially ionized hydrogen-helium plasma in the single-fluid MHD approximation and include realistic wave damping mechanisms that may operate in the chromosphere, namely Ohmic and ambipolar magnetic diffusion, viscosity, thermal conduction, and radiative losses. We perform an analytic local study in the limit of small amplitudes to approximately derive the lengthscales for critical damping and efficient dissipation of MHD wave energy. We find that the critical dissipation lengthscale for Alfv\'en waves depends strongly on the magnetic field strength and ranges from 10~m to 1~km for realistic field strengths. The damping of Alfv\'en waves is dominated by Ohmic diffusion for weak magnetic field and low heights in the chromosphere, and by ambipolar diffusion for strong magnetic field and medium/large heights in the chromosphere. Conversely, the damping of slow magnetoacoustic waves is less efficient, and spatial scales shorter than 10~m are required for critical damping. Thermal conduction and viscosity govern the damping of slow magnetoacoustic waves and play an equally important role at all heights. These results indicate that the spatial scales at which strong wave heating may work in the chromosphere are currently unresolved by observations.

Quiescent solar prominence fine structures are typically modelled as density enhancements, called threads, which occupy a fraction of a longer magnetic flux tube. The profile of the mass density along the magnetic field is however unknown and several arbitrary alternatives are employed in prominence wave studies. We present a comparison of theoretical models for the field-aligned density along prominence fine structures. We consider Lorentzian, Gaussian, and parabolic profiles. We compare their theoretical predictions for the period ratio between the fundamental transverse kink mode and the first overtone to obtain estimates for the ratio of densities between the central part of the tube and its foot-points and to assess which one would better explain observed period ratio data. Bayesian parameter inference and model comparison techniques are developed and applied. Parameter inference requires the computation of the posterior distribution for the density gradient parameter conditional on the observable period ratio. Model comparison involves the computation of the marginal likelihood as a function of the period ratio to obtain the plausibility of each density model and the computation of Bayes Factors to quantify the relative evidence for each model, given a period ratio observation. A Lorentzian density profile, with plasma density concentrated around the centre of the tube seems to offer the most plausible inversion result. A Gaussian profile would require unrealistically large values of the density gradient parameter and a parabolic density distribution does not enable us to obtain well constrained posterior probability distributions. However, our model comparison results indicate that the evidence points to the Gaussian and parabolic profiles for period ratios in between 2 and 3, while the Lorentzian profile is preferred for larger period ratio values.

Context. Recent observations of solar prominences show the presence of turbulent flows that may be caused by Kelvin-Helmholtz instabilites (KHI). However, the observed flow velocities are below the classical threshold for the onset of KHI in fully ionized plasmas. Aims. We investigate the effect of partial ionization on the onset of KHI in dense and cool cylindrical magnetic flux tubes surrounded by a hotter and lighter environment. Methods. The linearized governing equations of a partially ionized two-fluid plasma are used to describe the behavior of small-amplitude perturbations superimposed on a magnetic tube with longitudinal mass flow. A normal mode analysis is performed to obtain the dispersion relation for linear incompressible waves. We focus on the appearance of unstable solutions and study the dependence of their growth rates on various physical parameters. An analytical approximation of the KHI linear growth rate for slow flows and strong ion-neutral coupling is obtained. An application to solar prominence threads is given. Results. The presence of a neutral component in a plasma may contribute to the onset of the KHI even for sub-Alfv\'enic longitudinal shear flows. Collisions between ions and neutrals reduce the growth rates of the unstable perturbations but cannot completely suppress the instability. Conclusions. Turbulent flows in solar prominences with sub-Alfv\'enic flow velocities may be interpreted as consequences of KHI in partially ionized plasmas.

The Solar Simulations for the Atacama Large Millimeter Observatory Network (SSALMON) was initiated in 2014 in connection with two ALMA development studies. The Atacama Large Millimeter/submillimeter Array (ALMA) is a powerful new tool, which can also observe the Sun at high spatial, temporal, and spectral resolution. The international SSALMONetwork aims at coordinating the further development of solar observing modes for ALMA and at promoting scientific opportunities for solar physics with particular focus on numerical simulations, which can provide important constraints for the observing modes and can aid the interpretation of future observations. The radiation detected by ALMA originates mostly in the solar chromosphere - a complex and dynamic layer between the photosphere and corona, which plays an important role in the transport of energy and matter and the heating of the outer layers of the solar atmosphere. Potential targets include active regions, prominences, quiet Sun regions, flares. Here, we give a brief overview over the network and potential science cases for future solar observations with ALMA.

Cool and dense prominences found in the solar atmosphere are known to be partially ionized because of their relative low temperature. In this Letter, we address the long-standing problem of how the neutral component of the plasma in prominences is supported against gravity. Using the multiple fluid approach we solve the time-dependent equations in two dimensions considering the frictional coupling between the neutral and ionized components of the magnetized plasma representative of a solar prominence embedded in a hot coronal environment. We demonstrate that given an initial density enhancement in the two fluids, representing the body of the prominence, the system is able to relax in the vicinity of magnetic dips to a stationary state in which both neutrals and ionized species are dynamically suspended above the photosphere. Two different coupling processes are considered in this study, collisions between ions and neutrals and charge exchange interactions. We find that for realistic conditions ions are essentially static while neutrals have a very small downflow velocity. The coupling between ions and neutrals is so strong at the prominence body that the behavior is similar to that of a single fluid with an effective density equal to the sum of the ion and neutral species. We also find that the charge exchange mechanism is about three times more efficient sustaining neutrals than elastic scattering of ions with neutrals.

Magnetohydrodynamic (MHD) kink waves are ubiquitously observed in the solar atmosphere. The propagation and damping of these waves may play relevant roles for the transport and dissipation of energy in the solar atmospheric medium. However, in the atmospheric plasma dissipation of transverse MHD wave energy by viscosity or resistivity needs very small spatial scales to be efficient. Here, we theoretically investigate the generation of small scales in nonuniform solar magnetic flux tubes due to phase mixing of MHD kink waves. We go beyond the usual approach based on the existence of a global quasi-mode that is damped in time due to resonant absorption. Instead, we use a modal expansion to express the MHD kink wave as a superposition of Alfv\'en continuum modes that are phase mixed as time evolves. The comparison of the two techniques evidences that the modal analysis is more physically transparent and describes both the damping of global kink motions and the building up of small scales due to phase mixing. In addition, we discuss that the processes of resonant absorption and phase mixing are intimately linked. They represent two aspects of the same underlying physical mechanism: the energy cascade from large scales to small scales due to naturally occurring plasma and/or magnetic field inhomogeneities. This process may provide the necessary scenario for efficient dissipation of transverse MHD wave energy in the solar atmospheric plasma.

The ratio of the period of the fundamental mode to that of the first overtone of kink oscillations, from here on the "period ratio", is a seismology tool that can be used to infer information about the spatial variation of density along solar magnetic flux tubes. The period ratio is 2 in longitudinally homogeneous thin tubes, but it differs from 2 due to longitudinal inhomogeneity. In this paper we investigate the period ratio in longitudinally inhomogeneous prominence threads and explore its implications for prominence seismology. We numerically solve the two-dimensional eigenvalue problem of kink oscillations in a model of a prominence thread. We take into account three nonuniform density profiles along the thread. In agreement with previous works that used simple piecewise constant density profiles, we find that the period ratio is larger than 2 in prominence threads. When the ratio of the central density to that at the footpoints is fixed, the period ratio depends strongly on the form of the density profile along the thread. The more concentrated the dense prominence plasma near the center of the tube, the larger the period ratio. However, the period ratio is found to be independent of the specific density profile when the spatially averaged density in the thread is the same for all the profiles. An empirical fit of the dependence of the period ratio on the average density is given and its use for prominence seismology is discussed.

In this paper we present a numerical study of the time evolution of solar prominences embedded in sheared magnetic arcades. The prominence is represented by a density enhancement in a background stratified atmosphere and is connected to the photosphere through the magnetic field. By solving the ideal magnetohydrodynamic (MHD) equations in three dimensions we study the dynamics for a range of parameters representative of real prominences. Depending on the parameters considered, we find prominences that are suspended above the photosphere, i.e., detached prominences, but also configurations resembling curtain or hedgerow prominences whose material continuously connects to the photosphere. The plasma$-\beta$ is an important parameter that determines the shape of the structure. In many cases magnetic Rayleigh-Taylor (MRT) instabilities and oscillatory phenomena develop. Fingers and plumes are generated, affecting the whole prominence body and producing vertical structures in an essentially horizontal magnetic field. However, magnetic shear is able to reduce or even to suppress this instability.

Alfvenic waves are ubiquitous in the solar atmosphere and their dissipation may play an important role in atmospheric heating. In the partially ionized solar chromosphere, collisions between ions and neutrals are an efficient dissipative mechanism for Alfven waves with frequencies near the ion-neutral collision frequency. The collision frequency is proportional to the ion-neutral collision cross section for momentum transfer. Here, we investigate Alfven wave damping as a function of height in a simplified chromospheric model and compare the results for two sets of collision cross sections, namely those of the classic hard-sphere model and those based on recent quantum-mechanical computations. We find important differences between the results for the two sets of cross sections. There is a critical interval of wavelengths for which impulsively excited Alfven waves are overdamped as a result of the strong ion-neutral dissipation. The critical wavelengths are in the range from 1 km to 50 km for the hard-sphere cross sections, and from 1 m to 1 km for the quantum-mechanical cross sections. Equivalently, for periodically driven Alfven waves there is an optimal frequency for which the damping is most effective. The optimal frequency varies from 1 Hz to 10^2 Hz for the hard-sphere cross sections, and from 10^2 Hz to 10^4 Hz for the quantum-mechanical cross sections. Future observations at sufficiently high spatial or temporal resolution may show the importance of high-frequency Alfven waves for chromospheric heating. For instance, the Atacama Large Millimeter/submillimeter Array (ALMA) may be able to detect the critical wavelengths and optimal frequencies and so to test the effective collision cross section in the chromospheric plasma.

Seismology of coronal loops using observations of damped transverse oscillations in combination with results from theoretical models is a tool to indirectly infer physical parameters in the solar atmospheric plasma. Existing seismology schemes based on approximations to the period and damping time of kink oscillations are often used beyond their theoretical range of applicability. These approximations assume that the variation of density across the loop is confined to a nonuniform layer much thinner than the radius of the loop, but the results of the inversion problem often do not satisfy this preliminary hypothesis. Here, we determine the accuracy of the analytic approximations to the period and damping time, and its impact on seismology estimates, when largely nonuniform loops are considered. We find that the accuracy of the approximations when used beyond their range of applicability is strongly affected by the form of the density profile across the loop, that is observationally unknown and so must be arbitrarily imposed as part of the theoretical model. The error associated with the analytic approximations can be larger than 50% even for relatively thin nonuniform layers. This error directly affects the accuracy of approximate seismology estimates compared to actual numerical inversions. In addition, assuming different density profiles can produce noncoincident intervals of the seismic variables in inversions of the same event. The ignorance about the true shape of density variation across the loop is an important source of error that may dispute the reliability of parameters seismically inferred assuming an ad hoc density profile.

Observations of active regions and limb prominences often show cold, dense blobs descending with an acceleration smaller than that of free fall. The dynamics of these condensations falling in the solar corona is investigated in this paper using a simple fully ionised plasma model. We find that the presence of a heavy condensation gives rise to a dynamical rearrangement of the coronal pressure that results in the formation of a large pressure gradient that opposes gravity. Eventually this pressure gradient becomes so large that the blob acceleration vanishes or even points upwards. Then, the blob descent is characterised by an initial acceleration phase followed by an essentially constant velocity phase. These two stages can be identified in published time-distance diagrams of coronal rain events. Both the duration of the first stage and the velocity attained by the blob increase for larger values of the ratio of blob to coronal density, for larger blob mass, and for smaller coronal temperature. Dense blobs are characterised by a detectable density growth (up to 60% in our calculations) and by a steepening of the density in their lower part, that could lead to the formation of a shock. They also emit sound waves that could be detected as small intensity changes with periods of the order of 100 s and lasting between a few and about ten periods. Finally, the curvature of the falling path is only relevant when a very dense blob falls along inclined magnetic field lines.

Compressible disturbances propagate in a plasma in the form of magnetoacoustic waves driven by both gas pressure and magnetic forces. In partially ionized plasmas the dynamics of ionized and neutral species are coupled due to ion-neutral collisions. As a consequence, magnetoacoustic waves propagating through a partially ionized medium are affected by the ion-neutral coupling. The degree to which the behavior of the classic waves is modified depends on the physical properties of the various species and on the relative value of the wave frequency compared to the ion-neutral collision frequency. Here, we perform a comprehensive theoretical investigation of magnetoacoustic wave propagation in a partially ionized plasma using the two-fluid formalism. We consider an extensive range of values for the collision frequency, ionization ratio, and plasma $\beta$, so that the results are applicable to a wide variety of astrophysical plasmas. We determine the modification of the wave frequencies and study the frictional damping due to ion-neutral collisions. Approximate analytic expressions to the frequencies are given in the limit case of strongly coupled ions and neutrals, while numerically obtained dispersion diagrams are provided for arbitrary collision frequencies. In addition, we discuss the presence of cutoffs in the dispersion diagrams that constrain wave propagation for certain combinations of parameters. A specific application to propagation of compressible waves in the solar chromosphere is given.

Solar prominence models used so far in the analysis of MHD waves in such structures are quite elementary. In this work, we calculate numerically magnetohydrostatic models in two-dimensional configurations under the presence of gravity. Our interest is in models that connect the magnetic field to the photosphere and include an overlying arcade. The method used here is based on a relaxation process and requires solving the time-dependent nonlinear ideal MHD equations. Once a prominence model is obtained, we investigate the properties of MHD waves superimposed on the structure. We concentrate on motions purely two-dimensional neglecting propagation in the ignorable direction. We demonstrate how by using different numerical tools we can determine the period of oscillation of stable waves. We find that vertical oscillations, linked to fast MHD waves, are always stable and have periods in the 4-10 min range. Longitudinal oscillations, related to slow magnetoacoustic-gravity waves, have longer periods in the range of 28-40 min. These longitudinal oscillations are strongly influenced by the gravity force and become unstable for short magnetic arcades.

Magnetohydrodynamic (MHD) waves are ubiquitously observed in the solar atmosphere. Kink waves are a type of transverse MHD waves in magnetic flux tubes that are damped due to resonant absorption. The theoretical study of kink MHD waves in solar flux tubes is usually based on the simplification that the transverse variation of density is confined to a nonuniform layer much thinner than the radius of the tube, i.e., the so-called thin boundary approximation. Here, we develop a general analytic method to compute the dispersion relation and the eigenfunctions of ideal MHD waves in pressureless flux tubes with transversely nonuniform layers of arbitrary thickness. Results for kink waves are produced and are compared with fully numerical resistive MHD eigenvalue computations in the limit of small resistivity. We find that the frequency and resonant damping rate are the same in both ideal and resistive cases. The actual results for thick nonuniform layers deviate from the behavior predicted in the thin boundary approximation and strongly depend on the shape of the nonuniform layer. The eigenfunctions in ideal MHD are very different from those in resistive MHD. The ideal eigenfunctions display a global character regardless of the thickness of the nonuniform layer, while the resistive eigenfunctions are localized around the resonance and are indistinguishable from those of ordinary resistive Alfv\'en modes. Consequently, the spatial distribution of wave energy in the ideal and resistive cases is dramatically different. This poses a fundamental theoretical problem with clear observational consequences.

Ion-neutral collisions may lead to the damping of Alfven waves in chromospheric and prominence plasmas. Neutral helium atoms enhance the damping in certain temperature interval, where the ratio of neutral helium and neutral hydrogen atoms is increased. Therefore, the height-dependence of ionization degrees of hydrogen and helium may influence the damping rate of Alfven waves. We aim to study the effect of neutral helium in the damping of Alfven waves in stratified partially ionized plasma of the solar chromosphere. We consider a magnetic flux tube, which is expanded up to 1000 km height and then becomes vertical due to merging with neighboring tubes, and study the dynamics of linear torsional Alfven waves in the presence of neutral hydrogen and neutral helium atoms. We start with three-fluid description of plasma and consequently derive single-fluid magnetohydrodynamic (MHD) equations for torsional Alfven waves. Thin flux tube approximation allows to obtain the dispersion relation of the waves in the lower part of tubes, while the spatial dependence of steady-state Alfven waves is governed by Bessel type equation in the upper part of tubes. Consecutive derivation of single-fluid MHD equations results in a new Cowling diffusion coefficient in the presence of neutral helium which is different from previously used one. We found that shorter-period (< 5 s) torsional Alfven waves damp quickly in the chromospheric network due to ion-neutral collision. On the other hand, longer-period (> 5 s) waves do not reach the transition region as they become evanescent at lower heights in the network cores. Propagation of torsional Alfven waves through the chromosphere into the solar corona should be considered with caution: low-frequency waves are evanescent due to the stratification, while high-frequency waves are damped due to ion neutral collisions.

Alfv\'en waves are a particular class of magnetohydrodynamic waves relevant in many astrophysical and laboratory plasmas. In partially ionized plasmas the dynamics of Alfv\'en waves is affected by the interaction between ionized and neutral species. Here we study Alfv\'en waves in a partially ionized plasma from the theoretical point of view using the two-fluid description. We consider that the plasma is composed of an ion-electron fluid and a neutral fluid, which interact by means of particle collisions. To keep our investigation as general as possible we take the neutral-ion collision frequency and the ionization degree as free parameters. First, we perform a normal mode analysis. We find the modification due to neutral-ion collisions of the wave frequencies and study the temporal and spatial attenuation of the waves. In addition, we discuss the presence of cut-off values of the wavelength that constrain the existence of oscillatory standing waves in weakly ionized plasmas. Later, we go beyond the normal mode approach and solve the initial-value problem in order to study the time-dependent evolution of the wave perturbations in the two fluids. An application to Alfv\'en waves in the low solar atmospheric plasma is performed and the implication of partial ionization for the energy flux is discussed.

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 investigate theoretically the impact of partial ionization on MHD wave propagation in cylindrical flux tubes in the 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 take the neutral-ion collision frequency as an arbitrary parameter. Particular 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 realistic of the solar atmosphere ions and neutrals behave as a single fluid with an effective density corresponding to the sum of densities of both fluids and 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 taken into account for a consistent description of slow magnetoacoustic waves. The MHD waves are damped due to neutral-ion collisions. The damping is most efficient when the wave frequency and the collision frequency are of 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.

Observations show that propagating magnetohydrodynamic (MHD) waves are ubiquitous in the solar atmosphere. The technique of MHD seismology uses the wave observations combined with MHD wave theory to indirectly infer physical parameters of the solar atmospheric plasma and magnetic field. Here we present an analytical seismological inversion scheme for propagating MHD waves. This scheme uses in a consistent manner the observational information on wavelengths and damping lengths, along with observed values of periods or phase velocities, and is based on approximate asymptotic expressions for the theoretical values of wavelengths and damping lengths. The applicability of the inversion scheme is discussed and an example is given.

Transverse oscillations of thin threads in solar prominences are frequently reported in high-resolution observations. Two typical features of the observations are that the oscillations are damped in time and that simultaneous mass flows along the threads are detected. Flows cause the dense threads to move along the prominence magnetic structure while the threads are oscillating. The oscillations have been interpreted in terms of standing magnetohydrodynamic (MHD) kink waves of the magnetic flux tubes which support the threads. The damping is most likely due to resonant absorption caused by plasma inhomogeneity. The technique of seismology uses the observations combined with MHD wave theory to estimate prominence physical parameters. This paper presents a theoretical study of the joint effect of flow and resonant absorption on the amplitude of standing kink waves in prominence threads. We find that flow and resonant absorption can either be competing effects on the amplitude or both can contribute to damp the oscillations depending on the instantaneous position of the thread within the prominence magnetic structure. The amplitude profile deviates from the classic exponential profile of resonantly damped kink waves in static flux tubes. Flow also introduces a progressive shift of the oscillation period compared to the static case, although this effect is in general of minor importance. We test the robustness of seismological estimates by using synthetic data aiming to mimic real observations. The effect of the thread flow can significantly affect the estimation of the transverse inhomogeneity length scale. The presence of random background noise adds uncertainty to this estimation. Caution needs to be paid to the seismological estimates that do not take the influence of flow into account.

Magnetohydrodynamic (MHD) waves are ubiquitous in the solar atmosphere. Alfven waves and magneto-sonic waves are particular classes of MHD waves. These wave modes are clearly different and have pure properties in uniform plasmas of infinite extent only. Due to plasma non-uniformity MHD waves have mixed properties and cannot be classified as pure Alfven or magneto-sonic waves. However, vorticity is a quantity unequivocally related to Alfven waves as compression is for magneto-sonic waves. Here, we investigate MHD waves superimposed on a one-dimensional non-uniform straight cylinder with constant magnetic field. For a piecewise constant density profile we find that the fundamental radial modes of the non-axisymmetric waves have the same properties as surface Alfven waves at a true discontinuity in density. Contrary to the classic Alfveen waves in a uniform plasma of infinite extent, vorticity is zero everywhere except at the cylinder boundary. If the discontinuity in density is replaced with a continuous variation of density, vorticity is spread out over the whole interval with non-uniform density. The fundamental radial modes of the non-axisymmetric waves do not need compression to exist unlike the radial overtones. In thin magnetic cylinders the fundamental radial modes of the non-axisymmetric waves with phase velocities between the internal and the external Alfven velocities can be considered as surface Alfven waves. On the contrary, the radial overtones can be related to fast-like magneto-sonic modes.

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 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 the value of the collision frequency. In the compressible case, the domain of instability depends strongly on the plasma parameters, specially 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.

Context: Magnetohydrodynamic thermal modes may play an important role in the formation, plasma condensation, and evolution of solar prominences. Unstable thermal modes due to unbalance between radiative losses and heating can lead to rapid plasma cooling and condensation. An accurate description of the radiative loss function is therefore crucial for this process. Aims: We study the stability of thermal modes in unbounded and uniform plasmas with properties akin to those in solar prominences. Effects due to partial ionization are taken into account. Three different parametrizations of the radiative loss function are used. Methods: By means of a normal mode analysis, we investigate linear nonadiabatic perturbations superimposed on the equilibrium state. We find an approximate instability criterion for thermal modes, while the exact linear growth rate is obtained by numerically solving the general dispersion relation. The stability of thermal disturbances is compared for the three different loss functions considered. Results: Using up-to-date computations of radiative losses derived from the CHIANTI atomic database, we find that thermal modes may be unstable in prominences for lower temperatures than those predicted with previously existing loss functions. Thermal instability can take place for temperatures as low as 15,000 K, approximately. The obtained linear growth rates indicate that this instability might have an important impact on the dynamics and evolution of cool prominence condensations.

Given the difficulty in directly determining prominence physical parameters from observations, prominence seismology stands as an alternative method to probe the nature of these structures. We show recent examples of the application of magnetohydrodynamic (MHD) seismology techniques to infer physical parameters in prominence plasmas. They are based on the application of inversion techniques using observed periods, damping times, and plasma flow speeds of prominence thread oscillations. The contribution of Hinode to the subject has been of central importance. We show an example based on data obtained with Hinode's Solar Optical Telescope. Observations show an active region limb prominence, composed by a myriad of thin horizontal threads that flow following a path parallel to the photosphere and display synchronous vertical oscillations. The coexistence of waves and flows can be firmly established. By making use of an interpretation based on transverse MHD kink oscillations, a seismological analysis of this event is performed. It is shown that the combination of high quality Hinode observations and proper theoretical models allows flows and waves to become two useful characteristics for our understanding of the nature of solar prominences.

Prominences or filaments are cool clouds of partially ionized plasma living in the solar corona. Ground- and space-based observations have confirmed the presence of oscillatory motions in prominences and they have been interpreted in terms of magnetohydrodynamic (MHD) waves. Existing observational evidence points out that these oscillatory motions are damped in short spatial and temporal scales by some still not well known physical mechanism(s). Since prominences are partially ionized plasmas, a potential mechanism able to damp these oscillations could be ion-neutral collisions. Here, we will review the work done on the effects of partial ionization on MHD waves in prominence plasmas.

Context. Magnetohydrodynamic (MHD) waves are ubiquitous in the solar atmosphere. In magnetic waveguides resonant absorption due to plasma inhomogeneity naturally transfers wave energy from large-scale motions to small-scale motions. In the cooler parts of the solar atmosphere as, e.g., the chromosphere, effects due to partial ionization may be relevant for wave dynamics and heating. Aims. We study resonant Alfven waves in partially ionized plasmas. Methods. We use the multifluid equations in the cold plasma approximation. We investigate propagating resonant MHD waves in partially ionized flux tubes. We use approximate analytical theory based on normal modes in the thin tube and thin boundary approximations along with numerical eigenvalue computations. Results. We find that the jumps of the wave perturbations across the resonant layer are the same as in fully ionized plasmas. The damping length due to resonant absorption is inversely proportional to the frequency, while that due to ion-neutral collisions is inversely proportional to the square of the frequency. For observed frequencies in the solar atmosphere, the amplitude of MHD kink waves is more efficiently damped by resonant absorption than by ion-neutral collisions. Conclusions. Most of the energy carried by chromospheric kink waves is converted into localized azimuthal Alfven waves that can deposit energy in the coronal medium. The dissipation of wave energy in the chromosphere due to ion-neutral collisions is only effective for high-frequency waves. The chromosphere acts as a filter for kink waves with periods shorter than 10 s.

We analyse the oscillatory properties of resonantly damped transverse kink oscillations in two-dimensional prominence threads. The fine structures are modelled as cylindrically symmetric magnetic flux tubes with a dense central part with prominence plasma properties and an evacuated part, both surrounded by coronal plasma. The equilibrium density is allowed to vary non-uniformly in both the transverse and the longitudinal directions.We examine the influence of longitudinal density structuring on periods, damping times, and damping rates for transverse kink modes computed by numerically solving the linear resistive magnetohydrodynamic (MHD) equations. The relevant parameters are the length of the thread and the density in the evacuated part of the tube, two quantities that are difficult to directly estimate from observations. We find that both of them strongly influence the oscillatory periods and damping times, and to a lesser extent the damping ratios. The analysis of the spatial distribution of perturbations and of the energy flux into the resonances allows us to explain the obtained damping times. Implications for prominence seismology, the physics of resonantly damped kink modes in two-dimensional magnetic flux tubes, and the heating of prominence plasmas are discussed.

Recent observations by Hinode/SOT show that MHD waves and mass flows are simultaneously present in the fine structure of solar prominences. We investigate standing kink magnetohydrodynamic (MHD) waves in flowing prominence threads from a theoretical point of view. We model a prominence fine structure as a cylindrical magnetic tube embedded in the solar corona with its ends line-tied in the photosphere. The magnetic cylinder is composed of a region with dense prominence plasma, which is flowing along the magnetic tube, whereas the rest of the flux tube is occupied by coronal plasma. We use the WKB approximation to obtain analytical expressions for the period and the amplitude of the fundamental mode as functions of the flow velocity. In addition, we solve the full problem numerically by means of time-dependent simulations. We find that both the period and the amplitude of the standing MHD waves vary in time as the prominence thread flows along the magnetic structure. The fundamental kink mode is a good description for the time-dependent evolution of the oscillations, and the analytical expressions in the WKB approximation are in agreement with the full numerical results. The presence of flow modifies the period of the oscillations with respect to the static case. However, for realistic flow velocities this effect might fall within the error bars of the observations. The variation of the amplitude due to the flow leads to apparent damping or amplification of the oscillations, which could modify the real rate of attenuation caused by an additional damping mechanism.

It has been shown that resonant absorption is a robust physical mechanism to explain the observed damping of magnetohydrodynamic (MHD) kink waves in the solar atmosphere due to naturally occurring plasma inhomogeneity in the direction transverse to the direction of the magnetic field. Theoretical studies of this damping mechanism were greatly inspired by the first observations of post-flare standing kink modes in coronal loops using the Transition Region And Coronal Explorer (TRACE). More recently, these studies have been extended to explain the attenuation of propagating coronal kink waves observed by the Coronal Multi-Channel Polarimeter (CoMP). In the present study, for the first time we investigate the properties of propagating kink waves in solar waveguides including the effects of both longitudinal and transverse plasma inhomogeneity. Importantly, it is found that the wavelength is only dependent on the longitudinal stratification and the amplitude is simply a product of the two effects. In light of these results the advancement of solar atmospheric magnetoseismology by exploiting high spatial/temporal resolution observations of propagating kink waves in magnetic waveguides to determine the length scales of the plasma inhomogeneity along and transverse to the direction of the magnetic field is discussed.

Observations show the ubiquitous presence of propagating magnetohydrodynamic (MHD) kink waves in the solar atmosphere. Waves and flows are often observed simultaneously. Due to plasma inhomogeneity in the perpendicular direction to the magnetic field, kink waves are spatially damped by resonant absorption. The presence of flow may affect the wave spatial damping. Here, we investigate the effect of longitudinal background flow on the propagation and spatial damping of resonant kink waves in transversely nonuniform magnetic flux tubes. We combine approximate analytical theory with numerical investigation. The analytical theory uses the thin tube (TT) and thin boundary (TB) approximations to obtain expressions for the wavelength and the damping length. Numerically, we verify the previously obtained analytical expressions by means of the full solution of the resistive MHD eigenvalue problem beyond the TT and TB approximations. We find that the backward and forward propagating waves have different wavelengths and are damped on length scales that are inversely proportional to the frequency as in the static case. However, the factor of proportionality depends on the characteristics of the flow, so that the damping length differs from its static analogue. For slow, sub-Alfvenic flows the backward propagating wave gets damped on a shorter length scale than in the absence of flow, while for the forward propagating wave the damping length is longer. The different properties of the waves depending on their direction of propagation with respect to the background flow may be detected by the observations and may be relevant for seismological applications.

The fine structure of solar prominences and filaments appears as thin and long threads in high-resolution images. In H-alpha observations of filaments, some threads can be observed for only 5 - 20 minutes before they seem to fade and eventually disappear, suggesting that these threads may have very short lifetimes. The presence of an instability might be the cause of this quick disappearance. Here, we study the thermal instability of prominence threads as an explanation of their sudden disappearance from H-alpha observations. We model a prominence thread as a magnetic tube with prominence conditions embedded in a coronal environment. We assume a variation of the physical properties in the transverse direction, so that the temperature and density continuously change from internal to external values in an inhomogeneous transitional layer representing the particular prominence-corona transition region (PCTR) of the thread. We use the nonadiabatic and resistive magnetohydrodynamic equations, which include terms due to thermal conduction parallel and perpendicular to the magnetic field, radiative losses, heating, and magnetic diffusion. We combine both analytical and numerical methods to study linear perturbations from the equilibrium state, focusing on unstable thermal solutions. We find that thermal modes are unstable in the PCTR for temperatures higher than 80,000 K, approximately. These modes are related to temperature disturbances that can lead to changes in the equilibrium due to rapid plasma heating or cooling. For typical prominence parameters, the instability time scale is of the order of a few minutes and is independent of the form of the temperature profile within the PCTR of the thread. This result indicates that thermal instability may play an important role for the short lifetimes of threads in the observations.

Transverse oscillations and propagating waves are frequently observed in threads of solar prominences/filaments and have been interpreted as kink magnetohydrodynamic (MHD) modes. We investigate the spatial damping of propagating kink MHD waves in transversely nonuniform and partially ionized prominence threads. Resonant absorption and ion-neutral collisions (Cowling's diffusion) are the damping mechanisms taken into account. The dispersion relation of resonant kink waves in a partially ionized magnetic flux tube is numerically solved by considering prominence conditions. Analytical expressions of the wavelength and damping length as functions of the kink mode frequency are obtained in the Thin Tube and Thin Boundary approximations. For typically reported periods of thread oscillations, resonant absorption is an efficient mechanism for the kink mode spatial damping, while ion-neutral collisions have a minor role. Cowling's diffusion dominates both the propagation and damping for periods much shorter than those observed. Resonant absorption may explain the observed spatial damping of kink waves in prominence threads. The transverse inhomogeneity length scale of the threads can be estimated by comparing the observed wavelengths and damping lengths with the theoretically predicted values. However, the ignorance of the form of the density profile in the transversely nonuniform layer introduces inaccuracies in the determination of the inhomogeneity length scale.

We investigate standing kink magnetohydrodynamic (MHD) oscillations in a prominence fine structure modeled as a straight and cylindrical magnetic tube only partially filled with the prominence material, and with its ends fixed at two rigid walls representing the solar photosphere. The prominence plasma is partially ionized and a transverse inhomogeneous transitional layer is included between the prominence thread and the coronal medium. Thus, ion-neutral collisions and resonant absorption are the considered damping mechanisms. Approximate analytical expressions of the period, the damping time, and their ratio are derived for the fundamental mode in the thin tube and thin boundary approximations. We find that the dominant damping mechanism is resonant absorption, which provides damping ratios in agreement with the observations, whereas ion-neutral collisions are irrelevant for the damping. The values of the damping ratio are independent of both the prominence thread length and its position within the magnetic tube, and coincide with the values for a tube fully filled with the prominence plasma. The implications of our results in the context of the MHD seismology technique are discussed, pointing out that the reported short-period (2 - 10 min) and short-wavelength (700 - 8,000 km) thread oscillations may not be consistent with a standing mode interpretation and could be related to propagating waves. Finally, we show that the inversion of some prominence physical parameters, e.g., Alfv\'en speed, magnetic field strength, transverse inhomogeneity length-scale, etc., is possible using observationally determined values of the period and damping time of the oscillations along with the analytical approximations of these quantities.

Transverse oscillations of coronal loops are often observed and have been theoretically interpreted as kink magnetohydrodynamic (MHD) modes. Numerical simulations by Terradas et al. (2008, ApJ 687, L115) suggest that shear flows generated at the loop boundary during kink oscillations could give rise to a Kelvin-Helmholtz instability (KHI). Here, we investigate the linear stage of the KHI in a cylindrical magnetic flux tube in the presence of azimuthal shear motions. We consider the basic, linearized MHD equations in the beta = 0 approximation, and apply them to a straight and homogeneous cylindrical flux tube model embedded in a coronal environment. Azimuthal shear flows with a sharp jump of the velocity at the cylinder boundary are included in the model. We obtain an analytical expression for the dispersion relation of the unstable MHD modes supported by the configuration, and compute analytical approximations of the critical velocity shear and the KHI growth rate in the thin tube limit. A parametric study of the KHI growth rates is performed by numerically solving the full dispersion relation. We find that fluting-like modes can develop a KHI in time-scales comparable to the period of kink oscillations of the flux tube. The KHI growth rates increase with the value of the azimuthal wavenumber and decrease with the longitudinal wavenumber. However, the presence of a small azimuthal component of the magnetic field can suppress the KHI. Azimuthal motions related to kink oscillations of untwisted coronal loops may trigger a KHI, but this phenomenon has not been observed to date. We propose that the azimuthal component of the magnetic field is responsible for suppressing the KHI in a stable coronal loop. The required twist is small enough to prevent the development of the pinch instability.

Prominences are partially ionized, magnetized plasmas embedded in the solar corona. Damped oscillations and propagating waves are commonly observed. These oscillations have been interpreted in terms of magnetohydrodynamic (MHD) waves. Ion-neutral collisions and non-adiabatic effects (radiation losses and thermal conduction) have been proposed as damping mechanisms. We study the effect of the presence of helium on the time damping of non-adiabatic MHD waves in a plasma composed by electrons, protons, neutral hydrogen, neutral helium (He I), and singly ionized helium (He II) in the single-fluid approximation. The dispersion relation of linear non-adiabatic MHD waves in a homogeneous, unbounded, and partially ionized prominence medium is derived. The period and the damping time of Alfven, slow, fast, and thermal waves are computed. A parametric study of the ratio of the damping time to the period with respect to the helium abundance is performed. The efficiency of ion-neutral collisions as well as thermal conduction is increased by the presence of helium. However, if realistic abundances of helium in prominences (~10%) are considered, this effect has a minor influence on the wave damping. The presence of helium can be safely neglected in studies of MHD waves in partially ionized prominence plasmas.

Transverse oscillations of solar filament and prominence threads have been frequently reported. These oscillations have the common features of being of short period (2-10 min) and being damped after a few periods. Kink magnetohydrodynamic (MHD) wave modes have been proposed as responsible for the observed oscillations, whereas resonant absorption in the Alfven continuum and ion-neutral collisions are the best candidates to be the damping mechanisms. Here, we study both analytically and numerically the time damping of kink MHD waves in a cylindrical, partially ionized filament thread embedded in a coronal environment. The thread model is composed of a straight and thin, homogeneous filament plasma, with a transverse inhomogeneous transitional layer where the plasma physical properties vary continuously from filament to coronal conditions. The magnetic field is homogeneous and parallel to the thread axis. We find that the kink mode is efficiently damped by resonant absorption for typical wavelengths of filament oscillations, the damping times being compatible with the observations. Partial ionization does not affect the process of resonant absorption, and the filament plasma ionization degree is only important for the damping for wavelengths much shorter than those observed. To our knowledge, this is the first time that the phenomenon of resonant absorption is studied in a partially ionized plasma.

From recent high resolution observations obtained with the Swedish 1 m Solar Telescope in La Palma, we detect swaying motions of individual filament threads in the plane of the sky. The oscillatory character of these motions are comparable with oscillatory Doppler signals obtained from corresponding filament threads. Simultaneous recordings of motions in the line of sight and in the plane of the sky give information about the orientation of the oscillatory plane. These oscillations are interpreted in the context of the magnetohydrodynamic theory. Kink magnetohydrodynamic waves supported by the thread body are proposed as an explanation of the observed thread oscillations. On the basis of this interpretation and by means of seismological arguments, we give an estimation of the thread Alfv\'en speed and magnetic field strength by means of seismological arguments.

Oscillations and propagating waves are commonly seen in high-resolution observations of filament threads, i.e., the fine-structures of solar filaments/prominences. Since the temperature of prominences is typically of the order of 10^4 K, the prominence plasma is only partially ionized. In this paper, we study the effect of neutrals on the wave propagation in a filament thread modeled as a partially ionized homogeneous magnetic flux tube embedded in an homogeneous and fully ionized coronal plasma. Ohmic and ambipolar magnetic diffusion are considered in the basic resistive MHD equations. We numerically compute the eigenfrequencies of kink, slow, and Alfven linear MHD modes, and obtain analytical approximations in some cases. We find that the existence of propagating modes is constrained by the presence of critical values of the longitudinal wavenumber. In particular, the lower and upper frequency cut-offs of kink and Alfven waves owe their existence to magnetic diffusion parallel and perpendicular to magnetic field lines, respectively. The slow mode only has a lower frequency cut-off, which is caused by perpendicular magnetic diffusion and is significantly affected by the ionization degree. In addition, ion-neutral collisions is the most efficient damping mechanism for short wavelengths while ohmic diffusion dominates in the long-wavelength regime.

Transverse oscillations of small amplitude are commonly seen in high-resolution observations of filament threads, i.e. the fine-structures of solar filaments/prominences, and are typically damped in a few periods. Kink wave modes supported by the thread body offer a consistent explanation of these observed oscillations. Among the proposed mechanisms to explain the kink mode damping, resonant absorption in the Alfven continuum seems to be the most efficient as it produces damping times of about 3 periods. However, for a nonzero-beta plasma and typical prominence conditions, the kink mode is also resonantly coupled to slow (or cusp) continuum modes, which could further reduce the damping time. In this Letter, we explore for the first time both analytically and numerically the effect of the slow continuum on the damping of transverse thread oscillations. The thread model is composed of a homogeneous and straight cylindrical plasma, an inhomogeneous transitional layer, and the homogeneous coronal plasma. We find that the damping of the kink mode due to the slow resonance is much less efficient than that due to the Alfven resonance.

High resolution observations of solar filaments suggest the presence of groups of prominence threads, i.e. the fine-structures of prominences, which oscillate coherently (in phase). In addition, mass flows along threads have been often observed. Here, we investigate the effect of mass flows on the collective fast and slow nonadiabatic magnetoacoustic wave modes supported by systems of prominence threads. Prominence fine-structures are modeled as parallel, homogeneous and infinite cylinders embedded in a coronal environment. The magnetic field is uniform and parallel to the axis of threads. Configurations of identical and nonidentical threads are both explored. We apply the T-matrix theory of acoustic scattering to obtain the oscillatory frequency and the eigenfunctions of linear magnetosonic disturbances. We find that the existence of wave modes with a collective dynamics, i.e. those that produce significant perturbations in all threads, is only possible when the Doppler-shifted individual frequencies of threads are very similar. This can be only achieved for very particular values of the plasma physical conditions and flow velocities within threads.

Small-amplitude prominence oscillations are usually damped after a few periods. We study the attenuation of non-adiabatic magnetoacoustic waves in a slab prominence embedded in the coronal medium. We assume an equilibrium configuration with a transverse magnetic field to the slab axis and investigate wave damping by thermal conduction and radiative losses. The differential MHD equations that govern linear slow and fast modes are numerically solved to obtain the complex oscillatory frequency and the corresponding eigenfunctions. We find that coronal thermal conduction and radiative losses from the prominence plasma reveal as the most relevant damping mechanisms. Both mechanisms govern together the attenuation of hybrid modes, whereas prominence radiation is responsible for the damping of internal modes and coronal conduction essentially dominates the attenuation of external modes. In addition, the energy transfer between the prominence and the corona caused by thermal conduction has a noticeable effect on the wave stability, radiative losses from the prominence plasma being of paramount importance for the thermal stability of fast modes. We conclude that slow modes are efficiently damped, with damping times compatible with observations. On the contrary, fast modes are less attenuated by non-adiabatic effects and their damping times are several orders of magnitude larger than those observed. The presence of the corona causes a decrease of the damping times with respect to those of an isolated prominence slab, but its effect is still insufficient to obtain damping times of the order of the period in the case of fast modes.

High-resolution observations show that oscillations and waves in prominence threads are common and that they are attenuated in a few periods. In addition, observers have also reported the presence of material flows in such prominence fine-structures. Here we investigate the time damping of non-leaky oscillations supported by a homogeneous cylindrical prominence thread embedded in an unbounded corona and with a steady mass flow. Thermal conduction and radiative losses are taken into account as damping mechanisms, and the effect of these non-ideal effects and the steady flow on the attenuation of oscillations is assessed. We solve the general dispersion relation for linear, non-adiabatic magnetoacoustic and thermal waves supported by the model, and find that slow and thermal modes are efficiently attenuated by non-adiabatic mechanisms. On the contrary, fast kink modes are much less affected and their damping times are much larger than those observed. The presence of flow has no effect on the damping of slow and thermal waves, whereas fast kink waves are more (less) attenuated when they propagate parallel (anti-parallel) to the flow direction. Although the presence of steady mass flows improves the efficiency of non-adiabatic mechanisms on the attenuation of transverse, kink oscillations for parallel propagation to the flow, its effect is still not enough to obtain damping times compatible with observations.

Oscillations of magnetic structures in the solar corona have often been interpreted in terms of magnetohydrodynamic waves. We study the adiabatic magnetoacoustic modes of a prominence plasma slab with a uniform longitudinal magnetic field, surrounded by a prominence-corona transition region (PCTR) and a coronal medium. Considering linear small-amplitude oscillations, the dispersion relation for the magnetoacoustic slow and fast modes is deduced assuming evanescent-like perturbations in the coronal medium. In the system without PCTR, a classification of the oscillatory modes according to the polarisation of their eigenfunctions is made in order to distinguish modes with fast-like or slow-like properties. Internal and external slow modes are governed by the prominence and coronal properties respectively, and fast modes are mostly dominated by prominence conditions for the observed wavelengths. In addition, the inclusion of an isothermal PCTR does not substantially influence the mode frequencies, but new solutions (PCTR slow modes) are present.

Context. One of the typical features shown by observations of solar prominence oscillations is that they are damped in time and that the values of the damping times are usually between one and three times the corresponding oscillatory period. However, the mechanism responsible for the attenuation is still not well-known. Aims. Thermal conduction, optically thin or thick radiation and heating are taken into account in the energy equation, and their role on the attenuation of prominence oscillations is evaluated. Methods. The dispersion relation for linear non-adiabatic magnetoacoustic waves is derived considering an equilibrium made of a prominence plasma slab embedded in an unbounded corona. The magnetic field is orientated along the direction parallel to the slab axis and has the same strength in all regions. By solving the dispersion relation for a fixed wavenumber, a complex oscillatory frequency is obtained, and the period and the damping time are computed. Results. The effect of conduction and radiation losses is different for each magnetoacoustic mode and depends on the wavenumber. In the observed range of wavelengths the internal slow mode is attenuated by radiation from the prominence plasma, the fast mode by the combination of prominence radiation and coronal conduction and the external slow mode by coronal conduction. The consideration of the external corona is of paramount importance in the case of the fast and external slow modes, whereas it does not affect the internal slow modes at all. Conclusions. Non-adiabatic effects are efficient damping mechanisms for magnetoacoustic modes, and the values of the obtained damping times are compatible with those observed.