Source author record

J. de la Cruz Rodríguez

J. de la Cruz Rodríguez 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

17works
3topics
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

17 published item(s)

preprint2025arXiv

The making of robust and highly performing imaging spectropolarimeters for large solar telescopes

We discuss the requirements, concepts, simulations, implementation, and calibration of two dual Fabry-Perot based imaging spectropolarimeters, CRISP and CHROMIS, at the Swedish 1-meter Solar Telescope, and CRISP2 that is under construction. These instruments use a combination of a high-resolution and a low-resolution etalon together with an order-sorting prefilter to define the bandpass. The overall design is made robust and stable by tailoring the low-resolution etalon reflectivity to accommodate expected cavity errors from both etalons, and by using a compact optical design that eliminates the need for folding mirrors. By using a telecentric design based on lenses rather than mirrors, image degradation by the FPI system is negligible, as shown in a previous publication, and the throughput of the system is maximised. Initial alignment, and maintaining that alignment over time, is greatly simplified. The telecentric design allows full calibration and/or modelling of essential system parameters to be carried out without interfering with the optical setup. We also discuss briefly the polarimeters developed for CRISP and CHROMIS. The high performance of CRISP and CHROMIS has been demonstrated in an earlier publication through measurements of the granulation contrast and comparisons with similar measurements simultaneously made through broadband continuum filters. Here, we focus on the aspects of the design that are central to enabling high performance and robustness, but also discuss the calibration and processing of the data, and use a few examples of processed data to demonstrate the achievable image and data quality. We put forward a proposal for a similar conceptual design for the European Solar Telescope and conclude by discussing potential problems of the proposed approach to designs of this type. Some aspects of these FPI systems may be of interest also outside the solar community.

preprint2022arXiv

Heating of the solar chromosphere through current dissipation

The solar chromosphere is heated to temperatures higher than predicted by radiative equilibrium. This excess heating is greater in active regions where the magnetic field is stronger. We aim to investigate the magnetic topology associated with an area of enhanced millimeter (mm) brightness temperatures in a solar active region mapped by the Atacama Large Millimeter/submillimeter Array (ALMA) using spectropolarimetric co-observations with the 1-m Swedish Solar Telescope (SST). We used Milne-Eddington inversions, nonlocal thermodynamic equilibrium (non-LTE) inversions, and a magnetohydrostatic extrapolation to obtain constraints on the three-dimensional stratification of temperature, magnetic field, and radiative energy losses. We compared the observations to a snapshot of a magnetohydrodynamics simulation and investigate the formation of the thermal continuum at 3 mm using contribution functions. We find enhanced heating rates in the upper chromosphere of up to $\sim 5\rm\,kW\,m^{-2}$, where small-scale emerging loops interact with the overlying magnetic canopy leading to current sheets as shown by the magnetic field extrapolation. Our estimates are about a factor of two higher than canonical values, but they are limited by the ALMA spatial resolution ($\sim 1.2^{\prime\prime}$). Band 3 brightness temperatures reach about $\sim10^{4}\,$K in the region, and the transverse magnetic field strength inferred from the non-LTE inversions is on the order of $\sim 500\,$G in the chromosphere. We are able to quantitatively reproduce many of the observed features, including the integrated radiative losses in our numerical simulation. We conclude that the heating is caused by dissipation in current sheets. However, the simulation shows a complex stratification in the flux emergence region where distinct layers may contribute significantly to the emission in the mm continuum.

preprint2022arXiv

Radiative losses in the chromosphere during a C-class flare

Solar flares release an enormous amount of energy into the corona. A substantial fraction of this energy is transported to the lower atmosphere, which results in chromospheric heating. The mechanisms that transport energy to the lower solar atmosphere during a flare are still not fully understood. We aim to estimate the temporal evolution of the radiative losses in the chromosphere at the footpoints of a C-class flare, in order to set observational constraints on the electron beam parameters of a RADYN flare simulation. We estimated the radiative losses from hydrogen, and singly ionized Ca and Mg using semi-empirical model atmospheres. To estimate the integrated radiative losses in the chromosphere the net cooling rates were integrated between the temperature minimum and the height where the temperature reaches 10 kK. The stratification of the net cooling rate suggests that the Ca IR triplet lines are responsible for most of the radiative losses in the flaring atmosphere. During the flare peak time, the contribution from Ca II H & K and Mg II h & k lines are strong and comparable to the Ca IR triplet ($\sim$32 kW m$^{-2}$). Since our flare is a relatively weak event the chromosphere is not heated above 11 kK, which in turn yields a subdued Lyα contribution ($\sim$7 kW m$^{-2}$). The temporal evolution of total integrated radiative losses exhibits sharply-rising losses (0.4 kW m$^{-2}$ s$^{-1}$) and a relatively slow decay (0.23 kW~m$^{-2}$ s$^{-1}$). The maximum value of total radiative losses is reached around the flare peak time, and can go up to 175 kW m$^{-2}$ for a single pixel located at footpoint. After a small parameter study, we find the best model-data consistency in terms of the amplitude of radiative losses and the overall atmospheric structure with a RADYN flare simulation in the injected energy flux of $5\times10^{10}$ erg s$^{-1}$ cm$^{-2}$.

