Source author record

M. Collados

M. Collados 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

38works
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

38 published item(s)

preprint2022arXiv

First light of the Integral Field Unit of GRIS on the GREGOR solar telescope

An Integral Field Unit (IFU) based on image-slicers has been added to the GREGOR Infrared Spectrograph (GRIS). This upgrade to the instrument makes possible 2D spectropolarimetry in the near-infrared by simultaneously recording the full Stokes profiles of spectral lines (in a given spectral interval) at all the points in the field of view. It provides high-cadence spectropolarimetric observations at the instrument's high spatial resolution and high polarization sensitivity at the GREGOR solar telescope. The IFU is ideal for observing the polarized spectrum of fast-evolving solar features at high spatial and spectral resolutions. The high observing cadence opens the possibility of time-series observations. The analysis of observations to this level of accuracy is essential for understanding the complex dynamics and interactions of solar plasma and magnetic fields. The image slicer of the IFU has eight slices of width 100 micron, covering a total field of view of 6" x 3". It was designed and built within the framework of the European projects SOLARNET and GREST, as a prototype for future instruments of the European Solar Telescope (EST) and was integrated into GRIS. After two commissioning campaigns in 2017 and 2018, the IFU was finally installed at the end of September 2018 and offered to all observers who use the telescope.

preprint2022arXiv

Generalized Fluid Models of the Braginskii Type

Several generalizations of the well-known fluid model of Braginskii (Rev. of Plasma Phys., 1965) are considered. We use the Landau collisional operator and the moment method of Grad. We focus on the 21-moment model that is analogous to the Braginskii model, and we also consider a 22-moment model. Both models are formulated for general multi-species plasmas with arbitrary masses and temperatures, where all the fluid moments are described by their evolution equations. The 21-moment model contains two "heat flux vectors" (3rd and 5th-order moments) and two "viscosity-tensors" (2nd and 4th-order moments). The Braginskii model is then obtained as a particular case of a one ion-electron plasma with similar temperatures, with de-coupled heat fluxes and viscosity-tensors expressed in a quasi-static approximation. We provide all the numerical values of the Braginskii model in a fully analytic form (together with the 4th and 5th-order moments). For multi-species plasmas, the model makes calculation of transport coefficients straightforward. Formulation in fluid moments (instead of Hermite moments) is also suitable for implementation into existing numerical codes. It is emphasized that it is the quasi-static approximation which makes some Braginskii coefficients divergent in a weakly-collisional regime. Importantly, we show that the heat fluxes and viscosity-tensors are coupled even in the linear approximation, and that the fully contracted (scalar) perturbations of the 4th-order moment, which are accounted for in the 22-moment model, modify the energy exchange rates. We also provide several Appendices, which can be useful as a guide for deriving the Braginskii model with the moment method of Grad.

preprint2020arXiv

Evolution of Stokes V area asymmetry related to a quiet Sun cancelation observed with GRIS/IFU

A quiet Sun magnetic flux cancelation event at the disk centre was recorded using the Integral Field Unit (IFU) mounted on the GREGOR Infrared Spectrograph (GRIS). GRIS sampled the event in the photospheric Si I 10827 Å spectral line. The cancelation is preceded by a significant rise in line core intensity and excitation temperature, which is inferred from Stokes inversions under local thermodynamic equilibrium (LTE). The opposite polarity features seem to undergo reconnection above the photosphere. We also found that the border pixels neighboring the polarity inversion line of one of the polarities exhibit a systematic variation of area asymmetry. Area asymmetry peaks right after the line core intensity enhancement and gradually declines thereafter. Analyzing Stokes profiles recorded from either side of the polarity inversion line could therefore potentially provide additional information on the reconnection process related to magnetic flux cancelation. Further analysis without assuming LTE will be required to fully characterize this event.

preprint2020arXiv

Influence of ambipolar and Hall effects on vorticity in 3D simulations of magneto-convection

This paper presents the results of the analysis of 3D simulations of solar magneto-convection that include the joint action of the ambipolar diffusion and the Hall effect. Three simulation-runs are compared: one including both ambipolar diffusion and Hall effect; one including only ambipolar diffusion; and one without any of these two effects. The magnetic field is amplified from initial field to saturation level by the action of turbulent local dynamo. In each of these cases, we study 2 hours of simulated solar time after the local dynamo reaches the saturation regime. We analyze the power spectra of vorticity, of magnetic field fluctuations and of the different components of the magnetic Poynting flux responsible for the transport of vertical or horizontal perturbations. Our preliminary results show that the ambipolar diffusion produces a strong reduction of vorticity in the upper chromospheric layers and that it dissipates the vortical perturbations converting them into thermal energy. The Hall effect acts in the opposite way, strongly enhancing the vorticity. When the Hall effect is included, the magnetic field in the simulations becomes, on average, more vertical and long-lived flux tube-like structures are produced. We trace a single magnetic structure to study its evolution pattern and the magnetic field intensification, and their possible relation to the Hall effect.

preprint2020arXiv

Photospheric magnetic topology of a north polar region

