Showing 1–50 of 67 results for author: Terradas, J

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

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)

Three-dimensional solar active region magnetohydrostatic models and their stability using Euler potentials

Authors: Jaume Terradas, Thomas Neukirch

Abstract: Active regions (ARs) are typical magnetic structures found in the solar atmosphere. We calculate several magnetohydrostatic (MHS) equilibrium models that include the effect of a finite plasma-$β$ and gravity and that are representative of these structures in three dimensions. The construction of the models is based on the use of two Euler potentials, $α$ and $β$, that represent the magnetic field… ▽ More Active regions (ARs) are typical magnetic structures found in the solar atmosphere. We calculate several magnetohydrostatic (MHS) equilibrium models that include the effect of a finite plasma-$β$ and gravity and that are representative of these structures in three dimensions. The construction of the models is based on the use of two Euler potentials, $α$ and $β$, that represent the magnetic field as ${\bf B}=\nabla α\times \nabla β$. The ideal MHS nonlinear partial differential equations are solved numerically using finite elements in a fixed 3D rectangular domain. The boundary conditions are initially chosen to correspond to a potential magnetic field (current-free) with known analytical expressions for the corresponding Euler potentials. The distinctive feature is that we incorporate the effect of shear by progressively deforming the initial potential magnetic field. This procedure is quite generic and allows us to generate a vast variety of MHS models. The thermal structure of the ARs is incorporated through the dependence of gas pressure and temperature on the Euler potentials. Using this method we achieve the characteristic hot and over-dense plasma found in ARs, but we demonstrate that the method can also be applied to study configurations with open magnetic field lines. Furthermore, we investigate basic topologies that include neutral lines. Our focus is on the force balance of the structures and we do not consider the energy balance in the constructed models. In addition, we address the difficult question of the stability of the calculated 3D models. We find that if the plasma is convectively stable, then the system is not prone in general to develop magnetic Rayleigh-Taylor instabilities. △ Less

Submitted 9 December, 2022; originally announced December 2022.

Comments: Accepted for publication in A&A

Journal ref: A&A 671, A31 (2023)

Extension and validation of the pendulum model for longitudinal solar prominence oscillations

Authors: M. Luna, J. Terradas, J. Karpen, J. L. Ballester

Abstract: Longitudinal oscillations in prominences are common phenomena on the Sun. These oscillations can be used to infer the geometry and intensity of the filament magnetic field. Previous theoretical studies of longitudinal oscillations made two simplifying assumptions: uniform gravity and semi-circular dips on the supporting flux tubes. However, the gravity is not uniform and realistic dips are not sem… ▽ More Longitudinal oscillations in prominences are common phenomena on the Sun. These oscillations can be used to infer the geometry and intensity of the filament magnetic field. Previous theoretical studies of longitudinal oscillations made two simplifying assumptions: uniform gravity and semi-circular dips on the supporting flux tubes. However, the gravity is not uniform and realistic dips are not semi-circular. To understand the effects of including the nonuniform solar gravity on longitudinal oscillations, and explore the validity of the pendulum model with different flux-tube geometries. We first derive the equation describing the motion of the plasma along the flux tube including the effects of nonuniform gravity, yielding corrections to the original pendulum model. We also compute the full numerical solutions for the normal modes, and compare them with the new pendulum approximation. We have found that the nonuniform gravity introduces a significant modification in the pendulum model. We have also found a cut-off period, i.e. the longitudinal oscillations cannot have a period longer than 167 minutes. In addition, considering different tube geometries, the period depends almost exclusively on the radius of curvature at the bottom of the dip. We conclude that nonuniform gravity significantly modifies the pendulum model. These corrections are important for prominence seismology, because the inferred values of the radius of curvature and minimum magnetic-field strength differ substantially from those of the old model. However, we find that the corrected pendulum model is quite robust and is still valid for non-circular dips. △ Less

