A 3D radiative transfer framework: XI. multi-level NLTE

Abstract

Multi-level non-local thermodynamic equilibrium (NLTE) radiation transfer calculations have become standard throughout the stellar atmospheres community and are applied to all types of stars as well as dynamical systems such as novae and supernovae. Even today spherically symmetric 1D calculations with full physics are computationally intensive. We show that full NLTE calculations can be done with fully 3 dimensional (3D) radiative transfer. With modern computational techniques and current massive parallel computational resources, full detailed solution of the multi-level NLTE problem coupled to the solution of the radiative transfer scattering problem can be solved without sacrificing the micro physics description. We extend the use of a rate operator developed to solve the coupled NLTE problem in spherically symmetric 1D systems. In order to spread memory among processors we have implemented the NLTE/3D module with a hierarchical domain decomposition method that distributes the NLTE levels, radiative rates, and rate operator data over a group of processes so that each process only holds the data for a fraction of the voxels. Each process in a group holds all the relevant data to participate in the solution of the 3DRT problem so that the 3DRT solution is parallelized within a domain decomposition group. We solve a spherically symmetric system in 3D spherical coordinates in order to directly compare our well-tested 1D code to the 3D case. We compare three levels of tests: a) a simple H+He test calculation, b) H+He+CNO+Mg, c) H+He+Fe. The last test is computationally large and shows that realistic astrophysical problems are solvable now, but they do require significant computational resources. With presently available computational resources it is possible to solve the full 3D multi-level problem with the same detailed micro-physics as included in 1D modeling.

20 pages, 14 figures, A&A, in press

References in corpus (14)

• NLTE line formation of Fe for late-type stars. I. Standard stars with 1D and <3D> model atmospheres
• Analysis of the Flux and Polarization Spectra of the Type Ia Supernova SN 2001el: Exploring the Geometry of the High-velocity Ejecta
• Could There Be A Hole In Type Ia Supernovae?
• A 3D radiative transfer framework: I. non-local operator splitting and continuum scattering problems
• 3D LTE spectral line formation with scattering in red giant stars
• A 3D radiative transfer framework: II. line transfer problems
• A 3D radiative transfer framework: VI. PHOENIX/3D example applications
• Spectral Consequences of Deviation from Spherical Composition Symmetry in Type Ia Supernovae
• A 3D radiative transfer framework: V. Homologous Flows
• A 3D radiative transfer framework: IV. spherical & cylindrical coordinate systems
• A 3D radiative transfer framework: VII. Arbitrary velocity fields in the Eulerian frame
• A 3D radiative transfer framework: III. periodic boundary conditions
• A 3D radiative transfer framework IX. Time dependence
• A 3D radiative transfer framework: X. Arbitrary Velocity Fields in the Co-moving Frame

Cited by in corpus (9)

• Non-LTE line formation of Fe in late-type stars - III. 3D non-LTE analysis of metal-poor stars
• Solar science with the Atacama Large Millimeter/submillimeter Array - A new view of our Sun
• The Broad-Lined Type Ic SN 2012ap and the Nature of Relativistic Supernovae Lacking a Gamma-ray Burst Detection
• Supernova 1987A: neutrino-driven explosions in three dimensions and light curves
• Coupling hydrodynamics with comoving frame radiative transfer: II. Stellar wind stratification in the high-mass X-ray binary Vela X-1
• The stellar and wind parameters of six prototypical HMXBs and their evolutionary status
• MESAS: Measuring the Emission of Stellar Atmospheres at Submm/mm wavelengths
• Atomic line radiative transfer with MCFOST I. Code description and benchmarking
• Modeling the H$\alpha$ Emission Surrounding Spica using the Lyman Continuum from a Gravity-darkened Central Star