Showing 1–50 of 72 results for author: Soler, R

The influence of different coronal hole geometries on simulations of coronal wave -- coronal hole interaction

Authors: I. Piantschitsch, J. Terradas, E. Soubrie, S. G. Heinemann, S. J. Hofmeister, R. Soler, M. Temmer

Abstract: The geometry of a coronal hole (CH) affects the density profile of the reflected part of an incoming global coronal wave (CW). In this study, we perform for the first time magnetohydrodynamic (MHD) simulations of fast-mode MHD waves interacting with CHs of different geometries, such as circular, elliptic, convex, and concave shapes. We analyse the influence these geometries have on the density pro… ▽ More The geometry of a coronal hole (CH) affects the density profile of the reflected part of an incoming global coronal wave (CW). In this study, we perform for the first time magnetohydrodynamic (MHD) simulations of fast-mode MHD waves interacting with CHs of different geometries, such as circular, elliptic, convex, and concave shapes. We analyse the influence these geometries have on the density profiles of the reflected waves and we generate the corresponding simulation-based time-distance plots. Within these time-distance plots we determine regions that exhibit specific density features, such as large reflected density amplitudes. In a further step, these interaction features can be compared to actual observed CW-CH interaction events which makes it possible to explain interaction parameters from the observed interaction events, such as the density structure of the reflected wave, which are usually difficult to comprehensively understand by only analysing the measurements. Moreover, we show that the interaction between a concave shaped CH and CWs, whose density profile include an enhanced as well as a depleted wave, can lead to reflected density amplitudes that are more than two times larger than the incoming ones. Another effect of the interplay between the constructive and destructive interference of the reflected wave parts is a strongly depleted region in the middle of the CW-CH interaction process. In addition, we show how important the choice of the path is that is used to generate the time-distance plots and how this choice affects the interpretation of the CW-CH interaction results. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Comments: arXiv admin note: text overlap with arXiv:2308.08928

Numerical simulations of turbulence in prominence threads induced by torsional oscillations

Authors: Sergio Díaz-Suárez, Roberto Soler

Abstract: Threads are the main constituents of prominences and are subjected to oscillations that might be interpreted as MHD waves. Moreover, the Kelvin-Helmholtz instability (KHI) has been reported in prominences. Both waves and KHI may affect the thermodynamic state of the threads. We investigate the triggering of turbulence in a thread caused by the nonlinear evolution of standing torsional Alfvén waves… ▽ More Threads are the main constituents of prominences and are subjected to oscillations that might be interpreted as MHD waves. Moreover, the Kelvin-Helmholtz instability (KHI) has been reported in prominences. Both waves and KHI may affect the thermodynamic state of the threads. We investigate the triggering of turbulence in a thread caused by the nonlinear evolution of standing torsional Alfvén waves as well as possible observational signatures of this dynamics and the plasma heating. We modeled the thread as a radially and longitudinally nonuniform cylindrical flux tube with a uniform axial magnetic field embedded in a coronal environment. We perturbed the flux tube with the longitudinally fundamental mode of standing torsional Alfvén waves and numerically solved the 3D MHD equations to study the temporal evolution in both ideal and dissipative scenarios. We also performed forward modeling to calculate the synthetic Hα imaging. Standing torsional Alfvén waves undergo phase-mixing owing to the radially nonuniform density. The phase-mixing generates azimuthal shear flows that eventually trigger the KHI and, later, turbulence. If nonideal effects are included, plasma heating is localized in an annulus region at the thread boundary and does not increase the temperature in the cool core. Instead, the average temperature in the thread decreases owing to the mixing of internal and external plasmas. In the synthetic observations, first we find periodic pulsations in the Hα intensity caused by the integration of the phase-mixing flows along the line of sight. Later, we find fine strands that may be associated with the KHI vortices. Turbulence can be generated by standing torsional Alfvén waves in a thread after the onset of KHI, but this mechanism is not enough to heat globally the structure. The dynamics could be seen in high-resolution Hα observations. △ Less

Submitted 17 January, 2024; originally announced January 2024.

Comments: Accepted in Astronomy & Astrophysics

The role of initial density profiles in simulations of coronal wave - coronal hole interaction

Authors: I. Piantschitsch, J. Terradas, E. Soubrie, S. G. Heinemann, S. J. Hofmeister, R. Soler, M. Temmer

Abstract: Interactions between global coronal waves (CWs) and coronal holes (CHs) reveal many interesting features of reflected waves and coronal hole boundaries (CHB) but have fairly been studied so far. Magnetohydrodynamic (MHD) simulations can help us to better understand what is happening during these interaction events, and therefore, to achieve a broader understanding of the parameters involved. In th… ▽ More Interactions between global coronal waves (CWs) and coronal holes (CHs) reveal many interesting features of reflected waves and coronal hole boundaries (CHB) but have fairly been studied so far. Magnetohydrodynamic (MHD) simulations can help us to better understand what is happening during these interaction events, and therefore, to achieve a broader understanding of the parameters involved. In this study, we perform for the first time 2D MHD simulations of a CW-CH interaction including a realistic initial wave density profile that consists of an enhanced as well as a depleted wave part. We vary several initial parameters, such as the initial density amplitudes of the incoming wave, the CH density, and the CHB width, which are all based on actual measurements. We analyse the effects of different incident angles on the interaction features and we use the corresponding time-distance plots to detect specific features of the incoming and the reflected wave. We found that a particular combination of a small CH density, a realistic initial density profile and a sufficiently small incident angle leads to remarkable interaction features, such as a large density amplitude of the reflected wave with respect to the incoming one. The parameter studies in this paper provide a tool to compare time-distance plots based on observational measurements to those created from simulations and therefore enable us to derive interaction parameters from observed CW-CH interaction events that usually cannot be obtained directly. The simulation results in this study are augmented by analytical expressions for the reflection coefficient of the CW-CH interaction which allows us to verify the simulations results in an additional way. This work is the first of a series of studies aiming to finally reconstruct actual observed CW-CH interaction events by means of MHD-simulations. △ Less

