Draft version July 2, 2021 Typeset using LATEX twocolumn style in AASTeX63 Magnetar Models of Superluminous Supernovae from the Dark Energy Survey: Exploring Redshift Evolution Brian Hsu,1 Griffin Hosseinzadeh,1 and Edo Berger1 1Center for Astrophysics j Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138-1516, USA ABSTRACT Superluminous supernovae (SLSNe) are luminous transients that can be detected to high redshifts with upcoming optical time-domain surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). An interesting open question is whether the properties of SLSNe evolve through cosmic time. To address this question, in this paper we model the multi-color light curves of all 21 Type I SLSNe from the Dark Energy Survey (DES) with a magnetar spin-down engine, implemented in the Modular Open Source Fitter for Transients (MOSFiT). With redshifts up to z ≈ 2, this sample includes some of the highest-redshift SLSNe. We find that the DES SLSNe span a similar range of ejecta and magnetar engine parameters to previous samples of mostly lower-redshift SLSNe (spin 14 period P ≈ 0:79 − 13:61 ms, magnetic field B ≈ (0:03 − 7:33) × 10 G, ejecta mass Mej ≈ 1:54 − 30:32 4 −1 M , and ejecta velocity vej ≈ (0:55 − 1:45) × 10 km s ). The DES SLSN sample by itself exhibits the previously found negative correlation between Mej and P , with a pronounced absence of SLSNe with low ejecta mass and rapid spin. Combining our results for the DES SLSNe with 60 previous SLSNe modeled in the same way, we find no evidence for redshift evolution in any of the key physical parameters. Keywords: Supernovae (1668) 1. INTRODUCTION CCSNe. Interaction with a hydrogen-poor circumstellar Type I superluminous supernovae (hereafter, SLSNe) medium could in principle provide sufficient radiative are a rare subclass of core-collapse supernovae (CCSNe) energy, but there is no clear spectroscopic evidence for that have been discovered by wide-field time-domain op- such interaction (Jerkstrand et al. 2017; Liu et al. 2017). tical surveys over the past decade (Chomiuk et al. 2011; Instead, both the diverse light curves and the spectral Quimby et al. 2011). They were originally defined to evolution of SLSNe can be explained well with a cen- have a peak absolute magnitude of M < −21 (Gal-Yam tral engine, namely a rapidly spinning (∼ few ms) and 14 2012), but are now defined spectroscopically by an ab- highly magnetized (∼ 10 G) neutron star (a \mag- sence of hydrogen features, blue continua, and, usually, netar"; Kasen & Bildsten 2010; Woosley 2010; Dessart unique early-time \W"-shaped O II absorption lines at et al. 2012; Metzger et al. 2015; Nicholl et al. 2017c). ∼ 3600 − 4600 A˚ (e.g. Lunnan et al. 2013; Mazzali et al. This model accounts for the broad range of peak lumi- 2016; Quimby et al. 2018). nosities and timescales (e.g., Nicholl et al. 2017c; Blan- SLSNe radiate 10 − 100 times more energy in the chard et al. 2020), for the early UV/optical spectra (e.g., arXiv:2104.09639v2 [astro-ph.HE] 1 Jul 2021 UV/optical compared to normal CCSNe, and generally Nicholl et al. 2017b), for the nebular phase spectra (e.g., exhibit longer durations (e.g. Nicholl et al. 2015; Inserra Nicholl et al. 2016b, 2019; Jerkstrand et al. 2017), and, et al. 2017; Lunnan et al. 2018; De Cia et al. 2018). As most recently, for the power law decline rate observed 3 a result of their light curve behavior and spectral prop- in SN 2015bn at ∼ 10 d (Nicholl et al. 2018). Addi- erties, it has become clear that SLSNe are not powered tional support for a magnetar engine comes from the low by radioactive decay of 56Ni as in normal hydrogen-poor metallicity host galaxies of SLSNe, which most closely resemble the hosts of long gamma-ray bursts, another rare population of CCSNe that are powered by a central Corresponding author: Brian Hsu engine (e.g., Lunnan et al. 2014; Perley et al. 2016, but [email protected] cf. Schulze et al. 2018). 2 Hsu, Hosseinzadeh, & Berger 20 Systematic analyses of some SLSN samples in the con- DES 15 Blanchard et al. (2020) text of a magnetar engine model have been carried out 10 in recent years by Nicholl et al.(2017c), Villar et al. 5 (2018), and Blanchard et al.(2020) using the semi- 0 e 22.5 analytical Modular Open-Source Fitter for Transients d u t (MOSFiT; Guillochon et al. 2018). Recently, Angus et al. i 22.0 n g (2019) published a sample of 21 SLSNe at z ≈ 0:2 − 2 a M 21.5 discovered by the Dark Energy Survey (DES), presented e t u their general observed properties, and derived physical l 21.0 o s DES13S2cmm DES15X1noe properties from magnetar model fits to Gaussian pro- b DES14C1fi DES15X3hm A 20.5 DES14C1rhg DES16C2aix cess interpolated pseudo-bolometric light curves. Here, d DES14E2slp DES16C2nm n DES14S2qri DES16C3cv a DES14X2byo DES16C3dmp we use MOSFiT to model the multi-color light curves of B 20.0 - DES14X3taz DES16C3ggu g these DES SLSNe with a magnetar engine for the first DES15C3hav DES17C3gyp k DES15E2mlf DES17X1amf a 19.5 DES15S1nog DES17X1blv time to systematically compare them to the previously- e DES15S2nr P published SLSN samples. Taking advantage of the broad 0.0 0.5 1.0 1.5 2.0 0 2 4 6 8 redshift range of the DES sample, we also use the com- Redshift bined sample of about 80 SLSNe to explore any redshift Figure 1. Rest-frame peak g-band absolute magnitude ver- evolution in the engine and SN ejecta parameters. sus redshift for the DES SLSNe. The top panel shows the The paper is structured as follows. The DES SLSN redshift distribution of the DES SLSNe (green) compared to sample is summarized in §2. The model fitting proce- the SLSN compilation recently analyzed by Blanchard et al. dure is described in §3. The resulting fits and model (2020). Square symbols indicate observed pre-peak bumps in parameters are presented in §4, while in §5 we place the light curve via visual inspection. Solid and open mark- these results in the context of previous SLSN samples ers indicate the Gold and Silver spectroscopic classification, respectively, of Angus et al.(2019). analyzed in the same manner. Throughout the paper, we assume a flat ΛCDM cosmology with Ωm = 0:308 −1 et al.(2019) for Milky Way extinction using Schlafly & and H0 = 67:8 km s Mpc, based on the Planck 2015 results (Planck Collaboration 2016). Finkbeiner(2011). In Figure1 we plot the observed peak absolute magni- tudes versus redshift for the DES sample, as well as the 2. DATA SET redshift distribution in comparison to the SLSN compli- We model the DES SLSN sample presented in Angus ation from Blanchard et al.(2020). Because we do not et al.(2019). DES (Dark Energy Survey Collaboration have a complete set of spectra for the DES SLSNe, we 2005, 2016) was an optical imaging survey covering 5000 correct our observed peak magnitudes to a single rest- deg2 of the southern sky, using the 570 megapixel Dark frame filter by calculating synthetic photometry in the Energy Camera (DECam; Flaugher et al. 2015) on the rest-frame g-band and the closest observed filter using 4-m Blanco Telescope at Cerro Tololo Inter-American the modified blackbody SED returned by MOSFiT (see Observatory in Chile. The DES-SN component of the §3) at the epoch corresponding to the observed peak survey (Bernstein et al. 2012), accounting for ≈ 10% magnitude. See Table1 for a complete list of redshift, of the total survey time, covered ten 3-deg2 individual observer-frame filter used, K correction, and rest-frame fields in g, r, i, z with a roughly weekly cadence to typ- peak g-band absolute magnitude for each DES SLSN. ical limiting magnitudes of ≈ 23:5 (8 fields) and ≈ 24:5 The peak absolute magnitudes span a wide range of (2 fields). See Kessler et al.(2015) for a more detailed ≈ −19:5 to ≈ −22:5, with an apparent gap1 at about description of the observing strategy, and Morganson −20:5 to −21:5. As expected, the lower luminosity et al.(2018) for the DES image processing pipeline. events are mostly confined to lower redshift (z . 0:7) Angus et al.(2019) classified transient events as while the more luminous events extend over z ≈ 0:6 − 2. SLSNe using SUPERFIT (Howell et al. 2005) with the Angus et al.(2019) suspected that the intrinsic broad SLSN spectral template library of Quimby et al.(2018). Their final sample consists of 21 SLSNe with a redshift 1 range of z = 0:220 − 1:998, separated into two classi- There is one event, DES16C2aix, that has a peak rest-frame g- band absolute magnitude in this region. However, we note that fication \standards" of Gold (16 events) and Silver (5 the peak magnitude for DES16C2aix has large uncertainties and events) based on their confidence in the spectroscopic the next brightest point in rest-frame g-band resides on the lower classification, listed in Table1. In the following analy- bound of the plotted peak magnitude. sis we correct the photometric measurements in Angus DES SLSNe 3 Table 1. Peak Magnitudes of DES SLSNe in the Rest-frame g-band λeff of Closest Observed Peak Mag. Rest-frame SLSN Name Redshift Rest-frame Observed Peak MJD K correction g-band Peak Standardy g-band (A)˚ Band Mag. (J2000) Abs. Mag. DES13S2cmm 0.663 7934 i 22.42 56567.18 −0:56 −20:08 Gold DES14C1fi 1.302 10982 z 22.09 56956.35 −1:13 −21:65 Silver DES14C1rhg 0.481 7066 r 22.30 57012.08 −0:55 −19:37 Gold DES14E2slp 0.57 7490 i 21.60 57042.09 −0:61 −20:46 Silver DES14S2qri 1.50 11927 z 22.31 57018.05 −1:18 −21:76 Gold DES14X2byo 0.868 8912 z 21.37 56956.21 −0:84 −21:57 Gold DES14X3taz 0.608 7671 i 20.60 57078.02 −0:58 −21:66 Gold DES15C3hav 0.392 6641 r 21.84 57339.14 −0:38 −19:48 Gold DES15E2mlf