Probing the Final-stage Progenitor Evolution for Type IIP Supernova 2017eaw in NGC 6946

TL;DR

This study analyzes the progenitor of supernova 2017eaw, identifying it as a red supergiant with initial mass around 12 solar masses, and suggests recent violent mass loss episodes prior to explosion.

Contribution

It provides detailed observational evidence linking the progenitor's properties and recent mass loss to supernova characteristics, offering insights into the final evolution of Type IIP supernova progenitors.

Findings

Progenitor identified as an M4-type red supergiant with about 12 solar masses.

Detected circumstellar material indicating recent violent mass loss.

Observed luminosity decline inconsistent with standard RSG evolution models.

Abstract

We presented a detailed analysis of progenitor properties of type IIP supernova 2017eaw in NGC 6946, based on the pre-explosion images and early-time observations obtained immediately after the explosion. An unusually red star, with M$_{F814W}$ = $-$6.9 mag and m$_{F606W}$$-$ m$_{F814W}=$2.9$\pm$0.2 mag, can be identified at the SN position in the pre-discovery Hubble Space Telescope(HST) images taken in 2016. The observed spectral energy distribution of this star, covering the wavelength of 0.6-2.0$\ \mathrm{\mu m}$, matches that of an M4-type red supergiant (RSG) with a temperature of about 3550 K. These results suggest that SN 2017eaw has a RSG progenitor with an initial mass of 12$\pm$2 M$_\odot$. The absolute F814W-band magnitude of this progenitor star is found to evolve from $-$7.2 mag in 2004 to $-$6.9 mag in 2016. Such a dimming effect is, however, unpredicted for a RSG in its neon/oxygen burning phase when its luminosity should modestly increase. The spectrum of SN 2017eaw taken a few hours after discovery clearly shows a narrow H$_\alpha$ emission feature blueshifted by $\sim$160 km s$^{-1}$. This narrow component disappeared in the spectrum taken two days later, suggesting the presence of a circumstellar material (CSM) shell (i.e., at a distance of $<$2.1-4.3$\times$10$^{14}$ cm). Combining the inferred distance with the expansion velocity of the CSM, we suggest that the progenitor of SN 2017eaw should have experienced violent mass loss at about 1-2 years prior to explosion, perhaps invoked by pulsational envelop ejection. This mechanism may help explain its luminosity decline in 2016 as well as the lack of detections of RSGs with initial mass in the range of 17 M$_\odot<$ M $<$ 25 M$_\odot$ as progenitors of SNe IIP.

Figures

Tables

Table 1. Table 1: Log of the pre- and post-explosion HST images covering the explosion site of SN 2017eaw

Instrument	Filter	Obs date	Exposure time	Prospsal ID	PI	Drizzled imagea
ACS/WFC	F814W	2016-10-26	$600+590+690\times 2$	14786	Williams	jd9g03020_drc.fits
ACS/WFC	F606W	2016-10-26	$530+520+690\times 2$	14786	Williams	jd9g03010_drc.fits
ACS/WFC	F814W	2004-07-29	120	9788	Ho	hst_9788_d8_acs_wfc_f814w_drz.fits
ACS/WFC	F658N	2004-07-29	$350\times 2$	9788	Ho	j8mxd8010_drz.fits
WFC3/IR	F164N	2016-10-24	$599\times 4$	14638	Long	id6o03010_drz.fits
WFC3/IR	F160W	2016-10-24	$99\times 4$	14638	Long	id6o03020_drz.fits
WFC3/IR	F128N	2016-02-09	$353\times 4$	14156	Leroy	icyqe2030_drz.fits
WFC3/IR	F110W	2016-02-09	$228\times 2$	14156	Leroy	icyqf2010_drz.fits & icyqe2040_drz.fits
WFC3/UVIS	F814W	2017-05-29	$27\times 10$	14645	van Dyk	id9702010_drz.fits

Table 2. Table 2: Transformed SN positions on pre-explosion HST images using post-explosion TNT image