preprint2022arXiv

Spatio-temporal analysis of chromospheric heating in a plage region

Our knowledge of the heating mechanisms that are at work in the chromosphere of plage regions remains highly unconstrained from observational studies. The purpose of our study is to estimate the chromospheric heating terms from a plage dataset, characterize their spatio-temporal distribution and set constraints on the heating processes that are at work. We make use of NLTE inversions to infer a model of the photosphere and chromosphere of a plage dataset acquired with the Swedish 1-m Solar Telescope. We use this model atmosphere to calculate the chromospheric radiative losses from H i, Ca ii and Mg ii atoms. We approximate the chromospheric heating terms by the net radiative losses predicted by the inverted model. In order to make the analysis of time-series over a large field-of-view computationally tractable, we make use of a neural network. In the lower chromosphere, the contribution from the Ca ii lines is dominant and located in the surroundings of the photospheric footpoints. In the upper chromosphere, the H i contribution is dominant. Radiative losses in the upper chromosphere form an homogeneous patch that covers the plage region. The net radiative losses can be split in a periodic component with an average amplitude of ampQ = 7.6 kW m^{-2} and a static (or very slowly evolving) component with a mean value of -26.1 kW m^{-2}. Our interpretation is that in the lower chromosphere, the radiative losses are tracing the sharp lower edge of the hot magnetic canopy, where the electric current is expected to be large. In the upper chromosphere, both the magnetic field and the distribution of net radiative losses are room-filling, whereas the amplitude of the periodic component is largest. Our results suggest that acoustic wave heating may be responsible for one third of the energy deposition in the upper chromosphere, whereas other heating mechanisms must be responsible for the rest.

preprint2021arXiv

Active region chromospheric magnetic fields

Context. A proper estimate of the chromospheric magnetic fields is believed to improve modelling of both active region and coronal mass ejection evolution. Aims. We investigate the similarity between the chromospheric magnetic field inferred from observations and the field obtained from a magnetohydrostatic (MHS) extrapolation. Methods. Based Fe i 6173 Å and Ca ii 8542 Å observations of NOAA active region 12723, we employed the spatially-regularised weak-field approximation (WFA) to derive the vector magnetic field in the chromosphere from Ca ii, as well as non-LTE inversions of Fe i and Ca ii to infer a model atmosphere for selected regions. Milne-Eddington inversions of Fe i serve as photospheric boundary for the MHS model that delivers the three-dimensional field, gas pressure and density. Results. For the line-of-sight component, the MHS chromospheric field generally agrees with the non-LTE inversions and WFA, but tends to be weaker than those when larger in magnitude than 300 G. The observationally inferred transverse component is stronger, especially in magnetically weaker regions, yet the qualitative distribution with height is similar to the MHS results. For either field component the MHS chromospheric field lacks the fine structure derived from the inversions. Furthermore, the MHS model does not recover the magnetic imprint from a set of high fibrils connecting the main polarities. Conclusions. The MHS extrapolation and WFA provide a qualitatively similar chromospheric field, where the azimuth of the former is better aligned with Ca ii 8542 Å fibrils than that of the WFA, especially outside strong-field concentrations. The amount of structure as well as the transverse field strengths are underestimated by the MHS extrapolation. This underscores the importance of considering a chromospheric magnetic field constraint in data-driven modelling of active regions.

preprint2021arXiv

Bayesian Stokes inversion with Normalizing flows

