D*Dρ and B*Bρ strong couplings in light-cone sum rules

We present an improved calculation of the strong coupling constants and in light-cone sum rules, including one-loop QCD corrections of leading power with meson distribution amplitudes. We further compute subleading-power corrections from two-particle and three-particle higher-twist contributions at leading order up to twist-4 accuracy. The next-to-leading order corrections to the leading power contribution numerically offset the subleading-power corrections to a certain extent, and our numerical results are consistent with those of previous studies on sum rules. A comparison between our results and existing model-dependent estimations is also made.


Introduction
This study aims to provide a more precise determination of and coupling in the framework of light-cone sum rules (LCSR). These couplings describe the low-energy interaction between heavy and light mesons, which are important in understanding QCD longdistance dynamics. The coupling is a fundamental parameter of the effective Lagrangian in heavy meson chiral perturbative theory (HM PT) [1][2][3][4], which plays an important role in the study of heavy meson physics. Phenomenologically, it describes the strength of the final state interactions [5], which are important in the generation of the strong phase in B decays [6,7]. Moreover, the coupling relates the pole residue of to form factors at large momentum transfer by the dispersion relation, which is helpful to gain a better understanding of the behavior of form factors.

D(B)
ρ Various theoretical approaches have been suggested to evaluate the coupling. First, with the to form factors obtained in a certain region of the momentum transfer from LCSR [8][9][10] and the lattice QCD [11], the corresponding pole residue that relates to the coupling is extracted with appropriate extrapolation of the form factors. The simplest approach to perform this extrapola-χ tion is the vector meson dominance (VMD) hypothesis, which neglects continuum spectral characteristics. In addition to the VMD approximation, certain other modified parameterizations for the form factors were proposed in [12][13][14][15]. In addition to the above method, the coupling can also be estimated from HM PT with the VMD approximation, as performed in [16][17][18][19]. Another strategy to obtain the strong coupling is calculation from first principles of QCD. We study the strong coupling constants and in LCSR using the double dispersion relation. The LCSR was proposed in [20][21][22] based on the light-cone operator-product-expansion (OPE) relative to the conventional QCD sum rules (QCDSR) method. There have been several studies addressing the couplings and , starting from [23] with the inclusion of two-particle DAs corrections up to twist-3 at leading order (LO). Years later, [24][25][26][27] improved the formal calculation by considering the twoparticle twist-4 corrections [28] at LO. Meanwhile, is also calculated with three-point QCDSR by taking into account the dimension-5 quark-gluon condensate corrections under flavor symmetry [29]. In previous studies, coupling estimations showed a wide range of values. The results from LCSR studies [24][25][26][27] yield values that are smaller than those in other estimations; hence, we O(α s ) must perform a more accurate calculation to confirm whether this difference arises from radiative and power corrections. Moreover, the inclusion of NLO corrections can also decrease the scale-dependence. In this work, we present a calculation including corrections to the leading power contribution with the resummation of large logarithm to next-to-leading logarithmic accuracy (NLL). For the subleading-power corrections, our results also include the three-particle twist-4 corrections at LO.
The paper is organized as follows: in Section 2, we calculate the leading power contributions up to NLO. The procedures for analytic continuation and continuum subtraction are similar to the ones presented in [30,31]. Subleading-power corrections, including two-particle and three-particle corrections up to twist-4 at LO, are calculated in Section 3. Section 4 presents our numerical results and the phenomenological discussion. We summarize this work in the last section.