Submitted 16 February, 2022; originally announced February 2022.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 660, A54 (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

Excitation of ion-acoustic waves by non-linear finite-amplitude standing Alfvén waves

Authors: Jaume Terradas, Adolfo F. - Viñas, Jaime A. Araneda

Abstract: We investigate, using a multi-fluid approach, the main properties of standing ion-acoustic modes driven by nonlinear standing Alfvén waves. The standing character of the Alfvénic pump is because we study the superposition of two identical circularly polarised counter-propagating waves. We consider parallel propagation along the constant magnetic field and we find that left and right-handed modes g… ▽ More We investigate, using a multi-fluid approach, the main properties of standing ion-acoustic modes driven by nonlinear standing Alfvén waves. The standing character of the Alfvénic pump is because we study the superposition of two identical circularly polarised counter-propagating waves. We consider parallel propagation along the constant magnetic field and we find that left and right-handed modes generate via ponderomotive forces the second harmonic of standing ion-acoustic waves. We demonstrate that parametric instabilities are not relevant in the present problem and the secondary ion-acoustic waves attenuate by Landau damping in the absence of any other dissipative process. Kinetic effects are included in our model where ions are considered as particles and electrons as a massless fluid, and hybrid simulations are used to complement the theoretical results. Analytical expressions are obtained for the time evolution of the different physical variables in the absence of Landau damping. From the hybrid simulations we find that the attenuation of the generated ion-acoustic waves follows the theoretical predictions even under the presence of a driver Alfvénic pump. Due to the nonlinear induced ion-acoustic waves the system develops density cavities and an electric field parallel to the magnetic field. Theoretical expressions for this density and electric field fluctuations are derived. The implications of these results in the context of standing slow mode oscillations in coronal loops is discussed. △ Less

Submitted 24 December, 2021; originally announced December 2021.

Journal ref: A&A 660, A24 (2022)

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)

Kink oscillations of coronal loops

Authors: V. M. Nakariakov, S. A. Anfinogentov, P. Antolin, R. Jain, D. Y. Kolotkov, E. G. Kupriyanova, D. Li, N. Magyar, G. Nistico, D. J. Pascoe, A. K. Srivastava, J. Terradas, S. Vashegani Farahani, G. Verth, D. Yuan, I. V. Zimovets

Abstract: Kink oscillations of coronal loops, i.e., standing kink waves, is one of the most studied dynamic phenomena in the solar corona. The oscillations are excited by impulsive energy releases, such as low coronal eruptions. Typical periods of the oscillations are from a few to several minutes, and are found to increase linearly with the increase in the major radius of the oscillating loops. It clearly… ▽ More Kink oscillations of coronal loops, i.e., standing kink waves, is one of the most studied dynamic phenomena in the solar corona. The oscillations are excited by impulsive energy releases, such as low coronal eruptions. Typical periods of the oscillations are from a few to several minutes, and are found to increase linearly with the increase in the major radius of the oscillating loops. It clearly demonstrates that kink oscillations are natural modes of the loops, and can be described as standing fast magnetoacoustic waves with the wavelength determined by the length of the loop. Kink oscillations are observed in two different regimes. In the rapidly decaying regime, the apparent displacement amplitude reaches several minor radii of the loop. The damping time which is about several oscillation periods decreases with the increase in the oscillation amplitude, suggesting a nonlinear nature of the damping. In the decayless regime, the amplitudes are smaller than a minor radius, and the driver is still debated. The review summarises major findings obtained during the last decade, and covers both observational and theoretical results. Observational results include creation and analysis of comprehensive catalogues of the oscillation events, and detection of kink oscillations with imaging and spectral instruments in the EUV and microwave bands. Theoretical results include various approaches to modelling in terms of the magnetohydrodynamic wave theory. Properties of kink oscillations are found to depend on parameters of the oscillating loop, such as the magnetic twist, stratification, steady flows, temperature variations and so on, which make kink oscillations a natural probe of these parameters by the method of magnetohydrodynamic seismology. △ Less

Submitted 23 September, 2021; originally announced September 2021.

Journal ref: Space Science Reviews, 2021, Volume 217, Issue 6, article id.73

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)

Gravitational instability of solar prominence threads I. Curved magnetic fields without dips

Authors: A. Adrover-González, J. Terradas, R. Oliver, M. Carbonell

