REVIEW 2 major objections 8 minor 22 references
Towards more accurate $B_{(s)}\rightarrow\pi(K)$ and $D_{(s)}\rightarrow\pi(K)$ form factors
T0 review · 2 major / 8 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read The paper claims that the heavy-HISQ lattice QCD approach, with heavy quark masses from charm to near bottom, covers nearly the full kinematic range for the four decays $B\to\pi$, $B_s\to K$, $D\to\pi$, and $D_s\to K$, and that a chained…
desk verdict A competent, honest status report from a group extending a proven method to new channels; the main open question is whether the chained-fit split preserves all needed correlations. 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 central machinery is the heavy-HISQ formalism: fully relativistic highly improved staggered quarks for the heavy quark, with local scalar, vector, and tensor currents. The analysis relies on global correlator fits with Bayesian priors, implemented with gvar, lsqfit, and corrfitter, and uses a chained fitting procedure to handle the large amount of correlator data, splitting parameters into subsets based on their statistical correlations. The final extrapolation to the continuum, physical quark masses, and full $q^2$ range uses the modified $z$-expansion.
What would settle it
Compute the same form factor matrix elements on the f-5 ensemble using an independent renormalization procedure, such as the conserved vector current or a momentum-subtraction scheme, and compare with the results using the adopted $Z$ factors; a difference larger than the quoted uncertainties would indicate the renormalization assumption fails. Alternatively, test the kinematic identity $f_+(0)=f_0(0)$ at $q^2=0$ after the $z$-expansion, which must hold exactly; a violation beyond errors would signal a problem with the fitting or renormalization.
Extended reading notes
Core claim
The paper advances the claim that the heavy-HISQ approach, combined with HISQ valence quarks on MILC ensembles, can determine the form factors for $H\to\pi$ and $H_s\to K$ over nearly the entire kinematic range, where $H$ is a generic heavy meson whose heavy quark mass is tuned from charm to near bottom. The key evidence is the successful chained global fit to two- and three-point correlators, which models non-oscillating and oscillating states, and the resulting preliminary $f_0$, $f_+$, and $f_T$ results for $H\to\pi$ on the f-5 ensemble. The paper argues that this method, already demonstrated to give precise $B\to K$ and $D\to K$ form factors, will extend the same precision to the pion and kaon final states needed for $|V_{ub}|$ and $|V_{cd}|$.
Load-bearing premise
The renormalization factors for the scalar, vector, and tensor currents, taken from earlier calculations, are assumed to be correct for all heavy-quark masses and for the pion and kaon final states; if they are not, every form factor carries an unknown systematic shift.
Editorial extensions
If this is right
- If the chained fits succeed on all ensembles, the project will produce continuum-limit, physical-mass form factors for $B\to\pi$, $B_s\to K$, $D\to\pi$, and $D_s\to K$ over the full $q^2$ range.
- These form factors, combined with experimental decay rates, will yield exclusive determinations of $|V_{ub}|$ and $|V_{cd}|$.
- The nearly full kinematic coverage reduces the model dependence of the $z$-expansion extrapolation compared to calculations that only cover a limited $q^2$ range.
- The success of the chained fitting procedure demonstrates a practical way to handle the large correlator datasets required for multi-channel, multi-mass lattice calculations.
Reading between the lines
- If the anticipated precision materializes, the exclusive $|V_{ub}|$ derived here could sharpen the comparison with inclusive determinations and shed light on the long-standing inclusive-exclusive tension.
- The observed correlation grouping of the four currents suggests that scalar and temporal-vector parameters share much of their statistical fluctuation, so the renormalization uncertainties of these currents may be partially correlated; a combined analysis could exploit this.
- The method is naturally extensible to other light final states, such as $K^*$ or heavier resonances, as long as the correlator signals remain controllable, potentially broadening the set of decays accessible to the same formalism.
- The physical-mass ensembles included in the ongoing work will provide a direct check of the light-quark mass dependence, testing whether the 5:1 strange-to-light ratio results extrapolate consistently.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This proceedings paper reports work in progress toward lattice QCD calculations of scalar, vector, and tensor form factors for B→π, B_s→K, D→π, and D_s→K. The calculation uses HISQ valence quarks on MILC ensembles at several lattice spacings, including 0.044 fm, with heavy-HISQ heavy quarks ranging from charm to near bottom. The paper describes the two- and three-point correlator fit ansatz, the prior construction procedure, and a chained-fit strategy in which the data are split by current type. It shows a representative reconstruction check and preliminary f-5 ensemble form-factor plots. Section 3 explicitly states that the global fitting methodology is not yet finalized, so the numerical results are presented as preliminary.
Significance. If the program succeeds, it would provide form factors over nearly the full q^2 range for four phenomenologically important decays and could improve determinations of |V_ub| and |V_cd|. The paper's strengths are the broad ensemble coverage (including physical and 5:1 light masses, and lattice spacings down to 0.044 fm), the systematic use of the established heavy-HISQ framework of refs. [5,6], the explicit description of prior construction including Gaussian Bayes factor optimization, and the transparent presentation of preliminary reconstruction checks. Because the fitting methodology is not finalized (Section 3) and no continuum or physical-mass extrapolation is presented, the contribution should be read as a status report rather than as a source of final numerical values.
major comments (2)
- [Section 2.2, Fig. 1] The split of the chained fit into two groups (scalar plus temporal-vector, spatial-vector plus tensor) is justified only by the representative correlation matrix in Fig. 1. That figure shows that within-group correlations are high, but it does not demonstrate that cross-link correlations are negligible. Because all four currents are connected through the same H and pi energies and amplitudes in Eq. (3), the two groups are not obviously independent. The statement that this procedure 'preserves correlations between the various subsets of data' is therefore not established. Please quantify the cross-link correlations, and compare the chained-fit posterior with a full simultaneous fit on at least one ensemble (or a reduced data set), reporting differences in J^{nn}_{00} and its uncertainty. Without such a check, the reconstruction in Fig. 4 and the preliminary form factors in Fig. 5 cannot be taken as evidence that the method is unbiased.
- [Section 2.1, Eqs. (4)-(7)] The renormalization factors Z_disc, Z_V, and Z_T(μ) are taken from refs. [12-15], which were computed in other heavy-HISQ calculations. The manuscript does not state whether these factors are applicable to H→π and H_s→K for the full range of heavy masses used here, nor does it give an estimate of the systematic error if they are not. Since these factors enter every matrix element in Eq. (4), please state their provenance for these channels and masses, or provide an estimate of the induced uncertainty. This may be acceptable for a proceedings paper if the final publication will address it, but it should be made explicit.
minor comments (8)
- [Section 2, first paragraph] The phrase 'global (one fit per ensemble)' appears to conflict with the chained-fit procedure described in Section 2.2; please clarify whether 'global' refers to all data in an ensemble or to a single link.
- [Section 2.2, Fig. 1 caption] The caption lists only two twists (θ=0.0, 0.4281) while the text mentions five twists; please explain why only these are shown.
- [Figure 5 caption] The axes and the definition of the error bands are not given in the caption; please add them so the plot can be interpreted independently.
- [Section 2.3] The text 'denote the prior central value and uncertainty ... as \tilde{P}_i + \tilde{\sigma}_{P_i}' should read '\tilde{P}_i \pm \tilde{\sigma}_{P_i}' or 'with uncertainty \tilde{\sigma}_{P_i}'.
- [Throughout] There are several typos, including 'York Uiversity' in the author affiliation, 'apprach' in Section 3, and 'spacial' in the discussion of Eq. (6); these should be corrected.
- [Table 1] The column header for the separation widths is only 'T'; please define the unit (e.g., lattice time slices or physical time) in the table caption or a footnote.
- [Section 2.1, Eqs. (5)-(7)] The symbols M_H, M_π, p^μ, and the renormalization scale μ are not all defined with their lattice versus continuum conventions; a short definition sentence would improve readability.
- [Section 3] The explicit limitation that the global fitting methodology is not finalized is stated only near the end; the abstract and introduction should carry a similar caveat so that the preliminary nature of the results is clear from the outset.
Circularity Check
No circularity found: the form-factor extraction is a direct fit to lattice correlator data with broad priors and externally computed renormalization constants.
full rationale
The paper is a progress report whose central quantitative outputs (Figs. 4 and 5) are obtained by fitting two- and three-point correlators to the spectral decompositions in Eqs. (2) and (3), then converting the fitted ground-state current amplitude J_nn^00 into form factors via Eqs. (4)-(7). The fit parameters are determined by the correlator data themselves; the only external inputs are the renormalization factors Z_disc, Z_V, and Z_T(mu) taken from refs. [12-15], which are independent published quantities not derived in this paper and not computed by the present authors. The prior methodology in Sec. 2.3 is explicitly designed not to constrain the answers: prior widths are set to at least ten times the corresponding posterior widths, and the physics-based priors concern pion/kaon masses and dispersion relations rather than the target form factors. The chained-fit grouping in Sec. 2.2 is a statistical approximation based on a correlation-matrix inspection; it could affect the reliability of uncertainties if cross-link correlations are non-negligible, but it does not reduce any predicted quantity to an input by construction. The abstract's claim of 'nearly full coverage of the kinematic range' is a direct property of the chosen heavy-quark masses and simulated daughter momenta, not a fitted prediction. No uniqueness theorem, ansatz, or known result is imported from the authors' prior work in a load-bearing way; refs. [5,6] are ordinary methodological citations. Accordingly, no circular step is exhibited.
Assumptions & free parameters
free parameters (2)
- Prior width inflation factor =
at least 10x posterior uncertainty
- Number of exponential states N_exp =
4
assumptions (4)
- domain assumption HISQ action on MILC Nf=2+1+1 ensembles correctly describes QCD at the simulated lattice spacings.
- domain assumption The renormalization factors Z_disc, Z_V, Z_T from refs [12-15] are valid for the local currents used with these heavy quark masses and channels.
- standard math The modified z-expansion (ref [8]) is model-independent and will control the extrapolation to the full kinematic range.
- domain assumption The chained fit method [10] with the chosen grouping preserves the relevant statistical correlations.
Cite this review
Pith. "Pith review of Towards more accurate $B_{(s)}\rightarrow\pi(K)$ and $D_{(s)}\rightarrow\pi(K)$ form factors." pith.science (2026). https://pith.science/paper/5YDQ7BBF
@misc{pith2026250118586,
author = {Pith},
title = {Pith review of: Towards more accurate $B_(s)\rightarrow\pi(K)$ and $D_(s)\rightarrow\pi(K)$ form factors},
year = {2026},
howpublished = {\url{https://pith.science/paper/5YDQ7BBF}},
note = {Machine review of arXiv:2501.18586}
}
abstract
We present progress on the calculation of scalar, vector, and tensor form factors for the following meson decays: $B\rightarrow\pi$, $B_s\rightarrow K$, $D\rightarrow\pi$ and $D_s\rightarrow K$. This calculation uses the MILC HISQ gluon field ensembles with HISQ valence quarks. We generate ensembles of correlator data with varying lattice spacings, some as small as 0.044 fm. Some ensembles have a strange-to-light quark mass ratio of 5:1 and others use the physical light quark mass. The fully-relativistic, heavy-HISQ approach is used for the heavy quark, with simulation masses ranging from the charm to near the bottom. This heavy-HISQ approach provides nearly full coverage of the kinematic range.
Figures
Figures from the paper (2 more)
Reference graph
Works this paper leans on
-
[1]
B Meson Semileptonic Form Factors from Unquenched Lattice QCD
HPQCDCollaboration, E. Gulez, A. Gray, M. Wingate, C. T. H. Davies, G. P. Lepage, and J. Shigemitsu,𝐵 meson semileptonic form factors from unquenched lattice QCD, Phys. Rev. D 73(2006) 074502, [hep-lat/0601021]
work page Pith review arXiv 2006
-
[2]
E. Follana, Q. Mason, C. Davies, K. Hornbostel, G. P. Lepage, J. Shigemitsu, H. Trottier, and K. Wong,Highly improved staggered quarks on the lattice with applications to charm physics, Phys. Rev. D75 (2007) 054502, [hep-lat/0610092]
arXiv 2007
-
[3]
A. Bazavov, C. Bernard, C. DeTar, W. Freeman, S. Gottlieb, U. M. Heller, J. E. Hetrick, J.Laiho,L.Levkova, M.Oktay,J.Osborn,R.L.Sugar, D.Toussaint,andR.S.VandeWater, Scaling studies of QCD with the dynamical highly improved staggered quark action, Phys. Rev. D82(2010) 074501, [1004.0342]
arXiv 2010
-
[4]
A. Bazavov, C. Bernard, J. Komijani, C. DeTar, L. Levkova, W. Freeman, S. Gottlieb, R. Zhou, U. M. Heller, J. E. Hetrick, J. Laiho, J. Osborn, R. L. Sugar, D. Toussaint, and R. S. Van de Water,Lattice QCD ensembles with four flavors of highly improved staggered quarks, Phys. Rev. D87(2013) 054505, [1212.4768]
arXiv 2013
-
[5]
W. G. Parrott, C. Bouchard, C. T. H. Davies, and D. Hatton,Toward accurate form factors for 𝐵-to-light meson decay from lattice QCD, Phys. Rev. D103 (2021) 094506, [2010.07980]
arXiv 2021
-
[6]
W. G. Parrott, C. Bouchard, and C. T. H. Davies,𝐵→𝐾 and𝐷→𝐾 form factors from fully relativistic lattice QCD, Phys. Rev. D107 (2023) 014510, [2207.12468]
arXiv 2023
-
[7]
Navas et al.,Review of particle physics,Phys
Particle Data GroupCollaboration, S. Navas et al.,Review of particle physics,Phys. Rev. D 110 (2024) 030001
work page 2024
-
[8]
C. Bourrely, L. Lellouch, and I. Caprini,Model-independent description of𝐵→𝜋𝑙𝜈 decays and a determination of|𝑉𝑢𝑏|,Phys. Rev. D79(2009) 013008, [0807.2722]
arXiv 2009
Show all 22 references
-
[9]
Lepage, C
P. Lepage, C. Gohlke, and D. Hackett,gplepage/gvar: gvar version 13.1.1, 2024
2024
-
[10]
Lepage and C
P. Lepage and C. Gohlke,gplepage/lsqfit: lsqfit version 13.2.3, 2024
2024
-
[11]
Lepage,gplepage/corrfitter: corrfitter version 8.2, 2021
P. Lepage,gplepage/corrfitter: corrfitter version 8.2, 2021
2021
-
[12]
Monahan, J
C. Monahan, J. Shigemitsu, and R. Horgan,Matching lattice and continuum axial-vector and vector currents with nonrelativistic QCD and highly improved staggered quarks,Phys. Rev. D 87(2013) 034017, [1211.6966]
2013 arXiv
-
[13]
H. Na, C. T. H. Davies, E. Follana, G. P. Lepage, and J. Shigemitsu,𝐷→𝐾,ℓ𝜈 semileptonic decay scalar form factor and𝑉𝑐𝑠 from lattice QCD, Phys. Rev. D82 (2010) 114506, [1008.4562]
2010 arXiv
-
[14]
Koponen, C
J. Koponen, C. T. H. Davies, G. C. Donald, E. Follana, G. P. Lepage, H. Na, and J. Shigemitsu,The shape of the𝐷→𝐾 semileptonic form factor from full lattice QCD and 𝑉𝑐𝑠, 1305.1462. 9 Towards more accurate𝐵(𝑠)→𝜋(𝐾) and𝐷(𝑠)→𝜋(𝐾) form factors Logan Roberts
-
[15]
Hatton, C
HPQCDCollaboration, D. Hatton, C. T. H. Davies, G. P. Lepage, and A. T. Lytle, Renormalization of the tensor current in lattice QCD and the𝐽/𝜓 tensor decay constant, Phys. Rev. D102(2020) 094509, [2008.02024]
2020 arXiv
-
[16]
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]
2020 arXiv
-
[17]
Lepage, B
G. Lepage, B. Clark, C. Davies, K. Hornbostel, P. Mackenzie, C. Morningstar, and H. Trottier,Constrained curve fitting,Nucl. Phys. B Proc. Suppl.106-107(2002) 12, [hep-lat/0110175]
2002 arXiv
-
[18]
Sivia and J
D. Sivia and J. Skilling,Data Analysis: A Bayesian Tutorial. Oxford science publications. OUP Oxford, 2006
2006
-
[19]
Carlin and T
B. Carlin and T. Louis,Bayes and Empirical Bayes Methods for Data Analysis, Second Edition. Chapman & Hall/CRC Texts in Statistical Science. Taylor & Francis, 2010
2010
-
[20]
R. E. Kass and A. E. Raftery,Bayes factors, Journal of the American Statistical Association 90 (1995) 773
1995
-
[21]
McLean, C
HPQCDCollaboration, E. McLean, C. T. H. Davies, J. Koponen, and A. T. Lytle, 𝐵𝑠→𝐷𝑠ℓ𝜈 form factors for the full𝑞2 range from lattice QCD with nonperturbatively normalized currents, Phys. Rev. D101 (2020) 074513, [1906.00701]
2020 arXiv
-
[22]
Harrison, C
J. Harrison, C. T. Davies, and A. Lytle,𝐵𝑐→𝐽/𝜓 form factors for the full𝑞2 range from lattice QCD, Phys. Rev. D102 (2020) 094518, [2007.06957]. 10
2020 arXiv
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.