{"id":"02f90483-1560-4ee5-91f4-6fe0f8990b16","arxiv_id":"2504.18018","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"The QCD B-meson shape function in the endpoint region equals the HQET shape function multiplied by a perturbative matching coefficient whose one-loop form is derived in this paper.","lead":"This paper derives a one-loop factorization formula that relates the B-meson shape function defined in full QCD to the shape function defined in heavy quark effective theory. The formula separates the heavy quark mass scale from the nonperturbative hadronic scale, which is a step toward lattice QCD calculations of the shape function and improved inclusive B decay predictions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The one-loop matching coefficient (Eq. 22) behind the central formula (Eq. 23) is asserted but not derived: the cancellation of plus distributions, Lambda-dependence, and residual-momentum dependence is not shown, and the tail Z_tail (Eq. 19) omits the perturbative HQET tail.","rationale":"The paper's goal is clear: a leading-power, one-loop connection S_QCD(x, mu) = [1 + (alpha_s C_F/2pi) C(mu, m_b)] S_HQET(omega, mu) at omega v_+ = x m_B v_+ - m_b v_+ in the peak region, plus the tail coefficient Z_tail (Eq. 19). I read the derivation in good faith. The tree-level matching (Eqs. 8-9) is internally consistent. The local-vertex combination is dimensionally consistent with the single-log term: -S_QCD^(c) + S_HQET^(f) = -2 + ln(m_b^3/((-2v.k)^2 mu)) - ln(mu^2/(-2v.k)^2) = -2 + 3 ln(m_b/mu), matching the -3 ln(mu/m_b) single-log coefficient in Eq. (22). The tail-region sum of Eqs. (13b)+(14b) arithmetically reproduces Eq. (19) as the free-quark QCD tail. So the presented arithmetic is not obviously wrong, and the paper honestly acknowledges the uncontrolled transition region (Sec. IV) and the leading-power/one-loop scope (Sec. V). The load-bearing weakness is that the passage from Eqs. (13a)-(18) to Eq. (22), i.e., the cancellation of plus distributions, the Lambda cutoff, the v.k regulators, and all k_+-dependence, is asserted rather than shown; no code, data, or independent check is provided. This is the single most important point because if that combination retains any residual IR or regulator dependence, the central product formula (Eq. 23) fails. The second concern is the tail branch: the paper's own statement that S_HQET has a radiative tail at large |omega| implies that at x ~ 0 the convolution form of the factorization receives a non-negligible S_HQET^(1)(omega = x m_B - m_b) contribution of order alpha_s/m_b, so Eq. (19) as written is the bare QCD tail rather than a two-sided matching coefficient. The reader's CONDITIONAL verdict remains appropriate: the central formula is plausible and physically expected (Sudakov double log plus standard single logs), but the verification gap and the tail/transition incompleteness justify conditional acceptance rather than acceptance or rejection on the present evidence. My read does not change the verdict; it identifies the specific step that an independent re-derivation should settle.","tokens_in":15673,"tokens_out":32690,"duration_ms":316359,"concrete_test":"Independently recombine Eqs. (13a)-(18) in the matching equation (11) in d = 4 - 2epsilon with MS renormalization, keeping generic k_+, Lambda, and v.k: evaluate the delta-function subtraction integrals after the substitution omega v_+ = x m_B v_+ - m_b v_+, and verify term-by-term that the sum equals Eq. (22) with coefficient (1/2)ln^2(mu^2/m_b^2) - (3/2)ln(mu^2/m_b^2) + pi^2/12 - 2, with no residual dependence on k_+, Lambda, or v.k. If the cancellation fails, Eq. (23) is invalid. Separately, evaluate S_HQET^(1)(omega) from Eqs. (16)-(18) at omega v_+ = x m_B v_+ - m_b v_+ for x -> 0 and check whether S_QCD^(1)|tail - S_HQET^(1)|tail reproduces Z_tail^(1) of Eq. (19); if the HQET tail does not vanish, Eq. (19) requires the missing subtraction.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim, Eq. (23), follows only if the combination of the six one-loop expressions in Eqs. (13a)-(18) via Eq. (11) yields the local, k_+-independent coefficient in Eq. (22). This is the load-bearing step, and it is not demonstrated in the paper. After the delta functions impose omega v_+ = x m_B v_+ - m_b v_+ and omega v_+ = k_+, the delta-function subtraction pieces in Eqs. (13a), (14a), (16) and (17) combine into integrals whose limits involve k_+ and m_b (e.g., an integral from 2k_+ to m_b v_+ + 2k_+ over a logarithmic integrand, plus box/local terms). Such integrals are, on their face, dependent on the external residual momentum k_+, the IR cutoff Lambda, and the regulators v.k; the paper states that these cancel and that 'all dependence on x and omega resides in the delta function' (Sec. III.D), but the intermediate algebra is absent. No code or independent numerical check is provided. A second, closely related gap: Eq. (19) sets Z_tail^(1) equal to the free-quark QCD tail (sum of Eqs. (13b) and (14b)), although the paper itself notes S_HQET has a radiative tail at large |omega|; at x ~ 0 the convolution would evaluate S_HQET at omega = x m_B - m_b ~ -m_b, where S_HQET^(1) is nonzero at order 1/m_b. A genuine two-sided matching would give Z_tail = S_QCD|tail - S_HQET|tail, so Eq. (19) is either incomplete or the tail branch of Eq. (5) is a pure QCD OPE that does not connect to S_HQET. Both points are concrete and checkable; the peak-region cancellation is the more load-bearing because Eq. (23) fails if it does not go through exactly.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives a factorization formula relating the B-meson shape function defined in QCD, S_QCD(x, μ), to the shape function defined in HQET, S_HQET(ω, μ). The claimed relation is region-dependent: in the peak region x ~ 1 - Λ_QCD/m_b the QCD shape function equals a multiplicative one-loop matching coefficient times the HQET shape function evaluated at ωv_+ = x m_B v_+ - m_b v_+, while in the tail region x ~ 0 it is given by a purely perturbative coefficient. The authors compute the one-loop corrections to both shape functions using free-quark external states, list the individual diagram contributions for the QCD sail, box, and local vertex graphs and their HQET counterparts, and state the resulting one-loop matching coefficients Z_peak^(1) in Eq. (22) and Z_tail^(1) in Eq. (19). They then use a phenomenological model for S_HQET to plot the QCD shape function and discuss applications to lattice QCD via two-step factorization.","tokens_in":16118,"tokens_out":10684,"duration_ms":100628,"significance":"If correct, Eq. (23) provides a simple one-loop conversion between the two standard definitions of the B-meson shape function, enabling resummation of logarithms of m_b/Λ_QCD and serving as the second step in the recently proposed lattice two-step factorization scheme. The individual one-loop graph results are presented in substantial detail, including explicit IR regulators v·k and the plus-distribution definition, and the claimed cancellation of the IR regulator between QCD and HQET is physically expected and plausible. The numerical section uses model parameters taken from prior work, so the derivation is not circular. The main weakness is that the step from the individual graph results to the central matching coefficient is not demonstrated, and the treatment of the tail region is not fully specified. These issues are potentially fixable but are load-bearing for the paper's central claim.","major_comments":[{"comment":"The central result, Eq. (22), is asserted as the outcome of inserting Eqs. (20) and (21) into Eq. (11), but the intermediate algebra is not shown. The paper claims that the plus distributions, the v·k dependence, the Lambda dependence, and the residual-momentum dependence all cancel, leaving a local coefficient, yet the presented graph results in Eqs. (13a)-(18) contain terms that, after the substitution ωv_+ = x m_B v_+ - m_b v_+, generate nontrivial integrals: for example, the delta-function subtraction in Eq. (13a) involves an integral from x m_B v_+ + k_+ to Lambda, and Eq. (16) involves an integral from 2 ω v_+ to Lambda over logarithmic integrands. I could not verify the claimed cancellation from the text. Because Eq. (23) depends entirely on this step, the authors should provide the full combination of the six one-loop expressions, at least in an appendix, or supply an independent numerical check of Eq. (22).","section":"Section III.D, Eq. (22)"},{"comment":"The tail matching coefficient is written as the sum of the QCD tail contributions in Eqs. (13b) and (14b), with no subtraction of the perturbative HQET tail. The paper itself notes in this section that S_HQET has a calculable radiative tail at large |ω|, and at x ~ 0 the peak-branch variable would correspond to ω ~ -m_b, where that tail is not negligible at one-loop order. A genuine two-sided matching would give Z_tail = S_QCD|tail - S_HQET|tail. As written, Eq. (19) is either incomplete or the tail branch of Eq. (5) is intended to be a pure QCD OPE that does not connect to S_HQET. The manuscript should state explicitly which interpretation is intended and, if a subtraction is required, include it and assess its numerical effect on Fig. 2.","section":"Section III.D, Eq. (19)"},{"comment":"The factorization in Eq. (5) assumes a clean leading-power separation between the peak region x ~ 1 - Lambda_QCD/m_b and the tail region x ~ 0, but the transition region is not quantified. The paper acknowledges in Sec. IV that the shaded intermediate region requires higher-order (1-x) corrections, yet it does not specify a power-counting criterion for where Eq. (23) ceases to be valid. Please define the criterion used for the separation (for example, the size of Lambda_QCD / (m_b (1-x))) and indicate the corresponding boundary in Fig. 2.","section":"Section IV and Fig. 2"}],"minor_comments":[{"comment":"The argument of Z_tail^(0) is written as (x, ω, μ), but the tail coefficient has no ω dependence; the notation should be corrected to Z_tail^(0)(x, μ).","section":"Section II.B, Eq. (7b)"},{"comment":"The sentence \"We now proceed to the matching in the peak region, characterized by momentum fractions x ~ 0\" should read \"x ~ 1\" (or \"x ~ 1 - λ\"), since the immediately preceding paragraph discusses the tail region with x ~ 0.","section":"Section III.D, first paragraph"},{"comment":"The logarithm in the local vertex result is ambiguous as printed; it should be written with explicit parentheses, e.g., ln( m_b^3 / ( (-2 v·k)^2 μ ) ).","section":"Eqs. (15a) and (A11)"},{"comment":"There is a typographical error: \"Ferimi motion\" should be \"Fermi motion\".","section":"Introduction"},{"comment":"\"Identifying momentums\" should be \"Identifying momenta\".","section":"Section II.B"},{"comment":"The statement that the tail of S_HQET can be determined perturbatively is not used anywhere in the paper; if the tail branch is intended as a pure QCD OPE, this sentence is misleading and should be removed or clarified.","section":"Section III.C"}],"recommendation":"major_revision","confidential_remarks":"The paper is timely and from a recognized group, and the individual one-loop computations appear careful. The main risk is whether the stated matching coefficient Eq. (22) is actually correct; the missing combination algebra is the key gap. The tail-region issue in Eq. (19) is also important to resolve. Both are fillable within the scope of a revision, so I recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague],\n\nYou should know two things about Wang et al. (arXiv:2504.18018). First, the one-loop matching coefficient in Eq. (22), the basis for the multiplicative formula in Eq. (23), is asserted rather than demonstrated. Second, the paper is still worth a careful read: the matching between the QCD and HQET B-meson shape functions is new, the individual diagram calculations are provided in an appendix, and the claimed cancellations are plausible. My verdict is conditional, matching the reader's.\n\nWhat is actually new: the tree-level and one-loop matching functions connecting the two shape functions. This is the shape-function analogue of the LCDA factorization of Refs. [36,37], and I did not find it in the cited literature. The derivation is standard matching technology applied to a new object, and the numerical section is honestly labeled as a model illustration. The matching coefficient does not depend on the fitted model parameters, so there is no circularity. The citation pattern looks fine; no priority problem visible.\n\nThe main soft spot is that Eq. (22) appears with no intermediate algebra. The paper states that the v·k, Lambda, and plus-distribution dependences cancel, and that all x and omega dependence goes into the delta function, but the reader never sees the combination of Eqs. (20)–(21). Since the whole result hinges on this cancellation, a referee needs to see it. I do not think the result is wrong—similar cancellations are standard in matching calculations—but \"plausible\" is not \"checked.\" The appendix gives the individual graphs in good detail, so the fix is a few lines of algebra or an ancillary file.\n\nThe tail-region formula in Eq. (19) is terse. The stress-test note worries that S_HQET has a radiative tail at omega ~ -m_b, so the matching should be two-sided. On reading the paper, that concern is partly misplaced: the tail branch of Eq. (5) does not send S_QCD through S_HQET; it directly assigns the perturbative QCD tail. The authors do mention that both tails are perturbative. Still, calling Z_tail a \"matching function\" invites the question, and the paper should clarify whether the tail branch is a matching or an OPE result.\n\nThe transition region between peak and tail is explicitly left to future work. For a paper whose selling point is a factorization formula, that is a real limitation, though not fatal at leading power.\n\nWho benefits: people working on inclusive |Vub|, lattice QCD for light-cone quantities, and SCET matchings. The paper deserves peer review, but the referee should require the cancellation algebra before publication. I would not desk-reject.","headline":"New one-loop matching between QCD and HQET B-meson shape functions, but the paper skips the load-bearing cancellation algebra and leaves the tail region fuzzy; worth refereeing, not desk-rejecting.","tokens_in":16649,"tokens_out":8149,"would_cite":true,"duration_ms":74061,"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":"The paper derives a one-loop factorization formula connecting the QCD and HQET definitions of the B-meson shape function, with a multiplicative peak-region matching coefficient and a perturbative tail function.","keywords":["B-meson shape function","heavy quark effective theory","factorization","inclusive B decays","|Vub|","lattice QCD","one-loop matching","lightcone distribution"],"falsifier":"A two-loop computation of the matching coefficient would settle whether the multiplicative form of Eq. (23) survives beyond one loop; alternatively, a lattice extraction of the quasishape function matched to $S_{\\rm QCD}$ could be compared with Eq. (23) across the peak-to-tail transition, where a growing discrepancy would show the region separation is not clean.","tokens_in":15442,"feed_emoji":"⚛️","tokens_out":7873,"duration_ms":67167,"temperature":0.7,"pith_summary":"This paper derives a factorization formula connecting the two standard definitions of the B-meson shape function: the QCD definition, which captures the b quark's lightcone momentum fraction and contains physics at both $m_b$ and $\\Lambda_{\\rm QCD}$, and the HQET definition, which describes the heavy quark's residual momentum in the infinite-mass limit. The central claim is that in the peak region, $x \\sim 1 - \\Lambda_{\\rm QCD}/m_b$, the QCD shape function equals the HQET shape function multiplied by a simple one-loop coefficient, with the residual momentum shifted by $m_b$. In the tail region the two are connected by a purely perturbative matching function. The result matters because the HQET shape function is a major uncertainty in inclusive $B$ decays and the $|V_{ub}|$ determination, while the QCD version is more naturally accessible to lattice computation; the formula provides a bridge between the two.","feed_headline":"One-loop formula connects the two B-meson shape function definitions","feed_subtitle":"It bridges perturbative heavy-quark physics and lattice-computable input for inclusive B decays and |Vub|.","key_machinery":"The load-bearing object is the region-separated factorization formula Eq. (5), which splits $S_{\\rm QCD}(x,\\mu)$ into a peak region $x\\sim 1-\\Lambda_{\\rm QCD}/m_b$, matched to $S_{\\rm HQET}(\\omega,\\mu)$ by a kernel $Z_{\\rm peak}(x,\\omega,\\mu)$, and a tail region $x\\sim 0$, described by a perturbative matching function $Z_{\\rm tail}(x,\\mu)$. The argument is carried by explicit one-loop calculations of three Feynman diagrams in each theory (heavy-quark sail, box, and local vertex). The cancellation of the modified plus distribution of Eq. (12) and the infrared regulator $v\\cdot k$ in the matching step is what produces the simple multiplicative peak-region coefficient, Eq. (22).","core_discovery":"The central discovery is the one-loop factorization formula of Eq. (5), which in the peak region takes the explicit multiplicative form $S_{\\rm QCD}(x,\\mu) = \\left[1+\\frac{\\alpha_s C_F}{2\\pi}\\left(\\frac12\\ln^2\\frac{\\mu^2}{m_b^2}-\\frac32\\ln\\frac{\\mu^2}{m_b^2}+\\frac{\\pi^2}{12}-2\\right)\\right] S_{\\rm HQET}(\\omega,\\mu)$, with $\\omega v_+ = x m_B v_+ - m_b v_+$, together with the tail-region matching function $Z_{\\rm tail}^{(1)}(x,\\mu) = \\frac{1}{m_b v_+}\\frac{1+x^2}{1-x}\\left[-1+\\ln\\frac{\\mu^2}{(1-x)^2 m_b^2}\\right]$ from Eq. (19). The apparent convolution in the factorization reduces to multiplication because the plus distributions and the infrared regulator $v\\cdot k$ cancel between the QCD and HQET one-loop amplitudes, leaving a matching coefficient that depends only on $\\mu$ and $m_b$. This establishes, at leading power in $\\Lambda_{\\rm QCD}/m_b$ and one-loop order in $\\alpha_s$, that the two shape functions have identical infrared behavior and differ only by calculable short-distance physics.","pith_inferences":["Not derived in the paper: if the one-loop multiplicative coefficient exponentiates under renormalization-group evolution, the product form would allow an all-orders resummation of $\\ln(m_b/\\Lambda_{\\rm QCD})$; checking this requires computing the anomalous dimension of the matching coefficient.","The collapse of the convolution to multiplication at one loop may be an accident of this order; at two loops the matching could develop a genuinely nonlocal dependence on $\\omega$, which would change the simple form of Eq. (23).","A direct lattice determination of the quasishape function, combined with the first-step factorization, would test the predicted peak-region relation without invoking any model of the HQET shape function; this is an extension the paper leaves to future work."],"forward_implications":["The QCD shape function obtained from lattice simulations can be converted into the HQET shape function used in inclusive $B$-decay analyses, and the conversion is exact at leading power and one-loop order.","Separating the $m_b$ scale from $\\Lambda_{\\rm QCD}$ in this way makes the large logarithms $\\ln(m_b/\\Lambda_{\\rm QCD})$ resummable through the matching coefficient.","Using a phenomenological model for the HQET shape function, the paper constructs the corresponding QCD shape function and identifies the intermediate transition region where higher-order $(1-x)$ corrections are required.","Matching at a different heavy-quark mass would connect $B$-meson and $D$-meson shape functions, with the caveat that $\\Lambda_{\\rm QCD}/m_c$ power corrections become sizable in the charm case."],"supporting_citations":[{"why":"Defines the HQET B-meson shape function that the factorization connects to.","marker":"[19]"},{"why":"Supplies the asymptotic tail behavior, the phenomenological model used in Sec. IV, and the plus-distribution machinery adapted here.","marker":"[20]"},{"why":"Introduces the equivalent distribution regularization for the shape function endpoint.","marker":"[23]"},{"why":"Defines the modified plus distribution in Eq. (12) used in the one-loop calculations.","marker":"[35]"},{"why":"Proposes the two-step lattice scheme whose second step the shape-function factorization supplies.","marker":"[29]"},{"why":"Provides the analogous QCD-to-HQET factorization for B-meson light-cone distribution amplitudes that motivates the two-step lattice scheme.","marker":"[36]"},{"why":"Derives the corresponding LCDA matching that this shape-function factorization parallels.","marker":"[37]"}],"fun_headline_variants":["One-loop bridge unifies B-meson shape functions","QCD-HQET shape functions connected by exact formula","Factorization formula links B-meson shape functions","B-meson shape functions related at one loop","New formula enables lattice QCD for inclusive B decays"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The matching assumes a clean leading-power separation between the peak region and the tail region, with the two shape functions sharing identical infrared behavior; the paper does not control the transition region between them, so if that region is broad the stated one-loop formulas would be incomplete.","fun_headline_variants_meta":{"raw":{"variants":["One-loop bridge unifies B-meson shape functions","QCD-HQET shape functions connected by exact formula","Factorization formula links B-meson shape functions","B-meson shape functions related at one loop","New formula enables lattice QCD for inclusive B decays"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000208,"raw_usage":{"total_tokens":1443,"prompt_tokens":1024,"completion_tokens":419,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":640,"completion_tokens_details":{"reasoning_tokens":345}},"tokens_in":640,"tokens_out":419,"duration_ms":4737,"temperature":1.0,"reasoning_tokens":345,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T10:27:07.923932+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A two-loop computation of the matching coefficient would settle whether the multiplicative form of Eq. (23) survives beyond one loop; alternatively, a lattice extraction of the quasishape function matched to $S_{\\rm QCD}$ could be compared with Eq. (23) across the peak-to-tail transition, where a growing discrepancy would show the region separation is not clean.","supporting_citations":[{"cited_title":"Relation between Light Cone Distribution Amplitudes and Shape Function in B mesons","cited_arxiv_id":"0707.3027","evidence_quote":"Proposes the two-step lattice scheme whose second step the shape-function factorization supplies."}],"review_version":1}