Abstract: Prominence threads are dense and cold structures lying on curved magnetic fields that can be suspended in the solar atmosphere against gravity. The gravitational stability of threads, in the absence of non-ideal effects, is comprehensively investigated in the present work by means of an elementary but effective model. Based on purely hydrodynamic equations in one spatial dimension and applying lin… ▽ More Prominence threads are dense and cold structures lying on curved magnetic fields that can be suspended in the solar atmosphere against gravity. The gravitational stability of threads, in the absence of non-ideal effects, is comprehensively investigated in the present work by means of an elementary but effective model. Based on purely hydrodynamic equations in one spatial dimension and applying line-tying conditions at the footpoints of the magnetic field lines, we derive analytical expressions for the different feasible equilibria and the corresponding frequencies of oscillation. We find that the system allows for stable and unstable equilibrium solutions subject to the initial position of the thread, its density contrast and length, and the total length of the magnetic field lines. The transition between the two types of solutions is produced at specific bifurcation points that have been determined analytically in some particular cases. When the thread is initially at the top of the concave magnetic field, that is at the apex, we find a supercritical pitchfork bifurcation, while for a shifted initial thread position with respect to this point the symmetry is broken and the system is characterised by an S-shaped bifurcation. The plain results presented in this paper shed new light on the behaviour of threads in curved magnetic fields under the presence of gravity and help to interpret more complex numerical magnetohydrodynamics (MHD) simulations about similar structures. △ Less

Submitted 22 March, 2021; originally announced March 2021.

Comments: 13 pages, 9 figures

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

Geometrical properties of the interaction between oblique incoming coronal waves and coronal holes

Authors: I. Piantschitsch, J. Terradas

Abstract: Observations of coronal waves (CWs) interacting with coronal holes (CHs) show the formation of typical wave-like features, such as reflected, refracted and transmitted waves (collectively, secondary waves). In accordance with these observations, numerical evidence for the wave characteristics of CWs is given by simulations which demonstrate effects of deflection and reflection when a CW interacts… ▽ More Observations of coronal waves (CWs) interacting with coronal holes (CHs) show the formation of typical wave-like features, such as reflected, refracted and transmitted waves (collectively, secondary waves). In accordance with these observations, numerical evidence for the wave characteristics of CWs is given by simulations which demonstrate effects of deflection and reflection when a CW interacts with regions exhibiting a sudden density drop, such as CHs. However, secondary waves are usually weak in their signal and simulations are limited in the way the according idealisations have to be chosen. Hence, several properties of the secondary waves during a CW-CH interaction are unclear or ambiguous and might lead to misinterpretations. In this study we follow a theoretical approach and focus in particular on the geometrical properties of secondary waves caused by the interaction between oblique incoming CWs and CHs. Based on a linear theory, we derive analytical expressions for reflection and transmission coefficients, which tell us how strongly the amplitudes of the secondary waves increase and decrease with respect to the incoming wave, respectively. Additionally, we provide analytical terms for crucial incidence angles that are capable of giving information about the energy flux, the phase and the reflection properties of the secondary waves. These novel expressions provide a supplementary tool for estimating CW properties in a fast and straightforward way and, therefore, might have relevant consequences for a possible new interpretation of already studied CW-CH interaction events and the clarification of ambiguous observational data. △ Less

Submitted 22 March, 2021; originally announced March 2021.

Comments: 15 pages, 10 figures

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)

A new method for estimating global coronal wave properties from their interaction with solar coronal holes

Authors: I. Piantschitsch, J. Terradas, M. Temmer

Abstract: Global coronal waves (CWs) and their interaction with coronal holes (CHs) result, among other effects, in the formation of reflected and transmitted waves. Observations of such events provide us with measurements of different CW parameters, such as phase speed and intensity amplitudes. However, several of these parameters are provided with only intermediate observational quality, other parameters,… ▽ More Global coronal waves (CWs) and their interaction with coronal holes (CHs) result, among other effects, in the formation of reflected and transmitted waves. Observations of such events provide us with measurements of different CW parameters, such as phase speed and intensity amplitudes. However, several of these parameters are provided with only intermediate observational quality, other parameters, such as the phase speed of transmitted waves, can hardly be observed in general. We present a new method to estimate crucial CW parameters, such as density and phase speed of reflected as well as transmitted waves, Mach numbers and density values of the CH's interior, by using analytical expressions in combination with basic and most accessible observational measurements. The transmission and reflection coefficients are derived from linear theory and subsequently used to calculate estimations for phase speeds of incoming, reflected and transmitted waves. The obtained analytical expressions are validated by performing numerical simulations of CWs interacting with CHs. This new method enables to determine in a fast and straightforward way reliable CW and CH parameters from basic observational measurements which provides a powerful tool to better understand the observed interaction effects between CWs and CHs. △ Less

