Pith. sign in

REVIEW 4 major objections 5 minor 26 references

N$^3$LO soft-gluon corrections for $Z\gamma$ production

T0 review · 4 major / 5 minor · reviewed 2026-08-16 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read N$^3$LO soft-gluon corrections set LHC $Z\gamma$ rates

desk verdict 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. read the letter →

arxiv 2608.12556 v1 pith:MFTV4SZY submitted 2026-08-12 hep-ph

Nikolaos Kidonakis , Kaan Şimşek This is my paper · ORCID
classification hep-ph
keywords soft-gluonresummationthresholdlogarithmsZgammaproductionN3LOQCDone-particle-inclusivekinematicsLHCphotontransversemomentumanomalousneutraltriplegaugecouplings
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

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.

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 (4)
  1. [Section 3, paragraph 4; Section 1, paragraph 4] 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.
  2. [Appendix, after Eq. (A.8) and Eq. (A.12)] 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.
  3. [Section 3, paragraph after Eq. (3.2) and Table 1] 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.
  4. [Eq. (2.5) and Section 3, Eq. (3.4)] 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.
minor comments (5)
  1. [Section 3, first paragraph] 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.
  2. [Section 2, Eq. (2.9)] The equation appears garbled in the submitted text; please check that the δ-function argument and the overall normalization are typeset correctly.
  3. [Table 1] 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.
  4. [Section 4] 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.
  5. [References] Reference [2] is dated two years after the submission date of this manuscript; please verify that the arXiv number and year are correct.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the aNNLO/aN3LO predictions are constructed from threshold resummation with no fitting to the target Z gamma observables.

full rationale

The paper's derivation chain is self-contained against the target measurement: the N3LO soft-gluon corrections are obtained by expanding the resummed cross section in Eq. (2.12), with the exponent functions taken from the authors' previously published framework (Eq. (2.13) and the text 'We use the same notation and expressions for the exponents as in Ref. [21]'). This is a reuse of established, process-independent anomalous-dimension results rather than a definition of the predicted quantity in terms of itself. The aNNLO and aN3LO predictions are formed by adding these analytically derived corrections to full NLO results from MadGraph5_aMC@NLO; no parameter is fitted to Z gamma data, and the 'about 90%' soft-dominance statement in Section 3 is an empirical characterization of the NLO calculation, not an input assumption that forces the higher-order result. The convergence claim ('The difference between the aNNLO and aN3LO results is small... shows good convergence') is conditional on the threshold approximation being numerically dominant, which is a validity or correctness concern, not a circularity concern. The self-citations to Refs. [8-21] are legitimate reuse of a published resummation formalism; nothing in the paper's central claim reduces by construction to an input or to a self-citation chain.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The calculation introduces no fitted parameters in this paper. It relies on a previously established resummation framework, standard input constants and PDFs, and the assumed dominance of soft corrections. The main unverified burden is the availability and correctness of the omitted N3LO coefficients.

assumptions (4)
  • domain assumption Factorization and refactorization of the cross section into hard, soft, and collinear functions in 1PI kinematics.
    Invoked in Eqs. (2.9) to (2.12) and taken from Refs. [8-21]; no proof is given in this paper.
  • domain assumption The soft anomalous dimension and D_q exponents for q qbar to Z gamma are the same as in the generic q qbar to color-singlet case and are reused from Ref. [21].
    Section 2 states: "We use the same notation and expressions for the exponents as in Ref. [21]."
  • ad hoc to paper Non-soft contributions beyond NLO are negligible for aNNLO and aN3LO accuracy.
    The paper asserts soft corrections are "overwhelmingly dominant" in Section 3 but does not quantify the omitted non-soft NNLO and N3LO terms.
  • domain assumption MSHT20 nnlo and MSHT20 an3lo PDF sets are valid inputs at all perturbative orders considered.
    All numerical results in Section 3 use these PDFs, with the PDF-order choice held fixed across orders to isolate QCD corrections.

how reviews work

0 comments
Cite this review

Pith. "Pith review of N$^3$LO soft-gluon corrections for $Z\gamma$ production." pith.science (2026). https://pith.science/paper/MFTV4SZY

@misc{pith2026260812556,
  author       = {Pith},
  title        = {Pith review of: N$^3$LO soft-gluon corrections for $Z\gamma$ production},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/MFTV4SZY}},
  note         = {Machine review of arXiv:2608.12556}
}
abstract

We calculate soft-gluon corrections for $Z\gamma$ production through next-to-next-to-next-to-leading-order (N$^3$LO) in one-particle-inclusive (1PI) kinematics. We use these results to produce approximate N$^3$LO (aN$^3$LO) QCD predictions for $Z\gamma$ production cross sections in proton-proton collisions at LHC energies. Differential cross sections in photon transverse momentum $p_T^\gamma$ and rapidity $y^\gamma$ are computed at 13 TeV and 13.6 TeV collider energies. Higher-order $K$-factors are examined systematically, and scale and PDF uncertainties are quantified at each perturbative order. Soft-gluon corrections are found to be significant in the high-$p_T$ region, with good perturbative convergence observed at aN$^3$LO.