image	F128N	F110W(a)∗	F110W(b)∗∗	F814W(2004)	F658N
$s_{\mathrm{pre}}/{\left(\mathrm{arcsec}\cdot {\mathrm{pixel}}^{-1}\right)}^a$	0.13	0.13	0.13	0.05	0.05
$N_{\mathrm{stars}}^b$	25	24	26	9	14
$s_{\mathrm{post}}/{\left(\mathrm{arcsec}\cdot {\mathrm{pixel}}^{-1}\right)}^c$	0.52	0.52	0.52	0.52	0.52
${\sigma }_{\mathrm{total}}/{\mathrm{mas}}^d$	37/41	36/38	41/40	26/25	31/33
${(x,y)}_{\mathrm{SN}}^e$	(558.28, 650.24)	(569.27, 656.08)	(539.49, 633.24)	(3360.32, 3978.51)	(877.62, 267.01)
${(x,y)}^{\ast f}$	(558.23, 650.44)	(569.28, 656.50)	(539.65, 633.47)	(3360.48, 3978.36)	$\cdots$

Table 3. Table 3: Photometry of the Progenitor Star Identified for SN 2017eaw.

Telescope	Instrument	Filter	Exposure time(s)	ma
HST	ACS/WFC	F606W	530+520+690 $\times$ 2	26.419(053)
HST	ACS/WFC	F658N	350$\times$2	$>$23.7
HST	ACS/WFC	F814W(2004)	120	22.558(035)
HST	ACS/WFC	F814W(2016)	530+520+690 $\times$ 2	22.845(008)
HST	WFC3/IR	F110W	228 $\times$ 2	20.344(005)
HST	WFC3/IR	F128N	353 $\times$ 4	19.749(019)
HST	WFC3/IR	F160W	99 $\times$ 4	19.351(007)
HST	WFC3/IR	F164N	599 $\times$ 4	19.023(012)
Spitzer	IRAC	$3.6\mu m$	2144	18.02(070)
Spitzer	IRAC	$4.5\mu m$	2144	17.75(100)
Spitzer	IRAC	$5.8\mu m$	2144	$>$17.6
Spitzer	IRAC	$8.0\mu m$	2144	$>$16.9
Spitzer	MIPS	$24\mu m$	2144	$>$12.6

Full text

††pubyear: XXX††pagerange: Probing the Final-stage Progenitor Evolution for Type IIP Supernova 2017eaw in NGC 6946–Probing the Final-stage Progenitor Evolution for Type IIP Supernova 2017eaw in NGC 6946

1 Introduction

Massive stars with initial masses ranging from 8 to 25 M⊙ usually end their lives as type II core-collapse supernovae (Heger et al., 2003). Among them, type IIP supernovae (SNe IIP) are the dominant subclass, which are characterized by prominent Balmer lines in the spectra and a plateau phase lasting for about 100 days in the light curves (Barbon et al., 1979; Filippenko, 1997; Branch & Wheeler, 2017). These hydrogen-rich SNe are identified to arise from core-collapse of red supergiants (RSGs) with initial mass lying between 8 $-$ 16.5 M*⊙* (Li et al., 2007; Maund et al., 2008; Smartt, 2009, and reference therein). However, the observation of RSGs in Milky Way, Large Magellanic Cloud(LMC) and Small Magellanic Cloud(SMC) indicate a mass range of RSGs from 9 $-$ 25M⊙ (Levesque et al., 2005; Levesque et al., 2006, 2007). The lack of RSGs with mass in the range of 17 $-$ 25 M⊙ as progenitors of SNe IIP challenges current theory of stellar evolution. One possible explanation is that some massive RSGs experience prominent mass loss during its final stage evolution towards the explosion, due to the stellar wind, binary interaction, or pulsational eruption. The accumulated circumstellar materials (CSM) may obscure the RSGs and lead to the underestimation of the progenitor luminosity and its initial mass. The condensed dust shell will leave an imprint on the SN spectra via photoionization or interaction with the SN ejecta (e.g. Yaron et al., 2017; Hosseinzadeh et al., 2018).

