Pith. sign in

REVIEW 3 major objections 5 minor 2 cited by

Neutrinos from Primordial Black Hole Evaporation

T0 review · 3 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A neutrino telescope search homes in on the final 1,000 seconds of evaporating primordial black holes, using one year of track-like data to reach a local burst-rate-density sensitivity of…

desk verdict A legitimate proof-of-concept for an IceCube PBH-neutrino search, but the headline burst-rate sensitivity is computed from a dimensionally inconsistent formula and needs correction. read the letter →

arxiv 1908.05403 v1 pith:LFUC3YYA submitted 2019-08-15 astro-ph.HE

Pranav Dave , Ignacio Taboada (for the IceCube Collaboration) This is my paper · ORCID
classification astro-ph.HE
keywords primordialblackholesholeevaporationneutrinoburststime-dependentpointsourcesearchunbinnedlikelihoodtrack-likeeventsTeVastronomydarkmattercandidates
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Primordial black holes near $10^{15}$ g should be completing their evaporation around the present time, and their final moments should produce a burst of high-energy neutrinos. This proceedings paper adapts a neutrino telescope's standard time-dependent point-source flare search so that it can look for those bursts. It claims that with one year of track-like data the search is most sensitive to the last $10^3$ seconds of a black hole's life, corresponding to emission temperatures around $0.77$ TeV, and reaches a local burst-rate-density sensitivity of $3.8\times10^{6}\,\mathrm{pc}^{-3}\,\mathrm{yr}^{-1}$. Because the simulation software can only model time-independent spectra, the intrinsically time-dependent evaporation spectrum is approximated by three time bins; the paper states that this is known to be not optimal. This is the first neutrino-based search for evaporating primordial black holes and can be extended to ten years of data.

What carries the argument

The machinery is the three-bin time-integrated neutrino fluence, built from the instantaneous emission rate of the evaporation spectrum integrated over 0–10 s, 10–100 s, and 100–1000 s to evaporation at a reference distance. Each bin is treated as a time-independent broken power law, with spectral breaks set by the black-hole temperature at the bin boundary; subtracting fluences between times gives the fluence of any interval. This feeds an unbinned time-dependent likelihood with a signal PDF composed of directional, energy, and Gaussian temporal terms, and a test statistic whose 90% threshold defines the sensitivity distance. The three-bin approximation is the load-bearing adaptation, since the detector simulation software accepts only time-independent spectra.

What would settle it

Inject simulated bursts generated with a continuously evolving evaporation spectrum (or with many time bins) into the same one-year scrambled data and compare the 90% sensitivity distance and 5σ discovery distance with the three-bin result; if the distances differ by more than the statistical uncertainty of the pseudo-experiments, the reported $3.8\times10^{6}\,\mathrm{pc}^{-3}\,\mathrm{yr}^{-1}$ sensitivity and the claim that the last $10^3$ s dominate would be refuted.

Watch

Extended reading notes

Core claim

The paper's central claim is that a standard unbinned time-dependent flare analysis, in which signal events are modeled by a Gaussian time profile, can be repurposed to search for evaporating primordial black holes provided the time-dependent emission spectrum is binned into three constant-spectrum time intervals: 0–10 s, 10–100 s, and 100–1000 s before evaporation. The paper finds that the final $10^3$ s dominate the sensitivity, with the black-hole temperature around $0.77$ TeV, and that one year of track-like events yields a 90% sensitivity distance that converts to a local burst-rate-density sensitivity of $3.8\times10^{6}\,\mathrm{pc}^{-3}\,\mathrm{yr}^{-1}$. This is less sensitive than the previously published gamma-ray limits, but it is the first neutrino search at these energies and the paper argues it has a clear path to improvement with ten years of data.

Load-bearing premise

The reported sensitivity assumes that the intrinsically time-dependent evaporation spectrum is faithfully represented by three constant-spectrum time bins and a Gaussian flare time profile; the paper itself notes that the three-bin approximation is not optimal.

Editorial extensions

