Charmed and ϕ meson decay constants from 2+1-flavor lattice QCD

On a lattice with 2+1-flavor dynamical domain-wall fermions at the physical pion mass, we calculate the decay constants of , , and . The lattice size is , which corresponds to a spatial extension of fm, with a lattice spacing of fm. For the valence light, strange, and charm quarks, we use overlap fermions at several mass points close to their physical values. Our results at the physical point are MeV, MeV, MeV, MeV, and MeV. The couplings of and to the tensor current ( ) can be derived from ratios and , respectively, which are the first lattice quantum chromodynamics (QCD) results. We also obtain ratios and , which reflect the size of heavy quark symmetry breaking in charmed mesons. Ratios and can be taken as a measure of SU(3) flavor symmetry breaking.


I. INTRODUCTION
f P f P Meson decay constants are important nonperturbative quantities for the study of meson leptonic decays, and their results from lattice quantum chromodynamics (QCD) have received considerable attention. The pseudoscalar meson decay constants ( ) can be effectively used to determine the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, if combined with experimental measurements of the corresponding leptonic decays. The newest lattice QCD average of can be found in the review by the Flavor Lattice Averaging Group (FLAG) [1].
In principle, vector meson decay constants, , can also be used to determine the CKM matrix elements; however, experimental measurements of the leptonic decays of vector mesons are much more difficult than those of pseudoscalar mesons due to small branching ratios. With increasing statistics, the leptonic decay of may be expected to be measured via BES-III or Belle II for the f D * s first time in the near future for a vector meson [2]. Then, the comparison of from experimental and theoretical calculations can be used to study the low energy properties of QCD.
Furthermore, the decay constants of heavy-light vector mesons can be used to test the accuracy of the heavy quark effective theory (HQET). By neglecting the terms of , where is the heavy quark mass, one can obtain [3] from the leading order QCD calculation, which implies that the ratio approaches one since the strong coupling constant vanishes in the infinite heavy quark mass limit. We can obtain the corrections from the higher order terms in charmed mesons through from lattice QCD calculations. Moreover, the ratios for charmed mesons are input parameters for QCD factorization studies of charmed nonleptonic B meson decays [4,5]. Another important quantity, , is the coupling of a vector meson to f T V / f V the tensor current. The nonperturbative determination of ratio is important in the light cone QCD sum rule (LCSR) calculations of form factors in B to vector meson semileptonic decays (see discussions in [6][7][8]).
In this paper, we present a lattice calculation of , , and meson decay constants in a lattice setup with chiral fermions, which are usually expected to be important when light flavors are involved, since chiral symmetry is a fundamental property of QCD. We use overlap fermions for valence quarks and carry out the calculation on 2+1-flavor domain wall fermion gauge configurations generated by RBC-UKQCD Collaborations. The lattice size is big enough ( fm) to prevent large finite volume effects. The light sea quark mass is almost at the physical point. There have been four lattice QCD calculations of in the literature so far. Two of them were performed on 2-flavor gauge ensembles [9,10]. The other two were performed on 2+1-flavor ensembles [2] and 2+1+1-flavor ensembles [11], respectively. An unexpected large quenching effect of the strange quark was observed in and from the 2-flavor result [9] (confirmed in [12] but with a reduced effect), whereas the 2-flavor result from [10] demonstrated a much less pronounced effect. In this study, we develop an independent 2+1-flavor calculation for for comparison with the aforementioned calculations.
The rest of this paper is organized as follows. In Sec. II, we present our framework of the calculation, including the definitions of the decay constants and the lattice setup. Sec. III presents the details of the analyses and the numerical results and discussions. Finally, we summarize our work in Sec. IV.

