Pith. sign in

REVIEW 4 major objections 4 minor 12 references

$B\rightarrow D^{(*)}$ decays from $N_f=2+1+1$ highly improved staggered quarks and clover $b$-quark in the Fermilab interpretation

T0 review · 4 major / 4 minor · reviewed 2026-08-09 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read Preliminary lattice-QCD results show stable form factors for B→D and B→D* decays at non-zero recoil across four lattice spacings, a step toward a precise determination of the CKM element |Vcb|.

desk verdict A competent, transparent progress report from FNAL-MILC on 2+1+1 B->D(*) form factors; the new data are real but the conclusion overclaims systematic control, so the paper deserves referee time with a light touch. read the letter →

arxiv 2501.18969 v1 pith:DSBGWOU6 submitted 2025-01-31 hep-lat

Pietro Butti , Carleton DeTar , Aida El-Khadra , Elvira Gamiz , Steven Gottlieb , William Jay , Hwancheol Jeong , Andreas Kronfeld
show 4 more authors
Jack Laiho Andrew Lytle Pablo Trujillo Alejandro Vaquero
This is my paper · ORCID
classification hep-lat
keywords latticeQCDsemileptonicBdecaysB→D(*)formfactorsFermilabinterpretationHISQactioncloverCKMmatrixelementVcbnon-zerorecoil
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper reports preliminary lattice-QCD results for the semileptonic form factors of $B\to D$ and $B\to D^*$ decays at non-zero recoil, computed on $2+1+1$-flavor ensembles that use highly improved staggered quarks for the sea and light valence quarks and a clover bottom quark in the Fermilab interpretation, a heavy-quark lattice formalism. The data span four lattice spacings, from about $0.15$ fm to $0.06$ fm, and several light-quark masses. The authors aim to show that their two- and three-point correlation-function analysis produces stable form factors with controlled systematics; if true, the same pipeline can be carried through renormalization and continuum extrapolation to support a precise determination of the CKM matrix element $|V_{cb}|$. Because the current renormalization factors were not yet available, the plotted form factors are bare and blinded, and the paper stops short of final renormalized values.

What carries the argument

The load-bearing object is the ratio construction built from three-point correlation functions: ratios such as $R_+$ remove most overlap and kinematic factors and directly produce bare form-factor combinations for $B\to D$, while ratios such as $R_{A_1}$, $X_V$, $X_0$, and $X_1$ do the same for $B\to D^*$. Before fitting, the ratio data are smoothed over adjacent time slices to cancel oscillating staggered-parity contributions. The central fit model is $R(t,T) = F_0 + A(\mathbf{p})e^{-\Delta E_{\rm source} t} + B(\mathbf{p})e^{-\Delta E_{\rm sink}(T-t)}$, with $\Delta E_{D^{(*)}}$ and $\Delta M_B$ shared across momenta; $F_0$ is the ground-state ratio that feeds the form factors. Stability is assessed by varying $t_{\min}$, $t_{\max}$, and the number of states, with Takeuchi Information Criterion weights used for model averaging.

What would settle it

A concrete test: redo the global ratio fits with an additional excited state, for example a second exponential at the source or sink, and with larger source-sink separations; if the ground-state ratio $F_0$ shifts by more than the quoted statistical uncertainty while the fits remain acceptable, the single-exponential model is insufficient and the claimed stability does not guarantee unbiased form factors.

Watch

Extended reading notes

Core claim

The central claim is that the established ratio-based extraction pipeline, now applied to $2+1+1$-flavor ensembles, yields stable, bare $B\to D^{(*)}$ form factors at non-zero recoil. The extraction forms ratios of three-point functions that isolate $h_+$, $h_-$, $h_{A_1}$, $h_{A_2}$, $h_{A_3}$, and $h_V$, then fits the ratio data with a model containing one excited state at the source and one at the sink, with energy splittings shared across momenta. Stability is checked by scanning the lower fit window and the number of states, judging fits by their $p$-values, reduced $\chi^2$, and Takeuchi Information Criterion weights. The paper concludes that the resulting form factors are stable and that systematic effects are well controlled; final renormalization and continuum and physical-mass extrapolations are left to future work.