Aims. We aim to characterise the magnetism of a large fraction of the north polar region close to a maximum of activity, when the polar regions are reversing their dominant polarity. Methods. We make use of full spectropolarimetric data from the CRisp Imaging Spectro-Polarimeter installed at the Swedish Solar Telescope. The data consist of a photospheric spectral line, which is used to infer the various physical parameters of different quiet Sun regions by means of the solution of the radiative transfer equation. We focus our analysis on the properties found for the north polar region and their comparison to the same analysis applied to data taken at disc centre and low-latitude quiet Sun regions for reference. We also analyse the spatial distribution of magnetic structures throughout the north polar region. Results. We find that the physical properties of the polar region (line-of-sight velocity, magnetic flux, magnetic inclination and magnetic azimuth) are compatible with those found for the quiet Sun at disc centre and are similar to the ones found at low latitudes close to the limb. Specifically, the polar region magnetism presents no specific features. The structures for which the transformation from a line-of-sight to a local reference frame was possible harbour large magnetic fluxes ($>10^{17}$ Mx) and are in polarity imbalance with a dominant positive polarity, the largest ones ($>10^{19}$ Mx) being located below $73^{\circ}$ latitude.

preprint2020arXiv

Tracking downflows from the chromosphere to the photosphere in a solar arch filament system

We study the dynamics of plasma along the legs of an arch filament system (AFS) from the chromosphere to the photosphere, observed with high-cadence spectroscopic data from two ground-based solar telescopes: the GREGOR telescope (Tenerife) using the GREGOR Infrarred Spectrograph (GRIS) in the He I 10830 Å range and the Swedish Solar Telescope (La Palma) using the CRisp Imaging Spectro-Polarimeter to observe the Ca II 8542 Å and Fe I 6173 Å spectral lines. The temporal evolution of the draining of the plasma was followed along the legs of a single arch filament from the chromosphere to the photosphere. The average Doppler velocities inferred at the upper chromosphere from the He I 10830 Å triplet reach velocities up to 20-24~km~s$^{-1}$, in the lower chromosphere and upper photosphere the Doppler velocities reach up to 11~km~s$^{-1}$ and 1.5~km~s$^{-1}$ in the case of the Ca II 8542 Å and Si I 10827 Å spectral lines, respectively. The evolution of the Doppler velocities at different layers of the solar atmosphere (chromosphere and upper photosphere) shows that they follow the same LOS velocity pattern, which confirm the observational evidence that the plasma drains towards the photosphere as proposed in models of AFSs. The Doppler velocity maps inferred from the lower photospheric Ca I 10839 Å or Fe I 6173 Å spectral lines do not show the same LOS velocity pattern. Thus, there is no evidence that the plasma reaches the lower photosphere. The observations and the nonlinear force-free field extrapolations demonstrate that the magnetic field loops of the AFS rise with time. We found flow asymmetries at different footpoints of the AFS. The NLFFF values of the magnetic field strength give us a clue to explain these flow asymmetries.

preprint2019arXiv

An introductory guide to fluid models with anisotropic temperatures Part 1 -- CGL description and collisionless fluid hierarchy

We present a detailed guide to advanced collisionless fluid models that incorporate kinetic effects into the fluid framework, and that are much closer to the collisionless kinetic description than traditional magnetohydrodynamics. Such fluid models are directly applicable to modeling turbulent evolution of a vast array of astrophysical plasmas, such as the solar corona and the solar wind, the interstellar medium, as well as accretion disks and galaxy clusters. The text can be viewed as a detailed guide to Landau fluid models and it is divided into two parts. Part 1 is dedicated to fluid models that are obtained by closing the fluid hierarchy with simple (non Landau fluid) closures. Part 2 is dedicated to Landau fluid closures. Here in Part 1, we discuss the CGL fluid model in great detail, together with fluid models that contain dispersive effects introduced by the Hall term and by the finite Larmor radius (FLR) corrections to the pressure tensor. We consider dispersive effects introduced by the non-gyrotropic heat flux vectors. We investigate the parallel and oblique firehose instability, and show that the non-gyrotropic heat flux strongly influences the maximum growth rate of these instabilities. Furthermore, we discuss fluid models that contain evolution equations for the gyrotropic heat flux fluctuations and that are closed at the 4th-moment level by prescribing a specific form for the distribution function. For the bi-Maxwellian distribution, such a closure is known as the "normal" closure. We also discuss a fluid closure for the bi-kappa distribution. Finally, by considering one-dimensional Maxwellian fluid closures at higher-order moments, we show that such fluid models are always unstable. The last possible non Landau fluid closure is therefore the "normal" closure, and beyond the 4th-order moment, Landau fluid closures are required.

preprint2019arXiv

An introductory guide to fluid models with anisotropic temperatures Part 2 -- Kinetic theory, Padé approximants and Landau fluid closures