Submitted 17 August, 2023; originally announced August 2023.

Self-consistent equilibrium models of prominence thin threads heated by Alfvén waves propagating from the photosphere

Authors: Llorenç Melis, Roberto Soler, Jaume Terradas

Abstract: The fine structure of solar prominences is made by thin threads that outline the magnetic field lines. Observations show that transverse waves of Alfvénic nature are ubiquitous in prominence threads. These waves are driven at the photosphere and propagate to prominences suspended in the corona. Heating due to Alfvén wave dissipation could be a relevant mechanism in the cool and partially ionized p… ▽ More The fine structure of solar prominences is made by thin threads that outline the magnetic field lines. Observations show that transverse waves of Alfvénic nature are ubiquitous in prominence threads. These waves are driven at the photosphere and propagate to prominences suspended in the corona. Heating due to Alfvén wave dissipation could be a relevant mechanism in the cool and partially ionized prominence plasma. We explore the construction of 1D equilibrium models of prominence thin threads that satisfy energy balance between radiative losses, thermal conduction, and Alfvén wave heating. We assume the presence of a broadband driver at the photosphere that launches Alfvén waves towards the prominence. An iterative method is implemented, in which the energy balance equation and the Alfvén wave equation are consecutively solved. From the energy balance equation and considering no wave heating initially, we compute the equilibrium profiles along the thread of the temperature, density, ionisation fraction. We use the Alfvén wave equation to compute the wave heating rate, which is then put back in the energy balance equation to obtain new equilibrium profiles. The process is repeated until convergence to a self-consistent thread model heated by Alfvén waves is achieved. We have obtained equilibrium models composed of a cold and dense thread, a extremely thin PCTR, and an extended coronal region. The length of the cold thread decreases with the temperature at the prominence core and increases with the Alfvén wave energy flux. Equilibrium models are not possible for sufficiently large wave energy fluxes when the wave heating rate inside the cold thread becomes larger than radiative losses. The maximum value of the wave energy flux that allows an equilibrium depends on the prominence core temperature. This constrains the existence of equilibria in realistic conditions. △ Less

Submitted 23 June, 2023; originally announced June 2023.

Comments: 12 pages, 11 figures

Journal ref: A&A 676, A25 (2023)

Transition to turbulence in nonuniform coronal loops driven by torsional Alfvén waves. II. Extended analysis and effect of magnetic twist

Authors: Sergio Díaz-Suárez, Roberto Soler

Abstract: It has been shown in a previous work that torsional Alfvén waves can drive turbulence in nonuniform coronal loops with a purely axial magnetic field. Here we explore the role of the magnetic twist. We model a coronal loop as a transversely nonuniform straight flux tube, anchored in the photosphere, and embedded in a uniform coronal environment. We consider that the magnetic field is twisted and co… ▽ More It has been shown in a previous work that torsional Alfvén waves can drive turbulence in nonuniform coronal loops with a purely axial magnetic field. Here we explore the role of the magnetic twist. We model a coronal loop as a transversely nonuniform straight flux tube, anchored in the photosphere, and embedded in a uniform coronal environment. We consider that the magnetic field is twisted and control the strength of magnetic twist by a free parameter of the model. We excite the longitudinally fundamental mode of standing torsional Alfvén waves, whose temporal evolution is obtained by means of high-resolution three-dimensional ideal magnetohydrodynamic numerical simulations. We find that phase mixing of torsional Alfvén waves creates velocity shear in the direction perpendicular to the magnetic field lines. The velocity shear eventually triggers the Kelvin-Helmholtz instability (KHi). In weakly twisted magnetic tubes, the KHi is able to grow nonlinearly and, subsequently, turbulence is driven in the coronal loop in a similar manner as in the untwisted case. Provided that magnetic twist remains weak, the effect of magnetic twist is to delay the onset of the KHi and to slow down the development of turbulence. In contrast, magnetic tension can suppress the nonlinear growth of the KHi when magnetic twist is strong enough, even if the KHi has locally been excited by the phase-mixing shear. Thus, turbulence is not generated in strongly twisted loops △ Less

Submitted 19 September, 2022; v1 submitted 13 July, 2022; originally announced July 2022.

Comments: Published in A&A. 12 pages. 12 figures

Journal ref: A&A 665, A113 (2022)

Quasimodes in the cusp continuum in nonuniform magnetic flux tubes

Authors: Michaël Geeraerts, Pieter Vanmechelen, Tom Van Doorsselaere, Roberto Soler

