{"id":"f102d6d7-4a1c-4336-b317-e3ec681c4ed6","arxiv_id":"2504.20943","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Using the WI2easy code, the Standard Model gluon-coupled warm inflation model remains viable, with a Hubble-exit dissipation ratio Q* between 0.0076 and 30 and inflaton thermalization for Q* above about 0.08.","lead":"A model where the early universe's inflation field interacts with ordinary Standard Model particles was rechecked with sharper numerical calculations. The model still works, but the allowed parameter range and observable predictions differ from the original analysis.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The viability claim is controlled by Eq. (2.12), the effective dissipation coefficient adopted from Ref. [1] with Ntilde_f=5; this input is neither derived nor validated here, so the precision results inherit its uncertainty.","rationale":"The reader's weakest_assumption correctly identifies Eq. (2.12) as the load-bearing input: it is neither derived nor re-examined in this paper, and the paper's own text flags that fermion effects can strongly suppress or even eliminate sphaleron dissipation. I agree with this assessment. Other potential concerns—such as the use of the authors' own WI2easy code without an independent cross-check, or the approximate thermalization estimate—are secondary because they would only matter if the microphysical dissipation coefficient is reliable. The code is publicly available and the comparison with the old fit is a useful check, but the central physical input remains untested. The paper does not overclaim: it explicitly states the compatibility with Ref. [1] and the regions of deviation. Therefore, the CONDITIONAL verdict is appropriate: the claim is plausible but contingent on the validity of Eq. (2.12). My proposed test directly targets this contingency by comparing Eq. (2.12) with the explicit fermion-suppressed calculation and, if needed, recomputing the observables. No change to the reader's verdict is required; the same conditionality is retained.","tokens_in":12712,"tokens_out":10556,"duration_ms":113132,"concrete_test":"Compute the fermion-suppressed effective dissipation coefficient for the SM at T ∈ [6×10^-7, 5×10^-3] M_Pl using the explicit formulas of Ref. [29] (Drewes & Zell) with N_f = 6, including the top-quark Yukawa chirality-flip rate; compare it against Eq. (2.12) used with Ntilde_f = 5. If the ratio deviates from 1 by more than ~20% over the region giving Q* ∈ [0.0076, 30], re-run the WI2easy background and perturbation code with the corrected U_eff and check whether the allowed Q* interval and fa range survive within the Planck/BICEP/Keck bounds.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central result—the observationally allowed window Q* ∈ [0.0076, 30] and the corresponding fa range—is computed entirely through the effective dissipation coefficient Eq. (2.12), Υeff = Υsph / (1 + 4 Ntilde_f fa^2 Υsph / (H N_c T^2)) with Ntilde_f = 5. This formula is taken from Ref. [1] without derivation or independent re-examination. The paper itself notes (Sec. II, around Eqs. (2.11)–(2.12)) that fermions suppress sphaleron dissipation and 'can even make it completely vanish for massless fermions' [29,37,38]. At the temperatures of interest, T ∈ [6×10^-7, 5×10^-3] M_Pl, all six SM quark flavors are relativistic (m_q/T ≪ 1), so the active flavor content is N_f = 6, not 5. If the true suppression is stronger than Eq. (2.12) assumes—e.g., Ntilde_f should be 6, or the massless-fermion limit alters the functional form—then the mapping between fa, Q*, and the observables (ns, r, T/H) shifts, and the viability window could shrink or disappear. Because WI2easy is a precision tool but its input microphysics is an unverified cited formula, the headline robustness claim is only as strong as Eq. (2.12).","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper revisits the warm-inflation model of Berghaus, Drewes, and Zell (Ref. [1]), in which an axion-like inflaton couples to Standard Model gluons through the interaction φ G μν G̃ μν, and re-evaluates its viability for a quartic potential using the WI2easy code. The main quantitative claim is that for the dissipation ratio at Hubble exit Q* in the interval [0.0076, 30], the spectral tilt n_s and tensor-to-scalar ratio r remain within Planck/BICEP/Keck bounds, which translates into an axion decay constant f_a ∈ [4.98×10^10, 1.23×10^13] GeV. The paper also demonstrates that the approximate power-spectrum fit of Mirbabayi and Gruzinov used in Ref. [1] fails for Q* ≲ 0.6 and Q* ≳ 15, and that inflaton perturbations thermalize for Q* ≳ 0.08, based on a hard-thermal-loop scattering rate. The analysis is framed as a precision re-evaluation of Ref. [1] using a dedicated numerical code.","tokens_in":12999,"tokens_out":11572,"duration_ms":108004,"significance":"If the results hold, the paper provides a valuable precision confirmation of the minimal SM warm-inflation scenario and extends Ref. [1] by mapping the full observationally allowed parameter space, identifying the regime where the approximate fitting formula breaks down, and estimating the thermalization threshold. The work is reproducible in principle: the WI2easy code is public and the dataset is deposited on Zenodo. The central falsifiable prediction is a specific window of axion decay constants (≈10^10–10^13 GeV) and dissipation ratios for quartic-potential warm inflation. However, the quantitative results are controlled by the effective dissipation coefficient of Eq. (2.12), which is adopted from Ref. [1] without independent derivation, and the numerical code is authored by the present authors without an independent cross-check; these facts place a significant burden of robustness on the central viability claim.","major_comments":[{"comment":"The central quantitative results, including the allowed interval Q* ∈ [0.0076, 30] and the corresponding f_a range, are controlled by the effective dissipation coefficient Υ_eff of Eq. (2.12), which is adopted from Ref. [1] with Ñ_f = 5 without derivation or re-examination in this paper. At the temperatures reported in Fig. 5(c) (T ∈ [6×10^-7, 5×10^-3] M_Pl), all six SM quark flavors have m_q ≪ T, so the choice Ñ_f = 5 is not obvious; the paper's own text after Eq. (2.11) notes that fermions 'can even make it completely vanish for massless fermions.' Because this input controls the mapping between f_a, Q*, and the observables, the authors should either derive Eq. (2.12) and the value Ñ_f = 5 from the microphysics, or demonstrate that the viability window is robust to the uncertainty in this formula (for example, by varying Ñ_f between 4 and 6 and by considering the massless-fermion limit). The caveat in footnote 3 about UV completions further indicates that this is an input assumption rather than a tested prediction.","section":"Sec. II, Eq. (2.12)"},{"comment":"The paper presents point predictions for n_s, r, and f_a without any estimate of theoretical uncertainty, and the conclusion that the model is viable relies on the assumption that inflaton perturbations are thermalized for Q* ≳ 0.08, as inferred from a single hard-thermal-loop scattering rate (Eq. (2.17)). Since the axion thermalization literature is cited as debated (Refs. [43,44]), the threshold Q* ≈ 0.08 could shift, and with it the appropriate power-spectrum statistics in the weak dissipative regime. I ask for (i) a sensitivity test of the allowed Q* interval to the thermalization rate, for example using the estimates of Refs. [43,44] as an alternative, and (ii) an explicit statement of which curve (n* = 0, n* = n_BE, or a piecewise combination with threshold Q* = 0.08) is used to derive the quoted Q* ∈ [0.0076, 30] interval.","section":"Sec. III, Figs. 3–4"},{"comment":"The claim that the Mirbabayi-Gruzinov approximation (Eq. (2.15)) becomes unreliable for Q* ≲ 0.6 and Q* ≳ 15 is based entirely on the G(Q) function computed with WI2easy, a code authored by the present authors (Ref. [2]). To rule out that the observed deviations from Ref. [1] are artifacts of the code, the paper should validate G(Q) against an independent implementation, for example WarmSPy (Ref. [30]), or against an analytic limit for a simple dissipation coefficient such as Υ ∝ T^3. This is especially important because the deviations from Ref. [1] are the paper's main new quantitative content.","section":"Sec. III, Fig. 2 and G(Q)"}],"minor_comments":[{"comment":"The text says 'consistent with present bonds on the QCD axion'; 'bonds' should be 'bounds'.","section":"Sec. III, near Fig. 5"},{"comment":"The relationship between the expression for κ in Eq. (2.6) and the explicit Lambert-function result in Eq. (2.9) is not stated; Eq. (2.8) is self-referential, and the way Eq. (2.9) solves it is not shown.","section":"Sec. II, Eqs. (2.6)–(2.9)"},{"comment":"The legends and text do not specify which value of g* (106.75 or 106.75+1) is used for the thermalized-inflaton case; this should be stated explicitly.","section":"Sec. III, Figs. 4–5"},{"comment":"A table reporting n_s, r, T/H, and f_a at representative values of Q* for both n* = 0 and n* = n_BE would increase clarity and reproducibility beyond the figures.","section":"Sec. III"}],"recommendation":"major_revision","confidential_remarks":"The manuscript relies heavily on companion work by the same group (WI2easy, Ref. [2], and the ACT warm-inflation paper, Ref. [49]) and on Ref. [1] for the dissipation coefficient. This is not improper, but the lack of an independent check of the central dissipation input and of the numerical code should be weighed when assessing the robustness of the 'viability' claim. The paper is within the scope of a hep-ph phenomenology journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is a numerical reanalysis of the Berghaus-Drewes-Zell warm inflation model, not a new model. The paper's own contribution is a precision scan with WI2easy: it revises the allowed dissipation ratio at Hubble exit to Q* in [0.0076, 30], maps that to fa in [4.98e10, 1.23e13] GeV, argues inflaton perturbations thermalize for Q* above roughly 0.08, and shows that the analytic G(Q) fit used in the original paper fails for Q* below about 0.6 and above about 15. These are concrete, honest, reproducible outputs. The code and data are public, and the authors are explicit about what is inherited and what is new. Credit where due: they do not oversell.\n\nThe soft spot is exactly what the reader flagged: Eq. (2.12), the effective dissipation coefficient, is taken from Ref. [1] with Ntilde_f = 5 and not re-derived. Everything downstream--Q*, ns, r, fa--scales with that input. The paper notes fermion suppression can shut off sphaleron heating for massless fermions, and at the temperatures of interest all six SM quarks are relativistic, so Ntilde_f = 5 looks like a number worth interrogating rather than a parameter to adopt. It could be that the original authors have a subtle reason for five, but this paper does not supply it. That means the precision numbers carry the uncertainty of the cited formula. Also, the thermalization threshold comes from a single hard thermal loop rate, and the predictions have no error bars. None of this kills the paper; it does mean the headline viability window is conditional.\n\nOn the citation pattern: self-citations are present but not abusive. WI2easy is this group's tool, so citing it is natural. The companion ACT paper is self-cited, but the discussion is measured and actually notes it goes in the opposite direction here.\n\nWho is this for? Warm inflation practitioners who want a careful update on this specific SM-based model. A general cosmologist can skip it. It deserves peer review: a serious referee can pressure-test Eq. (2.12) and ask for a sensitivity study on Ntilde_f. If I were the editor, I would send it out.\n\nRegards.","headline":"A transparent precision reanalysis of a known warm inflation model; the central viability window is plausible but inherits a cited, unreviewed dissipation formula.","tokens_in":13559,"tokens_out":2602,"would_cite":true,"duration_ms":25962,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["98.80.Cq"],"model":"deepseek-v4-flash","headline":"Precision re-analysis keeps warm inflation built on Standard Model gluons inside Planck and BICEP/Keck bounds for $Q_*$ from 0.0076 to 30, with inflaton perturbations already thermal above $Q_*\\simeq 0.08$.","keywords":["warm inflation","standard model","axionlike inflaton","sphaleron dissipation","dissipation coefficient","primordial power spectrum","spectral tilt","tensor-to-scalar ratio"],"falsifier":"A first-principles computation of the sphaleron heating rate in the Standard Model quark-gluon plasma at the temperatures this model lives in ($T \\approx 10^{12}{-}10^{16}$ GeV) would settle the matter: if the five-flavor suppression and the Hubble-dependent correction of Eq. (2.12) are not reproduced, the window $Q_* \\in [0.0076, 30]$ and the associated axion decay constants are not real. On the observational side, the model traces a specific curve in the ($n_s$, $r$) plane, so a future CMB measurement that lands off that curve — a tensor-to-scalar ratio above 0.036, or a shift in $n_s$ beyond current error bars — would exclude it.","tokens_in":12443,"feed_emoji":"🌌","tokens_out":16845,"duration_ms":140443,"temperature":0.7,"pith_summary":"This paper asks whether the most minimal warm-inflation model on the market — an axionlike inflaton coupled to Standard Model gluons through $\\phi G\\tilde G$, with no other beyond-Standard-Model fields — still works when its dynamics are computed precisely instead of with fitting formulas. Using the WI2easy solver, it finds that for the quartic potential $V = \\lambda\\phi^4$ the model matches the Planck and BICEP/Keck bounds on the spectral tilt $n_s$ and the tensor-to-scalar ratio $r$ for dissipation ratios at horizon exit $Q_*$ between 0.0076 and 30, which corresponds to axion decay constants $f_a$ from about $4.98\\times 10^{10}$ to $1.23\\times 10^{13}$ GeV. It also finds that inflaton perturbations have already thermalized for $Q_* \\gtrsim 0.08$, and that the analytic power-spectrum fit used in the original analysis is reliable only for $Q_*$ roughly between 0.6 and 15. The payoff is that the most minimal warm-inflation scenario remains observationally alive, and its viable region is now mapped by a precision numerical tool instead of by approximations.","feed_headline":"Warm inflation built on Standard Model gluons still fits CMB data","feed_subtitle":"A precision rerun finds the quartic potential matches Planck and BICEP/Keck for Q* from 0.0076 to 30.","key_machinery":"The central object is the effective dissipation coefficient of Eq. (2.12), $\\Upsilon_{\\rm eff} = \\Upsilon_{\\rm sph}/(1 + 4\\tilde N_f f_a^2 \\Upsilon_{\\rm sph}/(H N_c T^2))$ with $\\tilde N_f = 5$, which controls how strongly the gluon bath drains energy from the inflaton: the sphaleron heating rate $\\Upsilon_{\\rm sph} = \\Gamma_{\\rm sphal}/(2T f_a^2)$ (built from the SU(3) sphaleron rate of Eq. (2.5) and the running strong coupling $\\alpha_g(T)$) is suppressed both by quark back-reaction and by the Hubble expansion. The second piece is WI2easy, a code that replaces stochastic sampling of warm-inflation perturbations with a deterministic Fokker–Planck formalism and returns the numerical factor $G(Q)$ multiplying the scalar power spectrum of Eq. (2.14); the paper compares this $G(Q)$ against the closed-form fit $F(Q)$ of Ref. [41]. The third is the thermalization check: the hard-thermal-loop scattering rate $\\Gamma_{\\rm scat}$ of Eq. (2.17) is compared with $H$, giving the threshold $Q_* \\approx 0.08$ above which the statistical factor $n_* = n_{\\rm BE}$ must be used.","core_discovery":"Re-evaluating the minimal Standard Model warm-inflation model of Ref. [1] with the WI2easy code, the paper finds the model observationally viable: for the quartic potential $V = \\lambda\\phi^4$, every dissipation ratio at Hubble exit in the interval $Q_* \\in [0.0076, 30]$ keeps the spectral tilt $n_s$ and the tensor-to-scalar ratio $r$ inside the Planck/BICEP/Keck bounds, which pins the axion decay constant to $f_a \\in [4.98\\times 10^{10}, 1.23\\times 10^{13}]$ GeV. Outside that window the model fails sharply: $Q_* < 0.0076$ drives $r$ above the bound $r < 0.036$, and $Q_* \\gtrsim 30$ pushes $n_s$ out of the two-$\\sigma$ region. The paper further claims that the inflaton-perturbation distribution reaches thermal (Bose–Einstein) form for $Q_* \\gtrsim 0.08$, and that the analytic power-spectrum fit used in Ref. [1] agrees with the full numerical result only for $Q_*$ between about 0.6 and 15; below 0.6 the fit is unreliable because it assumed vacuum inflaton perturbations, and above 15 because the temperature power of $\\Upsilon_{\\rm eff}$ drops below three, making $n_s$ redder than the fit predicted.","pith_inferences":["Editorial inference: if Eq. (2.12) is accurate, the model effectively predicts that any axionlike particle responsible for inflation must have a decay constant near $10^{11}{-}10^{13}$ GeV, a band that axion dark matter experiments are beginning to touch; the paper itself only notes compatibility with existing bounds.","Editorial inference: the thermalization threshold $Q_* \\approx 0.08$ relies on a hard-thermal-loop estimate of axion–gluon scattering; a full thermal QCD treatment of axion production at $T \\approx 10^{12}{-}10^{16}$ GeV could move this boundary, and with it the region in which the simpler vacuum-fluctuation formula applies.","Editorial inference: the demonstration that the Ref. [41] power-spectrum fit fails outside $Q_* \\in [0.6, 15]$ is a warning for other warm-inflation models that used the same fit; rerunning them with WI2easy-style numerics could change their claimed parameter windows as well.","Editorial inference: the authors note that ACT's newer preference for slightly larger $n_s$ pushes against this model, in the opposite direction to other warm-inflation models; if that trend is confirmed, the SM-only construction could lose its current two-sigma acceptance even though Planck/BICEP/Keck alone allow it."],"forward_implications":["If the paper is right, the minimal Standard-Model warm-inflation model is observationally viable for the quartic potential, with a dissipation window $Q_* \\in [0.0076, 30]$ that spans more than three orders of magnitude.","The axion decay constant inside the viable window is pinned to $f_a \\in [4.98\\times 10^{10}, 1.23\\times 10^{13}]$ GeV, a range that overlaps existing bounds on the QCD axion.","The old analytic fit to the warm-inflation power spectrum is reliable only for $Q_*$ roughly between 0.6 and 15; analyses that use it outside this range need the full numerical $G(Q)$, and above $Q_* \\approx 15$ the temperature dependence of the dissipation coefficient (effective power $c < 3$) makes the spectrum redder than the fit claimed.","Because inflaton perturbations thermalize already at $Q_* \\gtrsim 0.08$, the non-thermalized version of the power spectrum applies only to a sliver of the viable window; the rest must be evaluated with a Bose–Einstein statistical factor $n_* = n_{\\rm BE}$.","The failure modes are sharp: $Q_* < 0.0076$ violates the bound $r < 0.036$, and $Q_* \\gtrsim 30$ pushes $n_s$ outside the two-sigma region of the Planck/BICEP/Keck data."],"supporting_citations":[{"why":"The model under test: an axionlike inflaton coupled to gluons via phi G Gtilde, including the effective dissipation coefficient with Nf = 5 that the paper takes over.","marker":"[1]"},{"why":"WI2easy, the precision Fokker-Planck solver whose treatment of background and perturbations produces the paper's new results.","marker":"[2]"},{"why":"Origin of sphaleron-based warm inflation from the dimension-five axionlike coupling; supplies the relation between the sphaleron rate and the dissipation coefficient.","marker":"[14]"},{"why":"Provides the SU(N) sphaleron rate used to build the dissipation coefficient, including the coefficient kappa and the Debye-mass running.","marker":"[25]"},{"why":"Underlies the fermion suppression of sphaleron heating that enters the denominator of Eq. (2.12).","marker":"[29]"},{"why":"Gives the one-loop temperature running of the strong coupling used in the numerical evaluation.","marker":"[35]"},{"why":"Also cited for the connection between the sphaleron rate and the warm-inflation dissipation coefficient.","marker":"[36]"},{"why":"The analytic power-spectrum fit F(Q) whose validity the paper tests against WI2easy's numerical G(Q).","marker":"[41]"},{"why":"Supplies the hard-thermal-loop axion scattering rate used to set the thermalization threshold Q* about 0.08.","marker":"[45]"},{"why":"Sets the observational bounds (r < 0.036 and the n_s contours) that define the viable Q* window.","marker":"[47]"}],"fun_headline_variants":["SM gluon warm inflation survives precision CMB test","Warm inflation via Standard Model gluons still fits data","Quartic warm inflation passes CMB bounds for Q* 0.0076-30","Precision rerun keeps minimal SM warm inflation viable","Warm inflation with gluons: broad CMB-safe dissipation window"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole viable window sits on top of one formula the paper does not re-derive: the effective dissipation coefficient of Eq. (2.12), with the fermion suppression factor $1 + 4\\tilde N_f f_a^2 \\Upsilon_{\\rm sph}/(H N_c T^2)$ and the effective number of quark flavors $\\tilde N_f = 5$, is taken over from Ref. [1], and the paper's own footnote 3 concedes that heavy UV-completion fermions could open extra decay channels unless their coupling rates are tuned to stay far below the Hubble rate. If that dissipation law is inaccurate, the computed $Q_*$ window, and with it the claimed ranges for $n_s$, $r$, and $f_a$, all shift and the viability could disappear.","fun_headline_variants_meta":{"raw":{"variants":["SM gluon warm inflation survives precision CMB test","Warm inflation via Standard Model gluons still fits data","Quartic warm inflation passes CMB bounds for Q* 0.0076-30","Precision rerun keeps minimal SM warm inflation viable","Warm inflation with gluons: broad CMB-safe dissipation window"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000798,"raw_usage":{"total_tokens":3534,"prompt_tokens":995,"completion_tokens":2539,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":611,"completion_tokens_details":{"reasoning_tokens":2451}},"tokens_in":611,"tokens_out":2539,"duration_ms":18717,"temperature":1.0,"reasoning_tokens":2451,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:15:56.064160+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A first-principles computation of the sphaleron heating rate in the Standard Model quark-gluon plasma at the temperatures this model lives in ($T \\approx 10^{12}{-}10^{16}$ GeV) would settle the matter: if the five-flavor suppression and the Hubble-dependent correction of Eq. (2.12) are not reproduced, the window $Q_* \\in [0.0076, 30]$ and the associated axion decay constants are not real. On the observational side, the model traces a specific curve in the ($n_s$, $r$) plane, so a future CMB measurement that lands off that curve — a tensor-to-scalar ratio above 0.036, or a shift in $n_s$ beyond current error bars — would exclude it.","supporting_citations":[],"review_version":1}