Source author record

M. Luna

M. Luna appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2022arXiv

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

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.

preprint2016arXiv

The effects of magnetic-field geometry on longitudinal oscillations of solar prominences: Cross-sectional area variation for thin tubes

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.

preprint2015arXiv

Charge exchange in fluid description of partially ionized plasmas

The effects of charge exchange on waves propagating in weakly ionized plasmas are discussed. It is shown that for low-frequency processes, ions and neutrals should be treated as a single fluid with some effective charge on all of them. We have derived a new momentum equation which should be used in such an environment. As a result, the low-frequency magnetic waves can propagate even if particles are not magnetized, which is entirely due to the charge exchange and the fact that it is not possible to separate particles into two different populations as charged and neutral species. So there can be no friction force between ions and neutrals in the usual sense. The mean force per particle is proportional to the ionization ratio $n_i/(n_i+ n_n)$. Regarding the application of the theory to the Alfven wave propagation in the lower solar atmosphere, the results predict that the plane of displacement of the fluid must change by 90 degrees when an Alfven wave propagates from the area where particles are un-magnetized (photosphere) to the area where they are magnetized (chromosphere). Because of the most accurate cross sections which we have here, it is possible to very accurately determine altitudes at which such rotation of the Alfven wave takes place.

preprint2015arXiv

On the nature of transverse coronal waves revealed by wavefront dislocations

Coronal waves are an important aspect of the dynamics of the plasma in the corona. Wavefront dislocations are topological features of most waves in nature and also of magnetohydrodynamic waves. Are there dislocations in coronal waves? The finding and explanation of dislocations may shed light on the nature and characteristics of the propagating waves, their interaction in the corona and in general on the plasma dynamics. We positively identify dislocations in coronal waves observed by the Coronal Multi-channel Polarimeter (CoMP) as singularities in the Doppler shifts of emission coronal lines. We study the possible singularities that can be expected in coronal waves and try to reproduce the observed dislocations in terms of localization and frequency of appearance. The observed dislocations can only be explained by the interference of a kink and a sausage wave modes propagating with different frequencies along the coronal magnetic field. In the plane transverse to the propagation, the cross-section of the oscillating plasma must be smaller than the spatial resolution, and the two waves result in net longitudinal and transverse velocity components that are mixed through projection onto the line of sight. Alfvén waves can be responsible of the kink mode, but a magnetoacoustic sausage mode is necessary in all cases. Higher (flute) modes are excluded. The kink mode has a pressure amplitude that is smaller than the pressure amplitude of the sausage mode, though its observed velocity is larger. This concentrates dislocations on the top of the loop. To explain dislocations, any model of coronal waves must include the simultaneous propagation and interference of kink and sausage wave modes of comparable but different frequencies, with a sausage wave amplitude much smaller than the kink one.

preprint2015arXiv

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

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.

preprint2015arXiv

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

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.

preprint2014arXiv

Morphology and dynamics of solar prominences from 3D MHD simulations

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.

preprint2014arXiv

Observations and Implications of Large-Amplitude Longitudinal Oscillations in a Solar Filament

On 20 August 2010 an energetic disturbance triggered large-amplitude longitudinal oscillations in a nearby filament. The triggering mechanism appears to be episodic jets connecting the energetic event with the filament threads. In the present work we analyze this periodic motion in a large fraction of the filament to characterize the underlying physics of the oscillation as well as the filament properties. The results support our previous theoretical conclusions that the restoring force of large-amplitude longitudinal oscillations is solar gravity, and the damping mechanism is the ongoing accumulation of mass onto the oscillating threads. Based on our previous work, we used the fitted parameters to determine the magnitude and radius of curvature of the dipped magnetic field along the filament, as well as the mass accretion rate onto the filament threads. These derived properties are nearly uniform along the filament, indicating a remarkable degree of cohesiveness throughout the filament channel. Moreover, the estimated mass accretion rate implies that the footpoint heating responsible for the thread formation, according to the thermal nonequilibrium model, agrees with previous coronal heating estimates. We estimate the magnitude of the energy released in the nearby event by studying the dynamic response of the filament threads, and discuss the implications of our study for filament structure and heating.

preprint2014arXiv

Rayleigh-Taylor instability in prominences from numerical simulations including partial ionization effects

