Pith. sign in

REVIEW 2 major objections 4 minor 55 references

Particle decay as asymptotic narrow parametric resonance

T0 review · 2 major / 4 minor · reviewed 2026-08-15 · deepseek-v4-flash

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

Pith's one-line read A Gaussian-smoothed Boltzmann decay equation reproduces the asymptotic narrow parametric resonance growth rate exactly at the peak, and matches the full number density to within 20 percent.

desk verdict A clean derivation of the NPR rate from the Boltzmann equation, but the 20% number-density agreement is an artifact of the Gaussian ansatz rather than a prediction. read the letter →

arxiv 2505.03160 v2 pith:OXBZXH7G submitted 2025-05-06 hep-ph

Shao-Ping Li This is my paper · ORCID
classification hep-ph
keywords narrowparametricresonanceBoltzmannequationstimulateddecayGaussianmomentumspreadMathieuscalarfieldaxionphotonproductiondarkmatter
topics Dark Matter
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper claims that the ordinary Boltzmann equation for a trilinear decay process $\phi\to 2\chi$, with the energy-conserving delta function replaced by a Gaussian of width $\sigma = q m_\phi$, reproduces the asymptotic exponential growth of narrow parametric resonance (NPR). On the resonance peak the Boltzmann distribution becomes $f_\chi(z)=\frac{1}{2}e^{Gqz}-\frac{1}{2}$ with $G=1$, so the growth rate is $q m_\phi/2$, exactly the Mathieu-equation rate. After the occupation number exceeds one, the integrated number density from the Boltzmann equation tracks the Mathieu result to a constant ratio $n^M_\chi/n^B_\chi\approx 0.8$, nearly independent of $q$ for $q\lesssim 0.1$. If true, the simple particle-decay picture gives a quantitative and analytic handle on explosive boson production from oscillating scalar fields, replacing earlier Boltzmann treatments that mispredicted the density by orders of magnitude.

What carries the argument

The load-bearing object is the Gaussian replacement of the energy-conserving delta function, $\delta(m_\phi-2k_1)\approx (a/\sigma)\exp\left(-(k_1-m_\phi/2)^2/(b\sigma^2)\right)$, with the width $\sigma = q m_\phi$ fixed by the momentum spread inherited from nonrelativistic $\phi$ decay. This ansatz turns a singular decay kernel into a finite resonance band of width $q m_\phi$ centered at $k_1=m_\phi/2$, making the Boltzmann equation solvable in closed form and producing the factor $G(z)$ that matches the Mathieu band shape. The coefficients $a=2/\pi$, $b=\pi/16$ are set by normalization, and the decay channel $j=1$ dominates because higher-$j$ amplitudes are suppressed by powers of $(q/4)^{2j}$ for $q\ll 1$.

What would settle it

Take a fixed narrow-resonance value such as $q=0.01$ and solve the Mathieu equation numerically; if the late-time peak profile of $n_k$ is not well fitted by a Gaussian of width $q m_\phi$, or if the ratio $n^M_\chi/n^B_\chi$ does not approach $0.8$ as $z\to\infty$ (within, say, 10 percent), the central claim fails. A simpler variant: compute $n^B_\chi$ with a Lorentzian or boxcar line shape of the same width and check whether the asymptotic ratio still tends to a constant.

Watch

Extended reading notes

Core claim

The central claim is that stimulated decay in the Boltzmann equation, supplemented by a Gaussian model of the decay-produced momentum spread, is quantitatively equivalent to narrow parametric resonance in its asymptotic regime. The paper works with the trilinear coupling $\mathcal{L}=-\frac{1}{2}\mu\phi\chi^2$ and uses the Mathieu equation for the $\chi$ mode as the exact reference. It simulates the delta function $\delta(m_\phi-2k_1)$ by a Gaussian of width $\sigma = q m_\phi$, with $q\equiv 2\mu\bar\phi/m_\phi^2$, and chooses normalization coefficients $a=2/\pi$, $b=\pi/16$ so that the Boltzmann solution is $f_\chi(z)=\frac{1}{2}e^{Gqz}-\frac{1}{2}$, where $G=\exp\left(-16(k_1-m_\phi/2)^2/(\pi q^2 m_\phi^2)\right)$. At the peak $G=1$, giving the NPR growth rate $q m_\phi/2$, and the integrated number density $n^B_\chi$ matches the Mathieu density $n^M_\chi$ to a constant ratio $0.8$ that is essentially independent of $q$ in the narrow regime. The paper therefore concludes that the Boltzmann particle picture, already used for late-time thermalization, can be trusted for the explosive production phase.

