{"id":"0ecdd9ad-ac1c-4df8-81dc-89484d9fd11c","arxiv_id":"2505.04692","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"AMS-02 antiproton data set the tightest limits yet on Galactic primordial black holes with lognormal mass distributions, and cap the expected antideuteron flux below the reach of upcoming detectors.","lead":"Primordial black holes, if they exist, would evaporate and spray out antimatter particles; this paper works out the expected antiproton and antideuteron fluxes and compares them with AMS-02 measurements. The result is a tight upper limit on how many such black holes can be in the Galaxy, plus a clear statement that any detected antideuterons would imply physics beyond standard cosmic-ray production.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quoted PBH density bounds assume standard Hawking evaporation; memory-burden backreaction could suppress the M^2 tail and weaken the limits by orders of magnitude.","rationale":"Stress-test pass: The paper's central claim is the set of AMS-02 upper limits on the local PBH density from antiproton data, plus the prediction that any antideuteron detection would be only partly attributable to PBH evaporation. The calculations are carefully executed within the standard Hawking picture: mass evolution, fragmentation functions, propagation, and likelihood are all state-of-the-art, and the covariance treatment follows established work (Refs. [91,96]). The 'universal spectral shape' is a useful simplification and the authors report percent-level checks. However, the entire chain depends on the unmodified Hawking evaporation rate. The paper itself flags the memory-burden caveat (Sec. II A) but does not quantify how the bounds would change if evaporation is slowed. Since the masses probed (mu_c ~ M* = 5e14 g) are precisely those whose lifetime equals the age of the Universe, a slowdown due to memory burden would shift the evaporating population to heavier masses and suppress the M^2 tail that produces antinuclei; the quoted bounds could weaken by orders of magnitude. This is the single most load-bearing assumption. The reader's weakest_assumption identifies the same point, so agreement is 'agree'. The concern does not invalidate the paper's internal consistency, but it justifies the CONDITIONAL verdict; no change from the reader's verdict is needed.","tokens_in":22913,"tokens_out":11451,"duration_ms":116179,"concrete_test":"Recompute the antiproton source spectra and the 95% CL upper limits of Sec. IV B using a memory-burden-modified evaporation rate, e.g., dM/dt = -alpha(M)/M^2 / [1 + (M/M_b)^n] with M_b = q M_in, q in [0.5, 1], n in [1, 2], as in Thoss et al. (2024) arXiv:2402.17823, for mu_c = 5e14 g and sigma = 0.1, 0.5, 1, 2. Use the same USINE propagation and likelihood framework. If any of the four headline limits weakens by more than an order of magnitude, the central claim is conditional on standard Hawking evaporation and the paper should state the memory-burden dependence explicitly in the abstract and conclusions.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Sec. II A (after Eq. (6)) assumes the semi-classical picture holds up to Planck scale, with dM/dt = -alpha(M)/M^2 throughout the PBH lifetime. This underlies the mass evolution Eq. (9), the source spectra Eq. (12), and all upper limits in Sec. IV B, including the headline numbers for mu_c = M* = 5e14 g (3.5e-12, 2.2e-11, 4.4e-11, 8.1e-11 GeV cm^-3 for sigma = 0.1, 0.5, 1, 2). The paper explicitly flags that memory-burden backreaction (Refs. [76-82]) could slow evaporation and extend PBH lifetimes. If memory burden is operative, PBHs of initial mass ~5e14 g would not be completing evaporation today; the universal M^2 low-mass tail in the evolved distribution (Fig. 1) that drives antinuclei production would be suppressed, and the derived bounds would weaken or disappear. Because the 'universal spectral shape' argument in Sec. II B and the antideuteron forecast in Sec. V both rely on this tail, the central claim's robustness rests on an assumption the authors have not quantified. This is a correctness risk rather than an internal inconsistency: within standard Hawking evaporation the calculations are coherent.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes Galactic antiproton and antideuteron fluxes from Hawking evaporation of primordial black holes (PBHs) with a lognormal initial mass distribution, using public codes (BlackHawk, CosmiXs, USINE) and state-of-the-art propagation models. It then compares the predicted antiproton flux to AMS-02 data using the likelihood framework of Ref. [91] and derives 95% CL upper limits on the local PBH density as a function of the lognormal peak mass and width. The authors find that AMS-02 antiproton data exclude PBH densities far below those probed by other messengers, translate the bounds into fractions of local dark matter and into a local PBH explosion rate, and argue that the primary PBH antinuclei fluxes share a universal spectral shape. They also compute the maximal antideuteron flux allowed by the antiproton bounds and conclude that a future AMS-02 or GAPS antideuteron detection would be a new-physics signal that could only partly be explained by PBH evaporation.","tokens_in":23192,"tokens_out":12340,"duration_ms":134671,"significance":"If the results hold, the paper provides competitive and robust constraints on PBH abundances from a single well-measured cosmic-ray channel, and it produces a falsifiable prediction: any detected antideuteron flux above the secondary background would require non-PBH new physics. The analysis is built on public, reproducible tools and follows an established statistical pipeline, which is a genuine strength. The claimed universal spectral shape of PBH antinuclei fluxes is a useful simplification that connects the antiproton bounds to antideuteron forecasts in a transparent way. The central claims are not circular: the PBH density limits come from AMS-02 data, and the antideuteron prediction imports those limits as an input.","major_comments":[{"comment":"The headline constraints, including the values quoted in Sec. VI (rho_PBH < 3.5e-12, 2.2e-11, 4.4e-11, 8.1e-11 GeV/cm3 for sigma = 0.1, 0.5, 1, 2 at mu_c = M*), are derived from standard Hawking evolution dM/dt = -alpha(M)/M^2 assumed to hold for the entire PBH lifetime, and from the resulting universal M^2 low-mass tail shown in Fig. 1. The paper explicitly notes that memory-burden backreaction (Refs. [76-82]) could slow evaporation and extend PBH lifetimes, but it does not quantify the effect on Eq. (9), Eq. (12), the UL curves in Fig. 4, the explosion-rate bound in Eq. (21), or the antideuteron forecast in Sec. V. If memory burden is operative before PBHs of initial mass ~5e14 g complete evaporation, the low-mass tail that dominates antinuclei production would be suppressed and the quoted limits could be substantially weakened or shifted. I recommend adding a quantitative estimate, for example using a minimal memory-burden parametrization with a transition mass or occupation number as in Refs. [79-82], and/or stating prominently that the constraints apply only within standard Hawking evaporation and are not robust to memory-burden scenarios.","section":"Sec. II A (after Eq. (6)) and Sec. IV B"},{"comment":"The likelihood uses a covariance matrix C taken from the secondary-antiproton analysis of Ref. [96], with the justification that the PBH primary flux is subdominant. At the 95% upper limit, however, the primary flux is by construction large enough to be constrained, and its normalization inherits the modeling choices of footnote 4 (division by 2 for the single-jet treatment) as well as uncertainties in the Hawking spectra and fragmentation functions. A factor-2 error in the single-jet approximation would directly rescale all the upper limits in Fig. 4 by a factor of 2. Please quantify this systematic, or state the resulting uncertainty band on the quoted ULs, so that the comparison with other messengers in Fig. 5 is not read as more precise than the underlying source modeling allows.","section":"Sec. IV A, Eq. (18) and footnote 4"}],"minor_comments":[{"comment":"The text says that adding simulated GAPS data would 'increase the upper bounds by a factor 2', while the conclusions say GAPS 'should improve the above limits by a factor 2'; please clarify whether the combined AMS-02 + GAPS upper limit is larger or smaller.","section":"Sec. IV B, paragraph after Fig. 4"},{"comment":"The sentence about 'fusion modeling of the antinuclei to for and antideuteron' appears garbled and should be rephrased.","section":"Sec. V, paragraph after Fig. 6"},{"comment":"The phrase 'we assume the second one not to escape the PBH' is unclear; it should be stated more precisely that one of the two jets is effectively reabsorbed and only a single jet contributes to the hadronization spectrum.","section":"Footnote 4"},{"comment":"The text states the strongest f_PBH bound is at mu_c = M*, but the left panel of Fig. 5 starts at 1e15 g; please make the displayed range and the quoted value consistent.","section":"Sec. V and Fig. 5"},{"comment":"The symbol L is used both for the halo half-height and for the log-likelihood function, which makes the notation in Eq. (18) confusing; consider using a calligraphic symbol for the likelihood.","section":"Eq. (18)"}],"recommendation":"major_revision","confidential_remarks":"This is a careful, well-executed phenomenological paper that uses established public tools and a validated statistical framework. My main concern is the gap between the strong, precisely quoted constraints and the acknowledged memory-burden caveat: the authors cite the relevant literature but do not engage with it quantitatively. A revision that adds a memory-burden robustness check, or that clearly reframes all central claims as conditional on standard Hawking evaporation, would make the paper acceptable. The paper is well within the scope of the journal, and I see no citation or novelty concerns."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a careful, credible update of PBH antinuclei constraints. The genuinely new pieces are the lognormal mass distribution for antinuclei, the use of the 2021 AMS-02 antiproton sample with a covariance-based likelihood, and the explicit demonstration that all primary PBH antinuclei fluxes collapse to one spectral shape up to a normalization. That last point is more than a curiosity — it lets the AMS-02 antiproton bound be rescaled directly to antideuterons, which is why the forecast in Sec. V is as clean as this business gets. The central claim holds: within standard Hawking evaporation, AMS-02 antiprotons exclude PBHs as all of the local dark matter for the lognormal parameter space, and a future antideuteron detection would be mostly new physics beyond PBH evaporation.\n\nThe machinery is assembled honestly. BlackHawk, CosmiXs, and USINE are public; the propagation setup follows the recent AMS-02 antiproton analysis of Ref. [91]; the covariance matrix carries the main systematics. The paper flags its own caveats: no memory burden, the single-jet division for CosmiXs fragmentation functions, and the force-field solar modulation assumption. That is the right scientific posture.\n\nSoft spots, in proportion. The memory-burden point is the main one. The low-mass tail that actually produces antinuclei is generated under dM/dt ∝ M^{-2} up to the Planck scale. If memory burden slows evaporation, that tail changes and the bounds could weaken by orders of magnitude. The paper mentions this but does not quantify it. Not fatal — the assumption is standard and clearly stated — but a referee should ask how the headline upper limits respond to a slowdown. The factor-two difference between the 2016 and 2021 AMS-02 datasets also deserves more comment; the reader should know which dataset is doing the heavy lifting. The custom convolution code is not released. Everything else is public, so this is minor, but the authors should be asked to share it.\n\nCitation pattern is clean. Self-citations point to the actual propagation and antinuclei tools used, and the comparison with other messengers is properly qualified. The paper overstates nothing.\n\nWho this is for: anyone in PBH dark matter, cosmic-ray antinuclei, or GAPS/AMS-02 science. It deserves a serious referee — send it out.","headline":"Solid, transparent update of PBH antinuclei constraints; the memory-burden caveat is the one thing a referee should push on.","tokens_in":23743,"tokens_out":3909,"would_cite":true,"duration_ms":35042,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["98.70.Sa","95.35.+d","04.70.Dy"],"model":"deepseek-v4-flash","headline":"AMS-02 antiproton data set upper bounds on the local density of evaporating primordial black holes, and any future antideuteron detection would be a new-physics signal only partly attributable to PBH evaporation.","keywords":["primordial black holes","Hawking evaporation","antiprotons","antideuterons","cosmic rays","dark matter","lognormal mass distribution","AMS-02"],"falsifier":"Detect antideuterons: the paper predicts a maximum antideuteron flux obtained by saturating the antiproton bound, so an observed rate above that maximum excludes PBH evaporation as the source, while a clear non-detection at the predicted peak energy would tighten the density bound. A second check is whether PBHs near $5\\times10^{14}$ g still evaporate on the Universe's timescale; if memory-burden effects make them live longer, the predicted fluxes disappear.","tokens_in":22704,"feed_emoji":"🕳️","tokens_out":12152,"duration_ms":101795,"temperature":0.7,"pith_summary":"This paper asks how much of the Milky Way's dark matter could be made of light primordial black holes (PBHs) that are evaporating today, and answers with a tighter number using AMS-02 antiproton measurements. It computes the antiproton and antideuteron fluxes from Hawking evaporation for a lognormal PBH mass distribution, propagates them through Galactic transport models, and compares the result with the measured antiproton flux. The central result is an upper bound on the local PBH density: for a critical mass $\\mu_c \\simeq 5\\times 10^{14}$ g, the allowed density ranges from about $3.5\\times10^{-12}$ to $8.1\\times10^{-11}$ GeV cm$^{-3}$, depending on the width $\\sigma$ of the mass distribution. A second result is that all primary PBH antinuclei fluxes share a universal spectral shape, so a future antideuteron detection would not be attributable to PBH evaporation alone; it would be a signal of new physics whose PBH component is limited by the antiproton bounds.","feed_headline":"Antiproton data cap the density of evaporating black holes","feed_subtitle":"The AMS-02 flux limits leave no room for PBHs as dark matter and imply any antideuteron sighting is new physics.","key_machinery":"The central object is the evolved mass distribution of evaporating primordial black holes. The paper assumes semi-classical Hawking evaporation, $dM/dt = -\\alpha(M)/M^2$, holds for the entire PBH lifetime, and shows that any extended distribution develops a common low-mass tail proportional to $M^2$ below about $10^{13}$ g; that tail, made of PBHs evaporating now, dominates antinucleus production. The source term factorizes as $Q(r,z,E) = q(E)\\,\\tilde{\\rho}(r,z)$, separating the energy spectrum from the spatial density profile. Antideuteron formation is handled with a coalescence model based on the Wigner formalism, calibrated on collider data. The upper bounds come from a likelihood-ratio analysis of the antiproton data that treats the Galactic halo height as a nuisance parameter and uses a covariance matrix encoding experimental and theoretical uncertainties.","core_discovery":"The paper's central claim is that AMS-02 antiproton data set strong upper bounds on the local density of Galactic primordial black holes with a lognormal mass distribution, and that these bounds are comparable to or slightly stronger than constraints from other messengers. For the critical mass $M^* \\simeq 5\\times 10^{14}$ g, the mass of a PBH evaporating completely today, the maximum allowed density is $3.5\\times10^{-12}$, $2.2\\times10^{-11}$, $4.4\\times10^{-11}$, and $8.1\\times10^{-11}$ GeV cm$^{-3}$ for widths $\\sigma = 0.1, 0.5, 1, 2$, respectively. Translated into a fraction of the local dark-matter density, the strongest limit is $f_{\\mathrm{PBH}} < 10^{-11}$ at $\\sigma = 0.1$ and $\\mu_c = M^*$. The same calculation yields a universal spectral shape for every primary PBH antiproton and antideuteron flux, so the antiproton bound fixes the normalization of the antideuteron flux as $\\mu_c$ and $\\sigma$ vary. From this the authors conclude that a future antideuteron detection would clearly be a signal of new physics, but one that PBH evaporation alone cannot fully explain.","pith_inferences":["The universal spectral shape implies that a template search could extract the PBH component from low-energy antiproton data without fitting the full lognormal parameters, because only the normalization depends on the mass-distribution combination.","The same universality predicts exactly where to look: the antiproton peak sits near 1-2 GeV and the antideuteron peak near a few hundred MeV/n, so future sub-GeV measurements at solar minimum would either sharpen the bound or reveal residuals.","If the memory-burden effect slows evaporation as recent proposals suggest, the semi-classical assumption would break before the final stage; the predicted fluxes and bounds would shift, and an antideuteron detection would need a different interpretation."],"forward_implications":["At $\\mu_c = M^* \\simeq 5\\times10^{14}$ g and $\\sigma=0.1$, the AMS-02 antiproton data exclude PBHs as the whole of local dark matter, with $f_{\\mathrm{PBH}} < 10^{-11}$.","Because all primary fluxes share one spectral shape, the antiproton bounds translate directly into a maximum antideuteron flux: it exceeds the secondary background below a few GeV/n but stays near the sensitivity of upcoming detectors.","If an antideuteron is measured, the paper's conclusion is that it is a clear new-physics signal whose PBH component cannot account for the full event rate.","Fixing the lognormal distribution, the same bounds imply a local PBH explosion rate of at most about $7\\times10^{-5}$ pc$^{-3}$ yr$^{-1}$, an improvement of more than two orders of magnitude over an earlier estimate.","Broader lognormal widths weaken the density bound but extend the reach to larger critical masses, probing the asteroid-mass region where PBHs could conceivably be all of the dark matter."],"supporting_citations":[{"why":"Hawking's original prediction that black holes radiate thermally is the physical basis of the evaporation calculation.","marker":"[16]"},{"why":"It establishes particle creation by black holes and the spectrum used for the emitted fundamental particles.","marker":"[17]"},{"why":"It provides the quark and gluon jet emission from PBHs, whose hadronization yields the antinuclei source spectra.","marker":"[56]"},{"why":"It supplies the instantaneous spectra of directly emitted particles from PBHs that feed the convolution in the source term.","marker":"[57]"},{"why":"It is the earlier antideuteron-from-exotic-sources calculation that this study updates with newer data and propagation models.","marker":"[45]"},{"why":"It provides the first AMS-02 antiproton data release used in the statistical analysis and in comparisons.","marker":"[48]"},{"why":"It provides the larger, more recent AMS-02 antiproton dataset that sets the main upper bounds on the PBH density.","marker":"[49]"},{"why":"It supplies the likelihood-ratio analysis, propagation setup, and halo-size nuisance treatment adapted here to PBHs.","marker":"[91]"},{"why":"It defines the secondary antiproton background and its covariance matrix against which the PBH contribution is constrained.","marker":"[96]"},{"why":"It provides the Wigner-coalescence model for antideuteron formation used in the antideuteron flux predictions.","marker":"[89]"}],"fun_headline_variants":["Antiproton flux tightly bounds primordial black hole density","AMS-02 antiprotons cap dark matter from evaporating black holes","Antideuteron sighting would mean new physics, not just PBHs","No room for black hole dark matter in antiproton data"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that black holes evaporate exactly as in the semi-classical Hawking calculation, with no quantum memory-burden back-reaction slowing the final stage; if evaporation is slower, the mass evolution, source spectra, and all derived density bounds would change.","fun_headline_variants_meta":{"raw":{"variants":["Antiproton flux tightly bounds primordial black hole density","AMS-02 antiprotons cap dark matter from evaporating black holes","Antideuteron sighting would mean new physics, not just PBHs","No room for black hole dark matter in antiproton data"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000711,"raw_usage":{"total_tokens":3256,"prompt_tokens":1056,"completion_tokens":2200,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":672,"completion_tokens_details":{"reasoning_tokens":2126}},"tokens_in":672,"tokens_out":2200,"duration_ms":14550,"temperature":1.0,"reasoning_tokens":2126,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:23:35.076003+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Detect antideuterons: the paper predicts a maximum antideuteron flux obtained by saturating the antiproton bound, so an observed rate above that maximum excludes PBH evaporation as the source, while a clear non-detection at the predicted peak energy would tighten the density bound. A second check is whether PBHs near $5\\times10^{14}$ g still evaporate on the Universe's timescale; if memory-burden effects make them live longer, the predicted fluxes disappear.","supporting_citations":[{"cited_title":"Antideuterons in cosmic rays: sources and discovery potential","cited_arxiv_id":"1610.00699","evidence_quote":"It is the earlier antideuteron-from-exotic-sources calculation that this study updates with newer data and propagation models."}],"review_version":1}