If this is right

  • With the full ten-year sample, the same three-bin analysis is expected to push the local burst-rate-density sensitivity well below the one-year value of $3.8\times10^{6}\,\mathrm{pc}^{-3}\,\mathrm{yr}^{-1}$.
  • The most sensitive window is the final $10^3$ seconds of evaporation, at black-hole temperatures around $0.77$ TeV, which is a different phase of the burst than the windows probed by gamma-ray searches.
  • A detected burst would yield a neutrino-based measurement of the local rate of evaporating primordial black holes, the first at these energies.
  • The analysis shows that an existing all-sky flare search can be reused for a new source class with only a spectral-model modification, which keeps a decade-long search computationally tractable.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • My inference: switching from three time bins to a fine time grid would better capture the steeply rising spectrum and could either improve the discovery potential or shift the optimal search window; this is a simulation-driven extension the paper does not carry out.
  • My inference: the burst-rate-density limit can be translated into an upper bound on the present-day abundance of primordial black holes in the $10^{15}$ g range once a mass-to-rate mapping is assumed; the paper stops at the rate density.
  • My inference: since neutrinos traverse material that stops gamma rays, applying the same search near the Galactic plane could probe a region where gamma-ray burst searches are blind, provided atmospheric-background modeling stays under control.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper presents a proof-of-concept search for neutrino bursts from evaporating primordial black holes (PBHs) with initial masses around 10^15 g, using one year of IceCube track-like events. The Hawking spectrum is time-dependent, so the authors approximate the neutrino fluence over the last 1000 s of the PBH lifetime by three time bins with constant broken power-law spectra. They inject these signals into scrambled background data and use IceCube's unbinned time-dependent point-source likelihood to compute the number of signal events n_sens required for a 90% sensitivity and 5σ discovery potential as a function of burst duration τ. From n_sens they derive a maximum detectable distance d_sens and convert this to a local burst rate density, quoting a sensitivity of 3.8×10^6 pc^-3 yr^-1 at τ ≈ 10^3 s. They compare this with gamma-ray limits and propose that a future analysis with 10 years of IceCube data could improve the sensitivity.

Significance. The paper is significant as the first attempt to search for evaporating PBHs in neutrinos, extending the search to the highest energies and complementing gamma-ray instruments. The adaptation of the standard IceCube time-dependent search is described transparently, and the use of a broken power-law fluence approximation makes the PBH signal model computationally tractable. The qualitative finding that IceCube is most sensitive to the final ~10^3 seconds of evaporation is physically plausible. However, the central numerical result, the local burst rate density, is based on a dimensionally inconsistent conversion in Sec. 3 and is not currently supported; the single-declination treatment and the unquantified three-bin approximation add further uncertainty. If the conversion is corrected and the systematic effects quantified, the method would provide a valuable benchmark for future PBH searches, but the present manuscript needs revision before its headline number can be accepted.

major comments (3)
  1. [Sec. 3] The unnumbered equation after the distance conversion, ρ̇ = 3/(4π) n_sens / d_sens^3, is dimensionally inconsistent with the claimed units of pc^-3 yr^-1: the right-hand side has units pc^-3 (n_sens is dimensionless and d_sens is in pc) and no observation time appears. For a uniform, Poisson-distributed burst population with rate density ρ̇, the expected number of detectable bursts within a sphere of radius d_sens during an exposure T is N = (4π/3) d_sens^3 ρ̇ T, so the correct inversion is ρ̇ = N_th / [(4π/3) d_sens^3 T], with N_th ≈ 2.3 for a 90% confidence upper limit and T = 1 yr. Using n_sens, the number of signal events contributed by a single burst, in place of N_th is not a valid conversion. Because the abstract's headline value 3.8×10^6 pc^-3 yr^-1 is derived from this equation, the central quantitative claim is currently unsupported and must be recalculated.
  2. [Secs. 2.2 and 3] The sensitivity (n_sens, d_sens) is computed for a single declination, δ = 16°, but the paper then quotes a 'local density rate sensitivity' assuming a uniform PBH distribution. For an all-sky rate density, the conversion must be averaged over declination because the effective area and atmospheric neutrino background vary strongly with δ. The manuscript does not describe such a sky-average, so the quoted value 3.8×10^6 pc^-3 yr^-1 is at best a single-declination sensitivity and should not be presented as an all-sky local rate density without further justification.
  3. [Secs. 2.1 and 2.2] The three-bin approximation of the Hawking time dependence is explicitly acknowledged to be non-optimal, but the paper does not quantify how this discretization biases the computed sensitivity. In addition, the signal is injected using the three-bin time structure while the likelihood fits a Gaussian time profile; the claim that the Gaussian hypothesis is 'effective' is stated without any injection-recovery efficiency or bias test. The manuscript should provide a quantitative test (e.g., the fraction of injected PBH bursts correctly identified, or a comparison of n_sens obtained from the three-bin versus a finer time binning) so that the systematic uncertainty of the central sensitivity number is known.