Load-bearing premise

The Gaussian line shape with width $\sigma = q m_\phi$ is an ansatz: the momentum spread of produced particles is asserted rather than derived, and the quantitative match, including the 20 percent density ratio, holds only if the true resonance shape is this Gaussian.

Editorial extensions

If this is right

  • For $q\lesssim 0.1$, the $j=1$ decay channel alone reproduces the asymptotic NPR distribution; higher-order annihilation channels are negligible in the narrow regime.
  • The integrated particle number density from the Boltzmann equation is within 20 percent of the Mathieu result, so density-based observables can be computed analytically rather than by solving the Mathieu equation.
  • Explosive production works only when the condition $q^2\gg H_{\rm osc}/m_\phi$ holds, and with fast scattering the stronger condition $q^2\gg \xi H_{\rm osc}/m_\phi$ is needed; otherwise redshift and scattering erase the bandwidth before growth develops.
  • For oscillations beginning after MeV temperatures, an eV-scale mass readily satisfies these conditions, so axion-photon conversion and dark-matter production can be studied with the Boltzmann approach and connected to BBN and CMB observables.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • One could test whether the Gaussian width $\sigma = q m_\phi$ is exact by computing the one-loop decay line shape of a nonrelativistic condensate; if the true line shape is Lorentzian or asymmetric, the 20 percent density statement would need revision.
  • The constant ratio $0.8$ may be derivable analytically as a shape factor between the Gaussian envelope and the actual Floquet mode profile; such a derivation would sharpen the approximation and predict where it degrades as $q$ approaches the edge of the first stability band.
  • The same Gaussian Boltzmann construction could be extended to vector or fermionic decay products, where the Bose-enhancement factor changes, to see whether the 20 percent agreement persists.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 4 minor

Summary. The paper studies the relation between narrow parametric resonance (NPR) in a trilinear scalar interaction and the Boltzmann equation for stimulated decay ϕ → 2χ. The author solves the Mathieu equation for the χ mode functions and compares the occupation number n_k with the distribution function obtained from the Boltzmann equation, in which the energy-conserving delta function is replaced by a Gaussian of width σ = q m_ϕ (Eqs. (29)–(36)). The main results are: (i) the Boltzmann distribution function has the same peak growth rate q m_ϕ/2 as NPR; (ii) the integrated number density ratio n^M_χ/n^B_χ approaches a constant ≈ 0.8 independent of q; and (iii) the Boltzmann approach can be used in cosmological scenarios where expansion, backreaction, and scattering are negligible. The paper concludes that the Boltzmann particle-decay picture is quantitatively valid in the asymptotic NPR regime, unlike the order-of-magnitude discrepancy reported in Ref. [24].

Significance. If substantiated, the result would give a simple analytic Boltzmann description of explosive particle production from oscillating fields and would reconcile the S-matrix-based Boltzmann equation with Floquet theory in the narrow-resonance limit. The paper's explicit solution for f_χ, the peak-growth-rate matching, and the verification of the Feynman amplitudes with FeynArts/FeynCalc are clean and useful. However, the central quantitative claim rests on an ad hoc Gaussian line shape; the constant 0.8 is not a prediction of the Boltzmann collision term. The paper therefore does not yet establish the quantitative equivalence it advertises, though the approach is promising and the needed test is well defined.