Abstract: The study of MHD waves is important both for understanding heating in the solar atmosphere and for solar atmospheric seismology. The analytical investigation of wave mode properties in a cylinder is of particular interest in this domain, as many atmospheric structures can be modeled as such in a first approximation. We use linearized ideal MHD to study quasimodes (global modes that are damped thro… ▽ More The study of MHD waves is important both for understanding heating in the solar atmosphere and for solar atmospheric seismology. The analytical investigation of wave mode properties in a cylinder is of particular interest in this domain, as many atmospheric structures can be modeled as such in a first approximation. We use linearized ideal MHD to study quasimodes (global modes that are damped through resonant absorption) with a frequency in the cusp continuum, in a straight cylinder with a circular base and an inhomogeneous layer at its boundary which separates two homogeneous plasma regions inside and outside. We are in particular interested in the damping of these modes, and shall hence try to determine their frequency as a function of background parameters. After linearizing the ideal MHD equations, we find solutions to the second-order differential equation for the perturbed total pressure in the inhomogeneous layer in the form of Frobenius series around the regular singular points that are the Alfvén and cusp resonant positions, as well as power series around regular points. By connecting these solutions appropriately through the inhomogeneous layer and with the solutions of the homogeneous regions inside and outside the cylinder, we derive a dispersion relation for the frequency of the eigenmodes of the system. From the dispersion relation, it is also possible to find the frequency of quasimodes even though they are not eigenmodes. As an example, we find the frequency of the slow surface sausage quasimode as a function of the inhomogeneous layer's width, for values of the longitudinal wavenumber relevant for photospheric conditions. The results were found to match well the results found in another paper which studied the resistive slow surface sausage eigenmode. We also discuss the perturbation profiles of the quasimode and the eigenfunctions of continuum modes. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Journal ref: A&A 661, A100 (2022)

Construction of coronal hole and active region magnetohydrostatic solutions in two dimensions: Force and energy balance

Authors: J. Terradas, R. Soler, R. Oliver, P. Antolin, I. Arregui, M. Luna, I. Piantschitsch, E. Soubrié, J. L., Ballester

Abstract: Coronal holes and active regions are typical magnetic structures found in the solar atmosphere. We propose several magnetohydrostatic equilibrium solutions that are representative of these structures in two dimensions. Our models include the effect of a finite plasma-$β$ and gravity, but the distinctive feature is that we incorporate a thermal structure with properties similar to those reported by… ▽ More Coronal holes and active regions are typical magnetic structures found in the solar atmosphere. We propose several magnetohydrostatic equilibrium solutions that are representative of these structures in two dimensions. Our models include the effect of a finite plasma-$β$ and gravity, but the distinctive feature is that we incorporate a thermal structure with properties similar to those reported by observations. We developed a semi-analytical method to compute the equilibrium configuration. Using this method, we obtain cold and under-dense plasma structures in open magnetic fields representing coronal holes, while in closed magnetic configurations, we achieve the characteristic hot and over-dense plasma arrangements of active regions. Although coronal holes and active regions seem to be antagonistic structures, we find that they can be described using a common thermal structure that depends on the flux function. In addition to the force balance, the energy balance is included in the constructed models using an a posteriori approach. From the two-dimensional computation of thermal conduction and radiative losses in our models, we infer the required heating function to achieve energy equilibrium. We find that the temperature dependence on height is an important parameter that may prevent the system from accomplishing thermal balance at certain spatial locations. The implications of these results are discussed in detail. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: Accepted for publication in A&A

Transverse kink oscillations of inhomogeneous prominence threads: numerical analysis and H$α$ forward modelling

Authors: David Martínez-Gómez, Roberto Soler, Jaume Terradas, Elena Khomenko

Abstract: Prominence threads are very long and thin flux tubes which are partially filled with cold plasma. Observations have shown that transverse oscillations are frequent in these solar structures. The observations are usually interpreted as the fundamental kink mode, while the detection of the first harmonic remains elusive. Here, we aim to study how the density inhomogeneity in the longitudinal and rad… ▽ More Prominence threads are very long and thin flux tubes which are partially filled with cold plasma. Observations have shown that transverse oscillations are frequent in these solar structures. The observations are usually interpreted as the fundamental kink mode, while the detection of the first harmonic remains elusive. Here, we aim to study how the density inhomogeneity in the longitudinal and radial directions modify the periods and damping times of kink oscillations, and how this effect would be reflected in observations. We solve the ideal magnetohydrodynamics equations through two different methods: a) performing 3D numerical simulations, and b) solving a 2D generalised eigenvalue problem. We study the dependence of the periods, damping times and amplitudes of transverse kink oscillations on the ratio between the densities at the centre and at the ends of the tube, and on the average density. We apply forward modelling on our 3D simulations to compute synthetic H$α$ profiles. We confirm that the ratio of the period of the fundamental oscillation mode to the period of the first harmonic increases as the ratio of the central density to the footpoint density is increased or as the averaged density of the tube is decreased. We find that the damping times due to resonant absorption decrease as the central to footpoint density ratio increases. Contrary to the case of longitudinally homogeneous tubes, we find that the damping time to period ratio also increases as the density ratio is increased or the average density is reduced. We present snapshots and time-distance diagrams of the emission in the H$α$ line. The results presented here have implications for the field of prominence seismology. While the H$α$ emission can be used to detect the fundamental mode, the first harmonic is barely detectable in H$α$. This may explain the lack of detections of the first harmonic. △ Less

