{"id":"9001a55f-9979-4fc5-9513-a33cc9217289","arxiv_id":"2507.02821","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Combining relativistic accretion with memory-burdened evaporation widens the parameter space for primordial black holes as dark matter and changes dark matter and dark radiation emission predictions.","lead":"This paper computes how two processes, feeding mass into primordial black holes and slowing down their evaporation, work together. The results shift the mass range in which black holes could survive as dark matter and alter the expected dark radiation signals.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (3) and Eq. (4) contradict each other: the printed 16π accretion coefficient is ~2500 times too large to yield the Macc used in every later result.","rationale":"Reading the paper in good faith, the intended model is a parameter map combining two established effects, and the paper is transparent about its phenomenological parameters and sudden-transition assumption. The most load-bearing threat is not the sequential approximation (which for the quoted masses is numerically well separated: tacc is many orders of magnitude smaller than tev), but an internal algebraic inconsistency between the stated accretion law and the solution used everywhere. Retaining Eq. (3) as printed makes accretion about 2.5×10³ times stronger than Eq. (4) assumes, so the saturating mass Macc, which is the pivot for every later formula, would not exist. The most plausible resolution is a typographical factor in Eq. (3): under the reduced-Planck-mass convention, the Bondi coefficient should be 1/(16π), not 16π, in which case the paper's numbers follow and the scientific content is largely preserved. Because the authors never flag this factor and the paper's central claims are numerical, the safe verdict is CONDITIONAL: the paper should be accepted only after the authors confirm or correct Eq. (3) and recompute any affected tables. This matches the reader's CONDITIONAL verdict, so no verdict change is required; the reasoning differs, however, since the decisive test is algebraic rather than about overlapping accretion and evaporation.","tokens_in":25028,"tokens_out":22047,"duration_ms":244116,"concrete_test":"Take Eq. (3) verbatim and substitute Eq. (4), using ρ∞ = 4M_Pl²/[3(1+w)²t²] and Min = 6πγ(1+w)M_Pl²tin; the identity fails by a factor (16π)² ≈ 2526. As a numerical cross-check, recompute the w = 1/3, λc = λc^R case with the coefficient 1/(16π) in Eq. (3) and verify that M(t) again saturates at Macc ≈ 4 Min and that the quoted Min^0 = 1400 g survival bound is recovered. If the printed 16π coefficient is retained, the accretion term drives a divergence before evaporation begins, contradicting Fig. 1.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Direct integration of the printed accretion law, dM/dt = λc 16π M²ρ∞/M_Pl⁴ with ρ∞ = 4M_Pl²/[3(1+w)²t²] and Min = 6πγ(1+w)M_Pl² tin, gives M(t) = Min [1 − B(1 − tin/t)]^{-1} with B = 128π² λc γ/(1+w), not the paper's B = λc γ/[2(1+w)]. The two values differ by 256π² ≈ 2.5×10³. For the quoted λc^R at w = 1/3, the correct B is ≈ 1.9×10³, so the denominator in Eq. (4) crosses zero almost immediately and no finite Macc exists; the printed Eq. (3) cannot yield the Macc ≈ 4 Min used in Eqs. (5), (9)–(12), (17), (22), (26)–(35), and (43). If instead Eq. (4) is intended as the solution, Eq. (3) must carry 1/(16π) rather than 16π, which is the coefficient that follows from Bondi accretion written with the reduced Planck mass. Either way, the central quantitative claims—survival masses 5600/3200/1400 g, the βcr values, and the DM/ΔNeff contours—inherit an unstated factor that has not been validated. This is the most load-bearing point because every new result in the paper is a monotone function of the same combination (1 − λcγ/[2(1+w)])^{-1}; an order-unity error there shifts the Min bounds by a factor of about 4, while the printed error destroys the saturation assumption altogether.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the joint evolution of primordial black holes (PBHs) under relativistic or nonrelativistic accretion and memory-burdened Hawking evaporation. It derives analytic expressions for the PBH mass history, the minimum initial mass for survival until today, the maximum initial mass for evaporation before BBN, the critical abundance for PBH domination, the relic abundance of dark matter emitted during evaporation, the stable-PBH dark-matter fraction f_PBH, and the contribution to ΔN_eff. The central quantitative illustrations take κ=2, q=1/2 and show that accretion, especially relativistic accretion, shifts the survival threshold from about 5600 g (no accretion) to 3200 g (nonrelativistic) and 1400 g (relativistic) for w=1/3, with analogous shifts in the DM and ΔN_eff parameter spaces.","tokens_in":25349,"tokens_out":25871,"duration_ms":287952,"significance":"If the quantitative results are correct, the paper extends the memory-burden PBH program by demonstrating that the same mechanism that allows sub-10^15 g PBHs to survive also gains importance from early-universe accretion, lowering the survival boundary by a factor of a few and modifying the DM and dark-radiation constraints. The analytic treatment covers two equations of state and provides falsifiable predictions, e.g., the ΔN_eff reach of CMB-S4/CMB-HD for given PBH masses. The paper is not a fit to data: β, κ, q, and w are scanned parameters, and the analysis transparently rests on prior memory-burden literature and on the Bondi-type accretion formalism. The main weaknesses are the internal inconsistency of the printed accretion equation, the unquantified timescale-separation assumption, and the breakdown of the stated κ=0 limit in several general formulas.","major_comments":[{"comment":"Equation (3) as printed is inconsistent with Eq. (4). Substituting ρ∞ = 4M_Pl^2/[3(1+w)^2t^2] and t_in = M_in/[6πγ(1+w)M_Pl^2] into Eq. (3) gives dM/dt = [128π^2 λ_c γ/(1+w)] M^2/(M_in t^2), whose solution is Eq. (4) with B = 128π^2 λ_c γ/(1+w), not the printed B = λ_c γ/[2(1+w)]. Since 256π^2 ≈ 2.5×10^3, for λ_c^R at w=1/3 the denominator in Eq. (4) crosses zero almost immediately and no finite M_acc exists, whereas Eq. (5) and all subsequent results use M_acc/M_in ≈ 4. If the intended accretion law has the Bondi coefficient 1/(16π) rather than 16π, the authors must correct Eq. (3) and state this normalization explicitly, because every quantitative result in Eqs. (10)–(14), (17), (22), (26)–(35), and (43) is a monotone function of (1 − λ_cγ/[2(1+w)])^{-1}.","section":"II.A, Eqs. (3)–(5)"},{"comment":"The sequential treatment of accretion and evaporation is asserted rather than established. Equation (9) introduces t_acc without giving an expression for it, and the text only states that the interval [t_in, t_acc] is “typically much shorter” than [t_acc, t_ev]. Since the survival thresholds in Eqs. (13)–(14) and the β_cr values are computed by replacing M_in with M_acc as the effective initial mass, the authors should derive t_acc and demonstrate t_acc ≪ t_q and t_acc ≪ t_κ^ev for the full ranges of M_in and κ shown in Figs. 1–6. If accretion and evaporation overlap significantly, the quoted 5600/3200/1400 g thresholds and the corresponding DM and ΔN_eff contours would change, even though the qualitative direction of the accretion effect is robust.","section":"II.C, Eqs. (8)–(14)"},{"comment":"The claimed limit κ=0, q=1 is not recovered by the general formulas. In Eq. (12), t_κ^ev ∝ [2κ(3+2κ)]^{-1} diverges as κ→0 instead of reducing to the Hawking time M_in^3/(3ϵM_Pl^4); Eq. (14) then gives M_in ≥ 0 at κ=0; and Eqs. (26)–(27) vanish at κ=0. Since κ=0 is shown as a reference case in Figs. 2–4 and 6, the authors should state that the memory-burdened formulas apply for κ>0 and that the κ=0 curves are computed from the separate semiclassical expressions.","section":"II.B, Eqs. (7), (12)–(14); III, Eqs. (26)–(27)"},{"comment":"The factor q in Eq. (16) accounts for the mass qM_acc at memory-burden onset, but the accretion factor M_acc/M_in appears to be missing. For λ_c≠0, the late-time PBH mass is qM_acc = qM_in(1 − λ_cγ/[2(1+w)])^{-1}, so Eq. (16) should read t_BH = t_in[qβ(1 − λ_cγ/[2(1+w)])^{-1}]^{-(1+w)/(2w)}. Consequently β_cr in Eq. (17) should contain an extra factor (1 − λ_cγ/[2(1+w)]) relative to the printed expression; for relativistic accretion at w=1/3 this changes the quoted β_cr by M_acc/M_in ≈ 4. Please justify the present form or correct it.","section":"II.C, Eqs. (16)–(17)"}],"minor_comments":[{"comment":"The exponent in γ^{2w/(1+e)} should be 2w/(1+w); the “e” appears to be a typographical error.","section":"Eq. (43)"},{"comment":"The notation “m_j < t_acc^BH” and “m_j < T_acc^BH” mixes dimensions; t_acc should be T_BH^acc, the black-hole temperature at the start of the memory-burdened phase.","section":"Sec. III, Eq. (21) and surrounding text"},{"comment":"The statement that “setting κ=0 and q=1 recovers the standard case of Hawking evaporation” should be qualified, because Eq. (12) is singular at κ=0 and the first branch of Eq. (7) is the relevant one in that limit.","section":"II.B, Eq. (12)"},{"comment":"There are several typographical errors, e.g., “organisedd” in the Introduction and “lattitude” in the description of Eq. (45); these should be corrected in a proofreading pass.","section":"I and IV"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of the journal and the central idea is interesting, but the printed accretion law in Eq. (3) appears to have a factor-16π normalization error that is not propagated into the rest of the paper. This is likely fixable, but the authors must re-derive the quoted numbers after correcting the normalization. The κ=0 limit issue and the missing M_acc/M_in factor in β_cr should also be addressed before the paper can be accepted."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the paper is a useful parameter-space study, but as written its central numbers rest on an internal contradiction. I checked Eqs. (3)-(5) by hand. Integrating dM/dt = lambda_c 16pi M^2 rho_inf / M_Pl^4 with rho_inf = 4 M_Pl^2 / [3(1+w)^2 t^2] and t_in = M_in / [6pi gamma (1+w) M_Pl^2] gives M(t) = M_in [1 - (128 pi^2 lambda_c gamma / (1+w)) (1 - t_in/t)]^{-1}, not the paper's B = lambda_c gamma / [2(1+w)]. The printed coefficient is too large by a factor of 256 pi^2. With lambda_c^R at w = 1/3, B is about 1.9e3, so the denominator crosses zero almost immediately and no finite M_acc ~ 4 M_in exists. If instead Eq. (3) should read 1/(16 pi), then every later formula--survival masses 5600/3200/1400 g, beta_cr, Delta N_eff--goes through as written. This looks like a typo, but it is load-bearing and has to be fixed before the numbers can be trusted.\n\nWhat the paper does well: it cleanly combines two established effects into one mass evolution law, Eq. (9), and reduces to prior memory-burden results when lambda_c = 0 and to standard Hawking evaporation when kappa = 0, q = 1. The parameter maps for emitted DM, f_PBH, and Delta N_eff are laid out honestly, with kappa, q, beta, and w treated as free parameters rather than fitted. The author acknowledges the sudden-onset memory-burden and sequential accretion-then-evaporation idealizations. The citation pattern is appropriate, including the author's own prior derivations; self-citation here is not a problem.\n\nSofter spots: the sequential approximation is asserted rather than quantitatively demonstrated, and the sudden memory-burden transition is a known simplification. Those are secondary compared with the coefficient issue.\n\nBottom line: this deserves a serious referee, but not acceptance in its current form. The referee should ask the author to correct Eq. (3), re-run the numerics, and check the accretion coefficient against ref. [98]. Once that is done, I would cite it as a useful map of the joint parameter space. As it stands, I would not.","headline":"Useful PBH parameter map combining relativistic accretion and memory burden, but the printed accretion coefficient is off by ~2500 and every quantitative result inherits the typo.","tokens_in":25934,"tokens_out":5082,"would_cite":false,"duration_ms":51668,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Combining relativistic accretion with memory-burdened evaporation lowers the minimum initial mass for primordial black holes surviving today to roughly $10^3$ g, and shifts the resulting dark-matter and dark-radiation constraints.","keywords":["primordial black holes","memory burden","black hole evaporation","relativistic accretion","dark matter relic abundance","dark radiation","effective number of relativistic degrees of freedom","PBH dark matter"],"falsifier":"Integrate the full system $\\dot{M}=\\lambda_c M^2 \\rho_\\infty/(16\\pi M_{\\rm Pl}^4)-\\epsilon M_{\\rm Pl}^4 S^{-\\kappa}\\Theta(t-t_q)/M^2$ numerically for representative values of $\\kappa$, $q$, $M_{\\rm in}$, and $w$ without the sequential split, and compare the resulting evaporation times and survival masses with the paper's piecewise formulas; a difference of an order of magnitude in the survival threshold would show the sequential approximation fails.","tokens_in":24755,"feed_emoji":"🕳️","tokens_out":13495,"duration_ms":141290,"temperature":0.7,"pith_summary":"This paper tries to establish that two effects usually studied separately—accretion, which grows primordial black holes, and memory-burdened evaporation, which slows their mass loss—combine to change the fate of low-mass primordial black holes. In the combined picture, a black hole first accretes until its mass reaches a saturation value $M_{\\rm acc}$, then evaporates as if it had started at that larger mass, with evaporation further delayed once quantum memory-burden backreaction sets in. The result is an explicit survival threshold: with memory-burden parameters $\\kappa=2$, $q=1/2$ and a radiation-dominated background, initial masses as low as $1400$ g (for relativistic accretion) can survive to the present, versus $5600$ g with no accretion and about $10^{15}$ g in the standard semiclassical case. A sympathetic reader should care because this reopens a low-mass window for primordial black holes as dark matter and moves the parameter space for dark matter and dark radiation emitted during evaporation.","feed_headline":"Accretion lets kilogram-scale black holes survive to today","feed_subtitle":"Adding relativistic accretion and memory burden shifts the primordial black hole survival threshold down by twelve orders of magnitude.","key_machinery":"The load-bearing object is the piecewise mass-evolution law of Eq. (9): an accretion phase in which the mass follows the relativistic spherical accretion growth toward the saturation value $M_{\\rm acc}$, followed by a semiclassical evaporation phase and then a memory-burdened phase with evaporation rate $-(\\epsilon M_{\\rm Pl}^4/M^2)S^{-\\kappa}$. The memory burden is the suppression of evaporation once the black hole's entropy $S=M^2/(2M_{\\rm Pl}^2)$ has dropped, with $q$ fixing the transition mass and $\\kappa$ the suppression power. This machinery replaces a two-process dynamics by an effective single-process evolution starting from $M_{\\rm acc}$, which is what makes the survival threshold and the subsequent abundance formulas analytically tractable.","core_discovery":"The central claim is that relativistic accretion and memory-burdened evaporation are not independent corrections but combine through a simple mass rescaling. During the early dense phase, accretion dominates and the mass approaches $M_{\\rm acc}=M_{\\rm in}[1-\\lambda_c\\gamma/(2(1+w))]^{-1}$, which is about $4M_{\\rm in}$ in a radiation-dominated universe for relativistic accretion; evaporation takes over afterward with this saturated mass as the effective initial condition, and once the black hole has radiated down to a fraction $q$ of that mass, the memory burden suppresses further emission by the entropy power $S^{-\\kappa}$. The paper derives the survival threshold and finds, for $\\kappa=2$, $q=1/2$, $w=1/3$, minimum present-day survival masses of $5600$ g (no accretion), $3200$ g (nonrelativistic accretion), and $1400$ g (relativistic accretion). It further shows that accretion lowers the maximum initial mass that evaporates before big bang nucleosynthesis, shrinks the allowed $(m_j, M_{\\rm in})$ plane for dark matter produced by evaporation, shifts the $f_{\\rm PBH}$ constraints in the new mass window, and alters the predicted $\\Delta N_{\\rm eff}$ from dark radiation.","pith_inferences":["Beyond the paper: if the sequential split between accretion and evaporation is correct, the same 'replace $M_{\\rm in}$ by $M_{\\rm acc}$' rescaling should carry over to spinning black holes, where running existing Kerr evaporation codes with the rescaled initial mass would test whether the survival-window widening survives rotation.","Beyond the paper: the paper treats the transition into the memory-burdened phase as sudden; a gradual transition would smooth the mass evolution and could shift the survival threshold by an order of magnitude at fixed $\\kappa$, a direct extension of the present calculation.","Beyond the paper: because accretion changes the black-hole lifetime, the redshift and amplitude of secondary gravitational waves from evaporation would change as well; inserting the rescaled masses into existing stochastic-gravitational-wave computations gives a testable signature of the accretion boost."],"forward_implications":["Primordial black holes with initial masses as low as roughly $10^3$ g can survive to today, so the $f_{\\rm PBH}$ dark-matter window opens at masses far below the standard $10^{15}$ g threshold.","The maximum initial mass that evaporates completely before big bang nucleosynthesis drops to a few grams under relativistic accretion, so black holes that would have evaporated before BBN instead persist past it.","The allowed $(m_j, M_{\\rm in})$ parameter space for dark matter emitted during evaporation shrinks as the memory-burden parameter $\\kappa$ grows, and shrinks further when accretion, especially relativistic accretion, is included.","The critical initial abundance $\\beta_{\\rm cr}$ for primordial black holes to dominate the Universe before evaporation decreases by orders of magnitude with accretion, so a much smaller initial PBH fraction can alter cosmology.","The predicted $\\Delta N_{\\rm eff}$ contribution from dark radiation shifts with the accretion treatment and with the memory-burden parameters, and the shifts are large enough that future CMB measurements with sensitivity around $\\Delta N_{\\rm eff}\\sim 0.06$ could distinguish the cases."],"supporting_citations":[{"why":"supplies the classical steady-state spherical accretion rate that serves as the baseline for the nonrelativistic accretion term.","marker":"[84]"},{"why":"provides the relativistic accretion framework and the coefficient used for $\\lambda_c^R$ in the accretion rate.","marker":"[97]"},{"why":"derives the mass-growth solution that the paper takes as the accretion phase of the piecewise evolution.","marker":"[98]"},{"why":"introduces the memory-burden evaporation timescale that sets the delayed evaporation phase.","marker":"[51]"},{"why":"establishes the memory-burden survival-mass estimate and observational constraints that the paper extends to include accretion.","marker":"[52]"},{"why":"shows that the total number of emitted particles is independent of the emission rate, the basis for the dark-matter abundance formulas.","marker":"[53]"},{"why":"supplies the semiclassical evaporation rate that underlies the first evaporation phase.","marker":"[25]"},{"why":"provides the gamma-ray observations used to draw the $f_{\\rm PBH}$ exclusion curves in the new low-mass window.","marker":"[110]"}],"fun_headline_variants":["Accretion plus memory burden lets kg-scale black holes survive","Relativistic accretion lowers black hole survival mass to 1.4 kg","Memory-burdened evaporation and accretion shift PBH survival window","Joint accretion and memory effects rescue small primordial black holes","Accretion + memory burden: black holes as light as 1.4 kg survive"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that a black hole finishes growing by accretion and saturates before it starts losing much mass to evaporation, so the two processes can be treated one after the other; if they overlap substantially, the numerical thresholds shift.","fun_headline_variants_meta":{"raw":{"variants":["Accretion plus memory burden lets kg-scale black holes survive","Relativistic accretion lowers black hole survival mass to 1.4 kg","Memory-burdened evaporation and accretion shift PBH survival window","Joint accretion and memory effects rescue small primordial black holes","Accretion + memory burden: black holes as light as 1.4 kg survive"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000681,"raw_usage":{"total_tokens":3131,"prompt_tokens":1022,"completion_tokens":2109,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":638,"completion_tokens_details":{"reasoning_tokens":2016}},"tokens_in":638,"tokens_out":2109,"duration_ms":18163,"temperature":1.0,"reasoning_tokens":2016,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:20:50.809572+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Integrate the full system $\\dot{M}=\\lambda_c M^2 \\rho_\\infty/(16\\pi M_{\\rm Pl}^4)-\\epsilon M_{\\rm Pl}^4 S^{-\\kappa}\\Theta(t-t_q)/M^2$ numerically for representative values of $\\kappa$, $q$, $M_{\\rm in}$, and $w$ without the sequential split, and compare the resulting evaporation times and survival masses with the paper's piecewise formulas; a difference of an order of magnitude in the survival threshold would show the sequential approximation fails.","supporting_citations":[],"review_version":1}