B. Ruiz Cobo

An Open Source, Massively Parallel Code for Non-LTE Synthesis and Inversion of Spectral Lines and Zeeman-induced Stokes Profiles

With the advent of a new generation of solar telescopes and instrumentation, interpreting chromospheric observations (in particular, spectropolarimetry) requires new, suitable diagnostic tools. This paper describes a new code, NICOLE, that has been designed for Stokes non-LTE radiative transfer, for synthesis and inversion of spectral lines and Zeeman-induced polarization profiles, spanning a wide range of atmospheric heights from the photosphere to the chromosphere. The code features a number of unique features and capabilities and has been built from scratch with a powerful parallelization scheme that makes it suitable for application on massive datasets using large supercomputers. The source code is written entirely in Fortran 90/2003 and complies strictly with the ANSI standards to ensure maximum compatibility and portability. It is being publicly released, with the idea of facilitating future branching by other groups to augment its capabilities.

Twist, Writhe, and Helicity in the Inner Penumbra of a Sunspot

The aim of this work is the determination of the twist, writhe, and self-magnetic helicity of penumbral filaments located in an inner sunspot penumbra. For this purpose, we inverted data taken with the spectropolarimeter on board Hinode with the SIR (Stokes Inversion based on Response function) code. For the construction of a three-dimensional geometrical model we applied a genetic algorithm minimizing the divergence of and the net magnetohydrodynamic force, consequently a force-free solution would be reached if possible. We estimated two proxies to the magnetic helicity frequently used in literature: the force-free parameter ? z and the current helicity term . We show that both proxies are only qualitative indicators of the local twist as the magnetic field in the area under study significantly departs from a force-free configuration. The local twist shows significant values only at the borders of bright penumbral filaments with opposite signs on each side. These locations are precisely correlated to large electric currents. The average twist (and writhe) of penumbral structures is very small. The spines (dark filaments in the background) show a nearly zero writhe. The writhe per unit length of the intraspines diminishes with increasing length of the tube axes. Thus, the axes of tubes related to intraspines are less wrung when the tubes are more horizontal. As the writhe of the spines is very small, we can conclude that the writhe reaches only significant values when the tube includes the border of an intraspine.

Spectropolarimetric evidence for a siphon flow along an emerging magnetic flux tube

This work has been partially funded by the Spanish Ministerio de Economia y Competitividad, through Projects No. ESP2013-47349-C6 and ESP2014-56169-C6, including a percentage from European FEDER funds. The German contribution to sunrise and its reflight was funded by the Max Planck Foundation, the Strategic Innovations Fund of the President of the Max Planck Society (MPG), DLR, and private donations by supporting members of the Max Planck Society, which is gratefully acknowledged. The National Solar Observatory (NSO) is operated by the Association of Universities for Research in Astronomy (AURA) Inc. under a cooperative agreement with the National Science Foundation. The HAO contribution was partly funded through NASA grant number NNX13AE95G. This work was partly supported by the BK21 plus program through the National Research Foundation (NRF) funded by the Ministry of Education of Korea.

Analysis of a Spatially Deconvolved Solar Pore

Solar pores are active regions with large magnetic field strengths and apparent simple magnetic configurations. Their properties resemble the ones found for the sunspot umbra although pores do not show penumbra. Therefore, solar pores present themselves as an intriguing phenomenon that is not completely understood. We examine in this work a solar pore observed with Hinode/SP using two state of the art techniques. The first one is the spatial deconvolution of the spectropolarimetric data that allows removing the stray light contamination induced by the spatial point spread function of the telescope. The second one is the inversion of the Stokes profiles assuming local thermodynamic equilibrium that let us to infer the atmospheric physical parameters. After applying these techniques, we found that the spatial deconvolution method does not introduce artefacts, even at the edges of the magnetic structure, where large horizontal gradients are detected on the atmospheric parameters. Moreover, we also describe the physical properties of the magnetic structure at different heights finding that, in the inner part of the solar pore, the temperature is lower than outside, the magnetic field strength is larger than 2 kG and unipolar, and the LOS velocity is almost null. At neighbouring pixels, we found low magnetic field strengths of same polarity and strong downward motions that only occur at the low photosphere, below the continuum optical depth

Detection of emission in the Si i 1082.7 nm line core in sunspot umbrae

\log \tau=-1

Anomalous Polarization Profiles in Sunspots: Possible Origin of Umbral Flashes

. Finally, we studied the spatial relation between different atmospheric parameters at different heights corroborating the physical properties described before.

A GEOMETRICAL HEIGHT SCALE FOR SUNSPOT PENUMBRAE

Context. Determining empirical atmospheric models for the solar chromosphere is difficult since it requires the observation and analysis of spectral lines that are affected by non-local thermodynamic equilibrium (NLTE) effects. This task is especially difficult in sunspot umbrae because of lower continuum intensity values in these regions with respect to the surrounding brighter granulation. Umbral data is therefore more strongly affected by the noise and by the so-called scattered light, among other effects. Aims. The purpose of this study is to analyze spectropolarimetric sunspot umbra observations taken in the near-infrared Si I 1082.7 nm line taking NLTE effects into account. Interestingly, we detected emission features at the line core of the Si I 1082.7 nm line in the sunspot umbra. Here we analyze the data in detail and offer a possible explanation for the Si I 1082.7 nm line emission. Methods. Full Stokes measurements of a sunspot near disk center in the near-infrared spectral range were obtained with the GRIS instrument installed at the German GREGOR telescope. A point spread function (PSF) including the effects of the telescope, the Earth’s atmospheric seeing, and the scattered light was constructed using prior Mercury observations with GRIS and the information provided by the adaptive optics system of the GREGOR telescope during the observations. The data were then deconvolved from the PSF using a principal component analysis deconvolution method and were analyzed via the NICOLE inversion code, which accounts for NLTE effects in the Si I 1082.7 nm line. The information of the vector magnetic field was included in the inversion process. Results. The Si I 1082.7 nm line seems to be in emission in the umbra of the observed sunspot after the effects of scattered light (stray light coming from wide angles) are removed. We show how the spectral line shape of umbral profiles changes dramatically with the amount of scattered light. Indeed, the continuum levels range, on average, from 44% of the quiet Sun continuum intensity to about 20%. Although very low, the inferred levels are in line with current model predictions and empirical umbral models. The Si I 1082.7 nm line is in emission after adding more that 30% of scattered light so that it is very sensitive to a proper determination of the PSF. Additionally, we have thoroughly investigated whether the emission is a byproduct of the particular deconvolution technique but have not found any evidence to the contrary. Only the circular polarization signals seem to be more sensitive to the deconvolution strategy because of the larger amount of noise in the umbra. Interestingly, current umbral empirical models are not able to reproduce the emission in the deconvolved umbral Stokes profiles. The results of the NLTE inversions suggests that to obtain the emission in the Si I 1082.7 nm line, the temperature stratification should first have a hump located at about log  τ = −2 and start rising at lower heights when moving into the transition region. Conclusions. This is, to our knowledge, the first time the Si I 1082.7 nm line is seen in emission in sunspot umbrae. The results show that the temperature stratification of current umbral models may be more complex than expected with the transition region located at lower heights above sunspot umbrae. Our finding might provide insights into understanding why the sunspot umbra emission in the millimeter spectral range is less than that predicted by current empirical umbral models.