Pith. sign in

REVIEW 3 major objections 4 minor 5 cited by

The paper argues that in a quartic T-model inflaton potential, reheating through inflaton decays into fermions alone is effectively ruled out: for Yukawa couplings in the perturbative range, parametric resonance, kinematic blocking, and Pau

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-03 15:36 UTC pith:N43YTBCQ

load-bearing objection A careful non-perturbative study of fermion reheating in a quartic T-model, but the post-fragmentation impossibility claim rests on an uncomputed extrapolation near the BBN threshold. the 3 major comments →

arxiv 2512.16203 v2 pith:N43YTBCQ submitted 2025-12-18 hep-ph

Fermion (non)reheating with a quartic inflaton potential

Nabeen Bhusal , Ernesto Ch\'avez M. , Marcos A. G. Garcia , Adriana G. Menkara , Mathias Pierre This is my paper
classification hep-ph
keywords fermionic reheatinginflaton fragmentationquartic T-modelpreheatingPauli blockingparametric resonanceYukawa couplingBBN bound
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.

The paper tries to establish that reheating the universe after inflation via fermionic decays does not work in a quartic T-model inflaton potential. It shows that for Yukawa couplings y ≳ 10⁻⁷, which are intuitively perturbative, non-perturbative effects – parametric resonance, kinematic blocking, and Pauli suppression – unavoidably suppress fermion production, making the standard Boltzmann approximation fail. To equalize the inflaton and fermion energy densities before the inflaton condensate fragments requires enormous couplings y ≳ 0.2, for which backreaction and radiative corrections are not under control. After fragmentation, even accounting for full Fermi–Dirac statistics, fermionic decays of inflaton quanta cannot raise the reheating temperature above the big bang nucleosynthesis bound for couplings up to about 3×10⁻⁷. The paper concludes that successful reheating in this setup requires coupling the inflaton to a bosonic (integer-spin or spin-0) particle, not just to fermions.

Core claim

The central claim is that non-perturbative effects are inevitable for fermionic decay products in a quartic (λφ⁴) inflaton potential. Using the Heisenberg/Bogoliubov solution of the Dirac equation in the oscillating condensate, the authors find that for y, y₅ ≳ 10⁻⁷ the fermion phase-space distribution saturates the Pauli bound fψ = 1 over a wide range of momenta, and the energy density of produced fermions is suppressed by orders of magnitude relative to the perturbative Boltzmann prediction. Equalizing ρψ and ρφ prior to fragmentation requires y ≳ 0.2, a regime where backreaction on the inflaton dynamics is significant and was neglected. Treating post-fragmentation decays with the full Bol

What carries the argument

The central object is the fermion mode function in a time-dependent inflaton background, governed by the Dirac equation with effective mass matrix mψ,eff = mψ + yφ + i y₅φγ₅. Occupation numbers n_p are computed via Bogoliubov transformations, giving the fermion energy density ρψ = (2/a⁴) Σₛ ∫ d³p/(2π)³ ω_p n_p. For the inflaton sector, the Floquet analysis of Hill's equation for inflaton fluctuations (with Jacobi elliptic sn time dependence) determines the onset of fragmentation, and lattice data provide the inflaton phase-space distribution after fragmentation. The post-fragmentation analysis uses the full Boltzmann equation with Fermi–Dirac blocking factors, whose collision term includes P

Load-bearing premise

The impossibility conclusion rests on assuming that the produced fermions do not thermalize with the Standard Model on a timescale shorter than inflaton fragmentation; the paper states that this timescale is unknown, and if thermalization were fast, Pauli blocking would be lifted and fermionic reheating could proceed at lower couplings.

What would settle it

Perform a lattice simulation that includes the inflaton and the fermion field with full backreaction (energy-momentum feedback on the inflaton) for a coupling such as y = 0.1; if the fermion energy density approaches the inflaton energy density before fragmentation, the claimed y ≳ 0.2 threshold for pre-fragmentation reheating is falsified. Alternatively, measure the fermion energy density in a controlled 1+1D analog simulation where thermalization is artificially accelerated; if T_reh rises above the BBN bound, the post-fragmentation no-reheating conclusion would be evaded.

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

