Local-duality QCD sum rules for strong isospin breaking in the decay constants of heavy–light mesons

We discuss the leptonic decay constants of heavy–light mesons by means of Borel QCD sum rules in the local-duality (LD) limit of infinitely large Borel mass parameter. In this limit, for an appropriate choice of the invariant structures in the QCD correlation functions, all vacuum-condensate contributions vanish and all nonperturbative effects are contained in only one quantity, the effective threshold. We study properties of the LD effective thresholds in the limits of large heavy-quark mass mQ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_Q$$\end{document} and small light-quark mass mq\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_q$$\end{document}. In the heavy-quark limit, we clarify the role played by the radiative corrections in the effective threshold for reproducing the pQCD expansion of the decay constants of pseudoscalar and vector mesons. We show that the dependence of the meson decay constants on mq\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_q$$\end{document} arises predominantly (at the level of 70–80%) from the calculable mq\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_q$$\end{document}-dependence of the perturbative spectral densities. Making use of the lattice QCD results for the decay constants of nonstrange and strange pseudoscalar and vector heavy mesons, we obtain solid predictions for the decay constants of heavy–light mesons as functions of mq\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_q$$\end{document} in the range from a few to 100 MeV and evaluate the corresponding strong isospin-breaking effects: fD+-fD0=(0.96±0.09)MeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{D^+} - f_{D^0}=(0.96 \pm 0.09) \ \mathrm{MeV}$$\end{document}, fD∗+-fD∗0=(1.18±0.35)MeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{D^{*+}} - f_{D^{*0}}= (1.18 \pm 0.35) \ \mathrm{MeV}$$\end{document}, fB0-fB+=(1.01±0.10)MeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{B^0} - f_{B^+}=(1.01 \pm 0.10) \ \mathrm{MeV}$$\end{document}, fB∗0-fB∗+=(0.89±0.30)MeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$f_{B^{*0}} - f_{B^{*+}}=(0.89 \pm 0.30) \ \mathrm{MeV}$$\end{document}.


