pith. sign in

arxiv: 2604.07326 · v2 · submitted 2026-04-08 · ✦ hep-ph · astro-ph.CO

Analytic Approximations for Fermionic Preheating

Heather E. Logan , Daniel Stolarski , Fazlul Yasin This is my paper

Pith reviewed 2026-05-10 17:49 UTC · model grok-4.3

classification ✦ hep-ph astro-ph.CO
keywords fermionic preheatingnon-perturbative fermion productionparametric resonanceanalytic approximationsnumber density scalingresonance peaksdark matter boundsinflation
0
0 comments X

The pith

Analytic approximations show fermion number density during preheating scales as the square root of coupling q for small values and q to the three-quarters power for large values.

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

This paper develops analytic approximations for non-perturbative fermion production in the preheating phase after λφ^4 inflation. The total number density of produced fermions follows power-law scalings with the coupling parameter q: proportional to q to the power 1/2 for q less than or equal to 0.01 and to q to the power 3/4 for q greater than or equal to 10. A simple relation is given to locate the resonance peaks in the fermion momentum spectrum for any q. For larger q the main contribution shifts to an approximately half-filled Fermi sphere from non-adiabatic effects at low momenta. If these fermions account for all dark matter the results yield lower bounds on their mass.

Core claim

The authors establish that resonance peak momenta in the fermion spectrum obey a simple predictive relation valid for arbitrary q. They derive analytic approximations for the total fermion number density that scale as q to the 1/2 power when q is at most 0.01 and as q to the 3/4 power when q is at least 10. In the regime q greater than or equal to 0.01 the dominant production mechanism is non-adiabatic filling of roughly half a Fermi sphere at low momenta, while for smaller q discrete high-momentum resonance peaks provide the main contribution.

What carries the argument

The resonance peaks in the fermion momentum spectrum together with the power-law approximations to the integrated number density, both derived as functions of the coupling parameter q during inflaton oscillations.

If this is right

  • The total fermion abundance can be estimated directly from q without integrating the full momentum spectrum in the two regimes.
  • Lower bounds on the mass of these fermions follow immediately if they constitute all dark matter.
  • Production is dominated by high-momentum resonances rather than low-momentum non-adiabatic effects when q is small.
  • The transition between resonance-dominated and Fermi-sphere-dominated regimes occurs near q equal to 0.01.

Where Pith is reading between the lines

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

  • The resonance-peak relation could be tested or extended in numerical simulations of other inflationary potentials beyond λφ^4.
  • The distinct momentum distributions in the two q regimes would produce different free-streaming lengths if the fermions are dark matter, affecting small-scale structure.
  • An interpolation formula bridging the small-q and large-q power laws might be constructed for intermediate couplings around 0.1 to 1.

Load-bearing premise

The analytic approximations remain accurate across the stated ranges of q and that the produced fermions can make up the entire dark matter density when mass lower bounds are derived.

What would settle it

A high-precision numerical integration of the fermion production equations for q equal to 0.001 that yields a number density deviating from the predicted square-root scaling, or for q equal to 100 deviating from the three-quarter-power scaling, would falsify the approximations.

Figures

