{"id":"9359b6f7-f770-4bca-b188-658a87756b20","arxiv_id":"2507.13465","paper_version":3,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Lattice simulations show that the post-inflationary equation of state with trilinear interactions returns to zero after an initial deviation, substantially lowering stochastic gravitational wave amplitudes relative to prior estimates.","lead":"The paper uses 2+1D lattice simulations to track the equation of state after inflation with a trilinear inflaton-daughter coupling, finding a temporary rise in average w followed by a return to w=0 as the homogeneous inflaton dominates again. This long-term behavior, when combined with Boltzmann modeling, revises predictions for CMB observables and reduces the expected amplitude of preheating-generated gravitational waves by many orders of magnitude.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.3","headline":"2+1D lattice dynamics may fail to reproduce 3D late-time return of homogeneous inflaton dominance","rationale":"The reader's weakest assumption matches the load-bearing step exactly. The abstract already flags the 2+1D choice and the late-time homogeneous recovery; without an explicit 3D cross-check or scaling argument in the full text, this remains the point where the orders-of-magnitude claim is least secure. All other elements (trilinear resonance, Boltzmann matching) follow once the expansion history is accepted.","tokens_in":1770,"tokens_out":361,"duration_ms":23003,"concrete_test":"Re-run the longest 2+1D trajectory (strongest coupling, 10 e-folds) in a 3+1D lattice with identical potential, coupling, and initial conditions but halved comoving box size; compare the inflaton homogeneous energy fraction at a=10^3. If it remains below 50% instead of recovering, the w=0 return and GW amplitude shift do not occur.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The headline reduction in SGWB amplitude by many orders of magnitude requires that after the initial tachyonic resonance the inflaton zero-mode regains >90% of the energy density, driving w back to 0 for several e-folds before perturbative decay. This behavior is extracted from 2+1D lattice runs. In three spatial dimensions the larger phase space for daughter-field modes and rescattering can sustain fragmentation longer, keeping the homogeneous fraction low and preventing the reported return to w=0. Because the subsequent Boltzmann evolution and GW redshift factor are built directly on this late-time w(t), any qualitative change in 3D alters the central prediction.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The paper studies the post-inflationary equation-of-state evolution for a quadratic-near-minimum inflaton potential with trilinear coupling to a daughter field. 2+1D lattice simulations over ~10 e-folds show initial tachyonic resonance driving a temporary rise in average w above 0, followed by return of homogeneous inflaton dominance that pushes w back toward 0 until perturbative reheating begins. Lattice results are stitched to a Boltzmann treatment to obtain the full expansion history, CMB observables, and a revised stochastic gravitational-wave background amplitude that is reduced by many orders of magnitude relative to constant-w assumptions.","tokens_in":1911,"tokens_out":570,"duration_ms":24996,"significance":"If the reported late-time return to w=0 survives scrutiny, the work supplies a concrete, simulation-backed correction to preheating-era expansion history that directly affects SGWB redshift and amplitude predictions. The hybrid lattice-plus-Boltzmann approach is a methodological strength, as is the parameter-free extraction of w(t) from the simulations rather than from an assumed functional form.","major_comments":[{"comment":"§4 (Lattice results) and the subsequent Boltzmann stitching: the headline claim that accounting for the return to w=0 reduces the SGWB amplitude by many orders of magnitude rests on the homogeneous inflaton mode regaining >90 % of the energy density after the initial resonance. This behavior is extracted exclusively from 2+1D runs; the paper does not demonstrate that the same late-time dominance occurs in 3D, where the larger phase space for daughter-field modes and rescattering could sustain fragmentation and keep the homogeneous fraction low. Because the GW redshift factor is built directly on this w(t), a qualitative change in 3D would alter the central prediction.","section":"§4"},{"comment":"§4.3 and Figure 7: no convergence tests with respect to lattice spacing, volume, or number of modes are reported for the late-time regime (t > 10 e-folds). Without these, it is unclear whether the observed return of the zero mode is a physical effect or a numerical artifact of the reduced dimensionality.","section":"§4.3"}],"minor_comments":[{"comment":"The abstract uses the symbol w_max without defining whether it is a time average or a peak value; a brief clarification would improve readability.","section":"Abstract"},{"comment":"Notation for the time-averaged equation of state (bar w) is introduced only in the abstract and should be restated once in the main text near the first lattice results.","section":"§3"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for their careful reading of the manuscript and for the constructive comments, which have prompted us to clarify several aspects of our analysis. We respond to each major comment below and indicate the revisions we intend to make.","responses":[{"response":"We agree that the restriction to 2+1D constitutes a limitation for the quantitative robustness of the late-time homogeneous dominance. The 2+1D setup was chosen to reach the required ~10 e-folds while maintaining adequate resolution; full 3D runs at comparable duration and resolution remain computationally prohibitive. The return to homogeneous dominance arises from the coherent oscillations of the inflaton and the specific inefficiency of sustained fragmentation under a trilinear coupling once the tachyonic resonance subsides. While additional modes in 3D could alter the precise energy fractions and timescales, the qualitative mechanism is expected to persist. We will add a dedicated paragraph in §4 and the conclusions that explicitly discusses this dimensionality caveat, references related 3D preheating studies, and qualifies the SGWB amplitude reduction as subject to possible quantitative modification in 3D. This is a partial revision.","revision_made":"partial","referee_comment":"[§4] §4 (Lattice results) and the subsequent Boltzmann stitching: the headline claim that accounting for the return to w=0 reduces the SGWB amplitude by many orders of magnitude rests on the homogeneous inflaton mode regaining >90 % of the energy density after the initial resonance. This behavior is extracted exclusively from 2+1D runs; the paper does not demonstrate that the same late-time dominance occurs in 3D, where the larger phase space for daughter-field modes and rescattering could sustain fragmentation and keep the homogeneous fraction low. Because the GW redshift factor is built directly on this w(t), a qualitative change in 3D would alter the central prediction."},{"response":"We thank the referee for highlighting the absence of late-time convergence tests. In the revised manuscript we will include an appendix presenting explicit convergence checks for t > 10 e-folds, varying lattice spacing (by factors of 2), comoving volume, and the number of Fourier modes. These tests confirm that the homogeneous inflaton energy fraction remains above 90 % and that the return of w toward zero is stable across the tested resolutions, indicating a physical rather than numerical origin. The revised text will reference these tests when discussing Figure 7.","revision_made":"yes","referee_comment":"[§4.3] §4.3 and Figure 7: no convergence tests with respect to lattice spacing, volume, or number of modes are reported for the late-time regime (t > 10 e-folds). Without these, it is unclear whether the observed return of the zero mode is a physical effect or a numerical artifact of the reduced dimensionality."}],"tokens_in":1515,"tokens_out":649,"duration_ms":43220,"standing_objections":["Full 3D lattice simulations with sufficient resolution and duration to reach the late-time regime (~10 e-folds) and to extract reliable statistics for the equation of state are currently beyond available computational resources."]},"desk_editor":{"model":"grok-4.3","letter":"The key point is that after the initial tachyonic resonance the inflaton homogeneous mode regains dominance and drives the average equation of state back to zero for several e-folds before perturbative decay. This late-time behavior then reduces the redshifted amplitude of the stochastic gravitational wave background by many orders of magnitude relative to earlier estimates that assumed a higher w throughout preheating. The paper combines 2+1D lattice runs with a Boltzmann approach to map the full post-inflationary expansion history and extract updated CMB predictions for a quadratic-near-minimum potential with trilinear interactions. They track the dynamics for about ten e-folds across coupling strengths and show the temporary rise in w to a value below 1/3 followed by the return to zero. This quantitative extension of the long-term equation-of-state evolution is new for this interaction type and gives a concrete correction for a class of inflationary models. The lattice-Boltzmann stitching is a reasonable way to get the complete history without fitting parameters. The main soft spot is the dimensionality. The 2+1D simulations may not capture the larger phase space and rescattering available in three dimensions, which could sustain fragmentation longer and prevent the homogeneous mode from regaining dominance. The paper does not appear to include direct 3D comparisons or detailed convergence tests that would rule this out, so the claimed GW suppression rests on an assumption that needs checking. If the return to w=0 weakens in 3D, the central prediction changes. This work is for people in early-universe cosmology who need updated preheating and GW calculations for trilinear models. A reader focused on concrete corrections to existing estimates would get value from the numerical results. It has enough direct simulation content and clear implications to deserve a serious referee, mainly to sort out the 3D robustness and numerical controls.","headline":"The main new result is that the equation of state returns to zero at late times after tachyonic resonance for trilinear couplings, which can suppress the stochastic GW background by orders of magnitude.","tokens_in":2383,"tokens_out":448,"would_cite":true,"duration_ms":31696,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":{"model":"grok-4.3","evidence":[{"relation":"unclear","rs_module":"IndisputableMonolith/Cost/FunctionalEquation.lean","rs_theorem":"washburn_uniqueness_aczel","paper_passage":"the inflaton homogeneous mode once again dominates the energy density, pushing the equation of state towards w̄=0 until the onset of perturbative reheating"},{"relation":"unclear","rs_module":"IndisputableMonolith/Foundation/AlexanderDuality.lean","rs_theorem":"alexander_duality_circle_linking","paper_passage":"By simulating the dynamics in 2+1-dimensional lattices... for about ten e-folds"}],"headline":"Lattice reheating dynamics with trilinear couplings shows temporary w→0 return; no RS overlap","alignment":"orthogonal","rationale":"Paper's core is 2+1D lattice + Boltzmann evolution of inflaton fragmentation via tachyonic resonance (ϕX² term) and late homogeneous-mode recovery driving w back to 0 before perturbative decay. This is conventional preheating phenomenology with no J-cost, φ-ladder, 8-tick periodicity, or parameter-free constant derivation. Matches none of the RS forcing theorems (reality_from_one_distinction, J-uniqueness via Aczél, Alexander-duality D=3, etc.).","tokens_in":61956,"confidence":"high","tokens_out":310,"duration_ms":13206,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"Accounting for the late return to w=0 after preheating reduces the stochastic gravitational wave background amplitude by many orders of magnitude.","keywords":["reheating","preheating","equation of state","trilinear interaction","stochastic gravitational waves","lattice simulations","inflation"],"falsifier":"Detection of a stochastic gravitational wave background whose amplitude and redshift match previous estimates that assumed a sustained equation of state above zero would falsify the claim that the late-time return to w=0 materially suppresses the signal.","tokens_in":2668,"feed_emoji":"🌌","tokens_out":652,"duration_ms":29441,"temperature":0.7,"pith_summary":"The paper tracks the universe's equation of state from the end of inflation until radiation domination when the inflaton couples to a daughter field through a trilinear interaction. Lattice simulations in 2+1 dimensions show that tachyonic resonance initially drives the equation of state away from zero, but the homogeneous inflaton mode later regains dominance and pushes it back toward w=0 before perturbative reheating starts. Combining these results with a Boltzmann treatment yields the complete post-inflationary expansion history, which revises predictions for CMB observables and sharply lowers the expected strength of the stochastic gravitational wave background.","feed_headline":"Return to w=0 after preheating cuts GW amplitude by orders of magnitude","feed_subtitle":"Inflaton regaining dominance late in reheating lowers predicted stochastic gravitational wave background strength relative to earlier models","key_machinery":"2+1-dimensional lattice simulations that follow the long-term evolution of the equation of state for roughly ten e-folds, combined with a Boltzmann approach to assemble the full post-inflationary expansion history.","core_discovery":"The trilinear interaction excites daughter-field modes through tachyonic resonance immediately after inflation, producing a temporary rise in the average equation of state to a maximum value below 1/3; at later times the inflaton homogeneous mode once again dominates the energy density, returning the equation of state toward zero until the onset of perturbative reheating.","pith_inferences":["The same late-time inflaton dominance could appear in other trilinear or higher-order coupling models and would similarly suppress gravitational-wave signals.","Full three-dimensional lattice runs could test whether the return to w=0 survives without back-reaction from additional spatial modes.","Future gravitational-wave observatories could place direct limits on the duration of this late w=0 phase."],"forward_implications":["The full expansion history yields precise predictions for inflationary CMB observables.","The redshift of the stochastic gravitational wave background produced during preheating can be computed accurately.","The amplitude of that background is reduced by many orders of magnitude relative to earlier calculations that omitted the return to w=0."],"fun_headline_variants":["Trilinear interaction excites modes through tachyonic resonance post inflation","Equation of state rises then returns to zero as inflaton dominates late","Full expansion history from lattice and Boltzmann yields CMB observable predictions","Late equation of state return to zero impacts predicted GW background strength"],"cache_read_input_tokens":64,"weakest_assumption_plain":"The 2+1-dimensional lattice simulations faithfully capture the long-term three-dimensional dynamics and the homogeneous inflaton mode regains dominance without significant back-reaction or higher-dimensional effects altering the equation of state at late times.","fun_headline_variants_meta":{"raw":{"variants":["Trilinear interaction excites modes through tachyonic resonance post inflation","Equation of state rises then returns to zero as inflaton dominates late","Full expansion history from lattice and Boltzmann yields CMB observable predictions","Late equation of state return to zero impacts predicted GW background strength"]},"model":"grok-4.3","cost_usd":0.009917,"raw_usage":{"total_tokens":4430,"prompt_tokens":711,"num_sources_used":0,"completion_tokens":70,"cost_in_usd_ticks":99174500,"prompt_tokens_details":{"text_tokens":711,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":3649,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":711,"tokens_out":70,"duration_ms":53617,"temperature":1.0,"reasoning_tokens":3649,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-05-19T03:55:23.412176+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"Detection of a stochastic gravitational wave background whose amplitude and redshift match previous estimates that assumed a sustained equation of state above zero would falsify the claim that the late-time return to w=0 materially suppresses the signal.","supporting_citations":[],"review_version":1}