{"id":"68ee0c0f-969f-4654-81e6-9a3edf154da7","arxiv_id":"2412.13261","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":1,"one_line_summary":"The paper computes the O(alpha_s) soft-limit constants that interpolate between massive and massless scheme resummations, raising the accuracy of heavy-flavour fragmentation predictions from NLL to NLL'.","lead":"This paper extends a framework that resums both quark-mass and soft-gluon logarithms in heavy-quark production to the NLL' level by computing the finite terms in the soft limit. The improved predictions for charm fragmentation in electron-positron collisions change by up to about 15 percent at large moments, suggesting sizeable effects beyond the previous NLL accuracy.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The NLL' matching formula hinges on an unshown subtraction of logarithmic terms from the Mellin transform; the limit checks are reassuring, but an independent re-derivation is needed to confirm δC0 (4.7).","rationale":"The paper is a careful analytic extension of an established NLL framework to NLL' accuracy. The key limit checks in Eqs. (4.8)-(4.9) are nontrivial and provide strong internal support for the construction: the massless limit and the large-N limit are both reproduced with the correct O(αS) constants. I therefore do not believe the central claim is wrong. However, the derivation of δC0 from Eqs. (4.4)-(4.5) involves a subtraction of logarithmic terms from the exponent that is not shown, and the paper itself notes that F1 and F2 are neglected. While F1 and F2 are formally subleading in N, the subtraction step is the load-bearing link between the momentum-space calculation and the Mellin-space formula. If that subtraction were incomplete, δC0 would be altered at intermediate y and the interpolation would be incorrect, even though the asymptotic limits would remain unchanged. This is a genuine soft spot that warrants an independent check, but it does not rise to rejection. The reader's verdict of CONDITIONAL with moderate confidence is appropriate; my concern overlaps with the reader's weakest assumption but sharpens it by focusing on the subtraction rather than on F1 and F2. The concrete numerical test described above would settle whether Eq. (4.7) is exactly the correct NLL' coefficient function.","tokens_in":17684,"tokens_out":17955,"duration_ms":159188,"concrete_test":"Compute the exact Mellin transform of Eq. (3.22) for several (N, ξ) pairs (e.g., N = 5, 10, 50 and ξ = 10^-2, 10^-4) using high-precision numerical integration of the plus-distribution definition ∫_0^1 (x^{N-1}-1) f(x) dx. Then subtract the O(αS) expansion of Exp[j+bar j] obtained from the explicit Appendix A integrals. Verify that the remainder is independent of log\\bar N and log ξ and equals C0^(1) + D0^(1) + δC0(\\bar N ξ) with δC0 from Eq. (4.7), to within numerical tolerance. Repeat at y ≈ 0.1, 1, 10 to cover both Θ regions of Eq. (4.7).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is Eq. (4.6) with δC0 in Eq. (4.7). The derivation of δC0 proceeds by taking Mellin moments of Eq. (3.22) using identity (4.2), obtaining Eqs. (4.4) and (4.5), then 'subtracting the logarithmic terms already present in the resummed exponent'. This subtraction step is not shown explicitly. If the O(αS) expansion of Exp[j+bar j] contains terms of the same functional form as the Li2 and log(1+y) pieces in Eq. (4.7), then the residual that defines δC0 would be different, and the smooth interpolation between the massless and massive schemes would fail at intermediate y = \\bar N ξ even though the limits y→0 and y→∞ (Eqs. 4.8-4.9) would still be reproduced. The neglect of F1 and F2 is not the weak point: they are O(1/N) and vanish as ξ→0, so they are beyond NLL' accuracy. The real fragility is that the subtraction is asserted rather than demonstrated, so an independent check is required before the NLL' claim can be considered secure. The unfixed scale μ in Eq. (4.6) is a secondary ambiguity affecting numerics, not the formal accuracy claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a formalism to extend soft-gluon resummation for heavy-quark fragmentation in e+e- collisions from NLL to NLL' accuracy by including the O(alpha_s) constant term delta C_0^(1) in the Mellin-space resummed cross section. The constant is obtained from a momentum-space calculation of the quasi-collinear virtual, collinear, and anti-collinear contributions, and is designed to interpolate between the massless and massive schemes. The authors check that the resulting expression reduces to the known massless result as xi -> 0 and to the massive-scheme constant as N -> infinity, and they illustrate the impact on the charm-ratio observable.","tokens_in":18015,"tokens_out":14790,"duration_ms":128460,"significance":"The significance of this work, if the derivation is confirmed, is that it provides a scheme-consistent NLL' resummation for heavy-flavour production, which is a necessary ingredient for precision extraction of fragmentation functions and for comparisons with charm data. The calculation is first-principles and is validated by reproducing both limiting schemes. The numerical study indicates that the new constant has a sizeable effect (up to ~15% at N=50), which strengthens the motivation for a full NNLL calculation. The paper is clearly written and the technical framework is well-structured.","major_comments":[{"comment":"The derivation of the key constant delta C_0^(1) is not shown. The paper presents the Mellin transforms (4.4) and (4.5) without derivation and states that Eq. (4.7) follows by subtracting 'the logarithmic terms already present in the resummed exponent' from them, but this subtraction is not carried out or documented. Since this step determines the NLL' constant, the central claim is not fully verifiable as written. Please provide an explicit calculation of the O(alpha_s) expansion of Exp[j+jbar], or demonstrate that the O(alpha_s) expansion of Eq. (4.6) reproduces the Mellin transform of Eq. (3.22) up to O(1/N) terms.","section":"Sec. 4.1, Eqs. (4.4)-(4.7)"}],"minor_comments":[{"comment":"The scale mu in Eq. (4.6) is not fixed by NLL' accuracy, and the two ad hoc choices used for the numerics lead to different predictions in Fig. 2. Please quantify the sensitivity of the charm ratio to this ambiguity, for example by showing both mu choices with their associated scale-variation bands.","section":"Sec. 4.2, Eq. (4.11)"},{"comment":"The assertion that the limit 1-x << xi2 recovers Eq. (3.3) is not demonstrated. Given the non-commutativity of the soft and mass limits discussed throughout the paper, a short derivation or at least a reference to the relevant calculation would be helpful.","section":"Sec. 3, Eq. (3.18)"},{"comment":"The functions F1 and F2 are stated to vanish as 1/N and as xi2 -> 0, but their numerical size at the finite-N values used in Fig. 2 (10-50) is not discussed. Please comment on whether the neglected terms could affect the numerical comparison, or state that they are expected to be negligible.","section":"Sec. 4.1, Eqs. (4.4)-(4.5)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a natural extension of the authors' earlier NLL work [18] and the calculation appears technically sound. The main issue is the unshown subtraction in Mellin space; I recommend major revision. The paper is likely within scope for EPJC. The authors should also be encouraged to make their numerical code or tables available if possible, though this is not required."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this paper is the real thing. It takes the NLL merging framework from ref [18] and extends it to NLL' by computing the O(alpha_s) finite soft-limit constants that interpolate between the massive and massless schemes. The new object, delta C0^(1) in Eq. (4.7), comes from a first-principles momentum-space calculation in Sec. 3, not from fitting to the target observable. The derivation passes the obvious consistency checks: the xi2->0 limit reproduces the massless scheme result, the N->infinity limit gives the massive-scheme constant, and the IR poles cancel. There is no circularity; the reliance on their own earlier work is legitimate and the validation against externally known results is real.\n\nWhat I like: the paper is honest about its limits. It explicitly flags that the scale mu in Eq. (4.6) is not fixed by NLL' accuracy, that F1 and F2 are neglected (with justification that they are O(1/N) and vanish as xi2->0), and that the 'nNLL' curves are an estimate, not a controlled approximation. The charm-ratio numerics are presented as exploratory, and the uncertainty band is not oversold.\n\nThe soft spot, and the one I'd push on as a referee, is the step from Eqs. (4.4)-(4.5) to Eq. (4.7). The text says the logarithmic terms already present in the resummed exponent are subtracted, but the subtraction is not shown. The stress-test worry is legitimate: if the O(alpha_s) expansion of Exp[j+jbar] contains pieces of the same functional form as the Li2 and log(1+y) terms, the residual defining delta C0 could be different in the intermediate y region even though the y->0 and y->infinity limits are preserved. I don't think that is actually what happens—the exponent's logarithmic terms are known from ref [18], and the limit checks are very reassuring—but this is exactly the kind of step that should be displayed explicitly. An expanded appendix or a short derivation would remove the last bit of doubt.\n\nBottom line: this is a paper for the resummation/QCD fragmentation community, and it deserves a serious referee. I would accept it with a request for the explicit subtraction step. Numerics should be treated as indicative, not final, given the mu ambiguity. I'd cite it in my own work on heavy-quark fragmentation.","headline":"Genuine NLL' extension of the mass/massless soft-gluon resummation framework, with a solid momentum-space derivation; one matching step is under-documented but the limit checks hold.","tokens_in":18533,"tokens_out":3851,"would_cite":true,"duration_ms":34067,"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":"Heavy-quark fragmentation in e+e− collisions can now be resummed at NLL′ accuracy through a single O(α_s) constant that interpolates between massive and massless schemes.","keywords":["soft-gluon resummation","heavy-quark fragmentation","NLL′ accuracy","massive and massless schemes","quasi-collinear limit","Mellin transform","charm ratio","e+e− annihilation"],"falsifier":"Take the exact O(α_S) massive-scheme cross section for e+e− → h + anti-h + X, compute its Mellin transform numerically without the large-N expansion, keep the functions F_1 and F_2 in Eqs. (4.4) and (4.5), and compare the resulting N → ∞ constant with $C_0^{{(1)}}$ + $δC_0^{{(1)}}$ + $D_0^{{(1)}}$ = C_F($π^{2}$/2 − 1); an O(α_S) mismatch beyond power-suppressed 1/N terms would show the claimed constant is incomplete.","tokens_in":17409,"feed_emoji":"⚛️","tokens_out":7320,"duration_ms":65529,"temperature":0.7,"pith_summary":"The paper claims that soft-gluon resummation for heavy-quark fragmentation in e+e− collisions can be promoted from NLL to NLL′ accuracy by adding one O(α_S) constant, $δC_0^{{(1)}}$, that interpolates between the massive and massless schemes. The previous NLL formalism left the constant parts of the two schemes unequal, so it could not claim NLL′ accuracy; this work computes the missing constant from an O(α_S) momentum-space calculation in the quasi-collinear limit and inserts it into the Mellin-space resummed expression. If correct, the result gives a single resummed formula that smoothly reduces to the massless resummation as the quark mass goes to zero and to the massive-scheme resummation as N goes to infinity. A sympathetic reader would care because the matched formula removes an ambiguity in existing heavy-flavour predictions and gives a sharper basis for comparing charm and bottom fragmentation data with experiment.","feed_headline":"New constant bridges heavy- and massless-quark resummations","feed_subtitle":"A single α_s correction reproduces both the massless limit and the massive constant in one formula.","key_machinery":"The load-bearing object is the correction constant $δC_0^{{(1)}}$ of Eq. (4.7), a function of y = \\bar{N}\\xi built from logarithms and dilogarithm functions that interpolates between the massive and massless constants. The mechanism that produces it is the anti-collinear (recoiling-quark) sector of the O(α_S) cross section, evaluated in the quasi-collinear limit, together with the formal Mellin identity (4.2), $x^{{N-1}}$ → -$e^{{\\Gamma(1-\\partial/\\partial \\log \\bar{N}}$)}\\Theta(1-x-1/\\bar{N}), which converts plus distributions with ξ_2 dependence into constants plus logarithmically enhanced terms. Subtracting the logarithmic terms from the resummed exponent isolates $δC_0^{{(1)}}$; the functions F_1 and F_2 left over are asserted to be power-suppressed in N and are neglected. This is what carries the argument from momentum space to the Mellin-space NLL′ formula.","core_discovery":"The central claim is Eq. (4.6): the NLL′-accurate resummed cross section has the form \\tilde{\\$\\sigma$}^{(NLL')} = \\tilde{E}^{(sub)}(1 + \\frac{\\alpha_S}{\\pi}($C_0^{{(1)}}$ + \\delta $C_0^{{(1)}}$))(1 + \\frac{\\alpha_S}{\\pi} $D_0^{{(1)}}$) \\, \\mathrm{Exp}[j + \\bar{j}], where j and \\bar{j} are the quasi-collinear jet exponents and $δC_0^{{(1)}}$, defined in Eq. (4.7), is a new O(α_S) function of y = \\bar{N}\\xi whose two limits close the scheme gap. In the massless limit $δC_0^{{(1)}}$ → 0, so the formula collapses to the standard massless resummation; as N → ∞, $C_0^{{(1)}}$ + $δC_0^{{(1)}}$ + $D_0^{{(1)}}$ → C_F(\\$pi^{2}$/2 - 1), the known massive-scheme constant. The paper establishes this by computing the O(α_S) cross section from virtual, soft, collinear and anti-collinear sectors, using the anti-collinear sector in the quasi-collinear limit and plus-distribution identities to make the ξ → 0 limit well defined, then Mellin transforming with the large-N identity (4.2) and subtracting the logarithms already in the exponent.","pith_inferences":["The same anti-collinear subtraction and plus-distribution identities should also fix the O(α_S) constant in other heavy-flavour observables whose soft limit is sensitive to the recoiling massive quark, such as heavy-quark jet substructure, where the scale choice μ^2 = max(q^2/\\bar{N}, m_c^2) already appears.","A natural next test is to compare the NLL′ formula against exact NNLO massive-scheme results at moderate x: the predictions should agree up to power-suppressed terms, which would confirm that δC_0^{(1)} is complete.","If a future full NNLL calculation becomes available, the scale μ left unfixed by NLL′ accuracy can be pinned down, and the dynamic-scale approximate-NNLL estimate in the paper can be checked against the exact running-coupling terms."],"forward_implications":["Heavy-quark fragmentation in e+e− annihilation can now be resummed at NLL′ accuracy while treating the mass exactly in the constant terms, not just in the logarithms.","The new formula reproduces the massless scheme smoothly as ξ → 0, with δC_0^{(1)} → 0 and the jet exponents reducing to their massless forms, and the massive scheme at large N with the constant C_F(π^2/2 − 1).","For the charm ratio R(N, q_A, q_B), the NLL′ correction changes the resummed prediction by up to about 15% at large N, while the scale-uncertainty bands are not significantly reduced.","The numerical size of these corrections indicates that a full NNLL resummation of mass and soft logarithms is needed to improve precision and address the current tendency of theory to overshoot charm data.","Because the framework is general, the same matching constant can be incorporated into future phenomenological studies of b- and c-quark fragmentation energy spectra."],"supporting_citations":[{"why":"Supplies the NLL resummation formalism for mass and soft logarithms that this paper extends to NLL′.","marker":"[18]"},{"why":"Provides the massless-scheme resummed result and the constants C_0^{(1)} and D_0^{(1)} that enter the matched formula.","marker":"[21]"},{"why":"Documents the non-commutativity of the soft and massless limits that the new constant must resolve.","marker":"[10]"},{"why":"Gives the massive-scheme resummed expression and the soft anomalous dimension used in the interpolation.","marker":"[25]"},{"why":"Provides the large-N Mellin-space identity used to extract the constant terms from the momentum-space calculation.","marker":"[34]"},{"why":"Supplies the same Mellin-space identity in a form adapted to soft-gluon resummation for hard processes.","marker":"[35]"},{"why":"Gives the charm-ratio comparison at NLL that the numerical section extends with NLL′ and approximate-NNLL predictions.","marker":"[30]"}],"fun_headline_variants":["A single α_s constant closes the heavy-quark gap","NLL' accuracy unifies massive and massless resummations","One α_s correction reproduces both mass limits","Bridging soft and mass limits with a single α_s"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole construction rides on the assumption that the Mellin-space identity used to extract constants, together with dropping the two leftover functions F_1 and F_2, captures every constant piece needed at order α_S; if those neglected terms are not truly power-suppressed, the new constant is wrong and the smooth massless-to-massive interpolation fails.","fun_headline_variants_meta":{"raw":{"variants":["A single α_s constant closes the heavy-quark gap","NLL' accuracy unifies massive and massless resummations","One α_s correction reproduces both mass limits","Bridging soft and mass limits with a single α_s"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00087,"raw_usage":{"total_tokens":3793,"prompt_tokens":994,"completion_tokens":2799,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":610,"completion_tokens_details":{"reasoning_tokens":2731}},"tokens_in":610,"tokens_out":2799,"duration_ms":20108,"temperature":1.0,"reasoning_tokens":2731,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T13:18:37.418794+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the exact O(α_S) massive-scheme cross section for e+e− → h + anti-h + X, compute its Mellin transform numerically without the large-N expansion, keep the functions F_1 and F_2 in Eqs. (4.4) and (4.5), and compare the resulting N → ∞ constant with $C_0^{{(1)}}$ + $δC_0^{{(1)}}$ + $D_0^{{(1)}}$ = C_F($π^{2}$/2 − 1); an O(α_S) mismatch beyond power-suppressed 1/N terms would show the claimed constant is incomplete.","supporting_citations":[{"cited_title":"Resummation of Threshold Corrections for Single-Particle Inclusive Cross Sections","cited_arxiv_id":"hep-ph/9806467","evidence_quote":"Gives the massive-scheme resummed expression and the soft anomalous dimension used in the interpolation."},{"cited_title":"Catani and L","cited_arxiv_id":null,"evidence_quote":"Provides the large-N Mellin-space identity used to extract the constant terms from the momentum-space calculation."}],"review_version":1}