# Near-infrared Spectroscopy of Supernova 2017eaw in 2017: Carbon Monoxide and Dust Formation in a Type II-P Supernova

J. Rho, T. Geballe, D. Banerjee, L. Dessart, A. Evans, V. Joshi

► To cite this version:

J. Rho, T. Geballe, D. Banerjee, L. Dessart, A. Evans, et al. Near-infrared Spectroscopy of Supernova 2017eaw in 2017: Carbon Monoxide and Dust Formation in a Type II-P Supernova. The Astrophysical Journal Letters, 2018, 864 (1), pp.L20. 10.3847/2041-8213/aad77f. hal-02351311

HAL Id: hal-02351311

https://hal.science/hal-02351311v1

Submitted on 30 Aug 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Near-infrared Spectroscopy of Supernova 2017eaw in 2017: Carbon Monoxide and Dust Formation in a Type II-P Supernova

J. Rho¹,², T. R. Geballe³, D. P. K. Banerjee⁴, L. Dessart⁵, A. Evans⁶, and V. Joshi⁴

¹SETI Institute, 189 N. Bernardo Avenue, Mountain View, CA 94043, USA; jrho@seti.org

²Stratospheric Observatory for Infrared Astronomy, NASA Ames Research Center, MS 211-3, Moffett Field, CA 94035, USA
³Gemini Observatory, 670 North Aohoku Place, Hilo, HI 96720, USA

⁴Physical Research Laboratory, Navrangpura, Ahmedabad, Gujarat 380009, India

⁵Unidad Mixta Internacional Franco-Chilena de Astronomia (CNRS UMI 3386), Departamento de Astronomia, Universidad de Chile, Camino El Observatorio 1515, Las Condes, Santiago, Chile

⁶Astrophysics Group, Keele University, Keele, Staffordshire, ST5 5BG, UK

Received 2018 June 21; revised 2018 July 26; accepted 2018 July 31; published 2018 August 30

## Abstract

The origin of dust in the early universe has been the subject of considerable debate. Core-collapse supernovae (ccSNe), which occur several million years after their massive progenitors form, could be a major source of that dust, as in the local universe several ccSNe have been observed to be copious dust producers. Here we report nine near-infrared (0.8–2.5 µm) spectra, obtained with GNIRS on Gemini North, of the Type II-P supernova (SN) 2017eaw in NGC 6946, spanning the time interval 22–205 days after discovery. The spectra show the onset of CO formation and continuum emission at wavelengths greater than 2.1 µm from newly formed hot dust, in addition to numerous lines of hydrogen and metals, which reveal the change in ionization as the density of much of the ejecta decreases. The observed CO masses estimated from a local thermodynamic equilibrium (LTE) model are typically $10^{-4} M_{\odot}$ during days 124–205, but could be an order of magnitude larger if non-LTE conditions are present in the emitting region. The timing of the appearance of CO is remarkably consistent with the chemically controlled dust models of Sarangi & Cherchneff.

**Key words:** dust, extinction – ISM: molecules – molecular processes – supernovae: individual (SN 2017eaw)

## 1. Introduction

Whether core-collapse supernovae (ccSNe) are a significant source of dust in the universe is a long-standing question. The large quantities of dust observed in high-redshift galaxies (e.g., Isaak et al. 2002; Watson et al. 2015) raise a fundamental question as to the cause of its early presence, because stars, which are thought to produce most of the interstellar dust when on the asymptotic giant branch (AGB) in the modern Milky Way and nearby galaxies, could not have evolved to the dust-producing stage in very high-redshift galaxies. In contrast, supernovae (SNe) from massive stars occur just millions of years after they form. Michalowski (2015) concluded that the early appearance of dust could be explained if ccSNe do not subsequently destroy the dust that they create. However, dust forms and grows in the interstellar medium (Spitzer 1978; Zhukovska et al. 2018), and the dust in very high-z systems might be produced primarily there, seeded by heavy elements and/or surviving dust from SNe (Draine 2009). Depending on the progenitor mass, dust masses of 0.1–1.0 $M_{\odot}$ per ccSN (Todini & Ferrara 2001; Nozawa et al. 2003) can account for the amount of dust observed in very high-redshift galaxies. Evidence that supernova-produced dust is the primary source of dust in the early universe is still scarce (see discussion in Gall et al. 2011; Rho et al. 2018, and references therein) and requires further study.