Introduction
The method of QCD sum rules [1], based on the exploitation of Wilson's operator product expansion (OPE) in the study of properties of individual hadrons, has been extena e-mail: simula@roma3.infn.it sively applied to the decay constants of heavy mesons [2][3][4]. An important finding of these analyses was the observation of the strong sensitivity of the decay constants to the precise values of the input OPE parameters and to the algorithm used for fixing the effective threshold [5][6][7][8][9]. For any given approximation of the hadronic spectral density based on quark-hadron duality, the effective threshold determines to a large extent the numerical prediction for the decay constants inferred from QCD sum rules: even if the parameters of the truncated OPE are known with high precision, the decay constants may be predicted with only a limited accuracy, which represents their systematic uncertainty. In a series of papers [10][11][12][13][14], we proposed a new algorithm for fixing the effective threshold within the Borel QCD sum rules which allowed us to obtain realistic estimates of the systematic uncertainties. Our procedure opened the possibility to get predictions for the decay constants with a controlled accuracy [15][16][17][18] and thus allowed us to address subtle effects that call for a profound accurate treatment, such as the ratios of the decay constants of heavy vector and pseudoscalar mesons [19,20] or the strong isospin-breaking (IB) effects in the decay constants of heavy-light mesons [21], generated by the mass difference (m d − m u ) between up and down quarks.
Here, we discuss the application of another variant of QCD sum rules to the evaluation of the strong IB effects in the decay constants of heavy-light pseudoscalar and vector mesons. Our analysis takes advantage of the fact that the OPE provides the analytic dependence of the correlation functions on the quark masses; this allows us to study, e.g., the impact of the light-quark mass on heavy-meson decay constants, thus providing access to the strong IB effects. The approach we describe in this work seems quite promising for studying the dependence of a generic hadron observable on quark masses. 1.1 QCD sum rule in the local-duality limit A typical Borel QCD sum rule for the decay constant f H of a heavy (pseudoscalar or vector)Qq meson H of mass M H , consisting of a heavy quark Q with mass m Q and a light quark q with mass m q , has the form Here, τ is the Borel parameter, N is an integer that depends on the Lorentz structure in the correlation function chosen for the sum rule and on the number of subtractions in the corresponding dispersion representation, and s eff is the effective threshold such that √ s eff lies between the mass of the groundstate and the first excited state [1], namely s eff = (M H +z eff ) 2 with z eff 0.4-0.5 GeV.
Nonperturbative effects appear on the r.h.s. of (1) at two places: as power corrections given in terms of vacuum condensates and in the effective threshold s (N ) eff . Depending on the chosen value of N , nonperturbative effects are distributed in a different way between power corrections and the effective threshold. Perturbative effects are encoded in the spectral density ρ pert , in the effective threshold s (N ) eff and in the power corrections Π (N ) power . Recall that Eq. (1) is based on modelling the hadron continuum as the effective continuum, i.e., on the substitution ρ cont (s) = θ(s − s eff )ρ pert (s). This relation is fulfilled pointwise at large values of s above some s up , but is a "weak" relation and requires an appropriate smearing for s in the midenergy region above the physical hadron continuum threshold s th . An appropriate smearing is reached by performing the Borel transform For nonzero τ , the contribution of the hadron continuum given by (2) is exponentially suppressed compared to the ground-state contribution (1). Therefore, in the conventional use of QCD sum rules one works in some window of nonzero values of τ . However, one may ask whether or not Eq. (1) may be extended down to τ = 0, the so-called local-duality (LD) limit. 1 Obviously, at τ = 0 an appropriate smearing in (2) is guaranteed by the integration; on the other hand, the excited states are not suppressed and one can doubt that 1 The LD limit in Borel sum rules was introduced and discussed in [22][23][24][25] in connection with the pion and nucleon elastic form factors, and later applied to the analysis of meson transition form factors in [26][27][28][29]. A specific feature of this limit is the vanishing of power corrections in the two-and three-point Borelized correlation functions of axial-vector and vector currents of light quarks. for N = 1 [32][33][34] modelling the hadron continuum as the effective continuum remains a good approximation at small τ . First, note that power corrections contain singular terms of the form τ 2−N log(τ ). Therefore, the limit τ → 0 cannot easily be taken in the sum rule (1) for N ≥ 2. For N = 0 and N = 1, the limit τ → 0 in (1) is mathematically well defined. To demonstrate that this limit is also physically meaningful, one needs to show that the corresponding s eff indeed lies in the expected range. Figure 1a presents the effective threshold obtained in the case of the vector B * -meson by solving Eq. (1) for N = 0 and N = 1, using in its l.h.s. the results of recent lattice QCD simulations for f B * [30] and the experimental value for M B * [31]. Figure 1b shows the truncated series of power corrections including operators of dimension up to 6 for both N = 0 and N = 1. Note that power corrections at τ = 0 vanish for N = 0 and take a finite value for N = 1. The vanishing of N = 0 power corrections at τ = 0 is related to the absence in QCD of a dimension-2 condensate. Obviously, the truncated power corrections for N = 1 remain under control in a rather broad range of τ , but for N = 0 explode relatively soon as τ increases.
It is clear that for the N = 1 case the OPE is under good control in a broad range of τ and therefore the lower boundary of the Borel window can be safely extended down to τ = 0. For N = 0 only at relatively small τ the OPE is under control and the approximation of a τ -independent effective threshold may work well. The relevance of the unknown higher-order power corrections is reflected by the sharp rise of the effective threshold visible in Fig. 1a. Obviously, the standard QCD sum-rule analysis based on a stability window with a constant effective threshold may be problematic. In this respect an alternative approach based on a τ -dependent effective threshold seems to be more appropriate, but this issue goes well beyond the scope of the present paper. Here, the only important property is that the effective threshold at τ = 0 has the value expected on the basis of the standard considerations [1], i.e., it is around (M B * + z eff ) 2 with z eff ∼ 0.4-0.5 GeV. Therefore, for the correlators with N < 2, modelling the hadron continuum as an effective continuum remains a valid and equally accurate approximation as τ → 0, which does not represent a point of discontinuity.
In this work we show that the sum rule (1) for N = 0 can be of particular interest. Obviously, considering the sum rule at only one point, τ = 0, does not allow for the use of the usual sum-rule stability criteria [1] for determining s eff . Consequently, the decay constants cannot be determined entirely within the QCD sum rules; some "external" inputs are needed to determine s eff . Nevertheless, in this work we will show that, besides the reduction of the uncertainties related to the absence of the condensates, the LD sum rules represent an efficient tool to investigate the dependence of the pseudoscalar and vector meson decay constants on quark masses and their perturbative behaviour in QCD. Moreover, the LD sum rules turn out to be particularly suitable for the analysis of the strong IB effects in the two-point functions when implemented with only few "external" inputs, e.g., from experiment or lattice QCD.
1.2 Strong isospin breaking from a QCD sum rule in the LD limit We are interested in the dependence of the decay constants of heavy-light mesons on the quark masses, in particular, in the strong IB effects in the decay constants (i.e., the difference between the decay constants ofQd andQu mesons induced by the small mass difference δm = m d − m u ). We therefore need to properly take into account all effects depending on the light-quark flavour q in the correlation function of the appropriateQq interpolating currents. Clearly, the m q -dependence on the l.h.s. of (1) is encoded both in the decay constant f H and in the meson mass M H . On the r.h.s, the IB effects come from several sources: the m qdependence of ρ pert (s, m Q , m q , α s ), the m q -dependence of the effective threshold s eff , the m q -dependence of the power corrections, and the flavour dependence of the quark condensates, in particular, of qq . In general, all these effects mix together, which renders the goal of isolating the IB effects in f H a complicated task. A careful analysis has been carried out recently in [21], following the standard choice N = 2 for pseudoscalar and N = 1 for vector mesons.
There is, however, a special case which makes the sum rule (1) particularly suitable for the analysis of IB effects. As discussed above, for N = 0 and N = 1 power corrections are regular functions at τ = 0. Moreover, for N = 0 power corrections at τ = 0 vanish (power corrections for N = 1 are nonzero at τ = 0) and Eq. (1) is reduced to On the l.h.s., the M H contribution has dropped out, thus opening a direct access to the m q dependence of the decay constant f H . Since power corrections do not contribute to the sum rule in the LD limit, all nonperturbative effects enter now through a single quantity -the effective threshold. The functional dependence of the perturbative spectral density on the quark masses m Q and m q can be calculated to the necessary accuracy. The functional dependence of s eff on the quark masses may be determined from the general properties of the decay constants of heavy-light mesons in QCD. Namely, its dependence on m q may be parameterized by a polynomial formula in m q plus chiral logs, which can be determined by matching to heavy-meson chiral perturbation theory. The numerical coefficients in the polynomial function are not known but may be determined using only few results on the decay constants from lattice QCD, e.g., for nonstrange and strange heavy mesons. Having at our disposal the explicit m q dependence of the effective threshold and of the spectral densities opens direct access to the strong IB effects related to the small difference of the u-and d-quark masses in QCD. We will demonstrate that the main m q dependence of the decay constants originates from the calculable m q dependence of the perturbative spectral densities. Therefore, the LD limit opens the possibility of a reliable analysis of the m q dependence and the strong IB effects in the decay constants of heavy-light mesons (and, in principle, also in other quantities).
This paper is organized as follows: in Sect. 2, we recall the spectral densities of the QCD correlation functions relevant for our LD sum-rule analysis. In Sect. 3, we study the m Q -and m q -dependences of the effective thresholds by making use of an appropriate mass scheme (pole mass and running mass) for the heavy quarks. In Sect. 4, we perform the numerical analysis of the decay constants of heavy-light pseudoscalar and vector mesons and obtain predictions for strong IB effects in the decay constants. Section 5 gives our conclusions. The Appendix collects some details of treating the IB effects within the OPE, which, in our opinion, deserve to be presented.

