BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

J. Leenaarts

J. Leenaarts contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.SRastro-ph.IM

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
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

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Chromospheric Extension of the MURaM Code

Detailed numerical models of chromosphere and corona are required to understand the heating of the solar atmosphere. An accurate treatment of the solar chromosphere is complicated by the effects arising from Non Local Thermodynamic Equilibrium (NLTE) radiative transfer. A small number of strong, highly scattering lines dominate the cooling and heating in the chromosphere. Additionally, the recombination times of ionised hydrogen are longer than the dynamical timescales, requiring a non-equilibrium (NE) treatment of hydrogen ionisation. The MURaM code is extended to include the physical process required for accurate simulation of the solar chromosphere, as implemented in the Bifrost code. This includes a time-dependent treatment of hydrogen ionisation, a scattering multi-group radiation transfer scheme and approximations for NLTE radiative cooling. The inclusion of NE and NLTE physics has a large impact on the structure of the chromosphere; the NE treatment of hydrogen ionisation leads to a higher ionisation fraction and enhanced populations in the first excited state throughout cold inter-shock regions of the chromosphere. Additionally this prevents hydrogen ioniation from buffering energy fluctuations, leading to hotter shocks and cooler inter-shock regions. The hydrogen populations in the ground and first excited state are enhanced by $10^2-10^3$ in the upper chromosphere and up to $10^9$ near the transition region. Including the necessary NLTE physics leads to significant differences in chromospheric structure and dynamics. The thermodynamics and hydrogen populations calculated using the extended version of the MURaM code are consistent with previous non-equilibrium simulations. The electron number and temperature calculated using the non-equilibrium treatment of the chromosphere are required to accurately synthesise chromospheric spectral lines.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Quiet Sun Center to Limb Variation of the Linear Polarization Observed by CLASP2 Across the Mg II h & k Lines

The CLASP2 (Chromospheric LAyer SpectroPolarimeter 2) sounding rocket mission was launched on 2019 April 11. CLASP2 measured the four Stokes parameters of the Mg II h & k spectral region around 2800 Angstroms along a 200 arcsecond slit at three locations on the solar disk, achieving the first spatially and spectrally resolved observations of the solar polarization in this near ultraviolet region. The focus of the work presented here is the center-to-limb variation of the linear polarization across these resonance lines, which is produced by the scattering of anisotropic radiation in the solar atmosphere. The linear polarization signals of the Mg II h & k lines are sensitive to the magnetic field from the low to the upper chromosphere through the Hanle and magneto-optical effects. We compare the observations to theoretical predictions from radiative transfer calculations in unmagnetized semi-empirical models, arguing that magnetic fields and horizontal inhomogeneities are needed to explain the observed polarization signals and spatial variations. This comparison is an important step in both validating and refining our understanding of the physical origin of these polarization signatures, and also in paving the way toward future space telescopes for probing the magnetic fields of the solar upper atmosphere via ultraviolet spectropolarimetry.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Subarcsecond imaging of a solar active region filament with ALMA and IRIS

Quiescent filaments appear as absorption features on the solar disk when observed in chromospheric lines and at continuum wavelengths in the millimeter (mm) range. Active region (AR) filaments are their small-scale, low-altitude analogues, but they could not be resolved in previous mm observations. This spectral diagnostic can provide insight into the details of the formation and physical properties of their fine threads, which are still not fully understood. Here, we shed light on the thermal structure of an AR filament using high-resolution brightness temperature ($T_{\rm b}$) maps taken with ALMA Band 6 complemented by simultaneous IRIS near-UV spectra, Hinode/SOT photospheric magnetograms, and SDO/AIA extreme-UV images. Some of the dark threads visible in the AIA 304 Å passband and in the core of Mg II resonance lines have dark ($T_{\rm b}<5000$K) counterparts in the 1.25 mm maps, but their visibility significantly varies across the filament spine and in time. These opacity changes are possibly related to variations in temperature and electron density in filament fine structures. The coolest $T_{\rm b}$ values ($<$5000 K) coincide with regions of low integrated intensity in the Mg II h and k lines. ALMA Band 3 maps taken after the Band 6 ones do not clearly show the filament structure, contrary to the expectation that the contrast should increase at longer wavelengths based on previous observations of quiescent filaments. The ALMA maps are not consistent with isothermal conditions, but the temporal evolution of the filament may partly account for this.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The formation and heating of chromospheric fibrils in a radiation-MHD simulation