Spitzer mapping of the young supernova remnant (YSNR) Cas A revealed a remarkable similarity between the dust and CO distributions, confirming that dust forms and is preserved in at least some ccSNe ejecta (Rho et al. 2008, 2012). Four YSNRs, including SN 1987A (Matsuura et al. 2011, 2015), the Crab Nebula (Gomez et al. 2012; Temim et al. 2012; Temim & Dwek 2013), Cas A (Rho et al. 2008; De Looze et al. 2017),

and G54.1 + 0.3 (Temim et al. 2017; Rho et al. 2018), have current dust masses of 0.1–0.9 $M_{\odot}$ that are consistent with theoretical models (Todini & Ferrara 2001; Nozawa et al. 2003). The results suggest that some ccSNe are dust factories, and thus that they could be very important at high z, when other significant dust production processes may not exist. However, the fraction of ccSNe that produce this much dust is not known. Observing dust in additional local ccSNe and following its production and evolution in individual SNe over an extended period of time will help provide the answer.

Current understanding of dust chemistry and composition in ccSNe is limited. Some models (e.g., Nozawa et al. 2003) suggest that dust forms 350–900 days after a ccSN explosion, with carbon dust being one of the first condensates. Those models assume that little carbon is locked up in CO and almost all of it is available for dust formation. CO is one of the most powerful coolants in the ejecta of type II SNe and, once it forms at several thousand K in some portions of the ejecta, is believed to be responsible in large part for cooling those ejecta to temperatures $\lesssim 1700$ K, at which dust can form. Dust evolution models indicate that dust production depends significantly on the C/O ratio in the gas (Morgan & Edmunds 2003; Dwek 1998), as well as on the rate of destruction of CO, by impact both with the energetic electrons produced by the radioactive decay of $^{56}$Co and with $He^+(Cherchneff & Dwek 2009)$. Measurements of newly formed CO can shed light on the dust production rate.

Direct physical association of dust in high-z galaxies with high-z ccSNe is currently not possible, and thus to test the idea one must study dust production in local ccSNe. On 2017 May 14, a bright ccSN, 2017eaw, was discovered in the starburst galaxy NGC 6946. It has been classified as Type II-P (ATel
<table><caption>**Table 1**
Summary of Observations and Results of CO Modeling</caption><thead><tr><th>(yyyy mm dd)</th><th>Days</th><th>Exposure (s)</th><th>CO Mass (10<sup>-4</sup> M<sub>☉</sub>)</th><th>T<sup>a</sup> (K)</th><th>FWHM<sup>b</sup> (km s<sup>-1</sup>)</th><th>J, H, Ks<sup>c</sup> (mag)</th><th>Scale<sup>d</sup> Factor</th></tr></thead><tbody><tr><td>2017 Jun 05</td><td>22</td><td>240</td><td>...</td><td>...</td><td>...</td><td>...</td><td>1.000</td></tr><tr><td>2017 Jun 16</td><td>33</td><td>360</td><td>...</td><td>...</td><td>...</td><td>11.61 ± 0.04, 11.39 ± 0.05, 11.16 ± 0.04</td><td>0.505</td></tr><tr><td>2017 Jun 26</td><td>43</td><td>360</td><td>...</td><td>...</td><td>...</td><td>...</td><td>0.218</td></tr><tr><td>2017 Jul 06</td><td>53</td><td>360</td><td>...</td><td>...</td><td>...</td><td>...</td><td>0.112</td></tr><tr><td>2017 Aug 29</td><td>107</td><td>600</td><td>≤0.4</td><td>...</td><td>...</td><td>...</td><td>0.070</td></tr><tr><td>2017 Sep 15</td><td>124</td><td>600</td><td>0.6–1.6</td><td>3000</td><td>2800</td><td>...</td><td>0.051</td></tr><tr><td>2017 Oct 01</td><td>140</td><td>1200</td><td>1.0–1.9</td><td>3300</td><td>2850</td><td>...</td><td>0.012</td></tr><tr><td>2017 Oct 30</td><td>169</td><td>1200</td><td>1.6–2.2</td><td>3000</td><td>2850</td><td>13.65 ± 0.06, 13.45 ± 0.08, 13.31 ± 0.10</td><td>0.005</td></tr><tr><td>2017 Dec 05</td><td>205</td><td>1080</td><td>1.9–2.2</td><td>2700</td><td>2750</td><td>14.08 ± 0.07, 14.02 ± 0.10, 13.98 ± 0.12<sup>e</sup></td><td>0.002</td></tr></tbody></table>

**Notes.**

<sup>a</sup> σ ~ 200 K.

<sup>b</sup> σ ~ 200 km s<sup>-1</sup>.

<sup>c</sup> Photometry obtained by one of us (VJ) at 1.2 m Mount Abu Infrared Telescope, India.

<sup>d</sup> Spectra in Figure 1 have been multiplied by these values.

<sup>e</sup> The JHK imaging observations were taken on 2017 November 30.

10376; Xiang et al. 2017). The progenitor is known to be a red supergiant (RSG; ATel 10378; van Dyk et al. 2017; Kilpatrick & Foley 2018). Because of its proximity and location in the sky, SN 2017eaw is ideal for studying the onset and evolution of both molecule and dust formation from a northern hemisphere telescope. Here we present and describe nine Gemini near-infrared (NIR) spectra obtained during 2017 June–December and discuss the evolution of the spectrum, the detection of hot CO and an estimate of its mass, and emission from newly formed dust. ATel 10765 (Rho et al. 2017) is a brief report on the first six of these spectra. Estimates of the distance to NGC 6946 vary from 4 to 12.7 Mpc in the NASA/IPAC Extragalactic Database. We adopt 7.72 ± 0.32 Mpc from Anand et al. (2018).

