Free Access
Issue
A&A
Volume 626, June 2019
Article Number A117
Number of page(s) 18
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/201935171
Published online 24 June 2019

© ESO 2019

1. Introduction

High-energy neutrinos are produced in hadronic interactions in the most violent sources in the Universe and hence provide a smoking-gun signature for the origin of cosmic rays. Furthermore, neutrinos are capable of escaping even the densest environments such as stellar cores and black hole accretion discs. Along with gravitational waves (Abbott et al. 2016, 2017), they carry information of these environments that neither photons or charged particles can provide. After the first detection of a diffuse flux of high-energy astrophysical neutrinos (Aartsen et al. 2013a), the most pressing question in the field of neutrino astronomy is the astrophysical origin of the particles. No significant cluster in space or time has been found in un-triggered searches of neutrino data yet (Aartsen et al. 2017a) and the isotropic distribution of the neutrinos points to an extragalactic origin. Multi-wavelength observations are the key to identify neutrino sources through the detection of an electromagnetic counterpart. Recently, the high-energy 290 TeV neutrino event IceCube-170922A has been associated with a flare in the blazar TXS 0506+056 at 3σ significance (Aartsen et al. 2018a,b). Photohadronic processes in the cores of active galaxies have previously been proposed as the sources of these high-energy neutrinos (Padovani & Resconi 2014; Petropoulou et al. 2015) and the IceCube-170922A association is an intriguing development in the field. The only other identified extragalactic source of neutrinos is the core-collapse supernova (CCSN) 1987A. Undisputed detection of 20 neutrinos were made by two experiments in the energy range of 7.5–40 MeV (Hirata et al. 1987; Bionta et al. 1987), corresponding to the thermal regime in the collapsing core of the massive star that produced SN 1987A.

Gamma-ray bursts (GRB) have also been suggested as source candidates for the highest-energy cosmic rays and high-energy neutrinos (Waxman 1995; Waxman & Bahcall 1997). Aartsen et al. (2016) suggest that the bulk of the known population of bright GRBs can be excluded as the main contributor to IceCube’s diffuse neutrino flux. However, a large population of low-luminosity GRBs might contribute significantly to the observed flux (Murase et al. 2006; Murase & Ioka 2013). While highly-relativistic jets can explain gamma-ray bright bursts, “choked” jets within the star may explain relativistic SNe and low-luminosity GRBs, giving a unified picture of GRBs and GRB-SNe (Meszaros & Waxman 2001; Senno et al. 2016). In choked-jet SNe, the jet fails to penetrate the surface of the progenitor star. Neutrinos and gamma rays can be produced by hadronic processes in the jet. While neutrinos can leave the dense environment, gamma rays are absorbed. This scenario could be physically probed by the detection of high-energy neutrinos in coincidence with SNe that could harbour relativistic jets (Kowalski & Mohr 2007; Esmaili & Murase 2018). Such neutrino emission is expected in a relatively short time window (∼100 s) after core collapse. The signature of a choked-jet SN could therefore be a high-energy neutrino signal followed by a young CCSN. Explosions of CCSNe originate from stars more massive than M ≈ 8 M and are classified based on their spectroscopic features (Filippenko 1997). Type II SNe are H-rich and Type Ib/c SNe H-poor. Furthermore, Type Ib SNe show He in their spectra whereas Type Ic SNe lack these features. Type Ic SNe are thought to originate from massive stars that have lost their outer envelopes (Wolf-Rayet stars, e.g. Crowther 2007) or less massive stars in binary systems (e.g. Podsiadlowski et al. 1992). High-energy neutrino production in Type Ic SNe has been postulated to be possible, as a result of diffusive shock acceleration of protons in relativistic jets powered by a central engine (Razzaque et al. 2005; Ando & Beacom 2005; Tamborra & Ando 2016; Senno et al. 2016). Those protons can produce pions in interactions with ambient photon or matter fields. Charged pions produce neutrinos in their decay chain. The mechanism follows the GRB-fireball phenomenology (Piran 1999), but unlike conventional long GRBs with Γ factors of 102 − 103, the culprit SNe are assumed to harbour fewer relativistic jets. Such jets can be either choked or emerge, thereby producing low-luminosity or long-duration GRBs (potentially undetected by current facilities). Furthermore, it is expected that Type IIn SNe may also produce high-energy neutrinos arising from the shock interaction of the expanding SN ejecta and circumstellar medium. A Type IIn SN PTF12csy was discovered within the confidence region of two high-energy neutrinos detected by IceCube on 2012 March 30. However, the Pan-STARRS1 data showed the SN to be at least 158 days old and therefore likely unrelated (Aartsen et al. 2015). Contrary to CCSNe, Type Ia SNe are thermonuclear explosions of white dwarf stars and are not expected to produce high-energy neutrinos.

In the following sections, we briefly introduce the IceCube realtime neutrino search in Sect. 2 followed by a detailed discussion of the neutrino event IceCube-160427A in Sect. 3. Section 4 presents the observations and interpretation of a young SN PS16cgx discovered in the confidence region of the neutrino event. Summary of Pan-STARRS1 observations of other IceCube alerts is provided in Sect. 5. The probability to associate neutrino emitters with Pan-STARRS1 is discussed in Sect. 6 and the results of the paper are summarised in Sect. 7. In all cases we adopt a calibration to standard cosmology with H0 = 70 km s−1 Mpc−1, Ωm = 0.3, and ΩΛ = 0.7.

2. The IceCube realtime neutrino search

The IceCube neutrino Observatory is a 1 km3 water Cherenkov detector located in the clear glacial ice at the geographic South Pole. 5160 digital optical modules (DOMs) detect light emission in the ice at a depth ranging from 1450 to 2450 m (Abbasi et al. 2009; Aartsen et al. 2013b, 2017c). Charged-current interactions of muon neutrinos in the ice produce muons, which are detected using emitted Cherenkov light, and which – at TeV energies – travel long distances of several kilometres in the Antarctic ice (Chirkin & Rhode 2004). The resulting long lever arm of those “track events” gives a good reconstruction performance with median angular resolution of ΔΨ ≲ 1°. Charged-current interactions of electron or tau neutrinos, as well as neutral-current interactions of any neutrino type, produce “shower events” with almost spherically symmetric light emission resulting in an angular resolution of only about 10°.

The IceCube realtime programme selects neutrinos with a high probability of being of astrophysical origin, reconstructs their direction and notifies ground and space-based observatories with a median latency of 33 s to perform multi-wavelengths follow-up observations aiming for the detection of an electro-magnetic counterpart (see Aartsen et al. 2017b, for a detailed description).

Starting in late 2008, IceCube implemented an optical follow-up programme on multiplets. In order to suppress the background from atmospheric neutrinos, a multiplet of at least two neutrinos within 100 s and angular separation of 3 . ° $ \overset{\circ }{.} $5 or less was required to trigger an alert (Abbasi et al. 2012; Aartsen et al. 2017b). A number of alerts were sent to the Palomar Transient Factory (PTF) and Swift (Evans 2015) and which has not resulted in a compelling coincidence with an electromagnetic source (Aartsen et al. 2015). In 2016, IceCube began releasing single neutrinos in two neutrino candidate streams (Aartsen et al. 2017b) called High-Energy Starting Events (HESE) and Extremely High-energy Events (EHE), which are the focus of this paper.

The HESE sample consists of events which have their interaction vertices inside the detector volume and are selected by using the outer layers of the detector as a veto to suppress the background of atmospheric muons, which are produced in cosmic-ray air showers and enter the detector from above (Aartsen et al. 2014). The veto is only efficient at high energies and is therefore combined with a requirement on the charge deposited in the detector of Q >  6000 photoelectrons, which roughly corresponds to an energy deposited in the detector of 60 TeV. On average, four HESE tracks are detected per year, of which one is expected to be of astrophysical origin. Since April 2016, these have been made public in realtime using the Astrophysical Multimessenger Observatory Network (AMON; Smith et al. 2013; Aartsen et al. 2017b) and via the Gamma-Ray Coordinates Network1 (GCN).

EHE events, which do not have to start inside the detector but can be through-going, are selected by a combination of energy and zenith angle (see Aartsen et al. 2017b, for a detailed description). These neutrinos typically have energies > 200 TeV. The public EHE stream was activated in July 2016 and is expected to produce four alerts per year with an average signal purity of 50%. In this paper we introduce the neutrino detection electromagnetic counterpart search with Pan-STARRS1 with a focus on the first publicly released event IceCube-160427A.

3. IceCube-160427A

The first HESE event after the start of the programme was detected on 2016 April 27 at 05:52:32 UT (MJD 57 505.24481) and released to the public via a GCN notice2. A summary was published as a GCN circular a few days later (Blaufuss 2016).

Several instruments followed-up the IceCube event searching for a counterpart in various wavelengths. In gamma rays, observations were performed by the First G-APD Cherenkov Telescope (FACT; Biland & Dorner 2016), the Interplanetary Network (IPN3; Svinkin et al. 2016), the Very Energetic Radiation Imaging Telescope Array System (VERITAS; Mukherjee 2016), the High-Altitude Water Cherenkov Gamma-Ray Observatory (HAWC; Taboada 2016), the Fermi Large Area Telescope (Fermi-LAT; Vianello et al. 2016) and the Fermi Gamma-Ray Burst Monitor (Fermi-GBM; Blackburn et al. 2016), and in optical wavelength by the Intermediate Palomar Transient Factory (iPTF; Singer et al. 2016a), the Panoramic Survey Telescope and Rapid Response System 1 (Pan-STARRS1; Smartt et al. 2016a) and the Mobile Astronomical System of Telescope-Robots (MASTER; Lipunov et al. 2016a). Except a young SN discovered by Pan-STARRS1, no compelling counterpart detection was reported.

3.1. Direction

An angular reconstruction is applied to the HESE events at the South Pole, which is used to send a first notification. More time consuming reconstructions are applied, once the event data arrive in the north and provide a more accurate estimate on the neutrino direction within a few hours. A full likelihood scan is applied on a narrow grid with about 0 . ° $ \overset{\circ }{.} $06 distance between the grid points. The resulting map of the likelihood landscape allows us to find the global minimum and the error contours at a given confidence level. The calculation of the updated direction was not fully automated yet when the first public HESE track was detected, which led to an additional delay. Table 1 shows the initially published direction 81 s after the detection of the neutrino compared to the updated position 17.5 h later. The position moved by 2 . ° $ \overset{\circ }{.} $6.

Table 1.

Coordinate revisions of IceCube-160427A.

The likelihood scan was repeated after cleaning out late photomultiplier hits from afterpulsing, an effect not currently described by the track likelihood hypothesis. This improves the likelihood description and leads to a better fit. The final best fit position and the 50 and 90% confidence regions are shown in Fig. 1. The positions of the contour lines were determined using re-simulations of similar events while also varying the allowed models of the optical properties of the deep glacial ice in order to obtain a conservative range. The minimum yields our best estimate of the event direction: α J 2000 = 240 . ° 3 $ \alpha_{\mathrm{J2000}} = 240{{\overset{\circ}{.}}}3 $ and δ J 2000 = 9 . ° 7 $ \delta_{\mathrm{J2000}} = 9{{\overset{\circ}{.}}}7 $, with a 90% error ellipse contour with semi major (minor) axis of 0 . ° $ \overset{\circ }{.} $54 (0 . ° $ \overset{\circ }{.} $45).

thumbnail Fig. 1.