In Part 2 of our guide to collisionless fluid models, we concentrate on Landau fluid closures. These closures were pioneered by Hammett and Perkins and allow for the rigorous incorporation of collisionless Landau damping into a fluid framework. It is Landau damping that sharply separates traditional fluid models and collisionless kinetic theory, and is the main reason why the usual fluid models do not converge to the kinetic description, even in the long-wavelength low-frequency limit. We start with a brief introduction to kinetic theory, where we discuss in detail the plasma dispersion function $Z(ζ)$, and the associated plasma response function $R(ζ)=1+ζZ(ζ)=-Z'(ζ)/2$. We then consider a 1D (electrostatic) geometry and make a significant effort to map all possible Landau fluid closures that can be constructed at the 4th-order moment level. These closures for parallel moments have general validity from the largest astrophysical scales down to the Debye length, and we verify their validity by considering examples of the (proton and electron) Landau damping of the ion-acoustic mode, and the electron Landau damping of the Langmuir mode. We proceed by considering 1D closures at higher-order moments than the 4th-order, and as was concluded in Part 1, this is not possible without Landau fluid closures. We show that it is possible to reproduce linear Landau damping in the fluid framework to any desired precision, thus showing the convergence of the fluid and collisionless kinetic descriptions. We then consider a 3D (electromagnetic) geometry in the gyrotropic (long-wavelength low-frequency) limit and map all closures that are available at the 4th-order moment level. In the Appendix A, we provide comprehensive tables with Padé approximants of $R(ζ)$ up to the 8th-pole order, with many given in an analytic form.

preprint2016arXiv

Active region fine structure observed at 0.08 arcsec resolution

