Properties of the $\eta_q$ leading-twist distribution amplitude and its effects to the $B/D^+ \to\eta^{(\prime)}\ell^+ \nu_\ell$ decays

The $\eta^{(\prime)}$-mesons in the quark-flavor basis are mixtures of two mesonic states $|\eta_{q}\rangle=|\bar u u+\bar d d\rangle/\sqrt 2$ and $|\eta_{s}\rangle=|\bar s s\rangle$. In the previous work, we have made a detailed study on the $\eta_{s}$ leading-twist distribution amplitude. As a sequential work, in the present paper, we fix the $\eta_q$ leading-twist distribution amplitude by using the light-cone harmonic oscillator model for its wave function and by using the QCD sum rules within the QCD background field to calculate its moments. The input parameters of $\eta_q$ leading-twist distribution amplitude $\phi_{2;\eta_q}$ at an initial scale $\mu_0\sim 1$ GeV are then fixed by using those moments. The sum rules for the $0_{\rm th}$-order moment can also be used to fix the magnitude of $\eta_q$ decay constant, which gives $f_{\eta_q}=0.141\pm0.005$ GeV. As an application of the present derived $\phi_{2;\eta_q}$, we calculate the transition form factors $B(D)^+ \to\eta^{(\prime)}$ by using the QCD light-cone sum rules up to twist-4 accuracy and by including the next-to-leading order QCD corrections to the twist-2 part, and then fix the related CKM matrix element and the decay width for the semi-leptonic decays $B(D)^+ \to\eta^{(\prime)}\ell^+ \nu_\ell$.


I. INTRODUCTION
The mixing of η and η ′ mesons is essential to disentangle the standard model (SM) hadronic uncertainties with the new physics beyond the SM.It involves the dynamics and structure of the pseudoscalar mesons that has two mixing modes η − η ′ and η − η ′ − G, both of which have important theoretical significance.These mixings are caused by the QCD anomalies and are related to the breaking of chiral symmetry.However, since the matrix element of the exception operator is mainly nonperturbative, it still has not been calculated reliably.One may turn to phenomenological studies to obtain useful information on the non-perturbative QCD theory [1][2][3].At present, the η − η ′ − G mixing mode has been studied in detail in Refs [4][5][6][7][8][9][10].As for the η − η ′ mixing model, one can investigate it by using two distinct schemes, namely the singlet-octet (SO) scheme and the quark-flavor (QF) scheme.These two schemes reflect different understandings of the essential physics and they are related with a proper rotation of an ideal mixing angle [3].Practically, a dramatic simplification can be achieved by adopting the QF scheme [11][12][13][14], especially, the decay constants in the quark-flavor basis simply follow the same pattern of the state mixing due to the OZI-rule.In the present paper, we shall adopt the QF scheme to do our analysis and to achieve a better understanding of the mixing mechanism between η and η ′ .
As key components of B(D) → η (′) semileptonic decays, the B(D) → η (′) transition form factors (TFFs) need to be precisely calculated, whose main contribution comes from the |η q -component (the |η s -component gives negligible contribution here, but will have sizable contribution for B s (D s ) decays [23]).By further assuming the SU F (3) symmetry, the TFFs f B(D)→η (′)  + in the QF scheme satisfy the relations [23,33] f where φ is the mixing angle between the |η q -component and the |η s -component.
The TFFs of the heavy-to-light transitions at large and intermediate momentum transfers are among the most important applications of the light-cone sum rules (LCSR) approach.Using the LCSR approach, a twopoint correlation function (correlator) will be introduced and expanded near the light cone x 2 → 0, whose transition matrix elements are then parameterized as the light meson's light-cone distribution amplitudes (LCDAs) of increasing twists [34][35][36][37].It is thus important to know the properties of the LCDAs.
In present paper, we adopt the light cone harmonic oscillator (LCHO) model for the η q leading-twist LCDA φ 2;ηq .The LCHO model is based on the well-known Brodsky-Huang-Lepage (BHL) prescription [38,39] 1 for the light-cone wavefunction (LCWF), which is composed of the spin-space LCWF and the spatial one.The LCDA can be obtained by integrating over the transverse momentum from the LCWF.The parameters of φ 2;ηq at an initial scale will be fixed by using the derived moments of the LCDA, which will then be run to any scale region via proper evolution equation.Its moments will be calculated by using the QCD sum rules within the framework of the background field theory (BFTSR) [40,41].The QCD sum rules approach suggests to use the nonvanishing vacuum condensates to represent the nonperturbative effects [42].The QCD background field approach provides a simple physical picture for those vacuum condensates from the viewpoint of field theory [43][44][45][46].It assumes that the quark and gluon fields are composed of background fields and the quantum fluctuations around them.And the vacuum expectation values of the background fields describe the non-perturbative effects, while the quantum fluctuations represent the calculable perturbative effects.As a combination, the BFTSR approach provides a clean physical picture for separating the perturbative and non-perturbative properties of the QCD theory and provides a systematic way to derive the QCD sum rules for hadron phenomenology, which greatly simplifies the calculation due to its capability of adopting different gauges for quantum fluctuations and background fields.Till now, the BFTSR approach has been applied for dealing with the LCDAs of various mesons, some recent examples can be found in Refs.[47][48][49][50][51].
The remaining parts of the paper are organized as follows.In Sec.II, we give the calculation technology for deriving the moments of η q leading-twist LCDA φ 2;ηq by using the BFTSR approach, give a brief introduction of the LCHO model of φ 2;ηq , and then give the LCSRs for the TFFs of the semi-leptonic decay B(D) + → η q ℓ + ν ℓ .In Sec.III, we first determine the parameters of φ 2;ηq , and then the TFFs, the decay width and the CKM matrix element of the semi-leptonic decay B(D) + → η (′) ℓ + ν ℓ will be discussed.We will also compare our results with the experimental data and other theoretical predictions.Sec.IV is reserved for a summary.