Black (red) line: 50% (90%) confidence region. The confidence region likelihood ratios were determined using Monte Carlo re-simulations of the event, taking into account uncertainties in the optical properties of the deep glacial ice. The best fit position is shown as a × mark and the position of PS16cgx is indicated with a star symbol. Furthermore, the other faint transients discovered within the final 90% probability containment by Pan-STARRS1 are also shown and contextually classified based on their host or precursor information as labeled in the legend.

3.2. Energy

The total charge deposited in the detector amounts to 1.9 × 104 photoelectrons, corresponding to a deposited energy of 134 TeV. Using Monte Carlo simulations, we estimate a lower bound on the true energy of the primary neutrino. In the following we assume a spectral index of −2.58 for the astrophysical flux (Kopper et al. 2015), which is the best fit value obtained from the 4 yr HESE sample. We then simulate events with the same starting position in the detector (within ±30 m horizontally and ±15 m vertically), direction (within 2°) and deposited energy (within ±20%). By comparison of IceCube-160427A to the simulated events we estimate the true neutrino energy to be larger than 130 TeV (230 TeV) at 90% (50%) confidence (which, as expected, is similar to the deposited energy of the event). The result is only weakly dependent on the chosen spectral index (see Aartsen et al. 2018b).

3.3. Signal purity

On average 25% of the events in the HESE realtime stream are expected to be astrophysical neutrinos. The HESE events are characterised by the signal-trackness parameter, which is a number between 0 and 1, characterising how likely the event is both astrophysical as well as track-like, opposed to shower-like (see Aartsen et al. 2017b, for more details). For IceCube-160427A, the signal-trackness parameter is 0.92. We estimate the fraction of astrophysical signal events in the alert sample for events with a similar or larger deposited charge and which enter the detector from a similar zenith direction. We assume a power-law signal spectrum. We find a fraction of 46−58% depending on the assumed spectral index ranging from γ = −2.19, found in the analysis of the diffuse muon neutrino flux (Haack et al. 2017), to γ = −2.92, obtained in an analysis of high-energy starting-tracks (Kopper & IceCube Collaboration 2017). The fraction is similar to that of IceCube-170922A, which had a higher reconstructed energy but passed the EHE instead of the HESE stream selection criteria. The HESE event selection uses the outer detector layers as veto to select track events starting inside the instrumented detector volume.

3.4. Offline neutrino search

We investigate whether the HESE event was accompanied by an excess of low-energy neutrino events. We applied a maximum likelihood analysis following the procedure in Meagher & IceCube Collaboration (2017) to search for neutrino emission from a point source in the direction of PS16cgx on various time scales of 10, 100, 1000, 10 000 and 100 000 s centred on the time of the HESE event. The chosen time windows cover the shock propagation and shock breakout times scales estimated by Kistler et al. (2013). The HESE event was removed from the likelihood analysis to avoid a bias in the p-value calculation. No significant excess was found in any of the tested time windows leading to a 90% energy flux upper limit of 3.7 × 10−5 TeV cm−2 over the 100 000 s time window in an approximate energy range of 1 TeV to 1 PeV assuming a neutrino spectrum following an E−2 dependence. The non-detection of additional neutrinos is not surprising, since we expect sources of two detectable neutrinos at a median redshift of 0.06 (see Fig. 8 of Aartsen et al. 2017d) assuming the sources follow redshift evolution of the star-formation rate.

4. Pan-STARRS1 observations of IceCube EHE and HESE alerts

The 1.8 m Pan-STARRS1 telescope (Chambers et al. 2016) is located at Haleakala Observatory, Hawaii, USA. It is equipped with a 1.4 gigapixel mosaic CCD camera, called Gigapixel Camera 1 (GPC1). The field of view of GPC1 is ∼7 deg2 and has a pixel size of 0.25″ pix−1. The Pan-STARRS1 data are processed automatically by the real-time Pan-STARRS1 Image Processing Pipeline (IPP; Magnier et al. 2016a). The photometric calibration of the Pan-STARRS1 images is described in Magnier et al. (2016b). Each image has a reference sky subtracted from it to produce a calibrated difference image. After this reduction, calibration and image subtraction in Hawaii, the catalogues of all detections are transferred to Queen’s University Belfast where they are ingested into a MySQL database, assimilated into objects, and undergo a series of filters to produce a clean stream of real astrophysical transients (Wright et al. 2015). They are cross-matched with a large custom-built database of star and galaxy catalogues and annotated with associations to known sources. The processing is described in more detail in Smartt et al. (2016b) for searches for optical signatures of gravitational wave sources. Similar techniques were employed here. The diameter of the field of view of Pan-STARRS1 is 3° (Chambers et al. 2016) and hence one pointing typically encloses the full 90% confidence area estimated by IceCube, which tend to have radii of around ∼0 . ° $ \overset{\circ }{.} $5−1°. We have typically employed a strategy of multiple dithered exposures in iP1 or rP1 filters as were done in the Pan-STARRS1 Science Consortium Medium Deep Fields (for details see Rest et al. 2014; McCrum et al. 2015; Lunnan et al. 2018; Scolnic et al. 2018). These dithered exposures are usually 8 × 240 s long. They are stacked together to make one 1920 s exposure although we also process the individual 8 exposures. We carried out a concerted campaign to observe the error boxes of EHE and HESE events in 2016 and 2017 with Pan-STARRS1, and observed five in total (IceCube-160427A, IceCube-160731A, IceCube-161103A, IceCube-170321A, and IceCube-171106A). We focus here on details of the first event released (IceCube-160427A) since that was the only event in which we found a young, rising extragalactic transient with a possible temporal coincidence with the neutrino detection time. The results from the rest of the events are summarised in Sect. 5.

The Pan-STARRS1 telescope began observing the field of IceCube-160427A in the iP1-band on MJD 57 508.5 (2016 April 30.5 UT), 3.3 d after the IceCube detection, and the observations were repeated in the iP1-band on 4 subsequent nights. Individual dithered images were stacked together and the 3π Steradian Survey reference sky (see Chambers et al. 2016) was subtracted from these nightly stacked images. The 3π reference sky is constructed from the median combination of multiple individual exposures that are of length 43 (in filter gP1), 40 (rP1), 45 (iP1), 40 (zP1) and 30 (yP1) s. Typically there are around 12 individual exposures in each filter at each point in the sky that are combined to form the reference sky, which were taken over the period between 2009−2014. However, the number of exposures per point on the sky is a distribution covering the range between 5 and 25 (see Fig. 16 of Chambers et al. 2016). For follow-up of IceCube neutrino alerts, the iP1 filter was typically used and, in most cases, the number of exposures in the iP1 reference sky is around 16. The typical depth of the reference sky is miP1 ≈ 23.1 mag (see Table 11 and Fig. 17 of Chambers et al. 2016).

After subtraction, the stacked images typically had a 5σ point-source image depth of miP1 ≈ 22.5 mag. The initial results, based on the first revision, of the IceCube-160427A coordinates were reported by Smartt et al. (2016a). In this paper, we revise these results with the final IceCube-160427A position. As an important null result, we report that no persistent bright (mi ≲ 20 mag) transients were found in the IceCube-160427A field over a period of +3.3–45 days. While we did not cover the first 3 days, we rule out any relatively bright SN-like transient that exploded at a time coincident with IceCube-160427A. Other optical surveys by our team and others have found a similar null result for IceCube events, in that there are no bright candidates with SN-like light curves in the field (e.g. Lipunov et al. 2016a,b; Singer et al. 2016a,b).

For IceCube-160427A, we found 19 faint optical transients within the 90% probability containment of the IceCube-160427A position: 4 variable known quasi-stellar objects (QSOs), that is, extremely luminous active galactic nuclei (AGNs), based on the Sloan Digital Sky Survey (SDSS) Data Release 12 (DR12; Alam et al. 2015) archival spectra; 6 likely QSOs based on the host colour; 6 SN candidates; and 3 faint sources that are likely either faint Galactic variable stars or QSO variability, see Table 2. The SN candidates are transients that are either offset from or coincident with an apparent host galaxy which is not a known QSO or AGN. Five of the SN candidates show flat or declining light curves, consistent with SNe at peak or post maximum. However, one of the candidates, PS16cgx, shows a rising light curve of 0.4 mag over 2 days suggesting that it is a young SN and therefore possibly having an explosion epoch that is temporally consistent with IceCube-160427A. PS16cgx was located at αJ2000 = 240 . ° $ \overset{\circ }{.} $32755 and δJ2000 = +9 . ° $ \overset{\circ }{.} $86478 (16h01m18 . s $ \overset{\text{ s}}{.} $61 +09 ° 51 53 . 1 $ ^\circ 51\prime 53\mathop .\limits^{\prime \prime } 1 $), 10.0′ from the best-fit IceCube-160427A direction and consistent with the neutrino position within the uncertainties at 50% confidence.

Table 2.

19 transients discovered within the final 90% probability containment of the IceCube-160427A neutrino detection.

4.1. Follow-up of PS16cgx

The field of PS16cgx is shown in Fig. 2. We continued to observe the field with Pan-STARRS1 until MJD = 57 550, or 45 days after discovery of the neutrino. These follow-up images were primarily carried out in the iP1-band and the nightly stacked images of the field typically reached 22 ≲ miP1 ≲ 22.5 (AB mag, 5σ). One epoch of rP1 and two epochs of zP1-band imaging were also obtained. Furthermore, two epochs of i-band late-time imaging were obtained with the Canada France Hawaii Telescope (CFHT), using the MegaCam instrument (Boulade et al. 2003). The automatically calibrated level 2 data products were downloaded for the analysis. The CFHT observations were obtained at MJD 57 575.3 (2016 July 6.3 UT) suggesting mi = 23.649 ± 0.153 mag, and MJD 57 637.2 (2016 September 6.2 UT), with the latter epoch used as a reference template image for the first CFHT epoch. However, it is possible that the latter epoch contains some contamination from the SN. Assuming that PS16cgx is a Type Ia SN, we estimate the magnitude of the SN to be mi ≈ 24.5 mag during this epoch. Adding this contribution mathematically to the template subtracted magnitude at the MJD 57 575.3 epoch results in a magnitude of mi = 23.2 mag. The resulting photometry of PS16cgx is reported in Table 3.

thumbnail Fig. 2.

Pan-STARRS1 riz-band false-colour 1′×1′ subsection of the field of PS16cgx at MJD 57 525.5. North is up, east is left.

Table 3.

Photometry of PS16cgx.

Spectra of PS16cgx (see Table 4) were obtained with the following instruments: the 2.2 m University of Hawaii (UH) telescope with SuperNovae Integral Field Spectrograph (SNIFS; Lantz et al. 2004), the 8.1 m Gemini-North Telescope with Gemini Multi-Object Spectrograph (GMOS; Hook et al. 2004), and the 8.2 m Very Large Telescope (VLT) with FOcal Reducer and low dispersion Spectrograph 2 (FORS2; Appenzeller et al. 1998). The SNIFS spectrum was reduced using a dedicated calibration pipeline for the instrument (see also Smartt et al. 2016b, Sect. 2 for details of Pan-STARRS follow-up with SNIFS). The GMOS spectrum was reduced using the gemini IRAF4 package. The FORS2 spectra were processed using standard IRAF tasks and reduction steps, and combined for higher signal-to-noise ratio. The redshift of PS16cgx is constrained by the narrow host galaxy emission lines of Hα and [S II] 6716.4,6730.8 Å to z = 0.2895 ± 0.0001. This yields a luminosity distance of Dl = 1490 Mpc (μ = 40.87). The Galactic extinction towards PS16cgx is Ai = 0.073 mag based on the dust map calibration of Schlafly & Finkbeiner (2011), which is also the calibration system we adopt for all comparisons. Unfortunately, the Na I D feature at the rest wavelength of PS16cgx is redshifted in our observer frame to overlap with the heavily contaminating telluric A band around 7615 Å, and cannot be used to estimate the host galaxy reddening. For the analysis no host galaxy extinction is adopted. The K-corrections for the PS16cgx photometry are calculated from the GMOS spectrum with the SNAKE code (Inserra et al. 2018) yielding Krg = −0.32, Kir = −0.28, and Kzi = −0.29 mag.

