Pith. sign in

REVIEW 3 major objections 5 minor 18 references

Real photon emissions in leptonic decays

T0 review · 3 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

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

Pith's one-line read A lattice QCD calculation extracts the structure-dependent form factors FA and FV that govern radiative leptonic decays of pseudoscalar mesons, cleanly separating them from the point-like photon contribution.

desk verdict Useful progress report: first lattice extraction of D/Ds structure-dependent form factors, but the 'good precision' claim needs error bars and caveats on systematics. read the letter →

arxiv 1908.10160 v1 pith:SXXUH5KT submitted 2019-08-27 hep-lat hep-ph

G.M. de Divitiis , A. Desiderio , M. Di Carlo , R. Frezzotti , M. Garofalo , D. Giusti , M. Hansen , V. Lubicz
show 6 more authors
F. Mazzetti G. Martinelli C.T. Sachrajda F. Sanfilippo S. Simula N. Tantalo
This is my paper · ORCID
classification hep-lathep-ph
keywords radiativeleptonicdecayslatticeQCDstructure-dependentformfactorspseudoscalarmesonselectromagneticcorrectionsCKMmatrixtwistedboundaryconditionsinfrared-safeobservables
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

The paper reports a first-principles lattice QCD calculation of the two structure-dependent form factors $F_A$ and $F_V$ that enter the amplitude for radiative leptonic decays $P\to \ell\bar\nu_\ell\gamma$ of pseudoscalar mesons. It shows that, even with moderate statistics, Euclidean correlation functions can be arranged into plateau-forming ratios from which $F_A$ and $F_V$ are extracted with good precision for kaons and $D$ mesons, and that their momentum dependence can be mapped over the kinematic range. The method separates unambiguously the point-like, infrared-divergent part of the amplitude from the infrared-safe structure-dependent part, so the latter no longer has to be estimated by chiral perturbation theory or model assumptions. If the result holds, it opens the way to $O(\alpha_{\rm em})$ predictions of leptonic decay rates from the pion to the $B$ meson, improving tests of the Standard Model and the determination of the Cabibbo-Kobayashi-Maskawa matrix elements.

What carries the argument

The machinery is the ratio estimator: $R_A(t)$ projects out $F_A(x_\gamma)+2f_P/(m_P x_\gamma)$ and $R_V(t)$ projects out $F_V(x_\gamma)$ from a Euclidean time-ordered product of the hadronic weak current and the conserved electromagnetic current. Twisted spatial boundary conditions on the quark fields allow the two propagators attached to the electromagnetic current to carry different phases, so arbitrary meson and photon spatial momenta can be selected on a fixed lattice; the point-like term is later subtracted using $f_P$ from a standard two-point function. The argument carries because the amplitude decomposition contains only the four scalar form factors $H_1$, $H_2$, $F_A$, and $F_V$, and the point-like term saturates the Ward identity, making the separation between infrared-divergent and infrared-safe contributions unambiguous.

What would settle it

One concrete check is to compute the same two-current correlator on the same ensembles with the disconnected sea-quark contribution to the electromagnetic current included, and compare the resulting $F_A(x_\gamma)$ and $F_V(x_\gamma)$ with the electro-quenched values; a shift larger than the statistical errors would falsify the present extraction.

Watch

Extended reading notes

Core claim

The central claim is that the form factors needed for $P\to \ell\bar\nu_\ell\gamma$ can be extracted directly from lattice QCD. The radiative amplitude is decomposed into a point-like term fixed by the decay constant $f_P$, which saturates the Ward identity and produces the infrared divergence, and two structure-dependent form factors $F_A(x_\gamma)$ and $F_V(x_\gamma)$. The authors compute the two-current Euclidean correlator and build ratios $R_A(t)$ and $R_V(t)$ that exhibit plateaus at intermediate times; subtracting the known point-like piece using the independently measured $f_P$ isolates $F_A(x_\gamma)$. With a preliminary ensemble at one lattice spacing, they obtain non-zero signals for both form factors for kaons and $D_s$ mesons, cover the full physical $x_\gamma$ range for kaons, and compare the kaon results with chiral perturbation theory predictions. The conclusion is that, with moderate statistics, both light and heavy meson structure-dependent form factors can be extracted precisely enough to study their momentum dependence, and that a fully non-perturbative $O(a)$-improved calculation of inclusive leptonic decay rates is within reach.