II. CALCULATION TECHNOLOGY
A. Determination of the moments ξ n 2;ηq of the ηq leading-twist LCDA using the BFTSR For the QF scheme, the physical meson states |η and |η ′ are related to the QF basis |η q = |ūu + dd / √ 2 and |η s = |ss by an orthogonal transformation [14], where φ is the mixing angle.For the QF basis, one has two independent types of axial vector currents J q µ5 (q = u, d) and J s µ5 , e.g.
Their corresponding decay constants are where p is the momentum of η (′) .By extending the matrix elements (5) as the non-local operators over the light cone, one can achieve the definition of the corresponding LCDA.The LCDA of the valence quark momentum fraction distribution of η (′) meson can be defined similarly as those of other mesons by spreading out the non-local operators on the light cone via a way with increasing twists.In Ref. [52], the J s µ5 has been adopted to study the properties of η s , and at present, we will focus on J q µ5 to study the properties of η q .
To determine the properties of LCDA, one can firstly calculate its moments.The η (′) meson leading-twist LCDA are defined as [9,10] where Ψ = (u, d, s) represents the triplet of light-quark fields in flavour space, z µ stands for the light-like vector, [z, −z] is the path-ordered gauge connection which ensures the gauge invariance of the operator, and φ 2;η (′) are leading-twist LCDAs of η (′) mesons with respect to the current whose flavour content is given by C j with j = (q, s), respectively.We have which are derived in SO scheme [9], where λ 8 is the standard Gell-Mann matrix and 1 is 3 × 3 unit matrix.
By doing the series expansion near z 2 → 0 on both sides of Eq.( 6), one will get 0| Ψ(0)C j / zγ 5 (iz where the n th -order moment of η (′) leading-twist LCDA φ 2;η (′) j (x) has been defined as As mentioned above, the η (′) meson has two distinct components, |η q and |η s .The |η s -component has been studied by using B s (D s ) → η (′) semi-leptonic decays.Similarly, the |η q -component can also be studied by using B(D) → η (′) semi-leptonic decays.Eqs. ( .By further comparing theoretical predictions with the possible data on the B(D) → η (′) semi-leptonic decays, we can inversely achieve useful information on the η q leading-twist LCDA.We will calculate the moments of the η q leading-twist LCDA φ 2;ηq (x), which have been defined in Eq.( 9), within the framework of the BFTSR approach.
Secondly, one can apply the operator product expansion (OPE) to deal with the correlator in the deep Euclidean region.In deep Euclidean region q 2 ≪ 0, after applying the OPE, the correlator (10) becomes in which Tr[C q C q ] = 1, (iz (x, 0) represent the u-and d-quark propagators move from 0 → x and x → 0, respectively.The right-hand side of the correlator is perturbatively calculable within the framework of BFTSR.Following the standard procedures and by using the MS scheme to deal with the infrared divergences, the correlator can be expressed as a expansion series over the basic vacuum condensates with increasing dimensions.Since the current quark masses of u and d quarks are quite small, contributions from the u-and d-quark mass terms can be safely neglected in the calculation.
Thirdly, the correlator can also be calculated by inserting a complete set of the intermediate hadronic states in the physical region.By using the conventional quarkhadron duality [42], the hadronic expression of the correlator can be written as ImI (n,0) 2;ηq ,Had (q 2 ) = πδ(q 2 − m2 ηq )f 2 ηq ξ n 2;ηq | µ ξ 0 2;ηq | µ where µ represents the initial scale at which the η q leading-twist LCDAs have been defined.Because of the SU (3) F flavour symmetry, here mηq represents the η q effective mass [53], f ηq is the decay constant of η q and s ηq stands for the continuum threshold.By matching the hadronic expression with the OPE results with the help of the dispersion relation, one then obtains the required sum rules.Applying the Borel transformation on both sides, we can further suppress the uncertainties caused by the unwanted contributions from both the higher-order dimensional vacuum condensates and the continuum states, and our final sum rules for the moments of η q leading-twist LCDA becomes where the parameter κ = ss / q q , which comes from the use of relation g 2 s g s ψψ 2 = (2 + κ 2 ) g 2 s qq 2 with (ψ = u, d, s) in the OPE expansion.It has been shown that due to the anomalous dimension of the n th -order moment grows with the increment of n, contributions from the much higher moments at the large momentum transfer will be highly suppressed [54].Calculating the first few moments is sufficient, thus avoiding the need for further calculations.Specifically, the sum rules of the 0 th -order moment gives The effective mass mηq is taken as ∼ 370 MeV [53].To be self-consistent, we will adopt the relation ξ n 2;ηq | µ = ξ n 2;ηq | µ ξ 0 2;ηq | µ / ( ξ 0 2;ηq | µ ) 2 to calculate the n th -order moment [55].The decay constant is an important input for the B(D) → η (′) TFFs, which has been calculated under different methods such as the LCSR [56], the QCD sum rules (QCD SR) [57,58], the light-front quark model (LFQM) [59][60][61][62], the lattice QCD (LQCD) [63][64][65][66], the Bethe-Salpeter (BS) model [67][68][69], the relativistic quark model (RQM) [70][71][72], the non-relativistic quark model (NRQM) [73], and etc..As for the decay constant f ηq , those studies show f ηq is within a broader range [0.130, 0.168] GeV.At present, the sum rules of the η q decay constant can be inversely obtained by using Eq.( 14).The ξ 0 2;ηq | µ should be normalized in a suitable Borel window, which will be treated as an important criteria for determining the η q decay constant.

B. The LCHO model for ηq leading-twist LCDA
The meson's LCDA can be derived from its light-cone wave-function (LCWF) by integrating its transverse components.It is helpful to first construct the η q leadingtwist LCWF and then get its LCDA [74,75].Practically, the η q LCWF can be constructed by using the BHL prescription [38,39], and the LCHO model takes the follow-ing form [55]: where k ⊥ is the η q transverse momentum, χ 2;ηq (x, k ⊥ ) stands for the spin-space WF that comes from the Wigner-Melosh rotation and the spatial WF ψ R 2;ηq (x, k ⊥ ) comes from the approximate bound-state solution in the quark model for η q .Some more explanation on the LCWF construction can be found in Ref [55].Using the following relationship between the LCDA and LCWF, (16) and by integrating over the transverse momentum k ⊥ , one then obtains the leading-twist LCDA φ 2;ηq (x, µ), i.e.
where m q = m u = m d is the constituent light-quark mass, which is around one third of the proton mass and some typical choices are 250 MeV [76], 330 GeV [77,78] and 300 MeV [79,80], respectively.The overall parameter A 2;ηq and the transverse parameter β 2;ηq that dominates the LCWF's transverse behavior, can be fixed according to the following two constraints.One is the normalization condition, which is the same as the pionic case, e.g.
Another is the probability of finding the q q Fock state in a meson should be not larger than 1, We adopt P ηq ≈ 0.3 to carry out the following calculation, which is the same as that of pion LCWF [75].Equivalently, one can replace the constraint ( 19) by the quark transverse momentum k 2 ⊥ ηq , which is measurable and defined as [75] where the gamma function Γ[s, x] = x 0 t (s−1) e −t dt.The function ϕ 2;ηq (x) determines the dominant longitudinal behavior of φ 2;ηq , which can be expanded as a Gegenbauer series as For self-consistency, the parameters B n have been observed to closely approximate their corresponding Gegenbauer moments, i.e.B n ∼ a n , especially for the first few ones [80][81][82].The η q meson Gegenbauer moments at the scale µ can be calculated by the following way Then the Gegenbauer moments a n 2;ηq (µ) and the moments ξ n 2;ηq | µ satisfy the following relations, Using the sum rules (13) of ξ n 2;ηq | µ , one can determine the values of a n 2;ηq (µ), which will be used to fix B n .In the following we will adopt the given two Gegenbauer moments a 2,4 2;ηq to fix the parameters B 2,4 .
C. The B(D) + → ηqℓ + ν ℓ TFFs using the LCSR The LCSR approach is an effective tool in determining the non-perturbative properties of hadronic states.Here and after, we use the symbol "H" to indicate the B(D)meson for convenience.
Following the LCSR approach, one should first construct a correlator with the weak current and a current with the quantum numbers of H that are sandwiched between the vacuum and η q state.More explicitly, for H → η q , we need to calculate the correlator where the current j H = (m Q Qiγ 5 d) with Q = (b, c)quark for (B, D) meson, respectively.The LCSR calculation for the B(D) + → η q TFFs is similar to the case of B s (D s ) → η s , which has been done in Ref. [52].
In the following, we will give the main procedures for self-consistency, and the interesting reader may turn to Ref. [52] for more detail.
The dual property of the correlator ( 24) is used to connect the two different representations in different momentum transfer regions.In the time-like region, one can insert a complete set of the intermediate hadronic states in the correlator and obtain its hadronic representation by isolating out the pole term of the lowest meson state, i.e.
where the superscript "had" and "H" stand for the hadronic expression of the correlator and the continuum states of heavy meson, respectively.Here, the decay constant of B(D)-meson is defined via the equation, H| Qiγ 5 q|0 = m 2 H f H /m Q , and by using the hadronic dispersion relations in the virtuality (p + q) 2 of the current in the B(D) channel, we can relate the correlator to the H → η q matrix element [9] Due to chiral suppression, only the first term contributes to the semileptonic decay of H → η q with massless leptons in the final state.Then, the hadronic expression for the invariant amplitude can be written as where s 0 is the continuum threshold parameter, ρ H is the hadronic spectral density.
In the space-like region, the correlator can be calculated by using the operator production expansion (OPE).The OPE near the light cone x 2 ≈ 0 leads to a convolution of perturbatively calculable hard-scattering amplitudes and universal soft LCDAs.The contributions of the three-particle part being negligible [52], we solely focus on calculating the two-particle part here, and the corresponding matrix element is [83] The light-cone expansion for q 2 , (p+q) 2 ≪ m 2 b (or m 2 c ), the correlator Π OPE can be written in the general form In the above equation, the first term is the leading-order (LO) for all the LCDAs' contributions, and the second term stands for the gluon radiative corrections to the dominant leading-twist parts.
After an analytic continuation of the light-cone expansion to physical momenta using the dispersion relation, one can equate the above two representations by the assumption of quark-hadron duality.To improve the precision of the LCSR, we also apply the Borel transformation, which results in where F 0 (F 1 ) represents the leading-order (LO) or the next-to-leading order (NLO) contributions, respectively.Our final LCSR for the H → η q TFF is where The NLO invariant amplitude F 1 (q 2 , M 2 , s 0 ) can be found in Ref. [52], which is given as a factorized form of the convolutions.As will be shown below, the high-twist terms will be power suppressed and have quite small contributions to compare with those of the leading-twist terms, thus we will not discuss the uncertainties caused by the different choices of the high-twist LCDAs.For convenience, we take the η q twist-3 LCDAs φ p 3;ηq (u), φ σ 3;ηq (u), and the twist-4 LCDAs ψ 4;ηq (u), φ 4;ηq (u), together with their parameters, as those of Ref. [10].
Using the resultant B(D) → η (′) TFFs, one can further extract the CKM matrix element |V cd | or |V ub | by comparing with the predictions with the experimental data, i.e. via the following equation [84] where τ (H) is H-meson lifetime, and the maximum of squared momentum transfer q 2 max = (m H − m η (′) ) 2 .
When doing the numerical calculation, each vacuum condensates and current quark masses should be run from their initial values at an initial scale (µ 0 = 1 GeV) to the required scale by applying the renormalization group equations (RGEs) [55].
B. The ηq decay constant and the moments ξ n

2;ηq
The continuum threshold parameter (s 0 ) and the Borel parameter M 2 are two important parameters for the sum rules analysis.When calculating the decay constant f ηq , one may set its continuum threshold to be close to the squared mass of the η ′ meson, i.e. s 0 = 0.95 ± 0.1GeV 2 [57].To determine the allowable M 2 range, e.g. the Borel window, for the η q decay constant, we adopt the following criteria, • The continuum contribution is less than 30%; • The contributions of the six-dimensional condensates are no more than 5%; • The value of f ηq is stable in the Borel window; • The ξ 0 2;ηq | µ0 at the initial scale µ 0 = 1 GeV is normalized to 1, e.g.ξ 0 2;ηq | µ0 = 1.We put the decay constant f ηq versus the Borel parameter M 2 in Fig. 1, where the shaded band indicates the uncertainties from the errors of all the mentioned input parameters.The decay constant is flat in the allowable Borel window, which confirms the third criterion.Using the above four criteria and the chosen continuum threshold parameter, we put the numerical results of f ηq in Table I.As a comparison, we also present the predictions using the QCDSR and LQCD approaches.Our predictions are in good agreement with the QCDSR 2000 [57] and the LQCD 2021 predictions within errors [66].The reason why we are slightly different from QCDSR 2000 is that their calculation only includes the contributions up to five dimensional operators, and our present one includes the dimension-6 vacuum condensation terms.Using the determined f ηq , we then determine the moments of the leading-twist LCDA.Similarly, several important conditions need to be satisfied before the moments of η q LCDA can be determined [52].
The Borel window is one of the important parameters to determine the moments.When determining the Borel window, it is necessary to ensure that the contributions of both the continuum state and the dimensionsix condensates' contributions are sufficiently small.The lower limit of the Borel window is typically defined by the dimension-six condensates' contributions, while the upper limit is determined by the contribution of the continuum state.To find suitable Borel window for the moments, we adopt the dimension-six condensates' contributions to be no more than 5% and the continuum contribution to be no more than 40%.More explicitly, to fix the Borel window for the first two LCDA moments ξ n 2;ηq with n = (2, 4), we set the continuum contributions to be less than 35% and 40%, respectively.We find that the allowable Borel windows for the two moments ξ 2,4
Combining the normalization condition (18), the probability formula for q q Fock state P ηq ≈ 0.3, and the moments ξ (2,4) 2;ηq | µ0 given in Eqs.(34,35), the determined LCDA parameters are shown in Table II and their corresponding LCDA φ 2;ηq is given in Fig. 2. Its behavior of one peak with two humps is caused by a 2 2;ηq (µ 0 ) = 0.156 ± 0.042 and a 4 2;ηq (µ 0 ) = 0.055 ± 0.005, which is given by using their relations (22) to the moments ξ n 2;ηq that can be calculated by using the sum rules (13).In this paper, we take m q = 300 MeV to do the following calculation and use ∆m q = ±50 MeV to estimate its uncertainty.Table II shows that the parameters B 2 and B 4 and the quark transverse momentum k 2 ⊥ ηq increase with the increment of constituent quark mass, but the harmonious parameter β 2;ηq decreases gradually.Experimentally the average quark transverse momentum of pion, k 2 ⊥ π , is of the order (300 MeV) 2 [88].It is reasonable to require that k 2 ⊥ ηq have the value of about a few hundreds MeV [75].For the case of m q = 300 ± 50 MeV, we numerically obtain k 2 ⊥ ηq ≈ (351 +4 −3 MeV) 2 , which is reasonable and in some sense indicates the inner consistency of all the LCHO model parameters.Moreover, by using the RGE, one can get the φ 2;ηq (x, µ) at any scale µ [55].Fig. 3 shows the LCDA φ 2;ηq at several typical scales with m q = 300 MeV.At low scale, it shows double humped behavior and when the scale µ increases, the shape of φ 2;ηq becomes narrower; and when µ → ∞, it will tends to single-peak asymptotic behavior for the light mesons,
D. The TFFs and observable for the semileptonic decay B(D) One of the most important applications of the η qmeson LCDAs is the semileptonic decay H + → η (′) ℓ + ν ℓ , whose main contribution in the QF scheme comes from the |η q -component.Here H + stands for B + or D + , respectively.And to derive the required H + → η (′) TFFs, we take the mixing angle φ = (41.2+0.05 −0.06 ) • [52].The continuum threshold s H→η (′) 0 and Borel parameters M 2 are two important parameters for the LCSR of the TFFs.As usual choice of treating the heavy-tolight TFFs, we set the continuum threshold as the one near the squared mass of the first excited state of D or B-meson, accordingly.And to fix the Borel window for the TFFs, we require the contribution of the continuum states to be less than 30%.The determined values agree with Refs.[10,91], and we will take the following values to do our discussion B→η ′ = 18.0 ± 2.0 GeV.
Using Eqs.(1, 2) together with the LCSR (31) for the TFF f H→ηq + (q 2 ), we then get the results for f H→η (′) + (q 2 ), where H represents B or D, respectively.Fig. 5 shows 2 Since the leading-twist parts dominant the TFFs, this consistency also explains why our following LCSR predictions for the TFFs are close in shape with those of Ref. [10].
This work (LCSR) 0.145 LCSR 2013 [92] 0.552(51) 0.458 (105) how the total TFFs f H→η (′) + (q 2 ) change with the increment of q 2 , in which the twist-2 up to NLO QCD corrections, the twist-3 and the twist-4 contributions have been presented separately.The non-local operator matrix elements in LCSR can factorization into the universal hadron distribution amplitude, the latter term is suppressed by powers of a small parameter 1/M 2 as compared with the previous term, thereby suppressing the contribution from higher-order twist.Fig. 5 shows that the twist-2 terms dominant the TFFs.We also find that the NLO QCD corrections to the twist-2 terms are sizable and should be taken into consideration for a sound prediction.For examples, at the large recoil point, the twist-2 NLO terms give about 15.8% (17.6%) and 6.4% (7.2%) contributions to the total TFFs f D→η (′) + (0) and f B→η (′) + (0), respectively.Table III gives our present LCSR predictions for the TFFs f D→η (′) + (0) and f B→η (′) + (0).As a comparison, we have also presented the results derived from various theoretical approaches and experimental data in Table III, including the LCSR approach [9,10,92], the pQCD approach [93], the covariant light front (CLF) approach [94], the light front quark model (LFQM) approach [95], the covariant confining quark mode (CCQM) approach [96], and the BES-III Collaboration [28].The uncertainties of the + (q 2 ) with H = B + or D + in the allowable q 2 range, where the contributions from the twist-2, twist-3, twist-4 LCDAs are given separately.The twist-2 terms are given up to NLO QCD corrections.
Here the second equations show the squared averages of the errors from all the mentioned error sources.The errors are mainly caused by f ηq and m Q f H , and error caused by the mixing angle φ is quite small.As a common limitation of phenomenological quark models, accurately quantifying the theoretical uncertainty of predictions in LCSR is challenging.The net errors of these parameters are about 10% − 20%.+ (q 2 ), where all input parameters are set to be their central values.∆ is the measure of the quality of extrapolation.
f D→η The physically allowable ranges of the above four heavy-to-light TFFs are 67 GeV 2 , respectively.For light leptons, we have m 2 ℓ ≃ 0. The LCSR approach is applicable in low and intermediate q 2 region, which however can be extended to whole q 2 region via proper extrapolation approaches.In the present paper, we adopt the converging simplified series expansion (SSE) proposed in Refs.[97,98] to do the extrapolation, which suggest a simple parameterization for the heavy-to-light TFF, e.g.
where m R * = m B * = 5.325GeV (m D * = 2.010GeV) [10] are vector meson resonances, z(t, t 0 ) is a function Here is a free parameter.The free parameter b k can be fixed by requiring ∆ < 1%, where the parameter ∆ is used to measure the quality of extrapolation and it is defined as where t ∈ [0, 1  40 , • • • , 40  40 ] × 13.0(1.0)GeV for the case of η-meson, t ∈ [0, 1  40 , • • • , 40  40 ] × 11.2(0.5)GeV for the case of η ′ -meson.The two coefficients b 1,2 with all input parameters are set as their central values are listed in Table IV.The qualities of extrapolation parameter ∆ are less than ∼ 0.8%.It is noted that the fitting parameters b 1 and b 2 were obtained by rigorous fitting LCSR data.In the decay of B-meson, the difference between b 1 and b 2 is small, whereas in the decay of D-meson, there exists a significant disparity between the two.Except for the fact that the TFFs f B→η + (q 2 ) and f B→η ′ + (q 2 ) are closer in shape than the case of f D→η (′) + (q 2 ), there are other two reasons for this discrepancy, e.g.+ (q 2 ) is obtained using different parameterizations throughout the entire q 2 -region.+ (q 2 ) in whole q 2 -region, where the solid line is the central value and the shaded band shows its uncertainty.The darker part of the shaded band is the LCSR prediction, and the remaining part is the extrapolated result.As a comparison, predictions using different theoretical approaches and the experimental data, such as CCQM [96], LFQM [95], LCSR [10], pQCD [93] and BESIII collaboration [28], are also presented.
1) The whole physical region q 2 corresponding to them is quite different.The larger interval can be fitted with two similar parameters, whereas the smaller interval cannot be fitted by two similar parameters.
2) According to the SSE method employed, if two similar parameters, namely b 1 ≈ b 2 , are utilized in the decay process of D-meson, the curve fitted exhibits significant disparities compared to the outcomes derived from our LCSR calculation.And the larger the disparity between b 1 and b 2 within the same q 2 region, the greater the value of ∆, resulting in a more pronounced steepness of the fitted curve.
In addition to extrapolating with the SSE, there are other fitting methods, which have also been suggested in the literature, e.g.

1)
In early studies of beauty and charm semileptonic decays, the form of TFFs were commonly assumed to adhere to the simple pole model (also referred to as nearest pole dominance) [99].
Here m H and m H * represent the B(D)-meson mass and the corresponding vector meson resonance, respectively.This model oversimplifies the actual dy-namics, and the two free parameters (f H→η ( ′ ) + (0) and m H * ) do not fit the experimental data well when q 2 range is large.
2) Becirevic-Kaidalov (BK) parameterization [100]: where α BK is a free parameter.This parameterization of TFFs for heavy-to-light decay conforms to the heavy quark scaling law and avoids the introduction of explicit "dipole" forms, which has been used in the analyses of systematically lattice data and experimental studies of the semileptonic TFFs.
This parameterization is valid in the entire kinematical range of semileptonic decays and is consistent with vector-meson dominance at large momentum transfer.Besides, the parameters α and β are rather sensitive to the chosen range for q 2 in the actual calculation.ℓ + ν ℓ in whole q 2 -region, where the solid line is the central value and the shaded band shows its uncertainty.As a comparison, the predictions using different theoretical approaches and the experimental data, such as CCQM [96], LFQM [95], LCSR [10] and BESIII collaboration [27,28], are also presented.
4) The Boyd-Grinstein-Lebed (BGL) parametrization is based on the dispersion relation to describe the heavy mesons semileptonic decay form factors, independent of heavy quark symmetry [102].
where P (t) is Blaschke factor that depends on the masses of the sub-threshold resonances.The coefficients a k are unknown constants constrained to obey The BGL parameterization is often referred to as the series expansion (SE), whose starting point is to extend the TFFs defined in the physical range (from q 2 = 0 to (m H − m η ( ′ ) ) 2 ) to analytic functions throughout the complex t = q 2 plane.By selecting an appropriate normalization function φ(t), a simple dispersion bounds for SE coefficients can be obtained [98].
All the aforementioned parameterization methods can be utilized for the extrapolation of the TFFs, each with their own advantages and disadvantages [103,104].
The central value of the TFFs f H→η (′) + (q 2 ) are obtained through various parameterization methods is illustrated in Fig. 6.The TFFs f H→η ′ + (q 2 ) obtained by different parameterization methods exhibit consistency in across the entire q 2 region, but show discrepancy in f H→η + (q 2 ) at the larger q 2 region.Particularly for f D→η + (q 2 ), the behaviors derived from the four fitting methods differ significantly, while the trends of SSE and SE are relatively similar.The simple pole model is deemed overly simplistic, thus this portion of the graph is omitted.
Opting for SSE parameterization have a significant advantage as it effectively translates the near-threshold behavior of the TFFs into a constraining condition on the expansion coefficient.The extrapolated TFFs in whole q 2 -region are given in Fig. 7, where some typical theoretical and experimental results are presented as a comparison, such as CCQM [96], LFQM [95], LCSR 2015 [10], pQCD [93] and BESIII 2020 [28].The solid lines in Fig. 7 denote the center values of the LCSR predictions, where the shaded areas are theoretical uncertainties from all the mentioned error sources.The thicker shaded bands represent the LCSR predictions, which have been extrapolated to physically allowable q 2 -region.Fig. 7 indicates that: 1) Our present LCSR prediction of f D→η + (q 2 ) is in good agreement with BESIII data [28]; 2) Our present LCSR prediction of f D→η ′ + (q 2 ) is consistent with the LFQM prediction [95] and the LCSR 2015 [10] predictions within errors; 3) Our present LCSR predictions of f B→η (′) + (q 2 ) are close to the LCSR 2015 prediction [10], and their values at q 2 = 0 are consistent with the pQCD prediction [93] within errors.In Fig. 7, the upward trend for the f D→η + (q 2 ) TFF exhibits a significant difference in the large q 2 region with other theoretical results and experimental data, which there are two reasons for this difference, the main reason is because we get a smaller f D→η + (q 2 ) TFF at q 2 = 0.In order to obtain a smooth fitting curve, there is an overall upward trend in the large q 2 region.And the secondary reason is the utilization of distinct methods for extrapolation.The other theoretical groups in Fig 7 employed the double-pole fitting method to obtain the form factor across the entire physical region, whereas we opted for the simpler SSE approach.Fig. 8 presents the differential decay widthes for B(D) + → η (′) ℓ + ν ℓ without CKM matrix elements.As a comparison, the predictions using different theoretical approaches and the experimental data, such as CCQM [96], LFQM [95], LCSR [10] and BESIII collaboration [27,28], are also presented.The differential decay width dΓ/|V cd |dq 2 (D + → ηℓ + ν ℓ ) agrees with the BESIII 2018 [27] and BESIII 2020 [28] within errors.
By matching the branching fractions and the decay lifetimes given by the PDG with the decay widthes predicted by Eq.( 32), one may derive the CKM matrix elements |V ub | and |V cd |.We put our results in Table V, where the errors are caused by all the mentioned error sources and the PDG errors for the branching fractions and the decay lifetimes.Some typical measured values of |V ub | and |V cd | are also given in Table V.The predicted |V cd | is within the error range of experimental result BESIII 2020.Using the fixed CKM matrix ele-

