pith. sign in

arxiv: 2602.07829 · v2 · pith:NSLMV44Mnew · submitted 2026-02-08 · ✦ hep-lat · hep-ph

Radiative decay of heavy-light mesons from lattice QCD

Wen-Zheng Hou , Nan Wang , Long-Cheng Gui , Jun Hua , Jian Liang , Jun Shi , Yu Meng This is my paper

Pith reviewed 2026-05-21 14:33 UTC · model grok-4.3

classification ✦ hep-lat hep-ph
keywords lattice QCDcharmed mesonsradiative decaysdecay widthsheavy-light mesonscoupling constantselectromagnetic transitionsCLQCD ensembles
0
0 comments X

The pith

Lattice QCD yields precise first-principles values for the radiative decay couplings and widths of charmed mesons.

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

The paper computes the electromagnetic transition couplings for D, D*, and D_s mesons directly from lattice simulations rather than from models. It employs 2+1-flavor clover fermion ensembles that include a physical-pion-mass point and a fine lattice spacing near 0.05 fm. Systematic effects from matrix-element fits, momentum extrapolations, and chiral-continuum limits are folded into the final uncertainties. A sympathetic reader cares because these results supply controlled, non-perturbative benchmarks that can be compared with experiment and used to test descriptions of heavy-light meson structure.

Core claim

Using 2+1-flavor clover gauge ensembles, the authors extract the coupling constants g_{D^{*+} D^+ gamma} = -0.204(22) GeV^{-1}, g_{D^{*0} D^0 gamma} = 1.73(37) GeV^{-1}, and g_{D_s^{*+} D_s^+ gamma} = -0.120(14) GeV^{-1}. These values produce the decay widths Gamma_{D^{*+} -> D^+ gamma} = 0.253(55) keV, Gamma_{D^{*0} -> D^0 gamma} = 18.2(7.8) keV, and Gamma_{D_s^{*+} -> D_s^+ gamma} = 0.094(22) keV, with all quoted uncertainties incorporating the effects of fit forms, momentum-transfer extrapolations, and chiral-continuum extrapolations.

What carries the argument

Three-point lattice correlation functions that isolate the electromagnetic transition matrix elements, followed by simultaneous extrapolations in quark mass, lattice spacing, and momentum transfer.

If this is right

  • The computed widths supply direct first-principles inputs for phenomenological analyses of charm radiative transitions.
  • The method can be applied to other heavy-light systems once similar ensembles become available.
  • The difference between charged and neutral modes quantifies isospin-breaking effects in the electromagnetic decays.
  • The results test the consistency of lattice QCD with the expected hierarchy of decay rates for vector to pseudoscalar transitions.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Confirmation by experiment would strengthen in lattice techniques for computing electromagnetic matrix elements involving heavy quarks.
  • The same framework could be extended to bottom mesons to predict their radiative widths before experimental data arrive.
  • Discrepancies with older quark-model estimates would highlight the importance of sea-quark and discretization effects that only lattice calculations capture.
  • The quoted uncertainties already allow these widths to serve as constraints in global fits of heavy-meson parameters.

Load-bearing premise

The chiral, continuum, and momentum-transfer extrapolations together with the chosen fit forms for the matrix elements introduce no bias larger than the total uncertainties that are reported.

What would settle it

A high-precision experimental measurement of the D^{*0} -> D^0 gamma width that lies well outside the interval 10.4-25.9 keV would falsify the central result.

Figures

Figures reproduced from arXiv: 2602.07829 by Jian Liang, Jun Hua, Jun Shi, Long-Cheng Gui, Nan Wang, Wen-Zheng Hou, Yu Meng.