Table 4.

Spectroscopy of PS16cgx.

4.2. Nature of PS16cgx

Due to the faint observed brightness of PS16cgx around 21–22 mag in i-band, it is challenging to obtain high signal-to-noise spectra of such an event even with an 8 m class telescope. Due to this, the nature of PS16cgx is ambiguous between a Type Ia and a Type Ic SN based solely on spectroscopic classification (Fig. 3). A Type Ic origin of PS16cgx would be extremely interesting, since both GRB associated Type Ic SNe and also some normal Type Ic SNe (e.g. Paragi et al. 2010) are associated with relativistic jets immediately after the core collapse. These jets are expected to produce a short (∼10 s) burst of neutrinos; a detection of such neutrinos could produce valuable insights into the physics of these jets (e.g. Ando & Beacom 2005) and identify the class of choked-jet SNe as neutrino emitters.

thumbnail Fig. 3.

Panel a: spectra of PS16cgx compared to a selection of Type Ic SNe. The comparison events include the GRB associated Type Ic-BL event SN 1998bw (Patat et al. 2001); non-GRB associated Type Ic-BL event SN 2002ap (Modjaz et al. 2014); and normal Type Ic events, SN 2007gr (Hunter et al. 2009), and SN 2004aw (Taubenberger et al. 2006). The strong host galaxy emission lines of PS16cgx are labelled. Panel b: spectra of PS16cgx compared to a selection of Type Ia SNe. The comparison events include SN 2005eq, SN 2005na (Blondin et al. 2012), SN 2006py (Folatelli et al. 2013), and SN 2011fe (Pereira et al. 2013). In both panels the spectra have been corrected for line-of-sight reddening and to the rest-frame wavelength. The wavelengths of the main telluric bands in the spectra of PS16cgx are indicated with a ⊕ symbol. Comparison data were obtained via WISeREP (Yaron & Gal-Yam 2012).

The lack of very broad features in the spectra of PS16cgx suggest that it is not a broad-line Type Ic (Type Ic-BL) SN such as SN1998bw (Patat et al. 2001). While the spectra of PS16cgx do not exactly resemble those of the normal Type Ic SN 2007gr (Hunter et al. 2009) either, there is some similarity for example with the Type Ic SN 2004aw (Taubenberger et al. 2006). The main difference in comparison of PS16cgx to normal Type Ia SNe is the shallow absorption feature that could be associated with the Si II line. Nonetheless, the early spectrum of PS16cgx shares significant similarity to that of Type Ia SNe, such as SN 2006py (Folatelli et al. 2013).

Based on a third order polynomial fit to the Pan-STARRS1 iP1-band observations near maximum light, we estimate that the PS16cgx light curve peaked at 11 rest-frame days after the neutrino detection (corresponding to MJD 57 519.4) and at Mr ≈ −19.3 mag (Fig. 4a). Such an approach is often used to estimate the peak magnitudes of a SN (e.g. Suntzeff et al. 1999). Compared to the absolute peak magnitude ranges (corrected for host galaxy line-of-sight extinction) yielded by Drout et al. (2011), PS16cgx is brighter than normal Type Ic SNe with a peak magnitude of MR = −18.1 ± 0.6 mag, but consistent with Type Ic-BL SNe peaking at MR = −18.8 ± 1.1 mag (R-band magnitudes scaled to AB system). For comparison, Taddia et al. (2018) obtain an average peak magnitude of Mr = −17.66 ± 0.21 mag for Type Ic SNe. This is also reflected in our comparison of the absolute light curve of PS16cgx to a variety of nearby Type Ic SNe with well estimated explosion dates based on observational constraints or modelling (Fig. 4b). The comparison events include the GRB associated Type Ic-BL events SN 2006aj (Bianco et al. 2014) and SN 1998bw (Galama et al. 1998); non-GRB associated Type Ic-BL events SN 2016coi (Prentice et al. 2018), PTF10vgv (Corsi et al. 2012), SN 2009bb (Pignata et al. 2011), and SN 2002ap (Foley et al. 2003); and normal Type Ic events LSQ14efd (Barbarino et al. 2017), SN 2007gr (Hunter et al. 2009), and SN 2004aw (Taubenberger et al. 2006). The Johnson-Cousins light curves of SN 2007gr are also converted to the SDSS system using the transformations of Jester et al. (2005) in the comparison. Virgo and Great Attractor corrected host luminosity distances (via the NASA/IPAC Extragalactic Database) are used in comparisons if relevant.

thumbnail Fig. 4.

Panel a: observed light curves of PS16cgx. Zero epoch is set to the IceCube-160427A detection. A low-order polynomial light-curve fit is shown with a solid curve. Spectroscopic GMOS (30 days) and FORS2 (42 and 43 days) observations are indicated with vertical arrows. Panel b: PS16cgx r-band absolute light-curve comparison to a selection of nearby Type Ic SNe including the GRB associated Type Ic-BL events (SNe 2006aj and 1998bw); non-GRB associated Type Ic-BL events (PTF10vgv, and SNe 2016coi, 2009bb, and 2002ap); and normal Type Ic events (LSQ14efd, and SNe 2007gr and 2004aw). Panel c: selection of rapidly rising Type Ic SNe arbitrarily shifted as indicated in the figure legend. The comparison indicates an inconsistent spectrophotometric evolution of PS16cgx with Type Ic SNe which cannot be explained with extinction effects. Panel d: PS16cgx gri absolute light-curve comparison to a selection of similar Type Ia SNe. In addition, the slightly fainter r-band light curve of the canonical Type Ia SN 2011fe (Guillochon et al. 2017) is also included in the comparison. Furthermore, as a difference to the previous panels, the light curves of PS16cgx have been arbitrarily corrected for a modest host galaxy extinction of AV = 0.2 mag with RV = 1.27 derived by Burns et al. (2014) as an average value for their Type Ia sample. Assuming that the CFHT template reference image contains contamination from a Type Ia SN at mi = 24.5 mag, the final photometry point is also shown corrected for this effect with a square symbol in all panels.

Furthermore, most of the Type Ic SNe in our comparisons have slower rise times than PS16cgx would have to have to be associated with the IceCube-160427A event (assuming that neutrino emission is temporally coincident with core collapse or shock breakout). While some of the normal Type Ic SN light curves could be arbitrarily shifted to match the overall shape of the rest-frame r-band evolution of PS16cgx, the r − i colour evolution of events like SN 2006aj (Bianco et al. 2014) and SN 2007gr (Hunter et al. 2009) do not appear to be consistent with that of PS16cgx (Fig. 4c).

Contrary to the Type Ic SN comparisons, the absolute light curve evolution of PS16cgx is similar to that of many normal Type Ia SNe (Fig. 4d), for example SNe 2005eq, 2005hc, 2005na, 2006ax, and 2006py from the Carnegie Supernova Project (CSP; Contreras et al. 2010; Burns et al. 2014; Krisciunas et al. 2017). In particular, this is the case, if a small host galaxy extinction is assumed. Another approach to illustrate this is to compare the rest-frame colour curves of Type Ia and Ic SNe to the observations of PS16cgx (Fig. 5); this comparison suggests that the colours of PS16cgx are consistent Type Ia SNe, but bluer than normal Type Ic SNe. This discrepancy with Type Ic events cannot be explained with an underestimated host galaxy extinction of PS16cgx. Furthermore, while multi-band observations of PS16cgx are consistent with Type Ia SNe in the literature within errors, the match would be improved assuming a moderate host galaxy extinction of around AV ≈ 0.2 mag.

thumbnail Fig. 5.

Panel a: rest-frame Mg − Mr and Panel b: Mr − Mi colours of PS16cgx compared to the evolution of Type Ia SNe (solid curves; Burns et al. 2014; Krisciunas et al. 2017) including both normal and 91T-like events, and normal Type Ic SNe (dotted curves; Stritzinger et al. 2018a,b; Taddia et al. 2018) from CSP. Colour points assuming no host galaxy extinction are shown with open symbols, and assuming an arbitrary host extinction of AV = 0.2 mag in grey symbols.

To further quantify the compatibility of PS16cgx with the Type Ia SN hypothesis, we fit the riz-band light curve with the 4-parameter SALT2 template (Guy et al. 2005, 2007), assuming a Galactic extinction of E(B − V) = 0.054 mag derived from Finkbeiner et al. (1999) and a host redshift of z = 0.2895. The best-fit parameters were t0 = 57 517.55 ± 0.59 (MJD, i.e. 2016 May 9.15 UT) for peak, x0 = (4.39 ± 0.71)×10−5, x1 = −0.48 ± 0.61, and c = −0.04 ± 0.13. The allowed range for the x0 parameter corresponds to a distance modulus of μ = 40.89  ±  0.18, perfectly compatible with the ΛCDM prediction of 40.89 at z = 0.2895. Similarly, the x1 (stretch) and c (colour) parameters are compatible with the mean values (0) for Type Ia SNe. Furthermore, we studied the impact of colour information by testing the agreement between the data and template in each band. Since 19 of the 22 observations were taken in the i band, this band anchors the fit. However, the photometric points in both the r and z band fall slightly below the model expectation, which is an effect that cannot be explained with a simple reddening law. Quantitatively, the χ2/d.o.f. for the fit is 19.7/18 (p = 0.35), and the riz light curve can be considered compatible with the template to within photometric errors. To summarize, we conclude that there are no specific arguments to support a Type Ic classification of PS16cgx based on the available data. However the observations are in reasonable agreement with the expected Type Ia SN light curve and colour evolution. While some theoretical models (Wright et al. 2016, 2017) predict that MeV neutrinos from a nearby Type Ia SN would be detectable by IceCube, no neutrino emission is expected at TeV energies and above. Therefore, we find no plausible physical scenario where the neutrino event could be related to PS16cgx if it were indeed a Type Ia SN and conclude that in that case the two detections must be unrelated.

The host galaxy SDSS J160118.56+095153.7 of PS16cgx is an extended source in our images, however, no clear structure is resolved. Based on the Pan-STARRS1 DR1, the host has K-corrected absolute magnitudes of Mg ≈ −20.7 and Mr ≈ −21.2 mag. The luminosity-metallicity relation of Tremonti et al. (2004) suggests 12 + log(O/H)≈9. The metallicity using Pettini & Pagel (2004) N2 method is 12 + log(O/H)≈8.7 yielded from the VLT spectrum, consistent with the mass-luminosity relation. Therefore, the host has relatively normal, solar-like, metallicity.

4.3. Probabilities