Stokes inversion techniques are very powerful methods for obtaining information on the thermodynamic and magnetic properties of solar and stellar atmospheres. In recent years, very sophisticated inversion codes have been developed that are now routinely applied to spectro-polarimetric observations. Most of these inversion codes are designed for finding an optimum solution to the nonlinear inverse problem. However, to obtain the location of potentially multimodal cases (ambiguities), the degeneracies, and the uncertainties of each parameter inferred from the inversions, algorithms such as Markov chain Monte Carlo (MCMC), require to evaluate the likelihood of the model thousand of times and are computationally costly. Variational methods are a quick alternative to Monte Carlo methods by approximating the posterior distribution by a parametrized distribution. In this study, we introduce a highly flexible variational inference method to perform fast Bayesian inference, known as normalizing flows. Normalizing flows are a set of invertible, differentiable, and parametric transformations that convert a simple distribution into an approximation of any other complex distribution. If the transformations are conditioned on observations, the normalizing flows can be trained to return Bayesian posterior probability estimates for any observation. We illustrate the ability of the method using a simple Milne-Eddington model and a complex non-LTE inversion. However, the method is extremely general and other more complex forward models can be applied. The training procedure need only be performed once for a given prior parameter space and the resulting network can then generate samples describing the posterior distribution several orders of magnitude faster than existing techniques.

preprint2021arXiv

Design and Performance Analysis of a Highly Efficient Polychromatic Full-Stokes Polarization Modulator for the CRISP Imaging Spectrometer

We present the design and performance of a polychromatic polarization modulator for the CRisp Imaging SpectroPolarimeter (CRISP) Fabry-Perot tunable narrow-band imaging spectropolarimer at the Swedish 1-m Solar Telescope (SST). We discuss the design process in depth, compare two possible modulator designs through a tolerance analysis, and investigate thermal sensitivity of the selected design. The trade-offs and procedures described in this paper are generally applicable in the development of broadband polarization modulators. The modulator was built and has been operational since 2015. Its measured performance is close to optimal between 500 and 900~nm, and differences between the design and as-built modulator are largely understood. We show some example data, and briefly review scientific work that used data from SST/CRISP and this modulator.

preprint2021arXiv

Downflowing umbral flashes as an evidence of standing waves in sunspot umbrae

Umbral flashes are sudden brightenings commonly visible in the core of chromospheric lines. Theoretical and numerical modeling suggest that they are produced by the propagation of shock waves. According to these models and early observations, umbral flashes are associated with upflows. However, recent studies have reported umbral flashes in downflowing atmospheres. We aim to understand the origin of downflowing umbral flashes. We explore how the existence of standing waves in the umbral chromosphere impacts the generation of flashed profiles. We performed numerical simulations of wave propagation in a sunspot umbra with the code MANCHA. The Stokes profiles of the Ca II 8542 Å line were synthesized with NICOLE. For freely-propagating waves, the chromospheric temperature enhancements of the oscillations are in phase with velocity upflows. In this case, the intensity core of the Ca II 8542 Å atmosphere is heated during the upflowing stage of the oscillation. If we consider a different scenario with a resonant cavity, the wave reflections at the sharp temperature gradient of the transition region lead to standing oscillations. In this situation, temperature fluctuations are shifted backward and temperature enhancements partially coincide with the downflowing stage of the oscillation. In umbral flashes produced by standing oscillations, the reversal of the emission feature is produced when the oscillation is downflowing. The chromospheric temperature keeps increasing while the atmosphere is changing from a downflow to an upflow. During the appearance of flashed Ca II 8542 Å cores, the atmosphere is upflowing most of the time, and only 38\% of the flashed profiles are associated with downflows. We find a scenario that remarkably explains the recent empirical findings of downflowing umbral flashes as a natural consequence of the presence of standing oscillations above sunspot umbrae.

preprint2020arXiv

Inference of the chromospheric magnetic field configuration of solar plage using the Ca II 8542 Å line

It has so far proven impossible to reproduce all aspects of the solar plage chromosphere in quasi-realistic numerical models. The magnetic field configuration in the lower atmosphere is one of the few free parameters in such simulations. The literature only offers proxy-based estimates of the field strength, as it is difficult to obtain observational constraints in this region. Sufficiently sensitive spectro-polarimetric measurements require a high signal-to-noise ratio, spectral resolution, and cadence, which are at the limit of current capabilities. We use critically sampled spectro-polarimetric observations of the \cair line obtained with the CRISP instrument of the Swedish 1-m Solar Telescope to study the strength and inclination of the chromospheric magnetic field of a plage region. This will provide direct physics-based estimates of these values, which could aid modelers to put constraints on plage models. We increased the signal-to-noise ratio of the data by applying several methods including deep learning and PCA. We estimated the noise level to be $1\cdot10^{-3} I_c$. We then used STiC, a non-local thermodynamic equilibrium (NLTE) inversion code to infer the atmospheric structure and magnetic field pixel by pixel. We are able to infer the magnetic field strength and inclination for a plage region and for fibrils in the surrounding canopy. In the plage we report an absolute field strength of $|B| =440 \pm 90$ G, with an inclination of $10^\circ \pm 16^\circ$ with respect to the local vertical. This value for $|B|$ is roughly double of what was reported previously, while the inclination matches previous studies done in the photosphere. In the fibrillar region we found $|B| = 300 \pm 50$ G, with an inclination of $50^\circ \pm 13^\circ$.