Figures reproduced from arXiv: 2604.07326 by Daniel Stolarski, Fazlul Yasin, Heather E. Logan.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Numerically evaluated function [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p012_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p014_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Feynman diagrams describing the perturbative process [PITH_FULL_IMAGE:figures/full_fig_p015_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: General form of Feynman diagram describing the perturbative process [PITH_FULL_IMAGE:figures/full_fig_p016_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Contours corresponding to [PITH_FULL_IMAGE:figures/full_fig_p017_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: Contributions of different components as fractions of the total integral [PITH_FULL_IMAGE:figures/full_fig_p019_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10: Contributions of different components and the total integral for 10 [PITH_FULL_IMAGE:figures/full_fig_p020_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11 [PITH_FULL_IMAGE:figures/full_fig_p022_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12: Contributions of different components as fractions of the total integral [PITH_FULL_IMAGE:figures/full_fig_p024_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13: Different integrals as a function of [PITH_FULL_IMAGE:figures/full_fig_p025_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14: Total integral [PITH_FULL_IMAGE:figures/full_fig_p026_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: FIG. 15: Contributions of different components for 0 [PITH_FULL_IMAGE:figures/full_fig_p027_15.png] view at source ↗
read the original abstract

Non-thermal fermions can be produced non-perturbatively in the early universe during coherent oscillations of a scalar field. We explore fermion production in $\lambda\phi^{4}$ inflation through this mechanism and analyze the momentum spectrum of the fermions produced, which depends on a coupling parameter $q$. For $q \gtrsim 0.01$, the main contribution to the total number density comes from an approximately half-filled Fermi sphere as a result of non-adiabaticity. For $q\lesssim 0.01$, we find that the major contributions instead come from resonance peaks at higher momentum values. We find a simple relation to predict the momentum values corresponding to resonance peaks for any $q$. We also obtain analytic power-law approximations for the total number density of fermions and find that it is proportional to $q^{1/2}$ for $q\lesssim 0.01$ and proportional to $q^{3/4}$ for $q\gtrsim 10$. If fermions produced by this mechanism make up the entirety of dark matter, we estimate lower bounds on their mass.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 3 minor

Summary. The paper explores non-thermal production of fermions in λφ⁴ inflation via preheating, focusing on the momentum spectrum's dependence on the coupling q. It identifies two regimes: for q ≳ 0.01, production is dominated by an approximately half-filled Fermi sphere from non-adiabaticity, and for q ≲ 0.01, by resonance peaks. A simple predictive relation for resonance peak momenta is provided for any q. Analytic approximations yield total number density n ∝ q^{1/2} for q ≲ 0.01 and n ∝ q^{3/4} for q ≳ 10. Lower bounds on fermion mass are estimated assuming they comprise all dark matter.

Significance. Should the approximations prove accurate upon numerical validation, this work would provide practical analytic tools for calculating fermion yields in preheating, reducing reliance on computationally intensive simulations. The power-law relations and peak prediction formula could be particularly useful in exploring parameter spaces for fermionic dark matter models in the early universe.

major comments (3)
  1. [analytic approximations for total number density] The q^{1/2} scaling for q≲0.01 is obtained by approximating the sum over resonance peaks (as described in the regime analysis), but the manuscript lacks a direct comparison to numerical integration of the mode equations or error estimates, making it difficult to confirm the exponent remains accurate near the quoted boundaries.
  2. [Fermi sphere regime analysis] For q≳10, the n ∝ q^{3/4} scaling relies on modeling the low-momentum band as a half-filled Fermi sphere whose radius scales with q; the derivation does not quantify how expansion or back-reaction might alter the filling factor, which could invalidate the exponent if the adiabaticity assumption fails inside the regime.
  3. [resonance peaks section] The simple relation for resonance peak momenta is presented as valid for any q, yet the underlying approximations (including on the adiabaticity parameter) and tests at intermediate q values are not shown explicitly, leaving the generality of the relation unverified.
minor comments (3)
  1. The abstract states the regime boundaries at q~0.01 and q~10 without indicating how these thresholds were determined from the spectra; a brief justification in the main text would improve clarity.
  2. Adding a table or figure directly comparing the analytic n(q) to numerical results for several benchmark q values would strengthen the power-law claims.
  3. Notation for the mode functions and occupation numbers could be defined more explicitly at first use to aid readability.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments, which help improve the clarity and robustness of our analytic approximations. We address each major comment below and indicate the revisions planned for the manuscript.

read point-by-point responses

  1. Referee: [analytic approximations for total number density] The q^{1/2} scaling for q≲0.01 is obtained by approximating the sum over resonance peaks (as described in the regime analysis), but the manuscript lacks a direct comparison to numerical integration of the mode equations or error estimates, making it difficult to confirm the exponent remains accurate near the quoted boundaries.

    Authors: We agree that explicit numerical validation would strengthen the result. The q^{1/2} scaling follows directly from summing the analytically derived contributions of the resonance peaks, whose locations, widths, and heights are set by the resonance condition. In the revised manuscript we will add a direct comparison between the analytic sum and numerical integration of the mode equations for several representative values of q near the q=0.01 boundary, together with relative-error estimates to confirm that the exponent remains accurate throughout the quoted regime. revision: yes

  2. Referee: [Fermi sphere regime analysis] For q≳10, the n ∝ q^{3/4} scaling relies on modeling the low-momentum band as a half-filled Fermi sphere whose radius scales with q; the derivation does not quantify how expansion or back-reaction might alter the filling factor, which could invalidate the exponent if the adiabaticity assumption fails inside the regime.

    Authors: The half-filled Fermi-sphere model follows from the non-adiabatic production mechanism, which drives the occupation number to 1/2 for all modes below a q-dependent cutoff. For q≳10 the production occurs on a timescale much shorter than the Hubble time, so expansion during the resonance window is negligible; back-reaction remains small while the fermion energy density is still sub-dominant. We acknowledge, however, that a quantitative bound on possible deviations of the filling factor would be useful. In the revision we will add an explicit estimate of the adiabaticity parameter across the q≳10 regime and a brief discussion of the expected size of expansion and back-reaction corrections, together with the regime of validity of the q^{3/4} scaling. revision: partial

  3. Referee: [resonance peaks section] The simple relation for resonance peak momenta is presented as valid for any q, yet the underlying approximations (including on the adiabaticity parameter) and tests at intermediate q values are not shown explicitly, leaving the generality of the relation unverified.

    Authors: The relation is obtained by setting the adiabaticity parameter to O(1), which locates the centers of the resonance bands independently of the specific value of q. The underlying approximation is therefore general. To make this explicit, the revised manuscript will include the explicit expression for the adiabaticity parameter, together with numerical checks of the predicted peak locations at intermediate couplings (e.g., q=1) to verify that the formula remains accurate outside the two limiting regimes already discussed. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivations of resonance-peak relation and power-law scalings are self-contained analytic approximations

full rationale

The paper derives a predictive relation for resonance-peak momenta and the two power-law regimes for integrated fermion number density by approximating the mode functions, occupation-number integral, and resonance structure directly from the equations of motion in the λφ⁴ model. No step reduces a claimed prediction to a fitted parameter by construction, nor does any load-bearing premise rest on a self-citation chain or imported uniqueness theorem. The q-regime boundaries and adiabaticity assumptions are stated explicitly as conditions of validity rather than being smuggled in via prior work by the same authors. The central results therefore remain independent of the inputs they are applied to.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on the standard assumptions of inflationary cosmology and non-perturbative fermion production in an oscillating scalar background; no new entities or fitted constants beyond the input coupling q are introduced.

axioms (2)
  • domain assumption The background evolution is given by coherent oscillations of a scalar field with λφ⁴ potential during preheating.
    The entire analysis is performed inside this specific inflationary model.
  • standard math Fermion production is governed by the standard Dirac equation in a time-dependent scalar background.
    This is the usual QFT setup for preheating calculations.

pith-pipeline@v0.9.0 · 5492 in / 1459 out tokens · 30824 ms · 2026-05-10T17:49:06.636421+00:00 · methodology

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

41 extracted references · 41 canonical work pages · 2 internal anchors

  1. [1]

    Significant gravitational waves may also be generated as a consequence of the preheating mechanism [16]

    and leptogenesis [15] scenarios to explain the observed baryon asymmetry of the universe. Significant gravitational waves may also be generated as a consequence of the preheating mechanism [16]. Preheating can also lead to the production of particles out of thermal equilibrium that may act as a viable candidate for dark matter [17, 18]. For a general revi...

    work page

  2. [2]

    (4.9) for the maximumκgiven by the inequality on the right of Eq

    We usef ∗ from Eq. (4.9) for the maximumκgiven by the inequality on the right of Eq. (4.8) and see for the corresponding valueκ ∗: qf∗ ˙f∗ = κ2 ∗ +qf 2 ∗ 3/2 =⇒κ ∗ = 2q1/2 ˙f∗ 3 √ 3 2 27 1/4 q1/4 ≈0.52q 1/4.(4.10) Thisκ ∗ ≈0.52q 1/4 implies that the boundary of the bulk is∝q 1/4. If the major contribution to the fermion number density comes from the bulk ...

    work page

  3. [3]

    Reheating after Inflation

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

    work page Pith review arXiv 1994

  4. [4]

    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, 172 (1990)

    work page 1990

  5. [5]

    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, 2491 (1990)

    work page 1990

  6. [6]

    Universe Reheating after Inflation

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

    work page Pith review arXiv 1995

  7. [7]

    Preheating of fermions,

    P. B. Greene and L. Kofman, “Preheating of fermions,” Phys. Lett. B448, 6 (1999), hep- ph/9807339. 5 This case corresponds to considering bosons. 38

    work page arXiv 1999

  8. [8]

    Baryon Asymmetry in Inflationary Universe,

    A. D. Dolgov and A. D. Linde, “Baryon Asymmetry in Inflationary Universe,” Phys. Lett. B 116, 329 (1982)

    work page 1982

  9. [9]

    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, 3258 (1997), hep-ph/9704452

    work page Pith review arXiv 1997

  10. [10]

    A New method in the theory of superconductivity,

    N. N. Bogolyubov, V. V. Tolmachev, and D. V. Shirkov, “A New method in the theory of superconductivity,” Fortsch. Phys.6, 605 (1958)

    work page 1958

  11. [11]

    Particle creation from vacuum by homogeneous electric field with a periodical time dependence,

    V. M. Mostepanenko and V. M. Frolov, “Particle creation from vacuum by homogeneous electric field with a periodical time dependence,” Yad. Fiz.19, 885 (1974)

    work page 1974

  12. [12]

    Kolb and A.J

    E. W. Kolb and A. J. Long, “Cosmological gravitational particle production and its implica- tions for cosmological relics,” Rev. Mod. Phys.96, 045005 (2024), 2312.09042

    work page arXiv 2024

  13. [13]

    Fermion production during and after axion inflation

    P. Adshead and E. I. Sfakianakis, “Fermion production during and after axion inflation,” JCAP11, 021 (2015), 1508.00891

    work page Pith review arXiv 2015

  14. [14]

    Adshead, L

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

    work page arXiv 2018

  15. [15]

    Inflation Dynamics and Reheating

    B. A. Bassett, S. Tsujikawa, and D. Wands, “Inflation dynamics and reheating,” Rev. Mod. Phys.78, 537 (2006), astro-ph/0507632

    work page Pith review arXiv 2006

  16. [16]

    GUT baryogenesis after preheating,

    E. W. Kolb, A. D. Linde, and A. Riotto, “GUT baryogenesis after preheating,” Phys. Rev. Lett.77, 4290 (1996), hep-ph/9606260

    work page arXiv 1996

  17. [17]

    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, 014 (1999), hep-ph/9905242

    work page arXiv 1999

  18. [18]

    Cosmological Backgrounds of Gravitational Waves

    C. Caprini and D. G. Figueroa, “Cosmological Backgrounds of Gravitational Waves,” Class. Quant. Grav.35, 163001 (2018), 1801.04268

    work page Pith review arXiv 2018

  19. [19]

    Nonthermal supermassive dark matter,

    D. J. H. Chung, E. W. Kolb, and A. Riotto, “Nonthermal supermassive dark matter,” Phys. Rev. Lett.81, 4048 (1998), hep-ph/9805473

    work page arXiv 1998

  20. [20]

    Carena, N

    M. Carena, N. M. Coyle, Y.-Y. Li, S. D. McDermott, and Y. Tsai, “Cosmologically degenerate fermions,” Phys. Rev. D106, 083016 (2022), 2108.02785

    work page arXiv 2022

  21. [21]

    Reheating in Inflationary Cosmology: Theory and Applications

    R. Allahverdi, R. Brandenberger, F.-Y. Cyr-Racine, and A. Mazumdar, “Reheating in In- flationary Cosmology: Theory and Applications,” Ann. Rev. Nucl. Part. Sci.60, 27 (2010), 1001.2600

    work page Pith review arXiv 2010

  22. [22]

    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, 1530003 (2014), 1410.3808

    work page Pith review arXiv 2014

  23. [23]

    Chaotic Inflation,

    A. D. Linde, “Chaotic Inflation,” Phys. Lett. B129, 177 (1983)

    work page 1983

  24. [24]

    On the Theory of Fermionic Preheating

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

    work page Pith review arXiv 2000

  25. [25]

    Aspects of Preheating After Inflation,

    P. B. Greene, “Aspects of Preheating After Inflation,” Ph.D. thesis, Toronto U. (2002), URL http://hdl.handle.net/1807/121431

    work page 2002

  26. [26]

    Planck 2018 results. X. Constraints on inflation

    Y. Akrami et al. (Planck), “Planck 2018 results. X. Constraints on inflation,” Astron. Astro- phys.641, A10 (2020), 1807.06211

    work page internal anchor Pith review arXiv 2018

  27. [27]

    Fermion production during preheating after hybrid inflation,

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

    work page arXiv 2000

  28. [28]

    On initial conditions for the Hot Big Bang

    F. Bezrukov, D. Gorbunov, and M. Shaposhnikov, “On initial conditions for the Hot Big Bang,” JCAP06, 029 (2009), 0812.3622

    work page Pith review arXiv 2009

  29. [29]

    Garcia, K

    M. A. G. Garcia, K. Kaneta, Y. Mambrini, K. A. Olive, and S. Verner, “Freeze-in from preheating,” JCAP03, 016 (2022), 2109.13280

    work page arXiv 2022

  30. [30]

    Non-perturbative production of fermionic dark matter from fast preheating,

    J. Klaric, A. Shkerin, and G. Vacalis, “Non-perturbative production of fermionic dark matter from fast preheating,” JCAP02, 034 (2023), 2209.02668. 39

    work page arXiv 2023

  31. [31]

    WIMPzillas!,

    E. W. Kolb, D. J. H. Chung, and A. Riotto, “WIMPzillas!,” AIP Conf. Proc.484, 91 (1999), hep-ph/9810361

    work page internal anchor Pith review arXiv 1999

  32. [32]

    Dark Matter Decaying into a Fermi Sea of Neutrinos,

    O. E. Bjaelde and S. Das, “Dark Matter Decaying into a Fermi Sea of Neutrinos,” Phys. Rev. D82, 043504 (2010), 1002.1306

    work page arXiv 2010

  33. [33]

    Degenerate fermion dark matter from a broken U(1) B−L gauge symmetry,

    G. Choi, M. Suzuki, and T. T. Yanagida, “Degenerate fermion dark matter from a broken U(1) B−L gauge symmetry,” Phys. Rev. D102, 035022 (2020), 2004.07863

    work page arXiv 2020

  34. [34]

    Degenerate Sub-keV Fermion Dark Matter from a Solution to the Hubble Tension,

    G. Choi, M. Suzuki, and T. T. Yanagida, “Degenerate Sub-keV Fermion Dark Matter from a Solution to the Hubble Tension,” Phys. Rev. D101, 075031 (2020), 2002.00036

    work page arXiv 2020

  35. [35]

    E. W. Kolb and M. S. Turner,The Early Universe, vol. 69 (Taylor and Francis, 2019), ISBN 978-0-429-49286-0, 978-0-201-62674-2

    work page 2019

  36. [36]

    JacobiCN,

    Wolfram Research, “JacobiCN,”,https://reference.wolfram.com/language/ref/ JacobiCN.html(1988)

    work page 1988

  37. [37]

    Abramowitz and I

    M. Abramowitz and I. A. Stegun,Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables(Dover, New York City, 1964), ninth Dover printing, tenth GPO printing ed

    work page 1964

  38. [38]

    Hybrid Inflation

    A. D. Linde, “Hybrid inflation,” Phys. Rev. D49, 748 (1994), astro-ph/9307002

    work page Pith review arXiv 1994

  39. [39]

    Cosmology of the Invisible Axion,

    J. Preskill, M. B. Wise, and F. Wilczek, “Cosmology of the Invisible Axion,” Phys. Lett. B 120, 127 (1983)

    work page 1983

  40. [40]

    A Cosmological Bound on the Invisible Axion,

    L. F. Abbott and P. Sikivie, “A Cosmological Bound on the Invisible Axion,” Phys. Lett. B 120, 133 (1983)

    work page 1983

  41. [41]

    Thomson,Modern particle physics(Cambridge University Press, New York, 2013), ISBN 978-1-107-03426-6, 978-1-139-52536-7

    M. Thomson,Modern particle physics(Cambridge University Press, New York, 2013), ISBN 978-1-107-03426-6, 978-1-139-52536-7

    work page 2013