{"id":"7b247b2e-6187-4366-a7b3-3c84db478205","arxiv_id":"2412.13254","paper_version":2,"verdict":"REJECT","confidence":"HIGH","novelty_score":3.0,"correctness_risk":"high","formal_verification":"none","parameter_count":8,"one_line_summary":"A parameter-space scan shows that very large curvature-current couplings can fit the baryon asymmetry and dark matter abundance, while the electric charge of the black hole plays no role in the mechanism.","lead":"This paper studies the formation of baryon asymmetry and dark matter from the asymmetric Hawking evaporation of light, possibly charged, primordial black holes. It maps the coupling values that would simultaneously match the observed baryon abundance and dark matter relic density, with and without memory-burden effects.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper's central expansion fails internally: for its benchmark couplings, the curvature-induced chemical potential exceeds the Hawking temperature by 5–12 orders of magnitude at the epoch that dominates the asymmetry integral, invalidating Eqs. (2.15)–(2.33).","rationale":"The paper is a competent extension of known gravitational baryogenesis mechanisms to a dark sector and to memory-burdened PBHs, and it uses the public FRISBHEE code for the numerical scans. The analytical structure largely follows the prior literature, and the semiclassical formulas are cleanly derived under their stated assumptions. The Reader's weakest-assumption identification, the EFT cutoff falling below the Hawking temperature, is a serious and valid objection to the presented parameter space. However, an even more immediate problem is internal: the linear chemical-potential expansion used to derive the asymmetry rate is invalid for every benchmark coupling that fits the observed abundances. The condition µ_i ≪ T_BH, stated in Eqs. (2.12) and (2.15), fails by five to twelve orders of magnitude at the time that dominates the Q_B-L integral, because Q_B-L ∝ ∫ dt H^3 is dominated by t_in. This is not a matter of interpreting the coupling; it is a self-contained inconsistency in the semiclassical calculation itself. A corrected treatment using the exact distribution functions might open a different viable region, possibly with smaller couplings, but that region is not what the paper computes or claims. Until such a corrected calculation is provided, the central claim that Eqs. (2.33) and (3.6) define a viable parameter space is unsupported. I therefore concur with rejection, although the primary reason here is the internal breakdown of the linear approximation rather than the EFT cutoff alone.","tokens_in":20142,"tokens_out":15424,"duration_ms":145574,"concrete_test":"For benchmark BP1 (Min=1 g, β=10^-5, λ=1.5×10^18), compute t_in from Eq. (2.25), H=1/(2t_in), ˙R from Eq. (2.17), and µ_i=λ q_i^bl ˙R/(8π M_P^2). Evaluate the exact Fermi-Dirac charge-loss integral in Eq. (2.12) at t=t_in and compare it with the linearized expression in Eq. (2.15). If the ratio dQ_exact/dt / dQ_lin/dt exceeds O(1) (expected ≈10^5), then Eqs. (2.19)–(2.33) and the contours in Figs. 3–5 do not follow from the stated assumptions. Repeat the same comparison at t=10 t_in and t=100 t_in to map where the linear approximation becomes valid.","verdict_should_be":"REJECT","load_bearing_attack":"The asymmetry calculation in Eqs. (2.15)–(2.20) explicitly assumes the linear chemical-potential expansion µ_i/T_BH ≪ 1, and Eq. (2.12) states this before dropping µ_i^3 terms. For the couplings required by Eq. (2.33), this assumption is violated by orders of magnitude precisely when the asymmetry is generated. Using Eq. (2.17) and benchmark BP1 (Min=1 g, β=10^-5, λ=1.5×10^18), with t_in from Eq. (2.25), one finds H(t_in)≈2.7×10^13 GeV and T_BH≈1.1×10^13 GeV; with the SM value |1-3w|≈8×10^-3, µ_q = λ (1/3)|˙R|/(8π M_P^2) ≈ 6×10^18 GeV, so µ_q/T_BH ≈ 6×10^5. For BP2 and the memory-burdened benchmarks the ratio is even larger, reaching ~10^10 or more. Since dQ_B-L/dt ∝ H^3 ∝ t^-3, the integral for Q_B-L is dominated by the earliest times, exactly where the linear expansion fails the most. The exact Fermi-Dirac integral in Eq. (2.12) saturates when µ ≫ T, so Eqs. (2.15), (2.19), (2.33), and all derived contours and benchmark couplings are not reliable. This is an internal inconsistency, independent of any UV interpretation of Eq. (1.1). The Reader's EFT-cutoff objection is also valid, but this µ/T failure is more direct: it invalidates the semiclassical calculation even if one accepts the operator at face value.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that a curvature-coupled B-L current, Eq. (1.1), induces a chemical potential during the Hawking evaporation of light primordial black holes, thereby generating the observed baryon asymmetry, and that an analogous dark-sector coupling generates the dark matter abundance through asymmetric dark matter production. The authors derive analytic yield formulas for PBH evaporation in both radiation-dominated and PBH-dominated regimes, extend the treatment to memory-burdened PBHs, and present parameter-space scans (including FRISBHEE checks) with benchmark points in Table 1. The central claim is that a simultaneous explanation of baryon asymmetry and dark matter abundance is possible for the indicated couplings while satisfying BBN, CMB, gravitational-wave, and warm-dark-matter constraints.","tokens_in":20602,"tokens_out":23449,"duration_ms":214294,"significance":"If the calculation were correct, the paper would offer a purely gravitational cogenesis mechanism connecting the baryon asymmetry and dark matter abundance, with the memory-burden extension being a timely addition to the PBH evaporation literature. The analytic results are organized in a clear way, and the use of FRISBHEE to track PBH energy densities is a positive feature. However, the quantitative conclusions are not reliable: the central derivation rests on a linear chemical-potential expansion that is violated by many orders of magnitude at the epochs that dominate the asymmetry integrals, and the effective operator is used at energies far above its cutoff for the benchmark parameters. The claimed viable parameter space and benchmark couplings therefore do not follow from the paper's own equations.","major_comments":[{"comment":"The derivation assumes µ_i/T_BH ≪ 1 (stated before Eq. (2.12) and before Eq. (2.15)), but the benchmark couplings violate this condition by orders of magnitude at the epoch that dominates the asymmetry integral. For BP1 (M_in = 1 g, β = 10^-5, λ = 1.5×10^18), Eq. (2.25) gives t_in = M_in/(8πγ M_P^2) and thus H(t_in) = 4πγ M_P^2/M_in ≈ 2.7×10^13 GeV, while T_BH ≈ 1.1×10^13 GeV. Using Eq. (2.17) with the SM value |1-3w| ≈ 8×10^-3, the quark chemical potential is µ_q = λ (1/3)|˙R|/(8π M_P^2) ≈ 6×10^18 GeV, so µ_q/T_BH ≈ 6×10^5. For BP2 and the memory-burdened benchmarks the ratio is even larger. Since dQ_B-L/dt ∝ H^3 ∝ t^-3, the integrated charge is dominated by t ≈ t_in, exactly where the expansion fails. For µ ≫ T the Fermi-Dirac integral in Eq. (2.12) saturates and grows as µ^3, not as µ T^2, so Eqs. (2.15), (2.19), (2.20), (2.33), the contours in Figs. 3-5, and the benchmark points in Table 1 are not reliable. This is an internal inconsistency independent of any UV interpretation of Eq. (1.1).","section":"Sec. 2.1, Eqs. (2.15)-(2.20), Eq. (2.33), Table 1"},{"comment":"The effective operator ∂_µ R j^µ/M_P^2 has a cutoff Λ_cut ≈ M_P/√λ. For the benchmark couplings λ ~ 10^18-10^29, Λ_cut ranges from about 2×10^9 GeV down to about 7×10^3 GeV, while the Hawking temperatures for M_in = 1-100 g are T_BH ≈ 10^13-10^11 GeV. For example, BP2 (λ = 1.2×10^27) gives Λ_cut ≈ 7×10^4 GeV, well below T_BH ≈ 10^11 GeV. Thus the semiclassical emission calculation invokes the effective interaction for particle momenta p ~ T_BH far above its cutoff. Even if the µ/T problem were repaired by using the full Fermi-Dirac integral, the calculation would still be outside the domain of validity of Eq. (1.1).","section":"Eq. (1.1) and Table 1"},{"comment":"For β > β_c the PBHs form in a radiation-dominated universe and only later come to dominate, but the analytic Q_B-L(τ) used in the second line of Eq. (2.32) is the 'MD' expression from Eq. (2.20), which assumes a matter-dominated universe from t_in to τ. The early RD epoch contributes to Q_B-L through the trace-anomaly term (1-3w ≈ 8×10^-3), and because dQ_B-L/dt ∝ H^3 the integral is dominated by the earliest times, so the RD contribution can dominate over the PBH-dominated contribution. The correct treatment should split the integral at the PBH-radiation equality time. As written, the β > β_c contours in Figs. 3-5 and the benchmark points BP2 and BP2MB do not follow from the stated assumptions.","section":"Sec. 2.1, Eqs. (2.19)-(2.32), β>β_c branch"}],"minor_comments":[{"comment":"The DM mass is written as '105 g' where it should be 10^5 GeV; likewise, 'mDM = {105, 1010} GeV' in the text should read 10^5 and 10^10 GeV.","section":"Sec. 5, Fig. 3 caption"},{"comment":"There are several typos: 'remians' after Eq. (2.9) should be 'remains'; 'In order to to saturate' in Sec. 3 contains a duplicated 'to'; 'we find a the free parameter' in Sec. 5 should read 'we find that the free parameter'.","section":"Sec. 2 and Sec. 3"},{"comment":"The phrase 'provide observed baryon asymmetry as well as right DM abundance' should read 'provide the observed baryon asymmetry and the correct DM abundance'.","section":"Fig. 2 caption"}],"recommendation":"reject","confidential_remarks":"The rejection is based on internal technical errors in the central derivation, not on the novelty of the proposal. The µ/T ≪ 1 assumption fails by many orders of magnitude for the fitted couplings, and the β>β_c branch ignores the radiation-dominated epoch before PBH domination; these issues invalidate the benchmarks and parameter-space plots. The EFT-cutoff problem additionally suggests that the mechanism as formulated may not have a valid regime. The authors could potentially revise the calculation using the full Fermi-Dirac integrals and a piecewise cosmological evolution, but the present manuscript does not support its quantitative claims."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper's central calculation doesn't hold up: the curvature-induced chemical potential μ from Eq. (2.16) is orders of magnitude larger than the Hawking temperature T_BH at the earliest times, which dominate the asymmetry integral. For BP1 (Min=1g, β=1e-5, λ=1.5e18), μ/T_BH ~ 10^6 at t_in, and it gets worse for the memory-burdened benchmarks. The authors explicitly assume μ/T ≪ 1 when dropping μ^3 terms in Eqs. (2.12) and (2.15), so the linear expression for dQ_{B-L}/dt and the resulting λ in Eq. (2.33) are not valid. The exact Fermi-Dirac integral saturates for μ ≫ T, so the asymmetry is no longer linear in λ. This is an internal inconsistency, independent of any UV interpretation of the operator in Eq. (1.1).\n\nWhat the paper does well: it gives a systematic treatment of PBH evaporation with a dark sector and memory burden, checks analytical estimates against the FRISBHEE code, and maps a wide parameter space subject to BBN, CMB, GW, and WDM constraints. The memory-burden formalism follows the recent literature, and the figures are clear. The authors also cite the relevant prior work, including the Hook mechanism and earlier PBH cogenesis papers.\n\nBut the soft spots are serious. First, the advertised new ingredient, the PBH electric charge, decays exponentially and is gone long before memory burden matters (Fig. 1), so the title's \"charged\" aspect plays no role in the results. The mechanism is essentially the same as earlier papers, one by two of the same authors. Second, the β>β_c branch of Eq. (2.32) uses the matter-dominated expression for the asymmetry from t_in, but at formation the universe is radiation dominated by definition; the early-time integral should use the RD coefficient, so that branch is also mistreated. Third, λ and Λ_DM are solved from the observed baryon asymmetry and DM abundance, so the benchmarks are fits, not predictions.\n\nIn sum, the paper is a competent scan of a known idea, but the load-bearing assumption fails on its own terms. Without a regime where μ/T ≪ 1 and the asymmetry is still large enough, the parameter space shown is not reliable. I would not send this to peer review; it needs a fundamental reworking before it merits referee time.","headline":"The central calculation is internally inconsistent: the curvature-induced chemical potential exceeds the Hawking temperature by orders of magnitude at the epoch that determines the asymmetry, so the benchmarks and parameter space are not reliable.","tokens_in":21153,"tokens_out":5430,"would_cite":false,"duration_ms":50863,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Asymmetric Hawking radiation from light, charged primordial black holes can simultaneously produce the observed baryon asymmetry and the dark matter abundance.","keywords":["primordial black holes","Hawking radiation","baryogenesis","asymmetric dark matter","memory burden","gravitational baryogenesis","Reissner-Nordström black holes","dark matter abundance"],"falsifier":"Run the Boltzmann evolution with an explicit momentum cutoff $M_P/\\sqrt{\\lambda}$ imposed on the emission rates in Eqs. (2.16) and (3.2); if the asymmetry yield is suppressed whenever the Hawking temperature exceeds that cutoff, the joint fit to $Y_B^0 \\simeq 8.7\\times10^{-11}$ and $\\Omega_{\\rm DM}h^2 \\simeq 0.12$ fails and the central claim collapses.","tokens_in":19907,"feed_emoji":"🕳️","tokens_out":15646,"duration_ms":128797,"temperature":0.7,"pith_summary":"This paper tries to show that a single population of light primordial black holes, evaporating before big bang nucleosynthesis, can be the common source of both the observed matter–antimatter asymmetry and the dark matter relic density. The asymmetry arises because a curvature-dependent coupling between the Ricci scalar and the baryon-minus-lepton current, together with an analogous coupling to a dark-sector current, generates a chemical potential that makes Hawking radiation preferentially emit matter over antimatter and dark particles over their antiparticles. Working with Reissner-Nordström black holes, the authors find that the two measured abundances can be reproduced with couplings fixed by their Eqs. (2.33) and (3.6), within a region that respects bounds from BBN, CMB, gravitational waves, and warm dark matter. They further show that including the memory-burden effect, which slows evaporation once a black hole has radiated half its mass, moves the viable region toward lighter black holes and larger couplings. If this is right, baryogenesis and dark matter need no new TeV-scale particle physics, only the unavoidable gravitational interaction and an early population of light black holes.","feed_headline":"Charged light black holes can create both matter and dark matter","feed_subtitle":"A curvature-induced chemical potential biases Hawking radiation toward matter and sets the dark matter abundance.","key_machinery":"The load-bearing object is the curvature–current interaction $S = \\int d^4x\\sqrt{-g}\\,\\lambda\\,(\\partial_\\mu R)/M_P^2\\, j^\\mu_{B-L}$, with an identical dark-sector operator $S_{\\rm DM} = \\int d^4x\\sqrt{-g}\\,\\lambda_{\\rm DM}\\,(\\partial_\\mu R)/M_P^2\\, j^\\mu_{\\rm DM}$. In the expanding universe the trace anomaly gives $1-3w\\neq 0$, so $\\dot R\\neq 0$ and the operator acts as a chemical potential $\\mu_i = \\lambda q_i \\dot R/(8\\pi M_P^2)$ that biases Hawking emission toward particles over antiparticles. The second element is the memory-burden effect, a conjectured quantum backreaction that slows evaporation once the black hole radiates down to $qM_{\\rm in}$; it replaces the semiclassical decay rate with $dM_{\\rm BH}/dt = -\\epsilon [S(M_{\\rm BH})]^k M_P^4/M_{\\rm BH}^2$, prolonging the PBH lifetime. Integrated with the Friedmann–Boltzmann system for $\\rho_{\\rm BH}$ and $\\rho_R$, these two inputs produce the $B-L$ and dark charge yields that the paper matches to the observed abundances.","core_discovery":"The paper's central claim is that asymmetric Hawking radiation from a charged (Reissner-Nordström) primordial black hole can solve two cosmological puzzles at once, with no new particles beyond a dark-sector species. An effective operator of the form $S_{B-L} = \\int d^4x\\, \\sqrt{-g}\\, \\lambda\\, (\\partial_\\mu R)/M_P^2\\, j^\\mu_{B-L}$, motivated by gravitational baryogenesis, produces a chemical potential $\\mu_i = \\lambda q_i \\dot R/(8\\pi M_P^2)$ at the horizon; because the trace anomaly makes $1-3w$ nonzero during radiation domination, a net $B-L$ charge is radiated. The same mechanism in the dark sector, with coupling $\\lambda_{\\rm DM}$, generates an asymmetry in a dark matter species. Requiring $Y_B^0 \\simeq 8.7\\times10^{-11}$ and $\\Omega_{\\rm DM}h^2 \\simeq 0.12$ fixes $\\lambda$ and $\\Lambda_{\\rm DM}$ in ranges around $10^{17}$–$10^{27}$ and $10^8$–$10^{24}$ respectively, for PBH masses between roughly $0.1$ g and $10^5$ g. Including the memory-burden effect — with the semiclassical description breaking down at $q=0.5$ of the initial mass and a suppression index $k$ — preserves the same contours but shifts them toward lighter PBHs and larger couplings, while BBN, CMB, gravitational-wave, and warm-dark-matter bounds narrow the allowed region.","pith_inferences":["The paper leaves implicit that the fitted couplings ($\\lambda\\sim 10^{17}$–$10^{27}$) place the curvature-current operator's cutoff $M_P/\\sqrt{\\lambda}$ at or below $10^5$ GeV while the Hawking temperatures are $10^8$–$10^{13}$ GeV; a UV completion that preserves the operator's form at those temperatures is therefore needed for the mechanism to be self-consistent.","A natural extension is to compute the trace-anomaly coefficient $1-3w$ at two loops; because the yields are linear in this small number, such a correction would shift the required couplings and could open or close parts of the allowed region.","If the visible and dark couplings were related by a symmetry, the dark matter mass would become a prediction of the baryon-to-dark-matter ratio rather than an input; the paper scans $m_{\\rm DM}$ as a free parameter and could be adapted to test such a relation."],"forward_implications":["If the mechanism is correct, the baryon asymmetry and the dark matter abundance are tied to the same parameter set $\\{M_{\\rm in},\\beta,\\lambda\\}$ (plus $\\Lambda_{\\rm DM}$), so independent measurements of the PBH mass function and initial abundance would fix the curvature-current couplings.","Because evaporation must finish above the electroweak scale for sphaleron processes to act, the initial PBH mass is bounded by $M_{\\rm in}\\lesssim 10^5$ g, placing the scenario in the light-PBH window that evades stellar-mass PBH searches.","Memory burden prolongs PBH lifetimes, so the same final abundances can be produced with lighter PBHs; for fixed $\\lambda$, a larger memory-burden index $k$ shifts the viable region to smaller $M_{\\rm in}$.","The model produces an induced gravitational-wave background from PBH density fluctuations whose amplitude bounds $\\beta$; future gravitational-wave searches can therefore probe the allowed parameter region directly."],"supporting_citations":[{"why":"introduces the curvature–B−L current operator that is the starting action of the paper.","marker":"[85]"},{"why":"shows that a curvature-induced chemical potential at the horizon makes Hawking radiation asymmetric, the core mechanism used here.","marker":"[97]"},{"why":"applies asymmetric Hawking radiation to PBHs as dark matter and requires very large couplings, the baseline this paper extends.","marker":"[98]"},{"why":"demonstrates gravitational baryogenesis and dark matter from light black holes with a related curvature-scalar operator.","marker":"[55]"},{"why":"propose the memory-burden effect that stabilizes the black hole after half its mass is radiated, which the paper incorporates.","marker":"[69, 70]"},{"why":"establish the memory-burden mass window and the modified evaporation law used for nonzero $k$.","marker":"[75, 77]"},{"why":"supply the numerical machinery used to evolve PBH mass, charge, and the background bath.","marker":"[38, 39]"},{"why":"gives the gravitational-wave bound on the initial PBH fraction that closes the high-$\\beta$ region.","marker":"[107]"},{"why":"supplies the observed dark matter relic density target that the model must reproduce.","marker":"[100]"},{"why":"supplies the observed baryon asymmetry target used to fix the baryon yield.","marker":"[106]"}],"fun_headline_variants":["Charged black holes seed both matter and dark matter","Memory-burdened charged PBHs explain matter and dark matter","Asymmetric Hawking radiation from charged PBHs solves two puzzles","Gravitational chemical potential from charged PBHs sets baryon and DM asymmetry","Charged PBHs: one mechanism for baryon asymmetry and dark matter"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The curvature-current interaction in Eq. (1.1) is assumed to keep generating the asymmetry at black-hole temperatures of $10^8$–$10^{13}$ GeV, even though the fitted couplings put its natural cutoff near or below $10^5$ GeV.","fun_headline_variants_meta":{"raw":{"variants":["Charged black holes seed both matter and dark matter","Memory-burdened charged PBHs explain matter and dark matter","Asymmetric Hawking radiation from charged PBHs solves two puzzles","Gravitational chemical potential from charged PBHs sets baryon and DM asymmetry","Charged PBHs: one mechanism for baryon asymmetry and dark matter"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000856,"raw_usage":{"total_tokens":3746,"prompt_tokens":998,"completion_tokens":2748,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":614,"completion_tokens_details":{"reasoning_tokens":2671}},"tokens_in":614,"tokens_out":2748,"duration_ms":16980,"temperature":1.0,"reasoning_tokens":2671,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T13:19:35.302713+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the Boltzmann evolution with an explicit momentum cutoff $M_P/\\sqrt{\\lambda}$ imposed on the emission rates in Eqs. (2.16) and (3.2); if the asymmetry yield is suppressed whenever the Hawking temperature exceeds that cutoff, the joint fit to $Y_B^0 \\simeq 8.7\\times10^{-11}$ and $\\Omega_{\\rm DM}h^2 \\simeq 0.12$ fails and the central claim collapses.","supporting_citations":[],"review_version":1}