Signatures for the presence of nearby CSMs have been reported for SNe IIP/IIL/IIb. Early spectra of type IIb SN 2013cu (Gal-Yam et al., 2014) and type IIP SN 2013fs (Yaron et al., 2017) reveal flash-ionized emission lines from dense circumstellar wind. The narrow emission lines are also detected in the early-time spectra taken in few hours or days after explosion, i.e., SN 2016bkv (Hosseinzadeh et al., 2018), SN 2006bp (Quimby et al., 2007) and some young type IIP/IIL SNe reported by Khazov et al. (2016). The narrow interacting features could be generated by recombination of the CSMs ionized by ultraviolet (UV) radiation emitted during the shock breakout or the shock-cooling phase from the supernovae. However, the mechanism driving the mass loss of RSGs, i.e., whether the mass loss is stable or unstable, is not fully understood. Dense CSM surrounding a typical type IIP supernova is found to have a expansion velocity that is apparently higher than typical wind velocity of RSGs (Yaron et al., 2017) (i.e., 130 km s*-1* versus 10 km s*-1*), suggestive of eruptive mass loss during the late phase evolution of RSGs. Moreover, Yoon & Cantiello (2010) found that strong pulsation induced by partial ionization of hydrogen in the envelope is expected to enhance the mass loss during the RSG phase.

SN 2017eaw provides another rare opportunity to study the late-time evolution of the red supergiant with relatively low initial mass. This SN was discovered by Patrick Wiggins on 2017 May 14.238 (UT dates) in NGC 6946 at a distance of 5.5 Mpc, with an unfiltered CCD magnitude of 12.8 mag (Wiggins, 2017), and it was immediately classified as a young type IIP supernova (Xiang et al., 2017b). Note that NGC 6946 is a well-known nearby Scd galaxy that has recorded nine SNe before SN 2017eaw, which makes it one of the most prolific SN factories known in the local universe. Follow-up observations of SN 2017eaw were thus initiated less than 0.6 d after the discovery, including multi-color ultraviolet/optical photometry and rapid high-cadence spectroscopy. The UBVRI-band photometric evolution covering the first 200 days has been presented by Tsvetkov et al. (2017), which shows a long plateau characterized by a normal type IIP SN. Kilpatrick & Foley (2018) studied the progenitor properties of SN 2017eaw, suggesting that it has a red supergiant progenitor with a 13 M*⊙*. Based on the historical Spitzer data, their studies also suggested that the progenitor star might have experienced dust-shell formation a few years before the explosion.

In this paper, we also attempt to constrain the progenitor properties of type IIP supernova 2017eaw using photometric and spectroscopic observations obtained immediately after explosion and the archival images from HST and Spitzer. The description of early-time observation and brief data reduction process are presented in Section 2. In Section 3, we describe the identification of the progenitor candidate from HST and Spitzer in detail. We discuss the progenitor properties and explore its possible final-stage evolution before explosion in Section 4. Our conclusions are given in Section 5.

2 Early-time Followup Observations

The UBVRI photometry observations of SN 2017eaw were triggered immediately after the SN discovery using the Tsinghua-NAOC 80-cm telescope(TNT) (Huang et al., 2012) at Xinglong Observatory of NAOC (Located in Hebei, China), with the first observations obtained at t$\sim$0.6 days from discovery. The SN was observed by the Swift UVOT (Gehrels et al., 2004) even earlier in three ultraviolet (UV) (uvw2, uvm2, uvw1) and three optical filters (U, B, V), starting at t$\sim$0.4 days after discovery. Based on the explosion time deduced below, these two phases correspond to 1.6 days and 1.4 days after explosion, respectively.

For data obtained by TNT, the data were first processed with IRAF111IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities fo Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation (NSF)., this includes bias subtraction and flat field correction. The instrumental magnitudes were measured with point-spread function (PSF) photometry and photometric calibration was determined with 10 stars in the field of SN 2017eaw observed on photometric night. We adopted aperture photometry to reduce all Swift images using HEASoft(the High Energy Astrophysics Software)222http://www.swift.ac.uk/analysis/software.php and detailed data descriptions are addressed in Rui et al. (2019; in prep).

Figure 1 displays the evolution of the UV and optical photometry of SN 2017eaw up to $\sim$45 days after the discovery. The unfiltered light curve derived from the observations by Patrick Wiggins is overplotted. The first nondetection can be traced back to 2017 May 12.20 UT, with an unfiltered CCD magnitude upper limit of about 19.0 mag. The explosion time (or first light) can be estimated as JD 2457886.72$\pm$1.01 by taking the mean point between the first detection and the last non-detection. The rising phase was well captured in the $UBVRI$ bands, while the UV data already started to decline from the first data points. We fit the TNT $R$-band light curve during the first week after explosion with a $f\propto(t-t_{0})^{n}$ function, and find that the first (unfiltered) discovery point is brighter than the fitted curve by about 0.2 mag. After normalizing the unfiltered magnitudes to the R-band values in the plateau phase, this difference is about 0.16 mag, which could be related to the emission of shock breakout cooling tail. Such post-shock cooling decline before reaching the optical peak has been reported for SN 2014cx (Huang et al., 2016), SN 2016X (Huang et al., 2018), SN 2006bp (Quimby et al., 2007), KSN 2011d (Garnavich et al., 2016). We caution, however, that the evidence for shock cooling detection in SN 2017eaw should be less significant due to that the difference between the unfiltered and R-band magnitudes change with time as a result of evolution of the spectral energy distribution of SNe IIP in the early phase.