The various mechanisms of magneto-convective energy transport determines the structure of sunspots and active regions. We characterise the appearance of light bridges and other fine-structure details and elaborate on their magneto-convective nature. We present speckle-reconstructed images taken with the broad-band imager at the 1.5 m GREGOR telescope in the 486nm and 589nm bands. We estimate the spatial resolution from the noise characteristics of the image bursts and obtain 0.08" at 589nm. We describe structure details in individual best images as well as the temporal evolution of selected features. We find branched dark lanes extending along thin (~1") light bridges in sunspots at various heliocentric angles. In thick (~2") light bridges the branches are disconnected from the central lane and have a Y shape with a bright grain toward the umbra. The images reveal that light bridges exist on varying intensity levels and that their small-scale features evolve on timescales of minutes. Faint light bridges show dark lanes outlined by the surrounding bright features. Dark lanes are very common and are also found in the boundary of pores. They have a characteristic width of 0.1" or smaller. Intergranular dark lanes of that width are seen in active region granulation. While central dark lanes in thin light bridges are elevated and associated with an density increase above upflows, the dark lane branches correspond to locations of downflows and are depressed relative to the adjacent bright plasma. Thick light bridges with central dark lanes show no projection effect. They have a flat elevated plateau that falls off steeply at the umbral boundary. There, Y-shaped filaments form as they do in the inner penumbra. This indicates the presence of inclined magnetic fields, meaning that the umbral magnetic field is wrapped around the convective light bridge.

preprint2016arXiv

Deep probing of the photospheric sunspot penumbra: no evidence for magnetic field-free gaps

Some models for the topology of the magnetic field in sunspot penumbrae predict the existence of field-free or dynamically weak-field regions in the deep Photosphere. To confirm or rule out the existence of weak-field regions in the deepest photospheric layers of the penumbra. The magnetic field at $\logτ_5=0$ is investigated by means of inversions of spectropolarimetric data of two different sunspots located very close to disk center with a spatial resolution of approximately 0.4-0.45 arcsec. The data have been recorded using the GRIS instrument attached to the 1.5-meters GREGOR solar telescope at El Teide observatory. It includes three Fe I lines around 1565 nm, whose sensitivity to the magnetic field peaks at half a pressure-scale-height deeper than the sensitivity of the widely used Fe I spectral line pair at 630 nm. Prior to the inversion, the data is corrected for the effects of scattered light using a deconvolution method with several point spread functions. At $\logτ_5=0$ we find no evidence for the existence of regions with dynamically weak ($B<500$ Gauss) magnetic fields in sunspot penumbrae. This result is much more reliable than previous investigations done with Fe I lines at 630 nm. Moreover, the result is independent of the number of nodes employed in the inversion, and also independent of the point spread function used to deconvolve the data, and does not depend on the amount of straylight (i.e. wide-angle scattered light) considered.

preprint2016arXiv

Fitting peculiar spectral profiles in He I 10830 Å absorption features

The new generation of solar instruments provides better spectral, spatial, and temporal resolution for a better understanding of the physical processes that take place on the Sun. Multiple-component profiles are more commonly observed with these instruments. Particularly, the He I 10830 Å triplet presents such peculiar spectral profiles, which give information on the velocity and magnetic fine structure of the upper chromosphere. The purpose of this investigation is to describe a technique to efficiently fit the two blended components of the He I 10830 Å triplet, which are commonly observed when two atmospheric components are located within the same resolution element. The observations used in this study were taken on 2015 April 17 with the very fast spectroscopic mode of the GREGOR Infrared Spectrograph (GRIS) attached to the 1.5-meter GREGOR solar telescope, located at the Observatorio del Teide, Tenerife, Spain. We apply a double-Lorentzian fitting technique using Levenberg-Marquardt least-squares minimization. This technique is very simple and much faster than inversion codes. Line-of-sight Doppler velocities can be inferred for a whole map of pixels within just a few minutes. Our results show sub- and supersonic downflow velocities of up to 32 km/s for the fast component in the vicinity of footpoints of filamentary structures. The slow component presents velocities close to rest.

preprint2016arXiv

Flows in and around active region NOAA12118 observed with the GREGOR solar telescope and SDO/HMI

Accurate measurements of magnetic and velocity fields in and around solar active regions are key to unlocking the mysteries of the formation and the decay of sunspots. High spatial resolution image and spectral sequences with a high cadence obtained with the GREGOR solar telescope give us an opportunity to scrutinize 3-D flow fields with local correlation tracking and imaging spectroscopy. We present GREGOR early science data acquired in 2014 July - August with the GREGOR Fabry-Pérot Interferometer and the Blue Imaging Channel. Time-series of blue continuum (? 450.6 nm) images of the small active region NOAA 12118 were restored with the speckle masking technique to derive horizontal proper motions and to track the evolution of morphological changes. In addition, high-resolution observations are discussed in the context of synoptic data from the Solar Dynamics Observatory.

preprint2016arXiv

Horizontal flow fields in and around a small active region-- The transition period between flux emergence and decay

Aims. Combining high-resolution and synoptic observations aims to provide a comprehensive description of flux emergence at photospheric level and of the growth process that eventually leads to a mature active region. Methods. Small active region NOAA 12118 was observed on 2014 July 18 with the 1.5-meter GREGOR solar telescope on 2014 July 18. High-resolution time-series of blue continuum and G-band images acquired in the blue imaging channel (BIC) of the GREGOR Fabry-Pérot Interferometer (GFPI) were complemented by LOS magnetograms and continuum images obtained with the HMI onboard the SDO. Horizontal proper motions and horizontal plasma velocities were computed with local correlation tracking (LCT) and the differential affine velocity estimator, respectively. Morphological image processing was employed to measure the photometric/magnetic area, magnetic flux, and the separation profile of the EFR during its evolution. Results. The computed growth rates for photometric area, magnetic area, and magnetic flux are about twice as high as the respective decay rates. The space-time diagram using HMI magnetograms of five days traces a leaf-like structure, which is determined by the initial separation of the two polarities, a rapid expansion phase, a time when the spread stalls, and a period when the region slowly shrinks again. The separation rate of 0.26 km\s is highest in the initial stage, and it decreases when the separation comes to a halt. Horizontal plasma velocities computed at four evolutionary stages indicate a changing pattern of inflows. In LCT maps we find persistent flow patterns such as outward motions in the outer part of the two major pores, a diverging feature near the trailing pore marking the site of upwelling plasma and flux emergence, and low velocities in the interior of pores. We detected many elongated rapidly expanding granules between the two major polarities.

preprint2016arXiv

Magnetic fields of opposite polarity in sunspot penumbrae

Context. A significant part of the penumbral magnetic field returns below the surface in the very deep photosphere. For lines in the visible, a large portion of this return field can only be detected indirectly by studying its imprints on strongly asymmetric and three-lobed Stokes V profiles. Infrared lines probe a narrow layer in the very deep photosphere, providing the possibility of directly measuring the orientation of magnetic fields close to the solar surface. Aims. We study the topology of the penumbral magnetic field in the lower photosphere, focusing on regions where it returns below the surface. Methods. We analyzed 71 spectropolarimetric datasets from Hinode and from the GREGOR infrared spectrograph. We inferred the quality and polarimetric accuracy of the infrared data after applying several reduction steps. Techniques of spectral inversion and forward synthesis were used to test the detection algorithm. We compared the morphology and the fractional penumbral area covered by reversed-polarity and three-lobed Stokes V profiles for sunspots at disk center. We determined the amount of reversed-polarity and three-lobed Stokes V profiles in visible and infrared data of sunspots at various heliocentric angles. From the results, we computed center-to-limb variation curves, which were interpreted in the context of existing penumbral models. Results. Observations in visible and near-infrared spectral lines yield a significant difference in the penumbral area covered by magnetic fields of opposite polarity. In the infrared, the number of reversed-polarity Stokes V profiles is smaller by a factor of two than in the visible. For three-lobed Stokes V profiles the numbers differ by up to an order of magnitude.

preprint2016arXiv

On the magnetism and dynamics of prominence legs hosting tornadoes

Solar tornadoes are dark vertical filamentary structures observed in the extreme ultraviolet associated with prominence legs and filament barbs. Their true nature and relationship to prominences requires understanding their magnetic structure and dynamic properties. Recently, a controversy has arisen: is the magnetic field organized forming vertical, helical structures or is it dominantly horizontal? And concerning their dynamics, are tornadoes really rotating or is it just a visual illusion? Here, we analyze four consecutive spectropolarimetric scans of a prominence hosting tornadoes on its legs which help us shed some light on their magnetic and dynamical properties. We show that the magnetic field is very smooth in all the prominence, probably an intrinsic property of the coronal field. The prominence legs have vertical helical fields that show slow temporal variation probably related to the motion of the fibrils. Concerning the dynamics, we argue that 1) if rotation exists, it is intermittent, lasting no more than one hour, and 2) the observed velocity pattern is also consistent with an oscillatory velocity pattern (waves).

preprint2016arXiv

Three-dimensional structure of a sunspot light bridge

