Seismology of Standing Kink Oscillations of Solar Prominence Fine Structures

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.Comment: Accepted for publication in Ap

Nonlinear Instability of kink oscillations due to shear motions

First results from a high-resolution three-dimensional nonlinear numerical study of the kink oscillation are presented. We show in detail the development of a shear instability in an untwisted line-tied magnetic flux tube. The instability produces significant deformations of the tube boundary. An extended transition layer may naturally evolve as a result of the shear instability at a sharp transition between the flux tube and the external medium. We also discuss the possible effects of the instability on the process of resonant absorption when an inhomogeneous layer is included in the model. One of the implications of these results is that the azimuthal component of the magnetic field of a stable flux tube in the solar corona, needed to prevent the shear instability, is probably constrained to be in a very specific range

Three-Dimensional Propagation of Magnetohydrodynamic Waves in Solar Coronal Arcades

We numerically investigate the excitation and temporal evolution of oscillations in a two-dimensional coronal arcade by including the three-dimensional propagation of perturbations. The time evolution of impulsively generated perturbations is studied by solving the linear, ideal magnetohydrodynamic (MHD) equations in the zero-beta approximation. As we neglect gas pressure the slow mode is absent and therefore only coupled MHD fast and Alfven modes remain. Two types of numerical experiments are performed. First, the resonant wave energy transfer between a fast normal mode of the system and local Alfven waves is analyzed. It is seen how, because of resonant coupling, the fast wave with global character transfers its energy to Alfvenic oscillations localized around a particular magnetic surface within the arcade, thus producing the damping of the initial fast MHD mode. Second, the time evolution of a localized impulsive excitation, trying to mimic a nearby coronal disturbance, is considered. In this case, the generated fast wavefront leaves its energy on several magnetic surfaces within the arcade. The system is therefore able to trap energy in the form of Alfvenic oscillations, even in the absence of a density enhancement such as that of a coronal loop. These local oscillations are subsequently phase-mixed to smaller spatial scales. The amount of wave energy trapped by the system via wave energy conversion strongly depends on the wavelength of perturbations in the perpendicular direction, but is almost independent from the ratio of the magnetic to density scale heights.Comment: 27 pages, 11 figure

Magnetohydrodynamic kink waves in two-dimensional non-uniform prominence threads

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.Comment: 12 pages, 9 figures, A&A accepte

The role of Rayleigh-Taylor instabilities in filament threads

Many solar filaments and prominences show short-lived horizontal threads lying parallel to the photosphere. In this work the possible link between Rayleigh-Taylor instabilities and thread lifetimes is investigated. This is done by calculating the eigenmodes of a thread modelled as a Cartesian slab under the presence of gravity. An analytical dispersion relation is derived using the incompressible assumption for the magnetohydrodynamic (MHD) perturbations. The system allows a mode that is always stable, independently of the value of the Alfv\'en speed in the thread. The character of this mode varies from being localised at the upper interface of the slab when the magnetic field is weak, to having a global nature and resembling the transverse kink mode when the magnetic field is strong. On the contrary, the slab model permits another mode that is unstable and localised at the lower interface when the magnetic field is weak. The growth rates of this mode can be very short, of the order of minutes for typical thread conditions. This Rayleigh-Taylor unstable mode becomes stable when the magnetic field is increased, and in the limit of strong magnetic field it is essentially a sausage magnetic mode. The gravity force might have a strong effect on the modes of oscillation of threads, depending on the value of the Alfv\'en speed. In the case of threads in quiescent filaments, where the Alfv\'en speed is presumably low, very short lifetimes are expected according to the slab model. In active region prominences, the stabilising effect of the magnetic tension might be enough to suppress the Rayleigh-Taylor instability for a wide range of wavelengths

On the nature of kink MHD waves in magnetic flux tubes