Local-duality sum rules for f P and f V
Let us consider two-point QCD sum-rules for decay constants of pseudoscalar (P) and vector (V ) mesons built up of one massive quark Q with mass m Q and one light quark q with mass m q . We consider the axial-vector current and the vector current as interpolating currents for the pseudoscalar and vector mesons, respectively. The corresponding correlation functions involve two Lorentz structures, the transverse structure g μν p 2 − p μ p ν and the longitudinal structure p μ p ν : For Π 5 μν ( p), we study the longitudinal structure p μ p ν , as it contains the ground-state pseudoscalar-meson contribution with For Π μν ( p), we study the transverse structure g μν p 2 − p μ p ν , which contains the ground-state vector-meson contribution with As already noted, power corrections for dimension-2 correlation functions vanish in the LD limit τ = 0. We now present their explicit form. The leading power correction to Π 5 μν ( p) is given by the quark condensate and easily derived: The Borel transform By changing the sign of the light-quark mass, the power corrections for the vector correlator Obviously, the Borelized power corrections to both the pseudoscalar and the vector correlators vanish in the limit τ = 0.
Since the power corrections do not contribute to the LD sum rule under consideration, we need to consider only the perturbative contributions. After applying the duality cuts at s eff , separately in the pseudoscalar and the vector channels, performing the Borel transform p 2 → τ , and setting τ → 0, the corresponding sum rules take the form The functions ρ pert P and ρ pert V in (14) are the spectral densities of the invariant functions Π 5 L ( p 2 ) and Π T ( p 2 ), respectively. Let us emphasize that in (14) both the full spectral densities and the decay constants are scale independent quantities. Therefore, the effective thresholds are scale independent objects, too. In perturbation theory, the spectral densities are calculated as power expansions in a ≡ α s (μ)/π , α s (μ) the strong coupling in the MS-scheme at scale μ: with i = P, V . In practice, one adopts truncated expansions of the spectral densities; this leads to a scale dependence of the spectral densities. As the result, the effective thresholds will also depend on the scale, to compensate the scale dependence of the spectral densities emerging in the course of trun-cation. Explicitly, the leading-order (LO) spectral densities read Obviously, the lower integration limit in (14) is determined by the threshold in the correlation functions. In Eq. (15), we may employ different definitions of the quark masses: The most advanced calculation of the pseudoscalar and vector spectral densities including order-O(a 2 ) terms was performed [35,36], for a massless light quark, in terms of the heavy-quark pole mass. The expansion in terms of the heavy-quark pole mass is appropriate for considering the heavy-quark limit, which we address in Sect. 3.1.
However, the pole-mass expansion leads to a rather slow convergence of the perturbative expansion for the decay constants [15][16][17][18]37]. The convergence improves considerably when one rearranges the perturbative expansion in terms of the running MS masses. Therefore, for the practical analysis of the m q -dependences of the meson decay constants in Sect. 4, we make use of the perturbative expansion in terms of the running MS masses of the light and the heavy quarks. The corresponding NLO and NNLO functions ρ (i) (15) necessary for such an analysis are found from the spectral densities of the pseudoscalar correlation function given in [37] by multiplying them by 1/s 2 . Similarly, the transverse spectral densities ρ (15) are found from the spectral densities of [32][33][34] by multiplying them by 1/s. In our analysis, we make use of the exact LO perturbative spectral density given by (16), at the NLO we keep the terms O(am 0 q ) and O(am 1 q ), and in the NNLO we keep only the known terms of order O(a 2 m 0 q ). We would like to emphasize that the perturbative spectral density (15) does not generate terms of order m q log(m q ) in the dual correlator (14). This observation will be crucial for discussing properties of the effective thresholds in the next section.