Active regions are the most prominent manifestations of solar magnetic fields; their generation and dissipation are fundamental problems in solar physics. Light bridges are commonly present during sunspot decay, but a comprehensive picture of their role in the removal of photospheric magnetic field is still missing. We study the three dimensional configuration of a sunspot and in particular its light bridge during one of the last stages of its decay. We present the magnetic and thermodynamical stratification inferred from full Stokes inversions of the photospheric Si I 10827 Å and Ca I 10839 Å lines obtained with the GREGOR Infrared Spectrograph of the GREGOR telescope at Observatorio del Teide, Tenerife, Spain. The analysis is complemented by a study of continuum images covering the disk passage of the active region, which are provided by the Helioseismic and Magnetic Imager on board the Solar Dynamics Observatory. The sunspot shows a light bridge with penumbral continuum intensity that separates the central umbra from a smaller umbra. We find that in this region the magnetic field lines form a canopy with lower magnetic field strength in the inner part. The photospheric light bridge is dominated by gas pressure (high-$β$), as opposed to the surrounding umbra where the magnetic pressure is higher. A convective flow is observed in the light bridge. This flow is able to bend the magnetic field lines and to produce field reversals. The field lines close above the light bridge and become as vertical and strong as in the surrounding umbra. We conclude that it develops because of two highly magnetized regions which come closer during the sunspot evolution.

preprint2016arXiv

Upper Chromospheric Magnetic Field of a Sunspot Penumbra: Observations of Fine Structure

The fine-structure of magnetic field of a sunspot penumbra in the upper chromosphere is to be explored and compared to that in the photosphere. High spatial resolution spectropolarimetric observations were recorded with the 1.5-meter GREGOR telescope using the GREGOR Infrared Spectrograph (GRIS). The observed spectral domain includes the upper chromospheric He I triplet at 1083.0 nm and the photospheric Si I 1082.7 nm and Ca I 1083.3 nm spectral lines. The upper chromospheric magnetic field is obtained by inverting the He I triplet assuming a Milne-Eddington type model atmosphere. A height dependent inversion was applied to the Si I 1082.7 nm and Ca I 1083.3 nm lines to obtain the photospheric magnetic field. We find that the inclination of the magnetic field shows variations in the azimuthal direction both in the photosphere, but also in the upper chromosphere. The chromospheric variations remarkably well coincide with the variations in the inclination of the photospheric field and resemble the well-known spine and inter-spine structure in the photospheric layers of penumbrae. The typical peak-to-peak variations in the inclination of the magnetic field in the upper chromosphere is found to be 10-15 degree, i.e., roughly half the variation in the photosphere. In contrast, the magnetic field strength of the observed penumbra does not show variations on small spatial scales in the upper chromosphere. Thanks to the high spatial resolution observations possible with the GREGOR telescope at 1.08 microns, we find that the prominent small-scale fluctuations in the magnetic field inclination, which are a salient part of the property of sunspot penumbral photospheres, also persist in the chromosphere, although at somewhat reduced amplitudes. Such a complex magnetic configuration may facilitate penumbral chromospheric dynamic phenomena, such as penumbral micro-jets or transient bright dots.

preprint2015arXiv

Evershed flow observed in neutral and singly ionized iron lines

The amplitudes of the Evershed flow are measured using pairs of carefully selected FeI and FeII spectral lines located close in wavelength and registered simultaneously. A sunspot belonging to the NOAA 11582 group was scanned using the spectrograph of the German Vacuum Tower Telescope (Observatorio del Teide, Tenerife). Velocities were extracted from intensity profiles using the lambda-meter technique. The formation heights of the observed spectral lines were calculated using semi-empirical models of a bright and dark penumbral filament taking into account the sunspot location at the limb. Our objective is to compare azimuthally averaged amplitudes of the Evershed flow extracted from neutral and ion lines. We find measurable differences in the radial component of the flow. All five pairs of lines show the same tendency, with a few hundred m/s larger amplitude of the flow measured from FeI lines compared to FeII lines. This tendency is preserved at all photospheric heights and radial distances in the penumbra. We discuss the possible origin of this effect.

preprint2015arXiv

Magnetic and dynamical photospheric disturbances observed during an M3.2 solar flare

This letter reports on a set of full-Stokes spectropolarimetric observations in the near infrared He I 10830 A spectral region covering the pre-, flare, and post-flare phases of an M3.2 class solar flare. The flare originated on 2013 May 17 and belonged to active region NOAA 11748. We detected strong He I 10830 A emission in the flare. The red component of the He I triplet peaks at an intensity ratio to the continuum of about 1.86. During the flare, He I Stokes V is substantially larger and appears reversed compared to the usually larger Si I Stokes V profile. The photospheric Si I inversions of the four Stokes profiles reveal the following: (1) the magnetic field strength in the photosphere decreases or is even absent during the flare phase, as compared to the pre-flare phase. However, this decrease is not permanent. After the flare the magnetic field recovers its pre-flare configuration in a short time (i.e., in 30 minutes after the flare). (2) In the photosphere, the line-of-sight velocities show a regular granular up- and down-flow pattern before the flare erupts. During the flare, upflows (blueshifts) dominate the area where the flare is produced. Evaporation rates of ~ $10^{-3}$ and ~ $10^{-4}$ g cm$^{-2}$ s$^{-1}$ have been derived in the deep and high photosphere, respectively, capable of increasing the chromospheric density by a factor of two in about 400 seconds.

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.

preprint2014arXiv

Rayleigh-Taylor instability in partially ionized compressible plasmas: one fluid approach

