{"id":"ac2907c1-ab45-452c-beb9-27d3e34de8fc","arxiv_id":"1908.00985","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":8.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"The authors derive a generalized threshold factorization theorem for color-singlet hadron-collider processes, valid when only one parton is at large x, and show it outperforms soft threshold factorization.","lead":"This paper derives a new way to approximate particle-collision cross sections that works when only one of the two colliding quarks or gluons carries nearly all of its momentum, instead of requiring both to do so. The result captures more of the full quantum-chromodynamics prediction and fixes an inconsistency in earlier rapidity-dependent approximations.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The ultra-soft/Glauber cancellation asserted after Eq. (11) is load-bearing and not proven; it is inherited from inclusive collinear factorization, but Eq. (12) is qT-differential at a perturbative scale, so the gap deserves a dedicated check.","rationale":"The reader's weakest-assumption analysis identifies exactly the same load-bearing point: the unproven cancellation of ultra-soft and Glauber modes. I agree that this is the most delicate step in the derivation. However, the assumption is standard in collinear factorization, and the paper supplies substantial independent support: the analytic NLO agreement with Refs. [63,64], the NNLO numerical validation against Vrap in Fig. 3 and Fig. S2, the analytic consistency check of Eq. (20) to O(α_s^2), and the explicit mode analysis in Table S1. The qT-differential factorization Eq. (12) is not directly validated against an independent calculation, which is why I propose that specific check. The unprinted two-loop finite terms of the modified beam function are a reproducibility concern rather than a correctness risk, because the NNLO comparison against Vrap would fail if those terms were wrong. The critique of earlier rapidity-dependent soft threshold results in Sec. D is well supported by the explicit toy-model counterexamples in Eqs. (S29)-(S35). Overall, the central claim is credible and the identified concern does not rise to the level of changing the verdict; it warrants a further check rather than a rejection.","tokens_in":27840,"tokens_out":34374,"duration_ms":357688,"concrete_test":"Compute the O(α_s^2) contribution to dσ/dx_a dx_b d^2q_T for Drell-Yan in the region 1-x_a ~ q_T^2/Q^2 with generic x_b (e.g., x_b=0.1) using the factorized expression Eq. (12) with the double-differential beam functions of Refs. [55,56], and compare it to an independent NNLO qT-subtraction code. If the normalized difference does not vanish like a power of 1-x_a at fixed q_T^2/[Q^2(1-x_a)], a leading-power Glauber/usoft term is missing, invalidating Eq. (17).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The factorization theorem in Eq. (17) rests on the claim, made immediately after Eq. (11), that ultra-soft and Glauber modes cancel because the measurement is fully inclusive over perpendicular momenta at the scale λ_QCD Q. This cancellation is not demonstrated in the paper; it is asserted to follow from the standard collinear factorization theorem [1–3]. The concern is load-bearing because if any leading-power Glauber/usoft contribution survives in the generalized threshold limit, the separation into independent p̄n, P̄n, Pn, and Ps sectors in Eqs. (9) and (12) is spoiled, and Eq. (17) with it. The standard theorem applies to the fully inclusive (qT-integrated) Drell-Yan cross section, whereas Eq. (12) is explicitly differential in qT at the perturbative scale λQ. Although the supplement correctly notes that there is only a single collinear sector at the λQ scale, which excludes perturbative Glauber modes with pT ~ λQ, the cancellation of the nonperturbative Glauber modes with pT ~ λ_QCD Q is argued by scale separation rather than shown. The fixed-order NNLO validation of the qT-integrated cross section in Fig. 3 is strong and does exercise the leading singular terms, but it does not directly validate the qT-differential factorization against an independent calculation. This is the weakest link in an otherwise well-supported derivation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives a generalized threshold factorization theorem for color-singlet production, valid in the limit x_a → 1 for generic x_b (and symmetrically x_b → 1 for generic x_a). The central result, Eq. (17), expresses the cross section as a hard function times a sum of a threshold-PDF convolved with a (modified) beam function, with a soft-subtraction term for the overlap. The authors derive the theorem in SCET, introduce a new modified beam function, compute its matching coefficients through O(α_s^2) for quarks and O(α_s) for gluons, perform an analytic NLO check against published results in (z, y) variables for all partonic channels, and an NNLO numerical check against Vrap using flat PDFs. They use the theorem to predict a set of N3LO logarithmic terms and to identify the correct soft-threshold rapidity factorization among competing results in the literature.","tokens_in":28128,"tokens_out":18474,"duration_ms":182630,"significance":"If correct, this is the weakest known kinematic limit in which the process-dependent hard function factorizes, and it materially extends the domain of threshold resummation: it captures the full singular dependence in one momentum fraction while keeping the other exact, including flavor-nondiagonal channels at leading power, and it resolves conflicting soft-threshold results for rapidity spectra. The paper supplies strong evidence: the SCET derivation is standard and detailed in the supplemental material; NLO analytic agreement is shown in all partonic channels; the NNLO numerical validation against Vrap with flat PDFs shows the expected power-law vanishing of the difference; and the N3LO logarithmic predictions follow from RGE evolution rather than fitted constants. The new modified beam function is a useful ingredient in its own right. The only caveat, which I do not regard as blocking, is that the qT-differential factorization in Eq. (12) is not directly validated against an independent calculation; however, the main qT-integrated result in Eq.","major_comments":[],"minor_comments":[{"comment":"The assertion that ultra-soft and Glauber modes cancel after Eq. (11) is brief for a theorem that is qT-differential in Eq. (12); please add one or two sentences outlining the scale-separation argument (the measurement is inclusive over perpendicular momenta at the Glauber scale λ_QCD Q) and pointing to the supplemental discussion.","section":"Section II, after Eq. (11)"},{"comment":"The NNLO validation is for the qT-integrated factorization Eq. (14), not for the qT-differential theorem in Eq. (12); stating this explicitly would help the reader understand the scope of the numerical check.","section":"Section III, Fig. 3"},{"comment":"The ellipsis in Eq. (22) hides the lower logarithmic powers L0(1−z_a), L1(1−z_a), and L2(1−z_a); consider replacing it with an explicit placeholder or a note in the text that these terms are collected in the supplemental material.","section":"Equation (22)"},{"comment":"The two-loop modified beam function expressions are described as 'available from the authors upon request'; for a fully reproducible supplement, including them as auxiliary files would be preferable.","section":"Supplement, Sec. E.3"}],"recommendation":"minor_revision","confidential_remarks":"This is a well-executed and important paper with strong validation. The only point I would ask the authors to address is the ultra-soft/Glauber cancellation in the qT-differential step (Eq. (12)), which is currently supported mainly by scale separation and by the validation of the qT-integrated result; a short explicit argument or a citation to a theorem covering this case would remove the last doubt. This is a local clarification rather than a change to the central claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth your time. The paper proves a factorization for color-singlet production in the limit where only one momentum fraction goes to 1, keeping the other generic. That is genuinely new: it captures flavor-nondiagonal channels at leading power, and it contains the soft-threshold limit as a special case. As a byproduct, the authors show that several published rapidity-dependent soft-threshold results (Refs. [19-23]) miss a leading-power term. The diagnosis in Supplemental D is convincing, particularly the flaw in the Fourier-transform argument: if u ~ 1-z, then the conjugate M must be counted like 1/(1-z), not O(1), so replacing the distribution by delta(u) is not justified. That section alone is a useful contribution.\n\nThe paper earns its keep on evidence. The derivation is standard SCET but carefully laid out, with the mode table in the supplement. The NLO check is analytic for all partonic channels, and the NNLO validation against Vrap uses flat PDFs, which is exactly the right way to stress the singular terms; the difference vanishes as a power of 1-x_a, channel by channel in Fig. S2. Note that the NNLO check exercises the two-loop finite terms of the new modified beam function, not just the RGE-predicted logarithms, so the new object is genuinely tested. The N3LO predictions are a reasonable bonus, and the implementation ships with SCETlib.\n\nOn the stress-test concern about ultra-soft/Glauber cancellation after Eq. (11): it is real but not load-bearing. The assertion is inherited from the inclusive collinear factorization theorems [1-3], and the paper's scale argument (lambda_QCD << lambda, measurement fully inclusive over perpendicular momenta at the lambda_QCD Q scale) is plausible as stated. The worry that the NNLO check is qT-integrated does not actually land on the main theorem, because Eq. (17) is itself qT-integrated and therefore sits on firmer inclusive ground. And the NNLO agreement gives indirect empirical support: a missing leading-power Glauber contribution would almost certainly show up as a non-vanishing difference in Fig. 3. The qT-differential intermediate Eq. (12) is the piece most exposed; a dedicated check there would be honest, but it is not the step supporting the central claims. Minor irritants: the two-loop finite terms of the modified beam function are only partially printed (\"available from the authors upon request\"), and the axial singlet corrections are dropped as in Vrap. Both are minor. Citation pattern is fair; there is heavy self-citation, but it is to the beam-function machinery this paper builds on.\n\nYes, send it to referees. A serious referee will ask for the Glauber step to be sharpened and for more of the two-loop boundary terms to appear in print, but the result is important, the derivation is coherent, and the validation is strong.","headline":"A genuinely new factorization theorem for color-singlet production, with strong NLO and NNLO checks; the main soft spot (Glauber/usoft cancellation) is inherited from standard inclusive factorization and does not sink the central claim.","tokens_in":28661,"tokens_out":5743,"would_cite":true,"duration_ms":60564,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper establishes a generalized threshold factorization for color-singlet production that holds when only one momentum fraction tends to 1 and encodes the full singular structure to all orders.","keywords":["threshold factorization","soft-collinear effective theory","beam function","rapidity spectrum","Drell-Yan","Higgs production","large-x resummation","N3LO predictions"],"falsifier":"Take the exact $O(\\alpha_s^2)$ Drell-Yan cross section at fixed $x_b=10^{-2}$ and subtract the right-hand side of Eq. (14) evaluated with flat test PDFs; the theorem requires the residual to vanish as a power of $1-x_a$ with no logarithmically enhanced residue in any partonic channel. A sharper check is the $O(\\alpha_s^3)$ coefficient of $L_5(1-z_a)\\,\\delta(1-z_b)$, which Eq. (22) fixes in terms of known anomalous dimensions: an independent full $\\mathrm{N}^3\\mathrm{LO}$ calculation that disagrees with that coefficient would disprove the factorization.","tokens_in":27645,"feed_emoji":"⚛️","tokens_out":16813,"duration_ms":145173,"temperature":0.7,"pith_summary":"Soft-threshold factorization is a standard tool for hadron-collider predictions: it is valid when both incoming partons carry almost all of the parent proton's momentum, $x_a,x_b\\to 1$. This paper establishes the much weaker statement that it is enough for one momentum fraction to approach 1, say $x_a\\to 1$ with $x_b$ arbitrary, which is the large-rapidity edge of the spectrum for a produced color singlet such as a $Z$ boson or Higgs. In that limit the cross section factorizes into a process-dependent hard function, a threshold PDF on the large-$x$ side, and a modified beam function describing the other proton's collinear radiation, with the doubly soft overlap between the two one-sided limits subtracted once. The factorized expression contains the complete soft and/or collinear singular structure in both momentum fractions to all orders in perturbation theory, including flavor-off-diagonal channels. If correct, this gives the weakest known limit in which the hard function factorizes and supplies the large-$x$ resummation needed for rapidity-dependent predictions and PDF fits.","feed_headline":"One parton at threshold is enough to factorize the cross section","feed_subtitle":"Generalized threshold theorem captures large-rapidity singularities and predicts missing N3LO terms.","key_machinery":"The load-bearing object is the modified beam function $\\tilde B_i(\\tilde t,x)$, defined as the $\\vec{k}_T$-integrated projection of the double-differential beam function $B_i(t,\\vec{k}_T,x)$ with $\\tilde t=t+|\\vec{k}_T|^2/2$. It obeys the same renormalization-group evolution as the ordinary beam function, so its logarithmic structure is fixed, but it has new finite matching coefficients $\\tilde I_{ij}(\\tilde t,z)$; the paper computes these to $O(\\alpha_s^2)$ for quarks and $O(\\alpha_s)$ for gluons. The proof works by decomposing the effective theory into collinear, anti-collinear, soft, and threshold-PDF sectors, with the final-state radiation forced to be anti-collinear when $x_a\\to1$. The consistency relation Eq. (16) fixes the overlap: when $x_b\\to1$, the beam function must collapse to a convolution of the soft function with a threshold PDF, which is precisely the term subtracted in Eq. (17).","core_discovery":"The paper's central claim is the generalized threshold factorization theorem\n$$\\frac{\\mathrm{d}\\$\\sigma$}{\\mathrm{d}x_a\\,\\mathrm{d}x_b}=H_{ij}\\left[$f_i^{{\\mathrm{thr}}$}\\otimes\\tilde B_j+\\tilde B_i\\otimes $f_j^{{\\mathrm{thr}}$}-S\\otimes $f_i^{{\\mathrm{thr}}$}$f_j^{{\\mathrm{thr}}$}\\right],$$\nwhich holds at leading power in $1-x_a$ for generic $x_b$ (and symmetrically with $x_a\\leftrightarrow x_b$). Here $H_{ij}$ is the same hard function that appears in the classic soft-threshold limit, $f_i^{\\mathrm{thr}}$ is the threshold parton distribution function, $\\tilde B_j$ is a modified beam function obtained by integrating the double-differential $N$-jettiness beam function over transverse momentum, and $S$ is the soft function; the last term subtracts the overlap of the two one-sided limits. All convolutions are with respect to the appropriate momentum fractions and flavor indices are summed. The theorem encodes the complete leading-power singular behavior of the partonic cross section as either $z_a\\to1$ or $z_b\\to1$, including nondiagonal channels such as $qg\\to Lq$, and it reduces to the standard soft-threshold factorization when both $x_a,x_b\\to1$. The authors verify it analytically at NLO and numerically at NNLO and use it to predict a set of $\\mathrm{N}^3\\mathrm{LO}$ threshold logarithms.","pith_inferences":["The same one-sided factorization should carry over to other color-singlet-like final states, such as associated production with a jet or an identified hadron, with the modified beam function replaced by the appropriate generalized jet or fragmentation function.","The sharpest test of the theorem will be the fully differential $\\mathrm{N}^3\\mathrm{LO}$ result for Drell-Yan and $gg\\to H$: if Eq. (22)'s predicted $L_5$ and $L_4$ coefficients hold, resummation to $\\mathrm{N}^3\\mathrm{LL}$ in the generalized threshold limit becomes a practical default for rapidity spectra.","The paper's reliance on full inclusiveness over transverse momenta suggests that a $q_T$-vetoed version of the same observable is the natural place to look for corrections: there the ultra-soft and Glauber cancellation is less protected, and residual power corrections may become numerically visible."],"forward_implications":["The generalized threshold approximation reproduces the full Drell-Yan and $gg\\to H$ rapidity spectra at $O(\\alpha_s)$ and $O(\\alpha_s^2)$ in all partonic channels, whereas the standard soft-threshold approximation fails for gluon-initiated channels and gives only a poor approximation for the $q\\bar q$ channel.","At the partonic level, the leading-power generalized expansion contains the entire next-to-leading-power soft expansion, so a single order of the generalized expansion captures two orders of the soft expansion.","Equation (17) provides the resummation of large-$x$ logarithms in rapidity-dependent observables when only one PDF is probed near $x\\to 1$, the situation relevant for resummation-improved PDF fits.","The theorem predicts a rich set of $\\mathrm{N}^3\\mathrm{LO}$ threshold logarithms for any color-singlet process, including the full dependence on the non-threshold momentum fraction $z_b$, starting with the $L_5(1-z_a)\\delta(1-z_b)$ coefficient given in Eq. (22).","The two-dimensional convolution structure of Eq. (17) shows that rapidity-sensitive soft-threshold resummation cannot be reduced to a one-dimensional convolution; the paper identifies which published soft-threshold forms for rapidity spectra are incorrect and why."],"supporting_citations":[{"why":"Underpins the cancellation of ultra-soft and Glauber modes that lets the cross section separate into independent sectors.","marker":"[1–3]"},{"why":"Original two-variable threshold factorization in Mellin space, which Eq. (17) contains as the double limit $x_a,x_b\\to1$.","marker":"[5]"},{"why":"Defines the inclusive beam function and its matching onto PDFs, the key ingredient from which the modified beam function is built.","marker":"[39]"},{"why":"Supplies the beam-function renormalization-group equation and matching coefficients used to solve for the modified beam function's logarithmic structure.","marker":"[41]"},{"why":"Provides the soft-collinear effective theory whose mode decomposition is the derivation's mechanism.","marker":"[47–51]"},{"why":"Identifies the threshold PDF as the combined collinear-plus-soft matrix element, fixing the $f_i^{\\mathrm{thr}}$ factor in Eq. (17).","marker":"[53, 54]"},{"why":"Computes double-differential beam function matching coefficients from which the modified beam function's finite terms are obtained.","marker":"[55, 56]"},{"why":"Supplies the NLO partonic results in $(z,y)$ variables used for the analytic validation of Eq. (14).","marker":"[63, 64]"},{"why":"Provides the full fixed-order NNLO rapidity spectrum used as the numerical benchmark in the validation.","marker":"[68]"},{"why":"Contains the detailed derivation, the two-loop modified beam function coefficients, and the N3LO prediction.","marker":"[52]"}],"fun_headline_variants":["One parton at threshold suffices for full factorization","Generalized threshold theorem captures large-rapidity singularities","Single-sided threshold: new factorization for rapidity spectra","Factorization holds with just one parton at threshold","Rapidity logs resummed via generalized threshold factorization"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument would collapse if very soft gluon exchanges connecting the two incoming protons, the ultra-soft and Glauber modes, contribute at leading power; the paper assumes, following standard collinear factorization, that full inclusiveness over perpendicular momenta at the QCD scale makes them cancel.","fun_headline_variants_meta":{"raw":{"variants":["One parton at threshold suffices for full factorization","Generalized threshold theorem captures large-rapidity singularities","Single-sided threshold: new factorization for rapidity spectra","Factorization holds with just one parton at threshold","Rapidity logs resummed via generalized threshold factorization"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000378,"raw_usage":{"total_tokens":2087,"prompt_tokens":1100,"completion_tokens":987,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":716,"completion_tokens_details":{"reasoning_tokens":910}},"tokens_in":716,"tokens_out":987,"duration_ms":8695,"temperature":1.0,"reasoning_tokens":910,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:27:14.104422+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the exact $O(\\alpha_s^2)$ Drell-Yan cross section at fixed $x_b=10^{-2}$ and subtract the right-hand side of Eq. (14) evaluated with flat test PDFs; the theorem requires the residual to vanish as a power of $1-x_a$ with no logarithmically enhanced residue in any partonic channel. A sharper check is the $O(\\alpha_s^3)$ coefficient of $L_5(1-z_a)\\,\\delta(1-z_b)$, which Eq. (22) fixes in terms of known anomalous dimensions: an independent full $\\mathrm{N}^3\\mathrm{LO}$ calculation that disagrees with that coefficient would disprove the factorization.","supporting_citations":[{"cited_title":"Catani and L","cited_arxiv_id":null,"evidence_quote":"Original two-variable threshold factorization in Mellin space, which Eq. (17) contains as the double limit $x_a,x_b\\to1$."},{"cited_title":"Anastasiou, L","cited_arxiv_id":null,"evidence_quote":"Provides the full fixed-order NNLO rapidity spectrum used as the numerical benchmark in the validation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Contains the detailed derivation, the two-loop modified beam function coefficients, and the N3LO prediction."}],"review_version":1}