A comprehensive photoionization model of IC 418: physical conditions, ionization structure, and chemical abundances

Abstract

We constructed a self-consistent stellar+nebular photoionization model of IC 418 using the Cloudy code (v23.01) in 1-D spherical geometry. The central star is approximated by a 36,700 K blackbody. Nebular elemental abundances are set to values typical for a carbon-rich, Type II planetary nebula. Under these assumptions the model reproduces the observed emission and yields the expected ionization structure: an inner He++ zone, surrounded by a He+ shell, and an extended H+ region. The electron temperature rises to approximately 9000 K in high-ionization zones and falls to approximately 7000–8000 K at large radii, consistent with photoionization heating balanced by line cooling. The solution is in photoionization equilibrium (heating = cooling) by construction. Heavy-element abundances (Mg, Si, S, etc.) must be reduced relative to solar (approximately 0.5–3 dex) to match collisionally-excited line strengths, implying those elements are locked in grains. We include dust in the model and find carbonaceous grains reproduce the IR continuum; however, to avoid over-heating the grains one needs either larger grain sizes or confinement of dust to the neutral (PDR) zone. In summary, our Cloudy model (blackbody star, spherical gas, representative abundances) self-consistently matches IC 418’s emission, yielding the expected H+, He+, He++ ionization layers, realistic radial temperature gradients, and insights into its dust content and radiative transfer.