## 2. Observations and Data Reduction

We obtained 0.81–2.52 µm spectra of SN 2017eaw using the Gemini Near-Infrared Spectrograph (GNIRS) on the 8.1 m Frederick C. Gillett Gemini-North telescope, for programs GN-2017A-DD-8 and GN-2017B-DD-11. The observing dates and exposure times are listed in Table 1. GNIRS was configured in its cross-dispersed mode, using its 32 line mm⁻¹ grating and a 0.45 arcsec-wide slit to provide a resolving power of ~1200 (250 km s⁻¹). Observations were made in the standard stare/nod-along-slit mode with a nod angle of 3 arcsec. Nearby early A-type dwarfs, observed either just before or just after SN 2017eaw to minimize differences in airmass, served as telluric and flux standards. Although the key questions that we wished to address required only K-band spectra, and in particular the longer wavelength portion of that band, where hot newly formed dust and CO are both expected to emit, the full spectral coverage provides a wealth of information regarding the evolution of the kinematics, ionization state, and elemental abundances.

Data reduction utilized both the GNIRS cross-dispersed reduction pipeline and manual, order-by-order reduction for the shortest wavelength orders and for instances in other orders where the pipeline appeared to produce spurious results. The manual reduction used standard IRAF<fnref n="7" /> and Figaro tools for

flatfielding, spike removal, rectification of spectral images, extraction, wavelength calibration, and removal of hydrogen absorption lines in the spectra of the standard stars. When order-by-order reduction was used, spectral segments covering different orders were stitched together after small scaling factors were applied, to produce final continuous spectra between 0.81 and 2.52 µm. The reliability of the spectra in the two intervals from 1.35 to 1.45 µm and from 1.80 to 1.95 µm is low because of low atmospheric transmission; one should view the spectra at these wavelengths with caution. In several of the spectra that were reduced manually, the H I lines in the telluric standard in the 1.80–1.95 µm interval were not removed prior to ratioing. At most other wavelengths, although the signal-to-noise ratio varies considerably, it is very high judging by the repetitiveness of the weaker features in spectra obtained on adjacent dates. Because some narrow-slit spectra were obtained in conditions of poor seeing and slightly non-photometric skies, flux densities could be uncertain by ±30%.

## 3. Overview of Spectra

Figure 1 shows the nine 0.8–2.5 µm spectra of SN 2017eaw obtained during 2017. The earliest spectrum was obtained 22 days after discovery, the latest at 205 days. These data form the most extensive set of NIR spectra of a Type II-P supernova obtained to date. We are making them available to the astronomy community via WISEREP.<fnref n="8" />