Submitted 17 November, 2021; originally announced November 2021.

Comments: 17 pages, 16 figures, 4 movies. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 658, A106 (2022)

One dimensional prominence threads: I. Equilibrium models

Authors: J. Terradas, M. Luna, R. Soler, R. Oliver, M. Carbonell, J. L Ballester

Abstract: Threads are the building blocks of solar prominences and very often show longitudinal oscillatory motions that are strongly attenuated with time. The damping mechanism responsible for the reported oscillations is not fully understood yet. To understand the oscillations and damping of prominence threads it is mandatory to investigate first the nature of the equilibrium solutions that arise under st… ▽ More Threads are the building blocks of solar prominences and very often show longitudinal oscillatory motions that are strongly attenuated with time. The damping mechanism responsible for the reported oscillations is not fully understood yet. To understand the oscillations and damping of prominence threads it is mandatory to investigate first the nature of the equilibrium solutions that arise under static conditions and under the presence of radiative losses, thermal conduction and background heating. This provides the basis to calculate the eigenmodes of the thread models. The nonlinear ordinary differential equations for hydrostatic and thermal equilibrium under the presence of gravity are solved using standard numerical techniques and simple analytical expressions are derived under certain approximations. The solutions to the equations represent a prominence thread, i.e., a dense and cold plasma region of a certain length that connects with the corona through a prominence corona transition region (PCTR). The solutions can also match with a chromospheric-like layer if a spatially dependent heating function localised around the footpoints is considered. We have obtained static solutions representing prominence threads and have investigated in detail the dependence of these solutions on the different parameters of the model. Among other results, we have shown that multiple condensations along a magnetic field line are possible, and that the effect of partial ionisation in the model can significantly modify the thermal balance in the thread and therefore their length. This last parameter is also shown to be comparable to that reported in the observations when the radiative losses are reduced for typical thread temperatures. △ Less

Submitted 11 June, 2021; originally announced June 2021.

Comments: Accepted for publication in A&A

Journal ref: A&A 653, A95 (2021)

Alfven wave heating in partially ionized thin threads of solar prominences

Authors: Llorenc Melis, Roberto Soler, Jose Luis Ballester

Abstract: There is observational evidence of the presence of small-amplitude transverse magnetohydrodynamic (MHD) waves with a wide range of frequencies in the threads of solar prominences. It is believed that the waves are driven at the photosphere and propagate along the magnetic field lines up to prominences suspended in the corona. The dissipation of MHD wave energy in the partially ionized prominence p… ▽ More There is observational evidence of the presence of small-amplitude transverse magnetohydrodynamic (MHD) waves with a wide range of frequencies in the threads of solar prominences. It is believed that the waves are driven at the photosphere and propagate along the magnetic field lines up to prominences suspended in the corona. The dissipation of MHD wave energy in the partially ionized prominence plasma is a heating mechanism whose relevance needs to be explored. Here we consider a simple 1D model for a non-uniform thin thread and investigate the heating associated with dissipation of Alfven waves. The model assumes an ad hoc density profile and a uniform pressure, while the temperature and ionization degree are self-consistently computed considering either LTE or non-LTE approximations for the hydrogen ionization. A broadband driver for Alfven waves is placed at one end of the magnetic field line, representing photospheric excitation. The Alfvenic perturbations along the thread are obtained by solving the linearized MHD equations for a partially ionized plasma in the single-fluid approximation.We find that wave heating in the partially ionized part of the thread is significant enough to compensate for energy losses due to radiative cooling. A greater amount of heating is found in the LTE case because the ionization degree for core prominence temperatures is lower than that in the non-LTE approximation. This results in a greater level of dissipation due to ambipolar diffusion in the LTE case. Conversely, in the hot coronal part of the model, the plasma is fully ionized and wave heating is negligible. The results of this simple model suggest that MHD wave heating can be relevant for the energy balance in prominences. Further studies based on more elaborate models are required. △ Less

Submitted 30 March, 2021; originally announced March 2021.

Journal ref: A&A 650, A45 (2021)

A novel approach to identify resonant MHD wave modes in solar pores and sunspot umbrae: $B-ω$ analysis

Authors: M. Stangalini, D. B. Jess, G. Verth, V. Fedun, B. Fleck, S. Jafarzadeh, P. H. Keys, M. Murabito, D. Calchetti, A. A. Aldhafeeri, F. Berrilli, D. Del Moro, S. M. Jefferies, J. Terradas, R. Soler