preprint2019arXiv

The multi-thermal chromosphere: inversions of ALMA and IRIS data

Numerical simulations of the solar chromosphere predict a diverse thermal structure with both hot and cool regions. Observations of plage regions in particular feature broader and brighter chromospheric lines, which suggest that they are formed in hotter and denser conditions than in the quiet-Sun, but also implies a non-thermal component whose source is unclear. We revisit the problem of the stratification of temperature and microturbulence in plage now adding millimeter continuum observations provided by ALMA to inversions of near-ultraviolet IRIS spectra as a powerful new diagnostic to disentangle the two parameters. We fit cool chromospheric holes and track the fast evolution of compact mm brightenings in the plage region. We use the STiC non-LTE inversion code to simultaneously fit real ultraviolet and millimeter spectra in order to infer the thermodynamic parameters of the plasma. We confirm the anticipated constraining potential of ALMA in non-LTE inversions of the solar chromosphere. We find significant differences between the inversion results of IRIS data alone compared to the results of a combination with the mm data: the IRIS+ALMA inversions have increased contrast and temperature range, and tend to favor lower values of microturbulence in the chromosphere of plage. The average brightness temperature of the plage region at 1.25 mm is 8500 K, but the ALMA maps also show much cooler ($\sim3000$ K) and hotter ($\sim11\,000$ K) evolving features partially seen in other diagnostics. To explain the former, the inversions require the existence of localized, low temperature regions in the chromosphere where molecules such as CO could form. The hot features could sustain such high temperatures due to non-equilibrium hydrogen ionization effects in a shocked chromosphere - a scenario that is supported by low-frequency shock wave patterns found in the MgII lines probed by IRIS.

preprint2016arXiv

Spectropolarimetric capabilities of Ca II 8542 A line

The next generation of space and ground-based solar missions aim to study the magnetic properties of the solar chromosphere using the infrared Ca II lines and the He I 10830 Å line. The former seem to be the best candidates to study the stratification of magnetic fields in the solar chromosphere and their relation to the other thermodynamical properties underlying the chromospheric plasma. The purpose of this work is to provide a detailed analysis of the diagnostic capabilities of the Ca II 8542 Å line, anticipating forthcoming observational facilities. We study the sensitivity of the Ca II 8542 Å line to perturbations applied to the physical parameters of reference semi-empirical 1D model atmospheres using response functions and we make use of 3D MHD simulations to examine the expected polarization signals for moderate magnetic field strengths. Our results indicate that the Ca II 8542 Å line is mostly sensitive to the layers enclosed between $\log$ $τ=[0,-5.5]$, under the physical conditions that are present in our model atmospheres. In addition, the simulated magnetic flux tube generates strong longitudinal signals in its centre and moderate transversal signals, due to the vertical expansion of magnetic field lines, in its edge. Thus, observing the Ca II 8542 Å line we will be able to infer the 3D geometry of moderate magnetic field regions.

preprint2014arXiv

CRISPRED: A data pipeline for the CRISP imaging spectropolarimeter

The production of science-ready data from major solar telescopes requires expertise beyond that of the typical observer. This is a consequence of the increasing complexity of instruments and observing sequences, which require calibrations and corrections for instrumental and seeing effects that are not only difficult to measure, but are also coupled in ways that require careful analysis in the design of the correction procedures. Modern space-based telescopes have data-processing pipelines capable of routinely producing well-characterized data products. High-resolution imaging spectropolarimeters at ground-based telescopes need similar data pipelines. We present new methods for flat-fielding spectropolarimetric data acquired with telecentric Fabry-Perot instruments and a new approach for accurate camera co-alignment for image restoration. We document a procedure that forms the basis of current state of the art processing of data from the CRISP imaging spectropolarimeter at the Swedish 1-m Solar Telescope (SST). By collecting, implementing, and testing a suite of computer programs, we have defined a data reduction pipeline for this instrument. This pipeline, CRISPRED, streamlines the process of making science-ready data. It is implemented and operated in IDL, with time-consuming steps delegated to C. CRISPRED will be the basis also for the data pipeline of the forthcoming CHROMIS instrument.