If this is right

  • The perturbative Boltzmann estimate of the reheating temperature is not simply a small correction; for y ≳ 10⁻⁷ it overestimates fermionic energy density by orders of magnitude, so naive perturbative reheating calculations for quartic potentials are unreliable.
  • To reheat a quartic T-model inflaton to temperatures above BBN, a bosonic decay channel (scalar or vector) is required; fermion-only models essentially close the window for successful reheating.
  • For pre-fragmentation reheating to succeed, Yukawa couplings y ≳ 0.2 are needed, which is far outside the perturbative regime and demands a full backreaction treatment and control of radiative corrections.
  • In the post-fragmentation regime, the Pauli-suppressed Boltzmann treatment shows that even couplings as large as y = 3×10⁻⁷ yield T_reh ≲ 4 MeV, so the current BBN bound is only marginally avoided at the very edge of the explored parameter space.
  • If fermionic reheating is excluded, then alternative inflaton decay channels or a thermalization mechanism for the fermions with the Standard Model must be invoked.

Where Pith is reading between the lines

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

  • If the produced fermions thermalize with Standard Model degrees of freedom on a timescale shorter than the onset of fragmentation, Pauli blocking would be lifted, potentially allowing successful reheating at much smaller Yukawa couplings; the paper explicitly leaves this thermalization timescale unknown, so this is a testable loophole rather than a contradiction.
  • The constant kinematic ratio R for quartic potentials suggests the impossibility result may extend beyond T-models to any λφ⁴ minimum, since the non-adiabaticity parameter does not decrease in time as it does for quadratic minima.
  • A direct numerical simulation that couples the Dirac field to the inflaton including backreaction (e.g., a lattice code with fermionic degrees of freedom) for y around 0.1 could decide whether the pre-fragmentation threshold really is as high as 0.2 or whether backreaction changes the picture.
  • The suppression mechanism here is a concrete illustration of a general principle: Pauli blocking sets a finite capacity for energy transfer from a coherent scalar condensate to fermions, so fermionic reheating is generically less efficient than naive decay-rate estimates.

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 / 4 minor

Summary. This paper studies reheating of a quartic T-model inflaton into a Dirac fermion via Yukawa and pseudoscalar couplings, combining nonperturbative Dirac-equation integration before fragmentation with a full Boltzmann treatment after fragmentation. The authors find that for y ≳ 10^-7–10^-8, resonance and Pauli effects invalidate the perturbative approximation; pre-fragmentation energy equalization requires very large couplings, y ≳ 0.2; and post-fragmentation, the computed fermion energy density is suppressed by more than an order of magnitude relative to the no-statistics approximation. They conclude that fermionic reheating in this model generally cannot reach BBN temperatures, so a bosonic decay channel or fast thermalization is needed.

Significance. If fully established, this is a valuable negative result: it sharpens earlier perturbative analyses [19] by showing that nonperturbative effects make fermionic reheating even harder than previously estimated. The paper's strengths include the self-contained Dirac-equation computation with Bunch-Davies initial conditions, explicit phase-space distributions, lattice-informed inflaton fragmentation input, and solution of the full Boltzmann equation including Pauli blocking. The authors are also honest about their main caveats. However, the post-fragmentation impossibility claim is currently supported only up to y = 3×10^-7, more than three orders of magnitude below the perturbative BBN threshold; the central negative claim therefore needs either a threshold computation or a more conditional formulation.

major comments (3)
  1. [§4, Eq. (4.2), Figs. 8–10] The post-fragmentation impossibility claim rests on an extrapolation. The full Boltzmann computation (4.6) is presented for y = 10^-9, 10^-8, and 3×10^-7, far below the perturbative BBN threshold y ≈ 2.7×10^-4 from Eq. (4.2). The suppression seen at y = 3×10^-7 (more than an order of magnitude in ρψ) does not determine the behavior at y ≈ 3×10^-4, where the y^2 growth of the amplitude and a wider initially saturated region compete with the Pauli-blocking factors (1−fψ−fψ′). The statement that 'the trend is clear' (end of §4, p. 18) is an extrapolation, not a derivation. Please either extend the numerical integration to couplings at or above 2.7×10^-4, or explicitly restrict the conclusion to the computed range and soften the abstract's general 'reheating is not possible' claim.
  2. [§3.2, Fig. 7 and threshold y ≳ 0.2] The pre-fragmentation threshold y ≳ 0.2 is asserted without a shown coupling scan: Fig. 7 displays only y = 0.5, and the text jumps directly to the threshold value. In addition, this computation neglects backreaction of ψ on the inflaton dynamics (acknowledged in §3.2); the plotted ρφ, ρψ evolution does not conserve the total energy. Because the Dirac-equation background is fixed, the threshold is not a quantitative prediction. Please either show a systematic scan and estimate backreaction effects (e.g., from energy conservation or a simplified dissipative treatment), or label y ≳ 0.2 clearly as an order-of-magnitude guideline rather than a central quantitative result.
  3. [§5, Conclusions] The paper acknowledges that the impossibility conclusion assumes the absence of rapid thermalization of ψ with Standard Model degrees of freedom, and states that the thermalization timescale is unknown. Since fast thermalization would lift Pauli blocking and could restore a viable reheating channel, the central claim is conditional on an unquantified assumption. This conditionality should be stated in the abstract and in the final summary, not only in the concluding section.
