skip to main content


Title: Carbon monoxide formation and cooling in supernovae
Context. The inclusion of molecular physics is an important piece that tends to be missing from the puzzle when modeling the spectra of supernovae (SNe). Molecules have both a direct impact on the spectra, particularly in the infrared, and an indirect one as a result of their influence on certain physical conditions, such as temperature. Aims. In this paper, we aim to investigate molecular formation and non-local thermodynamic equilibrium (NLTE) cooling, with a particular focus on CO, the most commonly detected molecule in supernovae. We also aim to determine the dependency of supernova chemistry on physical parameters and the relative sensitivity to rate uncertainties. Methods. We implemented a chemical kinetic description of the destruction and formation of molecules into the SN spectral synthesis code SUMO . In addition, selected molecules were coupled into the full NLTE level population framework and, thus, we incorporated molecular NLTE cooling into the temperature equation. We produced a test model of the CO formation in SN 1987A between 150 and 600 days and investigated the sensitivity of the resulting molecular masses to the input parameters. Results. We find that there is a close inter-dependency between the thermal evolution and the amount of CO formed, mainly through an important temperature-sensitive CO destruction process with O + . After a few hundred days, CO completely dominates the cooling of the oxygen-carbon zone of the supernova which, therefore, contributes little optical emission. The uncertainty of the calculated CO mass scales approximately linearly with the typical uncertainty factor for individual rates. We demonstrate how molecular masses can potentially be used to constrain various physical parameters of the supernova.  more » « less
Award ID(s):
1927130
NSF-PAR ID:
10296883
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Astronomy & Astrophysics
Volume:
642
ISSN:
0004-6361
Page Range / eLocation ID:
A135
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We present the spectroscopic and photometric study of five intermediate-luminosity red transients (ILRTs), namely AT 2010dn, AT 2012jc, AT 2013la, AT 2013lb, and AT 2018aes. They share common observational properties and belong to a family of objects similar to the prototypical ILRT SN 2008S. These events have a rise time that is less than 15 days and absolute peak magnitudes of between −11.5 and −14.5 mag. Their pseudo-bolometric light curves peak in the range 0.5–9.0 × 10 40  erg s −1 and their total radiated energies are on the order of (0.3–3) × 10 47 erg. After maximum brightness, the light curves show a monotonic decline or a plateau, resembling those of faint supernovae IIL or IIP, respectively. At late phases, the light curves flatten, roughly following the slope of the 56 Co decay. If the late-time power source is indeed radioactive decay, these transients produce 56 Ni masses on the order of 10 −4 to 10 −3   M ⊙ . The spectral energy distribution of our ILRT sample, extending from the optical to the mid-infrared (MIR) domain, reveals a clear IR excess soon after explosion and non-negligible MIR emission at very late phases. The spectra show prominent H lines in emission with a typical velocity of a few hundred km s −1 , along with Ca II features. In particular, the [Ca  II ] λ 7291,7324 doublet is visible at all times, which is a characteristic feature for this family of transients. The identified progenitor of SN 2008S, which is luminous in archival Spitzer MIR images, suggests an intermediate-mass precursor star embedded in a dusty cocoon. We propose the explosion of a super-asymptotic giant branch star forming an electron-capture supernova as a plausible explanation for these events. 
    more » « less
  2. Supernova LSQ13abf was discovered soon after explosion by the La Silla-QUEST Survey and then followed by the Carnegie Supernova Project II at its optical and near-IR wavelengths. Our analysis indicates that LSQ13abf was discovered within two days of explosion and its first ≈10 days of evolution reveal a B -band light curve with an abrupt drop in luminosity. Contemporaneously, the V -band light curve exhibits a rise towards a first peak and the r - and i -band light curves show no early peak. The early light-curve evolution of LSQ13abf is reminiscent of the post-explosion cooling phase observed in the Type Ib SN 2008D, and the similarity between the two objects extends over weeks. Spectroscopically, LSQ13abf also resembles SN 2008D, with P Cygni He  I features that strengthen over several weeks. Spectral energy distributions are constructed from the broad-bandphotometry, a UVOIR light curve is constructed by fitting black-body (BB) functions, and the underlying BB-temperature and BB-radius profiles are estimated. Explosion parameters are estimated by simultaneously fitting an Arnett model to the UVOIR light curve and the velocity evolution derived from spectral features, and an in addition to a post-shock breakout cooling model to the first two epochs of the bolometric evolution. This combined model suggests an explosion energy of 1.27 ± 0.23 × 10 51 ergs, in addition to a relatively high ejecta mass of 5.94 ± 1.10 M ⊙ , a 56 Ni mass of 0.16 ± 0.02 M ⊙ , and a progenitor-star radius of 28.0 ± 7.5 R ⊙ . The ejecta mass suggests the origins of LSQ13abf lie with a > 25  M ⊙ zero-age-main-sequence mass progenitor and its estimated radius is three times larger compared to the result obtained from the same analysis applied to observations of SN 2008D, and nine times larger compared to SN 1999ex. Alternatively, a comparison of hydrodynamical simulations of ≳20−25 M ⊙ zero-age-main-sequence progenitors that evolve to pre-supernova envelope masses of ≲10 M ⊙ and extended (∼100 R ⊙ ) envelopes also broadly match the observations of LSQ13abf. 
    more » « less
  3. We present optical photometry and spectroscopy of the Type II supernova ASASSN-14jb, together with Very Large Telescope (VLT) Multi Unit Spectroscopic Explorer (MUSE) integral field observations of its host galaxy and a nebular-phase spectrum. This supernova, in the nearby galaxy ESO 467-G051 ( z  = 0.006), was discovered and followed-up by the all-sky automated survey for supernovae (ASAS-SN). We obtained well-sampled las cumbres network (LCOGTN) B V g r i and Swift w 2 m 1 w 1 u b v optical, near-UV/optical light curves, and several optical spectra in the early photospheric phases. The transient ASASSN-14jb exploded ∼2 kpc above the star-forming disk of ESO 467-G051, an edge-on disk galaxy. The large projected distance from the disk of the supernova position and the non-detection of any H II region in a 1.4 kpc radius in projection are in conflict with the standard environment of core-collapse supernova progenitors and suggests the possible scenario that the progenitor received a kick in a binary interaction. We present analysis of the optical light curves and spectra, from which we derived a distance of 25 ± 2 Mpc using state-of-the-art empirical methods for Type II SNe, physical properties of the SN explosion ( 56 Ni mass, explosion energy, and ejected mass), and properties of the progenitor; namely the progenitor radius, mass, and metallicity. Our analysis yields a 56 Ni mass of 0.0210  ±  0.0025  M ⊙ , an explosion energy of ≈0.25 × 10 51 ergs, and an ejected mass of ≈6  M ⊙ . We also constrained the progenitor radius to be R *  = 580  ±  28  R ⊙ which seems to be consistent with the sub-Solar metallicity of 0.3  ±  0.1  Z ⊙ derived from the supernova Fe II λ 5018 line. The nebular spectrum constrains strongly the progenitor mass to be in the range 10–12 M ⊙ . From the Spitzer data archive we detect ASASSN-14jb ≈330 days past explosion and we derived a total dust mass of 10 −4   M ⊙ from the 3.6 μ m and 4.5 μ m photometry. Using the F U V , N U V , B V g r i , K s , 3.6 μ m, and 4.5 μ m total magnitudes for the host galaxy, we fit stellar population synthesis models, which give an estimate of M *  ≈ 1 × 10 9   M ⊙ , an age of 3.2 Gyr, and a SFR ≈0.07  M ⊙ yr −1 . We also discuss the low oxygen abundance of the host galaxy derived from the MUSE data, having an average of 12 + log(O/H) = 8.27 +0.16 −0.20 using the O 3 N 2 diagnostic with strong line methods. We compared it with the supernova spectra, which is also consistent with a sub-Solar metallicity progenitor. Following recent observations of extraplanar H II regions in nearby edge-on galaxies, we derived the metallicity offset from the disk, being positive, but consistent with zero at 2 σ , suggesting enrichment from disk outflows. We finally discuss the possible scenarios for the unusual environment for ASASSN-14jb and conclude that either the in-situ star formation or runaway scenario would imply a low-mass progenitor, agreeing with our estimate from the supernova nebular spectrum. Regardless of the true origin of ASASSN-14jb, we show that the detailed study of the environment roughly agree with the stronger constraints from the observation of the transient. 
    more » « less
  4. null (Ed.)
    Context. The excitation of the filamentary gas structures surrounding giant elliptical galaxies at the center of cool-core clusters, also known as brightest cluster galaxies (BCGs), is key to our understanding of active galactic nucleus (AGN) feedback, and of the impact of environmental and local effects on star formation. Aims. We investigate the contribution of thermal radiation from the cooling flow surrounding BCGs to the excitation of the filaments. We explore the effects of small levels of extra heating (turbulence), and of metallicity, on the optical and infrared lines. Methods. Using the C LOUDY code, we modeled the photoionization and photodissociation of a slab of gas of optical depth A V  ≤ 30 mag at constant pressure in order to calculate self-consistently all of the gas phases, from ionized gas to molecular gas. The ionizing source is the extreme ultraviolet (EUV) and soft X-ray radiation emitted by the cooling gas. We tested these models comparing their predictions to the rich multi-wavelength observations from optical to submillimeter, now achieved in cool core clusters. Results. Such models of self-irradiated clouds, when reaching sufficiently large A V , lead to a cloud structure with ionized, atomic, and molecular gas phases. These models reproduce most of the multi-wavelength spectra observed in the nebulae surrounding the BCGs, not only the low-ionization nuclear emission region like optical diagnostics, [O  III ] λ 5007 Å/H β , [N  II ] λ 6583 Å/H α , and ([S  II ] λ 6716 Å+[S  II ] λ 6731 Å)/H α , but also the infrared emission lines from the atomic gas. [O  I ] λ 6300 Å/H α , instead, is overestimated across the full parameter space, except for very low A V . The modeled ro-vibrational H 2 lines also match observations, which indicates that near- and mid-infrared H 2 lines are mostly excited by collisions between H 2 molecules and secondary electrons produced naturally inside the cloud by the interaction between the X-rays and the cold gas in the filament. However, there is still some tension between ionized and molecular line tracers (i.e., CO), which requires optimization of the cloud structure and the density of the molecular zone. The limited range of parameters over which predictions match observations allows us to constrain, in spite of degeneracies in the parameter space, the intensity of X-ray radiation bathing filaments, as well as some of their physical properties like A V or the level of turbulent heating rate. Conclusions. The reprocessing of the EUV and X-ray radiation from the plasma cooling is an important powering source of line emission from filaments surrounding BCGs. C LOUDY self-irradiated X-ray excitation models coupled with a small level of turbulent heating manage to simultaneously reproduce a large number of optical-to-infrared line ratios when all the gas phases (from ionized to molecular) are modeled self-consistently. Releasing some of the simplifications of our model, like the constant pressure, or adding the radiation fields from the AGN and stars, as well as a combination of matter- and radiation-bounded cloud distribution, should improve the predictions of line emission from the different gas phases. 
    more » « less
  5. null (Ed.)
    ABSTRACT Carbon-enhanced metal-poor (CEMP) stars are the living fossils holding records of chemical enrichment from early generations of stars. In this work, we perform a set of numerical simulations of the enrichment from a supernova (SN) of a first generation of metal-free (Pop III) star and the gravitational collapse of the enriched cloud, considering all relevant cooling/heating processes and chemical reactions as well as the growth of dust grains. We adopt faint SN models for the first time with progenitor masses MPopIII = 13–$80 \ {\rm M_{\bigodot }}$, which yield C-enhanced abundance patterns ([C/Fe] = 4.57–4.75) through mixing and fallback of innermost layers of the ejecta. This model also considers the formation and destruction of dust grains. We find that the metals ejected by the SN can be partly re-accreted by the same dark matter minihalo, and carbon abundance of the enriched cloud A(C) = 3.80–5.06 is lower than the abundance range of observed CEMP stars (A(C) ≳ 6) because the mass of the metals ejected by faint SNe is smaller than normal core-collapse SNe due to extensive fallback. We also find that cloud fragmentation is induced by gas cooling from carbonaceous grains for $M_{\rm Pop III}= 13 \ {\rm M_{\bigodot }}$ even with the lowest iron abundance [Fe/H] ∼ −9. This leads to the formation of low-mass stars, and these ‘giga metal-poor’ stars can survive until the present-day Universe and may be found by future observations. 
    more » « less