Pith. sign in

REVIEW 3 major objections 5 minor 49 references

A derivative coupling between an axion-like inflaton and a heavy Majorana fermion can produce the right-handed neutrinos that source leptogenesis, even when the reheating temperature is low.

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 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-02 04:44 UTC pith:7LQLTLZI

load-bearing objection Solid phenomenological study of axion-inflaton-induced leptogenesis, but the quantitative baryon yield is softer than claimed because the delayed-decay production ignores decay on the saturation timescale and the intermediate lifetime window is left open. the 3 major comments →

arxiv 2607.13592 v1 pith:7LQLTLZI submitted 2026-07-15 hep-ph astro-ph.CO

Right-Handed Neutrino Production by an Axion-like Inflaton: Implications for Leptogenesis

Weiyi Deng , Chengcheng Han , Wuzhou Yin , Tong Ju This is my paper
classification hep-ph astro-ph.CO
keywords right-handed neutrinoaxion-like inflatonleptogenesispreheatingnon-thermal productionhelicity asymmetrybaryon asymmetryseesaw mechanism
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper argues that a derivative coupling between an axion-like inflaton and a heavy Majorana fermion can supply the right-handed neutrinos needed for leptogenesis, even when the reheating temperature is far below the neutrino mass scale. During inflation the rolling inflaton generates a helicity-asymmetric spectrum of these fermions; during preheating the oscillating inflaton produces repeated non-adiabatic bursts with a helicity-dependent momentum cutoff. The final abundance depends qualitatively on whether the produced fermions survive many oscillations (delayed decay, Pauli-blocked saturation) or decay within one oscillation (prompt decay, stage-by-stage sum). Applying the production history to non-thermal leptogenesis, the paper finds viable parameter regions in both regimes and shows the prompt-decay regime requires smaller CP violation because Pauli suppression does not accumulate.

Core claim

The central claim is that a dimension-five derivative interaction between an axion-like inflaton and a Majorana fermion—interpreted as the lightest right-handed neutrino—can, by itself, generate the observed baryon asymmetry through non-thermal leptogenesis. The rolling inflaton acts as a time-dependent axial background that shifts the effective momentum of each helicity mode by a term proportional to the inflaton velocity; this both imprints a helicity asymmetry during inflation and drives repeated non-adiabatic production whenever the shifted momentum crosses zero during preheating. The paper computes the resulting abundance in the two opposite lifetime limits and shows that both can match

What carries the argument

The load-bearing object is the helicity-dependent effective momentum k~r = r k/a − θdot (with θ = φ/f), whose zero crossings mark the non-adiabatic bursts in preheating and whose size sets the production cutoff. The calculation is carried out in the chirally rotated basis where the derivative coupling becomes a time-dependent complex mass, diagonalizing the instantaneous Hamiltonian to define the occupation number; for the prompt-decay regime, production is treated event-by-event by resetting the mode-mixing coefficients to vacuum at the start of each stage and truncating each stage at a cutoff proportional to the extremum of a|θdot|.

Load-bearing premise

The calculation hinges on the assumption that, in the prompt-decay regime, each stage starts from a truly empty fermion state because the produced neutrino decays within one inflaton oscillation and its decay products do not repopulate the modes or push back on the inflaton.

What would settle it