The first four spectra, obtained during the photospheric phase, are dominated by hydrogen recombination line emission and the Ca II triplet at 0.85 µm, all superimposed on a continuum that decreases monotonically with increasing wavelength. The spectrum at day 107 marks a transition, as part of the ejecta begins to turn nebular in density and numerous neutral atomic lines begin to appear across the entire wavelength range. According to the optical photometry of Tsvetkov et al. (2018), SN 2017eaw turned nebular (i.e., the nebular brightness follows the 0.01 mag day⁻¹ decay rate for ⁵⁶Co) at an age of ~125 days. At 124 days the atomic lines in the NIR spectra have become much more prominent, largely because the continuum has greatly decreased.

<footnote n="7">⁷ IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation.</footnote>

<footnote n="8">⁸ https://wiserep.weizmann.ac.il/</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Gemini GNIRS 0.81–2.52 µm spectra (dereddened) of SN 2017eaw obtained in 2017, in time order (top to bottom). Plotted spectra are multiplied by values in Table 1. Prominent and important lines are labeled. The wavelength interval with overtone CO emission is shaded.</figcaption></figure>

Spectral signatures of several CO overtone bands at 2.3–2.5 µm are evident in emission in the 124 days spectrum and afterward (see Figures 1 and 2), although the cut-on of the 2–0 CO emission near ~2.29 µm is not yet obvious until 169 days. There is marginal evidence for weak CO emission earlier, at day 107. In addition, beginning on day 124 there is a flattening of the continuum at $\lambda > 2.1$ µm, not seen in the earlier spectra, which we interpret as emission from newly formed hot dust. The flattening of the continuum beyond 2.1 µm is also present in the final three spectra.

# 4. Analysis and Discussion

## 4.1. Evolution of the Atomic Line Spectrum

During the photospheric phase, in addition to the H I and Ca II lines mentioned above, there are contributions from C I, Mg I, Mg II, and later S I (see Figure 1). Initially, the H I lines have P Cygni profiles and full extents of many thousands of km s⁻¹ but the latter decrease monotonically with time as shown in Figure 3. The extremely broad and shallow blueshifted absorption troughs at the H I Pa γ, Pa β, and Br γ lines in the early spectra may be the result of sustained absorption at high velocity by He I, which results from an ionization freeze-out of He in the outer ejecta (Dessart & Hillier 2008). Absorption by He I lines at 1.083, 1.253, 1.278, 2.113, and 2.165 µm may contribute to these extended troughs.

In the observed nebular phase, the SN luminosity has dropped by two magnitudes since the photospheric phase, and a much higher fraction of the ⁵⁶Co decay power emerges in lines. The nebular spectrum forms in the regions that absorb this decay power, located in what used to be the He core and the base of the H-rich envelope. While the lines present in the

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Zoomed-in CO first overtone band emission superposed on local thermodynamic equilibrium (LTE) model fits (red lines) at 124 and 205 days. Observed spectra are in black. Green and blue lines in uppermost panel are possible continua. Parameters of fits are in Table 1.</figcaption></figure>

photospheric phase (i.e., from H I, He I, Ca II) persist, additional lines from intermediate mass elements (O I, Na I, Mg I, Si I, S I, and Ca I), and iron-group elements (primarily Fe I and Fe II). Sr II 1.033 µm are observed (see Figure 1), especially around the transition to the nebular phase. The line profiles in the nebular phase are narrower than earlier and the peaks in some of lines are blueshifted by several hundreds of km s⁻¹ (Figure 3). As in the case of SN2010ji, which has similar velocity shifts (Gall et al. 2014), early formation of dust
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Upper panel: velocity profiles of hydrogen Pa β and Br γ and blended H/He Pa γ/1.083 µm on day 22. Doppler velocity relative to H I lines. Lower panel: evolution of the Pa β and Br γ line profiles.</figcaption></figure>

could be the cause of this; two other possible explanations are continuum opacity of the ejected gas and asymmetric ejecta.

We have generated a selection of simulations from Dessart et al. (2013) and Dessart et al. (2018) for representative Type II SNe from $15 M_{\odot}$ progenitors, exploded to yield ejecta masses of $12.5 M_{\odot}$ with $1.2 \times 10^{51}$ erg. For the photospheric phase the models are good matches to the observations. For example, the spectrum at 43 days corresponds to a model of a RSG explosion (model m15mlt3 from Dessart et al. 2013) with a photospheric temperature of 5450 K and a photospheric velocity of $4900 \text{ km s}^{-1}$. In the nebular phase, the overall agreement of models of blue supergiant explosion (model Bsm from Dessart et al. 2018) endowed with $0.084 M_{\odot}$ of $^{56}$Ni and the spectra is satisfactory, although the models overestimate the strengths of Mg I lines and underestimate the strengths of C I lines around $1 \mu\text{m}$. The apparent shallow continuum “bump” observed at 1.40–1.75 µm (see Figure 1) is due at least in part to numerous overlapping lines, largely from Fe I and Fe II.

