{"id":"f80f77ce-7256-4346-a3c0-3ed81a20ea8d","arxiv_id":"2608.12556","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The authors derive N3LO soft-gluon corrections for Z gamma production and present approximate N3LO LHC predictions for photon pT and rapidity distributions.","lead":"This paper calculates the next-to-next-to-next-to-leading-order soft-gluon corrections for Z gamma production at the LHC and combines them with NLO results to give approximate N3LO predictions. It provides new differential cross sections for the photon transverse momentum and rapidity at 13 and 13.6 TeV, which can serve as a Standard Model baseline for searches for new physics.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Soft-gluon dominance is asserted but never validated against the exact NNLO results that already exist for Zγ; the convergence claim is therefore conditional.","rationale":"Agree with the reader that the threshold-dominance assumption is the weakest link. The exact NNLO benchmark is the decisive test because it isolates exactly what the approximation adds: full NNLO results include both the threshold logarithms and the non-threshold terms that the aNNLO construction omits. The missing coefficients matter for reproducibility but are secondary, since they could be supplied. The comparison to exact NNLO is the one check that would make the central claim land. No change to the CONDITIONAL verdict is needed; the condition should be an explicit validation against exact NNLO and provision of the missing coefficients.","tokens_in":14798,"tokens_out":11821,"duration_ms":128603,"concrete_test":"Compute exact NNLO distributions for pp→Zγ with the same fiducial setup (p_T^γ>20 GeV, Frixione isolation R0=0.4, n=1, ε_γ=0.5, μR=μF=mZ) using a public NNLO code such as MATRIX or the implementation of Refs. [5,6] with MSHT20nnlo PDFs, and compare bin-by-bin with the aNNLO curves from Figures 9-12. If (aNNLO−NNLO_exact) exceeds the quoted aNNLO scale uncertainty or the aNNLO→aN3LO shift in any pT or rapidity bin, the threshold-only approximation fails and the convergence claim is unsupported; if it agrees within uncertainties, the approximation is validated at NNLO and the aN3LO results gain credibility.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that threshold/soft-gluon terms are overwhelmingly dominant in Zγ production, so adding them to full NLO gives reliable aNNLO/aN3LO cross sections and the small aNNLO-to-aN3LO shift demonstrates convergence. This is the load-bearing assumption. The only quantitative support is the statement in Section 3 that soft NLO terms are 'about 90%' of the total NLO corrections, for which no derivation or table entry is given, and the same dominance is assumed at NNLO/N3LO without a check. Because exact NNLO predictions for pp→Zγ are available (Refs. [5,6]) and use the same qT-subtraction formalism, the aNNLO approximation can be compared directly; the paper does not do so. Without that comparison, the apparent convergence could be an artifact of the threshold approximation rather than a property of QCD. In addition, C(2)_δ, C(3)_1, and C(3)_0, which are needed to reproduce the aNNLO/aN3LO numbers, are relegated to supplementary material, so the numerical central results are not independently reproducible from the submission.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes soft-gluon (threshold) corrections for Zγ production through N3LO in one-particle-inclusive (1PI) kinematics, using the Kidonakis–Sterman resummation framework. These corrections are added to fixed-order NLO results from MadGraph5 aMC@NLO to construct approximate NNLO (aNNLO) and approximate N3LO (aN3LO) predictions for the photon transverse momentum and rapidity distributions at 13 and 13.6 TeV, with MSHT20 nnlo and an3lo PDF sets. The paper reports integrated and differential cross sections, K-factors, and scale/PDF uncertainties, and concludes that the difference between aNNLO and aN3LO is small, indicating good perturbative convergence.","tokens_in":15053,"tokens_out":2802,"duration_ms":26541,"significance":"If the threshold-dominance assumption holds, this paper provides a new set of N3LO soft-gluon corrections for Zγ production, which could be valuable for interpreting LHC measurements and for searches for anomalous neutral triple gauge couplings. The work follows a well-established published framework, gives explicit analytic expressions for many of the NLO and NNLO coefficients, and presents results for two LHC energies and two PDF sets, which is useful for phenomenology. The main limitation is that the central approximation—that soft-gluon terms dominate the full higher-order corrections across the entire considered pT range—is asserted but not validated against existing exact NNLO results for this process. The convergence claim therefore remains conditional on an unquantified approximation.","major_comments":[{"comment":"The paper claims that soft-gluon corrections account for 'about 90% of the total NLO corrections' and are 'overwhelmingly dominant numerically', but no derivation, table, or quantitative comparison is provided. This dominance is then assumed to hold at NNLO and N3LO. Since exact NNLO predictions for pp→Zγ exist (Refs. [5,6]) and use the same qT-subtraction formalism, the aNNLO approximation can be compared directly with the exact NNLO result. The authors should add such a comparison (e.g., a bin-by-bin ratio of aNNLO to full NNLO for both dσ/dp_T^γ and dσ/d|y^γ|) to substantiate the load-bearing assumption. Without this validation, the small aNNLO-to-aN3LO difference in the conclusions could be an artifact of the threshold approximation rather than evidence of perturbative convergence.","section":"Section 3, paragraph 4; Section 1, paragraph 4"},{"comment":"The coefficients C(2)_δ, C(3)_1, and C(3)_0 are required for the numerical aNNLO and aN3LO predictions, but their expressions are only described as 'quite long' and relegated to supplementary material. Consequently, the central numerical results are not independently reproducible from the arXiv submission. Please include these expressions in the appendix, or at least provide the numerical values used for each coefficient at the specific scales and kinematic points that enter the calculation.","section":"Appendix, after Eq. (A.8) and Eq. (A.12)"},{"comment":"The scale uncertainties for aNNLO and aN3LO are described as resulting from 'quadrature addition of scale uncertainties from each successive soft-gluon contribution'. This procedure is not standard and is not fully documented. It is unclear whether this is equivalent to a 3-point variation of the complete approximate cross section or a separate propagation of uncertainties on individual soft-gluon terms. The authors should specify the exact algorithm used and justify its statistical interpretation; otherwise the quoted aNNLO and aN3LO uncertainties cannot be assessed.","section":"Section 3, paragraph after Eq. (3.2) and Table 1"},{"comment":"The threshold approximation drops the O((1-x_a)(1-x_b)) term in Eq. (2.5), yet the calculation is applied over the full pT range including the lowest bins (30 GeV) where partonic momentum fractions are not necessarily close to 1. The inaccuracy introduced by this approximation is never quantified. A concrete test would be to compute the NLO result both in the exact kinematics and in the threshold approximation, and to report the ratio across the pT bins presented in Eq. (3.4); this would indicate whether the approximation is reliable beyond the high-pT tail.","section":"Eq. (2.5) and Section 3, Eq. (3.4)"}],"minor_comments":[{"comment":"The statement 'the soft-gluon corrections at NLO account for about 90% of the total NLO corrections' is not supported by any table, figure, or reference; it should either be backed by a quantitative comparison or removed.","section":"Section 3, first paragraph"},{"comment":"The equation appears garbled in the submitted text; please check that the δ-function argument and the overall normalization are typeset correctly.","section":"Section 2, Eq. (2.9)"},{"comment":"The table would be easier to read if the rows were grouped by collider energy and the PDF-order labels (nnlo vs an3lo) were explained in the caption or a footnote.","section":"Table 1"},{"comment":"The conclusion states that the aN3LO scale uncertainty 'stays well below the LO uncertainty'; this is true for the integrated cross section, but for the high-pT bins the scale uncertainty at aN3LO (26–27%) exceeds the LO per-bin uncertainty (2–9%). The statement should be qualified accordingly.","section":"Section 4"},{"comment":"Reference [2] is dated two years after the submission date of this manuscript; please verify that the arXiv number and year are correct.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript follows the authors' well-established framework and is likely technically sound in its internal derivation. The central issue is the lack of validation of the threshold-dominance assumption against existing exact NNLO results, which is a fixable omission rather than a fatal error. The missing coefficients in the appendix also hamper reproducibility and should be addressed before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The new thing here is concrete: the N3LO soft-gluon coefficients for Zγ production in 1PI kinematics, and the aN3LO differential predictions that come from them. The formalism is not new—it is the same Kidonakis-Sterman framework used in the Wγ paper [21]—but the process-specific application is, and the numerical results (fiducial cross sections, pT and rapidity distributions, K-factors) are the sort of thing people in the LHC Zγ analysis community will want to use as a baseline. That is a fair and useful contribution.\n\nThe paper is also honestly presented. The setup is clear, the PDF choices and scale variation procedures are standard, and the authors do not oversell aN3LO as exact. The analytic coefficients that do appear in the appendix look plausible given the known results for other processes, and I saw no sign of circular reasoning or fitting to data.\n\nThat said, there are two soft spots, and one of them matters. First, three coefficients needed to reproduce the central numerical results—C(2)_δ, C(3)_1, C(3)_0—are relegated to supplementary material. That makes the numbers in the tables and figures not independently checkable from the submission. This is a real reproducibility issue, but a minor one if the supplementary material is actually provided.\n\nThe second issue is more substantive. The paper asserts that soft-gluon terms are “overwhelmingly dominant” and that they account for “about 90%” of the NLO corrections, but no derivation or table entry is given for that number, and the same dominance is assumed at NNLO and N3LO. The authors could have checked: exact NNLO predictions for Zγ already exist in Refs. [5,6]. Comparing the aNNLO result against that exact NNLO result would test whether the threshold approximation is actually doing what the paper claims. Without that comparison, the small aNNLO-to-aN3LO shift, which is presented as evidence of good convergence, is weaker evidence than the authors imply. It shows the approximate series is internally stable, but not necessarily that the approximation is close to the true QCD result.\n\nThe growth of scale uncertainties at aN3LO is explained by quadrature addition and is not a red flag by itself, though it tempers the “precision” claim somewhat.\n\nWho gets value from this: phenomenologists working on Zγ measurements, especially those using the process to constrain anomalous triple gauge couplings. They will get an updated, state-of-the-art approximate N3LO baseline, provided they keep the caveats in mind.\n\nRecommendation: this paper deserves peer review—not a desk rejection. The calculation is legitimate, the topic is relevant, and the missing pieces are fixable. I would send it to a referee and ask for three things before publication: the missing coefficients in the text or a clearly available supplement, a quantitative justification for the 90% dominance claim, and ideally a comparison of aNNLO with the exact NNLO results from Grazzini et al. If those are added, the paper would be a solid contribution to the Zγ literature.","headline":"A workmanlike application of well-established threshold resummation to Zγ, with a genuinely new aN3LO result; the main weaknesses are missing coefficients and an unvalidated dominance claim that the authors can fix.","tokens_in":15523,"tokens_out":2072,"would_cite":true,"duration_ms":23128,"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":"N$^3$LO soft-gluon corrections set LHC $Z\\gamma$ rates","keywords":["soft-gluon resummation","threshold logarithms","Z gamma production","N3LO QCD","one-particle-inclusive kinematics","LHC","photon transverse momentum","anomalous neutral triple gauge couplings"],"falsifier":"Compare the aNNLO predictions bin by bin with the exact NNLO calculation of $Z\\gamma$ production using identical cuts, binning, and PDFs; if the differences exceed the quoted scale uncertainties in the low-$p_T^\\gamma$ bins, the threshold-dominance assumption fails. A future complete N$^3$LO calculation that differs from aN$^3$LO by more than the quoted uncertainties would also falsify the approximation.","tokens_in":14639,"feed_emoji":"⚛️","tokens_out":7288,"duration_ms":62500,"temperature":0.7,"pith_summary":"This paper sets out to show that soft-gluon (threshold) corrections dominate the higher-order QCD series for $Z\\gamma$ production, and that their resummation through N$^3$LO yields concrete predictions for LHC observables. Working in one-particle-inclusive (1PI) kinematics, the authors compute the N$^3$LO soft-gluon coefficients, add them to the NLO result, and present approximate NNLO and approximate N$^3$LO differential cross sections in photon transverse momentum and rapidity at 13 and 13.6 TeV. The integrated fiducial cross sections for $p_T^\\gamma>20$ GeV come out around 74--79 pb, with the aN$^3$LO/NLO $K$-factor reaching roughly 1.33 in the high-$p_T^\\gamma$ tail and about 1.11 in the rapidity distribution. The small difference between the aNNLO and aN$^3$LO results is presented as evidence of perturbative convergence, which matters because precise $Z\\gamma$ predictions are a baseline for tests of the $ZZ\\gamma$ and $Z\\gamma\\gamma$ vertices and for new-physics searches.","feed_headline":"N3LO soft-gluon corrections set LHC Z-gamma rates","feed_subtitle":"Threshold logs through N3LO lift photon-$p_T$ and rapidity predictions while keeping scale uncertainty under control.","key_machinery":"The load-bearing object is the soft-gluon resummation formula for the partonic cross section in 1PI kinematics, meaning the photon is the observed final-state particle and the $Z$ boson is integrated over. The cross section is refactorized into a hard function, a soft function describing wide-angle soft gluons emitted from Wilson lines, and collinear functions, and each piece evolves under renormalization-group equations; the lightlike cusp anomalous dimension $A_q$ and the soft exponent $D_q$ generate the threshold logarithms. Expanding the resummed expression in $\\alpha_s$ through N$^3$LO and inverting the Laplace transform produces explicit coefficients $C_k^{(n)}$ for the soft-plus-virtual terms, including delta-function pieces, which are then added order by order to fixed-order results. The exponents and anomalous dimensions are taken from the existing single-particle-inclusive resummation framework, so the new work is the application to the $Z\\gamma$ hard function and the derivation of the N$^3$LO coefficients for this process.","core_discovery":"The central claim is that threshold logarithms of the form $\\ln^k(s_4/m_Z^2)/s_4$, where $s_4=s+t+u-m_Z^2$ is the partonic threshold variable, control the perturbative expansion of $pp\\to Z\\gamma$, and that these logarithms can be resummed through N$^3$LO in 1PI kinematics with the photon as the observed particle. The authors factorize the cross section into hard, soft, and collinear functions under a Laplace transform, expand the resummed expression to order $\\alpha_s^3$, invert back to momentum space, and add the resulting soft-plus-virtual terms to an NLO calculation. They find that soft contributions account for about 90% of the total NLO corrections, that the aNNLO/NLO and aN$^3$LO/NLO $K$-factors rise from roughly 1.09 and 1.11 at low $p_T^\\gamma$ to about 1.33 in the high-$p_T^\\gamma$ tail, and that the rapidity distribution receives almost uniform corrections around 1.09--1.11. The closeness of the aNNLO and aN$^3$LO curves is used to argue that the perturbative series is converging and that the predictions are stable under PDF and scale choices.","pith_inferences":["Because the threshold approximation is most accurate when $s_4$ is small, the low-$p_T^\\gamma$ bins are where the aNNLO and aN$^3$LO results are least protected; a bin-by-bin comparison with the existing exact NNLO calculation would quantify the error of the approximation.","The same 1PI resummation machinery applied here to $Z\\gamma$ should carry over to other two-body final states with a photon or heavy boson, so the N$^3$LO coefficients are a template for related processes.","If future 13.6 TeV data match the aN$^3$LO high-$p_T^\\gamma$ shape, that would corroborate the claim of threshold-logarithm dominance; if not, non-threshold contributions would need a full N$^3$LO calculation.","The growth of scale uncertainty with each added soft term suggests the approximate series has an intrinsic flattening point, so extending the same calculation to N$^4$LO would test whether the apparent convergence continues."],"forward_implications":["The aN$^3$LO differential distributions provide a fuller Standard Model baseline for $Z\\gamma$ searches for anomalous neutral triple gauge couplings.","In the high-$p_T^\\gamma$ tail, where threshold logarithms grow, the predicted $K$-factor of about 1.33 relative to NLO means measurements at 13 TeV should see noticeably higher rates than plain NLO would suggest.","The scale uncertainty pattern, which drops sharply at NLO and then grows again at aNNLO and aN$^3$LO, reflects the approximate nature of the added terms rather than a complete fixed-order calculation.","At 13.6 TeV the cross sections are about 5--6% larger than at 13 TeV at every order, so the energy dependence of the $K$-factors is mild."],"supporting_citations":[{"why":"Introduces the summation of large threshold logarithms on which the whole resummation rests.","marker":"[8]"},{"why":"Develops the treatment of subleading logarithms in QCD hard-scattering resummation.","marker":"[9]"},{"why":"Establishes the factorization of the cross section into hard, soft, and collinear functions.","marker":"[10]"},{"why":"Provides the single-particle-inclusive threshold resummation formalism used here.","marker":"[12]"},{"why":"Supplies the N$^3$LO soft-gluon exponents for 1PI kinematics that the paper extends to $Z\\gamma$.","marker":"[19]"},{"why":"Gives the notation and exponent expressions for soft and virtual corrections used in this paper.","marker":"[21]"},{"why":"Generates the LO and NLO fixed-order cross sections that receive the soft-gluon additions.","marker":"[22]"},{"why":"Defines the smooth-cone photon isolation criterion applied at NLO.","marker":"[24]"},{"why":"Provides the central parton distribution set used for the numerical predictions.","marker":"[25]"},{"why":"Provides the approximate N$^3$LO parton distribution set used to check PDF dependence.","marker":"[26]"}],"fun_headline_variants":["N3LO soft-gluon corrections sharpen Z-gamma high-pT predictions","N3LO threshold logs refine Z-gamma production at LHC","Resummed soft-gluon terms to N3LO improve Z-gamma cross sections","N3LO soft corrections cut high-pT uncertainty in Z-gamma","Threshold logs at N3LO refine Z-gamma predictions for LHC"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole approximation rests on the assumption that soft, low-energy gluon emissions near the kinematic limit dominate the higher-order QCD corrections at every photon transverse momentum and rapidity studied, so that adding just those threshold terms to the NLO result is enough to mimic the full NNLO and N$^3$LO answers.","fun_headline_variants_meta":{"raw":{"variants":["N3LO soft-gluon corrections sharpen Z-gamma high-pT predictions","N3LO threshold logs refine Z-gamma production at LHC","Resummed soft-gluon terms to N3LO improve Z-gamma cross sections","N3LO soft corrections cut high-pT uncertainty in Z-gamma","Threshold logs at N3LO refine Z-gamma predictions for LHC"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001078,"raw_usage":{"total_tokens":4527,"prompt_tokens":978,"completion_tokens":3549,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":594,"completion_tokens_details":{"reasoning_tokens":3448}},"tokens_in":594,"tokens_out":3549,"duration_ms":22505,"temperature":1.0,"reasoning_tokens":3448,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:05:17.153304+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compare the aNNLO predictions bin by bin with the exact NNLO calculation of $Z\\gamma$ production using identical cuts, binning, and PDFs; if the differences exceed the quoted scale uncertainties in the low-$p_T^\\gamma$ bins, the threshold-dominance assumption fails. A future complete N$^3$LO calculation that differs from aN$^3$LO by more than the quoted uncertainties would also falsify the approximation.","supporting_citations":[{"cited_title":"Sterman, Summation of large corrections to short-distance hadronic cross sections, Nucl","cited_arxiv_id":null,"evidence_quote":"Introduces the summation of large threshold logarithms on which the whole resummation rests."},{"cited_title":"Kidonakis and G","cited_arxiv_id":null,"evidence_quote":"Develops the treatment of subleading logarithms in QCD hard-scattering resummation."},{"cited_title":"N$^3$LO soft-gluon corrections in single-particle-inclusive kinematics and $H^+ H^-$ production","cited_arxiv_id":"2404.00089","evidence_quote":"Supplies the N$^3$LO soft-gluon exponents for 1PI kinematics that the paper extends to $Z\\gamma$."},{"cited_title":"Kidonakis and A","cited_arxiv_id":null,"evidence_quote":"Gives the notation and exponent expressions for soft and virtual corrections used in this paper."}],"review_version":1}