{"id":"589dc5b2-0812-4421-ae20-b9043b8baa33","arxiv_id":"2508.08238","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Primordial black holes heavier than about 10^11 solar masses generate isocurvature perturbations that violate Planck CMB and BAO bounds, ruling them out as a significant dark matter component.","lead":"This paper shows that stupendously large black holes with masses above about 10^11 times the Sun would leave isocurvature imprints in the cosmic microwave background that exceed current Planck and BAO limits, so they cannot be a significant part of dark matter. Generalists might read it because it closes one of the few remaining proposed dark-matter mass windows and constrains exotic explanations for JWST galaxy excesses and pulsar timing signals.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: Poisson-isocurvature CMB bounds robustly exclude SLABs as significant dark matter; residual form-factor/transfer-function issues affect only boundary details.","rationale":"The strongest claim is qualitative and does not depend on exact normalization. A quick scale estimate: at k=0.05 Mpc^-1, Δ_iso^2 ≈ 2×10^-6 f (M/10^11 M_sun), against curvature power 2.1×10^-9; even with a factor ~30 transfer-function ambiguity or a suppressed high-k window, f_PBH≳10^-3 is excluded for the SLAB range, far below 'significant DM'. The reader's conditional verdict is justified by the paper's technical sloppiness, but the identified weakest assumption (T_iso^2 in Eq. 18) is not the most dangerous point; the text resolves it. I therefore agree with the conditional outcome while partially disagreeing with the stated weakest assumption. The missing square in the window function is a more concrete internal defect, but it also does not move the qualitative conclusion.","tokens_in":9623,"tokens_out":34508,"duration_ms":381749,"concrete_test":"Recompute the pink CMB+BAO boundary in Fig. 2 using the squared window function [3j1(k/k_f)/(k/k_f)]^2 instead of the linear window in Eq. (15), and separately with Eq. (18) evaluated at a=a_f versus a_rec; if the f_PBH boundary at M=10^14, 10^17 and 10^19 M_sun shifts by more than an order of magnitude or rises above f_PBH≈10^-2, the high-mass part of the claim becomes normalization-dependent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim survives scrutiny. The exclusion rests on the shot-noise isocurvature amplitude at CMB scales, P_iso ∝ f_PBH M/ρ_DM. For M≥10^11 M_sun the computed dimensionless power at k≈0.05 Mpc^-1 exceeds the Planck/CMB+BAO bounds for f_PBH values many orders of magnitude below any plausible definition of 'significant dark matter'; the k≈0.002 Mpc^-1 bin independently excludes high masses near 10^19 M_sun. The reader's Eq. (18) normalization concern is a presentation ambiguity: the text explicitly compares '[Eq. (18)] without the transfer function' to Ref. [26], which is the correct primordial normalization. The form-factor issue is real (3j1(x)/x turns negative beyond x≈4.49), and Eq. (15) should arguably contain the squared window, but this only reshapes the high-k tail beyond k_f. Since the low-k, superhorizon-at-formation modes already drive the exclusion, neither defect can rescue a significant SLAB dark matter fraction. No load-bearing flaw in the headline claim was identified.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives the isocurvature power spectrum sourced by a population of Poisson-distributed, monochromatic, non-spinning primordial black holes (PBHs) with masses in the stupendously large black hole (SLAB) range, M >~ 1e11 solar masses. The central calculation models the PBH fluid as discrete objects with finite physical radius, obtains a shot-noise isocurvature spectrum, and compares its amplitude, after applying one of two UV cutoffs, with CMB+BAO upper bounds on the primordial isocurvature power spectrum taken from Ref. [26] and with the Planck isocurvature fractions of Eq. (30). The result is an exclusion region in the (M, f_PBH) plane that rules out SLABs as a significant dark matter component. The paper also places the new bounds in the context of existing PBH constraints and discusses implications for PTA gravitational-wave backgrounds and JWST early galaxy candidates.","tokens_in":9821,"tokens_out":10193,"duration_ms":119773,"significance":"If the result holds, this is a valuable model-independent exclusion: CMB isocurvature constraints are robust and do not rely on uncertain accretion or dynamical modeling, unlike several existing SLAB bounds. The qualitative conclusion is solid because the low-k, superhorizon-at-formation part of the Poisson spectrum is standard and the comparison with external CMB/BAO constraints is natural. The paper is also careful in presenting two UV cutoff prescriptions and in using the more conservative one as the primary bound. The main technical weakness is in the finite-size form factor, which enters the written spectrum with the wrong (unsquared) functional form; this affects the high-k tail and some boundary details but does not rescue a significant SLAB dark-matter fraction.","major_comments":[{"comment":"Equation (15) as written is not a valid power spectrum. The same-PBH contribution in Eq. (11) should involve the self-convolution of the PBH density profile, not a top-hat step, and the final Fourier-space expression should contain |W(k)|^2 = [3j1(k/kf)/(k/kf)]^2 rather than a single factor 3j1(k/kf)/(k/kf). As written, the spectrum becomes negative for k/kf above about 4.49, which is unphysical and contradicts the statement that the Bessel factor 'implements a cutoff'. The low-k limit is unaffected, so the central exclusion is robust, but the high-k tail and the quantitative boundaries in Fig. 2 for the largest masses, where k/kf at CMB scales can exceed unity, should be recomputed with the squared window.","section":"Eq. (15) and surrounding derivation"},{"comment":"The normalization epoch of the compared spectrum is ambiguous. Equation (18) defines <I_DM^2>(a) with an explicit T_iso^2(a) factor, but the text says the comparison to Ref. [26] uses '[Eq. (18)] without the transfer function'. The authors should define explicitly the primordial spectrum P_II(k) = f_PBH^2 <(delta_iso^f)^2> Theta(k-k_UV) with T_iso = 1, and state that the Aiso constraints of Ref. [26] refer to the primordial isocurvature amplitude. As it stands, a reader cannot tell whether Fig. 1 and Fig. 2 include a growth factor; if T_iso were accidentally included, the bounds would shift by a factor T_iso^2(a), which is order ten at recombination for the masses considered.","section":"Eq. (18) and 'New Isocurvature Constraints'"}],"minor_comments":[{"comment":"The effective degrees of freedom are written as g_star(k) and g_star,s(k), but the arguments are temperatures; the notation should be g_star(T_f) and g_star,s(T_f).","section":"Eqs. (17) and (29)"},{"comment":"The note added states that Ref. [61] already closed the range 1e14-1e16 solar masses with Lyman-alpha data; this should be integrated into the main text and figure discussion rather than left as a post-submission note, because it changes the novelty statement for part of the claimed range.","section":"Note added and Introduction"},{"comment":"The caption for Fig. 1 does not state whether the black curves include the transfer function T_iso(a) or represent the primordial spectrum with T_iso set to unity; this should be stated explicitly, especially given the ambiguity in Eq. (18).","section":"Figure 1 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper makes a robust qualitative point and is suitable for publication after the form-factor error is corrected and the normalization convention is stated precisely. I do not expect the corrected boundaries to overturn the central conclusion that SLABs cannot be a significant dark matter component."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis letter does one useful thing: applies the known Poisson shot-noise isocurvature calculation (Inman & Ali-Haïmoud) to the SLAB mass range and compares against the recent broad-k CMB+BAO bounds from Buckley et al. The result — SLABs cannot be a significant DM component — is solid. The calculation is standard but correctly executed, the UV cut-off choices are conservative and clearly stated, and the comparison with the external constraints is the right thing to do. The paper also connects to the PTA and JWST debates, which is fair context but not load-bearing.\n\nThe soft spots are real but minor. Eq. (15) writes the isocurvature spectrum with a single j1 factor, and that expression turns negative above k/kf ~ 4.5, which is unphysical. What they presumably want is the squared window or the appropriate small-k limit; the exclusion, however, rides on low-k superhorizon modes, so the high-k sign flip does not affect the headline. The text says that the comparison uses Eq. (18) \"without the transfer function,\" which is the correct primordial normalization; the reader's worry about an ambiguous mapping is a presentation issue, not a physics error, but the sentence could be clearer. No code or data are included, which is a minor annoyance given the simplicity of the calculation.\n\nThe citation practice is fine. They use external constraints from Ref. [26] and Planck, note the concurrent Lyman-alpha closure of part of the range, and do not lean on self-citation except where the result is theirs.\n\nBottom line: the qualitative conclusion is robust and the paper is worth a serious referee. I would accept it for review; it is a clean, useful addition to the PBH constraints literature, and the boundary-level technical issues are fixable with a clarifying note and a corrected form factor. Not groundbreaking, but not sloppy.","headline":"Robust exclusion of SLABs as dark matter via Poisson isocurvature CMB bounds; residual technical quirks affect boundary details, not the headline.","tokens_in":10411,"tokens_out":3851,"would_cite":true,"duration_ms":32241,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["98.80.-k","95.35.+d"],"model":"deepseek-v4-flash","headline":"Stupendously large black holes are excluded as a significant dark matter component because their isocurvature perturbations exceed CMB limits.","keywords":["primordial black holes","stupendously large black holes","SLAB","isocurvature perturbations","cosmic microwave background","dark matter","Poisson shot noise","baryon acoustic oscillations"],"falsifier":"Recompute the excluded region using $\\langle I_{\\mathrm{DM}}^2\\rangle(a_{\\mathrm{rec}}) = T_{\\mathrm{iso}}^2(a_{\\mathrm{rec}}) f_{\\mathrm{PBH}}^2 \\langle(\\delta_{\\mathrm{iso}}^f)^2\\rangle$ as the quantity compared with the 95% CL limits from Ref. [26]; if the curve for $M_{\\mathrm{PBH}} = 10^{11}\\,M_\\odot$ and $f_{\\mathrm{PBH}} = 10^{-2}$ falls below the limit, the central exclusion claim would be falsified at that normalization.","tokens_in":9401,"feed_emoji":"🌌","tokens_out":14677,"duration_ms":124834,"temperature":0.7,"pith_summary":"Primordial black holes with masses $M \\gtrsim 10^{11}\\,M_\\odot$ — the 'stupendously large' range — are too massive to reside within galaxies and so cannot make up all of the dark matter, yet earlier work suggested they could still be a significant fraction. This paper shows that such black holes cannot be significant: their rareness and discreteness generate an isocurvature component in the dark-matter density field whose power spectrum, scaled by $f_{\\mathrm{PBH}}^2$, exceeds the 95% CMB+BAO upper limits across $M_{\\mathrm{PBH}} \\in [10^{8}, 10^{19}]\\,M_\\odot$ unless $f_{\\mathrm{PBH}}$ is very small. The resulting exclusion curve, shown in pink in Fig. 2, closes the mass window that earlier work had left open. A sympathetic reader would take the paper's central claim to be that the CMB has already settled this part of the dark-matter question.","feed_headline":"CMB limits rule out stupendously large black holes as dark matter","feed_subtitle":"The discreteness of such black holes leaves isocurvature ripples that CMB data already exclude.","key_machinery":"The carrying object is the Poisson (shot-noise) isocurvature power spectrum of the discrete PBH population. For randomly placed black holes, the formation-epoch density contrast obeys $\\langle(\\delta_{\\mathrm{iso}}^f)^2\\rangle = \\frac{2}{3\\pi}\\left(\\frac{k}{k_{\\mathrm{PBH}}}\\right)^3 \\frac{3j_1(k/k_f)}{k/k_f}$, where $k_{\\mathrm{PBH}}$ is the inverse comoving interspacing of the black holes and the spherical Bessel factor $j_1$ imposes a cutoff at the comoving PBH radius at formation. This spectrum, multiplied by $f_{\\mathrm{PBH}}^2$ and evolved with the standard isocurvature transfer function $T_{\\mathrm{iso}}(a)$, is the quantity bounded by the CMB and BAO data.","core_discovery":"The paper's central claim is that the Poisson shot noise of a monochromatic, non-spinning PBH population with $M_{\\mathrm{PBH}} \\in [10^{8}, 10^{19}]\\,M_\\odot$ produces a dark-matter isocurvature perturbation whose primordial power spectrum, $f_{\\mathrm{PBH}}^2 \\langle(\\delta_{\\mathrm{iso}}^f)^2\\rangle$, exceeds the 95% CL CMB+BAO constraints on the isocurvature power spectrum derived in Ref. [26]. The full perturbation is $\\langle I_{\\mathrm{DM}}^2\\rangle(a) \\simeq T_{\\mathrm{iso}}^2(a) f_{\\mathrm{PBH}}^2 \\langle(\\delta_{\\mathrm{iso}}^f)^2\\rangle$, with the linear transfer function $T_{\\mathrm{iso}}(a) = (2+3y)/(2+3y_f)$, and the spectrum that is compared with the constraints is Eq. (18) evaluated without that transfer function. Under both ultraviolet cut-off prescriptions considered, the exclusion region in the $(M_{\\mathrm{PBH}}, f_{\\mathrm{PBH}})$ plane rules out SLABs as a significant dark matter component; using the three-point $\\beta_{\\mathrm{iso}}$ limits from the Planck satellite gives a slightly weaker but qualitatively identical bound.","pith_inferences":["Editorial extension: the same Poisson-shot-noise argument is not specific to SLABs; applied to other discrete dark-matter candidates, such as intermediate-mass black holes or MACHOs in unconstrained mass windows, it would yield analogous CMB isocurvature bounds that could be tested against current data.","Editorial extension: if the transfer-function normalization ambiguity goes the other way, the exact $f_{\\mathrm{PBH}}$ upper limits in Fig. 2 could move by up to an order of magnitude, but the qualitative exclusion of SLABs as a dominant component would likely survive because the excluded region spans many decades.","Editorial extension: generalizing to an extended PBH mass function would make the isocurvature spectrum scale-dependent; future CMB lensing or 21-cm observations could discriminate between monochromatic and extended cases."],"forward_implications":["SLABs in the mass range $[10^8, 10^{19}]\\,M_\\odot$ cannot make up a significant part of dark matter, closing the window that earlier SLAB studies had proposed.","The nanohertz gravitational-wave background detected by pulsar timing arrays is very unlikely to be dominated by SLAB mergers, since a population large enough to explain it would violate the CMB isocurvature bounds.","Proposals that PBH-induced density fluctuations seed the early massive galaxies seen by JWST are in tension with these limits.","Future CMB experiments, e.g., LiteBIRD and the Simons Observatory, should tighten the same bounds, pushing the allowed $f_{\\mathrm{PBH}}$ to even smaller values.","If isocurvature and adiabatic modes were fully anticorrelated, the constraint on the isocurvature fraction would strengthen to $\\mathcal{O}(10^{-3})$, tightening the abundance limits further."],"supporting_citations":[{"why":"Supplies the 95% CL CMB+BAO upper limits on the primordial isocurvature power spectrum that the PBH-induced spectrum is compared against; the central target of the exclusion.","marker":"[26]"},{"why":"Provides the linear perturbation theory for a PBH-PDM mixture, including the isocurvature transfer function and the nonlinear-scale estimates used for the ultraviolet cut-offs.","marker":"[14]"},{"why":"Gives the Planck isocurvature constraints and the three-point beta_iso upper limits used for the thinner exclusion curve and the anticorrelation discussion.","marker":"[16]"},{"why":"Defines the SLAB mass range and the earlier claim that SLABs could be a significant dark matter fraction for M in [10^14, 10^19] M_sun, the target of this work.","marker":"[12]"},{"why":"Provides the earlier large-scale-structure and seed-effect constraints and the extended-mass-function framework that this paper complements and compares against.","marker":"[25]"},{"why":"Supplies the definition of the comoving PBH interspacing that sets the scale of the Poisson isocurvature spectrum.","marker":"[18]"},{"why":"The standard solution for isocurvature mode growth that underlies the transfer function used to evolve the perturbation to later epochs.","marker":"[24]"}],"fun_headline_variants":["SLABs ruled out as dark matter by Planck","Isocurvature from SLABs clashes with Planck","SLABs' isocurvature noise breaks Planck limits","No room for SLABs as dark matter, CMB says","Stupendously large black holes fail as dark matter"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The exclusion boundaries assume the formation-epoch spectrum $f_{\\mathrm{PBH}}^2 \\langle(\\delta_{\\mathrm{iso}}^f)^2\\rangle$ is the quantity to compare with the CMB/BAO limits; if the transfer function $T_{\\mathrm{iso}}^2(a_{\\mathrm{rec}})$ should be included because the CMB constrains the later-time amplitude, the boundaries shift by up to an order of magnitude.","fun_headline_variants_meta":{"raw":{"variants":["SLABs ruled out as dark matter by Planck","Isocurvature from SLABs clashes with Planck","SLABs' isocurvature noise breaks Planck limits","No room for SLABs as dark matter, CMB says","Stupendously large black holes fail as dark matter"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000912,"raw_usage":{"total_tokens":3876,"prompt_tokens":859,"completion_tokens":3017,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":475,"completion_tokens_details":{"reasoning_tokens":2934}},"tokens_in":475,"tokens_out":3017,"duration_ms":18492,"temperature":1.0,"reasoning_tokens":2934,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T17:40:20.364677+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the excluded region using $\\langle I_{\\mathrm{DM}}^2\\rangle(a_{\\mathrm{rec}}) = T_{\\mathrm{iso}}^2(a_{\\mathrm{rec}}) f_{\\mathrm{PBH}}^2 \\langle(\\delta_{\\mathrm{iso}}^f)^2\\rangle$ as the quantity compared with the 95% CL limits from Ref. [26]; if the curve for $M_{\\mathrm{PBH}} = 10^{11}\\,M_\\odot$ and $f_{\\mathrm{PBH}} = 10^{-2}$ falls below the limit, the central exclusion claim would be falsified at that normalization.","supporting_citations":[],"review_version":1}