The probability of a chance detection of PS16cgx with the IceCube-160427A event can be estimated following a similar approach to that of Aartsen et al. (2015) and Smartt et al. (2016c). Both cases of a Type Ia or Type Ic SN are considered. Assuming that a young and very distant SN can be discovered and that it can be identified to be in rise, the detection would likely have to take place roughly 0.5 mag below the peak. The brightest events we consider are the Type Ic-BL SNe which could peak roughly at M = −22.1 mag assuming a 3σ deviation from an average magnitude (Drout et al. 2011). Therefore, with the Pan-STARRS1 survey limit of ∼22.5 mag, such events could be discovered at rise up to a luminosity distance of Dl ≈ 6610 Mpc (or z ≈ 1.00). However, the majority of the events within this volume are too faint to be discovered. The review on the cosmic star formation rate by Madau & Dickinson (2014) yields a CCSN rate of RCCSN = 2.0  ×  10−4 Mpc−3 yr−1 scaled to the redshift of PS16cgx, which very roughly corresponds also to a comoving mid-volume distance from within an average magnitude Type Ic-BL or a normal Type Ia SN could be discovered with this survey. We adopt this CCSN rate as a rough average for the calculations. An intrinsic number of exploding CCSNe can be roughly estimated from

N SN = 4 3 π D m 3 R SN A 41253 deg 2 t 365 d yr , $$ \begin{aligned} N_{\mathrm{SN} } = \frac{4}{3} \pi D_{\rm m}^{3} R_{\mathrm{SN} } \frac{A}{41253\,\mathrm{deg} ^{2}} \frac{t}{365\,\mathrm{d} }\,\mathrm{yr} , \end{aligned} $$(1)

where Dm is the comoving radial distance Dm = Dl/(1 + z). With a sky area given by the neutrino error region of A = 0.78 deg2 and a survey time window of t = 5 d, this suggests NCCSN = 7.8. For the fractions of normal Type Ic and Type Ic-BL SNe from all CCSNe in a volume-limited sample, Li et al. (2011) derive fIc = 14% and fIc − BL = 2%, respectively (see also Smith et al. 2011), suggesting NIc = 1.1 and NIc − BL = 0.16. To estimate which fraction of the events are bright enough to be discovered with the survey, two effects have to be taken into account. Firstly, Mattila et al. (2012) estimates that only 77% of CCSNe at the redshift of PS16cgx are observable by optical surveys. This is due to a missing fraction of CCSNe with very high localised host galaxy extinctions way above the ‘normal’ extinction distribution in high star formation rate environments. Such correction reduces the discrepancy between the directly observed CCSN rates and the predicted star formation rates (Horiuchi et al. 2011), and has been applied by CCSN rate studies (e.g. Dahlen et al. 2012; Melinder et al. 2012). Secondly, a transient that is possibly associated with a neutrino detection is a young SN discovered at > 0.5 mag below the following peak magnitudes: MIc = −18.1 mag and 1σ deviation σIc = 0.6 mag, or MIc − BL = −18.8 mag and σIc − BL = 1.1 mag (Drout et al. 2011). Furthermore, we assume a normal rest-frame r-band host extinction distribution for SNe based by Riello & Patat (2005) and scaled as in Mattila et al. (2012). Using the aforementioned limiting magnitude of ∼22.5 mag, this yields a detection of only 1.4% of normal Type Ic SNe within the comoving volume. Combining the detection and missing fraction corrections results in an observable NIc = 0.012, corresponding to a Poisson probability p = 1 − eN of pIc = 1.2%. Similarly, a normal detection efficiency of 4.7% can be yielded for Type Ic-BL SNe, which results in an observable NIc − BL = 0.006 and pIc − BL = 0.6% for a chance detection.

The Type Ia SN rate at the redshift of PS16cgx is RIa = 0.4 × 10−4 Mpc−3 yr−1 (Dilday et al. 2010). For the aforementioned large comoving volume this suggests NIa = 1.6. The host galaxy extinction distribution for Type Ia SNe is adopted from (Dilday et al. 2010). No additional missing fraction is assumed for Type Ia SNe that are rising from old stellar population. Luminosity functions at peak for Type Ia SNe are adopted from Melinder et al. (2011) assuming BR =0 mag and shifted to the AB system (see Table 5). These result in detection efficiencies at > 0.5 mag below peak of 4.7, 6.1, 0.9, and 3.4% for normal, 91T-like, 91bg-like, and faint Type Ia SNe, respectively. This yields a total detection efficiency of 4.4% within the volume and thus a number of detectable Type Ia SNe of NIa = 0.070 corresponding to a chance detection of pIa = 7%. Therefore, a chance detection of a young and unrelated Type Ia SN in the Pan-STARRS1 follow up of an IceCube neutrino alert sky region is not that unlikely. The adopted and yielded values of all the above probability calculations including different Type Ia subtypes are shown for clarity in Table 5.

Table 5.

Adopted and yielded values of the probability calculations of a chance discovery for the IceCube-160427A error region.

5. Summary of other Pan-STARRS1 observations of IceCube alerts

Over the two-year period 2016 April 1–2018 March 31, IceCube released 13 alerts of either EHE or HESE events. Of these thirteen, ten were visible in the night time sky (the other three were in solar conjunction) and with Pan-STARRS1 we observed the field of 5 of them (see Table 6 for a summary). Since usually the initially released neutrino direction changes slightly after application of more time consuming reconstructions, we will use the latest available reconstruction results as our final set of coordinates. We typically observed with a similar observing strategy as discussed here in Sect. 4 with dithered exposures stacked together to make a night stack (as for the Pan-STARRS1 Medium Deep Fields). The size of the Pan-STARRS1 camera field of view and typical sky localisation regions from the EHE and HESE tracks from IceCube are well matched and we typically covered the neutrino error region to a depth of miP1 ≈ 22 mag. Observations of the field started between a few hours to a few days after the IceCube alert and continued for 3 to 30 days afterwards. We were motivated to look both for rapidly evolving transients (timescales of order 24−48 h) and SN-like transients (timescales of order a few to 10 days). As listed in Table 6, we typically found 10−20 transients per field, with magnitude ranges and light curves broadly similar to those detailed for IceCube-160427A in Sect. 4. However, we found no clear candidate for a SN-like transient with a rising light curve (or indeed any other characteristic) that would point towards a temporal coincidence. PS16cgx was the best candidate for a SN with an explosion epoch coincident with the neutrino time, and we consider that to be disfavoured due to its likely thermonuclear origin. The numbers of transients that we found are entirely consistent with the typical SN field population in a Pan-STARRS1 footprint at this depth (the vast majority are old SNe) and variable QSOs and AGNs (see Appendix A).

Table 6.

IceCube HESE and EHE events and Pan-STARRS1 observations from 2016−2017.

We cannot make a definitive statement that all transients discovered can be discounted as counterparts. We can only say that none of them show any unusual characteristics in their light curves or absolute magnitudes (when a host galaxy redshift exists) or in spectral characteristics (when a spectrum has been taken, see Appendix A) that suggest they are candidate sources for the counterparts of any of the HESE or EHE neutrinos. In other words, we are confusion limited within the ∼1 − 2 deg2 localisation regions of IceCube. Given the number of sources detected, we have only the option of searching for an unusual observation characteristic or implied temporal coincidence. Either of these significantly impact the probability calculation of coincidence as shown in Sect. 4.3.

We can make a qualitative statement that we find no evidence of transients above miP1 ≈ 22 mag in the error regions. However the purity of the neutrino alerts is between 25−50% and we have not been able to cover all five events from 1 h to 10 days uniformly. Therefore we cannot rule out fast transients of timescales ∼min to 24 h or fainter transients with longer timescales.

6. Discussion

The distance distribution of CCSNe follow the star-formation rate. If we assume those are the sources of the diffuse neutrino signal, the median redshift of sources of high-energy neutrino singlets is at z ∼ 1. For example, with Pan-STARRS1 image depth of ∼22.5 mag typical normal Type Ic and Type Ic-BL SNe can be detected in rise up to z = 0.2 − 0.4. However, superluminous Type Ic events could be detected up to z ≈ 1. Out of the source, which produce a single detectable high-energy neutrino alert in IceCube, only 13% are closer than redshift of z = 0.29 (see Fig. 6). Finding a distant neutrino source is therefore not surprising. Due to the Eddington bias we cannot estimate the neutrino flux from the detection of a single high-energy neutrino event by using Poisson statistics. As a consequence of the Eddington bias, the detection of a single high-energy neutrino implies an expectation value much smaller than one (Strotjohann et al. 2019). Unlike previous searches for neutrino multiplets (Abbasi et al. 2012), the observation of single neutrinos do not constrain the distance or luminosity of the source very much. Therefore a deep optical follow-up is necessary to catch the counterparts of distant sources. In the following we calculate the probability to find an optical counterpart in the follow-up of single high-energy neutrinos as a function of the survey depth.

thumbnail Fig. 6.

Cumulative fraction, fd, of astrophysical neutrino sources with a single high-energy neutrino detection in IceCube, accumulated as a function of source redshift, assuming they follow the SFR (Madau & Dickinson 2014) given standard cosmology. The redshift of PS16cgx is marked with a vertical line. The maximum distance ranges where a normal Type Ic (squares) and Type Ic-BL (circles) SNe can be detected in rise are indicated assuming negligible line-of-sight extinction, ignoring K-corrections, and adopting 1σ distribution of the SN peak magnitudes, for image depths of 22.5 and 25 mag (white and grey symbols, respectively).

We focus on HESE and EHE track-like alerts for which we can assume that their location in the sky will be confined to within 1 deg2 uncertainty. The rate of HESE/EHE neutrinos produced in a redshift shell is nearly independent of the redshift up to z ∼ 0.6 if we assume that the redshift distribution of the sources (SN of some sort) follow the star-formation rate (SFR). Figure 6 shows the accumulated fraction of sources as a function of redshift. We expect about four HESE and another four EHE events per year, of which on average about one and two, respectively, should be of true astrophysical origin. With a rate of Rν = 3 yr−1 astrophysical track-like events we can expect to detect a true association with a rate Rν ⋅ fd ∼ 0.4 yr−1 up to a redshift 0.3 assuming that all IceCube alerts can be followed up.

Accordingly, a deeper search for optical counterparts is well motivated5. Extending the redshift range, for instance from 0.3 to 0.6 will nearly double the number of source identifications.

7. Conclusion

We presented the properties of the first IceCube high-energy starting event distributed to the public, IceCube-160427A. The neutrino direction was observed by various instruments in gamma-ray and optical wavelengths. The only compelling optical counterpart found was PS16cgx. Our detailed investigations show some ambiguity in the classification of this SN. While the optical spectra could be explained by both Type Ia and Ic scenarios, the light curve and colour evolution points to PS16cgx being a Type Ia SNe and therefore likely unrelated to IceCube-160427A. The higher chance probability for Type Ia SNe supports this scenario. Therefore, we report a Pan-STARRS1 follow-up non-detection of an optical counterpart for IceCube-160427A within the 90% confidence region down to a 5σ limit of miP1 ≈ 22.5 mag.

Our findings show that spectroscopic follow-up and multi-colour photometric coverage of the light curves of candidate neutrino sources are crucial. Especially, spectral coverage during the SNe peak would have made a classification easier, since at late times Type Ic and Type Ia spectra look somewhat similar.

We also find no other plausible candidate for an optical transient coincident with four IceCube EHE/HESE triggers down to the 5σ limit of miP1 ≈ 22 mag. A large fraction of potential CCSN counterparts of high-energy neutrino detections are not accessible by shallow surveys. Assuming that high-energy neutrinos are produced by transients that trace the star-formation rate and have an optical counterpart, we suggest to extend follow-up observations to surveys reaching redshifts of up to 0.6 to ensure that the volume encloses a high probability of detection of at least one counterpart per year.


