{"id":"36d11e0a-d564-40bf-aa01-2fac633da377","arxiv_id":"2501.01867","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"A tuned Einstein-Gauss-Bonnet inflation model can create primordial black holes from asteroid-sized to tens of solar masses and secondary gravitational waves, with abundances that change dramatically in a stiff post-inflation era.","lead":"The paper shows that adding an Einstein-Gauss-Bonnet coupling to a specific inflation model can briefly speed up the field's roll, boosting early-universe density ripples enough to form primordial black holes across a wide mass range. A generalist might care because the same mechanism also produces gravitational waves in the frequency bands of LISA, DECIGO, SKA, and pulsar timing arrays, connecting inflation to dark matter and future observatories.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"One-loop correction to the USR-enhanced curvature power spectrum (Kristiano-Yokoyama, ref. [44]) is cited but not evaluated; if it is O(1) at the peak, the Pζ ~ O(0.01) claims and all PBH and GW predictions built on it would be invalidated.","rationale":"The reader's weakest assumption identifies the same load-bearing concern: the linear, Bunch-Davies initial-condition calculation of Pζ through the USR phase is used without evaluating the one-loop correction, even though ref. [44] is cited. I agree this is the most consequential risk to the central claim. The paper otherwise follows standard numerical and analytic practice: it solves the Mukhanov-Sasaki equation, matches the pivot-scale amplitude to Planck, and uses publicly available SIGWfast code for the induced gravitational waves. Those parts are not the problem. The potential loop-level failure would invalidate the specific peak height and hence every PBH and secondary-GW prediction, while leaving the background dynamics and the existence of the USR phase intact. Because the one-loop question is not settled in the literature and the paper provides no estimate, a conditional verdict is appropriate. I therefore do not change the reader's conditional verdict; the missing one-loop check should be added before the quantitative claims are accepted.","tokens_in":29496,"tokens_out":3525,"duration_ms":39650,"concrete_test":"Evaluate the one-loop correction to the curvature power spectrum at the peak scale k_peak using the in-in formalism with the cubic action of this Einstein-Gauss-Bonnet model, using the numerically computed mode functions and background; if the ratio ΔPζ_one-loop/Pζ_tree is not small (≳ 0.1), the tree-level peak and all derived PBH and gravitational wave predictions are unreliable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central quantitative output is the curvature power spectrum peak of order O(0.01) computed from the linear Mukhanov-Sasaki equation (Eqs. 25–27) through an ultra-slow-roll phase where the slow-roll parameters violate the standard assumptions (e.g., eps2 exceeds one near the fixed point). The authors acknowledge in Section II that Eq. (12) cannot be used during USR and that a full numerical solution is required, and they do solve the linear equation. However, the same USR phase that produces the large tree-level peak can also enhance interaction terms, so the one-loop correction to Pζ may be comparable to the tree-level peak frequency by frequency. This is precisely the concern raised by Kristiano and Yokoyama (ref. [44]) for single-field USR models; the paper lists this reference but never applies its criterion. Since the PBH mass function and the induced gravitational wave spectrum are both computed from the same Pζ(k) via Eqs. (33), (39), and (48), an uncontrolled O(1) loop correction would change the abundance predictions by many orders of magnitude and could erase the claimed PBH formation windows. The concern is about the internal validity of the perturbation-theory tool, not about the existence of the USR phase itself, and it is distinct from the separate issue that several f_PBH entries in Tables III and IV exceed unity.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies primordial black hole formation and secondary gravitational waves in the Mutated Hilltop inflation model coupled to an Einstein-Gauss-Bonnet term. A tanh coupling is chosen so that the field passes near a fixed point and enters an ultra-slow-roll phase, and the authors solve the background equations and the Mukhanov-Sasaki equation numerically. They report a curvature power spectrum peak of order O(0.01) on small scales while keeping the CMB pivot amplitude at the Planck value, then use the Press-Schechter formalism generalized to an equation of state 1/3 < ω ≤ 1 to compute PBH masses and abundances for two values of the potential parameter and four benchmark sets. The same power spectra are fed into the SIGWfast code to produce induced gravitational wave spectra. The headline claims are a PBH mass range from O(10^-14) to O(10) solar masses, compatibility with LIGO-Virgo and OGLE constraints, and induced GW signals within the reach of LISA, DECIGO, SKA, and pulsar timing arrays.","tokens_in":29884,"tokens_out":5743,"duration_ms":61134,"significance":"If the numerical power spectrum survives scrutiny, the paper demonstrates a useful mechanism in EGB gravity: a small-scale enhancement of Pζ without adding ad hoc features to the potential, obtained through a coupling-induced USR phase, with a full numerical treatment of the background and perturbations. Strengths include solving the Mukhanov-Sasaki equation rather than relying on slow-roll formulae, reporting explicit benchmark parameter sets, and using the public SIGWfast code for the GW spectra. The significance is currently limited by three issues: the unphysical f_PBH > 1 entries in the non-standard epoch tables, an unevaluated one-loop correction during USR, and the fact that the advertised mass range is largely set by the input choices of φc.","major_comments":[{"comment":"The abundance tables contain f_PBH values far above unity: for α = 1, set-1, ω = 0.9 the value is 7575.07, and for α = 2, set-1, ω = 0.9 it is 3402.79. Since f_PBH is defined in Eq. (40) as the fraction of dark matter in PBHs, values greater than one are internally inconsistent. The text in §V.B says only that this \"seems unphysical,\" but the same numbers are then presented as predictions and plotted. The paper must either restrict the analysis to parameter and ω combinations with f_PBH ≤ 1, clearly label the remaining entries as excluded by self-consistency, or modify the abundance computation to account for the breakdown of the Press-Schechter/linear treatment in those regions. As written, the non-standard thermal history results are not quantitatively reliable.","section":"§V, Tables III and IV"},{"comment":"The central output is the tree-level linear spectrum through the USR phase. The paper cites Kristiano and Yokoyama (ref. [44]) but never computes or bounds the one-loop correction to Pζ. Because the same USR phase that produces the O(0.01) peak can also enhance interaction terms, an O(1) loop correction would change every PBH and induced-GW prediction built from Eqs. (33), (39), and (48). The authors should either evaluate the one-loop correction in this EGB model, derive a criterion showing it is subdominant, or explicitly state that all PBH and GW predictions are conditional on that correction being negligible. Simply citing the concern without engaging it is insufficient for the load-bearing claim.","section":"§II, Eqs. (25)-(27), and ref. [44]"},{"comment":"The advertised wide mass range from O(10^-14) M⊙ to O(10) M⊙ is a direct consequence of choosing different φc values, which set k_peak in Tables I and II, followed by the mapping k → M in Eq. (33). The predicted mass is therefore a re-parameterization of the input rather than an independent model prediction. The authors should reframe the claim as a demonstration that the model can accommodate PBH masses across this range by tuning parameters, and they should discuss whether the required φc choices are natural or otherwise constrained by the CMB-scale fit and the fixed-point condition (24).","section":"Abstract and §IV, Tables I and II"},{"comment":"The text says the scalar power spectrum is obtained by \"substituting the numerical solution of the Mukhanov-Sasaki equation in Eq. (12)\" in both subsections. Eq. (12) is the slow-roll approximation and is explicitly invalid during USR, as the authors themselves note in §II.A. If Eq. (12) was actually used, the peak amplitudes in Figs. 2 and 4 are suspect; if Eq. (27) was used, the text must say so. This clarification is essential because every PBH and GW result in the paper is derived from those spectra.","section":"§IV.A and §IV.B"}],"minor_comments":[{"comment":"The set-1 GW peak frequency is given as 10^-2 Hz in §VI but as 10^-1 Hz in the conclusion; please make the two statements consistent.","section":"§VI and §VII"},{"comment":"\"Mukhanov-Sasski\" appears in §IV; the correct spelling is \"Mukhanov-Sasaki.\"","section":"Throughout"},{"comment":"The transfer function is written with the notation \"Sin( l .√ω)\" and the definition of l appears only in the same sentence; please define l = q/k explicitly and use standard mathematical notation.","section":"Eq. (36)"},{"comment":"There is a typo \"fromO(10^-14)\" in the abstract; it should read \"from O(10^-14)\".","section":"Abstract"},{"comment":"The constraint labeled \"K\" is referred to as \"Kepler(k)\" in the body text but is not expanded in the caption; please make the caption self-contained.","section":"Fig. 9 caption"},{"comment":"The quantity M_PBH ψ(M_PBH) is used before the mass function ψ(M) is explicitly defined; please define the notation at first use near Eq. (37).","section":"§V.A"}],"recommendation":"major_revision","confidential_remarks":"The paper would benefit from releasing the numerical code or providing more detailed implementation information, since the parameter space is large and the central results are numerical. The citation of ref. [44] without engaging its criterion is a particular editorial risk, as readers in this community will expect a quantitative statement about the loop correction before accepting the O(0.01) peak at face value. The text also tends to overstate detectability when curves merely touch sensitivity boundaries; the GW claims should be phrased more cautiously."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nRead the EGB-PBH paper by Yogesh and Mohammadi. My quick take: the central numerical mechanism—an ultra-slow-roll phase generated by a tanh EGB coupling, producing a P_zeta peak at the 0.01 level—is plausible and worth checking. The paper is honest about some of its own problems. But it overreaches in the abstract and ignores a known loop-level threat, so it needs revision before its claims can be trusted.\n\nWhat's new: applying the existing EGB-USR trick (already in the cited literature) to the Mutated Hilltop potential, then running the general-omega PBH and induced-GW machinery. That's a legitimate extension, though not a new mechanism. The authors solve the Mukhanov-Sasaki equation numerically, check Planck-2018 at the pivot scale, and compute mass functions and GW curves for eight parameter sets. The GW code is public (SIGWfast). That's real, reproducible work.\n\nSoft spots, in order:\n\nFirst, the one-loop concern. The paper lists Kristiano-Yokoyama in the references but never engages with it. For a USR phase that pushes P_zeta to 0.01, the linear Bunch-Davies computation may not be the whole story. If the one-loop correction is O(1) at the peak, every PBH and GW prediction built on that peak is affected. This omission is real, even though the debate is unresolved.\n\nSecond, the f_PBH > 1 entries in Tables III and IV. The text admits these are unphysical, then presents them as results. That's the wrong frame. Those parameter points should be marked as excluded, not listed as predictions.\n\nThird, the 'wide mass range from 10^-14 to 10 solar masses' in the abstract is a parameter fit, not a prediction. The position of k_peak is set by phi_c, and the mass follows from that by construction.\n\nFourth, minor: Eq. (12) is cited where Eq. (27) is meant in the power-spectrum section.\n\nWho is this for? People working on EGB inflation and PBH phenomenology. It's a useful example but not a breakthrough. I'd send it to peer review: the numerics deserve scrutiny, and the issues are fixable with a serious referee. I would not cite it in its current form.","headline":"A credible EGB-inflation PBH mechanism that overstates its predictive power and ignores a known loop-level threat; the numerics deserve scrutiny but the paper needs revision before its claims can be accepted.","tokens_in":30434,"tokens_out":4807,"would_cite":false,"duration_ms":47117,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83F05","83C35","83D05"],"pacs":[],"model":"deepseek-v4-flash","headline":"A single Einstein-Gauss-Bonnet inflation model can form primordial black holes from $10^{-14}$ to $10$ solar masses and source gravitational waves across current and future detector bands.","keywords":["primordial black holes","Einstein-Gauss-Bonnet gravity","ultra-slow-roll inflation","scalar-induced gravitational waves","non-standard thermal history","Mutated Hilltop potential","curvature power spectrum enhancement"],"falsifier":"A decisive check is to compute the one-loop correction to the curvature power spectrum for the same parameter sets; if it suppresses the $\\mathcal{O}(0.01)$ peak by a factor of ten or more, the black-hole masses, abundances, and gravitational-wave amplitudes built on that peak would no longer follow. The paper cites this concern without evaluating it, so the calculation is the clearest way to settle whether the predicted $\\Omega_{\\rm GW}h^2\\sim10^{-8}$ signals are real. A null detection of that stochastic background by pulsar timing arrays or space-based interferometers would also put the benchmark peak in tension.","tokens_in":29264,"feed_emoji":"🕳️","tokens_out":14048,"duration_ms":115718,"temperature":0.7,"pith_summary":"The paper argues that in Einstein-Gauss-Bonnet gravity, the Mutated Hilltop inflation potential paired with a hyperbolic-tangent coupling between the scalar field and the Gauss-Bonnet term can generate an ultra-slow-roll phase without adding extra features to the potential. That phase amplifies the small-scale curvature power spectrum to about $10^{-2}$, enough to form primordial black holes with masses from roughly $10^{-14}$ to $10$ solar masses. The paper then shows how a non-standard post-inflationary epoch, with stiff equation of state $1/3<\\omega\\le1$, shifts both the masses and the abundances of those black holes, in some cases producing $f_{\\rm PBH}>1$ and thereby excluding those parameter choices. It also derives the scalar-induced gravitational-wave spectrum, with peaks near $10^{-8}$ in $\\Omega_{\\rm GW}h^2$, and shows the signals fall in the reach of pulsar timing arrays and space-based interferometers. If correct, a single modified-gravity inflation model connects CMB observations, dark-matter candidates, stellar-mass black-hole merger events, and future gravitational-wave detectors.","feed_headline":"One inflation model makes black holes from asteroid to stellar mass","feed_subtitle":"A Gauss-Bonnet coupling boosts inflationary perturbations enough to seed PBHs and the gravitational waves they betray.","key_machinery":"The load-bearing ingredient is the Gauss-Bonnet coupling function $\\xi(\\phi)=\\xi_0\\tanh[\\xi_1(\\phi-\\phi_c)]$: near $\\phi=\\phi_c$ the combination $V_{,\\phi}+4\\xi_{,\\phi}V^2/(3M_p^4)$ vanishes, creating a fixed point where $\\dot\\phi$, $\\ddot\\phi$, and $\\dot H$ nearly vanish. The field stalls there for several e-folds, $\\epsilon_1=-\\dot H/H^2$ drops to about $10^{-7}$, and the slow-roll approximation breaks down, so the power spectrum must be obtained numerically from the Mukhanov-Sasaki equation $\\nu_k''+(k^2-z_s''/z_s)\\nu_k=0$ with Bunch-Davies initial conditions; Eq. (27) turns its solution into $P_\\zeta(k)=k^3|\\zeta_k|^2/(2\\pi^2)$. This $P_\\zeta(k)$ is the single output that feeds every primordial-black-hole and gravitational-wave prediction in the paper. For the non-standard thermal history the analysis is carried by the mass-scale relation of Eq. (33), the Gaussian density variance of Eq. (36), and the equation-of-state-dependent collapse threshold of Eq. (42); for the gravitational waves it is carried by the second-order source kernel of Eq. (48).","core_discovery":"On its own terms, the central discovery is that the fixed point of the Gauss-Bonnet coupling function $\\xi(\\phi)=\\xi_0\\tanh[\\xi_1(\\phi-\\phi_c)]$ creates an ultra-slow-roll phase inside an otherwise standard Mutated Hilltop inflation, stalling the field and driving the first slow-roll parameter $\\epsilon_1$ down to about $10^{-7}$. The curvature power spectrum computed from the Mukhanov-Sasaki equation then rises to $\\mathcal{O}(0.01)$ at a scale set by $\\phi_c$, while remaining about $2.1\\times10^{-9}$ at the CMB pivot scale with $n_s=0.973$ and $r=0.004$. Feeding that spectrum into the $\\omega$-dependent Gaussian collapse estimate gives primordial black hole masses spanning about $10^{-14}$ to $10$ solar masses; one benchmark set reaches $f_{\\rm PBH}\\simeq1$ for asteroid-mass black holes, so those black holes could be all of the dark matter. The same spectrum sources second-order gravitational waves whose peak sits near $\\Omega_{\\rm GW}h^2\\sim10^{-8}$ and lands in different frequency bands depending on $\\phi_c$. Including a stiff post-inflationary epoch $1/3<\\omega\\le1$ changes the black-hole mass-abundance relation and can push $f_{\\rm PBH}$ above unity, which the paper reads as ruling out those combined parameter choices.","pith_inferences":["Because the quoted peak power spectrum is $\\mathcal{O}(0.01)$, the perturbations are not weakly coupled, so the Gaussian collapse estimate used for the abundances is only a first approximation; the same peak would generate sizeable non-Gaussianity, which can change $f_{\\rm PBH}$ by orders of magnitude and is not included in Eqs. (34)-(42).","The $f_{\\rm PBH}>1$ results in stiff epochs can be inverted: instead of discarding those epochs, one could lower the power-spectrum peak until $f_{\\rm PBH}=1$, turning the black-hole abundance into a constraint on the equation of state and the reheating temperature for each peak scale.","Because $\\phi_c$ moves the peak scale while $\\xi_0$ and $\\xi_1$ set its height and width, future null detections by gravitational-wave observatories would directly constrain these three coupling parameters, independently of black-hole abundance bounds."],"forward_implications":["The same Einstein-Gauss-Bonnet setup that matches CMB pivot-scale observations can also produce primordial black holes from about $10^{-14}$ to $10$ solar masses, so a single inflation model can cover both stellar-mass merger events and lighter dark-matter candidates.","For the benchmark sets with the earliest peaks, asteroid-mass primordial black holes have $f_{\\rm PBH}\\simeq1$, meaning they could account for all of the dark matter without adding a bump or dip to the inflation potential.","A stiff post-inflationary epoch with $1/3<\\omega\\le1$ changes both the masses and the abundances of the black holes; several combinations give $f_{\\rm PBH}>1$, which the model must count as excluded even though the same parameters are allowed in a radiation-dominated epoch.","The scalar-induced gravitational-wave spectrum peaks near $\\Omega_{\\rm GW}h^2\\sim10^{-8}$, and different parameter sets place that peak in the reach of pulsar timing arrays or of space-based interferometers.","Increasing $\\omega$ in the stiff epoch slightly enhances the gravitational-wave peak, so the non-standard thermal history is testable in the amplitude and shape of the stochastic background."],"supporting_citations":[{"why":"This supplies the general-equation-of-state primordial-black-hole formation formalism, including the mass function and variance, that Section V applies.","marker":"[75]"},{"why":"This fixes the pivot-scale power spectrum value $P_\\zeta(k_\\star)=2.098\\times10^{-9}$ and the spectral-index and tensor-to-scalar-ratio bounds the model must satisfy.","marker":"[76]"},{"why":"This establishes the earlier framework for primordial-black-hole formation in non-standard post-inflationary epochs that motivates the stiff-equation-of-state generalization.","marker":"[73]"},{"why":"This provides the analytic equation-of-state-dependent critical collapse threshold used in the Gaussian collapse abundance calculation.","marker":"[187]"},{"why":"This supplies the second-order scalar-induced gravitational-wave formalism that yields the spectral integral used for the predicted signals.","marker":"[79]"},{"why":"This provides the Einstein-Gauss-Bonnet coupling mechanism for enhancing the power spectrum that this paper adapts to the Mutated Hilltop potential.","marker":"[88]"},{"why":"This demonstrates the same family of Einstein-Gauss-Bonnet coupling functions producing primordial black holes, the direct precedent for the chosen coupling.","marker":"[91]"},{"why":"This shows another Einstein-Gauss-Bonnet inflation model with a hyperbolic-tangent coupling in which black-hole formation is achieved, supporting the choice of coupling.","marker":"[92]"},{"why":"This provides the public numerical code used to produce the induced-gravitational-wave curves in Fig. 10.","marker":"[232]"}],"fun_headline_variants":["Gauss-Bonnet inflation yields PBHs that could be all dark matter","Stiff post-inflation epoch reshapes PBH mass and abundance","Gauss-Bonnet coupling seeds PBHs and secondary GWs","Non-standard epoch alters PBH formation in Gauss-Bonnet inflation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire primordial-black-hole and gravitational-wave prediction rests on assuming that the usual linear perturbation calculation with Bunch-Davies initial conditions stays valid through the ultra-slow-roll phase, where the slow-roll approximation is known to break down.","fun_headline_variants_meta":{"raw":{"variants":["Gauss-Bonnet inflation yields PBHs that could be all dark matter","Stiff post-inflation epoch reshapes PBH mass and abundance","Gauss-Bonnet coupling seeds PBHs and secondary GWs","Non-standard epoch alters PBH formation in Gauss-Bonnet inflation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001545,"raw_usage":{"total_tokens":6304,"prompt_tokens":1193,"completion_tokens":5111,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":809,"completion_tokens_details":{"reasoning_tokens":5034}},"tokens_in":809,"tokens_out":5111,"duration_ms":34405,"temperature":1.0,"reasoning_tokens":5034,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T22:20:23.430420+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive check is to compute the one-loop correction to the curvature power spectrum for the same parameter sets; if it suppresses the $\\mathcal{O}(0.01)$ peak by a factor of ten or more, the black-hole masses, abundances, and gravitational-wave amplitudes built on that peak would no longer follow. The paper cites this concern without evaluating it, so the calculation is the clearest way to settle whether the predicted $\\Omega_{\\rm GW}h^2\\sim10^{-8}$ signals are real. A null detection of that stochastic background by pulsar timing arrays or space-based interferometers would also put the benchmark peak in tension.","supporting_citations":[],"review_version":1}