REVIEW
Direct renormalization of flavor-changing currents in vNRQCD for $B_c^{(*)}$ decays
High Energy Physics - Phenomenology · Submitted 22 Aug 2026 · 33 references
Pith's verdict: A careful vNRQCD derivation of the two-loop B_c/B_c* current anomalous dimensions that matches known results, but the direct-matching claim rests on an asserted rather than proven completeness of the potential-mode diagrams.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
At the matching scale $\nu=1$, the two-loop anomalous dimension of the leading flavor-changing currents in vNRQCD is one spin-unified formula, $\gamma_S^{(2)}(x)$, valid for any charm-to-bottom mass ratio and agreeing with full-QCD…
The reading
What carries the argument
The central machinery is the potential-mode evaluation of the two-loop current vertex in vNRQCD: after closing the energy contour, each diagram becomes a spatial integral in $d=3-2\epsilon$ dimensions with distinct unequal-mass denominators $D(q)=E-q^2/(2m_b)-q^2/(2m_c)$, and the logarithmic ultraviolet poles are isolated by the master integrals $I_0$ and $I_1$. The counterterm $\delta\tilde Z_S$ is minus the pole part, and the anomalous dimension follows by differentiating with respect to $\ln\mu_S$ at fixed bare coefficients, with the explicit powers of $\mu_S$ in the two-potential, kinetic, and $1/|\boldsymbol{k}|$ insertions producing the relative factors that combine into Eq. (49). Spin enters only through $S(S+1)$ via the spin-spin potential, so one formula serves both the pseudoscalar and vector channels.
Extended reading notes
Core claim
The paper claims that at the matching scale, where the soft and ultrasoft scales coincide ($\nu=1$, $\mu_\Lambda=\mu_S=\mu_U=m_b$), the two-loop current anomalous dimension is $\gamma_S^{(2)}(x)=-\pi^2\,[\,C_F C_A/2+(x^2+6x+1-\frac{8x}{3}S(S+1))/(2(1+x)^2)\,C_F^2\,]$, with $S=0$ for $B_c$ and $S=1$ for $B_c^*$. The counterterms are computed directly as the ultraviolet poles of effective-theory current matrix elements from potential-loop topologies, including Coulomb-potential insertions, spin and unequal-mass potentials, relativistic kinetic insertions, and a single $1/|\boldsymbol{k}|$ potential whose coefficient is imported from QCD matching. The claimed consequence is that these counterterms exactly cancel the residual infrared poles of the full-QCD short-distance coefficients for arbitrary $x$. In the equal-mass limit $x\to 1$ the formula reduces to known quarkonium anomalous dimensions, and it is symmetric under $x\leftrightarrow 1/x$, as required by exchanging charm and bottom.
Load-bearing premise
The load-bearing premise is that at the matching scale $\nu=1$ the only ultraviolet poles that renormalize the leading currents come from the potential-mode diagrams, with no soft- or ultrasoft-mode loop contributing a $1/\epsilon$ pole at that point—a completeness statement the paper asserts rather than proves.
Editorial extensions
If this is right
- The pseudoscalar and vector decay constants $f_{B_c}$ and $f_{B_c^*}$ can both be obtained from the same two-loop matching-scale renormalization constant, with spin dependence entering only through $S(S+1)$.
- Combined with the running of the long-distance matrix element given in Eq. (22), the product of short-distance coefficient and matrix element is scale independent at two-loop order, so the reported ratio $R_H$ isolates the current renormalization effect.
- The formula agrees with full-QCD matching results for both channels, so vNRQCD's mode separation is sufficient to control current renormalization at the matching point.
- Ultrasoft-induced mixing operators have Wilson coefficients that vanish at $\nu=1$, so they do not alter the fixed-order anomalous dimensions even though they are essential for velocity running below $\nu=1$.
- In the equal-mass limit the formula reduces to the known quarkonium expressions, smoothly connecting the $B_c$ sector to charmonium and bottomonium effective-theory results.
Reading between the lines
- If Eq. (49) is correct, the same direct-counterterm strategy could be extended to $O(\alpha_s^2 v^2)$ relativistic corrections, where the $S(S+1)$ structure would keep the pseudoscalar and vector channels unified.
- The explicit mass-ratio dependence and $x\leftrightarrow 1/x$ symmetry suggest the formula could serve as a perturbative anchor for interpolation between heavy-light and quarkonium limits, a use the paper does not itself explore.
- Once the full unequal-mass velocity running of the potentials is included, combining it with this matching-scale result would yield a complete VRG prediction for $B_c$ and $B_c^*$ decay-constant matrix elements that lattice QCD could test after current-normalization matching.
Original abstract
We present a direct vNRQCD calculation of the two-loop anomalous dimensions of the leading flavor-changing nonrelativistic currents that govern the leptonic decays of the $B_c$ and $B_c^*$ mesons. The calculation is performed at the matching scale $\nu=1$, where the anomalous dimensions are the ultraviolet counterterms required to cancel the residual infrared poles of the full-QCD short-distance coefficients. In contrast to the conventional extraction from those infrared poles, we determine the current renormalization constants directly in the effective theory. The pseudoscalar and vector channels are obtained from one spin-unified formula, with the spin dependence entering through $S(S+1)$. We also clarify the role of ultrasoft-induced mixing operators in the unequal-mass theory: they are essential for the full velocity running below $\nu=1$, but their Wilson coefficients vanish at the matching point and do not alter the fixed-order current anomalous dimensions reported here. The results agree with the known matching calculations and reduce smoothly to the equal-mass quarkonium limits.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
We present a direct vNRQCD calculation of the two-loop anomalous dimensions of the leading flavor-changing nonrelativistic currents that govern the leptonic decays of the $B_c$ and $B_c^*$ mesons. The calculation is performed at the matching scale $\nu=1$, where the anomalous dimensions are the ultraviolet counterterms required to cancel the residual infrared poles of the full-QCD short-distance coefficients. In contrast to the conventional extraction from those infrared poles, we determine the current renormalization constants directly in the effective theory. The pseudoscalar and vector channels are obtained from one spin-unified formula, with the spin dependence entering through $S(S+1)$. We also clarify the role of ultrasoft-induced mixing operators in the unequal-mass theory: they are essential for the full velocity running below $\nu=1$, but their Wilson coefficients vanish at the matching point and do not alter the fixed-order current anomalous dimensions reported here. The results agree with the known matching calculations and reduce smoothly to the equal-mass quarkonium limits.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | We present a direct vNRQCD calculation of the two-loop anomalous dimensions of the leading flavor-changing nonrelativistic currents that govern the leptonic decays of the $B_c$ and $B_c^*$ mesons. | contribution 1.00 | 0.97 | 0.48 | asserted 1.00 |
| s1 | The calculation is performed at the matching scale $\nu=1$, where the anomalous dimensions are the ultraviolet counterterms required to cancel the residual infrared poles of the full-QCD short-distance coefficients. | method 0.97 | 0.84 | 0.08 | asserted 1.00 |
| s2 | In contrast to the conventional extraction from those infrared poles, we determine the current renormalization constants directly in the effective theory. | method 0.70 | 0.94 | 0.17 | asserted 1.00 |
| s3 | The pseudoscalar and vector channels are obtained from one spin-unified formula, with the spin dependence entering through $S(S+1)$. | method 0.81 | 0.89 | 0.07 | asserted 1.00 |
| s4 | We also clarify the role of ultrasoft-induced mixing operators in the unequal-mass theory: they are essential for the full velocity running below $\nu=1$, but their Wilson coefficients vanish at the matching point and do not alter the fixed-order current anomalous dimensions reported here. | interpretation 0.67 | 0.73 | 0.12 | asserted 1.00 |
| s5 | The results agree with the known matching calculations and reduce smoothly to the equal-mass quarkonium limits. | result 0.81 | 0.82 | 0.21 | asserted 1.00 |
What was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
Figures reproduced from arXiv:2608.21771 by the authors.
References
How the literature uses it
What this paper does with the arXiv papers it cites
What it cites · 46 sentences without a confident label
Uses the method or tool of 1 paper
-
Near threshold, the homogeneous momentum regions are [ref] [equation] The familiar notation $p∼mv$ and $E∼mv^2$ is only a shorthand for these scalings.
Asymptotic expansion of Feynman integrals near threshold · Threshold expansion · 0.85
Uses data or values from 5 papers
-
Equations [ref] and [ref] therefore give [equation] Only after this differentiation do we insert the unequal-mass matching coefficient [equation] taken from Ref. [ref] .
Potential NRQCD for unequal masses and the $B_c$ spectrum at NNNLO · Reduction of the potential loops · 0.98
-
The potential sector has the form [equation] For a $cb̅$ system with unequal masses, the potential in the color basis used below can be written as [ref] [equation] where $𝐤=𝐩'-𝐩$, $m_r=m_bm_c/(m_b+m_c)$, and $𝐒=(σ_1+σ_2)/2$.
The QCD heavy-quark potential to order v^2: one loop matching conditions · Label fields and effective Lagrangian · 0.82
-
The potential sector has the form [equation] For a $cb̅$ system with unequal masses, the potential in the color basis used below can be written as [ref] [equation] where $𝐤=𝐩'-𝐩$, $m_r=m_bm_c/(m_b+m_c)$, and $𝐒=(σ_1+σ_2)/2$.
Ultrasoft Renormalization in Non-Relativistic QCD · Label fields and effective Lagrangian · 0.82
Compares with 6 papers
-
Equation [ref] reproduces the two-loop anomalous dimension of the pseudoscalar current inferred from the $B_c$ leptonic decay amplitude [ref] , while Eq. [ref] reproduces the corresponding vector-current result [ref] .
Three-loop QCD corrections to the decay constant of $B_c$ · Reduction of the potential loops · 0.98
-
Equation [ref] reproduces the two-loop anomalous dimension of the pseudoscalar current inferred from the $B_c$ leptonic decay amplitude [ref] , while Eq. [ref] reproduces the corresponding vector-current result [ref] .
Two-loop QCD Corrections to $B_c$ Meson Leptonic Decays · Reduction of the potential loops · 0.98
-
Equation [ref] reproduces the two-loop anomalous dimension of the pseudoscalar current inferred from the $B_c$ leptonic decay amplitude [ref] , while Eq. [ref] reproduces the corresponding vector-current result [ref] .
Next-to-next-to-leading order matching of beauty-charmed meson $B_{c}$ and $B^*_{c}$ decay constants · Reduction of the potential loops · 0.98
Cites as support 1 paper
-
It also provides a natural intermediate scale between the bottom and charm masses, consistent with the HPQCD observation that the physical $B_c$ system exhibits features of both heavy-heavy and heavy-light dynamics [ref] .
Heavy meson masses and decay constants from relativistic heavy quarks in full lattice QCD · Illustrative running of the long-distance matrix elements · 0.98
Cites as background 8 papers
-
For the vector channel, the NNLO matching coefficient was obtained in Ref. [ref] , and the $B_c^*$ decay constant was subsequently evaluated through $O(α_s^3)$ [ref] .
2210.02979 · Introduction · 0.82
-
For the vector channel, the NNLO matching coefficient was obtained in Ref. [ref] , and the $B_c^*$ decay constant was subsequently evaluated through $O(α_s^3)$ [ref] .
Next-to-next-to-leading order matching of beauty-charmed meson $B_{c}$ and $B^*_{c}$ decay constants · Introduction · 0.82
-
The one-loop QCD correction to $f_B_c$ was obtained by Braaten and Fleming [ref] .
QCD Radiative Corrections to the Leptonic Decay Rate of the B_c Meson · Introduction · 0.87
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.