We study the modification of the classical criterion for the linear onset and growth rate of the Rayleigh-Taylor instability (RTI) in a partially ionized (PI) plasma in the one-fluid description, considering a generalized induction equation. The governing linear equations and appropriate boundary conditions, including gravitational terms, are derived and applied to the case of the RTI in a single interface between two partially ionized plasmas. The boundary conditions lead to an equation for the frequencies in which some of them have positive complex parts, marking the appearance of the RTI. We study the ambipolar term alone first, extending the result to the full induction equation later. We find that the configuration is always unstable because of the presence of a neutral species. In the classical stability regime the growth rate is small, since the collisions prevent the neutral fluid to fully develop the RTI. For parameters in the classical instability regime the growth rate is lowered, but for the considered theoretical values of the collision frequencies and diffusion coefficients for solar prominences the differences with the compressible MHD case are small. We conclude that PI modifies some aspects of the linear RTI instability, since it takes into account that neutrals do not feel the stabilizing effect of the magnetic field. For the set of parameters representative for solar prominences, our model gives the resulting timescale comparable with observed lifetimes of RTI plumes.

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

Dislocations in magnetohydrodynamic waves in a stellar atmosphere

We describe the presence of wavefront dislocations in magnetohydrodynamic waves in stratified stellar atmospheres. Scalar dislocations such as edges and vortices can appear in Alfvén waves, as well as in general magneto-acoustic waves. We detect those dislocations in observations of magnetohydrodynamic waves in sunspots in the solar chromosphere. Through the measured charge of all the dislocations observed, we can give for the first time estimates of the modal contribution in the waves propagating along magnetic fields in solar sunspots.

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

A retrospective of the GREGOR solar telescope in scientific literature

In this review, we look back upon the literature, which had the GREGOR solar telescope project as its subject including science cases, telescope subsystems, and post-focus instruments. The articles date back to the year 2000, when the initial concepts for a new solar telescope on Tenerife were first presented at scientific meetings. This comprehensive bibliography contains literature until the year 2012, i.e., the final stages of commissioning and science verification. Taking stock of the various publications in peer-reviewed journals and conference proceedings also provides the "historical" context for the reference articles in this special issue of Astronomische Nachrichten/Astronomical Notes.

preprint2012arXiv

The GREGOR Fabry-Pérot Interferometer

The GREGOR Fabry-Pérot Interferometer (GFPI) is one of three first-light instruments of the German 1.5-meter GREGOR solar telescope at the Observatorio del Teide, Tenerife, Spain. The GFPI uses two tunable etalons in collimated mounting. Thanks to its large-format, high-cadence CCD detectors with sophisticated computer hard- and software it is capable of scanning spectral lines with a cadence that is sufficient to capture the dynamic evolution of the solar atmosphere. The field-of-view (FOV) of 50" x 38" is well suited for quiet Sun and sunspot observations. However, in the vector spectropolarimetric mode the FOV reduces to 25" x 38". The spectral coverage in the spectroscopic mode extends from 530-860 nm with a theoretical spectral resolution R of about 250,000, whereas in the vector spectropolarimetric mode the wavelength range is at present limited to 580-660 nm. The combination of fast narrow-band imaging and post-factum image restoration has the potential for discovery science concerning the dynamic Sun and its magnetic field at spatial scales down to about 50 km on the solar surface.

preprint2011arXiv

Magneto-acoustic wave energy from numerical simulations of an observed sunspot umbra

We aim at reproducing the height dependence of sunspot wave signatures obtained from spectropolarimetric observations through 3D MHD numerical simulations. A magneto-static sunspot model based on the properties of the observed sunspot is constructed and perturbed at the photosphere introducing the fluctuations measured with the \SiI\ $λ$ 10827 Å line. The results of the simulations are compared with the oscillations observed simultaneously at different heights from the \HeI\ $λ$ 10830 Å line, the \CaIIH\ core and the \FeI\ blends in the wings of the \CaIIH\ line. The simulations show a remarkable agreement with the observations. They reproduce the velocity maps and power spectra at the formation heights of the observed lines, as well as the phase and amplification spectra between several pair of lines. We find that the stronger shocks at the chromosphere are accompanied with a delay between the observed signal and the simulated one at the corresponding height, indicating that shocks shift the formation height of the chromospheric lines to higher layers. Since the simulated wave propagation matches very well the properties of the observed one, we are able to use the numerical calculations to quantify the energy contribution of the magneto-acoustic waves to the chromospheric heating in sunspots. Our findings indicate that the energy supplied by these waves is too low to balance the chromospheric radiative losses. The energy contained at the formation height of the lowermost \SiI\ $λ$ 10827 Å line in the form of slow magneto-acoustic waves is already insufficient to heat the higher layers, and the acoustic energy which reaches the chromosphere is around 3-9 times lower than the required amount of energy. The contribution of the magnetic energy is even lower.

preprint2010arXiv

Magneto-acoustic waves in sunspots from observations and numerical simulations