Load-bearing premise

The analysis assumes that contamination from excited states in the three-point correlators is fully captured by one extra state at the source and one at the sink, with the same energy splittings across momenta; if higher excited states contribute inside the chosen fit windows, the extracted form factors would be biased.

Editorial extensions

If this is right

  • If the stability holds, the same pipeline can be taken through renormalization and continuum extrapolation to produce competitive $B\to D^{(*)}$ form factors at non-zero recoil from $2+1+1$-flavor ensembles.
  • The non-zero recoil data constrain the shape of all six form factors, not only the zero-recoil normalization, which matters for comparing with the measured $w$-dependence of the decay rate.
  • With stable extraction and controlled systematics, the dominant remaining uncertainty in a future $|V_{cb}|$ determination is expected to come from renormalization and the continuum and physical-mass extrapolation rather than from excited-state contamination in the chosen fit windows.
  • The stability checks shown here provide a template for validating three-point correlation-function analyses of other semileptonic heavy-meson decays.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Beyond the paper: the single-exponential excited-state model in the fit could be self-consistent without being complete; a decisive test would be to redo the global fits with a second excited state at source or sink and check whether $F_0$ moves by more than the statistical error.
  • Beyond the paper: once renormalized and extrapolated, these data can be combined with new experimental differential decay-rate measurements to produce a $|V_{cb}|$ value that can be compared with inclusive determinations, sharpening the current exclusive-inclusive tension.
  • Beyond the paper: the ratio construction and stability-scan procedure transfer directly to other channels such as $B_s\to D_s^{(*)}$ decays, so the methodological conclusions here have a wider reach than the specific form factors shown.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. This manuscript is a LATTICE2024 proceedings contribution reporting progress on the FNAL-MILC computation of B→D(*) semileptonic form factors at nonzero recoil using N_f = 2+1+1 HISQ ensembles and a clover bottom quark in the Fermilab interpretation. The authors describe two-point fits (Sec. 3.1) with up to 3+3 states and a stability/model-averaging analysis, three-point ratio fits (Eq. (8), Sec. 3.2) performed simultaneously across momenta and source-sink separations, and the construction of bare form factors (Sec. 3.3, Figs. 4–5). They explicitly state that renormalization factors are not yet included and that no continuum or physical-mass extrapolation has been performed. The conclusion claims stable form factors and good control over systematic effects.

Significance. If completed, this analysis would extend the earlier 2+1-flavor FNAL-MILC form-factor work to the 2+1+1 ensembles, covering several lattice spacings down to about 0.06 fm and including physical light-quark masses, thereby setting up a future high-precision determination of |V_cb|. Strengths of the manuscript are its explicit use of an established pipeline, the detailed two-point stability analysis with TIC model averaging, the blinded presentation of bare form factors, and a clear statement of the quantities that remain to be computed. The current version, however, contains no renormalized form factors, no continuum limit, and no physical-mass extrapolation, and the excited-state control for the three-point fits is substantially less developed than for the two-point fits. It is therefore a status report rather than a final physics result.