minor comments (5)
  1. [Sec. 2.2] The sentence 'The test-statistic defined this way would follow a χ2 distribution where the degrees of freedom reflect the number of fit-parameters' is not correct in general for a likelihood ratio with a parameter at its physical boundary (n_s = 0); since the analysis uses scrambled-data and injection distributions, this theoretical claim should be softened or removed.
  2. [Fig. 3b and Sec. 3] The y-axis label in Fig. 3b gives units of yr^-1 pc^-3, but the rate-density equation in Sec. 3 has units pc^-3; the two should be made consistent.
  3. [Abstract and Discussion vs. Fig. 3b] The quoted sensitivity 3.8×10^6 pc^-3 yr^-1 lies below the lower limit (10^7 yr^-1 pc^-3) of the y-axis in Fig. 3b; the discrepancy between the quoted number and the plotted curve should be resolved.
  4. [Sec. 2.1] The broken power-law approximation for the time-integrated fluence is asserted without a reference or a validation plot; a comparison with a full numerical integration of Eq. (2.1) would support the approximation.
  5. [Sec. 1 and reference list] The citation 'Halzen and Zas [10]' refers to a three-authored paper (Halzen, Keszthelyi, and Zas); the in-text citation should list all authors or use 'et al.'

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the PBH sensitivity is computed from an external Hawking-radiation signal model and Monte Carlo injection, not from a fitted parameter or a self-citation chain.

full rationale

The paper's derivation chain is self-contained. The signal model starts from Hawking's formula (Eq. 2.1) and the time-integrated fluence is computed from that physics, not from IceCube data. The sensitivity is obtained by injecting simulated PBH-like signals into scrambled IceCube data, fitting the unbinned likelihood (Eq. 2.2), and extracting the required event count nsens for a 90% median upper-limit criterion. The distance scale dsens is then converted from nsens via nref, which depends only on the fluence and the IceCube effective area. The central numerical claim (3.8e6 pc^-3 yr^-1) therefore does not reduce to a fit of the target quantity or to a previous IceCube result. Citations to the standard time-dependent point-source search [14] and to IceCube software [12,13] are methodological references, not load-bearing self-references that substitute for the present calculation. The paper does acknowledge the 3-bin approximation in Sec. 2.1, but this is a modeling limitation, not circularity. If there is a weakness, it is the rate-density conversion in Sec. 3, which is dimensionally suspect (the right-hand side of rho_dot = 3 nsens/(4 pi dsens^3) has units of pc^-3 rather than pc^-3 yr^-1 and has no explicit exposure time), but that is a units/correctness concern, not an instance of a prediction being equivalent to its input by construction. Under the given hard rules, a correctness flaw is not circularity, so the circularity score is 0.

Assumptions & free parameters 1 free parameters · 6 assumptions · 0 invented entities

The central sensitivity claim rests on standard Hawking radiation physics and a series of modeling choices. The free parameter is the hand-chosen time binning; the axioms are mostly domain assumptions from prior literature, with two analysis-specific approximations (broken power-law binning and Gaussian time profile) that are acknowledged or asserted without quantitative support.

free parameters (1)
  • Time bin edges for signal simulation = 0-10 s, 10-100 s, 100-1000 s
    The continuous Hawking spectrum during the last 1000 s is approximated by three constant-spectrum bins because IceCube simulation software requires time-independent spectra (Section 2.1). The binning is acknowledged as non-optimal and affects the sensitivity estimate.
assumptions (6)
  • domain assumption Hawking radiation formula (Eq. 2.1) correctly describes particle production from an evaporating black hole
    The entire neutrino flux model is built on Hawking's 1974 result; the paper does not derive or validate it.
  • domain assumption The primordial black hole is neutral and non-spinning
    Section 2.1 explicitly assumes a neutral, non-spinning BH; charge or spin would change the emission spectrum and the sensitivity.
  • domain assumption Averaged absorption cross-sections from Halzen et al. 1995 are valid
    The paper uses these cross-sections for Γs(E,M) without re-derivation; errors here propagate to the fluence.
  • ad hoc to paper Time-integrated fluence is well approximated by a broken power law and by three time bins
    Section 2.1 says the fluence is 'well approximated' by a broken power law and then binned into three time intervals; this is an acknowledged, non-optimal modeling choice specific to this analysis.
  • domain assumption PBHs are uniformly distributed in the local universe
    Section 3 uses a uniform distribution to convert the distance sensitivity dsens into a burst rate density.
  • ad hoc to paper A Gaussian time profile is an effective model for PBH evaporation bursts
    Section 2.2 uses the standard Gaussian flare PDF and asserts it works for PBH bursts, but provides no quantitative demonstration.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Neutrinos from Primordial Black Hole Evaporation." pith.science (2026). https://pith.science/paper/LFUC3YYA

@misc{pith2026190805403,
  author       = {Pith},
  title        = {Pith review of: Neutrinos from Primordial Black Hole Evaporation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/LFUC3YYA}},
  note         = {Machine review of arXiv:1908.05403}
}
abstract