We study the propagation of waves from the photosphere to the chromosphere of sunspots. From time series of cospatial Ca II H (including its line blends) intensity spectra and polarimetric spectra of Si I 1082.7 nm and He I 1083.0 nm we retrieve the line-of-sight velocity at several heights. The analysis of the phase difference and amplification spectra shows standing waves for frequencies below 4 mHz and propagating waves for higher frequencies, and allows us to infer the temperature and height where the lines are formed. Using these observational data, we have constructed a model of sunspot, and we have introduced the velocity measured with the photospheric Si I 1082.7 nm line as a driver. The numerically propagated wave pattern fits reasonably well with the observed using the lines formed at higher layers, and the simulations reproduce many of the observed features. The observed waves are slow MHD waves propagating longitudinally along field lines.

preprint2010arXiv

Magneto-acoustic waves in sunspots: first results from a new 3D nonlinear magnetohydrodynamic code

Waves observed in the photosphere and chromosphere of sunspots show complex dynamics and spatial patterns. The interpretation of high-resolution sunspot wave observations requires modeling of three-dimensional non-linear wave propagation and mode transformation in the sunspot upper layers in realistic spot model atmospheres. Here we present the first results of such modeling. We have developed a 3D non-linear numerical code specially designed to calculate the response of magnetic structures in equilibrium to an arbitrary perturbation. The code solves the 3D nonlinear MHD equations for perturbations; it is stabilized by hyper-diffusivity terms and is fully parallelized. The robustness of the code is demonstrated by a number of standard tests. We analyze several simulations of a sunspot perturbed by pulses of different periods at subphotospheric level, from short periods, introduced for academic purposes, to longer and realistic periods of three and five minutes. We present a detailed description of the three-dimensional mode transformation in a non-trivial sunspot-like magnetic field configuration, including the conversion between fast and slow magneto-acoustic waves and the Alfvén wave, by calculation of the wave energy fluxes. Our main findings are the following: (1) the conversion from acoustic to the Alfvén mode is only observed if the the driving pulse is located out of the sunspot axis, but this conversion is energetically inefficient; (2) as a consequence of the cut-off effects and refraction of the fast magneto-acoustic mode, the energy of the evanescent waves with periods around 5 minutes remains almost completely below the level beta=1; (3) waves with frequencies above the cut-off propagate field-aligned to the chromosphere and their power becomes dominating over that of evanescent 5-minute oscillations, in agreement with observations.

preprint2010arXiv

Mode transformation and frequency change with height in 3D numerical simulations of magneto-acoustic wave propagation in sunspots

Three-dimensional numerical simulations of magnetoacoustic wave propagation are performed in a sunspot atmosphere with a computational domain covering from the photosphere to the chromosphere. The wave source, with properties resembling the solar spectrum, is located at different distances from the axis of the sunspot for each simulation. These results are compared with the theory of mode transformation and also with observational features. Simulations show that the dominant oscillation frequency in the chromosphere decreases with the radial distance from the sunspot axis. The energy flux of the different wave modes involved, including de Alfvén mode, is evaluated and discussed.

preprint2010arXiv

Multi-layer study of wave propagation in sunspots

We analyze the propagation of waves in sunspots from the photosphere to the chromosphere using time series of co-spatial Ca II H intensity spectra (including its line blends) and polarimetric spectra of Si I 10827 and the He I 10830 multiplet. From the Doppler shifts of these lines we retrieve the variation of the velocity along the line-of-sight at several heights. Phase spectra are used to obtain the relation between the oscillatory signals. Our analysis reveals standing waves at frequencies lower than 4 mHz and a continuous propagation of waves at higher frequencies, which steepen into shocks in the chromosphere when approaching the formation height of the Ca II H core. The observed non-linearities are weaker in Ca II H than in He I lines. Our analysis suggests that the Ca II H core forms at a lower height than the He I 10830 line: a time delay of about 20 s is measured between the Doppler signal detected at both wavelengths. We fit a model of linear slow magnetoacoustic wave propagation in a stratified atmosphere with radiative losses according to Newton's cooling law to the phase spectra and derive the difference in the formation height of the spectral lines. We show that the linear model describes well the wave propagation up to the formation height of Ca II H, where non-linearities start to become very important.

preprint2010arXiv

Numerical Simulation of Excitation and Propagation of Helioseismic MHD Waves in Magnetostatic Models of Sunspots

We present comparison of numerical simulations of propagation of MHD waves,excited by subphotospheric perturbations, in two different ("deep" and "shallow") magnetostatic models of the sunspots. The "deep" sunspot model distorts both the shape of the wavefront and its amplitude stronger than the "shallow" model. For both sunspot models, the surface gravity waves (f-mode) are affected by the sunspots stronger than the acoustic p-modes. The wave amplitude inside the sunspot depends on the photospheric strength of the magnetic field and the distance of the source from the sunspot axis. For the source located at 9 Mm from the center of the sunspot, the wave amplitude increases when the wavefront passes through the central part of the sunspot. For the source distance of 12 Mm, the wave amplitude inside the sunspot is always smaller than outside. For the same source distance from the sunspot center but for the models with different strength of the magnetic field, the wave amplitude inside the sunspot increases with the strength of the magnetic field. The simulations show that unlike the case of the uniform inclined background magnetic field, the p- and f-mode waves are not spatially separated inside the sunspot where the magnetic field is strongly non-uniform. These properties have to be taken into account for interpretation of observations of MHD waves traveling through sunspot regions.

preprint2010arXiv

The Imaging Magnetograph eXperiment (IMaX) for the Sunrise balloon-borne solar observatory