IV. SUMMARY
In this paper, we have suggested the LCHO model (17) for the η q -meson leading-twist LCDA φ 2;ηq (x, µ), whose moments have been calculated by using the QCD sum rules based on the QCD background field.To compare with the conventional Gegenbauer expansion for the LCDA, the LCHO model usually has better end-point behavior due to the BHL-prescription, which will be helpful to suppress the end-point singularity for the heavy-tolight meson decays.The QCD sum rules for the 0 th -order moment has been used to fix the η q decay constant, and we obtain f ηq = 0.141 ± 0.005 GeV, which is slightly larger than the conventional value of the pion decay constant f π ≃ 0.130 GeV [24].As an explicit application of φ 2;ηq , we have calculated the TFFs B(D) + → η (′) under the QF scheme for the η − η ′ mixing and by using the QCD light-cone sum rules up to twist-4 accuracy and by including the next-to-leading order QCD corrections to the dominant leading-twist part.Our LCSR prediction of TFFs are consistent with most of theoretical predictions and the recent BESIII data within errors.By applying those TFFs, we get the decay widths of B(D) + → η (′) ℓ + ν ℓ .The magnitudes of the CKM matrix elements |V ub | and |V cd | have also been discussed by inversely using the PDG values for the branching fractions and the decay lifetimes.The future more precise data at the high luminosity Belle II experiment [106] and super tau-charm factory [107] shall be helpful to test all those results.

20 FIG. 1 :
FIG.1:(Color online) The ηq decay constant fη q versus the Borel parameter M 2 , where the shaded band indicates the uncertainties from the input parameters.

TABLE I :
The decay constant fη q using the BFTSR approach.As a comparison, typical results derived from QCDSR and LQCD approaches have also been presented.

TABLE III :
Typical theoretical predictions on the TFFs f H→η(′)