EDP Sciences
Free access
Volume 559, November 2013
Article Number A68
Number of page(s) 5
Section Astrophysical processes
DOI http://dx.doi.org/10.1051/0004-6361/201322307
Published online 15 November 2013

© ESO, 2013

1. Introduction

Solving non-local thermodynamic equilibrium (NLTE) polarized radiative transfer problems in multidimensional media is a challenge from a numerical point of view. In principle, it requires much more memory and computing time than the unpolarized problem, because the source function for polarized radiation depends on both the spatial coordinates and the two angles (θ,ϕ), which define the propagation direction of the radiation. Furthermore, a polarized radiation is described by four measurable quantities, the Stokes parameters (instead of one, in the unpolarized scalar case), which are coupled by the physics of the absorption and emission processes.

In this paper, we limit ourselves to problems where the polarization arises from scattering processes and where the absorption of radiation does not introduce any coupling between the Stokes parameters. This is the case for line scattering problems when the Zeeman sublevels of the lower-level are evenly populated, and for Rayleigh or Thomson scattering problems in continua. We consider only linearly polarized radiation, as scattering processes do not couple linear and circular polarization. The Hanle effect due to the presence of weak magnetic fields is accounted for in the line scattering case.

The common point of these problems from a mathematical point of view is that the scattering can be described by a scattering matrix, which depends on both the directions and frequencies of the incoming and outgoing photons. It also depends on the magnetic field strength and direction in the presence of Hanle effect. A significant step toward the efficient solution of scattering problems in 3D media has been achieved thanks to the work of Anusha & Nagendra (2011a,b) and the subsequent series of papers by the same authors. Their work is based on the previous paper of Frisch (2007), who showed that the scattering phase matrix can be written as the sum of terms where the ingoing and outgoing directions are factorized. This factorization applies when the frequency redistribution due to scattering is independent of the directions of the photons. This is the case when the line is formed with complete frequency redistribution, whereas it relies on the assumption that the angle-averaged frequency redistribution function may be used instead of the full angle-dependent one in the presence of partial frequency redistribution. In this paper, we work in the framework of this assumption, but it could be relaxed following the same approach as in Supriya et al. (2013).

The NLTE polarized radiative transfer equation for the three coupled Stokes parameters I, Q, and U requires the computation of three source terms, depending on the position, the frequency, and two angles that define the radiation propagation direction. We see in the following that the reduced formalism allows us to deal with six reduced source terms, which do not depend on the propagation direction. In multidimensional media, the radiation field may be locally highly anisotropic and the anisotropy is an important source of polarization; fine angular grids are thus required for accurate numerical solutions of the polarized radiative transfer equation. The reduced formalism then drastically reduces the memory and computing time requirements of numerical solutions.

2. Reduced formalism for polarized radiative transfer

Let us summarize the main steps of the derivation leading to the reduced polarized radiative transfer equation.

2.1. Reduced radiative transfer equation

We assume that the scattering phase matrix may be written as a sum of factorized terms where the frequency and angular variables are separated, (1)where rII(x,x′) and rIII(x,x′) are the standard angle-averaged redistribution functions that correspond to coherent scattering and to complete frequency redistribution in the atomic rest frame, respectively. The matrices and are given in the Appendix C of Anusha & Nagendra (2011b), according to the results established in Bommier (1997a,b). They are related to the Rayleigh and Hanle phase matrices, which may be further factorized in terms of the irreducible spherical tensors as first introduced in Landi Degl’Innocenti (1984) (see also Landi Degl’Innocenti & Landolfi 2004).

The Hanle phase matrix is given by (2)where Ω and Ω′ denote the propagation directions of the outgoing and ingoing photons, respectively, and x and x′ are their frequencies measured with respect to the line central frequency in units of Doppler width, and B is the magnetic field vector. The index i = 0,1,2,3 refers to the four Stokes parameters. The matrix is given in Landi Degl’Innocenti & Landolfi (2004).

In the absence of a magnetic field, the Rayleigh phase matrix reduces to (3)The irreducible spherical tensors are well suited to handling the angular dependence of Rayleigh-type scattering phase matrices with K = 0, 1, 2 and −K ≤ Q ≤ K. In these problems, the circular and linear polarization decouple.

