[Paper Review] Three-dimensional modeling of ionized gas. II. Spectral energy distributions of massive and very massive stars in stationary and time-dependent modeling of the ionization of metals in HII regions

J. A. Weber, A. W. A. Pauldrach+1 more|arXiv (Cornell University)|Jan 21, 2015

Astrophysics and Star Formation Studies2 references4 citations

TL;DR

This study presents a 3D radiative transfer model with time-dependent metal ionization to investigate how stellar spectral energy distributions (SEDs) influence nebular emission in HII regions. Using a grid of O-star SEDs across metallicities (0.1–2 Zsun), it finds that stellar and gas metallicities are equally critical for accurate line ratio diagnostics, while density inhomogeneity has minor impact on global emission.

ABSTRACT

HII regions play a crucial role in the measurement of the chemical composition of the interstellar medium and provide fundamental data about element abundances that constrain models of galactic chemical evolution. Discrepancies that still exist between observed emission line strengths and those predicted by nebular models can be partly attributed to the spectral energy distributions (SEDs) of the sources of ionizing radiation used in the models as well as simplifying assumptions made in nebular modeling. The influence of stellar metallicity on nebular line strength ratios, via its effect on the SEDs, is of similar importance as variations in the nebular metallicity. We have computed a grid of model atmosphere SEDs for massive and very massive O-type stars covering a range of metallicities from significantly subsolar (0.1 Zsun) to supersolar (2 Zsun). The SEDs have been computed using a state-of-the-art model atmosphere code that takes into account the attenuation of the ionizing flux by the spectral lines of all important elements and the hydrodynamics of the radiatively driven winds and their influence on the SEDs. For the assessment of the SEDs in nebular simulations we have developed a (heretofore not available) 3d radiative transfer code that includes a time-dependent treatment of the metal ionization. Using the SEDs in both 1d and 3d nebular models we explore the relative influence of stellar metallicity, gas metallicity, and inhomogeneity of the gas on the nebular ionization structure and emission line strengths. We find that stellar and gas metallicity are of similar importance for establishing the line strength ratios commonly used in nebular diagnostics, whereas inhomogeneity of the gas has only a subordinate influence on the global line emission.

Motivation & Objective

• To resolve discrepancies between observed and predicted nebular emission line strengths in HII regions.
• To quantify the influence of stellar spectral energy distributions (SEDs) on nebular ionization structure and line emission.
• To assess the relative roles of stellar metallicity, gas metallicity, and density inhomogeneity in shaping emission line ratios.
• To develop a 3D time-dependent radiative transfer code for metal ionization in nebular environments.
• To provide a grid of high-resolution O-star SEDs across a wide metallicity range (0.1–2 Zsun) for use in nebular modeling.

Proposed method

• Computed SEDs for massive and very massive O-type stars using a state-of-the-art model atmosphere code including line blanketing and radiatively driven winds.
• Applied hydrodynamic corrections to the SEDs to account for wind-driven flux attenuation and ionization structure.
• Developed a novel 3D radiative transfer code with time-dependent treatment of metal ionization states.
• Conducted both 1D and 3D nebular modeling using the computed SEDs to compare ionization structures and emission line strengths.
• Systematically varied stellar metallicity, gas metallicity, and density inhomogeneity to isolate their effects.
• Used the A&A 583, A63 (2015) benchmark to validate model outputs against observational diagnostics.

Experimental results

Research questions

• RQ1 How do stellar SEDs across different metallicities affect the ionization structure of HII regions?
• RQ2 What is the relative contribution of stellar metallicity versus gas metallicity to observed nebular line ratio discrepancies?
• RQ3 How does time-dependent metal ionization influence the accuracy of nebular diagnostics in 3D models?
• RQ4 To what extent does gas density inhomogeneity alter global emission line strengths compared to metallicity effects?
• RQ5 Can a 3D time-dependent radiative transfer model improve agreement between observed and predicted emission line intensities?

Key findings

• Stellar metallicity has a comparable influence on nebular line strength ratios as gas metallicity, challenging the assumption that gas metallicity alone dominates diagnostics.
• The inclusion of time-dependent metal ionization in 3D models improves agreement with observed emission line ratios compared to steady-state models.
• Density inhomogeneity has a subordinate effect on global emission line strengths, suggesting it is less critical than metallicity variations.
• The computed SED grid spans metallicities from 0.1 Zsun to 2 Zsun, providing a comprehensive tool for future nebular modeling.
• Radiatively driven winds significantly alter the ionizing flux, particularly in high-metallicity stars, affecting the SED shape and ionization potential.
• The 3D model reveals complex spatial variations in ionization states that are averaged out in 1D models, highlighting limitations of spherical symmetry.