The Imaging Magnetograph eXperiment (IMaX) is a spectropolarimeter built by four institutions in Spain that flew on board the Sunrise balloon-borne telesocope in June 2009 for almost six days over the Arctic Circle. As a polarimeter IMaX uses fast polarization modulation (based on the use of two liquid crystal retarders), real-time image accumulation, and dual beam polarimetry to reach polarization sensitivities of 0.1%. As a spectrograph, the instrument uses a LiNbO3 etalon in double pass and a narrow band pre-filter to achieve a spectral resolution of 85 mAA. IMaX uses the high Zeeman sensitive line of Fe I at 5250.2 AA and observes all four Stokes parameters at various points inside the spectral line. This allows vector magnetograms, Dopplergrams, and intensity frames to be produced that, after reconstruction, reach spatial resolutions in the 0.15-0.18 arcsec range over a 50x50 arcsec FOV. Time cadences vary between ten and 33 seconds, although the shortest one only includes longitudinal polarimetry. The spectral line is sampled in various ways depending on the applied observing mode, from just two points inside the line to 11 of them. All observing modes include one extra wavelength point in the nearby continuum. Gauss equivalent sensitivities are four Gauss for longitudinal fields and 80 Gauss for transverse fields per wavelength sample. The LOS velocities are estimated with statistical errors of the order of 5-40 m/s. The design, calibration and integration phases of the instrument, together with the implemented data reduction scheme are described in some detail.

preprint2009arXiv

Sunspot seismic halos generated by fast MHD wave refraction

We suggest an explanation for the high-frequency power excess surrounding active regions known as seismic halos. The idea is based on numerical simulations of magneto-acoustic waves propagation in sunspots. We propose that such an excess can be caused by the additional energy injected by fast mode waves refracted in the higher atmosphere due to the rapid increase of the Alfven speed. Our model qualitatively explains the magnitude of the halo and allows to make some predictions of its behavior that can be checked in future observations.

preprint2009arXiv

The energy of waves in the photosphere and lower chromosphere: 1. Velocity statistics

Acoustic waves are one of the primary suspects besides magnetic fields for the chromospheric heating process to temperatures above radiative equilibrium (RE). We derived the mechanical wave energy as seen in line-core velocities to obtain a measure of mechanical energy flux with height for a comparison with the energy requirements in a semi-empirical atmosphere model. We analyzed a 1-hour time series and a large-area map of Ca II H spectra on the traces of propagating waves. We analyzed the velocity statistics of several spectral lines in the wing of Ca II H, and the line-core velocity of Ca II H. We converted the velocity amplitudes into volume and mass energy densities. For comparison, we used the increase of internal energy necessary to lift a RE atmosphere to the HSRA temperature stratification. We find that the velocity amplitude grows in agreement with linear wave theory and thus slower with height than predicted from energy conservation. The mechanical energy of the waves above around z~500 km is insufficient to maintain the chromospheric temperature rise in the semi-empirical HSRA model. The intensity variations of the Ca line core (z~1000 km) can be traced back to the velocity variations of the lowermost forming spectral line considered (z~ 250 km). The chromospheric intensity, and hence, (radiation) temperature variations are seen to be induced by passing waves originating in the photosphere.

preprint2008arXiv

Theoretical modeling of propagation of magneto-acoustic waves in magnetic regions below sunspots

We use 2D numerical simulations and eikonal approximation, to study properties of MHD waves traveling below the solar surface through the magnetic structure of sunspots. We consider a series of magnetostatic models of sunspots of different magnetic field strengths, from 10 Mm below the photosphere to the low chromosphere. The purpose of these studies is to quantify the effect of the magnetic field on local helioseismology measurements by modeling waves excited by sub-photospheric sources. Time-distance propagation diagrams and wave travel times are calculated for models of various field strength and compared to the non-magnetic case. The results clearly indicate that the observed time-distance helioseismology signals in sunspot regions correspond to fast MHD waves. The slow MHD waves form a distinctly different pattern in the time-distance diagram, which has not been detected in observations. The numerical results are in good agreement with the solution in the short-wavelength (eikonal) approximation, providing its validation. The frequency dependence of the travel times is in a good qualitative agreement with observations.

preprint2006arXiv

A near-IR line of Mn I as a diagnostic tool of the average magnetic energy in the solar photosphere

We report on spectropolarimetric observations of a near-IR line of Mn I located at 15262.702 A whose intensity and polarization profiles are very sensitive to the presence of hyperfine structure. A theoretical investigation of the magnetic sensitivity of this line to the magnetic field uncovers several interesting properties. The most important one is that the presence of strong Paschen-Back perturbations due to the hyperfine structure produces an intensity line profile whose shape changes according to the absolute value of the magnetic field strength. A line ratio technique is developed from the intrinsic variations of the line profile. This line ratio technique is applied to spectropolarimetric observations of the quiet solar photosphere in order to explore the probability distribution function of the magnetic field strength. Particular attention is given to the quietest area of the observed field of view, which was encircled by an enhanced network region. A detailed theoretical investigation shows that the inferred distribution yields information on the average magnetic field strength and the spatial scale at which the magnetic field is organized. A first estimation gives ~250 G for the mean field strength and a tentative value of ~0.45" for the spatial scale at which the observed magnetic field is horizontally organized.