## 4.2. CO and Dust

Carbon monoxide is first securely detected at 124 days. We have estimated the CO mass at days 124, 140, 169, and 205

using the LTE model developed by Das et al. (2009) and used by Banerjee et al. (2016), Joshi et al. (2017), and others. Two examples are shown in Figure 2. The model assumes pure $^{12}$C$^{16}$O and optically thin, LTE emission. The presence of CO in each of the spectra is clear, based on the close match of peaks in the spectrum to the wavelengths of the 3–1, 4–2, and 5–3 band heads (2.323, 2.351, and 2.383 µm, respectively), the fair match to the 2–0 band head at 2.294 µm (which may be contaminated by an atomic line, probably C I 2.291 µm) and, especially in the final two spectra, the large rise in flux density near the wavelength of the 2–0 band head. However, the strength of the CO emission and the CO mass derived from it are uncertain, due to the uncertainty in the location of the K-band continuum, especially on days 124 and 140. For self-consistency, in each of the final four 2017 spectra we initially adopted a flat continuum (due to the photosphere and presumed dust emission) passing through the bottom of the trough at 2.23 µm. In doing so we essentially assumed that atomic line emission accounts for the elevated “pseudo-continuum” at 2.0–2.3 µm on those two dates, but does not contribute significantly at longer wavelengths where the CO emits. However, it is quite possible that these minima are P Cygni troughs. If so, the above choice results in significant over-estimates of the CO mass on days 124 and 140, but not on days 169 and 205, as on those dates the “pseudo-continuum” is much weaker than the CO emission. Thus, the upper limits to the derived CO masses listed in Table 1 for 124 and 140 days should be regarded with caution. For example, if the adopted continuum in the top panel of Figure 2 is the blue line, the CO emission strength on that date would be ~3 times lower than that limit. The range of values of the CO mass in Table 1 reflect this uncertainty.

Best fit model spectra were determined by varying the CO mass, temperature, and velocity dispersion with the goal of minimizing the reduced $\chi^2$ value over the region 2.30 to 2.46 µm (thus excluding the 2–0 band head and the strong [Mg I] line at 2.48 µm). The fit results are summarized in columns 4–6 of Table 1. We obtain CO masses of ~$1 \times 10^{-4} M_{\odot}$ on all four dates from the LTE model and the assumed continua. The CO mass depends on the distance, $M_{7.72}$ ($d/d_{7.72}$)$^2$ where $M_{7.72}$ (in Table 1) and $d_{7.72}$ are the mass and distance assuming a distance of 7.72 Mpc, and d is an actual distance to SN 2017eaw. The masses on the first two dates may be considerably overestimated, as discussed above. On the other hand, if the molecular gas is not in LTE, the CO masses could be an order of magnitude larger, in analogy to the models of Liu et al. (1992) for SN 1987A. The best-fit CO temperatures and velocity dispersions for all four spectra are typically 3000 K and 2800 km s$^{-1}$, respectively.

Observations of CO in SNe are fairly rare; detections from about a dozen SNe, including a Type Ic (Hunter et al. 2009), are known and summarized in Sarangi et al. (2018). The evolution of the CO emission in ccSNe has only been studied over a lengthy time interval for SN 1987A (SNII-pec, a rare class but close to the Type II-P SN). In Figure 4 the estimated CO masses of SN 2017eaw are compared to those of SN 1987A and to the CO mass as a function of time predicted by Sarangi & Cherchneff (2013). In both objects overtone CO emission is detected for the first time shortly after 100 days: 112 days for SN 1987A (Spyromilio et al. 1988) and at 124 days for SN 2017eaw (implying that significant CO production began between 107 and 124 days). These dates
<figure><img src="image_6.png" /><figcaption>**Figure 4.** CO masses of SN 2017eaw (black, filled rectangles) plotted on model of CO produced by a 15 $M_{\odot}$ progenitor (Sarangi & Cherchneff 2013), for ejecta zones from innermost (1A) to outermost (4B). The vertical extents of black rectangles denote ranges of estimated CO masses, which depend on choice of 2.3–2.5 µm continuum. The zones include Si/S/Fe-rich (zone 1), O/Si/S-rich (zone 2), O/Mg/Si-rich (zone 3), and O/C/Mg-rich (zone 4) (Cherchneff & Dwek 2009). Zones 1 and 4 are subdivided into A and B subzones (see Table 1 of Sarangi & Cherchneff 2013). The gray line represents the mass summed over all zones. The largest amounts of CO primarily form in zones 4A and 4B. CO masses for SN 1987A derived for LTE (small triangles) and non-LTE (small squares) from models of Liu et al. (1992) and small circles for thermal chemical model of Liu & Dalgarno (1995) are also plotted.</figcaption></figure>