preprint2013arXiv

Heating of the magnetic chromosphere: observational constraints from Ca II 8542 spectra

The heating of the Sun's chromosphere remains poorly understood. While progress has been made on understanding what drives the quiet Sun internetwork chromosphere, chromospheric heating in strong magnetic field regions continues to present a difficult challenge, mostly because of a lack of observational constraints. We use high-resolution spectropolarimetric data from the Swedish 1-m Solar Telescope to identify the location and spatio-temporal properties of heating in the magnetic chromosphere. In particular, we report the existence of raised-core spectral line profiles in the Ca II 8542 line. These profiles are characterized by the absence of an absorption line core, showing a quasi-flat profile between +/- 0.5 Å, and are abundant close to magnetic bright-points and plage. Comparison with 3D MHD simulations indicates that such profiles occur when the line-of-sight goes through an "elevated temperature canopy" associated with the expansion with height of the magnetic field of flux concentrations. This temperature canopy in the simulations is caused by ohmic dissipation where there are strong magnetic field gradients. The raised-core profiles are thus indicators of locations of increased chromospheric heating. We characterize the location and temporal and spatial properties of such profiles in our observations, thus providing much stricter constraints on theoretical models of chromospheric heating mechanisms than before.

preprint2013arXiv

Opposite polarity field with convective downflow and its relation to magnetic spines in a sunspot penumbra

We discuss NICOLE inversions of Fe I 630.15 nm and 630.25 nm Stokes spectra from a sunspot penumbra recorded with the CRISP imaging spectropolarimeter on the Swedish 1-m Solar Telescope at a spatial resolution close to 0.15". We report on narrow radially extended lanes of opposite polarity field, located at the boundaries between areas of relatively horizontal magnetic field (the intra-spines) and much more vertical field (the spines). These lanes harbor convective downflows of about 1 km/s. The locations of these downflows close to the spines agree with predictions from the convective gap model (the "gappy penumbra") proposed six years ago, and more recent 3D MHD simulations. We also confirm the existence of strong convective flows throughout the entire penumbra, showing the expected correlation between temperature and vertical velocity, and having vertical RMS velocities of about 1.2 km/s.

preprint2013arXiv

Physical Properties of a Sunspot Chromosphere with Umbral Flashes

We present new high-resolution spectro-polarimetric Ca II 8542 observations of umbral flashes in sunspots. At nearly 0.18", and spanning about one hour of continuous observation, this is the most detailed dataset published thus far. Our study involves both LTE and non-LTE inversions (but includes also a weak field analysis as a sanity check) to quantify temperatures, mass flows and the full magnetic field vector geometry. We confirm earlier reports that UFs have very fine structure with hot and cool material intermixed at sub-arcsecond scales. The shock front is roughly 1000 K hotter than the surrounding material. We do not observe significant fluctuations of the field in the umbra. In the penumbra, however, the passage of the running penumbral waves alter the magnetic field strength by some 200 G (peak-to-peak amplitude) but it does not change the field orientation (at least not significantly within our sensitivity of a few degrees). We find a trend of decreasing high-frequency modulation power for more inclined magnetic fields for the line-of-sight velocity and magnetic field strength. In the outer penumbra we find an absence of high frequency power while there is increasingly more power at high frequencies towards the umbra.

preprint2012arXiv

Chromospheric magnetic fields. Observations, simulations and their interpretation

The magnetic field of the quiet-Sun chromosphere remains a mystery for solar physicists. The reduced number of chromospheric lines are intrinsically hard to model and only a few of them are magnetically sensitive. In this work, we use a 3D numerical simulation of the outer layers of the solar atmosphere, to asses the reliability of non-LTE inversions, in this case applied to the Ca II 8542 Å line. We show that NLTE inversions provide realistic estimates of physical quantities from synthetic observations.

preprint2011arXiv

Are solar chromospheric fibrils tracing the magnetic field?

Fibrils are thin elongated features visible in the solar chromosphere in and around magnetized regions. Because of their visual appearance, they have been traditionally considered a tracer of the magnetic field lines. For the first time, we challenge that notion, by comparing their orientation to that of the magnetic field, obtained via high-resolution spectropolarimetric observations of Ca ii lines. The short answer to the question posed in the title is that mostly yes, but not always.