{"id":"b7e99518-8b52-41fa-a08f-90d193ea6070","arxiv_id":"2412.03958","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Neutrino and antineutrino interaction rates with a QED plasma are defined and computed at NLO, with NLO corrections shifting N_eff only at the fourth decimal in the massless-electron limit.","lead":"This paper computes the full next-to-leading-order QED rates for neutrino and antineutrino production, annihilation, and scattering in the hot plasma of the early universe. These are the microscopic inputs needed to pin down the standard-model prediction for the effective number of neutrino species, N_eff, at the fourth decimal place.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The N_eff shift is obtained by inserting NLO Wightman rates into the leading-order Boltzmann collision form with simple Pauli blocking, a step the paper itself says cannot be extended to NLO; the fourth-decimal claim is an estimate, not a controlled prediction.","rationale":"The paper's central constructive result is the definition and NLO evaluation of double-differential rates. That part is internally coherent, reproduces the earlier neutrino interaction rate, and is supported by a public tabulation and interpolation routine, which deserves credit. The weakest point is not the rate calculation but the bridge to N_eff. The reader identified the same issue: the paper explicitly disclaims any direct NLO extension of Boltzmann equations, yet the final numerical N_eff values use the NLO rates inside precisely the leading-order Boltzmann collision structure. This is not a mere presentational slip; it is the step that converts the reliably computed Wightman rates into the headline claim that NLO effects affect the fourth decimal. The paper's own conclusion acknowledges that a full NLO estimate requires solving non-averaged kinetic equations, reinforcing the conditional reading. Electron-mass effects at LO are an order of magnitude larger than the quoted NLO shift, further supporting the need to treat the N_eff statement as an estimate within stated approximations. I therefore keep the reader's CONDITIONAL verdict unchanged and recommend no adjustment to the verdict.","tokens_in":25154,"tokens_out":10649,"duration_ms":116365,"concrete_test":"Construct the full NLO neutrino quantum kinetic equation from the Wightman self-energy underlying eqs. (3.25) and (3.29), keeping the p0 integration off-shell, then perform a gradient expansion and compare the resulting collision integral with the on-shell factorized form of eqs. (4.2)-(4.4). If the neglected off-shell and multi-particle terms change the energy-transfer rate at a relative level above about 10^-4, the Boltzmann insertion assumption controls the fourth decimal of N_eff, and the quoted NLO shift is not robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3 defines Psi, Psi-tilde, and Theta directly from thermal Wightman functions; this part is coherent, gauge-consistent, and validated by reproducing the neutrino interaction rate of ref. [6]. The load-bearing step appears in Secs. 4 and 6. There, the quantum-mechanical rates are inserted as gain and loss terms in the averaged Boltzmann equations of refs. [8,9], with the ordinary factors (1-f)(1-f) and f f. The paper states in Sec. 2 that Boltzmann equations themselves cannot be directly extended to the NLO level, yet eqs. (4.2)-(4.4) assume exactly this factorization: a fully inclusive transition rate computed from an initial state with no neutrino is multiplied by final-state Pauli factors as if the collision term at NLO had the same functional dependence on f_nu. NLO corrections generically involve off-shell propagation, multi-particle channels, and medium-modified self-energies, none of which are shown to reduce to this simple gain-loss structure. The quoted N_eff shift, 3.04859 to 3.04867, is much smaller than the LO electron-mass shift reported in the same paper (3.04859 to 3.04510), and the computation also uses a single averaged neutrino temperature. The central rate computation survives, but the N_eff number should be read as a model-dependent estimate within the averaged kinetic scheme, not as a complete NLO Standard Model prediction.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper defines double-differential production, annihilation, and scattering rates for neutrinos and antineutrinos in a hot QED plasma, deriving them from thermal Wightman functions rather than from a perturbative expansion of Boltzmann equations. The central results are the NLO expressions for the double-differential rates in Eqs. (3.25), (3.26), and (3.29), the energy transfer rates built from them in Eqs. (4.2)-(4.4), and a tabulation/interpolation package for the momentum-dependent rates. The NLO spectral functions are taken from previous two-loop computations, and the integrated neutrino interaction rate is shown to reproduce the earlier result of Ref. [6]. In the massless-electron limit the NLO corrections to energy transfer rates are about one percent, and when inserted into the averaged kinetic scheme of Refs. [8,9] the quoted N_eff changes from 3.04859 (LO) to 3.04867 (NLO). A LO treatment of finite electron mass, in Appendix B, changes N_eff to 3.04510.","tokens_in":25397,"tokens_out":2527,"duration_ms":25992,"significance":"If the double-differential NLO rates are correct, the paper provides a valuable and more complete object than the single integrated interaction rate of Ref. [6]: it supplies the momentum-resolved integrands needed for non-averaged kinetic equations, and it establishes a direct connection between the earlier rate computation and the energy transfer rates used in averaged codes. The cross-check against Ref. [6], the use of previously published NLO spectral functions from Refs. [12,13], and the public tabulation and interpolation routine are concrete strengths. The quoted N_eff shift, however, should be read as an estimate within a specific averaged kinetic scheme, not as a complete NLO Standard Model prediction; the paper itself says in Sec. 2 that Boltzmann equations cannot be directly extended to NLO, while Eqs. (4.2)-(4.4) adopt exactly the LO collision-term factorization.","major_comments":[{"comment":"The paper states in Sec. 2 that \"Boltzmann equations themselves cannot be directly extended to the NLO level\", yet Eqs. (4.2)-(4.4) insert the NLO Wightman-based rates into the leading-order gain-loss collision terms with simple Pauli blocking factors. The step from the quantum-mechanical rate, computed from an initial state with no neutrino population, to a Boltzmann collision term with factorized (1-f)(1-f) or f f dependence is an additional assumption at NLO, not a derived consequence of the QKE. This makes the N_eff numbers in Sec. 6 model-dependent estimates within the averaged kinetic scheme rather than a complete NLO prediction. The central rate computation is not undermined, but the N_eff claim needs to be explicitly labeled as such, or a derivation/justification of the factorization at NLO must be provided.","section":"Sec. 2 and Eqs. (4.2)-(4.4)"},{"comment":"The NLO rates are computed only for massless electrons, while the paper itself shows that the LO electron-mass correction changes N_eff in the third decimal (3.04859 -> 3.04510, Sec. 6). The quoted NLO-to-LO shift (3.04859 -> 3.04867) is therefore a massless-electron result, and the conclusion that \"NLO effects influence the fourth decimal\" is not a complete Standard Model statement once O(me/T) effects are admitted to be larger at LO. This should be stated more prominently in the abstract and conclusions, since otherwise readers may take the N_eff value as the full NLO SM result.","section":"Sec. 4.4 and Sec. 6"},{"comment":"The finite-temperature domain is restricted to T_gamma >= m_e, and the NLO expressions rely on the factorization of the hadronic uncertainty through the coefficient C_a. The uncertainty estimate in Eq. (4.28) is given as fixed percentages, but the overall reliability of the NLO interpolation is tested only against the integrated Q_LO and against the neutrino interaction rate of Ref. [6]. It would be useful to state explicitly whether the <0.05% accuracy quoted in Sec. 5 applies to all integrated quantities over the full T_nu/T_gamma range used in the N_eff estimate, or only to the specific test cases shown.","section":"Eq. (4.27) and Sec. 5"}],"minor_comments":[{"comment":"The footnote says the electron mass could have been kept non-zero, but the paper then proceeds with massless electrons for all NLO results; a brief statement near Eq. (3.13) clarifying that Appendix B only covers the LO mass dependence would help avoid confusion.","section":"Footnote 2, Sec. 3.2"},{"comment":"The loss-term curves labelled \"MB-LO\" are multiplied by 10^{-1} as stated in the caption, but the corresponding legend entry does not show this rescaling; please make the rescaling explicit in the legend or in the axis annotation.","section":"Fig. 2"},{"comment":"The mX/pX notation is explained, but the tables would be easier to read if the power-of-ten exponent were separated from the mantissa, or if the entries were written in scientific notation consistently.","section":"Tables 1 and 2"},{"comment":"The interpolation grid is described for the first quadrant, and the text says the shaded region is obtained by reflecting p0 -> -p0; it would be useful to state explicitly that the tabulated data therefore covers all domains of Fig. 3 used in the integrations.","section":"Sec. 5"},{"comment":"The sentence \"we find N_eff = 3.04858 with the improved Maxwell-Boltzmann approximation\" uses a different number of digits than the later values; please harmonize the significant digits in the quoted N_eff values.","section":"Sec. 6"}],"recommendation":"major_revision","confidential_remarks":"The manuscript's core derivation of the double-differential rates is convincing and the cross-checks are reassuring. The main editorial issue is that the N_eff claim is presented with a precision that is not supported by the NLO collision-term construction: the quoted fourth-decimal shift relies on inserting NLO rates into a LO Boltzmann structure. A revised version that clearly frames the N_eff value as an estimate within the averaged scheme, and that separates the massless-electron NLO statement from the LO electron-mass sensitivity, would bring the claims in line with the actual scope of the calculation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: the double-differential rates are the real result; the N_eff shift is a conditional estimate, and the paper mostly says so itself. The central QFT derivation is in good shape. Jackson and Laine show that the integrand of the previously computed neutrino interaction rate is a physical, gauge-invariant double-differential rate, and they compute production, annihilation, and scattering rates at NLO in QED for massless electrons. They verify the construction by recovering the known interaction rate of ref. [6], and they ship a tabulation and interpolation routine. That is a solid, reproducible piece of work, and it clears up the earlier confusion between the partial NLO estimates in refs. [5] and [7].\n\nThe soft spots are real but mostly located where the authors point. The NLO calculation is done for me=0; the LO electron-mass effect on N_eff is about ten times larger than the NLO shift, so the quoted 3.04867 is not a complete SM prediction. More importantly, the N_eff numbers come from inserting the quantum-mechanical rates into the averaged Boltzmann equations with the usual (1-f)(1-f) and f f factors, while the paper itself notes in Sec. 2 that Boltzmann equations cannot be directly extended to NLO. That step is an assumption, not a derivation from the full quantum kinetic equation. The N_eff shift should therefore be read as an estimate within the averaged single-temperature scheme, not as a controlled NLO prediction. The paper is transparent about this, so I would not call it a fatal flaw, but it could be stated more prominently in the abstract or conclusions.\n\nThe tabulation is a strong asset: the authors make the full momentum-dependent rates available, which is exactly what full kinetic-equation solvers need. The reproduction of the ref. [6] interaction rate and the agreement with ref. [7] on the sign and rough size of the NLO correction are good sanity checks.\n\nWho is this for: anyone working on precision N_eff, neutrino decoupling, or thermal field theory rate computations. It deserves a serious referee. The main referee asks should be whether the Wightman-function definition of the rates is fully rigorous with respect to the Boltzmann-factorization step, and whether the published version makes the approximate nature of the N_eff estimate unmistakable. I would engage with it.","headline":"Strong double-differential NLO rate computation with a transparent but real caveat on the N_eff estimate.","tokens_in":25978,"tokens_out":2254,"would_cite":true,"duration_ms":20638,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.20.-m","95.30.Cq","98.80.-k"],"model":"deepseek-v4-flash","headline":"Neutrino production, annihilation, and scattering rates in a hot QED plasma can be computed at NLO in QED, shifting N_eff only at the fourth decimal.","keywords":["neutrino decoupling","effective number of neutrino species","QED plasma","next-to-leading order","double-differential rates","thermal field theory","energy transfer rates","Boltzmann equations"],"falsifier":"A direct solution of the full quantum kinetic equation for neutrinos at $T_\\gamma\\sim1\\text{--}3\\,{\\rm MeV}$, using the same NLO spectral data but without the factorized Pauli-blocking approximation of eqs. (4.2)–(4.4), would settle the central claim: if the resulting $N_{\\rm eff}$ differs from $3.04867$ already in the third decimal, the factorization assumption fails at NLO.","tokens_in":24893,"feed_emoji":"🌌","tokens_out":13024,"duration_ms":107046,"temperature":0.7,"pith_summary":"This paper shows that the double-differential rates for neutrino–antineutrino pair production, annihilation, and scattering in a QED plasma at MeV temperatures can be defined in full thermal field theory and evaluated at next-to-leading order in QED. These rates are the momentum-resolved building blocks behind the collision terms that determine how neutrinos decouple in the early universe, and hence they feed the Standard Model prediction for the effective neutrino number $N_{\\rm eff}$. In the massless-electron limit the NLO corrections are about one percent, shifting $N_{\\rm eff}$ from $3.04859$ at leading order to $3.04867$ at NLO. The electron mass, treated at leading order, is a larger effect: it lowers $N_{\\rm eff}$ to $3.04510$. The paper also tabulates all double-differential rates with a fast interpolation routine so they can be used in non-averaged kinetic equations.","feed_headline":"NLO QED corrections move N_eff only at the fourth decimal","feed_subtitle":"Precise rates for neutrinos in the MeV plasma shift the decoupling parameter from 3.04859 to 3.04867.","key_machinery":"The object that carries the argument is the function $F(k;p_0,p)$ in eqs. (3.25), (3.26) and (3.29). It packages the dynamical information—transverse and longitudinal photon spectral functions $\\rho_T,\\rho_L$ and their NLO corrections, vertex corrections $\\chi_{T,L}$, and HTL-resummed photon propagators $R^*_{T,L}$—together with flavour and coupling factors, and it determines all three rates because the neutrino momentum $k$ enters only as a kinematic weight. The paper factorizes $F$ into coefficients $A_i,B_i,C_i,D_i$ that are precomputed on a wedge-shaped grid in the $(p_+,p_-)$ plane and linearly interpolated, making the otherwise costly two-loop spectral integrals available inside a kinetic-equation solver. The Wightman-correlator definition enforces detailed balance between production and annihilation, and reduces to the known NLO neutrino interaction rate in the literature when integrated with the appropriate measure.","core_discovery":"The central claim is that the integrand of a previously computed NLO neutrino interaction rate is itself a physical quantity: the double-differential production rate $\\Psi(\\mathbf{k}_\\nu,\\mathbf{q}_{\\bar\\nu})$, the annihilation rate $\\tilde{\\Psi}$, and the scattering rate $\\Theta(\\mathbf{q}_\\nu\\to\\mathbf{k}_\\nu)$ can be defined directly from Wightman correlators of the weak current and computed to $O(\\alpha_{\\rm em})$ for $T_\\gamma \\ge m_e$. Weighted integrals of these rates reproduce the energy transfer rates used in averaged kinetic equations, and in the massless limit the NLO corrections to those transfer rates are of order $\\sim 1\\%$. Inserting the NLO rates into the averaged decoupling equations gives $N_{\\rm eff}=3.04867$, against $3.04859$ at full LO and $3.04858$ in the improved Maxwell–Boltzmann approximation; the finite electron mass at LO gives $3.04510$, a third-decimal effect. The full momentum-dependent rates are tabulated so that future non-averaged kinetic equations can carry the same NLO accuracy.","pith_inferences":["Because the paper itself notes that Boltzmann equations cannot be directly extended to NLO, the largest remaining uncertainty may be the factorization of NLO rates into gain and loss terms with simple Pauli blocking; the tabulated rates are exactly what would be needed to test this in a full quantum kinetic equation.","The same Wightman-correlator construction should apply to other light weakly-interacting species coupled to a thermal electromagnetic current, provided the relevant current correlators and spectral decompositions are known.","The electron-mass effect at LO exceeds the NLO QED shift, so a finite-mass NLO computation is the natural next step; until it is done, the third decimal of the Standard Model $N_{\\rm eff}$ prediction rests on a leading-order mass treatment.","The single-temperature average used for the paper's $N_{\\rm eff}$ numbers may wash out momentum-dependent NLO effects; feeding the tabulated rates into a momentum-resolved solver is a direct way to check whether the fourth-decimal conclusion survives."],"forward_implications":["A full NLO computation of $N_{\\rm eff}$ from momentum-dependent kinetic equations is now in principle possible, since the rates that parametrize those equations have been evaluated at NLO without kinematic approximations.","In the massless-electron limit, NLO QED corrections change $N_{\\rm eff}$ only in the fourth decimal, from $3.04859$ to $3.04867$, placing the third decimal beyond the reach of this particular correction.","A finite electron mass moves the third decimal at leading order, $N_{\\rm eff}=3.04510$, so the mass dependence is a larger systematic than the NLO QED correction.","The released tabulation reproduces the full leading-order energy transfer rates to better than $0.05\\%$ and can also reconstruct the previously computed NLO neutrino interaction rate, giving independent cross-checks.","NLO results agree with one earlier partial estimate and disagree with another on the size of the shift, with the complete set of diagrams keeping the NLO effect small."],"supporting_citations":[{"why":"Supplies the previously computed NLO neutrino interaction rate; this paper shows its integrand is the physical double-differential rate and uses it for validation.","marker":"[6]"},{"why":"A partial NLO computation within the Maxwell-Boltzmann approximation that this paper compares against and disagrees with on the size of the N_eff shift.","marker":"[5]"},{"why":"An estimate of one NLO diagram in the collision integral; the full NLO result here agrees with its conclusion that the effect is at the fourth decimal.","marker":"[7]"},{"why":"Origin of the improved Maxwell-Boltzmann averaged kinetic equations and the analytic energy transfer rates reproduced in eq. (4.16).","marker":"[8]"},{"why":"Provides the averaged kinetic-equation code and the LO electron-mass N_eff value used to judge the impact of the new rates.","marker":"[9]"},{"why":"Supplies the two-loop thermal spectral functions needed for the NLO transverse and longitudinal photon spectral functions.","marker":"[12]"},{"why":"Supplies the numerical master integrals used to evaluate those two-loop spectral functions.","marker":"[13]"},{"why":"Supplies the low-energy hadronic input that sets the quoted uncertainty from the matching coefficient C_a.","marker":"[14]"}],"fun_headline_variants":["NLO QED rates shift N_eff by just 0.00008","N_eff stays at 3.04867 with NLO neutrino rates","Fourth-decimal N_eff change from NLO QED","MeV plasma neutrino rates at NLO: N_eff shift 0.00008","NLO QED nudges N_eff to 3.04867"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that NLO rates computed from equilibrium Wightman correlators can be inserted into Boltzmann equations as simple gain and loss terms multiplied by Pauli blocking factors, even though the paper states that Boltzmann equations themselves cannot be directly extended to NLO.","fun_headline_variants_meta":{"raw":{"variants":["NLO QED rates shift N_eff by just 0.00008","N_eff stays at 3.04867 with NLO neutrino rates","Fourth-decimal N_eff change from NLO QED","MeV plasma neutrino rates at NLO: N_eff shift 0.00008","NLO QED nudges N_eff to 3.04867"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000644,"raw_usage":{"total_tokens":2968,"prompt_tokens":963,"completion_tokens":2005,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":579,"completion_tokens_details":{"reasoning_tokens":1903}},"tokens_in":579,"tokens_out":2005,"duration_ms":13350,"temperature":1.0,"reasoning_tokens":1903,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T21:53:50.331555+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct solution of the full quantum kinetic equation for neutrinos at $T_\\gamma\\sim1\\text{--}3\\,{\\rm MeV}$, using the same NLO spectral data but without the factorized Pauli-blocking approximation of eqs. (4.2)–(4.4), would settle the central claim: if the resulting $N_{\\rm eff}$ differs from $3.04867$ already in the third decimal, the factorization assumption fails at NLO.","supporting_citations":[{"cited_title":"Jackson, Numerical code for master integrals for thermal spectral fu nctions, https://doi.org/10.5281/zenodo.3478143 (2019)","cited_arxiv_id":null,"evidence_quote":"Supplies the numerical master integrals used to evaluate those two-loop spectral functions."}],"review_version":1}