Aims: We examine the movements of mass elements within dense fibrils using passive tracer particles (corks) in order to understand fibril creation and destruction processes. Methods: Simulated fibrils were selected at times when they were visible in an H$α$ image proxy. The corks were selected within fibril H$α$ formation regions. From this set, a cork was selected, and the field line passing through it was constructed. Other fibrilar corks close to this fieldline were also selected. Pathlines were constructed, revealing the locations of the mass elements forward and backward in time. The forces acting on these mass elements were analysed. Results: The main process of fibrilar loading in the simulation is different to the mass loading scenario in which waves steepen into shocks and push material upwards along the fieldlines from near their footpoints. Twisted low lying fieldlines were destabilised and then they untwisted, lifting the material trapped above their apexes via the Lorentz force. Subsequently, the majority of the mass drained down the fieldlines towards one or both footpoints under gravity. Material with large horizontal velocities could also elevated in rising fieldlines, creating somewhat parabolic motions, but material was not generally moving upward along a stationary magnetic fieldline during loading. Conclusions: The processes observed in the simulation are plausible additional scenarios. Criteria for observing such events are described. It is desirable that our simulations can also form more densely-packed fibrils from material fed from the base of field footpoints. Experimental parameters required to achieve this are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Non-LTE inversions of a confined X2.2 flare: I. Vector magnetic field in the photosphere and chromosphere

Obtaining the magnetic field vector accurately in the solar atmosphere is essential for studying changes in field topology during flares and to reliably model space weather. We tackle this problem by applying various inversion methods to a confined X2.2 flare in NOAA AR 12673 on September 6, 2017, comparing the photospheric and chromospheric magnetic field vector with those from two numerical models of this event. We obtain the photospheric field from Milne-Eddington (ME) and (non-)local thermal equilibrium (non-LTE) inversions of Hinode SOT/SP Fe I 6301.5Å and 6302.5Å. The chromospheric field is obtained from a spatially-regularised weak field approximation (WFA) and non-LTE inversions of Ca II 8542Å observed with CRISP at the Swedish 1-m Solar Telescope. The LTE- and non-LTE-inferred photospheric field components are strongly correlated throughout the atmosphere, with stronger field and higher temperatures in the non-LTE inversions. For the chromospheric field, the non-LTE inversions correlate well with the spatially-regularised WFA. We find strong-field patches of over 4.5 kG in the photosphere, co-located with similar concentrations exceeding 3 kG in the chromosphere. The obtained field strengths are up to 2-3 times higher than in the numerical models, with more concentrated and structured photosphere-to-chromosphere shear close to the polarity inversion line. The LTE and non-LTE Fe I inversions yield essentially the same photospheric field, while ME inversions fail to reproduce the field vector orientation where Fe I is in emission. Our inversions confirm the locations of flux rope footpoints that are predicted by numerical models. However, pre-processing and lower spatial resolution lead to weaker and smoother field in the models than what the data indicate. This emphasises the need for higher spatial resolution in the models to better constrain pre-eruptive flux ropes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Observational constraints on the origin of the elements II. 3D non-LTE formation of Ba ii lines in the solar atmosphere

Context. The pursuit of more realistic spectroscopic modelling and consistent abundances has led us to begin a new series of papers designed to improve current solar and stellar abundances of various atomic species. To achieve this, we have began updating the three-dimensional (3D) non-local thermodynamic equilibrium (non-LTE) radiative transfer code, Multi3D, and the equivalent one-dimensional (1D) non-LTE radiative transfer code, MULTI. Aims. We examine our improvements to these codes by redetermining the solar barium abundance. Barium was chosen for this test as it is an important diagnostic element of the s-process in the context of galactic chemical evolution. New Ba II + H collisional data for excitation and charge exchange reactions computed from first principles had recently become available and were included in the model atom. The atom also includes the effects of isotopic line shifts and hyperfine splitting. Method. A grid of 1D LTE barium lines were constructed with MULTI and fit to the four Ba II lines available to us in the optical region of the solar spectrum. Abundance corrections were then determined in 1D non-LTE, 3D LTE, and 3D non-LTE. A new 3D non-LTE solar barium abundance was computed from these corrections. Results. We present for the first time the full 3D non-LTE barium abundance of $A({\rm Ba})=2.27\pm0.02\pm0.01$, which was derived from four individual fully consistent barium lines. Errors here represent the systematic and random errors, respectively.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access