Dependence of the effective thresholds on the quark masses
Let us now consider the dependences of the effective threshold on the quark masses m Q and m q .
3.1 Heavy-quark limit in the pole-mass scheme We start with the heavy-quark limit of the decay constants, originally discussed in Refs. [38,39] within the Heavy Quark Effective Theory (HQET). In what follows, however, we do not consider the static decay constants and we work in full QCD.
For the sake of argument we consider first a massless light quark: m q = 0. We can make use of any scheme for the heavy-quark mass, but we start with the pole-mass scheme, which leads to a more transparent behaviour of the effective threshold. We first isolate the pole mass, which we denote M Q , in the effective threshold: For the decay constant of pseudoscalar and vectorQq mesons, using results from [35,36], we obtain in the limit Hereafter, we denoteā ≡ᾱ s (M Q )/π ,ᾱ s (M Q ) the running strong coupling in the MS scheme at scale M Q , and we use the standard notations C F = (N 2 c − 1)/(2N c ), C A = N c , T = 1/2, and n l the number of massless quarks [35,36].
Since only the near-threshold behaviour of the spectral densities is relevant for the leading behaviour in the large-M Q limit, we may obtain also the O(ā 2 ) terms in the dual correlation functions [i.e., the r.h.s. of (19)] from the analytical expressions for these spectral densities given by Eqs. (30) and (31) of [35,36].
In the limit M Q → ∞, the dual correlation functions, expressed in terms of z pole eff (M Q ), do not contain corrections of orderā n M Q (this property will not hold in the runningmass scheme) but still contain log(M Q ) terms of the type (ā log(M Q /z pole eff )) n ,ā(ā log(M Q /z pole eff )) n−1 , etc. The terms (ā log(M Q /z pole eff )) n , although formally of orderā n , remain unsuppressed in the limit M Q → ∞. To treat all terms containing log(M Q ), it is important to emphasize that they are exactly the same in the vector and the pseudoscalar sum rules. Therefore, they may be resummed by introducing a properly defined effective threshold z HQ eff , one and the same in the pseudoscalar and the vector channels. The explicit relation between z pole eff (M Q ) and z HQ eff , including also the a 2 terms, reads The new quantity z HQ eff , which has the meaning of the effective threshold in HQET, absorbs all log(M Q ) terms on the r.h.s. of the sum rules (19); the latter assume a form in which the HQ limit may easily be taken: We did also calculate the O(ā 2 ) contributions but do not present their explicit expressions here. The expressions (21) immediately lead to the ratio of the decay constants in the heavy-quark limit [40,41]. Including also O(ā 2 ) corrections, we obtain 2 with ζ 3 1.202. The O(ā 2 ) term in (22) reproduces the result first presented in Eq (3.12) of [41].
Notice that for finite M Q , z pole eff contains not only the logarithmic corrections, which are the same in the pseudoscalar and the vector channels, but also the 1/M Q corrections, 2 The second-order pseudoscalar and vector spectral densities near the threshold, Eqs. (30) and (31) in [35,36], contain three unknown constants,c F F ,c F A , andc F L , which cancel in the ratio f V / f P .
which are different for the thresholds in the pseudoscalar and the vector sum rules.