minor comments (4)
  1. [Eq. (2.31)] The Hill equation is missing a factor of X_k on the right-hand side; it should read d²X_k/dz² + [(k/m_end)² + sn²(z/√6,−1)] X_k ≃ 0.
  2. [§2.2, first paragraph] Typo: 'homogeneities' should be 'inhomogeneities'.
  3. [§4, first paragraph] Typo: 'rapidly redfshifting' should be 'rapidly redshifting'.
  4. [§5] Typo: 'Heinsenberg equations' should be 'Heisenberg equations'.

Circularity Check

0 steps flagged

No significant circularity: the non-perturbative fermion production calculation is self-contained; the post-fragmentation impossibility claim is an extrapolation from computed low-coupling suppression, not a construction fitted to its own conclusion.

full rationale

The pre-fragmentation derivation chain is self-contained: the fermion mode functions are obtained by integrating the Dirac equations (2.21) with Bunch-Davies initial conditions (2.22), and the PSD and energy density are extracted from those mode functions via (2.27) and (2.26). No parameter is fitted to the y≳0.2 pre-fragmentation threshold; that threshold is read off from the computed ρψ versus ρφ crossing, with the paper explicitly noting that backreaction is not accounted for (Sec. 3.2). The post-fragmentation Boltzmann equation (4.6) is likewise solved with the non-perturbative PSD as initial condition and the lattice inflaton PSD as input; there is no inverse-engineering of the 'impossible' conclusion. Equation (4.2) from Ref. [19] is used as a perturbative benchmark, and the new low-coupling runs (Figs. 8–10) show Pauli suppression relative to that benchmark. The extension to couplings near and above y≃2.7×10^{-4} is an explicitly acknowledged extrapolation ('the trend is clear'), which is a robustness gap rather than a circular reduction. The same-author references [19,20,38] provide prior lattice constants, decay-rate formulas, and standard formalism; they are not an unverified uniqueness theorem or an ansatz that defines the present result. The paper's own stated limitations—unknown ψ thermalization (Sec. 5) and omitted backreaction in the pre-fragmentation calculation—also support reading the conclusion as a conservative guideline rather than a derivation that reduces to its inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 7 axioms · 0 invented entities

The paper introduces no new free parameters fitted here: λ is fixed by CMB normalization and the Yukawa couplings y, y5 are scanned model inputs. No new particles or forces are postulated. The central claim rests on standard background cosmology plus the domain assumptions of a single fermionic decay channel, absence of backreaction, and absence of prompt thermalization.

axioms (7)
  • domain assumption Quartic T-model potential V = λM_P^4 (√6 tanh(ϕ/√6M_P))^4 ≈ λϕ^4 near minimum, with λ ≃ 3.3×10^-12 fixed by CMB normalization.
    Section 2, Eqs. (2.2)-(2.4). The model choice supplies the λϕ^4 oscillatory regime and background evolution; the authors state qualitative robustness away from the minimum.
  • domain assumption Inflaton decay into a single Dirac fermion ψ through Yukawa couplings yϕψ̄ψ + iy5ϕψ̄γ5ψ is the only reheating channel; mψ ≪ mφ.
    Section 2, Eq. (2.8) and §3.1. The negative conclusion is conditional on absent bosonic decay channels and negligible fermion bare mass.
  • domain assumption Before fragmentation the inflaton condensate is spatially homogeneous; inhomogeneities enter only through δϕ and Floquet resonance.
    Section 2.1 Eqs. (2.9)-(2.10), Section 2.2 Eq. (2.28). Used for the background evolution and initial fermion production.
  • standard math Initial fermion modes are in the Bunch-Davies vacuum; normal-ordered Bogoliubov occupation numbers define ρψ.
    Section 2.1, Eqs. (2.22), (2.26)-(2.27); standard QFT-in-curved-spacetime setup.
  • domain assumption No backreaction of fermions on the inflaton in the pre-fragmentation non-perturbative computation.
    Section 3.2: "we have not accounted for the backreaction of ψ-production into the inflaton dynamics." This underpins the y≳0.2 threshold and is acknowledged by the authors.
  • domain assumption No prompt thermalization of the produced fermions; the Pauli-saturated spectrum persists.
    Conclusions: rapid thermalization would 'alleviate the saturation of the fermion energy levels', but "the time-scale for such thermalization is however unknown." If false, the central impossibility conclusion fails.
  • domain assumption After fragmentation, inflaton quanta decay as free particles with Γδϕ = y²mϕ/(8π); the lattice-derived inflaton PSD is reliable.
    Section 4, Eq. (4.1) and Fig. 2 from CosmoLattice; the decay is treated with a Boltzmann equation including quantum statistics.