4

IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation.

5

We note that his argument holds for sources with positive evolution. In the rare case of negative evolution, for example for tidal disruption events, a deeper search might not be beneficial.

Acknowledgments

Pan-STARRS: Pan-STARRS is supported by the University of Hawaii and the NASA’s Planetary Defense Office under Grant No NNX14AM74G and the Pan-STARRS Multi-messenger science collaboration is supported by Queen’s University Belfast. The Pan-STARRS1 Surveys (PS1) have been made possible through contributions of the Institute for Astronomy, the University of Hawaii, the Pan-STARRS Project Office, the Max-Planck Society and its participating institutes, the Max Planck Institute for Astronomy, Heidelberg and the Max Planck Institute for Extraterrestrial Physics, Garching, The Johns Hopkins University, Durham University, the University of Edinburgh, Queen’s University Belfast, the Harvard-Smithsonian Center for Astrophysics, the Las Cumbres Observatory Global Telescope Network Incorporated, the National Central University of Taiwan, the Space Telescope Science Institute, the National Aeronautics and Space Administration under Grant No. NNX08AR22G issued through the Planetary Science Division of the NASA Science Mission Directorate, the National Science Foundation under Grant No. AST-1238877, the University of Maryland, and Eotvos Lorand University (ELTE). IceCube: The IceCube collaboration acknowledges the significant contributions to this manuscript from Anna Franckowiak, Claudio Kopper and Jakob van Santen. The IceCube collaboration gratefully acknowledge the support from the following agencies: USA – U.S. National Science Foundation-Office of Polar Programs, U.S. National Science Foundation-Physics Division, Wisconsin Alumni Research Foundation, Center for High Throughput Computing (CHTC) at the University of Wisconsin-Madison, Open Science Grid (OSG), Extreme Science and Engineering Discovery Environment (XSEDE), U.S. Department of Energy-National Energy Research Scientific Computing Center, Particle astrophysics research computing center at the University of Maryland, Institute for Cyber-Enabled Research at Michigan State University, and Astroparticle physics computational facility at Marquette University; Belgium – Funds for Scientific Research (FRS-FNRS and FWO), FWO Odysseus and Big Science programmes, and Belgian Federal Science Policy Office (Belspo); Germany – Bundesministerium für Bildung und Forschung (BMBF), Deutsche Forschungsgemeinschaft (DFG), Helmholtz Alliance for Astroparticle Physics (HAP), Initiative and Networking Fund of the Helmholtz Association, Deutsches Elektronen Synchrotron (DESY), and High Performance Computing cluster of the RWTH Aachen; Sweden – Swedish Research Council, Swedish Polar Research Secretariat, Swedish National Infrastructure for Computing (SNIC), and Knut and Alice Wallenberg Foundation; Australia – Australian Research Council; Canada – Natural Sciences and Engineering Research Council of Canada, Calcul Québec, Compute Ontario, Canada Foundation for Innovation, WestGrid, and Compute Canada; Denmark – Villum Fonden, Danish National Research Foundation (DNRF), Carlsberg Foundation; New Zealand – Marsden Fund; Japan – Japan Society for Promotion of Science (JSPS) and Institute for Global Prominent Research (IGPR) of Chiba University; Korea – National Research Foundation of Korea (NRF); Switzerland – Swiss National Science Foundation (SNSF). Others: We thank the anonymous referee for useful comments. EK, SJS and KS acknowledge support from the Science and Technology Facilities Council (STFC; ST/P000312/1, ST/N002520/1). SJS acknowledges funding from the European Research Council under the European Union’s Seventh Framework Programme (FP7/2007-2013)/ERC Grant agreement no [291222]. TWC acknowledgments the funding provided by the Alexander von Humboldt Foundation. AF thanks Irene Tamborra for fruitful discussions on neutrino production in supernovae. AF is supported by the Initiative and Networking Fund of the Helmholtz Association. This work is based (in part) on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere, Chile as part of PESSTO, (the Public ESO Spectroscopic Survey for Transient Objects Survey) ESO program 188.D-3003 and the DDT programme 297.D-5024. Based on observations obtained at the Gemini Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the National Science Foundation (United States), the National Research Council (Canada), CONICYT (Chile), Ministerio de Ciencia, Tecnología e Innovación Productiva (Argentina), and Ministério da Ciência, Tecnologia e Inovação (Brazil). Observations were carried out under programme GN-2016A-DD-9. Based on observations obtained with MegaPrime/MegaCam, a joint project of CFHT and CEA/IRFU, at the Canada-France-Hawaii Telescope (CFHT) which is operated by the National Research Council (NRC) of Canada, the Institut National des Science de l’Univers of the Centre National de la Recherche Scientifique (CNRS) of France, and the University of Hawaii. This work is based in part on data products produced at Terapix available at the Canadian Astronomy Data Centre as part of the Canada-France-Hawaii Telescope Legacy Survey, a collaborative project of NRC and CNRS. We have made use of the Weizmann interactive supernova data repository – (http://wiserep.weizmann.ac.il). This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. Funding for SDSS-III has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, and the U.S. Department of Energy Office of Science. The SDSS-III web site is (http://www.sdss3.org/). SDSS-III is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSS-III Collaboration including the University of Arizona, the Brazilian Participation Group, Brookhaven National Laboratory, Carnegie Mellon University, University of Florida, the French Participation Group, the German Participation Group, Harvard University, the Instituto de Astrofisica de Canarias, the Michigan State/Notre Dame/JINA Participation Group, Johns Hopkins University, Lawrence Berkeley National Laboratory, Max Planck Institute for Astrophysics, Max Planck Institute for Extraterrestrial Physics, New Mexico State University, New York University, Ohio State University, Pennsylvania State University, University of Portsmouth, Princeton University, the Spanish Participation Group, University of Tokyo, University of Utah, Vanderbilt University, University of Virginia, University of Washington, and Yale University.

References

  1. Aartsen, M. G., Abbasi, R., Abdou, Y., et al. 2013a, Science, 342, 1242856 [CrossRef] [PubMed] [Google Scholar]
  2. Aartsen, M. G., Abbasi, R., Abdou, Y., et al. 2013b, Nucl. Instrum. Meth. A, 711, 73 [Google Scholar]
  3. Aartsen, M. G., Ackermann, M., Adams, J., et al. 2014, Phys. Rev. Lett., 113, 101101 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  4. Aartsen, M. G., Abraham, K., Ackermann, M., et al. 2015, ApJ, 811, 52 [NASA ADS] [CrossRef] [Google Scholar]
  5. Aartsen, M. G., Abraham, K., Ackermann, M., et al. 2016, ApJ, 824, 115 [CrossRef] [Google Scholar]
  6. Aartsen, M. G., Abraham, K., Ackermann, M., et al. 2017a, ApJ, 835, 151 [NASA ADS] [CrossRef] [Google Scholar]
  7. Aartsen, M. G., Ackermann, M., Adams, J., et al. 2017b, Astropart. Phys., 92, 30 [NASA ADS] [CrossRef] [Google Scholar]
  8. Aartsen, M. G., Ackermann, M., Adams, J., et al. 2017c, ArXiv e-prints [arXiv:astro-ph.HE/1710.01191] [Google Scholar]
  9. Aartsen, M. G., Ackermann, M., Adams, J., et al. 2017d, A&A, 607, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Aartsen, M. G., Ackermann, M., Adams, J., et al. 2018a, Science, 361, 147 [NASA ADS] [Google Scholar]
  11. Aartsen, M. G., Ackermann, M., Adams, J., et al. 2018b, Science, 361, eaat1378 [NASA ADS] [Google Scholar]
  12. Abbasi, R., Ackermann, M., Adams, J., et al. 2009, Nucl. Instrum. Meth. A, 601, 294 [Google Scholar]
  13. Abbasi, R., Abdou, Y., Abu-Zayyad, T., et al. 2012, A&A, 539, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2016, Phys. Rev. Lett., 116, 061102 [Google Scholar]
  15. Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2017, ApJ, 848, L12 [Google Scholar]
  16. Alam, S., Albareti, F. D., Allende Prieto, C., et al. 2015, ApJS, 219, 12 [NASA ADS] [CrossRef] [Google Scholar]
  17. Ando, S., & Beacom, J. F. 2005, Phys. Rev. Lett., 95, 061103 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  18. Appenzeller, I., Fricke, K., Fürtig, W., et al. 1998, The Messenger, 94, 1 [NASA ADS] [Google Scholar]
  19. Barbarino, C., Botticella, M. T., Dall’Ora, M., et al. 2017, MNRAS, 471, 2463 [NASA ADS] [CrossRef] [Google Scholar]
  20. Bianco, F. B., Modjaz, M., Hicken, M., et al. 2014, ApJS, 213, 19 [NASA ADS] [CrossRef] [Google Scholar]
  21. Biland, A., & Dorner, D. 2016, GRB Coord. Netw., 19427, 1 [Google Scholar]
  22. Bionta, R. M., Blewitt, G., Bratton, C. B., et al. 1987, Phys. Rev. Lett., 58, 1494 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  23. Blackburn, L., Briggs, M. S., Burns, E., et al. 2016, GRB Coord. Netw., 19364, 1 [Google Scholar]
  24. Blaufuss, E. 2016, GRB Coord. Netw., 19363, 1 [Google Scholar]
  25. Blondin, S., Matheson, T., Kirshner, R. P., et al. 2012, AJ, 143, 126 [NASA ADS] [CrossRef] [Google Scholar]
  26. Botticella, M. T., Trundle, C., Pastorello, A., et al. 2010, ApJ, 717, L52 [NASA ADS] [CrossRef] [Google Scholar]
  27. Boulade, O., Charlot, X., Abbon, P., et al. 2003, Proc. SPIE, 4841, 72 [NASA ADS] [CrossRef] [Google Scholar]
  28. Burns, C. R., Stritzinger, M., Phillips, M. M., et al. 2014, ApJ, 789, 32 [NASA ADS] [CrossRef] [Google Scholar]
  29. Chambers, K. C., Magnier, E. A., Metcalfe, N., et al. 2016, ArXiv e-prints [arXiv:1612.05560] [Google Scholar]
  30. Chirkin, D., & Rhode, W. 2004, ArXiv e-prints [arXiv:hep-ph/0407075] [Google Scholar]
  31. Contreras, C., Hamuy, M., Phillips, M. M., et al. 2010, AJ, 139, 519 [NASA ADS] [CrossRef] [Google Scholar]
  32. Corsi, A., Ofek, E. O., Gal-Yam, A., et al. 2012, ApJ, 747, L5 [NASA ADS] [CrossRef] [Google Scholar]
  33. Crowther, P. A. 2007, ARA&A, 45, 177 [NASA ADS] [CrossRef] [Google Scholar]
  34. Dahlen, T., Strolger, L.-G., Riess, A. G., et al. 2012, ApJ, 757, 70 [NASA ADS] [CrossRef] [Google Scholar]
  35. Dilday, B., Smith, M., Bassett, B., et al. 2010, ApJ, 713, 1026 [NASA ADS] [CrossRef] [Google Scholar]
  36. Drout, M. R., Soderberg, A. M., Gal-Yam, A., et al. 2011, ApJ, 741, 97 [NASA ADS] [CrossRef] [Google Scholar]
  37. Esmaili, A., & Murase, K. 2018, JCAP, 1812, 008 [NASA ADS] [CrossRef] [Google Scholar]
  38. Evans, P. A., Osborne, J. P., Kennea, J. A., et al. 2015, MNRAS, 448, 2210 [NASA ADS] [CrossRef] [Google Scholar]
  39. Filippenko, A. V. 1997, ARA&A, 35, 309 [NASA ADS] [CrossRef] [Google Scholar]
  40. Flesch, E. W. 2015, PASA, 32, e010 [NASA ADS] [CrossRef] [Google Scholar]
  41. Finkbeiner, D. P., Davis, M., & Schlegel, D. J. 1999, ApJ, 524, 867 [NASA ADS] [CrossRef] [Google Scholar]
  42. Folatelli, G., Morrell, N., Phillips, M. M., et al. 2013, ApJ, 773, 53 [NASA ADS] [CrossRef] [Google Scholar]
  43. Foley, R. J., Papenkova, M. S., Swift, B. J., et al. 2003, PASP, 115, 1220 [NASA ADS] [CrossRef] [Google Scholar]
  44. Galama, T. J., Vreeswijk, P. M., van Paradijs, J., et al. 1998, Nature, 395, 670 [NASA ADS] [CrossRef] [Google Scholar]
  45. Guillochon, J., Parrent, J., Kelley, L. Z., & Margutti, R. 2017, ApJ, 835, 64 [NASA ADS] [CrossRef] [Google Scholar]
  46. Guy, J., Astier, P., Nobili, S., Regnault, N., & Pain, R. 2005, A&A, 443, 781 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Guy, J., Astier, P., Baumont, S., et al. 2007, A&A, 466, 11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Haack, C., Wiebusch, C., & IceCube Collaboration. 2017, Int. Cosmic Ray Conf., 301, 1005 [NASA ADS] [CrossRef] [Google Scholar]
  49. Hirata, K., Kajita, T., Koshiba, M., et al. 1987, Phys. Rev. Lett., 58, 1490 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  50. Hook, I. M., Jørgensen, I., Allington-Smith, J. R., et al. 2004, PASP, 116, 425 [NASA ADS] [CrossRef] [Google Scholar]
  51. Horiuchi, S., Beacom, J. F., Kochanek, C. S., et al. 2011, ApJ, 738, 154 [NASA ADS] [CrossRef] [Google Scholar]
  52. Hunter, D. J., Valenti, S., Kotak, R., et al. 2009, A&A, 508, 371 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. IceCube Collaboration 2017a, GRB Coord. Netw., 22016, 1 [Google Scholar]
  54. IceCube Collaboration 2017b, GRB Coord. Netw., 22065, 1 [Google Scholar]
  55. IceCube Collaboration 2017c, GRB Coord. Netw., 22105, 1 [Google Scholar]
  56. Inserra, C., Smartt, S. J., Gall, E. E. E., et al. 2018, MNRAS, 475, 1046 [NASA ADS] [CrossRef] [Google Scholar]
  57. Jester, S., Schneider, D. P., Richards, G. T., et al. 2005, AJ, 130, 873 [NASA ADS] [CrossRef] [Google Scholar]
  58. Kistler, M. D., Haxton, W. C., & Yuksel, H. 2013, ApJ, 778, 81 [NASA ADS] [CrossRef] [Google Scholar]
  59. Kopper, C., & IceCube Collaboration 2017, Int. Cosmic Ray Conf., 301, 981 [NASA ADS] [CrossRef] [Google Scholar]
  60. Kopper, C., Giang, W., & Kurahashi, N. 2015, 34th Int. Cosmic Ray Conf. (ICRC2015), 34, 1081 [NASA ADS] [Google Scholar]
  61. Kowalski, M., & Mohr, A. 2007, Astropart. Phys., 27, 533 [NASA ADS] [CrossRef] [Google Scholar]
  62. Krisciunas, K., Contreras, C., Burns, C. R., et al. 2017, AJ, 154, 211 [NASA ADS] [CrossRef] [Google Scholar]
  63. Lantz, B., Aldering, G., Antilogus, P., et al. 2004, Proc. SPIE, 5249, 146 [NASA ADS] [CrossRef] [Google Scholar]
  64. Li, W., Leaman, J., Chornock, R., et al. 2011, MNRAS, 412, 1441 [NASA ADS] [CrossRef] [Google Scholar]
  65. Lipunov, V., Gorbovskoy, E., Kornilov, V., et al. 2016a, GRB Coord. Netw., 19362, 1 [Google Scholar]
  66. Lipunov, V., Gorbovskoy, E. S., Tyurina, N., et al. 2016b, GRB Coord. Netw., 19748, 1 [Google Scholar]
  67. Lunnan, R., Chornock, R., Berger, E., et al. 2018, ApJ, 852, 81 [NASA ADS] [CrossRef] [Google Scholar]
  68. Madau, P., & Dickinson, M. 2014, ARA&A, 52, 415 [NASA ADS] [CrossRef] [Google Scholar]
  69. Magnier, E. A., Chambers, K. C., Flewelling, H. A., et al. 2016a, ArXiv e-prints [arXiv:1612.05240] [Google Scholar]
  70. Magnier, E. A., Schlafly, E. F., Finkbeiner, D. P., et al. 2016b, ArXiv e-prints [arXiv:1612.05242] [Google Scholar]
  71. Mattila, S., Dahlen, T., Efstathiou, A., et al. 2012, ApJ, 756, 111 [NASA ADS] [CrossRef] [Google Scholar]
  72. McBrien, O., Smith, K. W., Smartt, S. J., et al. 2017, GRB Coord. Netw., 22173, 1 [Google Scholar]
  73. McCrum, M., Smartt, S. J., Rest, A., et al. 2015, MNRAS, 448, 1206 [NASA ADS] [CrossRef] [Google Scholar]
  74. Meagher, K., & IceCube Collaboration 2017, Int. Cosmic Ray Conf. Proc., 35, 1007 [NASA ADS] [CrossRef] [Google Scholar]
  75. Melinder, J., Dahlen, T., Mencía-Trinchant, L., et al. 2011, A&A, 532, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  76. Melinder, J., Dahlen, T., Mencía Trinchant, L., et al. 2012, A&A, 545, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  77. Meszaros, P., & Waxman, E. 2001, Phys. Rev. Lett., 87, 171102 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  78. Modjaz, M., Blondin, S., Kirshner, R. P., et al. 2014, AJ, 147, 99 [NASA ADS] [CrossRef] [Google Scholar]
  79. Mukherjee, R. 2016, GRB Coord. Netw., 19377, 1 [Google Scholar]
  80. Murase, K., & Ioka, K. 2013, Phys. Rev. Lett., 111, 121102 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  81. Murase, K., Ioka, K., Nagataki, S., & Nakamura, T. 2006, ApJ, 651, L5 [NASA ADS] [CrossRef] [Google Scholar]
  82. Padovani, P., & Resconi, E. 2014, MNRAS, 443, 474 [NASA ADS] [CrossRef] [Google Scholar]
  83. Paragi, Z., Taylor, G. B., Kouveliotou, C., et al. 2010, Nature, 463, 516 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  84. Patat, F., Cappellaro, E., Danziger, J., et al. 2001, ApJ, 555, 900 [NASA ADS] [CrossRef] [Google Scholar]
  85. Pereira, R., Thomas, R. C., Aldering, G., et al. 2013, A&A, 554, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  86. Petropoulou, M., Dimitrakoudis, S., Padovani, P., Mastichiadis, A., & Resconi, E. 2015, MNRAS, 448, 2412 [NASA ADS] [CrossRef] [Google Scholar]
  87. Pettini, M., & Pagel, B. E. J. 2004, MNRAS, 348, L59 [NASA ADS] [CrossRef] [Google Scholar]
  88. Pignata, G., Stritzinger, M., Soderberg, A., et al. 2011, ApJ, 728, 14 [NASA ADS] [CrossRef] [Google Scholar]
  89. Piran, T. 1999, Phys. Rep., 314, 575 [NASA ADS] [CrossRef] [Google Scholar]
  90. Podsiadlowski, P., Joss, P. C., & Hsu, J. J. L. 1992, ApJ, 391, 246 [NASA ADS] [CrossRef] [Google Scholar]
  91. Prentice, S. J., Ashall, C., Mazzali, P. A., et al. 2018, MNRAS, 478, 4162 [NASA ADS] [CrossRef] [Google Scholar]
  92. Razzaque, S., Mészáros, P., & Waxman, E. 2005, Mod. Phys. Lett. A, 20, 2351 [NASA ADS] [CrossRef] [Google Scholar]
  93. Rest, A., Scolnic, D., Foley, R. J., et al. 2014, ApJ, 795, 44 [NASA ADS] [CrossRef] [Google Scholar]
  94. Reynolds, T., Fraser, M., Mattila, S., et al. 2016, ATel, 9337, 1 [NASA ADS] [Google Scholar]
  95. Richards, G. T., Myers, A. D., Gray, A. G., et al. 2009, ApJS, 180, 67 [NASA ADS] [CrossRef] [Google Scholar]
  96. Riello, M., & Patat, F. 2005, MNRAS, 362, 671 [NASA ADS] [CrossRef] [Google Scholar]
  97. Scalzo, R. A., Yuan, F., Childress, M. J., et al. 2017, PASA, 34, e030 [NASA ADS] [CrossRef] [Google Scholar]
  98. Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 [NASA ADS] [CrossRef] [Google Scholar]
  99. Scolnic, D. M., Jones, D. O., Rest, A., et al. 2018, ApJ, 859, 101 [NASA ADS] [CrossRef] [Google Scholar]
  100. Senno, N., Murase, K., & Mészáros, P. 2016, Phys. Rev. D, 93, 083003 [NASA ADS] [CrossRef] [Google Scholar]
  101. Singer, L. P., Kasliwal, M. M., Cenko, S. B., et al. 2016a, GRB Coord. Netw., 19392, 1 [Google Scholar]
  102. Singer, L. P., Kasliwal, M. M., Bhalerao, V., et al. 2016b, GRB Coord. Netw., 19760, 1 [Google Scholar]
  103. Smartt, S. J., Valenti, S., Fraser, M., et al. 2015, A&A, 579, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  104. Smartt, S., Chambers, K. C., Smith, K., et al. 2016a, GRB Coord. Netw., 19381, 1 [Google Scholar]
  105. Smartt, S. J., Chambers, K. C., Smith, K. W., et al. 2016b, MNRAS, 462, 4094 [NASA ADS] [CrossRef] [Google Scholar]
  106. Smartt, S. J., Chambers, K. C., Smith, K. W., et al. 2016c, ApJ, 827, L40 [NASA ADS] [CrossRef] [Google Scholar]
  107. Smith, N., Li, W., Filippenko, A. V., & Chornock, R. 2011, MNRAS, 412, 1522 [NASA ADS] [CrossRef] [Google Scholar]
  108. Smith, M. W. E., Fox, D. B., Cowen, D. F., et al. 2013, Astropart. Phys., 45, 56 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  109. Stritzinger, M. D., Anderson, J. P., Contreras, C., et al. 2018a, A&A, 609, A134 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  110. Stritzinger, M. D., Taddia, F., Burns, C. R., et al. 2018b, A&A, 609, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  111. Strotjohann, N., Kowalski, M., Burns, C. R., et al. 2019, A&A, 622, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  112. Suntzeff, N. B., Phillips, M. M., Covarrubias, R., et al. 1999, AJ, 117, 1175 [NASA ADS] [CrossRef] [Google Scholar]
  113. Svinkin, D., Golenetskii, S., Aptekar, R., et al. 2016, GRB Coord. Netw., 19426, 1 [Google Scholar]
  114. Taboada, I. 2016, GRB Coord. Netw., 19361, 1 [Google Scholar]
  115. Taddia, F., Stritzinger, M. D., Bersten, M., et al. 2018, A&A, 609, A136 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  116. Tamborra, I., & Ando, S. 2016, Phys. Rev. D, 93, 053010 [NASA ADS] [CrossRef] [Google Scholar]
  117. Taubenberger, S., Pastorello, A., Mazzali, P. A., et al. 2006, MNRAS, 371, 1459 [NASA ADS] [CrossRef] [Google Scholar]
  118. Tremonti, C. A., Heckman, T. M., Kauffmann, G., et al. 2004, ApJ, 613, 898 [NASA ADS] [CrossRef] [Google Scholar]
  119. Vianello, G., Magill, J. D., Omodei, N., et al. 2016, GRB Coord. Netw., 19360, 1 [Google Scholar]
  120. Waxman, E. 1995, Phys. Rev. Lett., 75, 386 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  121. Waxman, E., & Bahcall, J. 1997, Phys. Rev. Lett., 78, 2292 [CrossRef] [Google Scholar]
  122. Wright, D. E., Smartt, S. J., Smith, K. W., et al. 2015, MNRAS, 449, 451 [NASA ADS] [CrossRef] [Google Scholar]
  123. Wright, W. P., Nagaraj, G., Kneller, J. P., et al. 2016, Phys. Rev. D, 94, 025026 [NASA ADS] [CrossRef] [Google Scholar]
  124. Wright, W. P., Kneller, J. P., Ohlmann, S. T., et al. 2017, Phys. Rev. D, 95, 043006 [NASA ADS] [CrossRef] [Google Scholar]
  125. Yaron, O., & Gal-Yam, A. 2012, PASP, 124, 668 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Details of Pan-STARRS1 observations of other IceCube Events

In the tables that follow, we provide the information on the most probable type for each transient found:

  • Stellar: the transient is spatially coincident with a faint source likely to be a star.

  • QSO: the transient is spatially coincident with the core of a galaxy that has previously been catalogued as an AGN or QSO.

  • SN: the transient is associated with a host galaxy at a plausible physical or angular separation.

  • NT: nuclear transient. The transient is coincident with the core of a galaxy that has not been previously catalogued as an AGN or QSO. These could be AGN, or nuclear SNe or even tidal disruption events.

  • Orphan: there is no plausible host galaxy or stellar counterpart detected either in the Pan-STARRS reference images or any known catalogue.

There were no transients that we found that had characteristics that made them peculiar from any others in the field and would lead to further investigation as to their association with the neutrinos. Typically, the IceCube sky locations have a 90% confidence region given by a likelihood scan with a radius of 0 . ° $ \overset{\circ }{.} $5−1 . ° $ \overset{\circ }{.} $2. However, these values depend on the sophistication of the IceCube track analysis. Therefore, for consistency and completeness, we report all transients in the tables below that were found within the Pan-STARRS1 camera footprint of around 1 . ° $ \overset{\circ }{.} $4 radius.

A.1. IceCube-160427A

The SN PS16cgx has been extensively discussed in the main manuscript. The rest of the transients discovered within the 90% confidence region of IceCube-160427A are shortly summarised below.

PS16cfv is associated with SDSS J160133.99+094759.0. This is a stellar object in SDSS.

PS16cgw is offset from the host SDSS J160039.70+093919.6 by 1 . $ \overset{\prime \prime }{.} $6 N and 0 . $ \overset{\prime \prime }{.} $2 E.

PS16chy is associated with a stellar object SDSS J160044.22+095040.6.

PS16cfy is matched with SDSS J160205.21+094630.3, which is a spectroscopic QSO at z = 2.344.

PS16dow is offset from the host SDSS J160202.18+094907.6 core by 0 . $ \overset{\prime \prime }{.} $1 which is Milliquas source at a photometric redshift of z = 1.325.

PS16dov is coincident with the spectroscopic QSO SDSS J160153.91+095244.3 at z = 2.282.

PS16cfu is offset by 0 . $ \overset{\prime \prime }{.} $1 N and 0 . $ \overset{\prime \prime }{.} $2 E from the host SDSS J160115.65+092504.8.

PS16cgb is coincident with a host SDSS J160219.12+093450.2 at a photometric redshift of z = 0.755. The transient light curve is flat; it could be an old SN or AGN variability.

PS16cfz is closely matched with the nucleus of SDSS J160211.95+095407.7 at a photometric redshift of z = 0.169. It is either a SN or AGN variability with a flat light curve.

PS16fns is coincident with SDSS J160124.77+100201.5; an AGN in Milliquas at a photometric redshift of z = 1.955.

PS16cfw is matched with the spectroscopic QSO SDSS J160139.09+092159.5 at z = 0.862.

PS16fnt is matched with the Milliquas AGN SDSS J160128.08+100250.9 at a photometric redshift z = 0.975.

PS16fnu is consistent with SDSS J155941.20+094619.4. This is a stellar object in SDSS.

PS16fnv is associated with SDSS J155941.54+094749.8. This is a source in the Milliquas catalogue at a photometric redshift of z = 1.575.

PS16fnw is coincident with SDSS J160153.37+100521.7, a spectroscopic QSO at z = 2.285.

PS16fnx is located 0 . $ \overset{\prime \prime }{.} $3 N and 0 . $ \overset{\prime \prime }{.} $2 W from SDSS J155928.45+095151.4. This is a faint galaxy in SDSS at a photometric redshift of z = 0.579.

PS16fny host SDSS J160227.74+100309.0 is a source in Milliquas at a photometric redshift of z = 0.375.

PS16fnz is coincident with SDSS J160119.77+101218.7, which is a Milliquas source at a photometric redshift of z = 0.225.

A.2. IceCube-160731A

IceCube-160731A was discovered on 31 July 2016 01:55:04.00 UT (MJD = 57 600.07990). Pan-STARRS1 began observing the field 31 July 2016 07:41:39 UT (MJD = 57 600.32) and continued to observe every night until 7 August 2016, and then took another 6 epochs up until 26 August 2016. This gave an observing baseline of 26 days. Table A.1 contains the candidate counterparts identified by the Pan-STARRS1 follow-up campaign, out to a 1 . ° $ \overset{\circ }{.} $45 radius. Three transient were spectroscopically classified.

Table A.1.

31 transients discovered within the local vicinity of the IceCube-160731A neutrino detection (MJD 57 600.07990).

PS16dpm was discovered on the first epoch of observing and was observed to have a rising light curve. A spectrum was taken with the University of Hawaii SNIFS spectrometer on 3 August 2016, showing it to be a convincing Type Ia at z = 0.12 around 6 days before maximum light. A PESSTO spectrum confirmed this result (Smartt 2015; Reynolds et al. 2016). The Pan-STARRS1 light curve and these two spectra imply a peak of MJD = 57 610 and an explosion epoch of MJD = 57 592 ± 2, or 8 days before IceCube-160731A.

PS16dpl was also discovered on the first epoch of observing and the light curve stayed flat for 7 days in the iP1 band. The host galaxy SDSS J142303.54+000552.1 has a spectroscopic redshift of z = 0.141. It was observed by SNIFS on 7 August 2016 when it was still miP1 = 20.8 ± 0.1 mag and shows features consistent with a Type Ia at this redshift.

PS16dpi was a known, old SN discovered by the SkyMapper transient survey as SMT16ojf (Scalzo 2017), it was an old Type Ia at z = 0.085115 (Reynolds et al. 2016) which exploded about 30 days before IceCube-160731A.

The rest of the transients are summarised below and none had any characteristic to make it appear possibly related to IceCube-160731A.

PS16cxe is not matched to any known or catalogued source.

PS16dpu has a SDSS host which is stellar.

6H3Ghhv is associated with SDSS J141721.29−003050.5. This is a faint and stellar object in SDSS.

6H3Ghhn is a marginal candidate, being faint and low signal-to-noise in zP1 and absent in iP1.

6H3Ghik is associated with SDSS J141718.86−003927.2 and is offset from the host core by 0 . $ \overset{\prime \prime }{.} $2.

PS16dqg is not matched to any known or catalogued source.

6H3Ghhq is associated with SDSS J141843.68−005204.3. This is a faint and stellar object in SDSS.

PS16dpe is not matched to any known or catalogued source, but appears to be a SN at peak.

PS16dpc is not matched to any known or catalogued source, but appears to be a genuine orphaned SN.

PS16dpv is not matched to any known source and has a flat light curve.

PS16dpd is coincident with a stellar like source that was classified in SDSS DR12 as a star, but Milliquas catalogue (Flesch 2015) suggests it is a QSO. In any event, it is not a SN candidate.

PS16dpf is not matched to any known source and has a flat light curve. There appears to be a faint host in the Pan-STARRS1 reference image making this a potential SN candidate.

6H3Ghhx is associated with SDSS J141445.76−002739.7 but has a projected offset of only 0.3 kpc, at the host redshift, implying it may be a nuclear transient.

PS16cyu is of offset from SDSS J141953.14−010041.5 and shows a declining light curve. It was previously discovered by Pan-STARRS on 28 June 2016 and is therefore an old SN.

6H3Ghhk is matched to SDSS J142101.10−004537.9. This is a faint and stellar object in SDSS.

PS16dpk, although not matched to a known source, showed a rising light curve after discovery marking it as a plausible SN candidate.

6H3Ghis is associated with SDSS J141506.06+001806.3, a stellar source at a photometric redshift of z ∼ 0.3.

PS16dpt has a SDSS host which is stellar.

6H3Ghhu is not matched to any known or catalogued source.

6H3Ghia is associated with SDSS J141417.75+000207.8 and is offset from the host by 0 . $ \overset{\prime \prime }{.} $10 N and 0 . $ \overset{\prime \prime }{.} $20 W

PS16dps is not matched to any known source, but does appear near an extended source in the Pan-STARRS1 footprint. This is a plausible SN candidate.

6H3Ghet is matched to SDSS J142152.68−004016.9 and is located 0 . $ \overset{\prime \prime }{.} $03 S and 0 . $ \overset{\prime \prime }{.} $08 W of the host core. At a redshift of 0.072, this gives a projected offset of approximately 10 kpc.

6H3Ghjd is associated with SDSS J141600.56−011920.2 and is offset from the host by 0 . $ \overset{\prime \prime }{.} $2, or a projected offset of 0.5 kpc at a redshift of 0.171.

PS16dpj is not matched to any catalogued source and has a relatively flat light curve, making discerning the nature of this transient difficult.

6H3Ghib is associated with SDSS J141318.12−002527.8 and is offset from the host by 0 . $ \overset{\prime \prime }{.} $16 S and 0 . $ \overset{\prime \prime }{.} $09 W.

6H3Ghij is possibly associated with SDSS J141442.38+004325.7, a stellar source, or a QSO in Milliquas at a redshift of z = 0.5.

PS16dqf appears to be a nuclear transient.

PS16dpr is associated with a host catalogued in Milliquas, making it a likely QSO at a photometric redshift of z = 2.043.

A.3. IceCube-161103A

IceCube-161103A was discovered on 3 November 2016 09:07:31.12 (MJD 57 695.38022). Pan-STARRS1 follow-up observations began the same night on 3 November 2016 11:09:55 UT (MJD 57 695.47) and proceeded for a subsequent 9 nights. Of these, 2 nights were lost due to weather. A total of 7 transients were found while searching the confidence region in this time, detailed in Table A.2. None of these events were spectroscopically classified.

Table A.2.

7 transients discovered within the local vicinity of the IceCube-161103A neutrino detection (MJD 57 695.38022).

PS16exc is not associated with any known or catalogued source and is faint, at ∼22 mag in iP1. The light curve remains relatively flat for almost 25 days.

PS16exd is not associated with any known or catalogued sources. The transient was observed for several days in the iP1 band before fading and it is a potential SN candidate at maximum. The host is a faint object, either a compact galaxy or star.

PS16exh has no associated catalogued host but a faint point source is visible in the Pan-STARRS1 reference stacks. This is a potential stellar source but appears more as a variable QSO.

PS16exi has no associated catalogued host but a faint point source is visible in the Pan-STARRS1 reference stacks. This is a potential stellar source but appears more as a variable QSO.

PS16exf has no associated catalogued host but a faint point source is visible in the Pan-STARRS1 reference stacks. This is likely a stellar source or variable QSO.

PS16exg is not associated with any known or catalogued sources. The Pan-STARRS1 exposures show a point source, coupled with a relatively flat light curve (declining < 0.5 mag in 44 days) suggests this may be an AGN/QSO or a stellar source.

6K3Ghjo is not associated with any known or catalogued source. It was observed initially in iP1 but also in zP1 from +40 days from discovery. Typically, these zP1 band observations were one magnitude brighter than the iP1 band data, but both sets were relatively flat, implying a stellar source. There is a clear point source at the transient position.

PS16fhi was found outside the 90% radius (at αJ2000 = 02h38m35 . s $ {{\overset{\text{ s}}{.}}} $74, δJ2000 = +13 ° 05 06 . 2 $ ^\circ 05\prime 06\mathop .\limits^{\prime \prime } 2 $) and is not listed in Table A.2. However it initially caught our attention as it began quickly rising. A Gemini GMOS spectrum and its rising light curve showed it to be a bright Type II at z = 0.085, resembling SN 2009kf (Botticella et al. 2010). We rule it out as being related as a combination of these data implied it exploded significantly after (up to 10−12 days) IceCube-161103A.

A.4. IceCube-170321A

IceCube-170321A was discovered on 21 March 2017 07:32:20.69 UT (MJD 57 833.31413). Pan-STARRS1 began following this event the next night on 22 March 2017 05:44:56 UT (MJD 57 834.24) and proceeded for a subsequent 3 nights. From the night of 27 March 2017 to 31 March 2017, observations were performed on alternate nights. Following this, a final night was dedicated to observing on 6 April 2017. Table A.3 contains the candidate counterparts identified by the Pan-STARRS1 follow-up campaign while searching the confidence region. A total of 7 transients were identified, though none were observed spectroscopically.

Table A.3.

7 transients discovered within the local vicinity of the IceCube-170321A neutrino detection (MJD 57 805.31413).

PS17ceb is not matched to any known or catalogued source, but does appear to be within a faint, extended source. This is most likely a galaxy, making this transient a SN candidate.

PS17cdy appears as a point source but is not matched to any known source. It is likely an uncatalogued stellar source.

PS17cdz is not matched to any known or catalogued source.

7C3Gjfw is not matched to any known or catalogued sources.

PS17cea is not matched to any known or catalogued source.

PS17byo is not matched to any known sources, but does show a faint potential host in the reference images. As such, this is likely a SN candidate.

PS17byn is not matched to any known or catalogued source, but its light curve shows evidence of some decline. It is possible that this is a SN candidate.

A.5. IceCube-171106A

The field of IceCube-171106A was observed starting from 8 November 2017 08:09:36.0 UT (MJD 58 065.34), some 37.44 hours following the detection and proceeded for the subsequent 10 nights. Observations were made using the iP1 and zP1 bands, with difference imaging being performed with respect to the Pan-STARRS1 3π stacked reference sky. 5σ limiting magnitudes of miP1 ≈ 22.5 mag were achieved. Table A.4 contains the candidate counterparts identified by the Pan-STARRS1 follow-up campaign. The initial results of the IceCube-171106A Pan-STARRS1 follow-up were reported by McBrien et al. (2017).

PS17eym is a SN candidate associated with SDSS J223922.20+071638.0. It is located 1 . $ \overset{\prime \prime }{.} $80 S and 0 . $ \overset{\prime \prime }{.} $97 W from the galaxy core, or a projected offset of 10.7 kpc.

Table A.4.

7 transients discovered within the local vicinity of the IceCube-171106A neutrino detection (MJD 58 063.77754).

PS17fcc is a SN candidate that was found associated with the galaxy SDSS J224158.66+074653.2, being located 0 . $ \overset{\prime \prime }{.} $59 N and 0 . $ \overset{\prime \prime }{.} $28 W from the galaxy core, or a projected offset of 1.9 kpc.

PS17eyn is a SN candidate associated with the galaxy SDSS J224224.58+065502.1. It is located 0 . $ \overset{\prime \prime }{.} $58 N and 0 . $ \overset{\prime \prime }{.} $84 E from the galaxy core, or a projected offset of 1.8 kpc.

PS17fem is a SN candidate associated with the galaxy SDSS J224219.29+064237.5. It is offset from the galaxy core by 0 . $ \overset{\prime \prime }{.} $57 N and 0 . $ \overset{\prime \prime }{.} $41 E or a projected offset of 2.1 kpc.

PS17eyo is a SN candidate that is not associated with any known or catalogued host.

PS17eyl is a nuclear transient, being associated with the galaxy SDSS J223902.90+062526.6. It is offset from the galaxy core by 0 . $ \overset{\prime \prime }{.} $28 S and 0 . $ \overset{\prime \prime }{.} $04 E, a projected offset of 1.1 kpc.

PS17fcd is a nuclear transient located in the galaxy SDSS J224401.85+081327.5, 0 . $ \overset{\prime \prime }{.} $37 N and 0 . $ \overset{\prime \prime }{.} $13 W from the core. This is a projected offset of 2.1 kpc. This transient was first discovered in the zP1 band on MJD 58 066.20450 at magnitude 20.53.

All Tables

Table 1.

Coordinate revisions of IceCube-160427A.

Table 2.

19 transients discovered within the final 90% probability containment of the IceCube-160427A neutrino detection.

Table 3.

Photometry of PS16cgx.

Table 4.

Spectroscopy of PS16cgx.

Table 5.

Adopted and yielded values of the probability calculations of a chance discovery for the IceCube-160427A error region.

Table 6.

IceCube HESE and EHE events and Pan-STARRS1 observations from 2016−2017.

Table A.1.

31 transients discovered within the local vicinity of the IceCube-160731A neutrino detection (MJD 57 600.07990).

Table A.2.

7 transients discovered within the local vicinity of the IceCube-161103A neutrino detection (MJD 57 695.38022).

Table A.3.

7 transients discovered within the local vicinity of the IceCube-170321A neutrino detection (MJD 57 805.31413).

Table A.4.

7 transients discovered within the local vicinity of the IceCube-171106A neutrino detection (MJD 58 063.77754).

All Figures

thumbnail Fig. 1.

Black (red) line: 50% (90%) confidence region. The confidence region likelihood ratios were determined using Monte Carlo re-simulations of the event, taking into account uncertainties in the optical properties of the deep glacial ice. The best fit position is shown as a × mark and the position of PS16cgx is indicated with a star symbol. Furthermore, the other faint transients discovered within the final 90% probability containment by Pan-STARRS1 are also shown and contextually classified based on their host or precursor information as labeled in the legend.

In the text
thumbnail Fig. 2.

Pan-STARRS1 riz-band false-colour 1′×1′ subsection of the field of PS16cgx at MJD 57 525.5. North is up, east is left.

In the text
thumbnail Fig. 3.

Panel a: spectra of PS16cgx compared to a selection of Type Ic SNe. The comparison events include the GRB associated Type Ic-BL event SN 1998bw (Patat et al. 2001); non-GRB associated Type Ic-BL event SN 2002ap (Modjaz et al. 2014); and normal Type Ic events, SN 2007gr (Hunter et al. 2009), and SN 2004aw (Taubenberger et al. 2006). The strong host galaxy emission lines of PS16cgx are labelled. Panel b: spectra of PS16cgx compared to a selection of Type Ia SNe. The comparison events include SN 2005eq, SN 2005na (Blondin et al. 2012), SN 2006py (Folatelli et al. 2013), and SN 2011fe (Pereira et al. 2013). In both panels the spectra have been corrected for line-of-sight reddening and to the rest-frame wavelength. The wavelengths of the main telluric bands in the spectra of PS16cgx are indicated with a ⊕ symbol. Comparison data were obtained via WISeREP (Yaron & Gal-Yam 2012).

In the text
thumbnail Fig. 4.

Panel a: observed light curves of PS16cgx. Zero epoch is set to the IceCube-160427A detection. A low-order polynomial light-curve fit is shown with a solid curve. Spectroscopic GMOS (30 days) and FORS2 (42 and 43 days) observations are indicated with vertical arrows. Panel b: PS16cgx r-band absolute light-curve comparison to a selection of nearby Type Ic SNe including the GRB associated Type Ic-BL events (SNe 2006aj and 1998bw); non-GRB associated Type Ic-BL events (PTF10vgv, and SNe 2016coi, 2009bb, and 2002ap); and normal Type Ic events (LSQ14efd, and SNe 2007gr and 2004aw). Panel c: selection of rapidly rising Type Ic SNe arbitrarily shifted as indicated in the figure legend. The comparison indicates an inconsistent spectrophotometric evolution of PS16cgx with Type Ic SNe which cannot be explained with extinction effects. Panel d: PS16cgx gri absolute light-curve comparison to a selection of similar Type Ia SNe. In addition, the slightly fainter r-band light curve of the canonical Type Ia SN 2011fe (Guillochon et al. 2017) is also included in the comparison. Furthermore, as a difference to the previous panels, the light curves of PS16cgx have been arbitrarily corrected for a modest host galaxy extinction of AV = 0.2 mag with RV = 1.27 derived by Burns et al. (2014) as an average value for their Type Ia sample. Assuming that the CFHT template reference image contains contamination from a Type Ia SN at mi = 24.5 mag, the final photometry point is also shown corrected for this effect with a square symbol in all panels.

In the text
thumbnail Fig. 5.

Panel a: rest-frame Mg − Mr and Panel b: Mr − Mi colours of PS16cgx compared to the evolution of Type Ia SNe (solid curves; Burns et al. 2014; Krisciunas et al. 2017) including both normal and 91T-like events, and normal Type Ic SNe (dotted curves; Stritzinger et al. 2018a,b; Taddia et al. 2018) from CSP. Colour points assuming no host galaxy extinction are shown with open symbols, and assuming an arbitrary host extinction of AV = 0.2 mag in grey symbols.

In the text
thumbnail Fig. 6.

Cumulative fraction, fd, of astrophysical neutrino sources with a single high-energy neutrino detection in IceCube, accumulated as a function of source redshift, assuming they follow the SFR (Madau & Dickinson 2014) given standard cosmology. The redshift of PS16cgx is marked with a vertical line. The maximum distance ranges where a normal Type Ic (squares) and Type Ic-BL (circles) SNe can be detected in rise are indicated assuming negligible line-of-sight extinction, ignoring K-corrections, and adopting 1σ distribution of the SN peak magnitudes, for image depths of 22.5 and 25 mag (white and grey symbols, respectively).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.