Hρ
The strong coupling constant is defined by . Here, we choose and H to represent the meson and meson, respectively, and is . and are the polarization vectors of and mesons, respectively. We use the conventions and . The couplings of different charge states are related by isospin symmetry, for example, We construct the following correlation function at the starting point , Q is the heavy quark field. We introduce two light-cone vectors and with . For the momentum, we choose a coordinate to have . Within the frame of the LCSR, the meson is on the light-cone, and its momentum is chosen as On the hadronic level, taking advantage of the following definitions for decay constants the correlation function (4) can be written as where the second term counts the contributions from higher resonances and continuum states. The ellipses denote the terms that vanish after the double Borel transformation.
After the double Borel transformation, we obtain the hadronic representation of the correlation function For the boundary of the integral , we take with as the threshold of excited and continuum states. The Borel parameters associated with and are quite similar in magnitude; hence, we set the same value .
On the quark level, the leading-twist tree diagram is displayed in Fig. 1. The correlation function reads ρ Using the definition of the leading-twist DA [28] in Appendix A, we obtain the leading twist tree level factorization formula The result (9) can be written in the form of double dispersion relation where the double dispersion density is defined as The expression of wave function in terms of Gegenbauer polynomials can be written as where is the function of Gegenbauer moments . The spectral density can be obtained as follows [30] We defined two variables in Eq. (13), which are similar to those in our previous study [32].
Equating (6) and (10) and applying the double Borel transformation, then subtracting the continuum states using the quark-hadronic duality, we obtain the leadingtwist strong coupling constant at LO

Next-to-leading order corrections
The one-loop diagrams are shown in Fig. 2. The loop calculations are similar to those in [32]. We borrow the results from [32] and take the asymptotic form of wave function in Eq. (12) as for feasibility of calculation. Subsequently, we adopt a similar procedure to derive the double dispersion spectral densities and apply the double Borel transformation and continuum subtraction to obtain the result. Here, we omit the detailed procedures for simplicity, and the explicit operations with final results are provided in our previous study [32].
Combing the NLO result from [32] with Eq. (14), we obtain the leading-twist sum rules up to , and the coupling constant reads where and For the scale dependence, we use and we obtain d d ln µ µ It is evident that our result is independent of the factorization scale at the one-loop level.

Subleading-power corrections at LO
We perform the subleading power corrections to coupling constants, which are involved with two-particle and three-particle corrections up to twist-4. The higher twist (up to twist-4) Q-quark propagator in the background field is adopted [33]: ρ For the two-particle corrections, inserting the first term of Eq. (20) into the correlation function (4) and using the DAs of meson, we gain the higher twist twoparticle corrections The twist-3 corrections are suppressed by , and twist-4 corrections are suppressed by compared with the leading-twist corrections (10). To maintain consistency, we must also consider the subleading-power corrections of leading-twist. According to Eq. (8), the correlation function reads Similar to the leading-twist LO case, after applying double Borel transformations and subtracting the continuum states using quark-hadronic duality, we obtain the strong coupling constant for two-particle subleadingpower corrections at LO Considering three-particle corrections at the tree qqg level, the correlation function of the three-particle corrections can be written as Employing the conformal expansion of the three-particle DAs in Appendix B, we obtain Collecting Eqs. (15), (23), and (25), the final LCSR reads F LT, (1) where is defined in Eq. (16).

MS ρ
The masses of quarks in the scheme and the values of decay constants are listed in Table 1. The parameters in DAs are listed in Table 2, which are extracted from the method of QCDSR [9,28].
The solution to the two-loop evolution of the Gegenbauer moment is where the explicit expressions of and can be found in [39]. The scale evolution of other non-perturbative parameters to leading logarithmic accuracy is listed in Appendix B.

