Interpreting spectropolarimetric observations of the solar atmosphere takes much longer than the acquiring the data. The most important reason for this is that the model fitting, or "inversion", used to infer physical quantities from the observations is extremely slow, because the underlying models are numerically demanding. We aim to improve the speed of the inference by using a neural network that relates input polarized spectra to the output physical parameters. We first select a subset of the data to be interpreted and infer physical quantities from corresponding spectra using a standard minimization-based inversion code. Taking these results as reliable and representative of the whole data set, we train a convolutional neural network to connect the input polarized spectra to the output physical parameters (nodes, in context of spectropolarimetric inversion). We then apply the neural network to the various other data, previously unseen to the network. As a check, we apply the referent inversion code to the unseen data and compare the fit quality and the maps of the inferred parameters between the two inversions. The physical parameters inferred by the neural network show excellent agreement with the results from the inversion, and are obtained in a factor of $10^5$ less time. Additionally, substituting the results of the neural network back in the forward model, shows excellent agreement between inferred and original spectra. The method we present here is very simple for implementation and extremely fast. It only requires a training data set, which can be obtained by inverting a representative subset of the observed data. Applying these (and similar) machine learning techniques will yield orders of magnitude acceleration in the routine interpretation of spectropolarimetric data.

Studying the properties of the solar convection using high-resolution spectropolarimetry began in the early 90's with the focus on observations in the visible wavelength regions. Its extension to the infrared (IR) remains largely unexplored. The IR iron lines around 15600\,$\rm{\AA}$, most commonly known for their high magnetic sensitivity, also have a non-zero response to line-of-sight velocity below $\log (\tau)=0.0$. In this paper we aim to tap this potential to explore the possibility of using them to measure sub-surface convective velocities. By assuming a snapshot of a three-dimensional magnetohydrodynamic simulation to represent the quiet Sun, we investigate how well the iron IR lines can reproduce the LOS velocity in the cube and up to what depth. We use the recently developed spectropolarimetric inversion code SNAPI and discuss the optimal node placements for the retrieval of reliable results from these spectral lines. We find that the IR iron lines can measure the convective velocities down to $\log (\tau)=0.5$, below the photosphere, not only at original resolution of the cube but also when degraded with a reasonable spectral and spatial PSF and stray light. Meanwhile, the commonly used Fe~{\sc i} 6300\,\AA{} line pair performs significantly worse. Our investigation reveals that the IR iron lines can probe the subsurface convection in the solar photosphere. This paper is a first step towards exploiting this diagnostic potential.

Inversion codes are computer programs that fit a model atmosphere to the observed Stokes spectra, thus retrieving the relevant atmospheric parameters. The rising interest in the solar chromosphere, where spectral lines are formed by scattering, requires developing, testing, and comparing new non-local thermal equilibrium (NLTE) inversion codes. We present a new NLTE inversion code that is based on the analytical computation of the response functions. We named the code SNAPI, which is short for spectropolarimetic NLTE analytically powered inversion. SNAPI inverts full Stokes spectrum in order to obtain a depth-dependent stratification of the temperature, velocity, and the magnetic field vector. It is based on the so-called node approach, where atmospheric parameters are free to vary in several fixed points in the atmosphere, and are assumed to behave as splines in between. We describe the inversion approach in general and the specific choices we have made in the implementation. We test the performance on one academic problem and on two interesting NLTE examples, the Ca\,II\,8542 and Na\,I\,D spectral lines. The code is found to have excellent convergence properties and outperforms a finite-difference based code in this specific implementation by at least a factor of three. We invert synthetic observations of Na lines from a small part of a simulated solar atmosphere and conclude that the Na lines reliably retrieve the magnetic field and velocity in the range $-3<\log \tau < -0.5$.

Most of the quantitative information about the magnetic field vector in solar prominences comes from the analysis of the Hanle effect acting on lines formed by scattering. As these lines can be of non-negligible optical thickness, it is of interest to study the line formation process further. We investigate the multidimensional effects on the interpretation of spectropolarimetric observations, particularly on the inference of the magnetic field vector. We do this by analyzing the differences between multidimensional models, which involve fully self-consistent radiative transfer computations in the presence of spatial inhomogeneities and velocity fields, and those which rely on simple one-dimensional geometry. We study the formation of a prototype line in ad hoc inhomogeneous, isothermal 2D prominence models. We solve the NLTE polarized line formation problem in the presence of a large-scale oriented magnetic field. The resulting polarized line profiles are then interpreted (i.e. inverted) assuming a simple 1D slab model. We find that differences between input and the inferred magnetic field vector are non-negligible. Namely, we almost universally find that the inferred field is weaker and more horizontal than the input field. Spatial inhomogeneities and radiative transfer have a strong effect on scattering line polarization in the optically thick lines. In real-life situations, ignoring these effects could lead to a serious misinterpretation of spectropolarimetric observations of chromospheric objects such as prominences.