pith-pipeline@v1.3.0-alltime-deepseek · 21072 in / 14984 out tokens · 146098 ms · 2026-08-03T15:36:04.740522+00:00 · methodology

0 comments
read the original abstract

Any viable inflationary model must account for reheating of the universe prior to the onset of primordial nucleosynthesis. In this work, we study the (p)reheating mechanism for an inflaton field with a quartic minimum of the T-model kind with coupling $\lambda$, prior to and post fragmentation, making a clear distinction between the two regimes. We assume that the main particle production channel corresponds to the decay into a pair of spin 1/2 fermions via Yukawa-like interactions. On top of its decays, we also consider the self-interaction of the inflaton, which sources the resonant growth of inflaton inhomogeneities, possibly leading to its eventual fragmentation. By means of a combination of non-perturbative (Heisenberg/Bogoliubov) and perturbative (Boltzmann) methods, we find that for Yukawa couplings that seemed to be intuitively perturbative, such as $y\gtrsim 10^{-8}$ ($y^2/\lambda\gtrsim 3\times10^{-5}$), parametric resonance, kinematic blocking, and Pauli suppression effects cannot be ignored. Additionally, we show that achieving $\rho_\phi \sim \rho_\psi$ prior to fragmentation requires large couplings, $y\gtrsim 0.2$ ($y^2/\lambda\gtrsim 10^{10}$), which needs a detailed study of backreaction and radiative corrections. Thus the rest of our work constitutes studying post-fragmentation fermion production where we conclude that, in general, reheating in this setup is not possible and thus we conclude that in order to successfully reheat, one must invoke a coupling to a integer- and/or 0-spin particle like a scalar boson.

discussion (0)

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

Forward citations

Cited by 5 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Leptogenesis and Low Reheating Temperatures

    hep-ph 2026-07 conditional novelty 7.0

    Standard non-thermal leptogenesis works at arbitrarily low reheating temperatures above the BBN bound when the inflaton potential has a quartic minimum (k≥4), because the inflaton's evolving mass kinematically shuts o...

  2. Asymmetric Reheating of Dark QED

    hep-ph 2026-05 unverdicted novelty 6.0

    Asymmetric reheating in Dark QED produces dark matter via a new channel where DM particles annihilate while still being created by inflaton decay, with the hidden-to-visible temperature ratio tied to the square root o...

  3. Graviton Production from Inflaton Condensate: Boltzmann vs Bogoliubov

    hep-ph 2026-04 unverdicted novelty 6.0

    For quadratic inflaton potentials Boltzmann and Bogoliubov spectra agree at short wavelengths, but for steeper potentials non-adiabatic transition effects captured only by Bogoliubov are sizable across a broad momentum range.

  4. Graviton Production from Inflaton Condensate: Boltzmann vs Bogoliubov

    hep-ph 2026-04 accept novelty 6.0

    For n=2 inflaton potentials Boltzmann and Bogoliubov agree on short-wavelength gravitons; for n>2 the non-adiabatic transition dominates and requires the Bogoliubov formalism.

  5. Thermal effects on Dark Matter production during cosmic reheating

    hep-ph 2026-04 unverdicted novelty 4.0

    Thermal corrections to reheating and freeze-in DM production rates are generally small in the computable regime but can be large in constructed counter-examples.

Reference graph

Works this paper leans on