Submitted 12 June, 2020; originally announced June 2020.

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

3D numerical simulations of oscillations in solar prominences

Authors: A. Adrover-González, J. Terradas

Abstract: Oscillations in solar prominences are a frequent phenomenon, and they have been the subject of many studies. A full understanding of the mechanisms that drive them and their attenuation has not been reached yet. We numerically investigate the periodicity and damping of transverse and longitudinal oscillations in a 3D model of a curtain-shaped prominence. We carried out a set of numerical simulatio… ▽ More Oscillations in solar prominences are a frequent phenomenon, and they have been the subject of many studies. A full understanding of the mechanisms that drive them and their attenuation has not been reached yet. We numerically investigate the periodicity and damping of transverse and longitudinal oscillations in a 3D model of a curtain-shaped prominence. We carried out a set of numerical simulations of vertical, transverse and longitudinal oscillations with the high-order finite-difference Pencil Code. We solved the ideal magnetohydrodynamic (MHD) equations for a wide range of parameters, including the width and density of the prominence, and the magnetic field strength (B) of the solar corona. We studied the periodicity and attenuation of the induced oscillations. We found that longitudinal oscillations can be fit with the pendulum model, whose restoring force is the field aligned component of gravity, but other mechanisms such as pressure gradients may contribute to the movement. On the other hand, transverse oscillations are subject to magnetic forces. The analysis of the parametric survey shows, in agreement with observational studies, that the oscillation period (P) increases with the prominence width. For transverse oscillations we obtained that P increases with density and decreases with B. For longitudinal oscillations we also found that P increases with density, but there are no variations with B. The attenuation of transverse oscillations was investigated by analysing the velocity distribution and computing the Alfvén continuum modes. We conclude that resonant absorption is the mean cause. Damping of longitudinal oscillations is due to some kind of shear numerical viscosity. △ Less

Submitted 9 December, 2019; originally announced December 2019.

Journal ref: A&A 633, A113 (2020)

Normal modes of transverse coronal loop oscillations from numerical simulations: I. Method and test case

Authors: S. Rial, I. Arregui, R. Oliver, J. Terradas

Abstract: The purpose of this work is to develop a procedure to obtain the normal modes of a coronal loop from time-dependent numerical simulations with the aim of better understanding observed transverse loop oscillations. To achieve this goal, in this paper we present a new method and test its performance with a problem for which the normal modes can be computed analytically. In a follow-up paper, the app… ▽ More The purpose of this work is to develop a procedure to obtain the normal modes of a coronal loop from time-dependent numerical simulations with the aim of better understanding observed transverse loop oscillations. To achieve this goal, in this paper we present a new method and test its performance with a problem for which the normal modes can be computed analytically. In a follow-up paper, the application to the simulations of \citet{rial2013} is tackled. The method proceeds iteratively and at each step consists of (i) a time-dependent numerical simulation followed by (ii) the Complex Empirical Orthogonal Function (CEOF) analysis of the simulation results. The CEOF analysis provides an approximation to the normal mode eigenfunctions that can be used to set up the initial conditions for the numerical simulation of the following iteration, in which an improved normal mode approximation is obtained. The iterative process is stopped once the global difference between successive approximate eigenfunctions is below a prescribed threshold. The equilibrium used in this paper contains material discontinuities that result in one eigenfunction with a jump across these discontinuities and two eigenfunctions whose normal derivatives are discontinuous there. After 6 iterations, the approximation to the frequency and eigenfunctions are accurate to $\lesssim 0.7$\% except for the eigenfunction with discontinuities, which displays a much larger error at these positions. △ Less

Submitted 13 June, 2019; v1 submitted 1 February, 2019; originally announced February 2019.

Journal ref: Rial et al. (2019) ApJ, 876, 86

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

GONG Catalog of Solar Filament Oscillations Near Solar Maximum

Authors: Manuel Luna, Judith Karpen, José Luís Ballester, Karin Muglach, Jaume Terradas, Therese Kucera, Holly Gilbert

