• Of the two solar lines, K I $D_1$ and $D_2$, almost all attention so far has been devoted to the $D_1$ line, as $D_2$ is severely affected by an O$_2$ atmospheric band. This, however, makes the latter appealing for balloon and space observations from above (most of) the Earth's atmosphere. We estimate the residual effect of the O$_2$ band on the K I $D_2$ line at altitudes typical for stratospheric balloons. Our aim is to study the feasibility of observing the 770 nm window. Specifically, this paper serves as a preparation for the third flight of the Sunrise balloon-borne observatory. The results indicate that the absorption by O$_2$ is still present, albeit much weaker, at the expected balloon altitude. We applied the obtained O$_2$ transmittance to K I $D_2$ synthetic polarimetric spectra and found that in the absence of line-of-sight motions, the residual O$_2$ has a negligible effect on the K I $D_2$ line. On the other hand, for Doppler-shifted K I $D_2$ data, the residual O$_2$ might alter the shape of the Stokes profiles. However, the residual O$_2$ absorption is sufficiently weak at stratospheric levels that it can be divided out if appropriate measurements are made, something that is impossible at ground level. Therefore, for the first time with Sunrise III, we will be able to perform polarimetric observations of the K I $D_2$ line and, consequently, we will have improved access to the thermodynamics and magnetic properties of the upper photosphere from observations of the K I lines.
  • We analyze spectropolarimetric sunspot umbra observations taken in the near-infrared Si i 1082.7 nm line taking NLTE effects into account. The data were obtained with the GRIS instrument installed at the German GREGOR telescope. A point spread function (PSF) was constructed using prior Mercury observations with GRIS and the information provided by the adaptive optics system of the GREGOR telescope. The data were then deconvolved from the PSF using a principal component analysis deconvolution method and were analyzed via the NICOLE inversion code. 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 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%. The inferred levels are in line with current model predictions and empirical umbral models. 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 tau -2 and start rising at lower heights when moving into the transition region. 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.
  • In this publication we continue the work started in Quintero Noda et al. (2017) examining this time a numerical simulation of a magnetic flux tube concentration. Our goal is to study if the physical phenomena that take place in it, in particular, the magnetic pumping, leaves a specific imprint on the examined spectral lines. We find that the profiles from the interior of the flux tube are periodically dopplershifted following an oscillation pattern that is also reflected in the amplitude of the circular polarization signals. In addition, we analyse the properties of the Stokes profiles at the edges of the flux tube discovering the presence of linear polarization signals for the Ca II lines, although they are weak with an amplitude around 0.5% of the continuum intensity. Finally, we compute the response functions to perturbations in the longitudinal field and we estimate the field strength using the weak field approximation. Our results indicate that the height of formation of the spectral lines changes during the magnetic pumping process which makes the interpretation of the inferred magnetic field strength and its evolution more difficult. These results complement those from previous works demonstrating the capabilities and limitations of the 850 nm spectrum for chromospheric Zeeman polarimetry in a very dynamic and complex atmosphere.
  • We characterize the K I D1 & D2 lines in order to determine whether they could complement the 850 nm window, containing the Ca II infrared triplet lines and several Zeeman sensitive photospheric lines, that was studied previously. We investigate the effect of partial redistribution on the intensity profiles, their sensitivity to changes in different atmospheric parameters, and the spatial distribution of Zeeman polarization signals employing a realistic magnetohydrodynamic simulation. The results show that these lines form in the upper photosphere at around 500 km and that they are sensitive to the line of sight velocity and magnetic field strength at heights where neither the photospheric lines nor the Ca II infrared lines are. However, at the same time, we found that their sensitivity to the temperature essentially comes from the photosphere. Then, we conclude that the K I lines provide a complement to the lines in the 850 nm window for the determination of atmospheric parameters in the upper photosphere, especially for the line of sight velocity and the magnetic field.
  • We study the dynamics and topology of an emerging magnetic flux concentration using high spatial resolution spectropolarimetric data acquired with the Imaging Magnetograph eXperiment on board the Sunrise balloon-borne solar observatory. We obtain the full vector magnetic field and the line-of-sight (LOS) velocity through inversions of the Fe I line at 525.02 nm with the SPINOR code. The derived vector magnetic field is used to trace magnetic field lines. Two magnetic flux concentrations with different polarity and LOS velocities are found to be connected by a group of arch-shaped magnetic field lines. The positive polarity footpoint is weaker (1100 G) and displays an upflow, while the negative polarity footpoint is stronger (2200 G) and shows a downflow. This configuration is naturally interpreted as a siphon flow along an arched magnetic flux tube.
  • 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 $\log \tau=-1$. Finally, we studied the spatial relation between different atmospheric parameters at different heights corroborating the physical properties described before.
  • Polar faculae are bright features that can be detected in solar limb observations and they are related to magnetic field concentrations. Although there is a large number of works studying them, some questions about their nature as their magnetic properties at different heights are still open. Thus, we aim to improve the understanding of solar polar faculae. In that sense, we infer the vertical stratification of the temperature, gas pressure, line of sight velocity and magnetic field vector of polar faculae regions. We performed inversions of the Stokes profiles observed with Hinode/SP after removing the stray light contamination produced by the spatial point spread function of the telescope. Moreover, after solving the azimuth ambiguity, we transform the magnetic field vector to local solar coordinates. The obtained results reveal that the polar faculae are constituted by hot plasma with low line of sight velocities and single polarity magnetic fields in the kilogauss range that are nearly perpendicular to the solar surface. We also found that the spatial location of these magnetic fields is slightly shifted respect to the continuum observations towards the disc centre. We believe that this is due to the hot wall effect that allows detecting photons that come from deeper layers located closer to the solar limb.
  • Observations of the Sun from the Earth are always limited by the presence of the atmosphere, which strongly disturbs the images. A solution to this problem is to place the telescopes in space satellites, which produce observations without any (or limited) atmospheric aberrations. However, even though the images from space are not affected by atmospheric seeing, the optical properties of the instruments still limit the observations. In the case of diffraction limited observations, the PSF establishes the maximum allowed spatial resolution, defined as the distance between two nearby structures that can be properly distinguished. In addition, the shape of the PSF induce a dispersion of the light from different parts of the image, leading to what is commonly termed as stray light or dispersed light. This effect produces that light observed in a spatial location at the focal plane is a combination of the light emitted in the object at relatively distant spatial locations. We aim to correct the effect produced by the telescope's PSF using a deconvolution method, and we decided to apply the code on Hinode/SP quiet Sun observations. We analyze the validity of the deconvolution process with noisy data and we infer the physical properties of quiet Sun magnetic elements after the deconvolution process.
  • With the advent of a new generation of solar telescopes and instrumentation, the interpretation of chromospheric observations (in particular, spectro-polarimetry) requires new, suitable diagnostic tools. This paper describes a new code, NICOLE, that has been designed for Stokes non-LTE radiative transfer, both 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 fosters 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 being publicly released, with the idea of facilitating future branching by other groups to augment its capabilities.
  • We have examined the spatial and polarimetric properties of these events using a variety of data from the Hinode spacecraft. We have also inferred the atmospheric stratification of the physical parameters by means of the inversion of the observed Stokes profiles employing the Stokes Inversion based on Response functions (SIR) code. Finally, we analyzed their evolution using a time series from the same instrument. Blue-shifted events tend to appear over bright regions at the edge of granules, while red-shifted events are seen predominantly over dark regions on intergranular lanes. Large linear polarization signals can be seen in the region that connects them. The magnetic structure inferred from the time series revealed that the structure corresponds to a $\Omega$-loop, with one footpoint always over the edge of a granule and the other inside an intergranular lane. The physical parameters obtained from the inversions of the observed Stokes profiles in both events show an increase with respect to the Harvard-Smithonian reference atmosphere in the temperature at $\log\tau_{500} \in (-1, -3)$ and a strong magnetic field, $B \ge 1$ kG, at the bottom of the atmosphere that quickly decreases upward until vanishing at $\log\tau_{500} \approx -2$. In the blue-shifted events, the line of sight velocities change from upflows at the bottom to downflows at the top of the atmosphere. Red-shifted events display the opposite velocity stratification. The change of sign in line of sight velocity happens at the same optical depth in which the magnetic field becomes zero. The physical mechanism that best explains the inferred magnetic field configuration and flow motions is a siphon flow along an arched magnetic flux tube. Further investigation is required however, as the expected features of a siphon flow cannot be unequivocally identified.
  • We study the presence of reversed polarity magnetic flux in sunspot penumbra. We applied a new regularized method to deconvolve spectropolarimetric data observed with the spectropolarimeter SP onboard Hinode. The new regularization is based on a principal component decomposition of the Stokes profiles. The resulting Stokes profiles were inverted to infer the magnetic field vector using SIR. We find, for the first time, reversed polarity fields at the border of many bright penumbral filaments in the whole penumbra.
  • A 50 min time series of one-dimensional slit-spectrograms, taken in quiet sun at disk center, observed at the German Vacuum Tower Telescope (Observatorio del Teide), was used to study the global and spatial variations of different line parameters. In order to determine the vertical structure of the photosphere two lines with well separated formation heights have been considered. The data have been filtered of p-modes to isolate the pure convective phenomenon. From our studies of global correlation coefficients and coherence and phase shift analyzes between the several line parameters, the following results can be reported. The convective velocity pattern preserves structures larger than 1.0" up to the highest layers of the photosphere (~ 435 km). However, at these layers, in the intensity pattern only structures larger than 2.0" are still connected with those at the continuum level although showing inverted brightness contrast. This confirms an inversion of temperature that we have found at a height of ~140 km. A possible evidence of gravity waves superimposed to the convective motions is derived from the phase shift analysis. We interpret the behavior of the full width at half maximum and the equivalent width as a function of the distance to the granular borders, as a consequence of enhanced turbulence and/or strong velocity gradients in the intergranular lanes.
  • From the inversion of a time series of high resolution slit spectrograms obtained from the quiet sun, the spatial and temporal distribution of the thermodynamical quantities and the vertical flow velocity is derived as a function of logarithmic optical depth and geometrical height. Spatial coherence and phase shift analyzes between temperature and vertical velocity depict the height variation of these physical quantities for structures of different size. An average granular cell model is presented, showing the granule-intergranular lane stratification of temperature, vertical velocity, gas pressure and density as a function of logarithmic optical depth and geometrical height. Studies of a specific small and a specific large granular cell complement these results. A strong decay of the temperature fluctuations with increasing height together with a less efficient penetration of smaller cells is revealed. The T -T coherence at all granular scales is broken already at log tau =-1 or z~170 km. At the layers beyond, an inversion of the temperature contrast is revealed. Vertical velocities are in phase throughout the photosphere and penetrate into the highest layers under study.
  • We determine the entire electrical current density vector in a geometrical 3D volume of the inner penumbra of a sunspot from an inversion of spectropolarimetric data obtained with Hinode/SP. Significant currents are seen to wrap around the hotter, more elevated regions with lower and more horizontal magnetic field that harbor strong upflows and radial outflows (the intraspines). The horizontal component of the current density vector is 3-4 times larger than the vertical; nearly all previous studies only obtain the vertical component and thus strongly underestimate the current density. The current density and the magnetic field vectors form an angle of about 20 degrees. The plasma beta at the 0 km level is larger than 1 in the intraspines and is one order of magnitude lower in the background component of the penumbra (spines). At the 200 km level, the plasma beta is below 0.3 nearly everywhere. The plasma beta surface as well as the surface optical depth unity are very corrugated. At the borders of intraspines and inside, the magnetic field is not force-free at deeper layers and nearly force free at the top layers. The magnetic field of the spines is close to being potential everywhere. The dissipated ohmic energy is five orders of magnitudes smaller than the solar energy flux and thus negligible for the energy balance of the penumbra.
  • Inversions of spectropolarimetric observations of penumbral filaments deliver the stratification of different physical quantities in an optical depth scale. However, without establishing a geometrical height scale their three-dimensional geometrical structure can not be derived. This is crucial in understanding the correct spatial variation of physical properties in the penumbral atmosphere and to provide insights into the mechanism capable of explaining the observed penumbral brightness. The aim of this work is to determine a global geometrical height scale in the penumbra by minimizing the divergence of the magnetic field vector and the deviations from static equilibrium as imposed by a force balance equation that includes pressure gradients, gravity and the Lorentz force. Optical depth models are derived from the SIR inversion of spectropolarimetric data of an active region observed with SOT on-board the Hinode satellite. We use a genetic algorithm to determine the boundary condition for the inference of geometrical heights. The retrieved geometrical height scale permits the evaluation of the Wilson depression at each pixel and the correlation of physical quantities at each height. Our results fit into the uncombed penumbral scenario, i.e., a penumbra composed of flux tubes with channelled mass flow and with a weaker and more horizontal magnetic field as compared with the background field. The ascending material is hotter and denser than their surroundings. We do not find evidence of overturning convection or field free regions in the inner penumbral area analyzed. The penumbral brightness can be explained by the energy transfer of the ascending mass carried by the Evershed flow, if the physical quantities below z=-75km are extrapolated from the results of the inversion.
  • We report on the discovery of mostly horizontal field channels just outside sunspot penumbrae (in the so-called `moat' region) that are seen to sustain supersonic flows (line-of-sight component of 6 km s{-1}). The spectral signature of these supersonic flows corresponds to circular polarization profiles with an additional, satellite, third lobe of the same sign as the parent sunspot' Stokes V blue lobe, for both downflows and upflows. This is consistent with an outward directed flow that we interpret as the continuation of the magnetized Evershed flow outside sunspots at supersonic speeds. In Stokes Q and U, a clear signature of a transverse field connecting the two flow streams is observed. Such an easily detectable spectral signature should allow for a clear identification of these horizontal field channels in other spectropolarimetric sunspot data. For the spot analyzed in this paper, a total of 5 channels with this spectral signature have been unambiguously found.
  • We analyse spectropolarimetric observations of the penumbra of the NOAA AR 10953 at high spatial resolution (0.3"). The full Stokes profiles of the Fe I lines at 630.1 nm and 630.2 nm have been obtained with the Solar Optical Telescope (SOT) on board the Hinode satellite. The data have been inverted by means of the SIR code, deriving the stratifications of temperature, line of sight velocity, and the components of the magnetic field vector in optical depth. In order to evaluate the gas pressure and to obtain an adequate geometrical height scale, the motion equation has been integrated for each pixel taking into account the terms of the Lorentz force. To establish the boundary condition, a genetic algorithm has been applied. The final resulting magnetic field has a divergence compatible with 0 inside its uncertainties. First analyses of the correlation of the Wilson depression with velocity, temperature, magnetic field strength, and field inclination strongly support the uncombed penumbral model proposed by Solanki & Montavon (1993).
  • High-resolution observations of sunspots have revealed the existence of dark cores inside the bright filaments of the penumbra. Here we present the stationary solution of the heat transfer equation in a stratified penumbra consisting of nearly horizontal magnetic flux tubes embedded in a stronger and more vertical field. The tubes and the external medium are in horizontal mechanical equilibrium. This model produces bright filaments with dark cores as a consequence of the higher density of the plasma inside the flux tube, which shifts the surface of optical depth unity toward higher (cooler) layers. Our results suggest that the surplus brightness of the penumbra is a natural consequence of the Evershed flow, and that magnetic flux tubes about 250 km in diameter can explain the morphology of sunspot penumbra.
  • Context: Observations at 0.1" have revealed the existence of dark cores in the bright filaments of sunspot penumbrae. Expectations are high that such dark-cored filaments are the basic building blocks of the penumbra, but their nature remains unknown. Aims: We investigate the origin of dark cores in penumbral filaments and the surplus brightness of the penumbra. To that end we use an uncombed penumbral model. Methods: The 2D stationary heat transfer equation is solved in a stratified atmosphere consisting of nearly horizontal magnetic flux tubes embedded in a stronger and more vertical field. The tubes carry an Evershed flow of hot plasma. Results: This model produces bright filaments with dark cores as a consequence of the higher density of the plasma inside the tubes, which shifts the surface of optical depth unity toward higher (cooler) layers. Our calculations suggest that the surplus brightness of the penumbra is a natural consequence of the Evershed flow, and that magnetic flux tubes about 250 km in diameter can explain the morphology of sunspot penumbrae.
  • We study the contradictory magnetic field strength distributions retrieved from independent analyses of spectropolarimetric observations in the near-infrared (1.56 micron) and in the visible (630 nm) at internetwork regions. In order to solve this apparent controversy, we present simultaneous and co-spatial 1.56 micron and 630 nm observations of an internetwork area. The properties of the circular and linear polarization signals, as well as the Stokes V area and amplitude asymmetries, are discussed. As a complement, inversion techniques are also used to infer the physical parameters of the solar atmosphere. As a first step, the infrared and visible observations are analysed separately to check their compatibility. Finally, the simultaneous inversion of the two data sets is performed. The magnetic flux densities retrieved from the individual analysis of the infrared and visible data sets are strongly correlated. The polarity of the Stokes V profiles is the same at co-spatial pixels in both wavelength ranges. This indicates that both 1.56 micron and 630 nm observations trace the same magnetic structures on the solar surface. The simultaneous inversion of the two pairs of lines reveals an internetwork full of sub-kG structures that fill only 2 % of the resolution element. A correlation is found between the magnetic field strength and the continuum intensity: equipartition fields (B ~500 G) tend to be located in dark intergranular lanes, whereas weaker field structures are found inside granules. The most probable unsigned magnetic flux density is 10 Mx/cm2. The net magnetic flux density in the whole field of view is nearly zero. This means that both polarities cancel out almost exactly in our observed internetwork area.
  • The reliability of quiet Sun magnetic field diagnostics based on the \ion{Fe}{1} lines at 6302 \AA has been questioned by recent work. We present here the results of a thorough study of high-resolution multi-line observations taken with the new spectro-polarimeter SPINOR, comprising the 5250 and 6302 \AA spectral domains. The observations were analyzed using several inversion algorithms, including Milne-Eddington, LTE with 1 and 2 components, and MISMA codes. We find that the line-ratio technique applied to the 5250 \AA lines is not sufficiently reliable to provide a direct magnetic diagnostic in the presence of thermal fluctuations and variable line broadening. In general, one needs to resort to inversion algorithms, ideally with realistic magneto-hydrodynamical constrains. When this is done, the 5250 \AA lines do not seem to provide any significant advantage over those at 6302 \AA . In fact, our results point towards a better performance with the latter (in the presence of turbulent line broadening). In any case, for very weak flux concentrations, neither spectral region alone provides sufficient constraints to fully disentangle the intrinsic field strengths. Instead, we advocate for a combined analysis of both spectral ranges, which yields a better determination of the quiet Sun magnetic properties. Finally, we propose the use of two other \ion{Fe}{1} lines (at 4122 and 9000 \AA) with identical line opacities that seem to work much better than the others.
  • The aim of this work is to study the structure of the magnetic field vector in the internetwork and search for the presence of small-scale loops. We invert 1.56 micron spectropolarimetric observations of internetwork regions at disc centre by applying the SIR code. This allows us to recover the atmospheric parameters that play a role in the formation of these spectral lines. We are mainly interested in the structure of the magnetic field vector. We find that many opposite polarity elements of the internetwork are connected by short (2-6''), low-lying (photospheric) loops. These loops connect at least the 10-20 % of the internetwork flux visible in our data. Also we have some evidence that points towards a dynamic scenario which can be produced by the emergence of internetwork magnetic flux.
  • The purpose of this work is to analyze the reliability of the magnetic field strengths inferred from the 630 nm pair of Fe I lines at internetwork quiet Sun regions. Some numerical experiments have been performed that demonstrate the inability of these lines to recover the magnetic field strength in such low flux solar regions. It is shown how different model atmospheres, with magnetic field strengths ranging from few hundred Gauss to kiloGauss, give rise to Stokes profiles that can not be distinguished. The reasons for this degeneracy are discussed.
  • An inversion technique has been developed to recover LTE one-dimensional model photospheres for late-type stars from very high-resolution high signal-to-noise stellar line profiles. It is successfully applied to the Sun using a set of unblended Ti I, Ca I, Cr I and Fe I lines with accurate transition probabilities. Temperature stratification, continuum flux, centre-to-limb variation and wings of strong metal lines obtained from the resulting model are compared with those from other well-known theoretical and empirical solar models and show the reliability of the procedure.
  • First results of the inversion of Stokes I and V profiles from plage regions near disk center are presented. Both low and high spatial resolution spectra of FeI 6301.5 and FeI 6302.5 A obtained with the Advanced Stokes Polarimeter (ASP) have been considered for analysis. The thin flux tube approximation, implemented in an LTE inversion code based on response functions, is used to describe unresolved magnetic elements. The code allows the simultaneous and consistent inference of all atmospheric quantities determining the radiative transfer with the sole assumption of hydrostatic equilibrium. By considering velocity gradients within the tubes we are able to match the full ASP Stokes profiles. The magnetic atmospheres derived from the inversion are characterized by the absence of significant motions in high layers and strong velocity gradients in deeper layers. These are essential to reproduce the asymmetries of the observed profiles. Our scenario predicts a shift of the Stokes V zero-crossing wavelengths which is indeed present in observations made with the Fourier Transform Spectrometer.