All-Flavor Time-Dependent Search For Transient Neutrino Sources
All-Flavor Time-Dependent Search For Transient Neutrino Sources
Neutrino Sources
arXiv:2507.08775v1 [astro-ph.HE] 11 Jul 2025
Transient sources are among the preferred candidates for the sources of high-energy neutrino
emission. Intriguing examples so far include blazar flares and tidal disruption events coinci-
dent with IceCube neutrinos. Here, we report the first all-flavor, all-sky time-dependent search
for neutrino sources by combining IceCube throughgoing tracks, starting tracks and cascades.
Throughgoing tracks provide the best sensitivity in the Northern Sky, while cascades have worse
angular resolution but yield better sensitivity in the Southern Sky than tracks. The relatively new
starting tracks sample has reduced contamination from atmospheric muons. This analysis takes
advantage of the strengths of each of the datasets, combining them for increased statistics and
obtaining the best accessible all-sky sensitivity for transient searches. In this search, we look for
unbound 𝐸 −𝛾 power-law sources, as well as 𝐸 −2 sources with low and high-energy exponential
cutoffs, optimizing the sensitivity for the duration of the flares.
∗ Presenter
© Copyright owned by the author(s) under the terms of the Creative Commons
Attribution-NonCommercial-NoDerivatives 4.0 International License (CC BY-NC-ND 4.0). https://pos.sissa.it/
All-flavor Time-dependent Search for Transient Neutrino Sources
1. Introduction
Transient neutrino sources, such as blazar flares, tidal disruption events, and gamma-ray bursts,
are among the primary candidates for the origin of high-energy neutrino emission [1]. In 2017,
IceCube detected a 290 TeV neutrino spatially coincident with the blazar TXS 0506+056 [2].
Time-integrated searches with 10 years of neutrino tracks showed the transient TXS 0506+056 as
the second most significant source in the Northern Sky [3]. However, the neutrino signal is not
necessarily accompanied by an electromagnetic counterpart, leaving neutrinos as the only way to
observe these sources. These would appear in neutrino data as a clustering of events in the time
domain.
Time-dependent searches provide an advantage over the time-integrated ones, for flares that
last for Δ𝑇 ≲ 100 days, as it reduces the atmospheric neutrino background, which is relatively
uniform in time. Previously, the IceCube collaboration has performed a transient search [4] and
multi-flare search [5], but did not find any new significant excesses.
In this study, we aim for improving upon previous searches by utilizing all channels of obser-
vations to obtain the best accessible sensitivity for IceCube.
Furthermore, neutrino point source searches have usually looked for generic power-law fluxes
−𝛾
𝐸 . However, the expected neutrino flux generally has features beyond a simple power-law. The
IceCube collaboration has used alternative parametrizations for the diffuse neutrino flux (see e.g.,
[6]). More recently, point source searches have also used more sophisticated signal energy spectra,
such as Seyfert flux models to search for neutrino emission from Seyfert galaxies [7]. Here, we
implement a power-law with low and high-energy cutoffs in the search that could capture the
predicted features in prominent transient sources of neutrinos.
2. Method
For the analysis we will use the likelihood (LLH) function L [8]
∑︁ h 𝑛𝑠 𝑛𝑠 i
ln L = ln S𝑖 + 1 − B𝑖 , (1)
𝑖
𝑁 𝑁
where S𝑖 and B𝑖 are the signal and background probability density functions (pdfs), respectively,
for each event 𝑖. The signal pdf can be factorized into a spatial pdf term, which describes the spatial
clustering of the events around the source, an energy pdf term, describing their energy distribution,
and a temporal pdf term, describing the arrival time distribution of the events.
This time-dependent analysis assumes that the temporal pdf follows a Gaussian distribution,
centered at 𝑇0 and a standard deviation 𝜎𝑇 . The spatial pdfs are two-dimensional Gaussians, where
the reconstructed direction has a standard deviation 𝜎𝑖 . The NT sample was later enhanced by using
a Kernel Density Estimate (KDE) so that the spectral shape of the flux is also taken into account,
making the spatial pdf more realistic (see e.g., supplemental material in [9] for a description of the
KDE method).
The energy pdf is the probability that an event is observed with a given energy at a fixed
declination 𝛿 and for a given spectral hypothesis. In the case of 𝐸 −𝛾 power-law sources, this
introduces an additional parameter, the power-law index 𝛾.
2
All-flavor Time-dependent Search for Transient Neutrino Sources
δ = 0 ◦ , E −2 source 0.10
δ = 0 ◦ , E −2 source
DNN NT + DNN + ESTES NT KDE 90% sensitivity
∆T E 2 dN/dE| 1TeV [GeV cm −2 ]
10 1 0.04
Combined/NT KDE
1.0 ICECUBE PRELIMINARY 1.0 ICECUBE PRELIMINARY
Combined/NT
0.9 0.9
0.8 3 2 1 0.8
10 10 10 100 101 102 10 3
10 2
10 1
100 101 102
σT [days] σT [days]
Figure 1: Top left panel: Per-flavor 90% CL sensitivity (solid lines) and 3𝜎 evidence (dashed lines) fluences
for an 𝐸 −2 source at declination 𝛿 = 0◦ , for different flare half-widths 𝜎𝑇 and 𝐸 −𝛾 signal flux hypothesis.
The blue, red, orange and black lines correspond to DNNCascade, NT, ESTES and combined datasets,
respectively. Here, the NT sample assumes Gaussian spatial pdfs. Top right panel: Same as top left panel,
but the NT dataset is using KDEs in its spatial and energy pdfs. Bottom left panel: Combined-to-NT fluence
ratios. Bottom right panel: Combined-to-NT KDE fluence ratios
ˆ 𝑇ˆ0 , 𝜎
𝜎ˆ 𝑇 L ( 𝑛ˆ 𝑠 , 𝛾, ˆ 𝑇)
TS = 2 ln , (2)
𝑇live L (𝑛𝑠 = 0)
3
All-flavor Time-dependent Search for Transient Neutrino Sources
102
δ = − 30 ◦ , σT = 1 day
DNN 90% sensitivity
∆T E 2 dN/dE [GeV cm −2 ] ESTES
DNN+ESTES
3σ evidence potential
101
ICECUBE PRELIMINARY
100
Figure 2: Per-flavor 90% CL sensitivity (solid lines) and 3𝜎 evidence potential (dashed lines) fluences to a
neutrino flare of half-width 𝜎𝑇 . The blue, orange and black lines correspond to DNNCascade, ESTES and
combined datasets, respectively.
3. Analysis performance
4
All-flavor Time-dependent Search for Transient Neutrino Sources
where 𝐸 𝐿 and 𝐸 𝐻 are the low and high-energy cutoffs, respectively. We will also refer to this flux
as a two-sided cutoff. The inclusion of 𝐸 𝐻 is motivated by the presence of a maximum cosmic ray
energy, while 𝐸 𝐿 accounts for the minimum proton energy required to initiate 𝑝𝛾 or 𝑝 𝑝 interactions
in the source and lead to neutrino production. The use of a more representative signal hypothesis
allows for a more accurate sensitivity. Based on the differential sensitivities, having a two-sided
cutoff flux that peaks around 100 TeV could potentially be missed in a power-law fit.
With this new flux hypothesis, the power-law index 𝛾 = 2 is fixed when calculating TS. We
restrict 𝐸 𝐿 ∈ [100 GeV, 10 TeV] and 𝐸 𝐻 ∈ [10 TeV, 100 PeV]. When fitting for two-sided cutoffs,
the fit parameters are 𝑛𝑠 , 𝐸 𝐿 , 𝐸 𝐻 , 𝑇0 and 𝜎𝑇 .
We proceed to compare the performance between flux hypotheses. We assume that the true
signal flux is given by Eq. (3), with 𝐸 𝐿 = 1 TeV , 𝐸 𝐻 = 100 TeV. We then compare the 3𝜎 evidence
potential for this injected flux under a two-sided cutoff and a power-law hypothesis separately. First,
we identify the 3𝜎 evidence potential for the two-sided cutoff hypothesis. We then inject this 3𝜎
two-sided cutoff flux, as the true signal flux, and perform a power-law fit instead. Note that by
swapping to a different signal hypothesis, the signal pdfs going into equation (2) have changed, and
thus the TS adopts a different value even if the events in a given trial are the same. Hence, we have to
create new signal and background TS distributions. In a power-law fit, the trials for the injected 3𝜎
two-sided cutoff flux will no longer generate a TS which is 50% of the time above the 3𝜎 threshold
of the new background TS distribution (i.e. it is not a 3𝜎 flux under this new hypothesis). Under
the power-law hypothesis, we can compute the new number of sigma associated with this two-sided
cutoff flux, which we call the recovered significance. In this sense, we say that the 3𝜎 two-sided
cutoff flux has been recovered with a different significance.
We show the significance recovery in Figure 3. For this example, the left (right) panel uses
the DNNCascade (NT) sample only. The NT sample in the right panel does not use KDEs. We
find that a two-sided cutoff flux is not recovered at the 3𝜎 level in power-law searches. In fact, we
only recover it at 1𝜎 for the shortest time windows. As the flare duration increases, the recovery
improves. We found that the TS distribution is not significantly different between signal hypotheses.
5
All-flavor Time-dependent Search for Transient Neutrino Sources
Figure 3: Recovered significance in power-law fits after injecting the 3𝜎 evidence flux for two-sided cutoffs,
for different flare half-widths 𝜎𝑇 . Left (right) panel corresponds to the DNNCascade (NT) sample.
Hence, the loss in significance is mostly tied to the reconstruction of the signal pdf parameters. The
power-law hypothesis tends to fit for 𝛾 > 2, which causes an overfitting of the signal parameter
𝑛𝑠 . In the case of flares with 𝜎𝑇 ≲ 1 day, the injection of a 3𝜎 two-sided cutoff flux amounts to
injecting ≈ 3 − 4 signal events. For these values of 𝜎𝑇 , the injected number of events is too low for
the LLH to fully take advantage of the background-free regime and reduce the 𝑛𝑠 bias.
4. Conclusions
This is the first multi-flavor time-dependent all-sky point source search. Here we covered the
difference in performance between datasets and its improvement upon combining them. We report
a ≈ 15% improvement in sensitivity and evidence potential for 𝐸 −2 sources at the horizon, when
using the combined dataset. We also find that a power-law fit in our LLH analysis might miss a
two-sided cutoff signal, particularly for flares with 𝜎𝑇 ≲ 1 day.
References
[1] K. Murase and I. Bartos, Ann. Rev. Nucl. Part. Sci. 69 (2019) 477–506.
[2] IceCube Collaboration, M. G. Aartsen et al., Science 361 no. 6398, (2018) 147–151.
[3] IceCube Collaboration, M. G. Aartsen et al., Phys. Rev. Lett. 124 no. 5, (2020) 051103.
[4] IceCube Collaboration, R. Abbasi et al., Astrophys. J. 967 no. 1, (2024) 48.
[5] IceCube Collaboration, R. Abbasi et al., Astrophys. J. Lett. 920 no. 2, (2021) L45.
6
All-flavor Time-dependent Search for Transient Neutrino Sources
[8] J. Braun, M. Baker, J. Dumm, C. Finley, A. Karle, and T. Montaruli, Astropart. Phys. 33
(2010) 175–181.
[9] IceCube Collaboration, R. Abbasi et al., Science 378 no. 6619, (2022) 538–543.
7
All-flavor Time-dependent Search for Transient Neutrino Sources
8
All-flavor Time-dependent Search for Transient Neutrino Sources
USA
40 Institute of Physics, University of Mainz, Staudinger Weg 7, D-55099 Mainz, Germany
41 Department of Physics, Marquette University, Milwaukee, WI 53201, USA
42 Institut für Kernphysik, Universität Münster, D-48149 Münster, Germany
43 Bartol Research Institute and Dept. of Physics and Astronomy, University of Delaware, Newark, DE 19716, USA
44 Dept. of Physics, Yale University, New Haven, CT 06520, USA
45 Columbia Astrophysics and Nevis Laboratories, Columbia University, New York, NY 10027, USA
46 Dept. of Physics, University of Oxford, Parks Road, Oxford OX1 3PU, United Kingdom
47 Dipartimento di Fisica e Astronomia Galileo Galilei, Università Degli Studi di Padova, I-35122 Padova PD, Italy
48 Dept. of Physics, Drexel University, 3141 Chestnut Street, Philadelphia, PA 19104, USA
49 Physics Department, South Dakota School of Mines and Technology, Rapid City, SD 57701, USA
50 Dept. of Physics, University of Wisconsin, River Falls, WI 54022, USA
51 Dept. of Physics and Astronomy, University of Rochester, Rochester, NY 14627, USA
52 Department of Physics and Astronomy, University of Utah, Salt Lake City, UT 84112, USA
53 Dept. of Physics, Chung-Ang University, Seoul 06974, Republic of Korea
54 Oskar Klein Centre and Dept. of Physics, Stockholm University, SE-10691 Stockholm, Sweden
55 Dept. of Physics and Astronomy, Stony Brook University, Stony Brook, NY 11794-3800, USA
56 Dept. of Physics, Sungkyunkwan University, Suwon 16419, Republic of Korea
57 Institute of Physics, Academia Sinica, Taipei, 11529, Taiwan
58 Dept. of Physics and Astronomy, University of Alabama, Tuscaloosa, AL 35487, USA
59 Dept. of Astronomy and Astrophysics, Pennsylvania State University, University Park, PA 16802, USA
60 Dept. of Physics, Pennsylvania State University, University Park, PA 16802, USA
61 Dept. of Physics and Astronomy, Uppsala University, Box 516, SE-75120 Uppsala, Sweden
62 Dept. of Physics, University of Wuppertal, D-42119 Wuppertal, Germany
63 Deutsches Elektronen-Synchrotron DESY, Platanenallee 6, D-15738 Zeuthen, Germany
a also at Institute of Physics, Sachivalaya Marg, Sainik School Post, Bhubaneswar 751005, India
b also at Department of Space, Earth and Environment, Chalmers University of Technology, 412 96 Gothenburg, Sweden
c also at INFN Padova, I-35131 Padova, Italy
d also at Earthquake Research Institute, University of Tokyo, Bunkyo, Tokyo 113-0032, Japan
e now at INFN Padova, I-35131 Padova, Italy
9
All-flavor Time-dependent Search for Transient Neutrino Sources
Acknowledgments
The authors gratefully acknowledge the support from the following agencies and institutions: USA – U.S. National Science Foundation-
Office of Polar Programs, U.S. National Science Foundation-Physics Division, U.S. National Science Foundation-EPSCoR, U.S. National
Science Foundation-Office of Advanced Cyberinfrastructure, Wisconsin Alumni Research Foundation, Center for High Throughput
Computing (CHTC) at the University of Wisconsin–Madison, Open Science Grid (OSG), Partnership to Advance Throughput Computing
(PATh), Advanced Cyberinfrastructure Coordination Ecosystem: Services & Support (ACCESS), Frontera and Ranch computing
project at the Texas Advanced Computing Center, 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, Astroparticle physics computational facility at Marquette University, NVIDIA Corporation, and Google Cloud Platform;
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 Forschung, Technologie und Raumfahrt (BMFTR), 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; European Union – EGI Advanced Computing for research; Australia – Australian Research Council;
Canada – Natural Sciences and Engineering Research Council of Canada, Calcul Québec, Compute Ontario, Canada Foundation for
Innovation, WestGrid, and Digital Research Alliance of Canada; Denmark – Villum Fonden, Carlsberg Foundation, and European
Commission; 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).
10