Abstract: The umbral regions of sunspots and pores in the solar photosphere are generally dominated by $3$~mHz oscillations, which are due to $p$-modes penetrating the magnetic region. In these locations, wave power is also significantly reduced with respect to the quiet Sun. However, here we study a pore where the power of the oscillations in the umbra is not only comparable, or even larger than that of th… ▽ More The umbral regions of sunspots and pores in the solar photosphere are generally dominated by $3$~mHz oscillations, which are due to $p$-modes penetrating the magnetic region. In these locations, wave power is also significantly reduced with respect to the quiet Sun. However, here we study a pore where the power of the oscillations in the umbra is not only comparable, or even larger than that of the quiet Sun, but the main dominant frequency is not $3$~mHz as expected, but instead $5$~mHz. By combining Doppler velocities and spectropolarimetry and analysing the relationship between magnetic field strength and frequency, the resultant $B-ω$ diagram reveals distinct ridges which are remarkably clear signatures of resonant MHD oscillations confined within the pore umbra. In addition to velocity oscillations, we demonstrate that these modes are also accompanied by magnetic oscillations, as predicted from MHD theory. The novel technique of $B-ω$ analysis, proposed in this Letter, opens an exciting new avenue for identifying MHD wave modes in the umbral regions of both pores and sunspots. △ Less

Submitted 22 March, 2021; originally announced March 2021.

Comments: Accepted in A&A

Journal ref: A&A 649, A169 (2021)

Transition to turbulence in nonuniform coronal loops driven by torsional Alfven waves

Authors: S. Diaz-Suarez, R. Soler

Abstract: Both observations and numerical simulations suggest that Alfvenic waves may carry sufficient energy to sustain the hot temperatures of the solar atmospheric plasma. However, the thermalization of wave energy is inefficient unless very short spatial scales are considered. Phase mixing is a mechanism that can take energy down to dissipation lengths, but it operates over too long a timescale. Here, w… ▽ More Both observations and numerical simulations suggest that Alfvenic waves may carry sufficient energy to sustain the hot temperatures of the solar atmospheric plasma. However, the thermalization of wave energy is inefficient unless very short spatial scales are considered. Phase mixing is a mechanism that can take energy down to dissipation lengths, but it operates over too long a timescale. Here, we study how turbulence, driven by the nonlinear evolution of phase-mixed torsional Alfven waves in coronal loops, is able to take wave energy down to the dissipative scales much faster than the theory of linear phase mixing predicts. We consider a simple model of a transversely nonuniform cylindrical flux tube with a constant axial magnetic field. The flux tube is perturbed by the fundamental mode of standing torsional Alfven waves. We solved the three-dimensional (3D) ideal magnetohydrodynamics equations numerically to study the temporal evolution. Initially, torsional Alfven waves undergo the process of phase mixing because of the transverse variation of density. After only few periods of torsional waves, azimuthal shear flows generated by phase mixing eventually trigger the Kelvin-Helmholtz instability (KHi), and the flux tube is subsequently driven to a turbulent state. Turbulence is very anisotropic and develops transversely only to the background magnetic field. After the onset of turbulence, the effective Reynolds number decreases in the flux tube much faster than in the initial linear stage governed by phase mixing alone. We conclude that the nonlinear evolution of torsional Alfven waves, and the associated KHi, is a viable mechanism for the onset of turbulence in coronal loops. △ Less

Submitted 31 March, 2021; v1 submitted 12 February, 2021; originally announced February 2021.

Comments: Accepted in A&A

Journal ref: A&A 648, A22 (2021)

Resonant absorption: transformation of compressive motions into vortical motions

Authors: M. Goossens, I. Arregui, R. Soler, T. Van Doorsselaere

Abstract: This paper investigates the changes in spatial properties when magnetohydrodynamic (MHD) waves undergo resonant damping in the Alfvén continuum. The analysis is carried out for a 1D cylindrical pressure-less plasma with a straight magnetic field. The effect of the damping on the spatial wave variables is determined by using complex frequencies that arise as a result of the resonant damping. Compre… ▽ More This paper investigates the changes in spatial properties when magnetohydrodynamic (MHD) waves undergo resonant damping in the Alfvén continuum. The analysis is carried out for a 1D cylindrical pressure-less plasma with a straight magnetic field. The effect of the damping on the spatial wave variables is determined by using complex frequencies that arise as a result of the resonant damping. Compression and vorticity are used to characterise the spatial evolution of the MHD wave. The most striking result is the huge spatial variation in the vorticity component parallel to the magnetic field. Parallel vorticity vanishes in the uniform part of the equilibrium. However, when the MHD wave moves into the non-uniform part, parallel vorticity explodes to values that are orders of magnitude higher than those attained by the transverse components in planes normal to the straight magnetic field. In the non-uniform part of the equilibrium plasma, the MHD wave is controlled by parallel vorticity and resembles an Alfvén wave, with the unfamiliar property that it has pressure variations even in the linear regime. △ Less

Submitted 17 September, 2020; originally announced September 2020.

Journal ref: A&A, 641, A106 (2020)

Phase Mixing of Kink MHD Waves in the Solar Corona: Viscous Dissipation and Heating

Authors: Zanyar Ebrahimi, Roberto Soler, Kayoomars Karami

Abstract: 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 wa… ▽ More 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. △ Less

Submitted 9 May, 2020; originally announced May 2020.

Comments: 21 pages, 8 figures

Report number: ApJ 893(2), (2020), 157

Transverse waves in coronal flux tubes with thick boundaries: The effect of longitudinal flows

Authors: Roberto Soler

Abstract: 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én 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 t… ▽ More 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én 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én resonance, while the numerical part of the method consists of solving iteratively the transcendental dispersion relation together with the equation for the Alfvén 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. △ Less

Submitted 30 January, 2019; originally announced January 2019.

Comments: Accepted in A&A