Theory predictions
The coupling constants for and are listed in Table 3. It is apparent that the perturbative QCD corrections decrease the tree level leading power contributions of coupling constants by 10% and 20% for and , respectively. Meanwhile, contributions from tree level sub-leading power corrections have evident compensation effects, which increase the leading twist contributions by approximately 20% and 30% for and , respectively. We also note that three-particle sub-leading power corrections have more minor effects compared with two-particle sub-leading power corrections, particularly for the coupling of , which is consistent with the validity of OPE.
The uncertainties are also included for separate terms of coupling constants in Table 3. The uncertainties are estimated by varying independent input parameters, and the individual uncertainties are presented in Table 4. To estimate total uncertainties, we add them in quadrature to obtain the final results. From Table 4, we can see that the primary uncertainties stem from the Borel mass and parameter both for and . The second Gegenbauer moment refers to transversely polarized meso twist-4 LCDA, which is indicated in Appendix A, B. The quasi-distribution amplitudes (quasi-DAs), which are based on the large-momentum effective theory (LaMET) [43,44], along with the simulation from lattice QCD [45] are also capable of extracting meson LCDA [46,47] other than QCDSR. In the future, a more precise determination of the meson LCDA will be helpful to decrease the uncertainties of couplings.
Next, we compare our results of coupling constants with those obtained in other studies using sum rules, as shown in Table 5. Among them, the LCSR studies [24,26] did not take the perturbative QCD corrections and 3-particle corrections into account, and the QCDSR study considers the dimension-5 quark-gluon condensate cor- Table 1. Heavy quark masses [34,35] and decay constants [9,[36][37][38], with scale-dependent quantity given at GeV.

Central value
Next, we extract the coupling from form factors. The form factors and , which are related to , are defined as From the dispersion relation of the form factor the strong coupling relates to the pole of the form factors at the unphysical point . Then, we have the relation where is the tensor coupling of the meson, and it is defined as . (32) Using Eq. (31) and choosing the size , we extract several numerical values of the coupling from recent LCSR studies [9,10], which we list in Table 6.
As shown, our central value is smaller than the extrapolations from LCSR form factors. To explain this discrepancy, on the one hand, uncertainties from the parameterizations of form factors to obtain the un-physical singularity may be underestimated. On the other hand, for the double dispersion relation, the Borel suppression is not sufficient to neglect the isolated excitation contributions, according to the discussion in [48]. Apart from , the and form factors can also be used to extract the coupling. To the best of our knowledge, there are no relevant studies that address corresponding form factors in LCSR. A future study could verify whether different choices of form factors are consistent across the coupling extractions.
Finally, we analyze our study and compare it with other model-dependent studies. The HM PT effective Lagrangian to parametrize the coupling can be written as [1] where is a matrix for light meson nonet, and the heavy mesons are represented by the doublet field with the conventional normalization. The parameter . In the chiral and heavy quark limits, we have the following relation for the coupling We find the coupling at leading power from our value of . In Table 7, we compare this value with other model estimations. Similar to the estimations from form factors, our result is smaller than the model predictions.
One possibility for this discrepancy may be that the model predictions have potentially larger errors. Consequently, in principle, the LCSR is more credible. However, in addition to the possible influences of excitation contributions mentioned above, the values of NLO corrections to higher twists and the sub-sub-leading power contributions are unknown, which may likewise represent sources of the discrepancy. Improvements in the calculations of both the sum rules and models may help better understand such discrepancies in the future.

Conclusion
We compute the and strong couplings to subleading-power in LCSR. Long-distance dynamics are incorporated in the DAs. We calculate the corrections to leading power of the sum rules. The subleadingpower corrections are calculated at LO by accounting the two-and three-particle wave functions up to twist-4. The analytical results of the double spectral densities are obtained, and in accordance with the previous study, we also perform a continuum subtraction for the higher-twist corrections. The LO results are independent of the choice of the duality region because of the special form of the double spectral density.
The NLO corrections reduce the numerical tree-level results by 10% and 20% for and , respectively. However, the subleading-power corrections at LO can give rise to values about 20% and 30% for and , respectively. Summing all the contributions, our values are consistent with the predictions from previous sum rules studies. Moreover, we also predict the coup-λ χ ling in HM PT at leading power. The central value of our result is smaller than the existing model-dependent estimations. Possible explanations for the discrepancy are potentially large errors within models and influences of excitation contributions in LCSR. Moreover, radiative corrections to higher twists and sub-sub-leading power contributions may also be sources of error. A better understanding of this discrepancy will prove beneficial to shedding light on long-distance QCD dynamics in the future.
We are grateful to Prof. Cai-dian Lü and Prof. Yuelong Shen for helpful discussions.
126 55 For the twist-4 three-particle DAs, The eight parameters in three-particle DAs can be written as The scale evolution of non-perturbative parameters to leading logarithmic accuracy is as follows where and , being the number of flavors involved. The scale-dependence of Γ ω where is given by . (B7)