{"id":"8292b444-0676-4007-9230-af68e0bf6b82","arxiv_id":"2608.10977","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":9,"one_line_summary":"In mixed dark matter models combining Planck-mass black hole relics, Hawking-evaporated particles, and a WIMP, freeze-in, or axion component, the total abundance's fine-tuning is dominated by the primordial power spectrum dependence of black hole formation.","lead":"This paper calculates how finely tuned mixed dark matter models with evaporating primordial black holes must be to produce the observed dark matter abundance. It finds the required tuning is dominated by the exponential sensitivity of black hole formation to primordial fluctuations, not by the particle dark matter candidate.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Axion extension of the central claim is asserted, not demonstrated: Section 4.1.1 states the PBH-power-spectrum sensitivity dominates for QCD axions but provides no numbers or figure, so the claimed interplay cannot be checked.","rationale":"I identify the unquantified axion extension as the most load-bearing concern rather than the Press-Schechter assumption, because the exponential sensitivity of PBH formation to P_ζ is a generic feature of rare-event collapse and survives non-Gaussian corrections, whereas the paper's advertised scope includes axions and yet Section 4.1.1 provides no numbers to support the claim. The reader's verdict already flags this as a red flag; my concern agrees with that red flag but differs from the reader's formal weakest_assumption. The paper's own limitations, such as the local nature of the BG measure, are acknowledged and do not undermine the qualitative conclusion. The missing axion calculation is easily supplied and would settle the issue.","tokens_in":18,"tokens_out":16269,"duration_ms":201051,"concrete_test":"Compute the BG measures for the tripartite QCD-axion model of Section 4.1.1: set θ_i = π/√3, use Eq. (4.1) for Ω_a^sc, Eq. (2.2) for f_relic with the RD branch, and the same benchmark choices as Fig. 3 (e.g., f_relic = f_axion = 0.5, M_PBH from 1 g to 10^6 g). Evaluate Δ_Pζ, Δ_ma, and Δ_θ_i; verify that Δ_Pζ ≫ max(Δ_ma, Δ_θ_i) over the full mass range. If the hierarchy fails for any allowed benchmark, the conclusion 'robust across different DM candidates' is falsified; if it holds, the claimed extension is confirmed quantitatively.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract and conclusion claim that the fine-tuning hierarchy is robust for WIMP freeze-out, freeze-in, and QCD axion misalignment (e.g., 'this conclusion remains unchanged for all DM candidates considered, including the QCD axion'). However, Section 4.1.1 only asserts 'We found very similar results' and refers to Fig. 3 without presenting any quantitative Barbieri-Giudice values for the tripartite axion model. Appendix A2 tabulates axion-only sensitivities in non-standard cosmologies (Δ_ma, Δ_Tend of order O(1)), but not the Δ_Pζ of the full tripartite model. Because the total DM in the axion case is f_relic + f_axion (with evaporated axions contributing to ΔNeff rather than DM), the fraction of the total abundance carried by relics directly controls Δ_Pζ = (Ω_relic/Ω_total) Δ_Pζ,relic. Without specifying this benchmark and showing the resulting hierarchy, the paper's key claim of particle-candidate independence is unsupported for one of the three candidates it advertises. This is a missing support, not an internal contradiction, but it is load-bearing because the central claim explicitly includes the axion.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the sensitivity of the total dark-matter abundance in a tripartite model consisting of (i) Planck-mass relics left after PBH evaporation, (ii) particles produced by Hawking evaporation, and (iii) an independently produced component, which is taken to be WIMP freeze-out, freeze-in, or QCD axion misalignment. The authors apply the Barbieri–Giudice measure to parameters such as P_ζ, M_PBH, m_χ, and ⟨σv⟩, and find that in radiation domination the dominant fine-tuning is Δ_Pζ ∼ 10–100 because the Press–Schechter relation β = Erfc(δ_c/√(2P_ζ)) makes the PBH abundance exponentially sensitive to the curvature power spectrum. They argue that an early PBH-dominated phase removes the local sensitivity to β while still requiring a tuned inflationary amplitude, and that alternative formation channels (supercooled phase transitions, domain-wall collapse) transfer the exponential sensitivity to different parameters. The paper concludes that a natural order-unity PBH relic abundance is hard to motivate.","tokens_in":22635,"tokens_out":8749,"duration_ms":83325,"significance":"The central derivation is transparent and the qualitative hierarchy is robust: the functional dependence of β on P_ζ is exponential, and the paper correctly emphasizes that local Barbieri–Giudice measures can miss the global tuning of the inflationary sector. The manuscript is also honest about the known limitations of the BG measure and about the model dependence of the FOPT and domain-wall estimates. However, the quantitative claims for the QCD axion are not backed by any numbers or figure, and several quoted values (Δ_Pζ and Δ_M_PBH) are not reproducible from the text. These issues are fixable and do not overturn the qualitative conclusion.","major_comments":[{"comment":"The abstract and the conclusion state that the fine-tuning hierarchy is robust for all three DM candidates considered, 'including the QCD axion.' Section 4.1.1, however, only asserts that 'we found very similar results' and refers to Fig. 3, which shows the thermal-WIMP case; no Barbieri–Giudice values, benchmark parameters, or a figure for the tripartite axion model are provided. Appendix A2 tabulates axion-only sensitivities in non-standard cosmologies but does not compute Δ_Pζ for the full model. Because the total DM fraction in the axion case is f_relic + f_χ + f_axion, the relative size of f_relic directly sets Δ_Pζ, and the paper does not specify this benchmark. Please add a quantitative axion benchmark (e.g., post-inflationary θ_i = π/√3, a chosen m_a, and the relic/evaporation/axion split) and show the resulting Δ_Pζ and the competing measures, either in a table or in a figure analogous to Fig. 3.","section":"4.1.1"},{"comment":"Equation (5.2) (and the simplified scaling in Eq. (5.4)) gives f_relic ∝ M_PBH^(−3/2) P_coll, so if P_coll is independent of M_PBH, the Barbieri–Giudice measure for M_PBH is 3/2. The text quotes Δ_M_PBH ≃ 0.72, which is not consistent with this scaling. Either the quoted value comes from a different parameterization (e.g., varying M_PBH while keeping the total DM fixed by some relation to the phase-transition temperature) or it is a typo; please state the exact derivation and, if the 3/2 value is the correct one, correct the number. This does not change the qualitative conclusion that the exponential β/H dependence dominates, but the quantitative comparison should be internally consistent.","section":"5.1"},{"comment":"The quantitative fine-tuning values quoted in the text—Δ_Pζ ≃ 10–100 in §3.1, Δ_Pζ ≃ 19–32 in §3.3, and Δ_β/H ≃ 32–69 in §5.1—are given as ranges with no benchmark table. The text fixes m_χ = 100 GeV and f_χ = f_relic, but it does not state the corresponding values of P_ζ, β_i, and the fractional components for the different M_PBH values, nor does it specify how the constraint Ω_total = Ω_DM is imposed. Please provide a table of benchmark points and, where possible, an estimate of the sensitivity of the quoted ranges to δ_c and to the choice of the Press–Schechter threshold; without this the numerical hierarchy cannot be reproduced or compared across the three production mechanisms.","section":"3.1"}],"minor_comments":[{"comment":"'regime tof Hawking evaporation' should be 'regime of Hawking evaporation'.","section":"2.2.2"},{"comment":"'inequivelant' should be 'inequivalent'.","section":"3.2"},{"comment":"'Appendex' should be 'Appendix'.","section":"Table of Contents"},{"comment":"The coefficient 4×10^(−26)/g in Eq. (7.5) should be reconciled with the 4×10^(−28) quoted in Eq. (1.2); the text should state explicitly that g ≃ 108 is used in the numerical coefficient.","section":"Equations (1.2) and (7.5)"},{"comment":"'Pζ dominates the dominant fine-tuning' is awkward and should be rephrased, e.g., 'Δ_Pζ is the dominant fine-tuning measure'.","section":"Figure 3 caption"},{"comment":"The remark that bubble collisions give 'a single exponential suppression … rather than the double exponential dependence' of density-perturbation collapse is misleading; β(P_ζ) is a single exponential in 1/P_ζ, and the double-exponential statement should be clarified or removed.","section":"5.1"}],"recommendation":"major_revision","confidential_remarks":"The central relic abundance formula is attributed to the authors' own Ref. [33], and the fine-tuning results inherit this formula. This is not a circularity problem because the dominant Δ_Pζ comes from the independent Press–Schechter relation, but the manuscript would be strengthened by a self-contained derivation or a clearer statement of which parts depend on Ref. [33]. The main concern for the editor is the gap between the advertised candidate-independence claim and the missing quantitative axion analysis."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper applies the Barbieri–Giudice fine-tuning measure to a tripartite dark matter model — Planck-mass relics from evaporating ultra-light PBHs, DM particles emitted by the same PBHs, and one independent particle candidate (WIMP, freeze-in, or QCD axion). The central result is that for the first two cases, the fine-tuning is dominated by the exponential sensitivity of PBH formation to the primordial curvature power spectrum P_ζ, with Δ_Pζ of order 10–100, largely independent of particle-physics parameters. That conclusion is well supported by the Press–Schechter erfc relation and the standard abundance formulas. The paper also checks alternative PBH formation channels and shows they replace the inflationary tuning with a different exponential sensitivity, and it gives an honest treatment of the ePBHd attractor, noting that Δ=0 locally does not mean the model is natural globally.\n\nThe main derivation is clean and the qualitative hierarchy is robust. The authors are transparent about the limitations of the BG measure, the local nature of the criterion, and the ongoing debate about FOPT PBH production. The use of Ref. [33] for the relic abundance is standard and does not create a circularity problem.\n\nThe soft spots are real but not fatal. The QCD axion case in Section 4.1.1 is asserted, not demonstrated: no numbers or figure are given for the tripartite axion model, yet the abstract and conclusion claim robustness \"for all DM candidates considered, including the QCD axion.\" Without specifying the benchmark fractions and showing the resulting hierarchy, the claim of particle-candidate independence is unsupported for one of the three advertised candidates. The quantitative Δ values (10–100, 19–32, 32–69) also come without error budgets or sensitivity to the chosen thresholds (δ_c, masses, cross sections, order-one coefficients). These do not overturn the qualitative insight, but they make the numerical hierarchy less crisp than the paper suggests.\n\nThe paper is for readers working on PBH dark matter, naturalness, and mixed DM scenarios. It deserves to go through peer review: the central claim is likely correct, and the missing axion analysis can be fixed by a table or figure together with a sensitivity scan. I would recommend acceptance after major revision, with the axion numbers and uncertainty estimates as mandatory additions.","headline":"Fine-tuning in mixed DM with PBH relics is dominated by P_ζ sensitivity as claimed for WIMP and freeze-in, but the axion extension is asserted without numbers; worth publishing after major revision.","tokens_in":23298,"tokens_out":4257,"would_cite":true,"duration_ms":38700,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The tuning of PBH-relic dark matter lives in the black hole formation rate, not in the particle dark matter candidate.","keywords":["primordial black holes","dark matter","fine-tuning","Hawking evaporation","Planck-mass relics","primordial curvature power spectrum","QCD axion","freeze-in"],"falsifier":"A direct calculation of the PBH formation fraction with non-Gaussian primordial perturbations, or a numerical-relativity determination of the collapse threshold as a function of peak height, would settle the quantitative claim: if $\\beta(P_\\zeta)$ is not erfc-like over the relevant amplitude range, the reported hierarchy $\\Delta_{P_\\zeta}\\sim10$–$100$ versus $O(1)$ particle parameters would not survive.","tokens_in":22097,"feed_emoji":"🕳️","tokens_out":8752,"duration_ms":85140,"temperature":0.7,"pith_summary":"The paper asks how much the early universe must be tuned for a three-component dark-matter sector—Planck-mass relics left when tiny primordial black holes evaporate, particles emitted during that evaporation, and an independently produced dark-matter component—to add up to the observed $\\Omega_{\\rm DM}$. It finds that the dominant fine-tuning always sits in the black-hole formation step: the initial PBH fraction $\\beta(P_\\zeta)=\\mathrm{Erfc}(\\delta_c/\\sqrt{2P_\\zeta})$ depends exponentially on the primordial curvature power spectrum $P_\\zeta$, so the Barbieri–Giudice sensitivity $\\Delta_{P_\\zeta}\\sim 10$–$100$ outweighs sensitivities to particle masses, cross sections, or evaporation details. This holds whether the third component is a thermal WIMP, a freeze-in particle, or a QCD axion. An early PBH-dominated era dilutes pre-existing abundances and locally erases the $\\beta_i$ dependence, but the universe still needs a large tuned perturbation amplitude to reach that era. Supercooled phase transitions and collapsing domain walls replace the inflationary exponential with their own exponentials, leading the paper to conclude that an order-unity PBH relic abundance is hard to motivate naturally.","feed_headline":"PBH relic dark matter faces exponential fine-tuning at formation","feed_subtitle":"WIMPs, freeze-in and axions all point to the same source: the primordial power spectrum, not the particle physics.","key_machinery":"The central object is the Press–Schechter formation fraction $\\beta(P_\\zeta)=\\mathrm{Erfc}(\\delta_c/\\sqrt{2P_\\zeta})$ with collapse threshold $\\delta_c\\simeq0.45$, which converts the primordial curvature power spectrum into an initial PBH abundance. Because $\\beta$ depends exponentially on $P_\\zeta$, the Barbieri–Giudice measure $\\Delta_{P_\\zeta}=|\\partial\\ln\\Omega/\\partial\\ln P_\\zeta|$ becomes the dominant sensitivity in the tripartite DM abundance, while parameters entering through power laws, such as $m_\\chi$, $\\langle\\sigma v\\rangle$, and $M_{\\rm PBH}$, give $O(1)$ or smaller measures. The early-PBH-dominated branch is governed by formulas in which final relic and evaporated yields become independent of $\\beta_i$ after entropy injection, so the local measure gives $\\Delta_{\\beta_i}=0$ there, but the same large $P_\\zeta$ amplitude is still needed to trigger the domination epoch.","core_discovery":"The paper's central claim is that in a tripartite dark-matter model with Planck-mass PBH relics, the total abundance's parameter sensitivity is dominated by the PBH formation mechanism rather than by the particle candidate or evaporation physics. Across WIMP freeze-out, freeze-in, and QCD axion misalignment, the qualitative hierarchy is unchanged: the primordial curvature power spectrum $P_\\zeta$ gives the largest Barbieri–Giudice measure because the formation fraction $\\beta(P_\\zeta)=\\mathrm{Erfc}(\\delta_c/\\sqrt{2P_\\zeta})$ is exponentially sensitive to $P_\\zeta$. Even when an early PBH-dominated era makes the final relic and evaporated abundances independent of the initial PBH fraction $\\beta_i$, the required initial perturbation amplitude remains a tuned input. Alternative formation channels such as supercooled phase transitions and domain-wall collapse do not remove the exponential; they transfer it to parameters like $\\beta/H$ or $\\alpha_{\\rm ann}$. The paper therefore concludes that a natural realization of an order-unity PBH relic abundance is difficult to motivate.","pith_inferences":["Editorial extension: the paper itself notes that the Barbieri–Giudice measure is local and not reparameterization invariant; a global, prior-based naturalness measure would likely make the early-domination branch look even more tuned than the local $\\Delta=0$ suggests, because most of the logarithmic prior on $P_\\zeta$ produces no PBH domination at all.","Editorial extension: the same machinery could be applied to any scenario where the PBH abundance is exponentially sensitive to an underlying amplitude, such as power spectra generated by spectator fields or cosmic defects, to compare their naturalness on a common footing.","Editorial extension: a concrete test of the hierarchy would be to compute $\\beta(P_\\zeta)$ including primordial non-Gaussianity; if the functional form changes significantly, the reported values of $\\Delta_{P_\\zeta}$ would need revision even if a large sensitivity remains.","Editorial extension: future gravitational-wave constraints on the small-scale curvature power spectrum, for example from pulsar timing arrays or space-based interferometers, could locate the required peak amplitude and directly test whether it falls in the fine-tuned regime this paper identifies."],"forward_implications":["The observed DM abundance in these models is most sensitive to the amplitude of primordial curvature perturbations, so a measurement or bound on $P_\\zeta$ at PBH scales directly constrains how tuned the model must be.","Choosing a different particle DM candidate—WIMP freeze-out, freeze-in, or QCD axion—does not change the qualitative hierarchy of fine-tuning parameters.","A PBH-dominated era makes the final relic abundance insensitive to the initial PBH fraction, but the large initial perturbation amplitude needed to reach that era remains a required tuned input.","If PBH formation instead proceeds through supercooled phase transitions or domain-wall collapse, the exponential sensitivity reappears in $\\beta/H$ or $\\alpha_{\\rm ann}$, so the naturalness problem is shifted rather than solved.","A natural model of PBH relics must either produce a large curvature perturbation peak without exponential parameter dependence or invoke a formation channel whose exponential is controlled by a parameter that can be order-one without tuning."],"supporting_citations":[{"why":"Establishes the exponential dependence of the PBH mass fraction on density fluctuations, central to the $P_\\zeta$ sensitivity.","marker":"[47]"},{"why":"Supplies the Press–Schechter formalism giving the erfc form of the formation fraction used throughout.","marker":"[69]"},{"why":"Supplies the collapse threshold $\\delta_c\\simeq0.45$ used in the Press–Schechter relation.","marker":"[70]"},{"why":"Previous work by the same authors whose formulas for $f_{\\rm relic}$ and $\\beta_c$ define the relic and early-domination branches.","marker":"[33]"},{"why":"Provides the final DM yield $Y_\\chi$ expressions for evaporated particles in both radiation-dominated and PBH-dominated branches.","marker":"[64]"},{"why":"Supplies the simplified freeze-in abundance expression and the critique that local Barbieri–Giudice measures are not reparameterization invariant.","marker":"[54]"},{"why":"Presents the alternative Wilson-style fine-tuning criterion against which the local BG result is contrasted in the PBH-dominated era.","marker":"[53]"},{"why":"Gives the semi-analytic collapse probability for PBH formation from supercooled first-order phase transitions, used to estimate $\\Delta_{\\beta/H}$.","marker":"[59]"},{"why":"Provides the domain-wall collapse suppression factor $F(\\alpha_{\\rm ann})$ used for the alternative PBH formation channel.","marker":"[81]"}],"fun_headline_variants":["Exponential PBH formation tuning dominates mixed DM relics","PBH relic DM: particle physics not the source of fine-tuning","Dark matter with PBH relics: power spectrum drives the tuning","Planck-mass PBH relics shift tuning to inflation, not particles","Natural PBH relic abundance hard to motivate, paper finds"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quantitative hierarchy assumes the Gaussian Press–Schechter relation $\\beta(P_\\zeta)=\\mathrm{Erfc}(\\delta_c/\\sqrt{2P_\\zeta})$ with $\\delta_c\\simeq0.45$ and that the local Barbieri–Giudice sensitivity is the right way to judge naturalness; if non-Gaussianities change the functional form, or a global prior-based measure is used, the numerical conclusions can shift.","fun_headline_variants_meta":{"raw":{"variants":["Exponential PBH formation tuning dominates mixed DM relics","PBH relic DM: particle physics not the source of fine-tuning","Dark matter with PBH relics: power spectrum drives the tuning","Planck-mass PBH relics shift tuning to inflation, not particles","Natural PBH relic abundance hard to motivate, paper finds"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000232,"raw_usage":{"total_tokens":1513,"prompt_tokens":990,"completion_tokens":523,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":606,"completion_tokens_details":{"reasoning_tokens":438}},"tokens_in":606,"tokens_out":523,"duration_ms":5932,"temperature":1.0,"reasoning_tokens":438,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:58:50.976846+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct calculation of the PBH formation fraction with non-Gaussian primordial perturbations, or a numerical-relativity determination of the collapse threshold as a function of peak height, would settle the quantitative claim: if $\\beta(P_\\zeta)$ is not erfc-like over the relevant amplitude range, the reported hierarchy $\\Delta_{P_\\zeta}\\sim10$–$100$ versus $O(1)$ particle parameters would not survive.","supporting_citations":[{"cited_title":"Formation of Galaxies and Clusters of Galaxies by Self- Similar Gravitational Condensation","cited_arxiv_id":null,"evidence_quote":"Supplies the Press–Schechter formalism giving the erfc form of the formation fraction used throughout."},{"cited_title":"Comprehensively Constraining Ultra-Light Primordial Black Holes Through Relic Formation and Early Mergers","cited_arxiv_id":"2506.16154","evidence_quote":"Previous work by the same authors whose formulas for $f_{\\rm relic}$ and $\\beta_c$ define the relic and early-domination branches."}],"review_version":1}