{"id":"613f6a91-f2e6-45f0-910c-4d9735d5ee1d","arxiv_id":"2504.14185","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":11,"one_line_summary":"PBH-emitted ALPs converting in cosmic magnetic fields, and PBH electrons scattering off ALP dark-matter halos, could create gamma-ray signals detectable by AMEGO, e-ASTROGAM, and MAST, with projected coupling reach near 1e-13 GeV^-1.","lead":"Scientists calculate whether future gamma-ray telescopes could spot axion-like particles or primordial black holes through the particles the black holes evaporate into. The paper estimates new detection reach for three proposed satellites and finds the low-energy AMEGO telescope could improve current axion-photon coupling limits by roughly an order of magnitude.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The conversion probabilities that set the gaγγ forecast assume every extragalactic PBH sits in a 10 μG cluster core and every galactic PBH follows one GC-to-Sun path; realistic line-of-sight averaging could shift the projected limits by a large factor.","rationale":"The paper does the standard thing: BlackHawk spectra, standard ALP-photon mixing equations, and Fisher forecasts with profiled nuisance parameters; these are not the problem. The problem is that every forecast limit is proportional to a conversion probability that is effectively assumed, not computed from a realistic distribution of PBH sources and magnetic fields. Eq. (16) multiplies the galactic and extragalactic PBH fluxes by single P_gal and P_exgal values. The text in Sec. III.A explicitly says galactic PBHs are treated with one GC-to-Sun path and extragalactic PBHs are assumed to sit at cluster centers. Those assumptions are in tension with the actual source density: NFW galactic PBHs are spread along the line of sight, and cosmological PBHs trace the cosmic web, where cluster cores occupy a small volume and mass fraction. For MPBH=3e16 g, both flux terms are comparable, so an order-of-magnitude error in P_exgal translates almost linearly into the gaγγ limit. The additional choice of the median over random field orientations, followed by Gaussian smoothing, is also not the ensemble average needed for an incoherent sum of many ALP paths. The reader's conditional verdict is appropriate; it should remain conditional pending this check. If the check shows the limit is stable to within a factor of about 2, I would be more comfortable accepting the headline claim as a valid forecast.","tokens_in":18153,"tokens_out":11052,"duration_ms":111387,"concrete_test":"Recompute the evaporation-conversion signal (Eq. (16)) with a Monte Carlo over source positions and magnetic environments: for galactic PBHs, propagate each source to Earth through the full 3D JF12 field (including halo/X components) instead of the single GC-Sun path; for extragalactic PBHs, draw source environments from a realistic cluster/filament/void occupancy (e.g., cluster mass fraction ~10-20% with the adopted B0=10 μG profile, filaments with B~0.1-1 μG, voids with nG fields), and use the arithmetic mean of P over realizations, not the median, without Gaussian smoothing. Redo the Fisher forecast of Fig. 6 (right panel) at MPBH=3×10^16 g, a*=0; if the projected gaγγ limit at ma<1e-10 eV shifts by more than a factor of 3, the claimed order-of-magnitude improvement is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Sec. III.A, the signal in Eq. (16) is built from P_gal computed on a single Galactic Center-to-Sun path (8.2 kpc, eight spiral-arm domains, Eq. (15) with ψ=0) and P_exgal computed assuming every extragalactic PBH is at a galaxy-cluster center, propagates through 100 domains of a 10 μG turbulent field, and contributes the median of 1000 realizations, Gaussian-smoothed (Eq. (14) and surrounding text). This is not the physical per-source line-of-sight integral. Most extragalactic PBHs do not reside in 10 μG cluster cores; the cluster mass fraction is small, while the EBL contribution is dropped as negligible, so P_exgal is likely overestimated for most sources. For MPBH=3×10^16 g, the extragalactic column is comparable to the Galactic-center column, so the AMEGO limit in Fig. 6 is directly affected. The same smoothing removes real oscillations; using the mean rather than the median of P over random field realizations is what a large number of incoherent sources actually requires. If these choices change the effective conversion probability by even a factor of 3, the claimed improvement over current astrophysical bounds at gaγγ~1e-13 GeV^-1 no longer holds as stated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies two gamma-ray production mechanisms involving primordial black hole (PBH) evaporation and axion-like particles (ALPs). In the first, \"evaporation-conversion\" scenario, PBH-emitted ALPs propagate through cosmic magnetic fields and convert into photons; in the second, \"evaporation-scattering\" scenario, PBH-emitted relativistic electrons scatter off a Galactic ALP DM halo. Using BlackHawk for Hawking spectra, a transfer-matrix treatment for ALP-photon conversion, and a Fisher forecasting procedure, the authors derive projected 95% C.L. sensitivities to the PBH DM fraction f_PBH and the ALP-photon coupling g_aγγ for the future AMEGO, e-ASTROGAM, and MAST telescopes. The central quantitative claims are that in the conversion channel AMEGO can probe g_aγγ down to about 1e-13 GeV^-1 for m_a < 1e-10 eV with M_PBH = 3e16 g, improving on current astrophysical bounds by one order of magnitude, while MAST gives complementary coverage at higher masses; and that in the scattering channel PBHs of about 1e14 g can give limits comparable to the strongest cosmic-ray-scattering constraints, though the corresponding parameter space is already excluded.","tokens_in":18540,"tokens_out":8270,"duration_ms":79949,"significance":"If the projected sensitivities are robust, the paper offers a genuinely new window onto light ALPs: Hawking radiation from asteroid-mass PBHs provides an ALP source whose subsequent conversion in magnetic fields can be probed by future MeV-GeV satellites, and the same signal carries information about f_PBH. The use of established tools (BlackHawk, the Raffelt-Stodolsky transfer matrix, and publicly available detector responses) is a strength, as is the explicit, falsifiable nature of the forecasts. The main caveat is that the conversion probability that sets the entire g_aγγ reach is computed under strong idealized assumptions about the astrophysical magnetic-field environments and source locations; if those assumptions are relaxed, the magnitude and even the existence of the claimed one-order-of-magnitude improvement is uncertain. The paper is therefore interesting and timely, but its headline projection is not yet demonstrated to be robust.","major_comments":[{"comment":"","section":"Sec. III.A, Eqs. (14)-(16), Fig. 6"},{"comment":"","section":"Sec. III.B, Eq. (20), Fig. 7"},{"comment":"","section":"Sec. III.A, Eq. (14) and Fig. 2"}],"minor_comments":[{"comment":"","section":"Fig. 5 (right)"},{"comment":"","section":"Sec. III.A"},{"comment":"","section":"Sec. IV, Eq. (24)"},{"comment":"","section":"Fig. 6 (right)"},{"comment":"","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a timely and interesting question, and the central mechanism (PBH Hawking emission of ALPs followed by conversion in cosmic magnetic fields) is a plausible new probe. However, the headline claim of a one-order-of-magnitude improvement over current constraints rests on conversion probabilities that are computed under very specific, idealized assumptions about source locations and magnetic-field geometries. The authors are transparent about these assumptions, but they do not quantify the associated uncertainty or provide a conservative estimate. I would encourage the editor to request a revision that includes a realistic treatment of the extragalactic source population and a robustness study of the magnetic-field parameters, or that lowers the claimed reach accordingly."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does something genuinely new: it puts PBH-emitted ALPs through magnetic-field conversion and PBH-emitted electrons through inverse Primakoff scattering off ALP DM, then projects what AMEGO, e-ASTROGAM, and MAST could see. The two channels are not in the cited literature, and the quantitative forecasts for fPBH and gaγγ are a natural extension of the authors' earlier PBH-ALP work and of cosmic-ray electron scattering studies. Credit where due: the core equations are standard and correctly assembled, the Hawking fluxes come from BlackHawk, the Fisher procedure is sensible, and the paper does not invent new entities or fit the target result. The figures are readable and the comparison to existing bounds is helpful.\n\nThe soft spots are real and concentrated in the conversion geometry. The Galactic signal uses one fixed GC-to-Sun path for every PBH in the Milky Way, regardless of position. The extragalactic signal assumes every PBH sits at the center of a galaxy cluster with a 100-domain, 10 µG turbulent field, takes the median over 1000 random realizations, and then Gaussian-smooths it. Most extragalactic PBHs are not in cluster cores, and the EBL contribution is dropped as negligible, so this likely overestimates the conversion probability for the typical source. If the effective P_a→γ changes by a factor of a few, the headline claim of an order-of-magnitude improvement over current astrophysical bounds at gaγγ ~ 1e-13 GeV^-1 no longer holds as stated. The electron-cooling neglect in the scattering channel is a smaller issue, but it is not quantified, and a code release would help others test the geometry dependence.\n\nThe stress-test concern lands. It is not a manufactured flaw; it is the load-bearing approximation of the forecast. I would not call the paper fatally flawed, because the two channels are still worth exploring and the equations are sound. But the projected limits should be treated as conditional on a more realistic line-of-sight or Monte Carlo treatment of the magnetic field, and on energy-loss estimates for the scattering channel.\n\nWho is this for? People planning MeV/GeV gamma-ray missions and anyone working on PBH or ALP dark matter. The paper deserves a serious referee, but the referee should push for actual line-of-sight averaging, a sensitivity check on the cluster-field assumptions, and a public code release. With those changes, the projections could become a solid reference; without them, the quantitative claims are too fragile to cite as-is.","headline":"A useful but geometry-sensitive forecast paper: the new channels are real, the projected one-order improvement in gaγγ rests on optimistic conversion assumptions that need line-of-sight averaging before the numbers can be trusted.","tokens_in":19030,"tokens_out":1245,"would_cite":true,"duration_ms":14040,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["95.35.+d","14.80.Va","04.70.Dy","95.85.Pw"],"model":"deepseek-v4-flash","headline":"A MeV satellite could spot axions from black-hole evaporation","keywords":["axion-like particles","primordial black holes","Hawking radiation","ALP-photon conversion","gamma-ray astronomy","Fisher forecast","dark matter","satellite gamma-ray telescopes"],"falsifier":"Compute the conversion signal using the actual PBH spatial distribution and a modern Bayesian galactic magnetic-field model, assigning a separate line-of-sight conversion probability to every PBH rather than one shared path; if the flux averaged over the $|l|\\le5^\\circ$, $|b|\\le5^\\circ$ region differs from the paper's single-path value by more than about a factor of two, the projected $g_{a\\gamma\\gamma}$ and $f_{\\mathrm{PBH}}$ limits in Section IV are not robust.","tokens_in":17920,"feed_emoji":"🔭","tokens_out":7525,"duration_ms":64280,"temperature":0.7,"pith_summary":"This paper asks whether evaporating primordial black holes, one dark-matter candidate, can be turned into a source of axion-like particles, another dark-matter candidate. It works out two gamma-ray production paths: ALPs emitted by PBHs convert into photons in cosmic magnetic fields, and PBH-emitted electrons scatter off an ALP dark-matter halo, producing gamma rays through the ALP-photon coupling. Forecasting one year of exposure for the proposed AMEGO, e-ASTROGAM, and MAST satellites, it projects sensitivity to the ALP-photon coupling $g_{a\\gamma\\gamma}$ down to about $10^{-13}\\,\\mathrm{GeV}^{-1}$ for ALP masses below $10^{-10}$ eV, roughly an order of magnitude better than current astrophysical limits. If the forecast holds, the same satellite data would also constrain the PBH dark-matter fraction $f_{\\mathrm{PBH}}$, joining two dark-matter programmes in a single signal.","feed_headline":"A MeV satellite could spot axions from black-hole evaporation","feed_subtitle":"PBH-emitted axions turning into gamma rays would beat current bounds by one order.","key_machinery":"The argument runs on the propagation-matrix treatment of ALP-photon conversion in a turbulent magnetized medium. A Schr\\\"odinger-like mixing equation for the photon and ALP amplitudes is integrated domain by domain: 100 randomly oriented coherent cells for a galaxy cluster, and eight spiral-arm segments along a fixed Galactic Center-to-Sun path in a simplified JF12 galactic-field model. The cluster conversion probability is the median of 1000 Monte Carlo realizations, and both galactic and extragalactic probabilities are Gaussian-smoothed before being folded into the signal. Hawking emission spectra are generated with the BlackHawk v2.1 code, and the projected limits come from a profiled Fisher information matrix that treats astrophysical-background parameters as nuisance variables. For the scattering channel, the load-bearing object is the differential cross section of the inverse Primakoff process $e^- + a \\to e^- + \\gamma$, integrated over the PBH electron flux and the ALP DM halo density.","core_discovery":"The central claim is that current and near-future gamma-ray satellites can test two dark-matter candidates at once through Hawking radiation from primordial black holes. In the evaporation-conversion channel, ALPs of mass $m_a$ emitted by PBHs of mass $M_{\\mathrm{PBH}}$ oscillate into photons in the Milky Way's magnetic field and in galaxy-cluster magnetic fields; the paper finds that AMEGO can reach $g_{a\\gamma\\gamma}\\sim 1\\times 10^{-13}\\,\\mathrm{GeV}^{-1}$ for $m_a<10^{-10}$ eV at $M_{\\mathrm{PBH}}=3\\times10^{16}$ g and zero PBH spin, improving on existing astrophysical constraints by an order of magnitude, while MAST extends coverage to $m_a>10^{-10}$ eV. In the evaporation-scattering channel, relativistic electrons from PBHs near $10^{14}$ g scatter on non-relativistic ALP dark matter and produce gamma rays, yielding projected constraints on $g_{a\\gamma\\gamma}$ comparable to the strongest cosmic-ray-scattering limits, although the parameter space for heavier PBHs is already excluded by other bounds. The paper also derives projected 95% confidence upper limits on the PBH fraction $f_{\\mathrm{PBH}}$ for both channels.","pith_inferences":["Inference: because the conversion signal and the direct PBH gamma-ray background both scale with $f_{\\mathrm{PBH}}$, a simultaneous spectral fit could separate the conversion bump from the Hawking continuum; this separation could be tested with existing detector simulations before launch.","Inference: a multi-line-of-sight treatment that assigns each galactic PBH its own conversion path would likely smooth the oscillatory probability and could shift the projected limits by a factor of a few, so the order-of-magnitude improvement is a target for re-analysis rather than a fixed promise.","Inference: if future surveys tighten the allowed PBH abundance, the same AMEGO and MAST data would automatically convert those bounds into stronger ALP limits, linking PBH searches with axion searches in a way the paper only partially exploits."],"forward_implications":["A non-detection by AMEGO in its 150 keV to 5 MeV window would exclude ALP-photon couplings above roughly $10^{-13}\\,\\mathrm{GeV}^{-1}$ for ultralight ALPs from $3\\times10^{16}$ g PBHs, a decade beyond current astrophysical bounds.","MAST's large effective area at 100 MeV to 3 GeV gives complementary sensitivity for ALP masses above $10^{-10}$ eV, where AMEGO's reach declines.","The evaporation-scattering channel would place PBH-electron constraints on $g_{a\\gamma\\gamma}$ comparable to cosmic-ray-scattering limits for PBHs near $10^{14}$ g, effectively converting a PBH-abundance bound into an ALP-coupling probe.","Both channels turn one observing campaign into joint limits on $f_{\\mathrm{PBH}}$ and $g_{a\\gamma\\gamma}$, so a single future gamma-ray telescope could constrain two dark-matter candidates simultaneously."],"supporting_citations":[{"why":"Establishes the ALP-photon mixing Hamiltonian and propagation framework that underlies the conversion probability.","marker":"[84]"},{"why":"Supplies the JF12 Galactic magnetic-field model whose simplified eight-segment version sets the Milky Way conversion probability.","marker":"[88]"},{"why":"Provides the turbulent multi-domain cluster model used for the 100-cell extragalactic conversion calculation.","marker":"[92]"},{"why":"Sets the current extragalactic gamma-ray and CMB bounds on fPBH that anchor the benchmark values and comparisons.","marker":"[36]"},{"why":"Supplies the BlackHawk v2.1 code that generates the Hawking-emission spectra of ALPs, electrons, and photons.","marker":"[78]"},{"why":"Updates the BlackHawk code used for the evaporation spectra in the numerical analysis.","marker":"[79]"},{"why":"Defines AMEGO's effective area and energy window used in the low-energy Fisher forecast.","marker":"[67]"},{"why":"Defines MAST's large effective area and high-energy window used for the complementary forecast.","marker":"[70]"},{"why":"Provides the electron-ALP scattering cross section and kinematics on which the evaporation-scattering channel is built.","marker":"[64]"},{"why":"Supplies the current cosmic-ray-scattering limits against which the projected scattering-channel constraints are compared.","marker":"[65]"},{"why":"Provides the Fisher forecasting formalism, including nuisance profiling, used to project the limits.","marker":"[106]"}],"fun_headline_variants":["AMEGO could detect axions from PBH evaporation, beating limits by 10x","Two dark-matter candidates linked: PBH evaporations could illuminate axions","MeV gamma-ray satellites might reveal axions from evaporating black holes","PBH evaporation could expose axion dark matter through gamma-ray signals"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Every galactic PBH is assumed to lie on one fixed Galactic Center-to-Sun path and every extragalactic PBH is assumed to sit at the center of a galaxy cluster with a prescribed 100-domain turbulent field, so if real sightlines or field configurations change the conversion probability by even a factor of a few, the projected limits shift.","fun_headline_variants_meta":{"raw":{"variants":["AMEGO could detect axions from PBH evaporation, beating limits by 10x","Two dark-matter candidates linked: PBH evaporations could illuminate axions","MeV gamma-ray satellites might reveal axions from evaporating black holes","PBH evaporation could expose axion dark matter through gamma-ray signals"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00091,"raw_usage":{"total_tokens":3940,"prompt_tokens":1004,"completion_tokens":2936,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":620,"completion_tokens_details":{"reasoning_tokens":2853}},"tokens_in":620,"tokens_out":2936,"duration_ms":20252,"temperature":1.0,"reasoning_tokens":2853,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:55:16.707222+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the conversion signal using the actual PBH spatial distribution and a modern Bayesian galactic magnetic-field model, assigning a separate line-of-sight conversion probability to every PBH rather than one shared path; if the flux averaged over the $|l|\\le5^\\circ$, $|b|\\le5^\\circ$ region differs from the paper's single-path value by more than about a factor of two, the projected $g_{a\\gamma\\gamma}$ and $f_{\\mathrm{PBH}}$ limits in Section IV are not robust.","supporting_citations":[{"cited_title":"Probing axion-like particles through the gamma-ray production from cosmic-ray scattering in the Milky Way dark matter halo","cited_arxiv_id":"2501.11691","evidence_quote":"Supplies the current cosmic-ray-scattering limits against which the projected scattering-channel constraints are compared."}],"review_version":1}