Using these factorizations, Anusha & Nagendra (2011a,b) showed that the polarized radiative transfer equation along the ray for the Stokes vector I = (I,Q,U)T, namely (4)may be cast into a set of six formally independent transfer equations for the six-component reduced Stokes vector , where the corresponding six-component reduced source vector is angle-independent. That is (5)As usual, ktot(r,x) = kl(r)φ(x) + kc(r), where kl and kc denote, respectively, the line-integrated and the continuum absorption coefficients, and φ(x) is the normalized line absorption profile. The coordinate s is the path length along the ray.

The reason for this remarkable property is that the source vector in Eq. (4) may be written in the irreducible spherical tensors basis as (6)Here the irreducible components of the total source vector are independent of the propagation direction. This is a direct consequence of the factorization of the scattering phase matrix.

As already mentioned in Frisch (2007), the reduced intensity components defined as the solutions of the reduced radiative transfer equations are different from the tensors , defined in Landi Degl’Innocenti & Landolfi (2004), whose angle-averages appear in the statistical equilibrium equations for the element of the atomic density matrix.

The components of the total source function are simply related to the reduced components of the line and continuum source functions and by (7)The reduced line source function is given in Eqs. (18)–(19) of Anusha & Nagendra (2011b) in terms of the reduced six-component radiation vector ℐ, so we may write it as (8)where the 6 × 6 matrix is given by (9)In Eq. (8), the second term of the right-hand side denotes the reduced components of the primary source of line photons, and is a generic notation for a (6 × 6) angle-averaged frequency redistribution matrix. Its analytical expressions, which vary according to different frequency domains, can be derived from the formulae given in Appendices B and C of Anusha & Nagendra (2011b).

The reduced continuum source function may be cast into a form similar to Eq. (8), replacing by the primary source of photons in the continuum and the frequency redistribution matrix by where is the identity matrix.

The reduced source function and intensity vector ℐ are complex quantities, to deal with real quantities Frisch (2007) introduced the following notations: (10)She showed that obeys a radiative transfer equation analogous to Eq. (5), where the corresponding source term can be derived from expressions analogous to Eq. (8) when the matrix is replaced by the matrix , as given in Appendix D of Anusha & Nagendra (2011b). The derivation of the transfer equation for the real reduced components of the radiation field relies on the symmetry property .

The reduced formalism for polarized transfer that we have summarized here is well suited to an efficient numerical solution, such as Polarized Accelerated Lambda Iterations (PALI) or Jacobi methods, which may be applied in 3D media. Once the reduced components of the radiation field have been obtained, one can easily compute the Stokes parameters of the radiation field using the relations given in the Appendix B of Frisch (2007).

2.2. Boundary conditions for the reduced radiative transfer equation

In the papers quoted above, the radiative transfer problem was solved for self-emitting media in the absence of incident radiation. In Milić (2013), the same formalism was used but for computing line-scattering polarization emerging from a circumstellar disk illuminated by its central star. To self-consistently solve the reduced transfer equation, one needs to cast the boundary conditions, i.e. the known incident radiation field, into the same reduced basis. However, there is no unique way of deriving the reduced components from a known radiation field described in terms of the Stokes parameters (only the reversed operation is well defined). The reason is probably that the reduced formalism is well suited to describe the angular dependence of a radiation field that arises from scattering processes, or from unpolarized and isotropic primary creation of photons, but it is not clear how to apply it to any arbitrary radiation. For problems where the medium is illuminated by an arbitrary radiation field, which can be polarized or unpolarized, we have to find a way of writing the boundary conditions for Eq. (5).

In the following, we show that a consistent way of solving the problem is to write the radiation field as the sum of a transmitted radiation and a diffuse part. The diffuse radiation obeys a transfer equation like Eq. (5) with an additional primary creation term that arises from the first scattering of the incident radiation in the medium. By definition, the diffuse radiation transfer problem has no incident radiation. The additional primary source term is easily cast into reduced components because it originates in a scattering process. The diffuse radiation may then be computed in the reduced formalism presented above and its Stokes parameters can then be derived. The full radiation is finally obtained in the Stokes formalism by adding the Stokes components of the transmitted radiation and the Stokes components of the diffuse radiation.

3. Transmitted and diffuse radiation fields

