Context. Prominences show a surprising amount of fine structure and it is widely believed that their threads, as seen in Hα observations, provide indirect information concerning magnetic field topology. Both prominence and coronal rain condensations most likely originate from thermal instabilities in the solar corona. It is still not understood how non-linear instability evolution shapes the observed fine structure of prominences. Investigating this requires multidimensional, high-resolution simulations to resolve all emerging substructure in great detail. Aims. We investigate the spontaneous emergence and evolution of fine structure in high-density condensations formed through the process of thermal instability under typical solar coronal conditions. Our study reveals intricate multidimensional processes that occur through in situ condensations in a representative coronal volume in a low plasma beta regime. Methods. We quantified slow wave eigenfunctions used as perturbations and discuss under which conditions the thermal mode is unstable when anisotropic thermal conduction effects are included. We performed 2D and 3D numerical simulations of interacting slow magnetohydrodynamic (MHD) wave modes when all relevant non-adiabatic effects are included. Multiple levels of adaptive mesh refinement achieve a high resolution near regions with high density, thereby resolving any emerging fine structure automatically. Our study employs a local periodic coronal region traversed by damped slow waves inspired by the presence of such waves observed in actual coronal magnetic structures. Results. We show that the interaction of multiple slow MHD wave modes in a regime unstable to the thermal mode leads to thermal instability. This initially forms pancake-like structures almost orthogonal to the local magnetic field, while low-pressure induced inflows of matter generate rebound shocks. This is succeeded by the rapid disruption of these pancake sheets through thin-shell instabilities evolving naturally from minute ram pressure imbalances. This eventually creates high-density blobs accompanied by thread-like features from shear flow effects. The further evolution of the blobs follows the magnetic field lines, such that a dynamical realignment with the background magnetic field appears. However, the emerging thread-like features are not at all field-aligned, implying only a very weak link between fine structure orientation and magnetic field topology. Conclusions. As seen in our synthetic Hα views, threads formed by non-linear thermal instability evolution do not strictly outline magnetic field structure and this finding has far-reaching implications for field topology interpretations based on Hα observations.

Context. Thermal instabilities give rise to condensations in the solar corona, and are the most probable scenario for coronal rain and prominence formation. We revisit the original theoretical treatment done by Field (1965, ApJ, 142, 531) in a homogeneous plasma with heat-loss effects and combine this with state-of-the-art numerical simulations to verify growth-rate predictions and address the long-term non-linear regime. We especially investigate interaction between multiple magnetohydrodynamic (MHD) wave modes and how they in turn trigger thermal mode development. Aims. We assess how well the numerical MHD simulations retrieve the analytically predicted growth rates. We complete the original theory with quantifications of the eigenfunctions, calculated to consistently excite each wave mode. Thermal growth rates are fitted also in the non-linear regime of multiple wave-wave interaction setups, at the onset of thermal instability formation. Methods. We performed 2D numerical MHD simulations, including an additional (radiative) heat-loss term and a constant heating term to the energy equation. We mainly focus on the thermal (i.e. entropy) and slow MHD wave modes and use the wave amplitude as a function of time to make a comparison to predicted growth rates. Results. It is shown that the numerical MHD simulations retrieve analytically predicted growth rates for all modes, of thermal and slow or fast MHD type. In typical coronal conditions, the latter are damped due to radiative losses, but their interaction can cause slowly changing equilibrium conditions which ultimately trigger thermal mode development. Even in these non-linear wave-wave interaction setups, the growth rate of the thermal instability agrees with the exponential profile predicted by linear theory. The nonlinear evolutions show systematic field-guided motions of condensations with fairly complex morphologies, resulting from thermal modes excited through damped slow MHD waves. These results are of direct interest to the study of solar coronal rain and prominence fine structure. Our wave-wave interaction setups are relevant for coronal loop sections which are known to host slow wave modes, and hence provide a new route to explain the sudden onset of coronal condensation.

