Search

Massive stars end their life as core-collapse supernovae, amongst which some extremes are Type Ic broad-lined supernovae associated with long-duration gamma-ray bursts (LGRBs) having powerful relativistic jets. Their less-extreme brethren make unsuccessful jets that are choked inside the stars, appearing as X-ray flashes or low-luminosity GRBs. On the other hand, there exists a population of extragalactic fast X-ray transients (EFXTs) with timescales ranging from seconds to thousands of seconds, whose origins remain obscure. Known sources that contribute to the observed EFXT population include the softer analogs of LGRBs, shock breakouts of supernovae, or unsuccessful jets. Here, we report the discovery of the bright X-ray transient EP240414a detected by the Einstein Probe (EP), which is associated with the Type Ic supernova SN 2024gsa at a redshift of 0.401. The X-ray emission evolution is characterised by a very soft energy spectrum peaking at < 1.3 keV, which makes it distinct from known LGRBs, X-ray flashes, or low-luminosity GRBs. Follow-up observations at optical and radio bands revealed the existence of a weak relativistic jet that interacts with an extended shell surrounding the progenitor star. Located on the outskirts of a massive galaxy, this event reveals a new population of explosions of Wolf-Rayet stars characterised by a less powerful engine that drives a successful but weak jet, possibly owing to a progenitor star with a smaller core angular momentum than in traditional LGRB progenitors., Comment: 43 pages, 9 figures, 4 tables, submitted. Comments are welcome

Hydrogen-rich superluminous supernovae (SLSNe II) are rare. The exact mechanism producing their extreme light curve peaks is not understood. Analysis of single events and small samples suggest that CSM interaction is the main responsible for their features. However, other mechanisms can not be discarded. Large sample analysis can provide clarification. We aim to characterize the light curves of a sample of 107 SLSNe II to provide valuable information that can be used to validate theoretical models. We analyze the gri light curves of SLSNe II obtained through ZTF. We study peak absolute magnitudes and characteristic timescales. When possible we compute g-r colors, pseudo-bolometric light curves, and estimate lower limits for their total radiated energy. We also study the luminosity distribution of our sample and estimate the percentage of them that would be observable by the LSST. Finally, we compare our sample to other H-rich SNe and to H-poor SLSNe I. SLSNe II are heterogeneous. Their median peak absolute magnitude is -20.3 mag in optical bands. Their rise can take from two weeks to over three months, and their decline from twenty days to over a year. We found no significant correlations between peak magnitude and timescales. SLSNe II tend to show fainter peaks, longer declines and redder colors than SLSNe I. We present the largest sample of SLSNe II light curves to date, comprising of 107 events. Their diversity could be explained by considering different CSM morphologies. Although, theoretical analysis is needed to explore alternative scenarios. Other luminous transients, such as Active Galactic Nuclei, Tidal Disruption Events or SNe Ia-CSM, can easily become contaminants. Thus, good multi-wavelength light curve coverage becomes paramount. LSST could miss 30 percent of the ZTF events in the its footprint in gri bands. Redder bands become important to construct complete samples., Comment: 35 pages. 27 figures. 9 tables