Evolve the same mode equations with a finite neutrino decay width between bursts instead of resetting the coefficients to vacuum for Γ_N1 between m_phi/100 and m_phi; if the resulting baryon asymmetry does not interpolate between the delayed and prompt results, the binary regime classification fails.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Non-thermal leptogenesis works even when the reheating temperature is orders of magnitude below the seesaw scale, as long as the inflaton–neutrino coupling is active during inflation and preheating.
  • In both delayed-decay regimes the final baryon yield is controlled by the same saturated preheating abundance, so the required CP asymmetry is identical whether the neutrino decays before or after reheating; only the allowed neutrino-mass/effective-mass window differs.
  • In the prompt-decay regime the same momentum modes can be repopulated in successive bursts, so the total abundance can exceed the Pauli-blocked value and the required CP violation is smaller.
  • The produced spectrum has a hard, helicity-dependent momentum cutoff that scales roughly linearly with the coupling M_Pl/f, giving a concrete prediction for the momentum range of the non-thermal neutrinos.
  • The physically correct production rate vanishes as the fermion mass goes to zero, a property that is obscured if the occupation number is defined in the original derivative-coupling basis.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the intermediate-lifetime region (m_phi/100 ≲ Γ_N1 ≲ m_phi) were solved with a coupled production–decay evolution, the baryon yield would likely interpolate between the delayed and prompt estimates, which could either enlarge or shrink the viable parameter space shown here.
  • The helicity asymmetry of the produced right-handed neutrinos is a new handle: since only one helicity is efficiently produced, the subsequent decays into left-handed leptons and Higgs may carry a chiral bias that standard leptogenesis calculations average over; exploring this could modify the effective CP-asymmetry efficiency.
  • The same derivative-coupling mechanism could apply to other fermions in the early Universe, such as dark-matter candidates or additional sterile states; because production vanishes for massless fermions, this framework predicts negligible production of exactly massless states and a strong relation between fermion mass and coupling.
  • A direct test would be to compute the gravitational-wave or non-Gaussianity signatures of the helicity-asymmetric fermion bursts; the paper does not do this, but the time-dependent effective momentum leaves a characteristic scale that might be probed by future observations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper develops a unified treatment of Majorana-fermion production from a derivative coupling to an axion-like inflaton, covering both the slow-roll inflationary stage and the subsequent preheating stage. The authors solve the mode equations in a Hamiltonian-diagonal ψ basis, obtain a helicity-asymmetric spectrum during inflation, and identify repeated non-adiabatic bursts during preheating with a helicity-dependent momentum cutoff. They distinguish a delayed-decay regime, where Pauli blocking leads to a saturated non-thermal abundance, from a prompt-decay regime, where each burst is treated as an independent source. These results are then applied to non-thermal leptogenesis by identifying the produced fermion with the lightest right-handed neutrino N1. The paper derives baryon-yield formulas in three lifetime regimes and claims that both the delayed and prompt regimes can reproduce the observed baryon asymmetry in viable regions of parameter space, with the prompt regime requiring a smaller CP asymmetry.

Significance. If the quantitative claims hold, the paper provides a useful connection between axion-like inflation and non-thermal leptogenesis, showing that the observed baryon asymmetry can be reproduced even when the reheating temperature is far below the seesaw scale. The analysis has several genuine strengths: the analytic slow-roll solution is a nontrivial derivation; the discussion of why the Y-basis occupation number is unphysical in the small-mass limit and why the ψ-basis should be used is clear and important; the preheating calculation includes Pauli blocking in the delayed regime and a stage-by-stage summation in the prompt regime; and the leptogenesis section makes use of standard, named constraints such as the Davidson-Ibarra bound and ΔL=2 washout. The paper is also transparent about several limitations, notably the unanalyzed intermediate-lifetime regime. The central claims are plausible but, for the reasons detailed below, the quantitative parameter-space boundaries and the delayed/prompt dichotomy need further work before the results can be accepted at face value.