Primordial Black Holes (PBHs) are candidates for dark matter as well as ultra-high energy cosmic rays. PBHs are speculated to exist over a large range of masses, from below $10^{15}$ g to $10^3$ M$_\odot$. Here we search for PBHs with an initial mass of $\sim 10^{15}$ g. Hawking radiation by black holes of this initial mass predicts their evaporation at present time. PBHs are expected to produce copious amounts of high-energy neutrinos and gamma rays right before evaporating. Gamma-ray instruments such as Fermi, VERITAS, HAWC, HESS, and Milagro have conducted searches for evaporating PBHs during their last second to a year of existence. They are able to detect bursts from PBHs in a range of $10^{-3}$ to $0.1$ pc. We present sensitivity to PBH evaporation using one year of neutrino data by IceCube. In these proceedings, we detail the changes to adapt IceCube's standard neutrino flare search, aka time-dependent point source search, into one that is appropriate for evaporating BHs. These proceedings serve as proof of concept for a first-ever search for evaporating PBHs using neutrinos that can use 10 years of IceCube data.

Figures

Figures reproduced from arXiv: 1908.05403 by the authors.

Figure 1
Figure 1. The approximated time-integrated neutrino fluence as a function of neutrino energy over [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. The number of signal events required for sensitivity and 5 [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. The left panel (3a) shows the maximum distance to a PBH burst required to achieve 90% median upper limit sensitivity and 5σ discovery potential. The right panel (3b) converts the distances in the left panel to burst rate density. Furthermore, we assume a uniform PBH distribution in the local universe and calculate the burst rate density, ρ˙, using nsens and ndisc in the following relation: ρ˙ = 3 4π nsens d 3 sens 4… view at source ↗

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Hawking emission of massive vector fields by Kerr black holes

    gr-qc 2026-07 conditional novelty 7.0 of 10

    First computation of massive vector (Proca) Hawking emission spectra from Kerr black holes, including polarization-dependent greybody factors, Page functions, and mass-enhanced superradiance up to ~7%.

  2. High-Energy Neutrinos from Black Hole Evaporation in Neutron Stars

    hep-ph 2026-07 conditional novelty 6.0 of 10

    Repeated collapse of asymmetric dark matter inside neutron stars into evaporating microscopic black holes can produce a Galactic-Center-concentrated high-energy neutrino flux at the 10^-12 GeV cm^-2 s^-1 level, subdom...

Reference graph

Works this paper leans on

14 extracted references · 9 canonical work pages · cited by 2 Pith papers

  1. [1]

    S. W. Hawking, Nature 248 (Mar., 1974) 30–31

    work page 1974

  2. [2]

    Carr, arXiv e-prints (Jan, 2019) arXiv:1901.07803

    B. Carr, arXiv e-prints (Jan, 2019) arXiv:1901.07803

    arXiv 2019

  3. [3]

    LIGO Scientific Collaboration and Virgo Collaboration Collaboration, B. P. Abbott et al., Phys. Rev. Lett. 116 (Feb, 2016) 061102

    work page 2016

  4. [4]

    Ackermann et al., Astrophys

    Fermi-LAT Collaboration, M. Ackermann et al., Astrophys. J. 857 (2018) 49. 6 Neutrinos from PBH evaporation Pranav Dave

    work page 2018

  5. [5]

    A. A. Abdo et al., Astropart. Phys. 64 (2015) 4–12

    work page 2015

  6. [6]

    G. T. and, Journal of Physics: Conference Series 375 (jul, 2012) 052024

    work page 2012

  7. [7]

    Glicenstein, A

    J.-F. Glicenstein, A. Barnacka, M. Vivier, T. Herr, and for the H. E. S. S. Collaboration, arXiv e-prints (Jul, 2013) arXiv:1307.4898

    arXiv 2013

  8. [8]

    IceCube Collaboration, M. G. Aartsen et al., JINST 12 (2017) P03012

    2017

Show all 14 references
  1. [9]

    Williams, PoS(ICRC2019)16 (these proceedings)

    IceCube Collaboration, D. Williams, PoS(ICRC2019)16 (these proceedings)

  2. [10]

    Halzen, B

    F. Halzen, B. Keszthelyi, and E. Zas, Phys. Rev. D52 (1995) 3239–3247

    1995

  3. [11]

    IceCube Collaboration, M. G. Aartsen et al., Astropart. Phys. 92 (2017) 30–41

    2017

  4. [12]

    Braun, J

    J. Braun, J. Dumm, F. De Palma, C. Finley, A. Karle, and T. Montaruli, Astroparticle Physics 29 (May, 2008) 299–305

    2008

  5. [13]

    Braun, M

    J. Braun, M. Baker, J. Dumm, C. Finley, A. Karle, and T. Montaruli, Astroparticle Physics 33 (Apr.,

  6. [14]

    IceCube Collaboration, M. G. Aartsen et al., Astrophys. J. 807 (2015) 46. 7

    2015

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.