{"id":"f0f6e7b0-0a7b-4aeb-92cc-173d477a9d83","arxiv_id":"2411.16634","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"The authors derive and solve the renormalization-group equation for the subleading shape function g17(ω,ω1;μ) in inclusive B decays, and simplify the RG evolution of a related exclusive soft function.","lead":"This paper computes the one-loop 'shape function' that controls a large theoretical uncertainty in rare B-meson decays, and solves its scale-evolution equation. It is the first step toward next-to-leading-order predictions for the resolved-photon contribution in Bbar→Xsγ and Bbar→Xsℓ+ℓ−, whose current errors are among the largest in these modes.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The factorized anomalous dimension in Eq. (3.3) hinges on dropping ΔZ17; the NLO verification is deferred and only shown in the massless limit, so the finite-charm analyticity assumption is unchecked.","rationale":"The manuscript is careful, honestly scoped, and contains real cross-checks: the numerical discretisation in Figs. 4 and 6 supports the analytic solutions, and the comparison with the independent anomalous dimension of [40] in Section 4 provides non-trivial confirmation. The single most load-bearing point is nevertheless the removal of ΔZ17. The one-loop calculation of Z17 is presented in detail, but the reason ΔZ17 is 'irrelevant' is an analyticity statement about the full jet functions; the only explicit example is tree level. The all-order validity is not derivable from the material shown, and the NLO check that would settle it is deferred to a forthcoming publication and performed in the massless limit. Since the physical process is dominated by the charm-quark loop, the smoothness of m_c→0 is essential but not demonstrated. This does not make the central one-loop result wrong; it makes the generalized factorization and the reduced evolution equation conditional on an unshown calculation. The reader's conditional verdict is appropriate.","tokens_in":32767,"tokens_out":9368,"duration_ms":93002,"concrete_test":"Compute the two-loop anti-hard-collinear jet function at finite charm mass (or at least repeat the calculation with m_c≠0) and explicitly evaluate the convolution of ΔZ17 in Eq. (2.33) with the cut NLO jet functions. Verify that all 1/ε singularities cancel without ΔZ17 and that the convolutions vanish. If finite-m_c jet functions have branch cuts in both half-planes, Eq. (3.3) and the solution (3.31) require modification; if the check passes, the factorized RG equation is validated beyond the one-loop example.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result Eq. (3.3) is obtained by dropping ΔZ17 from the one-loop Z-factor in Eq. (2.33). This is the load-bearing step: if ΔZ17 does not vanish after convolution with the relevant jet functions, the RG equation does not factorize and the closed-form solution (3.31) is not the evolution of g17. The authors argue by contour deformation, relying on the claim that the jet functions have all poles and branch cuts in one half-plane (Sections 2.4 and 4.1), but only the leading-order propagator (ω+i0)^{-1} is exhibited. For the physical Q1^c-Q7γ contribution the relevant two-loop anti-hard-collinear jet function involves a charm-quark loop; the authors state that they verified the 1/ε cancellation, but the calculation is deferred to [56] and performed in the massless limit m_c→m_u=0. If finite-m_c corrections move branch cuts into both half-planes, or if the m_c→0 limit is not smooth, the convolution with ΔZ17 need not vanish and Eq. (3.3) would need a non-factorized correction. The one-loop cut-diagram argument in Section 2.4 reduces but does not eliminate this concern, because the decisive NLO check is not shown.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives the one-loop renormalisation of the HQET operator O17 whose hadronic matrix element defines the subleading shape function g17(ω,ω1;μ) in the factorisation of the resolved-photon Q1^c–Q7γ interference in B̄→Xsγ and B̄→Xsℓ+ℓ−. The one-loop Z-factor is computed with explicit Feynman rules, modified plus-distributions, and dimensional regularisation. The authors find that, after removing an ostensibly non-factorising contribution ΔZ17, the anomalous dimension separates into two kernels acting on the two light-cone momentum variables ω and ω1, Eq. (3.3). They then solve the resulting integro-differential RG equation analytically by Mellin transformations, expressing the solution in terms of evolution functions involving Meijer-G and hypergeometric functions, Eq. (3.31), and provide numerical checks against direct discretisation. The same analyticity argument is applied to the amplitude-level soft function ΦG appearing in exclusive B̄d,s→γγ decays, for which a simplified 'reduced' RG equation is solved and compared with the full result of [40] through a double inverse moment, Eq. (4.28).","tokens_in":32929,"tokens_out":4005,"duration_ms":42740,"significance":"If the central factorisation result Eq. (3.3) holds, the paper supplies the first ingredient for a systematic next-to-leading-order analysis of one of the largest hadronic uncertainties in inclusive penguin B decays, and it provides a closed-form RG evolution valid beyond leading power. The calculation is performed with unusual transparency: the operator Feynman rules, the treatment of modified plus-distributions, and the analytic continuation prescriptions are given in enough detail to be checked independently. The Abelian part correctly reproduces the known evolution of the leading shape function [36], and the non-Abelian kernel is consistent with the independent calculation of [40] for the exclusive counterpart. The proposed simplification of multi-light-cone soft-function renormalisation, if confirmed, has broader applicability to other subleading-power and QED-generalised soft functions. The analytic solutions and the convergence properties of the evolved moments are valuable, and the numerical consistency checks in Figs. 4 and 6 support the algebra of the solution.","major_comments":[{"comment":"The factorised form of the anomalous dimension, Eq. (3.3), is obtained by dropping ΔZ17 from the one-loop Z-factor. The paper's justification is that ΔZ17 vanishes after convolution with jet functions because all jet-function singularities lie in one half-plane, but the only explicit example given is the leading-order propagator (ω+i0)^{-1}. The decisive next-to-leading-order check is stated to have been performed but is deferred to a forthcoming publication [56], and it is carried out in the massless limit mc→mu=0. This is a load-bearing point: if the finite-charm jet function develops branch cuts in both half-planes, or if the mc→0 limit is not smooth, Eq. (3.3) would need a non-factorised correction and the closed-form solution (3.31) would not describe the evolution of g17. The paper itself flags this limitation in Section 2.4 ('we have explicitly checked... in the massless case') and in the Conclusions ('will be discussed in a forthcoming publication'). The manuscript should either include the NLO cancellation check, or clearly state the central result as conditional on the analyticity assumption and on the smoothness of the massless limit.","section":"Section 2.4, Eqs. (2.33), (3.3), (3.31)"},{"comment":"The same analyticity argument is used to replace ΔΓG by a local imaginary term, Eq. (4.8), and to solve the 'reduced' RG equation of ΦG, Eq. (4.18). The identification H+ = H− on the space of jet-function test functions is plausible and consistent with the leading-order analytic structure, but the paper does not provide a proof for all orders or for the finite-charm case; the statement is supported only by the leading-order propagator example and by the assertion that the NLO 1/ε cancellation has been checked in the massless limit. The cross-check in Eq. (4.28) compares a single double inverse moment of the reduced-evolution result with the full result of [40], which is a necessary but weak test of the claimed equivalence of the evolution kernels. Since the simplified RG solution is a central by-product of the paper, the reduction should either be proven more rigorously or the scope of the claim should be restricted accordingly.","section":"Section 4.1, Eqs. (4.8), (4.9), (4.18)"}],"minor_comments":[{"comment":"The text twice refers to 'Meier-G functions'; the correct name is Meijer-G functions.","section":"Section 4.2, after Eq. (4.21)"},{"comment":"The agreement between the analytic evolution and the discretised RG equation is described only as 'good agreement'; adding a residual plot or a quantitative error measure would make the verification more convincing.","section":"Figs. 4 and 6"},{"comment":"The definitions of the modified plus-distributions [⋯]⊕ and [⋯]⊖ are essential for the central calculation, but the notation is dense; a short worked example showing how a test function is evaluated at ω or at ±ω would improve readability.","section":"Eq. (2.12) and surrounding text"},{"comment":"The phenomenological discussion explicitly leaves a serious numerical estimate to future work; this should be stated already in the introduction or section heading so that the reader does not expect a quantitative estimate of the scale ambiguity reduction.","section":"Section 3.3"}],"recommendation":"major_revision","confidential_remarks":"The paper is methodologically sound and the one-loop calculation is presented with a high degree of transparency. The main concern is that the central factorisation claim depends on an analyticity property whose decisive NLO verification is relegated to a forthcoming publication [56] and is only demonstrated in the massless limit. This is exactly the kind of situation where a referee should ask for the missing check to be included or for the claims to be explicitly conditioned on it. The manuscript is otherwise well within the scope of a high-energy physics journal and the results, if confirmed, would be a useful contribution to subleading-power SCET phenomenology."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my read. The paper supplies the missing soft-scale ingredient for NLO corrections to the resolved-photon Q1c–Q7gamma interference in B->Xs gamma and B->Xs l+l-. That alone makes it worth attention. The actual calculation looks careful: explicit Feynman rules, distribution identities, regulator prescriptions, and the Abelian part reproduces the known leading shape-function evolution. The non-Abelian piece matches the independent exclusive soft-function calculation in [40]. The analytic momentum-space solution of the RG equation, with the factorized evolution functions, is new and technically neat. The 'reduced' RG equation for the exclusive Phi_G is a useful simplification.\n\nThe soft spot is the treatment of Delta Z17. The factorized anomalous dimension in Eq. (3.3) depends on dropping this term, and the justification is a contour-deformation argument: the jet functions' singularities are claimed to lie in one half-plane. The paper shows this explicitly only for the leading-order propagator; the NLO check is deferred to a forthcoming publication and is computed in the massless limit m_c -> m_u = 0. If finite charm mass shifts branch cuts, or if the massless limit is not smooth, the factorization would need modification. This is a real caveat, and it is load-bearing. But the authors flag it themselves, they have performed the NLO cancellation check in the massless case, and the one-loop cut-diagram argument is plausible. So it is a reason to demand the follow-up calculation, not a reason to distrust the central one-loop result.\n\nThe paper is honest about its scope: no serious numerical estimate of evolution effects is given, and the phenomenological section is illustrative. I don't see circularity. The citation pattern is fine; the overlap with [47] is natural given shared authorship, and they generalize rather than repackage that work.\n\nI'd send this to peer review. A good referee should verify the one-loop integrals and press for the NLO check (or at least a clear statement of the massless-limit caveat) before fully endorsing the all-order claim. I'd cite it for the one-loop anomalous dimension, and I'd happily discuss it in a reading group focused on soft functions at subleading power.","headline":"Careful one-loop calculation of the g17 anomalous dimension with a plausible but not yet fully demonstrated argument for dropping the non-factorizing term.","tokens_in":33609,"tokens_out":6094,"would_cite":true,"duration_ms":50977,"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 derives the one-loop anomalous dimension of the subleading B-meson shape function $g_{17}(\\omega,\\omega_1;\\mu)$, shows that the naively mixed terms are irrelevant in physical convolutions, and solves the resulting factorised RG…","keywords":["subleading power factorization","shape function g17","soft-collinear effective theory","heavy-quark effective theory","resolved photon contributions","B to X_s gamma decay","B to X_s l+l- decay","renormalisation group evolution"],"falsifier":"Evaluate the two-loop ($O(\\alpha_s^2)$) anti-hard-collinear jet function in the massless-quark limit and compute its convolution with $\\Delta Z_{17}$ from Eq. (2.33); if the integral over $\\omega$ and $\\omega_1$ does not vanish on the physical jet-function space, the factorised anomalous dimension in Eq. (3.3) is not the correct RG kernel for $g_{17}$. For the exclusive case, one can instead test whether the full solution of the anomalous dimension in Ref. [40] and the 'reduced' solution in Eq. (4.18) agree for a low-scale model of $\\Phi_G$ that is allowed negative support.","tokens_in":32445,"feed_emoji":"💥","tokens_out":15886,"duration_ms":123118,"temperature":0.7,"pith_summary":"This paper establishes that the scale dependence of the subleading shape function $g_{17}(\\omega,\\omega_1;\\mu)$ --- the non-perturbative soft function controlling the resolved-photon $Q_1^c$--$Q_{7\\gamma}$ interference in $\\bar{B}\\to X_s\\gamma$ and $\\bar{B}\\to X_s\\ell^+\\ell^-$ --- is governed by a one-loop anomalous dimension that factorises into two independent kernels, one acting on each light-cone momentum variable. Because of this factorisation, the renormalisation-group equation separates and can be solved analytically in momentum space, leading to the closed-form solution in Eq. (3.31). The authors argue that the terms in the $Z$-factor that appear to mix the two collinear sectors are irrelevant: they vanish when convolved with jet functions whose singularities all lie on one side of the real axis, and they are eliminated by the correct cut (closed-time-path) definition of the soft operator. The same reasoning simplifies the renormalisation of an amplitude-level soft function $\\Phi_G$ relevant to exclusive $\\bar{B}_{d,s}\\to\\gamma\\gamma$, yielding a 'reduced' RG equation that is solved explicitly. This is the first ingredient needed for an NLO RG-improved calculation of the largest resolved-photon uncertainty in the inclusive penguin modes.","feed_headline":"Subleading B-decay shape function evolution solved in closed form","feed_subtitle":"Resolves the largest soft-scale uncertainty in B to X_s gamma and B to X_s l+l- at subleading power.","key_machinery":"The central object is the subleading shape function $g_{17}(\\omega,\\omega_1;\\mu)$, the Fourier-transformed forward matrix element of an HQET operator $O_{17}$ that contains two heavy-quark fields separated along the $n$ light-cone and a soft gluon field-strength smeared along the $\\bar n$ light-cone. The machinery carrying the argument is the one-loop $Z$-factor and its anomalous dimension, computed both with ordinary time-ordered Feynman rules and with cut (closed-time-path) rules that distinguish fields in the amplitude from fields in the complex-conjugate amplitude. The decisive mechanism is analyticity: the relevant jet functions have their singularities on one side of the real axis ($\\omega+i0$ for the inclusive cut process, $\\omega-i0$ for the exclusive time-ordered process), so the mixed convolution terms $\\Delta Z_{17}$ and $\\Delta\\Gamma_G$ can be removed by contour deformation. Once removed, the anomalous dimension decouples into separate kernels for $\\omega$ and $\\omega_1$; each sector is then solved by a Mellin transformation, and the momentum-space evolution functions are expressed with Meijer-G functions.","core_discovery":"The central claim is that the one-loop $\\overline{\\rm MS}$ anomalous dimension of the operator $O_{17}$, after the contribution $\\Delta Z_{17}$ in Eq. (2.33) is dropped, takes the factorised form\n$$\\gamma_{17}(\\omega,\\omega_1,\\omega',\\omega'_1;\\mu)=\\frac{\\alpha_s}{\\pi}\\left\\{C_F\\,\\delta(\\omega_1-\\omega'_1)\\,\\gamma_n(\\omega,\\omega';\\mu)+\\frac{C_A}{2}\\,\\delta(\\omega-\\omega')\\,\\gamma_{\\bar n}(\\omega_1,\\omega'_1;\\mu)\\right\\}.$$\nThe Abelian kernel $\\gamma_n$ acts only on the soft variable associated with the $n$ direction and is identical in structure to the leading shape-function kernel; the non-Abelian kernel $\\gamma_{\\bar n}$ acts only on the variable associated with the $\\bar n$ direction. The mixed terms $\\Delta Z_{17}$ are not needed for physical quantities because the time-ordered jet functions have poles and branch cuts only in the lower half-plane, and the cut prescription in Eq. (2.37) enforces $t>0$, so their convolutions vanish. For the exclusive amplitude-level soft function $\\Phi_G$, an analogous set of mixed terms in the anomalous dimension becomes irrelevant after convolution with jet functions whose singularities lie in the upper half-plane; there the identity $H_+=H_-$ holds on the relevant test-function space, reducing the anomaly to the factorised form in Eq. (4.9). The paper checks the cancellation of all $1/\\varepsilon$ singularities between hard, jet and soft loops at NLO in the massless-quark limit, using a two-loop anti-hard-collinear jet-function calculation that is reported in a forthcoming publication.","pith_inferences":["If the one-sided analyticity property holds beyond one loop, the irrelevance of mixed two-light-cone operator-mixing terms is likely a general structural property of multi-light-cone soft functions; a testable extension would be to apply the same reduction to soft-quark functions in Drell-Yan or $gg\\to h$ at next-to-leading power.","Appendix A's requirement that matrix elements of $O_{17}$ be integrated over the full real axis in $\\omega$ implies that practical models for $g_{17}$ with a hard cut-off in $\\omega$ are not consistent with renormalisation; model builders may need to work with full-line support and control the radiative tail by evolution.","The claim that reduced evolution introduces no complex phases in inverse moments is phenomenologically consequential for CP asymmetries in $\\bar{B}\\to X_s\\gamma$: it would mean strong phases in the resolved contribution come only from the initial soft function at $\\mu_0$, not from resummation. This could be tested once a complete model of $g_{17}$ at the low scale is available.","The identity $H_+=H_-$ on the jet-function test space suggests that an anomalous dimension of a soft operator is defined only modulo terms with one-sided analyticity; classifying such 'irrelevant' terms could simplify future subleading-power RG analyses."],"forward_implications":["The resolved-photon contribution from the $Q_1^c$--$Q_{7\\gamma}$ interference in $\\bar{B}\\to X_s\\gamma$ and $\\bar{B}\\to X_s\\ell^+\\ell^-$ can now be RG-improved at next-to-leading order, removing the large scale ambiguity of the leading-order result.","The scale evolution of $g_{17}$ is completely fixed by its initial condition at the low scale through Eq. (3.31); the function remains real, so no strong phases are generated by resummation.","For the exclusive $\\bar{B}_{d,s}\\to\\gamma\\gamma$ soft function $\\Phi_G$, the reduced RG equation factorises in $\\omega$ and $\\omega_1$, a positively-supported initial function stays positively supported, and inverse moments acquire no complex phases, simplifying the resummation of large logarithms.","The cancellation of $1/\\varepsilon$ singularities between hard, jet and soft loops at NLO, checked in the massless-quark limit with the two-loop $\\bar n$-jet function, supports the consistency of the underlying factorisation formula."],"supporting_citations":[{"why":"defines the shape function g17 and supplies the leading-order factorization formula (1.2) and the identity (3.35).","marker":"[16]"},{"why":"provides the leading-power shape-function renormalisation and the Abelian solution that the paper generalises.","marker":"[36]"},{"why":"introduces the modified plus-distributions H±, the Mellin-space solution method, and the concept of irrelevant anomalous-dimension terms.","marker":"[47]"},{"why":"derives the full anomalous dimension and RG solution for the amplitude-level soft function ΦG that the paper confirms and reduces.","marker":"[40]"},{"why":"establishes the non-local resolved-photon power corrections in which g17 enters.","marker":"[10]"},{"why":"supplies the model for h17(ω1;μ0) and moment bounds used in the numerical illustration.","marker":"[28]"},{"why":"quantifies the leading-order resolved-photon uncertainty and scale ambiguity that motivates the NLO analysis.","marker":"[29]"},{"why":"provides the refactorisation of the subleading Bbar→Xsγ formula used in this work.","marker":"[31]"}],"fun_headline_variants":["Subleading soft function RG evolution solved for B decays","Shape function g17 RG equation solved at subleading power","Taming the largest soft uncertainty in B to X_s gamma","Closed-form RG evolution for subleading B-decay shape function"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument hinges on the analytic behaviour of the jet functions: if at higher orders these functions develop singularities or branch cuts on both sides of the real axis, or do not fall off fast enough for contour deformation, the troublesome mixed terms in the anomalous dimension would not vanish, and the factorised RG equation would have to be modified.","fun_headline_variants_meta":{"raw":{"variants":["Subleading soft function RG evolution solved for B decays","Shape function g17 RG equation solved at subleading power","Taming the largest soft uncertainty in B to X_s gamma","Closed-form RG evolution for subleading B-decay shape function"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000135,"raw_usage":{"total_tokens":1642,"prompt_tokens":1192,"completion_tokens":450,"prompt_tokens_details":{"cached_tokens":1152},"prompt_cache_hit_tokens":1152,"prompt_cache_miss_tokens":40,"completion_tokens_details":{"reasoning_tokens":382}},"tokens_in":40,"tokens_out":450,"duration_ms":37034,"temperature":1.0,"reasoning_tokens":382,"cache_read_input_tokens":1152,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:53:49.993039+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the two-loop ($O(\\alpha_s^2)$) anti-hard-collinear jet function in the massless-quark limit and compute its convolution with $\\Delta Z_{17}$ from Eq. (2.33); if the integral over $\\omega$ and $\\omega_1$ does not vanish on the physical jet-function space, the factorised anomalous dimension in Eq. (3.3) is not the correct RG kernel for $g_{17}$. For the exclusive case, one can instead test whether the full solution of the anomalous dimension in Ref. [40] and the 'reduced' solution in Eq. (4.18) agree for a low-scale model of $\\Phi_G$ that is allowed negative support.","supporting_citations":[],"review_version":1}