II. DEFINITIONS AND LATTICE SETUP
A. Decay constants of pseudoscalar and vector mesons f P The decay constant, , of a pseudoscalar meson P is defined as being the momentum of the meson. Using the partially conserved axial vector current (PCAC) relation, we can obtain from the matrix element of the pseudoscalar density is the polarization vector of meson with helicity . We use the local vector current on the lattice to compute the above matrix element for convenience; however, we must also calculate the finite renormalization constant for the local current, which was obtained nonperturbatively in Ref. [13] for our lattice setup.
In addition to , vector mesons have another decay constant, , which is defined through the following matrix element of the tensor current Here, in the tensor current, . Since the tensor current has a nonzero anomalous dimension, we will give values of in the commonly used scheme and at a scale of GeV. The matching factor from the lattice to the continuum scheme for the tensor current was presented in Ref. [13].

B. Lattice setup
Our calculation is carried out on the gauge configurations of domain wall fermions generated by RBC-UKQCD Collaborations [14]. We use the gauge ensemble named 48I, with a lattice size of and a pion mass of MeV from the sea quarks. The lattice spacing was determined to be [14]; thus, the spatial extension of the lattice is approximately . The parameters of the configurations are given in Table 1.
We use overlap fermions for valence quarks (as in [13]) to perform a mixed action study. The mismatch of the mixed valence and sea pion masses between the domain-wall fermion and the overlap fermion, measured us- am phy c Table 1. Parameters of the gauge configurations used in this work. are the valence quark mass parameters in lattice units, and the corresponding pion masses (in MeV) are from Ref. [15]. The physical charm quark mass, , is estimated to be approximately 0.73 (see below). ing , is GeV 4 [16], which is very small, reflecting a small partial quenching effect. The multimass algorithm of overlap fermions [17] permits calculations of multiple quark propagators at a reasonable cost. We calculate propagators with a range of masses from the light to charm quark on 45 configurations. The valence quark masses, , in lattice units are given in Table 1. The deflation algorithm is adopted to accelerate the inversion by projecting the 1000 low eigenvectors (including zero modes) of the overlap Dirac operator, which are calculated explicitly beforehand.
We use four mass parameters, , (as listed in Table 1) for the light valence quarks for chiral interpolation. The corresponding pion masses range from 114 MeV to 208 MeV [15]. Two strange quark mass parameters are used to extrapolate to the physical strange quark mass point. The bare charm quark masses that we use are approximately in lattice units, which are not small. For chiral lattice fermions, although the discretization error due to the heavy quark mass starts at , it could still be large. Thus, we shall try to estimate the finite lattice spacing effects in our results for D-mesons.

C. Two-point correlators
The matrix elements in Eqs. (1), (3), and (4), using which the decay constants are defined, can be derived directly from the related two-point functions, with the currents being the sink operators. Since the mesons involved in this study are all ground state hadrons, for the matrix elements to be determined precisely, it is desirable that the two-point functions be dominated by the contribution from the ground states. In this work, we adopt the Coulomb wall-source technique. In other words, we perform the Coulomb gauge fixing to the gauge configurations first and then, calculate the two-point functions using the following wall-source operators, which are obviously gauge dependent, where , and for pseudoscalar mesons and ( ) for vector mesons. Based on our experience [15], in addition to the suppression of excited states, the choice of wall source can also suppress the Pwave scattering states in the vector channels.
For the sink operators, we use spatially extended operators, , by splitting the quark and anti-quark field with spatial displacement , namely, . The operators with the same spatial separation are averaged to guarantee the correct quantum number and also to increase the statistics as a by-product. Thus, the two-point functions that we calculate are and where is the number of values with the same . Two-point functions with different r can be calculated simultaneously without expensive extra inversions. After the insertion of the intermediate states, the spectral expression of a two-point function is where is proportional to the Bethe-Salpeter amplitude for the n-th state. Since the r dependences of are different for different states in each channel, a proper linear combination of several values with different r may yield an optimal twopoint function , which is dominated by the ground state.
Obviously, the parameterization of Eq. (9) reveals that spectral weight is proportional to the matrix element that defines the decay constant of a specific meson state. However, to obtain the decay constant, we need to remove factor , which is the matrix element of the wall-source operator, , between the vacuum and the meson state and can be derived from the wall-to-wall correlation function

