In three-dimensional (3D) radiative transfer (RT) problems, the tensor product quadratures are generally not optimal in terms of the number of discrete ray directions needed for a given accuracy of the angular integration of the radiation field. In this paper, we derive a new set of angular quadrature rules that are more suitable for solving 3D RT problems with the short- and long-characteristics formal solvers. These quadratures are more suitable than the currently used ones for the numerical calculation of the radiation field tensors that are relevant in the problem of the generation and transfer of polarised radiation without assuming local thermodynamical equilibrium (non-LTE). We show that our new quadratures can save up to about 30\,\% of computing time with respect to the Gaussian-trapezoidal product quadratures with the same accuracy.

We investigated the off-limb spicules observed in the Mg II h and k lines by IRIS in a solar polar coronal hole. We analyzed the large dataset of obtained spectra to extract quantitative information about the line intensities, shifts, and widths. The observed Mg II line profiles are broad and double-peaked at lower altitudes, broad but flat-topped at middle altitudes, and narrow and single-peaked with the largest Doppler shifts at higher altitudes. We use 1D non-LTE vertical slab models (i.e. models which consider departures from Local Thermodynamic Equilibrium) in single-slab and multi-slab configurations to interpret the observations and to investigate how a superposition of spicules along the line of sight (LOS) affects the synthetic Mg II line profiles. The used multi-slab models are either static, i.e. without any LOS velocities, or assume randomly assigned LOS velocities of individual slabs, representing the spicule dynamics. We conducted such single-slab and multi-slab modeling for a broad set of model input parameters and showed the dependence of the Mg II line profiles on these parameters. We demonstrated that the observed line widths of the h and k line profiles are strongly affected by the presence of multiple spicules along the LOS. We later showed that the profiles obtained at higher altitudes can be reproduced by single-slab models representing individual spicules. We found that the multi-slab model with a random distribution of the LOS velocities ranging from -25 to 25 km s$^{-1}$ can well reproduce the width and the shape of Mg II profiles observed at middle altitudes.

Interpreting the Stokes profiles observed in quiet regions of the solar chromosphere is a challenging task. The Stokes Q and U profiles are dominated by the scattering polarisation and the Hanle effect, and these processes can only be correctly quantified if 3D radiative transfer effects are taken into account. Forward-modelling of the intensity and polarisation of spectral lines using a 3D model atmosphere is a suitable approach in order to statistically compare the theoretical and observed line profiles. Our aim is to present novel observations of the Ca 8542 line profiles in a quiet region at the centre of the solar disc and to quantitatively compare them with the theoretical Stokes profiles. We aim at estimating the reliability of the 3D model atmosphere using not only the line intensity but the full vector of Stokes parameters. We used data obtained with the ZIMPOL instrument at the IRSOL and compared the observations with the theoretical profiles computed with the PORTA radiative transfer code, using as solar model atmosphere a 3D snapshot taken from a radiation-magnetohydrodynamics simulation. The synthetic profiles were degraded to match the instrument and observing conditions. The degraded theoretical profiles of the Ca 8542 line are qualitatively similar to the observed ones. We confirm that there is a fundamental difference in the widths of all Stokes profiles: the observed lines are wider than the theoretical lines. We find that the amplitudes of the observed profiles are larger than those of the theoretical ones, which suggests that the symmetry breaking effects in the solar chromosphere are stronger than in the model atmosphere. This means that the isosurfaces of temperature, velocity, and magnetic field strength and orientation are more corrugated in the solar chromosphere than in the currently available 3D radiation-magnetohydrodynamics simulation.

On 3rd September 2015, the Chromospheric Lyman-Alpha SpectroPolarimeter (CLASP) successfully measured the linear polarization produced by scattering processes in the hydrogen Lyman-$\alpha$ line of the solar disk radiation, revealing conspicuous spatial variations in the $Q/I$ and $U/I$ signals. Via the Hanle effect the line-center $Q/I$ and $U/I$ amplitudes encode information on the magnetic field of the chromosphere-corona transition region (TR), but they are also sensitive to the three-dimensional structure of this corrugated interface region. With the help of a simple line formation model, here we propose a statistical inference method for interpreting the Lyman-$\alpha$ line-center polarization observed by CLASP.