Figure 1
Figure 1. Figure 1: shows the effective energies and the corresponding fit￾ting results of the D meson on the F32P21 ensemble. We start the fits at time slice t/a = 12, and the correlations between time slices are taken into consideration. The fitting results are shown as bands, with uncertainties that are comparable to those of the chosen starting data points. The subsequent data points almost all lie within one sigma of the… view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Dispersion relation results for the [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Results for [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Comparison of the two di [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: The upper panel illustrates the distribution of [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: The momentum transfer extrapolation of [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: The first three rows display the combined continuum and chiral extrapolation for [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Momentum transfer extrapolations for [PITH_FULL_IMAGE:figures/full_fig_p011_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: The same as Fig [PITH_FULL_IMAGE:figures/full_fig_p012_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10: The same as Fig [PITH_FULL_IMAGE:figures/full_fig_p012_10.png] view at source ↗
read the original abstract

We present the first systematic study of the radiative decays of charmed mesons using $2+1$-flavor clover fermion gauge ensembles generated by the CLQCD collaboration. One of the ensembles is at the physical pion mass, and one has a fine lattice spacing $a\sim 0.05 ~\text{fm}$. We determine the coupling constants to be $g_{D^{\ast+} D^+ \gamma} = -0.204(22)$ GeV$^{-1}$, $g_{D^{\ast0} D^0 \gamma} = 1.73(37)$ GeV$^{-1}$, and $g_{D_s^{\ast+} D_s^+ \gamma} =-0.120(14)$ GeV$^{-1}$, respectively. Compared with previous studies, our results demonstrate significant improvements in precision. In particular, we carefully estimate the systematic uncertainty arising from matrix element fits, momentum transfer extrapolations, and chiral and continuum limit extrapolations, which are included in the reported total uncertainties. These couplings yield the following predictions of decay widths: $\Gamma_{D^{\ast+} \rightarrow D^+ \gamma} = 0.253(55)$ keV, $\Gamma_{D^{\ast0} \rightarrow D^0 \gamma} = 18.2(7.8)$ keV, and $\Gamma_{D_s^{\ast+}\rightarrow D_s^+ \gamma} = 0.094(22)$ keV. This work establishes first-principles results of the charmed meson radiative transitions and provides inputs for understanding the structure and properties of heavy-light mesons.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript presents the first systematic lattice QCD study of radiative decays of charmed heavy-light mesons (D*, D_s*) on 2+1-flavor clover fermion ensembles generated by CLQCD. One ensemble is at the physical pion mass and one has a fine spacing a≈0.05 fm. The authors extract the couplings g_{D^{*+} D^+ gamma} = -0.204(22) GeV^{-1}, g_{D^{*0} D^0 gamma} = 1.73(37) GeV^{-1}, and g_{D_s^{*+} D_s^+ gamma} = -0.120(14) GeV^{-1}, then derive the decay widths Gamma_{D^{*+} -> D^+ gamma} = 0.253(55) keV, Gamma_{D^{*0} -> D^0 gamma} = 18.2(7.8) keV, and Gamma_{D_s^{*+} -> D_s^+ gamma} = 0.094(22) keV, with total uncertainties stated to incorporate matrix-element fit, q^2 extrapolation, and chiral-continuum systematics.

Significance. If the central values and uncertainties prove robust, the work supplies valuable first-principles inputs for the radiative transitions of charmed mesons that can be used in phenomenological analyses of heavy-light meson structure. The presence of a physical-pion-mass ensemble and a fine a≈0.05 fm ensemble, together with the explicit attempt to fold multiple sources of systematic uncertainty into the quoted errors, are clear strengths that advance beyond earlier lattice studies.

major comments (1)
  1. [Sections describing the fit forms and extrapolation procedure (likely §4–5)] The multi-stage extrapolation chain (finite-q^2 matrix-element fits, q^2→0 extrapolation, followed by simultaneous chiral/continuum extrapolation) is load-bearing for the final central values and error bars. With only one physical-pion-mass ensemble and a modest number of total ensembles, the stability of the chosen functional forms (polynomial orders in m_π² and q², possible omission of heavy-meson chiral logarithms) must be demonstrated explicitly; otherwise the claim that all systematic uncertainties are absorbed into the reported totals remains under-supported.
minor comments (2)
  1. Clarify the precise number of gauge ensembles, the range of pion masses, and the values of the lattice spacings in a single table for quick reference.
  2. Add a brief discussion of how the quoted total uncertainties were constructed (e.g., quadrature sum, or more conservative envelope) to aid reproducibility.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comment regarding the extrapolation procedure. We appreciate the positive assessment of the overall approach and the recognition of the strengths provided by the physical-pion-mass ensemble and the fine-spacing ensemble. We address the major comment below.

read point-by-point responses

  1. Referee: [Sections describing the fit forms and extrapolation procedure (likely §4–5)] The multi-stage extrapolation chain (finite-q^2 matrix-element fits, q^2→0 extrapolation, followed by simultaneous chiral/continuum extrapolation) is load-bearing for the final central values and error bars. With only one physical-pion-mass ensemble and a modest number of total ensembles, the stability of the chosen functional forms (polynomial orders in m_π² and q², possible omission of heavy-meson chiral logarithms) must be demonstrated explicitly; otherwise the claim that all systematic uncertainties are absorbed into the reported totals remains under-supported.

    Authors: We agree that explicit demonstration of stability is important given the modest number of ensembles. In the revised manuscript we have added an appendix that varies the polynomial orders in m_π² and q² and tests the effect of including heavy-meson chiral logarithms in the fit forms. The central values and uncertainties remain stable under these variations, with any differences folded into the systematic error budget. This additional material supports the claim that the quoted total uncertainties encompass the relevant extrapolation systematics. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results follow from direct lattice matrix-element evaluation on external ensembles

full rationale

The derivation computes the radiative couplings g from three-point correlation functions on CLQCD 2+1-flavor clover ensembles (including a physical-pion-mass point and a fine a≈0.05 fm ensemble), followed by standard q², chiral, and continuum extrapolations whose functional forms are chosen and whose uncertainties are folded into the quoted total errors. The decay widths are obtained from these g values via the usual tree-level kinematic formula relating width to coupling and phase space; this step is a fixed conversion, not a fit to the target observable. No self-definitional loop, no fitted parameter renamed as a prediction, and no load-bearing self-citation chain appears in the abstract or described procedure. The calculation is therefore self-contained against external gauge configurations and does not reduce to its own inputs by construction.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The calculation rests on standard lattice QCD with clover fermions and on the external CLQCD gauge ensembles; no new particles or forces are postulated and the only adjustable elements are the usual lattice parameters already fixed by the collaboration.

free parameters (1)
  • lattice spacing and quark-mass parameters of the CLQCD ensembles
    These are inputs taken from the generating collaboration rather than fitted inside the present analysis.
axioms (1)
  • domain assumption Clover fermion action with 2+1 flavors reproduces the correct low-energy QCD dynamics on the lattice
    Invoked by the choice of the CLQCD gauge ensembles used for all correlation functions.

pith-pipeline@v0.9.0 · 5836 in / 1390 out tokens · 74411 ms · 2026-05-21T14:33:38.512708+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

24 extracted references · 24 canonical work pages · 13 internal anchors

  1. [1]

    J. E. Barteltet al.(CLEO), Phys. Rev. Lett.80, 3919 (1998), arXiv:hep-ex/9711011

    work page internal anchor Pith review Pith/arXiv arXiv 1998

  2. [2]

    Precision measurement of the $D^{*0}$ decay branching fractions

    M. Ablikimet al.(BESIII), Phys. Rev. D91, 031101 (2015), arXiv:1412.4566 [hep-ex]

    work page internal anchor Pith review Pith/arXiv arXiv 2015

  3. [3]

    Ablikimet al.(BESIII), Phys

    M. Ablikimet al.(BESIII), Phys. Rev. D107, 032011 (2023), arXiv:2212.13361 [hep-ex]

    work page arXiv 2023

  4. [4]

    Li, C.-D

    H.-D. Li, C.-D. L ¨u, C. Wang, Y .-M. Wang, and Y .-B. Wei, JHEP 04, 023, arXiv:2002.03825 [hep-ph]

    work page arXiv 2002

  5. [5]

    Pullin and R

    B. Pullin and R. Zwicky, JHEP09, 023, arXiv:2106.13617 [hep-ph]

    work page arXiv

  6. [6]

    J. F. Amundson, C. G. Boyd, E. E. Jenkins, M. E. Luke, A. V . Manohar, J. L. Rosner, M. J. Savage, and M. B. Wise, Phys. Lett. B296, 415 (1992), arXiv:hep-ph/9209241

    work page internal anchor Pith review Pith/arXiv arXiv 1992

  7. [7]

    C.-T. Tran, M. A. Ivanov, P. Santorelli, and Q.-C. V o, Chin. Phys. C48, 023103 (2024), arXiv:2311.15248 [hep-ph]

    work page arXiv 2024

  8. [8]

    Strong and radiative decays of heavy mesons in a covariant model

    C.-Y . Cheung and C.-W. Hwang, JHEP04, 177, arXiv:1401.3917 [hep-ph]

    work page internal anchor Pith review Pith/arXiv arXiv

  9. [9]

    A. H. Orsland and H. Hogaasen, Eur. Phys. J. C9, 503 (1999), arXiv:hep-ph/9812347

    work page internal anchor Pith review Pith/arXiv arXiv 1999

  10. [10]

    Jaus, Phys

    W. Jaus, Phys. Rev. D53, 1349 (1996), [Erratum: Phys.Rev.D 54, 5904 (1996)]

    work page 1996

  11. [11]

    J. L. Goity and W. Roberts, Phys. Rev. D64, 094007 (2001), arXiv:hep-ph/0012314

    work page internal anchor Pith review Pith/arXiv arXiv 2001

  12. [12]

    T. M. Aliev, E. Iltan, and N. K. Pak, Phys. Lett. B334, 169 (1994)

    work page 1994

  13. [13]

    D*-->Dpi and D*-->Dgamma decays: Axial coupling and Magnetic moment of D* meson

    D. Becirevic and B. Haas, Eur. Phys. J. C71, 1734 (2011), arXiv:0903.2407 [hep-lat]

    work page internal anchor Pith review Pith/arXiv arXiv 2011

  14. [14]

    G. C. Donald, C. T. H. Davies, J. Koponen, and G. P. Lepage, Phys. Rev. Lett.112, 212002 (2014), arXiv:1312.5264 [hep- lat]

    work page internal anchor Pith review Pith/arXiv arXiv 2014

  15. [15]

    Meng, J.-L

    Y . Meng, J.-L. Dang, C. Liu, Z. Liu, T. Shen, H. Yan, and K.- L. Zhang, Phys. Rev. D109, 074511 (2024), arXiv:2401.13475 [hep-lat]

    work page arXiv 2024

  16. [16]

    Frezzotti, G

    R. Frezzotti, N. Tantalo, G. Gagliardi, F. Sanfilippo, S. Simula, V . Lubicz, F. Mazzetti, G. Martinelli, and C. T. Sachrajda, Phys. Rev. D108, 074505 (2023), arXiv:2306.05904 [hep-lat]

    work page arXiv 2023

  17. [17]

    Huet al.(CLQCD), Phys

    Z.-C. Huet al.(CLQCD), Phys. Rev. D109, 054507 (2024), arXiv:2310.00814 [hep-lat]

    work page arXiv 2024

  18. [18]

    J. P. Leeset al.(BaBar), Phys. Rev. Lett.111, 111801 (2013), arXiv:1304.5657 [hep-ex]

    work page internal anchor Pith review Pith/arXiv arXiv 2013

  19. [19]

    Hadron Spectrum with Wilson fermions

    T. Bhattacharya, R. Gupta, G. Kilcup, and S. R. Sharpe, Phys. Rev. D53, 6486 (1996), arXiv:hep-lat/9512021

    work page internal anchor Pith review Pith/arXiv arXiv 1996

  20. [20]

    Lattice QCD estimate of the $\eta_{c}(2S)\to J/\psi\gamma$ decay rate

    D. Be ˇcirevi´c, M. Kruse, and F. Sanfilippo, JHEP05, 014, arXiv:1411.6426 [hep-lat]

    work page internal anchor Pith review Pith/arXiv arXiv

  21. [21]

    R. J. Hill and G. Paz, Phys. Rev. D82, 113005 (2010), arXiv:1008.4619 [hep-ph]

    work page internal anchor Pith review Pith/arXiv arXiv 2010

  22. [22]

    Sint and P

    S. Sint and P. Weisz, Nucl. Phys. B502, 251 (1997), arXiv:hep- lat/9704001

    work page arXiv 1997

  23. [23]

    Non-perturbative renormalization and O$(a)$-improvement of the non-singlet vector current with $N_{\mathrm{f}}=2+1$ Wilson fermions and tree-level Symanzik improved gauge action

    A. Gerardin, T. Harris, and H. B. Meyer, Phys. Rev. D99, 014519 (2019), arXiv:1811.08209 [hep-lat]

    work page internal anchor Pith review Pith/arXiv arXiv 2019

  24. [24]

    Value without systematic

    J. Heitger and F. Joswig (ALPHA), Eur. Phys. J. C81, 254 (2021), arXiv:2010.09539 [hep-lat]. Appendix A: Additional figures and tables 10 Ensemble F32P21 C32P23 H48P32 F32P30 C24P29 C48P14 gD∗+ D+γ Value without systematic (GeV−1) -0.291(12) -0.397(14) -0.1761(71) -0.2499(86) -0.361(10) -0.443(17) Momentum transfer extrapolation (GeV−1) 0.015 0.00002 0.01...

    work page arXiv 2021