A. Meson masses
To extract the meson masses, we apply two fitting strategies. One strategy is applying correlated simultaneous fittings to the correlation functions with different r using one (for vector mesons) or two (for pseudoscalar Charmed and ϕ meson decay constants from 2+1-flavor lattice QCD Chin. Phys. C 45, 023109 (2021) 023109-3 mesons) mass terms. The function form used in the simultaneous fits is , and the second term in the brackets on the right hand side is obtained from the propagation of the correlator in the negative time direction. and are fitted with the minimum method. We vary the number of mass terms to two or three and check the stability of the fitting results. Within statistical uncertainties, the fitted ground state mass, , does not depend on the number of mass terms. The upper limit of fitting range is chosen using the following criteria. For the pseudoscalar channel, is fixed to the maximum value where the relative errors of correlators satisfy 5%. For the vector channel, is chosen by requiring 10%. The lower limit of the fitting range is varied in a wide range during fitting, and we check the stability of the results. Then, among all the fittings that have and yield a consistent ground state mass, we choose the earliest to obtain our final results. The uncertainties are obtained from Jackknife analyses to take into account the correlations among the data as we repeat the fitting for each Jackknife ensemble.
In the left panel of Fig. 1, we show fitted ground state mass in lattice units as a function of . Here, we finally choose the fitting range [11,18] for the D meson. In the right panel of Fig. 1, the obtained ground state mass, , (the band in the graph) is compared with the corresponding effective masses, , from various correlators with different r. The data points in the magenta squares are from the correlator with ( and permutations averaged). The ones in the blue triangles are from the local sink correlator with . The ones in the black circles are the effective masses from a combination of two correlators where we can tune parameter and use various to make the effective mass plateau from appear as early as possible. This leads to our second fitting strategy. Different states with the same quantum number contribute differently to correlators , and these contributions vary as r varies. Thus, it is possible to find a large r such that the contribution of the lowest excited state to cancels that to , and is dominated by the ground state.
In the right panel of Fig. 1, the black circles present a mass plateau that starts much earlier than that from correlator or . Therefore, we can fit combined correlator easily with a single exponential term. We check that this fitting yields stable and consistent ground state mass as we vary parameter . We also confirm that the results from the above two fitting strategies are consistent. The fitting results of and from the two strategies with are summarized in Table 2 for comparison. The second strategy yields smaller statistical uncertainties, since the mass plateau from the combined correlator appears earlier; thus, data points with less errors are used in fittings. Similar advantages of the second strategy are observed in the analyses of other meson masses; therefore, we adopt strategy II to obtain meson masses in the following.
The results of the pion and kaon masses are shown in Table 3. The pion mass and the combination are used to fix the physical up (degenerate with the down quark) and strange quark mass, respectively. From Table 3, we can observe that is independent of the pion mass (or equivalently, the up/down quark mass) within the statistical uncertainties. This is exactly what we expect from the lowest-order analysis of the chiral perturbation theory and is the reason we use this combination. The results of the meson masses and decay constants will be interpolated/extrapolated to the physical point where and by using MeV and MeV [18]. Here, the uncertainties come from the error of the lattice spacing. Since these uncertainties are much smaller than our statistical error or the discretization error, as we will observe later, we ignore them in our estimate of the systematic uncertainty. Table 4, we collect the masses of and at our valence quark masses. From the data, we can observe that the mass of barely depends on the light quark mass with our current statistical uncertainties. To obtain at the physical point, we use the following interpolation/extrapolation form: and . This is the Taylor expansion around the physical and strange quark masses, and we keep only the lowest order, i.e., the linear terms, since our quark masses are close to their physical values. Then, we obtain m phys K * = 895(10) MeV, where the error includes the statistical/fit uncertainty and the uncertainty of the lattice spacing. Parameter from the fitting is consistent with zero within uncertainty, as expected from the raw data.  Table 4. For the corresponding at each of the two strange quark masses, we use the average of the four values in the last column of Table  3. This extrapolation yields  Vector mesons can decay to two pseudoscalar mesons through P-wave. On our lattice, the minimal nonzero momentum is 226 MeV, which is not small. The thresholds of P-wave decays for , , and mesons are not open on our lattice, but can decay to on our lattice. We observed mass plateaus for the meson, but not for the scattering states of , which we believe are suppressed by the usage of the Coulomb gauge wall source while calculating the 2-point functions [15]. The agreement of and (from our interpolation/extrapolation) with their experimental values indicates that it is safe to ignore the threshold effects at our current precision.  Table 5. We use the experimental value of (together with in the above) to set the physical charm (and strange) quark mass. With our lattice spacing, we obtain by using MeV from the Particle Data Group (PDG2018) [18]. We use the following function, similar to Eq. (13), to interpolate/extrapolate to the physical strange and charm quark mass point:     (7) 0.0030 0.503 (7) 0.0060 0.504 (6) 0.0650 0.579(5) 0.127 (7) 0.0017 0.514 (7) 0.0024 0.512 (7) 0.0030 0.511 (7) 0.0060 0.512 (7) Charmed and ϕ meson decay constants from 2+1-flavor lattice QCD Chin. Phys. C 45, 023109 (2021) 023109-5 where , and is another free parameter. From this, we obtain which agrees with the experimental value of GeV [18]. The interpolation/extrapolation is shown in Fig. 2. The function in Eq. (16) can describe the data very well. The dependence of on the strange quark mass is relatively small. Therefore, the slope of the straight lines in the left plot of Fig. 2 is small. This is also why the two lines in the right plot are very close to each other. The dependence on the charm quark mass is apparent. The position of the physical point in the left plot indicates that the physical charm quark mass is approximately .
The masses of D and mesons are listed in Table 6. The following ansatz is used to interpolate/extrapolate our numerical results to the physical quark mass point: Here, the term appears because our lattice results, , are not calculated at the physical strange quark mass, and is used to set the physical charm quark mass. We obtain for the two mesons, respectively, after the interpolations/ extrapolations. Our D meson mass agrees with the PDG2018 value of GeV within . However, our meson mass is heavier than the PDG2018 value of GeV by approximately 1%. Thus, we estimate the discretization error associated with the large charm quark mass to be approximately 1% in our results for the charmed meson masses.