major comments (4)
  1. [Sec. 3.2, Eq. (8)] The central fit model, Eq. (8), includes exactly one excited state on each side with momentum-independent splittings ΔE_source and ΔE_sink. The two-point analysis in Sec. 3.1 finds that up to 3+3 states are needed to stabilize the ground state, and the paper presents no equivalent state-count variation for the three-point ratio fits. If higher excited states contribute in the chosen windows (t_min ≈ 0.6 fm and the two source-sink separations T), or if the splittings vary with momentum, the extracted F0 values would be biased in a momentum-dependent way that directly propagates into the form factors. The claim of "good control over systematic effects" in the Conclusion therefore requires either explicit tests of Eq. (8) against models with more excited states or a quantitative demonstration that the one-state truncation is sufficient for the chosen windows.
  2. [Sec. 3.2, footnote 4] The assumption that ΔE_D(*) and ΔM_B can be shared across all momenta is justified in footnote 4 by the limits p^2 << M^2 and ΔE_D(*) ≈ ΔE_B, but this is not demonstrated quantitatively. At the largest momenta shown, p ≈ 1 GeV, while M_D ≈ 1.87 GeV, so the momentum dependence of sqrt(M_1^2+p^2) - sqrt(M_0^2+p^2) is not negligible. Please provide a numerical check of the momentum dependence of the energy splittings across the fitted momentum range, or allow momentum-dependent splittings and show that the resulting F0 values are unchanged.
  3. [Sec. 3.3, Figs. 4 and 5; Sec. 4] The form factors shown are explicitly bare and unrenormalized, and no continuum or physical-mass extrapolation is performed. The Conclusion's statement that "the results show stable form factors and a good control over systematic effects" therefore overstates the current status: several dominant systematic effects, including current renormalization, discretization, and chiral/continuum extrapolation, are not yet controlled in this analysis. I suggest either restricting the systematic-control claim to the correlation-function fit systematics or presenting a more explicit list of the systematics that have and have not been addressed.
  4. [Sec. 3.2, Figs. 2 and 3] The captions of Figs. 2 and 3 state that the fitted results are insensitive to the fitting window, but this appears to be based on visual inspection; no p-values, reduced chi-squared values, or TIC weights are shown for the three-point fits, in contrast to the two-point analysis in Fig. 1. Reporting these diagnostics for the ratio fits would make the stability claim quantitative and would allow the reader to compare the three-point fits with the two-point fits on equal footing.
minor comments (4)
  1. [Sec. 3.2, Eq. (7)] The smoothing formula R(t,T) = 1/2 R(t,T) + 1/4 R(t,T+1) + 1/2 R(t+1,T+1) has coefficients summing to 1.25; please clarify whether this is an intentional weighted average or a typographical error, since an unbiased smoothing would normally have coefficients summing to one.
  2. [Fig. 5 caption] The caption says "for the ensemble indicated in the legend," but the legend appears to list several ensembles; please specify which ensemble or ensembles are actually shown.
  3. [Sec. 1 and Sec. 4] The abstract uses the future tense ("We will present") while the conclusion uses the past tense ("We presented"); the tense should be made consistent throughout.
  4. [Sec. 2, Eqs. (3) and (4)] The notation for the barred ratios is introduced only in passing; a sentence explicitly defining the bar notation as the value obtained after the ratio fit on the lattice would improve readability.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the form factors are extracted from new 2+1+1 lattice correlator data using a framework inherited from the collaboration's own earlier data-driven analyses, which does not make the present results equivalent to their inputs.

full rationale

The paper does not claim to predict form factors from first principles; it extracts them from two- and three-point correlation functions evaluated on new $N_f=2+1+1$ ensembles. The ratio model in Eq. (8) is a truncated spectral decomposition in which $F_0$ is a fit parameter determined by the same data; presenting the fitted value as a preliminary form factor is a measurement, not a prediction forced by an input. The workflow is taken from Refs. [3,5,12], all prior FNAL-MILC lattice calculations based on independent ensembles and independent data; citing one's own earlier analysis framework is not circularity when the cited results are data-driven and do not assume the present 2+1+1 form factors. The procedure is self-contained in the sense that the final numbers come from the correlator fits on the new ensembles, not from the cited papers' numerical values. The skeptical concern about excited-state contamination in Eq. (8) — one excited state per side with momentum-independent shared splittings — is a legitimate systematic-error or model-validity risk, but it is not circularity: even if the model were biased, the output is not equivalent to the input by construction. No fitted parameter is relabeled as a prediction, no uniqueness theorem is imported from the authors' prior work, no known result is renamed in new coordinates, and the paper makes no external benchmark claim that would be circular. Therefore the circularity score is 0.

Assumptions & free parameters 3 free parameters · 4 assumptions · 0 invented entities

This proceedings paper is an incremental update of the collaboration's B to D(*) form factor analysis. The central claim (preliminary form factor extraction) rests on standard lattice QCD methodology from earlier collaboration papers, plus the numerical fits described here. No new physics entities or fundamental constants are introduced; the fitted quantities are the bare form factors themselves and the nuisance parameters of the correlator fits.