Combined heavy-quark and chiral limits in the pole-mass scheme
The results (19) are obtained for a massless light quark.
Switching on a small light-quark mass m q , the leading corrections generated by integration of the perturbative spectral densities are proportional to m q . As already noted in [21], no chiral logs of the kind m q log(m q ) arise from integrating the spectral densities. Therefore, chiral logs in the decay constants may be generated only by chiral logs in the effective threshold. Moreover, in order to study the chiral logs in the decay constants, it is sufficient to make use of the perturbative spectral densities for m q = 0. On the other hand, heavy-meson chiral perturbation theory (HMChPT) (S. R. Sharpe, private communication; see the appendix in [21]) [42], requires the appearance of chiral logs, which we denote as z HQ L in the chiral expansion of the decay constants in the heavy-quark limit. Since the only source of such terms is the effective threshold, we write where the dots denote linear and higher-order terms in the light-quark mass m q . The coefficient z HQ L can now be fixed by matching to HMChPT (S. R. Sharpe, private communication; see the appendix in [21]) [42], which provides the explicit chiral logs R χ (m q ) in the ratio f H q (m q )/ f H ud , with H ud a heavy meson with a light valence quark of the average mass m ud ≡ (m u + m d )/2. Finally, we obtain The explicit expression for R χ (m q ) was derived in (S. R. Sharpe, private communication; see the appendix in [21]) [42] and presented by Eq. (A.3) of [21].