By fitting a low-order polynomial to the data around maximum light, the supernova reached the $U$- and $B$-band peak of 12.39 mag and 13.16 mag at t$\sim$+5.0 days and t$\sim$+5.1 days from the discovery respectively. We note that the rise time in $B$ band is almost the shortest among a large sample of SNe II (González-Gaitán et al., 2015).

For the spectrum of SN 2017eaw displayed in Fig. 2, the first two were obtained with BFOSC and the third one was obtained with OMR on the Xinglong 2.16-m telescope. All spectra were reduced using standard IRAF routines. Flux of these spectra were calibrated using nearby spectrophotometric standard stars with similar airmass. Atmospheric extinction was corrected basing on the extinction curve of Xinglong observatory.

The earliest spectrum, obtained at $\sim$1.4d after explosion, shows a blue continuum and weak features of Balmer lines that are characteristic of young type II supernovae (Xiang et al., 2017a). These Balmer lines become dominant features in the following spectra, similar to those seen in SNe IIP(Fig. 2). Broad absorption feature of He I $\lambda$ 5876 become visible in the $+$3.4d spectra. The spectrum of SN 2017eaw at t$=+$8d is similar to other well-studied type IIP supernova(especially SN 2004et), while showing weaker and broader profiles of Balmer and He I lines compared to SN 1999em. Noticeable narrow interstellar Na I D absorption features can be seen at around 5892.5Å in the spectra, with an equivalent width (EW) of 1.6$\pm$0.4Å. Owing to the very close distance of the SN and the low resolution of our spectra, it is difficult to separate the absorption component due to the Milky way from that due to NGC 6946. The reddening towards SN 2017eaw can be estimated with the strength of this interstellar absorption line. Based on the existing empirical relation between the EW of Na ID absorption and dust extinction (Barbon et al., 1990; Faran et al., 2014), the total extinction of AV=1.24$\pm$0.30 or 2.41$\pm$0.62 mag can be inferred for SN 2017eaw. Due to this large discrepancy, we adopted an average value of AV=1.83$\pm$0.59 mag for the extinction correction of the progenitor. Given the Galactic extinction of AV =0.94 (Schlafly & Finkbeiner, 2011) towards SN 2017eaw, the host-galaxy extinction is estimated as AV=0.89$\pm$0.59 mag and it is adopted in the progenitor analysis.

3 Identification of the Progenitor Candidate

3.1 Images from Hubble Space Telescope

Pre-explosion images of SN 2017eaw are available from Mikulski Archive for Space Telescopes (MAST) and the Hubble Legacy Archive (HLA). By searching the publicly available archive images, we found there are F128N- and F110W-band (WFC3/IR) images that were taken on 9 February 2016 (PI Leroy), F164N- and F160W-band images taken on 24 October 2016 (PI Long), and F606W- and F814W-band (ACS/WFC) images taken on 26 October 2016 (PI Williams). Besides many images taken in 2016, there are also F658N- and F814W-band (ACS/WFC) images taken on 29 July 2004 (PI Ho). Table 1 list the observation log of these images.

To measure the SN position on these HST images, we used $I$-band image of SN 2017eaw taken by the Tsinghua-NAOC 0.8-m telescope (TNT) on 17 May, which has a typical full-width at half maximum (FWHM) of $1.6^{\prime\prime}$. The position of SN 2017eaw was determined by averaging the results from six different methods (centroid, Gaussian and ofilter center algorithm of IRAF task phot, IRAF task imexamine, daofind, and SExtractor) as $(735.154\pm 0.007,644.454\pm 0.011)$. The uncertainty of this position is estimated as the standard deviation of the mean value. Several stars commonly seen on the HST images and the $I$-band TNT image are used to establish the transformation function (a 2nd-order polynomial) converting the coordinates of the $I$-band TNT image to that of the HST images. The resulting positions and uncertainties are listed in Table 2. The uncertainties include errors in SN position and geometric transformation. The reference stars used in position transformation varies in different HST wavebands, and the exact number is listed as $N_{\mathrm{stars}}$ in Table 2.

Figure 3 shows the region centering on the SN site with the transformed position marked. A point-like source can be detected near the position of the SN on the F128N-, F110W- and F814W-band images, while no detection on the F658N-band image. Difference between the positions of SN and this point-like source is generally less than the uncertainty of the transformed positions.

On 29 May 2017, an F814W-band image of the SN was taken again by the HST with WFC3/UVIS (PI Van Dyk), which became publicly available immediately. We retrieved the image and measured the SN position, with the X-pixel being at 547.90$\pm$0.02 and Y-pixel being at 576.79$\pm$0.04. This pixel position is then transformed at the SN positions on the pre-explosion HST images (see also Figure 3), which has higher precision but is overall consistent with those derived from the TNT images. Based on these two results, we concluded that the point-like source marked in Figure 3 is the candidate of the progenitor of SN 2017eaw.

Photometry of the pre-explosion HST images is conducted with the DOLPHOT 2.0333http://americano.dolphinsim.com/dolphot/ package on the bias-subtracted, flat-corrected, so-called FLT FITS images that are retrieved from the MAST. Photometry is simultaneously measured using all ACS/WFC images and WFC3/IR images, following the routines that include the selection of a deep drizzled image as reference image, masking of bad pixels, calculations of sky background, alignment of images to the reference, and adoption of the recommended photometry parameters as described in the User’s Guide of DOLPHOT. Magnitudes and their uncertainties of the progenitor candidate are extracted from the output of DOLPHOT. The Magnitudes are PSF magnitudes, and all in Vega system. An upper limit of magnitude is given for the non-detection F658N image. The measured magnitudes of the SN progenitor in different wavebands and phases are reported in Table 3, with F606W-band and F814W-band magnitudes being 26.419 mag and 22.845 mag. The magnitudes of the SN progenitor were also independently measured by Kilpatrick & Foley (2018) using the published HST images. In F606W, F814W, F160W bands, our results are consistent with theirs within 0.05 mag, while in F110W and F164N their results are fainter than ours by 0.4 mag and 0.1 mag, respectively.

3.2 Images from the Spitzer Telescope

The position calibration step was skipped for the Spitzer images, as the retrieved WCS information is well consistent with that from the GSC-2.3 standard catalogue. We find the counterpart of the progenitor candidate at the expected coordinates on $3.6\ \mathrm{\mu m}$ and $4.5\ \mathrm{\mu m}$ images (Figure 4). On the $5.8\ \mathrm{\mu m}$, $8.0\ \mathrm{\mu m}$ and $24\ \mathrm{\mu m}$ images, there are no clear point-like sources at the SN site. Photometry of the progenitor on Spitzer images is directly taken from (Khan, 2017a), and the limiting magnitudes of non-detection images are estimated using the same method. It should be pointed out the flux at the SN site of these images could be seriously contaminated by nearby sources that are not related to the progenitor of the SN(Khan, 2017b) due to the low angular resolution of the Spitzer images. For example, the photometric aperture is adopted as 2.4" for the 3.6um/4.5um images(Khan, 2017a), which corresponds to 48 pixels on the HST ACS images and 19 pixels on the WFC3/IR images, respectively. Adopting this large aperture for the HST images would result in an overestimate of the F814W-band magnitude by about 1.0 mag (i.e., $\sim$21.8 mag). Thus, we caution the use of the Spitzer magnitudes to constrain the progenitor properties of SN 2017eaw.

4 Analysis

4.1 Constraining Properties of the Progenitor

SN 2017eaw provides a good opportunity to study the progenitor of SNe IIP, as the pre-explosion Hubble Space Telescope (HST) images are available in both optical and near-infrared bands. Inspecting the HST WFC3 images taken in 2016 reveals that a point source can be identified at the SN position in the F606W-, F814W-, F110W-, F128N-, F160W-, and F164N-bands (see Fig. 3 and Sec. 3.1 for details of the progenitor identification). This point source tends to become progressively brighter at longer wavelengths, consistent with a very red star (see Table 1 and Fig.3). In the F814W band, the SN progenitor is measured to have a magnitude of 22.845$\pm$0.008 mag, which corresponds to an absolute magnitude of $-$6.93$\pm$0.51 mag after adopting a distance of 5.5 Mpc ($\mu$=28.67$\pm$0.43) and AV = 1.83$\pm$0.59 mag. This brightness is comparable to that of a red supergiant star. Such identification is also favored by the detection of a luminous and red source in the archival Spitzer images at $3.6\ \mathrm{\mu m}$ and $4.5\ \mathrm{\mu m}$ wavelengths (see Fig.4). The point source has an absolute magnitude of $-$10.65$\pm$0.38 mag at IRAC [3.6$\mu$ m] and a color of 0.27$\pm$0.06 in IRAC [3.6$\mu$ m]-[4.5$\mu$ m], which are consistent with the typical values of RSGs(Szczygieł et al., 2010). Note, however, that the nearby red stars could contaminate the Spitzer magnitudes due to the low angular resolution of the Spitzer images.

It is interesting that the region including the site of SN 2017eaw was also observed in 2004 by the HST Advanced Camera for Surveys (ACS). These observations were made in the F658N and F814W bands, where the progenitor can be clearly detected in F814W with mF814W = 22.558$\pm$0.005 mag. In comparison with this value, the SN progenitor appears dimmer by $\sim$0.29 mag in 2016, which is likely related to circumstellar dust formed around the progenitor shortly before its explosion (see discussions below). Note that the non-detection in the F658N band is due to the relatively short exposure time and narrow bandwidth of the filter. The measured magnitudes of the SN progenitor in different wavebands and phases are reported in Table 3.

To determine the spectral type of the progenitor of SN 2017eaw, we use 85 observed RSGs with effective temperatures ranging from 3325$K$ to 4200$K$ to fit the observed spectral energy distribution (SED). This sample consists of 74 Galactic RSGs(Levesque et al., 2005), 7 RSGs in the Large Magellanic Cloud (LMC), and 4 RSGs in the Small Magellanic Cloud (SMC)(Levesque et al., 2007). A second-order polynomial is used to fit the $BVRIJH$-band magnitudes of the observed RSGs and then derive their SED curves. These curves are then compared with the reddening-corrected SED of the progenitor. The comparison shows that the progenitor matches well with an M4-type red supergiant with an effective temperature of about 3550$\pm$100 K (see the upper panel of Fig.5). In comparison, Kilpatrick & Foley (2018) derived the effective temperature of 3350${}^{+450}_{-250}$ K by fitting the SED of progenitor with stellar SED and dust models. Based on the derived absolute K-band magnitude (i.e., $-$10.33$\pm$0.20 mag) and the bolometric correction inferred from this band (Levesque et al., 2006), the bolometric magnitude or luminosity of the progenitor is derived as $-$7.45$\pm$0.50 mag or log($L_{bol}/L_{\odot}$) = 4.88$\pm$0.20. With the derived temperature and bolometric luminosity, the progenitor star can be well located in the Hertzsprung-Russel diagram. Inspection of stellar evolutionary tracks(Ekström et al., 2012) shown in the lower panel of Fig.5, one can see that the progenitor of SN 2017eaw is likely to be a RSG star with an initial mass of 10 M⊙ to 14 M⊙ and a radius of about 700 R*⊙*.

The radius of progenitor can be estimated from the evolution of photospheric temperature in the early stage of the explosion, which is primarily determined by the radius, opacity, and density profile of the progenitor (Rabinak & Waxman, 2011). This early-time temperature evolution can be derived using the ultraviolet (UV) and optical photometry obtained within the first week after the explosion. Adopting the typical density profile as fρ = 0.13 and Thomson scattering opacity as 0.34 cm2 g*-1*, the progenitor radius of SN 2017eaw can be determined as 636$\pm$155R*⊙* (Fig.6). Moreover, the rise time of the light curve can be also used to estimate the radius of the progenitor of type IIP supernova through the relation between the progenitor radius and the rise time of the light curve, i.e., log R[R⊙] = (1.225$\pm$0.178) log trise [day] + (1.692$\pm$0.490) (Morozova et al., 2016). For SN 2017eaw, the rise time is estimated to be trise = 6.81$\pm$0.85 day in the g band, which gives an estimate of the radius as 515$\pm$180 R*⊙*. The above two estimates of the progenitor radius agree within 1-$\sigma$ error.