Probing the magnetism of the upper solar chromosphere requires measuring and modeling the scattering polarization produced by anisotropic radiation pumping in UV spectral lines. Here we apply PORTA (a novel radiative transfer code) to investigate the hydrogen Ly$\alpha$ line in a 3D model of the solar atmosphere resulting from a state of the art MHD simulation. At full spatial resolution the linear polarization signals are very significant all over the solar disk, with a large fraction of the field of view showing line-center amplitudes well above the 1% level. Via the Hanle effect the line-center polarization signals are sensitive to the magnetic field of the model's transition region, even when its mean field strength is only 15 G. The breaking of the axial symmetry of the radiation field produces significant forward-scattering polarization in Ly$\alpha$, without the need of an inclined magnetic field. Interestingly, the Hanle effect tends to decrease such forward-scattering polarization signals in most of the points of the field of view. When the spatial resolution is degraded, the line-center polarization of Ly$\alpha$ drops below the 1% level, reaching values similar to those previously found in 1D semi-empirical models (i.e., up to about 0.5%). The center to limb variation of the spatially-averaged polarization signals is qualitatively similar to that found in 1D models, with the largest line-center amplitudes at $\mu=\cos\theta\approx 0.4$ ($\theta$ being the heliocentric angle). These results are important, both for designing the needed space-based instrumentation and for a reliable interpretation of future observations of the Ly$\alpha$ polarization.

Magnetic field measurements in the upper chromosphere and above, where the gas-to-magnetic pressure ratio $\beta$ is lower than unity, are essential for understanding the thermal structure and dynamical activity of the solar atmosphere. Recent developments in the theory and numerical modeling of polarization in spectral lines have suggested that information on the magnetic field of the chromosphere-corona transition region could be obtained by measuring the linear polarization of the solar disk radiation at the core of the hydrogen Lyman-$\alpha$ line at 121.6~nm, which is produced by scattering processes and the Hanle effect. The Chromospheric Lyman-$\alpha$ Spectropolarimeter (CLASP) sounding rocket experiment aims to measure the intensity (Stokes $I$) and the linear polarization profiles ($Q/I$ and $U/I$) of the hydrogen Lyman-$\alpha$ line. In this paper we clarify the information that the Hanle effect can provide by applying a Stokes inversion technique based on a database search. The database contains all theoretical $Q/I$ and $U/I$ profiles calculated in a one-dimensional semi-empirical model of the solar atmosphere for all possible values of the strength, inclination, and azimuth of the magnetic field vector, though this atmospheric region is highly inhomogeneous and dynamic. We focus on understanding the sensitivity of the inversion results to the noise and spectral resolution of the synthetic observations as well as the ambiguities and limitation inherent to the Hanle effect when only the hydrogen Lyman-$\alpha$ is used. We conclude that spectropolarimetric observations with CLASP can indeed be a suitable diagnostic tool for probing the magnetism of the transition region, especially when complemented with information on the magnetic field azimuth that can be obtained from other instruments.

There is an ongoing debate about the origin and even the very existence of a high degree of linear polarization of some chromospheric spectral lines observed in solar flares. The standard explanation of these measurements is in terms of the impact polarization caused by non-thermal proton and/or electron beams. In this work, we study the possible role of resonance line polarization due to radiation anisotropy in the inhomogeneous medium of the flare ribbons. We consider a simple two-dimensional model of the flaring chromosphere and we solve self-consistently the non-LTE problem taking into account the role of resonant scattering polarization and of the Hanle effect. Our calculations show that the horizontal plasma inhomogeneities at the boundary of the flare ribbons can lead to a significant radiation anisotropy in the line formation region and, consequently, to a fractional linear polarization of the emergent radiation of the order of several percent. Neglecting the effects of impact polarization, our model can provide a clue for resolving some of the common observational findings, namely: (1) why a high degree of polarization appears mainly at the edges of the flare ribbons; (2) why polarization can also be observed during the gradual phase of a flare; (3) why polarization is mostly radial or tangential. We conclude that the radiation transfer in the realistic multi-dimensional models of solar flares needs to be considered as an essential ingredient for understanding the observed spectral line polarization.

The interpretation of the intensity and polarization of the spectral line radiation produced in the atmosphere of the Sun and of other stars requires solving a radiative transfer problem that can be very complex, especially when the main interest lies in modeling the spectral line polarization produced by scattering processes and the Hanle and Zeeman effects. One of the difficulties is that the plasma of a stellar atmosphere can be highly inhomogeneous and dynamic, which implies the need to solve the non-equilibrium problem of the generation and transfer of polarized radiation in realistic three-dimensional (3D) stellar atmospheric models. Here we present PORTA, an efficient multilevel radiative transfer code we have developed for the simulation of the spectral line polarization caused by scattering processes and the Hanle and Zeeman effects in 3D models of stellar atmospheres. The numerical method of solution is based on the non-linear multigrid iterative method and on a novel short-characteristics formal solver of the Stokes-vector transfer equation which uses monotonic B\'ezier interpolation. Therefore, with PORTA the computing time needed to obtain at each spatial grid point the self-consistent values of the atomic density matrix (which quantifies the excitation state of the atomic system) scales linearly with the total number of grid points. Another crucial feature of PORTA is its parallelization strategy, which allows us to speed up the numerical solution of complicated 3D problems by several orders of magnitude with respect to sequential radiative transfer approaches, given its excellent linear scaling with the number of available processors. The PORTA code can also be conveniently applied to solve the simpler 3D radiative transfer problem of unpolarized radiation in multilevel systems.