Magnetohydrodynamic (MHD) waves are often reported in the solar atmosphere and usually classified as slow, fast, or Alfv\'en. The possibility that these waves have mixed properties is often ignored. The goal of this work is to study and determine the nature of MHD kink waves. This is done by calculating the frequency, the damping rate and the eigenfunctions of MHD kink waves for three widely different MHD waves cases: a compressible pressure-less plasma, an incompressible plasma and a compressible plasma with non-zero plasma pressure which allows for MHD radiation. In all three cases the frequency and the damping rate are for practical purposes the same as they differ at most by terms proportional to $(k_z R)^2$. In the magnetic flux tube the kink waves are in all three cases, to a high degree of accuracy incompressible waves with negligible pressure perturbations and with mainly horizontal motions. The main restoring force of kink waves in the magnetised flux tube is the magnetic tension force. The total pressure gradient force cannot be neglected except when the frequency of the kink wave is equal or slightly differs from the local Alfv\'{e}n frequency, i.e. in the resonant layer. Kink waves are very robust and do not care about the details of the MHD wave environment. The adjective fast is not the correct adjective to characterise kink waves. If an adjective is to be used it should be Alfv\'{e}nic. However, it is better to realize that kink waves have mixed properties and cannot be put in one single box

Analytic approximate seismology of transversely oscillating coronal loops

We present an analytic approximate seismic inversion scheme for damped transverse coronal loop oscillations based on the thin tube and thin boundary approximation for computing the period and the damping time. Asymptotic expressions for the period and damping rate are used to illustrate the process of seismological inversion in a simple and easy to follow manner. The inversion procedure is formulated in terms of two simple functions, which are given by simple closed expressions. The analytic seismic inversion shows that an infinite amount of 1-dimensional equilibrium models can reproduce the observed periods and damping times. It predicts a specific range of allowable values for the Alfven travel time and lower bounds for the density contrast and the inhomogeneity length scale. When the results of the present analytic seismic inversion are compared with those of a previous numerical inversion, excellent agreement is found up to the point that the analytic seismic inversion emerges as a tool for validating results of numerical inversions. Actually it helped us to identify and correct inaccuracies in a previous numerical investigation.Comment: 7 pages, 1 figure, A&A, accepte

Mixed Properties of MHD Waves in Non-uniform Plasmas

This paper investigates the mixed properties of MHD waves in a non-uniform plasma. It starts with a short revision of MHD waves in a uniform plasma of infinite extent. In that case the MHD waves do not have mixed properties. They can be separated in Alfvén waves and magneto-sonic waves. The Alfvén waves propagate parallel vorticity and are incompressible. In addition they have no parallel displacement component. The magneto-sonic waves are compressible and in general do have a parallel component of displacement but do not propagate parallel vorticity. This clear separation has been the reason why there has been a strong inclination in the literature to use this classification in the study of MHD waves in non-uniform plasmas. The main part of this paper is concerned with MHD waves in a non-uniform plasma. It is shown that the MHD waves in that situation in general propagate both vorticity and compression and hence have mixed properties. Finally, the close connection between resonant absorption and MHD waves with mixed properties is discussed

MCMC-driven importance samplers

Monte Carlo sampling methods are the standard procedure for approximating complicated integrals of multidimensional posterior distributions in Bayesian inference. In this work, we focus on the class of layered adaptive importance sampling algorithms, which is a family of adaptive importance samplers where Markov chain Monte Carlo algorithms are employed to drive an underlying multiple importance sampling scheme. The modular nature of the layered adaptive importance sampling scheme allows for different possible implementations, yielding a variety of different performances and computational costs. In this work, we propose different enhancements of the classical layered adaptive importance sampling setting in order to increase the efficiency and reduce the computational cost, of both upper and lower layers. The different variants address computational challenges arising in real-world applications, for instance with highly concentrated posterior distributions. Furthermore, we introduce different strategies for designing cheaper schemes, for instance, recycling samples generated in the upper layer and using them in the final estimators in the lower layer. Different numerical experiments show the benefits of the proposed schemes, comparing with benchmark methods presented in the literature, and in several challenging scenarios