Quark-mass dependences of the effective threshold in the running-mass scheme
For practical sum-rule analyses of decay constants, one prefers the MS running-mass scheme since it entails a better convergence of the perturbative expansion [15][16][17][18]37]. It is not difficult to perform the limit m Q (μ) → ∞ also for the running-mass correlation function. Also therein one can write The effective thresholdz eff (μ) in the MS scheme is related to z pole eff introduced in the pole-mass scheme through an obvious relation which just expresses the fact that the upper integration limit s eff is a scheme-independent quantity: In particular, for μ = m Q , taking into account that one finds Since z pole eff does not contain terms scaling as M Q in the limit M Q → ∞,z eff (μ) should contain terms which diverge as powers of a n M Q in this limit. This is, of course, no obstacle for usingz eff (μ) in the analysis of the decay constants of charmed or beauty mesons but makes this quantity not particularly convenient for studying the heavy-quark limit of the sum rules. The terms inz eff (μ) divergent as m Q → ∞, however, do not lead to divergent terms in the decay constants; also, the behaviour of the spectral densities in the MS scheme is a bit more tricky than in the pole-mass scheme. The dual correlator is determined by the end-point behaviour of the spectral densities; as already mentioned in [37], the higherorder spectral densities in the MS scheme do not vanish at the threshold but behave as Finally, when the MS spectral densities are used and the duality cut is expressed via z HQ eff , all terms containing powers of m Q -those coming from the integrals of the spectral densities and those contained in z eff (m Q ) -cancel each other, yielding a sum rule for f 2 H which can also be obtained just by expressing M Q via m Q in (19), e.g., Let us now switch on a small light-quark mass m q . The spectral densities are now treated as functions of m Q (μ) and m q (μ). Taking into account that the effective threshold depends on the scale μ only because of the truncation of the perturbative series, and that the chiral logs have been fixed in the pole-mass scheme, it is convenient to work with the following parameterisation for s eff : The pole mass M (2) Q here is understood as being expressed via the running mass m Q (μ) (e.g., [43]) at O(a 2 ) accuracy, the available accuracy of the correlation function. We can rewrite this expression in a form similar to (26) Let us recall that the chiral logs z HQ L have been calculated in the heavy-quark limit; at finite values of m Q , chiral logs receive corrections which are unknown. So we take into account only the known leading effect of chiral logs, to study whether or not their impact on the IB is crucial. Two other parameters of the effective thresholdz pole eff andz 1 (μ) if one makes use of the parameterisation (31), orz 0 (μ) andz 1 (μ) if one works with (32) -are unknown and will be fixed by using some external benchmark results for the decay constants from lattice QCD. The inclusion of higher-order terms in the light-quark mass has no impact on the decay constants; thus, such terms are not considered.

Numerical analysis of the sum rules
Now, we turn to the numerical estimates. For the relevant OPE parameters, we adopt the following numerical input: We work with the effective threshold in the form (32) and consider the following three Ansätze: 1. "Constant" threshold: thez 1 (μ) term in the effective threshold (32) and the chiral logs z HQ L are neglected; the only unknown parameterz 0 (μ) is fixed from the lattice results for the decay constants of the isospin-symmetric heavy mesons, with m q = m ud . Table 1 Parameters of the effective thresholds and resulting IB in the decay constants of heavy pseudoscalar and vector mesons. The parameter z L in the effective threshold for the "linear + log" ansatz is fixed by ChHQET in the heavy-quark limit

Meson
Threshold 2. "Linear" threshold: the chiral logs z HQ L are neglected and the parametersz 0 (μ) andz 1 (μ) are fixed by the lattice QCD results for the decay constants at two m q values, for the isospin-symmetric and the strange heavy mesons. 3. "Linear + log" threshold: the known leading chiral logs represented by z HQ L are included; the parametersz 0 (μ) andz 1 (μ) are fixed from the lattice QCD results for the decay constants at two m q values, for isospin-symmetric and strange heavy mesons.
As we have already noted, because of the truncation of the perturbative expansion, the truncated spectral densities depend on the scale μ. Obviously, the parametersz 0 andz 1 are also μ-dependent.
For fixing the parameters of the effective thresholds, we make use of the following results from lattice QCD: f D s f D = 1.1716 ± 0.0032 [24], In these formulae, f H denotes the decay constant of the isospin-averaged heavy-light mesons with the light-quark mass m ud , whereas f H s denotes the decay constant of the heavy strange mesons. Table 1 summarizes the effective thresholds corresponding to our three Ansätze and presents estimates of the strong IB effect. For our final estimates, we perform a bootstrap analysis of the uncertainties assuming that the OPE parameters in (33) have Gaussian distributions with corresponding Gaussian errors, whereas the scale μ has a flat distribu-tion in the range 1 < μ (GeV) < 3 for charmed mesons and 3 < μ (GeV) < 6 for beauty mesons.
As soon as the effective thresholds are known, we readily get the decay constants f H q as a function of the scale independent ratio (m q − m ud )/(m s − m ud ). The results for the ratios of the decay constants f H q / f H ud are shown in Figs. 2 and 3.
Notice that the results corresponding to a constant effective threshold [ansatz (1)] are quite close to those obtained including the m q -dependence [ansatz (2)] and to the results of Ref. [21], which contain effects in the decay constants at any order in the light-quark mass. So, an important conclusion to be drawn from our results is that effects at order O(m 2 q ) in the effective threshold are not crucial for describing the m q -dependence of the decay constants and for estimating the slope of the IB effect at the physical value of the light-quark mass: the latter are both determined to a large extent by the known m q -dependence of the spectral densities and can thus be reliably controlled in our approach.