are remarkably consistent with the chemically controlled dust models by Sarangi & Cherchneff (2013, 2015), in which the first appearance of CO is predicted to be at ~100 days. The CO mass observed in SN 2017eaw (~$10^{-4} M_{\odot}$) is also roughly comparable to amounts in the models, in which most of the CO is produced in zones 4A and 4B (Sarangi & Cherchneff 2013). Individually each represents less than 0.1 % of the available C and O masses ($0.1 M_{\odot}$ and $1.0 M_{\odot}$, respectively, based on the models described above). However, CO may continue to form during the first few or a few ten years after the explosion and thus may result in a vastly larger amount of CO than observed to be emitting in the first overtone band at any one time.

We have interpreted the flattening of the K-band continuum starting at 2.1 µm, in the four GNIRS spectra obtained from 124 to 205 days, to be due to emission from hot newly condensed dust. When we fit the continua of the spectra above 2.1 µm using carbon dust, dust temperatures of ~1400 K and ~1200 K are consistent with the spectra at 124 days and 205 days, respectively. We suggest that the carbon dust formed in SN 2017eaw is graphitic carbon because it condenses at 1150–1690 K, whereas amorphous carbon condenses at ~1100 K (Ebel 2000; Fedkin et al. 2010). Dust formation at 124 days is much earlier than predicted by some models (Todini & Ferrara 2001; Nozawa et al. 2003). However, the models of Schneider et al. (2004), Sarangi & Cherchneff (2015), and Lazzati & Heger (2016) suggest that dust formation begins at ~150 days and continues until 2000 days. Recent observations show that dust formed at during 40–868 days in SN2010j1 (Gall et al. 2014) and continues to be present in SN 1987A as late as 25 years after outburst (Matsuura et al. 2011; Indebetouw et al. 2014; Zanardo et al. 2014).

The early sign of dust formation in SN 2017eaw appears to be supported by an excess of 3.6 µm emission seen on day 193 in Spitzer IRAC observations (S. Tinyanont et al. 2018, private

communication). This work serves as an early alert for further monitoring of the formation/evolution of the dust and CO in this object, particularly in the infrared.

As of the date of the final spectrum (UT 2017 December 5), the rate of CO production remained strong. Because of this, we plan to obtain additional NIR measurements of SN 2017eaw during 2018.

This paper is based on observations obtained at the Gemini Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the National Science Foundation (United States), the National Research Council (Canada), CONICYT (Chile), Ministerio de Ciencia, Tecnología e Innovación Productiva (Argentina), and Ministrio da Cincia, Tecnologia e Inovao (Brazil). J.R. acknowledges support in part for this work from NASA ADAP grant (NNX12AG97G). Research at the Physical Research Laboratory is funded by the Department of Space, Government of India. We thank Dr. Arka Sarangi for providing his models, and the anonymous referee for insightful comments.

## References

* Anand, G. S., Rizzi, L., & Tully, R. B. 2018, AJ, 156, 105

* Banerjee, D. P. K., Srivastava, M. K., Ashok, N. M., & Venkataraman, V. 2016, MNRAS, 455, L109

* Cherchneff, I., & Dwek, E. 2009, ApJ, 703, 642

* Das, R. K., Banerjee, D. P. K., & Ashok, N. M. 2009, MNRAS, 398, 375

* De Looze, I., Barlow, M. J., Swinyard, B. M., et al. 2017, MNRAS, 465, 3309

* Dessart, L., & Hillier, D. J. 2008, MNRAS, 383, 57

