skip to main content

Title: Neutrino follow-up with the Zwicky transient facility: results from the first 24 campaigns
ABSTRACT

The Zwicky Transient Facility (ZTF) performs a systematic neutrino follow-up programme, searching for optical counterparts to high-energy neutrinos with dedicated Target-of-Opportunity (ToO) observations. Since first light in March 2018, ZTF has taken prompt observations for 24 high-quality neutrino alerts from the IceCube Neutrino Observatory, with a median latency of 12.2 h from initial neutrino detection. From two of these campaigns, we have already reported tidal disruption event (TDE) AT 2019dsg and likely TDE AT 2019fdr as probable counterparts, suggesting that TDEs contribute >7.8 per cent of the astrophysical neutrino flux. We here present the full results of our programme through to December 2021. No additional candidate neutrino sources were identified by our programme, allowing us to place the first constraints on the underlying optical luminosity function of astrophysical neutrino sources. Transients with optical absolutes magnitudes brighter that −21 can contribute no more than 87 per cent of the total, while transients brighter than −22 can contribute no more than 58 per cent of the total, neglecting the effect of extinction and assuming they follow the star formation rate. These are the first observational constraints on the neutrino emission of bright populations such as superluminous supernovae. None of the neutrinos were coincident with bright more » optical AGN flares comparable to that observed for TXS 0506+056/IC170922A, with such optical blazar flares producing no more than 26 per cent of the total neutrino flux. We highlight the outlook for electromagnetic neutrino follow-up programmes, including the expected potential for the Rubin Observatory.

