Astronomy and Astrophysics – Astronomy
Scientific paper
May 1999
adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=1999apj...516..436t&link_type=abstract
The Astrophysical Journal, Volume 516, Issue 1, pp. 436-450.
Astronomy and Astrophysics
Astronomy
37
Line: Profiles, Methods: Numerical, Polarization, Radiative Transfer, Scattering, Stars: Atmospheres, Sun: Atmosphere
Scientific paper
This paper shows how to generalize to non-LTE polarization transfer some operator splitting methods that were originally developed for solving unpolarized transfer problems. These are the Jacobi-based accelerated Lambda-iteration (ALI) method of Olson, Auer, & Buchler and the iterative schemes based on Gauss-Seidel and successive overrelaxation (SOR) iteration of Trujillo Bueno and Fabiani Bendicho. The theoretical framework chosen for the formulation of polarization transfer problems is the quantum electrodynamics (QED) theory of Landi Degl'Innocenti, which specifies the excitation state of the atoms in terms of the irreducible tensor components of the atomic density matrix. This first paper establishes the grounds of our numerical approach to non-LTE polarization transfer by concentrating on the standard case of scattering line polarization in a gas of two-level atoms, including the Hanle effect due to a weak microturbulent and isotropic magnetic field. We begin demonstrating that the well-known Lambda-iteration method leads to the self-consistent solution of this type of problem if one initializes using the ``exact'' solution corresponding to the unpolarized case. We show then how the above-mentioned splitting methods can be easily derived from this simple Lambda-iteration scheme. We show that our SOR method is 10 times faster than the Jacobi-based ALI method, while our implementation of the Gauss-Seidel method is 4 times faster. These iterative schemes lead to the self-consistent solution independently of the chosen initialization. The convergence rate of these iterative methods is very high; they do not require either the construction or the inversion of any matrix, and the computing time per iteration is similar to that of the Lambda-iteration method.
Manso Sainz Rafael
Trujillo Bueno Javier
No associations
LandOfFree
Iterative Methods for the Non-LTE Transfer of Polarized Radiation: Resonance Line Polarization in One-dimensional Atmospheres does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with Iterative Methods for the Non-LTE Transfer of Polarized Radiation: Resonance Line Polarization in One-dimensional Atmospheres, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Iterative Methods for the Non-LTE Transfer of Polarized Radiation: Resonance Line Polarization in One-dimensional Atmospheres will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-1123938