Abstract: We have catalogued 196 filament oscillations from the GONG $Hα$ network data during several months near the maximum of solar cycle 24 (January - June 2014). Selected examples from the catalog are described in detail, along with our statistical analyses of all events. Oscillations were classified according to their velocity amplitude: 106 small-amplitude oscillations (SAOs), with velocities… ▽ More We have catalogued 196 filament oscillations from the GONG $Hα$ network data during several months near the maximum of solar cycle 24 (January - June 2014). Selected examples from the catalog are described in detail, along with our statistical analyses of all events. Oscillations were classified according to their velocity amplitude: 106 small-amplitude oscillations (SAOs), with velocities $<10\mathrm{\, km \; s^{-1}}$, and 90 large-amplitude oscillations (LAOs), with velocities $>10\mathrm{\, km \; s^{-1}}$. Both SAOs and LAOs are common, with one event of each class every two days on the visible side of the Sun. For nearly half of the events we identified their apparent trigger. The period distribution has a mean value of 58$\pm$15 min for both types of oscillations. The distribution of the damping time per period peaks at $τ/P=1.75$ and $1.25$ for SAOs and LAOs respectively. We confirmed that LAO damping rates depend nonlinearly on the oscillation velocity. The angle between the direction of motion and the filament spine has a distribution centered at $27^\circ$ for all filament types. This angle agrees with the observed direction of filament-channel magnetic fields, indicating that most of the catalogued events are longitudinal (i.e., undergo field-aligned motions). We applied seismology to determine the average radius of curvature in the magnetic dips, $R\approx89$ Mm, and the average minimum magnetic-field strength, $B\approx16$ G. The catalog is available to the community online, and is intended to be expanded to cover at least 1 solar cycle. △ Less

Submitted 10 April, 2018; originally announced April 2018.

Comments: Accepted for publication in The Astrophysical Journal Supplement Series

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

Effect of magnetic twist on nonlinear transverse kink oscillations

Authors: J. Terradas, N. Magyar, T. Van Doorsselaere

Abstract: Magnetic twist is thought to play an important role in many structures of the solar atmosphere. One of the effects of twist is to modify the properties of the eigenmodes of magnetic tubes. In the linear regime standing kink solutions are characterized by a change in polarization of the transverse displacement along the twisted tube. In the nonlinear regime, magnetic twist affects the development o… ▽ More Magnetic twist is thought to play an important role in many structures of the solar atmosphere. One of the effects of twist is to modify the properties of the eigenmodes of magnetic tubes. In the linear regime standing kink solutions are characterized by a change in polarization of the transverse displacement along the twisted tube. In the nonlinear regime, magnetic twist affects the development of shear instabilities that appear at the tube boundary when it is oscillating laterally. These Kelvin-Helmholtz instabilities (KHI) are produced either by the jump in the azimuthal component of the velocity at the edge of the sharp boundary between the internal and external part of the tube, or either by the continuous small length-scales produced by phase-mixing when there is a smooth inhomogeneous layer. In this work the effect of twist is consistently investigated by solving the time-dependent problem including the process of energy transfer to the inhomogeneous layer. It is found that twist always delays the appearance of the shear instability but for tubes with thin inhomogeneous layers the effect is relatively small for moderate values of twist. On the contrary, for tubes with thick layers, the effect of twist is much stronger. This can have some important implications regarding observations of transverse kink modes and the KHI itself. △ Less

Submitted 19 December, 2017; originally announced December 2017.

Comments: Accepted for publication in ApJ (11-22-2017)

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 effects of magnetic-field geometry on longitudinal oscillations of solar prominences: Cross-sectional area variation for thin tubes

Authors: M. Luna, A. J. Diaz, R. Oliver, J. Terradas, J. Karpen

Abstract: Solar prominences are subject to both field-aligned (longitudinal) and transverse oscillatory motions, as evidenced by an increasing number of observations. Large-amplitude longitudinal motions provide valuable information on the geometry of the filament-channel magnetic structure that supports the cool prominence plasma against gravity. Our pendulum model, in which the restoring force is the grav… ▽ More Solar prominences are subject to both field-aligned (longitudinal) and transverse oscillatory motions, as evidenced by an increasing number of observations. Large-amplitude longitudinal motions provide valuable information on the geometry of the filament-channel magnetic structure that supports the cool prominence plasma against gravity. Our pendulum model, in which the restoring force is the gravity projected along the dipped field lines of the magnetic structure, best explains these oscillations. However, several factors can influence the longitudinal oscillations, potentially invalidating the pendulum model. The aim of this work is to study the influence of large-scale variations in the magnetic field strength along the field lines, i.e., variations of the cross-sectional area along the flux tubes supporting prominence threads. We studied the normal modes of several flux tube configurations, using linear perturbation analysis, to assess the influence of different geometrical parameters on the oscillation properties. We found that the influence of the symmetric and asymmetric expansion factors on longitudinal oscillations is small.}{We conclude that the longitudinal oscillations are not significantly influenced by variations of the cross-section of the flux tubes, validating the pendulum model in this context. △ Less

Submitted 11 July, 2016; originally announced July 2016.

Comments: Accepted for publication in Astronomy & Astrophysics

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)

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