In order to understand the heating of the solar corona it is crucial to obtain empirical information on the magnetic field in its lower boundary (the transition region). To this end, we need to measure and model the linear polarization produced by scattering processes in strong UV lines, such as the hydrogen Ly$\alpha$ line. The interpretation of the observed Stokes profiles will require taking into account that the outer solar atmosphere is highly structured and dynamic, and that the height of the transition region may well vary from one place in the atmosphere to another. Here we report on the Ly$\alpha$ scattering polarization signals we have calculated in a realistic model of an enhanced network region, resulting from a state-of-the-art radiation MHD simulation. This model is characterized by spatially complex variations of the physical quantities at transition region heights. The results of our investigation lead us to emphasize that scattering processes in the upper solar chromosphere should indeed produce measurable linear polarization in Ly$\alpha$. More importantly, we show that via the Hanle effect the model's magnetic field produces significant changes in the emergent $Q/I$ and $U/I$ profiles. Therefore, we argue that by measuring the polarization signals produced by scattering processes and the Hanle effect in Ly$\alpha$ and contrasting them with those computed in increasingly realistic atmospheric models, we should be able to decipher the magnetic, thermal and dynamic structure of the upper chromosphere and transition region of the Sun.

Recent theoretical investigations have pointed out that the cores of the Ly-alpha lines of H I and He II should show measurable scattering polarization signals when observing the solar disk, and that the magnetic sensitivity, through the Hanle effect, of such linear polarization signals is suitable for exploring the magnetism of the solar transition region. Such investigations were carried out in the limit of complete frequency redistribution (CRD) and neglecting quantum interference between the two upper J-levels of each line. Here we relax both approximations and show that the joint action of partial frequency redistribution (PRD) and J-state interference produces much more complex fractional linear polarization (Q/I) profiles, with large amplitudes in their wings. Such wing polarization signals turn out to be very sensitive to the temperature structure of the atmospheric model, so that they can be exploited for constraining the thermal properties of the solar chromosphere. Finally, we show that the approximation of CRD without J-state interference is however suitable for estimating the amplitude of the linear polarization signals in the core of the lines, where the Hanle effect operates.

In order to obtain empirical information on the magnetism of the solar transition region we need to measure and interpret the linear polarization produced by scattering processes in FUV and EUV spectral lines. Via the Hanle effect such linear polarization signals are sensitive to the magnetic fields expected for the quiet and active regions of the outer solar atmosphere. For example, the Ly$\alpha$ line of H\,{\sc i} at 1216\,\AA\ is mainly sensitive to magnetic strengths between 10 and 100 G. The interpretation of the observed spectral line polarization requires the development of suitable modeling tools. To this end, we have developed a three-dimensional (3D), non-LTE multilevel radiative transfer code for modeling the intensity and linear polarization produced by scattering processes in spectral lines and its modification by the Hanle effect.

The Ly-alpha line of He II at 304 \AA\ is one of the spectral lines of choice for EUV channels of narrow-band imagers on board space telescopes, which provide spectacular intensity images of the outer solar atmosphere. Since the magnetic field information is encoded in the polarization of the spectral line radiation, it is important to investigate whether the He II line radiation from the solar disk can be polarized, along with its magnetic sensitivity. Here we report some theoretical predictions concerning the linear polarization signals produced by scattering processes in this strong emission line of the solar transition region, taking into account radiative transfer and the Hanle effect caused by the presence of organized and random magnetic fields. We find that the fractional polarization amplitudes are significant (${\sim}1%$), even when considering the wavelength-integrated signals. Interestingly, the scattering polarization of the Ly-alpha line of He II starts to be sensitive to the Hanle effect for magnetic strengths B>100 G (i.e., for magnetic strengths of the order of and larger than the Hanle saturation field of the hydrogen Ly-alpha line at 1216 \AA). We therefore propose simultaneous observations of the scattering polarization in both Ly-alpha lines to facilitate magnetic field measurements in the upper solar chromosphere. Even the development of a narrow-band imaging polarimeter for the He II 304 \AA\ line alone would be already of great diagnostic value for probing the solar transition region.