major comments (2)
  1. [Section III, Eqs. (29)–(36) and Fig. 2] The asymptotic ratio n^M_χ/n^B_χ ≈ 0.8 is fixed by the Gaussian line-shape ansatz and is not a property of the Boltzmann collision term. For qz >> 1, the Mathieu occupation near the first band grows as exp(qz(1 − 8 r^2)) with r = (k − m_ϕ/2)/(q m_ϕ), whereas the Gaussian ansatz of Eq. (36) gives exp(qz(1 − (16/π) r^2)). The saddle-point ratio of the two momentum integrals is sqrt((16/π)/8) = sqrt(2/π) ≈ 0.798, independent of q and z, in agreement with the plateau in Fig. 2. After imposing the normalization (30) and the peak-rate matching, the combination b σ^2 is fixed, so the solution (35) is independent of the stated value of σ; the only free input is the Gaussian form in Eq. (29). Thus the 20% offset measures the difference between the assumed Gaussian curvature and the actual Mathieu band curvature, not the accuracy of the Boltzmann equation. The authors should either derive the Gaussian line shape from the QFT/Mathieu solution or insert the exact Floquet profile into the Boltzmann equation and check whether the ratio approaches 1.
  2. [Section III, Eq. (33)] The identification of the momentum spread is not supported. From Eq. (9), the first instability band is |A_k − 1| < q, i.e., |k − m_ϕ/2| < q m_ϕ/4, whereas Eq. (32) gives σ = q m_ϕ, a factor-of-four difference from the Mathieu half-width. Moreover, as noted above, σ cancels from the final distribution function once the normalization and peak-rate conditions are imposed. The paper should state clearly that the Gaussian is an assumption and address how the line shape could be fixed by the underlying theory, rather than presenting it as a 'proper simulation' of the decay kinematics.
minor comments (4)
  1. [Section III, Eq. (38) and Fig. 2] The value of the cutoff n is not specified; the authors should state the value used and demonstrate that the plotted ratio is insensitive to the choice for the range of z shown.
  2. [References] Reference [35] appears in the bibliography but is not cited in the text; please add the citation or remove the entry.
  3. [Section II, Eq. (8)] The sign in Eq. (8) appears inconsistent with the Lagrangian (4): for L = −1/2 μϕχ^2, the equation of motion gives ω_k^2 = k^2 + μϕ, not k^2 − μϕ. This can be absorbed by a shift of the oscillation phase, but the text should be made consistent.
  4. [Fig. 1 caption] The caption should specify the time (or z) at which the top-panel momentum profiles are plotted, and the bottom-panel label 'rationk/fχ' should read 'ratio n_k/f_χ'.

Circularity Check

2 steps flagged · score 6.0 of 10

The 20% number-density offset is fixed by the hand-set Gaussian curvature b=π/16 and the width σ=qmφ is imported from the Mathieu bandwidth; the exponential rate qmφ/2 itself is not circular.

  1. self definitional [Section III, Eq. (33) (Gaussian width from Mathieu equation)]
    "Therefore, the momentum has a spread σ≡δk1=µ¯ϕ/k1=qmϕ, (33), which is the same as that derived from Eq. (9) by using e.g., the harmonic balancing [19]."

    The Gaussian width inserted into the Boltzmann equation is not obtained from the S-matrix decay kinematics of the collision term, but is read off from Eq. (9), the Mathieu equation of the very narrow parametric resonance whose asymptotic regime the Boltzmann treatment is meant to reproduce. The momentum scale and band location of the produced particles are therefore inputs from the target calculation, not independent outputs. The peak growth rate and the occupation value are not forced by this step alone, so the circularity is partial.

  2. fitted input called prediction [Section III, Eqs. (29)-(36) and Fig. 2]
    "Substituting Eq. (29) into Eq. (27) with the following simulation coefficients a = 2/π, b = π/16, we arrive at fχ(z) = 1/2 e^{Gqz} − 1/2 ... We see that after the ratio nk/fχ reaches 1, the ratio nMχ/nBχ will reach a constant at 0.8, basically independent of q."

    With σ=qmφ and r≡(k−mφ/2)/(q mφ), Eq. (36) gives G=exp(−16r^2/π), so nB is a saddle integral over exp(qz exp(−16r^2/π)). The first Mathieu band grows as exp(qz√(1−16r^2)) ≈ exp(qz(1−8r^2)) near the peak. For qz≫1, nM/nB → √((16/π)/8)=√(2/π)≈0.798, independent of q and z. Thus the announced 20% offset is the ratio of the Gaussian curvature set by the unexplained parameter b=π/16 to the Mathieu curvature; it is an output of the chosen line shape, not of the Boltzmann collision integral. The normalization condition (31) fixes only a√b, leaving b free.

full rationale