Load-bearing premise

The load-bearing approximation is the electro-quenched one, which neglects the disconnected sea-quark contribution to the electromagnetic current without quantifying its effect; if that contribution is not negligible, the extracted $F_A$ and $F_V$ would be biased.

Editorial extensions

If this is right

  • Together with non-perturbative virtual-photon corrections, the extracted form factors enable $O(\alpha_{\rm em})$ predictions for leptonic decay rates of pseudoscalar mesons from the pion to the $B$ meson.
  • Structure-dependent corrections no longer have to be borrowed from chiral perturbation theory or from model assumptions, removing a systematic uncertainty from radiative decay phenomenology.
  • More precise leptonic decay rates translate into more precise determinations of the CKM matrix elements, tightening Standard Model consistency tests.
  • Hard-photon emission rates for heavy $D$ and $B$ mesons, previously only available from model-dependent predictions, become calculable from first principles and directly comparable with experiment.

Reading between the lines

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

  • A direct extension is to repeat the same ratio extraction with the disconnected sea-quark diagram of the electromagnetic current included, which would quantify the electro-quenched bias left unestimated by this paper.
  • The same boundary-condition ratio technique should transfer to $B$ mesons and to other radiative processes with one final-state photon, where analogous structure-dependent form factors appear.
  • Experimental hard-photon spectra for kaon and $D_s$ radiative leptonic decays could serve as a direct cross-check: if the lattice form factors reproduce the measured photon-energy distribution, the method is validated; if not, the discrepancy would point to missing contributions or new physics.
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

3 major / 5 minor

Summary. This Lattice 2019 proceedings paper reports a non-perturbative lattice QCD calculation of the structure-dependent form factors F_A and F_V entering the radiative leptonic decay amplitude P -> l nu_l gamma. The authors extend their earlier strategy for O(alpha_em) corrections to leptonic decays by extracting, from Euclidean correlator ratios, the infrared-safe parts of the hadronic tensor after subtracting the point-like contribution. Twisted boundary conditions are used to tune arbitrary meson and photon three-momenta. Numerical results are presented for kaons and D_s mesons from one lattice ensemble at a = 0.0619 fm with an unphysical light-quark mass (m_pi = 255 MeV), including plateau plots and comparisons with chiral perturbation theory predictions. The paper concludes that, with moderate statistics, the form factors can be extracted with good precision and their momentum dependence studied.

Significance. If the quantitative uncertainties can be established, this method is genuinely valuable: it would give first-principles, model-independent predictions for radiative leptonic decays of light and heavy mesons, with direct impact on CKM determinations and on tests of the Standard Model, particularly for D and B decays where current predictions are model-dependent. The technical construction is a clear strength: the decomposition of the amplitude in Eq. (2.2) separates the infrared-divergent point-like term from the finite structure-dependent form factors, and the use of twisted boundary conditions to scan the kinematical range is elegant. The comparison with tree-level chiral perturbation theory is used as an external benchmark rather than as an input, so the extraction is not circular. The main weakness is that the reported data are explicitly preliminary and lack statistical uncertainties on the extracted form factors, so the central precision claim is not yet supported.

