REVIEW 3 major objections 5 minor 31 references
Update of HPQCD $B_c\to J/\psi$ Form Factors
T0 review · 3 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read This lattice QCD update determines the Bc to J/psi semileptonic form factors across the full kinematical range with roughly twice the previous precision, and adds tensor form factors for the first time.
desk verdict Solid proceedings update with genuinely new tensor form factors and an ultra-fine ensemble, but the factor-of-two precision gain rests on a heuristic mass interpolation that carries no assigned systematic error. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is the full dispersive $z$-expansion parameterisation of the form factors, $F_Y(z)=[P_Y(z)\phi_Y(z)]^{-1}\sum_n a^Y_n p_n(z)$, where $P_Y$ are Blaschke factors built from the masses of single-particle $\bar h c$ states below the pair-production threshold, $\phi_Y$ are outer functions depending on continuum susceptibilities, and $p_n$ are orthonormal polynomials on the unit circle. Around this sits the physical-continuum extrapolation model: heavy-quark mass dependence enters through heuristic linear interpolation formulas for masses not computed on the lattice, $M_{\mathrm{pole},i}=M_{\mathrm{pole},i}^{\mathrm{phys}}+M_{H_c}^{\mathrm{latt}}-M_{B_c}^{\mathrm{phys}}$ and $M_H=M_B^{\mathrm{phys}}+M_{H_c}^{\mathrm{latt}}-M_{B_c}^{\mathrm{phys}}$, with the $z$-expansion coefficients truncated at $O(z^4)$ under uniform priors, together with quark-mass mistuning and discretisation terms. The same machinery converts the lattice matrix elements into the QCD-basis form factors $V,A_0,A_1,A_{12},T_1,T_2,T_{23}$ and carries them to the physical continuum.
What would settle it
Repeat the physical-continuum extrapolation with an alternative heavy-quark mass dependence, for example adding a curvature term to the interpolation formulas or replacing $M_H$ by $1/M_H$ scaling, and check whether the physical-point form factors move by more than the quoted uncertainties; a particularly direct test is whether the new $a\approx 0.03\,\mathrm{fm}$, $am_h\approx 0.625$ ensemble point, which sits near the physical bottom mass, is reproduced by the fitted $m_h$ trend within errors.
Extended reading notes
Core claim
The central claim is that lattice QCD, using the Highly Improved Staggered Quark action for all valence quarks and nonperturbatively renormalised current operators, now determines the $B_c\to J/\psi$ vector and axial-vector form factors and, newly, the tensor form factors across the full $q^2$ range with roughly twice the precision of the previous calculation and consistent with it. The improvement comes from adding a physical-light-quark ensemble and a very fine ensemble that reaches the physical bottom quark mass, so the chiral-continuum extrapolation is anchored directly at the physical heavy-quark point. The extrapolation uses the full dispersive parameterisation in the QCD basis, with Blaschke factors built from sub-threshold pole masses and outer functions determined from continuum susceptibilities computed on the lattice. The tensor form factors are presented as essential for constraining possible new physics in this decay mode.
Load-bearing premise
The load-bearing assumption is that the heavy-quark mass dependence of the meson masses not computed on the lattice is captured by the simple linear interpolation formulas $M_{\mathrm{pole},i}=M_{\mathrm{pole},i}^{\mathrm{phys}}+M_{H_c}^{\mathrm{latt}}-M_{B_c}^{\mathrm{phys}}$ and $M_H=M_B^{\mathrm{phys}}+M_{H_c}^{\mathrm{latt}}-M_{B_c}^{\mathrm{phys}}$; if that linear dependence is wrong, the pole and outer-function factors in the $z$-expansion are wrong and the extrapolated form factors are biased.
Editorial extensions
If this is right
- The Standard Model vector and axial-vector form factors are consistent with the earlier calculation and roughly twice as precise, so $B_c\to J/\psi\ell\nu$ rate predictions inherit a smaller uncertainty.
- The new tensor form factors complete the lattice description of the decay's operator basis, giving a direct handle on new-physics contributions in $B_c\to J/\psi\tau^+\nu_\tau$ versus $B_c\to J/\psi\mu^+\nu_\mu$.
- With data on two ensembles spanning the full kinematical range, the analysis can test how much the choice of parameterisation and truncation order changes integrated observables, addressing a known source of underestimated uncertainty in $B\to D^*\ell\nu$ analyses.
- A sharper Standard Model prediction for $R(J/\psi)$ follows if the quoted precision is confirmed, making the existing experimental measurements a more discriminating lepton-flavour-universality test.
Reading between the lines
- The paper's heuristic mass-interpolation formulas could be checked by adding a third heavy-quark mass value on the fine ensemble and assigning a systematic uncertainty; until then, the factor-of-two gain in precision may understate the true model dependence.
- The same full dispersive parameterisation could be applied to other heavy-to-heavy exclusive decays, such as $B_s\to D_s^*\ell\nu$, where parameterisation-related tensions with $B\to D^*$ have been reported.
- If the observed disagreement between the lattice axial-tensor susceptibility and $O(\alpha_s)$ perturbation theory persists, the new tensor form factors offer a clean observable in which to test whether that discrepancy is perturbative or physical.
- A testable extension is to compare the new $z$-expansion shapes for $T_1,T_2,T_{23}$ against future differential decay-rate measurements in $B_c\to J/\psi\ell\nu$, checking whether truncation at $O(z^4)$ is adequate.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This proceedings paper reports an update of the lattice QCD calculation of the B_c -> J/psi semileptonic form factors (vector, axial-vector, and new tensor) using the HISQ action on the MILC n_f=2+1+1 ensembles. Two new ensembles are used: one with near-physical light quarks at a≈0.06 fm and one at a≈0.03 fm where the physical b quark can be simulated. The correlation functions are fit with a spectral decomposition, the currents are renormalised nonperturbatively (with tensor renormalisation in RI-SMOM matched to MS), and the physical-continuum extrapolation uses a full dispersive parameterisation with Blaschke factors, outer functions, and lattice-computed susceptibilities. The paper claims that the new SM form factors are consistent with the previous HPQCD result and roughly a factor of two more precise.
Significance. If substantiated, this update would provide a more precise determination of the B_c -> J/psi form factors across the full kinematic range, relevant for R(J/psi) and new-physics searches. The addition of the physical-light-quark ensemble and the fine ensemble reaching the physical b mass is a clear strength, as is the use of nonperturbative renormalisation and the full dispersive parameterisation with lattice susceptibilities. The new tensor form factors are a valuable addition. However, the paper's quantitative claims are not backed by numerical tables or a systematic error budget, and the treatment of the heavy-quark mass interpolation in the extrapolation is heuristic. The significance is therefore conditional on the details promised for the full paper.
major comments (3)
- [Sec. 2.2, mass interpolation formulas] The heuristic relations M_pole,i = M_phys_pole,i + M_latt_Hc - M_phys_Bc and M_H = M_phys_B + M_latt_Hc - M_phys_Bc are introduced without derivation or an assigned systematic uncertainty. These masses enter the Blaschke factors P_Y(z,t+,t0), the outer functions phi_Y, and the conformal variable z through t+ and t0. Because the lattice heavy-quark masses differ from the physical b mass by several hundred MeV on the coarser ensembles, the omitted O(Lambda_QCD/m_h) corrections to the linear shift could bias the central values at the percent level, comparable to the claimed factor-of-two improvement in precision. The authors should either justify the heuristic with a lattice calculation of the relevant pole masses or assign and propagate a systematic uncertainty.
- [Sec. 3, Conclusions] The claim that the new SM form-factor results are 'consistent with our previous calculation, though roughly a factor of 2 more precise' is not supported by any numerical evidence in the paper. No table of the physical-continuum form factors with uncertainties is provided, and the only comparison with the previous result [3] is the visual overlay in Fig. 1. To substantiate the precision claim and to make the update usable, the paper should include a table of the final form factors (or the z-expansion coefficients) with total uncertainties, along with a quantitative comparison of the new uncertainties with those of [3] (e.g., the ratio of uncertainty bands at selected q^2 values).
- [Sec. 2.2, O(z^4) truncation and priors] The z-expansion is truncated at O(z^4) and the coefficients are assigned uniform priors |a_Y^n|<=1, but no variation of the truncation order or prior width is reported. The conclusions acknowledge that 'This will allow for a detailed study of systematic uncertainties associated with the choice of kinematical parameterisation,' which indicates the present paper does not yet assess that model dependence. Since the claimed improvement in precision is the central message, the authors should demonstrate stability of the extrapolated form factors under, for example, increasing the truncation to O(z^5) or changing the prior range, and include the spread as a systematic uncertainty.
minor comments (5)
- [Eq. (11)] The symbol N_n is used before it is defined in Eq. (13); the definition should be moved forward or the notation should be made consistent (e.g., N_Y(s)_n).
- [Sec. 2.2, after Eq. (10)] The text says the coefficients 'satisfy the bounds sum_{Y->Gamma} sum_n |a_Y^n|^2' but does not state the actual bound value; the later imposition of the weaker bound |a_Y^n|<=1 should be introduced more explicitly as a choice.
- [Fig. 1 caption] The red dashed lines are described as the '±1σ confidence interval' of the previous result; in a Bayesian context this is usually a credible interval, so the wording should be harmonised with the statistical framework used.
- [General] The paper would benefit from a short summary of the fit quality (e.g., chi^2/dof) for both the correlator fits and the physical-continuum extrapolation fit, since Fig. 1 alone does not quantify the goodness of fit.
- [Abstract and Sec. 2] The abstract states the a≈0.03 fm ensemble 'reaches the physical bottom quark mass', while Sec. 2 specifies that this ensemble has heavier-than-physical light quarks; this is consistent but should be made explicit in the abstract to avoid confusion.
Circularity Check
No significant circularity: the form factors are extracted from independent lattice three-point correlator data, and the cited prior results (renormalization factors, susceptibilities) are independent external inputs rather than fitted outputs.
full rationale
The paper's derivation chain starts from lattice two- and three-point correlation functions (Sec. 2, Eqs. (1)-(3)) and extracts matrix elements that are then converted to form factors via Eqs. (4)-(6). The physical-continuum extrapolation in Sec. 2.2 uses the dispersive z-expansion of Eq. (10), with Blaschke factors and outer functions. The outer functions depend on susceptibilities taken from Ref. [30], a separate lattice QCD calculation of current-current correlators by the same author; this is an independent input, not the B_c -> J/psi form-factor data being predicted, and it is externally falsifiable. The renormalization factors from Refs. [27,23] are likewise independent calculations. The heuristic mass formulas M_pole,i = M_phys_pole,i + M_latt_Hc - M_phys_Bc and M_H = M_phys_B + M_latt_Hc - M_phys_Bc are admittedly not derived and carry no explicit systematic uncertainty, but they are modeling assumptions used to set inputs, not quantities derived from the fit; they do not make any output equal to an input by construction. The consistency statement comparing to the previous HPQCD result [3] is a comparison, not a fitted input. No equation in the paper reduces to another by definition, and no fitted parameter is renamed as a prediction. The claimed factor-of-two improvement in precision is therefore a genuine result of new lattice data and the stated fit procedure, subject to the unquantified model dependence noted for the heavy-quark mass interpolation, which is a correctness risk rather than circularity.
Assumptions & free parameters
free parameters (3)
- aY_n z-expansion coefficients =
not reported
- mistuning coefficients AY_n, BY_n, CY_n, DY_n =
not reported
- discretization coefficients c(nu,Gamma),jkl_n =
not reported
assumptions (6)
- domain assumption The HISQ action on the second-generation MILC n_f=2+1+1 gluon field configurations is a valid discretization of QCD, and the continuum limit is described by the extrapolation ansatz in Eq. (17).
- standard math The dispersive parameterisation of Ref. [26], including Blaschke factors, outer functions, and unitarity bounds, correctly describes the continuum form factors over the full kinematic range.
- ad hoc to paper The heuristic heavy-quark mass interpolation formulas M_pole,i = M_phys_pole,i + M_latt_Hc - M_phys_Bc and M_H = M_phys_B + M_latt_Hc - M_phys_Bc correctly model the m_h dependence of masses not computed on the lattice.
- ad hoc to paper Truncating the z-expansion at O(z^4) with uniform priors |aY_n| <= 1 is sufficient for the full kinematic range.
- domain assumption The RI-SMOM renormalization factors for the vector, axial-vector, and tensor currents are correct, and the 3-loop running of the tensor factors to m_h(m_h) is valid.
- domain assumption The heavy-light susceptibilities computed in Ref. [30] are correct enough to normalize the outer functions in Eq. (10), despite the noted disagreement with O(alpha_s) perturbation theory for the tensor susceptibility.
Cite this review
Pith. "Pith review of Update of HPQCD $B_c\to J/\psi$ Form Factors." pith.science (2026). https://pith.science/paper/4QNZQVY3
@misc{pith2026250115180,
author = {Pith},
title = {Pith review of: Update of HPQCD $B_c\to J/\psi$ Form Factors},
year = {2026},
howpublished = {\url{https://pith.science/paper/4QNZQVY3}},
note = {Machine review of arXiv:2501.15180}
}
abstract
We present an update of our lattice QCD determination of the $B_c\to J/\psi$ vector and axial-vector form factors, including new results for the tensor form factors. We use the Highly Improved Staggered Quark action for all valence quarks, together with the second generation MILC $n_f=2+1+1$ HISQ gluon field configurations. This calculation includes two additional ensembles, one with physically light up and down quarks and $a\approx 0.06 \mathrm{fm}$ and one with $a\approx 0.03\mathrm{fm}$ on which we are able to reach the physical bottom quark mass. Our calculation uses nonperturbatively renormalised current operators and covers the full kinematical range of the decay. Our physical-continuum extrapolation utilises the full dispersive parameterisation for $B_c\to J/\psi$.
Figures
Reference graph
Works this paper leans on
-
[3]
J. Harrison, C.T.H. Davies and A. Lytle,𝐵𝑐→𝐽/𝜓 form factors for the full𝑞2 range from lattice QCD, Phys. Rev. D102 (2020) 094518 [2007.06957]
arXiv 2020
- [1]
-
[2]
L.J. Cooper, C.T. Davies, J. Harrison, J. Komijani and M. Wingate,𝐵𝑐→𝐵𝑠(𝑑) form factors from lattice QCD, Phys. Rev. D102 (2020) 014513 [2003.00914]
arXiv 2020
-
[4]
J. Harrison and C.T.H. Davies,𝐵𝑠→𝐷∗ 𝑠 form factors for the full𝑞2 range from lattice QCD, Phys. Rev. D105 (2022) 094506 [2105.11433]
arXiv 2022
-
[5]
A. Bazavov et al.,Fermilab Lattice, MILC collaboration, Semileptonic form factors for 𝐵→𝐷∗ℓ𝜈 at nonzero recoil from2+ 1-flavor lattice QCD,Eur. Phys. J. C82(2022) 1141 [2105.14019]
arXiv 2022
-
[6]
L.J. Cooper, C.T.H. Davies and M. Wingate,Form factors for the processes𝐵+ 𝑐→𝐷0ℓ+𝜈ℓ and𝐵+ 𝑐→𝐷+ 𝑠ℓ+ℓ+(𝜈 ¯𝜈) from lattice QCD, Phys. Rev. D105 (2022) 014503 [2108.11242]
work page Pith review arXiv 2022
-
[7]
B. Colquhoun, S. Hashimoto, T. Kaneko and J. Koponen,JLQCD collaboration,Form factors of B→𝜋ℓ𝜈 and a determination of |Vub| with Möbius domain-wall fermions, Phys. Rev. D106 (2022) 054502 [2203.04938]
arXiv 2022
-
[8]
W.G. Parrott, C. Bouchard and C.T.H. Davies,B→K and D→K form factors from fully relativistic lattice QCD, Phys. Rev. D107 (2023) 014510 [2207.12468]
arXiv 2023
Show all 31 references
-
[9]
Flynn, R.C
J.M. Flynn, R.C. Hill, A. Jüttner, A. Soni, J.T. Tsang and O. Witzel,RBC/UKQCD collaboration, Exclusive semileptonic Bs→Kℓ𝜈 decays on the lattice,Phys. Rev. D107 (2023) 114512 [2303.11280]
2023 arXiv
-
[10]
Harrison and C.T.H
J. Harrison and C.T.H. Davies,𝐵→𝐷∗ and𝐵𝑠→𝐷∗ 𝑠 vector, axial-vector and tensor form factors for the full q2 range from lattice QCD,Phys. Rev. D109 (2024) 094515 [2304.03137]
2024 arXiv
-
[11]
Y. Aoki, B. Colquhoun, H. Fukaya, S. Hashimoto, T. Kaneko, R. Kellermann et al.,JLQCD collaboration, B→D*ℓ𝜈ℓ semileptonic form factors from lattice QCD with Möbius domain-wall quarks,Phys. Rev. D109 (2024) 074503 [2306.05657]
2024 arXiv
-
[12]
Aoki et al.,Flavour Lattice A veraging Group (FLAG)collaboration,FLAG Review 2024, 2411.04268
Y. Aoki et al.,Flavour Lattice A veraging Group (FLAG)collaboration,FLAG Review 2024, 2411.04268
2024 arXiv
-
[13]
Bordone and A
M. Bordone and A. Juttner,New strategies for probing𝐵→𝐷∗ℓ ¯𝜈ℓ lattice and experimental data, 2406.10074
-
[14]
C. Boyd, B. Grinstein and R.F. Lebed,Precision corrections to dispersive bounds on form-factors,Phys. Rev. D56 (1997) 6895 [hep-ph/9705252]. 9 Update of HPQCD𝐵𝑐→𝐽/𝜓 Form Factors Judd Harrison
1997 arXiv
-
[15]
Aaij et al.,LHCb collaboration, Measurement of the ratio of branching fractions B(𝐵+ 𝑐 → 𝐽/𝜓𝜏+𝜈𝜏)/B(𝐵+ 𝑐 → 𝐽/𝜓𝜇+𝜈𝜇), Phys
R. Aaij et al.,LHCb collaboration, Measurement of the ratio of branching fractions B(𝐵+ 𝑐 → 𝐽/𝜓𝜏+𝜈𝜏)/B(𝐵+ 𝑐 → 𝐽/𝜓𝜇+𝜈𝜇), Phys. Rev. Lett.120 (2018) 121801 [1711.05623]
2018 arXiv
-
[16]
CMScollaboration, measurement of R(Jpsi) in the hadronic tau decay channel, CMS-PAS-BPH-23-001(2024)
2024
-
[17]
Bazavov et al.,MILCcollaboration, Lattice QCD Ensembles with Four Flavors of Highly Improved Staggered Quarks, Phys
A. Bazavov et al.,MILCcollaboration, Lattice QCD Ensembles with Four Flavors of Highly Improved Staggered Quarks, Phys. Rev. D87(2013) 054505 [1212.4768]
2013 arXiv
-
[18]
Bazavov et al.,MILCcollaboration, Scaling studies of QCD with the dynamical HISQ action,Phys
A. Bazavov et al.,MILCcollaboration, Scaling studies of QCD with the dynamical HISQ action,Phys. Rev. D82(2010) 074501 [1004.0342]
2010 arXiv
-
[19]
Borsanyi et al.,High-precision scale setting in lattice QCD,JHEP 09(2012) 010 [1203.4469]
S. Borsanyi et al.,High-precision scale setting in lattice QCD,JHEP 09(2012) 010 [1203.4469]
2012 arXiv
-
[20]
Dowdall, C.T.H
R.J. Dowdall, C.T.H. Davies, G.P. Lepage and C. McNeile,Vus from pi and K decay constants in full lattice QCD with physical u, d, s and c quarks, Phys. Rev. D88 (2013) 074504 [1303.1670]
2013 arXiv
-
[21]
Chakraborty, C.T.H
B. Chakraborty, C.T.H. Davies, P.G. de Oliviera, J. Koponen, G.P. Lepage and R.S. Van de Water,The hadronic vacuum polarization contribution to𝑎𝜇 from full lattice QCD,Phys. Rev. D96(2017) 034516 [1601.03071]
2017 arXiv
-
[22]
Chakraborty, C.T.H
B. Chakraborty, C.T.H. Davies, B. Galloway, P. Knecht, J. Koponen, G.C. Donald et al., High-precision quark masses and QCD coupling from𝑛𝑓 = 4 lattice QCD,Phys. Rev. D91 (2015) 054508 [1408.4169]
2015 arXiv
-
[23]
Hatton, C.T.H
D. Hatton, C.T.H. Davies, B. Galloway, J. Koponen, G.P. Lepage and A.T. Lytle, Charmonium properties from lattice𝑄𝐶𝐷+QED : Hyperfine splitting,𝐽/𝜓 leptonic width, charm quark mass, and𝑎𝑐 𝜇, Phys. Rev. D102 (2020) 054511 [2005.01845]
2020 arXiv
-
[24]
Bazavov et al.,𝐵- and𝐷-meson leptonic decay constants from four-flavor lattice QCD, Phys
A. Bazavov et al.,𝐵- and𝐷-meson leptonic decay constants from four-flavor lattice QCD, Phys. Rev. D98 (2018) 074512 [1712.09262]
2018 arXiv
-
[25]
Lepage,corrfitter,corrfitter Version 8.0.2
G.P. Lepage,corrfitter,corrfitter Version 8.0.2
-
[26]
Gubernari, M
N. Gubernari, M. Reboud, D. van Dyk and J. Virto,Dispersive analysis of B→ K(∗) and B𝑠→𝜙 form factors,JHEP12 (2023) 153 [2305.06301]
2023 arXiv
-
[27]
Hatton, C
D. Hatton, C. Davies, G. Lepage and A. Lytle,Renormalizing vector currents in lattice QCD using momentum-subtraction schemes, Phys. Rev. D100 (2019) 114513 [1909.00756]
2019 arXiv
-
[28]
Sharpe and A
S.R. Sharpe and A. Patel,Perturbative corrections for staggered four fermion operators, Nucl. Phys. B417 (1994) 307 [hep-lat/9310004]
1994 arXiv
-
[29]
Gracey,Three loopMS tensor current anomalous dimension in QCD,Phys
J.A. Gracey,Three loopMS tensor current anomalous dimension in QCD,Phys. Lett. B488 (2000) 175 [hep-ph/0007171]. 10 Update of HPQCD𝐵𝑐→𝐽/𝜓 Form Factors Judd Harrison
2000 arXiv
-
[30]
Harrison, ¯𝑏𝑐 susceptibilities from fully relativistic lattice QCD,Phys
J. Harrison, ¯𝑏𝑐 susceptibilities from fully relativistic lattice QCD,Phys. Rev. D110 (2024) 054506 [2405.01390]
2024 arXiv
-
[31]
McLean, C.T.H
E. McLean, C.T.H. Davies, A.T. Lytle and J. Koponen,Lattice QCD form factor for 𝐵𝑠→𝐷∗ 𝑠𝑙𝜈 at zero recoil with non-perturbative current renormalisation, Phys. Rev. D99 (2019) 114512 [1904.02046]. 11
2019 arXiv
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.