Modeling of the nebular-phase spectral evolution of stripped-envelope supernovae

Abstract

We present an extended grid of multi-epoch 1D nonlocal thermodynamic equilibrium radiative transfer calculations for nebular-phase Type Ibc supernovae (SNe) from He-star explosions. Compared to our previous work, which was focused on a post-explosion epoch of 200 days, here we study the spectral evolution from 100 to about 450 days. We also augment the model set with progenitors that evolved without wind mass loss. Models with the same final, pre-SN mass have similar yields and produce essentially the same emergent spectra. Hence, the uncertain progenitor mass loss history compromises the inference of the initial, main sequence mass. This shortcoming does not affect Type IIb SNe in which mass-loss has left a small residual H-rich envelope in the progenitor star at core collapse and, hence, an intact He core. However, our 1D models with a different pre-SN mass tend to yield widely different spectra, as seen through variations in the strong emission lines due to [N II] λλ 6548, 6583, [O I] λλ 6300, 6364, [Ca II] λλ 7291, 7323, [Ni II] λ 7378, and the forest of Fe II lines below 5500 Å. At the lower mass end, the ejecta are He-rich, and at 100 days, they cool through He I, N II, Ca II, and Fe II lines, with N II and Fe II dominating at 450 days. These models, associated with He giants, stand in conflict to observed SNe Ib, which typically lack strong N II emission. Instead, they may lead to SNe Ibn or, because of additional stripping by a companion star, ultra-stripped SNe Ic. In contrast, for higher pre-SN masses, the ejecta are progressively He poor and cool at 100 days through O I, Ca II, and Fe II lines, with O I and Ca II dominating at 450 days. Non-uniform, aspherical, large-scale mixing is more likely to determine the SN type at intermediate pre-SN masses, rather than any compositional differences. Variations in clumping and mixing, as well as departures from spherical symmetry would increase the spectral diversity, but also introduce additional degeneracies. More robust predictions from spectral modeling thus require that careful attention be paid to the initial conditions by incorporating the salient features of physically consistent 3D explosion models.