REVIEW 2 major objections 5 minor 1 cited by
This paper argues that the reheating temperature is set by the seesaw sector, not by the inflaton, through a long-lived intermediate scalar.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-02 07:44 UTC pith:GKEFS46G
load-bearing objection A clean, plausible realization of mediator-dominated reheating in a seesaw context; the central T_RH ~ sqrt(Γ_S M_P) scaling holds up, but the advertised plateau and slope changes rest on approximations the paper itself defers. the 2 major comments →
Seesaw reheating
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that reheating does not have to be a direct inflaton-to-radiation process. In the proposed seesaw reheating chain, phi -> S S -> N_R N_R, the intermediate scalar S is the field responsible for spontaneous lepton-number breaking, and its decay into right-handed neutrinos reheats the Universe. The reheating temperature obeys T_RH proportional to sqrt(Gamma_S M_P), replacing the standard T_RH proportional to sqrt(Gamma_phi M_P), so T_RH probes the mediator rather than the inflaton. The paper further shows that the thermal history contains a distinctive plateau, with rho_R roughly constant while S is relativistic and continuously produced, followed by a gradual transition rh
What carries the argument
The central object is the long-lived intermediate scalar S, and the key identity is its effective decay rate Gamma_eff^S = Gamma_S / <gamma_S>, where <gamma_S> is the production-weighted average Lorentz factor of the continuously produced scalar population. During inflaton domination this average is approximated by <gamma_S> ~ (3/5) (m_phi / (2 m_S)). This single averaging expression encodes the relativistic time-dilation suppression that creates the low-temperature plateau and produces the sequence of slopes in rho_R: roughly constant, then a^{-1}, then a^{-3/2}.
Load-bearing premise
The argument assumes that the continuously produced scalar population can be summarized by a single production-weighted average Lorentz factor; if the true momentum distribution decays differently, the plateau duration, the slope changes, and the final reheating temperature would shift.
What would settle it
Perform a full phase-space Boltzmann calculation for the decay chain phi -> S S -> N_R N_R, tracking the scalar momentum distribution instead of using a single average Lorentz factor. The decisive check is whether this calculation reproduces the plateau rho_R roughly constant followed by the slopes rho_R ~ a^{-1} and rho_R ~ a^{-3/2}; the paper itself defers this computation.
If this is right
- The reheating temperature no longer constrains the inflaton sector; it measures the mediator's decay width, so the usual direct link between inflationary observables and T_RH is broken in this class of models.
- The extended low-temperature plateau means the maximum temperature during reheating can be much lower than in conventional instantaneous-decay reheating, with direct consequences for thermal dark-matter production and freeze-in scenarios.
- The relativistic-to-non-relativistic transition of the intermediate scalar adds a new timescale to the pre-radiation expansion history, which can imprint on inflationary observables and on the stochastic gravitational-wave background.
- Through the relation T_RH proportional to lambda_S^{1/4} sqrt(y_N M_P M_N), the reheating temperature becomes tied to the seesaw scale and to neutrino-mass parameters, so laboratory measurements of the neutrino sector could inform the early-Universe thermal history.
- The framework offers a well-motivated setting for sterile-neutrino dark matter and for leptogenesis, since the same scalar decays generate right-handed neutrinos.
- The paper states that this relies on prompt thermalization of the produced right-handed neutrinos; if that fails, the reheating history and the connection to neutrino masses would need revision.
Where Pith is reading between the lines
- Extension: The exact plateau shape is an artifact of using a single production-weighted Lorentz factor; a full momentum-space treatment would likely smear the plateau into a slow rise, although the qualitative separation of timescales should survive.
- Extension: The mechanism is not limited to a lepton-number-breaking scalar: any long-lived massive mediator that reheats through its decay would shift T_RH from the inflaton width to the mediator width, making this a general template for 'mediator reheating' scenarios.
- Extension: A genuinely low peak temperature suppresses the thermal production of heavy relics such as gravitinos and moduli, a potential benefit the paper leaves implicit and which could be tested by combining the plateau temperature with relic-abundance calculations.
- Extension: If the relation T_RH proportional to lambda_S^{1/4} sqrt(y_N M_P M_N) holds, then current and future bounds on neutrino masses and searches for lepton-number-breaking scalars could indirectly constrain the reheating temperature, and vice versa.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a new reheating scenario in which the inflaton first decays into a long-lived scalar S (identified with the B−L breaking field), which subsequently decays into right-handed neutrinos and reheats the Universe. Coupled Boltzmann equations are written for ρ_φ, ρ_S, and ρ_R, and piecewise analytical scalings are derived: during inflaton domination ρ_S ∝ a^{-3/2} and ρ_R remains roughly constant (a low-temperature plateau); after inflaton decay, relativistic S-domination gives ρ_R ∝ a^{-1}, and after S becomes non-relativistic ρ_R ∝ a^{-3/2}. The reheating temperature is argued to scale as T_RH ∝ sqrt(Γ_S M_P), replacing the standard T_RH ∝ sqrt(Γ_φ M_P), and, when S is the B−L scalar, T_RH is connected to neutrino parameters via T_RH ∝ λ_S^{1/4} sqrt(y_N M_P M_N). The paper is a proposal with analytical estimates and a benchmark numerical figure.
Significance. If the central scaling holds, the result is conceptually significant: the reheating temperature need not be a direct probe of the inflaton sector, but can instead be set by the lifetime of a seesaw mediator. The two advertised signatures — a low-temperature plateau and a double slope change in ρ_R from the relativistic-to-non-relativistic transition of S — are clear and, in principle, falsifiable through the thermal history. The paper builds on earlier work for the average energy of continuously produced scalars, cites it explicitly, and does not fit data. The main weakness is that the quantitative version of these signatures rests on a single-shell effective-Lorentz-factor approximation whose accuracy is not yet demonstrated; the authors themselves defer the full phase-space treatment. The asymptotic T_RH ∝ sqrt(Γ_S M_P) scaling is plausible and likely robust, but the quantitative reheating history is not yet fully established.
major comments (2)
- [Analytical understanding of Seesaw Reheating, Eqs. (4), (11), (14), (17), (20), (23)] The advertised thermal history is derived by replacing the S population with a single shell at the production-weighted mean energy. The actual decay term in Eq. (5) is Γ_S m_S n_S; writing it as Γ_eff^S ρ_S requires the energy-density-weighted mean Lorentz factor of the present population, which is not shown to equal the production-weighted average once low-γ particles have decayed away and once the source has turned off. The full phase-space treatment is explicitly deferred, yet Eqs. (14), (17), and (20) — the plateau and the double-slope signatures — and the numerical coefficient in Eq. (23) are controlled by this closure. The asymptotic scaling T_RH ∝ sqrt(Γ_S M_P) is plausibly robust, but the quantitative reheating history is not. Please provide the full phase-space solution for the benchmark, or at least a controlled moment-equation estimate of the error, and state how w_S and ⟨γ_S⟩
- [The Seesaw Reheating Mechanism, footnote 2 and Eqs. (22)–(26)] The identification of ρ_R with the thermal bath temperature assumes prompt thermalization of the right-handed neutrinos (Γ_th ≫ H). The paper defers a general treatment and states that for the benchmark thermalization becomes efficient, but no quantitative condition is given. Since T_RH in Eq. (22) is a temperature of the radiation bath, and the neutrino-mass connection in Eq. (26) uses this identification, the assumption is load-bearing. Please give the thermalization rate for the benchmark and state the parameter range in which prompt thermalization holds.
minor comments (5)
- [References] Several reference entries are incomplete, with missing years or volume numbers, e.g. Refs. [4], [10], [11], [17], [19], [20], [27]. Please standardize.
- [Fig. 1] The numerical integration used to produce Fig. 1 is not described. In particular, it is unclear whether w_S and ⟨γ_S⟩ are evolved via the single-shell approximation or via a full distribution. This should be stated in the caption or text.
- [Eq. (11)] The estimate ⟨E_S⟩ ≃ (3/5)(m_φ/2) should explicitly state the assumptions: matter-dominated expansion and neglect of S decays. This is important because the factor 3/5 follows from the production history and would change if the source or expansion history differs.
- [Eqs. (21)–(22)] The definition of reheating temperature via ρ_R(a_RH)=ρ_S(a_RH) differs from the common H=Γ definition. After ρ_R=ρ_S, S continues to decay and heat the bath, so the maximum temperature may occur later. The authors should clarify how their definition relates to the conventional one.
- [The Seesaw Reheating Mechanism] The time-dependent effective mass of S induced by μφS^2 is acknowledged but deferred. A brief statement of the parameter range in which this effect is negligible would help the reader assess the benchmark.
Circularity Check
No significant circularity: T_RH follows from a self-contained Boltzmann calculation; plateau and slope changes are scaling solutions, not fitted outputs. The one self-cited average-energy input is an independent published result.
full rationale
The derivation is not circular. Starting from the coupled Boltzmann equations (5) and the Friedmann equation (6), the paper solves the source-dominated regime: ρ_φ ∝ a^{-3}, H ∝ a^{-3/2}, and Eq. (8) gives ρ_S ∝ a^{-3/2}, so ρ_S/H is constant. With Γ_eff^S constant during continuous replenishment, Eq. (10) then gives ρ_R ≃ const, i.e. the plateau (14). The later slope changes ρ_R ∝ a^{-1} and ρ_R ∝ a^{-3/2} follow from the assumed scalings ⟨γ_S⟩ ∝ a^{-1} (relativistic redshift) and ⟨γ_S⟩ ≃ 1 (non-relativistic), respectively; each is a scaling consequence of the stated single-shell average, not an imposition of the desired answer. The central result T_RH ∝ sqrt(Γ_S M_P) (Eq. 23/28) is obtained by the standard reheating condition ρ_R(a_RH)=ρ_S(a_RH) with H~Γ_S; no data are fitted and no fitted parameter is relabeled as a prediction. The only external input is ⟨E_S⟩ ≈ (3/5)(m_φ/2) (Eq. 11), cited to Refs. [27,28]; Ref. [27] includes the author, but it is an independent published calculation of continuous non-thermal production, not a restatement of this paper's target. The paper explicitly defers a full phase-space treatment and finite thermalization to companion papers; these are unquantified approximations, but they are not circular reductions. There is no self-definition, no unique-theorem imported from the authors, and no renaming of a known result as a new prediction.
Axiom & Free-Parameter Ledger
free parameters (1)
- production-weighted average energy coefficient <E_S> ~ (3/5)(m_phi/2) =
3/5 m_phi/2
axioms (5)
- domain assumption Boltzmann equations (5) with Gamma_eff^S = Gamma_S/<gamma_S> capture the scalar decay dynamics
- domain assumption <E_S> ~ (3/5)(m_phi/2) for continuously produced scalars
- domain assumption Right-handed neutrinos thermalize promptly so their energy can be identified with rho_R
- domain assumption The S decay width uses its vacuum expression; oscillating-condensate-induced mass delay of phi -> SS is neglected
- domain assumption Standard seesaw relations M_N = y_N v_S/sqrt(2), m_S = sqrt(2 lambda_S) v_S, m_nu = y_nu^2 v^2/M_N
read the original abstract
We introduce the Seesaw Reheating scenario, in which the inflaton transfers its energy to a long-lived intermediate scalar associated with the spontaneous breaking of lepton number before the Universe is reheated through its decay into right-handed neutrinos. As a result, the reheating temperature is no longer determined by the inflaton decay width but by the dynamics of the seesaw sector. We derive analytical solutions describing the complete reheating history, revealing two characteristic features of this scenario: the relativistic time dilation of the intermediate scalar, which suppresses its decay and delays the transfer of energy to the thermal bath, and its subsequent transition from a relativistic to a non-relativistic regime, introducing a new characteristic timescale in the thermal history. Together, these effects lead to simple analytical expressions for the reheating temperature. When the intermediate scalar is identified with the field responsible for the spontaneous breaking of lepton number, the same framework naturally connects the reheating temperature to the origin of neutrino masses and provides a well-motivated setting for sterile-neutrino dark matter.
Figures
Forward citations
Cited by 1 Pith paper
-
Seesaw Cosmology
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.
Reference graph
Works this paper leans on
-
[1]
J. Martin and C. Ringeval, Phys. Rev. D82, 023511 (2010), arXiv:1004.5525 [astro-ph.CO]
Pith/arXiv arXiv 2010
-
[2]
L. Dai, M. Kamionkowski, and J. Wang, Phys. Rev. Lett. 113, 041302 (2014), arXiv:1404.6704 [astro-ph.CO]
Pith/arXiv arXiv 2014
-
[3]
M. A. G. Garcia, K. Kaneta, Y. Mambrini, and K. A. Olive, Phys. Rev. D101, 123507 (2020), arXiv:2004.08404 [hep-ph]
Pith/arXiv arXiv 2020
-
[4]
M. A. G. Garcia, K. Kaneta, Y. Mambrini, K. A. Olive, and S. Verner, JCAP03(03), 016, arXiv:2109.13280 [hep-ph]
-
[5]
S. E. Henrich, Y. Mambrini, and K. A. Olive, Phys. Rev. D111, 083501 (2025), arXiv:2412.13288 [hep-ph]
Pith/arXiv arXiv 2025
-
[6]
S. E. Henrich, M. Gross, Y. Mambrini, and K. A. Olive, Phys. Rev. D112, 103538 (2025), arXiv:2505.04703 [hep- ph]
Pith/arXiv arXiv 2025
-
[7]
S. E. Henrich, Y. Mambrini, and K. A. Olive, Phys. Rev. Lett.135, 221002 (2025), arXiv:2511.02117 [hep-ph]
arXiv 2025
-
[8]
X. Bertou, O. Deligny, M. Gross, Y. Mambrini, and I.- E. Mellouki, When direct detection constrains reheat- ing temperature: freeze-in with stronger couplings and inflaton-seeded freeze-in (2026), arXiv:2606.12408 [hep- ph]
Pith/arXiv arXiv 2026
-
[9]
M. Drewes and L. Ming, Phys. Rev. Lett.133, 031001 (2024), arXiv:2208.07609 [hep-ph]
Pith/arXiv arXiv 2024
-
[10]
J. L. Cook, E. Dimastrogiovanni, D. A. Easson, and L. M. Krauss, JCAP04, 047, arXiv:1502.04673 [astro-ph.CO]
-
[11]
C. Gross, O. Lebedev, and Y. Mambrini, JHEP08, 158, arXiv:1505.07480 [hep-ph]
-
[12]
Tenkanen, JHEP09, 049, arXiv:1607.01379 [hep-ph]
T. Tenkanen, JHEP09, 049, arXiv:1607.01379 [hep-ph]
-
[13]
T. Asaka, K. Hamaguchi, M. Kawasaki, and T. Yanagida, Phys. Rev. D61, 083512 (2000), arXiv:hep-ph/9907559
Pith/arXiv arXiv 2000
-
[14]
Lazarides and Q
G. Lazarides and Q. Shafi, Phys. Lett. B258, 305 (1991)
1991
-
[15]
W. Buchmuller, R. D. Peccei, and T. Yanagida, Ann. Rev. Nucl. Part. Sci.55, 311 (2005), arXiv:hep- ph/0502169
arXiv 2005
- [16]
-
[17]
C. Cosme, F. Costa, and O. Lebedev, JCAP06, 031, arXiv:2402.04743 [hep-ph]
- [18]
-
[19]
M. A. Amin and M. A. G. Garcia, JCAP08, 042, arXiv:1812.05559 [hep-ph]
- [20]
-
[21]
Y. Mambrini and K. A. Olive, Phys. Rev. D103, 115009 (2021), arXiv:2102.06214 [hep-ph]
Pith/arXiv arXiv 2021
-
[22]
N. Bernal, S. Cl´ ery, Y. Mambrini, and Y. Xu, JHEP06, 023, arXiv:2311.12694 [hep-ph]
-
[23]
S. Clery, Y. Mambrini, K. A. Olive, and S. Verner, Phys. Rev. D105, 075005 (2022), arXiv:2112.15214 [hep-ph]
Pith/arXiv arXiv 2022
-
[24]
R. T. Co, Y. Mambrini, and K. A. Olive, Phys. Rev. D 106, 075006 (2022), arXiv:2205.01689 [hep-ph]
Pith/arXiv arXiv 2022
-
[25]
M. A. G. Garcia, K. Kaneta, Y. Mambrini, and K. A. Olive, JCAP04, 012, arXiv:2012.10756 [hep-ph]
Pith/arXiv arXiv 2012
-
[26]
B. Barman, N. Bernal, and J. Rubio, Nucl. Phys. B1018, 116996 (2025), arXiv:2503.19980 [hep-ph]
Pith/arXiv arXiv 2025
-
[27]
M. A. G. Garcia, M. Gross, Y. Mambrini, K. A. Olive, M. Pierre, and J.-H. Yoon, JCAP12, 028, arXiv:2308.16231 [hep-ph]
-
[28]
M. A. G. Garcia, S. E. Henrich, W. Ke, and K. A. Olive, Leptogenesis and Low Reheating Temperatures (2026), arXiv:2607.08663 [hep-ph]
Pith/arXiv arXiv 2026
-
[29]
M. Gross, S. E. Henrich, and F. Koutroulis, Gravita- tional ultra-relativistic freeze-out during general reheat- ing (2026), arXiv:2606.25033 [hep-ph]
Pith/arXiv arXiv 2026
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.