skip to main content


Title: ASAS-SN follow-up of IceCube high-energy neutrino alerts
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 neutrino sources for a range of assumed evolution models.

 
more » « less
Award ID(s):
1908952
NSF-PAR ID:
10371056
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
516
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
p. 2455-2469
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 redshift 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. 
    more » « less
  2. 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 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.

     
    more » « less
  3. ABSTRACT

    The origin of cosmic high-energy neutrinos remains largely unexplained. For high-energy neutrino alerts from IceCube, a coincidence with time-variable emission has been seen for three different types of accreting black holes: (1) a gamma-ray flare from a blazar (TXS 0506+056), (2) an optical transient following a stellar tidal disruption event (TDE; AT2019dsg), and (3) an optical outburst from an active galactic nucleus (AGN; AT2019fdr). For the latter two sources, infrared follow-up observations revealed a powerful reverberation signal due to dust heated by the flare. This discovery motivates a systematic study of neutrino emission from all supermassive black hole with similar dust echoes. Because dust reprocessing is agnostic to the origin of the outburst, our work unifies TDEs and high-amplitude flares from AGN into a population that we dub accretion flares. Besides the two known events, we uncover a third flare that is coincident with a PeV-scale neutrino (AT2019aalc). Based solely on the optical and infrared properties, we estimate a significance of 3.6σ for this association of high-energy neutrinos with three accretion flares. Our results imply that at least ∼10 per cent of the IceCube high-energy neutrino alerts could be due to accretion flares. This is surprising because the sum of the fluence of these flares is at least three orders of magnitude lower compared to the total fluence of normal AGN. It thus appears that the efficiency of high-energy neutrino production in accretion flares is increased compared to non-flaring AGN. We speculate that this can be explained by the high Eddington ratio of the flares.

     
    more » « less
  4. Several tentative associations between high-energy neutrinos and astrophysical sources have been recently reported, but a conclusive identification of these potential neutrino emitters remains challenging. We explore the use of Monte Carlo simulations of source populations to gain deeper insight into the physical implications of proposed individual source–neutrino associations. In particular, we focus on the IC170922A–TXS 0506+056 observation. Assuming a null model, we find a 7.6% chance of mistakenly identifying coincidences between γ -ray flares from blazars and neutrino alerts in 10-year surveys. We confirm that a blazar–neutrino connection based on the γ -ray flux is required to find a low chance coincidence probability and, therefore, a significant IC170922A–TXS 0506+056 association. We then assume this blazar–neutrino connection for the whole population and find that the ratio of neutrino to γ -ray fluxes must be ≲10 −2 in order not to overproduce the total number of neutrino alerts seen by IceCube. For the IC170922A–TXS 0506+056 association to make sense, we must either accept this low flux ratio or suppose that only some rare sub-population of blazars is capable of high-energy neutrino production. For example, if we consider neutrino production only in blazar flares, we expect the flux ratio of between 10 −3 and 10 −1 to be consistent with a single coincident observation of a neutrino alert and flaring γ -ray blazar. These constraints should be interpreted in the context of the likelihood models used to find the IC170922A–TXS 0506+056 association, which assumes a fixed power-law neutrino spectrum of E −2.13 for all blazars. 
    more » « less
  5. ABSTRACT Using blazar light curves from the optical All-Sky Automated Survey for Supernovae (ASAS-SN) and the γ-ray Fermi-LAT telescope, we performed the most extensive statistical correlation study between both bands, using a sample of 1180 blazars. This is almost an order of magnitude larger than other recent studies. Blazars represent more than 98 per cent of the AGNs detected by Fermi-LAT and are the brightest γ-ray sources in the extragalactic sky. They are essential for studying the physical properties of astrophysical jets from central black holes. However, their γ-ray flare mechanism is not fully understood. Multiwavelength correlations help constrain the dominant mechanisms of blazar variability. We search for temporal relationships between optical and γ-ray bands. Using a Bayesian Block Decomposition, we detect 1414 optical and 510 γ-ray flares, we find a strong correlation between both bands. Among all the flares, we find 321 correlated flares from 133 blazars, and derive an average rest-frame time delay of only 1.1$_{-8.5}^{+7.1}$ d, with no difference between the flat-spectrum radio quasars, BL Lacertae-like objects or low, intermediate, and high-synchrotron peaked blazar classes. Our time-delay limit rules out the hadronic proton-synchrotron model as the driver for non-orphan flares and suggests a leptonic single-zone model. Limiting our search to well-defined light curves and removing 976 potential but unclear ‘orphan’ flares, we find 191 (13 per cent) and 115 (22 per cent) clear ‘orphan’ optical and γ-ray flares. The presence of ‘orphan’ flares in both bands challenges the standard one-zone blazar flare leptonic model and suggests multizone synchrotron sites or a hadronic model for some blazars. 
    more » « less