We present some theoretical predictions concerning the amplitude and magnetic sensitivity of the linear polarization signals produced by scattering processes in the hydrogen Ly-alpha line of the solar transition region. To this end, we have calculated the atomic level polarization (population imbalances and quantum coherences) induced by anisotropic radiation pumping in semi-empirical and hydrodynamical models of the solar atmosphere, taking into account radiative transfer and the Hanle effect caused by the presence of organized and random magnetic fields. The line-center amplitudes of the emergent linear polarization signals are found to vary typically between 0.1% and 1%, depending on the scattering geometry and the strength and orientation of the magnetic field. The results shown here encourage the development of UV polarimeters for sounding rockets and space telescopes with the aim of opening up a diagnostic window for magnetic field measurements in the upper chromosphere and transition region of the Sun.

The Hanle effect in strong spectral lines is the physical mechanism that should be increasingly exploited for quantitative explorations of the magnetism of the quiet solar chromospheric plasma. Here we show, by means of multilevel radiative transfer calculations and new spectropolarimetric observations, that the amplitude and shape of the scattering polarization profiles of the H-alpha line is very sensitive to the strength and structure of the chromospheric magnetic field. The interpretation of the observations in terms of one-dimensional radiative transfer modeling suggests that there is an abrupt magnetization in the upper chromosphere of the quiet Sun.

Although the spectral lines of hydrogen contain valuable information on the physical properties of a variety of astrophysical plasmas, including the upper solar chromosphere, relatively little is known about their scattering polarization signals whose modification via the Hanle effect may be exploited for magnetic field diagnostics. Here we report on a basic theoretical investigation of the linear polarization produced by scattering processes and the Hanle effect in Ly-a, Ly-b and H-a taking into account multilevel radiative transfer effects in an isothermal stellar atmosphere model, the fine-structure of the hydrogen levels, as well as the impact of collisions with electrons and protons. The main aim of this first paper is to elucidate the physical mechanisms that control the linear polarization in the three lines, as well as its sensitivity to the perturbers density and to the strength and structure of micro-structured and deterministic magnetic fields. To this end, we apply an efficient radiative transfer code we have developed for performing numerical simulations of the Hanle effect in multilevel systems with overlapping line transitions. For low density plasmas such as that of the upper solar chromosphere collisional depolarization is caused mainly by collisional transitions between the fine-structure levels of n=3, so that it is virtually insignificant for Ly-a but important for Ly-b and H-a. We show the impact of the Hanle effect on the three lines taking into account the radiative transfer coupling between the different hydrogen line transitions. For example, we demonstrate that the linear polarization profile of the H-a line is sensitive to the presence of magnetic field gradients in the line core formation region and that in solar-like chromospheres selective absorption of polarization components does not play any significant role on the emergent scattering polarization.

We report on a detailed radiative transfer modeling of the observed scattering polarization in the H-alpha line, which allows us to infer quantitative information on the magnetization of the quiet solar chromosphere. Our analysis suggests the presence of a magnetic complexity zone with a mean field strength <B> > 30 G lying just below the sudden transition region to the coronal temperatures. The chromospheric plasma directly underneath is very weakly magnetized, with <B> ~ 1 G. The possible existence of this abrupt change in the degree of magnetization of the upper chromosphere of the quiet Sun might have large significance for our understanding of chromospheric (and, therefore, coronal) heating.

A particularly interesting line for exploring the physical conditions of the quiet solar chromosphere is H$\alpha$, but its intensity profile is magnetically insensitive and the small circular polarization signatures produced by the longitudinal Zeeman effect come mainly from the underlying photosphere. Here we show that the Hanle effect in H$\alpha$ provides quantitative information on the magnetism of the quiet chromosphere. To this end, we calculate the response function of the emergent scattering polarization to perturbations in the magnetic field.

We will show that isotropic collisions of electrons and protons with neutral hydrogen can lead to creation of net orientation of the atomic levels in the presence of a magnetic field. Consequently, the emitted Stokes-V profile of the spectral lines can be almost symmetric in contrast to the typical antisymmetric signature of the Zeeman effect. Moreover, the amplitude of the symmetric lobe can be significantly higher than the amplitude of the antisymmetric components. This mechanism is caused by a $\pm{M}$ symmetry breaking of the collisional transitions between different Zeeman sublevels. We will show an example of our first results for the H$\alpha$ line. This new mechanism could perhaps explain the net circular polarization of spectral lines observed in some solar limb observations and which are currently not understood. However, our results are very preliminary and more developments are needed for going further on.