major comments (3)
  1. [Sec. 3.2.1, Sec. 4.1, Eqs. (4.41)–(4.55)] The delayed-decay abundance is computed by solving the Bogolyubov equations (3.39) with no decay term and evaluating the result at t_sat = 300 m_φ^{-1} (Eq. 4.43). The same saturated comoving abundance is then used for the entire delayed region Γ_N1 ≲ m_φ/100. However, for Γ_N1 = m_φ/100 the lifetime is 100 m_φ^{-1}, so a population produced near the start of the build-up has decayed by a factor e^{-3} ≈ 5% by t_sat. Decay during production is therefore not negligible at the upper boundary of the delayed region. The manuscript does not establish whether the correct quantity for the baryon yield is the stable-particle abundance, the number of N1 present at t_sat, or the cumulative number of decays including those before t_sat; these differ at the O(1) level or more when Γ_N1 is near m_φ/100. Consequently, the baryon-yield relations (4.55) and (4.68) and the required-ϵ curves in Figs. 9–10
  2. [Sec. 4, Eqs. (4.13), (4.36), (4.78), Fig. 8] The lifetime inequalities used to separate the delayed and prompt regimes compare the rest-frame width Γ_N1 with m_φ. But the N1 population produced during preheating is relativistic: the physical momenta extend up to k_phys ∼ (MPl/f) m_φ / a(t), so for the benchmarks M1 = 10 m_φ and MPl/f = O(10^3), the typical Lorentz factor is γ ∼ O(10–100). In the cosmic frame the decay rate is Γ_N1/γ, not Γ_N1. The prompt-decay condition (4.78) should therefore be Γ_N1/γ ≳ m_φ (up to O(1) factors from the inter-burst spacing), and the boundary Γ_N1 = m_φ in Fig. 8 is not the physical boundary for relativistic particles. This affects the classification of parameter points and the meaning of the prompt-regime abundance (4.89). The authors should either use a momentum-dependent decay rate, integrate over the produced spectrum, or restrict the benchmarks to parameter points where γ = O(1).
  3. [Sec. 3.2.2, Eqs. (3.55)–(3.59), Sec. 4.3] The prompt-decay calculation is an event-based prescription: the Bogolyubov coefficients are reset to vacuum at the start of each stage (Eq. 3.55) and each stage is truncated at k_{r,max,n} = C a(t_ext)|θdot(t_ext)| with C = 1.1 (Eq. 3.59). This assumes complete decay between bursts, no repopulation of the fermion modes by decay products, and no backreaction on the inflaton. These are reasonable limiting assumptions, but the manuscript explicitly leaves the intermediate regime m_φ/100 ≲ Γ_N1 ≲ m_φ unanalyzed. Since the delayed and prompt calculations use different definitions of the relevant abundance, and since the final baryon yield is a continuous function of Γ_N1, the binary delayed/prompt approximation is not sufficient to establish the claimed viable regions across the whole parameter space shown in Fig. 8. A quantitative version of the central claim needs either a continuous produ
minor comments (5)
  1. [Sec. 3.1] The text states that the analytic inflationary spectrum agrees with Ref. [42], but no comparison plot or explicit mapping is shown. Since this is an external validation of a central analytic result, a figure or a brief numerical comparison should be included.
  2. [Sec. 4, Eqs. (4.33)–(4.34)] The numerical EFT cutoff estimate is internally inconsistent. With a(t_max) ≈ 1.5 and k_{r,max} ≈ 0.9 (MPl/f) m_φ, the inequality m_φ MPl/(a f) ≲ 2π f gives f ≳ 3.4 × 10^15 GeV and MPl/f ≲ 3.6 × 10^3, not the quoted f ≳ 1.6 × 10^15 GeV and MPl/f < 1500. The plotted range MPl/f ≤ 1400 is safe, but the quoted numbers should be corrected.
  3. [Sec. 4.1, Sec. 4.3] The Davidson–Ibarra bound is used in Secs. 4.1 and 4.2 but is only introduced at the end of Sec. 4.3. It would be clearer to define it before the first use.
  4. [Eq. (2.2) and throughout] There is a typo in Eq. (2.2): 'FLR W' should be 'FLRW'. Also, the notation 'Ybasis' and 'ψbasis' is occasionally run together; consistent typesetting would improve readability.
  5. [Figs. 9–11] The required-ϵ curves are presented as functions of MPl/f for fixed T_RH and M1. It would be helpful to indicate on the plots (or in the captions) the corresponding saturated comoving abundance or (n_{N1}/ρ_φ)_pre, since this is the quantity that controls the normalization through Eq. (4.55).

Circularity Check

0 steps flagged

No significant circularity: the production spectra are obtained by solving mode equations and the leptogenesis analysis is a parameter-space scan against an external CP bound.

full rationale

The central derivation is self-contained. Fermion occupation numbers are computed by evolving the Bogolyubov-mode equations (3.38)-(3.39) derived from the Lagrangian (2.2), with the instantaneous number defined by diagonalizing the ψ-basis Hamiltonian (2.31)-(2.36); the final abundances are integrals over these numerically computed distributions (3.49), not fits to the baryon asymmetry. The inflationary spectrum is independently cross-checked against Ref. [42], and the preheating momentum cutoffs are tied to zero crossings of the effective momentum k̃_r (3.46)-(3.48). The prompt-decay coefficient C=1.1 in Eq. (3.59) is explicitly described as a practical numerical choice with only mild sensitivity, so it is not a load-bearing fitted parameter. In the leptogenesis section, the paper does not claim to predict n_B/s; it solves for the CP asymmetry ϵ required to match the observed n_B/s (e.g., Eq. (4.55)) and then tests whether that ϵ lies below the Davidson-Ibarra bound and other constraints. That is standard parameter-space mapping rather than circularity. The self-citations [33,34] appear only in an introductory list of prior perturbative inflaton-decay mechanisms and play no role in the derivation. The explicitly unanalyzed intermediate lifetime region (m_ϕ/100 ≲ Γ_N1 ≲ m_ϕ) and the reset-to-vacuum approximation (3.55) are acknowledged limitations and potential correctness risks, but they are approximations, not reductions of the output to the input. No load-bearing step is equivalent by construction to its own input, so the circularity score is 0.

