REVIEW 3 major objections 4 minor 22 references
Study on the $P$-wave form factors contributing to $ B_s $ to $D_s$ inclusive semileptonic decays from lattice simulations
T0 review · 3 major / 4 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read This pilot lattice QCD study extracts the P-wave form factors of $B_s$ to $D_s$ semileptonic decays from four-point correlators and finds that the ground-state doublet leaves a deficit in the Uraltsev sum rule, evidence that radial…
desk verdict A genuine method paper: four-point correlators can be used to extract exclusive P-wave form factors, but the Uraltsev deficit claim is not yet robust because single-exponential fits and external inputs are not fully controlled. 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 object is the four-point correlator $C_{J_\mu J_\nu}(\mathbf{q},t)$, related to the forward hadronic tensor by a Laplace transform; inserting a complete set of states casts it as a sum of exponentials with form-factor prefactors. The load-bearing identity is the Uraltsev sum rule, $\sum_n (|\tau_{3/2}^{(n)}(1)|^2 - |\tau_{1/2}^{(n)}(1)|^2) = 1/4$, which is the benchmark the extracted ground-state values are compared against. The conversion of the lattice numbers $g_+(1)$ and $f_{V1}(1)$ into $\tau_{1/2}$ and $\tau_{3/2}$ runs through first-order HQET relations involving the mass splittings $\Lambda_s^* - \Lambda_s$ and $\Lambda_s' - \Lambda_s$.
What would settle it
Compute the same zero-recoil four-point correlator decomposition on a finer lattice with physical quark masses, extracting the P-wave form factors and also determining the mass splittings $\Lambda_s^* - \Lambda_s$ and $\Lambda_s' - \Lambda_s$ from the lattice spectrum itself. If $\tau_{3/2}(1)^2 - \tau_{1/2}(1)^2$ then agrees with $1/4$ within errors, the radial-excitation interpretation would be falsified; alternatively, a direct three-point correlator measurement of $\tau_{1/2}(1)$ and $\tau_{3/2}(1)$ with properly constructed P-wave interpolating operators could settle the same question.
Extended reading notes
Core claim
The central discovery claim is that the Euclidean four-point correlator $C_{J_\mu J_\nu}(\mathbf{q},t)$, after inserting a complete set of final states, is a sum of exponentials whose prefactors are the exclusive hadronic matrix elements, so a multi-exponential fit can extract the P-wave form factors without constructing excited-state interpolating operators. At zero recoil the paper obtains $|g_+(1)| = 0.166 \pm 0.049$ and $|f_{V1}(1)| = 0.445 \pm 0.055$ from the $D_{s0}^*$ and $D_{s1}$ channels, and through the heavy-quark relations (19) and (20) converts them to $|\tau_{1/2}(1)| = 0.35 \pm 0.10$ and $|\tau_{3/2}(1)| = 0.423 \pm 0.052$. The resulting difference $\tau_{3/2}(1)^2 - \tau_{1/2}(1)^2 = 0.053 \pm 0.079$ is below the Uraltsev sum rule value $1/4$, and the paper takes this deficit as evidence that radial excitations contribute significantly to the sum rule, while noting the single-ensemble pilot nature of the study.
Load-bearing premise
The numerical values of the Isgur-Wise form factors, and hence the size of the Uraltsev deficit, depend on assumed heavy-quark-limit mass splittings (0.28 GeV and 0.41 GeV) and quark masses ($m_b = 4.8$ GeV, $m_c = 1.1$ GeV) taken from a phenomenological analysis, with no uncertainty propagated on these inputs.
Editorial extensions
If this is right
- The same four-point correlator data used for inclusive decay rates can be re-fitted to obtain exclusive excited-state form factors, providing a common lattice handle on both sides of the inclusive-exclusive $V_{cb}$ comparison.
- The extracted ground-state P-wave form factors do not saturate the Uraltsev sum rule, so radial excitations of the $D_s$ system must carry a non-negligible fraction of the rule's strength.
- The zero-recoil suppression of $g_{V1}$ relative to $f_{V1}$, predicted by the HQET factor $(\epsilon_c - 3\epsilon_b)$, is confirmed and justifies neglecting the $D_{s1}'$ contribution in the correlator decomposition.
- The method reproduces the known mass hierarchy of the excited $D_s$ states from the same single-exponential fits, with splittings consistent with experiment within errors.
- At non-zero recoil the approach can in principle constrain the slopes $\tau'$ and $\zeta'$ of the Isgur-Wise form factors, although the current errors ($\tau' = -1 \pm 17$, $\zeta' = -18 \pm 10$) are too large to be decisive.
Reading between the lines
- If the Uraltsev deficit survives on finer ensembles with physical quark masses and propagated input uncertainties, the $1/2$-versus-$3/2$ puzzle may be resolved by radial-excitation contributions within the ground-state doublet's own sector, rather than by $1/m_Q$ corrections.
- The same four-point correlator decomposition is directly transferable to other weak transitions, such as $B \to D^*$ or charmed-meson decays, where constructing orbitally excited interpolating operators is even harder.
- Because the conversion to $\tau$ values currently relies on external phenomenological mass splittings, a version of this analysis that extracts $\Lambda_s^* - \Lambda_s$ and $\Lambda_s' - \Lambda_s$ from the same lattice spectrum would turn the sum-rule test into a fully internal lattice calculation.
- Jointly analyzing the inclusive width and the exclusive form factors from one set of four-point correlators could provide a data-driven check of quark-hadron duality at finite Euclidean time, rather than assuming saturation channel by channel.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a pilot lattice study using B_s four-point correlators on a single 24^3 x 64 ensemble with a = 0.11 fm to extract P-wave form factors of semileptonic B_s -> D_s transitions. At zero recoil, single-exponential fits to C_A0A0 and C_V||V|| yield |g_+(1)| = 0.166 ± 0.049 and |f_V1(1)| = 0.445 ± 0.055, which are converted via HQET relations into Isgur-Wise form factors |tau_1/2(1)| = 0.35 ± 0.10 and |tau_3/2(1)| = 0.423 ± 0.052. The combination tau_3/2(1)^2 - tau_1/2(1)^2 = 0.053 ± 0.079 is compared with the Uraltsev sum rule value 1/4, and the deficit is interpreted as evidence for significant radial-excitation contributions. The paper also reports mass splittings consistent with PDG and a non-zero recoil analysis using "Approximation A," which yields only poorly constrained slope parameters.
Significance. If the results hold, the paper demonstrates a new route to exclusive P-wave form factors from four-point correlators, potentially connecting inclusive and exclusive determinations of V_cb and addressing the 1/2-versus-3/2 puzzle. The work is not circular: the tau values are computed from lattice form factors with external HQET inputs, and the comparison to 1/4 is an external benchmark, not a fit target. The zero-recoil formalism in Secs. 2-3 is clear, and the masses are consistent with experiment within large errors. However, the central quantitative claim rests on the ground-state identification of the fitted amplitudes and on unpropagated external inputs, so the significance of the specific numerical deficit is not yet established.
major comments (3)
- [Sec. 5, Eqs. (16)-(17), Fig. 1] The central claim that radial excitations contribute significantly to the Uraltsev sum rule rests on the zero-recoil amplitudes |g_+(1)| and |f_V1(1)|, extracted from single-exponential fits to C_A0A0 and C_V||V||. These correlators contain, in principle, an infinite sum over radial excitations with the same J^P, and the paper does not show a t_min scan, a two-exponential fit, or any other demonstration that the fitted amplitudes correspond to the n=0 states. The quoted P-wave mass splittings have large errors (e.g., M_Ds0* - M_Ds = 420 ± 160 MeV), so the fitted energies do not certify amplitude purity. If excited-state contamination biases |tau_1/2| upward, the deficit tau_3/2^2 - tau_1/2^2 would be artificially suppressed, undermining the interpretation in terms of radial-excitation contributions. This should be addressed before the deficit is used as evidence.
- [Sec. 5, Eqs. (19)-(20)] The conversion of the lattice results into tau_1/2(1) and tau_3/2(1) uses the external inputs Lambda*_s - Lambda_s = 0.28 GeV, Lambda'_s - Lambda_s = 0.41 GeV, m_b = 4.8 GeV, and m_c = 1.1 GeV, but the uncertainties on these inputs are not propagated. The quoted errors on tau_1/2(1), tau_3/2(1), and on the deficit tau_3/2(1)^2 - tau_1/2(1)^2 = 0.053 ± 0.079 are therefore not the full uncertainty, and the size of the deficit could shift appreciably under reasonable variations of these parameters. A systematic error estimate for these external inputs is needed for the main quantitative conclusion.
- [Sec. 5, Eqs. (17) and (21)] The approximation in Eq. (17) neglects the D'_s1 contribution to C_V||V|| based on the HQET expectation that g_V1 is suppressed. The numerical check in Eq. (21), however, yields |g_V1(1)| = 0.0282 ± 0.083, whose error is much larger than the central value and therefore does not quantitatively validate the neglect. Please provide an upper bound on the neglected contribution or include both 1+ states in the fit, since a non-negligible D'_s1 amplitude would bias the extracted |f_V1(1)| and hence tau_3/2(1).
minor comments (4)
- [Fig. 2 caption] The caption contains 'A^{D_{s9}}_{A∥A∥}', which should read 'A^{D^*_{s0}}_{A∥A∥}'.
- [Sec. 2 and Sec. 5] The text states that exclusive information can be extracted by multiple-exponential fits, but the zero-recoil analysis in Sec. 5 uses single-exponential fits; please state the fit ranges and any plateau analysis that justifies this choice.
- [Sec. 5, after Eq. (21)] The value |g_V1(1)| = 0.0282 ± 0.083 has an error larger than the central value; reporting an upper bound or confidence interval would be more informative.
- [Abstract and Sec. 6] With a deficit of 0.053 ± 0.079 (about 2.5 sigma from 1/4), the wording 'suggest significant contributions' is stronger than the present statistics warrant; consider 'tentatively suggest' or 'are consistent with significant contributions'.
Circularity Check
No significant circularity: the Isgur-Wise form factors are computed from lattice correlator amplitudes via HQET relations with external inputs, and the Uraltsev sum-rule comparison is an external benchmark, not a fit target.
full rationale
The derivation chain is self-contained with respect to the paper's central claim. The extracted zero-recoil amplitudes |g_+(1)| and |f_V1(1)| come from single-exponential fits to the lattice four-point correlators (Eqs. 16-17). The Isgur-Wise quantities |tau_1/2(1)| and |tau_3/2(1)| are then obtained by dividing these amplitudes by assumed mass splittings and quark masses (Eqs. 19-20), taken from the external phenomenological analysis [22]; they are not fitted to the Uraltsev sum rule. The comparison of tau_3/2(1)^2 - tau_1/2(1)^2 to 1/4 is an external benchmark, so the reported deficit is a measurement, not a construction. The nonzero-recoil prefactors labeled 'predicted by Approx. A' in Fig. 2 are consistency checks of the approximation against independently extracted amplitudes, not independent predictions used as load-bearing evidence; they are presented with large uncertainties and no decisive conclusion is drawn. Self-citations to Refs. [3-8] establish the four-point correlator formalism, but this formalism is a published, independently available method, and the new exclusive extraction is not defined in terms of the Uraltsev result. The main limitations - possible contamination from radial excitations in single-exponential zero-recoil fits, and the external mass-splitting inputs - are systematic and input-accuracy risks, not circular dependencies.
Assumptions & free parameters
free parameters (4)
- Lambda*_s - Lambda_s =
0.28 GeV
- Lambda'_s - Lambda_s =
0.41 GeV
- Bottom quark mass m_b =
4.8 GeV
- Charm quark mass m_c =
1.1 GeV
assumptions (5)
- domain assumption HQET leading-order relations in Eqs. (19)-(20) are valid at the lattice quark masses.
- standard math The four-point correlator is saturated by a sum over single-hadron states as in Eq. (4).
- domain assumption Contributions from D'_s1 and D_s2 can be neglected at the current precision.
- domain assumption Approximation A (treating w - 1 as order 1/m_Q) is valid for the small q^2 values used.
- domain assumption One coarse 24^3 x 64 ensemble with a = 0.11 fm and near-physical quark masses is sufficient for a pilot estimate.
Cite this review
Pith. "Pith review of Study on the $P$-wave form factors contributing to $ B_s $ to $D_s$ inclusive semileptonic decays from lattice simulations." pith.science (2026). https://pith.science/paper/QLCX5J6J
@misc{pith2026250119284,
author = {Pith},
title = {Pith review of: Study on the $P$-wave form factors contributing to $ B_s $ to $D_s$ inclusive semileptonic decays from lattice simulations},
year = {2026},
howpublished = {\url{https://pith.science/paper/QLCX5J6J}},
note = {Machine review of arXiv:2501.19284}
}
abstract
We present a pilot study on extracting the form factors of the semileptonic decay of a $ B_s $ meson to the $P$-wave $ D_s^{**} $ states from $B_s$ four-point correlators. With their inclusive nature, four-point correlators include contributions from all possible final states. From the extracted $ P $-wave form factors, we obtain numerical results for the corresponding Isgur-Wise form factors. The results suggest significant contributions from radial excitations to the Uraltsev sum rule at zero-recoil. In this pilot study, a coarse lattice of $ 24^3\times 64 $ with lattice spacing of $0.11\,\mathrm{fm}$ is used for the analysis.
Figures
Reference graph
Works this paper leans on
-
[1]
HFLA Vcollaboration, Averages of b-hadron, c-hadron, and𝜏-lepton properties as of 2021, Phys. Rev. D107(2023) 052008 [2206.07501]
arXiv 2023
-
[2]
I.I. Bigi, B. Blossier, A. Le Yaouanc, L. Oliver, O. Pene, J.C. Raynal et al.,Memorino on the ‘1/2 versus 3/2 puzzle’ in𝐵→𝑙𝜈𝑋𝑐 - a year later and a bit wiser,Eur. Phys. J. C52 (2007) 975 [0708.1621]
arXiv 2007
-
[3]
S. Hashimoto,Inclusive semi-leptonic B meson decay structure functions from lattice QCD, PTEP2017 (2017) 053B03 [1703.01881]
arXiv 2017
-
[4]
P. Gambino and S. Hashimoto,Inclusive Semileptonic Decays from Lattice QCD, Phys. Rev. Lett.125 (2020) 032001 [2005.13730]
arXiv 2020
- [5]
-
[6]
R. Kellermann, A. Barone, S. Hashimoto, A. Jüttnerc and T. Kanekoa,Studies on finite-volume effects in the inclusive semileptonic decays of charmed mesons,PoS LATTICE2023(2024) 272 [2312.16442]
arXiv 2024
- [7]
-
[8]
Updates on inclusive charmed and bottomed meson decays from the lattice
R. Kellermann, A. Barone, S. Hashimoto, A. Jüttner and T. Kaneko,Updates on inclusive charmed and bottomed meson decays from the lattice, in12th International Workshop on the CKM Unitarity Triangle2405.06152
Show all 22 references
-
[9]
Atoui, B
M. Atoui, B. Blossier, V. Morénas, O. Pène and K. Petrov,Semileptonic𝐵→𝐷∗∗ decays in Lattice QCD : a feasibility study and first results,Eur. Phys. J. C75(2015) 376 [1312.2914]
2015 arXiv
-
[10]
Leibovich, Z
A.K. Leibovich, Z. Ligeti, I.W. Stewart and M.B. Wise,Semileptonic B decays to excited charmed mesons, Phys. Rev. D57(1998) 308 [hep-ph/9705467]
1998 arXiv
-
[11]
Isgur and M.B
N. Isgur and M.B. Wise,Excited charm mesons in semileptonic anti-B decay and their contributions to a Bjorken sum rule,Phys. Rev. D43(1991) 819
1991
-
[12]
Uraltsev,New exact heavy quark sum rules,Phys
N. Uraltsev,New exact heavy quark sum rules,Phys. Lett. B501 (2001) 86 [hep-ph/0011124]
2001 arXiv
-
[13]
RBC/UKQCD collaboration, Exclusive semileptonic Bs→Kℓ𝜈 decays on the lattice, Phys. Rev. D107 (2023) 114512 [2303.11280]
2023 arXiv
-
[14]
Shamir,Chiral fermions from lattice boundaries, Nucl
Y. Shamir,Chiral fermions from lattice boundaries, Nucl. Phys. B406 (1993) 90 [hep-lat/9303005]. 9 𝑃-wave form factors of𝐵𝑠 to𝐷𝑠 semileptonic decays Z. Hu
1993 arXiv
-
[15]
Furman and Y
V. Furman and Y. Shamir,Axial symmetries in lattice QCD with Kaplan fermions, Nucl. Phys. B439(1995) 54 [hep-lat/9405004]
1995 arXiv
-
[16]
Brower, H
R.C. Brower, H. Neff and K. Orginos,The Möbius domain wall fermion algorithm, Comput. Phys. Commun.220 (2017) 1 [1206.5214]
2017 arXiv
-
[17]
Y.-G. Cho, S. Hashimoto, A. Jüttner, T. Kaneko, M. Marinkovic, J.-I. Noaki et al.,Improved lattice fermion action for heavy quarks,JHEP05 (2015) 072 [1504.01630]
2015 arXiv
-
[18]
Christ, M
N.H. Christ, M. Li and H.-W. Lin,Relativistic Heavy Quark Effective Action,Phys. Rev. D76 (2007) 074505 [hep-lat/0608006]
2007 arXiv
-
[19]
Lin and N
H.-W. Lin and N. Christ,Non-perturbatively Determined Relativistic Heavy Quark Action, Phys. Rev. D76 (2007) 074506 [hep-lat/0608005]
2007 arXiv
-
[20]
Particle Data Groupcollaboration, Review of particle physics,Phys. Rev. D110 (2024) 030001
2024
-
[21]
Flavour Lattice A veraging Group (FLAG)collaboration, FLAG Review 2024, 2411.04268
2024 arXiv
-
[22]
Bernlochner and Z
F.U. Bernlochner and Z. Ligeti,Semileptonic𝐵(𝑠) decays to excited charmed mesons with 𝑒,𝜇,𝜏 and searching for new physics with𝑅(𝐷∗∗),Phys. Rev. D95(2017) 014022 [1606.09300]. 10
2017 arXiv
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.