57 extracted references · 48 linked inside Pith · cited by 4 Pith papers

  1. [1]

    Inflation,

    K. A. Olive, “Inflation,”Phys. Rept.190(1990) 307–403

    1990

  2. [2]

    A. D. Linde,Particle physics and inflationary cosmology, vol. 5. 1990.arXiv:hep-th/0503203. – 19 –

    Pith/arXiv arXiv 1990

  3. [3]

    Particle physics models of inflation and the cosmological density perturbation,

    D. H. Lyth and A. Riotto, “Particle physics models of inflation and the cosmological density perturbation,”Phys. Rept.314(1999) 1–146,arXiv:hep-ph/9807278

    Pith/arXiv arXiv 1999

  4. [4]

    Inflationary cosmology,

    A. D. Linde, “Inflationary cosmology,”Phys. Rept.333(2000) 575–591

    2000

  5. [5]

    A New Type of Isotropic Cosmological Models Without Singularity,

    A. A. Starobinsky, “A New Type of Isotropic Cosmological Models Without Singularity,”Phys. Lett. B91(1980) 99–102

    1980

  6. [6]

    Universality Class in Conformal Inflation,

    R. Kallosh and A. Linde, “Universality Class in Conformal Inflation,”JCAP07(2013) 002,arXiv:1306.5220 [hep-th]. [7]PlanckCollaboration, Y. Akramiet al., “Planck 2018 results. X. Constraints on inflation,”Astron. Astrophys.641 (2020) A10,arXiv:1807.06211 [astro-ph.CO]. [8]BICEP , KeckCollaboration, P. A. R. Adeet al., “Improved Constraints on Primordial Gr...

    Pith/arXiv arXiv 2013

  7. [9]

    Improved limits on the tensor-to-scalar ratio using BICEP and Planck data,

    M. Tristramet al., “Improved limits on the tensor-to-scalar ratio using BICEP and Planck data,”Phys. Rev. D 105no. 8, (2022) 083524,arXiv:2112.07961 [astro-ph.CO]. [10]ACTCollaboration, T. Louiset al., “The Atacama Cosmology Telescope: DR6 Power Spectra, Likelihoods and ΛCDM Parameters,”arXiv:2503.14452 [astro-ph.CO]. [11]ACTCollaboration, E. Calabreseet ...

    Pith/arXiv arXiv 2022

  8. [13]

    On Particle Creation by a Time Dependent Scalar Field,

    A. D. Dolgov and D. P. Kirilova, “On Particle Creation by a Time Dependent Scalar Field,”Sov. J. Nucl. Phys.51 (1990) 172–177

    1990

  9. [14]

    Particle Production During Out-of-equilibrium Phase Transitions,

    J. H. Traschen and R. H. Brandenberger, “Particle Production During Out-of-equilibrium Phase Transitions,” Phys. Rev. D42(1990) 2491–2504

    1990

  10. [15]

    Universe reheating after inflation,

    Y. Shtanov, J. H. Traschen, and R. H. Brandenberger, “Universe reheating after inflation,”Phys. Rev. D51(1995) 5438–5455,arXiv:hep-ph/9407247

    Pith/arXiv arXiv 1995

  11. [16]

    Reheating the postinflationary universe,

    D. Boyanovsky, M. D’Attanasio, H. J. de Vega, R. Holman, D. S. Lee, and A. Singh, “Reheating the postinflationary universe,”arXiv:hep-ph/9505220

    Pith/arXiv arXiv

  12. [17]

    Catastrophic particle production under periodic perturbation,

    M. Yoshimura, “Catastrophic particle production under periodic perturbation,”Prog. Theor. Phys.94(1995) 873–898,arXiv:hep-th/9506176

    Pith/arXiv arXiv 1995

  13. [18]

    Towards the theory of reheating after inflation,

    L. Kofman, A. D. Linde, and A. A. Starobinsky, “Towards the theory of reheating after inflation,”Phys. Rev. D56 (1997) 3258–3295,arXiv:hep-ph/9704452

    Pith/arXiv arXiv 1997

  14. [19]

    Reheating after inflaton fragmentation,

    M. A. G. Garcia and M. Pierre, “Reheating after inflaton fragmentation,”JCAP11(2023) 004,arXiv:2306.08038 [hep-ph]

    Pith/arXiv arXiv 2023

  15. [20]

    Effects of fragmentation on post-inflationary reheating,

    M. A. G. Garcia, M. Gross, Y. Mambrini, K. A. Olive, M. Pierre, and J.-H. Yoon, “Effects of fragmentation on post-inflationary reheating,”JCAP12(2023) 028,arXiv:2308.16231 [hep-ph]

    Pith/arXiv arXiv 2023

  16. [21]

    Impact of dark sector preheating on CMB observables,

    M. A. G. Garcia and A. Pereyra-Flores, “Impact of dark sector preheating on CMB observables,”JCAP08(2024) 043,arXiv:2403.04848 [astro-ph.CO]

    Pith/arXiv arXiv 2024

  17. [22]

    Structure of resonance in preheating after inflation,

    P. B. Greene, L. Kofman, A. D. Linde, and A. A. Starobinsky, “Structure of resonance in preheating after inflation,”Phys. Rev. D56(1997) 6175–6192,arXiv:hep-ph/9705347

    Pith/arXiv arXiv 1997

  18. [23]

    Preheating in an expanding universe: Analytic results for the massless case,

    D. I. Kaiser, “Preheating in an expanding universe: Analytic results for the massless case,”Phys. Rev. D56(1997) 706–716,arXiv:hep-ph/9702244

    Pith/arXiv arXiv 1997

  19. [24]

    Nonlinear inflaton fragmentation after preheating,

    G. N. Felder and L. Kofman, “Nonlinear inflaton fragmentation after preheating,”Phys. Rev. D75(2007) 043518, arXiv:hep-ph/0606256. – 20 –

    Pith/arXiv arXiv 2007

  20. [25]

    Preheating in the Standard Model with the Higgs-Inflaton coupled to gravity,

    J. Garcia-Bellido, D. G. Figueroa, and J. Rubio, “Preheating in the Standard Model with the Higgs-Inflaton coupled to gravity,”Phys. Rev. D79(2009) 063531,arXiv:0812.4624 [hep-ph]

    Pith/arXiv arXiv 2009

  21. [26]

    Non-linear Dynamics and Primordial Curvature Perturbations from Preheating,

    A. V. Frolov, “Non-linear Dynamics and Primordial Curvature Perturbations from Preheating,”Class. Quant. Grav.27(2010) 124006,arXiv:1004.3559 [gr-qc]

    Pith/arXiv arXiv 2010

  22. [27]

    Oscillons After Inflation,

    M. A. Amin, R. Easther, H. Finkel, R. Flauger, and M. P. Hertzberg, “Oscillons After Inflation,”Phys. Rev. Lett. 108(2012) 241302,arXiv:1106.3335 [astro-ph.CO]

    Pith/arXiv arXiv 2012

  23. [28]

    Theory of self-resonance after inflation. I. Adiabatic and isocurvature Goldstone modes,

    M. P. Hertzberg, J. Karouby, W. G. Spitzer, J. C. Becerra, and L. Li, “Theory of self-resonance after inflation. I. Adiabatic and isocurvature Goldstone modes,”Phys. Rev. D90(2014) 123528,arXiv:1408.1396 [hep-th]

    Pith/arXiv arXiv 2014

  24. [29]

    Nonperturbative Dynamics Of Reheating After Inflation: A Review,

    M. A. Amin, M. P. Hertzberg, D. I. Kaiser, and J. Karouby, “Nonperturbative Dynamics Of Reheating After Inflation: A Review,”Int. J. Mod. Phys. D24(2014) 1530003,arXiv:1410.3808 [hep-ph]

    Pith/arXiv arXiv 2014

  25. [30]

    Equation of State and Duration to Radiation Domination after Inflation,

    K. D. Lozanov and M. A. Amin, “Equation of State and Duration to Radiation Domination after Inflation,”Phys. Rev. Lett.119no. 6, (2017) 061301,arXiv:1608.01213 [astro-ph.CO]

    Pith/arXiv arXiv 2017

  26. [31]

    Self-resonance after inflation: oscillons, transients and radiation domination,

    K. D. Lozanov and M. A. Amin, “Self-resonance after inflation: oscillons, transients and radiation domination,” Phys. Rev. D97no. 2, (2018) 023533,arXiv:1710.06851 [astro-ph.CO]

    Pith/arXiv arXiv 2018

  27. [32]

    Parametric resonance in the early Universe—a fitting analysis,

    D. G. Figueroa and F. Torrenti, “Parametric resonance in the early Universe—a fitting analysis,”JCAP02(2017) 001,arXiv:1609.05197 [astro-ph.CO]

    Pith/arXiv arXiv 2017

  28. [33]

    Production of gravitational waves during preheating with nonminimal coupling,

    C. Fu, P. Wu, and H. Yu, “Production of gravitational waves during preheating with nonminimal coupling,”Phys. Rev. D97no. 8, (2018) 081303,arXiv:1711.10888 [gr-qc]

    Pith/arXiv arXiv 2018

  29. [34]

    Characterizing the postinflationary reheating history: Single daughter field with quadratic-quadratic interaction,

    S. Antusch, D. G. Figueroa, K. Marschall, and F. Torrenti, “Characterizing the postinflationary reheating history: Single daughter field with quadratic-quadratic interaction,”Phys. Rev. D105no. 4, (2022) 043532, arXiv:2112.11280 [astro-ph.CO]

    Pith/arXiv arXiv 2022

  30. [35]

    Dark matter production via a non-minimal coupling to gravity,

    O. Lebedev, T. Solomko, and J.-H. Yoon, “Dark matter production via a non-minimal coupling to gravity,”JCAP 02(2023) 035,arXiv:2211.11773 [hep-ph]

    arXiv 2023

  31. [36]

    Freeze-in from preheating,

    M. A. G. Garcia, K. Kaneta, Y. Mambrini, K. A. Olive, and S. Verner, “Freeze-in from preheating,”JCAP03 no. 03, (2022) 016,arXiv:2109.13280 [hep-ph]. [37]PlanckCollaboration, N. Aghanimet al., “Planck 2018 results. VI. Cosmological parameters,”Astron. Astrophys. 641(2020) A6,arXiv:1807.06209 [astro-ph.CO]. [Erratum: Astron.Astrophys. 652, C4 (2021)]

    Pith/arXiv arXiv 2022

  32. [38]

    Inflaton Oscillations and Post-Inflationary Reheating,

    M. A. G. Garcia, K. Kaneta, Y. Mambrini, and K. A. Olive, “Inflaton Oscillations and Post-Inflationary Reheating,”JCAP04(2021) 012,arXiv:2012.10756 [hep-ph]

    Pith/arXiv arXiv 2021

  33. [39]

    Constraints on Attractor Models of Inflation and Reheating from Planck, BICEP/Keck, ACT DR6, and SPT-3G Data,

    J. Ellis, M. A. G. Garcia, K. A. Olive, and S. Verner, “Constraints on Attractor Models of Inflation and Reheating from Planck, BICEP/Keck, ACT DR6, and SPT-3G Data,”arXiv:2510.18656 [hep-ph]

    arXiv

  34. [40]

    Coherent Scalar Field Oscillations in an Expanding Universe,

    M. S. Turner, “Coherent Scalar Field Oscillations in an Expanding Universe,”Phys. Rev. D28(1983) 1243

    1983

  35. [41]

    Reheating and Post-inflationary Production of Dark Matter,

    M. A. G. Garcia, K. Kaneta, Y. Mambrini, and K. A. Olive, “Reheating and Post-inflationary Production of Dark Matter,”Phys. Rev. D101no. 12, (2020) 123507,arXiv:2004.08404 [hep-ph]

    Pith/arXiv arXiv 2020

  36. [42]

    Production of massive fermions at preheating and leptogenesis,

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

    Pith/arXiv arXiv 1999

  37. [43]

    Preheating of fermions,

    P. B. Greene and L. Kofman, “Preheating of fermions,”Phys. Lett. B448(1999) 6–12,arXiv:hep-ph/9807339

    Pith/arXiv arXiv 1999

  38. [44]

    Matter creation via vacuum fluctuations in the early universe and observed ultrahigh-energy cosmic ray events,

    V. Kuzmin and I. Tkachev, “Matter creation via vacuum fluctuations in the early universe and observed ultrahigh-energy cosmic ray events,”Phys. Rev. D59(1999) 123006,arXiv:hep-ph/9809547

    Pith/arXiv arXiv 1999

  39. [45]

    Fermion production during preheating after hybrid inflation,

    J. Garcia-Bellido, S. Mollerach, and E. Roulet, “Fermion production during preheating after hybrid inflation,” JHEP02(2000) 034,arXiv:hep-ph/0002076

    Pith/arXiv arXiv 2000

  40. [46]

    Preheating of massive fermions after inflation: Analytical results,

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

    Pith/arXiv arXiv 2000

  41. [47]

    Coupled fields in external background with application to nonthermal production of gravitinos,

    H. P. Nilles, M. Peloso, and L. Sorbo, “Coupled fields in external background with application to nonthermal production of gravitinos,”JHEP04(2001) 004,arXiv:hep-th/0103202. – 21 –

    Pith/arXiv arXiv 2001

  42. [48]

    Gravitational Fermion Production in Inflationary Cosmology,

    D. J. H. Chung, L. L. Everett, H. Yoo, and P. Zhou, “Gravitational Fermion Production in Inflationary Cosmology,”Phys. Lett. B712(2012) 147–154,arXiv:1109.2524 [astro-ph.CO]

    Pith/arXiv arXiv 2012

  43. [49]

    Phenomenology of fermion production during axion inflation,

    P. Adshead, L. Pearce, M. Peloso, M. A. Roberts, and L. Sorbo, “Phenomenology of fermion production during axion inflation,”JCAP06(2018) 020,arXiv:1803.04501 [astro-ph.CO]

    Pith/arXiv arXiv 2018

  44. [50]

    The Big Bang, CPT, and neutrino dark matter,

    L. Boyle, K. Finn, and N. Turok, “The Big Bang, CPT, and neutrino dark matter,”Annals Phys.438(2022) 168767,arXiv:1803.08930 [hep-ph]

    Pith/arXiv arXiv 2022

  45. [51]

    On energy and particle production in cosmology: the particular case of the gravitino,

    G. Casagrande, E. Dudas, and M. Peloso, “On energy and particle production in cosmology: the particular case of the gravitino,”JHEP06(2024) 003,arXiv:2310.14964 [hep-th]

    Pith/arXiv arXiv 2024

  46. [52]

    Cosmological gravitational particle production and its implications for cosmological relics,

    E. W. Kolb and A. J. Long, “Cosmological gravitational particle production and its implications for cosmological relics,”Rev. Mod. Phys.96no. 4, (2024) 045005,arXiv:2312.09042 [astro-ph.CO]

    Pith/arXiv arXiv 2024

  47. [53]

    On the theory of fermionic preheating,

    P. B. Greene and L. Kofman, “On the theory of fermionic preheating,”Phys. Rev. D62(2000) 123516, arXiv:hep-ph/0003018

    Pith/arXiv arXiv 2000

  48. [54]

    Magnus and S

    W. Magnus and S. Winkler,Hill’s Equation. Dover Books on Mathematics. Dover Publications, Mineola, NY, Feb., 2004

    2004

  49. [55]

    The art of simulating the early Universe – Part I,

    D. G. Figueroa, A. Florio, F. Torrenti, and W. Valkenburg, “The art of simulating the early Universe – Part I,” JCAP04(2021) 035,arXiv:2006.15122 [astro-ph.CO]

    Pith/arXiv arXiv 2021

  50. [56]

    CosmoLattice,

    D. G. Figueroa, A. Florio, F. Torrenti, and W. Valkenburg, “CosmoLattice,”arXiv:2102.01031 [astro-ph.CO]

    Pith/arXiv arXiv

  51. [57]

    Inflationary Imprints on Dark Matter,

    S. Nurmi, T. Tenkanen, and K. Tuominen, “Inflationary Imprints on Dark Matter,”JCAP11(2015) 001, arXiv:1506.04048 [astro-ph.CO]

    Pith/arXiv arXiv 2015

  52. [58]

    Scalar dark matter production from preheating and structure formation constraints,

    M. A. G. Garcia, M. Pierre, and S. Verner, “Scalar dark matter production from preheating and structure formation constraints,”Phys. Rev. D107no. 4, (2023) 043530,arXiv:2206.08940 [hep-ph]

    Pith/arXiv arXiv 2023

  53. [59]

    How warm are non-thermal relics? Lyman-αbounds on out-of-equilibrium dark matter,

    G. Ballesteros, M. A. G. Garcia, and M. Pierre, “How warm are non-thermal relics? Lyman-αbounds on out-of-equilibrium dark matter,”JCAP03(2021) 101,arXiv:2011.13458 [hep-ph]

    Pith/arXiv arXiv 2021

  54. [60]

    Production and backreaction of fermions from axion-SU(2) gauge fields during inflation,

    L. Mirzagholi, A. Maleknejad, and K. D. Lozanov, “Production and backreaction of fermions from axion-SU(2) gauge fields during inflation,”Phys. Rev. D101no. 8, (2020) 083528,arXiv:1905.09258 [hep-th]

    Pith/arXiv arXiv 2020

  55. [61]

    Production and backreaction of massive fermions during axion inflation with non-Abelian gauge fields,

    P. Adshead, A. Liu, and K. D. Lozanov, “Production and backreaction of massive fermions during axion inflation with non-Abelian gauge fields,”JCAP09(2022) 043,arXiv:2203.09370 [hep-ph]

    Pith/arXiv arXiv 2022

  56. [62]

    Nonequilibrium dynamics of fermions in a spatially homogeneous scalar background field,

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

    Pith/arXiv arXiv 1998

  57. [63]

    Quantum theory of fermion production after inflation,

    J. Berges, D. Gelfand, and J. Pruschke, “Quantum theory of fermion production after inflation,”Phys. Rev. Lett. 107(2011) 061301,arXiv:1012.4632 [hep-ph]. – 22 –

    Pith/arXiv arXiv 2011