We study the Rayleigh-Taylor instability (RTI) at a prominence-corona transition region in a non-linear regime. Our aim is to understand how the presence of neutral atoms in the prominence plasma influences the instability growth rate, and the evolution of velocity, magnetic field vector and thermodynamic parameters of turbulent drops. We perform 2.5D numerical simulations of the instability initiated by a multi-mode perturbation at the corona-prominence interface using a single-fluid MHD approach including a generalized Ohm's law. The initial equilibrium configuration is purely hydrostatic and contains a homogeneous horizontal magnetic field forming an angle with the direction in which the plasma is perturbed. We analyze simulations with two different orientations of the magnetic field. For each field orientation we compare two simulations, one for the pure MHD case, and one including the ambipolar diffusion in the Ohm's law (AD case). Other than that, both simulations for each field orientation are identical. The numerical results in the initial stage of the instability are compared with the analytical linear calculations. We find that the configuration is always unstable in the AD case. The growth rate of the small-scale modes in the non-linear regime is up to 50% larger in the AD case than in the purely MHD case and the average velocities of flows are a few percent larger. Significant drift momenta are found at the interface between the coronal and the prominence material at all stages of the instability, produced by the faster downward motion of the neutral component with respect to the ionized component. The differences in temperature of the bubbles between the ideal and non-ideal case are also significant, reaching 30%. There is an asymmetry between large rising bubbles and small-scale down flowing fingers, favoring the detection of upward velocities in observations.

preprint2013arXiv

Observational Study of Large Amplitude Longitudinal Oscillations in a Solar Filament

On 20 August 2010 an energetic disturbance triggered damped large-amplitude longitudinal (LAL) oscillations in almost an entire filament. In the present work we analyze this periodic motion in the filament to characterize the damping and restoring mechanism of the oscillation. Our method involves placing slits along the axis of the filament at different angles with respect to the spine of the filament, finding the angle at which the oscillation is clearest, and fitting the resulting oscillation pattern to decaying sinusoidal and Bessel functions. These functions represent the equations of motion of a pendulum damped by mass accretion. With this method we determine the period and the decaying time of the oscillation. Our preliminary results support the theory presented by Luna and Karpen (2012) that the restoring force of LAL oscillations is solar gravity in the tubes where the threads oscillate, and the damping mechanism is the ongoing accumulation of mass onto the oscillating threads. Following an earlier paper, we have determined the magnitude and radius of curvature of the dipped magnetic flux tubes hosting a thread along the filament, as well as the mass accretion rate of the filament threads, via the fitted parameters.

preprint2013arXiv

Propagating Waves Transverse to the Magnetic Field in a Solar Prominence

We report an unusual set of observations of waves in a large prominence pillar which consist of pulses propagating perpendicular to the prominence magnetic field. We observe a huge quiescent prominence with the Solar Dynamics Observatory (SDO) Atmospheric Imaging Assembly (AIA) in EUV on 2012 October 10 and only a part of it, the pillar, which is a foot or barb of the prominence, with the Hinode Solar Optical Telescope (SOT) (in Ca II and Hαlines), Sac Peak (in Hα, Hβ and Na-D lines), THEMIS ("Télescope Héliographique pour l' Etude du Magnétisme et des Instabilités Solaires") with the MTR (MulTi-Raies) spectropolarimeter (in He D_3 line). The THEMIS/MTR data indicates that the magnetic field in the pillar is essentially horizontal and the observations in the optical domain show a large number of horizontally aligned features on a much smaller scale than the pillar as a whole. The data is consistent with a model of cool prominence plasma trapped in the dips of horizontal field lines. The SOT and Sac Peak data over the 4 hour observing period show vertical oscillations appearing as wave pulses. These pulses, which include a Doppler signature, move vertically, perpendicular to the field direction, along thin quasi-vertical columns in the much broader pillar. The pulses have a velocity of propagation of about 10 km/s, a period about 300 sec, and a wavelength around 2000 km. We interpret these waves in terms of fast magneto-sonic waves and discuss possible wave drivers.

preprint2013arXiv

Rayleigh-Taylor instability in partially ionized prominence plasma

We study Rayleigh-Taylor instability (RTI) at the coronal-prominence boundary by means of 2.5D numerical simulations in a single-fluid MHD approach including a generalized Ohm's law. The initial configuration includes a homogeneous magnetic field forming an angle with the direction in which the plasma is perturbed. For each field inclination we compare two simulations, one for the pure MHD case, and one including the ambipolar diffusion in the Ohm's law, otherwise identical. We find that the configuration containing neutral atoms is always unstable. The growth rate of the small-scale modes in the non-linear regime is larger than in the purely MHD case.