Journal ref: A&A 623, A32 (2019)

Energy transport and heating by torsional Alfvén waves propagating from the photosphere to the corona in the quiet Sun

Authors: Roberto Soler, Jaume Terradas, Ramon Oliver, Jose Luis Ballester

Abstract: In the solar atmosphere, Alfvén 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én waves propagating in magnetic flux tubes that expand from the photosphere to the corona in… ▽ More In the solar atmosphere, Alfvén 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én 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. △ Less

Submitted 4 December, 2018; originally announced December 2018.

Comments: Accepted for publication in ApJ

Multi-fluid approach to high-frequency waves in plasmas. III. Nonlinear regime and plasma heating

Authors: David Martínez-Gómez, Roberto Soler, Jaume Terradas

Abstract: 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 k… ▽ More 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én 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. △ Less

Submitted 22 February, 2018; originally announced February 2018.

Comments: 26 pages, 15 figures, 3 movies. Accepted for publication in The Astrophysical Journal

The effect of a twisted magnetic field on the phase mixing of the kink magnetohydrodynamic waves in coronal loops

Authors: Zanyar Ebrahimi, Kayoomars Karami, Roberto Soler

Abstract: 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.… ▽ More 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é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. △ Less

Submitted 26 August, 2017; originally announced August 2017.

Comments: 25 pages, 22 figures

Journal ref: The Astrophysical Journal, 845:86 (20pp), 2017 August 10

Partially Ionized Plasmas in Astrophysics

Authors: Jose Luis Ballester, Igor Alexeev, Manuel Collados, Turlough Downes, Robert F. Pfaff, Holly Gilbert, Maxim Khodachenko, Elena Khomenko, Ildar F. Shaikhislamov, Roberto Soler, Enrique Vazquez-Semadeni, Teimuraz Zaqarashvili

Abstract: 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… ▽ More 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. △ Less

Submitted 26 July, 2017; v1 submitted 25 July, 2017; originally announced July 2017.

Comments: 163 pages, submitted to Space Science Reviews

Multi-fluid Approach to High-frequency Waves in Plasmas. II. Small-amplitude Regime in Partially Ionized Media

Authors: David Martínez-Gómez, Roberto Soler, Jaume Terradas

Abstract: 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én waves and alters their oscillation frequency and phase speed. When the collision frequencies are larger than the frequency of the waves, single-fluid magnet… ▽ More 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én 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. △ Less

Submitted 15 March, 2017; originally announced March 2017.

Comments: 24 pages, 15 figures. Revised version published in The Astrophysical Journal

Journal ref: ApJ, 837:80, 2017

Multi-fluid approach to high-frequency waves in plasmas: I. Small-amplitude regime in fully ionized medium

Authors: David Martínez-Gómez, Roberto Soler, Jaume Terradas

Abstract: Ideal MHD provides an accurate description of low-frequency Alfvén 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… ▽ More Ideal MHD provides an accurate description of low-frequency Alfvén 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. △ Less

Submitted 20 September, 2016; originally announced September 2016.

Comments: 29 pages, 15 figures. Accepted for publication in The Astrophysical Journal

The role of Alfvén wave heating in solar prominences

Authors: Roberto Soler, Jaume Terradas, Ramon Oliver, Jose Luis Ballester

Abstract: 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 disc… ▽ More 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én 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én 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. △ Less

Submitted 23 May, 2016; originally announced May 2016.

Comments: Accepted for publication in A&A

Journal ref: A&A 592, A28 (2016)

ALMA Observations of the Sun in Cycle 4 and Beyond

Authors: S. Wedemeyer, B. Fleck, M. Battaglia, N. Labrosse, G. Fleishman, H. Hudson, P. Antolin, C. Alissandrakis, T. Ayres, J. Ballester, T. Bastian, J. Black, A. Benz, R. Brajsa, M. Carlsson, J. Costa, B. DePontieu, G. Doyle, G. Gimenez de Castro, S. Gunár, G. Harper, S. Jafarzadeh, M. Loukitcheva, V. Nakariakov, R. Oliver , et al. (12 additional authors not shown)

Abstract: 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 interestin… ▽ More 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. △ Less

Submitted 11 April, 2016; v1 submitted 4 January, 2016; originally announced January 2016.

Comments: SSALMON White Paper with focus on potential solar science with ALMA in Cycle 4; 54 pages. Version 1.2, March 29th, 2016 (updated technical capabilities and observing plans)

Solar prominences embedded in flux ropes: morphological features and dynamics from 3D MHD simulations

Authors: J. Terradas, R. Soler, M. Luna, R. Oliver, J. L. Ballester, A. N. Wright

Abstract: 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, charact… ▽ More 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. △ Less

Submitted 22 December, 2015; originally announced December 2015.

Comments: Submitted to ApJ

Dynamics of coronal rain and descending plasma blobs in solar prominences: II. partially ionized case

Authors: R. Oliver, R. Soler, J. Terradas, T. V. Zaqarashvili

Abstract: 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 materia… ▽ More 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. △ Less

Submitted 28 October, 2015; originally announced October 2015.

Journal ref: 2016, ApJ, 818, 128

Apparent cross-field superslow propagation of magnetohydrodynamic waves in solar plasmas