On the robustness of the pendulum model for large-amplitude longitudinal oscillations in prominences

Authors: M. Luna, J. Terradas, E. Khomenko, M. Collados, A. de Vicente

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

Submitted 16 December, 2015; originally announced December 2015.

Comments: Accepted for publication in The Astrophysical Journal

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

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

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

Propagation and dispersion of sausage wave trains in magnetic flux tubes

Authors: R. Oliver, M. S. Ruderman, J. Terradas

Abstract: A localized perturbation of a magnetic flux tube produces a pair of wave trains that propagate in opposite directions along the tube. These wave packets disperse as they propagate, where the extent of dispersion depends on the physical properties of the magnetic structure, on the length of the initial excitation, and on its nature (e.g., transverse or axisymmetric). In Oliver et al. (2014) we cons… ▽ More A localized perturbation of a magnetic flux tube produces a pair of wave trains that propagate in opposite directions along the tube. These wave packets disperse as they propagate, where the extent of dispersion depends on the physical properties of the magnetic structure, on the length of the initial excitation, and on its nature (e.g., transverse or axisymmetric). In Oliver et al. (2014) we considered a transverse initial perturbation, whereas the temporal evolution of an axisymmetric one is examined here. In both papers we use a method based on Fourier integrals to solve the initial value problem. Previous studies on wave propagation in magnetic wave guides have emphasized that the wave train dispersion is influenced by the particular dependence of the group velocity on the longitudinal wavenumber. Here we also find that long initial perturbations result in low amplitude wave packets and that large values of the magnetic tube to environment density ratio yield longer wave trains. To test the detectability of propagating transverse or axisymmetric wave packets in magnetic tubes of the solar atmosphere (e.g., coronal loops, spicules, or prominence threads) a forward modelling of the perturbations must be carried out. This is left for a future work. △ Less

Submitted 4 February, 2015; originally announced February 2015.

Journal ref: 2015, ApJ, 806, 56

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

Rayleigh-Taylor instabilities with sheared magnetic fields

Authors: M. S. Ruderman, J. Terradas, J. L. Ballester

Abstract: Magnetic Rayleigh-Taylor (MRT) instabilities may play a relevant role in many astrophysical problems. In this work the effect of magnetic shear on the growth rate of the MRT instability is investigated. The eigenmodes of an interface and a slab model under the presence of gravity are analytically calculated assuming that the orientation of the magnetic field changes in the equilibrium, i.e., there… ▽ More Magnetic Rayleigh-Taylor (MRT) instabilities may play a relevant role in many astrophysical problems. In this work the effect of magnetic shear on the growth rate of the MRT instability is investigated. The eigenmodes of an interface and a slab model under the presence of gravity are analytically calculated assuming that the orientation of the magnetic field changes in the equilibrium, i.e., there is magnetic shear. We solve the linearised magnetohydrodynamic (MHD) equations in the incompressible regime.We find that the growth rate is bounded under the presence of magnetic shear. We have derived simple analytical expressions for the maximum growth rate, corresponding to the most unstable mode of the system. These expressions provide the explicit dependence of the growth rate on the various equilibrium parameters. For small angles the growth time is linearly proportional to the shear angle, and in this regime the single interface problem and the slab problem tend to the same result. On the contrary, in the limit of large angles and for the interface problem the growth time is essentially independent of the shear angle. In this regime we have also been able to calculate an approximate expression for the growth time for the slab configuration. Magnetic shear can have a strong effect on the growth rates of the instability. As an application of the results found in this paper we have indirectly determined the shear angle in solar prominence threads using their lifetimes and the estimation of the Alfvén speed of the structure. △ Less

Submitted 26 February, 2014; originally announced February 2014.

