{"id":"26da9f41-17fe-4d0a-a280-9bb3887be79c","arxiv_id":"2411.16004","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"State-of-the-art LHC predictions for W/Z transverse momentum spectra at N3LL' and approximate N4LL accuracy, with a rigorous reduction of nonperturbative TMD physics to one effective function per process.","lead":"This paper produces the most precise theoretical predictions to date for how often W and Z bosons are produced with a given sideways momentum at the LHC. It also proves that all nonperturbative quantum effects in these spectra can be compressed into one effective function, which sharpens precision measurements such as the W boson mass.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The ported two-parameter nonsingular model (eq. 2.23) is the load-bearing weak point: the W/fiducial-shape predictions and cumulative-cross-section claim inherit an unvalidated O(alpha_s^3) nonsingular fitted only for Z+jet ATLAS cuts and set to zero below q*_T = 8 GeV.","rationale":"After reading the paper in good faith, I find the main threat to the central claim is exactly the one identified by the reader: the O(alpha_s^3) nonsingular is not a full N3LO calculation but a two-parameter model fitted to Z+jet data at ATLAS cuts, with a hard cutoff at 8 GeV, and then ported to W and other settings. The paper is transparent about this, and the uncertainty bands are carefully built, but the quoted uncertainty on the nonsingular only reflects the fit errors c1,c2, not the process/cut universality of fnons or the zero-integral assumption below q*_T for W. Because the cumulative pT^Z cross-section claim (Sec. 5) and the W spectra (Sec. 6) are central advertised results at state-of-the-art precision, an unvalidated port is load-bearing. Other potential concerns, such as the nonperturbative effective-function construction in Sec. 3 or the approximate N4LL Sudakov ingredients, are either controlled by explicit power-counting statements or shown to have small numerical impact, and the paper's uncertainty decomposition gives them reasonable coverage. I therefore agree with the reader's weakest assumption, and the verdict should remain CONDITIONAL: the predictions are plausible and well-constructed, but the W and cumulative-cross-section claims should be qualified by, or validated with, a dedicated N3LO W+jet nonsingular calculation.","tokens_in":57283,"tokens_out":8504,"duration_ms":83014,"concrete_test":"Using the public N3LO W+jet implementation of Neumann and Campbell (ref. [40], also for Z with CMS cuts if available), with the fiducial cuts of eq. (6.1), compute the O(alpha_s^3) nonsingular coefficient dsigma_nons^(3)/dqT at a technical cutoff qcut_T = 0.447 GeV as in ref. [43]. Then check two numbers: (a) I_W = integral_0^{8 GeV} dqT dsigma_nons^(3)/dqT. If |I_W| is not compatible with zero at the level of the Z value (-0.34 +/- 0.76 pb) rescaled by the relative W/Z cross-section size, the N3LO normalization of the W spectrum is biased. (b) The ratio [dsigma_nons^(3)/dqT]/[dsigma_NLO1/dqT] for 8 < qT < 60 GeV, compared with the fitted fnons from eq. (2.26). If the deviation exceeds the factor-10-scaled uncertainty band used in the paper, the W predictions are overconfident. This single computational campaign settles whether the ported model is the bottleneck.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim of N3LL' + O(alpha_s^3) matched predictions for fiducial W and Z pT spectra requires the O(alpha_s^3) nonsingular contribution to be accurate where it matters. The only true O(alpha_s^3) fixed-order input is Z+jet production with ATLAS 13 TeV cuts from NNLOjet [43]. From it, eq. (2.23) defines a two-parameter model fnons(qT) for the ratio dsigma_nons^(3)/dsigma_NLO1, fitted for qT > q*_T = 8 GeV and set to zero below q*_T. Eq. (2.25) shows the integral of the Z nonsingular below 8 GeV is consistent with zero. This fnons is then taken to be universal: it is applied to W+ production and to CMS Z cuts, with the uncertainty only scaled by a factor 10 for W. No calculation verifies (i) that the sub-q*_T integral vanishes for W or for other fiducial cuts, or (ii) that the ratio fnons has the same shape for W in the 8-20 GeV region, where the nonsingular is a few-percent effect. The cumulative pT^Z cross section, advertised as a percent-level N3LO PDF constraint, inherits any bias in the total normalization from the missing or non-universal low-qT nonsingular. The quoted approximation uncertainty propagates only the c1,c2 fit errors (eq. 2.26); it does not include the systematic uncertainty from porting fnons across processes and cuts or from assuming the sub-q*_T integral vanishes. This is the weakest link in the headline claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript presents resummed and matched predictions for the transverse-momentum spectra of Z and W bosons in fiducial LHC phase space at N3LL' and approximate N4LL accuracy, as implemented in SCETlib. It develops an effective-function reduction of nonperturbative TMD physics, gives a detailed five-component perturbative uncertainty estimate (resummation, fixed-order, DGLAP, matching, and recoil), compares with ATLAS and CMS 13 TeV data, and proposes the cumulative fiducial pT^Z cross section as a constraint for approximate N3LO PDF determinations. The highest-order matched predictions use an approximate O(alpha_s^3) nonsingular cross section modeled from NNLOjet Z+jet data, which is then ported to W production and other fiducial cuts.","tokens_in":57861,"tokens_out":6341,"duration_ms":63296,"significance":"If the claims hold, this is among the most complete descriptions of the Drell-Yan pT spectrum currently available, combining a rigorous treatment of fiducial power corrections with a transparent and unusually comprehensive perturbative uncertainty budget. The paper provides concrete evidence of perturbative convergence, fixed-order recovery, and a novel effective description of nonperturbative TMD effects. The main caveat is that the headline accuracy claim rests in part on an approximate and ported O(alpha_s^3) nonsingular model, so the advertised 'matching to available O(alpha_s^3) fixed-order results' needs to be either strengthened or explicitly delimited.","major_comments":[{"comment":"The O(alpha_s^3) nonsingular cross section is obtained by fitting the two-parameter model fnons(qT) to NNLOjet Z+jet data for the ATLAS 13 TeV cuts only, with the coefficient set to zero below q*_T = 8 GeV using the cancellation in Eq. (2.25). The same fnons is then applied to W+- production and to CMS Z fiducial cuts, with the fit uncertainty scaled by factors of 3 and 10. The manuscript does not provide a calculation or independent validation that the shape of fnons, or the vanishing of the sub-q*_T integral, holds for these other processes and fiducial acceptances. Because this approximate coefficient enters the matched spectra at qT above about 20 GeV and contributes to the total normalization used in Section 5, the quoted approximation uncertainty from the c1 and c2 fit errors does not cover the systematic uncertainty of the portability assumption. Please either validate the ported model with independent O(alpha_s^3) or high-statistics data in the relevant qT range, or quantify and propagate a portability uncertainty, and adjust the abstract's phrasing about matching to available O(alpha_s^3) fixed-order results accordingly.","section":"2.4, Eqs. (2.23)-(2.26)"},{"comment":"The claim that the cumulative pT^Z cross section provides strong N3LO PDF constraining power relies on the total normalization being controlled at the sub-percent level. However, the cumulative cross section at qmax_T = 20-30 GeV integrates over the region qT < 8 GeV, where the O(alpha_s^3) nonsingular is set to zero using the ATLAS-Z cancellation of Eq. (2.25). That cancellation is not established for the CMS fiducial acceptance used in the data comparison in Fig. 14, and any bias from this model assumption enters the cumulative prediction directly. Please assess the sensitivity of the cumulative cross section to alternative assumptions for the sub-q*_T nonsingular and include that variation in the quoted uncertainty.","section":"5, Figs. 14-16"},{"comment":"The central values of the nonperturbative TMD parameters (c_nu, omega_nu,q and Omega_V) are chosen for illustration and are not fitted to data. The comparison in Fig. 8 shows residual differences in the first few bins that vary with these parameters and are not covered by the perturbative band. Since the paper's stated scope includes the nonperturbative region qT ~ Lambda_QCD, the quoted perturbative uncertainties do not represent the full model uncertainty of the low-qT predictions. The paper should either profile or fit these parameters against the data shown, or state more prominently in the abstract and conclusions that the central predictions below about qT = 10 GeV are illustrative and conditional on an unconstrained model choice.","section":"4.2, Eqs. (3.9), (3.28)-(3.29)"}],"minor_comments":[{"comment":"The envelope in Eq. (4.17) is labeled Delta_match, but it is the recoil uncertainty and should be labeled Delta_recoil to be consistent with Eq. (4.18).","section":"4.1, Eq. (4.17)"},{"comment":"The two elements of V_recoil are typeset as {GJ, GJ} and are visually indistinguishable; please use an explicit overline or subscript to distinguish the anti-Gottfried-Jackson frame from the Gottfried-Jackson frame in the text and in Eq. (4.16).","section":"4.1, Eq. (4.16)"},{"comment":"The word 'conributing' in the caption of Figure 17 should be 'contributing'.","section":"Figure 17 caption"},{"comment":"The optimal fit window is given as 8 GeV <= qT <= 77.5 GeV, while the reported predictions are restricted to qT <= 65 GeV (Z) and qT <= 60 GeV (W); a sentence explaining why the fit extends beyond the prediction range would be helpful.","section":"2.4, around Eq. (2.26)"}],"recommendation":"major_revision","confidential_remarks":"This is a substantial and careful paper, and the resummation machinery appears sound. The main risk is that the abstract and Section 5 overstate the reach of the approximate O(alpha_s^3) nonsingular model, which is ported from Z+jet ATLAS cuts to W and CMS Z predictions. I would not recommend rejection, but the claims should be brought into line with what is actually computed, or the portability uncertainty should be quantified. The SCETlib code is not yet public, which is not a blocker but would strengthen reproducibility once released."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: this is a serious precision QCD paper that builds on an active program and makes two genuinely new contributions — a complete uncertainty treatment for matching and recoil scheme choices, and a rigorous reduction of the flavor/x-dependent TMD nonperturbative function space to a single effective function per vector boson and collider. The numerical predictions are careful, self-critical, and mostly convincing.\n\nThe load-bearing approximation is the O(alpha_s^3) nonsingular: a two-parameter model fitted to NNLOjet Z+jet data with ATLAS cuts, then ported to W and CMS cuts. The stress-test note is right that this is the soft spot, but I'd phrase it more gently. The authors state plainly that their fit does not meet the standard of state-of-the-art nonsingular fits, they scale the fit error by 3 (Z) and 10 (W) to cover the porting and functional-form uncertainty, and they show the resulting approximation error is ~0.2% (Z) and ~0.6% (W) in the 10-20 GeV window against a ~1% perturbative band. That is a quantified, bounded approximation, not a hidden one. The cumulative-cross-section claim does inherit any bias in the total normalization, but the integral check below 8 GeV for Z is a real consistency test, and the cumulative curve is dominated by moderate qT where the nonsingular is a small effect.\n\nThe effective-function reduction in Section 3 is genuinely useful: it clarifies exactly what information can and cannot be extracted from single-boson pT spectra, and the pocket formula is practical. The N4LL RGE kernels in the appendix are concrete, and the exact analytic solutions are referenced to published work. The comparisons to ATLAS/CMS data and multiple PDF sets are honest — the finding that NNLO PDF sets overshoot at intermediate qT while aN3LO PDFs fix it is a concrete, checkable statement.\n\nMinor soft spots: the nonperturbative parameters (Omega_V, CS-kernel model) are illustrative, not fitted, which the authors admit; and SCETlib is not yet public, so the numerics can't be independently reproduced today. Neither is disqualifying, but both argue against using these numbers as a black box for mW extraction yet.\n\nWho it's for: LHC precision practitioners working on mW, alpha_s from pT, or N3LO PDF fits. It deserves a serious referee. I'd accept after the authors add an explicit W-side validation of the ported nonsingular and a clearer statement of which parameters are extracted versus illustrative.","headline":"Serious, careful N3LL' precision phenomenology with a genuinely new effective-function result; the ported O(alpha_s^3) nonsingular is a quantified approximation, not a hidden flaw, and the paper deserves refereeing.","tokens_in":58292,"tokens_out":2053,"would_cite":true,"duration_ms":21807,"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":"This paper claims that the full LHC W and Z transverse-momentum spectra can now be predicted from QCD with percent-level perturbative uncertainties, and that a single effective function captures all nonperturbative TMD effects for each…","keywords":["Drell-Yan","transverse momentum spectrum","Sudakov resummation","TMD factorization","Collins-Soper kernel","nonperturbative TMD model","parton distribution functions","fiducial cuts"],"falsifier":"Compute the $O(\\alpha_s^3)$ nonsingular Z+jet cross section differentially in $q_T$ from 8 to 65 GeV with a fully independent NNLO calculation and compare it to the fitted $f_{\\rm nons}$ model scaled by the NLO cross section; a deviation exceeding the red uncertainty band in figure 3 would invalidate the matching. Alternatively, measure the cumulative fiducial $p_T^Z$ cross section at $q_T^{\\rm max}$ between 20 and 30 GeV with a total uncertainty below 0.3%; if the aN$^3$LO-PDF prediction misses the measured value by more than the combined uncertainty, the claimed PDF-constraining power is falsified.","tokens_in":57097,"feed_emoji":"⚛️","tokens_out":6243,"duration_ms":56075,"temperature":0.7,"pith_summary":"This paper aims to show that the transverse-momentum spectra of W and Z bosons at the LHC can now be described end to end by QCD perturbation theory, from the nonperturbative few-GeV region to the fixed-order tail, with residual perturbative uncertainties near the percent level. The authors build resummed predictions at N$^3$LL$'$ and approximate N$^4$LL order, matched to the available $O(\\alpha_s^3)$ fixed-order results, and they include a comprehensive account of uncertainties from scale choices, matching, $\\alpha_s$, PDFs, and nonperturbative TMD effects. They also claim that for a single resonantly produced vector boson at a given collider, all the flavor and Bjorken-$x$ complexity of nonperturbative TMD physics collapses into one effective function of the transverse distance $b_T$, which simplifies both phenomenology and the interpretation of single-process fits. If these claims hold, the predictions give the most complete available theoretical description of the Drell-Yan $p_T$ spectrum, and the cumulative $p_T^Z$ cross section becomes a competitive new input for N$^3$LO parton distribution functions.","feed_headline":"W and Z momentum spectra predicted to percent-level precision","feed_subtitle":"New three-loop-order QCD predictions cover the whole spectrum and sharpen PDF constraints from the LHC.","key_machinery":"The machinery is the $b_T$-space factorization theorem for $q_T$ spectra, which splits the cross section into hard, beam, and soft functions whose renormalization-group and rapidity evolutions resum the Sudakov logarithms, together with additive matching to the fixed-order nonsingular cross section through hybrid profile scales that interpolate between the canonical $b_T$ scales and the fixed-order scales. The new organizing object is the effective nonperturbative function $F^{\\rm np}(\\Theta, b_T)$, defined as a flavor-, $x$-, and rapidity-averaged ratio of TMD boundary terms against perturbative weights, which captures all leading nonperturbative TMD effects in a single one-dimensional function for fiducial $q_T$ spectra. Around it, the paper also introduces a two-parameter model for the $O(\\alpha_s^3)$ nonsingular cross section fitted to existing Z+jet fixed-order data and ported to W production.","core_discovery":"The central discovery is a consistent, state-of-the-art prediction framework: the fiducial W and Z $p_T$ spectra are computed at N$^3$LL$'$ and approximate N$^4$LL resummed accuracy matched to $O(\\alpha_s^3)$ fixed order, with all known perturbative ingredients at four and five loops where needed, and with a quantitative uncertainty budget that decomposes resummation, fixed-order, DGLAP, matching, and recoil-scheme uncertainties. A second, conceptual result is that the complicated flavor- and $x$-dependent nonperturbative TMD boundary conditions reduce, for any single resonant vector boson at a fixed collider, to a single effective function of $b_T$ (or of $b_T$ and rapidity if the rapidity is resolved), formally generalizing earlier leading-power averaging results. The paper further establishes that NNLO PDF sets systematically overshoot the measured $p_T^Z$ spectrum in the intermediate region, while approximate N$^3$LO PDF sets remove most of this discrepancy, and that the additional N$^4$LL Sudakov effects beyond N$^3$LL$'$ are numerically marginal.","pith_inferences":["The effective-function reduction should carry over to other color-singlet processes that use the same factorization, such as fiducial Higgs $p_T$ spectra, where an analogous one-function model could simplify nonperturbative uncertainties.","If the cumulative $p_T^Z$ cross section is included in a global PDF fit at aN$^3$LO, the small nonsingular pieces can be held fixed at a reference PDF set while propagating the PDF through the resummed cross section, giving an inexpensive but genuinely three-loop-accurate fitting template; this strategy is implicit in the paper's discussion and looks testable.","The spread of lattice QCD determinations of the Collins-Soper kernel at large $b_T$ could be used as an external prior in the effective model, turning the illustrative parameter variations into a quantitative constraint; the paper leaves that step to future work.","The observed preference of the data for weaker nonperturbative effects at low $q_T$ and for aN$^3$LO PDFs in the tail suggests that a combined fit of $\\alpha_s$, PDFs, and the effective nonperturbative parameter to the $p_T^Z$ spectrum could resolve the current mild discrepancies, but such a fit is not performed here."],"forward_implications":["Fiducial Z and W $p_T$ spectra can be predicted with roughly one-percent perturbative uncertainty in the 5-30 GeV range, with uncertainties carefully split into resummation, fixed-order, DGLAP, matching, and recoil components.","The comparison with LHC data shows that NNLO PDF sets overshoot the intermediate-$q_T$ spectrum by a few percent and that approximate N$^3$LO PDF sets resolve most of the discrepancy, implying the spectrum carries real three-loop PDF information.","The cumulative fiducial $p_T^Z$ cross section is predicted at sub-percent perturbative accuracy and is nearly free of nonperturbative TMD effects for $q_T^{\\rm max}$ above about 20-25 GeV, making it a clean constraint for full N$^3$LO PDF fits.","For single-boson analyses such as an $m_W$ extraction, all relevant nonperturbative TMD effects can be varied with a single effective parameter (or a single function of rapidity), decoupling the analysis from detailed flavor assumptions.","The extra N$^4$LL Sudakov terms change the prediction by less than about 0.5% relative to N$^3$LL$'$ with the same aN$^3$LO PDFs, so the practical gain of N$^4$LL is currently small."],"supporting_citations":[{"why":"Supplies the $O(\\alpha_s^3)$ Z+jet fixed-order data used to fit the approximate nonsingular model and to validate the normalization against fixed-order results.","marker":"[43]"},{"why":"Provides the factorization framework with exact leptonic kinematics and fiducial power corrections on which the resummed predictions are built.","marker":"[22]"},{"why":"The default MSHT20nnlo collinear PDF set used for the central predictions and for the baseline PDF uncertainty estimates.","marker":"[147]"},{"why":"The approximate N$^3$LO MSHT PDF set whose use removes most of the intermediate-$q_T$ discrepancy with LHC data.","marker":"[27]"},{"why":"The NNPDF approximate N$^3$LO set that confirms the PDF-related shift seen with the MSHT set.","marker":"[28]"},{"why":"The ATLAS 13 TeV fiducial $p_T^Z$ measurement used for comparison, normalization, and the shape of the nonperturbative tests.","marker":"[5]"},{"why":"The CMS 13 TeV measurement used both for the fiducial spectrum comparison and for the cumulative cross-section comparisons.","marker":"[11]"},{"why":"The earlier leading-power effective-parameter result for $O(\\Lambda_{\\rm QCD}^2/q_T^2)$ effects that this paper generalizes to all powers in $\\Lambda_{\\rm QCD}$.","marker":"[152]"},{"why":"Provides the exact analytic renormalization-group solutions for the Sudakov kernels and strong coupling used in the N$^3$LL$'$ predictions.","marker":"[139]"}],"fun_headline_variants":["W and Z pT spectra at N3LL' with full uncertainty breakdown","One effective TMD function captures all W and Z flavor complexity","Percent-level W/Z pT predictions sharpen N3LO PDF constraints","N4LL effects negligible in high-precision W/Z pT spectra"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the true $O(\\alpha_s^3)$ nonsingular cross section, which is too expensive to compute for every process and fiducial set, is well represented by the two-parameter logarithmic model fitted to Z+jet data above 8 GeV and then reused for W production and other cuts; if that model is wrong at the level of the quoted uncertainties, the matched spectra and cumulative cross sections shift beyond what the bands show.","fun_headline_variants_meta":{"raw":{"variants":["W and Z pT spectra at N3LL' with full uncertainty breakdown","One effective TMD function captures all W and Z flavor complexity","Percent-level W/Z pT predictions sharpen N3LO PDF constraints","N4LL effects negligible in high-precision W/Z pT spectra"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000706,"raw_usage":{"total_tokens":3228,"prompt_tokens":1036,"completion_tokens":2192,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":652,"completion_tokens_details":{"reasoning_tokens":2115}},"tokens_in":652,"tokens_out":2192,"duration_ms":17336,"temperature":1.0,"reasoning_tokens":2115,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:39:00.110382+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the $O(\\alpha_s^3)$ nonsingular Z+jet cross section differentially in $q_T$ from 8 to 65 GeV with a fully independent NNLO calculation and compare it to the fitted $f_{\\rm nons}$ model scaled by the NLO cross section; a deviation exceeding the red uncertainty band in figure 3 would invalidate the matching. Alternatively, measure the cumulative fiducial $p_T^Z$ cross section at $q_T^{\\rm max}$ between 20 and 30 GeV with a total uncertainty below 0.3%; if the aN$^3$LO-PDF prediction misses the measured value by more than the combined uncertainty, the claimed PDF-constraining power is falsified.","supporting_citations":[{"cited_title":"Analytic results for Sudakov form factors in QCD","cited_arxiv_id":"2110.11360","evidence_quote":"Provides the exact analytic renormalization-group solutions for the Sudakov kernels and strong coupling used in the N$^3$LL$'$ predictions."}],"review_version":1}