Axiom & Free-Parameter Ledger

5 free parameters · 7 axioms · 0 invented entities

The central production calculation relies on standard QFT-in-curved-spacetime tools plus a set of benchmark/modeling assumptions; the leptogenesis estimate adds standard thermal formulas and a free CP asymmetry. No new particles, forces, or conserved quantities are introduced.

free parameters (5)
  • Inflaton mass m_phi = 10^13 GeV
    Benchmark choice for the quadratic preheating potential (Eq. 3.35); the leptogenesis contours are quoted for this value.
  • Preheating initial conditions phi(0), phidot(0) = M_Pl, -0.7 M_Pl m_phi
    Chosen initial conditions for the numerical background (Eq. 3.36); not derived from the inflationary phase.
  • Slow-roll velocity |phidot|/H_inf^2 = O(10^3)
    Benchmark stated in Sec. 3.1 for the analytic inflationary spectrum.
  • Stage cutoff fudge factor C = 1.1
    Introduced in Eq. (3.59) to include the spectral tail; the authors state it is a practical numerical choice.
  • CP asymmetry epsilon = Required values from matching nB/s (Figs. 9-11)
    Not computed from the model's Yukawa sector; solved from matching the observed baryon asymmetry and bounded by the Davidson-Ibarra bound. This is the parameter that makes the scenario 'consistent'.
axioms (7)
  • domain assumption Inflaton is a homogeneous classical background; metric perturbations and fermion backreaction are neglected
    Stated in Sec. 2.1 before Eq. (2.2).
  • domain assumption Bunch-Davies vacuum initial state; slow-roll with approximately constant phidot
    Sec. 3.1, Eqs. (3.3)-(3.4).
  • domain assumption Quadratic inflaton potential during preheating and Friedmann background evolution
    Sec. 3.2, Eqs. (3.32)-(3.34); this is not the full inflationary potential.
  • domain assumption Instantaneous reheating; no washout before reheating; no thermal regeneration for M1 > TRH
    Sec. 4, Eqs. (4.16)-(4.18) and the text before Eq. (4.14); if false, the final asymmetry could be reduced.
  • ad hoc to paper Prompt-decay stage independence: vacuum reset between bursts and no repopulation/backreaction
    Sec. 3.2.2, Eq. (3.55); this is the key modeling assumption distinguishing prompt from delayed regimes.
  • domain assumption Effective field theory cutoff Lambda = 2πf and derivative coupling as the leading operator
    Sec. 4, Eq. (4.33).
  • standard math Standard Bogolyubov formalism for fermions in time-dependent backgrounds
    Used throughout Secs. 2-3; not proved in the paper.

pith-pipeline@v1.3.0-alltime-deepseek · 33655 in / 17867 out tokens · 175531 ms · 2026-08-02T04:44:17.743532+00:00 · methodology

0 comments
read the original abstract

We study heavy right-handed neutrino production induced by a derivative coupling to an axion-like inflaton and its implications for non-thermal leptogenesis. We develop a unified treatment of Majorana-fermion production during inflation and preheating. During inflation, the rolling inflaton background generates a helicity-asymmetric spectrum that can be obtained analytically in the slow-roll regime, while during preheating the oscillating background drives repeated non-adiabatic production events with a helicity-dependent effective momentum and a natural momentum cutoff. We compute the resulting abundance in both the delayed-decay and prompt-decay limits, including the effects of Pauli blocking and successive production stages. We further clarify the relation between the fermion basis used to identify non-adiabatic production and the Hamiltonian-diagonal basis used to define instantaneous occupation numbers. Applying these results to non-thermal leptogenesis, we identify parameter regions consistent with the observed baryon asymmetry.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