major comments (3)
  1. [Sec. 4 and concluding paragraph] The central claim that "with moderate statistics, it is possible to extract with good precision" the form factors is not supported by the data as displayed: Figs. 2, 3, and 4 show plateau curves and extracted values without any statistical error bars or uncertainty bands, and no numerical values with errors are quoted in the text. Without uncertainties one cannot distinguish a clean signal from noise-dominated fluctuations, nor assess whether the momentum dependence visible in Fig. 4 is statistically significant. Please include statistical errors on every displayed quantity (ideally with a covariance estimate), or explicitly restrict the conclusion to a demonstration of feasibility.
  2. [Sec. 4, lattice setup] All numerical results are obtained at a single lattice spacing (a = 0.0619 fm) and a single unphysical light-quark mass (m_ud(2 GeV) = 11.7 MeV, M_pi = 255 MeV), with no continuum or chiral extrapolation. The paper itself labels the results as preliminary and postpones "a detailed analysis of all the systematics"; under these conditions the quantitative claim of good precision and the comparison to chiPT predictions in Fig. 4 cannot yet validate the conclusion beyond a proof of principle. Either add explicit estimates of discretization, finite-volume, and quark-mass effects, or rephrase the conclusion as a feasibility statement.
  3. [Sec. 3, Fig. 1 (right panel)] The electro-quenched approximation, namely the neglect of the disconnected sea-quark diagram in the electromagnetic-current correlator, is stated explicitly but its effect on F_A and F_V is never quantified. Since the paper aims at first-principles non-perturbative results at O(alpha_em), the size of this bias is load-bearing for the physical interpretation of the extracted form factors. Please provide at least a crude estimate of the disconnected contribution, for example by comparing with a stochastic or perturbative evaluation, or by quoting an order-of-magnitude bound from a partially connected calculation.
minor comments (5)
  1. [Sec. 4] The list of lattice spacings "a[fm] = 0.0085(36), 0.00815(30), 0.0619(18)" contains an apparent typographical error: the first entry has a relative uncertainty far larger than any reasonable lattice-spacing determination, and the first two entries appear inconsistent with the ensemble values cited from ref. [13]. Please correct the numbers.
  2. [Sec. 4] The sentence about the momentum combinations contains a typo ("assignements") and it is not clear whether the 100 combinations are per quark-mass choice or in total; please clarify.
  3. [Eq. (3.1)] The notation in the integrand, with expressions like "eEγty" and vectors written as "yyy", is non-standard and hard to read; the exponent should be written as exp(E_gamma t_y - i k . y + i p . x) with the integration variable named consistently.
  4. [Abstract and Introduction] The abstract states that the non-perturbative calculation will allow accurate predictions at O(alpha_em) for the inclusive leptonic decay rates, but this paper reports only form factors, not rates. Please clarify that the rates are the eventual target of the ongoing program and are not computed here.
  5. [Figure captions] The captions of Figs. 2-4 should state that the data are preliminary, at a = 0.0619 fm and unphysical quark masses, and should define whether error bars are present; currently the absence of uncertainties is not mentioned.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: FA/FV are extracted from Euclidean correlators with an algebraic fP subtraction; the chiPT comparison is an external benchmark.

full rationale

The paper derives FA and FV from lattice Euclidean correlators through the estimators RA(t) and RV(t) in Eqs. (3.3)-(3.4), followed by plateau fits. No form-factor parameter is tuned to reproduce the final FA/FV values; the only subtraction, 2 fP/(mP x_gamma), uses the decay constant fP measured independently from the pseudoscalar-axial two-point function. This is an algebraic removal of a known term, not a fit of FA to the point-like contribution. The comparison with chiral perturbation theory uses external low-energy constants from ref. [18] and serves as a benchmark, not as an input to the lattice extraction. Self-citations to refs. [8]-[13] supply the gauge ensembles and the overall strategy, but the numerical values of FA and FV do not reduce to those citations; the cited results are not used as a substitute for the computed form factors. The paper's preliminary status and the absence of displayed statistical uncertainties are limitations of the evidence, but they are not circularity.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

The paper introduces no new free parameters fitted to data for the central derivation and no invented entities. It relies on standard Lorentz decompositions, the analytic continuation of correlators, and the electro-quenched approximation, the latter being the most questionable assumption because its numerical effect is unquantified.

