in a Plane-Parallel Medium

The non-LTE source function in a multiple-scattering/ absorbing plane-parallel atmosphere is generally determined by the integral equation.
_{v} being the total optical depth), t the corresponding optical depth for the absorbers, S_{0}(τ) the initial source function, and the kernel function K_{1} gives the probability for a photon to traverse unscattered and unabsorbed between τ' and τ and be (isotropically) scattered at the latter level. In case of a continuous (i.e. frequency-independent) scattering/absorption cross section, K_{1} is simply determined by the usual exponential absorption law (i.e. for a plane-parallel atmosphere by the first exponential integral E_{1}), but for the case of a spectral line it is given by
^{-x2}, or in case of combined Doppler and natural broadening Φ(x)=H(a,x) (the Voigt function)).

This form for K_{1} implies that both the emission and the absorption profile are given by Φ(x) (complete frequency redistribution) which is furthermore assumed to be independent of τ here (the latter implies for instance that the atmosphere is isothermal). Both of these assumptions are usually reasonably good approximations if one is merely interested in the line-integrated intensities and source functions rather than the details of the spectral profile.
The method of choice here is an iterative solution of Eq.(1). First of all, Eq.(1) is discretized by assuming the source function to be constant within a given sub-interval j and integrating over τ'. This yields
_{2} the second exponential integral, and the optical depth of the absorbers assumed to change proportionally to the optical depth of the scatterers in each interval, i.e.
_{j}) (which also appears on the right hand side of the equation (for k=j)). This yields
_{j}), which can be solved by iteration, using the given initial source function S_{0}(τ_{j}) as starting values.

The source functions are of course not directly observable, but only the intensities, which for 'up-looking' directions (looking into directions of decreasing τ) under an angle η to the normal are here determined by (using μ=cos(η))

Analogously one obtains for the 'down-looking' intensities at level τ_{j}

On the page Radiative Transfer Code for Multiple Scattering in a Plane Parallel Atmosphere I have translated this formalism into a corresponding C/C++ program (this can also be used for the case of continuum- (monochromatic scattering) as the corresponding kernel functions are simply obtained from the above by leaving out the frequency integration and setting Φ(x)=1 (resulting thus in straight exponential functions and exponential integrals))

(note that contrary to the above formulation, the program for the intensity uses the line-of-sight optical depth variable rather than the vertical optical depth; this makes it suitable for non-plane-parallel geometries as well).

(1) S(τ) = S_{0}(τ) + _{0}∫^{τv} dτ' S(τ') ^{.} K_{1}(|τ'-τ|, |t'(τ')-t(τ)|)

(2) K_{1}(τ,t)= 1/2/√π^{.}_{-∞}∫^{∞}dx Φ^{2}(x)^{.}E_{1}(τ^{.}Φ(x)+t) ,

This form for K

_{N}

(3) S(τ_{j}) = S_{0}(τ_{j}) + Σ S(τ_{k}) ^{.}[ | K_{2}(|τ_{k+1}-τ_{j}| , |t_{k+1}-t_{j}| , α_{k}) - K_{2}(|τ_{k}-τ_{j}|, |t_{k}-t_{j}| , α_{k}) | ] ,

^{k=1}

(4) K_{2}(τ,t,α) = _{0}∫^{τ} dτ' K_{1}(τ',t'(τ')) = 1/2/√π^{.}_{-∞}∫^{∞}dx Φ^{2}(x)/(Φ(x)+α) ^{.}E_{2}(τ^{.}Φ(x)+t) .

(5) t'(τ') = t_{k} + α_{k}^{.}(τ'- τ_{k}) ; τ_{k}≤τ'≤τ_{k+1}

_{N}

(6) S(τ_{j}) = [ S_{0}(τ_{j}) + Σ S(τ_{k}) ^{.}dK_{2}(j,k) ] / [1-dK_{2}(j,j) ] ,

^{k=1(≠j)}

(7) dK_{2}(j,k) = | K_{2}(|τ_{k+1}-τ_{j}| , |t_{k+1}-t_{j}| , α_{k}) - K_{2}(|τ_{k}-τ_{j}|, |t_{k}-t_{j}| , α_{k}) | .

The source functions are of course not directly observable, but only the intensities, which for 'up-looking' directions (looking into directions of decreasing τ) under an angle η to the normal are here determined by (using μ=cos(η))

(8) I_{up}(τ_{j},μ) = 1/μ^{.} _{0}∫^{τj} dτ' S(τ') ^{.} T( (τ_{j}-τ')/μ , (t_{j}-t')/μ ) )

(9) T(τ,t) = 1/√π^{.}_{-∞}∫^{∞}dx Φ(x)^{.}e^{-(τ.Φ(x)+t)} :

_{j-1}

(10) I_{up}(τ_{j},μ) = Σ S(τ_{k}) ^{.}[ Z(τ_{j}-τ_{k})/μ , (t_{j}-t_{k})/μ , α_{k}) - Z(τ_{j}-τ_{k+1})/μ , (t_{j}-t_{k+1})/μ , α_{k}) ] ,

^{k=1}

(11) Z(τ,t,α) = _{0}∫^{τ} dτ' T(τ',t'(τ')) = 1/√π^{.}_{-∞}∫^{∞}dx Φ(x)/(Φ(x)+α)^{.}[e^{-(t-ατ)} -e^{-(τ.Φ(x)+t)}] .

Analogously one obtains for the 'down-looking' intensities at level τ

_{N}

(12) I_{down}(τ_{j},μ) = Σ S(τ_{k}) ^{.}[ Z(τ_{k+1}-τ_{j})/μ , (t_{k+1}-t_{j})/μ , α_{k}) - Z(τ_{k}-τ_{j})/μ , (t_{k}-t_{j})/μ , α_{k}) ] .

^{k=j}

On the page Radiative Transfer Code for Multiple Scattering in a Plane Parallel Atmosphere I have translated this formalism into a corresponding C/C++ program (this can also be used for the case of continuum- (monochromatic scattering) as the corresponding kernel functions are simply obtained from the above by leaving out the frequency integration and setting Φ(x)=1 (resulting thus in straight exponential functions and exponential integrals))

(note that contrary to the above formulation, the program for the intensity uses the line-of-sight optical depth variable rather than the vertical optical depth; this makes it suitable for non-plane-parallel geometries as well).

Thomas Smid (M.Sc. Physics, Ph.D. Astronomy)

email: