Asphericity, Interaction, and Dust in SN 2013ej

Asphericity, Interaction, and Dust in the Type II-P/II-L Supernova 2013ej in Messier 74

Jon C. Mauerhan11affiliation: Department of Astronomy, University of California, Berkeley, CA 94720-3411, USA 2 2affiliationmark: , Schuyler D. Van Dyk33affiliation: Infrared Processing and Analysis Center, California Institute of Technology, 1200 E. California Blvd., Pasadena, CA, 91125 USA , Joel Johansson44affiliation: Department of Particle Physics and Astrophysics, Weizmann Institute of Science, 234 Herzl St., Rehovot, Israel , Maokai Hu55affiliation: Department of Physics, Texas A&M University, College Station, TX 77843, USA , Ori D. Fox66affiliation: Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA , Lifan Wang55affiliation: Department of Physics, Texas A&M University, College Station, TX 77843, USA , Melissa L. Graham11affiliation: Department of Astronomy, University of California, Berkeley, CA 94720-3411, USA 77affiliation: Department of Astronomy, University of Washington, Box 351580, U.W., Seattle, WA 98195-1580, USA , Alexei V. Filippenko11affiliation: Department of Astronomy, University of California, Berkeley, CA 94720-3411, USA , and Isaac Shivvers11affiliation: Department of Astronomy, University of California, Berkeley, CA 94720-3411, USA ,

Supernova (SN) 2013ej is a well-studied core-collapse supernova (SN) that stemmed from a directly identified red supergiant (RSG) progenitor in galaxy M74. The source exhibits signs of substantial geometric asphericity, X-rays from persistent interaction with circumstellar material (CSM), thermal emission from warm dust, and a light curve that appears intermediate between supernovae of Types II-P and II-L. The proximity of this source motivates a close inspection of these physical characteristics and their potential interconnection. We present multiepoch spectropolarimetry of SN 2013ej during the first 107 days and deep optical spectroscopy and ultraviolet through infrared photometry past days. SN 2013ej exhibits the strongest and most persistent continuum and line polarization ever observed for a SN of its class during the recombination phase. Modeling indicates that the data are consistent with an oblate ellipsoidal photosphere, viewed nearly edge-on and probably augmented by optical scattering from circumstellar dust. We suggest that interaction with an equatorial distribution of CSM, perhaps the result of binary evolution, is responsible for generating the photospheric asphericity. Relatedly, our late-time optical imaging and spectroscopy show that asymmetric CSM interaction is ongoing, and the morphology of broad H emission from shock-excited ejecta provides additional evidence that the geometry of the interaction region is ellipsoidal. Alternatively, a prolate ellipsoidal geometry from an intrinsically bipolar explosion is also a plausible interpretation of the data but would probably require a ballistic jet of radioactive material capable of penetrating the hydrogen envelope early in the recombination phase. Finally, our latest space-based optical imaging confirms that the late interaction-powered light curve dropped below the stellar progenitor level, confirming the RSG star’s association with the explosion.

Subject headings:
supernovae: general — supernovae: individual (SN 2013ej) — galaxies: individual (NGC 628)
22affiliationtext: e-mail:

1. Introduction

Type II-P supernovae (SNe) are the most common class of core-collapse explosion, marking the deaths of red supergiant (RSG) stars with thick hydrogen envelopes. Direct detection of progenitors indicates that their initial stellar masses lie in the range 8–18 M (Smartt 2009). The namesake plateau light curves of SNe II-P are produced as the photospheric size of the expanding envelope is balanced by the inwardly propagating wave of recombination, releasing the shock-deposited energy over a period of days. Once the envelope fully recombines and the photosphere passes into the inner layers, the light curve drops steeply and the SN enters the nebular phase, the decline of which is governed by the emission of radioactive decay energy. However, the luminosity evolution during both the recombination and nebular phases can be augmented by interaction with circumstellar material (CSM) lost by the progenitor leading up to core collapse.

SN 2013ej is a bright core-collapse supernova (SN) in the nearby galaxy NGC 628 (Messier 74). The stellar progenitor was detected in archival imaging data from the Hubble Space Telescope (HST) (Van Dyk et al. 2013), and subsequently inferred to be an M-type supergiant with an estimated initial mass in the range 9.5–15.5 M (Fraser et al. 2014). The SN exhibited a peak absolute magnitude of , and a nebular-phase decline rate consistent with a small mass of radioactive Ni in the range 0.013–0.023 M (Dhungana et al. 2016; Yuan et al. 2016). Although SN 2013ej was initially classified as a SN II-P (Leonard et al. 2013; Valenti et al. 2014), subsequent studies showed that the relatively fast decline rate during the recombination phase ( mag in the first 50 days after rising up to the plateau; Valenti et al. 2016) and some of the spectroscopic features appear similar to those of SNe from the II-L class (Bose et al. 2015; Huang et al. 2015; Dhungana et al. 2016; Valenti et al. 2016). The distinction between SN II-L and II-P light curves is not always obvious (Anderson et al. 2014; Sanders et al. 2015); indeed, large samples have revealed a continuum of light-curve morphologies intermediate between sources classified as SNe II-P and SNe II-L. Thus, a simple separation criterion for the two classes is probably not valid (Valenti et al. 2016). It is generally clear, however, that both SN types usually stem from cool supergiant progenitors. The detailed evolution of SN II light curves (rise time, peak brightness, and plateau/linear recombination phases) is thought to be a combined function of progenitor envelope mass and radius (possibly reduced by pre-SN mass loss in the case of objects classified as SNe II-L) and explosion energy (Gall et al. 2015) — but it is also likely that the observational characterization of the SN II-P and II-L subclasses could be affected by CSM interaction, which can provide an added source of continuum luminosity, and/or aspherical explosion geometries that are viewed at different angles.

Aspherical supernova geometries are diagnosable with spectropolarimetric observations (see review by Wang & Wheeler 2008). For SNe II-P, in particular, most objects that have been well studied with this technique have exhibited only weak linear polarization (%) during their recombination phases, indicating that the outer envelopes of these explosions tend to be highly spherical (Leonard et al. 2006; Chornock et al. 2010; but also see Dessart & Hillier 2011). However, substantial increases in continuum and line polarization have been observed in many SNe II-P at the end of the plateau stage and into the nebular phase. It has been suggested that as the electron-scattering photosphere recedes into the ejecta, the more aspherical inner layers become exposed. Nonuniform absorption of the photosphere by ions can yield enhanced polarization of the associated spectral lines. A commonly invoked origin for explosion asphericities is an asymmetric distribution of radioactive Ni from the core, which could manifest itself in a variety of geometric configurations, from irregular plumes or clumps to axisymmetric lobes, tori, or jets (Chugai et al. 1992, 2006). It is important to note, however, that optical scattering of SN light by an aspherical distribution of circumstellar dust is another potential means of producing net continuum polarization (the light-echo effect; e.g., Wang & Wheeler 1996). This latter mechanism might be particularly important after the direct light from the recombination phase diminishes and the delayed scattering of that light off of dusty CSM continues.

In addition to spectropolarimetry, late-time spectroscopy can probe for the presence of explosion asphericities as well. If SN ejecta encounter CSM created by the progenitor’s pre-SN mass loss, broad-line emission resulting from the passage of ejecta through the reverse shock can be detected (Chevalier & Fransson 1994; Mauerhan & Smith 2012; Milisavljevic et al. 2012). The signature of an aspherical interaction zone, either from aspherical ejecta or CSM, can become imprinted on the emission-line profiles. However, the causes of emission-line asymmetries are not limited to the geometric configuration of the ionized material; the distribution of newly synthesized or pre-existing circumstellar dust, if present, can result in varying degrees of obscuration for the approaching and receding regions of the outflow and can thus have a substantial effect on the line-profile morphologies as well.

Leonard et al. (2013) were the first to report strong linear polarization from SN 2013ej, measured as early as day 7. Although no correction for the effects of interstellar polarization (ISP) were applied, substantial intrinsic asphericity was nonetheless convincingly indicated by the data; the implied asphericity for the outer hydrogen envelope is unusual in the context of previously studied SNe II-P. More recently, Kumar et al. (2016) presented -band imaging polarimetry of only the post-plateau nebular phase, reporting null polarization beginning at days that climbed to % by day ; however, no data during the earlier recombination phase were obtained in that study. In addition to indications of asphericity from polarimetry data, SN 2013ej exhibited signs of asphericity in its optical flux spectrum during its nebular phase, in the form of an asymmetric red shoulder on the H profile (Bose et al. 2015).

SN 2013ej also exhibited signs of weak CSM interaction in X-ray data. Chakraborti et al. (2016) reported the results of long-term X-ray monitoring observations between days 13 and 145, confirming persistent interaction with CSM from a pre-explosion steady wind having a density and mass-loss rate typical of an RSG (), which must have been blowing for at least 400 yr. Thermal dust emission is another documented feature of SN 2013ej. Tinyanont et al. (2016) present measurements of the SN out to day 438, showing the infrared (IR) decline to be broadly consistent with the behavior of most core-collapse SNe. SN 2013ej also exhibited some peculiar ionization characteristics at early phases. In particular, Si ii emission appeared during the first 2 weeks (first noted by Leonard et al. 2013).

In this paper, we investigate the nature of explosion asphericity, CSM interaction, and thermal dust emission in SN 2013ej. Multi-epoch spectropolarimetry of SN 2013ej during the recombination and early nebular phases is given in §2, and we interpret the data in the context of simple electron-scattering and dust-scattering models. In §3, we present deep late-time optical spectroscopy and ultraviolet (UV) through optical imaging to investigate the persistence of SN emission and CSM interaction. Additional epochs of late-time mid-IR photometry are given in §4 and compared with existing data, and we investigate the evolution of thermal emission from dust. In §5, we discuss potential interpretations of the data and the possible physical interconnections between the asphericity, CSM interaction, and thermal dust emission.

2. Spectropolarimetry

2.1. Kast Observations

Spectropolarimetry was performed at Lick Observatory using the the Shane 3 m reflector and Kast spectrograph (Miller et al. 1988; Miller & Stone 1993). SN 2013ej was observed on the following dates in 2013: Aug. 4, 8, 12, 30; Sep. 6; Oct. 1, 5, 26, Nov. 2, 8 (UTC is used throughout this manuscript). These dates correspond to epochs between 11.5 and 107.2 days, assuming the explosion occurred on July 23.95 (JD = 2,456,497.45; Valenti et al. 2014). Kast is a dual-beam spectropolarimeter that utilizes a rotatable semiachromatic half-waveplate to modulate the incident polarization and a Wollaston prism in the collimated beam to separate the two orthogonally polarized spectra onto the CCD detector. Only the red channel of Kast was used for spectropolarimetry; a GG455 order-blocking filter suppressed all second-order light at wavelengths shortward of 9000 Å.

Observations were made with the 300 line mm grating and a 3-wide slit, yielding a spectral resolution of  Å pixel and a full width at half-maximum intensity (FWHM) spectral resolution of  Å (as determined by the FWHM of calibration-lamp lines). The useful wavelength coverage with this setup is 4600–9000 Å. The orientation of the slit on the sky was always set to a position angle of (i.e., aligned north-south). In most cases the source was observed through transit of the meridian, and thus tracked through parallactic angle. However, some observations were performed with the slit at a position angle far enough from the parallactic angle that we experience a small amount of blue light loss. This does not affect the spectropolarimetry, owing to the dual-beam instrumentation, but does slightly affect the individual flux spectra we obtain from the data. Exposures of 900 s were obtained at each of four waveplate positions (, , , and ). In most cases, three waveplate sequences were performed. Flatfield and calibration-lamp spectra were obtained immediately after each sequence, without moving the telescope. Data were extracted and calibrated using generic IRAF routines and our own IDL functions. Spectropolarimetric analysis was also performed in IRAF and IDL following the methods described by Miller et al. (1988) and implemented by Leonard et al. (2001). All spectral data of the SN have been shifted to the rest-frame wavelength, assuming a host-galaxy redshift of .

For polarimetric calibrations, standard stars were selected from the sample of Schmidt et al. (1992a,b). Unpolarized standard stars HD 212311, BD32 3739, and HD 57702 were observed to verify the low instrumental polarization of the Kast spectrograph. We constrained the average fractional Stokes and values to %. By observing the unpolarized standard stars through a 100% polarizing filter, we determined that the polarimetric response is so close to 100% that no correction was necessary, and we obtained the instrumental polarization position-angle curve, which we used to correct the data. We observed the high-polarization stars HD 204827, HD 19820, and BD59 389 to obtain the zeropoint of the polarization position angle on the sky () and to determine the accuracy of polarimetric measurements, which were generally consistent with previously published values to within % and .

2.2. Polarimetry Formalism

Linear polarization is expressed as the quadratic sum of the and Stokes vectors, . The position angle on the sky is given by , taking into account the quadrant in the plane where the inverse tangent angle is located. Since is a positive-definite quantity, it is overestimated in situations where the signal-to-noise ratio (S/N) is low. It is thus typical to express the “debiased” (or bias-corrected) form of as , where and are the uncertainties in the and Stokes parameters. Note, however, that at low S/N, is also not a reliable function, as it has a peculiar probability distribution (Miller et al. 1988). Thus, for extracting statistically reliable values of polarization within a particular waveband, we have binned the calibrated and Stokes spectra separately over the wavelength range of interest before calculating and . All quoted and tabulated values in this paper were determined in this manner, while some figures present . For , if , then we set a 1 upper limit on of . In cases where , is essentially undetermined and is not graphically displayed.

Figure 1.— Spectral region near Na D interstellar absorption feature for SN 2013ej and SN 2002ap, both of which occurred in M 74.
UT Date MJD Epoch  (%)  (deg)  (%)  (deg)  (%)  (deg)
Observed Values
Aug. 04 56509.0 11.5
Aug. 08 56513.0 15.5
Aug. 12 56517.0 19.5
Aug. 30 56535.0 37.5
Sep. 6 56541.9 44.4
Oct. 1 56567.0 69.5
Oct. 5 56570.8 73.3
Oct. 26 56591.8 94.3
Nov. 2 56598.8 101.3
Nov. 8 56604.7 107.2
ISP-Subtracted Values
Aug. 04 56509.0 11.5
Aug. 08 56513.0 15.5
Aug. 12 56517.0 19.5
Aug. 30 56535.0 37.5
Sep. 6 56541.9 44.4
Oct. 1 56567.0 69.5
Oct. 5 56570.8 73.3
Oct. 26 56591.8 94.3
Nov. 2 56598.8 101.3
Nov. 8 56604.7 107.2

Uncertainties are statistical.

Epoch is given in days since the adopted explosion date (JD 2,456,497.45; Valenti et al. 2014).

-band values averaged over the wavelength range 5050–5950 Å.

-band values averaged over the wavelength range 5890–7270  Å.

Continuum sample region is 7800–8150 Å; some weak line contamination might be present.

% and values were adopted from Leonard et al. (2002).

Table 1Integrated Broadband and Continuum Polarization of SN 2013ej.

2.3. Interstellar Polarization (ISP)

As light traverses interstellar media, non-spherical dust grains aligned with the local magnetic field will preferentially scatter incident photons with electric vectors that are aligned with the grains (Hiltner 1949; Davis 1955). This results in the dichroic extinction of light, which polarizes the transmitted signal. The resulting polarization of an object is thus the vector sum of the interstellar and intrinsic source Stokes parameters. Therefore, to obtain a reliable measurement of the intrinsic polarization of an extragalactic source, the interstellar component of the polarization must be determined and subtracted from the observed Stokes vectors. This is often a very problematic element of polarimetric analysis, as an accurate assessment of the ISP from both the Milky Way (MW) and, in particular, the SN host galaxy, can be very difficult to obtain. Foreground stars in the MW can be used to gauge only the Galactic component of ISP.

Prior to SN 2013ej, M 74 was host to SN 2002ap, for which spectropolarimetry was obtained by Leonard et al. (2002), Kawabata et al. (2002), and Wang et al. (2003). Considering the results of Kawabata et al. (2002) along with their own, Leonard et al. (2002) concluded a MW-dominated ISP value of % and % (%, ) in the direction of SN 2002ap. Wang et al. (2003) arrived at a similar value of % and % (%, ) but concluded that the host ISP is stronger than the MW. To determine whether the ISP estimates for SN 2002ap can be reasonably applied to SN 2013ej, we must first confirm that the ISP is indeed MW dominated. Figure 1 compares the spectra of both SNe in the spectral region around the Na i D interstellar absorption lines (spectra reproduced from Foley et al. 2003 and Dhungana et al. 2016). No significant absorption features are detected at the redshift of the host, while the lines from the MW have comparable strengths. Since the extinction is clearly dominated by the MW, one could reasonably conclude that the ISP is probably dominated by the MW as well. We therefore adopt values % and from Leonard et al. (2002), and applied them to our spectropolarimetric data according to the functional form of Serkowski et al. (1975) and Whittet et al. (1992).

However, in their recent broadband imaging polarimetry study of SN 2013ej, Kumar et al. (2016) performed an analysis of foreground ISP-probe stars in the broad vicinity of SN 2013ej, while also considering the ISP values derived by Leonard et al. (2002), Kawabata et al. (2002), and Wang et al. (2003) for the case of SN 2002ap. Narrowing the sample to four ISP probe stars within 5 of the SN, and after removing those which exhibited evidence for polarimetric variability, Kumar et al. (2016) derived with in the band for SN 2013ej. If this value is scaled from the band to the entire spectrum according to Serkowski et al. (1975), then the resulting -band Stokes parameters are % and % (). This value of is only slightly higher than the ISP estimates for SN 2002ap discussed above, and at a position angle that differs by ; but it is also 0.1% larger than the maximum polarization expected from the limiting relation suggested by Serkowski et al., where , assuming  mag in the direction of M74 (Schlafly & Finkbeiner 2011). We therefore consider this to be an overestimate, and instead we move forward adopting the ISP values derived by Leonard et al. (2002).

A previous employed method to check the validity of ISP estimates is to assume that the peaks of the strongest emission-line features in the flux spectrum, such as H, are intrinsically unpolarized (e.g., see Tran et al. 1997); this method can be particularly useful when the host-galaxy ISP is also significant (note that this does not appear to be the case for SN 2013ej). The justification for this assumption is that strong emission lines form in the outer regions of the ejecta, well above the polarized continuum photons emerging from the electron-scattering photosphere. However, this assumption might be inappropriate; although substantial depressions in polarization are indeed normally observed at the wavelengths of broad emission components, radiative transfer calculations by Dessart & Hillier (2011) have demonstrated that variations in optical depth can result in non-null polarization of the H emission feature. Furthermore, weak features of He that overlap with broad H might also contribute polarized flux at these wavelengths (Maund et al. 2007). Moreover, the assumption of unpolarized emission-line peaks assumes that electron scattering is the sole source of the polarized flux. For example, if dust scattering in CSM is an important contributor to the net polarization, then the use of this method will obviously result in erroneous estimate of ISP. As we will demonstrate in the following sections, the ISP-subtracted data indeed show a changing amount of residual polarization remaining at the wavelength of the H flux peak for all epochs, which could be indicative of CSM dust scattering.

Figure 2.— Debiased polarization () and position angle () for days 11 through 107. The flux spectra are plotted in red for reference. Prominent polarized line features are marked (note that polarized peaks are always blueshifted with respect to the central wavelengths of the associated transitions). Wavelengths for which have been omitted from the curve.

2.4. Polarization of SN 2013ej

The observed and ISP-subtracted values of integrated and for SN 2013ej are presented in Table 1, and the polarized spectra are presented in Figure 2 (flux spectra are also included for reference). Hereafter all reference to the polarization data will imply the ISP-subtracted values. Spectral line identifications are adopted from the literature — in particular, from Bose et al. (2015) and Dhungana et al. (2016), whose line identifications are supported by their spectral modeling. Those line identifications are also corroborated by other spectroscopic studies of SN 2013ej (e.g., Valenti et al. 2014; Huang et al. 2015; Yuan et a. 2016). We refer the reader to these works for further detailed justification of line identifications.

Figure 3.— Left: Temporal evolution of the integrated broadband and for SN 2013ej, after subtracting our favored value for the ISP (see §2.3). The red, green, and black dots represent the R-band, V-band, and continuum regions, respectively. A scaled version of the light curve is shown in gray color for comparison. Right: Same data after subtracting the slightly different value of ISP adopted from Kumar et al. (2016), whose -band imaging polarimetry measurements of SN 2013ej are also shown for comparison (orange triangles).

From day 11, we observe substantial polarization of the continuum at a level of %. We chose the wavelength range of 7800–8150 Å to sample the continuum because this region remains mostly devoid of substantial line features throughout the evolution of the SN. There are modulations in and across H, H, and slightly for He/Na around 5875 Å. The level of polarization near the peak of the H emission component of the P-Cygni profile is at %, while rotates by between the blueshifted and the redshifted sides of the emission feature, straying from the continuum position angle of the continuum on both sides. The data on day 15 appear similar, but the modulation across H dips near null levels of and the rotation across the H feature increase to .

Figure 4.— Temporal evolution of the H (upper left), Na i-D (upper right), Ca ii IR (lower left), and O i (lower right) line polarization on days 16 (gray), 38 (black), 74 (magenta), and 94 (blue). For each plot, the upper-left panel shows the flux spectrum, the lower-left panel displays the polarization, and the right three panels show the behavior in the plane. The filled dots in the plane mark the integrated continuum polarization for that date; the day 16 data are shown in all panels for reference. For Na i, if the feature is instead associated with He i, then the inferred velocities are  km s less.

By day 19, enhanced polarization across the absorption components of line features starts becoming prominent. The strongest line feature at this stage appears to be associated with H, peaking near %. Weaker enhancements are seen for the Ca ii IR triplet and O i near 7650 Å. There is an absorption line just blueward of the H component, together appearing as a double-dipped feature. The bluer feature has been identified by Valenti et al. (2014) as Si ii 6355. We observe a slight enhancement in by % for this feature, relative to the polarization level on either side. Leonard et al. (2013) noted that the appearance of this potential Si ii feature at such an early phase is unusual, suggesting peculiar ionization conditions not standard for an RSG envelope; this issue will be revisited in §5. Leonard et al. (2013) also considered the alternative possibility that this feature is a separate high-velocity component of H; however, there is no comparable high-velocity feature associated with H. Therefore, we favor the Si ii identification, which is supported by the spectral modeling of Bose et al. (2015) and Dhungana et al. (2016), and is also the consensus among other studies of SN 2013ej in the literature (Valenti et al. 2014; Huang et al. 2015; Yuan et al. 2016).

By day 38, the line features in the flux spectrum have strengthened, with substantial deepening of their absorption components. H absorption begins exhibiting enhanced polarization, with clear depolarization across the emission component. Ca ii IR polarization dominates the spectrum, peaking near %, with complete depolarization of its emission component. Interestingly, it is at this epoch that the Ca ii absorption develops multi-component structure. O i is also enhanced at %. Less prominent line features of Na D 5890, 5896, the Sc ii multiplet near 5700 Å, and Fe ii 5535, 5169 also exhibit enhanced polarization. The overall characteristics remain similar on day 44, but with the Fe ii polarization climbing above 2.0%.

Some interesting changes occurred by day 70. Both the line and continuum polarization have further increased by a few tenths of a percent, with strong, broad modulations in polarization becoming prominent, particularly in the region between 5800 Å and 6500 Å. The global position angle also appears to have rotated slightly, as illustrated by Figure 3. Just 3 days later at day 73, highly polarized line features grow narrower in width. He i features become discernible across much of the optical spectrum. H line polarization is particularly striking at this epoch, with peak polarization exceeding % and a very narrow morphology; this strongly peaked feature appears to be coincident with the emergence of a asymmetric notch on the blue side of the absorption component, which will be shown in more detail below in Figure 4. It is remarkable that the large modulations in polarization, particularly in the wavelength range 6000–6500 Å, exhibit such constancy in position angle.

By day 94, H and H polarization have diminished in strength, while line features associated with Ba ii 6142 and the Sc ii multiplet near 5700 Å have increased to almost . The continuum polarization at this epoch has increased to its peak level for our coverage at %; however, we note that our chosen continuum sample region also appears to contain a small but non-negligible contribution from line polarization. At this epoch the H profile develops a blue bump, which is also shown in more detail in Figure 4.

On day 101 the continuum polarization has dropped back down to %, while only minor changes in line polarization are seen. Just 6 days later, on day 107, the Ba ii and Sc ii features around 6100–6200 Å are no longer discernible, but the continuum polarization has only marginally decreased to %. The line-blanketed region in the blue near  Å appears to drop to null values of polarization. H and the Sc ii multiplet near 5700 Å remain prominent, however.

Figure 3 (left panel) illustrates the temporal evolution of the integrated broadband polarization and the position angle. During the first 20 days there is a % fluctuation in at nearly constant . Subsequently, and through day 44, the continuum and -band remain within a range of 0.8% to 1.1%, while the band exhibits higher values in the range 1.2% and 1.3%. By day 70, the values of all bands approach convergence near , but quickly drop again all together by several tenths of a percent only a few days later; also makes a substantial shift at this point and continues to drift thereafter. Specifically, on day 94, when the light curve is near the edge of the plateau, the continuum and reach their largest excursions from the average, shifting by much more than the or bands. This implies that there might be a changing source continuum that is separate from the line . However, never deviates by more than at any point. As the light-curve plateau drops down to the radioactive decay tail, the continuum also drops by several tenths of a percent, while shifts back toward the values measured on day 73 and earlier.

Figure 3 (right panel) shows the same data after subtracting the slightly larger value of ISP derived by Kumar et al. (2016). As described in §2.3, we adopted their -band ISP value and scaled to the -band and continuum regions according to the functional form of Serkowski et al. (1975) and Whittet et al. (1992). The -band imaging polarimetry measurements from Kumar et al. (2016) are also included in the Figure 3 for comparison. We note that their first two epochs of polarimetry near day 100 are significantly below the persistently high polarization trend exhibited by our data, although the error bars associated with their second epoch overlap with our inferred trend. Their third epoch measurement is consistent with the trend exhibited by our data, but their associated value appears offset by . Their following two measurements past day 110, during the nebular phase, exhibit a significant rise in polarization beyond the temporal range of our coverage. Given the relatively large uncertainties in their measurements, only their first epoch might appear to be inconsistent with the trend inferred from our data. However, the apparent evolution that one might glean from their data alone — a weakly polarized recombination phase exhibiting a continuous increase in polarization after the plateau, similar to some previous studies of SNe II-P (e.g., Leonard et al. 2006; Chornock et al. 2010) — would be inconsistent with the persistently strong polarization we have measured for SN 2013ej throughout the entirety of the recombination phase.

2.4.1 Behavior in the Plane

The behavior of the most prominent polarized line features in the plane, and their strength as a function of radial velocity, are illustrated in Figure 4 for days 16, 38, 73, and 94. On day 16, H forms a cluster of points exhibiting slight elongation along the axis, with maximum deviation of % occurring at a radial velocity near  km s. On day 38 the feature is stronger (%) and clearly elongated along the symmetry axis defined by the continuum region. Peak polarization is near  km s. By day 73 the feature becomes a tight linear feature and is still aligned with the continuum symmetry axis. Strengthened polarization near 2.3% is seen out to velocities of  km s, but a prominent (%) narrow peak has developed at 8000 km s. This is the epoch where the high-velocity notch first mentioned in the previous section has appeared in the absorption component. By day 94, the strong feature has weakened only slightly and remains near the same position angle, but has broadened to lower radial velocities of 5000 km s, where the position angle shifts slightly, following the small drift of the continuum axis. At this epoch a bump has emerged on the blue side of the H profile.

For Ca ii IR, the behavior in the plane is quite different from that of H. On day 38, strongly enhanced polarization forms a relatively wide loop in the Stokes plane; maximum polarization still occurs near the position angle of the continuum symmetry axis, but polarization at lower velocities exhibits a departure from axial symmetry on the sky ( excursion in the Stokes plane). The maximum polarization of Ca ii also occurs at significantly higher radial velocities than H, peaking near  km s. By day 73, maximum polarization remains near the same axisymmetric angle, but the position angle at lower velocities exhibits more erratic behavior in the plane. But by day 94, the Ca ii feature forms a more elongated structure in the plane that extends into the same position angle as that of H and the other lines in Figure 4, but at significantly higher velocities.

For the Na feature, maximum polarization occurs on day 73, as it does for H. The line also forms a linear structure in the plane at this epoch. The position angle is comparable to that of H, but is slightly shifted by about on the sky. The velocities associated with the enhanced polarization are significantly lower than in the case of H, peaking near  km s (assuming association with Na i), and they likely originate at deeper layers in the homologous outflow.

O i appears to follow an evolutionary trend similar to that of other lines, and at velocities comparable to those of Na i. However, by day 94 the development of multiple absorption components in the flux spectrum suggests that contamination by other lines might be an issue for this later epoch.

2.5. Spectropolarimetric Modeling

2.5.1 Dust-Scattering Model

SN 2013ej exhibited X-ray evidence for interaction with tenuous CSM (Chakraborti et al. 2016), the density of which is consistent with the mass loss of a typical RSG wind. Moreover, as we will show in §3, CSM interaction persists at late times, implying the presence of a highly extended tenuous CSM envelope. RSG winds have long been known to contain dust (Zuckerman 1980), and high-resolution imaging and polarimetric observations of Betelgeuse and Antares have shown that such dust can be clumpy (Kervella et al. 2011, 2016; Ohnaka 2014); this material could provide an asymmetric scattering medium and, if viewed at a suitable orientation, a major source of polarization. Depending upon the composition and size of dust grains and their distance from the SN, circumstellar dust can potentially survive if the UV pulse of SN shock breakout fades faster than the sublimation timescale at a given radius. The polarizing effect of dust scattering should therefore be considered in the case of SN 2013ej.

Prior to Wang & Wheeler (1996), strong, broad modulations in polarization across line features were taken as evidence against circumstellar scattering; those authors showed that interpretation to be incorrect, following their re-analysis of the polarization data from SN 1987A. Indeed, the observed polarization depends on both the intrinsic SN light curve and the scattering process in a complicated way (see Wang & Wheeler 1996, their Equation 3). The net signal is a mix of direct light from the SN, with highly polarized scattered light from an earlier spectroscopic phase, echoed by distant CSM dust. In short, the net polarization develops wavelength dependence because the luminosity evolves differently for different wavelengths; this is especially true for P-Cygni features. Therefore, polarization modulations across lines are not exclusive to electron scattering in the SN envelope.

Figure 5.— The polarization spectrum from dust-scattering spectropolarimetry models (purple) compared to the data (black) at the observed epochs. We refer the reader to Figure 2 for line identifications.

We therefore constructed a model for which the polarized signal is the result of optical scattering by dust particles in the CSM, building upon the methodology of Wang & Wheeler (1996). The detected photons consist of those propagating directly from the SN mixed with those scattered toward the observer by the circumstellar dust. This mode of polarization does not require intrinsic asphericity of the SN. The basic elements of the model, which will be presented in more detail in a forthcoming paper (M. Hu et al., in preparation), are summarized as follows. The available flux spectra were used as the input source of light, and the luminosity evolution was calibrated using the observed light curve (for this we used all of the spectra and photometry originally presented by Huang et al. 2015). The asymmetric component of CSM was modeled as a single stationary sphere of dust with radius  cm, displaced from the SN by  cm in the plane of the sky ( day light-travel time from the SN), resulting in a scattering angle for the observer; this orientation is a reasonable approximation, since photons that scatter at right angles will dominate the polarized signal. The models were constructed by trial-and-error to provide the best overall visual matches to the multi-epoch data. In the end, the separation of  cm for the dust sphere was used because it provided the best match to the observed spectropolarimetric evolution. Model calculations were performed and are presented for , 73, 101, and 107 days. The optical depth to dust scattering is allowed to vary from epoch to epoch, and the best-matching models were obtained for optical depths that decreased continually from values of 0.1 to 0.05 between 44 and 107 days, respectively. The model does not distinguish between graphite and silicate grains, although we note that graphite grains are much more likely to survive the SN at radii within  cm for grain sizes m (see Wang & Wheeler 1996, their Figure 1).

The synthetic polarized spectra from our dust-scattering models are compared with the data in Figure 5. Considering its geometric simplicity, the model reproduces many of the observed features in the polarized spectrum reasonably well. For day 44, the model generates substantial continuum polarization and strongly polarized line features. In particular, Ca ii IR line polarization is closely matched in terms of strength and radial velocity. The model H line exhibits similar morphology and blueshift relative to the data, but does not fully reproduce the depth of the polarization trough associated with the emission component. The Fe ii feature near wavelength 5000 Å and H at the blue end of the spectrum are generated by the model, although at significantly higher strengths than we observed. The weak O i 7774 is also fairly well reproduced on this date in terms of strength and radial velocity.

On days 73 and 101, the observed 1% continuum level is well matched. However, the model fails to reproduce the polarized line profile of Ca ii IR, showing a relatively narrow trough in polarization where an enhancement is observed. The model does generate line features that correspond with most of those that are present in the data; however, the features are not matched well in strength. A particularly noteworthy aspect of the modeling at these epochs is the relative blueshifts of polarized line features. While the model Na i D, Sc ii, Fe ii, and Ba ii features match the data in radial velocity, the H, O i, and Ca ii features in the data exhibit substantially higher radial velocities. We speculate that this discrepancy might provide a means to distinguish between multiple components of line polarization (i.e., those intrinsic to the SN and those produced by scattering in CSM).

By day 107, the model match improves for some of the broad modulations across the polarized spectrum, while worsening for some of the narrow line features. Specifically, the depth of the broad H, O i, and Ca ii appears well matched, while the model exhibits a strong Na i D polarization feature that is no longer present in the data at this epoch. The model also predicts the vanishing of the narrow H feature altogether, yet this feature remains strong in the data on day 107 and thus appears to be better matched by the day 101 model. The blue end of the model spectrum exhibits significantly higher polarization than the data, although the downturning trend in the spectrum toward shorter wavelengths is reproduced.

Our basic scattering model has substantial limitations in reproducing some of the detailed features and evolution. Perhaps most importantly, the single dust-sphere configuration results in no scattered photons being observed until after the associated light-travel delay time of days, yet strong polarization is indeed observed from the earliest epoch on day 11. A more realistic model might involve multiple clumps in a more continuous distribution of circumstellar dust. Furthermore, pure dust scattering will always produce some polarized flux at the wavelengths of the strong emission peaks in the flux spectrum. In this regard, the model best matches the data from day 107, after the recombination phase has ended, and, to a lesser extent day 44. However, on day 73, in particular, we see that H and Ca ii emission exhibit nearly complete depolarization, which suggests that electron scattering is likely to be the dominant mode of scattering on these epochs.

2.5.2 Electron-Scattering Model

We experimented with the construction of electron-scattering models to match the flux and polarization spectra from day 73, which is the epoch where the line polarization is most pronounced and is not satisfactorily reproduced by the dust-scattering models presented above. This epoch also exhibits nearly complete depolarization of the H and Ca ii emission peaks, which is another indication that electron scattering denominates the polarization. We used the program synpol, which was utilized for the spectropolarimetric analysis of the Type Ia SN 2001el in the spectral vicinity of the Ca ii IR feature (see Kasen et al. 2003 for details); synpol is an extension of the program synow (e.g., Branch et al. 2007), which is used for modeling flux spectra of SNe. For polarization modeling, synpol allows for an ellipsoidal photosphere, geometrically represented as an inner unpolarized boundary surface surrounded by a pure electron-scattering envelope having a power-law electron-density profile of the form . The program synpol does not satisfactorily reproduce the strong, broad emission features in the flux spectra of SNe, because, like synow, the code assumes local thermodynamic equilibrium (LTE) and thus does not account for collisional excitation, which is important in the line-forming regions of SNe. Consequently, synpol is not able to fully reproduce the broad dips in polarization associated with the emission components of the flux spectrum, which result from the fact that the line-forming region is optically thin to electron scattering and adds a source of mostly unpolarized flux to the net signal.

Figure 6.— Our best-matching synpol model for the normalized flux (upper panel) and polarization spectrum (lower panel) of SN 2013ej on day 73. We refer the reader to Figure 2 for line identifications.

Uneven absorption of polarized photospheric light by ions in the overlying line-forming region can produce enhanced polarization features at the wavelengths of absorption components in the flux spectrum. From the observer’s perspective, the fastest radial velocity region on the approaching side of the homologous outflow will preferentially obscure the central portion of the underlying photosphere. Since photospheric light emerging from that central region is mostly forward scattered, it does not contribute substantially to the net polarization of the unresolved source. The blueshifted absorption thus removes the weakly polarized light from the net signal, while strongly polarized light from the photospheric limb (where radially diffusing photons will be scattered at angles near ) emerges. This effect thus increases the fractional polarization of the blueshifted side of the absorption profile, blueward of the flux minimum. This mode of line polarization will produce line polarization exhibiting a high degree of axisymmetry, and will produce linear features in the plane, such as those we see for H on day 73 (see Figure 4); however, this mode requires that the underlying photosphere is aspherical, in order to prevent net cancellation of all electric vectors emerging from the photospheric limb photons. The alternative case of patchy line absorption by clumps of ions over a spherical photosphere can also produce polarized lines, but in this case the features will be non-axisymmetric, and will produce circular loops in the plane — for example, more consistent with the relatively high-velocity Ca ii feature on day 38 in Figure 4.

The synpol program allows for various geometric configurations for the line-forming region, including ellipsoids, spherical clumps, and toroids. A homologous flow is assumed, whereby the radial extents of the photosphere and line-forming regions are defined by velocity coordinates, set by the user. For clumps of line absorption randomly distributed across an aspherical photosphere, one would expect to observe substantial modulations in across the wavelength range of the spectral absorption features, since both axisymmetric and non-axisymmetric components of the photosphere will be obscured over the associated range of radial velocities. For similar reasons, a toroidal distribution of line absorption will also exhibit substantial shifts in with respect to radial velocity, particularly if the toroid is viewed nearly edge-on. For SN 2013ej, the observed high degree of axisymmetry between the polarized absorption features and the continuum appears be more consistent with an ellipsoidal geometry that matches the ellipsoidal geometry of the underlying photosphere.

Figure 7.— Left panel: late-time Keck/LRIS spectra of SN 2013ej (black) from days 434 and 807. A comparison spectrum of the Type II-P SN 2004et (blue; Sahu et al. 2006) on day 428 is also shown in the left panel for comparison, which was obtained from the WISeREP database (Yaron & Gal-Yam 2012). Right panel: an expanded view of the region near H and O i. The boxy profile on day 807 appears to have been present on day 434 as well, underneath the strong nebular features.

We therefore constructed our synthetic polarization spectra using an ellipsoidal photosphere and line-forming region. The following atomic transitions were included in the model: H i, O i, Na,i, Mg i, Si ii, Ca ii, Fe ii, Ti i, Ti ii, and Ba ii. The radii and optical-depth profiles of the line-forming envelopes were also chosen by trial and error until a satisfactory match to the polarized features was obtained. The power-law exponent, , of the radial density profile, , of the line-forming region was set to have for all transitions. However, reproducing the strongly peaked polarized line profile of H required . The best matches to the flux and polarization spectrum were obtained with blackbody and excitation temperatures of 5500 K and a photospheric radius of 4100 km s (velocity coordinates). We note that matching the strongly peaked H line required setting the line optical depth, , to a value of 150, whereas most of the other line matches were obtained with optical depths in the range 5–10, with the exception of Ca ii, which required to get the best match.

Our synpol model for the flux and polarized spectrum is presented in Figure 6. The best overall match was found using an ellipsoidal photosphere with an axis ratio of 1.4, viewed at 90 inclination angle (edge-on). Slightly higher values of axis ratio also produced consistent matches if the inclination was slightly decreased. However, our best matches occurred for axis ratios in the range 1.3–1.5 at inclination angles within of edge-on. The strongly peaked polarized line features were matched reasonably well using line-forming regions that are attached to the photosphere (radii beginning at velocity coordinate 4100 km s) and with the same ellipsoidal geometry (spherical line-forming regions worked similarly well), with the exception of Ca ii IR, for which we obtained our best match using a detached line-forming region with inner radius at 5400 km s; this necessary detachment might be related to the non-axisymmetric geometry of this feature exhibited in Figure 4.

There are large discrepancies between the data and model emission components of the flux spectrum. This is to be expected, since, as stated earlier, our calculations assume LTE and do not account for collisional ionization, which contributes to the emission-line flux. The depolarizing effect of collisional ionization also results in deeper polarization troughs than our model produces, which if taken into account would probably improve the match.

The ellipsoidal model produces purely axisymmetric line polarization, roughly consistent with the high degree of axisymmetry we measure (i.e., the drift in ) across most of the polarized line. Realistically, small deviations from a perfectly aligned photosphere and line-forming region could likely account for the minor spread in observed across the spectrum. However, the erratic structure of Ca ii in the plane (see Figure 4) indicates non-axisymmetric asphericity for the absorption associated with that ion, which our preliminary model does not explain.

Overall, the synpol model provides a reasonable match to the day 73 data, and thus demonstrates that electron scattering in an ellipsoidal geometry viewed nearly edge-on could explain the continuum polarization, and the morphologies of most of the narrowly peaked polarized line features and their associated radial velocities. However, the ability for both electron-scattering and dust-scattering models to generate many of the general spectropolarimetric features of SN 2013ej indicates that both processes could potentially contribute to the observed characteristics. If the position angles of dust- and electron-scattering components are similar, the constructive interference of their associated Stokes vectors could lead to an overestimate of the axis ratio of the electron-scattering photosphere.

Figure 8.— Radial-velocity profiles of selected emission lines in the Keck/LRIS spectrum of SN 2013ej from day 807.
Figure 9.— Keck/DEIMOS spectrum of SN 2013ej from day 897 (black). The peak of line polarization from day 73 (red) is included to highlight the similar velocities of the polarization feature and the edge of the late-time flux profile. Arrows mark regions of interest, including possible discrete individual emission components within the multi-peaked flux profile.

3. Deep Late-Time Spectroscopy

Late-time optical spectroscopy of SN 2013ej was obtained on 2014 Oct. 2 (day 434) and 2016 Jan. 8 (day 897) at Keck Observatory using the Deep Imaging Multi-Object Spectrograph (DEIMOS; Faber et al. 2003), and on 2015 Oct. 10 (day 807) using the Low-Resolution Imaging Spectrometer (LRIS; Oke et al. 1995). The day 434 DEIMOS data were obtained with the 600ZD grating and 10 slit for a FWHM resolution of  Å; two exposures were obtained for a total integration time of 1200 s. The day 807 LRIS observations were performed using the 600/400 grism for the blue channel and the 400/8500 grating for the red channel; three exposures were obtained for a total integration time of 3600 s . The day 897 DEIMOS data were taken with the 1200 l mm grating and 08 slit for a FWHM resolution of 0.9–1.3 Å; five exposures were obtained for a total integration time of 7500 s. The DEIMOS and LRIS observations at all epochs were performed at an airmass of 1.0–1.1 and with the slit aligned with the parallactic angle (Filippenko 1982). All data were flat-fielded and wavelength calibrated using standard techniques. Flux calibration was applied using observations of standard stars at similar airmass, and telluric lines were removed.

The spectra from days 434 and 807 are presented in Figure 7. On day 434, SN 2013ej exhibits a late nebular-phase emission spectrum, dominated by strong lines of H, [O i], [O ii], [Ca ii], [Fe ii], Na i, and Ca ii IR. The spectrum closely resembles that of SN 2004et on day 428 of its evolution (Sahu et al. 2006), and is characteristic of a nebular-phase spectrum from a core-collapse SN.

By day 807 the spectrum of the SN has changed dramatically. H line emission exhibits a very broad and boxy profile, much broader than the nebular spectrum from a year prior, with a bumpy multi-peaked structure. There is also a broad asymmetric emission line of [O I] 6300, 6364 which overlaps with the blue side of H, in addition to a blend of [O ii] 7325 and [Ca ii] 7291, 7324, and emission from [Fe ii] 7155. Closer inspection of the earlier day 434 spectrum, facilitated by the right-hand panel of Figure 7, shows that the same broad, weak components of H and [O i] from day 807 were already present on day 434, underlying the stronger, narrower nebular line profiles. There is weak emission likely to be associated with Na D or He i 5876. Blueward of that feature the spectrum forms a pseudocontinuum resulting from a blend of multiple emission lines, including [Fe ii] transitions and prominent [O iii] 5007.

Figure 10.— Day 807 spectrum of SN 2013ej (black) compared to late-time spectra of SN 1993J (red; Matheson et al. 2000), SN 2004et (blue; Faran et al. 2014), and SN 2007od (Andrews et al. 2010). The spectra have been normalized by their average values (i.e., divided by a constant). The spectrum of SN 2007od was further scaled and offset by a constant for clarity.

All of the broad emission lines exhibit strongly blueshifted peaks. Figure 8 shows an expanded view of H, [O iii], [O i], and [Ca ii] plotted as a function of radial velocity (day 807). The H profile exhibits steep edges extending out to  km s on the blue side and  km s on the red side. There appears to be preferential suppression of flux on the red side of the profile, with the exception of enhanced emission at the far-red edge. The same overall line morphology is apparent for [O iii] 5007, but with stronger red-flux suppression. For this line there also appears to be a strong emission component on the red side near  km s, although this might be emission from a different line. A very narrow component of nebular [O iii] emission is seen at rest velocity, and is probably from the host H ii region. [O i] and [Ca ii] also exhibit line morphologies with blueshifted peaks, although the velocity widths appear to be slightly smaller than those of H and [O iii].

The higher-resolution DEIMOS spectrum of H from day 897 is shown in Figure 9, along with the polarized line profile on day 73, when the line polarization has reached maximum. The flux-profile morphology is roughly consistent with the LRIS spectrum from day 807, but the higher-resolution data exhibits more pronounced detail. In particular, there is possible indication of a triple-peaked structure, including a central intermediate-width (FWHM km s) component near rest velocity, bounded by two broader components. Interestingly, the blue edge of H extends to the same high value of radial velocity as the peak of the polarized line profile from day 73, which suggests that the polarized line feature is an intrinsic feature of SN originating in the outer ejecta; we discuss this further in §5.

Figure 11.— Keck/LRIS -band image (black and white background) of SN 2013ej and M74 on 2015 Feb. 9. The inset shows HST data from 2015 Dec. 29 used to extract the photometry presented in Table 2; the color images are a blue-green-red composite constructed from the respective F275W, F336W, F438W (left inset) and F438W, F555W, F814W images (right inset). North is up and east toward the left in all images.

The overall appearance of the day 807 spectrum resembles that of the Type IIb SN 1993J at day 975, shown in Figure 10. H and Na i/He i emission, with respect to the other spectral features, are much stronger in SN 1993J, but the widths of the lines are roughly the same and both exhibit red suppression. Their resemblance is particularly striking in the blue pseudocontinuum, including the broad [O iii]/H and blended [Fe ii] features. SN 2013ej also shares similarities with the Type II-P SN 2004et (day 819), although the latter does not exhibit the strong emission features in its blue continuum and the H profile is more asymmetric and less boxy in shape. The relative strength of the blueshifted component of emission from [O i] and [O i][Ca ii] near 7200 Å is most pronounced in SN 2013ej.

The broad, boxy line profiles observed on days 807 and 897 are the result of CSM interaction, generated as the SN ejecta become energized by passage of the reverse shock. The shocked ejecta experience collisional excitation while pre-shock ejecta can be photoionized by X-rays and UV emission from the hot plasma near the reverse shock (Chevalier & Fransson 1994). Both processes potentially contribute to the broad-line emission, although post-shock collisional excitation is probably the dominant contributor to H.

CSM interaction through day 897 probes mass loss from the progenitor going back at least  yr before core collapse, assuming shock and CSM-wind velocities of 9700 km s and 10 km s, respectively (Chakraborti et al. 2016). Meanwhile, the earlier presence of broad emission on day 434, seen underneath the stronger nebular line profiles, indicates that CSM interaction was occurring at this phase as well. This is consistent with long-term X-ray observations, which showed that SN 2013ej was indeed persistently interacting with the RSG wind for at least 145 days (Chakraborti et al. 2016). Such highly extended, tenuous, and asymmetric CSM is reminiscent of the environments surrounding the RSGs Antares and Betelgeuse (Kervella et al. 2011, 2016; Ohnaka et al. 2014), and also motivates our construction of dusty-CSM scattering models to investigation the spectropolarimetry (see §2).

4. Late-Time Ultraviolet–Infrared Photometry

4.1. Hst/wfc3

Band mag333The HST flight-system magnitudes were transformed to following Sirianni et al. (2005). Uncertainties are statistical. flux (Jy) Epoch (day)
F275W 888
F336W 888
F438W 888
F555W 888
F814W 888

Table 2Late-time HST Photometry of SN 2013ej.222MJD 57387 (day 888) adn MJD 57666 (day 1169).
Figure 12.— HST photometry from days 888 and 1169 (squares) shown with the earlier light-curve data from Dhungana et al. (2016), whose late-time decline fit of 0.011 mag day is illustrated by the green dashed line. The radioactive Co decay rate is represented by the solid gray line. The HST F814W (-band) progenitor level is illustrated by the solid orange horizontal line.
Figure 13.— RGB composite of the Spitzer/IRAC template images of M74 from 2004 (left), and the template-subtracted images of the region around SN 2013ej (tiled frames, right).

High-resolution imaging observations of SN 2013ej were performed with the HST and Wide-Field Camera 3 (WFC3) on 2015 Dec. 29 (day 888) under HST program GO-14116 (PI S. Van Dyk) and on 2016 Oct. 4 (day 1169) under HST Snapshot program GO-14668 (PI A. Filippenko). Exposures were obtained in the F275W (), F336W (), F438W (), F555W (), and F814W () filters on day 888, and only the latter two filters on day 1169. The point source is detected in all available bands from both of these late epochs, shown in the inset of Figure 11. The aim of the program was to detect a light echo around SN 2013ej.

Photometry of the source was extracted from the images using Dolphot. We tried two different approaches to estimate the background, including the use of an annulus region to measure the sky (FitSky=1) and, alternatively, measuring the sky within the point-spread function (PSF) aperture (FitSky=3; best to use when the field is very crowded). Our annulus-based background subtraction produced the most consistent results for all bands, although the results from each setting are within the respective uncertainties. The photometry is listed in Table 2.

Figure 12 shows the long-term light curve of SN 2013ej from the beginning of the explosion out to day 1169. Photometry of the cool supergiant progenitor of SN 2013ej, presented by Fraser et al. (2014), was contaminated by a neighboring blue source with (F435W) = 25.05 mag and (F555W) = 24.74 mag; the photocenter through those filters was offset from that of the SN. However, the photocenter of the F814W source in pre-explosion imaging, for which Fraser et al. (2014) measured m(F814W) =  mag, accurately matched that of the SN and has therefore been attributed to the progenitor star. We performed Dolphot runs on both of the available pre-explosion F814W epochs from 2003 Nov. 20 and 2005 Jun. 16, obtaining and , respectively (both values are averages from the results of Dolphot runs using FitSky=1 and FitSky=3 background estimation parameters), slightly brighter than the results from Fraser et al. (2014). Our photometry from days 888 and 1169 shows that the F814W flux of the SN remnant dropped below the faintest recorded state of the progenitor by  mag, confirming the RSG stars’s association with the SN.

Our day 888 and 1169 HST photometry shows that the SN remnant is substantially brighter than what is expected from both the earlier 0.011 mag day optical decline rate measured by Dhungana et al. (2016) and the expected rate of decay of Co to Fe. Moreover, the source has increased slightly in brightness between days 888 and 1169, by  mag and  mag in the and bands, respectively; the corresponding color change is a slight reddening by  mag. We attribute the slow rate of decline between the last ground-based coverage and the HST observations on day 888, and recent slight rise in flux between days 888 and 1169, to line emission from CSM interaction, which was revealed by our deep late-time spectroscopy presented in §3. The color change could be due to emission-line evolution or dust extinction (see §4.3). However, we note the possibility that the photometry in the bluer bands on day 888 might be slightly contaminated by the neighboring blue source identified by Fraser et al. (2014) in the pre-explosion data. The latest available data show that the SN remnant is still mag brighter than the blue neighboring source, so any contamination from the neighboring source is probably not substantial in our most recent measurements. Future epochs of HST photometry will help disentangle the SN emission from neighboring sources more accurately, once the emission from the SN and CSM interaction fades.

4.2. Keck/LRIS

We obtained a -band image of SN 2013ej on 2016 Feb. 9 with Keck/LRIS. Under 1 seeing conditions, we were unable to resolve the SN emission from neighboring sources. We thus rely on our higher-resolution HST images for our late-time optical photometry, but we show the SN position in the host galaxy with the -band LRIS image in Figure 11.

4.3. Warm Spitzer

SN 2013ej was observed between days 238 and 1006 during the Spitzer Space Telescope Warm Mission utilizing channels 1 (3.6 m) and 2 (4.5 m) of the Infrared Array Camera (IRAC; Fazio et al. 2004). We acquired fully coadded and calibrated data from the Spitzer Heritage Archive444, including data from Program IDs 11063, 11053, and 10046 (PIs Kasliwal, Fox, and Sanders, respectively). The first three epochs of photometry of SN 2013ej were reported by Tinyanont et al. (2016). We extracted photometry from all epochs independently, registering each SN image with earlier pre-SN images of the host galaxy, which were used as subtraction templates. The template-subtracted images are shown in Figure 13. We performed aperture photometry on the template-subtracted (PBCD / Level 2) images using a 6-pixel aperture radius and aperture corrections listed in Table 4.7 of the Spitzer IRAC Instrument Handbook555 The background and noise levels in the subtracted images were computed from 100 randomly placed apertures near the SN location.The photometry results are listed in Table 3.

The SN is detected at 4.5 m in all epochs up to day 1006 (see light curve in Figure 14). In the 3.6 m data, however, the SN is visible only until day 472, after which we conservatively list the flux as 3 upper limits in Table 3.

The mid-IR light curve of SN 2013ej is consistent with the general behavior of other core-collapse SNe (e.g., see Tinyanont et al. 2016). It most closely resembles the mid-IR light curves of SN 2004et, which was optically classified as a SN II-P (Li et al. 2005), and SN 2007od, which was also classified as a SN II-P (Andrews et al. 2010), but potentially classifiable as a SN II-L, based on the criterion of Anderson et al. (2014). Interestingly, both SNe exhibited evidence for late-time CSM interaction, similar to that of SN 2013ej (see Figure 10). Based on their optical to IR SEDs, the mid-IR counterparts of SN 2004et and SN 2007od have been shown to be consistent with emission from warm dust (Kotak et al. 2009; Andrews et al. 2010; Inserra et al. 2011). Late-time Spitzer observations of SN 2004et showed that the mid-IR flux increased dramatically after day 1000, interpreted as the result of dust formation within a dense post-shock shell that developed in response to the CSM interaction (Kotak et al. 2009). There is marginal evidence of a similar late-time increase developing in SN 2013ej; this is an interesting possibility, given the complementary evidence for CSM interaction at this late stage, and motivates continued monitoring.

MJD Epoch (day) 3.6 m (Jy)666Uncertainties are statistical. 4.5 m (Jy)
56735.0 237.6
56758.3 260.9
56936.6 439.2
56965.3 467.9
56970.1 472.7
57313.0 815.6
57320.5 823.1
57334.2 836.8
57474.9 977.5
57482.4 985.0
57503.6 1006.2
Table 3Spitzer Photometry of SN 2013ej.

The UV–IR spectral energy distribution (SED) of SN 2013ej is shown in Figure 15, combining the HST photometry from day 888 with the Spitzer photometry from day 1006, when the late 4.5 m flux has slightly rebrightened. The broadband optical photometry from day 888 is more or less consistent with the Keck/LRIS spectrum from day 807. The IR excess from the SN, like the mid-IR light curve, appears similar to that of SN 2004et and SN 2007od. The source of the excess is likely to be thermal, although it is also possible that the 4.5 m flux is also influenced by line emission.

Figure 14.— Long-term mid-IR light curve of SN 2013ej (black), and comparison objects SN 2004et (blue; Kotak et al. 2009), SN 2007od (green; Andrews et al. 2010), and SN 1993J (red; Tinyanont et al. 2016). Data for SN 2011dh (yellow; Helou et al. 2013; Ergon et al. 2015) and SN 1998S (magenta; Gerardy et al. 2002; Pozzo et al. 2004) are also included for comparison.
Figure 15.— Late-time optical through IR SED for the remnant of SN 2013ej, using the HST and Spitzer data listed in Tables 2 and 3. The day 807 Keck spectrum is included (gray curve). The SEDs of SN 2004et (Kotak et al. 2009) and SN 2007od (Andrews et al. 2010) are included for comparison, and their fluxes have been scaled to the distance of SN 2013ej.

4.3.1 Dust Emission Model

We analyze the Spitzer data using a simple dust model (Fox et al. 2011). We assume that the source of the mid-IR excess is warm dust, and fit the SED as a function of the dust mass () and temperature (). The flux is given by (Hildebrand 1983)


where  is the dust grain radius, is the dust absorption coefficient, and is the distance of the dust from the observer. We assume m and derive for graphite dust, following Fox et al. (2010; their Figure 4). For simplicity (and given the limited number of data points), we assume optically thin dust emitting at a single equilibrium temperature (e.g., Hildebrand 1983). We calculated the dust temperature and mass up to day 815; later epochs are excluded owing to low S/N. Table 4 summarizes the best-fitting parameters and limits. Our values are consistent with the results of Tinyanont et al. (2016), who performed similar calculations up to day 438. Including the longer range of dates presented here, there is only marginal indication of a decrease in dust temperature and mass, considering the uncertainties. Overall, the ranges of mass and temperature values are broadly consistent with the current sample of SNe II (Tinyanont et al. 2016).

Epoch (day) Temp. (K) Mass ()
Table 4Dust Emission Model Parameters.

5. Discussion

5.1. Origins of Continuum and Line Polarization

SN 2013ej exhibits the strongest and most persistent continuum and line polarization ever reported for a SN II-P or II-L during the recombination phase, which indicates substantial asphericity in the explosion and/or circumstellar environment. Modeling in §2.5 shows that the general features of strong continuum polarization and prominent modulations across line features can be mutually generated by electron scattering in an oblate ellipsoidal envelope, viewed nearly edge-on, or dust scattering in an asymmetric distribution of CSM. We suspect that both mechanisms are potentially involved in the strong and peculiar spectropolarimetric evolution of SN 2013ej. The relative contributions of each scattering mode with time are difficult to ascertain. However, the fact that the dust-scattering model does not fully reproduce the substantial depolarization of the emission lines during the recombination phase (for example, H and Ca ii on day 73), nor the narrowness and high radial velocity of the strongly polarized H feature, suggests that electron scattering is probably the dominant mode of polarization for most of the recombination phase.

By day 107, however, when recombination of the SN envelope is near complete, the increasing polarization level at the wavelengths of the emission peaks suggests that dust scattering becomes more influential as the optical depth to electron scattering drops. It is therefore interesting that on day 107, when the contribution of dust scattering might be at its maximum, the position angle remains within of the average value during the earlier recombination phase. The position angle also never deviates more than for the continuum at any time during the evolution of the SN. This implies that, if the spectropolarimetric evolution is indeed an evolving mix of dust scattering (CSM) and electron scattering (aspherical photosphere), then the two scattering geometries are probably similar. Below, we consider several physical scenarios that can potentially explain the spectropolarimetric data.

5.1.1 Polarization via Aspherical CSM Interaction

We suspect that the unusual polarization characteristics of SN 2013ej are the result of CSM interaction, which we already know was occurring during the recombination phase, based on the persistent X-ray emission (Chakraborti et al. 2016). We also know that the CSM interaction was geometrically asymmetric based on the multi-peaked morphology of broad H emission at late times (see Figure 9). The aspherical component of the scattering geometry could thus be governed by the aspherical CSM distribution. Specifically, as the SN ejecta expand into the CSM, the added heat from interaction could result in a radially extended electron-scattering photosphere along directions where the CSM is densest, increasing the Thomson optical depth and producing excess continuum luminosity. For CSM that is in an equatorially enhanced distribution, perhaps stemming from mass loss in a binary system (e.g., common-envelope evolution and/or non-conservative mass transfer), the photosphere of the outflow could develop a bulged equator as the ejecta photosphere overtakes and envelops the region of intense interaction, where the shock slows. This sort of buried interaction was proposed to explain the temporary diminishing of narrow-line emission from the peculiar interacting SN PTF11iqb (Smith et al. 2015), as well as its relatively high luminosity, relative to normal SNe II-P, and fast decline rate during the plateau phase. An equatorially enhanced photosphere will result in a strong polarization signal, if viewed at a suitable inclination angle; we suggest this to be the case for SN 2013ej. An illustration of this potential scenario is presented in Figure 16. Note, the detection of persistent X-rays by Chakraborti et al. (2016) implies CSM interaction outside the photosphere, with the RSG wind; the embedded equatorially enhanced interaction that potentially creates the ellipsoidal photosphere is a separate component. At late times, emission from the reverse shock-excited ejecta will be enhanced along the same equatorial direction, potentially creating the multi-peaked H morphology we observe at days.

Since it is the CSM distribution that fosters the aspherical component of the photosphere, any dust in the outer reaches of the same distribution of CSM could result in an added source of scattering and polarization in a geometry that is axisymmetric with the photospheric scattering. In the case of an ellipsoidal photosphere from an equatorial distribution of dusty CSM, the two sets of Stokes vectors from dust scattering and electron scattering could constructively interfere, resulting in persistently strong polarization with little drift in position angle as their relative contributions change with time. Meanwhile, the SN photosphere along the more weakly interacting poles will recombine at a faster rate, temporarily amplifying the effect of the equatorial component of the photosphere. This could produce the rise in polarization strength and a small change in position angle shown in Figure 3. As electron scattering diminishes, dust particles farther out in the same equatorially enhanced distribution of CSM could continue scattering and polarizing the SN light at a similar position angle.

Figure 16.— Illustration of our proposed scenario for producing an oblate ellipsoidal photosphere via CSM interaction (observer’s view). Geometric features are exaggerated for clarity, and not to scale.

The CSM-interaction scenario above could also explain the strong polarized H feature and its high degree of axisymmetry with respect to the continuum. Indeed, the fact that its high radial velocity matches that of the outermost ejecta (see Figure 9) implies that the region of aspherical interaction at late times and the asymmetric H absorption that is creating the polarized feature during the earlier recombination phase have a common origin in the outermost ejecta. As noted in §2.4.1, the rapid development of the narrow component of the polarized H feature coincided with the sudden appearance of the high-velocity “notch” in the absorption component of the flux line profile (see Figure 4, upper-left panel). Bose et al. (2015) previously interpreted the morphology of this absorption feature as the result of a separate high-velocity shell of H i excited by CSM interaction. Similar signs of weak CSM interaction have been noted in other SNe II-P, such as SN 2004dj (Chugai et al. 2007), SN 2009bw (Inserra et al. 2012), and SN 2012aw (Bose et al. 2013).

The process by which CSM interaction influences the morphology of the H absorption profile was investigated by Chugai et al. (2007); those authors showed that, upon encountering CSM, outer SN ejecta that would otherwise be optically thin are excited by X-rays from the forward and reverse shocks. This effect sustains high optical depth for H in the outer high-velocity ejecta for a longer period of time than in the absence of CSM interaction, resulting in an observable feature in the blue shoulder of the absorption profile (e.g., a “notch” or multiple minima). Since this added absorption occurs in the outermost ejecta, from the observer’s perspective the absorption will obscure only the weakly polarized forward-scattered light from the central region of the photosphere, thus amplifying the effect of the highly polarized light scattered from the limbs of the photosphere (near scattering angle), producing the strong and narrow polarization peak at a fast radial velocity and with a position angle that is axisymmetric with the continuum. Radial variations in CSM density could thus lead to a rapid evolution of the blue side of the absorption profile, which, in turn, will result in a rapidly evolving polarized line feature.

It is important to note that rapid changes in the polarized line spectrum can also occur in the dust-scattering (echo) scenario, since at any given time the total signal is a mix of direct SN light and highly polarized dust-scattered light from an earlier phase. Thus, if the SN suddenly develops enhanced line absorption, whether from enhanced CSM interaction or something else, the direct SN light will have a decreased contribution to the total flux. Since the more highly polarized echo spectrum has yet to reflect the additional absorption owing to larger light-travel time, there will be a temporary increase in the fraction of polarized light at the wavelength of the absorption enhancement, and a polarization feature will emerge and/or strengthen. We note that in the case of dust scattering, substantial deviations in position angle over line features are not expected.

In each of the scenarios postulated above, it is the CSM and interaction that is responsible for the continuum polarization and the spectral changes that lead to the rapid development of enhanced line polarization. We note that similarly strong continuum and line polarization was observed in the Type IIb SN 1993J, which also showed signs of asymmetric CSM interaction during its recombination phase and at late times (see Figure 10). This SN was also modeled as a nearly edge-on oblate ellipsoidal photosphere (Höflich et al. 1996), with potential contribution from circumstellar dust (Tran et al. 1997). It appears that SN 2013ej and SN 1993J could have exploded into similar asymmetric circumstellar environments of tenuous, dusty CSM.

CSM interaction also might also explain the unusually early appearance of polarized Si ii in SN 2013ej from day 11. Indeed, the presence of this feature, first noted by Leonard et al. (2013), indicates peculiar ionization characteristics for a young SN II, for which we suspect enhanced line optical depth from CSM interaction is potentially responsible. The lack of substantial changes in across this feature (see Figure 2), similar to the polarized H feature, is also consistent with a high degree of axisymmetry.

5.1.2 Polarization via Aspherical Explosion Geometry

We must also consider the possibility that the electron-scattering mode of polarization is the result of an intrinsically aspherical explosion geometry, irrespective of CSM influence. In this case, the early onset of polarization would imply that the electron-scattering photosphere was aspherical while in the hydrogen envelope. In previous spectropolarimetric studies of normal SNe II-P, strong continuum polarization has only been observed to develop at relatively late times, near the onset of the nebular phase, after the photosphere has moved into the helium layer. Asphericity that appears only in deep layers has been attributed to an asymmetric distribution of radioactive material from the core (Wang et al. 2003; Leonard et al. 2006; Chornock et al. 2010). The resulting asymmetric distribution of heating produces an aspherical photosphere. Since radioactive asphericities are core-born, they are not expected to be observable at earlier phases when the photosphere is far out in the thick hydrogen envelope, unless such structures can penetrate the envelope rapidly — for example, in the form of energetic jets of Ni. Simulations by Couch et al. (2009) showed that non-relativistic radioactive jets that have a high ratio of kinetic to thermal energy move ballistically, and can quickly punch through the thick RSG envelope and produce an aspherical SN photosphere at the earliest phases. We thus speculate that the polarization of SN 2013ej could be the signature of such a jet-driven bipolar explosion. It is noteworthy in this regard that a bipolar distribution of Ni was invoked to explain the nebular-phase line profile of H in SN 2013ej (Bose et al. 2015; Yuan et al. 2016). The spectropolarimetry of a bipolar explosion, which could be approximated by a prolate ellipsoid, might appear similar to that of the oblate ellipsoidal configuration that we modeled in §2.5.2. It is also plausible that such a bipolar explosion geometry results in the multi-peaked morphology of the broad H emission that we observed at days. The interpretation of a ballistic jet is also interesting in that the models from Couch et al. (2009) have shown that some of the the radioactive material overtakes the hydrogen envelope within the first week. We thus question whether the early appearance of Si ii line polarization could be an additional signature of high-velocity Ni from a jet-driven explosion.

5.2. The Source of Warm Dust: Ejecta or CSM?

The warm dust responsible for the mid-IR emission could involve multiple and physically distinct components. The light curve of SN 2013ej and the inferred dust temperature and mass are broadly consistent with the properties of most SNe II (Tinyanont et al. 2016), including several objects that exhibit evidence for interaction with tenuous CSM. Pre-existing dust in the circumstellar environment can be heated by the SN radiation and/or swept up and heated by the interaction shock, generating a source of mid-IR emission. Interestingly, our best-matching dust-scattering model for the spectropolarimetric data on SN 2013ej involved a dust component located at a distance of  km from the SN. Assuming a speed of 12,000 km s for the fastest ejecta, this material would encounter the same light-scattering CSM dust at days, which is near the time we see the mid-IR emission possibly beginning to rise (see Figure 14). In the relatable cases of SN 2004et and SN 2007od, both which exhibited late-time interaction-powered spectra similar to that of SN 2013ej (see §3), their thermal excesses were interpreted as a result of dust formation within a dense post-shock shell that developed in response to the CSM interaction (Kotak et al. 2009; Andrews et al. 2010; Inserra et al. 2011). This interpretation is therefore appealing for SN 2013ej as well, not only because of its spectroscopic and SED similarities with SN 2004et and SN 2007od, but because of the possible signs of scattering by pre-existing circumstellar dust indicated by the spectropolarimetry in §2.5.1. Interestingly, Leonard et al. (2009) noted that SN 2004et also exhibited substantial continuum polarization weeks into the photospheric phase, while preliminary spectropolarimetry of SN 2007od also showed persistent continuum polarization at the level of %–1.0% throughout its photospheric phase (R. Chornock, private communication). All of these similarities seem to promote a physical connection between polarization, CSM interaction, and circumstellar dust.

It is not clear, however, if the dust producing the mid-IR emission from SN 2013ej is the same dust responsible for the red-side flux suppression of the broad H profile on days 807 and 897. Indeed, this suggests that the obscuring dust is interior to the reverse shock, preferentially obscuring the receding ejecta on the far side of the outflow. If the obscuring dust were in a post-shock shell on the outside of the reverse shock, then the entire H profile should be evenly obscured. A clumpy distribution of dust in the post-shock medium, as opposed to a continuous shell, might provide a solution, as some blueshifted emission from the near side of the reverse shock could potentially avoid obscuration, thus resulting in the asymmetric line profiles we observe. Such a clumpy distribution of dusty material was suggested in the relatable case of SN 2007od (Andrews et al. 2010; Inserra et al. 2011).

Another possible source of thermal dust emission is the metal-rich SN ejecta, interior to the reverse shock. Such a configuration naturally results in a suppressed red side for the H line profile from shock-excited ejecta. This hypothesis might also be more consistent with the dust masses inferred from our calculations in §4.3.1 (Eq. 1), since such a pre-existing CSM dust mass would require a total CSM much more massive than expected around a normal RSG, and more on par with strongly interacting SNe IIn (e.g., see Fox et al. 2011) — unless the CSM had a different origin than an RSG wind (e.g., equatorial mass loss via binary interactions). Of course, it is also possible that the dust producing the mid-IR emission is not the same dust responsible for the uneven obscuration of the H emission. Indeed, as stated previously, multiple components of dust within the ejecta, CSM, and post-shock medium could contribute to the mid-IR evolution.

6. Summary and Concluding Remarks

We have presented spectropolarimetry, deep late-time spectroscopy, and space-based UV–IR photometry of the weakly interacting SN 2013ej. From day 11 through 107, the SN exhibits the strongest and most persistent continuum polarization ever reported for a SN II-P or II-L, and also unusually prominent and narrowly peaked line polarization for this type of SN, particularly after months. Preliminary modeling in §2.5.1 has shown that the continuum polarization and broad modulations across the line features can be mutually generated by an ellipsoidal photosphere, and by optical scattering off dust particles in an asymmetric distribution of CSM. The spectropolarimetric evolution of SN 2013ej could be governed by a complex time-dependent combination of Stokes vectors from both of these scattering processes. However, electron scattering in an ellipsoidal photosphere is likely to be the dominant source of polarization, particularly on day 73, when nearly full depolarization of the H and Ca ii emission-line peaks is observed. The dust-scattering mode appears to become more influential on our last epoch at the end of the recombination phase (day 107), based on the increased amount of residual polarization for the H emission peak, and on the reasonable match of our model. The continued increase in fractional polarization indicated by the late-time imaging polarimetry from Kumar et al. (2016) could indicate that scattering by circumstellar dust continues beyond day 130, or that the asphericity in the inner electron-scattering environment became even more pronounced.

Deep optical spectroscopy past 800 days reveals persistent and geometrically asymmetric CSM interaction into late phases, in the form of broad multi-peaked emission-line profiles and a blue pseudocontinuum. The broad spectral features are indicative of SN ejecta energized by the passage of the reverse shock, and exhibit similarities with other well-studied core-collapse SNe that have exhibited evidence for persistent interaction, substantial polarization, and late-time thermal excess from dust, including SN 1993J (Matheson et al. 2000; Höflich et al. 1996; Tran et al. 1997), SN 2004et (Kotak et al. 2009; polarization noted by Leonard et al. 2009), and SN 2007od (Andrews et al. 2010; polarization information from R. Chornock, private communication). The late time UV-optical photometry of SN 2013ej from HST, out to day 1169, which is much brighter than what is expected from radioactive decay, is also indicative of ongoing CSM interaction. The CSM responsible for this late-time interaction was probably created at least  yr before core collapse. The hint of triple-peaked line morphology is interesting in the context of the spectropolarimetry, as it appears to be consistent with an equatorially enhanced distribution of CSM, perhaps stemming from mass loss in a binary system.

We currently lean toward a scenario in which interaction with an equatorially enhanced distribution of CSM, viewed nearly edge-on, is responsible for the photospheric asphericity, the resulting polarization, and (perhaps) for providing some additional luminosity early in the recombination phase that may have influenced the relatively rapid decline rate. Dust in the same distribution of CSM could provide an additional source of scattering that is nearly axisymmetric with the interaction-influenced photosphere, constructively adding to the net polarization. The late-time mid-IR emission detected with Spitzer out to days could involve the same circumstellar dust that potentially contributes to the optical scattering at earlier phases, radiatively heated by the propagating interaction shock. This scenario is reminiscent of interpretations previously suggested to explain the spectropolarimetry of the Type IIb SN 1993J (Höflich et al. 1996), and of more strongly interacting SNe IIn, such as SN 2009ip (Mauerhan et al. 2014), SN 1998S (Leonard et al. 2000; Wang et al 2001), and SN 1997eg (Hoffman et al. 2008). It is noteworthy that, if viewed at a suitable inclination angle, weakly interacting SNe can potentially produce continuum and line polarization levels comparable to those of SNe IIn, even if the CSM densities and shock strengths are not high enough to generate the strong narrow emission lines that characterize the flux spectra of SNe IIn. Future studies involving a larger sample of objects for which there is high-quality spectropolarimetry, deep late-time spectroscopy, and mid-IR photometry will further elucidate the physical connections between CSM interaction, asphericity/polarization, light-curve evolution, and thermal dust emission.

We also considered the alternative scenario that the polarization is produced by an intrinsically aspherical explosion, perhaps the result of a ballistic jet that drives a bipolar outflow. If this interpretation is correct, it would be remarkable that such a ballistic jet could result from the collapse of a rather modest 9.5–15.5 M RSG, which is the mass estimated for the directly identified progenitor of SN 2013ej (Fraser et al. 2014).

We note that, although the physical interpretations we have considered are plausible and guided by theoretical expectations, none offer unique solution to the data on SN 2013ej. Models involving more detailed treatment of complex scattering geometries and their combined spectropolarimetric effects are required, and we look forward to their construction in the future.

Finally, the late-time UV–optical photometry from HST shows that the interaction-powered light curve faded (temporarily, at least) well below the level of the RSG progenitor, confirming its association with the explosion (Fraser et al. 2013; Van Dyk et al. 2013).


