{"id":"d4fe248a-ac7b-4c51-bf48-306f2b5e9254","arxiv_id":"2505.08611","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Using the renormalon model, the paper estimates that power corrections to the QCD-to-HQET light-cone distribution amplitude matching are around 22% for the D meson and 7% for the B meson in the peak region.","lead":"This paper estimates missing power corrections in the standard method for extracting heavy meson light-cone distribution amplitudes from lattice QCD. Using a renormalon model, it finds corrections around 22% for D mesons and 7% for B mesons, which matter for precision predictions of heavy quark decays.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 22% and 7% estimates rest on dropping the real theta(u-s) term of the vertex bubble diagram (from Eq. (29) vs. Eq. (A2)); no quantitative proof that it is power-suppressed in the peak region is given.","rationale":"The paper has real strengths: the two calculational routes (direct bubble-chain sum and Borel-transformed propagator) cross-check each other, and the pole locations follow transparently from the Gamma functions in Eq. (40). The concern is not about the renormalon framework or about the phenomenological LCDA models, both of which are legitimate modeling choices. The weak point is narrower: every route starts from Eq. (29), where the real-emission part of the same diagram is discarded. The full result in Appendix A includes that part, and the claim that it is absent in the peak region is asserted rather than derived. The reader's weakest_assumption identifies exactly this step, and my reading agrees. Because both the pole set and the residues feed directly into the 22% and 7% numbers, this is the most load-bearing assumption in the paper. I would keep the reader's CONDITIONAL verdict rather than move it: the calculation is a plausible renormalon-model estimate, but the headline numbers should be contingent on a proof of the virtual-part reduction or on an explicit computation retaining the real term.","tokens_in":17620,"tokens_out":12302,"duration_ms":128359,"concrete_test":"Recompute the matching kernel at fixed bubble order n from the full QCD-A result of Eq. (A2), retaining the theta(u-s) term through the Fourier transform and through the convolution (7) with an exponential HQET LCDA phi_+(omega) = (omega/omega_0^2) exp(-omega/omega_0), in the peak region u ~ Lambda_QCD/m_Q; then take the Borel transform of the resulting n-dependent coefficients and compare the residues at w = 1/2 and w = 1 with Eq. (40). If the retained residues change at leading order in lambda, the quoted 22% and 7% corrections are not supported; if the real term is genuinely suppressed, the present numbers stand.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central numerical result is not yet secure because Eq. (29) drops the real-emission term of the QCD bubble-chain vertex diagram, while the full bare result in Eq. (A2) contains a term proportional to theta(u-s) / (u-s)^(2(n+1)epsilon+1). The text says that only the virtual part of Fig. 1(a) contributes in the peak region and that this can be demonstrated by expansion by regions or by explicit calculation, but no demonstration is shown. In the convolution (7), s = omega/m_H, and the HQET LCDA phi_+(omega) peaks at omega of order Lambda_QCD, so for u ~ lambda the difference u-s is also of order lambda; the discarded term is not suppressed by a power of m_Q relative to the delta(u-s) term. In addition, in epsilon expansion the product theta(u-s)/(u-s)^(1+2(n+1)epsilon) contains a 1/epsilon delta(u-s) piece plus a plus-distribution, so it cannot be dismissed as a power-suppressed effect by power counting alone. Since the Borel poles and residues of Eqs. (40)-(41), and therefore the 22% (D) and 7% (B) corrections, are computed from the virtual-only expression, an unquantified contribution from the real term could change the headline numbers.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper estimates the power corrections to the factorization formula that matches heavy-meson QCD light-cone distribution amplitudes onto bHQET/HQET LCDAs, using the renormalon model. The authors compute the bubble-chain vertex diagram of Fig. 1(a) in coordinate space, renormalize the fermion-loop insertions, and derive a Borel transform of the matching kernel whose singularities they locate at w = n + 1/2 (n ∈ N) and w = 1. Identifying the renormalon ambiguity with the leading Λ_QCD/m_Q corrections, they then use phenomenological models for the D- and B-meson LCDAs and report power corrections of about 22% for the D meson and 7% for the B meson in the peak region.","tokens_in":17902,"tokens_out":7043,"duration_ms":74229,"significance":"If the claims are correct, the paper provides a concrete, quantitative warning that leading-power QCD-HQET matching is insufficient for charm physics and only marginally safe for bottom physics. The work has several strengths: the analytic bubble-chain calculation is presented in detail, the bare results for all four diagrams are collected in Appendix A, two independent calculational methods (direct bubble-chain evaluation and Borel-transformed propagator) are used as cross-checks, and the final numbers are concrete and falsifiable by future lattice determinations of HQET LCDAs. The main weakness is that the numerical results rest on an unproven truncation to the virtual part of the vertex diagram, so the headline percentages are not yet secure.","major_comments":[{"comment":"The central numerical results rely on dropping the theta(u-s) term in the full bare QCD vertex diagram. The text states that in the peak region only the virtual part of Fig. 1(a) contributes and that this follows from expansion by regions or explicit calculation, but no demonstration is given. In the convolution Eq. (7), s = ω/m_H and the HQET LCDA φ_+(ω) peaks at ω of order Λ_QCD, while u is of order λ = Λ_QCD/m_Q; hence u-s is generically of the same order as u and s, not suppressed by a power of m_Q. Moreover, the discarded term in Eq. (A2), proportional to θ(u-s)/(u-s)^{2(n+1)ε+1}, contains in its ε expansion a 1/ε δ(u-s) piece plus a plus-distribution, so it cannot be dismissed by simple power counting. Since the Borel pole positions and residues in Eqs. (40)-(41), and therefore the 22% and 7% estimates, are computed from the virtual-only expression, please either supply the promised expansion-by-regions proof or compute the real-emission term explicitly and show that its contribution to the final ratios is numerically negligible.","section":"Sec. II.B and Eq. (29) vs. Eq. (A2)"},{"comment":"The numerical pipeline from the Borel residues to the reported percentages is not shown. No explicit formula is given for δφ_+/φ_+ in terms of Res B[J](w), the set of poles included, the definition of Λ_QCD, the values of α_s, or the sum over singularities. The statement that the ratio is independent of the momentum fraction u is asserted but not derived; it presumably follows from a formula such as δφ_+/φ_+ = δJ/J in the peak region, but the reader cannot verify this from Eqs. (39)-(41) alone. Please present the complete expression used for Figs. 4 and 6 and list the numerical inputs, so that the 22% and 7% numbers can be reproduced.","section":"Sec. IV, Figs. 4 and 6"}],"minor_comments":[{"comment":"The index in the term f^{[-1]}_{i+1} appears to be a typo; it should presumably be f^{[-1]}_{n+1}. If not, the index i is undefined.","section":"Eq. (38)"},{"comment":"References [35] and [71] are identical (Han et al., Lattice Parton Collaboration, Phys. Rev. D 111, 034503 (2025)). Please remove the duplicate. In addition, references [33] and [73]-[78] do not appear to be cited in the text; please check the citation list.","section":"References"},{"comment":"The equation for f_H appears to have a typographical formatting issue: the right-hand side should show the tilde on the first decay constant, i.e. f_H = \\tilde f_H(μ)(1 - α_s C_F/(4π)(3/2 ln μ^2/m_Q^2 + 2) + ...). Please fix the notation so the matching relation is unambiguous.","section":"Eq. (6)"},{"comment":"The replacement μ^{2w} e^{-w/(β_0 a_s)} = Λ_{QCD}^{2w} implicitly fixes a scheme and running-order convention for Λ_QCD. Please state the convention (e.g., one-loop MS-bar with a specified number of active flavors and reference value) so the numerical estimates are reproducible.","section":"Eq. (41)"}],"recommendation":"major_revision","confidential_remarks":"The paper is competently written and the technical machinery is mostly transparent, but the missing quantitative treatment of the theta(u-s) term in Eq. (A2) is a genuine load-bearing gap. The fix is within scope: the authors can either prove the virtual-only dominance by an explicit expansion-by-regions calculation or include the real term and show its numerical effect. I would be willing to reconsider after that point is addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nRead the paper. What is actually new: this is the first bubble-chain renormalon calculation for the recent QCD-to-HQET LCDA matching kernel, and it produces a Borel pole set (w = n+1/2 and w = 1) that differs from the usual HQET renormalon pattern. The analytic part is genuinely worked out: coordinate-space Schwinger calculation, two independent methods (two-step renormalization and Borel-transformed propagator) that cross-check, and full bare results in Appendix A. That is real work, and the authors mostly avoid overclaiming—they call it an estimate and note model dependence.\n\nThe soft spot is serious, and it is exactly the one flagged in the stress test. Equation (29) drops the theta(u-s)/(u-s)^{2(n+1)epsilon+1} term from the bare QCD diagram A, with only the assertion that in the peak region only the virtual part contributes. No expansion-by-regions demonstration is given. The stress-test objection lands: in the convolution (7), both u and s are of order lambda in the peak region, so u-s is not a large scale, and the epsilon expansion of theta(x)/x^{1+a epsilon} contains delta(x)/epsilon. That term is not obviously power-suppressed relative to the delta(u-s) virtual term. Since the Borel residues in Eqs. (40)-(41), and therefore the 22%/7%, come from the virtual-only expression, this is a load-bearing gap rather than a cosmetic one. The claim may still be true—the real term could cancel against QCD diagram B or the HQET diagrams, or be genuinely subleading by regions—but the paper needs to show it.\n\nA second, smaller weakness is that the numerical pipeline is not fully auditable: no alpha_s or Lambda_QCD inputs, no explicit residue sum, no error bars on 22% and 7%. The scale-variation bands in Figs. 4 and 6 are suspiciously flat, which may be a ratio artifact but needs a sentence. And since Ref. [42] shares three authors and overlapping machinery, the comparison to that work should be explicit; the reader should not have to guess what is new.\n\nIf the virtual-only step is correct, the pole positions follow and the estimates are plausible. On the evidence in the text, I would treat the numerical results as conditional. But the question matters: a 22% D-meson correction would directly affect the lattice extraction of HQET LCDAs and exclusive B/D phenomenology. This deserves a serious referee. I would send it to review and ask for the missing demonstration of the virtual-only reduction and a transparent numerical input table.","headline":"A plausible first renormalon estimate of power corrections in the QCD-to-HQET LCDA matching, but the headline 22%/7% numbers rest on an unproven 'virtual part only' step that the authors need to demonstrate.","tokens_in":18482,"tokens_out":3642,"would_cite":false,"duration_ms":37970,"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":"Renormalon-model calculation puts the power corrections to the QCD-to-HQET light-cone distribution amplitude matching at roughly 22% for the D meson and 7% for the B meson.","keywords":["heavy meson LCDA","HQET","power corrections","renormalon model","Borel transform","bubble chain diagrams","D meson","B meson"],"falsifier":"Evaluate the Borel transform of the complete bare QCD bubble-chain result including the discarded real $\\theta(u-s)$ term, and compare its residues at $w=1/2$ and $w=1$ with the virtual-only residues used here; alternatively, a lattice extraction of the $D$-meson HQET LCDA that directly measures the size of the $1/m_c$ correction would settle whether 22% is the right order of magnitude.","tokens_in":17383,"feed_emoji":"⚛️","tokens_out":10363,"duration_ms":84977,"temperature":0.7,"pith_summary":"To extract the light-cone distribution amplitude (LCDA) of a heavy meson, one matches the QCD LCDA onto the simpler HQET LCDA at leading power. This paper estimates the power-suppressed corrections to that matching using the renormalon model, in which the divergent tail of the perturbative matching kernel parameterizes higher-twist contributions. The authors compute the vertex bubble-chain diagrams and find Borel-plane poles at $w=n+\\frac12$ for all $n$ and at $w=1$, with the dominant ambiguity coming from the virtual part of the vertex diagram. Using phenomenological models for both LCDAs, the correction is about 22% for the $D$ meson and 7% for the $\\overline{B}$ meson in the peak region. If those numbers are right, leading-power matching is not sufficient for high-precision charm physics and needs a few-percent correction even for bottom.","feed_headline":"Power corrections hit D-meson LCDAs at 22%, B at 7%","feed_subtitle":"Renormalon analysis: charm is not heavy enough for leading-power accuracy, and bottom needs ~7% correction.","key_machinery":"The central object is the Borel transform of the matching kernel $J_{\\rm peak}$ that connects the QCD and HQET LCDAs. In the renormalon model the perturbative series is made divergent by chains of fermion bubbles inserted into the gluon propagator; after renormalization the factorial growth appears as poles in the Borel plane, and the residue of each pole gives the renormalon ambiguity that stands in for a higher-twist power correction. The calculation either renormalizes the bubble chain step by step or uses the Borel-transformed bubble-chain propagator, reducing the vertex diagram to a residue formula involving $\\Gamma(-2w)\\Gamma(w)/\\Gamma(2-w)$ with $e^{5w/3}$. Evaluating the residues at $w=n+\\frac12$ and $w=1$ produces the 22% and 7% estimates.","core_discovery":"The paper claims that in the renormalon model the dominant power corrections to the QCD-to-HQET LCDA matching kernel come from the virtual part of the vertex bubble-chain diagram, whose Borel transform has poles at $w=n+\\frac12$ for every non-negative integer $n$ and at $w=1$. These poles make the Borel integral ambiguous at precisely the orders expected for higher-twist operators, $(\\Lambda_{\\rm QCD}/m_Q)^{2w}$. With a polynomial parametrization of the QCD LCDA and an exponential model for the HQET LCDA, the renormalon ambiguity divided by the one-loop matching result is approximately 22% for the $D$ meson and 7% for the $\\overline{B}$ meson at $\\mu\\sim m_Q$; because the matching kernel is independent of the light-cone momentum fraction, the correction is a constant shift of the peak-region amplitude rather than a distortion of its shape.","pith_inferences":["A natural next step is to compute the full bubble-chain result without dropping the real $\\theta(u-s)$ term; if its Borel residue is comparable, the saturation assumption breaks and the numerical estimates are not robust.","Extrapolating the mass dependence shown in the paper's Fig. 7, the renormalon estimate grows as the meson mass drops below the charm scale, so the model is least trustworthy in exactly the region where the correction is largest; explicit higher-twist matrix elements would be needed there.","The same machinery could be applied to other factorization theorems with a heavy-quark scale, such as $B$-meson form factors or $W$-decay-to-$B$ processes, to produce comparable power-correction estimates before lattice input is available."],"forward_implications":["Charm-quark extractions of the HQET LCDA from lattice QCD will need an explicit treatment of power corrections or a model for higher-twist operators, since a 22% normalization shift is larger than typical target precision.","For the bottom quark, leading-power matching carries an intrinsic uncertainty of order 7%, so precision comparisons with $B$-factory data should include this correction or estimate it independently.","Because the matching kernel is independent of the light-cone momentum fraction $u$, the renormalon power correction changes the overall size of the peak-region LCDA rather than its shape.","Heavy-quark spin symmetry extends the size of the estimated correction to vector mesons, so the same percent-level power corrections should be expected there."],"supporting_citations":[{"why":"Supplies the QCD-to-HQET factorization formula and the polynomial-moment input for the D and B meson LCDAs used in the numerics.","marker":"[31]"},{"why":"Supplies the two renormalization procedures (two-step renormalization and Borel-transformed propagator) used to compute the bubble chains.","marker":"[38]"},{"why":"Establishes the renormalon and ultraviolet-dominance logic that connects the factorial growth of the series to higher-twist power corrections.","marker":"[40]"},{"why":"Justifies naive nonabelianization, the replacement of fermion-loop insertions by the one-loop beta-function coefficient in the bubble-chain calculation.","marker":"[69]"},{"why":"Provides the phenomenological exponential model for the D-meson HQET LCDA used in the numerical estimate.","marker":"[71]"}],"fun_headline_variants":["Renormalon poles: 22% power correction to D-meson LCDA","Vertex bubble chain drives 22% (D) and 7% (B) corrections","Power corrections from renormalons: D 22%, B 7%","Renormalon ambiguity: 22% for D meson, 7% for B meson"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The estimate rests on the assumption that in the peak region the power correction is saturated by the virtual part of the vertex bubble-chain diagram; the real part proportional to $\\theta(u-s)$ is discarded from the bare QCD result (equation A2), and if that discarded piece contributes at the same order, the pole residues and the 22%/7% numbers would change.","fun_headline_variants_meta":{"raw":{"variants":["Renormalon poles: 22% power correction to D-meson LCDA","Vertex bubble chain drives 22% (D) and 7% (B) corrections","Power corrections from renormalons: D 22%, B 7%","Renormalon ambiguity: 22% for D meson, 7% for B meson"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000227,"raw_usage":{"total_tokens":1469,"prompt_tokens":939,"completion_tokens":530,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":555,"completion_tokens_details":{"reasoning_tokens":438}},"tokens_in":555,"tokens_out":530,"duration_ms":4376,"temperature":1.0,"reasoning_tokens":438,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:50:50.079008+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the Borel transform of the complete bare QCD bubble-chain result including the discarded real $\\theta(u-s)$ term, and compare its residues at $w=1/2$ and $w=1$ with the virtual-only residues used here; alternatively, a lattice extraction of the $D$-meson HQET LCDA that directly measures the size of the $1/m_c$ correction would settle whether 22% is the right order of magnitude.","supporting_citations":[{"cited_title":"Zichichi, doi:10.1007/978-1-4613-4208-3 13","cited_arxiv_id":null,"evidence_quote":"Provides the phenomenological exponential model for the D-meson HQET LCDA used in the numerical estimate."}],"review_version":1}