Renormalization constants
Before we delve into the data analyses for the meson decay constants, first, we present the renormalization constants (RCs) for the local vector current and the tensor current. The RCs of the quark bilinear operators for our lattice setup (overlap fermions on domain-wall fermion configurations) were calculated nonperturbatively in Refs. [13,19]. For the 48I ensemble used in this work, we employed both the RI/MOM and the RI/SMOM schemes to calculate those constants nonperturbatively [13]. The matching factors to the scheme for the local axial vector current, , and for the tensor current (at a scale of 2 GeV) are listed in Table 7. Because we use chiral fermions, we obtain , which was also confirmed numerically in Ref. [13].

and
To obtain decay constants and , we perform simultaneous fits to the wall-to-point ( ) and wall-to-wall ( ) correlators for a given meson, M. These fittings are with two exponentials, and the ground state mass is constrained within to its fitted result from the above strategy II as we determined the meson masses. After removing the matrix element of source operator from the spectral weight of , we obtain and then, decay constants by using Eqs. (2), (3) and the fitted meson mass. This fitting and calculation process is repeated for each Jackknife sample to obtain the statistical uncertainty of . For obtained from the local vector current, we need to multiply it with , as discussed in Section IIIB.1. In the following, we use a superscript "bare" to indicate decay constants obtained directly from the local vector current.  Table 4. The two center values are almost the same, and our statistical uncertainty is large (~6%). Thus, it is difficult to tell the strange quark mass dependence of . If we perform a constant fit to the two numbers, then we obtain or MeV after multiplying it with and . If we perform a linear extrapolation to physical strange quark mass point , then we obtain a value of MeV. We choose the value with a larger error from the constant fit as our result at the physical point. Therefore, we obtain We use to estimate our discretization error due to the large charm quark mass since we cannot extrapolate to the continuum limit with only one lattice spacing. The decay constant in lattice units for all charm and strange quark masses are given in Table 5. One can use the function form given in Eq. (16) (replacing with ) to extrapolate/interpolate our lattice results in Table  5 to the physical charm and strange quark mass point. We find the following: The difference in the center values of calculated in this work and in our previous work (254(2)(4) MeV) [20] is 5 MeV or 2%. Since our previous result was obtained in the continuum limit, we treat this 2% difference as an estimate of the discretization error and assign it to all our decay constants for the charmed mesons in this work. The vector meson decay constant, , from our lattice data is given in the sixth column of Table 5. Again, we use the function form given in Eq. (16) (replacing with ) to extrapolate/interpolate our lattice results to the physical charm and strange quark mass point. The fitting is shown on the left panel of Fig. 3. Compared with the case of the quark mass dependence of is difficult to see relative to the large statistical errors. . Multiplying this number with GeV and , we find MeV. Here, the uncertainty includes the errors from the statistics, extrapolation/interpolation, lattice spacing, and . If we assign a 2% discretization error, then we finally obtain  Table 7. Matching factors to the scheme for the local axial vector current and for the tensor current [13].

