A new nonlocal thermodynamical equilibrium radiative transfer method for cool stars Method and numerical implementation

Research output: Contribution to journalArticle


Context. The solution of the nonlocal thermodynamical equilibrium (non-LTE) radiative transfer equation usually relies on stationary iterative methods, which may falsely converge in some cases. Furthermore, these methods are often unable to handle large-scale systems, such as molecular spectra emerging from, for example, cool stellar atmospheres. Aims. Our objective is to develop a new method, which aims to circumvent these problems, using nonstationary numerical techniques and taking advantage of parallel computers. Methods. The technique we develop may be seen as a generalization of the coupled escape probability method. It solves the statistical equilibrium equations in all layers of a discretized model simultaneously. The numerical scheme adopted is based on the generalized minimum residual method. Results. The code has already been applied to the special case of the water spectrum in a red supergiant stellar atmosphere. This demonstrates the fast convergence of this method, and opens the way to a wide variety of astrophysical problems.


External organisations
  • University of Montpellier
  • Joseph Fourier University
Research areas and keywords

Subject classification (UKÄ) – MANDATORY

  • Astronomy, Astrophysics and Cosmology


  • radiative transfer, stars: atmospheres, methods: numerical
Original languageEnglish
Article numberA50
JournalAstronomy & Astrophysics
Publication statusPublished - 2015
Publication categoryResearch

Related projects

Nils Ryde & Julien Lambert


Project: ResearchIndividual research project

View all (2)