Propagation and dispersion of transverse wave trains in magnetic flux tubes

Authors: R. Oliver, M. S. Ruderman, J. Terradas

Abstract: The dispersion of small amplitude, impulsively excited wave trains propagating along a magnetic flux tube is investigated. The initial disturbance is a localized transverse displacement of the tube that excites a fast kink wave packet. The spatial and temporal evolution of the perturbed variables (density, plasma displacement, velocity, ...) is given by an analytical expression containing an integ… ▽ More The dispersion of small amplitude, impulsively excited wave trains propagating along a magnetic flux tube is investigated. The initial disturbance is a localized transverse displacement of the tube that excites a fast kink wave packet. The spatial and temporal evolution of the perturbed variables (density, plasma displacement, velocity, ...) is given by an analytical expression containing an integral that is computed numerically. We find that the dispersion of fast kink wave trains is more important for shorter initial disturbances (i.e. more concentrated in the longitudinal direction) and for larger density ratios (i.e. for larger contrasts of the tube density with respect to the environment density). This type of excitation generates a wave train whose signature at a fixed position along a coronal loop is a short event (duration ~ 20 s) in which the velocity and density oscillate very rapidly with typical periods of the order of a few seconds. The oscillatory period is not constant but gradually declines during the course of this event. Peak values of the velocity are of the order of 10 km/s and are accompanied by maximum density variations of the order of 10-15% the unperturbed loop density. △ Less

Submitted 18 February, 2014; originally announced February 2014.

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

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

Transverse kink oscillations in the presence of twist

Authors: J. Terradas, M. Goossens

Abstract: Magnetic twist is thought to play an important role in coronal loops. The effects of magnetic twist on stable magnetohydrodynamic (MHD) waves is poorly understood because they are seldom studied for relevant cases. The goal of this work is to study the fingerprints of magnetic twist on stable transverse kink oscillations. We numerically calculated the eigenmodes of propagating and standing MHD wav… ▽ More Magnetic twist is thought to play an important role in coronal loops. The effects of magnetic twist on stable magnetohydrodynamic (MHD) waves is poorly understood because they are seldom studied for relevant cases. The goal of this work is to study the fingerprints of magnetic twist on stable transverse kink oscillations. We numerically calculated the eigenmodes of propagating and standing MHD waves for a model of a loop with magnetic twist. The azimuthal component of the magnetic field was assumed to be small in comparison to the longitudinal component. We did not consider resonantly damped modes or kink instabilities in our analysis. For a nonconstant twist the frequencies of the MHD wave modes are split, which has important consequences for standing waves. This is different from the degenerated situation for equilibrium models with constant twist, which are characterised by an azimuthal component of the magnetic field that linearly increases with the radial coordinate. In the presence of twist standing kink solutions are characterised by a change in polarisation of the transverse displacement along the tube. For weak twist, and in the thin tube approximation, the frequency of standing modes is unaltered and the tube oscillates at the kink speed of the corresponding straight tube. The change in polarisation is linearly proportional to the degree of twist. This has implications with regard to observations of kink modes, since the detection of this variation in polarisation can be used as an indirect method to estimate the twist in oscillating loops. △ Less

Submitted 30 October, 2012; originally announced October 2012.

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

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

The role of Rayleigh-Taylor instabilities in filament threads

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

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

Submitted 20 March, 2012; originally announced March 2012.

Prominence seismology

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

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

Submitted 22 January, 2012; originally announced January 2012.

Comments: 8 pages, 3 figures

Journal ref: 4th Hinode Science Meeting: Unsolved Problems and New Insights ASP Conference Series, Vol. 455, pp. 211-218 (2012)

Wave Leakage and Resonant Absorption in a Loop Embedded in a Coronal Arcade

Authors: S. Rial, I. Arregui, J. Terradas, R. Oliver, J. L. Ballester