J.C.M. would like to thank Nathan Smith for informative scientific discussion and commentary on this manuscript. Support for HST programs GO-14116 and GO-14668 was provided by NASA through grants from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. We are grateful to the staffs at Lick and Keck Observatories for their excellent assistance. This work is based in part on observations from the Low Resolution Imaging Spectrometer at the Keck-1 telescope. The W. M. Keck Observatory is operated as a scientific partnership among the California Institute of Technology, the University of California, and NASA; it was made possible by the generous financial support of the W. M. Keck Foundation. We extend special thanks to those of Hawaiian ancestry on whose sacred mountain we are privileged to be guests. A.V.F.’s group at U.C. Berkeley is supported by NSF grant AST-1211916, Gary & Cynthia Bengier, the Richard & Rhoda Goldman Fund, the Christopher R. Redlich Fund, and the TABASGO Foundation. Research at Lick Observatory is partially supported by a generous gift from Google. This work is based in part on observations and in part on archival data obtained with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under a contract with NASA; support was provided by NASA through an award issued by JPL/Caltech.


  • Anderson et al. (2014) Anderson, J. P., González-Gaitán, S., Hamuy, M., et al. 2014, ApJ, 786, 67
  • Andrews et al. (2010) Andrews, J. E., Gallagher, J. S., Clayton, G. C., et al. 2010, ApJ, 715, 541
  • Bose et al. (2013) Bose, S., Kumar, B., Sutaria, F., et al. 2013, MNRAS, 433, 1871
  • Bose et al. (2015) Bose, S., Sutaria, F., Kumar, B., et al. 2015, ApJ, 806, 160
  • Branch et al. (2007) Branch, D., Parrent, J., Troxel, M. A., et al. “The Multicolored Landscape of Compact Objects and Their Explosive Origins”, 200, APIC Conf. Series, 924, 342, eds. di Salvo, T. and Israel, G. L. and Piersant, L. and Burderi, L. and Matt, G. and Tornambe, A. and Menna, M. T.
  • Chakraborti et al. (2016) Chakraborti, S., Ray, A., Smith, R., et al. 2016, ApJ, 817, 22
  • Chevalier & Fransson (1994) Chevalier, R. A., & Fransson, C. 1994, ApJ, 420, 268
  • Chornock et al. (2010) Chornock, R., Filippenko, A. V., Li, W., & Silverman, J. M. 2010, ApJ, 713, 1363
  • Chugai (1992) Chugai, N. N. 1992, Soviet Astronomy Letters, 18, 168
  • Chugai (2006) Chugai, N. N. 2006, Astronomy Letters, 32, 739
  • Chugai et al. (2007) Chugai, N. N., Chevalier, R. A., & Utrobin, V. P. 2007, ApJ, 662, 1136
  • Couch et al. (2009) Couch, S. M., Wheeler, J. C., & Milosavljević, M. 2009, ApJ, 696, 953
  • Davis (1955) Davis, L., Jr. 1955, Vistas in Astronomy, 1, 336
  • Dessart & Hillier (2011) Dessart, L., & Hillier, D. J. 2011, MNRAS, 415, 3497
  • Dhungana et al. (2016) Dhungana, G., Kehoe, R., Vinko, J., et al. 2016, ApJ, 822, 6
  • Ergon et al. (2015) Ergon, M., Jerkstrand, A., Sollerman, J., et al. 2015, A&A, 580, A142
  • Faber et al. (2003) Faber, S. M., Phillips, A. C., Kibrick, R. I., et al. 2003, SPIE, 4841, 1657
  • Faran et al. (2014) Faran, T., Poznanski, D., Filippenko, A. V., et al. 2014, MNRAS, 442, 844
  • Fazio et al. (2004) Fazio, G. G., Hora, J. L., Allen, L. E., et al. 2004, ApJS, 154, 10
  • Filippenko (1982) Filippenko, A. V. 1982, PASP, 94, 715
  • Foley et al. (2003) Foley, R. J., Papenkova, M. S., Swift, B. J., et al. 2003, PASP, 115, 1220
  • Fox et al. (2011) Fox, O. D., Chevalier, R. A., Skrutskie, M. F., et al. 2011, ApJ, 741, 7
  • Fraser et al. (2014) Fraser, M., Maund, J. R., Smartt, S. J., et al. 2014, MNRAS, 439, L56
  • Gall et al. (2015) Gall, E. E. E., Polshaw, J., Kotak, R., et al. 2015, A&A, 582, A3
  • Gerardy et al. (2002) Gerardy, C. L., Fesen, R. A., Nomoto, K., et al. 2002, ApJ, 575, 1007
  • Helou et al. (2013) Helou, G., Kasliwal, M. M., Ofek, E. O., et al. 2013, ApJL, 778, L19
  • Hildebrand (1983) Hildebrand, R. H. 1983, QJRAS, 24, 267
  • Hiltner (1949) Hiltner, W. A. 1949, ApJ, 109, 471
  • Hoffman et al. (2008) Hoffman, J. L., Leonard, D. C., Chornock, R., et al. 2008, ApJ, 688, 1186-1209
  • Höflich et al. (1996) Höflich, P., Wheeler, J. C., Hines, D. C., & Trammell, S. R. 1996, ApJ, 459, 307
  • Huang et al. (2015) Huang, F., Wang, X., Zhang, J., et al. 2015, ApJ, 807, 59
  • Indebetouw et al. (2005) Indebetouw, R., Mathis, J. S., Babler, B. L., et al. 2005, ApJ, 619, 931
  • Inserra et al. (2012) Inserra, C., Turatto, M., Pastorello, A., et al. 2012, MNRAS, 422, 1122
  • Inserra et al. (2011) Inserra, C., Turatto, M., Pastorello, A., et al. 2011, MNRAS, 417, 261
  • Kasen et al. (2003) Kasen, D., Nugent, P., Wang, L., et al. 2003, ApJ, 593, 788
  • Kawabata et al. (2002) Kawabata, K. S., Jeffery, D. J., Iye, M., et al. 2002, ApJL, 580, L39
  • Kervella et al. (2016) Kervella, P., Lagadec, E., Montargès, M., et al. 2016, A&A, 585, A28
  • Kervella et al. (2011) Kervella, P., Perrin, G., Chiavassa, A., et al. 2011, A&A, 531, A117
  • Kotak et al. (2009) Kotak, R., Meikle, W. P. S., Farrah, D., et al. 2009, ApJ, 704, 306
  • Kumar et al. (2016) Kumar, B., Pandey, S. B., Eswaraiah, C., & Kawabata, K. S. 2016, MNRAS, 456, 3157
  • Leonard et al. (2013) Leonard, D. C., Pignata, G., Dessart, L., et al. 2013, ATel, 5275.
  • Leonard et al. (2009) Leonard, D. C., Filippenko, A. V., Ganeshalingam, M., Li, W., et al. 2009, BAAS 41, 466
  • Leonard et al. (2006) Leonard, D. C., Filippenko, A. V., Ganeshalingam, M., et al. 2006, Nature, 440, 505
  • Leonard et al. (2002) Leonard, D. C., Filippenko, A. V., Chornock, R., & Foley, R. J. 2002, PASP, 114, 1333
  • Leonard et al. (2001) Leonard, D. C., Filippenko, A. V., Ardila, D. R., & Brotherton, M. S. 2001, ApJ, 553, 861
  • Leonard et al. (2000) Leonard, D. C., Filippenko, A. V., Barth, A. J., & Matheson, T. 2000, ApJ, 536, 239
  • Li et al. (2005) Li, W., Van Dyk, S. D., Filippenko, A. V., & Cuillandre, J.-C. 2005, PASP, 117, 121
  • Matheson et al. (2000) Matheson, T., Filippenko, A. V., Ho, L. C., Barth, A. J., & Leonard, D. C. 2000, AJ, 120, 1499
  • Mauerhan & Smith (2012) Mauerhan, J., & Smith, N. 2012, MNRAS, 424, 2659
  • Mauerhan et al. (2014) Mauerhan, J., Williams, G. G., Smith, N., et al. 2014, MNRAS, 442, 1166
  • Maund et al. (2007) Maund, J. R., Wheeler, J. C., Patat, F., et al. 2007, ApJ, 671, 1944
  • Milisavljevic et al. (2012) Milisavljevic, D., Fesen, R. A., Chevalier, R. A., et al. 2012, ApJ, 751, 25
  • Miller et al. (1988) Miller, J. S., Robinson, L. B., & Goodrich, R. W. 1988, in Instrumentation for Ground-Based Optical Astronomy, ed. L. B. Robinson (New York: Springer), 157
  • Miller & Stone (1993) Miller, J. S., & Stone, R. P. S. 1993, Lick Obs. Tech. Rep. 66 (Univ. California, Santa Cruz)
  • Ohnaka (2014) Ohnaka, K. 2014, A&A, 568, A17
  • Oke et al. (1995) Oke, J. B., Cohen, J. G., Carr, M., et al. 1995, PASP, 107, 375
  • Pozzo et al. (2004) Pozzo, M., Meikle, W. P. S., Fassia, A., et al. 2004, MNRAS, 352, 457
  • Sahu et al. (2006) Sahu, D. K., Anupama, G. C., Srividya, S., & Muneer, S. 2006, MNRAS, 372, 1315
  • Sanders et al. (2015) Sanders, N. E., Soderberg, A. M., Gezari, S., et al. 2015, ApJ, 799, 208
  • Schlafly & Finkbeiner (2011) Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103
  • Schmidt et al. (1992) Schmidt, G. D., Elston, R., & Lupie, O. L. 1992a, AJ, 104, 1563
  • Schmidt et al. (1992) Schmidt, G. D., Stockman, H. S., & Smith, P. S. 1992b, ApJL, 398, L57
  • Serkowski et al. (1975) Serkowski, K., Mathewson, D. S., & Ford, V. L. 1975, ApJ, 196, 261
  • Sirianni et al. (2005) Sirianni, M., Jee, M. J., Benítez, N., et al. 2005, PASP, 117, 1049
  • Smartt (2009) Smartt, S. J. 2009, ARA&A, 47, 63
  • Smith et al. (2015) Smith, N., Mauerhan, J. C., Cenko, S. B., et al. 2015, MNRAS§, 449, 1876
  • Tanaka et al. (2012) Tanaka, M., Kawabata, K. S., Hattori, T., et al. 2012, ApJ, 754, 63
  • Tanaka et al. (2008) Tanaka, M., Kawabata, K. S., Maeda, K., Hattori, T., & Nomoto, K. 2008, ApJ, 689, 1191
  • Tinyanont et al. (2016) Tinyanont, S., Kasliwal, M. M., Fox, O. D., et al. 2016, arXiv:1601.03440
  • Tran et al. (1997) Tran, H. D., Filippenko, A. V., Schmidt, G. D., et al. 1997, PASP, 109, 489
  • Valenti et al. (2016) Valenti, S., Howell, D. A., Stritzinger, M. D., et al. 2016, MNRAS, 459, 3939
  • Valenti et al. (2014) Valenti, S., Sand, D., Pastorello, A., et al. 2014, MNRAS, 438, L101
  • Van Dyk et al. (2013) Van Dyk, S. D., Petigura, E. A., Cenko, S. B., et al.  2013, ATel, 5229
  • Wang et al. (2003) Wang, L., Baade, D., Höflich, P., & Wheeler, J. C. 2003, ApJ, 592, 457
  • Wang & Wheeler (1996) Wang, L., & Wheeler, J. C. 1996, ApJL, 462, L27
  • Wang & Wheeler (2008) Wang, L., & Wheeler, J. C. 2008, ARAA, 46, 433
  • Whittet et al. (1992) Whittet, D. C. B., Martin, P. G., Hough, J. H., et al. 1992, ApJ, 386, 562
  • Yaron & Gal-Yam (2012) Yaron, O., & Gal-Yam, A. 2012, PASP, 124, 668
  • Yuan et al. (2016) Yuan, F., Jerkstrand, A., Valenti, S., et al. 2016, MNRAS, 461, 2003
  • Zuckerman (1980) Zuckerman, B. 1980, ARA&A, 18, 263
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description