49 extracted references · 34 linked inside Pith

  1. [5]

    Minkowski,µ→eγat a Rate of One Out of10 9 Muon Decays?,Phys

    P. Minkowski,µ→eγat a Rate of One Out of10 9 Muon Decays?,Phys. Lett. B67(1977) 421

    1977

  2. [6]

    Yanagida,Horizontal gauge symmetry and masses of neutrinos,Conf

    T. Yanagida,Horizontal gauge symmetry and masses of neutrinos,Conf. Proc. C7902131 (1979) 95

    1979

  3. [7]

    Yanagida,Horizontal Symmetry and Masses of Neutrinos,Prog

    T. Yanagida,Horizontal Symmetry and Masses of Neutrinos,Prog. Theor. Phys.64(1980) 1103

    1980

  4. [8]

    Gell-Mann, P

    M. Gell-Mann, P. Ramond and R. Slansky,Complex Spinors and Unified Theories,Conf. Proc. C790927(1979) 315 [1306.4669]

    Pith/arXiv arXiv 1979

  5. [9]

    Mohapatra and G

    R.N. Mohapatra and G. Senjanovic,Neutrino Mass and Spontaneous Parity Nonconservation,Phys. Rev. Lett.44(1980) 912. – 48 –

    1980

  6. [10]

    Mohapatra and G

    R.N. Mohapatra and G. Senjanovic,Neutrino Masses and Mixings in Gauge Models with Spontaneous Parity Violation,Phys. Rev. D23(1981) 165

    1981

  7. [11]

    Schechter and J.W.F

    J. Schechter and J.W.F. Valle,Neutrino Masses in SU(2) x U(1) Theories,Phys. Rev. D22 (1980) 2227

    1980

  8. [12]

    Schechter and J.W.F

    J. Schechter and J.W.F. Valle,Neutrino Decay and Spontaneous Violation of Lepton Number,Phys. Rev. D25(1982) 774

    1982

  9. [13]

    Sakharov,Violation of CP Invariance, C asymmetry, and baryon asymmetry of the universe,Pisma Zh

    A.D. Sakharov,Violation of CP Invariance, C asymmetry, and baryon asymmetry of the universe,Pisma Zh. Eksp. Teor. Fiz.5(1967) 32

    1967

  10. [14]

    Fukugita and T

    M. Fukugita and T. Yanagida,Baryogenesis Without Grand Unification,Phys. Lett. B174 (1986) 45

    1986

  11. [15]

    Kuzmin, V.A

    V.A. Kuzmin, V.A. Rubakov and M.E. Shaposhnikov,On the Anomalous Electroweak Baryon Number Nonconservation in the Early Universe,Phys. Lett. B155(1985) 36

    1985

  12. [16]

    Harvey and M.S

    J.A. Harvey and M.S. Turner,Cosmological Baryon and Lepton Number in the Presence of Electroweak Fermion Number Violation,Phys. Rev. D42(1990) 3344

    1990

  13. [17]

    Luty,Baryogenesis via leptogenesis,Phys

    M.A. Luty,Baryogenesis via leptogenesis,Phys. Rev. D45(1992) 455

    1992

  14. [18]

    Plumacher,Baryogenesis and lepton number violation,Z

    M. Plumacher,Baryogenesis and lepton number violation,Z. Phys. C74(1997) 549 [hep-ph/9604229]

    Pith/arXiv arXiv 1997

  15. [19]

    Barbieri, P

    R. Barbieri, P. Creminelli, A. Strumia and N. Tetradis,Baryogenesis through leptogenesis, Nucl. Phys. B575(2000) 61 [hep-ph/9911315]

    Pith/arXiv arXiv 2000

  16. [20]

    Giudice, A

    G.F. Giudice, A. Notari, M. Raidal, A. Riotto and A. Strumia,Towards a complete theory of thermal leptogenesis in the SM and MSSM,Nucl. Phys. B685(2004) 89 [hep-ph/0310123]

    Pith/arXiv arXiv 2004

  17. [21]

    Buchmuller, P

    W. Buchmuller, P. Di Bari and M. Plumacher,Leptogenesis for pedestrians,Annals Phys. 315(2005) 305 [hep-ph/0401240]

    Pith/arXiv arXiv 2005

  18. [22]

    Buchmuller, R.D

    W. Buchmuller, R.D. Peccei and T. Yanagida,Leptogenesis as the origin of matter,Ann. Rev. Nucl. Part. Sci.55(2005) 311 [hep-ph/0502169]

    Pith/arXiv arXiv 2005

  19. [23]

    Garbrecht,Why is there more matter than antimatter? Calculational methods for leptogenesis and electroweak baryogenesis,Prog

    B. Garbrecht,Why is there more matter than antimatter? Calculational methods for leptogenesis and electroweak baryogenesis,Prog. Part. Nucl. Phys.110(2020) 103727 [1812.02651]

    Pith/arXiv arXiv 2020

  20. [24]

    Lazarides and Q

    G. Lazarides and Q. Shafi,Origin of matter in the inflationary cosmology,Phys. Lett. B258 (1991) 305

    1991

  21. [25]

    Asaka, K

    T. Asaka, K. Hamaguchi, M. Kawasaki and T. Yanagida,Leptogenesis in inflationary universe,Phys. Rev. D61(2000) 083512 [hep-ph/9907559]

    Pith/arXiv arXiv 2000

  22. [26]

    Asaka, K

    T. Asaka, K. Hamaguchi, M. Kawasaki and T. Yanagida,Leptogenesis in inflaton decay, Phys. Lett. B464(1999) 12 [hep-ph/9906366]

    Pith/arXiv arXiv 1999

  23. [27]

    Asaka, H.B

    T. Asaka, H.B. Nielsen and Y. Takanishi,Nonthermal leptogenesis from the heavier Majorana neutrinos,Nucl. Phys. B647(2002) 252 [hep-ph/0207023]

    Pith/arXiv arXiv 2002

  24. [28]

    Senoguz and Q

    V.N. Senoguz and Q. Shafi,GUT scale inflation, nonthermal leptogenesis, and atmospheric neutrino oscillations,Phys. Lett. B582(2004) 6 [hep-ph/0309134]

    Pith/arXiv arXiv 2004

  25. [29]

    Hahn-Woernle and M

    F. Hahn-Woernle and M. Plumacher,Effects of reheating on leptogenesis,Nucl. Phys. B806 (2009) 68 [0801.3972]. – 49 –

    Pith/arXiv arXiv 2009

  26. [30]

    Buchmuller, V

    W. Buchmuller, V. Domcke, K. Kamada and K. Schmitz,A Minimal Supersymmetric Model of Particle Physics and the Early Universe,1309.7788

    Pith/arXiv arXiv

  27. [31]

    Barman, D

    B. Barman, D. Borah and R. Roshan,Nonthermal leptogenesis and UV freeze-in of dark matter: Impact of inflationary reheating,Phys. Rev. D104(2021) 035022 [2103.01675]

    Pith/arXiv arXiv 2021

  28. [32]

    Barman, A

    B. Barman, A. Basu, D. Borah, A. Chakraborty and R. Roshan,Testing leptogenesis and dark matter production during reheating with primordial gravitational waves,Phys. Rev. D 111(2025) 055016 [2410.19048]

    Pith/arXiv arXiv 2025

  29. [33]

    Han, H.-J

    C. Han, H.-J. He, L. Song and J. You,Cosmological signatures of neutrino seesaw mechanism,Phys. Rev. D112(2025) L081309 [2412.21045]

    arXiv 2025

  30. [34]

    J. You, L. Song, H.-J. He and C. Han,Cosmological non-Gaussianity from the neutrino seesaw mechanism,Phys. Rev. D112(2025) 083555 [2412.16033]

    arXiv 2025

  31. [35]

    Baacke, K

    J. Baacke, K. Heitmann and C. Patzold,Nonequilibrium dynamics of fermions in a spatially homogeneous scalar background field,Phys. Rev. D58(1998) 125013 [hep-ph/9806205]

    Pith/arXiv arXiv 1998

  32. [36]

    Felder, L

    G.N. Felder, L. Kofman and A.D. Linde,Instant preheating,Phys. Rev. D59(1999) 123523 [hep-ph/9812289]

    Pith/arXiv arXiv 1999

  33. [37]

    Greene and L

    P.B. Greene and L. Kofman,Preheating of fermions,Phys. Lett. B448(1999) 6 [hep-ph/9807339]

    Pith/arXiv arXiv 1999

  34. [38]

    Giudice, M

    G.F. Giudice, M. Peloso, A. Riotto and I. Tkachev,Production of massive fermions at preheating and leptogenesis,JHEP08(1999) 014 [hep-ph/9905242]

    Pith/arXiv arXiv 1999

  35. [39]

    Peloso and L

    M. Peloso and L. Sorbo,Preheating of massive fermions after inflation: Analytical results, JHEP05(2000) 016 [hep-ph/0003045]

    Pith/arXiv arXiv 2000

  36. [40]

    Adshead and E.I

    P. Adshead and E.I. Sfakianakis,Leptogenesis from left-handed neutrino production during axion inflation,Phys. Rev. Lett.116(2016) 091301 [1508.00881]

    Pith/arXiv arXiv 2016

  37. [41]

    Adshead and E.I

    P. Adshead and E.I. Sfakianakis,Fermion production during and after axion inflation,JCAP 11(2015) 021 [1508.00891]

    Pith/arXiv arXiv 2015

  38. [42]

    Adshead, L

    P. Adshead, L. Pearce, M. Peloso, M.A. Roberts and L. Sorbo,Phenomenology of fermion production during axion inflation,JCAP06(2018) 020 [1803.04501]

    Pith/arXiv arXiv 2018

  39. [43]

    Kulkarni and L

    A. Kulkarni and L. Sorbo,Spectra of fermions produced by a time-dependent axion in the radiation- and matter-dominated Universe,JCAP01(2025) 039 [2408.04550]

    Pith/arXiv arXiv 2025

  40. [44]

    Kaloper and L

    N. Kaloper and L. Sorbo,Where in the String Landscape is Quintessence,Phys. Rev. D79 (2009) 043528 [0810.5346]

    Pith/arXiv arXiv 2009

  41. [45]

    Kofman, A.D

    L. Kofman, A.D. Linde and A.A. Starobinsky,Reheating after inflation,Phys. Rev. Lett.73 (1994) 3195 [hep-th/9405187]

    Pith/arXiv arXiv 1994

  42. [46]

    Shtanov, J.H

    Y. Shtanov, J.H. Traschen and R.H. Brandenberger,Universe reheating after inflation,Phys. Rev. D51(1995) 5438 [hep-ph/9407247]

    Pith/arXiv arXiv 1995

  43. [47]

    Kitajima, T

    N. Kitajima, T. Sekiguchi and F. Takahashi,Cosmological abundance of the QCD axion coupled to hidden photons,Phys. Lett. B781(2018) 684 [1711.06590]

    Pith/arXiv arXiv 2018

  44. [48]

    Wang,Inflation, Cosmic Perturbations and Non-Gaussianities,Commun

    Y. Wang,Inflation, Cosmic Perturbations and Non-Gaussianities,Commun. Theor. Phys. 62(2014) 109 [1303.1523]

    Pith/arXiv arXiv 2014

  45. [49]

    Pilaftsis,Heavy Majorana neutrinos and baryogenesis,Int

    A. Pilaftsis,Heavy Majorana neutrinos and baryogenesis,Int. J. Mod. Phys. A14(1999) 1811 [hep-ph/9812256]. – 50 –

    Pith/arXiv arXiv 1999

  46. [50]

    Buchmuller and M

    W. Buchmuller and M. Plumacher,CP asymmetry in Majorana neutrino decays,Phys. Lett. B431(1998) 354 [hep-ph/9710460]

    Pith/arXiv arXiv 1998

  47. [51]

    Garbrecht and E

    B. Garbrecht and E. Wang,The neutrino mass bound from leptogenesis revisited,JHEP03 (2025) 008 [2411.09765]

    Pith/arXiv arXiv 2025

  48. [52]

    Davidson and A

    S. Davidson and A. Ibarra,A Lower bound on the right-handed neutrino mass from leptogenesis,Phys. Lett. B535(2002) 25 [hep-ph/0202239]

    Pith/arXiv arXiv 2002

  49. [53]

    Buchmuller, P

    W. Buchmuller, P. Di Bari and M. Plumacher,The Neutrino mass window for baryogenesis, Nucl. Phys. B665(2003) 445 [hep-ph/0302092]. – 51 –

    Pith/arXiv arXiv 2003