A statistical approach to identify superluminous supernovae and probe their diversity

We investigate the identification of hydrogen-poor superluminous supernovae (SLSNe I) using a photometric analysis, without including an arbitrary magnitude threshold. We assemble a homogeneous sample of previously classified SLSNe I from the literature, and fit their light curves using Gaussian processes. From the fits, we identify four photometric parameters that have a high statistical significance when correlated, and combine them in a parameter space that conveys information on their luminosity and color evolution. This parameter space presents a new definition for SLSNe I, which can be used to analyse existing and future transient datasets. We find that 90% of previously classified SLSNe I meet our new definition. We also examine the evidence for two subclasses of SLSNe I, combining their photometric evolution with spectroscopic information, namely the photospheric velocity and its gradient. A cluster analysis reveals the presence of two distinct groups. `Fast' SLSNe show fast light curves and color evolution, large velocities, and a large velocity gradient. `Slow' SLSNe show slow light curve and color evolution, small expansion velocities, and an almost non-existent velocity gradient. Finally, we discuss the impact of our analyses in the understanding of the powering engine of SLSNe, and their implementation as cosmological probes in current and future surveys.Comment: 16 pages, 9 figures, accepted by ApJ on 23/01/201

Euclid: Superluminous supernovae in the Deep Survey

Context. In the last decade, astronomers have found a new type of supernova called superluminous supernovae (SLSNe) due to their high peak luminosity and long light-curves. These hydrogen-free explosions (SLSNe-I) can be seen to z ~ 4 and therefore, offer the possibility of probing the distant Universe. Aims. We aim to investigate the possibility of detecting SLSNe-I using ESA’s Euclid satellite, scheduled for launch in 2020. In particular, we study the Euclid Deep Survey (EDS) which will provide a unique combination of area, depth and cadence over the mission. Methods. We estimated the redshift distribution of Euclid SLSNe-I using the latest information on their rates and spectral energy distribution, as well as known Euclid instrument and survey parameters, including the cadence and depth of the EDS. To estimate the uncertainties, we calculated their distribution with two different set-ups, namely optimistic and pessimistic, adopting different star formation densities and rates. We also applied a standardization method to the peak magnitudes to create a simulated Hubble diagram to explore possible cosmological constraints. Results. We show that Euclid should detect approximately 140 high-quality SLSNe-I to z ~ 3.5 over the first five years of the mission (with an additional 70 if we lower our photometric classification criteria). This sample could revolutionize the study of SLSNe-I at z > 1 and open up their use as probes of star-formation rates, galaxy populations, the interstellar and intergalactic medium. In addition, a sample of such SLSNe-I could improve constraints on a time-dependent dark energy equation-of-state, namely w(a), when combined with local SLSNe-I and the expected SN Ia sample from the Dark Energy Survey. Conclusions. We show that Euclid will observe hundreds of SLSNe-I for free. These luminous transients will be in the Euclid data-stream and we should prepare now to identify them as they offer a new probe of the high-redshift Universe for both astrophysics and cosmology.Acknowledgements. We thank the internal EC referees (P. Nugent and J. Brichmann) as well as the many comments from our EC colleagues and friends. C.I. thanks Chris Frohmaier and Szymon Prajs for useful discussions about supernova rates. C.I. and R.C.N. thank Mark Cropper for helpful information about the V IS instrument. C.I. thanks the organisers and participants of the Munich Institute for Astro- and Particle Physics (MIAPP) workshop “Superluminous supernovae in the next decade” for stimulating discussions and the provided online material. The Euclid Consortium acknowledges the European Space Agency and the support of a number of agencies and institutes that have supported the development of Euclid. A detailed complete list is available on the Euclid web site (http://www.euclid-ec.org). In particular the Agenzia Spaziale Italiana, the Centre National dEtudes Spatiales, the Deutsches Zentrum für Luft- and Raumfahrt, the Danish Space Research Institute, the Fundação para a Ciênca e a Tecnologia, the Ministerio de Economia y Competitividad, the National Aeronautics and Space Administration, The Netherlandse Onderzoekschool Voor Astronomie, the Norvegian Space Center, the Romanian Space Agency, the State Secretariat for Education, Research and Innovation (SERI) at the Swiss Space Office (SSO), the United Kingdom Space Agency, and the University of Helsinki. R.C.N. acknowledges partial support from the UK Space Agency. D.S. acknowledges the Faculty of Technology of the University of Portsmouth for support during his PhD studies. C.I. and S.J.S. acknowledge funding from the European Research Council under the European Union’s Seventh Framework Programme (FP7/2007-2013)/ERC Grant agreement No. [291222]. C.I. and M.S. acknowledge support from EU/FP7-ERC grant No. [615929]. E.C. acknowledge financial contribution from the agreement ASI/INAF/I/023/12/0. The work by KJ and others at MPIA on NISP was supported by the Deutsches Zentrum für Luft- und Raumfahrt e.V. (DLR) under grant 50QE1202. M.B. and S.C. acknowledge financial contribution from the agreement ASI/INAF I/023/12/1. R.T. acknowledges funding from the Spanish Ministerio de Economía y Competitividad under the grant ESP2015-69020-C2- 2-R. I.T. acknowledges support from Fundação para a Ciência e a Tecnologia (FCT) through the research grant UID/FIS/04434/2013 and IF/01518/2014. J.R. was supported by JPL, which is run under a contract for NASA by Caltech and by NASA ROSES grant 12-EUCLID12-0004

Using late-time optical and near-infrared spectra to constrain Type Ia supernova explosion properties

The late-time spectra of Type Ia supernovae (SNe Ia) are powerful probes of the underlying physics of their explosions. We investigate the late-time optical and near-infrared spectra of seven SNe Ia obtained at the VLT with XShooter at $>$200 d after explosion. At these epochs, the inner Fe-rich ejecta can be studied. We use a line-fitting analysis to determine the relative line fluxes, velocity shifts, and line widths of prominent features contributing to the spectra ([Fe II], [Ni II], and [Co III]). By focussing on [Fe II] and [Ni II] emission lines in the ~7000-7500 \AA\ region of the spectrum, we find that the ratio of stable [Ni II] to mainly radioactively-produced [Fe II] for most SNe Ia in the sample is consistent with Chandrasekhar-mass delayed-detonation explosion models, as well as sub-Chandrasekhar mass explosions that have metallicity values above solar. The mean measured Ni/Fe abundance of our sample is consistent with the solar value. The more highly ionised [Co III] emission lines are found to be more centrally located in the ejecta and have broader lines than the [Fe II] and [Ni II] features. Our analysis also strengthens previous results that SNe Ia with higher Si II velocities at maximum light preferentially display blueshifted [Fe II] 7155 \AA\ lines at late times. Our combined results lead us to speculate that the majority of normal SN Ia explosions produce ejecta distributions that deviate significantly from spherical symmetry.Comment: 17 pages, 12 figure, accepted for publication in MNRA

Moderately Luminous type II Supernovae

Core-collapse Supernovae (CC-SNe) descend from progenitors more massive than about 8 Msun. Because of the young age of the progenitors, the ejecta may eventually interact with the circumstellar medium (CSM) via highly energetic processes detectable in the radio, X-ray, ultraviolet (UV) and, sometimes, in the optical domains. In this paper we present ultraviolet, optical and near infrared observations of five type II SNe, namely SNe 2009dd, 2007pk, 2010aj, 1995ad, and 1996W. Together with few other SNe they form a group of moderately luminous type II events. We collected photometry and spectroscopy with several telescopes in order to construct well-sampled light curves and spectral evolutions from the photospheric to the nebular phases. Both photometry and spectroscopy indicate a degree of heterogeneity in this sample. The light curves have luminous peak magnitudes ($-16.95<M_{B}<-18.70$). The ejected masses of ^56\ni for three SNe span a wide range of values ($2.8\times10^{-2}$Msun$<$M(\ni)$<1.4\times10^{-1}$Msun), while for a fourth (SN2010aj) we could determine a stringent upper limit ($7\times10^{-3}$Msun). Clues of interaction, such as the presence of high velocity (HV) features of the Balmer lines, are visible in the photospheric spectra of SNe 2009dd and 1996W. For SN2007pk we observe a spectral transition from a type IIn to a standard type II SN. Modelling the observations of SNe 2009dd, 2010aj and 1995ad with radiation hydrodynamics codes, we infer kinetic plus thermal energies of about 0.2-0.5 foe, initial radii of 2-5$\times10^{13}$ cm and ejected masses of $\sim$5.0-9.5 Msun. These values suggest moderate-mass, super-asymptotic giant branch (SAGB) or red super-giants (RSG) stars as SN precursors, in analogy with other luminous type IIP SNe 2007od and 2009bw.Comment: 28 pages, 27 fig, accepted by A&A, 3 pages of online material, abstract abridged. revised significantly with respect to the previous versio

The Type IIP SN 2007od in UGC 12846: from a bright maximum to dust formation in the nebular phase

Ultraviolet (UV), optical and near infrared (NIR) observations of the type IIP supernova (SN) 2007od are presented, covering from the maximum light to the late phase, allowing to investigate in detail different physical phenomena in the expanding ejecta. These data turn this object into one of the most peculiar IIP ever studied. The early light curve of SN 2007od is similar to that of a bright IIPs with a short plateau, a bright peak (MV = -18 mag), but a very faint optical light curve at late time. However, with the inclusion of mid infrared (MIR) observations during the radioactive decay we have estimate a M(56Ni) ~ 2\times10^-2 M\odot. Modeling the bolometric light curve, ejecta expansion velocities and black-body temperature, we estimate a total ejected mass was 5 - 7.5 M\odot with a kinetic energy of at least 0.5 \times 10^51 erg. The early spectra reveal a boxy H{\alpha} profile and high velocities features of the Balmer series that suggest interaction between the ejecta and a close circum-stellar matter (CSM). SN 2007od may be, therefore, an intermediate case between a Type IIn SN and a typical Type IIP SN. Also late spectra show a clear evidence of CSM and the presence of dust formed inside the ejecta. The episodes of mass loss short before explosion, the bright plateau, along with the relatively small amount of 56Ni and the faint [O I] observed in the nebular spectra are consistent with a super-asympthotic giant branch (super-AGB) progenitor (M~9.7 - 11 M\odot).Comment: V2, some test added and three figures changed from the first version. 21 pages, 18 figures, accepted for publication in MNRAS on May 24, 201

The host galaxy and late-time evolution of the Super-Luminous Supernova PTF12dam

Super-luminous supernovae of type Ic have a tendency to occur in faint host galaxies which are likely to have low mass and low metallicity. PTF12dam is one of the closest and best studied super-luminous explosions that has a broad and slowly fading lightcurve similar to SN 2007bi. Here we present new photometry and spectroscopy for PTF12dam from 200-500 days (rest-frame) after peak and a detailed analysis of the host galaxy (SDSS J142446.21+461348.6 at z = 0.107). Using deep templates and image subtraction we show that the full lightcurve can be fit with a magnetar model if escape of high-energy gamma rays is taken into account. The full bolometric lightcurve from -53 to +399 days (with respect to peak) cannot be fit satisfactorily with the pair-instability models. An alternative model of interaction with a dense CSM produces a good fit to the data although this requires a very large mass (~ 13 M_sun) of hydrogen free CSM. The host galaxy is a compact dwarf (physical size ~ 1.9 kpc) and with M_g = -19.33 +/- 0.10, it is the brightest nearby SLSN Ic host discovered so far. The host is a low mass system (2.8 x 10^8 M_sun) with a star-formation rate (5.0 M_sun/year), which implies a very high specific star-formation rate (17.9 Gyr^-1). The remarkably strong nebular lines provide detections of the [O III] \lambda 4363 and [O II] \lambda\lambda 7320,7330 auroral lines and an accurate oxygen abundance of 12 + log(O/H) = 8.05 +/- 0.09. We show here that they are at the extreme end of the metallicity distribution of dwarf galaxies and propose that low metallicity is a requirement to produce these rare and peculiar supernovae.Comment: 20 pages, 12 figures, 8 tables, accepted for publication to MNRA

A comparative study of Type II-P and II-L supernova rise times as exemplified by the case of LSQ13cuw

We report on our findings based on the analysis of observations of the Type II-L supernova LSQ13cuw within the framework of currently accepted physical predictions of core-collapse supernova explosions. LSQ13cuw was discovered within a day of explosion, hitherto unprecedented for Type II-L supernovae. This motivated a comparative study of Type II-P and II-L supernovae with relatively well-constrained explosion epochs and rise times to maximum (optical) light. From our sample of twenty such events, we find evidence of a positive correlation between the duration of the rise and the peak brightness. On average, SNe II-L tend to have brighter peak magnitudes and longer rise times than SNe II-P. However, this difference is clearest only at the extreme ends of the rise time versus peak brightness relation. Using two different analytical models, we performed a parameter study to investigate the physical parameters that control the rise time behaviour. In general, the models qualitatively reproduce aspects of the observed trends. We find that the brightness of the optical peak increases for larger progenitor radii and explosion energies, and decreases for larger masses. The dependence of the rise time on mass and explosion energy is smaller than the dependence on the progenitor radius. We find no evidence that the progenitors of SNe II-L have significantly smaller radii than those of SNe II-P.Comment: 19 pages, 10 figures, accepted by A&

Observations of the GRB afterglow ATLAS17aeu and its possible association with GW170104

We report the discovery and multi-wavelength data analysis of the peculiar optical transient, ATLAS17aeu. This transient was identified in the skymap of the LIGO gravitational wave event GW170104 by our ATLAS and Pan-STARRS coverage. ATLAS17aeu was discovered 23.1hrs after GW170104 and rapidly faded over the next 3 nights, with a spectrum revealing a blue featureless continuum. The transient was also detected as a fading x-ray source by Swift and in the radio at 6 and 15 GHz. A gamma ray burst GRB170105A was detected by 3 satellites 19.04hrs after GW170104 and 4.10hrs before our first optical detection. We analyse the multi-wavelength fluxes in the context of the known GRB population and discuss the observed sky rates of GRBs and their afterglows. We find it statistically likely that ATLAS17aeu is an afterglow associated with GRB170105A, with a chance coincidence ruled out at the 99\% confidence or 2.6$\sigma$. A long, soft GRB within a redshift range of $1 \lesssim z \lesssim 2.9$ would be consistent with all the observed multi-wavelength data. The Poisson probability of a chance occurrence of GW170104 and ATLAS17aeu is $p=0.04$. This is the probability of a chance coincidence in 2D sky location and in time. These observations indicate that ATLAS17aeu is plausibly a normal GRB afterglow at significantly higher redshift than the distance constraint for GW170104 and therefore a chance coincidence. However if a redshift of the faint host were to place it within the GW170104 distance range, then physical association with GW170104 should be considered.Comment: 16 pages, 6 figures, accepted to Ap