« less
Authors:
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; « less
Award ID(s):
2034437
Publication Date:
NSF-PAR ID:
10405416
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
521
Issue:
4
Page Range or eLocation-ID:
p. 5046-5063
ISSN:
0035-8711
Publisher:
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    We report on the search for optical counterparts to IceCube neutrino alerts released between 2016 April and 2021 August with the All-Sky Automated Survey for SuperNovae (ASAS-SN). Despite the discovery of a diffuse astrophysical high-energy neutrino flux in 2013, the source of those neutrinos remains largely unknown. Since 2016, IceCube has published likely astrophysical neutrinos as public real-time alerts. Through a combination of normal survey and triggered target-of-opportunity observations, ASAS-SN obtained images within 1 h of the neutrino detection for 20 per cent (11) of all observable IceCube alerts and within one day for another 57 per cent (32). For all observable alerts, we obtained images within at least two weeks from the neutrino alert. ASAS-SN provides the only optical follow-up for about 17 per cent of IceCube’s neutrino alerts. We recover the two previously claimed counterparts to neutrino alerts, the flaring-blazar TXS 0506 + 056 and the tidal disruption event AT2019dsg. We investigate the light curves of previously detected transients in the alert footprints, but do not identify any further candidate neutrino sources. We also analysed the optical light curves of Fermi 4FGL sources coincident with high-energy neutrino alerts, but do not identify any contemporaneous flaring activity. Finally, we derive constraints on the luminosity functions of neutrinomore »sources for a range of assumed evolution models.

    « less
  2. In order to identify the sources of the observed diffuse high-energy neutrino flux, it is crucial to discover their electromagnetic counterparts. To increase the sensitivity of detecting counterparts of transient or variable sources by telescopes with a limited field of view, IceCube began releasing alerts for single high-energy ( E ν  >  60 TeV) neutrino detections with sky localisation regions of order 1° radius in 2016. We used Pan-STARRS1 to follow-up five of these alerts during 2016–2017 to search for any optical transients that may be related to the neutrinos. Typically 10–20 faint ( m i P1  ≲ 22.5 mag) extragalactic transients are found within the Pan-STARRS1 footprints and are generally consistent with being unrelated field supernovae (SNe) and AGN. We looked for unusual properties of the detected transients, such as temporal coincidence of explosion epoch with the IceCube timestamp, or other peculiar light curve and physical properties. We found only one transient that had properties worthy of a specific follow-up. In the Pan-STARRS1 imaging for IceCube-160427A (probability to be of astrophysical origin of ∼50%), we found a SN PS16cgx, located at 10.0′ from the nominal IceCube direction. Spectroscopic observations of PS16cgx showed that it was an H-poor SN at redshiftmore »z  = 0.2895 ± 0.0001. The spectra and light curve resemble some high-energy Type Ic SNe, raising the possibility of a jet driven SN with an explosion epoch temporally coincident with the neutrino detection. However, distinguishing Type Ia and Type Ic SNe at this redshift is notoriously difficult. Based on all available data we conclude that the transient is more likely to be a Type Ia with relatively weak Si  II absorption and a fairly normal rest-frame r -band light curve. If, as predicted, there is no high-energy neutrino emission from Type Ia SNe, then PS16cgx must be a random coincidence, and unrelated to the IceCube-160427A. We find no other plausible optical transient for any of the five IceCube events observed down to a 5 σ limiting magnitude of m i P1  ≈ 22 mag, between 1 day and 25 days after detection.« less
  3. ABSTRACT

    Searches for optical transients are usually performed with a cadence of days to weeks, optimized for supernova discovery. The optical fast transient sky is still largely unexplored, with only a few surveys to date having placed meaningful constraints on the detection of extragalactic transients evolving at sub-hour time-scales. Here, we present the results of deep searches for dim, minute-time-scale extragalactic fast transients using the Dark Energy Camera, a core facility of our all-wavelength and all-messenger Deeper, Wider, Faster programme. We used continuous 20 s exposures to systematically probe time-scales down to 1.17 min at magnitude limits g > 23 (AB), detecting hundreds of transient and variable sources. Nine candidates passed our strict criteria on duration and non-stellarity, all of which could be classified as flare stars based on deep multiband imaging. Searches for fast radio burst and gamma-ray counterparts during simultaneous multifacility observations yielded no counterparts to the optical transients. Also, no long-term variability was detected with pre-imaging and follow-up observations using the SkyMapper optical telescope. We place upper limits for minute-time-scale fast optical transient rates for a range of depths and time-scales. Finally, we demonstrate that optical g-band light-curve behaviour alone cannot discriminate between confirmed extragalactic fast transients such asmore »prompt GRB flashes and Galactic stellar flares.

    « less
  4. We report on SRG/eROSITA, ZTF, ASAS-SN, Las Cumbres, NEOWISE-R, and Swift XRT/UVOT observations of the unique ongoing event AT 2019avd, located in the nucleus of a previously inactive galaxy at z = 0.029. eROSITA first observed AT 2019avd on 2020-04-28 during its first all sky survey, when it was detected as an ultra-soft X-ray source ( kT ~ 85 eV) that was ≳90 times brighter in the 0.2−2 keV band than a previous 3 σ upper flux detection limit (with no archival X-ray detection at this position). The ZTF optical light curve in the ~450 days preceding the eROSITA detection is double peaked, and the eROSITA detection coincides with the rise of the second peak. Follow-up optical spectroscopy shows the emergence of a Bowen fluorescence feature and high-ionisation coronal lines ([Fe  X ] 6375 Å, [Fe  XIV ] 5303 Å), along with persistent broad Balmer emission lines ( FWHM ~ 1400 km s −1 ). Whilst the X-ray properties make AT 2019avd a promising tidal disruption event (TDE) candidate, the optical properties are atypical for optically selected TDEs. We discuss potential alternative origins that could explain the observed properties of AT 2019avd, such as a stellar binary TDE candidate, ormore »a TDE involving a super massive black hole binary.« less
  5. ABSTRACT

    Active galactic nuclei (AGNs) make up about 35 per cent of the more than 250 sources detected in very high-energy (VHE) gamma rays to date with the imaging atmospheric Cherenkov telescopes. Apart from four nearby radio galaxies and two AGNs of unknown type, all known VHE AGNs are blazars. Knowledge of the cosmological redshift of gamma-ray blazars is key to enabling the study of their intrinsic emission properties, as the interaction between gamma rays and the extragalactic background light (EBL) results in a spectral softening. Therefore, the redshift determination exercise is crucial to indirectly placing tight constraints on the EBL density, and to studying blazar population evolution across cosmic time. Due to the powerful relativistic jets in blazars, most of their host galaxies’ spectral features are outshined, and dedicated high signal-to-noise (S/N) spectroscopic observations are required. Deep medium- to high-resolution spectroscopy of 33 gamma-ray blazar optical counterparts was performed with the European Southern Observatory, New Technology Telescope, Keck II telescope, Shane 3-metre telescope, and the Southern African Large Telescope. From the sample, spectra from 25 objects display spectral features or are featureless and have high S/N. The other eight objects have low-quality featureless spectra. We systematically searched for absorptionmore »and emission features and estimated, when possible, the fractional host galaxy flux in the measured total flux. Our measurements yielded 14 firm spectroscopic redshifts, ranging from 0.0838 to 0.8125, one tentative redshift, and two lower limits: one at $z > 0.382$ and the other at z > 0.629.

    « less