One of the necessary parameters needed for the interpretation of the light curves of transiting exoplanets or eclipsing binaries, as well as interferometric measurements of a star or microlensing events is how the intensity and polarization of a light change from the center to the limb. Scattering and absorption processes in stellar atmosphere affect both the center-to limb variation of intensity (CLVI) and polarization (CLVP). In this paper, we present a study of the CLVI and CLVP in continuum spectra considering different contributions of scattering and absorption opacity for different spectral type stars with spherical atmospheres. We solve the polarized radiative transfer equation in the presence of continuum scattering, considering spherical stellar model atmospheres. We developed two independent codes based on Feautrier and short characteristics methods to cross-check our results. We calculate the CLVI and CLVP in continuum for the Phoenix grid of spherical stellar model atmospheres for a range of $T_{eff} = 4000 - 7000 \rm K$, $\log g = 1.0 - 5.5$ and $\lambda = 4000 - 7000 \rm \AA$, which are tabulated and available at the CDS. For sub-giant and dwarf stars ($\log g = 3.0 - 4.5$), lower $\log g$ and lower $T_{eff}$ of a star lead to higher limb polarization of the star. For giant and supergiant stars ($\log g = 1.0 - 2.5$), the highest effective temperature yields the largest polarization. By decreasing of the $T_{eff}$ of a star down to $4500 - 5500 \rm K$ (depending on $\log g$) the limb polarization decreases and reaches a local minimum. It increases again down to $T_{eff}$ of $4000 \rm K$. For the most compact dwarf stars ($\log g = 5.0 - 5.5$) the limb polarization degree shows a maximum for models with $T_{eff}$ in the range $4200 - 4600 \rm K$ (depending on $\log g$) and decreases toward higher and lower temperatures.

To interpret observations of astrophysical disks it is essential to understand the formation process of the emitted light. If the disk is optically thick, scattering dominated and permeated by a Keplerian velocity field, Non-Local Thermodynamic Equilibrium radiative transfer modeling must be done to compute the emergent spectrum from a given disk model. We investigate Non-local thermodynamic equilibrium polarized line formation in different simple disk models and aim to demonstrate the importance of both radiative transfer effects and scattering as well as the effects of velocity fields. We self-consistently solve the coupled equations of radiative transfer and statistical equilibrium for a two level atom model by means of Jacobi iteration. We compute scattering polarization, that is Q/I and U/I line profiles. The degree of scattering polarization is significantly influenced by the inclination of the disk with respect to observer, but also by the optical thickness of the disk and the presence of rotation. Stokes U shows double-lobed profiles with amplitude which increases with the disk rotation. Our results suggest that the line profiles, especially the polarized ones, emerging from gaseous disks differ significantly from the profiles predicted by simple approximations. The profiles are diverse in shape, but typically symmetric in Stokes Q and antisymmetric in Stokes U. A clear indicator of disk rotation is the presence of Stokes U, which might prove to be a useful diagnostic tool. We also demonstrate that, for moderate rotational velocities, an approximate treatment can be used, where non-local thermodynamic equilibrium radiative transfer is done in the velocity field-free approximation and Doppler shift is applied in the process of spatial integration over the whole emitting surface.

State-of-the-art methods in multidimensional NLTE radiative transfer are based on the use of local approximate lambda operator within either Jacobi or Gauss-Seidel iterative schemes. Here we propose another approach to the solution of 2D NLTE RT problems, Forth-and-Back Implicit Lambda Iteration (FBILI), developed earlier for 1D geometry. In order to present the method and examine its convergence properties we use the well-known instance of the two-level atom line formation with complete frequency redistribution. In the formal solution of the RT equation we employ short characteristics with two-point algorithm. Using an implicit representation of the source function in the computation of specific intensities, we compute and store the coefficients of the linear relations J = a + bS between the mean intensity J and the corresponding source function S. The use of iteration factors in the 'local' coefficients of these implicit relations in two 'inward' directions, along with the update of the source function in other two, 'outward', directions leads to four times faster solution than the Jacobi's one. Moreover, the update made in all four consecutive sweeps of the grid leads to an acceleration by a factor of 6-7 compared to the Jacobi iterative scheme.

Polarized NLTE radiative transfer in the presence of scattering in spectral lines and/or in continua may be cast in a so-called reduced form for six reduced components of the radiation field. In this formalism the six components of the reduced source function are angle-independent quantities. It thus reduces drastically the storage requirement of numerical codes. This approach encounters a fundamental problem when the medium is illuminated by a polarized incident radiation, because there is a priori no way of relating the known (and measurable) Stokes parameters of the incident radiation to boundary conditions for the reduced equations. The origin of this problem is that there is no unique way of deriving the radiation reduced components from its Stokes parameters (only the inverse operation is clearly defined). The method proposed here aims at enabling to work with arbitrary incident radiation field (polarized or unpolarized). In previous works an ad-hoc treatment of the boundary conditions, applying to cases where the incident radiation is unpolarized, has been used. In this note we show that it is possible to account for the incident radiation in a rigorous way, without any assumption on its properties, by expressing the radiation field as the sum of a directly transmitted radiation and of a so-called diffuse radiation. The diffuse radiation field obeys a transfer equation with no incident radiation that may be solved in the reduced form. The first scattering of the incident radiation introduces primary creation terms in the six components of the reduced source function. Once the reduced polarized transfer problem is solved for the diffuse radiation field, its Stokes parameters can be computed. We perform numerical computations of such cases, showing that the emergent line-polarization may be significantly affected by the polarization of the incident radiation.