assumptions (5)
  • domain assumption The Euclidean correlator can be analytically continued from Minkowski space and converges because no intermediate state is lighter than the pseudoscalar meson.
    Converts the Minkowski amplitude Eq (2.1) to the computable Euclidean correlator Eq (3.1); if this continuation fails, the extraction of H is invalid.
  • domain assumption Neglecting the disconnected sea-quark contribution to the electromagnetic current (electro-quenched approximation) is a valid approximation for these form factors.
    Explicitly stated as a limitation; its effect on FA and FV is not quantified, so the central claim rests on this approximation being small.
  • standard math The Lorentz decomposition of the amplitude Eq (2.2) into H1, H2, FV, FA plus the point-like term is complete, and the point-like term saturates the Ward identity.
    Standard decomposition taken from refs [8,15]; it is assumed without derivation and is the basis for separating structure-dependent terms.
  • domain assumption Twisted boundary conditions on the two quark propagators connected to the electromagnetic current yield arbitrary meson and photon momenta with only exponentially suppressed unitarity violations.
    Enables the momentum setup used in the calculation; relies on the formalism of refs [16,17] and the smallness of the induced symmetry violations.
  • domain assumption The photon emission from the final-state lepton can be treated perturbatively using the decay constant fP, so only the hadronic tensor needs a non-perturbative calculation.
    Underlies the whole separation at O(alpha_em); standard in the literature but not derived in this paper.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Real photon emissions in leptonic decays." pith.science (2026). https://pith.science/paper/SXXUH5KT

@misc{pith2026190810160,
  author       = {Pith},
  title        = {Pith review of: Real photon emissions in leptonic decays},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/SXXUH5KT}},
  note         = {Machine review of arXiv:1908.10160}
}
abstract

We present a non-perturbative calculation of the form factors which contribute to the amplitudes for the radiative decays $P\to \ell \bar \nu_\ell \gamma$, where $P$ is a pseudoscalar meson and $\ell$ is a charged lepton. Together with the non-perturbative determination of the virtual photon corrections to the processes $P\to \ell \bar \nu_\ell$, this will allow accurate predictions to be made at $O(\alpha_{em})$ for leptonic decay rates for pseudoscalar mesons ranging from the pion to the $B$ meson. We are able to separate unambiguously the point-like contribution, the square of which leads to the infrared divergence in the decay rate, from the structure dependent, infrared-safe, terms in the amplitude. The fully non-perturbative, $O(a)$ improved calculation of the inclusive leptonic decay rates will lead to significantly improved precision in the determination of the corresponding Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. Precise predictions for the emission of a hard photon are also very interesting, especially for the decays of heavy $D$ and $B$ mesons for which currently only model-dependent predictions are available to compare with existing experimental data.

Figures

Figures reproduced from arXiv: 1908.10160 by the authors.