preprint2012arXiv

Formation and evolution of a multi-threaded prominence

We investigate the process of formation and subsequent evolution of prominence plasma in a filament channel and its overlying arcade. We construct a three-dimensional time-dependent model of an intermediate quiescent prominence. We combine the magnetic field structure with one-dimensional independent simulations of many flux tubes, of a three-dimensional sheared double arcade, in which the thermal nonequilibrium process governs the plasma evolution. We have found that the condensations in the corona can be divided into two populations: threads and blobs. Threads are massive condensations that linger in the field line dips. Blobs are ubiquitous small condensations that are produced throughout the filament and overlying arcade magnetic structure, and rapidly fall to the chromosphere. The threads are the principal contributors to the total mass. The total prominence mass is in agreement with observations, assuming a reasonable filling factor. The motion of the threads is basically horizontal, while blobs move in all directions along the field. The peak velocities for both populations are comparable. We have generated synthetic images of the whole structure in an H$α$ proxy and in two EUV channels of the AIA instrument aboard SDO, thus showing the plasma at cool, warm, and hot temperatures. The predicted differential emission measure of our system agrees very well with observations. We conclude that the sheared-arcade magnetic structure and plasma behavior driven by thermal nonequilibrium fit well the abundant observational evidence for typical intermediate prominences.

preprint2012arXiv

The effects of magnetic-field geometry on longitudinal oscillations of solar prominences

We investigate the influence of the geometry of the solar filament magnetic structure on the large-amplitude longitudinal oscillations. A representative filament flux tube is modeled as composed of a cool thread centered in a dipped part with hot coronal regions on either side. We have found the normal modes of the system, and establish that the observed longitudinal oscillations are well described with the fundamental mode. For small and intermediate curvature radii and moderate to large density contrast between the prominence and the corona, the main restoring force is the solar gravity. In this full wave description of the oscillation a simple expression for the oscillation frequencies is derived in which the pressure-driven term introduces a small correction. We have also found that the normal modes are almost independent of the geometry of the hot regions of the tube. We conclude that observed large-amplitude longitudinal oscillations are driven by the projected gravity along the flux tubes, and are strongly influenced by the curvature of the dips of the magnetic field in which the threads reside.

preprint2010arXiv

Transverse oscillations of a multi-stranded loop

We investigate the transverse oscillations of a line-tied multi-stranded coronal loop composed of several parallel cylindrical strands. First, the collective fast normal modes of the loop are found with the T-matrix theory. There is a huge quantity of normal modes with very different frequencies and a complex structure of the associated magnetic pressure perturbation and velocity field. The modes can be classified as bottom, middle, and top according to their frequencies and spatial structure. Second, the temporal evolution of the velocity and magnetic pressure perturbation after an initial disturbance are analyzed. We find complex motions of the strands. The frequency analysis reveals that these motions are a combination of low and high frequency modes. The complexity of the strand motions produces a strong modulation of the whole tube movement. We conclude that the presumed internal fine structure of a loop influences its transverse oscillations and so its transverse dynamics cannot be properly described by those of an equivalent monolithic loop.

preprint2008arXiv

Transverse oscillations of systems of coronal loops

We study the collective kinklike normal modes of a system of several cylindrical loops using the T-matrix theory. Loops that have similar kink frequencies oscillate collectively with a frequency which is slightly different from that of the individual kink mode. On the other hand, if the kink frequency of a loop is different from that of the others, it oscillates individually with its own frequency. Since the individual kink frequency depends on the loop density but not on its radius for typical 1 MK coronal loops, a coupling between kink oscillations of neighboring loops take place when they have similar densities. The relevance of these results in the interpretation of the oscillations studied by \citet{schrijver2000} and \citet{verwichte2004}, in which transverse collective loop oscillations seem to be detected, is discussed. In the first case, two loops oscillating in antiphase are observed; interpreting this motion as a collective kink mode suggests that their densities are roughly equal. In the second case, there are almost three groups of tubes that oscillate with similar periods and therefore their dynamics can be collective, which again seems to indicate that the loops of each group share a similar density. All the other loops seem to oscillate individually and their densities can be different from the rest.