The use of iteration factors in the solution of the NLTE line transfer problem-II. Multilevel atom

The iteration factors method (IFM) developed in Paper I [1] (Atanacković-Vukmanović and Simonneau, 1994) to solve the NLTE line transfer problem for a two-level atom model, is extended here to deal with a multilevel atom case. At the beginning of each iteration step, for each line transition, angle and frequency averaged depth-dependent iteration factors are computed from the formal solution of radiative transfer (RT) equation and used to close the system of the RT equation moments, non-linearly coupled with the statistical equilibrium (SE) equations. Non-linear coupling of the atomic level populations and the corresponding line radiation field intensities is tackled in two ways. One is based on the linearization of the equations with respect to the relevant variables, and the other on the use of the old (known from the previous iteration) level populations in the line-opacity-like terms of the SE equations. In both cases the use of quasi-invariant iteration factors provided very fast and accurate solution. The properties of the proposed procedures are investigated in detail by applying them to the solution of the prototype multilevel RT problem of Avrett and Loeser [2], and compared with the properties of some other methods.