Researcher profile

J. de la Cruz Rodríguez

J. de la Cruz Rodríguez contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
2topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 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.