More than 3000 spectroscopically confirmed Type Ia supernovae (SNe Ia) are presented in the Zwicky Transient Facility SN Ia Data Release 2 (ZTF DR2). In this paper, we detail the spectral properties of 482 SNe Ia near maximum light, up to a redshift limit of $z$ $\leq$ 0.06. We measure the velocities and pseudo-equivalent widths (pEW) of key spectral features (Si II $\lambda$5972 and Si II $\lambda$6355) and investigate the relation between the properties of the spectral features and the photometric properties from the SALT2 light-curve parameters as a function of spectroscopic sub-class. We discuss the non-negligible impact of host galaxy contamination on SN Ia spectral classifications, as well as investigate the accuracy of spectral template matching of the ZTF DR2 sample. We define a new subclass of underluminous SNe Ia (`04gs-like') that lie spectroscopically between normal SNe Ia and transitional 86G-like SNe Ia (stronger Si II $\lambda$5972 than normal SNe Ia but significantly weaker Ti II features than `86G-like' SNe). We model these `04gs-like' SN Ia spectra using the radiative-transfer spectral synthesis code tardis and show that cooler temperatures alone are unable to explain their spectra; some changes in elemental abundances are also required. However, the broad continuity in spectral properties seen from bright (`91T-like') to faint normal SN Ia, including the transitional and 91bg-like SNe Ia, suggests that variations within a single explosion model may be able to explain their behaviour.

The unprecedented statistics of detected Type Ia supernovae (SNe Ia) brought by the Zwicky Transient Facility enables us to probe the impact of the Large-Scale Structure on the properties of these objects. The goal of this paper is to explore the possible impact of the under-dense part of the large-scale structure on the intrinsic SALT2 light curve properties of SNe Ia and uncover possible biases in SN Ia analyses. With a volume-limited selection of ZTF-Cosmo-DR2 Type Ia supernovae overlapping with the SDSS-DR7 survey footprint, we investigate the distribution of their properties with regard to voids detected in the SDSS-DR7 galaxy sample. We further use Voronoi volumes as proxy for local density environments within the large-scale structure. We find a moderate dependency of the stretch toward the localisation around the void centre and none when considering colour. The local Voronoi volumes mostly affect the fraction of low/high stretch supernovae. With the current statistics available, we consider that the impact of high or low local density environment can be considered as a proxy for the colour of the host galaxy. Under-dense environments should not cause any biases in supernova analyses., Comment: 10 pages, 8 figures. Submitted to A&A

To get distances, Type Ia Supernovae magnitudes are corrected for their correlation with lightcurve width and colour. Here we investigate how this standardisation is affected by the SN environment, with the aim to reduce scatter and improve standardisation. We first study the SN Ia stretch distribution, as well as its dependence on environment, as characterised by local and global (g-z) colour and stellar mass. We then look at the standardisation parameter $\alpha$, which accounts for the correlation between residuals and stretch, along with its environment dependence and linearity. We finally compute magnitude offsets between SNe in different astrophysical environments after colour and stretch standardisation, aka steps. This analysis is made possible due to the unprecedented statistics of the ZTF SN Ia DR2 volume-limited sample. The stretch distribution exhibits a bimodal behaviour, as previously found in literature. However, we find the distribution means to decrease with host stellar mass at a 9.0$\sigma$ significance. We demonstrate, at the 14.3$\sigma$ level, that the stretch-magnitude relation is non-linear, challenging the usual linear stretch-residuals relation. Fitting for a broken-$\alpha$ model, we indeed find two different slopes between stretch regimes ($x_1<-0.49\pm0.06$): $\alpha_{low}=0.28\pm0.01$ and $\alpha_{high}=0.09\pm0.01$, a $\Delta_{\alpha}=-0.19\pm0.01$ difference. As the relative proportion of SNe Ia in the high-/low-stretch modes evolves with redshift and environment, this implies that a linear $\alpha$ also evolves with redshift and environment. Concerning the environmental magnitude offset $\gamma$, we find it to be greater than 0.14 mag regardless of the considered environmental tracer used (local or global colour and stellar mass), all measured at the $\geq 6\sigma$ level, increased to $\sim0.18\pm0.01$ mag when accounting for the stretch-non linearity., Comment: 15 pages, 9 figures, submitted to Astronomy and Astrophysics

Obtaining spectroscopic observations of the progenitors of core-collapse supernovae is often unfeasible due to an inherent lack of knowledge as to which stars will go supernova and when they will explode. In this letter, we present photometric and spectroscopic observations of the progenitor activity of SN 2023fyq in the preceding 150 days before the He-rich progenitor exploded as a Type Ibn supernova. The progenitor of SN 2023fyq shows an exponential rise in flux prior to core-collapse. Complex He I emission line features are observed, with a P-Cygni like profile, as well as an evolving broad base with velocities on the order of 10,000 km/s, possibly due to electron scattering. The luminosity and evolution of SN 2023fyq are consistent with a faint Type Ibn, reaching a peak r-band magnitude of 18.1 mag, although there is some uncertainty in the distance to the host, NGC 4388, located in the Virgo cluster. We present additional evidence of asymmetric He-rich material being present prior to the explosion of SN 2023fyq, as well as after, suggesting this material has survived the ejecta-CSM interaction. Broad [O I] and the Ca II triplet lines are observed at late phases, confirming that SN 2023fyq was a genuine supernova rather than a non-terminal interacting transient. SN 2023fyq provides insight into the final moments of a massive star's life, highlighting that the progenitor is likely highly unstable before core-collapse., Comment: 7 Pages, 5 Figures, accepted to A&A Letters

Rapid identification of the optical counterparts of Neutron Star (NS) merger events discovered by gravitational wave detectors may require observing a large error region and sifting through a large number of transients to identify the object of interest. Given the expense of spectroscopic observations, a question arises: How can we utilize photometric observations for candidate prioritization, and what kinds of photometric observations are needed to achieve this goal? NS merger kilonova exhibits low ejecta mass (~5x10^-2 solar mass) and a rapidly evolving photospheric radius (with a velocity ~0.2c). As a consequence, these sources display rapid optical-flux evolution. Indeed, selection based on fast flux variations is commonly used for young supernovae and NS mergers. In this study, we leverage the best currently available flux-limited transient survey - the Zwicky Transient Facility Bright Transient Survey - to extend and quantify this approach. We focus on selecting transients detected in a 3-day cadence survey and observed at a one-day cadence. We explore their distribution in the phase space defined by g-r, g-dot, and r-dot. Our analysis demonstrates that for a significant portion of the time during the first week, the kilonova AT 2017gfo stands out in this phase space. It is important to note that this investigation is subject to various biases and challenges; nevertheless, it suggests that certain photometric observations can be leveraged to identify transients with the highest probability of being fast-evolving events. We also find that a large fraction (~0.75) of the transient candidates with |g-dot|>0.7 mag/day, are cataclysmic variables or active galactic nuclei with radio counterparts., Comment: 9 pages, 4 figures, MNRAS in press

Observing a supernova explosion shortly after it occurs can reveal important information about the physics of stellar explosions and the nature of the progenitor stars of supernovae (SNe). When a star with a well-defined edge explodes in vacuum, the first photons to escape from its surface appear as a brief shock-breakout flare. The duration of this flare can extend to at most a few hours even for nonspherical breakouts from supergiant stars, after which the explosion ejecta should expand and cool. Alternatively, for stars exploding within a distribution of sufficiently dense optically thick circumstellar material, the first photons escape from the material beyond the stellar edge, and the duration of the initial flare can extend to several days, during which the escaping emission indicates photospheric heating. The difficulty in detecting SN explosions promptly after the event has so far limited data regarding supergiant stellar explosions mostly to serendipitous observations that, owing to the lack of ultraviolet (UV) data, were unable to determine whether the early emission is heating or cooling, and hence the nature of the early explosion event. Here, we report observations of SN 2023ixf in the nearby galaxy M101, covering the early days of the event. Using UV spectroscopy from the Hubble Space Telescope (HST) as well as a comprehensive set of additional multiwavelength observations, we trace the photometric and spectroscopic evolution of the event and are able to temporally resolve the emergence and evolution of the SN emission.

We report on our study of supernova (SN) 2022xxf based on observations obtained during the first four months of its evolution. The light curves (LCs) display two humps of similar maximum brightness separated by 75 days, unprecedented for a broad-lined (BL) Type Ic supernova (SN IcBL). SN 2022xxf is the most nearby SN IcBL to date (in NGC 3705, $z = 0.0037$, at a distance of about 20 Mpc). Optical and near-infrared photometry and spectroscopy are used to identify the energy source powering the LC. Nearly 50 epochs of high signal-to-noise-ratio spectroscopy were obtained within 130 days, comprising an unparalleled dataset for a SN IcBL, and one of the best-sampled SN datasets to date. The global spectral appearance and evolution of SN 2022xxf points to typical SN Ic/IcBL, with broad features (up to $\sim14000$ km s$^{-1}$) and a gradual transition from the photospheric to the nebular phase. However, narrow emission lines (corresponding to $\sim1000-2500$ km s$^{-1}$) are present in the spectra from the time of the second rise, suggesting slower-moving circumstellar material (CSM). These lines are subtle, in comparison to the typical strong narrow lines of CSM-interacting SNe, for example, Type IIn, Ibn, and Icn, but some are readily noticeable at late times such as in Mg I $\lambda$5170 and [O I] $\lambda$5577. Unusually, the near-infrared spectra show narrow line peaks in a number of features formed by ions of O and Mg. We infer the presence of CSM that is free of H and He. We propose that the radiative energy from the ejecta-CSM interaction is a plausible explanation for the second LC hump. This interaction scenario is supported by the color evolution, which progresses to the blue as the light curve evolves along the second hump, and the slow second rise and subsequent rapid LC drop. (Abstract abridged), Comment: Accepted version

We present and analyze observations of the Type Ib supernova (SN) 2019odp (a.k.a ZTF19abqwtfu) covering epochs within days of the explosion to the nebular phase at 360 d post-explosion. We discuss them in the context of recombination cooling emission for the early excess emission and consider progenitor models based on the nebular phase spectra. Our observations include photometric observations mainly in the optical and low to medium-resolution spectroscopic observations covering the complete observable time-range. We expand on existing methods to derive oxygen mass estimates from nebular phase spectroscopy. Our spectroscopic observations confirm the presence of He in the SN ejecta and we thus (re)classify it as a Type Ib supernova. From the pseudo-bolometric lightcurve we estimate a high ejecta mass $M_\text{ej} \sim 4 - 7~M_\odot$. The high ejecta mass, large nebular [O I]/[Ca II] line flux ratio ($1.2 - 1.9$) and an oxygen mass above $\gtrapprox 0.5\, M_\odot$ point towards a progenitor with pre-explosion mass higher than $18\,M_\odot$. The compact nature of the progenitor ($\lesssim 10\,R_\odot$) suggests a Wolf-Rayet (WR) star as progenitor., Comment: Submitted to A&A

If high-mass stars (>20-25 Msun) are the progenitors of stripped-envelope (SE) supernovae (SNe), their massive ejecta should lead to broad, long-duration lightcurves (LCs). Instead, literature samples of SE~SNe have reported relatively narrow LCs with ejecta masses between 1-4 Msun that favor progenitors <20-25 Msun. Working with the untargeted sample of (i)PTF SNe to better constrain their rates, we search for SE~SNe with broad LCs. Using a simple LC stretch compared to a template to measure broadness, we identified eight significantly broader Type~Ibc SNe after applying quantitative sample selection criteria. The LCs, colors, and spectra of these SNe are found to evolve more slowly relative to typical Type~Ibc SNe, proportional with the stretch. Bolometric LC modeling and their nebular spectra indicate high ejecta and nickel masses, assuming radioactive decay powering. Additionally, these objects are preferentially located in low-metallicity host galaxies with high star-formation rates, which may account for their massive progenitors, as well as their relative absence from the literature. Our study thus supports the link between broad LCs (as measured by stretch) and high-mass progenitor stars in SE~SNe with independent evidence from bolometric LC modeling, nebular spectra, host environment properties, and photometric evolution. In the first systematic search of its kind using an untargeted sample, we use the stretch distribution to identify a higher than previously appreciated fraction of SE~SNe with broad LCs (~13%). Correcting for Malmquist and LC duration observational biases, we conservatively estimate that a minimum of ~6% of SE~SNe are consistent with high-mass progenitors. This result has implications for the progenitor channels of SE~SNe, late stages of massive stellar evolution, oxygen fraction in the universe, and formation channels for stellar-mass black holes.(Abridged), Comment: 38 pages, 24 figures, 12 tables. Submitted to A&A

SN\,2020qlb (ZTF20abobpcb) is a hydrogen-poor superluminous supernova (SLSN-I) that is among the most luminous (maximum M$_{g} = -22.25$ mag) and that has one of the longest rise times (77 days from explosion to maximum). We estimate the total radiated energy to be $>2.1\times10^{51}$ erg. SN\,2020qlb has a well-sampled light curve that exhibits clear near and post peak undulations, a phenomenon seen in other SLSNe, whose physical origin is still unknown. We discuss the potential power source of this immense explosion as well as the mechanisms behind its observed light curve undulations. We analyze photospheric spectra and compare them to other SLSNe-I. We constructed the bolometric light curve using photometry from a large data set of observations from the Zwicky Transient Facility (ZTF), Liverpool Telescope (LT), and Neil Gehrels Swift Observatory and compare it with radioactive, circumstellar interaction and magnetar models. Model residuals and light curve polynomial fit residuals are analyzed to estimate the undulation timescale and amplitude. We also determine host galaxy properties based on imaging and spectroscopy data, including a detection of the [O III]$\lambda$4363, auroral line, allowing for a direct metallicity measurement. We rule out the Arnett $^{56}$Ni decay model for SN\,2020qlb's light curve due to unphysical parameter results. Our most favored power source is the magnetic dipole spin-down energy deposition of a magnetar. Two to three near peak oscillations, intriguingly similar to those of SN\,2015bn, were found in the magnetar model residuals with a timescale of $32\pm6$ days and an amplitude of 6$\%$ of peak luminosity. We rule out centrally located undulation sources due to timescale considerations; and we favor the result of ejecta interactions with circumstellar material (CSM) density fluctuations as the source of the undulations., Comment: 22 pages, 25 figures, submitted to A&A

In the new era of time-domain surveys Type Ia supernovae are being caught sooner after explosion, which has exposed significant variation in their early light curves. Two driving factors for early time evolution are the distribution of nickel in the ejecta and the presence of flux excesses of various causes. We perform an analysis of the largest young SN Ia sample to date. We compare 115 SN Ia light curves from the Zwicky Transient Facility to the turtls model grid containing light curves of Chandrasekhar-mass explosions with a range of nickel masses, nickel distributions and explosion energies. We find that the majority of our observed light curves are well reproduced by Chandrasekhar-mass explosion models with a preference for highly extended nickel distributions. We identify six SNe Ia with an early-time flux excess in our g- and r-band data (four `blue' and two `red' flux excesses). We find an intrinsic rate of 18+/-11 per cent of early flux excesses in SNe Ia at z < 0.07, based on three detected flux excesses out of 30 (10 per cent) observed SNe Ia with a simulated efficiency of 57 per cent. This is comparable to rates of flux excesses in the literature but also accounts for detection efficiencies. Two of these events are mostly consistent with CSM interaction, while the other four have longer lifetimes in agreement with companion interaction and nickel-clump models. We find a higher frequency of flux excesses in 91T/99aa-like events (44+/-13 per cent).

Supernova (SN) siblings -- two or more SNe in the same parent galaxy -- are useful tools for exploring progenitor stellar populations as well as properties of the host galaxies such as distance, star formation rate, dust extinction, and metallicity. Since the average SN rate for a Milky Way-type galaxy is just one per century, a large imaging survey is required to discover an appreciable sample of SN siblings. From the wide-field Zwicky Transient Facility (ZTF) Bright Transient Survey (BTS; which aims for spectroscopic completeness for all transients which peak brighter than $r{<}$18.5 mag) we present 10 SN siblings in 5 parent galaxies. For each of these families we analyze the SN's location within the host and its underlying stellar population, finding agreement with expectations that SNe from more massive progenitors are found nearer to their host core and in regions of more active star formation. We also present an analysis of the relative rates of core collapse and thermonuclear SN siblings, finding a significantly lower ratio than past SN sibling samples due to the unbiased nature of the ZTF., Comment: 14 pages, 4 tables, 7 figures, accepted to MNRAS

Decades after their discovery, only four hydrogen-deficient carbon (HdC) stars were known to have no circumstellar dust shell. This is in complete contrast to the $\sim$130 known Galactic HdC stars that are notorious for being heavy dust producers, i.e. the R Coronae Borealis (RCB) stars. Together they form a rare class of supergiant stars that are thought to originate from the merger of CO/He white dwarf binary systems, otherwise known as the double-degenerate scenario. We searched for new dustless HdC (dLHdC) stars. We used primarily the 2MASS and GAIA eDR3 catalogues to select candidates that were followed-up spectroscopically. We discovered 27 new dLHdC stars, one new RCB star and two new EHe stars. Surprisingly, 20 of the new dLHdC stars share a characteristic of the known dLHdC star HD 148839, having lower atmospheric hydrogen deficiencies. The uncovered population of dLHdC stars exhibit a Bulge-like distribution, like the RCB stars, but show multiple differences from those that indicate they are a different population of HdC stars following its own evolutionary sequence with a fainter luminosity and also a narrow range of effective temperature. We found indication of a current low dust production activity for four of the new dLHdC stars which could be typical RCB stars passing through a transition time. We have evidence for the first time of a large range of absolute magnitudes in the overall population of HdC stars, spanning over 3 mag. In the favoured formation framework, this is explained by a large range in the initial total WD binary mass which leads to a series of evolutionary sequences with distinct maximum brightness and initial temperature. The cold Galactic RCB stars are also noticeably fainter than the Magellanic ones, possibly due to a difference in metallicity resulting in different WD mass ratio. In our Galaxy, there could be as many dLHdC stars as RCB stars., Comment: 22 pages, 18 figures, accepted in A&A. V3: minor corrections on Table 2 and 3 related to HD stars names and (V-I)0 colours (Fig.10 changed accordingly). V4: changed the colour to light blue for the two new EHe stars indicated in Fig.10, page 13

The final explosive fate of massive stars, and the nature of the compact remnants they leave behind (black holes and neutron stars), are major open questions in astrophysics. Many massive stars are stripped of their outer hydrogen envelopes as they evolve. Such Wolf-Rayet (W-R) stars emit strong and rapidly expanding (v_wind>1000 km/s) winds indicating a high escape velocity from the stellar surface. A fraction of this population is also helium depleted, with spectra dominated by highly-ionized emission lines of carbon and oxygen (Types WC/WO). Evidence indicates that the most commonly-observed supernova (SN) explosions that lack hydrogen and helium (Types Ib/Ic) cannot result from massive WC/WO stars, leading some to suggest that most such stars collapse directly into black holes without a visible supernova explosions. Here, we present observations of supernova SN 2019hgp, discovered about a day after explosion. The short rise time and rapid decline place it among an emerging population of rapidly-evolving transients (RETs). Spectroscopy reveals a rich set of emission lines indicating that the explosion occurred within a nebula composed of carbon, oxygen, and neon. Narrow absorption features show that this material is expanding at relatively high velocities (>1500 km/s) requiring a compact progenitor. Our observations are consistent with an explosion of a massive WC/WO star, and suggest that massive W-R stars may be the progenitors of some rapidly evolving transients., Comment: Unedited author version, Nature in press

We present 75 near-infrared (NIR; 0.8$-$2.5 $\mu$m) spectra of 34 stripped-envelope core-collapse supernovae (SESNe) obtained by the Carnegie Supernova Project-II (CSP-II), encompassing optical spectroscopic Types IIb, Ib, Ic, and Ic-BL. The spectra range in phase from pre-maximum to 80 days past maximum. This unique data set constitutes the largest NIR spectroscopic sample of SESNe to date. NIR spectroscopy provides observables with additional information that is not available in the optical. Specifically, the NIR contains the resonance lines of He I and allows a more detailed look at whether Type Ic supernovae are completely stripped of their outer He layer. The NIR spectra of SESNe have broad similarities, but closer examination through statistical means reveals a strong dichotomy between NIR "He-rich" and "He-poor" SNe. These NIR subgroups correspond almost perfectly to the optical IIb/Ib and Ic/Ic-BL types, respectively. The largest difference between the two groups is observed in the 2 $\mu$m region, near the He I $\lambda$2.0581 $\mu$m line. The division between the two groups is not an arbitrary one along a continuous sequence. Early spectra of He-rich SESNe show much stronger He I $\lambda$2.0581 $\mu$m absorption compared to the He-poor group, but with a wide range of profile shapes. The same line also provides evidence for trace amounts of He in half of our SNe in the He-poor group.

Type Ia supernovae (SNe~Ia) in the nearby Hubble flow are excellent distance indicators in cosmology. The Zwicky Transient Facility (ZTF) has observed a large sample of supernovae from an untargeted, rolling survey, reaching $20.8, 20.6, 20.3$ mag in $g$ $r$, and $i$-band, respectively. With a FoV of 47 sq.deg, ZTF discovered $>$ 3000 SNe~Ia in a little over 2.5 years. Here, we report on the sample of 761 spectroscopically classified SNe~Ia from the first year of operations (DR1). The sample has a median redshift $\bar z =$ 0.057, nearly a factor of two higher than the current low-$z$ sample. Our sample has a total of 934 spectra, of which 632 were obtained with the robotic SEDm on Palomar P60. We assess the potential for precision cosmology for a total of 305 SNe with redshifts from host galaxy spectra. The sample is already comparable in size to the entire combined literature low-$z$ anchor sample. The median first detection is 13.5 days before maximum light, about 10 days earlier than the median in the literature. Furthermore, six SNe from our sample are at $D_L < 80$ Mpc, for which host galaxy distances can be obtained in the JWST era, such that we have calibrator and Hubble flow SNe observed with the same instrument. In the entire duration of ZTF-I, we have observed nearly fifty SNe for which we can obtain calibrator distances, key for percent level distance scale measurements., Comment: 16 pages, 9 figures, accepted for publication in MNRAS. Data and associated information for the sample can be found at https://github.com/ZwickyTransientFacility/ztfcosmodr

We report on the discovery of AT2018lqh (ZTF18abfzgpl) -- a rapidly-evolving extra-galactic transient in a star-forming host at 242 Mpc. The transient g-band light curve's duration above half-maximum light is about 2.1 days, where 0.4/1.7 days are spent on the rise/decay, respectively. The estimated bolometric light curve of this object peaked at about 7x10^42 erg/s -- roughly seven times brighter than AT2017gfo. We show that this event can be explained by an explosion with a fast (v~0.08 c) low-mass (~0.07 Msun) ejecta, composed mostly of radioactive elements. For example, ejecta dominated by Ni-56 with a time scale of t_0=1.6 days for the ejecta to become optically thin for gamma-rays fits the data well. Such a scenario requires burning at densities that are typically found in the envelopes of neutron stars or the cores of white dwarfs. A combination of circumstellar material (CSM) interaction power at early times and shock cooling at late times is consistent with the photometric observations, but the observed spectrum of the event may pose some challenges for this scenario. The observations are not consistent with a shock breakout from a stellar envelope, while a model involving a low-mass ejecta ramming into low-mass CSM cannot explain both the early- and late-time observations., Comment: 10 pages, ApJ in press

We present a detailed analysis of SN 2020qmp, a nearby type IIP core-collapse supernova (CCSN), discovered by the Palomar Gattini-IR (PGIR) survey in the galaxy UGC07125. We illustrate how the multiwavelength study of this event helps our general understanding of stellar progenitors and circumstellar medium (CSM) interactions in CCSNe. We also highlight the importance of near-infrared (NIR) surveys for early detections of SNe in dusty environments. SN 2020qmp displays characteristic hydrogen lines in its optical spectra, as well as a plateau in its optical LC, hallmarks of a type IIP SN. We do not detect linear polarization during the plateau phase, with a 3 sigma upper limit of 0.78%. Through hydrodynamical LC modeling and an analysis of its nebular spectra, we estimate a progenitor mass of around 11 solar masses, and an explosion energy of around 0.8e51 erg. We find that the spectral energy distribution cannot be explained by a simple CSM interaction model, assuming a constant shock velocity and steady mass-loss rate, and the excess X-ray luminosity compared with the synchrotron radio luminosity suggests deviations from equipartition. Finally, we demonstrate the advantages of NIR surveys over optical surveys for the detection of dust-obscured CCSNe in the local universe. Specifically, our simulations show that the Wide-Field Infrared Transient Explorer will detect about 14 more CCSNe out of 75 expected in its footprint within 40 Mpc, over five years than an optical survey equivalent to the Zwicky Transient Facility would detect. We have determined or constrained the main properties of SN 2020qmp and of its progenitor, highlighting the value of multiwavelength follow-up observations of nearby CCSNe. We have also shown that forthcoming NIR surveys will finally enable us to do a nearly complete census of CCSNe in the local universe., Comment: 15 pages, 14 figures, Accepted in Astronomy & Astrophysics; Revised version after referee comments and language edits from journal

We present observations and analysis of SN 2020cxd, a Low luminous (LL), long-lived Type IIP SN. This object was a clear outlier in the magnitude-limited SN sample recently presented by the ZTF Bright Transient Survey. We demonstrate that SN 2020cxd is an additional member of the group of LL SNe, and discuss the rarity of LL SNe in the context of the ZTF survey, and how further studies of these faintest members of the CC SN family might help understand the underlying initial mass function for stars that explode., Comment: 22 pages, 16 figures. This is the version resubmitted to A&A with replies to referee report

We report a discovery of a new long-period X-ray pulsar SRGA J204318.2+443815/SRGe J204319.0+443820 in the Be binary system. The source was found in the second all-sky survey by the Mikhail Pavlinsky telescope on board the SRG mission. The follow-up observations with XMM-Newton, NICER and NuSTAR observatories allowed us to discover a strong coherent signal in the source light curve with the period of $\sim742$ s. The pulsed fraction was found to depend on the energy increasing from $\sim20$% in soft X-rays to $>50$% at high energies, as it is typical for X-ray pulsars. The source demonstrate a quite hard spectrum with an exponential cutoff at high energies and bolometric luminosity of $L_X \simeq 4\times10^{35}$ erg/s. Dedicated optical and infrared observations with the RTT-150, NOT, Keck and Palomar telescopes revealed a number of emission lines (H$_{\alpha}$, HeI, Pashen and Braket series) with the strongly absorbed continuum. All of above suggests that SRGAJ204318.2+443815/ SRGeJ204319.0+443820 is a new persistent low luminosity X-ray pulsar in a distant binary system with a Be-star of the B0-B2e class. Thus the SRG observatory allow us to unveil the hidden population of faint persistent objects including the population of slowly rotating X-ray pulsars in Be systems., Comment: Submitted to A&A for the Special Issue: The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission. 9 pages, 9 figures

We present optical and near-infrared (NIR, $YJH$-band) observations of 42 Type Ia supernovae (SNe Ia) discovered by the untargeted intermediate Palomar Transient Factory (iPTF) survey. This new data-set covers a broad range of redshifts and host galaxy stellar masses, compared to previous SN Ia efforts in the NIR. We construct a sample, using also literature data at optical and NIR wavelengths, to examine claimed correlations between the host stellar masses and the Hubble diagram residuals. The SN magnitudes are corrected for host galaxy extinction using either a global total-to-selective extinction ratio, $R_V$=2.0 for all SNe, or a best-fit $R_V$ for each SN individually. Unlike previous studies which were based on a narrower range in host stellar mass, we do not find evidence for a "mass-step", between the color- and stretch-corrected peak $J$ and $H$ magnitudes for galaxies below and above $\log(M_{*}/M_{\odot}) = 10$. However, the mass-step remains significant ($3\sigma$) at optical wavelengths ($g,r,i$) when using a global $R_V$, but vanishes when each SN is corrected using their individual best-fit $R_V$. Our study confirms the benefits of the NIR SN Ia distance estimates, as these are largely exempted from the empirical corrections dominating the systematic uncertainties in the optical., Comment: 20 pages, 13 figures. Submitted to ApJ. Comments are welcome!

Cosmological Gamma-Ray Bursts (GRBs) are known to arise from distinct progenitor channels: short GRBs mostly from neutron star mergers and long GRBs from a rare type of core-collapse supernova (CCSN) called collapsars. Highly magnetized neutron stars called magnetars also generate energetic, short-duration gamma-ray transients called Magnetar Giant Flares (MGFs). Three have been observed from the Milky Way and its satellite galaxies and they have long been suspected to contribute a third class of extragalactic GRBs. We report the unambiguous identification of a distinct population of 4 local ($<$5 Mpc) short GRBs, adding GRB 070222 to previously discussed events. While identified solely based on alignment to nearby star-forming galaxies, their rise time and isotropic energy release are independently inconsistent with the larger short GRB population at $>$99.9% confidence. These properties, the host galaxies, and non-detection in gravitational waves all point to an extragalactic MGF origin. Despite the small sample, the inferred volumetric rates for events above $4\times10^{44}$ erg of $R_{MGF}=3.8_{-3.1}^{+4.0}\times10^5$ Gpc$^{-3}$ yr$^{-1}$ place MGFs as the dominant gamma-ray transient detected from extragalactic sources. As previously suggested, these rates imply that some magnetars produce multiple MGFs, providing a source of repeating GRBs. The rates and host galaxies favor common CCSN as key progenitors of magnetars., Comment: Accepted for publication in ApJL. Updated versions fix typos in the table and updates citations to published versions

Early observations of Type Ia supernovae (SNe$\,$Ia) provide essential clues for understanding the progenitor system that gave rise to the terminal thermonuclear explosion. We present exquisite observations of SN$\,$2019yvq, the second observed SN$\,$Ia, after iPTF$\,$14atg, to display an early flash of emission in the ultraviolet (UV) and optical. Our analysis finds that SN$\,$2019yvq was unusual, even when ignoring the initial flash, in that it was moderately underluminous for an SN$\,$Ia ($M_g \approx -18.5\,$mag at peak) yet featured very high absorption velocities ($v \approx 15,000\,\mathrm{km\,s}^{-1}$ for Si II $\lambda$6355 at peak). We find that many of the observational features of SN$\,$2019yvq, aside from the flash, can be explained if the explosive yield of radioactive $^{56}\mathrm{Ni}$ is relatively low (we measure $M_{^{56}\mathrm{Ni}} = 0.31 \pm 0.05\,M_\odot$) and it and other iron-group elements are concentrated in the innermost layers of the ejecta. To explain both the UV/optical flash and peak properties of SN$\,$2019yvq we consider four different models: interaction between the SN ejecta and a nondegenerate companion, extended clumps of $^{56}\mathrm{Ni}$ in the outer ejecta, a double-detonation explosion, and the violent merger of two white dwarfs. Each of these models has shortcomings when compared to the observations; it is clear additional tuning is required to better match SN$\,$2019yvq. In closing, we predict that the nebular spectra of SN$\,$2019yvq will feature either H or He emission, if the ejecta collided with a companion, strong [Ca II] emission, if it was a double detonation, or narrow [O I] emission, if it was due to a violent merger., Comment: 27 pages, 14 figures, accepted in ApJ

We present optical, radio, and X-ray observations of SN2020bvc (=ASASSN20bs; ZTF20aalxlis), a nearby ($z=0.0252$; $d$=114 Mpc) broad-lined (BL) Type Ic supernova (SN). Our observations show that SN2020bvc shares several properties in common with the Ic-BL SN2006aj, which was associated with the low-luminosity gamma-ray burst (LLGRB) 060218. First, the 10 GHz radio light curve is on the faint end of LLGRB-SNe ($L_\mathrm{radio} \approx 10^{37}$erg/s): we model our VLA observations (spanning 13-43 d) as synchrotron emission from a mildly relativistic ($v \gtrsim 0.3c$) forward shock. Second, with Swift and Chandra we detect X-ray emission ($L_X \approx 10^{41}$erg/s) that is not naturally explained as inverse Compton emission or as part of the same synchrotron spectrum as the radio emission. Third, high-cadence ($6\times$/night) data from the Zwicky Transient Facility (ZTF) shows a double-peaked optical light curve, the first peak from shock-cooling emission from extended low-mass material (mass $M<10^{-2} M_\odot$ at radius $R>10^{12}$cm) and the second peak from the radioactive decay of Ni-56. SN2020bvc is the first confirmed double-peaked Ic-BL SN discovered without a GRB trigger, and shows X-ray and radio emission similar to LLGRB-SNe: this is consistent with models in which the same mechanism produces both the LLGRB and the shock-cooling emission. For four of the five other nearby ($z\lesssim0.05$) Ic-BL SNe with ZTF high-cadence data, we rule out a first peak like that seen in SN2006aj and SN2020bvc, i.e. that lasts $\approx 1$d and reaches a peak luminosity $M \approx -18$. X-ray and radio follow-up observations of future such events will establish whether double-peaked optical light curves are indeed predictive of LLGRB-like X-ray and radio emission., Comment: Accepted to ApJ. 39 pages, 17 figures. Comments welcome

We report the results from spectroscopic observations of the multiple images of the strongly lensed Type Ia supernova (SN Ia), iPTF16geu, obtained with ground based telescopes and the Hubble Space Telescope (HST). From a single epoch of slitless spectroscopy with HST, we can resolve spectra of individual lensed supernova images for the first time. This allows us to perform an independent measurement of the time-delay between the two brightest images, $\Delta t = 1.4 \pm 5.0$ days, which is consistent with the time-delay measured from the light-curves. We also present measurements of narrow emission and absorption lines characterizing the interstellar medium in the host galaxy at z=0.4087, as well as in the foreground lensing galaxy at z=0.2163. We detect strong Na ID absorption in the host galaxy, indicating that iPTF16geu belongs to a subclass of SNe Ia displaying "anomalously" large Na ID column densities in comparison to the amount of dust extinction derived from their light curves. For the deflecting galaxy, we refine the measurement of the velocity dispersion, $\sigma = 129 \pm 4$ km/s, which significantly constrains the lens model. Since the time-delay between the SN images is negligible, we can use unresolved ground based spectroscopy, boosted by a factor ~70 from lensing magnification, to study the properties of a high-z SN Ia with unprecedented signal-to-noise ratio. The spectral properties of the supernova, such as pseudo-Equivalent widths of several absorption features and velocities of the Si II-line indicate that iPTF16geu, besides being lensed, is a normal SN Ia, indistinguishable from well-studied ones in the local universe, providing support for the use of SNe Ia in precision cosmology. We do not detect any significant deviations of the SN spectral energy distribution from microlensing of the SN photosphere by stars and compact objects in the lensing galaxy., Comment: 10 pages, 10 figures. Comments are welcome!

AT2016dah and AT2017fyp are fairly typical Andromeda Galaxy (M31) classical novae. AT2016dah is an almost text book example of a 'very fast' declining, yet uncommon, Fe II'b' (broad-lined) nova, discovered during the rise to peak optical luminosity, and decaying with a smooth broken power-law light curve. AT2017fyp is classed as a 'fast' nova, unusually for M31, its early decline spectrum simultaneously shows properties of both Fe II and He/N spectral types - a 'hybrid'. Similarly, the light curve of AT2017fyp has a broken power-law decline but exhibits an extended flat-topped maximum. Both novae were followed in the UV and X-ray by the Neil Gehrels Swift Observatory, but no X-ray source was detected for either nova. The pair were followed photometrically and spectroscopically into their nebular phases. The progenitor systems were not visible in archival optical data, implying that the mass donors are main sequence stars. What makes AT2016dah and AT2017fyp particularly interesting is their position with respect to M31. The pair are close on the sky but are located far from the centre of M31, lying almost along the semi-minor axis of their host. Radial velocity measurements and simulations of the M31 nova population leads to the conclusion that both novae are members of the Andromeda Giant Stellar Stream (GSS). We find the probability of at least two M31 novae appearing coincident with the GSS by chance is ~1%. Therefore, we claim that these novae arose from the GSS progenitor, not M31 - the first confirmed novae discovered in a tidal steam., Comment: 22 pages, 14 figures, 4 tables. Accepted for publication in MNRAS

We present $81$ near-infrared (NIR) spectra of $30$ Type II supernovae (SNe II) from the Carnegie Supernova Project-II (CSP-II), the largest such dataset published to date. We identify a number of NIR features and characterize their evolution over time. The NIR spectroscopic properties of SNe II fall into two distinct groups. This classification is first based on the strength of the He I $\lambda1.083\,\mu$m absorption during the plateau phase; SNe II are either significantly above (spectroscopically strong) or below $50$ angstroms (spectroscopically weak) in pseudo equivalent width. However between the two groups, other properties, such as the timing of CO formation and the presence of Sr II, are also observed. Most surprisingly, the distinct weak and strong NIR spectroscopic classes correspond to SNe II with slow and fast declining light curves, respectively. These two photometric groups match the modern nomenclature of SNe IIP and IIL. Including NIR spectra previously published, 18 out of 19 SNe II follow this slow declining-spectroscopically weak and fast declining-spectroscopically strong correspondence. This is in apparent contradiction to the recent findings in the optical that slow and fast decliners show a continuous distribution of properties. The weak SNe II show a high-velocity component of helium that may be caused by a thermal excitation from a reverse-shock created by the outer ejecta interacting with the red supergiant wind, but the origin of the observed dichotomy is not understood. Further studies are crucial in determining whether the apparent differences in the NIR are due to distinct physical processes or a gap in the current data set., Comment: 32 pages, 17 figures, accepted to ApJ

We present photometry and spectroscopy of four hydrogen-poor luminous supernovae discovered during the two-month science commissioning and early operations of the Zwicky Transient Facility (ZTF) survey. Three of these objects, SN2018bym (ZTF18aapgrxo), SN2018avk (ZTF18aaisyyp) and SN2018bgv (ZTF18aavrmcg) resemble typical SLSN-I spectroscopically, while SN2018don (ZTF18aajqcue) may be an object similar to SN2007bi experiencing considerable host galaxy reddening, or an intrinsically long-lived, luminous and red SN Ic. We analyze the light curves, spectra, and host galaxy properties of these four objects and put them in context of the population of SLSN-I. SN2018bgv stands out as the fastest-rising SLSN-I observed to date, with a rest-frame g-band rise time of just 10 days from explosion to peak -- if it is powered by magnetar spin-down, the implied ejecta mass is only ~1 M$_{\odot}$. SN2018don also displays unusual properties -- in addition to its red colors and comparatively massive host galaxy, the light curve undergoes some of the strongest light curve undulations post-peak seen in a SLSN-I, which we speculate may be due to interaction with circumstellar material. We discuss the promises and challenges of finding SLSNe in large-scale surveys like ZTF given the observed diversity in the population., Comment: Updated to match published version. Expanded discussion following referee comments; conclusions unchanged

We present a sample of supernovae Type IIn (SNe IIn) from the untargeted, magnitude-limited surveys of the Palomar Transient Factory (PTF) and its successor, the intermediate PTF (iPTF). The SNe IIn found and followed by the PTF/iPTF were used to select a sample of 42 events with useful constraints on the rise times as well as with available post-peak photometry. The sample SNe were discovered in 2009-2016 and have at least one low-resolution classification spectrum, as well as photometry from the P48 and P60 telescopes at Palomar Observatory. We study the light-curve properties of these SNe IIn using spline fits (for the peak and the declining portion) and template matching (for the rising portion). We find that the typical rise times are divided into fast and slow risers at $20\pm6$ d and $50\pm11$ d, respectively. The decline rates are possibly divided into two clusters, but this division has weak statistical significance. We find no significant correlation between the peak luminosity of SNe IIn and their rise times, but the more luminous SNe IIn are generally found to be more long-lasting. Slowly rising SNe IIn are generally found to decline slowly. The SNe in our sample were hosted by galaxies of absolute magnitude $-22 \lesssim M_g \lesssim -13$ mag. The K-corrections at light-curve peak of the SNe IIn in our sample are found to be within 0.2 mag for the observer's frame $r$-band, for SNe at redshifts $z < 0.25$. By applying K-corrections and also including ostensibly "superluminous" SNe IIn, we find that the peak magnitudes are $M_{\rm peak}^{r} = -19.18\pm1.32$ mag. We conclude that the occurrence of conspicuous light-curve bumps in SNe IIn, such as in iPTF13z, are limited to $1.4^{+14.6}_{-1.0} \%$ of the SNe IIn. We also investigate a possible sub-type of SNe IIn with a fast rise to a $\gtrsim 50$ d plateau followed by a slow, linear decline., Comment: Matches journal version. Table 2 on CDS. Durations in Fig. 13 corrected, conclusions unchanged. Abstract abridged. 33 pages, 22 figures, 7 tables

We present the volumetric rate of normal type Ia supernovae (SNe Ia) discovered by the Palomar Transient Factory (PTF). Using strict data-quality cuts, and considering only periods when the PTF maintained a regular cadence, PTF discovered 90 SNe Ia at $z\le0.09$ in a well-controlled sample over three years of operation (2010-2012). We use this to calculate the volumetric rate of SN Ia events by comparing this sample to simulations of hundreds of millions of SN Ia light curves produced in statistically representative realisations of the PTF survey. This quantifies the recovery efficiency of each PTF SN Ia event, and thus the relative weighting of each event. From this, the volumetric SN Ia rate was found to be $r_v=2.43\pm0.29\,\text{(stat)}_{-0.19}^{+0.33}\text{(sys)}\times10^{-5}\,\text{SNe yr}^{-1}\,\text{Mpc}^{-3}\, h_{70}^{3}$. This represents the most precise local measurement of the SN Ia rate. We fit a simple SN Ia delay-time distribution model, $\propto\mathrm{t}^{-\beta}$, to our PTF rate measurement combined with a literature sample of rate measurements from surveys at higher-redshifts. We find $\beta{\sim}1$, consistent with a progenitor channel governed by the gravitational in-spiral of binary white dwarfs., Comment: 14 pages, 8 figures

SN2010jl was a luminous Type IIn supernova (SN), detected in radio, optical, X-ray and hard X-rays. Here we report on its six year R- and g-band light curves obtained using the Palomar Transient Factory. The light curve was generated using a pipeline based on the proper image subtraction method and we discuss the algorithm performances. As noted before, the R-band light curve, up to about 300 days after maximum light is well described by a power-law decline with a power-law index of about -0.5. Between day 300 and day 2300 after maximum light, it is consistent with a power-law decline, with a power-law index of about -3.4. The longevity of the light curve suggests that the massive circum-stellar material around the progenitor was ejected on time scales of at least tens of years prior to the progenitor explosion., Comment: 6 pages, PASP in press

We present a detailed study of a transient in the center of SDSS1115+0544 based on the extensive UV, optical, mid-IR light curves (LC) and spectra over 1200 days. The host galaxy is a quiescent early type galaxy at $z$ = 0.0899 with a blackhole mass of $2\times10^7M_\odot$. The transient underwent a 2.5 magnitude brightening over $\sim120$ days, reaching a peak $V$-band luminosity (extinction corrected) of $-20.9$ magnitude, then fading 0.5 magnitude over 200 days, settling into a plateau of $>600$ days. Following the optical brightening are the significant mid-IR flares at $3.4$ and $4.5\mu$m, with a peak time delay of $\sim180$ days. The mid-IR LCs are explained as the echo of UV photons by a dust medium with a radius of $5\times10^{17}$ cm, consistent with $\rm E(B-V)$ of 0.58 inferred from the spectra. This event is very energetic with an extinction corrected $L_{bol} \sim 4\times10^{44}$ erg s$^{-1}$. Optical spectra over 400 days in the plateau phase revealed newly formed broad H$\alpha, \beta$ emission with a FWHM of $\sim3750$ km s$^{-1}$ and narrow coronal lines such as [Fe VII], [Ne V]. This flare also has a steeply rising UV continuum, detected by multi-epoch $Swift$ data at $+700$ to $+900$ days post optical peak. The broad Balmer lines and the UV continuum do not show significant temporal variations. The slow evolving LCs over 1200 days, the constant Balmer lines and UV continuum at late-times rule out TDE and SN IIn as the physical model for this event. We propose that this event is a `turn-on' AGN, transitioning from a quiescent state to a type 1 AGN with a sub-Eddington accretion rate of $0.017M_\odot$/yr. This change occurred on a very short time scale of $\sim 120- 200$ days. The discovery of such a rapid `turn-on' AGN poses challenges to accretion disk theories and may indicate such event is not extremely rare., Comment: Comments are welcome. Emails to the first author. Accepted for publication in ApJ

We describe a dynamic science portal called the GROWTH Marshal that allows time-domain astronomers to define science programs, program filters to save sources from different discovery streams, co-ordinate follow-up with various robotic or classical telescopes, analyze the panchromatic follow-up data and generate summary tables for publication. The GROWTH marshal currently serves 137 scientists, 38 science programs and 67 telescopes. Every night, in real-time, several science programs apply various customized filters to the 10^5 nightly alerts from the Zwicky Transient Facility. Here, we describe the schematic and explain the functionality of the various components of this international collaborative platform., Comment: Published in PASP Focus Issue on the Zwicky Transient Facility (doi: 10.1088/1538-3873/aafbc2). 12 Pages, 6 Figures

We present a catalog of 417 luminous infrared variable stars with periods exceeding 250 days. These were identified in 20 nearby galaxies by the ongoing SPIRITS survey with the Spitzer Space Telescope. Of these, 359 variables have $M_{[4.5]}$ (phase-weighted mean magnitudes) fainter than $-12$ and periods and luminosities consistent with previously reported variables in the Large Magellanic Cloud. However, 58 variables are more luminous than $M_{[4.5]} = -12$, including 11 that are brighter than $M_{[4.5]} = -13$ with the brightest having $M_{[4.5]} = -15.51$. Most of these bright variable sources have quasi-periods longer than 1000 days, including four over 2000 days. We suggest that the fundamental period-luminosity relationship, previously measured for the Large Magellanic Cloud, extends to much higher luminosities and longer periods in this large galaxy sample. We posit that these variables include massive AGB stars (possibly super-AGB stars), red supergiants experiencing exceptionally high mass-loss rates, and interacting binaries. We also present 3.6, 4.5, 5.8 and 8.0 $\mu$m photometric catalogs for all sources in these 20 galaxies., Comment: 18 pages, 25 figures

We present the best 265 sampled R-band light curves of spectroscopically identified Type Ia supernovae (SNe) from the Palomar Transient Factory (PTF; 2009-2012) survey and the intermediate Palomar Transient Factory (iPTF; 2013-2017). A model-independent light curve template is built from our data-set with the purpose to investigate average properties and diversity in our sample. We searched for multiple populations in the light curve properties using machine learning tools. We also utilised the long history of our light curves, up to 4000 days, to exclude any significant pre- or post- supernova flares. From the shapes of light curves we found the average rise time in the R band to be $16.8^{+0.5}_{-0.6}$ days. Although PTF/iPTF were single-band surveys, by modelling the residuals of the SNe in the Hubble-Lema\^{i}tre diagram, we estimate the average colour excess of our sample to be $<$E$($B$-$V$)> \approx 0.05(2)$ mag and thus the mean corrected peak brightness to be $M_R = -19.02\pm0.02$ $+5 \log( {\rm H}_0 [{\rm km} \cdot{\rm s}^{-1} {\rm Mpc}^{-1}]/70)$ mag with only weakly dependent on light curve shape. The intrinsic scatter is found to be $\sigma_R = 0.186 \pm 0.033$ mag for the redshift range $0.05

We present new 0.6-10 GHz observations of the binary neutron star merger GW170817 covering the period up to 300 days post-merger, taken with the Karl G. Jansky Very Large Array, the Australia Telescope Compact Array, the Giant Metrewave Radio Telescope and the MeerKAT telescope. We use these data to precisely characterize the decay phase of the late-time radio light curve. We find that the temporal decay is consistent with a power-law slope of t^-2.2, and that the transition between the power-law rise and decay is relatively sharp. Such a slope cannot be produced by a quasi-isotropic (cocoon-dominated) outflow, but is instead the classic signature of a relativistic jet. This provides strong observational evidence that GW170817 produced a successful jet, and directly demonstrates the link between binary neutron star mergers and short-hard GRBs. Using simple analytical arguments, we derive constraints on the geometry and the jet opening angle of GW170817. These results are consistent with those from our companion Very Long Baseline Interferometry (VLBI) paper, reporting superluminal motion in GW170817., Comment: 11 pages, 3 figures, 3 tables. Accepted for publication in ApJ Letters

Shifting the focus of Type Ia supernova (SN Ia) cosmology to the near-infrared (NIR) is a promising way to significantly reduce the systematic errors, as the strategy minimizes our reliance on the empirical width-luminosity relation and uncertain dust laws. Observations in the NIR are also crucial for our understanding of the origins and evolution of these events, further improving their cosmological utility. Any future experiments in the rest-frame NIR will require knowledge of the SN Ia NIR spectroscopic diversity, which is currently based on a small sample of observed spectra. Along with the accompanying paper, Phillips et al. (2018), we introduce the Carnegie Supernova Project-II (CSP-II), to follow up nearby SNe Ia in both the optical and the NIR. In particular, this paper focuses on the CSP-II NIR spectroscopy program, describing the survey strategy, instrumental setups, data reduction, sample characteristics, and future analyses on the data set. In collaboration with the Harvard-Smithsonian Center for Astrophysics (CfA) Supernova Group, we obtained 661 NIR spectra of 157 SNe Ia. Within this sample, 451 NIR spectra of 90 SNe Ia have corresponding CSP-II follow-up light curves. Such a sample will allow detailed studies of the NIR spectroscopic properties of SNe Ia, providing a different perspective on the properties of the unburned material, radioactive and stable nickel produced, progenitor magnetic fields, and searches for possible signatures of companion stars., Comment: 20 pages, 7 figures, accepted for publication in PASP

Compact neutron star binary systems are produced from binary massive stars through stellar evolution involving up to two supernova explosions. The final stages in the formation of these systems have not been directly observed. We report the discovery of iPTF 14gqr (SN 2014ft), a Type Ic supernova with a fast evolving light curve indicating an extremely low ejecta mass ($\approx 0.2$ solar masses) and low kinetic energy ($\approx 2 \times 10^{50}$ ergs). Early photometry and spectroscopy reveal evidence of shock cooling of an extended He-rich envelope, likely ejected in an intense pre-explosion mass loss episode of the progenitor. Taken together, we interpret iPTF 14gqr as evidence for ultra-stripped supernovae that form neutron stars in compact binary systems., Comment: Published in Science on October 12

Hydrogen-poor superluminous supernovae (SLSN-I) are a class of rare and energetic explosions discovered in untargeted transient surveys in the past decade. The progenitor stars and the physical mechanism behind their large radiated energies ($\sim10^{51}$ erg) are both debated, with one class of models primarily requiring a large rotational energy, while the other requires very massive progenitors to either convert kinetic energy into radiation via interaction with circumstellar material (CSM), or engender a pair-instability explosion. Observing the structure of the CSM around SLSN-I offers a powerful test of some scenarios, though direct observations are scarce. Here, we present a series of spectroscopic observations of the SLSN-I iPTF16eh, which reveal both absorption and time- and frequency-variable emission in the Mg II resonance doublet. We show that these observations are naturally explained as a resonance scattering light echo from a circumstellar shell. Modeling the evolution of the emission, we find a shell radius of 0.1 pc and velocity of 3300 km s$^{-1}$, implying the shell was ejected three decades prior to the supernova explosion. These properties match theoretical predictions of pulsational pair-instability shell ejections, and imply the progenitor had a He core mass of $\sim 50-55~{\rm M}_{\odot}$, corresponding to an initial mass of $\sim 115~{\rm M}_{\odot}$., Comment: Accepted. Fixed typo in table header, otherwise unchanged from previous version

We present an analysis of 507 spectra of 173 stripped-envelope (SE) supernovae (SNe) discovered by the untargeted Palomar Transient Factory (PTF) and intermediate PTF (iPTF) surveys. Our sample contains 55 Type IIb SNe (SNe IIb), 45 Type Ib SNe (SNe Ib), 56 Type Ic SNe (SNe Ic), and 17 Type Ib/c SNe (SNe Ib/c). We compare the SE SN subtypes via measurements of the pseudo-equivalent widths (pEWs) and velocities of the He I $\lambda\lambda5876, 7065$ and O I $\lambda7774$ absorption lines. Consistent with previous work, we find that SNe Ic show higher pEWs and velocities in O I $\lambda7774$ compared to SNe IIb and Ib. The pEWs of the He I $\lambda\lambda5876, 7065$ lines are similar in SNe Ib and IIb after maximum light. The He I $\lambda\lambda5876, 7065$ velocities at maximum light are higher in SNe Ib compared to SNe IIb. We have identified an anticorrelation between the He I $\lambda7065$ pEW and O I $\lambda7774$ velocity among SNe IIb and Ib. This can be interpreted as a continuum in the amount of He present at the time of explosion. It has been suggested that SNe Ib and Ic have similar amounts of He, and that lower mixing could be responsible for hiding He in SNe Ic. However, our data contradict this mixing hypothesis. The observed difference in the expansion rate of the ejecta around maximum light of SNe Ic ($V_{\mathrm{m}}=\sqrt{2E_{\mathrm{k}}/M_{\mathrm{ej}}}\approx15,000$ km s$^{-1}$) and SNe Ib ($V_{\mathrm{m}}\approx9000$ km s$^{-1}$) would imply an average He mass difference of $\sim1.4$ $M_{\odot}$, if the other explosion parameters are assumed to be unchanged between the SE SN subtypes. We conclude that SNe Ic do not hide He but lose He due to envelope stripping., Comment: 18 pages, 6 figures, accepted for publication in Astronomy and Astrophysics

We apply a newly-developed On-the-Fly mosaicing technique on the NSF's Karl G. Jansky Very Large Array (VLA) at 3 GHz in order to carry out a sensitive search for an afterglow from the Advanced LIGO binary black hole merger event GW 151226. In three epochs between 1.5 and 6 months post-merger we observed a 100 sq. deg region, with more than 80% of the survey region having a RMS sensitivity of better than 150 uJy/beam, in the northern hemisphere having a merger containment probability of 10%. The data were processed in near-real-time, and analyzed to search for transients and variables. No transients were found but we have demonstrated the ability to conduct blind searches in a time-frequency phase space where the predicted afterglow signals are strongest. If the gravitational wave event is contained within our survey region, the upper limit on any late-time radio afterglow from the merger event at an assumed mean distance of 440 Mpc is about 1e29 erg/s/Hz. Approximately 1.5% of the radio sources in the field showed variability at a level of 30%, and can be attributed to normal activity from active galactic nuclei. The low rate of false positives in the radio sky suggests that wide-field imaging searches at a few Gigahertz can be an efficient and competitive search strategy. We discuss our search method in the context of the recent afterglow detection from GW 170817 and radio follow-up in future gravitational wave observing runs., Comment: 11 pages. 6 figures. 1 table. Accepted for publication in ApJ Letters

We present continued radio follow-up observations of PTF11qcj, a highly energetic broad-lined Type Ic supernova (SN), with a radio peak luminosity comparable to that of the $\gamma$-ray burst (GRB) associated SN 1998bw. The latest observations, carried out with the Karl G. Jansky Very Large Array (VLA), extend up to $\sim$5 years after the PTF11qcj optical discovery. The radio light curve shows a double-peak profile, possibly associated with density variations in the circumstellar medium (CSM), or with the presence of an off-axis GRB jet. Optical spectra of PTF11qcj taken during both peaks of the radio light curve do not show the broad H$\alpha$ features typically expected from H-rich circumstellar interaction. Modeling of the second radio peak within the CSM interaction scenario requires a flatter density profile and an enhanced progenitor mass-loss rate compared to those required to model the first peak. Although our radio data alone cannot rule out the alternative scenario of an off-axis GRB powering the second radio peak, the implied off-axis GRB parameters are unusual compared to typical values found for cosmological long GRBs. Deep X-ray observations carried out around the time of the second radio peak could have helped distinguish between the density variation and off-axis GRB scenarios. Future VLBA measurements of the PTF11qcj radio ejecta may unambiguously rule out the off-axis GRB jet scenario., Comment: 13 pages, 6 figures

There has been speculation of a class of relativistic explosions with an initial Lorentz factor smaller than that of classical Gamma-Ray Bursts (GRBs). These "dirty fireballs" would lack prompt GRB emission but could be pursued via their optical afterglow, appearing as transients that fade overnight. Here we report a search for such transients (transients that fade by 5-$\sigma$ in magnitude overnight) in four years of archival photometric data from the intermediate Palomar Transient Factory (iPTF). Our search criteria yielded 45 candidates. Of these, two were afterglows to GRBs that had been found in dedicated follow-up observations to triggers from the Fermi GRB Monitor (GBM). Another (iPTF14yb; Cenko et al. 2015) was a GRB afterglow discovered serendipitously. Two were spurious artifacts of reference image subtraction and one was an asteroid. The remaining 37 candidates have red stellar counterparts in external catalogs. The photometric and spectroscopic properties of the counterparts identify these transients as strong flares from M dwarfs of spectral type M3-M7 at distances of d ~ 0.15-2.1 kpc; two counterparts were already spectroscopically classified as late-type M stars. With iPTF14yb as the only confirmed relativistic outflow discovered independently of a high-energy trigger, we constrain the all-sky rate of transients that peak at m = 18 and fade by $\Delta$2 mag in $\Delta$3 hr to be 680 per year with a 68% confidence interval of 119-2236 per year. This implies that the rate of visible dirty fireballs is at most comparable to that of the known population of long-duration GRBs., Comment: accepted to ApJ Letters. 16 pages, 5 figures, 2 tables

GW170817 is the first gravitational wave detection of a binary neutron star merger. It was accompanied by radiation across the electromagnetic spectrum and localized to the galaxy NGC 4993 at a distance of 40 Mpc. It has been proposed that the observed gamma-ray, X-ray and radio emission is due to an ultra-relativistic jet launched during the merger, directed away from our line of sight. The presence of such a jet is predicted from models positing neutron star mergers as the central engines driving short-hard gamma-ray bursts (SGRBs). Here we show that the radio light curve of GW170817 has no direct signature of an off-axis jet afterglow. While we cannot rule out the existence of a jet pointing elsewhere, the observed gamma-rays could not have originated from such a jet. Instead, the radio data requires a mildly relativistic wide-angle outflow moving towards us. This outflow could be the high velocity tail of the neutron-rich material dynamically ejected during the merger or a cocoon of material that breaks out when a jet transfers its energy to the dynamical ejecta. The cocoon scenario can explain the radio light curve of GW170817 as well as the gamma-rays and X-rays (possibly also ultraviolet and optical emission), and hence is the model most consistent with the observational data. Cocoons may be a ubiquitous phenomenon produced in neutron star mergers, giving rise to a heretofore unidentified population of radio, ultraviolet, X-ray and gamma-ray transients in the local universe., Comment: 41 pages, 4 figures (main text), 4 figures (supplementary text), 1 table. This revised version incorporates comments from the referees

We performed a wide-area (2000 deg$^{2}$) g and I band experiment as part of a two month extension to the Intermediate Palomar Transient Factory. We discovered 36 extragalactic transients including iPTF17lf, a highly reddened local SN Ia, iPTF17bkj, a new member of the rare class of transitional Ibn/IIn supernovae, and iPTF17be, a candidate luminous blue variable outburst. We do not detect any luminous red novae and place an upper limit on their rate. We show that adding a slow-cadence I band component to upcoming surveys such as the Zwicky Transient Facility will improve the photometric selection of cool and dusty transients., Comment: 14 pages, 17 figures. Accepted to PASP

With the first direct detection of merging black holes in 2015, the era of gravitational wave (GW) astrophysics began. A complete picture of compact object mergers, however, requires the detection of an electromagnetic (EM) counterpart. We report ultraviolet (UV) and X-ray observations by Swift and the Nuclear Spectroscopic Telescope ARray (NuSTAR) of the EM counterpart of the binary neutron star merger GW170817. The bright, rapidly fading ultraviolet emission indicates a high mass ($\approx0.03$ solar masses) wind-driven outflow with moderate electron fraction ($Y_{e}\approx0.27$). Combined with the X-ray limits, we favor an observer viewing angle of $\approx 30^{\circ}$ away from the orbital rotation axis, which avoids both obscuration from the heaviest elements in the orbital plane and a direct view of any ultra-relativistic, highly collimated ejecta (a gamma-ray burst afterglow)., Comment: Science, in press; 56 pages, 12 figures

Merging neutron stars offer an exquisite laboratory for simultaneously studying strong-field gravity and matter in extreme environments. We establish the physical association of an electromagnetic counterpart EM170817 to gravitational waves (GW170817) detected from merging neutron stars. By synthesizing a panchromatic dataset, we demonstrate that merging neutron stars are a long-sought production site forging heavy elements by r-process nucleosynthesis. The weak gamma-rays seen in EM170817 are dissimilar to classical short gamma-ray bursts with ultra-relativistic jets. Instead, we suggest that breakout of a wide-angle, mildly-relativistic cocoon engulfing the jet elegantly explains the low-luminosity gamma-rays, the high-luminosity ultraviolet-optical-infrared and the delayed radio/X-ray emission. We posit that all merging neutron stars may lead to a wide-angle cocoon breakout; sometimes accompanied by a successful jet and sometimes a choked jet., Comment: Science, in press DOI 10.1126/science.aap9455, 83 pages, 3 tables, 16 figures