Figures

Figures reproduced from arXiv: 2608.12556 by the authors.

Figure 1
Figure 1. Single-differential p γ T distributions at 13 TeV with MSHT20 nnlo PDFs, showing scale and PDF uncertainties, at LO (top left), NLO QCD (top right), aNNLO QCD (bottom left), and aN3LO QCD (bottom right). 7 [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. The same as Figure 1 but at 13.6 TeV [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. The same as Figure 1 but with MSHT20 an3lo PDFs. [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figures from the paper (9 more)
Figure 4
Figure 4. Figure 4: The same as Figure 3 but at 13.6 TeV [PITH_FULL_IMAGE:figures/full_fig_p010_4.png]
Figure 5
Figure 5. Figure 5: Single-differential |y γ | distributions at 13 TeV with MSHT20 nnlo PDFs, showing scale and PDF uncertainties, at LO (top left), NLO QCD (top right), aNNLO QCD (bottom left), and aN3LO QCD (bottom right). 11 [PITH_FULL_IMAGE:figures/full_fig_p011_5.png]
Figure 6
Figure 6. Figure 6: The same as Figure 5 but at 13.6 TeV [PITH_FULL_IMAGE:figures/full_fig_p012_6.png]
Figure 7
Figure 7. Figure 7: The same as Figure 5 but with MSHT20 an3lo PDFs. [PITH_FULL_IMAGE:figures/full_fig_p013_7.png]
Figure 8
Figure 8. Figure 8: The same as Figure 7 but at 13.6 TeV [PITH_FULL_IMAGE:figures/full_fig_p014_8.png]
Figure 9
Figure 9. Figure 9: The p γ T (left) and |y γ | (right) distributions at 13 TeV with MSHT20 nnlo PDFs at LO (dotted), NLO QCD (dot-dashed), aNNLO QCD (dashed), and aN3LO QCD (solid). The inset shows the K-factors relative to NLO QCD, with aNNLO QCD (dashed) and aN3LO QCD (solid). 50 100 5…
Figure 10
Figure 10. Figure 10: The same as Figure 9 but at 13.6 TeV [PITH_FULL_IMAGE:figures/full_fig_p015_10.png]
Figure 11
Figure 11. Figure 11: The same as Figure 9 but with MSHT20 an3lo PDFs. [PITH_FULL_IMAGE:figures/full_fig_p016_11.png]
Figure 12
Figure 12. Figure 12: The same as Figure 11 but at 13.6 TeV [PITH_FULL_IMAGE:figures/full_fig_p016_12.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

26 extracted references · 11 canonical work pages

  1. [21]

    Kidonakis and A

    N. Kidonakis and A. Tonero, Higher-order soft and virtual corrections inpp→γWproduction at the LHC, Eur. Phys. J. C85, 1270 (2025) [arXiv:2506.01590]

    work page arXiv 2025

  2. [1]

    ATLAS Collaboration, Measurement of theZγ→ν¯νγproduction cross section inppcolli- sions at√s= 13 TeV with the ATLAS detector and limits on anomalous triple gauge-boson couplings, JHEP12, 010 (2018) [arXiv:1810.04995]

    arXiv 2018

  3. [2]

    CMS Collaboration, Measurement of theZγproduction cross section and search for anomalous neutral triple gauge couplings inppcollisions at √s= 13 TeV, arXiv:2601.14102

    work page arXiv

  4. [3]

    Degrande, A basis of dimension-eight operators for anomalous neutral triple gauge boson interactions, JHEP02, 101 (2014) [arXiv:1308.6323]

    C. Degrande, A basis of dimension-eight operators for anomalous neutral triple gauge boson interactions, JHEP02, 101 (2014) [arXiv:1308.6323]

    arXiv 2014

  5. [4]

    Ohnemus, Order-αs calculations of hadronicW ±γandZγproduction, Phys

    J. Ohnemus, Order-αs calculations of hadronicW ±γandZγproduction, Phys. Rev. D47, 940 (1993)

    work page 1993

  6. [5]

    Grazzini, S

    M. Grazzini, S. Kallweit, D. Rathlev, and A. Torre,Zγproduction at hadron colliders in NNLO QCD, Phys. Lett. B731, 204 (2014) [arXiv:1309.7000]

    arXiv 2014

  7. [6]

    Grazzini, S

    M. Grazzini, S. Kallweit, and D. Rathlev,WγandZγproduction at the LHC in NNLO QCD, JHEP07, 085 (2015) [arXiv:1504.01330]

    arXiv 2015

  8. [7]

    ATLAS Collaboration, Measurement of theZ(→l +l−)γproduction cross-section inppcolli- sions at√s= 13 TeV with the ATLAS detector, JHEP03, 054 (2020) [arXiv:1911.04813]

    work page Pith review arXiv 2020

Show all 26 references
  1. [8]

    Sterman, Summation of large corrections to short-distance hadronic cross sections, Nucl

    G. Sterman, Summation of large corrections to short-distance hadronic cross sections, Nucl. Phys. B281, 310 (1987)

    1987

  2. [9]

    Kidonakis and G

    N. Kidonakis and G. Sterman, Subleading logarithms in QCD hard scattering, Phys. Lett. B 387, 867 (1996)

    1996

  3. [10]

    Kidonakis and G

    N. Kidonakis and G. Sterman, Resummation for QCD hard scattering, Nucl. Phys. B505, 321 (1997) [arXiv:hep-ph/9705234]

    1997 arXiv

  4. [11]

    Kidonakis, G

    N. Kidonakis, G. Oderda, and G. Sterman, Evolution of color exchange in QCD hard scattering, Nucl. Phys. B531, 365 (1998) [arXiv:hep-ph/9803241]

    1998 arXiv

  5. [12]

    Laenen, G

    E. Laenen, G. Oderda, and G. Sterman, Resummation of threshold corrections for single- particle inclusive cross sections, Phys. Lett. B438, 173 (1998) [arXiv:hep-ph/9806467]

    1998 arXiv

  6. [13]

    Kidonakis, Collinear and soft gluon corrections to Higgs production at next-to-next-to-next- to-leading order, Phys

    N. Kidonakis, Collinear and soft gluon corrections to Higgs production at next-to-next-to-next- to-leading order, Phys. Rev. D77, 053008 (2008) [arXiv:0711.0142]

    2008 arXiv

  7. [14]

    Kidonakis, Two-loop soft anomalous dimensions for single top quark associated production with aW − orH − , Phys

    N. Kidonakis, Two-loop soft anomalous dimensions for single top quark associated production with aW − orH − , Phys. Rev. D82, 054018 (2010) [arXiv:1005.4451]

    2010 arXiv

  8. [15]

    Kidonakis and R.J

    N. Kidonakis and R.J. Gonsalves, NNLO soft-gluon corrections for theZboson andWboson transverse momentum distributions, Phys. Rev. D89, 094022 (2014) [arXiv:1404.4302]

    2014 arXiv

  9. [16]

    Kidonakis, Higher-order corrections fortZproduction via anomalous couplings, Phys

    N. Kidonakis, Higher-order corrections fortZproduction via anomalous couplings, Phys. Rev. D97, 034028 (2018) [arXiv:1712.01144]

    2018 arXiv

  10. [17]

    Forslund and N

    M. Forslund and N. Kidonakis, Resummation for 2→nprocesses in single-particle-inclusive kinematics, Phys. Rev. D102, 034006 (2020) [arXiv:2003.09021]. 20

    2020 arXiv

  11. [18]

    Kidonakis and N

    N. Kidonakis and N. Yamanaka, Soft-gluon corrections fortqZproduction, Phys. Lett. B838, 137708 (2023) [arXiv:2210.09542]

    2023 arXiv

  12. [19]

    Kidonakis and A

    N. Kidonakis and A. Tonero, N3LO soft-gluon corrections in single-particle-inclusive kinematics andH +H− production, JHEP06, 138 (2024) [arXiv:2404.00089]

    2024 arXiv

  13. [20]

    Kidonakis and C

    N. Kidonakis and C. Foster, Higher-order soft-gluon corrections fort ¯tZcross sections, Phys. Lett. B860, 139146 (2025) [arXiv:2410.01214]

    2025 arXiv

  14. [22]

    J. Alwallet al., The automated computation of tree-level and next-to-leading order differen- tial cross sections, and their matching to parton shower simulations, JHEP07, 079 (2014) [arXiv:1405.0301]

    2014 arXiv

  15. [23]

    Frederixet al., The automation of next-to-leading order electroweak calculations, JHEP07, 185 (2018) [Erratum: JHEP11, 085 (2021)] [arXiv:1804.10017]

    R. Frederixet al., The automation of next-to-leading order electroweak calculations, JHEP07, 185 (2018) [Erratum: JHEP11, 085 (2021)] [arXiv:1804.10017]

    2018 arXiv

  16. [24]

    Frixione, Isolated photons in perturbative QCD, Phys

    S. Frixione, Isolated photons in perturbative QCD, Phys. Lett. B429, 369 (1998) [arXiv:hep- ph/9801442]

    1998

  17. [25]

    Bailey, T

    S. Bailey, T. Cridge, L.A. Harland-Lang, A.D. Martin, and R.S. Thorne, Parton distributions from LHC, HERA, Tevatron and fixed target data: MSHT20 PDFs, Eur. Phys. J. C81, 341 (2021) [arXiv:2012.04684]

    2021 arXiv

  18. [26]

    McGowan, T

    J. McGowan, T. Cridge, L.A. Harland-Lang, and R.S. Thorne, Approximate N 3LO parton distribution functions with theoretical uncertainties: MSHT20aN 3LO PDFs, Eur. Phys. J. C 83, 185 (2023) [arXiv:2207.04739]. 21

    2023 arXiv

Pith tools

Reviewed August 16, 2026 · model on record in the stance chip above.