{"id":"f948eb91-8033-487e-a02b-ee4216e6c155","arxiv_id":"2505.03160","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A Boltzmann equation with Gaussian momentum spread reproduces the asymptotic growth rate and number density of narrow parametric resonance to within about 20 percent.","lead":"This paper shows that a simple particle-decay equation, the Boltzmann equation, can reproduce the explosive particle production that happens in narrow parametric resonance, as long as the decay momentum spread is modeled as a Gaussian. The result gives an analytic handle for estimating photon and dark matter production from oscillating fields in the early universe.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"0.8 number-density ratio is set by the Gaussian ansatz, not by the Boltzmann equation: the exact Mathieu band curvature differs and gives sqrt(2/pi)≈0.8; test with the actual Floquet line shape is needed.","rationale":"In good faith, the paper does useful work: it shows that a Gaussian-smoothed Boltzmann equation can reproduce the peak growth rate q m/2 and an order-unity number density, and the FeynArts/FeynCalc cross-check plus the q <= 0.1 numerics are genuine supporting evidence. However, the central quantitative conclusion, the 20% density deviation independent of q, is fixed before any Boltzmann physics enters. Once the peak rate and normalization are imposed, the Gaussian curvature is determined, and the saddle-point ratio to the exact Mathieu band is sqrt(2/pi) ≈ 0.8. The q-independence and the value 0.8 therefore do not validate the Boltzmann equation; they validate that a Gaussian with matching peak and normalization gives a 20% density error relative to the Mathieu solution. This sharpens the reader's weakest assumption: it is not merely that the true line shape might be non-Gaussian; the specific numerical claim 0.8 is exactly the width mismatch between the chosen Gaussian and the actual Floquet profile. A different normalized profile with the same width, such as a Lorentzian or top-hat, gives a different ratio, so the claim is not robust. The cutoff n flagged by the reader is secondary: for qz >> 1 the integrals are saturated well before the cutoff, so the asymptotic ratio is insensitive to n = O(1). No code was shipped, so an independent recomputation is the appropriate test. The paper's applications to delayed oscillations and MeV-scale observables inherit this uncertainty, but the central claim is where the issue lies. I therefore keep the conditional verdict and request the exact-line-shape test.","tokens_in":11411,"tokens_out":12659,"duration_ms":135805,"concrete_test":"Replace the Gaussian delta simulation in Eq. (29) by the exact Mathieu line-shape factor rho(Δ) = sqrt(q^2 − (4Δ/m)^2) for |Δ| ≤ q m/4, normalized so that ∫_0^∞ dk rho = 1/2, and recompute Fig. 2 from Eq. (20). If n_M/n_B moves to 1, or becomes q-dependent, rather than staying at 0.8, the claimed 20% Boltzmann-NPR agreement is a consequence of the Gaussian choice rather than of the Boltzmann equation. A secondary check with a Lorentzian or top-hat profile of the same width would confirm whether the 0.8 ratio is shape-dependent.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The quantitative core of the paper is Fig. 2: n_M/n_B → 0.8, independent of q, claimed as evidence that the Boltzmann equation reproduces narrow parametric resonance. This constant is not a prediction of the Boltzmann collision term; it is fixed by the Gaussian shape chosen in Eq. (29). In the q<<1 first instability band, the Mathieu occupation grows as n_k ∝ exp(qz sqrt(1 − 16 r^2)) with r ≡ (k − m/2)/(q m), so near the peak it behaves as exp(qz(1 − 8 r^2)). The Gaussian ansatz of Eqs. (29)-(36) gives f_chi ∝ exp(qz exp(−16 r^2/π)), i.e. exp(qz(1 − (16/π) r^2)) near the peak. Both are normalized to the same peak growth rate q/2, but their quadratic curvatures differ: 8 versus 16/π ≈ 5.09. For qz >> 1, both integrals are saturated near the peak, and their ratio is sqrt((16/π)/8) = sqrt(2/π) ≈ 0.798, precisely the 0.8 shown in Fig. 2, independent of q and z. Thus the headline 20% statement is the saddle-point width ratio of the two chosen line shapes, not a property of the Boltzmann equation. The width sigma = q m is imported from the Mathieu bandwidth, and the Gaussian coefficient b is then fixed by normalization and the peak rate, leaving no freedom to match the true band profile. If the exact Floquet shape were inserted, the comparison would actually test the Boltzmann equation; as it stands, the central quantitative claim is contingent on the ansatz shape.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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].","tokens_in":11746,"tokens_out":31632,"duration_ms":311900,"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":[{"comment":"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.","section":"Section III, Eqs. (29)–(36) and Fig. 2"},{"comment":"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.","section":"Section III, Eq. (33)"}],"minor_comments":[{"comment":"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.","section":"Section III, Eq. (38) and Fig. 2"},{"comment":"Reference [35] appears in the bibliography but is not cited in the text; please add the citation or remove the entry.","section":"References"},{"comment":"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.","section":"Section II, Eq. (8)"},{"comment":"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_χ'.","section":"Fig. 1 caption"}],"recommendation":"major_revision","confidential_remarks":"The Mathieu-side calculation is sound and the paper is clearly written. The obstacle is that the Boltzmann-side comparison is not a genuine test of the collision term because the Gaussian line shape fixes the 20% constant. I would advise requiring a direct comparison with the exact Floquet profile or a first-principles derivation of the line shape before acceptance; otherwise the abstract's quantitative claim should be substantially softened."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my read. The paper does one genuinely useful thing: from the Boltzmann equation with the j=1 decay channel, it gets the correct narrow parametric resonance growth rate q m/2, with the right prefactor 1/2, not the pi/4 of Ref. [23]. The peak distribution function f_chi = (1/2) e^{G q z} - 1/2 also matches the Mathieu result once the occupation is large. That part is clean and a real step beyond previous attempts.\n\nBut the headline claim—the number density ratio n^M/n^B approaching 0.8 independent of q—does not hold up under closer inspection. The 0.8 is fixed by the Gaussian ansatz in Eq. (29), not by the Boltzmann equation. In the first instability band, the Mathieu occupation near the peak grows as exp(q z (1 - 8 r^2)) with r ≡ (k - m/2)/(q m). The Gaussian ansatz of Eq. (35) behaves as exp(q z (1 - (16/pi) r^2)) near the peak. Both are normalized to the same peak growth rate q/2, but the band curvatures differ: 8 versus 16/pi ≈ 5.09. For q z >> 1 the integrals saturate near the peak, and their ratio is sqrt((16/pi)/8) = sqrt(2/pi) ≈ 0.798, which is exactly the 0.8 in Fig. 2. In other words, the 20% discrepancy is the saddle-point width ratio of the chosen line shape, and it would be different for any other line shape. The width sigma = q m is imported from the Mathieu bandwidth, so the comparison is partly circular.\n\nThere are also smaller issues: the momentum cutoff n in Eq. (38) is never specified, and the numerical check only covers q ≲ 0.1. None of these are fatal for the rate calculation, which survives. But the central quantitative claim of the paper, the q-independent 0.8, is not a prediction of the Boltzmann equation.\n\nThe right fix is to replace the Gaussian with the actual Floquet line shape and see whether the Boltzmann equation still gives the same number density. If it does, the paper would be genuinely important; if not, the rate matching remains a useful but much weaker statement.\n\nThis is a paper for people working on axion/ALP photon production and oscillating-field dark matter, especially those connecting to BBN/CMB Boltzmann codes. I'd send it to peer review—the issue I raised is precise and testable, and the rate derivation deserves referee time. But I would not cite the 20% number-density claim until the shape dependence is resolved.","headline":"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.","tokens_in":12298,"tokens_out":3175,"would_cite":false,"duration_ms":30496,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["narrow parametric resonance","Boltzmann equation","stimulated decay","Gaussian momentum spread","Mathieu equation","scalar field decay","axion photon production","dark matter production"],"falsifier":"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.","tokens_in":11151,"feed_emoji":"⚛️","tokens_out":6815,"duration_ms":54275,"temperature":0.7,"pith_summary":"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.","feed_headline":"Boltzmann decay matches narrow parametric resonance to 20 percent","feed_subtitle":"A Gaussian momentum width lets decay equations reproduce field-theory growth for axion and dark-matter studies.","key_machinery":"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$.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Defines narrow parametric resonance, the $q$ parameter, and the exponential growth rate $\\Gamma\\propto e^{qmt}$ that the paper aims to reproduce.","marker":"[2]"},{"why":"Supplies the Mathieu equation and its stability-band analysis, from which the resonance shape and the momentum spread $\\sigma = q m_\\phi$ are inferred.","marker":"[19]"},{"why":"Earlier Boltzmann treatment that, with a QFT-constructed distribution, overproduced the number density by orders of magnitude; this is the discrepancy the paper resolves.","marker":"[24]"},{"why":"Previous stimulated-decay estimate $\\tilde\\mu = q m_a\\pi/4$ that differs from the NPR rate $q m_a/2$ and motivates the quantitative matching.","marker":"[23]"},{"why":"Provides the N-body amplitude formulas $|M_j|^2$ used to justify keeping only the decay channel $j=1$ for $q\\ll 1$.","marker":"[32]"},{"why":"Earlier use of the uncertainty principle to model Bose-enhanced photon production from axions, giving the qualitative stimulated-decay picture.","marker":"[22]"},{"why":"Establishes that a scalar condensate can be described by the S-matrix decay rate in quantum field theory, grounding the Boltzmann approach.","marker":"[20]"}],"fun_headline_variants":["Boltzmann with Gaussian width hits narrow resonance to 20%","Decay equations reproduce parametric resonance growth","Gaussian momentum spread bridges decay and resonance","Boltzmann approximation captures explosive production phase"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Boltzmann with Gaussian width hits narrow resonance to 20%","Decay equations reproduce parametric resonance growth","Gaussian momentum spread bridges decay and resonance","Boltzmann approximation captures explosive production phase"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000488,"raw_usage":{"total_tokens":2409,"prompt_tokens":958,"completion_tokens":1451,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":574,"completion_tokens_details":{"reasoning_tokens":1406}},"tokens_in":574,"tokens_out":1451,"duration_ms":9390,"temperature":1.0,"reasoning_tokens":1406,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:59:03.240590+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Mathieu equation and its stability-band analysis, from which the resonance shape and the momentum spread $\\sigma = q m_\\phi$ are inferred."},{"cited_title":"Decay of Scalar Condensation in Quantum Field Theory","cited_arxiv_id":"0709.4338","evidence_quote":"Establishes that a scalar condensate can be described by the S-matrix decay rate in quantum field theory, grounding the Boltzmann approach."}],"review_version":1}