Authors: T. Kaneko, M. Goossens, R. Soler, J. Terradas, T. Van Doorsselaere, T. Yokoyama, A. N. Wright

Abstract: In this paper we show that the phase mixing of continuum Alfvé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 b… ▽ More In this paper we show that the phase mixing of continuum Alfvé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é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é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én/slow frequency across the magnetic surfaces and is slower than the Alfvé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. △ Less

Submitted 10 September, 2015; originally announced September 2015.

Comments: 41 pages, 7 figures, accepted for publication in ApJ

Model comparison for the density structure across solar coronal waveguides

Authors: I. Arregui, R. Soler, A. Asensio Ramos

Abstract: 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 compar… ▽ More 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én 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. △ Less

Submitted 8 September, 2015; originally announced September 2015.

Comments: 12 pages, 9 figures, 3 tables

Damped transverse oscillations of interacting coronal loops

Authors: Roberto Soler, Manuel Luna

Abstract: 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 observatio… ▽ More 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. △ Less

Submitted 4 September, 2015; originally announced September 2015.

Comments: Accepted for publication in A&A

Journal ref: A&A 582, A120 (2015)

On the spatial scales of wave heating in the solar chromosphere

Authors: Roberto Soler, Marc Carbonell, Jose Luis Ballester

Abstract: 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én waves and slow magnetoacoustic waves. We consider the governing… ▽ More 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én 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én waves depends strongly on the magnetic field strength and ranges from 10~m to 1~km for realistic field strengths. The damping of Alfvén 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. △ Less

Submitted 6 August, 2015; originally announced August 2015.

Comments: Accepted for publication in ApJ

Model comparison for the density structure along solar prominence threads

Authors: I. Arregui, R. Soler

Abstract: 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 pro… ▽ More 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. △ Less

Submitted 13 May, 2015; originally announced May 2015.

Comments: 9 pages, 5 figures, accepted in A&A

Solar science with the Atacama Large Millimeter/submillimeter Array - A new view of our Sun

Authors: S. Wedemeyer, T. Bastian, R. Brajsa, H. Hudson, G. Fleishman, M. Loukitcheva, B. Fleck, E. P. Kontar, B. De Pontieu, P. Yagoubov, S. K. Tiwari, R. Soler, J. H. Black, P. Antolin, E. Scullion, S. Gunar, N. Labrosse, H. -G. Ludwig, A. O. Benz, S. M. White, P. Hauschildt, J. G. Doyle, V. M. Nakariakov, T. Ayres, P. Heinzel , et al. (13 additional authors not shown)

Abstract: 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 play… ▽ More 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. △ Less

Submitted 4 January, 2016; v1 submitted 26 April, 2015; originally announced April 2015.

Comments: 73 pages, 21 figures ; Space Science Reviews (accepted December 10th, 2015); accepted version

Onset of the Kelvin-Helmholtz instability in partially ionized magnetic flux tubes

Authors: David Martínez-Gómez, Roberto Soler, Jaume Terradas

Abstract: 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 surrounde… ▽ More 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énic 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énic flow velocities may be interpreted as consequences of KHI in partially ionized plasmas. △ Less

Submitted 21 April, 2015; originally announced April 2015.

Comments: 8 pages, 4 figures. Accepted for publication in A&A

Journal ref: A&A 578, A104 (2015)

On the support of neutrals against gravity in solar prominences

Authors: J. Terradas, R. Soler, R. Oliver, J. L. Ballester

Abstract: 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… ▽ More 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. △ Less

Submitted 18 March, 2015; originally announced March 2015.

Comments: Accepted for publication in ApJL

SSALMON - The Solar Simulations for the Atacama Large Millimeter Observatory Network

Authors: S. Wedemeyer, T. Bastian, R. Brajsa, M. Barta, H. Hudson, G. Fleishman, M. Loukitcheva, B. Fleck, E. Kontar, B. De Pontieu, S. Tiwari, Y. Kato, R. Soler, P. Yagoubov, J. H. Black, P. Antolin, S. Gunar, N. Labrosse, A. O. Benz, A. Nindos, M. Steffen, E. Scullion, J. G. Doyle, T. Zaqarashvili, A. Hanslmeier , et al. (5 additional authors not shown)

Abstract: 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 s… ▽ More 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. △ Less

Submitted 6 April, 2015; v1 submitted 19 February, 2015; originally announced February 2015.

Comments: 15 pages, 3 figures, to appear in Advances in Solar Physics (special issue of Advances in Space Research, AdSR), proceedings of the 14th European Solar Physics Meeting (ESPM-14, Dublin, Ireland, September 2014); (resubmitted version after minor revision)

Magnetohydrodynamic kink waves in nonuniform solar flux tubes: phase mixing and energy cascade to small scales

Authors: Roberto Soler, Jaume Terradas

Abstract: 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 theor… ▽ More 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én 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. △ Less

Submitted 13 February, 2015; originally announced February 2015.

Comments: Accepted for publication in ApJ

Prominence seismology using the period ratio of transverse thread oscillations

Authors: R. Soler, M. Goossens, J. L. Ballester

Abstract: 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 w… ▽ More 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. △ Less

Submitted 21 January, 2015; originally announced January 2015.

Comments: Accepted for publication in A&A

Journal ref: A&A 575, A123 (2015)

