REVIEW 3 major objections 5 minor 41 references
Predicting the spectrum and decay constants of positive-parity heavy-strange mesons using domain-wall fermions
T0 review · 3 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read This lattice-QCD study computes masses and decay constants for the positive-parity heavy-strange mesons $D^*_{s0}$, $D_{s1}$, $B^*_{s0}$, and $B_{s1}$, with the bottom-strange decay constants appearing for the first time in a lattice…
desk verdict First lattice decay constants for B*_s0 and B_s1, but the central numbers come from a fit the authors themselves suspect is biased; a transparent, well-scoped work-in-progress report rather than an established result. 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 a correlation-matrix analysis of zero-momentum two-point functions built from a four-operator basis: a local quark-antiquark operator $\bar{Q}s$, a derivative operator $\bar{Q}\gamma_i\nabla_i s$, the renormalized current $J_{V_0}$ or $J_{A_i}$ that defines each decay constant, and a two-meson operator of the form $\Phi_K \Phi_H$ (for $0^+$) or $\Phi_K \Phi_{H^*}$ (for $1^+$) that couples to the $D^{(*)}K$ / $B^{(*)}K$ scattering states. The two-meson operator appears only at the source, reusing precomputed light and strange propagators through a sequential heavy-quark propagator; the fitted overlap $A_3$ of the current operator gives the decay constant via $f = A_3\sqrt{2/E}$. The central issue is whether this basis, with or without the two-meson operator, extracts the ground-state energy without bias; the paper's two analysis types are exactly this comparison.
What would settle it
Reanalyze the same correlation functions with two-meson operators at both source and sink, or extract the $DK$ and $BK$ bound states from a finite-volume quantization of the scattering phase shifts. If the resulting pole sits below the no-two-meson fit and above the two-meson-source fit, the two reported analyses bracket the truth; if it agrees with one, that analysis is the unbiased one.
Extended reading notes
Core claim
On its own terms, the paper's claim is that a lattice-QCD calculation with $2+1$ dynamical domain-wall fermions and a tuned anisotropic clover heavy-quark action reproduces the expected below-threshold positive-parity heavy-strange mesons and delivers the first lattice estimates of $f_{B^*_{s0}}$ and $f_{B_{s1}}$. The masses are reported as finite-volume binding energies $\Delta_H = E - E_{H^{(*)}} - E_K$, and the decay constants are obtained from the ground-state overlap of the renormalized vector and axial-vector currents. The paper explicitly labels the results preliminary and says that fits without two-meson operators at the source are more stable yet may overestimate the ground-state energies, while the alternative fits underestimate them; the final answer is therefore bracketed rather than stated.
Load-bearing premise
The load-bearing premise is that the fits without two-meson operators at the source recover the true ground-state energies; the paper itself concedes these fits may overestimate the energies, leaving the final values bracketed between two biased analyses.
Editorial extensions
If this is right
- The bottom-strange decay constants $f_{B^*_{s0}}$ and $f_{B_{s1}}$ would become the first lattice determinations in that sector, replacing QCD sum-rule estimates as the direct, systematically improvable input.
- The observed strong lattice-spacing dependence in the $B^*_{s0}$ and $B_{s1}$ decay constants, if physical and not a fit artifact, will demand a careful continuum extrapolation before these numbers can be used in phenomenology.
- The mass results, once combined with a finite-volume quantization analysis of the $DK$ and $BK$ systems, would convert the bracketed finite-volume energies into infinite-volume bound-state pole positions for comparison with experiment.
- The decay constants would enter the single-particle bound-state contributions in dispersive and unitarity bounds on $b\to s$ and $c\to s$ form factors, sharpening the parametrization of semileptonic decays used in new-physics searches.
Reading between the lines
- A direct extension the authors leave implicit: computing with two-meson operators at both source and sink, or at the sink only, would break the asymmetry that causes their false plateaus and could decide which analysis is biased, using the same correlators they already have.
- The bracketing structure suggests the true ground-state energies, and hence the true decay constants, lie between the two reported fit results; a weighted combination of the two analyses would already give a range, even before the planned chiral-continuum extrapolation.
- If the lattice-spacing dependence of the bottom-strange decay constants survives the removal of the two-meson bias, the anisotropic clover action's improvement would be the prime suspect, and tests on finer lattices or at fixed physical volume would be the natural check.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This proceedings contribution reports a lattice-QCD study of the masses and decay constants of the positive-parity heavy-strange mesons D*_{s0}, D_{s1}, B*_{s0}, and B_{s1}. The simulations use seven RBC/UKQCD domain-wall fermion ensembles with pion masses from 139 to 431 MeV and an anisotropic clover action for charm and bottom quarks. The paper describes two analysis strategies: fits that include two-meson operators only at the source, and fits that omit two-meson operators entirely. The latter are more stable but, as the authors state, may overestimate the ground-state energies. The authors present only the second analysis type, finding below-threshold ground states for all four mesons. They emphasize that the results are preliminary, that the largest systematic uncertainty is the missing two-meson-operator analysis, and that continuum extrapolation and a Luescher analysis are planned. The claimed novelty is the first lattice determination of the B*_{s0} and B_{s1} decay constants.
Significance. If the results survive the planned improvements, they would provide a useful addition to heavy-strange spectroscopy and to the parametrization of semileptonic form factors. The use of multiple ensembles including a near-physical pion mass, exact chiral symmetry for light quarks, and a heavy-quark action tuned to D/B spectroscopy are strengths. The authors are also commendably explicit about the limitations of the present fits. However, the current quantitative claims, especially the first-time values for f_{B*_{s0}} and f_{B_{s1}}, are not yet supported because the extraction omits two-meson operators at the sink and the B-sector decay constants show strong lattice-spacing dependence without a continuum extrapolation. At this stage the paper is an honest status report rather than a finalized calculation.
major comments (3)
- [Sec. 4, Eqs. (4)-(5)] The decay constants are extracted from the fits without two-meson operators via f = A_3 sqrt(2/E), while the paper itself states that excluding these operators 'may ... overestimation of the ground-state energy levels' (Sec. 4). Since A_3 and E come from the same fit, any bias in E propagates directly into f. Because D*_{s0}, D_{s1}, B*_{s0}, and B_{s1} lie close to the DK or BK thresholds, a basis containing only quark-antiquark and local current operators cannot reliably project the physical state. The claim in Sec. 1 that this is 'for the first time' the decay constants of B*_{s0} and B_{s1} are computed is therefore not supported by the presented numbers; either the two-meson-operator analysis at both source and sink must be included, or the claim must be weakened to an exploratory study.
- [Fig. 3 and Sec. 4] The B-sector decay constants in Fig. 3 show a strong dependence on lattice spacing, and the paper provides no continuum extrapolation. The text itself states that the fit results 'may still be unreliable due to the exclusion of the two-meson operators.' Without a controlled extrapolation or a Luescher analysis, the plotted values are not final lattice determinations. The abstract and figure captions should state this explicitly, and the phrase 'first time' should be reserved for a publication that includes the planned analyses (two-meson operators at the sink, Luescher method, chiral-continuum extrapolation).
- [Sec. 4] The paper presents two analysis types that bracket the true result: one underestimates the ground-state energy due to false plateaus (Refs. [38,39]) and one may overestimate it. The paper does not provide a combined estimate or a final central value. This is not only a caveat: the paper's central claim to 'present a lattice-QCD calculation of the masses and decay constants' is not quantitatively fulfilled until a controlled extraction is available. The manuscript should either present such an extraction or be reframed as a status report.
minor comments (5)
- [Sec. 3.2, Eq. (3)] The convention that 'Phi(3)' denotes the current is introduced only later in the text; the operator numbering in Eqs. (1)-(2) should be made explicit at first use.
- [Table 2] There is a formatting issue in the pion mass column ('0 .13917(35)'), and the caption should be grammatically corrected.
- [Figs. 2 and 3] The figure captions should state explicitly which bands are experimental PDG values, which are previous lattice results, and whether the plotted points include the systematic uncertainty from the t_min variation.
- [Sec. 3.2] The procedure for combining statistical uncertainties with the t_min shift is described only in Sec. 4; it would be clearer to mention it in the fitting section as well.
- [Throughout] The paper does not provide a table of numerical results, which limits usability; even a preliminary summary table would improve the proceedings contribution.
Circularity Check
No circularity: the spectrum and decay constants are extracted from standard multi-exponential fits to lattice correlation functions, with input parameters tuned to different (negative-parity) states and external benchmarks.
full rationale
The derivation chain is a standard lattice-QCD spectroscopy analysis. Energies and current amplitudes are obtained by fitting correlation functions C_lm(t) built from local quark-antiquark and (at the source only) two-meson interpolating operators; the ground-state decay constant is then computed from the fitted amplitude via f = A_3 sqrt(2/E). This is a normalization conversion, not a fitted input recycled as a prediction. The heavy-quark action parameters are taken from prior work (Refs. [26,29]), but they were tuned by matching the negative-parity D_s^{(*)} and B_s^{(*)} dispersion relations and hyperfine splittings, not the positive-parity states computed here; the target results are therefore not imposed by the input. The current renormalization factors come from nonperturbative determinations and one-loop matching, again independent of the final decay-constant values. The paper explicitly flags its main systematic caveat in Sec. 4: 'concerns remain that without the two-meson operators there may now be an overestimation of the ground-state energy levels,' and it notes the strong lattice-spacing dependence of the bottom-strange decay constants. That is an acknowledged operator-basis uncertainty, not a circular step. The self-citations (e.g., Ref. [26]) supply infrastructure and parameter values, but the central claim is checked against external PDG masses and prior lattice results, so the derivation is self-contained. No circularity is present.
Assumptions & free parameters
free parameters (2)
- Heavy quark action parameters (a m_Q, nu, c_B=c_E) =
from Ref. [26]
- Renormalization factors Z and O(a) improvement coefficients =
nonperturbative Z; 1-loop rho and improvement
assumptions (4)
- domain assumption The RBC/UKQCD ensembles with 2+1 dynamical domain-wall fermions correctly sample the QCD path integral.
- domain assumption The anisotropic clover action tuned per Ref. [26] correctly describes charm and bottom quarks at the lattice spacings used.
- ad hoc to paper Ground-state energies can be extracted from fits without two-meson operators.
- domain assumption Correlators with two-hadron operators at one end only produce negative bias (Refs. [38,39]), motivating the choice of analysis.
Cite this review
Pith. "Pith review of Predicting the spectrum and decay constants of positive-parity heavy-strange mesons using domain-wall fermions." pith.science (2026). https://pith.science/paper/C2UUHNYU
@misc{pith2026250117846,
author = {Pith},
title = {Pith review of: Predicting the spectrum and decay constants of positive-parity heavy-strange mesons using domain-wall fermions},
year = {2026},
howpublished = {\url{https://pith.science/paper/C2UUHNYU}},
note = {Machine review of arXiv:2501.17846}
}
abstract
We present a lattice-QCD calculation of the masses and decay constants of the positive-parity heavy-strange mesons $D^*_{s0}$, $D_{s1}$, $B^*_{s0}$, and $B_{s1}$. The calculations are performed with domain-wall fermions for the light and strange quarks and an anisotropic clover action for the charm and bottom quarks. We use seven different RBC/UKQCD ensembles with pion masses ranging from a near-physical 139 MeV up to 431 MeV. We consider two different analysis types, with or without two-meson operators at the source. We observe the expected below-threshold ground states. The fits without the two-meson operators appear to be more stable, but may overestimate the ground-state energies, while preliminary fits with two-meson operators at the source only appear to underestimate the ground-state energies.
Figures
Reference graph
Works this paper leans on
-
[1]
T. Barnes, F.E. Close and H.J. Lipkin,Implications of a DK molecule at 2.32-GeV,Phys. Rev. D68(2003) 054006 [hep-ph/0305025]
arXiv 2003
-
[2]
BaBarcollaboration, Observation of a narrow meson decaying to𝐷+ 𝑠𝜋0 at a mass of 2.32-GeV/c2,Phys. Rev. Lett.90 (2003) 242001 [hep-ex/0304021]
arXiv 2003
-
[3]
S. Godfrey and R. Kokoski,The Properties of p Wave Mesons with One Heavy Quark,Phys. Rev. D43(1991) 1679. 8 Spectrum and decay constants of positive-parity heavy-strange mesons Forrest Guyton
work page 1991
-
[4]
F.-K. Guo, C. Hanhart, U.-G. Meißner, Q. Wang, Q. Zhao and B.-S. Zou,Hadronic molecules,Rev. Mod. Phys.90(2018) 015004 [1705.00141]. [Erratum:Rev. Mod. Phys.94 (2022) 029901]
arXiv 2018
-
[5]
T. Barnes,The Status of molecules, in29th Rencontres de Moriond: QCD and High-energy Hadronic Interactions, hep-ph/9406215
-
[6]
Weinberg,Evidence That the Deuteron Is Not an Elementary Particle,Phys
S. Weinberg,Evidence That the Deuteron Is Not an Elementary Particle,Phys. Rev.137 (1965) B672
work page 1965
-
[7]
UKQCDcollaboration, Comprehensive spectroscopy of𝐷,𝐷𝑠 and charmonium systems from relativistic lattice QCD,Nucl. Phys. B Proc. Suppl.53(1997) 398
work page 1997
-
[8]
D. Mohler, C.B. Lang, L. Leskovec, S. Prelovsek and R.M. Woloshyn,𝐷∗ 𝑠0(2317) Meson and𝐷-Meson-Kaon Scattering from Lattice QCD, Phys. Rev. Lett.111 (2013) 222001 [1308.3175]
arXiv 2013
Show all 41 references
-
[9]
G.S. Bali, S. Collins, A. Cox and A. Schäfer,Masses and decay constants of the𝐷∗ 𝑠0(2317) and𝐷𝑠1(2460) from𝑁𝑓 = 2 lattice QCD close to the physical point, Phys. Rev. D96 (2017) 074501 [1706.01247]
2017 arXiv
-
[10]
Hadron Spectrumcollaboration,𝐷𝐾 𝐼= 0,𝐷𝐾 𝐼= 0, 1scattering and the𝐷∗ 𝑠0(2317) from lattice QCD, JHEP 02(2021) 100 [2008.06432]
2021 arXiv
-
[11]
Alexandrou, J
C. Alexandrou, J. Berlin, J. Finkenrath, T. Leontiou and M. Wagner,Tetraquark interpolating fields in a lattice QCD investigation of the𝐷∗ 𝑠0(2317) meson, Phys. Rev. D101 (2020) 034502 [1911.08435]
2020 arXiv
-
[12]
Burch, C
T. Burch, C. Hagen, C.B. Lang, M. Limmer and A. Schafer,Excitations of single-beauty hadrons,Phys. Rev. D79 (2009) 014504 [0809.1103]
2009 arXiv
-
[13]
C.B. Lang, D. Mohler, S. Prelovsek and R.M. Woloshyn,Predicting positive parity𝐵𝑠 mesons from lattice QCD,Phys. Lett. B750 (2015) 17 [1501.01646]
2015 arXiv
-
[14]
Hudspith and D
R.J. Hudspith and D. Mohler,Exotic tetraquark states with two¯𝑏 quarks and𝐽𝑃 = 0+ and 1+ 𝐵𝑠 states in a nonperturbatively tuned lattice NRQCD setup, Phys. Rev. D107 (2023) 114510 [2303.17295]
2023 arXiv
-
[15]
Capdevila, A
B. Capdevila, A. Crivellin and J. Matias,Review of semileptonic𝐵 anomalies, Eur. Phys. J. ST 1(2023) 20 [2309.01311]
2023 arXiv
-
[16]
Blake, S
T. Blake, S. Meinel, M. Rahimi and D. van Dyk,Dispersive bounds for local form factors in Λ𝑏→Λ transitions, Phys. Rev. D108 (2023) 094509 [2205.06041]
2023 arXiv
-
[17]
Flavour Lattice A veraging Group (FLAG)collaboration,FLAG Review 2024, 2411.04268
2024 arXiv
-
[18]
Narison,Decay Constants of Heavy-Light Mesons from QCD, Nucl
S. Narison,Decay Constants of Heavy-Light Mesons from QCD, Nucl. Part. Phys. Proc. 270-272(2016) 143 [1511.05903]. 9 Spectrum and decay constants of positive-parity heavy-strange mesons Forrest Guyton
2016 arXiv
-
[19]
Q. Li, T. Wang, Y. Jiang, G.-L. Wang and C.-H. Chang,Mixing angle and decay constants of 𝐽𝑃 = 1+ heavy-light mesons, Phys. Rev. D100 (2019) 076020 [1802.06351]
2019 arXiv
-
[20]
Particle Data Groupcollaboration, Review of particle physics, Phys. Rev. D110 (2024) 030001
2024
-
[21]
RBC, UKQCD collaboration, Domain wall QCD with physical quark masses, Phys. Rev. D 93 (2016) 074505 [1411.7017]
2016 arXiv
-
[22]
RBC/UKQCDcollaboration,SU(3)-breaking ratios for𝐷(𝑠) and𝐵(𝑠) mesons, 1812.08791
-
[23]
Kaplan,A Method forsimulating chiral fermions onthe lattice,Phys
D.B. Kaplan,A Method forsimulating chiral fermions onthe lattice,Phys. Lett. B288(1992) 342 [hep-lat/9206013]
1992 arXiv
-
[24]
Brower, H
R.C. Brower, H. Neff and K. Orginos,Mobius fermions, Nucl. Phys. B Proc. Suppl.153 (2006) 191 [hep-lat/0511031]
2006 arXiv
-
[25]
Shamir,Chiral fermions from lattice boundaries,Nucl
Y. Shamir,Chiral fermions from lattice boundaries,Nucl. Phys. B406 (1993) 90 [hep-lat/9303005]
1993 arXiv
-
[26]
Meinel,Status of next-generationΛ𝑏→ 𝑝,Λ,Λ𝑐 form-factor calculations,PoS LATTICE2023(2024) 275 [2309.01821]
S. Meinel,Status of next-generationΛ𝑏→ 𝑝,Λ,Λ𝑐 form-factor calculations,PoS LATTICE2023(2024) 275 [2309.01821]
2024 arXiv
-
[27]
Shintani, R
E. Shintani, R. Arthur, T. Blum, T. Izubuchi, C. Jung and C. Lehner,Covariant approximation averaging, Phys. Rev. D91(2015) 114511 [1402.0244]
2015 arXiv
-
[28]
El-Khadra, A.S
A.X. El-Khadra, A.S. Kronfeld and P.B. Mackenzie,Massive fermions in lattice gauge theory, Phys. Rev. D55(1997) 3933 [hep-lat/9604004]
1997 arXiv
-
[29]
RBC, UKQCD collaboration,Nonperturbative tuning of an improved relativistic heavy-quark action with application to bottom spectroscopy, Phys. Rev. D86 (2012) 116003 [1206.2554]
2012 arXiv
-
[30]
Hashimoto, A.X
S. Hashimoto, A.X. El-Khadra, A.S. Kronfeld, P.B. Mackenzie, S.M. Ryan and J.N. Simone, Lattice QCD calculation of¯𝐵→𝐷ℓ ¯𝜈 decay form-factors at zero recoil,Phys. Rev.D61 (1999) 014502 [hep-ph/9906376]
1999 arXiv
-
[31]
El-Khadra, A.S
A.X. El-Khadra, A.S. Kronfeld, P.B. Mackenzie, S.M. Ryan and J.N. Simone,The Semileptonic decays𝐵→𝜋ℓ𝜈 and𝐷→𝜋ℓ𝜈 from lattice QCD, Phys. Rev.D64 (2001) 014502 [hep-ph/0101023]
2001 arXiv
-
[32]
Boyle, L
P.A. Boyle, L. Del Debbio, A. Jüttner, A. Khamseh, F. Sanfilippo and J.T. Tsang,The decay constants fD and fDs in the continuum limit ofNf = 2+ 1 domain wall lattice QCD,JHEP 12 (2017) 008 [1701.02644]
2017 arXiv
-
[33]
Meinel,Λ𝑐→Λ𝑙+𝜈𝑙 form factors and decay rates from lattice QCD with physical quark masses,Phys
S. Meinel,Λ𝑐→Λ𝑙+𝜈𝑙 form factors and decay rates from lattice QCD with physical quark masses,Phys. Rev. Lett.118 (2017) 082001 [1611.09696]. 10 Spectrum and decay constants of positive-parity heavy-strange mesons Forrest Guyton
2017 arXiv
-
[34]
Meinel and G
S. Meinel and G. Rendon,Λ𝑏→Λ∗ 𝑐(2595, 2625)ℓ− ¯𝜈 form factors from lattice QCD, Phys. Rev. D103 (2021) 094516 [2103.08775]
2021 arXiv
-
[35]
Detmold, C
W. Detmold, C. Lehner and S. Meinel,Λ𝑏→ 𝑝ℓ− ¯𝜈ℓ andΛ𝑏→Λ𝑐ℓ− ¯𝜈ℓ form factors from lattice QCD with relativistic heavy quarks,Phys. Rev. D92(2015) 034503 [1503.01421]
2015 arXiv
-
[36]
Lehner,Automated lattice perturbation theory and relativistic heavy quarks in the Columbia formulation,PoS LATTICE2012(2012) 126 [1211.4013]
C. Lehner,Automated lattice perturbation theory and relativistic heavy quarks in the Columbia formulation,PoS LATTICE2012(2012) 126 [1211.4013]
2012 arXiv
-
[37]
Meinel and G
S. Meinel and G. Rendon,Λ𝑏→Λ∗(1520)ℓ+ℓ− form factors from lattice QCD,Phys. Rev. D 103 (2021) 074505 [2009.09313]
2021 arXiv
-
[38]
B. Hörz et al.,Two-nucleon S-wave interactions at the𝑆𝑈(3) flavor-symmetric point with 𝑚𝑢𝑑≃𝑚phys 𝑠 : A first lattice QCD calculation with the stochastic Laplacian Heaviside method,Phys. Rev. C103 (2021) 014003 [2009.11825]
2021 arXiv
-
[39]
Iritani et al.,Mirage in Temporal Correlation functions for Baryon-Baryon Interactions in Lattice QCD, JHEP 10(2016) 101 [1607.06371]
T. Iritani et al.,Mirage in Temporal Correlation functions for Baryon-Baryon Interactions in Lattice QCD, JHEP 10(2016) 101 [1607.06371]
2016 arXiv
-
[40]
Luscher,Two particle states on a torus and their relation to the scattering matrix,Nucl
M. Luscher,Two particle states on a torus and their relation to the scattering matrix,Nucl. Phys. B354(1991) 531
1991
-
[41]
Boerner, S
T.J. Boerner, S. Deems, T.R. Furlani, S.L. Knuth and J. Towns,ACCESS: Advancing Innovation: NSF’s Advanced Cyberinfrastructure Coordination Ecosystem: Services & Support, inPractice and Experience in Advanced Research Computing 2023: Computing for the Common Good, PEARC ’23, (...
2023
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.