We now deal with the polarized radiative transfer equation in the Stokes formalism, given by Eq. (4). Here we want to solve it the case where the medium is illuminated by an incident radiation field denoted by Iinc(r0,x) with known Stokes parameters , where i = 0,1,2,3. As in Chandrasekhar (1960) and Ivanov et al. (1997), we write the solution of Eq. (4) as the sum of two terms: (11)where τΩ,x(r) is the line optical depth from the boundary at the location r along the direction Ω at frequency x. The first term in the right-hand side of Eq. (11) is the transmitted radiation at depth r and the second term is the diffuse radiation field. Substituting Eq. (11) in Eq. (4), we easily derive the radiative transfer equation for Id(r,x), namely (12)The source term is the same as in Eq. (4); it may be written in terms of the diffuse and transmitted radiation by substituting Eq. (11) into the scattering integral. The scattering integral is now the sum of two terms, of which both may be expanded on the irreducible tensors . Thus, we can repeat for the diffuse radiation field the same procedure as described in the previous section and cast the radiative transfer Eq. (12) into a reduced form similar to Eq. (5), where the reduced source vector is given in Eq. (8).

Here, we write the line source vector on its real form in terms of the diffuse field . Namely, (13)where is the six-component additional source term, related to the real and imaginary parts of the complex quantities, (14)The coefficients WK are constants that depend on the quantum numbers of the line atomic levels. Analytical expressions of as a function of the Stokes parameters of the incident radiation are given in Appendix A. The continuum source function in the reduced form has an expression similar to Eq. (13) but for coherent Rayleigh scattering.

In non-magnetic media, the redistribution matrix in Eq. (14) reduces to (15)with (16)In the presence of Hanle effect, we recall that the redistribution matrix has different expressions, according to different frequency domains (see Appendix B in Anusha & Nagendra 2011b). For the lines formed in dilute media (such as prominences) where collisions are negligible, we have: (17)in the frequency domain where and . In the complementary domain, (18)The cutoff frequency νc between the line core and line wings depends on the value of the Voigt parameter a. It is usually on the order of a few Doppler widths.

The diffuse radiation field is zero at the boundaries of the medium. We can then solve its radiative transfer equation in reduced form and derive its Stokes parameters from the relations given in Appendix B of Frisch (2007). When the Stokes parameters of the diffuse field are computed, we obtain the Stokes parameters of the total radiation field through Eq. (11).

4. Numerical tests

The numerical solution of the reduced transfer problem for the diffuse radiation may be performed with the Jacobi iterative method introduced in Anusha & Nagendra (2011a; see also Milić 2013). The additional source terms and have to be computed first on long characteristics from the boundaries to each point of the spatial grid. We note that this significantly increases the computational work because of the use of the long characteristics, moreover integration must be performed on dense grids if high precision is to be achieved. For the numerical tests presented below, we use a 100-point Gaussian quadrature for integration over θ, while integration over azimuth is not needed because of the symmetry of the incident radiation. For this particular test case, the use of quadrature with 20 Gaussian points results in differences smaller than 1%.

In all the tests we have used a slightly modified second order short-characteristics (Olson & Kunasz 1987) formal solution. A 2D-Forth and Back Implicit Lambda Iterative procedure (FBILI) is used to obtain a self-consistent solution of the unpolarized NLTE problem1. We then perform several Λ iterations for the polarized part with the same formal solver to obtain the solution for the reduced source function.

4.1. Unpolarized incident radiation field

In Anusha & Nagendra (2011a,b) it is assumed that one may use as boundary conditions for the reduced radiation field, if the incident radiation is not polarized, the expression ℐinc = (Iinc,0,0,0,0,0)T. In that case one can verify that both this approach and the diffuse radiation field approach lead to identical expressions for the Stokes parameters of the total radiation field.

However, the diffuse field approach allows us to consider the anisotropy of the incident radiation in the computation of the polarized source function in a more accurate way, even in the case when the incident radiation is unpolarized, because of the following: If one was to solve the radiative transfer problem, one would have to specify ℐinc on a coarse angular grid. In the case of an anisotropic illumination, which also encounters abrupt changes (such as in the following examples), this would introduce errors in the computation of the reduced source vector. The accuracy is improved if one uses the averaging procedure described in Gouttebroze (2005) and subsequently in Milić (2013). However, the averaging process is constructed to accurately reproduce the component (identical to the scalar mean intensity J), while the computation of the other components is still not accurate. For externally illuminated objects of low and moderate optical thickness, the incident radiation is the dominant part of the reduced source function. This part is computed with a much higher precision with our approach because the averaging process is avoided, while the diffuse part of the computation can still be done with a coarser angular grid.

