Methods for high-precision determinations of radiative-leptonic decay form factors using lattice QCD

TL;DR

This paper develops and compares lattice-QCD methods to accurately determine form factors for radiative leptonic decays of pseudoscalar mesons, achieving 5% precision across the entire kinematic range.

Contribution

It introduces an optimized lattice-QCD approach combining 3d sequential propagators and infinite-volume approximation to control systematic errors and compute form factors comprehensively.

Findings

Achieved 5% statistical and fitting uncertainties in form factors.

Demonstrated the effectiveness of 3d sequential propagators and infinite-volume methods.

Provided numerical results for $D_s^+ o au^+ u \gamma$ decay form factors.

Abstract

We present a study of lattice-QCD methods to determine the relevant hadronic form factors for radiative leptonic decays of pseudoscalar mesons. We provide numerical results for $D_s^+ \to \ell^+ \nu \gamma$. Our calculation is performed using a domain-wall action for all quark flavors and on a single RBC/UKQCD lattice gauge-field ensemble. The first part of the study is how to best control two sources of systematic error inherent in the calculation, specifically the unwanted excited states created by the meson interpolating field, and unwanted exponentials in the sum over intermediate states. Using a 3d sequential propagator allows for better control over unwanted exponentials from intermediate states, while using a 4d sequential propagator allows for better control over excited states. We perform individual analyses of the 3d and 4d methods as well as a combined analysis using both methods, and find that the 3d sequential propagator offers good control over both sources of systematic uncertainties for the smallest number of propagator solves. From there, we further improve the use of a 3d sequential propagator by employing an infinite-volume approximation method, which allows us to calculate the relevant form factors over the entire allowed range of photon energies. We then study improvements gained by performing the calculation using a different three-point function, using ratios of three-point functions, averaging over positive and negative photon momentum, and using an improved method for extracting the structure-dependent part of the axial form factor. The optimal combination of methods yields results for the $D_s^+ \to \ell^+ \nu \gamma$ structure-dependent vector and axial form factors in the entire kinematic range with statistical plus fitting uncertainties of order 5%, using 25 gauge configurations with 64 samples per configuration.