The paper is not globally circular. The exponential rate qmφ/2 follows from the j=1 decay term of the Boltzmann equation (27) once q=2μ̄ϕ/m_φ^2 is identified, and the peak distribution comparison in Fig. 1 is an external numerical check against the Mathieu equation, not a self-citation loop. However, the headline quantitative statement—the constant 20% deviation of the integrated number density—is not an independent prediction. The Gaussian delta simulation (29) has two coefficients; the normalization condition (31) fixes only their product, and the choice b=π/16 in Eq. (34) is unexplained. This choice sets the curvature of G in Eq. (36), and the saddle-point ratio of the Mathieu to Boltzmann integrals is √((16/π)/8)=√(2/π)≈0.798, exactly the 0.8 reported in Fig. 2. Moreover, the width σ=qmφ is explicitly taken from the Mathieu equation (9), i.e., from the NPR calculation whose asymptotic regime is being reproduced. Hence the integrated-density claim reduces by construction to the Gaussian ansatz, while the exponential-rate claim retains independent content; this yields a partial circularity score of 6 rather than a higher score.

Assumptions & free parameters 2 free parameters · 5 assumptions · 0 invented entities

The central claim rests on a small set of standard domain assumptions plus one ad hoc shape choice: the Gaussian momentum spread. The width sigma is not fitted to the final growth rate, but it is taken from the Mathieu bandwidth, so it is not fully independent of the target calculation. The unspecified cutoff n is a second, smaller parameter. No new particles, forces, or entities are introduced.

free parameters (2)
  • Gaussian width sigma (momentum spread) = sigma = q m_phi
    The width of the Gaussian replacing the energy delta function is set equal to the Mathieu first-bandwidth (Eq. 33). It is an input imported from the NPR calculation, and the central match depends on it.
  • Integration cutoff n = O(1), unspecified
    Eq. (38) computes n^B_chi by integrating only |k - m_phi/2| < n sigma. The numerical value of n is never given; the claimed 20 percent density ratio may depend on it.
assumptions (5)
  • domain assumption The Mathieu equation (9) with q << 1 accurately describes the chi mode evolution.
    Standard in the NPR literature [2,19]; the paper uses the Mathieu solution as the reference truth for the comparison.
  • domain assumption Neglect of cosmic expansion, backreaction, and thermal scattering over the exponential-growth phase.
    Stated in Sections II and IV; the central equivalence is derived in their absence, so the Boltzmann match does not extend to realistic early-universe conditions unless (40) and (43) hold.
  • domain assumption Nonrelativistic phi condensate with occupation f_i >> 1 and replacement of the phase-space integral by (n_phi / 2 m_phi)^j.
    Used in Eq. (28) to convert the Boltzmann collision term into the compact form (26).
  • domain assumption The j = 1 decay channel dominates; higher-j channels are suppressed by (q/4)^{2j} for q << 1.
    Used to reduce Eq. (26) to Eq. (27); valid only in the narrow-resonance regime.
  • ad hoc to paper The Dirac delta may be simulated by a Gaussian with coefficients a and b fixed by normalization.
    Eqs. (29)-(34); no first-principles derivation of the Gaussian line shape is given. This is the main approximation controlling the agreement.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Particle decay as asymptotic narrow parametric resonance." pith.science (2026). https://pith.science/paper/OXBZXH7G

@misc{pith2026250503160,
  author       = {Pith},
  title        = {Pith review of: Particle decay as asymptotic narrow parametric resonance},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/OXBZXH7G}},
  note         = {Machine review of arXiv:2505.03160}
}
read the original abstract

Parametric resonance can produce particles from oscillating scalar field, where an exponential growth of the particle number density can be developed. While it has been noticed that stimulated decay in Boltzmann equations exhibits similar parametric dependence of the exponential growth, it is not quantitatively clear yet under what circumstance can the two phenomena reconcile. We demonstrate that trilinear particle interaction in the Boltzmann equation can provide a good approximation to describe the distribution function and the number density in the asymptotic regime of narrow parametric resonance. We find that the crucial treatment leading to the quantitative agreement is a proper Gaussian simulation of the momentum spread inherited from nonrelativistic particle decay. With the simple particle picture, the analytic Boltzmann approximation can be applied to explosive photon production from axions/axion-like particles and to dark matter production from oscillating fields.

Figures

Figures reproduced from arXiv: 2505.03160 by the authors.