Summary and conclusions
We addressed the local-duality (LD) limit, τ = 0, of the Borel QCD sum rules for the decay constants of heavy-light pseudoscalar and vector mesons. An invaluable feature of the LD limit is that for a proper choice of the correlation function, all vacuum-condensate contributions vanish and the full nonperturbative QCD dynamics is parameterized in terms of merely one quantity -the effective threshold. Our analysis demonstrates that the effective threshold has a nontrivial functional dependence on the masses of the heavy and the light quarks, m Q and m q , respectively. This dependence has been parameterized in the form suggested by the behaviour of the decay constants in the known limits: the chiral limit for m q and the heavy-quark limit for m Q . In the heavy-quark  Table 1 are displayed. We also show results from an alternative analysis based on Borel QCD sum rules [21] (a) (b) Fig. 3 The same as in Fig. 2 but for pseudoscalarbq (a) and vectorbq (b) mesons limit, we clarify the role played by the radiative corrections in the effective threshold for reproducing the pQCD expansion of the decay constants of pseudoscalar and vector mesons.
This paper elucidates the dependence of the decay constants on a light-quark mass m q in the range m ud < m q < m s . Fixing a few numerical parameters of the effective threshold by using the available accurate inputs from lattice QCD, we have derived the full analytic dependence of the decay constants f H (m q ) on the light-quark mass m q . The resulting dependence of the decay constants f H (m q ) on m q emerges from two sources: (i) from the m q -dependence of the QCD perturbative spectral densities known explicitly as expansion in powers of α s and (ii) from the m q -dependence of the effective threshold known approximately. An important outcome of our analysis is that the variation of the decay constants with respect to m q comes to a great extent (70-80% of the full effect) comes from the rigorously calculable dependence on m q of the perturbative spectral densities and is therefore under a good theoretical control.
Noteworthy, the known perturbative expansion of the correlation functions [32][33][34][35][36][37], where the sea-quark mass effects are neglected, limits the accuracy of the decay constants of the heavy-light mesons to O(m sā 2 ) accuracy,ā ∼ 0.1 at the appropriate renormalisation scales. Therefore the accuracy of the individual decay constants obtained from QCD sum rules does not exceed a few MeV. Nevertheless, we would like to emphasize that the IB difference of the decay constants, f M d − f M u , where the sea-quark contributions of order O(m s,u,dā 2 ) cancel each other, may be predicted with a much higher accuracy, O(δmā 2 ). Therefore, the proposed method can potentially provide a higher accuracy of the IB effects than other approaches.
As our final estimates of the IB, we take the average of the results corresponding to the linear and the linear + log effective thresholds in Table 1: Sizeably larger uncertainties of the IB in the decay constants of vector mesons compared to pseudoscalar mesons are related to larger uncertainties of the input lattice QCD results for the corresponding ratios f H s / f H ud . These estimates are in good agreement with the results of our recent analysis within a different version of QCD sum rules -the Borel sum rules with τ -dependent threshold [21]: The only exception is the D * case, where one observes tension between these two sets of the results; note, however, that also the uncertainties of these predictions are rather large.
Very recently [48] a new precise determination of the strong IB effect in the decay constants of D-and B-mesons has been carried out by the FNAL and MILC lattice collaborations.
In the charm sector their result is f D + − f D 0 = 1.13(15) MeV, which nicely agrees with our findings (35) and (39). As for the bottom sector, it is shown that the available HPQCD and RBC/UKQCD calculations [49][50][51] overestimate significantly the strong IB effect because of an inappropriate use of unitary lattice points (i.e. those having the same mass for valence and sea light-quarks). The FNAL/MILC result is f B 0 − f B + = 1.12 (15) MeV, which is in excellent agreement with our findings (37) and (41).
Thus, our sum-rule predictions are nicely confirmed quantitatively by lattice QCD both for the central values and the overall uncertainties. This is reassuring: the strong IB effect and its uncertainty in the decay constants of heavy-light mesons can be reliable and accurately estimated within the QCD sum-rule approach.
It should be emphasized that the present approach based on the combination of OPE and a few inputs from lattice QCD potentially has fewer theoretical uncertainties than other formulations of QCD sum rules: first, the condensate contributions, in particular, those of the quark condensate, which produced the main OPE error in the decay constants, vanish in the LD limit; second, the systematic uncertainty of the sum-rule method is now encoded in only one quantity -the effective threshold, which may be fixed to good accuracy due to the use of the few accurate lattice inputs.
Thus, QCD sum rules for the mass dimension-2 Borelized invariant amplitudes at τ = 0 (i.e., an infinitely large Borel mass parameter) provide an efficient tool for the analysis of the dependence of decay constants (and potentially of other hadron observables) on quark masses.
Finally, we want to mention that, besides the strong IB effect due to the up and down quark mass difference, there are other isospin violating effects due to electromagnetism, i.e. to the difference between the up and down quark electric charges. However, the inclusion of such electromagnetic corrections within a sum-rule approach is not a trivial task and it requires the development of new strategies going beyond the traditional QCD sum-rule approaches. In this respect it is worth mentioning a new lattice strategy [52] developed to deal with QCD + QED effects on quantities that require the cancellation of infrared divergences in the intermediate steps of the calculation, like, e.g., the decay rate of charged pseudoscalar mesons [53].
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .

Appendix: Isospin breaking in the OPE
The two-point Green function Π of interest is given by the functional integral where j 1 and j 2 are (gauge-invariant) operators constructed from quark and gluon fields, and Here, L (0) (x) is the SU (2)-symmetric Lagrangian describing two equal-mass quarks, with the quark-mass term m(dd +ūu), m ≡ 1 2 (m d + m u ).
After expanding Eq. (43) in powers of δm, one finds with the superscript "0" indicating that the full Green functions correspond to the SU (2)-symmetric QCD with two light quarks with degenerate mass m. Let us emphasize the appealing feature of the expansion (47): at each order in δm, one encounters the full Green function of the SU (2)symmetric QCD. One may expect that the power corrections in the OPE for the three-point Green function Γ of the SU (2)-symmetric QCD are the SU (2)-symmetric condensates, e.g., ūu = d d ≡ qq . However, in order to obtain the contributions of interest, we need to perform the limit q → 0. This step cannot be done easily: the OPE for Γ is only given by local SU (2)symmetric condensates if one keeps q 2 large and negative; a straightforward extension of the known power corrections to q → 0 leads to a wrong result: it is well known that if one naively extends power corrections in the vector three-point function to q = 0, then they do not satisfy the Ward identity (see, e.g., [26][27][28][29]).
On the other hand, one can proceed by expanding Π( p) in powers of the small quark mass; then the mass derivatives emerge. Translating the expression (47) into momentum space, to O(δm) accuracy we obtain where Π (0) ( p) is the full two-point function of SU (2)symmetric QCD, and Γ (0) ( p, q = 0) is the three-point function of the scalar currentqq at zero momentum transfer, also calculated in the full SU (2)-symmetric theory. Consequently, finding the leading-order SU (2)-breaking effects reduces to calculating the Green functions in SU (2) Using the well-known relation the three-point function at zero momentum transfer may be related to the mass derivative of the two-point function, which then leads to the appearance of the mass derivatives of the quark condensate.