We also try to fit the observed SED of the progenitor star from the prediscovery HST and Spitzer images with synthesized spectra calculated using the stellar evolution code MARCS444http:marcs.astro.uu.se/. However, the best-fit model spectrum has an effective temperature of only about 2500 K (A similar value of 2600 K was also suggested by Kilpatrick & Foley (2018)) with AV=1.83$\pm$0.59 mag by including the Spitzer data(Fig.7), much lower than that inferred for the HST data alone. Given the luminosity, the lower temperature would result in an effective radius of 1330 R⊙ for the progenitor, which is inconsistent with the estimate either from the rise time or temperature evolution of the shock cooling in the early phase. Despite an overall reasonable fit, the measurements in two Spitzer bands may suffer large contaminations of the nearby red sources due to large aperture photometry and/or the presence of circumstellar dust. In the latter case, the light of the progenitor star is scattered/absorbed by the surrounding dust and re-emitted at longer wavelengths (see also the analysis by (Kilpatrick & Foley, 2018)). We also derive the F606W-band magnitude in 2004 by assuming that the dimming behavior in F814W-band is totally due to the CSM dust. Taking into account this speculated F606W-band magnitude and the observed F814W-band magnitude obtained in 2004, the SED of the progenitor inferred in 2004 matches with the theoretical stellar spectra model characterized by an effective temperature of 3550K$\pm$150K and a radius of 636$\pm$155R⊙ (see Fig. 7). As an alternative, the dimming behavior in F814W band could be also due to the intrinsic variation of a RSG star, but this explanation is inconsistent with either the observation taken at one year to a few days before the explosion or the theoretical prediction for its final stage evolution (see discussions in Sec 4.3).

4.2 The Dust Around the Progenitor

Some recent studies reveal that the ion-flashed features can be detected in very early spectra of SN 2013fs (Yaron et al., 2017), and SN 2016bkv (Nakaoka et al., 2018), which proves the existence of CSM around SNe IIP. The detection of prominent asymmetric Hα emission feature in the t=$+$1.4 day spectrum provides strong evidence for the existence of CSM dust around SN 2017eaw (see Fig.8). This weak feature is also detected in a spectrum obtained an hour earlier (Cheng et al., 2017). The asymmetric Hα line profile seen in SN 2017eaw can be decomposed into two components, with the broad component formed in the SN ejecta and the weak narrow peak formed due to the ionization of the surrounding CSM. The weak emission feature is measured to have a wavelength of 6559.50$\pm$1.31Å , corresponding to an expansion velocity of about 163 km s*-1* for the CSM, as shown in Fig.9. This measurement of velocity is validated by carefully inspecting the position of sky emission lines [O I] $\lambda$5577 and [O I] $\lambda$6300 in the spectrum to ensure the accuracy of wavelength calibration (Fig.10). The stellar wind inferred from the ion-flashed signature in the spectrum of SN 2013fs shows a similar velocity (i.e., $\sim$100 km -1) (Yaron et al., 2017).

Similar asymmetric Hα line profile can also be seen in the earliest spectra of SN 2006bp and iPTF13dkk (Khazov et al., 2016), but only narrow emission features are shortly detected in the flash spectra of SN 2013fs, SN 2016bkv, and SN 2018zd. In the very early phase, the spectra of SN 2013fs and SN 2006bp also show prominent He II lines, while this feature is blended with C III/N III in iPTF13dkk, SN 2016bkv, and SN 2017eaw. Such an observed diversity of the flashed spectroscopy suggests that the CSM/dust environments around their progenitor stars are quite diverse.

The duration of narrow emission features also varies for different SNe II (see Fig.8). For SN 2017eaw, the weak narrow Hα emission lines seen at t$\sim$+1.4 days disappeared in the t$\sim$+3.4 spectrum, suggesting that the obscuring dust shell should be thin and close to the progenitor star. Such a quick variation can be explained by that the CSM was initially ionized by the supernova photons, and then it was captured and destroyed by the fast-expanding SN ejecta. Assuming that the outermost layer of the SN ejecta has a velocity of  15,000-20,000 km s*-1* as inferred from the blue-shifted absorption of Hα and it swept the CSM within two days, and we estimate that the CSM has a distance $<$2.1-4.3$\times$1014cm from the supernova. With this distance, one can speculate that the progenitor star ejected the CS materials at 1-2 years before the explosion with the derived stellar wind velocity of about 160 km s*-1*.

