{"id":"53536ee9-56ec-454a-b5d7-617c01885693","arxiv_id":"2507.01421","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A 0.1% dark matter fraction of 1e-19 to 1e-17 solar mass primordial black holes fits Coma's GeV gamma-ray emission, but evaporating PBHs contribute negligibly to the unresolved gamma-ray background and Draco X shows no gamma-ray excess.","lead":"This paper asks whether tiny primordial black holes evaporating via Hawking radiation could explain the gamma-ray glow of the Coma cluster, and whether cross-correlating gamma-ray maps with lensing can detect them. It finds that a 0.1 percent PBH fraction can be made to fit Coma's GeV emission, but PBH evaporation contributes negligibly to the unresolved gamma-ray background and the Draco X cluster concentration shows no gamma-ray counterpart.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Mass–energy mapping is off by ~100: stated 10^-19–10^-17 Msun PBHs radiate at 0.5–50 MeV, not GeV, so the claimed 5 GeV cutoff and Coma fit are physically inconsistent as written.","rationale":"The reader correctly identifies the assumed PBH IMF as a weak point: the predicted spectral slope and the location of the 5 GeV cutoff depend on n(m) ~ m^-1, and the paper gives no independent justification. However, the more load-bearing problem is that the stated mass range cannot produce the claimed GeV spectrum under standard Hawking physics at all. A 10^-19 M_sun PBH has Hawking temperature ~50 MeV and peak photon energy ~0.1–0.3 GeV; a 10^-18 M_sun PBH peaks near 15 MeV, not 5 GeV. The paper's own statements about a 5 GeV cutoff and a 5 GeV peak for these masses indicate a mass-to-energy unit or conversion error, shifting the predicted spectrum by roughly two orders of magnitude. Additionally, the claimed exhaustion of the low-mass end is inconsistent with the evaporation lifetime of 10^-19 M_sun PBHs, which exceeds the age of the Universe. Because the Coma fit to the Baghmanyan et al. (2022) data is the paper's central claim, this inconsistency invalidates the fit as written. I would therefore recommend rejection unless the authors correct the mass/energy mapping and refit; the UGRB non-detection and Draco X sections may survive, but the title claim does not.","tokens_in":10796,"tokens_out":17366,"duration_ms":202571,"concrete_test":"Run BlackHawk v2 to compute the primary photon spectrum from monochromatic PBH populations at M = 10^-19 M_sun and M = 10^-18 M_sun, normalized to f = 10^-3 of the Coma DM mass at 100 Mpc, and overlay them on the Baghmanyan et al. (2022) data. If the 5 GeV flux is <1% of the 0.2 GeV flux, as expected since E/T ~ 100, the Figure 2 fit and the Sec. 4 'peaks at ~5 GeV' claim are not reproducible, and the mass range must be shifted down by about 2 dex for the central claim to be tested.","verdict_should_be":"REJECT","load_bearing_attack":"The central fit is internally inconsistent at the level of basic Hawking evaporation physics. For a Schwarzschild PBH, T ~ 1.06e13/(M/1 g) GeV. The stated range 10^-19 <= M/M_sun <= 10^-17 corresponds to M ~ 2e14–2e16 g, hence T ~ 50 MeV to 0.5 MeV. The primary photon spectrum peaks near 0.1–0.3 GeV for the lightest end and is exponentially negligible at 5 GeV. Nevertheless, the Figure 2 caption says the cutoff 'just above 5 GeV' is due to exhaustion of lower-mass PBHs, and Sec. 4 states that at m = 10^-18 M_sun the intensity 'peaks at ~5 GeV'. Both assertions require M ~ 10^-21 M_sun, two orders of magnitude below the quoted range. A related inconsistency: a 10^-19 M_sun PBH has an evaporation lifetime ~2.1e10 yr, longer than the age of the Universe, so the low-mass end is not exhausted by z = 0. Consequently, the claimed fit to the Baghmanyan et al. (2022) 0.2–300 GeV data cannot be produced by the stated PBH population; the central claim fails as written unless the mass range is shifted or the energy scale is corrected.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper has three parts. First, it proposes that Hawking radiation from primordial black holes with masses 1e-19 to 1e-17 solar masses, at a fraction f=1e-3 of the Coma cluster dark matter, can explain the cluster's GeV emission measured by Baghmanyan et al. (2022); this fraction is presented as an upper limit because alternative emission from a giant radio halo may also explain the data. Second, it uses 12 years of Fermi-LAT data and DES Y3 weak-lensing measurements, together with a PBH window function, to argue that even with f=1 and a single PBH mass of 1e-18 solar masses, the predicted intensity is about five orders of magnitude below the unresolved gamma-ray background, so the cross-correlation approach cannot constrain PBH evaporation. Third, it identifies a cluster of X-ray clusters, Draco X, at z=0.12, estimates its mass and luminosity, and discusses the absence of corresponding GeV emission.","tokens_in":11133,"tokens_out":11479,"duration_ms":119380,"significance":"If the central calculation were sound, the Coma PBH interpretation would be an interesting application of Hawking radiation to cluster-scale dark matter, and the UGRB cross-correlation analysis would provide a useful null result. The manuscript is transparent that f is fitted and that the UGRB estimate uses f=1 as a conservative upper bound, and it makes use of real Fermi-LAT and DES data. However, the mass-energy mapping in the Coma fit is off by about two orders of magnitude, and the quoted UGRB peak energy is off by about three orders of magnitude. These errors invalidate the quantitative conclusions of Sections 3 and 4, so the paper's current significance is limited.","major_comments":[{"comment":"The central claim that PBHs with masses 1e-19 to 1e-17 solar masses can explain the Coma GeV emission is inconsistent with Hawking evaporation physics. For a Schwarzschild PBH, T = 1.06e13 (1 g/M) GeV, so the stated mass range corresponds to T between about 53 MeV and 0.53 MeV; the primary photon spectrum peaks below about 100 MeV and is exponentially suppressed at 5 GeV. The Figure 2 caption's cutoff 'just above 5 GeV' therefore cannot arise from the stated PBH population. The caption's explanation that the cutoff is due to exhaustion of lower-mass PBHs is also not viable: a 1e-19 solar mass PBH has an evaporation timescale of about 2e10 yr, longer than the age of the Universe. A PBH light enough to emit a 5 GeV peak (M about 1e-21 solar masses) would have evaporated in about 2e4 yr and would not be present at z=0. The mass range and the resulting fit need to be corrected before the Coma interpretation can be evaluated.","section":"Sec. 3, Fig. 2"},{"comment":"The statement after Eq. (2) that for f=1 and m=1e-18 solar masses the average PBH intensity 'peaks at ~5 GeV' is incorrect by about three orders of magnitude: the Hawking temperature for 1e-18 solar masses is about 5.3 MeV, so the primary photon spectrum peaks in the MeV range, not at 5 GeV. Because Eq. (1) uses this spectrum to compute the PBH window function, the quoted intensity of about 1e-13 cm^-2 s^-1 sr^-1 and the conclusion that PBHs contribute negligibly to the UGRB must be recalculated with the correct energy assignment.","section":"Sec. 4, Eq. (2)"},{"comment":"The Draco X mass estimate is internally inconsistent and lacks derivation. The text gives an integrated mass of about 7.2e16 solar masses with a 90% confidence interval ranging from 0.66e16 to 2.0e16 solar masses, but this interval does not contain the quoted central value, and no method is given for either number. This needs to be corrected or substantiated before Draco X can be used as a cosmological probe.","section":"Sec. 5"},{"comment":"The agreement with the Baghmanyan et al. (2022) spectrum is not a quantitative fit: f=1e-3 is chosen to match the observed flux, and no error bars or statistical measure are presented for the predicted spectral shape. Since the assumed initial mass function n(m) proportional to m^-1 determines both the slope and the cutoff location, the reported agreement cannot by itself constrain the PBH fraction or distinguish this IMF from other choices.","section":"Sec. 3"}],"minor_comments":[{"comment":"The quantities dt/dz and e^-tau in Eq. (1) are not clearly defined; tau is said to be integrated along the line of sight for the DES source galaxies, but gamma-ray attenuation should also account for extragalactic background light absorption, which is not discussed.","section":"Eq. (1)"},{"comment":"The y-axis label 'GeV/M /s' and the relation between the 'proportion of PBH fuel consumed' and the plotted curves are confusing; the color-bar labels 'Timestep' should be explained in the caption.","section":"Fig. 1"},{"comment":"The footnote spells the author name as 'Fusco-Fermiano', while the reference list has 'Fusco-Fermiano R.'; please check the correct spelling and make it consistent.","section":"Table 1"},{"comment":"The abstract states 'If 0.1% of the dark matter in the Coma cluster is PBH in this mass range, a fit to the cluster's GeV emission is obtained' as if it were a result; the conditional nature of the assumption should be explicit throughout the text.","section":"Abstract"},{"comment":"The phrase 'a lower limit on its X-ray luminosity, relative to Coma, is estimated to be 58 times greater' is grammatically ambiguous; please clarify whether the lower limit is 58 times the Coma luminosity.","section":"Sec. 5"}],"recommendation":"reject","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe central claim in this paper doesn't survive a basic Hawking-evaporation check. The authors say PBHs in the mass range 10^-19 to 10^-17 Msun produce GeV emission in Coma, with a cutoff just above 5 GeV. For that mass range, the Hawking temperature is 0.5–50 MeV, primary photons peak between roughly 0.1 and 0.3 GeV, and the flux at 5 GeV is exponentially small. The claimed cutoff from 'exhaustion of lower mass PBHs' is also wrong: a 10^-19 Msun PBH has an evaporation lifetime around 2×10^10 yr, longer than the age of the Universe. So the Coma fit, as written, cannot be produced by the stated PBH population. This is a load-bearing error.\n\nWhat is actually new: the specific Coma fit with f=10^-3, the identification of the Draco X overdensity, and a UGRB cross-correlation feasibility statement that applies an existing pipeline. The authors are honest that f is fitted, not predicted, and they treat it as an upper limit because of the alternative radio-halo explanation. Draco X looks like a real empirical structure from the MCXC-II catalog, and the window-function analysis in Sec. 4 is a reasonable application of Thakore et al.\n\nSoft spots, proportionally: the energy-scale error invalidates the headline result and also enters the UGRB calculation, where 10^-18 Msun is claimed to peak at ~5 GeV; the correct peak is about 30 MeV. The qualitative conclusion (f=1 PBHs contribute negligibly to the UGRB in this energy range) probably survives, but the numbers need redoing. The Draco X mass of 7.2×10^16 Msun is given without derivation, and the quoted 90% interval does not bracket the central value—likely a typo. The IMF slope n ∝ m^-1 is assumed from Mould (2025) with no independent justification, and the fit cannot distinguish between IMF choices. There is also no comparison with existing PBH evaporation constraints, which is a serious gap even for a short paper.\n\nWho gains: readers working on PBH constraints might find Draco X worth following, and the corrected UGRB cross-correlation argument could be a useful footnote. The Coma claim, as stated, is not credible.\n\nRecommendation: send it to referees, but ask for a quantitative re-derivation of the evaporation spectrum and mass range. If the authors shift the mass range to ~10^-21–10^-19 Msun (where the cutoff would actually fall near 5 GeV) or reframe it as MeV emission, the fit could become testable. As written, it should not be accepted.","headline":"The Coma fit rests on a factor-of-100 energy error—the claimed 5 GeV cutoff requires PBHs two orders of magnitude lighter than quoted—but the Draco X overdensity and the honest UGRB feasibility analysis give the paper salvageable pieces.","tokens_in":11689,"tokens_out":10320,"would_cite":false,"duration_ms":107789,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Hawking radiation from tiny primordial black holes at 0.1% of Coma's dark matter can explain the cluster's GeV glow.","keywords":["primordial black holes","Hawking radiation","Coma cluster","unresolved gamma-ray background","weak gravitational lensing","dark matter","Fermi-LAT","Draco X"],"falsifier":"A gamma-ray measurement of Coma that resolves the predicted spectral turnover just above 5 GeV would settle this: if the spectrum continues smoothly beyond 5 GeV instead of cutting off, the assumed one-over-mass distribution and the 0.1% abundance fit are ruled out.","tokens_in":10587,"feed_emoji":"🕳️","tokens_out":21481,"duration_ms":187972,"temperature":0.7,"pith_summary":"This paper argues that the GeV (giga-electronvolt) gamma-ray glow of the Coma cluster can be explained by Hawking radiation—the quantum evaporation of black holes—from primordial black holes (PBHs, black holes formed in the early universe) with masses between $10^{-19}$ and $10^{-17}$ solar masses, provided those PBHs make up about 0.1% of the cluster's dark matter. It treats that fraction as an upper limit, because the same GeV emission could instead come from a giant radio halo. The paper then tests whether evaporating PBHs can be seen statistically in cross-correlations between the unresolved gamma-ray background (the diffuse glow not tied to individual sources) and weak gravitational lensing (the subtle distortion of distant galaxy shapes by intervening mass), and concludes that the measured correlation is consistent with known astrophysical sources such as blazars (active galaxies whose jets point toward Earth), not with PBH evaporation. It also identifies a massive cluster-of-clusters, Draco X, that emits no detectable GeV radiation despite being a scaled-up analogue of Coma. If the central claim is right, galaxy clusters become a concrete site for testing this particular PBH dark-matter scenario, while unresolved-background cross-correlations remain too insensitive to see it.","feed_headline":"Black holes in 0.1% of Coma's dark matter could explain GeV glow","feed_subtitle":"Tiny black holes in Coma's dark matter fit its GeV glow and cap their own abundance.","key_machinery":"The machinery that carries the Coma argument is the Hawking evaporation rate $\\dot{m}=-\\alpha\\hbar c^4/(G^2m^2)$, equivalently a fractional mass-loss rate of about $3\\times10^{-23}\\,\\mathrm{s}^{-1}$, which makes lighter black holes evaporate faster and radiate at higher energies; when convolved with the assumed scale-invariant initial mass function $n(m)\\propto m^{-1}$, it sets both the downward slope of the GeV spectrum and the $\\sim5$ GeV cutoff. For the large-scale-structure test, the load-bearing object is the redshift-dependent PBH window function $W_{\\mathrm{PBH}}$, a line-of-sight integral over the evaporation rate with $(1+z)^{-4}$ dimming and attenuation, whose poor overlap with the weak-lensing shear window functions is the reason the cross-correlation approach cannot see PBH evaporation.","core_discovery":"The paper's central claim is that a specific, modest abundance of evaporating PBHs can account for Coma's GeV emission. Using the Hawking mass-loss rate and an assumed initial mass function $n(m)\\propto m^{-1}$, a dark-matter fraction $f=10^{-3}$ in the mass window $10^{-19}{-}10^{-17}\\,M_\\odot$ yields a spectrum that fits the cluster's measured GeV flux, with a turnover just above 5 GeV produced by the exhaustion of the lowest-mass black holes, whose evaporation rate scales as $m^{-2}$. Because a competing explanation (secondary emission from a giant radio halo) can also match the data, the paper frames $f=10^{-3}$ as an upper limit. Extrapolating the same evaporation model to cosmological scales with $f=1$ gives an unresolved gamma-ray background intensity roughly five orders of magnitude below the observed one, so the detected UGRB--lensing cross-correlation (signal-to-noise 8.9) is attributed to clustered astrophysical sources rather than PBHs, and the method is judged unable to constrain PBH evaporation. Separately, the paper reports the identification of Draco X, an overdensity of X-ray clusters at $z=0.12$ with mass around $7\\times10^{16}\\,M_\\odot$, which shows no corresponding excess of extended gamma-ray sources.","pith_inferences":["The $f=10^{-3}$ limit is conditional on the assumed $n(m)\\propto m^{-1}$ mass function; a different mass distribution would change the spectral slope and the inferred abundance, so the limit is only as secure as that external input.","Because the predicted UGRB contribution sits about five orders of magnitude below the measured background, more sensitive cross-correlation analyses are unlikely to reveal PBH evaporation in this mass range; the Coma spectral cutoff and independent abundance probes would be sharper tests.","Draco X's null GeV detection could mean that Coma's glow is produced by the radio halo rather than Hawking radiation, or that Draco X's dark-matter distribution differs from Coma's; matched X-ray and gamma-ray observations of the two systems would separate those options.","A targeted search for the predicted 5 GeV turnover in Coma, with good energy resolution, would provide a nearly parameter-free test of the $m^{-1}$ mass function and the $10^{-3}$ abundance."],"forward_implications":["If the Coma fit holds, the cluster's GeV glow can be powered by Hawking radiation from PBHs with masses $10^{-19}$ to $10^{-17}\\,M_\\odot$ at a dark-matter fraction $f=10^{-3}$, and that fraction is an upper limit if the giant-radio-halo model is correct.","The predicted spectrum declines with decreasing energy according to the assumed $m^{-1}$ mass function and cuts off just above 5 GeV, giving future gamma-ray measurements a sharp feature to confirm or reject.","Attributing a measurable part of the unresolved gamma-ray background to PBH evaporation would require PBH abundances far above what is plausible, so current UGRB--lensing cross-correlations cannot constrain evaporating PBHs in this mass range.","Because the PBH gamma-ray window and the shear window overlap only weakly in redshift, deeper lensing surveys under these assumptions will not substantially improve sensitivity to PBH evaporation.","A Coma-like Draco X should emit more than twice Coma's GeV flux, yet no extended gamma-ray source excess is seen there, implying that the GeV emission mechanism is not universal across massive clusters."],"supporting_citations":[{"why":"Introduces Hawking radiation, the quantum evaporation mechanism at the center of the paper's Coma interpretation.","marker":"Hawking (1974)"},{"why":"Supplies the evaporation mass-loss rate, including photon and particle channels, used for the spectra and window functions.","marker":"Mosbech & Picker (2022)"},{"why":"Provides the assumed initial mass function $n(m)\\propto m^{-1}$ for the PBH population.","marker":"Mould (2025)"},{"why":"Gives the Coma GeV flux measurements from 200 MeV to 300 GeV that the model is fit to.","marker":"Baghmanyan et al. (2022)"},{"why":"Proposes a giant radio halo as an alternative origin for the same GeV emission, making $f=10^{-3}$ an upper limit.","marker":"Kushnir et al. (2024)"},{"why":"Supplies the Planck 2018 cosmological parameters used in the window-function calculations.","marker":"Aghanim et al. (2020)"},{"why":"Provides the 12-year gamma-ray and weak-lensing cross-correlation pipeline and the 8.9-sigma UGRB-lensing detection.","marker":"Thakore et al. (2025)"},{"why":"Catalogs the X-ray clusters from which the Draco X concentration is identified.","marker":"Sadibekova et al. (2024)"},{"why":"Provides the source-galaxy redshift distributions used for the weak-lensing shear window functions.","marker":"Myles et al. (2021)"}],"fun_headline_variants":["Coma's GeV glow might come from primordial black holes","0.1% PBH fraction in Coma fits its gamma-ray emission","Primordial black holes in Coma could explain GeV signal","Evaporating black holes in Coma: a possible GeV source","Coma's GeV glow consistent with primordial black holes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes the black holes are distributed in mass as one-over-mass, a power law taken from an earlier paper, and everything about the predicted spectrum and the 0.1% abundance limit follows from that choice; if the true mass distribution is different, the fit changes, and the paper gives no independent evidence for this power law.","fun_headline_variants_meta":{"raw":{"variants":["Coma's GeV glow might come from primordial black holes","0.1% PBH fraction in Coma fits its gamma-ray emission","Primordial black holes in Coma could explain GeV signal","Evaporating black holes in Coma: a possible GeV source","Coma's GeV glow consistent with primordial black holes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000749,"raw_usage":{"total_tokens":3448,"prompt_tokens":1170,"completion_tokens":2278,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":786,"completion_tokens_details":{"reasoning_tokens":2192}},"tokens_in":786,"tokens_out":2278,"duration_ms":19591,"temperature":1.0,"reasoning_tokens":2192,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:52:59.763635+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A gamma-ray measurement of Coma that resolves the predicted spectral turnover just above 5 GeV would settle this: if the spectrum continues smoothly beyond 5 GeV instead of cutting off, the assumed one-over-mass distribution and the 0.1% abundance fit are ruled out.","supporting_citations":[{"cited_title":"W., 1974, Nature, 248, 30","cited_arxiv_id":null,"evidence_quote":"Introduces Hawking radiation, the quantum evaporation mechanism at the center of the paper's Coma interpretation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the evaporation mass-loss rate, including photon and particle channels, used for the spectra and window functions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the assumed initial mass function $n(m)\\propto m^{-1}$ for the PBH population."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the Coma GeV flux measurements from 200 MeV to 300 GeV that the model is fit to."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Planck 2018 cosmological parameters used in the window-function calculations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the 12-year gamma-ray and weak-lensing cross-correlation pipeline and the 8.9-sigma UGRB-lensing detection."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Catalogs the X-ray clusters from which the Draco X concentration is identified."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the source-galaxy redshift distributions used for the weak-lensing shear window functions."}],"review_version":1}