Aims. We investigate the effect of the electric return currents in solar flares on the profiles of hydrogen Balmer lines. We consider the monoenergetic approximation for the primary beam and runaway model of the neutralizing return current. Methods. Propagation of the 10 keV electron beam from a coronal reconnection site is considered for the semiempirical chromosphere model F1. We estimate the local number density of return current using two approximations for beam energy fluxes between $4\times 10^{11}$ and $1\times 10^{12} {\rm erg cm^{-2} s^{-1}}$. Inelastic collisions of beam and return-current electrons with hydrogen are included according to their energy distributions, and the hydrogen Balmer line intensities are computed using an NLTE radiative transfer approach. Results. In comparison to traditional NLTE models of solar flares that neglect the return-current effects, we found a significant increase emission in the Balmer line cores due to nonthermal excitation by return current. Contrary to the model without return current, the line shapes are sensitive to a beam flux. It is the result of variation in the return-current energy that is close to the hydrogen excitation thresholds and the density of return-current electrons.

Context. A derivation of a generalized sqrt(epsilon)-law for nonthermal collisional rates of excitation by charged perturbers is presented. Aims. Aim of this paper is to find a more general analytical expression for a surface value of the source function which can be used as an addtional tool for verification of the non-LTE radiative transfer codes. Methods. Under the impact approximation hypothesis, static, one-dimensional, plane-parallel atmosphere, constant magnetic field of arbitrary strength and direction, two-level atom model with unpolarized lower level and stimulated emission neglected, we introduce the unphysical terms into the equations of statistical equilibrium and solve the appropriate non-LTE integral equations. Results. We derive a new analytical condition for the surface values of the source function components expressed in the basis of irreducible spherical tensors.

We review the problem of proton beam bombardment of solar chromosphere considering the self-consistent NLTE polarized radiation transfer in hydrogen lines. Several observations indicate a linear polarization of the H-alpha line of the order of 5% or higher and preferentially in radial direction. This polarization is often explained as anisotropic collisional excitation of the n=3 level by vertical proton beams. Our calculations indicate that deceleration of the proton beam with initial power-law energy distribution together with increased electron and proton densities in the H-alpha forming layers lead to a negligible line polarization. Thus the proton beams seem not to be a good candidate for explanation of the observed polarization degree. On the other hand, the effect of electric return currents could perhaps provide a better explanation of the observed linear polarization. We report the new calculations of this effect.

Context. We present a theoretical review of the effect of impact polarization of a hydrogen H$\alpha$ line due to an expected proton beam bombardment in solar flares. Aims. Several observations indicate the presence of the linear polarization of the hydrogen H$\alpha$ line observed near the solar limb above 5% and preferentially in the radial direction. We theoretically review the problem of deceleration of the beam originating in the coronal reconnection site due to its interaction with the chromospheric plasma, and describe the formalism of the density matrix used in our description of the atomic processes and the treatment of collisional rates. Methods. We solve the self-consistent NLTE radiation transfer problem for the particular semiempirical chromosphere models for both intensity and linear polarization components of the radiation field. Results. In contrast to recent calculations, our results show that the energy distribution of the proton beam at H$\alpha$ formation levels and depolarizing collisions by background electrons and protons cause a significant reduction of the effect below 0.1%. The radiation transfer solution shows that tangential resonance-scattering polarization dominates over the impact polarization effect in all considered models. Conclusions. In the models studied, proton beams are unlikely to be a satisfying explanation for the observed linear polarization of the H$\alpha$ line.

A new iterative method for non-LTE multilevel polarized radiative transfer in hydrogen lines is presented. Iterative methods (such as the Jacobi method) tend to damp out high-frequency components of the error fast, but converges poorly due to slow reduction of low-frequency components. The idea is to use a set of differently coarsed grids to reduce both the short- and long-period errors. This leads to the so-called multigrid (MG) methods. For the grid of $N$ spatial points, the number of iterations required to solve a non-LTE transfer problem is of the order of O(N). This fact could be of great importance for problems with fine structure and for multi-dimensional models. The efficiency of the so-called standard MG iteration in comparison to Jacobi iteration is shown. The formalism of density matrix is applied to the demonstrative example of 1D, semi-infinite, non-magnetic, 3-principal level hydrogen atmospheric model. The effect of depolarizing collisions with thermal electrons is taken into account as well as general treatment of overlapping profiles.