4.2. Polarized incident radiation field

The formalism developed in Anusha & Nagendra (2011a) does not deal with the boundary conditions in cases where the incident radiation is polarized. In that case, we have to use the diffuse field approach. Here we consider slabs that are effectively 1D (computations are performed with a 2D code), placed horizontally or vertically above the solar surface and illuminated by a limb darkened continuum radiation given by (similar example as in: Landi Degl’Innocenti & Landolfi 2004) (19)We consider a case with the following parameters: there is no continuum opacity, the line absorption profile is given by a Doppler function, the line is formed by only scattering processes and W0 = W2 = 1 is assumed. The line-integrated optical depth in both cases is equal to unity. This creates some radiative transfer effects while still maintaining a large impact of the incident radiation. In both cases, the slab is standing 50   000 km above the solar surface. Note that it introduces an additional geometrical anisotropy since the incident radiation is within a limited solid angle between the solar surface and the slab. As a coarse computational grid, we use a grid with 12 × 12 angles where integrations over θ and ϕ are done with a Gaussian and trapezoid quadrature, respectively.

thumbnail Fig. 1

Q/I line profiles emergent from a vertical slab illuminated by a polarized incident radiation. The different curves correspond to different maximum degrees of polarization of the incident radiation field.

Open with DEXTER

Figure 1 shows the emergent Q/I line profiles in the direction θ = π/2, ϕ = 0 from a vertical slab (simple model of a solar prominence), which is illuminated by a polarized radiation with the following angular dependence: (20)We vary Pmax between 0% and 20%. It is clearly seen that the degree of polarization of the incident radiation significantly influences the polarization in the line. A similar effect can be seen in Fig. 2 in the case of the horizontal slab (simple model of a solar filament). Note that the emergent radiation is a combination of the transmitted and scattered radiations in the latter case.

thumbnail Fig. 2

Same as Fig. 1, except for a horizontal slab, μ = 0.6.

Open with DEXTER

5. Conclusion

We have demonstrated a consistent way of dealing with boundary conditions for polarized radiative transfer in the framework of the reduced intensity formalism. Our method is based on splitting the radiation field into a transmitted and a diffuse component. We then incorporate the incident radiation (known in the Stokes formalism) in the source function by performing high-precision angular integration and appropriate attenuation via long characteristics. This method is computationally slower but allows for higher precision, even in the unpolarized case and, more importantly, allows accounting for polarized incident radiation.

We have applied the method on simple 1D slabs, which mimic prominences or filaments illuminated by the solar continuum radiation. We have shown that the degree of polarization of the incoming radiation may significantly influence the emergent line polarization from such structures. It is possible that the interpretation of UV lines formed by scattering of polarized continuum radiation in the upper layers of the Sun might require such computational techniques.


Work in preparation by Milić & Atanacković.


I.M. is thankful to Eiffel Scholarship of French Government and to COST MP1104 for supporting his stay in Nice. O.A. is thankful to the University of Nice-Sophia Antipolis for supporting her stay in Nice. This research is also partially funded by the Serbian Ministry of Science and Education through the project 176004, “Stellar Physics”.


Appendix A: Explicit expressions for the real reduced primary source term due to the incident radiation

Here, we give the analytical expressions for the additional source term in the RT equation for the diffuse radiation in real form, .

In the non-magnetic case, In the presence of a magnetic field for lines formed in dilute media, the same magnetic matrix as given in Frisch (2007) can be applied to to account for the Hanle effect.

All Figures

thumbnail Fig. 1

Q/I line profiles emergent from a vertical slab illuminated by a polarized incident radiation. The different curves correspond to different maximum degrees of polarization of the incident radiation field.

Open with DEXTER
In the text
thumbnail Fig. 2

Same as Fig. 1, except for a horizontal slab, μ = 0.6.

Open with DEXTER
In the text