4.3 Variability of the progenitor

To understand the dimming behavior of the progenitor star in F814W-band over the period from 2004 to 2016, we examined in Figure 11 the magnitude variation of the stars within 200 pc from the progenitor. On a timescale of a decade, most of the field stars are found to have light variations less than  0.3 mag when considering the larger photometric uncertainty at the faint end (i.e., m${}_{F814W}>$24.0 mag; see left panel of Fig.11); while the SN progenitor and star G are the two ones that show noticeable light variations with a relatively high significant level. Although a significant magnitude change might be common for a red supergiant star during its evolution, becoming faint is not predicted by theory for a red supergiant in its final stage towards explosion. In 2004 and 2016, the progenitor of SN 2017eaw should be near the end of carbon burning and in the neon/oxygen burning phase, respectively (Woosley & Janka, 2005). In the neon/oxygen burning stage, the density wave energy is mainly deposited into the hydrogen envelope and consequently the surface luminosity will modestly increase instead of becoming faint (Fuller, 2017). Thus the dimming behavior observed for the progenitor of SN 2017eaw is unlikely due to the intrinsic variability of a RSG. Moreover, the Ks-band observations of the progenitor star with Keck/MOSFIRE and P200/WIRC, obtained at one year to a few days before the explosion, did not show significant light variation at a 6% level Tinyanont et al. (2019), which is consistent with our speculation that the light variation is not intrinsic to the red supergiant.

We further examined the distribution of the mF606W-mF814W color as measured from the 2016 HST images. After corrections for the reddening, the progenitor is found to have a mF606W-mF814W color of 2.94$\pm$0.20 mag, which is noticeably redder than the nearby RSGs. Its position is very close to the instability regime of supergiants with mass in a range of 10-15 M*⊙, in particular when considering the effect of CSM dust (see Fig.11). On the other hand, instability mass ejection is also needed to reproduce the coolest (or reddest) variables with masses lower than 10-15 M⊙* according to the census of luminous stellar variability in M51 (Conroy et al., 2018). In the supergiant instability regime, a pulsational superwind could lead to a dramatically enhanced mass loss at late phases, which may account for an ejection of the H-envelope of RSG stars just before their explosions (Yoon & Cantiello, 2010; Conroy et al., 2018).

5 Conclusion

From the above analysis, we conclude that the normal type IIP supernova 2017eaw has an M4-type red supergiant star with an initial mass 12$\pm$2 M*⊙, which further confirms the trend that most SNe IIP arise from core-collapse explosions of RSGs. However, the multi-epoch HST images obtained before the SN detection reveal that the RSG progenitor became faint by 30% one year before it exploded. Although this dimming phenomenon could be due to the intrinsic light variation of the RSG star, detailed analysis of the neighbouring stars and early time spectra suggest that it is more likely due to the obscuration of circumstellar dust shell newly-formed at a distance of $<$2.1-4.3$\times$1014* cm from the star. The fast-moving circumstellar materials (at a velocity of $\sim$ 160 km s*-1*) indicate that they were ejected from the star in a violent manner within a few years prior to the explosion. Thus, the observations of SN 2017eaw indicates that some low-mass RSGs could also experience a strong mass loss, perhaps invoked by pulsational instability, during its final-stage evolution towards explosion.

Acknowledgements

We thank Patrick Wiggins for useful unfiltered images. This work is supported by the National Natural Science Foundation of China (NSFC grants 11325313 and 11633002), and the National Program on Key Research and Development Project (grant no. 2016YFA0400803). J.-J. Zhang is supported by NSFC (grants 11403096, 11773067), the Youth Innovation Promotion Association of the CAS(grants 2018081), the Western Light Youth Project, and the Key Research Program of the CAS (Grant NO. KJZD-EW-M06). T.-M. Zhang is supported by the NSFC (grants 11203034). This work was also partially Supported by the Open Project Program of the Key Laboratory of Optical Astronomy, National Astronomical Observatories, Chinese Academy of Sciences. The research of JRM is supported through a Royal Society University Research Fellowship.