free parameters (3)
  • Energy splittings and amplitudes in three-point fits (Delta E_source, Delta E_sink, A(p), B(p)) = not quoted
    These are fitted to the ratio data in Eq. (8); the extracted ground-state value F0 is the bare form factor. The values are not reported in the paper.
  • Fit analysis hyperparameters (t_min, N_states) = t_min around 0.6 fm for D(*), N_states=3
    Chosen by hand and scanned to identify stability regions; the final results are taken from model averages, so these choices carry systematic uncertainty.
  • Overlap factors Z_a^(n) and energies E_n in two-point fits = not quoted
    Fitted to two-point correlators; they enter the ratios via the normalization of the three-point functions.
assumptions (4)
  • domain assumption The Fermilab interpretation provides a valid effective theory for the b-quark on these lattices.
    Invoked in Sec. 2. If the clover action with Fermilab interpretation fails to reproduce the heavy-quark dynamics, the extracted form factors would not correspond to physical B mesons.
  • domain assumption The spectral decomposition of Eqs. (5) and (6) can be truncated to a small number of states (N_states=3 in two-point, single excited state in three-point).
    Used throughout Sec. 3. The fits assume the truncation is sufficient for the chosen time windows; the paper tests stability across windows but does not prove the absence of further excited states.
  • domain assumption The smoothing procedure of Eq. (7) removes the oscillating staggered-fermion contributions.
    Adopted from Ref. [12] to define the ratios; this is a standard but non-trivial step in staggered analyses.
  • ad hoc to paper The renormalization factors for the currents are approximately one.
    Stated in Sec. 3.3 to justify displaying unrenormalized form factors; this assumption is used only for presentation, not for the final physics claim.

how reviews work

0 comments
Cite this review

Pith. "Pith review of $B\rightarrow D^{(*)}$ decays from $N_f=2+1+1$ highly improved staggered quarks and clover $b$-quark in the Fermilab interpretation." pith.science (2026). https://pith.science/paper/DSBGWOU6

@misc{pith2026250118969,
  author       = {Pith},
  title        = {Pith review of: $B\rightarrow D^(*)$ decays from $N_f=2+1+1$ highly improved staggered quarks and clover $b$-quark in the Fermilab interpretation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/DSBGWOU6}},
  note         = {Machine review of arXiv:2501.18969}
}
abstract

We present an update on the analysis of semileptonic $B\rightarrow D^{(*)}$ decays at non-zero recoil. Our computation employs $2 + 1 + 1$ FNAL-MILC ensembles with highly improved staggered quark (HISQ) action for sea and light valence quarks, while the bottom quark is treated using the clover action in the Fermilab interpretation. Simulations are performed across several lattice spacings, ranging approximately from $\sim 0.15$ fm to $\sim 0.06$ fm, and for various quark masses. We will present an overview of the analysis and show some preliminary results for the form factors.

Figures

Figures reproduced from arXiv: 2501.18969 by the authors.

