{"id":"a688f003-b846-44d4-ab84-519ac3be49c6","arxiv_id":"1908.02179","paper_version":4,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Resumming next-to-leading-power non-relativistic QCD corrections near the top-pair threshold raises the predicted cross section by about 9% and shifts the extracted top mass by about 1.5 GeV.","lead":"This paper calculates missing higher-order corrections to the production rate of top-quark pairs near the threshold energy and shows they raise the predicted rate by about 9%, bringing theory closer to CMS data. A reader should care because the size of this correction shifts the extracted top-quark mass by about 1.5 GeV, toward the directly measured value.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No-soft-function assertion in Eq. (6) is the load-bearing premise; it is asserted, not demonstrated, and if wrong the +9% enhancement and 1.5 GeV shift are misattributed.","rationale":"The central quantitative claim is a 9% enhancement of the NNLO cross section and a 1.5 GeV shift in m_t, which follows from the NLP resummation of Eq. (6). The structure of Eq. (6) is H×J, with no soft function. This is the least secure element: the argument given in the text addresses only ultrasoft initial-final exchanges, and the soft modes listed in Eq. (5) are not shown to be confined to the top-antitop pair. A soft function coupling initial and final states would break the product form and alter the size and z-dependence of the resummed correction; the 9% and 1.5 GeV numbers would change. The paper does provide useful internal checks: the β→0 limit of the NLP expansion is claimed to reproduce the exact NLO correction in [340,380] GeV (Fig. 2), and the width-insensitivity of the [300,380] integral is noted. However, the NLO hard functions are not given, so these checks cannot be reproduced. I therefore agree with the reader that the verdict should be conditional, with the condition made more specific: demonstrate the absence of a soft function and provide the analytic NLO hard functions. The most decisive single test is an independent derivation of H^(1) and a factorization test of the β→0 NLO cross section; if the product form survives, the concern is resolved.","tokens_in":9332,"tokens_out":14814,"duration_ms":176516,"concrete_test":"Compute the known NLO q qbar -> t tbar partonic cross section in the beta->0 limit with z fixed away from 1, and test whether the singular (1/β and ln β) terms factorize as H^(0)(z)·J^(1)(E) + H^(1)(z)·J^(0)(E). Concretely: for fixed z<1, Q_T=0, extract the β-singular coefficient; if it depends on z in a way that is not proportional to H^(0)(z) (i.e., requires a convolution or a soft subtraction), then a soft function is present at NLP. Recompute Table I with the resulting corrected formula; if the NNLO+NLP average changes by more than ~5%, the shift estimate fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing premise is the factorization dσ_NLP = H(z,Q_T,Y)·J(E) in Eq. (6), specifically the assertion immediately after Eq. (6) that \"there is no soft function in our factorization formula.\" The stated justification covers only ultrasoft (k~M_tt β^2) exchanges among initial-state partons and the pair. But the mode list Eq. (5) also contains soft modes with k~M_tt β, and J is defined to absorb only potential/soft/ultrasoft exchanges between the two heavy quarks, not exchanges connecting the colored initial-state partons to the pair. A soft gluon exchanged between an eikonal initial parton and the color-octet or color-singlet tbar pair has the same momentum scaling and is not power-suppressed by any argument given. If such a mode contributes at NLP, Eq. (6) must be replaced by H ⊗ S ⊗ J with a convolution in the soft momentum, entangling z and E; then the beta->0 limit of the exact NLO cross section would not have the product form H(z)·J(E) claimed implicitly. The entire resummation, the +9% enhancement in Table I, and the 1.5 GeV shift rest on this product form. Because the NLO hard functions H^(1) that appear in Eq. (6) are announced but not displayed (\"We have calculated... analytically\"), no reader can independently test the product form from the manuscript.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the ttbar invariant-mass distribution near the 2mt threshold at the LHC. It derives a factorization formula, Eq. (6), of the form d\\hat{\\sigma}_NLP ∼ H(z,QT,Y) × J(E), where H is a hard function and J is the pNRQCD potential/soft/ultrasoft Green function, and uses it to resum (αs/β)^n and ln^nβ corrections to next-to-leading power. The new NLO hard functions are announced but not displayed. The resummed result is matched to fixed-order NLO and NNLO predictions via Eq. (8), yielding NNLO+NLP = 1.434 pb/GeV averaged over M_tt in [300,380] GeV, compared with NNLO = 1.319 pb/GeV and CMS = 1.664 ± 0.166 pb/GeV. The authors claim a 9% enhancement and estimate a 1.5 GeV upward shift in the extracted top-quark mass.","tokens_in":9580,"tokens_out":4725,"duration_ms":51115,"significance":"If the factorization and the claimed absence of a soft function are correct, this work identifies a genuinely missing class of higher-order non-relativistic corrections that could resolve part of the longstanding discrepancy in the ttbar invariant-mass distribution, with direct consequences for top-quark mass extraction. The paper's numerical validation in Fig. 2—the β→0 limit reproducing the exact NLO correction in the [340,380] GeV window—is a strong check of the core mechanism, and Fig. 3 showing the fixed-order expansion diverging near threshold while the resummed result remains finite is compelling. The careful matching procedure and the use of external fixed-order results (MCFM, refs. [6,7]) add credibility. The main weakness is that the load-bearing 'no soft function' assertion in Eq. (6) is stated rather than derived, and the new NLO hard functions, which are essential for the product form and for the resummation, are not shown.","major_comments":[{"comment":"The assertion 'there is no soft function in our factorization formula' is not demonstrated for soft exchanges between the initial-state partons and the ttbar pair. The mode list in Eq. (5) includes soft modes with k∼Mtβ, and the potential function Jα(E) is defined to describe only exchanges between the top and anti-top quarks. A soft gluon of this scaling exchanged between an eikonal initial-state parton and the color-singlet or color-octet ttbar pair would have the same parametric scaling and is not shown to be power-suppressed at NLP. If such a contribution exists, Eq. (6) would need a soft function with a convolution entangling z and E, breaking the product form H×J on which the resummation, the +9% enhancement in Table I, and the 1.5 GeV mass shift all rest. The paper should either provide a power-counting derivation for why such exchanges vanish at this order, or include the soft function and redo the analysis.","section":"Section II, after Eq. (6)"},{"comment":"The text states 'We have calculated the NLO corrections to the hard functions analytically' and that they involve singular distributions in 1−z, QT, and Y, but the results are not displayed anywhere in the manuscript or in an ancillary file. Since these hard functions are the new ingredient required to implement the factorization at NLP, and since they are the only way a reader can verify the cancellation of infrared divergences and the validity of the product form, they should be provided explicitly (for example, in an appendix or as a supplementary file). Without them, the central derivation of Eq. (6) is not independently checkable.","section":"Section II, NLO hard functions"},{"comment":"The quoted central 9% enhancement and the resulting 1.5 GeV shift in mt depend entirely on the no-soft-function assumption. The uncertainty band quoted for NNLO+NLP (+0.014/−0.060 pb/GeV) accounts only for renormalization/factorization scale variation and does not include any uncertainty from the possible omitted soft contributions. At minimum, the authors should estimate the magnitude of a potential NLO soft contribution, for instance by inserting a model soft function and observing the change in the integrated cross section, and state how sensitive the 1.5 GeV shift is to this assumption.","section":"Section III, Table I and the 1.5 GeV estimate"}],"minor_comments":[{"comment":"The notation 'dσ(n)nLO' in Eq. (8) is confusing; the authors should clarify that the subtraction terms are the fixed-order expansions of the NLP resummed result to the same order as the corresponding fixed-order NLO and NNLO predictions, and align the labels with the 'niLO' notation defined in Section III.","section":"Section II, Eq. (8) and surrounding text"},{"comment":"The sentence explaining that ultrasoft gluon exchanges among initial-state partons and the ttbar pair do not contribute at NLP would benefit from a one-line power-counting justification, since the mode list in Eq. (5) contains both soft and ultrasoft modes and their different treatment is not obvious.","section":"Section II, after Eq. (6)"},{"comment":"The caption appears to read 'nLO LO' which is likely a typo for 'n0LO' and 'nLO'; please ensure the labels in the figure match the niLO notation defined in the text.","section":"Figure 3 caption"}],"recommendation":"major_revision","confidential_remarks":"This is a potentially important Letter whose central claim depends on the absence of a soft function in the NLP factorization. The numerical validation in Fig. 2 is encouraging, but the manuscript as written does not supply enough detail for a referee to verify the factorization of the NLO hard functions or the power counting that excludes soft exchanges with initial-state partons. I would be willing to accept a revised version that provides the NLO hard functions (in an appendix or supplementary material) and a rigorous argument for the missing soft function, or that includes such a soft function if it turns out to be present. The paper is well within the scope of the journal and the topic is timely."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nQuick take: this is a real contribution, not a housekeeping paper. The authors keep the full partonic sqrt(s) instead of forcing sqrt(s)->M_ttbar, which lets them absorb extra hard radiation into the hard function, and they carry out the resummation at NLP in the pNRQCD framework. That is a genuinely different setup from the earlier soft-threshold resummations. The NLO hard functions with full kinematic dependence are new (at least as stated), and the beta->0 limit of the EFT reproduces the exact NLO correction in the 340-380 GeV window (Fig. 2). The fixed-order expansion diverges approaching threshold; the resummed curve is finite and allows M_ttbar below 2mt via bound-state/width effects. Quantitatively they find a ~9% enhancement over NNLO and estimate a 1.5 GeV upward shift of the extracted top mass. The estimate is clearly labeled as such. That is honest and useful.\n\nCriticisms, in proportion.\n\nThe biggest weakness is presentation: the new hard functions are announced but not displayed, and no code or data are provided to regenerate Table I. As a letter, some omissions are expected, but for a result whose main selling point is a new NLO ingredient, at least the analytic expressions (or a link) should be available. Without them, the beta->0 validation is the only independent check, and that check is limited.\n\nThe second issue is the no-soft-function claim. After Eq. (6) the authors state that at NLP there is no soft function because ultrasoft exchanges among initial-state partons and the pair do not contribute. The stress-test note worries that the mode list includes soft modes at k~M beta, and that a soft gluon exchanged between an eikonal initial parton and the pair might not be power-suppressed. On reading the text, the argument is not fully spelled out. The claim may be correct — in the pNRQCD counting, such a contribution could be beyond NLP — but the paper does not demonstrate it. Since the entire +9% hangs on the product form H x J, this deserves a clear explanation in the paper or in a reply to referees.\n\nMinor: the 1.5 GeV shift is a rough estimate, not a full fit; the scale uncertainty on the shift is not given. The central scale choice HT/4 is fine, but the error band in Table I is asymmetric and driven by the resummed part.\n\nOverall, the core mechanism and numeric results are plausible. The paper should not be desk-rejected; it should go to referees, with a request for the hard functions and a precise statement on the absence of a soft function.\n\nBest.","headline":"Worth a serious referee: a genuinely new NLP threshold resummation for the ttbar invariant-mass spectrum, with real validation, but the unpublished hard functions and the asserted no-soft-function form need to be addressed.","tokens_in":10181,"tokens_out":3066,"would_cite":true,"duration_ms":34544,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["14.65.Ha","12.38.Bx","12.38.Cy"],"model":"deepseek-v4-flash","headline":"Resumming Coulomb and soft gluon exchanges near the top-pair threshold raises the NNLO differential cross section by about 9% and shifts the fitted top-quark mass by about 1.5 GeV toward the world average.","keywords":["top-quark pair production","invariant mass distribution","top-quark mass determination","threshold resummation","next-to-leading power","pNRQCD","Coulomb gluon resummation","NNLO matching"],"falsifier":"Compute the complete N3LO fixed-order differential cross section in the 300–380 GeV bin and compare it with the claimed $\\sim$9% enhancement; if the full third-order correction is much smaller, then the NLP resummation overestimates the missing contribution. Alternatively, repeat the kinematic top-mass fit using the NNLO+NLP prediction with the full experimental covariance of the measured distribution; if the fitted mass moves below the NLO-based value rather than up toward the world average, the claimed shift is refuted.","tokens_in":9085,"feed_emoji":"⚛️","tokens_out":10522,"duration_ms":101971,"temperature":0.7,"pith_summary":"The paper aims to close a persistent gap between theory and measurement for the top-antitop invariant-mass distribution just above twice the top mass, a region that strongly influences top-quark mass fits from kinematic distributions. It shows that higher-order non-relativistic corrections, which scale as powers of $\\alpha_s/\\beta$ with $\\beta$ the relative top velocity, are large near threshold and absent from current NNLO and NNLL$'$ predictions. By factorizing the cross section into a hard function times a potential function and resumming these Coulomb and soft-gluon effects to all orders at next-to-leading power, the paper obtains a finite threshold prediction that can be matched to NLO and NNLO fixed-order results. The matched NNLO+NLP prediction is about 9% higher than NNLO alone in the 300–380 GeV range and agrees much better with the measured distribution, implying an upward shift of about 1.5 GeV in the top-quark mass extracted from such fits, closer to the world average.","feed_headline":"Missing threshold corrections raise fitted top mass by 1.5 GeV","feed_subtitle":"A new all-order resummation of gluon exchanges makes theory match measured top-pair rates near threshold.","key_machinery":"The load-bearing object is the pNRQCD-based factorization formula $d\\hat\\sigma_{\\rm NLP}\\sim H\\times J$ near threshold. The hard functions $H_{ij,\\alpha}$ are Wilson coefficients of the effective theory describing exchanges and emissions of gluons with momenta of order $M_{t\\bar t}$, evaluated at NLO with full kinematic dependence so that the fixed-order matching is consistent with $H_T$-dependent scales. The potential functions $J_\\alpha(E)$ are related to the imaginary part of the pNRQCD Green function of the $t\\bar t$ pair at the origin and encode all potential, soft, and ultrasoft interactions; the top-quark width enters through the replacement $E\\to E+i\\Gamma_t$. The mode separation of Eq. (5), into hard, potential, soft, and ultrasoft momentum regions, is what lets the small-$\\beta$ singularities be identified and resummed, and the claimed absence of collinear modes and of a soft function at NLP is what fixes the simple $H\\times J$ form.","core_discovery":"The central discovery is that previously missing non-relativistic higher-order corrections constitute the dominant missing contribution in the threshold region. The paper derives a next-to-leading-power factorization $d\\hat\\sigma_{\\rm NLP}\\sim \\sum_\\alpha c_{ij,\\alpha}\\,H_{ij,\\alpha}\\,J_\\alpha$, in which the hard function $H$ contains all hard gluon exchanges and emissions, with explicit dependence on $z$, $M_{t\\bar t}$, $Q_T$, $Y$, and the renormalization and factorization scales, and the potential function $J_\\alpha(E)$ with $E=M_{t\\bar t}-2m_t$ resums Coulomb, soft, and ultrasoft exchanges between the slowly moving top and antitop. At this order the formula has no collinear modes and no soft function. Matching the resummed expression to NLO and NNLO results, the paper finds the threshold cross section is enhanced by about 9% in the 300–380 GeV window, making its central prediction $1.434\\,\\text{pb/GeV}$ compared with the NNLO value $1.319\\,\\text{pb/GeV}$ and a measured value near $1.664\\,\\text{pb/GeV}$; it estimates that the corresponding top-mass shift would be about $+1.5\\,\\text{GeV}$.","pith_inferences":["If the 9% enhancement is correct, analogous threshold resummations should improve other heavy-quark observables whose fixed-order predictions fall short of measured rates, and the same shift mechanism would appear in those fits.","The claim that no soft function exists at next-to-leading power is a power-counting statement; computing the first contribution from initial-state ultrasoft exchange would show whether the enhancement is stable or receives a large correction at the next order.","A direct test is to extend the factorization to the normalized $M_{t\\bar t}$ distribution, where scale and PDF uncertainties partially cancel; the predicted 1.5 GeV shift should become sharper and can be compared with a future higher-statistics measurement."],"forward_implications":["Near $M_{t\\bar t}=2m_t$, the resummed result stays finite where fixed-order expansions diverge: the n3LO curve goes to $+\\infty$ and the n4LO curve to $-\\infty$ as $\\beta\\to 0$.","The integrated cross section in the 300–380 GeV bin is insensitive to the top-quark width, even though the shape below $2m_t$ depends on it.","Using the NNLO+NLP prediction in the mass fit shifts $m_t$ by about $+1.5\\,\\text{GeV}$ relative to NLO-based fits, moving the kinematic value toward the direct-measurement world average.","The factorization can be combined with existing NNLL$'$ soft-gluon resummation and electroweak corrections to produce a precision prediction across the full phase space.","The same formalism extends to top-quark-pair-plus-jet production, another channel used for top-mass extraction."],"supporting_citations":[{"why":"The 13 TeV measurement of the $M_{t\\bar t}$ distribution whose excess over fixed-order theory motivates the paper.","marker":"[2]"},{"why":"The NNLO differential prediction that serves as the fixed-order baseline being matched and enhanced.","marker":"[6]"},{"why":"Earlier NNLL$'$ soft-gluon resummation representing the previous state of the art that misses the non-relativistic corrections.","marker":"[8]"},{"why":"The NNLO+NNLL$'$ prediction used in the comparison to data and as the state of the art that the paper extends.","marker":"[9]"},{"why":"The kinematic extraction of $m_t$ from top-pair distributions whose low central value the paper's correction shifts upward.","marker":"[15]"},{"why":"Earlier $H\\times J$ factorization for heavy-quark pair production that the paper generalizes with full kinematic scale dependence and subleading-$\\beta$ prefactors.","marker":"[19]"},{"why":"A similar potential-function resummation framework used as a comparison point and source of Green-function expressions.","marker":"[20]"},{"why":"The pNRQCD effective field theory that justifies the hard, potential, soft, and ultrasoft mode separation.","marker":"[22]"},{"why":"The next-to-leading-power pNRQCD Green function with finite top width used to build the potential functions $J_\\alpha$.","marker":"[24]"},{"why":"The Monte Carlo program used to compute exact NLO corrections for validating the small-$\\beta$ approximation and for the NLO matching.","marker":"[30]"}],"fun_headline_variants":["Missing threshold corrections lift top-quark mass by 1.5 GeV","Resummed gluon exchanges raise fitted top mass by 1.5 GeV","Threshold top-pair rate enhanced 9% by all-order gluon resummation","NLP resummation brings fitted top mass closer to direct value"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation hinges on the mode separation of Eq. (5) being complete at next-to-leading power: all hard radiation must be fully absorbed into $H$, all pair interactions into $J$, and no initial-state ultrasoft interactions may contribute at this order; if any omitted mode contributes, the resummed correction is misassigned and the 1.5 GeV shift would not follow.","fun_headline_variants_meta":{"raw":{"variants":["Missing threshold corrections lift top-quark mass by 1.5 GeV","Resummed gluon exchanges raise fitted top mass by 1.5 GeV","Threshold top-pair rate enhanced 9% by all-order gluon resummation","NLP resummation brings fitted top mass closer to direct value"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000638,"raw_usage":{"total_tokens":2951,"prompt_tokens":966,"completion_tokens":1985,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":582,"completion_tokens_details":{"reasoning_tokens":1902}},"tokens_in":582,"tokens_out":1985,"duration_ms":15866,"temperature":1.0,"reasoning_tokens":1902,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:51:58.763021+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the complete N3LO fixed-order differential cross section in the 300–380 GeV bin and compare it with the claimed $\\sim$9% enhancement; if the full third-order correction is much smaller, then the NLP resummation overestimates the missing contribution. Alternatively, repeat the kinematic top-mass fit using the NNLO+NLP prediction with the full experimental covariance of the measured distribution; if the fitted mass moves below the NLO-based value rather than up toward the world average, the claimed shift is refuted.","supporting_citations":[],"review_version":1}