Magnetohydrodynamic (MHD) spectroscopy is central to many astrophysical disciplines, ranging from helio-to asteroseismology, over solar coronal (loop) seismology, to the study of waves and instabilities in jets, accretion disks, or solar/stellar atmospheres. MHD spectroscopy quantifies all linear (standing or travelling) wave modes, including overstable (i.e. growing) or damped modes, for a given configuration that achieves force and thermodynamic balance. Here, we present Legolas a) , a novel, open-source numerical code to calculate the full MHD spectrum of one-dimensional equilibria with flow, that balance pressure gradients, Lorentz forces, centrifugal effects and gravity, enriched with non-adiabatic aspects like radiative losses, thermal conduction and resistivity. The governing equations use Fourier representations in the ignorable coordinates, and the set of linearised equations are discretised using Finite Elements in the important height or radial variation, handling Cartesian and cylindrical geometries using the same implementation. A weak Galerkin formulation results in a generalised (non-Hermitian) matrix eigenvalue problem, and linear algebraic algorithms calculate all eigenvalues and corresponding eigenvectors. We showcase a plethora of well-established results, ranging from p-and g-modes in magnetised, stratified atmospheres, over modes relevant for coronal loop seismology, thermal instabilities and discrete overstable Alfvén modes related to solar prominences, to stability studies for astrophysical jet flows. We encounter (quasi-)Parker, (quasi-)interchange, currentdriven and Kelvin-Helmholtz instabilities, as well as non-ideal quasi-modes, resistive tearing modes, up to magneto-thermal instabilities. The use of high resolution sheds new light on previously calculated spectra, revealing interesting spectral regions that have yet to be investigated.

The quantification of all possible waves and instabilities in a given system is of paramount importance, and knowledge of the full magnetohydrodynamic (MHD) spectrum allows one to predict the (in)stability of a given equilibrium state. This is highly relevant in many (astro)physical disciplines, and when applied to the solar atmosphere it may yield various new insights in processes like prominence formation and coronal loop oscillations. In this work we present a detailed, high-resolution spectroscopic study of the solar atmosphere, where we use our newly developed Legolas code to calculate the full spectrum with corresponding eigenfunctions of equilibrium configurations that are based on fully realistic solar atmospheric models, including gravity, optically thin radiative losses and thermal conduction. Special attention is given to thermal instabilities, known to be responsible for the formation of prominences, together with a new outlook on the thermal and slow continua and how they behave in different chromospheric and coronal regions. We show that thermal instabilities are unavoidable in our solar atmospheric models and that there exist certain regions where both the thermal, slow and fast modes all have unstable wave mode solutions. We also encounter regions where the slow and thermal continua become purely imaginary and merge on the imaginary axis. The spectra discussed in this work illustrate clearly that thermal instabilities (both discrete and continuum modes) and magneto-thermal overstable propagating modes are ubiquitous throughout the solar atmosphere, and may well be responsible for much of the observed fine-structuring and multi-thermal dynamics.

Many linear stability aspects in plasmas are heavily influenced by non-ideal effects beyond the basic ideal magnetohydrodynamics (MHD) description. Here, the extension of the modern open-source MHD spectroscopy code Legolas with viscosity and the Hall current is highlighted and benchmarked on a stringent set of historic and recent findings. The viscosity extension is demonstrated in a cylindrical set-up featuring Taylor–Couette flow and in a viscoresistive plasma slab with a tearing mode. For the Hall extension, we show how the full eigenmode spectrum relates to the analytic dispersion relation in an infinite homogeneous medium. We quantify the Hall term influence on the resistive tearing mode in a Harris current sheet, including the effect of compressibility, which is absent in earlier studies. Furthermore, we illustrate how Legolas mimics the incompressible limit easily to compare with literature results. Going beyond published findings, we emphasise the importance of computing the full eigenmode spectrum, and how elements of the spectrum are modified by compressibility. These extensions allow for future stability studies with Legolas that are relevant to ongoing dynamo experiments, protoplanetary disks or magnetic reconnection.