{"id":"9a937aeb-c8b9-4896-ad0d-dd06be166c42","arxiv_id":"2607.27325","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"In seesaw reheating, the post-inflation universe can pass through four alternating matter/radiation eras, with the Standard Model temperature falling as a^{-1/4} and then a^{-3/8}, which changes dark-matter production.","lead":"This paper analyzes a two-step reheating after inflation: the inflaton decays into heavy right-handed neutrinos, which later decay into ordinary particles. It derives new temperature scalings for the intermediate eras and shows they can boost or alter dark-matter production.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The four-stage history and the DM predictions rest on the unquantified assumption that RHNs are collisionless; if seesaw Yukawa scatterings or inverse decays thermalize the RHN population, the production-time-dependent scalings (T∝a^{-1/4}, a^{-3/8}; p=12,20) fail.","rationale":"The paper presents a coherent integro-differential system and the analytic solutions are internally consistent and match the numerical plots. After checking the heavy-tier derivations, I find no algebraic error; the stage equations (4.1)–(4.13) and the temperature history (4.14) follow. The most load-bearing (and indeed the reader's) concern is the collisionless-RHN assumption stated in Section 3. The entire novelty—the production-time distribution, the alternating ω_eff, the a^{-1/4}/a^{-3/8} scalings, and the p=12,20 freeze-in powers—depends on RHNs not scattering before decay. The paper does not quantify the scattering/inverse-decay rates, so a skeptical reader cannot tell whether the regime is actually realizable or whether it collapses into the prompt-thermalization picture of [29]. I agree with the reader's identification; however, because the Yukawa couplings are tiny in the seesaw parameter space and the paper's benchmark is plausibly far below the scattering horizon, I do not see this as a fatal flaw. The concrete check (scattering rate vs. H for the benchmark and the claimed region) would settle the matter. Thus the ACCEPT verdict stands, and I would not adjust it.","tokens_in":16689,"tokens_out":24795,"duration_ms":190538,"concrete_test":"Compute the dominant RHN thermalization rates from the seesaw Lagrangian (2.1) — elastic scattering N+ℓ↔N+ℓ and inverse decay ℓ+H↔N — using the Casas-Ibarra couplings (2.5), and compare with H(a) over the interval a_φ ≲ a ≲ a_rh for the benchmark of Fig. 1 (H_I=10^9 GeV, m_φ=3.3×10^9 GeV, m_N=1.7×10^6 GeV, Γ_φ=2.5×10^6 GeV, Γ_N=3.1×10^{-5} GeV). If Γ_scat/H ≥ 1 anywhere, the collisionless assumption fails and the T(a) scalings and p=12,20 predictions must be revised; if Γ_scat/H ≪ 1 throughout, the central claim is supported. Also scan the parameter region where the four-stage history is claimed to confirm the condition is generic.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3 (after Eq. 3.9) treats the nonthermal RHN population as collisionless between production and decay, neglecting elastic scattering, number-changing reactions, and inverse processes. This approximation is the load-bearing pillar of the entire paper: the production-time-resolved distribution δn_N(a,a') (Eq. 3.9) and the resulting energy/pressure integrals (3.13–3.14) are what generate the alternating ω_eff and the T(a) scalings in Eq. (1.3)/(4.14). If RHNs scatter or thermalize efficiently, the distribution collapses to a thermal one, and the history reduces to the prompt-thermalization treatment of [29]; the a^{-1/4}/a^{-3/8} scalings and the critical powers p=12,20 in Section 5.2 would not apply. The paper gives no quantitative estimate of the RHN scattering or inverse-decay rate relative to H in the parameter region where the four-stage sequence is claimed (Section 4; benchmark in Fig. 1). The Yukawa couplings that set Γ_N also control these scatterings; while they are small in the seesaw regime, the condition Γ_scat ≪ H during a_φ ≲ a ≲ a_rh is not derived or checked. Without that check, the central claim is contingent on an unverified assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a reheating history in which the inflaton decays exclusively into long-lived right-handed neutrinos (RHNs), which subsequently decay into the Standard Model. The authors solve the coupled Boltzmann system while retaining the production-time dependence of the nonthermal RHN distribution, including the relativistic-to-nonrelativistic transition and the time-dilated decay rate of each production cohort. They find an alternating effective equation of state (inflaton matter, RHN radiation, RHN matter, SM radiation) and derive piecewise temperature scalings T(a) proportional to a^0, a^{-1/4}, a^{-3/8}, and a^{-1} in the four eras. They apply this thermal history to dark matter from direct inflaton decay, obtaining an enhancement of order m_phi/(2 m_N) at fixed reheating temperature and branching fraction, and to UV freeze-in, where they identify critical powers p=12 and p=20 that control how much of the DM is produced in the intermediate RHN-dominated eras. The benchmark numerical evolution in Fig. 1 is consistent with the analytic estimates.","tokens_in":17149,"tokens_out":25597,"duration_ms":204188,"significance":"If the underlying assumptions hold, this is a qualitatively new pre-BBN thermal history that connects type-I seesaw neutrino masses to the expansion and temperature evolution before nucleosynthesis. The paper contains a genuinely useful technical step: the exact source identity in Eq. (3.11), which removes momentum dependence from the SM radiation source, is elegant and correctly implemented. The analytic scalings are derived, not fitted, and the benchmark numerics confirm the four-stage sequence. The DM implications are concrete and falsifiable: the p=12 and p=20 critical powers give sharply different predictions from standard reheating and from the near-simultaneous prompt-thermalization treatment in Ref. [29]. The main weakness is that the entire production-time-resolved picture is conditional on the RHNs being collisionless from production to decay, and this condition is asserted but not quantitatively established.","major_comments":[{"comment":"The paper states that RHNs are 'collisionless between their production and decay, neglecting elastic scattering, number-changing reactions, and inverse inflaton decays.' This assumption is load-bearing: the production-time-dependent distribution delta n_N(a,a') in Eq. (3.9), the energy/pressure integrals in Eqs. (3.13)-(3.14), and the resulting temperature scalings in Eq. (1.3)/(4.14) all rely on it. If RHN scattering or inverse processes are efficient enough to redistribute momenta or change the RHN number density, the distribution collapses toward the prompt-thermalization/characteristic-momentum treatment of Ref. [29] and the claimed departures disappear. The paper does not estimate any relevant rate (e.g., NN->lH, lH->N, or pair annihilation) relative to H in the region where the four-stage sequence is claimed. For the benchmark parameters of Fig. 1, the minimal Casas-Ibarra Yukawa c","section":"Section 3, Eq. (3.9) paragraph"},{"comment":"The text claims that inverse inflaton decay (production of inflatons from two RHNs) is 'always kinematically blocked.' This is stronger than what follows from the stated assumptions. For two RHNs produced from different inflaton decays, the invariant mass s of the pair is not fixed; after redshift, a continuum of pair energies is present, and pairs with s = m_phi^2 are kinematically allowed. The total 3-momentum of such a pair is generally nonzero, so the final inflaton has nonzero momentum, but the process is not forbidden. Whether it is negligible is a dynamical question controlled by the small inflaton-RHN coupling and the phase-space density of resonant pairs. The statement as written is not justified and should either be corrected or replaced by a quantitative bound on the rate of NN -> phi relative to H.","section":"Section 3, paragraph after Eq. (3.9)"}],"minor_comments":[{"comment":"The phrase 'with n′ referring to the N_i fulfilling...' is awkwardly typeset; please clarify that the sum is over RHN species with M_i < m_phi/2.","section":"Section 2, Eq. (2.2)"},{"comment":"The header 'T able 1' contains a formatting typo; should be 'Table 1'.","section":"Table 1"},{"comment":"The definition of a_nr in Eq. (4.9) uses the highest-momentum cohort produced at a_phi. This is a clean asymptotic criterion, but before a_nr some older cohorts are already nonrelativistic and after a_nr the highest-momentum cohort is not instantly cold. The analytic scalings are therefore only strictly valid in the limit R2 >> 1. The paper already requires well-separated regimes in Section 4, but it would be helpful to state this explicitly near Eq. (4.9) and to quantify the residual correction.","section":"Section 4, transition at a_nr"},{"comment":"In Eq. (5.20) and the surrounding text, g_star and g_star s should be evaluated at T_rh and taken constant across the eras; this is a standard approximation but should be stated once for clarity.","section":"Section 5.2"},{"comment":"The note comparing Eq. (4.16) with Ref. [29] is useful. It would strengthen the paper to state explicitly which predictions beyond T_rh differ from the prompt-thermalization treatment and to emphasize that the difference is controlled by the collisionless assumption.","section":"Note added"}],"recommendation":"major_revision","confidential_remarks":"The paper is physically interesting and technically sound conditional on the collisionless-RHN approximation. My main concern is that this approximation is central to every new result and is not quantified. I would be comfortable with acceptance after the authors add a dedicated estimate of the relevant RHN scattering and inverse-decay rates relative to the Hubble rate, specifically in the benchmark parameter region and with the Casas-Ibarra parametrization. The kinematic-blocking claim for NN -> phi should also be corrected. The overlap with Ref. [29] is honestly disclosed, but the distinguishing claim depends on the same assumption and therefore needs the same scrutiny."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. The production-time-resolved treatment of the nonthermal RHN distribution is the real novelty, and the T(a) scalings (plateau, a^{-1/4}, a^{-3/8}, a^{-1}) are clean results that follow from the stated equations. But the whole four-stage history and the DM predictions sit on the assumption that the RHNs are collisionless between production and decay, and the paper does not quantify that.\n\nThe stress-test note is on target: Section 3 just states the assumption, with no estimate of elastic scattering, number-changing reactions, or inverse decays relative to H. In the tiny-Yukawa seesaw regime those rates are probably small, but the paper should show it. I would add a second, more direct worry: the neglect of Pauli blocking in the inflaton decay itself. The paper says Γ_phi << H makes Pauli blocking negligible, but that is not the right criterion. A steady-state estimate for the benchmark gives a phase-space occupancy f ~ Γ_phi n_phi/(H p^3) >> 1 (many orders of magnitude). If that is right, the perturbative decay rate is suppressed by (1-f)^2 and the RHN energy density would not reach the values in Fig. 1. The authors should check this; it may force them to heavier inflaton masses or multiple RHN species. This issue is more serious than the scattering one, because it affects the validity of the background equations themselves.\n\nWhat is genuinely good: the integro-differential system in Eqs. (3.1)-(3.15) is coherent, and the source identity (3.11) is elegant. The analytic stage solutions match the numerical evolution in Fig. 1. The DM applications are carefully derived: the m_phi/(2m_N) enhancement for direct decay and the p=12,20 critical powers for UV freeze-in follow from the stated background. The note added discloses the overlap with [29] and the agreement in Eq. (4.16). The literature review is honest and the free parameters are clearly identified.\n\nMinor weaknesses: the abstract promises baryogenesis but no computation is given; the analytic stage solutions are stated without full derivations; no code is shipped. None of these are load-bearing.\n\nWho gets value: pre-BBN cosmologists and DM model-builders working on RHN-mediated reheating. The paper deserves a serious referee, but the referee should push for a quantitative analysis of the collisionless assumption and the Pauli-blocking consistency of the benchmark. As it stands, the central scenario is promising but resting on an unverified pillar.","headline":"Elegant two-step reheating formalism with clean scalings, but the load-bearing collisionless-RHN assumption is unchecked and a Pauli-blocking estimate calls the benchmark into question.","tokens_in":17596,"tokens_out":23082,"would_cite":true,"duration_ms":200143,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"When inflaton energy passes through long-lived seesaw neutrinos, the Universe cycles matter→radiation→matter→radiation and the Standard Model temperature cools with the unusual powers a^-1/4 and a^-3/8, reshaping dark-matter production.","keywords":["seesaw cosmology","right-handed neutrinos","reheating","nonthermal distribution","equation of state","dark matter production","UV freeze-in","pre-BBN thermal history"],"falsifier":"If one includes even a small elastic-scattering or number-changing interaction among RHNs (or inverse decays that repopulate the inflaton), the distribution in Eq. (3.9) no longer holds. A Boltzmann simulation with a nonzero scattering cross section would show whether the a^-1/4 and a^-3/8 scalings and the p=12,20 behavior persist; if RHNs thermalize before anr, the temperature history collapses to the characteristic-momentum/prompt-thermalization result and the claimed departures disappear.","tokens_in":16569,"feed_emoji":"🌡️","tokens_out":4384,"duration_ms":35847,"temperature":0.7,"pith_summary":"Seesaw cosmology is the claim that reheating can proceed in two steps: the inflaton first decays into right-handed neutrinos (RHNs), and only then do those RHNs decay into Standard Model particles. Because the RHNs are produced relativistically and live long enough to dominate the energy budget, the expansion of the Universe oscillates between matter-like and radiation-like behavior—0, 1/3, 0, 1/3—instead of going directly to radiation. The key result is a distinctive temperature history: after an early plateau, T falls as a^-1/4 during relativistic-RHN domination and as a^-3/8 during nonrelativistic-RHN domination, before the standard a^-1. This matters because pre-BBN cosmology is almost unconstrained, and any process that depends on the temperature during reheating—dark-matter production, baryogenesis, gravitational-wave backgrounds—will shift accordingly. The paper shows concretely that direct dark-matter production from inflaton decays is enhanced by m_phi/(2m_N), while ultraviolet freeze-in acquires critical powers p=12 and p=20.","feed_headline":"Reheating through seesaw neutrinos alternates matter and radiation eras","feed_subtitle":"The same particles that give neutrinos mass dictate a new pre-BBN temperature history, boosting some dark-matter production.","key_machinery":"The central object is the production-time-dependent differential RHN number density δn_N(a,a') da' (Eq. 3.9), which tracks how many RHNs created between a' and a'+da' survive until a, including the redshift of their momentum p_N(a,a') = (a'/a)√(m_phi^2/4 - m_N^2) and the time-dilated decay rate Γ~_N = (m_N/E_N) Γ_N. Integrating this distribution over production times yields the total RHN number density, energy density, pressure, and SM radiation source; the relativistic-to-nonrelativistic transition emerges from the distribution rather than from any single-momentum approximation. This machinery produces the a^-1/4 and a^-3/8 temperature scalings and the p=12,20 freeze-in critical powers.","core_discovery":"The central discovery is that the full production-time dependence of the nonthermal RHN distribution cannot be replaced by a single average momentum: RHNs born at different times have different redshifted momenta, time-dilated decay rates, and equations of state. Retaining that distribution, the paper derives an integro-differential Boltzmann system and analytic solutions in four asymptotic regimes. The effective equation of state alternates as 0 -> 1/3 -> 0 -> 1/3, and the SM temperature follows T ∝ a^0, a^-1/4, a^-3/8, a^-1. In the hierarchical limit m_phi >> 2m_N, the final DM yield from inflaton decays is controlled by m_N and enhanced by m_phi/(2m_N); for UV freeze-in, the critical expo","pith_inferences":["The collisionless assumption is likely to break when RHN number densities are high enough for elastic scattering or inverse decays; a natural testable extension is to include a finite thermalization rate and check whether the four-stage sequence survives when RHNs scatter even weakly.","The alternating equation of state 0→1/3→0→1/3 should leave a characteristic imprint on primordial gravitational-wave spectra, since the spectral tilt changes with each era; a dedicated computation could distinguish seesaw cosmology from single-momentum approximations.","The enhancement m_phi/(2m_N) suggests that in minimal seesaw scenarios with normal or inverted mass ordering, the dark-matter abundance and the duration of the RHN eras could be predicted from neutrino oscillation data plus the inflaton mass, defining target parameter regions for experiments.","The p=12 and p=20 critical powers single out operators of mass dimension D=8 and D=12 as having special freeze-in behavior; laboratory probes of such operators, combined with the measured DM abundance, could constrain the ratios R1 and R2 and hence the seesaw parameters."],"forward_implications":["If correct, the pre-BBN thermal history is not generic radiation domination but a four-stage sequence whose duration is set by seesaw parameters, making reheating more predictive because RHN decay widths are bounded by neutrino mass-squared differences.","Dark-matter yield from a subdominant inflaton decay channel is set by m_N rather than m_phi and can be boosted by m_phi/(2m_N) relative to standard reheating at fixed Trh and branching fraction.","UV freeze-in production is not dominated by the final reheating temperature for p ≥ 12: for p=12 there is a logarithmic enhancement, for p=20 a R^3 ln(R2) boost, and for p>20 a double power-law enhancement, so the RHN-dominated eras can dominate the DM abundance.","The maximum temperature Tmax depends on both Γ_phi and Γ_N, while Trh depends only on Γ_N; both are seeded by neutrino parameters, so the seesaw scenario links observable neutrino masses to the duration of each cosmological era."],"fun_headline_variants":["Seesaw neutrinos set the cosmic clock before BBN","Neutrino masses drive a new reheating history","Seesaw cosmology: four eras from one particle species","Right-handed neutrinos orchestrate early-universe eras","How seesaw neutrinos reshape dark matter production"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The whole four-stage history assumes the right-handed neutrinos do not scatter, annihilate, or inverse-decay between production and decay: they are collisionless, so their nonthermal momentum distribution is only redshifted and time-dilated.","fun_headline_variants_meta":{"raw":{"variants":["Seesaw neutrinos set the cosmic clock before BBN","Neutrino masses drive a new reheating history","Seesaw cosmology: four eras from one particle species","Right-handed neutrinos orchestrate early-universe eras","How seesaw neutrinos reshape dark matter production"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000823,"raw_usage":{"total_tokens":3464,"prompt_tokens":799,"completion_tokens":2665,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":543,"completion_tokens_details":{"reasoning_tokens":2587}},"tokens_in":543,"tokens_out":2665,"duration_ms":17104,"temperature":1.0,"reasoning_tokens":2587,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T09:28:22.144748+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"If one includes even a small elastic-scattering or number-changing interaction among RHNs (or inverse decays that repopulate the inflaton), the distribution in Eq. (3.9) no longer holds. A Boltzmann simulation with a nonzero scattering cross section would show whether the a^-1/4 and a^-3/8 scalings and the p=12,20 behavior persist; if RHNs thermalize before anr, the temperature history collapses to the characteristic-momentum/prompt-thermalization result and the claimed departures disappear.","supporting_citations":[],"review_version":1}