{"id":"ad0e2b8e-2ec0-49d1-a597-e2db128ce2a2","arxiv_id":"hep-ph/9407247","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Parametric resonance makes post-inflationary scalar production much more efficient than first-order Born estimates, by a factor of order [λ log(1/λ)]⁻¹ for a λφ⁴ inflaton.","lead":"A rapidly oscillating inflaton field after inflation drives parametric resonance, producing particles far faster than standard perturbative estimates suggest. This paper computes the enhancement for quartic and quadratic inflaton potentials and locates where the effect shuts off.","discovery_kind":"extension","skeptic_critique":{"model":"selfhosted/deepseek-v4.1-flash","headline":"The central claim is internally consistent within its stated perturbative scope. The load-bearing soft spot is the backreaction ansatz (Eqs. 67-70) that fixes the enhancement factor; the rescattering caveat the reader flags is acknowledged and plausibly subdominant for the quoted parameters.","rationale":"The paper's central claim is carefully scoped to the perturbative regime (Eq. 44, clarified by Eq. 53), and within that scope the derivation in Sections 4-5 and Appendices A-B is internally consistent: the Fourier coefficients (77)-(78) are correct, the sweep-rate algebra (69)-(72) checks out dimensionally, and the self-consistent inversion (74)-(76) reproduces the stated enhancement. I verified the arithmetic of Eq. (72) and the self-consistency of Eq. (76) against Eq. (74) and found no error. The reader's three reasons for CONDITIONAL are all legitimate, but they are not equally load-bearing. The abstract's claim to obtain both the lowest and the highest freeze-out bounds (abstract:s8) is an overclaim: the body attributes the highest bounds (136)-(137) to Ref. 13, and the reader correctly notes this. The rescattering caveat is acknowledged in the Introduction and, for the quoted parameters, a rough estimate of the χχ scattering rate against the resonance growth rate suggests it is subdominant during the buildup; the reader's phrasing that it could reduce the rates 'possibly by a large amount' is stronger than the numbers support, though the general point that it is unresolved stands. The backreaction closure (Eqs. 67-70) is the item that actually fixes the quantitative enhancement, and it is presented as an energy-conservation argument without an independent check against the coupled mode equations. That is the concern I would weight most heavily, and it is the one the reader lists first in the rationale but does not make the weakest_assumption. I therefore agree with the reader's verdict of CONDITIONAL but only partially with the identification of rescattering as the single weakest assumption; the backreaction ansatz is at least as load-bearing for the quoted numbers, and the qualitative conclusion is robust to O(1) errors in it. No change to the verdict is warranted: the paper's durable contribution, that parametric resonance dominates Born estimates in the narrow-resonance regime, survives the concerns raised, while the specific reheating temperatures carry the O(1) uncertainty the reader flags.","tokens_in":24483,"tokens_out":22087,"duration_ms":427663,"concrete_test":"Numerically integrate the coupled system for a representative point (λ=10⁻¹², h=10⁻¹³, φ₀=10⁻⁴M_P): the inflaton mode equation including the backreaction force 2hφ⟨χ²⟩, together with the χ mode equation (48) for a band of comoving k, starting from vacuum. Extract the effective Γ_χ^(res) from the late-time growth of ρ_χ and compare to Eq. (81). If the numerical rate differs from Eq. (81) by more than a factor of ~3, the quoted T_rh in Eq. (101) shifts by the corresponding factor and the O(1) uncertainty the reader flags is confirmed as material.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The enhancement factors (83)-(84) and the reheating temperature (101) are not read off the mode equation directly; they follow from the backreaction closure in Section 5. Eq. (67) posits that the inflaton loses energy only through χ production at rate Γ, Eq. (68) then forces φ̇₀/φ₀ = -(H + ¼Γ), and Eq. (69) converts this into the resonance sweep rate dΔω_k/dt = ¼Γω_res. Eq. (70) is the resulting energy budget, and Eq. (76) is its inversion. Three things are assumed rather than derived: (i) the frequency-amplitude relation ω = c√λφ₀ (Eq. 31) survives backreaction, i.e. the 2hφ⟨χ²⟩ force term in the inflaton equation is negligible; (ii) the single dominant band saturates the sum in (74); (iii) the O(1) coefficient inside the log, 32π²/c⁴λ, is correct. The paper never checks (i)-(iii) against a direct solution of the coupled inflaton-plus-χ mode equations. This matters because the log is the entire difference between the resonance rate and the Born rate: a factor-of-few error in the log argument propagates into Γ_χ^(res) and hence into T_rh. The qualitative statement (enhancement ~1/λ) is robust to such errors, but the quoted 10¹² GeV is not. The rescattering caveat the reader identifies is real but is acknowledged in the Introduction and, by a rough estimate (Γ_scatt/μ₊ ~ 10⁻¹⁴ for λ=10⁻¹², h=0.1λ), is not obviously dominant for the quoted parameters; the backreaction closure is the sharper unaddressed item.","agreement_with_reader":"partial"},"referee_report":{"model":"selfhosted/deepseek-v4.1-flash","summary":"The paper computes scalar particle production by a rapidly oscillating inflaton after chaotic inflation, using Bogoliubov transformations and paying particular attention to parametric resonance. For a quartic potential V = λφ⁴ the authors argue that, within the regime where perturbation theory in the inflaton couplings is valid (their Eq. 44), resonance dominates over first-order Born estimates by a factor of order [λ log(1/λ)]⁻¹ (Eqs. 83–84), giving a reheating temperature T_rh ~ 10⁻¹ M_P λ^{1/4} (h/λ)² log⁻¹(1/λ) (Eq. 101), about 10¹² GeV for λ ~ 10⁻¹² and h/λ ~ 10⁻¹. For a quadratic potential V = ½m²φ² the resonance is narrow and time-limited, and the paper derives lowest and highest possible freeze-out boundaries (Eqs. 127, 132, 136, 137) below which particle creation becomes inefficient and the inflaton may decouple and provide cold dark matter. Two appendices supply the Bogoliubov averaging for the first instability band and the time-dependent oscillator formalism.","tokens_in":24853,"tokens_out":7499,"duration_ms":149963,"significance":"This is a substantial and timely contribution. It is one of the papers that established parametric resonance as the dominant reheating channel and corrected the Born-approximation estimates that had been standard since Dolgov–Linde, Abbott–Farhi–Wise, and Albrecht–Steinhardt–Turner–Wilczek. The derivation is checkable: the mode equation (50), the resonance condition (54), and the growth rate (60) are given, and Appendices A and B make the Bogoliubov machinery explicit rather than asserted. The results are in principle falsifiable through the λ, h, σ dependence of T_rh and through the freeze-out boundaries. The principal physical caveat — rescattering and decay of the produced χ particles — is acknowledged in the Introduction (the paragraph beginning 'The opinion of the authors of Ref. 11…') and deferred to future work; the quantitative claims are contingent on it, and the authors say so.","major_comments":[{"comment":"The enhancement factor is not read off the mode equation directly; it follows from the backreaction closure in this section. Eq. (67) posits that the inflaton loses energy only through χ production at rate Γ, Eq. (68) then forces φ̇₀/φ₀ = −(H + ¼Γ), and Eq. (69) converts this into the resonance sweep rate dΔω_k/dt = ¼Γω_res; Eq. (70) is the resulting energy budget and Eq. (76) its inversion. This assumes (i) that the frequency–amplitude relation ω = c√λφ₀ (Eq. 31) survives backreaction, i.e. that the 2hφ⟨χ²⟩ force term in the inflaton equation is negligible compared with the Γ term, and (ii) that a single dominant band saturates the sum in Eq. (74). Neither assumption is checked against a direct solution of the coupled inflaton-plus-χ mode equations. Because the entire difference between the resonance rate and the Born rate is the logarithm in Eqs. (79)/(81), a factor-of-few error in the log argument propagates directly into Γ_χ^(res) and hence into T_rh (Eq. 101). Please justify (i) and (ii), or state explicitly that the quoted T_rh is an order-of-magnitude estimate.","section":"§5, Eqs. (67)–(70)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The perturbativity condition is printed as f², h, σφ₀ ≪ (ωφ₀)². With f and h dimensionless, σ of mass dimension one, and φ₀ of mass dimension one, the right-hand side has mass dimension four, so the inequality is not dimensionally consistent as written. The natural comparison for the hφ²χ² interaction is hφ₀² ≪ ω², which differs from the printed form by powers of φ₀. This matters because the perturbativity ranges (128) and (133) in §6 are derived from Eq. (44); if the condition is misstated, those ranges and the associated freeze-out bounds need rederivation. Please state the intended dimensionless combination and recheck (128) and (133) against it.","section":"§3, Eq. (44)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The O(1) coefficient inside the logarithm, 32π²/c⁴λ (and 512π²/c⁴λ), is not derived from first principles; it comes from the phase-space volume and the Gaussian integral in Eq. (72), and the prefactors in Eqs. (62) and (64) are likewise not tracked. The qualitative statement (enhancement ~1/λ) is robust to such errors, but the numerical value T_rh ~ 10¹² GeV (Eq. 102) is not. Please either derive the coefficient or present the result with an explicit uncertainty band, and state the expected accuracy of the phase-space estimate in Eq. (62).","section":"§5, Eqs. (72)–(76) and (83)–(84)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The freeze-out boundaries are obtained by substituting extreme values of H into the condition Γ > 4H: the lowest bound (127)/(132) uses the inflaton-only H from Eq. (12), while the highest (136)/(137) uses H ~ m at the end of inflation. The authors themselves note that the latter lies 'just on the boundary of the applicability of perturbation theory.' The bracketing logic is reasonable, but the paper should state clearly that (136)–(137) lie outside the perturbative regime in which the rest of the analysis is conducted, so that the bracket is not a controlled error bar. Please make this consequence explicit where the bounds are quoted.","section":"§6, Eqs. (126)–(137)","authors_acknowledge":null,"authors_acknowledge_p":null}],"minor_comments":[{"comment":"The inflaton field is written both as φ and as ϕ, sometimes within the same paragraph (e.g. the abstract versus §2). Please unify the notation.","section":"Throughout","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The text refers to 'q = 4 (quadric potential)'; this should read 'quartic potential.' Similar typographical slips occur elsewhere: 'oscllations' (§4), 'Wee see' (§5), 'not ony' (§7), 'coinsides' and 'procede' (§6), 'detailes' (Appendix B), and 'deser ves' (Introduction).","section":"§2, after Eq. (26)","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The authors state that Sections 2 and 3 closely follow Ref. 12 (Shtanov, Ukr. Fiz. Zh. 38, 1425, 1993, in Russian). Please state explicitly which results are new relative to that reference, so the reader can locate the novel contribution.","section":"§1 and §3","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The comparison with Ref. 13 (Kofman–Linde–Starobinsky) is qualitative and conducted in the notation of the present paper. A short table mapping the two sets of estimates (narrow versus broad resonance, the h ≃ 2λ/5 and h ≃ 0.22λ correspondences) would make the comparison much easier to follow.","section":"§7","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The abstract states that the reheating temperature 'is [λ log(1/λ)]⁻¹ times larger' than the Born estimate. Since the quoted T_rh carries order-unity prefactors (Eqs. 101, 115), please phrase this as an estimate or scaling rather than an exact ratio.","section":"Abstract, sentence [abstract:s5]","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The factor 1/2π² and the 1/2 in Eq. (62) are not explained. Please state the convention used for the phase-space volume so the reader can reproduce the coefficient.","section":"§4, Eq. (62)","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"major_revision","confidential_remarks":"The paper is a landmark contribution and the central qualitative claim — that parametric resonance dominates over the Born approximation in the perturbative regime, with an enhancement of order 1/λ — is well supported by the derivation and the two appendices. My recommendation of major revision is driven by two addressable items: the dimensional consistency of the perturbativity condition (Eq. 44), on which the ranges (128) and (133) depend, and the justification of the backreaction closure (Eqs. 67–70), on which the numerical enhancement factor and T_rh rest. I expect both can be handled in revision without altering the paper's conclusions. Note that the rescattering/decay caveat is explicitly acknowledged in the Introduction, so it should not be treated as an undisclosed limitation."},"author_rebuttal":null,"desk_editor":{"model":"selfhosted/deepseek-v4.1-flash","letter":"This is one of the papers that made preheating standard, and it earns that place. Traschen and Brandenberger (their Ref. 10) established parametric resonance in new inflation; this paper carries the analysis into chaotic inflation with λφ⁴ and ½m²φ², adds a backreaction treatment for the quartic case, and derives freeze-out boundaries for the quadratic case. The Bogoliubov machinery is set out properly in Appendices A and B — the first instability band, the growth rate μ₊, the adiabatic generalization to Eq. (B32) — so a reader can check the mode equation rather than take the result on faith. The scoping is honest: everything is claimed only where Eq. (44) holds, and the paper says plainly that perturbation theory breaks down at the onset of oscillations for the coupling ranges (128) and (133). It also discloses that the first version had errors corrected after comments from Kofman, Linde and Starobinsky. That is the right way to handle it. One correction to the record: the classifier reports no stated limitations, which is wrong. The Introduction explicitly flags rescattering and decay of the produced χ particles as an open question deferred to future work, and the Section 5 footnote says thermalization is assumed instantaneous. Both are acknowledged, not hidden. The sharper unaddressed item is the backreaction closure in Section 5. Eqs. (67)–(70) posit that the inflaton loses energy only through χ production at rate Γ, convert that into a resonance sweep rate dΔω_k/dt = ¼Γω_res, and invert the resulting energy budget in Eq. (76). Three things are assumed rather than checked: that ω = c√λφ₀ survives backreaction, that one dominant band saturates the sum in (74), and that the O(1) coefficient inside the log, 32π²/c⁴λ, is right. The qualitative enhancement ~1/λ is robust to all three; the quoted 10¹² GeV is not, since the log is the whole difference between the resonance and Born rates. Minor but real: the abstract claims both the lowest and highest freeze-out bounds, while the body attributes the highest (Eqs. 136–137) to Ref. 13. Who should read it: anyone working on reheating, preheating, or post-inflationary cosmology, and anyone who wants the perturbative-regime calculation done carefully. Send it to referees — the central claim holds as scoped, and the backreaction closure is exactly the kind of thing a good referee should press on.","headline":"A founding preheating paper: the parametric-resonance enhancement over Born estimates is real and well-scoped, but the backreaction closure that fixes the quoted numbers is an ansatz, not a derivation.","tokens_in":25592,"tokens_out":1027,"would_cite":true,"duration_ms":19238,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"pith_extraction":{"msc":[],"pacs":["98.80.Cq","98.80.Bp"],"model":"selfhosted/deepseek-v4.1-flash","headline":"Parametric resonance, not ordinary perturbative decay, controls how the universe reheats after inflation, and the correction can raise the reheating temperature by many orders of magnitude.","keywords":["parametric resonance","reheating","chaotic inflation","inflaton","Bogoliubov transformation","narrow resonance","inflaton freeze-out","inflaton dark matter"],"falsifier":"A coupled lattice or classical-statistical simulation of the inflaton together with the produced scalar, including rescattering and any $chi$ self-interaction, that tracks the occupation number of a mode as it passes through a resonance band: if $N_k$ saturates near unity instead of growing as $sinh^2(int mu_+ dt)$, the enhancement factors of Eqs. 83-84 and the reheating temperatures of Eq. 115 are wrong. Conversely, if the simulated $N_k$ follows the predicted exponential growth through the band, the central claim stands.","tokens_in":25761,"feed_emoji":"🔥","tokens_out":3066,"duration_ms":58867,"temperature":0.7,"pith_summary":"The paper asks how the cold, empty universe produced by inflation gets filled with particles. The standard answer computed the rate at which an oscillating inflaton field decays into fermions and scalars using first-order Born approximation. The authors show that for the scalar channel this misses the dominant physics: the oscillating inflaton acts as a periodic drive on the mode equation for the produced field $k$, and, exactly as in a parametrically driven oscillator, modes whose frequency sits near half-integer multiples of the inflaton frequency, $omega_k approx (n/2),omega$, grow exponentially rather than linearly in time. In the quartic case $V(varphi)=lambda varphi^4$ this makes reheating efficient even for tiny couplings, and the resulting temperature exceeds the Born estimate by a factor of order $[lambda log(1/lambda)]^{-1}$. In the quadratic case $V(varphi)=tfrac12 m^2varphi^2$ the resonance is narrow and shuts off as the universe expands, so the paper instead derives freeze-out boundaries below which particle creation becomes inefficient and the residual inflaton oscillations can survive as cold dark matter.","feed_headline":"Reheating after inflation runs on resonance, not decay","feed_subtitle":"Ignoring parametric resonance can underestimate the post-inflation reheating temperature by more than ten orders of magnitude.","key_machinery":"The load-bearing device is the reduction of the produced-field mode equation to a Hill/Matthiue-type oscillator, $ddot{Y}_k + (omega_k^2(t) + g(omega t))Y_k = 0$, with $g$ a periodic drive containing harmonics at $omega$ (from the linear coupling $sigma varphi$) and $2omega$ (from the quadratic coupling $h varphi^2$). Bogoliubov averaging in the small drive amplitude gives the first instability band, its width $|g_n|$, and the exponential rate $mu_+$; the occupation number follows from $N_k simeq sinh^2(int mu_+ dt)$ with the integral taken over the time a mode spends inside the band. The second device is the bookkeeping that turns this into a cosmological rate: in the quartic case, energy conservation ties the inflaton's amplitude loss to the drift speed $|dDeltaomega_k/dt|$ through the band, so the resonance rate and the back-reaction on the inflaton must be solved simultaneously. That closure is what converts a statement about instability bands into a number for $T_{rh}$.","core_discovery":"The central claim is that within the regime where perturbation theory in the inflaton couplings is legitimate, parametric resonance dominates scalar particle production and first-order Born estimates systematically underestimate it. The authors reduce the mode equation for the produced scalar $chi$ to a periodically driven oscillator, identify the instability bands $omega_k^2 - (nomega/2)^2 < |g_n|$, and compute the growth exponent $mu_+ = (1/nomega)sqrt{|g_n|^2 - Delta_n^2}$, so that the mean occupation number of a mode swept through a band is $N_k simeq sinh^2(int mu_+ dt)$. For $V=lambda varphi^4$ the resonance frequency redshifts at the same rate as the mode frequencies, so the authors close the system with a back-reaction equation: energy lost by the inflaton fixes the rate at which modes drift through the band, which in turn fixes $N_k$, giving an implicit equation for the effective production rate $Gamma_chi^{(res)}$ (their Eq. 74). The resulting enhancement over Born approximation is a factor $sim 32pi^2 c^{-2}lambda^{-1}log^{-1}(512pi^2 c^{-4}lambda^{-1})$ for the linear coupling and $sim pi^2 c^{-2}lambda^{-1}log^{-1}(32pi^2 c^{-4}lambda^{-1})$ for the quadratic coupling, which for $lambda sim 10^{-12}$ is enormous. The reheating temperature follows as $T_{rh} sim lambda^{1/4}phi_0$ with $phi_0$ at the onset of efficient reheating, giving $T_{rh} sim 10^{-1}M_P lambda^{1/4}(h/lambda)^2 log^{-1}(1/lambda)$, roughly $10^{12}$ GeV for the illustrative couplings. For the quadratic potential the frequency is constant, the drift rate is set by Hubble expansion rather than back-reaction, and the same machinery instead yields the amplitude at which production switches off: $(phi_0/m)^3 gtrsim 65(m/h^2 M_P)log(m/h^2 M_P)$ for quadratic coupling and $(phi_0/m) gtrsim (m^3/sigma^2 M_P)log(m^3/sigma^2 M_P)$ for linear coupling, with the minimum possible freeze-out energy density of order $m^4$.","pith_inferences":["If resonance really drives reheating in the quartic model, the produced particles arrive in narrow, highly occupied bands rather than in a thermal spectrum, which would leave a non-thermal transient and possibly a distinctive gravitational-wave or fluctuation signature; the paper does not compute one.","The freeze-out boundaries for the quadratic potential give a concrete handle on inflaton dark matter: the relic abundance depends on where between the lowest and highest boundary the inflaton actually freezes out, which is exactly the range the paper leaves undetermined.","The narrow-resonance treatment here and the broad-resonance treatment of the cited concurrent work should connect continuously as the coupling grows; a calculation that interpolates between them would test whether the enhancement factors quoted here are the small-coupling limit of a single curve.","Because the whole result rests on occupation numbers growing before anything depletes them, a coupled lattice simulation of the inflaton and the produced field would be the natural arbiter, and would also indicate how far into the reheating epoch the analytic rate survives."],"forward_implications":["For $V=lambda varphi^4$, reheating is efficient for arbitrarily small inflaton couplings, so the quartic chaotic model does not suffer the reheating problem that Born-approximation estimates suggested; the temperature scales as $lambda^{1/4}(h/lambda)^2 log^{-1}(1/lambda)$ rather than as a much smaller Born rate.","The reheating temperature inferred from any model with a quartic inflaton potential should be revised upward by the factor $[lambda log(1/lambda)]^{-1}$ relative to first-order calculations, which for $lambda sim 10^{-12}$ is a difference of many orders of magnitude.","For $V=tfrac12 m^2varphi^2$, there is a definite amplitude below which particle production stops competing with Hubble redshift; the inflaton then decouples and its residual oscillations can supply cold dark matter, with the freeze-out energy density bracketed between about $m^4$ and the highest possible boundary set by the Hubble scale at the end of inflation.","Fermions are unaffected: the resonance mechanism has no fermionic analogue, so the Born result $Gamma_psi = O(1) f^2 omega/(4pi)$ and the reheating condition $phi_0 lesssim 10^{-2}f^2 M_P$ stand as previously computed.","Which channel dominates reheating depends on whether $phi_0$ is above or below $sigma/h$, giving a sharp crossover in the predicted reheating temperature as a function of the two scalar couplings."],"supporting_citations":[{"why":"Traschen and Brandenberger established that parametric resonance, not Born-level decay, governs boson production by an oscillating inflaton in the new-inflation context; this paper extends that mechanism to chaotic inflation with power-law potentials.","marker":"[10]"},{"why":"Dolgov and Kirilova supply the result that fermion production has no resonance enhancement, which the paper adopts to keep the fermionic channel at Born order, and also the concern about rescattering and decay that the authors explicitly defer.","marker":"[11]"},{"why":"Kofman, Linde and Starobinsky's concurrent work defines the narrow-versus-broad resonance split and the inflaton-as-dark-matter scenario; the paper compares its narrow-resonance numbers and freeze-out boundaries against theirs.","marker":"[13]"},{"why":"Shtanov's earlier treatment supplies the Born-approximation production rates and inflaton dynamics that Sections 2 and 3 reproduce and then improve upon.","marker":"[12]"},{"why":"Bogoliubov and Mitropolsky's averaging method is the technical basis for the instability-band analysis in Appendix A that yields the growth exponent.","marker":"[14]"},{"why":"Linde's review supplies the chaotic inflation framework and the benchmark coupling $lambda sim 10^{-12}$ used to convert the enhancement factors into a numerical reheating temperature.","marker":"[1]"}],"fun_headline_variants":["Parametric resonance drives reheating when perturbation theory holds","Quartic inflaton reheating outpaces first-order Born estimates","Quadratic inflaton reheating depends on coupling type and strength","Suppressed couplings may let the inflaton freeze out as dark matter","Resonance bands, not simple decay, set scalar production rates"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The calculation assumes that nothing removes particles from a resonantly amplified mode before the mode leaves the resonance band; the authors state in the Introduction, in the paragraph beginning 'The opinion of the authors of Ref. 11 is that parametric resonance effects of scalar particle production are likely to be suppressed by the processes of scattering and decay,' that rescattering and decay of the produced $chi$ particles are not treated here and are left to future work, so the enhancement factors and reheating temperatures hold only if those processes are slower than the resonance.","fun_headline_variants_meta":{"raw":{"variants":["Parametric resonance drives reheating when perturbation theory holds","Quartic inflaton reheating outpaces first-order Born estimates","Quadratic inflaton reheating depends on coupling type and strength","Suppressed couplings may let the inflaton freeze out as dark matter","Resonance bands, not simple decay, set scalar production rates"]},"model":"selfhosted/deepseek-v4.1-flash","effort":"low","attempts":1,"cost_usd":0.0,"raw_usage":{"total_tokens":2330,"prompt_tokens":1335,"completion_tokens":995,"prompt_tokens_details":null,"completion_tokens_details":{"reasoning_tokens":922}},"tokens_in":1335,"tokens_out":995,"duration_ms":15750,"image_pages":[],"temperature":1.0,"reasoning_budget":50,"reasoning_tokens":922,"transport_provider":"selfhosted","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-10-29T19:38:11.477015+00:00","model_set":{"reader":"selfhosted/deepseek-v4.1-flash"},"falsifier":"A coupled lattice or classical-statistical simulation of the inflaton together with the produced scalar, including rescattering and any $chi$ self-interaction, that tracks the occupation number of a mode as it passes through a resonance band: if $N_k$ saturates near unity instead of growing as $sinh^2(int mu_+ dt)$, the enhancement factors of Eqs. 83-84 and the reheating temperatures of Eq. 115 are wrong. Conversely, if the simulated $N_k$ follows the predicted exponential growth through the band, the central claim stands.","supporting_citations":[{"cited_title":"Phys. Rev. D42, 2491","cited_arxiv_id":null,"evidence_quote":"Traschen and Brandenberger established that parametric resonance, not Born-level decay, governs boson production by an oscillating inflaton in the new-inflation context; this paper extends that mechanism to chaotic inflation with power-law potentials."},{"cited_title":"Sov. Nucl. Phys. 51, 273","cited_arxiv_id":null,"evidence_quote":"Dolgov and Kirilova supply the result that fermion production has no resonance enhancement, which the paper adopts to keep the fermionic channel at Born order, and also the concern about rescattering and decay that the authors explicitly defer."},{"cited_title":"UH-IfA-94/35, SU -ITP-94-13, YITP/U- 94-15, and HEP-TH-9405187, Submitted to: Phys. Rev. Lett","cited_arxiv_id":null,"evidence_quote":"Kofman, Linde and Starobinsky's concurrent work defines the narrow-versus-broad resonance split and the inflaton-as-dark-matter scenario; the paper compares its narrow-resonance numbers and freeze-out boundaries against theirs."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shtanov's earlier treatment supplies the Born-approximation production rates and inflaton dynamics that Sections 2 and 3 reproduce and then improve upon."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Bogoliubov and Mitropolsky's averaging method is the technical basis for the instability-band analysis in Appendix A that yields the growth exponent."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Linde's review supplies the chaotic inflation framework and the benchmark coupling $lambda sim 10^{-12}$ used to convert the enhancement factors into a numerical reheating temperature."}],"review_version":1}