Atomic line radiative transfer with MCFOST I. Code description and benchmarking
arXiv:2101.04722 · doi:10.1051/0004-6361/202039697
Abstract
Aims. We present MCFOST-art, a new non-local thermodynamic equilibrium radiative transfer solver for multilevel atomic systems. The code is embedded in the 3D radiative transfer code MCFOST and is compatible with most of the MCFOST modules. The code is versatile and designed to model the close environment of stars in 3D. Methods. The code solves for the statistical equilibrium and radiative transfer equations using the Multilevel Accelerated Lambda Iteration (MALI) method. We tested MCFOST-art on spherically symmetric models of stellar photospheres as well as on a standard model of the solar atmosphere. We computed atomic level populations and outgoing fluxes and compared these values with the results of the TURBOspectrum and RH codes. Calculations including expansion and rotation of the atmosphere were also performed. We tested both the pure local thermodynamic equilibrium and the out-of-equilibrium problems. Results. In all cases, the results from all codes agree within a few percent at all wavelengths and reach the sub-percent level between RH and MCFOST-art. We still note a few marginal discrepancies between MCFOST-art and TURBOspectrum as a result of different treatments of background opacities at some critical wavelength ranges.
References in corpus (19)
- HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
- A grid of MARCS model atmospheres for late-type stars I. Methods and general properties
- Monte Carlo radiative transfer in protoplanetary disks
- Benchmark problems for continuum radiative transfer. High optical depths, anisotropic scattering, and polarisation
- Circumbinary, not transitional: On the spiral arms, cavity, shadows, fast radial flows, streamers and horseshoe in the HD142527 disc
- STiC -- A multi-atom non-LTE PRD inversion code for full-Stokes solar observations
- The AMBRE project: A new synthetic grid of high-resolution FGKM stellar spectra
- Numerical methods for non-LTE line radiative transfer: Performance and convergence characteristics
- Multidimensional models of hydrogen and helium emission line profiles for classical T Tauri Stars: method, tests and examples
- Accretion dynamics in the classical T Tauri star V2129 Oph
- Excitation and charge transfer in low-energy hydrogen atom collisions with neutral atoms: Theory, comparisons, and application to Ca
- Ring formation and dust dynamics in wind-driven protoplanetary discs: global simulations
- Modeling the Halpha line emission around classical T Tauri stars using magnetospheric accretion and disk wind models
- Using 3D Voronoi grids in radiative transfer simulations
- IRIS: A Generic Three-Dimensional Radiative Transfer Code
- Magritte, a modern software library for 3D radiative transfer: I. Non-LTE atomic and molecular line modelling
- 3D radiative transfer: Continuum and line scattering in non-spherical winds from OB stars
- Near optimal angular quadratures for polarised radiative transfer
- A 3D radiative transfer framework: XI. multi-level NLTE
Cited by in corpus (4)
- Combined analysis of stellar and planetary absorption lines via global forward-transit simulations
- 3D simulations of AGB stellar winds -- II. Ray-tracer implementation and impact of radiation on the outflow morphology
- Spectroscopic and interferometric signatures of magnetospheric accretion in young stars
- Tetrahedral grids in Monte Carlo radiative transfer