{"id":"f986e4da-d8bb-4d1f-bb5e-26d48f412e18","arxiv_id":"2604.18007","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"CMB data limits the s-wave annihilation cross section of thermal dark matter particles to ≲ 10^{-30} cm³/s scaled by PBH fraction and mass for PBHs heavier than ~10^{-10} solar masses.","lead":"The paper uses a statistical analysis of cosmic microwave background data to set limits on mixed dark matter models where primordial black holes cluster self-annihilating particles into dense halos. This restricts the allowed annihilation rates for such particles when PBHs exceed a certain mass threshold.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.3","headline":"Central limits depend on unvalidated dense halo profiles from Parts I/II","rationale":"The reader's identification of the halo modeling assumption as the weakest link is correct and remains the single load-bearing step even after the full text becomes available; all subsequent CMB statistical steps inherit the uncertainty from that prior modeling. No other internal inconsistency (e.g., in the likelihood construction or parameter marginalization) rises to the same level of centrality.","tokens_in":1766,"tokens_out":428,"duration_ms":17409,"concrete_test":"Take the benchmark parameters f_BH=10^{-6}, m_PBH=10^{-9} M_⊙, mχ=100 GeV from the paper; recompute the halo density profile using only the analytic matching and turnaround conditions stated in §2 of Part I (without importing the numerical results); insert the resulting annihilation luminosity into the CMB Boltzmann solver used in §4 of Part III; if the Δχ² relative to Planck 2018 TT+EE+lowE shifts by more than 2σ from the published value, the quoted cross-section limit is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline constraint on s-wave <σv> ≲ 10^{-30} cm³/s (mχ/100 GeV) (f_BH/10^{-6})^{-3} is obtained by folding an enhanced annihilation rate—arising from the ρ(r) ∝ r^{-9/4} or steeper cusps built during radiation domination—into the standard CMB energy-injection calculation (ionization, heating, and distortion of the recombination history). This rate is taken directly from the halo models of the prior papers in the series; no independent derivation or robustness test of those profiles appears in Part III. If the assumed halo overdensity or the matching to the background DM density at the turnaround radius is off by even a factor of a few, the injected energy per PBH scales linearly with that factor and the derived upper bound on <σv> shifts proportionally, rendering the quoted numerical limit unreliable.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript, as Part III in a series, examines a mixed dark matter scenario where a fraction f_BH of primordial black holes (PBHs) coexist with thermally produced self-annihilating particles. It models dense halos of particles around PBHs formed during radiation domination and derives constraints via a full statistical analysis of CMB data. The central result limits the s-wave annihilation cross section to ≲ 10^{-30} cm³/s (m_χ/100 GeV) (f_BH/10^{-6})^{-3} for PBHs heavier than ∼10^{-10} M_⊙, which can also constrain PBH abundances in this range; lighter (asteroid-mass) PBHs are compatible with standard thermal relics. Implications for Subaru-HSC microlensing events are briefly discussed.","tokens_in":1990,"tokens_out":834,"duration_ms":43842,"significance":"If the halo clustering model holds, the work provides a valuable bridge between PBH and particle dark matter by showing how even tiny f_BH can exclude large regions of thermal WIMP parameter space through enhanced annihilation in dense cusps. The use of a full statistical CMB analysis (rather than approximate energy-injection estimates) is a clear strength, as is the explicit scaling of limits with f_BH and the timely link to Subaru-HSC data. This could meaningfully restrict viable mixed DM models and motivate further multi-probe studies.","major_comments":[{"comment":"The headline bound on <σv> is obtained by folding the enhanced annihilation rate from the assumed halo profiles (ρ(r) ∝ r^{-9/4} or steeper, built during radiation domination) into the standard CMB energy-injection calculation. This rate is taken directly from the halo models of Parts I and II with no independent derivation, robustness test, or sensitivity analysis to variations in overdensity or turnaround-radius matching appearing in Part III. Because the injected energy per PBH scales linearly with the central density, any systematic offset in those profiles shifts the quoted <σv> limit proportionally (see the scaling in the abstract and the results section).","section":"Introduction and results section"},{"comment":"The distinction that PBHs heavier than ∼10^{-10} M_⊙ impose strong limits while asteroid-mass or lighter PBHs 'live in perfect peace' with the particles rests on the halo formation and survival assumptions carried over from prior papers. A quantitative demonstration (e.g., via an explicit calculation or figure showing suppression of the annihilation rate below this mass threshold) is needed to support the claim that the effect vanishes for lighter PBHs.","section":"Results and discussion"},{"comment":"The full statistical CMB analysis is presented as the basis for the numerical limits, yet the manuscript does not detail the specific likelihoods (Planck or otherwise), foreground treatment, or how the PBH-induced ionization/heating term is implemented in the recombination code. Without these, it is difficult to assess whether the f_BH^{-3} scaling and the numerical prefactor are robust to reasonable variations in the analysis pipeline.","section":"CMB analysis section"}],"minor_comments":[{"comment":"The abstract states that the limits 'can also be turned into constraints on PBHs' but does not show the corresponding f_BH upper bounds as a function of mass; adding a brief table or plot would improve clarity.","section":"Abstract"},{"comment":"Notation for f_BH, m_χ, and M_⊙ is introduced without a dedicated symbols table or consistent definition on first use in the main text.","section":"Introduction"}],"recommendation":"major_revision","confidential_remarks":"The central results depend entirely on the authors' own prior halo-modeling papers; the editor should consider whether an independent validation or external cross-check of the ρ(r) profiles is required before publication. The arXiv identifier 2604.18007 appears anomalously dated relative to current numbering conventions."},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the positive assessment of our work and for the detailed comments that help improve the clarity and robustness of the manuscript. We address each major comment point-by-point below, indicating the revisions made to the manuscript.","responses":[{"response":"We thank the referee for highlighting this point. As this is Part III of the series, the detailed derivation of the halo profiles is presented in Parts I and II. In the current manuscript, we have added a new subsection in the introduction summarizing the key features of the halo models, including the ρ ∝ r^{-9/4} profile and the assumptions on overdensity and turnaround radius. Additionally, we have performed a sensitivity analysis by varying the central density by factors of 2 and 0.5, showing that the resulting limits on σv shift proportionally as expected, but the overall conclusions remain unchanged. This is now included in the revised results section.","revision_made":"partial","referee_comment":"[Introduction and results section] The headline bound on σv is obtained by folding the enhanced annihilation rate from the assumed halo profiles (ρ(r) ∝ r^{-9/4} or steeper, built during radiation domination) into the standard CMB energy-injection calculation. This rate is taken directly from the halo models of Parts I and II with no independent derivation, robustness test, or sensitivity analysis to variations in overdensity or turnaround-radius matching appearing in Part III. Because the injected energy per PBH scales linearly with the central density, any systematic offset in those profiles shifts the quoted σv limit proportionally (see the scaling in the abstract and the results section)."},{"response":"The mass threshold of ∼10^{-10} M_⊙ corresponds to the point where the PBH-induced halos form sufficiently early and densely during radiation domination to survive until recombination without significant disruption. For lighter PBHs, the smaller mass leads to earlier formation but also to halos that are more susceptible to tidal disruption or have lower peak densities due to the matching conditions at turnaround. To address the request for quantitative demonstration, we have added an explicit calculation of the annihilation rate as a function of PBH mass in the revised manuscript, including a new figure that shows the effective suppression factor for masses below 10^{-10} M_⊙, confirming that the annihilation signal becomes negligible compared to the standard thermal relic case.","revision_made":"yes","referee_comment":"[Results and discussion] The distinction that PBHs heavier than ∼10^{-10} M_⊙ impose strong limits while asteroid-mass or lighter PBHs 'live in perfect peace' with the particles rests on the halo formation and survival assumptions carried over from prior papers. A quantitative demonstration (e.g., via an explicit calculation or figure showing suppression of the annihilation rate below this mass threshold) is needed to support the claim that the effect vanishes for lighter PBHs."},{"response":"We apologize for the lack of detail in the original submission. In the revised version, we have expanded the CMB analysis section to specify that we use the Planck 2018 likelihoods (TT, TE, EE + lowE + lensing), with standard foreground marginalization as implemented in the Planck analysis. The PBH-induced energy injection is incorporated by modifying the ionization and heating rates in the recombination code (using a customized version of CosmoRec), where the additional term is proportional to the PBH number density times the per-PBH annihilation luminosity. We have verified that the f_BH^{-3} scaling arises directly from the combination of PBH density (∝ f_BH) and the halo annihilation rate (∝ f_BH^2 from the two-body process in the dense halo), and the numerical prefactor is robust to small variations in the pipeline as tested with alternative likelihood combinations.","revision_made":"yes","referee_comment":"[CMB analysis section] The full statistical CMB analysis is presented as the basis for the numerical limits, yet the manuscript does not detail the specific likelihoods (Planck or otherwise), foreground treatment, or how the PBH-induced ionization/heating term is implemented in the recombination code. Without these, it is difficult to assess whether the f_BH^{-3} scaling and the numerical prefactor are robust to reasonable variations in the analysis pipeline."}],"tokens_in":1650,"tokens_out":904,"duration_ms":41713,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The one thing to know is that this Part III folds the authors' earlier clustering calculation into a standard CMB energy-injection analysis and comes out with a concrete upper limit on the s-wave annihilation cross section of roughly 10^{-30} cm³/s times mass and f_BH factors, for PBHs above 10^{-10} solar masses. That specific numerical result is not in the prior literature they cite, so it is a genuine extension of their program. The paper also notes that asteroid-mass or lighter PBHs evade the constraint and briefly links the result to Subaru-HSC microlensing events. Those are the useful outputs. The analysis itself is a full statistical treatment of CMB data rather than a back-of-the-envelope estimate, which is the right approach for this kind of limit. The authors are clear about how the bound tightens the allowed space for thermal particle DM when even a small PBH fraction is present. That part is straightforward and worth having on record for people working in this narrow corner of mixed DM models. The soft spot is exactly the one flagged in the stress test. The enhanced annihilation rate comes entirely from the ρ(r) ∝ r^{-9/4} cusps and the matching to background density that were derived in the previous papers; Part III does not re-derive those profiles, test their sensitivity, or provide independent benchmarks. If the central overdensity is off by a factor of a few, the injected energy per PBH shifts proportionally and the quoted cross-section limit moves with it. That dependence is load-bearing for the headline number, not peripheral. The paper treats the halo model as an input rather than revisiting it, which is acceptable for a continuation but leaves the result conditional. This is work for specialists already following the series or looking for updated numbers in PBH-plus-annihilating-DM scenarios. A reader who wants to use the bound should first check the halo assumptions in Parts I and II. It shows honest engagement with the literature and a concrete calculation, so it deserves a serious referee who can examine the propagation of the halo uncertainties into the final limits.","headline":"The CMB bounds here are new numbers but scale linearly with the dense halo profiles taken from the authors' own Parts I and II.","tokens_in":2463,"tokens_out":493,"would_cite":false,"duration_ms":24072,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"Even a tiny fraction of primordial black holes can force co-existing thermal dark matter particles to annihilate far more slowly than usual, as shown by CMB data analysis.","keywords":["primordial black holes","dark matter","cosmic microwave background","annihilation cross section","mixed dark matter","dense halos","thermal relics"],"falsifier":"A CMB power spectrum measurement showing no excess energy injection from annihilation at the predicted level for a confirmed population of PBHs above 10^{-10} solar masses with fraction 10^{-6} and a larger annihilation cross section would falsify the derived limits.","tokens_in":2669,"feed_emoji":"🌌","tokens_out":797,"duration_ms":42074,"temperature":0.7,"pith_summary":"The paper examines mixed dark matter models containing both primordial black holes and self-annihilating particles. It focuses on the dense halos of particles that accumulate around the black holes during the early radiation-dominated era. These halos cause extra annihilation that releases energy into the plasma, altering the cosmic microwave background. A full statistical analysis of CMB observations then translates this effect into tight upper limits on the annihilation cross section when black holes exceed a mass threshold. Readers would care because the result shows that the presence of black holes and particle dark matter cannot be treated independently; their interaction through these halos shrinks the allowed parameter space for standard thermal dark matter candidates.","feed_headline":"Tiny PBH fraction caps DM annihilation at 10^{-30} cm³/s","feed_subtitle":"CMB analysis shows dense halos around heavy primordial black holes restrict thermal particle dark matter parameters.","key_machinery":"Dense halos of self-annihilating dark matter particles that form around primordial black holes and inject annihilation energy into the early universe, modifying CMB anisotropies.","core_discovery":"In a mixed dark matter scenario, primordial black holes are surrounded by extremely dense halos of self-annihilating particles built up during radiation domination. A full statistical analysis of CMB data shows that for PBHs heavier than about 10^{-10} solar masses, even a fraction f_BH as small as 10^{-6} restricts the s-wave annihilation cross section to values ≲ 10^{-30} cm³/s (m_χ/100 GeV) (f_BH/10^{-6})^{-3}. Lighter PBHs in the asteroid mass range or below impose no such restrictions and can coexist with the particles without tension. The analysis also briefly considers implications for tentative Subaru-HSC microlensing events interpreted as PBHs.","pith_inferences":["Future CMB experiments with higher sensitivity could extend these limits to smaller PBH fractions or even lower cross sections.","The same halo mechanism may produce observable gamma-ray signals from individual nearby PBHs, providing an independent test.","Model builders should treat PBH and particle dark matter production as potentially coupled rather than separate processes.","Direct searches for PBHs should prioritize mass ranges where the halo effect does or does not apply to avoid or exploit these constraints."],"forward_implications":["The allowed parameter space for thermal relic dark matter particles shrinks when PBHs heavier than 10^{-10} solar masses are present.","The abundance of such PBHs is itself constrained by the requirement that the annihilation cross section not violate CMB bounds.","Asteroid-mass or lighter PBHs remain fully compatible with standard thermal dark matter annihilation rates.","Any interpretation of microlensing events as PBHs in the relevant mass range would require strongly suppressed annihilation cross sections for co-existing particles."],"fun_headline_variants":["Heavy PBH fractions cap DM annihilation at 10^{-30} cm³/s","CMB data shows dense halos limit DM cross section to 10^{-30} cm³/s","Tiny PBH share restricts annihilation below 10^{-30} cm³/s for heavy PBHs","Asteroid mass PBHs compatible with DM particles per CMB analysis"],"cache_read_input_tokens":64,"weakest_assumption_plain":"That primordial black holes develop extremely dense halos of self-annihilating dark matter particles during radiation domination.","fun_headline_variants_meta":{"raw":{"variants":["Heavy PBH fractions cap DM annihilation at 10^{-30} cm³/s","CMB data shows dense halos limit DM cross section to 10^{-30} cm³/s","Tiny PBH share restricts annihilation below 10^{-30} cm³/s for heavy PBHs","Asteroid mass PBHs compatible with DM particles per CMB analysis"]},"model":"grok-4.3","cost_usd":0.012334,"raw_usage":{"total_tokens":5409,"prompt_tokens":736,"num_sources_used":0,"completion_tokens":88,"cost_in_usd_ticks":123337000,"prompt_tokens_details":{"text_tokens":736,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":4585,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":736,"tokens_out":88,"duration_ms":60713,"temperature":1.0,"reasoning_tokens":4585,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-05-10T04:07:10.789751+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A CMB power spectrum measurement showing no excess energy injection from annihilation at the predicted level for a confirmed population of PBHs above 10^{-10} solar masses with fraction 10^{-6} and a larger annihilation cross section would falsify the derived limits.","supporting_citations":[],"review_version":1}