{"id":"1f5e12c2-28ad-4a8d-9472-d87cf7b3e1c7","arxiv_id":"2502.00914","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":4,"one_line_summary":"Poisson shot noise from stellar-mass primordial black holes can boost early-universe ultradense dark matter halo formation, with predictions highly sensitive to the assumed shape of the small-scale isocurvature power spectrum.","lead":"This paper calculates how the random, discrete positions of stellar-mass primordial black holes add extra small-scale power to the dark matter density field, and how that extra power changes the predicted abundance of ultradense dark matter halos in the early universe. It finds that heavier black holes and a milder high-wavenumber cutoff strongly boost the number of these halos, with the preferred dark matter composition depending on the assumed cutoff shape.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed preference for multi- vs single-component dark matter rests on comparing UDMH mass fractions to PBH abundance upper limits, which are different quantities; the inference is unsupported as written.","rationale":"The paper's analytic framework is coherent, and the qualitative statement that PBH discreteness enhances small-scale power is plausible and consistent with prior work. The most load-bearing issue is not the phenomenological cutoff in Eq. (25) by itself, since the authors transparently treat n as a free parameter and scan it. The decisive move is the comparison of the predicted UDMH mass fraction to observational upper limits on the PBH abundance: df/dlogM is the fraction of dark matter in halos of mass M, while the quoted constraints bound the fraction of dark matter in PBHs of mass MPBH. Because M_halo/MPBH can be large, the plotted curves can exceed the PBH exclusion regions even when fPBH is small, and no extra dark matter component is implied. This is an internal mismatch between the predicted quantity and the observational bound used to evaluate it, not a disagreement with consensus. The reader's conditional verdict remains appropriate, but for this more specific reason; the formal weakest_assumption of the reader centered on the cutoff in Eq. (25), whereas the comparison issue is at least as important and was only noted in the reader's rationale.","tokens_in":15571,"tokens_out":28601,"duration_ms":310129,"concrete_test":"Replot the high-fPBH curves of Fig. 2 as the implied PBH fraction f_implied = (MPBH / M_halo) * df/dlogM, assuming one seed PBH per UDMH of mass M_halo, and overlay the same OGLE/LVK/UFD upper limits. If the converted f_implied lies below the limits for the curves currently shown as violating them, the claimed violation and the multi-component conclusion are an artifact of comparing df/dlogM with fPBH bounds. If the converted curves still cross the exclusion regions, the comparison can be salvaged.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"In Figs. 2-4 the differential mass fraction df/dlogM of UDMHs is plotted and directly compared with shaded upper limits on PBH abundance from OGLE, LVK, and UFD disruption. These constraints bound the fraction of dark matter in stellar-mass black holes, not the mass fraction in extended ultradense halos. In the model, a PBH can seed a UDMH whose mass M is much larger than MPBH; the plotted curves reach ~10^6 solar masses even for MPBH = 1-100 solar masses. Thus the predicted df/dlogM can lie above the PBH exclusion region while the actual fPBH is far below the bound, and exceeding the shaded regions does not by itself imply that a second dark matter component is needed. The paper's conclusion that low n favors multi-component dark matter while high n favors single-component dark matter is drawn precisely from where df/dlogM sits relative to these PBH-only limits. That comparison is apples-to-oranges, so the multi-component inference is not supported as stated. The excursion-set computation may remain useful as a model exploration, but the observational interpretation and the multi/single-component claim must be re-expressed in terms of the same constrained quantity.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the effect of Poisson noise from stellar-mass primordial black holes (PBHs) on the formation of ultradense dark matter halos (UDMHs) during the radiation-dominated era. It constructs a modified power spectrum that adds a PBH-induced isocurvature term with an exponential cutoff to the adiabatic CDM power spectrum, and then uses excursion-set theory with several mass-function prescriptions (PS, ST, DP1, DP2) to compute differential mass fractions of UDMHs. The authors vary the PBH mass, the PBH fraction, and a suppression parameter n, and compare the resulting mass functions with observational upper limits on PBH abundance from OGLE, LVK, and ultrafaint dwarf disruption. The central claims are that heavier PBHs shift UDMH mass functions to higher masses, lower n boosts UDMH abundance and favors multi-component dark matter, higher n aligns with single-component dark matter, and the DP1/DP2 mass functions predict higher abundances than Press-Schechter.","tokens_in":15768,"tokens_out":9667,"duration_ms":110108,"significance":"If the calculation were robust, the paper would be a useful exploration of an interesting and understudied effect: the discreteness of PBHs can strongly affect small-scale density fluctuations and thus the formation of compact dark matter structures. The paper is clearly written in its main equations and covers a reasonable parameter space, and the use of several halo mass functions is a helpful comparison. However, the advertised observational conclusions are currently not supported by the comparison made in Figs. 2-4, and the quantitative results are governed by an empirically unspecified cutoff in the isocurvature power spectrum. The paper would be strengthened by re-expressing the comparison in terms of a single constrained quantity and by treating the suppression parameter as a free knob in a conditional model exploration rather than as the basis for claims about which dark matter scenario is favored.","major_comments":[{"comment":"The central observational inference is based on comparing the differential mass fraction of UDMHs, df/dlogM, with shaded upper limits on the PBH abundance f_PBH from OGLE, LVK, and UFD disruption. These are different quantities: the shaded regions constrain the fraction of dark matter in stellar-mass black holes, while the curves describe the mass fraction in extended ultradense halos that can be much more massive than the seeding PBH (the plotted curves extend to ~10^6 M_sun for M_PBH = 1-100 M_sun). A single PBH can therefore produce a UDMH whose mass fraction lies above a PBH exclusion region while f_PBH itself is far below the bound. Consequently, the statement in Section 3 that the curves 'clearly violate' regions excluded by PBH constraints and thereby 'lend compelling support to a multi-component dark matter scenario' is unsupported as written. The same issue underlies the conclusion that low n favors multi-component and high n favors single-component dark matter. The comparison should be recast in terms of one constrained quantity, for example by converting the predicted UDMH abundances into the relevant observational observable or by explicitly stating that the curves are model predictions not directly bounded by PBH abundance limits.","section":"Section 3, Figs. 2-4"},{"comment":"The quantitative results are governed by the assumed exponential cutoff P_iso(k) = Q(k) exp(-[(k-k')/sigma]^n / 2) with k' ~ 10^2 Mpc^-1, sigma ~ 4.5 x 10^2 Mpc^-1, and free index n, together with the adopted amplitude A_iso = 3.2 x 10^-12 f_PBH (M_PBH/30 M_sun). The paper itself acknowledges that the microphysical origin of this suppression is model-dependent. Because P_iso is the dominant small-scale contribution, the conclusions that lower n boosts UDMH formation and that n = 2 or 3 brings the results into agreement with single-component dark matter are direct consequences of the chosen shape and parameters rather than robust predictions. The authors should either derive these parameters from a microphysical model, constrain them with independent data, show a sensitivity analysis over k' and sigma, or clearly frame the results as a conditional model exploration. Without such an addition, statements in the abstract and conclusions about favoring multi- versus single-component dark matter are overreaching.","section":"Section 2.2, Eqs. (25)-(27)"},{"comment":"There is a notational inconsistency that affects reproducibility: Eq. (23) defines P(k) as the power spectrum of primordial curvature perturbations, while Fig. 1 is labeled as the modified matter power spectrum and Eq. (7) requires the matter power spectrum P(k,a). The connection is made through Eq. (9), but the text and figure should use consistent notation so that a reader can reproduce sigma^2(M,a) without guessing whether Eq. (23) is P_zeta(k) or P_matter(k).","section":"Section 2.2 and Fig. 1"}],"minor_comments":[{"comment":"The caption lists f_PBH = 1, 0.1, 0.01, 0.0001, and 0.0001; the last value is presumably 10^-5 and should be corrected.","section":"Fig. 1 caption"},{"comment":"The lower-left panel appears to have a duplicated axis label 'MPBH=10M'; one of the two labels should be removed.","section":"Fig. 4"},{"comment":"The shading colors are described inconsistently: the main text refers to blue, red, and green shaded regions for OGLE, LVK, and UFD, while the captions list orange, wheat, and cyan with different hatching.","section":"Figs. 2-4 captions and text"},{"comment":"The LVK citation in the figures is malformed as '?R. e. a. Abbott et al. 2022' and should be fixed to a standard author-year reference.","section":"Figs. 2-4 captions"},{"comment":"Equation (12) appears to have a typographical error in its prefactor: as printed it reads as a product of terms rather than a normalized first-crossing distribution with a denominator.","section":"Eq. (12)"}],"recommendation":"major_revision","confidential_remarks":"The main advertised conclusion about multi-component versus single-component dark matter rests on comparing UDMH mass fractions to PBH abundance limits, which are different quantities; this needs to be fixed or substantially reframed. The suppression parameter n is effectively a free knob, and the paper's qualitative conclusions track it directly. I would encourage the editor to require a revised version that either derives or constrains the cutoff parameters, or explicitly presents the work as a conditional model study with the observational language softened."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the useful part: this paper does a careful analytic scan of how Poisson shot noise from stellar-mass PBHs alters the UDMH mass function, using the excursion-set formalism with several mass functions (PS, ST, DP1, DP2). It confirms the known result that heavier PBHs (lower number density) give stronger shot noise and push UDMH formation to higher masses, and it shows the suppression parameter n and f_PBH matter. The DP1/DP2 mass functions, which include angular momentum and dynamical friction, predict higher abundances than PS, which is expected and worth stating. The algebra is mostly internally consistent, and the paper is honest that the exponential cutoff in Eq. (25) is phenomenological.\n\nNow the problems. The central claim about multi- vs single-component dark matter is not supported. Figures 2-4 plot df/dlogM of UDMHs and compare it to shaded upper limits on f_PBH from OGLE, LVK, UFD disruption. Those bounds constrain the fraction of dark matter in stellar-mass black holes, not the mass fraction in extended UDMHs. A single PBH can seed a UDMH much more massive than M_PBH (up to ~10^6 M_sun), so the predicted df/dlogM can exceed the PBH exclusion region while actual f_PBH is far below the bound. The paper's conclusion that low n favors multi-component and high n favors single-component is drawn precisely from where df/dlogM sits relative to these PBH-only limits. That is apples-to-oranges. To make that inference, they need to compare with constraints on UDMH abundance or compact halo mass fraction, not PBH abundance.\n\nSecond, the quantitative results are controlled by the hand-chosen cutoff in Eq. (25): k' ~ 10^2 Mpc^-1, sigma ~ 4.5x10^2 Mpc^-1, and free index n. The paper scans n and reads off conclusions about which dark matter scenario is favored. Since the cutoff is not derived, the results are a model exploration, not a prediction. That's fine as a conditional contribution, but the abstract and conclusions overstate the robustness.\n\nThird, the novelty is modest. The Poisson-noise power spectrum, the UDMH mass function framework, and the DP mass functions all come from prior work, much of it by the same authors. The paper's contribution is a parameter scan combining them. I don't see a new mechanism or a new result that is independent of the assumptions.\n\nWho is this for? Someone working in PBH substructure who wants a quick map of how UDMH abundance depends on M_PBH, f_PBH, and the shape of the isocurvature cutoff. That reader will find the plots useful, but should treat the observational comparisons with caution. It deserves a serious referee because the topic is relevant and the formalism is sound enough to engage with, but the referee should ask for a re-expression of the constraints and a clear separation between derived and assumed inputs.\n\nRecommendation: send to peer review, but with a request for major revision.","headline":"A clean parameter scan confirming known Poisson-noise effects on UDMHs, but the multi-component dark matter claim rests on comparing UDMH mass fractions to PBH abundance bounds.","tokens_in":16383,"tokens_out":4818,"would_cite":false,"duration_ms":36941,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper argues that the discrete, Poisson-distributed positions of stellar-mass primordial black holes add an isocurvature term to the matter power spectrum that strongly boosts the formation of ultradense dark matter halos, with the…","keywords":["primordial black holes","Poisson noise","ultradense dark matter halos","matter power spectrum","halo mass function","dark matter","small-scale structure","isocurvature perturbations"],"falsifier":"A measurement of the small-scale matter power spectrum at $k\\sim10^2$–$10^4\\,\\mathrm{Mpc}^{-1}$ (for example from Lyman-$\\alpha$ forest or 21-cm observations) that shows no isocurvature bump of the predicted height, or a high-resolution N-body simulation of stellar-mass PBH shot noise whose UDMH mass function does not shift toward higher masses with increasing $M_{\\mathrm{PBH}}$, would settle against the central claim.","tokens_in":15283,"feed_emoji":"🕳️","tokens_out":6874,"duration_ms":58562,"temperature":0.7,"pith_summary":"This paper tries to establish that the discrete, random spatial distribution of stellar-mass primordial black holes—their Poisson noise—leaves a measurable imprint on the formation of ultradense dark matter halos in the radiation-dominated era. It constructs a modified matter power spectrum that adds an isocurvature term from PBH shot noise to the usual adiabatic spectrum, and shows that this term strongly enhances small-scale density fluctuations, especially for more massive PBHs and higher PBH fractions. The authors find that the resulting UDMH mass function shifts toward higher masses as the PBH mass grows, and that the steepness of the small-scale cutoff, parameterized by $n$, controls whether the predictions favor multi-component dark matter (low $n$) or a single-component stellar-mass-PBH dark matter (high $n$). A sympathetic reader would care because these ultradense structures could be observable probes of PBH dark matter that are distinct from the usual continuous-fluid approximation.","feed_headline":"Stellar-mass PBH noise boosts ultradense dark matter halos","feed_subtitle":"Heavier black holes shift halo masses up; steeper cutoffs favor a single-component dark matter picture.","key_machinery":"The central object is the modified linear matter power spectrum $P(k)=P_{\\mathrm{ad}}(k)+P_{\\mathrm{iso}}(k)$ of Eq. (23), where $P_{\\mathrm{ad}}$ is the nearly scale-invariant adiabatic spectrum and $P_{\\mathrm{iso}}$ is the PBH Poisson-noise (isocurvature) contribution with amplitude proportional to $f_{\\mathrm{PBH}}M_{\\mathrm{PBH}}$ and an exponential cutoff controlled by the parameter $n$. This spectrum converts into halo abundances through the excursion-set peak-height variable $\\nu=\\delta_c/\\sigma(M)$ and the multiplicity functions PS, ST, DP1, DP2; the DP-family barriers include angular momentum and dynamical friction, which raise the predicted UDMH abundance.","core_discovery":"The paper's central claim is that the shot noise from a discrete population of stellar-mass PBHs cannot be ignored: it contributes an isocurvature term $P_{\\mathrm{iso}}(k) = Q(k)\\exp[-(k-k')^n/(2\\sigma^n)]$ to the matter power spectrum, with amplitude $A_{\\mathrm{iso}} = 3.2\\times10^{-12} f_{\\mathrm{PBH}}(M_{\\mathrm{PBH}}/30M_\\odot)$ at the pivot scale. This Poisson-noise term enhances small-scale power by many orders of magnitude, and when fed through excursion-set halo mass functions it shifts the differential mass function of UDMHs toward higher masses as $M_{\\mathrm{PBH}}$ increases. The steepness index $n$ of the exponential cutoff controls how much of this enhancement survives: low $n$ leaves a broad small-scale boost that favors multi-component dark matter, while high $n$ suppresses low-mass UDMHs and brings predictions into line with a single-component PBH dark matter scenario. The paper also claims that mass functions including angular momentum and dynamical friction (DP1 and DP2) yield more UDMHs than the Press-Schechter form.","pith_inferences":["Because the exponential cutoff parameters ($k'$, $\\sigma$, $n$) are phenomenological, the same analytic machinery could be re-run with a physically motivated damping spectrum to see whether the multi-component versus single-component preference survives.","The Poisson-noise contribution is generic to any discrete compact-object dark matter component, so the predicted UDMH boost should apply to other massive compact halo objects, not only to PBHs.","If UDMHs are as abundant as claimed at low $n$, their annihilation or lensing signatures could be searched for in gamma-ray and microlensing data even when the PBH fraction is too small to be seen directly.","The degeneracy among $n$, $M_{\\mathrm{PBH}}$, and $f_{\\mathrm{PBH}}$ means that observational upper limits on UDMH abundance translate into joint constraints, not independent bounds on the PBH fraction."],"forward_implications":["For fixed $f_{\\mathrm{PBH}}$, increasing $M_{\\mathrm{PBH}}$ from $1\\,M_\\odot$ to $100\\,M_\\odot$ shifts the differential UDMH mass function to higher masses because fewer, heavier PBHs produce stronger shot noise.","Lower values of the suppression index $n$ leave more small-scale power and boost the UDMH abundance, a regime in which the predicted halo counts exceed PBH observational bounds and favor multi-component dark matter.","Higher $n$ values damp low-mass UDMH formation and bring the predicted distributions into closer agreement with OGLE, LVK, and ultra-faint-dwarf constraints, supporting a single-component PBH dark matter scenario for lighter PBHs.","The DP1 and DP2 mass functions, which include angular momentum and dynamical friction, consistently produce higher UDMH abundances than the Press-Schechter form, especially in the high-mass tail.","The modified power spectrum's small-scale amplitude is enhanced by at least seven orders of magnitude relative to large scales, making the Poisson-noise effect potentially accessible to small-scale structure probes."],"supporting_citations":[{"why":"Supplies the Poisson-noise power spectrum $P_p \\simeq n_{\\mathrm{PBH}}^{-1}$ for a discrete PBH distribution.","marker":"(N. Afshordi et al. 2003)"},{"why":"Provides the isocurvature amplitude scaling and the exponential-cutoff form adopted for $P_{\\mathrm{iso}}(k)$.","marker":"(J.-O. Gong & N. Kitajima 2017)"},{"why":"Establishes the UDMH formation picture and the first-crossing multiplicity function approximation used to compute halo abundances.","marker":"(M. S. Delos & J. Silk 2023)"},{"why":"Gives the ellipsoidal-collapse moving barrier and the Sheth-Tormen multiplicity function.","marker":"(R. K. Sheth et al. 2001)"},{"why":"Adds angular momentum and cosmological-constant effects, defining the DP1 mass function.","marker":"(A. Del Popolo 2006)"},{"why":"Adds dynamical friction to the barrier, defining the DP2 mass function that yields higher UDMH abundances.","marker":"(A. Del Popolo et al. 2017)"},{"why":"Sets the adiabatic spectrum amplitude $A_s$ and tilt $n_s$ used as the baseline.","marker":"(Planck Collaboration et al. 2020)"},{"why":"Provides the OGLE microlensing upper limits that the predicted UDMH abundances are compared against.","marker":"(P. Mr´oz & et al 2024)"}],"fun_headline_variants":["PBH Poisson noise seeds ultradense dark matter halos","Poisson noise from PBHs creates ultradense dark matter halos","PBH shot noise shifts ultradense halo masses upward","Low suppression index favors multi-component dark matter halos","Stellar-mass PBHs amplify small-scale dark matter structure"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the Poisson-noise contribution from PBHs is accurately described by an exponential cutoff with the chosen location, width, and steepness ($k'\\simeq10^2\\,\\mathrm{Mpc}^{-1}$, $\\sigma\\simeq4.5\\times10^2\\,\\mathrm{Mpc}^{-1}$, free index $n$); if the true small-scale damping or PBH clustering differs, the predicted UDMH abundances and the favored dark matter composition would change.","fun_headline_variants_meta":{"raw":{"variants":["PBH Poisson noise seeds ultradense dark matter halos","Poisson noise from PBHs creates ultradense dark matter halos","PBH shot noise shifts ultradense halo masses upward","Low suppression index favors multi-component dark matter halos","Stellar-mass PBHs amplify small-scale dark matter structure"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000673,"raw_usage":{"total_tokens":3088,"prompt_tokens":996,"completion_tokens":2092,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":612,"completion_tokens_details":{"reasoning_tokens":2008}},"tokens_in":612,"tokens_out":2092,"duration_ms":13875,"temperature":1.0,"reasoning_tokens":2008,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T17:16:15.462805+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A measurement of the small-scale matter power spectrum at $k\\sim10^2$–$10^4\\,\\mathrm{Mpc}^{-1}$ (for example from Lyman-$\\alpha$ forest or 21-cm observations) that shows no isocurvature bump of the predicted height, or a high-resolution N-body simulation of stellar-mass PBH shot noise whose UDMH mass function does not shift toward higher masses with increasing $M_{\\mathrm{PBH}}$, would settle against the central claim.","supporting_citations":[{"cited_title":"2017, JCAP, 2017, 017, doi: 10.1088/1475-7516/2017/08/017","cited_arxiv_id":null,"evidence_quote":"Provides the isocurvature amplitude scaling and the exponential-cutoff form adopted for $P_{\\mathrm{iso}}(k)$."},{"cited_title":"S., & Silk, J","cited_arxiv_id":null,"evidence_quote":"Establishes the UDMH formation picture and the first-crossing multiplicity function approximation used to compute halo abundances."}],"review_version":1}