{"id":"cbc38176-44d0-4eac-a10f-f96accaf723d","arxiv_id":"1908.04266","paper_version":3,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":7,"one_line_summary":"A modified potential can terminate constant-roll inflation and produce an oscillatory reheating phase, yielding a reheating temperature within a few orders of magnitude of the GUT scale.","lead":"A cosmology paper proposes a way to end a particular inflation model by changing its energy potential at late times, causing the early universe to reheat. It then uses Planck satellite data to estimate the temperature right after inflation, a quantity that affects later cosmology.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reheating-temperature formula is not derived: Eq. (64) inserts an exponent 3γ where Eq. (55) gives 1/(3γ), so Eq. (66) and Figs. 11-13 are not consequences of the model.","rationale":"I read the paper as attempting two things: (i) show that a modified constant-roll potential can end inflation and produce oscillations, and (ii) use the horizon-exit-to-present e-fold bookkeeping to convert Planck data into a reheating temperature. The first part is illustrated numerically and appears plausible as a proof of mechanism. The second part carries the title's quantitative claim, and it is where the argument breaks. I independently checked the algebra flagged by the reader. Eq. (55) gives N2 = -(1/(3γ)) ln(ρ_reh/ρ_end), so e^{N2} = (ρ_end/ρ_reh)^{1/(3γ)}. Eq. (64) instead uses (ρ_end/ρ_reh)^{3γ}. The resulting T_reh exponent in Eq. (66), 1/(12γ − 1), is not the inverse of the exponent (3γ − 4)/(3γ) that follows from the stated definitions. Figs. 11-13 therefore cannot be read as predictions of the model. The unconstrained γ reinforces the problem: the paper never fixes w_eff from interactions or data, so even a typographically corrected formula leaves T_reh a function of a free input. I agree with the reader's REJECT verdict, though I would phrase the weakest assumption as the algebraic mismatch plus the free γ rather than the free γ alone.","tokens_in":11515,"tokens_out":9124,"duration_ms":96499,"concrete_test":"Re-derive T_reh from Eqs. (34), (50), and (55) without using Eq. (64): substitute e^{N2} = (ρ_end/ρ_reh)^{1/(3γ)} into e^N = H_*/k_0, use ρ_end = 3M_P^2 H^2(φ_end) and ρ_reh = (g_reh/30)π^2 T_reh^4, and solve for T_reh. Then evaluate the ratio of this corrected T_reh to Eq. (66)'s value at γ = 2/3 with the paper's parameters (71). If the ratio differs by more than an order of magnitude, or if the sign of the exponent changes, Eq. (66) and Figs. 11-13 are invalid as printed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim—T_reh from Planck data via Eq. (66)—is unsupported because an algebraic inconsistency appears between the reheating e-fold definition and the assembled formula. Eq. (55) defines N2 = ln(a_reh/a_end) = -(1/(3γ)) ln(ρ_reh/ρ_end). Taking the exponential of both sides gives e^{N2} = (ρ_end/ρ_reh)^{1/(3γ)}. However Eq. (64) (and then Eq. (66)) uses (3M_P^2 H^2(φ_end)/((g_reh/30)π^2 T_reh^4))^{3γ}, i.e. the density ratio raised to 3γ rather than to 1/(3γ). Consequently the exponent 1/(12γ−1) in Eq. (66) is not the inverse of the T_reh exponent implied by the paper's own definitions; the correct solving exponent would be 3γ/(3γ−4), with a different sign for γ < 4/3. Figs. 11-13, including the stated 'few orders below GUT' result, therefore do not follow from the model. Separately, γ = w_eff + 1 is introduced in Eq. (56) as an average over the reheating period but is never fixed by the inflaton interactions or by data, so even after correcting the algebra the numerical temperatures remain conditional on a free parameter. The end-of-inflation position φ_end is likewise set by the ad hoc estimate λ(φ_end − φ_min)^q ≃ 1 rather than by the dynamics of the modified potential. None of these points attacks the idea that an f(φ) modification can end constant-roll inflation; they attack the paper's quantitative reheating-temperature prediction.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes a modification of the constant-roll inflaton potential, V_new(φ)=f(φ)V_old(φ) with f(φ)≈1 during inflation, so that constant-roll inflation ends and the inflaton oscillates around a minimum, reheating the Universe. The authors split the post-inflationary history into inflation, reheating, radiation, and recombination, and use the e-fold consistency condition N=ln(H_*/k_0) to solve for the reheating temperature in terms of Planck-normalized H_*, the constant-roll parameters β and M, the position φ_end, and an effective equation-of-state parameter w_eff. A numerical example with f(φ)=1−exp[−λ(φ−φ_min)^q] demonstrates exit and oscillation, and the paper plots T_reh versus γ and n_s.","tokens_in":11863,"tokens_out":9494,"duration_ms":93879,"significance":"If the quantitative derivation were correct, the paper would supply a useful complete inflation-reheating history for a constant-roll model and would show that T_reh is only weakly sensitive to n_s in this class of models. The mechanism for ending constant-roll inflation is coherent, the numerical example is explicit, and the calculation is not circular: H_* is fixed by the Planck-normalized spectrum and the constant-roll background comes from previous work. However, the central reheating-temperature formula contains an algebraic error that invalidates the reported numerical temperatures, and the result is also conditional on a free effective EoS parameter. With corrected algebra and a more transparent treatment of the EoS dependence, the paper could be publishable.","major_comments":[{"comment":"There is an inconsistency between the definition of the reheating e-fold number and the assembled formula. Eq. (55) defines N2 = −(1/(3γ)) ln(ρ_reh/ρ_end), and Eq. (59) is consistent with that, giving e^{N2} = (ρ_end/ρ_reh)^{1/(3γ)}. However, Eq. (64) contains the factor (3M_P^2 H^2(φ_end)/((g_reh/30)π^2 T_reh^4))^{3γ}, i.e. the density ratio raised to 3γ instead of 1/(3γ). Solving the correctly assembled equation for T_reh yields the exponent 3γ/(3γ−4), not 1/(12γ−1) as in Eq. (66); the sign of the correct exponent changes for γ<4/3. Since Figs. 11–13 and the conclusion that T_reh is a few orders of magnitude below the GUT scale are computed from Eq. (66), those quantitative results do not follow from the model as written.","section":"Sec. 3.4, Eqs. (55), (59), (64), (66)"},{"comment":"The position φ_end at which inflation ends is fixed by the order-of-magnitude criterion λ(φ_end−φ_min)^q ≈ 1 rather than by solving the background dynamics (2)–(4) with the modified potential (32) and (68). This matters because N1 in Eq. (37), H(φ_end) in Eq. (57), and therefore T_reh in Eq. (66) all depend on φ_end. The paper should either derive φ_end from the numerical solution or quantify how much T_reh changes under the uncertainty in this estimate.","section":"Sec. 4, Eq. (69)"},{"comment":"The effective EoS parameter w_eff (equivalently γ = w_eff + 1) is introduced as an average over the reheating era but is never determined from the inflaton's couplings or from an explicit preheating calculation. Because Eq. (66) depends on γ through a power law, all numerical temperatures reported in Figs. 11–13 are conditional on this free input. The paper acknowledges that a precise determination is not currently feasible, but the abstract and conclusion should not present T_reh as a model-independent estimate based on astrophysical data.","section":"Sec. 3.2, Eq. (56)"}],"minor_comments":[{"comment":"The sentence 'For β < 0, and for the potential (19)' appears to be a typo: Eq. (28) uses H = M coth(−βM t), which is the β<0 solution corresponding to potential (14), not potential (19).","section":"Sec. 2, after Eq. (28)"},{"comment":"The caption says 'The plot of (14) for β=0.015', but Eq. (14) is the β<0 potential; the plotted potential for β>0 is Eq. (19).","section":"Fig. 2 caption"},{"comment":"The text reads 'So NI depends on φ∗, φend, and β'; the subscript should be N1 for consistency with Eq. (34).","section":"Sec. 3.1, after Eq. (37)"},{"comment":"The statement 'we numerically solved (2) and (44)' likely refers to the Friedmann and Klein-Gordon equations (2) and (4), not Eq. (44), which is the power-spectrum amplitude.","section":"Sec. 4, final paragraph"},{"comment":"The quantities \tilde T_reh, \tilde T_CMB, and \tilde k_0 are dimensionless (in units of M_P), but this is never stated explicitly; the authors should define these reduced variables before Eq. (66).","section":"Eq. (66) and surrounding text"},{"comment":"Reference [35] contains a typo in the journal citation: 'JCAP 11 02, 021' should be 'JCAP 1102, 021'.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is not publishable in its current form because Eq. (66) is algebraically inconsistent with Eq. (55), and all quantitative reheating-temperature results depend on that formula. The inconsistency is local and fixable, so I recommend major revision rather than rejection; the authors should correct the exponent, rerun the numerical plots, and carefully state the dependence on the free EoS parameter."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the mechanism is real, the central formula is not. The authors modify the constant-roll potential (19) by multiplying with f(φ)=1−e^{−λ(φ−φ_min)^q}, so inflation exits into oscillations around a minimum. Their numerical example (β=0.015, λ=30M_P^{-2}, φ_min=6M_P) shows the exit, the oscillatory phase, and an attractor, and the e-fold count N1≈42 is reasonable. That part of the paper is coherent and worth publishing as a proof of concept.\n\nThe problem is the reheating temperature. Eq. (55) defines N2 = −(1/(3γ̄)) ln(ρ_reh/ρ_end). Exponentiating gives (ρ_end/ρ_reh)^{1/(3γ̄)}. But Eq. (64) inserts (ρ_end/ρ_reh)^{3γ̄} — the exponent is wrong by a factor of 9γ̄^2, including the sign. Eq. (66) inherits this, so the stated exponent 1/(12γ̄−1) and every number in Figs. 11–13 do not follow from the paper's own definitions. This is not a cosmetic typo: the load-bearing quantitative claim, T_reh from Planck data, is unsupported as printed.\n\nTwo softer issues. First, γ̄ is introduced in Eq. (56) as an average EoS parameter but is never fixed by the inflaton's interactions or by data; Figs. 11–13 are parametric plots, which is honest, but the 'reheating temperature a few orders below GUT' statement is a function of an unconstrained input. Second, φ_end is set by the ad hoc condition λ(φ_end−φ_min)^q≃1 (Eq. 69) rather than derived from the potential's dynamics; the numerical solution is consistent with it, but it is an estimate.\n\nThe citations are appropriate, including self-citations to the same e-fold method; there is no circularity in the argument. One small typo: in Section 2 the text says β=3 gives a flat potential, when the standard statement is β=−3 (ultra-slow-roll); the later analysis uses β>0, so it doesn't affect results.\n\nBottom line: the exit-and-oscillation mechanism is a legitimate extension of constant-roll and the numerical demonstration is convincing, but the reheating-temperature derivation has a load-bearing algebraic error. This deserves a serious referee — the flaw is fixable and the underlying idea is sound — but it should not be published until Eq. (64)/(66) are corrected and the numbers re-computed. My vote: major revision, not acceptance as-is.","headline":"The exit mechanism works, but the reheating-temperature formula has a load-bearing exponent error; don't trust the numbers, do trust the qualitative mechanism.","tokens_in":12464,"tokens_out":2841,"would_cite":false,"duration_ms":28011,"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":"The paper shows that multiplying the constant-roll potential by a function that switches on near a minimum lets constant-roll inflation end and reheat the Universe, and it derives the reheating temperature from the observed perturbation…","keywords":["constant-roll inflation","reheating temperature","inflation exit","scalar field potential","cosmological e-folds","CMB constraints","preheating equation of state"],"falsifier":"One concrete test is to compute the effective equation of state during reheating from first principles for the paper's example potential, $V_{\\rm new}(\\phi)=f(\\phi)V_{\\rm old}(\\phi)$ with $f=1-e^{-30(\\phi-6)^2}$ and $\\beta=0.015$ in Planck units, using a lattice preheating simulation; inserting the simulated $\\bar w_{\\rm eff}$ into eq. (66) gives a definite $T_{\\rm reh}$, and if that value falls outside the few-times-$10^{15}$ GeV band shown in figs. 11–13, the model's reheating prediction is wrong. A measurement of the primordial gravitational-wave background from preheating, or a cosmological bound on the reheating temperature, would provide the same test without needing the microphysics.","tokens_in":11244,"feed_emoji":"🌌","tokens_out":11100,"duration_ms":99500,"temperature":0.7,"pith_summary":"Constant-roll inflation—a scalar field moving with fixed acceleration relative to the Hubble drag, $\\ddot{\\phi}=\\beta H\\dot{\\phi}$—does not naturally end: the roll continues until a turnaround or forever. This paper proposes multiplying the constant-roll potential by a function $f(\\phi)$ that stays near $1$ during inflation and turns on near a minimum, so the field exits inflation and oscillates coherently, reheating the Universe. Counting e-folds from horizon exit to the present and matching to the observed curvature power spectrum and spectral index yields an explicit formula for the reheating temperature, eq. (66). In a worked example the reheating temperature comes out a few orders of magnitude below the grand-unification scale and is nearly insensitive to the spectral index. The result matters because it converts constant-roll inflation from an inflationary stage without an exit into a complete inflation-to-radiation history.","feed_headline":"Modified potential ends constant-roll inflation, sets reheating","feed_subtitle":"E-fold counting plus CMB data yields a GUT-scale reheating temperature in a worked example.","key_machinery":"The load-bearing object is the modified potential of eq. (32), $V_{\\rm new}(\\phi)=f(\\phi)V_{\\rm old}(\\phi)$, with $f\\simeq1$ in the inflationary epoch and $V_{\\rm new}\\simeq \\Lambda(\\phi-\\phi_{\\rm min})^q$, $q$ even, near the minimum. It carries the argument because $f$ converts the unending constant-roll trajectory into an oscillatory phase: $\\phi_{\\rm end}$ is set by the condition $\\lambda(\\phi_{\\rm end}-\\phi_{\\rm min})^q\\simeq1$, and the oscillating field behaves like matter with equation-of-state parameter $w=(q-2)/(q+2)$. Equation (66) is the working identity: it collects the e-fold budget $N=\\ln(H_*/k_0)$ from the four eras and inverts it for $T_{\\rm reh}$, with $H_*$ fixed through the power-spectrum normalization (50). The correction factor $\\xi=\\frac12(f'/f)^2$ stays small in the example, so the constant-roll predictions for $n_s$ and $r$ are preserved.","core_discovery":"The paper's central claim is that the apparent dead end of constant-roll inflation—the roll never stops, so there is no natural reheating—can be removed by extending the potential to $V_{\\rm new}(\\phi)=f(\\phi)V_{\\rm old}(\\phi)$ with $f\\simeq 1$ during inflation and $f$ dropping to zero at $\\phi_{\\rm min}$. The field then leaves inflation at a computable $\\phi_{\\rm end}$ and oscillates around $\\phi_{\\rm min}$; near the minimum the potential behaves as $V_{\\rm new}\\simeq \\Lambda(\\phi-\\phi_{\\rm min})^q$ with even $q$, and the oscillation reheats the Universe. Equating the total e-fold number from horizon exit to today with the sum of the constant-roll, reheating, radiation, and recombination contributions gives the reheating temperature (66), with $H_*$ fixed by the curvature power spectrum. In the explicit example $f(\\phi)=1-e^{-\\lambda(\\phi-\\phi_{\\rm min})^q}$ with $\\lambda=30M_P^{-2}$, $\\phi_{\\rm min}=6M_P$, $q=2$, $\\beta=0.015$, inflation ends after about 42 e-folds and the reheating temperature is within a few orders of magnitude of $10^{16}$ GeV, growing with $\\bar\\gamma=\\bar w_{\\rm eff}+1$ and changing by only about 0.1% across the observed range of $n_s$.","pith_inferences":["If the inflaton's decay products and couplings were specified, a lattice or analytic preheating computation could replace the free parameter $\\bar w_{\\rm eff}$ with a number; the temperature formula (66) would then become a sharp prediction rather than a curve family.","The same $f(\\phi)$ construction is not tied to the particular potential (19): any never-ending constant-roll or non-attractor stage could be ended by inserting a factor that vanishes at a nearby minimum before $V_{\\rm old}$ becomes negative.","A future detection of a stochastic gravitational-wave background from preheating would probe the oscillatory stage directly and could select which value of $\\bar\\gamma$ actually operated, testing the model against its own temperature curves.","The near-insensitivity of $T_{\\rm reh}$ to $n_s$ is specific to this model; applying the same e-fold counting to slow-roll models amplifies small changes in $n_s$, so the comparison offers a way to distinguish the two frameworks observationally."],"forward_implications":["The constant-roll model that fits the observed perturbation spectra can be evolved all the way past inflation, so exit and reheating need not be imposed by hand.","For the example parameters, $T_{\\rm reh}$ changes by only about 0.1% as $n_s$ runs across its observed 1-$\\sigma$ range, making the temperature a sharp prediction once the post-inflationary equation of state is fixed.","The predicted reheating temperature sits a few orders of magnitude below the grand-unification scale, in line with the usual assumption that reheating happens before or at the scale of grand unification.","The final state at $\\phi_{\\rm min}=6M_P$ is an attractor for all initial conditions with $0<\\phi(0)\\le6M_P$, so the reheating history is insensitive to the field's initial data.","Because the corrections to $n_s$ and $r$ are controlled by the small quantity $\\xi$, the observational success of the base constant-roll potential carries over to the extended model."],"supporting_citations":[{"why":"Defines the constant-roll condition and derives the $H(\\phi)$ solutions; the paper builds its inflation stage on this.","marker":"[19]"},{"why":"Provides the $\\beta>0$ potential and the slow-roll approximations for $n_s$ and $r$ used to constrain the model.","marker":"[20]"},{"why":"Supplies the observed values of the curvature power spectrum, spectral index, and tensor limit that fix $H_*$ and the allowed parameter region.","marker":"[40]"},{"why":"Establishes the coherent-oscillation reheating picture and the equation-of-state parameter $w=(q-2)/(q+2)$ near the minimum.","marker":"[8]"},{"why":"Gives the radiation energy density $\\rho_{\\rm reh}=g_{\\rm reh}\\pi^2T_{\\rm reh}^4/30$ and the thermal relations used in the e-fold counting.","marker":"[45]"},{"why":"Provides the standard decomposition of e-folds into the four eras that is assembled in eq. (34).","marker":"[44]"},{"why":"One of the template calculations for deriving the reheating temperature from the e-fold budget; the procedure here follows it.","marker":"[42]"},{"why":"Shows how the reheating temperature depends on inflation-model parameters in another setup, the method adapted for the extended potential.","marker":"[37]"},{"why":"Motivates the effective equation-of-state parameter by reviewing why perturbative reheating fails and preheating must be included.","marker":"[12]"}],"fun_headline_variants":["Potential fix ends constant-roll, sets reheating","Constant-roll inflation ends, reheats via potential tweak","Ending constant-roll inflation via potential extension","From constant-roll to GUT-scale reheating via potential change"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the whole reheating era can be represented by a single effective equation-of-state parameter $\\bar w_{\\rm eff}$ (eqs. 55–56) that is never derived from the inflaton's interactions or fixed by data; every numerical reheating temperature scales with this input.","fun_headline_variants_meta":{"raw":{"variants":["Potential fix ends constant-roll, sets reheating","Constant-roll inflation ends, reheats via potential tweak","Ending constant-roll inflation via potential extension","From constant-roll to GUT-scale reheating via potential change"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000938,"raw_usage":{"total_tokens":3972,"prompt_tokens":867,"completion_tokens":3105,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":483,"completion_tokens_details":{"reasoning_tokens":3038}},"tokens_in":483,"tokens_out":3105,"duration_ms":19954,"temperature":1.0,"reasoning_tokens":3038,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:57:21.794739+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"One concrete test is to compute the effective equation of state during reheating from first principles for the paper's example potential, $V_{\\rm new}(\\phi)=f(\\phi)V_{\\rm old}(\\phi)$ with $f=1-e^{-30(\\phi-6)^2}$ and $\\beta=0.015$ in Planck units, using a lattice preheating simulation; inserting the simulated $\\bar w_{\\rm eff}$ into eq. (66) gives a definite $T_{\\rm reh}$, and if that value falls outside the few-times-$10^{15}$ GeV band shown in figs. 11–13, the model's reheating prediction is wrong. A measurement of the primordial gravitational-wave background from preheating, or a cosmological bound on the reheating temperature, would provide the same test without needing the microphysics.","supporting_citations":[{"cited_title":"Mukhanov, Physical Foundations of Cosmology (Cambr idge Univer- sity Press 2005)","cited_arxiv_id":null,"evidence_quote":"Gives the radiation energy density $\\rho_{\\rm reh}=g_{\\rm reh}\\pi^2T_{\\rm reh}^4/30$ and the thermal relations used in the e-fold counting."},{"cited_title":"Mielczarek, Phys","cited_arxiv_id":null,"evidence_quote":"One of the template calculations for deriving the reheating temperature from the e-fold budget; the procedure here follows it."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows how the reheating temperature depends on inflation-model parameters in another setup, the method adapted for the extended potential."}],"review_version":1}