Figure 1
Figure 1. FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. The time evolution of the ratios [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. The evolution of [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

55 extracted references · 12 canonical work pages

  1. [23]

    Carenza, A

    P. Carenza, A. Mirizzi, and G. Sigl, Dynamical evolution of axion condensates under stimulated decays into photons, Phys. Rev. D 101 (2020), no. 10 103016, [arXiv:1911.07838]

    arXiv 2020

  2. [24]

    Moroi and W

    T. Moroi and W. Yin, Particle Production from Oscillating Scalar Field and Consistency of Boltzmann Equation, JHEP 03 (2021) 296, [arXiv:2011.12285]

    arXiv 2021

  3. [1]

    Kofman, A

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

    arXiv 1994

  4. [2]

    Kofman, A

    L. Kofman, A. D. Linde, and A. A. Starobinsky, Towards the theory of reheating after inflation, Phys. Rev. D 56 (1997) 3258–3295, [hep-ph/9704452]

    arXiv 1997

  5. [3]

    Agrawal, G

    P. Agrawal, G. Marques-Tavares, and W. Xue,Opening up the QCD axion window, JHEP 03 (2018) 049, [arXiv:1708.05008]

    arXiv 2018

  6. [4]

    R. T. Co, L. J. Hall, and K. Harigaya, QCD Axion Dark Matter with a Small Decay Constant, Phys. Rev. Lett. 120 (2018), no. 21 211602, [arXiv:1711.10486]

    arXiv 2018

  7. [5]

    J. A. Dror, K. Harigaya, and V . Narayan,Parametric Resonance Production of Ultralight Vector Dark Matter, Phys. Rev. D 99 (2019), no. 3 035036, [arXiv:1810.07195]

    arXiv 2019

  8. [6]

    I. I. Tkachev, Fast Radio Bursts and Axion Miniclusters, JETP Lett. 101 (2015), no. 1 1–6, [arXiv:1411.3900]

    arXiv 2015

Show all 55 references
  1. [7]

    M. P. Hertzberg and E. D. Schiappacasse, Dark Matter Axion Clump Resonance of Photons, JCAP 11 (2018) 004, [arXiv:1805.00430]

    2018 arXiv

  2. [8]

    Arza, Photon enhancement in a homogeneous axion dark matter background, Eur

    A. Arza, Photon enhancement in a homogeneous axion dark matter background, Eur. Phys. J. C 79 (2019), no. 3 250, [arXiv:1810.03722]

    2019 arXiv

  3. [9]

    Suyama, T

    T. Suyama, T. Tanaka, B. Bassett, and H. Kudoh, Black hole production in tachyonic preheating, JCAP 04 (2006) 001, 8 [hep-ph/0601108]

    2006 arXiv

  4. [10]

    Y .-F. Cai, X. Tong, D.-G. Wang, and S.-F. Yan,Primordial Black Holes from Sound Speed Resonance during Inflation, Phys. Rev. Lett. 121 (2018), no. 8 081306, [arXiv:1805.03639]

    2018 arXiv

  5. [11]

    Kawai and J

    S. Kawai and J. Kim, Gauss–Bonnet Chern–Simons gravitational wave leptogenesis, Phys. Lett. B 789 (2019) 145–149, [arXiv:1702.07689]

    2019 arXiv

  6. [12]

    Y .-F. Cai, C. Lin, B. Wang, and S.-F. Yan,Sound speed resonance of the stochastic gravitational wave background, Phys. Rev. Lett. 126 (2021), no. 7 071303, [arXiv:2009.09833]

    2021 arXiv

  7. [13]

    Fonseca, E

    N. Fonseca, E. Morgante, R. Sato, and G. Servant, Axion fragmentation, JHEP 04 (2020) 010, [arXiv:1911.08472]

    2020 arXiv

  8. [14]

    R. T. Co, L. J. Hall, K. Harigaya, K. A. Olive, and S. Verner, Axion Kinetic Misalignment and Parametric Resonance from Inflation, JCAP 08 (2020) 036, [arXiv:2004.00629]

    2020 arXiv

  9. [15]

    Kasuya and M

    S. Kasuya and M. Kawasaki, Restriction to parametric resonant decay after inflation, Phys. Lett. B 388 (1996) 686–691, [hep-ph/9603317]

    1996 arXiv

  10. [16]

    Shtanov, J

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

    1995 arXiv

  11. [17]

    Yoshimura, Catastrophic particle production under periodic perturbation, Prog

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

    1995 arXiv

  12. [18]

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

    2014 arXiv

  13. [19]

    R. R. Kovacic, I. and S. Mohamed Sah, Mathieu’s Equation and Its Generalizations: Overview of Stability Charts and Their Features, Appl. Mech. Rev. 70 (2018) 020802

    2018

  14. [20]

    Matsumoto and T

    S. Matsumoto and T. Moroi, Decay of scalar condensation in quantum field theory, Phys. Rev. D 77 (2008) 045014, [arXiv:0709.4338]

    2008 arXiv

  15. [21]

    Yoshimura, Decay rate of coherent field oscillation, in Frontiers in Quantum Field Theory in Honor of the 60th Birthday of Prof

    M. Yoshimura, Decay rate of coherent field oscillation, in Frontiers in Quantum Field Theory in Honor of the 60th Birthday of Prof. K. Kikkawa, pp. 394–400, 3, 1996. hep-ph/9603356

    1996 arXiv

  16. [22]

    Alonso-Álvarez, R

    G. Alonso-Álvarez, R. S. Gupta, J. Jaeckel, and M. Spannowsky, On the Wondrous Stability of ALP Dark Matter, JCAP 03 (2020) 052, [arXiv:1911.07885]

    2020 arXiv

  17. [25]

    Fujisaki, K

    H. Fujisaki, K. Kumekawa, M. Yamaguchi, and M. Yoshimura, Particle production and dissipative cosmic field, Phys. Rev. D 53 (1996) 6805–6812, [hep-ph/9508378]

    1996 arXiv

  18. [26]

    Fujisaki, K

    H. Fujisaki, K. Kumekawa, M. Yamaguchi, and M. Yoshimura, Particle production and gravitino abundance after inflation, Phys. Rev. D 54 (1996) 2494–2503, [hep-ph/9511381]

    1996 arXiv

  19. [27]

    S. Y . Khlebnikov and I. I. Tkachev,The Universe after inflation: The Wide resonance case, Phys. Lett. B 390 (1997) 80–86, [hep-ph/9608458]

    1997 arXiv

  20. [28]

    J. F. Dufaux, G. N. Felder, L. Kofman, M. Peloso, and D. Podolsky, Preheating with trilinear interactions: Tachyonic resonance, JCAP 07 (2006) 006, [hep-ph/0602144]

    2006 arXiv

  21. [29]

    Parker, Quantized fields and particle creation in expanding universes

    L. Parker, Quantized fields and particle creation in expanding universes. 1., Phys. Rev. 183 (1969) 1057–1068

    1969

  22. [30]

    Asaka and H

    T. Asaka and H. Nagao, Non-perturbative Corrections to Particle Production from Coherent Oscillation, Prog. Theor. Phys. 124 (2010) 293–314, [arXiv:1004.2125]

    2010 arXiv

  23. [31]

    Kaneta, S

    K. Kaneta, S. M. Lee, and K.-y. Oda, Boltzmann or Bogoliubov? Approaches compared in gravitational particle production, JCAP 09 (2022) 018, [arXiv:2206.10929]

    2022 arXiv

  24. [32]

    Li and K.-P

    S.-P. Li and K.-P. Xie, Photon proliferation from N-body dark matter annihilation, arXiv:2412.15749

  25. [33]

    Hahn, Generating Feynman diagrams and amplitudes with FeynArts 3, Comput

    T. Hahn, Generating Feynman diagrams and amplitudes with FeynArts 3, Comput. Phys. Commun. 140 (2001) 418–431, [hep-ph/0012260]

    2001 arXiv

  26. [34]

    Shtabovenko, R

    V . Shtabovenko, R. Mertig, and F. Orellana,FeynCalc 10: Do multiloop integrals dream of computer codes?, Comput. Phys. Commun. 306 (2025) 109357, [arXiv:2312.14089]

    2025 arXiv

  27. [35]

    A. Arza, T. Schwetz, and E. Todarello, How to suppress exponential growth—on the parametric resonance of photons in an axion background, JCAP 10 (2020) 013, [arXiv:2004.01669]

    2020 arXiv

  28. [36]

    Jaeckel and S

    J. Jaeckel and S. Schenk, Challenging the Stability of Light Millicharged Dark Matter, Phys. Rev. D 103 (2021), no. 10 103523, [arXiv:2102.08394]

    2021 arXiv

  29. [37]

    D. J. E. Marsh, Axion Cosmology, Phys. Rept. 643 (2016) 1–79, [arXiv:1510.07633]

    2016 arXiv

  30. [38]

    Di Luzio, M

    L. Di Luzio, M. Giannotti, E. Nardi, and L. Visinelli, The landscape of QCD axion models, Phys. Rept. 870 (2020) 1–117, [arXiv:2003.01100]

    2020 arXiv

  31. [39]

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

    1983

  32. [40]

    Dine and W

    M. Dine and W. Fischler, The Not So Harmless Axion, Phys. Lett. B 120 (1983) 137–141

    1983

  33. [41]

    Preskill, M

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

    1983

  34. [42]

    Lee and K.-W

    D.-S. Lee and K.-W. Ng, Photon production of axionic cold dark matter, Phys. Rev. D 61 (2000) 085003, [hep-ph/9909282]

    2000 arXiv

  35. [43]

    Yokoyama, Fate of oscillating scalar fields in the thermal bath and their cosmological implications, Phys

    J. Yokoyama, Fate of oscillating scalar fields in the thermal bath and their cosmological implications, Phys. Rev. D 70 (2004) 103511, [hep-ph/0406072]

    2004 arXiv

  36. [44]

    Yokoyama, Can oscillating scalar fields decay into particles with a large thermal mass?, Phys

    J. Yokoyama, Can oscillating scalar fields decay into particles with a large thermal mass?, Phys. Lett. B 635 (2006) 66–71, [hep-ph/0510091]

    2006 arXiv

  37. [45]

    R. T. Co, L. J. Hall, and K. Harigaya, Axion Kinetic Misalignment Mechanism, Phys. Rev. Lett. 124 (2020), no. 25 251802, [arXiv:1910.14152]

    2020 arXiv

  38. [46]

    R. T. Co and K. Harigaya, Axiogenesis, Phys. Rev. Lett. 124 (2020), no. 11 111602, [arXiv:1910.02080]

    2020 arXiv

  39. [47]

    Di Luzio, B

    L. Di Luzio, B. Gavela, P. Quilez, and A. Ringwald, Dark matter from an even lighter QCD axion: trapped misalignment, JCAP 10 (2021) 001, [arXiv:2102.01082]

    2021 arXiv

  40. [48]

    Di Luzio, B

    L. Di Luzio, B. Gavela, P. Quilez, and A. Ringwald, An even lighter QCD axion, JHEP 05 (2021) 184, [arXiv:2102.00012]

    2021 arXiv

  41. [49]

    Daido, N

    R. Daido, N. Kitajima, and F. Takahashi, Level crossing between the QCD axion and an axionlike particle, Phys. Rev. D 93 (2016), no. 7 075027, [arXiv:1510.06675]

    2016 arXiv

  42. [50]

    Daido, N

    R. Daido, N. Kitajima, and F. Takahashi, Domain Wall Formation from Level Crossing in the Axiverse, Phys. Rev. D 92 (2015), no. 6 063512, [arXiv:1505.07670]

    2015 arXiv

  43. [51]

    S. Y . Khlebnikov and I. I. Tkachev,Classical decay of inflaton, Phys. Rev. Lett. 77 (1996) 219–222, [hep-ph/9603378]

    1996 arXiv

  44. [52]

    S. Y . Khlebnikov and I. I. Tkachev,Resonant decay of Bose condensates, Phys. Rev. Lett. 79 (1997) 1607–1610, [hep-ph/9610477]. 9

    1997 arXiv

  45. [53]

    G. B. Gelmini and M. Roncadelli, Left-Handed Neutrino Mass Scale and Spontaneously Broken Lepton Number, Phys. Lett. B 99 (1981) 411–415

    1981

  46. [54]

    Pitrou, A

    C. Pitrou, A. Coc, J.-P. Uzan, and E. Vangioni, Precision big bang nucleosynthesis with improved Helium-4 predictions, Phys. Rept. 754 (2018) 1–66, [arXiv:1801.08023]

    2018 arXiv

  47. [55]

    Aghanim et al., Planck 2018 results

    Planck Collaboration, N. Aghanim et al., Planck 2018 results. VI. Cosmological parameters, Astron. Astrophys. 641 (2020) A6, [arXiv:1807.06209]. [Erratum: Astron.Astrophys. 652, C4 (2021)]

    2020 arXiv

Pith tools

Reviewed August 15, 2026 · model on record in the stance chip above.