Morphology and dynamics of solar prominences from 3D MHD simulations

Authors: J. Terradas, R. Soler, M. Luna, R. Oliver, J. L. Ballester

Abstract: 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 param… ▽ More 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$-β$ 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. △ Less

Submitted 23 December, 2014; originally announced December 2014.

Comments: Accepted for publication in ApJ

Overdamped Alfven waves due to ion-neutral collisions in the solar chromosphere

Authors: R. Soler, J. L. Ballester, T. V. Zaqarashvili

Abstract: 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 s… ▽ More 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. △ Less

Submitted 21 November, 2014; originally announced November 2014.

Comments: Accepted for publication in A&A

Journal ref: A&A 573, A79 (2015)

The behavior of transverse waves in nonuniform solar flux tubes. II. Implications for coronal loop seismology

Authors: Roberto Soler, Marcel Goossens, Jaume Terradas, Ramon Oliver

Abstract: 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 approximat… ▽ More 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. △ Less

Submitted 18 December, 2013; originally announced December 2013.

Comments: Accepted for publication in ApJ

Dynamics of coronal rain and descending plasma blobs in solar prominences: I. Fully ionised case

Authors: R. Oliver, R. Soler, J. Terradas, T. V. Zaqarashvili, M. L. Khodachenko

Abstract: 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… ▽ More 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. △ Less

Submitted 4 November, 2013; originally announced November 2013.

Journal ref: 2014, ApJ, 784, 21

Magnetoacoustic waves in a partially ionized two-fluid plasma

Authors: Roberto Soler, Marc Carbonell, Jose Luis Ballester

Abstract: 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 t… ▽ More 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 $β$, 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. △ Less

Submitted 27 September, 2013; originally announced September 2013.

Comments: Accepted for publication in ApJ Supplement Series

MHD waves in two-dimensional prominences embedded in coronal arcades

Authors: J. Terradas, R. Soler, A. J. Díaz, R. Oliver, J. L. Ballester

Abstract: 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… ▽ More 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. △ Less

Submitted 19 September, 2013; originally announced September 2013.

The behavior of transverse waves in nonuniform solar flux tubes. I. Comparison of ideal and resistive results

Authors: Roberto Soler, Marcel Goossens, Jaume Terradas, Ramon Oliver

Abstract: 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 o… ▽ More 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én 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. △ Less

Submitted 13 September, 2013; originally announced September 2013.

Comments: Accepted for publication in ApJ

Alfven waves in a partially ionized two-fluid plasma

Authors: R. Soler, M. Carbonell, J. L. Ballester, J. Terradas

Abstract: Alfvén waves are a particular class of magnetohydrodynamic waves relevant in many astrophysical and laboratory plasmas. In partially ionized plasmas the dynamics of Alfvén waves is affected by the interaction between ionized and neutral species. Here we study Alfvén waves in a partially ionized plasma from the theoretical point of view using the two-fluid description. We consider that the plasma i… ▽ More Alfvén waves are a particular class of magnetohydrodynamic waves relevant in many astrophysical and laboratory plasmas. In partially ionized plasmas the dynamics of Alfvén waves is affected by the interaction between ionized and neutral species. Here we study Alfvén 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én waves in the low solar atmospheric plasma is performed and the implication of partial ionization for the energy flux is discussed. △ Less

Submitted 18 March, 2013; originally announced March 2013.

Comments: Accepted for publication in ApJ

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

Authors: R. Soler, A. J. Diaz, J. L. Ballester, M. Goossens

Abstract: 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… ▽ More 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. △ Less

Submitted 22 January, 2013; originally announced January 2013.

Comments: Accepted for publication in A&A

Torsional Alfvén waves in solar partially ionized plasma: effects of neutral helium and stratification

Authors: T. V. Zaqarashvili, M. L. Khodachenko, R. Soler

Abstract: 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 e… ▽ More 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. △ Less

Submitted 10 May, 2013; v1 submitted 6 November, 2012; originally announced November 2012.

Comments: 9 pages, 7 figures (accepted in A&A)

Analytic approximate seismology of propagating MHD waves in the solar corona

Authors: M. Goossens, R. Soler, I. Arregui, J. Terradas

Abstract: 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 co… ▽ More 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. △ Less

Submitted 9 October, 2012; originally announced October 2012.

Comments: Accepted for publication in ApJ. arXiv admin note: text overlap with arXiv:0804.3877

Damped kink oscillations of flowing prominence threads

Authors: R. Soler, M. S. Ruderman, M. Goossens

Abstract: 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 h… ▽ More 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. △ Less

Submitted 15 September, 2012; originally announced September 2012.

Comments: Accepted for publication in A&A

Surface Alfven waves in solar flux tubes

Authors: M. Goossens, J. Andries, R. Soler, T. Van Doorsselaere, I. Arregui, J. Terradas

Abstract: 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… ▽ More 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. △ Less

Submitted 4 May, 2012; originally announced May 2012.

Comments: Accepted for publication in ApJ

Kelvin-Helmholtz instability in partially ionized compressible plasmas

Authors: R. Soler, A. J. Diaz, J. L. Ballester, M. Goossens

Abstract: 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… ▽ More 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. △ Less

Submitted 20 February, 2012; originally announced February 2012.

Comments: Accepted for publication in ApJ