Figure 1
Figure 1. Fit results for 𝐷-meson correlator at momentum (2, 0, 0) for 𝑎 = 0.09, 𝑚𝑙/𝑚𝑠 = 0.2 for different values of the lower end of the fitting window 𝑡min and different number of excited states considered in the fitting model. Each dot corresponds to a value extracted from a different fit, each 𝑥-coordinate corresponds to a different value of 𝑡min each colour to a different 𝑁states. Upper panel: The 𝑦-axis correspond to th… view at source ↗
Figure 2
Figure 2. Results of the global fit of the relevant ratios for the 𝐵 → 𝐷 process for the same ensemble of the example case in Fig.1. The dashed lines indicate the fitted functional forms projected to every case, while the colored bands represent the value of 𝐹0 with its error. We observe that the fit results are insensitive to the fitting window’s choice. Each color represents a different momentum as in the legend, except for… view at source ↗
Figure 3
Figure 3. [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: Bare, blinded form factors for the 𝐵 → 𝐷 process, defined in Eq. (3). 0.6 0.7 0.8 0.9 preliminary hA1 [not renorm.] a ∼ 0.12 fm, ml/ms ∼ 0.2 a ∼ 0.12 fm, ml/ms ∼ 0.1 a ∼ 0.12 fm, ml/ms ∼[phys] a ∼ 0.09 fm, ml/ms ∼ 0.1 0.9 1.0 1.1 1.2 1.3 1.4 preliminary hV [not renorm.…
Figure 5
Figure 5. Figure 5: Bare, blinded form factors for the 𝐵 → 𝐷 ∗ process, for the ensemble indicated in the legend, defined in Eq. (4) 7 [PITH_FULL_IMAGE:figures/full_fig_p007_5.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

12 extracted references · 5 canonical work pages

  1. [1]

    Flavour Lattice A veraging Group (FLAG)collaboration, FLAG Review 2024, 2411.04268

    arXiv 2024

  2. [2]

    $B\to D^{(\ast)}\ell\nu$ at non-zero recoil

    A. Vaquero,𝐵→𝐷(∗)ℓ𝜈 semileptonic decays at non-zero recoil, PoS LATTICE2022 (2022) 250 [2212.10217]

    work page Pith review arXiv 2022

  3. [3]

    Fermilab Lattice and MILC collaboration,Semileptonic form factors for𝐵→𝐷∗ℓ𝜈 at nonzero recoil from2+ 1-flavor lattice QCD: Fermilab Lattice and MILC Collaborations, Eur. Phys. J. C82(2022) 1141 [2105.14019]

    arXiv 2022

  4. [4]

    Fermilab Lattice and MILC collaboration, Lattice QCD Ensembles with Four Flavors of Highly Improved Staggered Quarks,Phys. Rev. D87(2013) 054505 [1212.4768]

    arXiv 2013

  5. [5]

    Fermilab Lattice and MILC collaboration, B→Dℓ𝜈 form factors at nonzero recoil and |V𝑐𝑏| from 2+1-flavor lattice QCD, Phys. Rev. D92(2015) 034506 [1503.07237]. 8 𝐵→𝐷(∗) decaysfrom 𝑁𝑓 = 2+ 1+ 1highlyimprovedstaggeredquarksandclover 𝑏-quarkintheFermilab interpretation Pietro Butti

    arXiv 2015

  6. [6]

    Richardson,The Heavy Quark Potential and the Upsilon, J/psi Systems,Phys

    J.L. Richardson,The Heavy Quark Potential and the Upsilon, J/psi Systems,Phys. Lett. B82 (1979) 272

    work page 1979

  7. [7]

    Fermilab Lattice, MILC collaboration, B- and D-meson decay constants from three-flavor lattice QCD, Phys. Rev. D85(2012) 114506 [1112.3051]

    arXiv 2012

  8. [8]

    Ledoit and M

    O. Ledoit and M. Wolf,Improved estimation of the covariance matrix of stock returns with an application to portfolio selection, Journal of Empirical Finance10 (2003) 603–621

    work page 2003

Show all 12 references
  1. [9]

    Toussaint and W

    D. Toussaint and W. Freeman,Sample size effects in multivariate fitting of correlated data, 0808.2211

    arXiv

  2. [10]

    Bruno and R

    M. Bruno and R. Sommer,On fits to correlated and auto-correlated data,Comput. Phys. Commun. 285 (2023) 108643 [2209.14188]

    2023 arXiv

  3. [11]

    Frison,Towards fully bayesian analyses in Lattice QCD, 2302.06550

    J. Frison,Towards fully bayesian analyses in Lattice QCD, 2302.06550

    arXiv

  4. [12]

    Fermilab Lattice and MILC Collaborations collaboration,𝐵→𝐷∗ℓ𝜈 form factor at zero recoil from three-flavor lattice qcd: A model independent determination of|𝑉𝑐𝑏|,Phys. Rev. D79(2009) 014506. 9

    2009

Pith tools

Reviewed August 9, 2026 · model on record in the stance chip above.