Abstract: We investigate the temporal evolution of impulsively generated perturbations in a potential coronal arcade with an embedded loop. As the initial configuration we consider a coronal loop, represented by a density enhancement, which is unbounded in the ignorable direction of the arcade. The linearized time-dependent magnetohydrodynamic equations have been numerically solved in field-aligned coordina… ▽ More We investigate the temporal evolution of impulsively generated perturbations in a potential coronal arcade with an embedded loop. As the initial configuration we consider a coronal loop, represented by a density enhancement, which is unbounded in the ignorable direction of the arcade. The linearized time-dependent magnetohydrodynamic equations have been numerically solved in field-aligned coordinates and the time evolution of the initial perturbations has been studied in the zero-beta approximation. For propagation constrained to the plane of the arcade, the considered initial perturbations do not excite trapped modes of the system. This weakness of the model is overcome by the inclusion of wave propagation in the ignorable direction. The inclusion of perpendicular propagation produces two main results. First, damping by wave leakage is less efficient because the loop is able to act as a wave trap of vertical oscillations. Second, the consideration of an inhomogeneous corona enables the resonant damping of vertical oscillations and the energy transfer from the interior of the loop to the external coronal medium. △ Less

Submitted 19 January, 2012; originally announced January 2012.

Resonantly Damped Propagating Kink Waves in Longitudinally Stratified Solar Waveguides

Authors: Roberto Soler, Jaume Terradas, Gary Verth, Marcel Goossens

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

Submitted 30 April, 2011; originally announced May 2011.

Comments: Accepted for publication in ApJ

Spatial Damping of Propagating Kink Waves Due to Resonant Absorption: Effect of Background Flow

Authors: Roberto Soler, Jaume Terradas, Marcel Goossens

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

Submitted 10 April, 2011; originally announced April 2011.

Comments: Accepted for publication in ApJ

Seismology of transversely oscillating coronal loops with siphon flows

Authors: J. Terradas, I. Arregui, G. Verth, M. Goossens

Abstract: There are ubiquitous flows observed in the solar atmosphere of sub-Alfvénic speeds, however after flaring and coronal mass ejection events flows can become Alfvénic. In this Letter, we derive an expression for the standing kink mode frequency due to siphon flow in coronal loops, valid for both low and high speed regimes. It is found that siphon flow introduces a linear spatially dependent phase sh… ▽ More There are ubiquitous flows observed in the solar atmosphere of sub-Alfvénic speeds, however after flaring and coronal mass ejection events flows can become Alfvénic. In this Letter, we derive an expression for the standing kink mode frequency due to siphon flow in coronal loops, valid for both low and high speed regimes. It is found that siphon flow introduces a linear spatially dependent phase shift along coronal loops and asymmetric eigenfunctions. We demonstrate how this theory can be used to determine the kink and flow speed of oscillating coronal loops with reference to an observational case study. It is shown that the presence of siphon flow can cause the underestimation of magnetic field strength in coronal loops using the traditional seismological methods. △ Less

Submitted 27 January, 2011; originally announced January 2011.

On the linear coupling between fast and slow MHD waves due to line-tying effects

Authors: J. Terradas, J. Andries, E. Verwichte

Abstract: Oscillations in coronal loops are usually interpreted in terms of uncoupled magnetohydrodynamic (MHD) waves. Examples of these waves are standing transverse motions, interpreted as the kink MHD modes, and propagating slow modes, commonly reported at the loop footpoints. Here we study a simple system in which fast and slow MHD waves are coupled. The goal is to understand the fingerprints of the cou… ▽ More Oscillations in coronal loops are usually interpreted in terms of uncoupled magnetohydrodynamic (MHD) waves. Examples of these waves are standing transverse motions, interpreted as the kink MHD modes, and propagating slow modes, commonly reported at the loop footpoints. Here we study a simple system in which fast and slow MHD waves are coupled. The goal is to understand the fingerprints of the coupling when boundary conditions are imposed in the model. The reflection problem of a fast and slow MHD wave interacting with a rigid boundary, representing the line-tying effect of the photosphere, is analytically investigated. Both propagating and standing waves are analysed and the time-dependent problem of the excitation of these waves is considered. An obliquely incident fast MHD wave on the photosphere inevitably generates a slow mode. The frequency of the generated slow mode at the photosphere is exactly the same as the frequency of the incident fast MHD mode, but its wavelength is much smaller, assuming that the sound speed is smaller than the Alfvén speed. The main signatures of the generated slow wave are density fluctuations at the loop footpoints. We have derived a simple formula that relates the velocity amplitude of the transverse standing mode with the density enhancements at the footpoints due to the driven slow modes. Using these results it is shown that there are possible evidences in the observations of the coupling between these two modes. △ Less

Submitted 3 November, 2010; originally announced November 2010.