1.055(31) MS
Then, we use Eq. (16) and Eq. (18) to interpolate/extrapolate our raw data to the physical quark mass point for and respectively. After multiplying the results with renormalization factor in the scheme and assigning a 2% discretization uncertainty, we find at the scale of 2 GeV. Here the first uncertainty includes the errors from statistics and interpolation/extrapolation and the error of and is dominated by the error of the renormalization factor. The second uncertainty is from the finite lattice spacing effect.

IV. SUMMARY
We calculated decay constants , , and of the charmed and light mesons, including and , by using 2+1-flavor domain wall fermion gauge configurations at one lattice spacing. The valence overlap fermion has 4, 2, and 4 mass values, respectively, for the light, strange, and charm quarks. We use the experimental values of , , and to set the physical light, strange, and charm quark masses. The masses of D, , , and at the physical point are found via interpolation/extrapolation using the lowest order of Taylor expansion (i.e., a linear interpolation/extrapolation) since our valence quark masses are close to their physical values.
Masses , , , and , obtained from our lattice calculation, are in good agreement with their experimental measurements. The mass we found is 1% higher than its experiment value. The center value of from this calculation is 2% away from our previous lattice QCD calculation extrapolated to the continuum limit [20]. Thus, we estimate the discretization uncertainty in this work to be around 2%.
The final results of this work for the decay constants are given in Eqs. (20), (22)-(26), (28). By quadratically adding the statistical/fitting uncertainty and the systematic uncertainty, we obtain the decay constants, shown in Table 9, and some of their ratios, shown in Table 10. For     In Fig. 5, we compare and the ratio from this work and other lattice QCD calculations [2,[9][10][11][12]. The values from the 2+1-flavor and 2+1+1-flavor simulations are consistent. There might be tension between the 2-flavor calculations and other calculations, including the dynamical strange quark. This may reflect an unexpected large quenching effect from the strange quark. However, the 2-flavor calculation of in [10] shows that this quenching effect is not as significant as that observed in [9]. The two calculations employ different lattice actions of the two-flavor theory. The computation in [12] is performed on the same 2-flavor gauge ensembles as those used in [9] and employs the same analysis method as that used in [11]. It yields a with a smaller strange quark quenching effect and, therefore, is more in agreement with [10]. Thus, more lattice QCD calculations, especially those with two dynamical flavors, are certainly welcome to clarify this situation.
The ratios of decay constants of charmed mesons in Table 10 show that the size of heavy quark symmetry breaking is approximately 10%, whereas the size of SU(3) flavor symmetry breaking is approximately 17%. ϕ To better control the systematic uncertainty from discretization effects in our work, we need to perform our calculation at more lattice spacings in the future. We also need to include the quark-line disconnected diagram for the meson two-point function. To accurately estimate the threshold effects of strong decays of the vector mesons, further studies on larger volumes are necessary. ter (NERSC) for providing HPC resources that have contributed to the research results reported within this paper.