{"id":"faab272e-749a-43b8-ba5a-ecca101bd95d","arxiv_id":"2508.21011","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Evaporating primordial black holes, biased by a new gravitational interaction, can reproduce the observed baryon asymmetry once entropy dilution and chemical-potential-dependent emission are included.","lead":"This paper shows that tiny black holes evaporating in the early universe could have produced the observed excess of matter over antimatter, if a hypothetical gravitational interaction biased their radiation. The authors compute this more carefully than previous work, including how the emitted radiation dilutes the asymmetry, and find working regions for several black hole mass distributions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Pure-gradient chemical potential: Appendix A's truncated A_μ introduces a spurious electric field and invalidates the claimed μ-dependent greybody factors.","rationale":"The reader's weakest_assumption identified the imported chemical potential and the semiclassical Hawking law as the least protected premises. My concern is more specific: the greybody-factor calculation in Appendix A does not follow from the operator (1.1) because it drops the spatial component of ∂_μ K, turning a pure-gradient background into a radial electric field. This is not merely a matter of unknown UV physics; it is an internal inconsistency that affects the central numerical yields. A pure-gradient vector potential is gauge-equivalent to the free Dirac equation in a static background, so the greybody factors must be independent of the chemical potential; the large μ-dependence claimed in Fig. 2 is an artifact of the truncation. If this test confirms the greybody factors are unchanged, then the paper's 'exact greybody factors' are not exact, and the quantitative predictions (Figs. 5–8) will shift to those of the geometric-optics/low-energy approximations, which the paper shows differ by up to a factor of 15 in the rate. The central mechanism might still yield the observed asymmetry, but the paper's specific numerical claims would need revision, and the parameter space might change. Therefore a conditional acceptance with a targeted verification is appropriate, rather than a rejection, because the mechanism itself is not disproven by this flaw.","tokens_in":18055,"tokens_out":32604,"duration_ms":327465,"concrete_test":"Reproduce the greybody factors for a massless spin-1/2 particle in the Schwarzschild background with the full vector field A_μ = c M⋆^{-4} ∂_μ K (including both A_t and A_r) for a fixed M. After the exact phase redefinition ψ → exp(i c K/M⋆^4)ψ, the radial equations should be identical to the μ=0 equations; numerically verify that the transmission coefficients agree to machine precision with the standard μ=0 greybody factors (e.g., Page 1976). If they agree, recompute the asymmetry rate with μ only in the occupation numbers of Eq. (2.15) and the μ-independent greybody factors; compare the resulting n_L/s for a benchmark point (e.g., M⋆ = 10^{-3}Mpl, M_in = 10^5 Mpl, β=10^{-5}) with the paper's Fig. 5. If the yields change by more than a factor of 2, the paper's quantitative claim is not robust.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The numerical results rest on the greybody factors computed in App. A, where the interaction of Eq. (1.1) is implemented as a vector potential A_μ = −(A_0(r),0,0,0) with A_0 = −α(r_S/r)^6. But the operator (1.1) gives A_μ = M⋆^{-4} ∂_μ K, a pure gradient with F_{μν}=0. The radial component A_r = M⋆^{-4} ∂_r K is not negligible (at r_S, A_r/A_0 ∼ [12π ε(M)]^{-1}(M/Mpl)^2, large for M ≫ Mpl) and is required for the gradient property. A pure-gradient background can be removed by the phase redefinition ψ → e^{iK/M⋆^4}ψ, so in a static background the fermion scattering (and hence the greybody factors) must coincide with the μ=0 case. The truncated potential used in App. A, by contrast, has F_{0r} ≠ 0 and generates a spurious 'electric' modification of the greybody factors. The 15% (and up to factor-15) deviations in Fig. 2 are therefore not consequences of Eq. (1.1). Since the published asymmetry yields (Figs. 5–8) use these greybody factors, the central claim is not yet supported by the calculation as written.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper revisits baryogenesis from primordial black hole evaporation, assuming a dimension-eight CP-violating operator (1.1) that couples the derivative of the Kretschmann scalar to a baryon/lepton-number-violating current. The authors argue that the time dependence of the black-hole mass generates a chemical potential at the horizon, biasing Hawking radiation. They solve the coupled Friedmann/Boltzmann equations for the PBH population and radiation bath, include entropy dilution from photon emission, and compute greybody factors with a chemical-potential dependence in Appendix A. They scan monochromatic, log-normal, critical-collapse, power-law and generalized-critical-collapse mass distributions and claim that the observed baryon asymmetry, n_B/s ~ 9 x 10^-11, can be reproduced in a viable parameter region.","tokens_in":18555,"tokens_out":15553,"duration_ms":171077,"significance":"If the calculation is correct, the paper would be a useful quantitative update of PBH baryogenesis, improving on Ref. [19] by treating entropy dilution, extended mass spectra, and the full coupled evolution. The authors are appropriately conservative in stopping the evolution before the EFT breaks down and at masses near the Planck scale, and they disclose the numerical zigzag artifact in Fig. 6. However, the central new ingredient—the chemical-potential-dependent greybody factors—is undermined by the treatment of the operator as a truncated vector potential in Appendix A. Until that issue is resolved, the quantitative results, including the claimed parameter regions and the enhancement factors in Fig. 2, are not supported by the calculation as written.","major_comments":[{"comment":"The implementation of operator (1.1) as A_0 = -alpha (r_S/r)^6 is incomplete. The operator gives A_mu = M_*^{-4} partial_mu K, which has a radial component A_r = M_*^{-4} partial_r K. At the horizon, r_S A_r / A_0 ~ (3/(64 pi^2 epsilon))(M/M_pl)^3, which is large for the masses considered. Dropping A_r while keeping A_0 produces a nonzero F_{0r}, i.e., a spurious electric field. Because the full A_mu is a pure gradient, the bulk Dirac equation is equivalent to the free equation under psi -> exp(i c K/M_*^4) psi; for the static background used in a greybody calculation this phase is regular at the horizon and tends to 1 at infinity, so the greybody factors must coincide with the mu=0 case. The differences shown in Fig. 2 (a factor 1.15 and up to 15) are therefore artifacts of the truncation, and the yields in Figs. 5-8 obtained through Eq. (4.13) are not consequences of Eq. (1.1). The aut","section":"Appendix A, Eq. (A.13) and Fig. 2"},{"comment":"The central physical input—that the operator (1.1) generates the chemical potential of Eq. (1.3) at the horizon—is imported from Ref. [19] without a derivation. This is load-bearing: the asymmetry rate in Eq. (3.1), the cutoff conditions in Eq. (4.4), and all numerical results depend on it. Please provide a self-contained derivation, or at least state the thermodynamic/boundary assumptions under which a pure-gradient bulk coupling yields a non-zero chemical potential. Without this, it is difficult to distinguish the mechanism from a gauge artifact.","section":"Section 1, Eq. (1.3)"},{"comment":"The abstract and conclusion state that the observed asymmetry is reproduced for power-law mass spectra, but the results section only says 'We also observe the same qualitative features with a power-law mass distribution' without showing a figure or giving a quantitative statement. If this claim is retained, the corresponding result should be presented or quantified; otherwise the claim should be softened.","section":"Section 5 and Abstract"}],"minor_comments":[{"comment":"The quantities Delta and Sigma in Eq. (4.11) are introduced with little motivation. A short explanation of the terms in the radiation-energy equation would improve readability.","section":"Section 4.2, Eq. (4.10)"},{"comment":"The sentence 'The large suppression at large values of omega comes from the exponential function' is only true for moderate mu/T_BH. For large positive mu/T_BH the exponential can become less suppressing for some momenta, as the text later notes. Consider rephrasing.","section":"Section 2.2, Eq. (2.15)"},{"comment":"There are several typos and minor notation inconsistencies: 'F riedmann' in the Section 4 header, 't e pseudo-Riemannian' in Appendix A, and 'MP' instead of 'M_pl' in the Fig. 6 caption.","section":"Throughout"},{"comment":"The transition from Eq. (4.6) to Eq. (4.7) is terse; spelling out the definition of T0/Tev and the origin of the (1+beta T0/Tev)^{-3/4} factor would help the reader follow the analytic approximation.","section":"Section 4.1, Eq. (4.7)"}],"recommendation":"major_revision","confidential_remarks":"The decisive issue is the pure-gradient treatment in Appendix A. If the authors can show that the operator (1.1) produces a real chemical potential through a horizon boundary term and then recompute the greybody factors consistently, the paper could be publishable. As written, the main numerical results rely on a truncated vector potential whose spurious electric field is the likely source of the claimed mu-dependent greybody factors. I would ask the editor to require the authors to address this head-on before any acceptance decision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe new content here is real: full coupled Boltzmann equations for PBH evaporation with entropy dilution, extended mass distributions, and a documented numerical setup. The comparison between analytic and numerical treatments in Sec. 4 is useful, and the conservative stopping prescriptions are honest. If the physical input were right, this would be a quantitative step beyond Hamada–Iso.\n\nThe soft spot is the greybody factor computation in App. A, and it is load-bearing. The operator (1.1) is ∂_α K J^α, which in the Dirac equation becomes a vector coupling A_μ = ∂_μ K / M*^4. That is a pure gradient. For fixed black hole mass, a pure-gradient background can be removed by the phase redefinition ψ → e^{iK/M*^4}ψ, so the fermion scattering (and greybody factors) should be identical to the μ=0 case. The paper instead keeps only A_0 = −α(r_S/r)^6 and discards A_r. That truncation has F_{0r} ≠ 0 — a spurious radial electric field — which is precisely what generates the claimed μ-dependent greybody factors. The radial term is not small: at the horizon A_r/A_0 ∼ [12πε(M)]^{-1}(M/M_pl)^2, i.e., large for the masses of interest. So the 15% (or factor-15) deviations in Fig. 2 are not consequences of Eq. (1.1).\n\nThis means the central result — the reproduced asymmetry in Figs. 5–8 — is not supported by the calculation as written. The chemical potential of Eq. (1.3) is imported from Ref. [19] without an independent derivation, so the flaw is not caught there. The minor issues (the disclosed zigzag, the choice of final mass) are just that — minor. The absence of code is a hindrance, but not the main problem.\n\nIf a referee can fix the operator treatment — e.g., by showing how a time-dependent gradient produces a genuine chemical potential in the emission spectrum, or by moving the asymmetry into a separate baryon-number-violating interaction — the numerical framework is reusable. As it stands, the paper's headline claim is not supported. That said, this deserves a serious referee: the physics is subtle, the objection needs airing, and the cosmology machinery is worth preserving.\n\nMy take: send to peer review, but the authors will need to rework the greybody factors or the mechanism. I would not cite it in its current form.\n\nBest.","headline":"The cosmological machinery is solid, but the μ-dependent greybody factors — the linchpin of the asymmetry — look like an artifact of truncating a pure-gradient vector potential.","tokens_in":18920,"tokens_out":7990,"would_cite":false,"duration_ms":89453,"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":"Evaporating primordial black holes can generate the observed matter–antimatter asymmetry.","keywords":["primordial black holes","baryogenesis","Hawking radiation","chemical potential","Kretschmann scalar","baryon asymmetry","greybody factors","entropy dilution"],"falsifier":"The cleanest check is to recompute the asymmetric emission with the chemical potential evolved self-consistently inside the greybody factors instead of treated as a constant background during a given mode; if the resulting nL/s moved below roughly 9 × 10^-11 across the whole (M⋆, M) plane, the central claim would fail. Observationally, tightening PBH-abundance bounds β so that the band where nL/s ≃ 9 × 10^-11 is excluded for every M⋆ would also settle it.","tokens_in":18031,"feed_emoji":"🕳️","tokens_out":8599,"duration_ms":77373,"temperature":0.7,"pith_summary":"This paper claims that evaporating primordial black holes can explain why the universe contains more matter than antimatter. The mechanism is a higher-dimensional, CP-violating operator that couples the changing curvature (the time derivative of the Kretschmann scalar) to a baryon- or lepton-number-violating current, producing an effective chemical potential at the black-hole horizon. That chemical potential biases Hawking radiation toward emitting more baryons (or leptons) than antibaryons (or antileptons). The authors' advance is to compute this effect with chemical-potential-dependent greybody factors, to include entropy dilution from the photons the black holes also emit, and to solve the full coupled evolution of a black-hole population and the radiation bath numerically. They find that the observed baryon-to-entropy ratio of about 9 × 10^-11 is reproduced for monochromatic, log-normal, critical-collapse, and power-law mass spectra within a viable region of parameter space.","feed_headline":"Primordial black holes can forge the universe's matter surplus","feed_subtitle":"A spacetime chemical potential biases Hawking emission, and full numerical evolution reproduces the observed baryon excess.","key_machinery":"The central object is the dimension-eight, CP-violating operator (1.1), (∂_αK)J^α/M⋆^4, where K is the Kretschmann scalar, the square of the Riemann curvature tensor, and J^α is a baryon/lepton-number-violating current. Because K depends on the black-hole mass, an evaporating hole has ∂_0K ≠ 0, which generates the effective chemical potential µ = ∂_0K/M⋆^4 at the horizon. That chemical potential biases the Hawking spectrum; the paper computes the bias with greybody factors—the factors that account for the gravitational and centrifugal barriers modifying the Hawking spectrum—that depend on µ, and feeds it into coupled Boltzmann equations for the PBH population and radiation bath.","core_discovery":"On the paper's own terms, the central discovery is that the observed baryon asymmetry can be produced without new CP-violating particle decays: gravity itself, through the time-varying Kretschmann scalar of an evaporating black hole, can supply the necessary bias. As a Schwarzschild black hole loses mass, the curvature invariant K = 3M(t)^2/(4π^2 M_pl^2 r^6) changes, so ∂_0K ≠ 0, and the dimension-eight operator (∂_αK)J^α/M⋆^4 becomes active, generating a chemical potential µ = ∂_0K/M⋆^4 at the horizon. The chemical potential grows as the black hole shrinks, making Hawking emission increasingly asymmetric. Solving the coupled Friedmann–Boltzmann system with full greybody factors, including p","pith_inferences":["The authors do not discuss what a fully time-dependent treatment of the chemical potential during emission would do; recomputing the greybody factors with A0(t) evolving inside them is the minimal next calculation that could shift the quoted yields.","Because the operator couples to a generic baryon/lepton-number-violating current, the same machinery could yield model-independent lower bounds on the cutoff M⋆ once PBH abundance constraints tighten.","The early PBH-dominated epochs shown in the paper modify the expansion history, so gravitational-wave and CMB spectral-distortion searches could indirectly probe the viable parameter region.","The analysis is restricted to non-rotating black holes; extending it to Kerr black holes would change both the Kretschmann scalar and the greybody factors and could plausibly enhance or suppress the asymmetry."],"forward_implications":["If the central claim is right, the observed baryon asymmetry can be produced with no new CP-violating particle decays; the emission bias comes from gravity through a single higher-dimension operator.","Entropy dilution from photon evaporation is a controlling effect: it shrinks the viable (M⋆, M) region compared with earlier analytic estimates, especially for heavier PBHs.","The claim holds across four mass spectra, so the conclusion does not depend on a single PBH formation mechanism.","Because sphalerons convert the lepton asymmetry to a baryon asymmetry before electroweak symmetry breaking, the mechanism is also a leptogenesis scenario; washout bounds require the asymmetry to be produced below T ~ 5 × 10^11 GeV.","The results are conservative in the paper's own reading: stopping the effective field theory at M⋆ omits asymmetry that would be generated beyond the cutoff, so a UV completion could widen the viable region."],"supporting_citations":[{"why":"Introduces the Kretschmann-scalar operator and the effective chemical potential at the horizon; the paper's central mechanism is inherited from it.","marker":"[19]"},{"why":"Supplies the coupled PBH–radiation Boltzmann treatment whose numerical solution is extended here to asymmetric emission.","marker":"[12]"},{"why":"Provides the entropy-dilution estimate from PBH evaporation that the paper incorporates and improves upon.","marker":"[13]"},{"why":"Defines the mass evaporation function ε(M) entering the mass-loss rate and the chemical-potential formula.","marker":"[25]"},{"why":"The Dirac-equation and greybody-factor method in Schwarzschild spacetime that Appendix A adapts to include the chemical potential.","marker":"[62]"},{"why":"Gives the comoving PBH density and temperature-evolution equations used for the fully coupled numerical solution.","marker":"[77]"},{"why":"Supplies the two-component-spinor Dirac equation in the Newman–Penrose formalism used for the numerical greybody factors.","marker":"[82]"},{"why":"Establishes black-hole evaporation as the emission process that the asymmetry mechanism biases.","marker":"[9]"}],"fun_headline_variants":["Black holes tilt the cosmic matter-antimatter scales","Curvature warps Hawking radiation to yield matter surplus","Primordial black holes forge the seen baryon excess","Spacetime twist biases black hole decay into matter","Gravity's grip skews black hole emission to make matter"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The calculation assumes the semiclassical Hawking emission law and the effective field theory describing the chemical potential both remain valid up to the cutoff scale and near the Planck mass; if either breaks down before the chemical potential becomes large, the reproduced asymmetry would not follow.","fun_headline_variants_meta":{"raw":{"variants":["Black holes tilt the cosmic matter-antimatter scales","Curvature warps Hawking radiation to yield matter surplus","Primordial black holes forge the seen baryon excess","Spacetime twist biases black hole decay into matter","Gravity's grip skews black hole emission to make matter"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000269,"raw_usage":{"total_tokens":1422,"prompt_tokens":675,"completion_tokens":747,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":419,"completion_tokens_details":{"reasoning_tokens":681}},"tokens_in":419,"tokens_out":747,"duration_ms":7769,"temperature":1.0,"reasoning_tokens":681,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T14:38:11.553782+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"The cleanest check is to recompute the asymmetric emission with the chemical potential evolved self-consistently inside the greybody factors instead of treated as a constant background during a given mode; if the resulting nL/s moved below roughly 9 × 10^-11 across the whole (M⋆, M) plane, the central claim would fail. Observationally, tightening PBH-abundance bounds β so that the band where nL/s ≃ 9 × 10^-11 is excluded for every M⋆ would also settle it.","supporting_citations":[{"cited_title":"Page, Particle Emission Rates from a Black Hole","cited_arxiv_id":null,"evidence_quote":"The Dirac-equation and greybody-factor method in Schwarzschild spacetime that Appendix A adapts to include the chemical potential."},{"cited_title":"Page, Dirac Equation Around a Charged, Rotating Black Hole, Phys","cited_arxiv_id":null,"evidence_quote":"Supplies the two-component-spinor Dirac equation in the Newman–Penrose formalism used for the numerical greybody factors."},{"cited_title":"Hawking, Black hole explosions, Nature248 (1974) 30","cited_arxiv_id":null,"evidence_quote":"Establishes black-hole evaporation as the emission process that the asymmetry mechanism biases."}],"review_version":1}