* Dessart, L., Hillier, D. J., Waldman, R., & Livne, E. 2013, MNRAS, 433, 1745

* Dessart, L., Hillier, D. J., & Wilk, K. D. 2018, A&A, accepted (arXiv:1808.08724)

* Draine, B. T. 2009, in ASP Conf. Ser. 414, Cosmic Dust—Near and Far, ed. T. Henning, E. Grün, & J. Steinacker (San Francisco, CA: ASP), 453

* Dwek, E. 1998, ApJ, 501, 643

* Ebel, D. S. 2000, JGR, 105, 10363

* Fedkin, A. V., Meyer, B. S., & Grossman, L. 2010, GeCoA, 74, 3642

* Gall, C., Hjorth, J., & Andersen, A. C. 2011, A&ARv, 19, 43

* Gall, C., Hjorth, J., Watson, D., et al. 2014, Natur, 511, 326

* Gomez, H. L., Krause, O., Barlow, M. J., et al. 2012, ApJ, 760, 96

* Hunter, D. J., Valenti, S., Kotak, R., et al. 2009, A&A, 508, 371

* Indebetouw, R., Matsuura, M., Dwek, E., et al. 2014, ApJL, 782, L2

* Isaak, K. G., Priddey, R. S., McMahon, R. G., et al. 2002, MNRAS, 329, 149

* Joshi, V., Banerjee, D. P. K., & Srivastava, M. 2017, ApJL, 851, L30

* Kilpatrick, C. D., & Foley, R. J. 2018, arXiv:1806.00348

* Lazzati, D., & Heger, A. 2016, ApJ, 817, 134

* Liu, W., & Dalgarno, A. 1995, ApJ, 454, 472

* Liu, W., Dalgarno, A., & Lepp, S. 1992, ApJ, 396, 679

* Matsuura, M., Dwek, E., Meixner, M., et al. 2011, Sci, 333, 1258

* Matsuura, M., Dwek, E., Barlow, M. J., et al. 2015, ApJ, 800, 50

* Michalowski, M. J. 2015, A&A, 577, A80

* Morgan, H. L., & Edmunds, M. G. 2003, MNRAS, 343, 427

* Nozawa, T., Kozasa, T., Umeda, H., Maeda, K., & Nomoto, K. 2003, ApJ, 598, 785

* Rho, J., Geballe, T. R., Banerjee, D. P. K., et al. 2017, ATel, 10765, 1

* Rho, J., Gomez, H. L., Boogert, A., et al. 2018, MNRAS, 479, 5101

* Rho, J., Kozasa, T., Reach, W. T., et al. 2008, ApJ, 673, 271

* Rho, J., Onaka, T., Cami, J., & Reach, W. T. 2012, ApJL, 747, L6

* Sarangi, A., & Cherchneff, I. 2013, ApJ, 776, 107

* Sarangi, A., & Cherchneff, I. 2015, A&A, 575, A95

* Sarangi, A., Dwek, E., & Arendt, R. G. 2018, ApJ, 859, 66

* Schneider, R., Ferrara, A., & Salvaterra, R. 2004, MNRAS, 351, 1379

* Spitzer, L. 1978, Physical Processes in the Interstellar Medium (New York: Wiley)

* Spyromilio, J., Meikle, W. P. S., Learner, R. C. M., & Allen, D. A. 1988, Natur, 334, 327

* Temim, T., & Dwek, E. 2013, ApJ, 774, 8

* Temim, T., Dwek, E., Arendt, R. G., et al. 2017, ApJ, 836, 129
Temim, T., Sonneborn, G., Dwek, E., et al. 2012, ApJ, 753, 72

Todini, P., & Ferrara, A. 2001, MNRAS, 325, 726

Tsvetkov, D. Y., Shugarov, S. Y., Volkov, I. M., et al. 2018, arXiv:1801.00340

van Dyk, S. D., Filippenko, A. V., Fox, O. D., et al. 2017, ATel, 10378, 1

Watson, D., Christensen, L., Knudsen, K. K., et al. 2015, Natur, 519, 327

Xiang, D., Rui, L., Wang, X., et al. 2017, ATel, 10376, 1

Zanardo, G., Staveley-Smith, L., Indebetouw, R., et al. 2014, ApJ, 796, 82

Zhukovska, S., Henning, T., & Dobbs, C. 2018, ApJ, 857, 94