Figure 1
Figure 1. The connected diagram on the left shows our choice of the spatial boundary conditions. By treating the two propagators attached to the electromagnetic current as two different flavours, with the same mass and electric charge but different boundary conditions, we may choose arbitrary values for the meson and photon spatial momenta.The diagram on the right represents the contribution associated with the emission of th… view at source ↗
Figure 2
Figure 2. Examples of plateaux fits for the ratios RA(t,T/2) (left) and RV (t,T/2) (right). For t 0 we get the following numerical estimators for the form-factors RA(t) = mP 4p · k ∑ r=1,2 ∑ j=1,2 R jr A (t; p, k) ε j r →  FA(xγ ) + 2 fP mPxγ  , RV (t) = mP 4 ∑ r=1,2 ∑ j=1,2 R jr V (t; p, k) i [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. The extracted value of RA(xγ ), Eq. (3.5), as a function of xγ for the K meson (left) and for the Ds meson (right). The (red) squares represent the point-like contribution given by 2 fP/(mPxγ ). 0 0.5 1 1.5 -0.2 -0.1 0 0.1 0.2 0.3 0 0.5 1 1.5 -0.1 -0.05 0 0.05 0.1 0.15 0.2 [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: The extracted value of the kaon form factors FA(xγ ) (left) and FV (xγ ) (right) as a function of xγ . The (red) lines correspond to the χPT predictions obtained by using the formulae discussed in the text. MD = 1933(50) MeV, MK = 535(14) MeV and Mπ = 255(7) MeV. Simil…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

18 extracted references · 5 canonical work pages

  1. [1]

    Tanabashi et al

    M. Tanabashi et al. [Particle Data Group], Phys. Rev. D 98 (2018) no.3, 030001

    2018

  2. [2]

    C. Kane, C. Lehner, S. Meinel and A. Soni, arXiv:1907.00279 [hep-lat]

    work page Pith review arXiv 1907

  3. [3]

    $\pi \to l\nu \gamma$ Form Factors at Two-loop

    J. Bijnens and P. Talavera, Nucl. Phys. B 489 (1997) 387 [hep-ph/9610269]

    work page Pith review arXiv 1997

  4. [4]

    C. Q. Geng, I. L. Ho and T. H. Wu, Nucl. Phys. B 684 (2004) 281 [hep-ph/0306165]

    arXiv 2004

  5. [5]

    Mateu and J

    V . Mateu and J. Portoles, Eur. Phys. J. C 52 (2007) 325 [arXiv:0706.1039 [hep-ph]]

    arXiv 2007

  6. [6]

    Unterdorfer and H

    R. Unterdorfer and H. Pichl, Eur. Phys. J. C 55 (2008) 273 [arXiv:0801.2482 [hep-ph]]

    arXiv 2008

  7. [7]

    Cirigliano, G

    V . Cirigliano, G. Ecker, H. Neufeld, A. Pich and J. Portoles, Rev. Mod. Phys.84 (2012) 399 [arXiv:1107.6001 [hep-ph]]

    arXiv 2012

  8. [8]

    Carrasco et al., Phys

    N. Carrasco et al., Phys. Rev. D 91 (2015) no.7, 074506 [arXiv:1502.00257 [hep-lat]]

    arXiv 2015

Show all 18 references
  1. [9]

    Lubicz et al

    V . Lubicz et al. , Phys. Rev. D95 (2017) no.3, 034504 [arXiv:1611.08497 [hep-lat]]

    2017 arXiv

  2. [10]

    Lubicz et al., PoS LATTICE 2016 (2016) 290 [arXiv:1610.09668 [hep-lat]]

    V . Lubicz et al., PoS LATTICE 2016 (2016) 290 [arXiv:1610.09668 [hep-lat]]

    2016 arXiv

  3. [11]

    Tantalo et al., arXiv:1612.00199 [hep-lat]

    N. Tantalo et al., arXiv:1612.00199 [hep-lat]

    arXiv

  4. [12]

    Giusti et al., Phys

    D. Giusti et al., Phys. Rev. Lett. 120 (2018) no.7, 072001 [arXiv:1711.06537 [hep-lat]]

    2018 arXiv

  5. [13]

    Di Carlo et al., Phys

    M. Di Carlo et al., Phys. Rev. D 100 (2019) no.3, 034514 [arXiv:1904.08731 [hep-lat]]

    2019 arXiv

  6. [14]

    Becirevic, B

    D. Becirevic, B. Haas and E. Kou, Phys. Lett. B 681 (2009) 257 [arXiv:0907.1845 [hep-ph]]

    2009 arXiv

  7. [15]

    Bijnens, G

    J. Bijnens, G. Ecker and J. Gasser, Nucl. Phys. B 396 (1993) 81 [hep-ph/9209261]

    1993 arXiv

  8. [16]

    G. M. de Divitiis, R. Petronzio and N. Tantalo, Phys. Lett. B 595 (2004) 408 [hep-lat/0405002]

    2004 arXiv

  9. [17]

    C. T. Sachrajda and G. Villadoro, Phys. Lett. B 609 (2005) 73 [hep-lat/0411033]

    2005 arXiv

  10. [18]

    Bijnens and G

    J. Bijnens and G. Ecker, Ann. Rev. Nucl. Part. Sci. 64 (2014) 149 [arXiv:1405.6488 [hep-ph]]. 6

    2014 arXiv

Pith tools

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