# Mass Spectra and Decay Constants of Heavy-Light Mesons: A Case Study of QCD Sum Rules and Quark Model

Authors: Halil Mutuk
Journal: Advances in High Energy Physics
 (2018)

---

## Abstract

We visited mass spectra and decay constants of pseudoscalar and vector heavy-light mesons (B, Bs, D, and DS) in the framework of QCD sum rule and quark model. The harmonic oscillator wave function was used in quark model while a simple interpolating current was used in QCD sum rule calculation. We obtained good results in accordance with the available experimental data and theoretical studies.


---

## Body

## 1. Introduction

The ultimate objective of particle physics is to investigate and examine the structure and the origin of matter. For this purpose many theoretical and experimental endeavors are made, and a resulting model was theorized, which we call the Standard Model of particle physics. Quark model which was proposed by Gell-Mann and Zweig in 1964 [1] is a part of the Standard Model, and interprets hadrons fairly compatible with the experimental data. According to the quark model, mesons are made of quark-antiquark pairs (qq¯) and baryons are made of three quarks qqq or antiquarks q¯q¯q¯. These quarks interact with each other via emitting and/or absorbing gluons. The resulting theory which explains these interactions is the Quantum Chromodynamics (QCD).The interaction of quarks is described by QCD, which is part of the Standard Model of particle physics. QCD is thought to be thetrue theory of strong interactions. QCD is a SU(3) gauge theory describing the interactions of six quarks which transform under the fundamental representation of SU(3) group via the exchange of gluons that transform under the adjoint representation. Although it has been more than 50 years that QCD has been proposed, a solution has been evaded. Contrary to electroweak theory, where it is possible to obtain precise results using perturbation theory, the order of precision obtained in QCD has been lower by orders of magnitude. The main reason for this is that the coupling constant (which should be the perturbation parameter) of QCD is of the order one in low energies; hence the truncation of the perturbative expansion cannot be carried out. However, it is an important subject to study the spectrum of particles predicted by QCD.Since perturbation theory is not applicable, a nonperturbative approach has to be used to study systems that involve strong interactions. Some of the nonperturbative approaches to strongly interacting systems are the QCD sum rules, quark models, and lattice QCD. The advantage of QCD sum rules and lattice QCD is that they are based on QCD itself, whereas, in quark models, one assumes a potential energy between the quarks and solves a Schrödinger-like equation. The advantage of quark models, on the other hand, is that it allows one to study also the excited states, whereas, in QCD sum rules and lattice QCD, only the ground state or in some exceptional cases the first excited state can be studied.In quark models one assumes a potential interaction among quarks which makes model as a nonrelativistic approach. Therefore, the systems that are best suited for study in quark models are the heavy quark system which containc or b quarks. The bare masses of u, d, and s quarks are 2 MeV, 4 MeV, and 96 MeV, respectively [2]. At a first look, quark model seems rather difficult to apply to light quarks. Capstick et al. presented reasonable explanations to link quark models including a minimal amount of relativity to the basics of QCD [3]. Although the pole masses of u, d and s quarks are very low and hence they are relativistic, in constituent quark models, instead of treating the physical u, d, and s quarks, one treats the so-called constituent quarks, which are nothing else than quarks dressed by gluons and other sea quarks inside the hadron. The masses of constituent quarks are around 300 MeV and hence they can also be treated in nonrelativistic quark models. Such an approach has been applied to light quark systems with a surprising success [4–6], leading to that model so-called Constituent Quark Model (CQM), which, based on the Gell Mann-Zweig idea, explains meson and baryon bound systems.A different situation is for heavy-light quark systems(Qq¯). For example an electron is more relativistic in the hydrogen atom (p,e-) than in the positronium atom (e+e-) [7]. Positronium can be taken as a naive model for quarkonium. The binding energy of the positronium is half of the hydrogen atom and is small compared to the electron mass. For this reason the positronium bound state can be described by nonrelativistic quantum mechanics. But the decay of the positronium resonance is a purely relativistic phenomenon. Nevertheless, we can attempt to apply the quark model to heavy-light mesons. The outcome of this attempt is not directly using of Heavy Quark Symmetry (HQS), but one aspect of it. Mesons are two particle systems and the reduced mass is dominated by the light quark mass, 1/μ=1/m+1/M≃1/M if M≫m. The spectra for (cq¯) and (bq¯) should be very similar under this assumption [7]. Indeed reasonable spectroscopy of D and B mesons can be obtained. There is a rich literature for the spectrum and dynamics of the heavy-light mesons, for example [8–25]. In [26], they studied semileptonic D and Ds decays based on the predictions of the relevant form factors from the covariant light-front quark model. In [27], the authors studied the Cabibbo-Kobayashi-Maskawa matrix element |Vub| which is not determined up to now in inclusive or exclusive B decays.Light quark physics is a key topic to understand the nature of QCD. They can be thought of a probe of the strong interactions by means of nonperturbative effects [28]. Heavy-light meson systems (Qq¯) is also central to enlighten the nature of QCD and strong interactions. Heavy-light meson spectroscopy has been the subject of both theoretical and experimental studies since the 2000s. Especially in the charm sector, new excited states were observed in D and Ds mesons [28–32].An important feature of B meson physics is that it is sensitive to New Physics (NP) Beyond the Standard Model (BSM) via rare decays. Furthermore hadronic decay channels of B mesons might have more systematic uncertainties due to the model indetermination, compared to the lepton/photon decay channels. Thus studyingB→lepton/photon decays present a play field for the search of NP. Besides that, B factory experiments BaBar and Belle were built to test the description of quark mixing in the Standard Model. The first theoretical description of quark mixing was proposed by Cabibo in 1963 [33]. One year later in 1964, Christenson et al. discovered CP violation in neutral kaon decays with a tiny friction [34]. This phenomenon is referred to as conclusion that matter and antimatter might behave differently. Kobayashi and Maskawa generalized Cabibbo’s idea by adjusting new quarks to the model [35]. In the framework of Standard Model, CP violation can be accommodated by introducing a complex phase in the 3×3 unitary Cabibo-Kobayashi-Maskawa (CKM) matrix. Indeed this phase can be measured in experiments. The cost of adding a parameter is to use a third generation of quarks. CP violation also occurs in B decays. The B factories were built to test for this purpose. B factories gave a substantial contribution to particle physics such as first observation of CP violation apart from the kaon, measurements of CKM matrix elements, measurements of purely leptonic B meson decays, and searches for new physics.In this work, we obtained mass spectrum and decay constants of theD and B mesons via QCD Sum Rule and a Quark Model potential. We also predicted decay constant for the Bs meson where there are no specific experimental data. Harmonic oscillator wave function is used in the quark model and a sufficiently trivial interpolating current is used in QCD Sum Rule calculations. We studied ground states since they are accessible in the framework of QCD Sum Rules.


## 2. QCD Sum Rule Formalism

In perturbation theory we assume that the eigenvalues and eigenfunctions can be expanded in a power series as follows:(1)En=En0+λEn1+λ2En2+⋯n=n0+λn1+λ2n2+⋯,where n is the principal quantum number and λ is a parameter. These series are in principle divergent, but they are asymptotic. This means that when the perturbation parameter is small, the first two or three terms are convergent so that the rest of the series can be ignored. In the case of QCD, due to the largeness of the parameter in lower energies, such a truncation, cannot be performed. The nonperturbative aspect of QCD makes it almost impossible to study bound states in terms of perturbation theory. For this reason, there is a need of nonperturbative methods to overwhelm this situation and study bound states. Among others such as Effective Field Theory and Lattice QCD, QCD Sum Rule is maybe the most popular nonperturbative method.QCD Sum Rule is first formulated by Shifman, Vansthein, and Zakharov for mesons in [36] and generalized to baryons by Iofe in [37]. The basic idea of the this formalism is to study bound state phenomena in QCD from the asymptotic freedom side, i.e., to start evaluation of correlation function at short distances, where quark-gluon dynamics are perturbative and move to larger distances where hadronization occurs, including nonperturbative effects and using some approximate procedure to get information on hadronic properties [38].To obtain physical observables from QCD sum rules, a correlator of two hadronic currents which is defined as follows(2)Π=i∫d4xeipx0Tjxj†00,is studied. Here p is momentum and j(x) is a current composed of quarks and gluon fields with the hadron’s quantum numbers. When this operator is applied to vacuum, it can create the hadron that we study. Equation (2) is known as correlation function. The fundamental assumption of the QCD sum rules is that there is a region of p where correlation function can be equivalently described at both quark and hadron sector. The former is known as QCD or OPE (Operator Product Expansion) side, and the latter is known as the phenomenological side. Matching these two sides of the sum rule, one can obtain information about hadron properties [38].Forp2>0, resolution of identity operator of hadron states can be written between the operators. This results in correlation function as follows:(3)Π=∑h0jhp1p2-mh2hpj†0+higherstates.It can be seen from (3) the poles in the correlation function, which indicates the presence of hadrons, created by operator j(x).For-p2≫ΛQCD2(p2<0), major contribution to the correlation function will come from the x~0 region [39]. In this case the product of two operators can be written in terms of OPE:(4)Tjxj†0=∑dCdxOd.Here Cd(x) are the coefficients, which can be calculated by the perturbation theory, and Od′s are the operators with the mass dimension d. If Fourier transformation applies to (4), correlation function can be written as follows:(5)Π=∑dCdfpOdpd,where 〈Od〉 are the vacuum condensates that cannot be calculated by perturbation theory except d=0. d=0 corresponds to unitary operator and can be calculated via perturbation theory. Other operators can be written as 〈q¯q〉 (d=3); mq〈q¯q〉 (d=4), 〈GμνGμν〉 (d=4), 〈q¯gσGq〉 (d=5). For d=1 and d=2 there exists no operator. As a result of this, the expansion converges quickly although it is an infinite summation.In order to get sum rules we must equate (3) and (5). But these two expressions are obtained in different regions of p. By using spectral density representation of correlation function, this matching can be made:(6)Πp2=∫0∞ρss-p2+polynomialsofp2.Spectral densityρ(s) can be acquired from (3). Inserting ρ(s) into (6), an expression of correlation function can be obtained from (3) for p2<0 region. If we denote ρphen(s) as spectral density from (3) and ρQCD(s) from (5), we get(7)∫0∞dsρphenss-p2+polynomials=∫0∞dsρQCDss-p2+polynomials.In order to extract physical properties from this expression, one must eliminate the polynomial terms, for example, by using derivatives. In principle, no one knows the polynomial degree and how many polynomials are. The correct procedure is then to use the Borel transformation, which contains infinite derivative:(8)BM2Πq2=lim-q2,n→∞,-q2/n=M2-q2n+1n!ddq2nΠq2.Here M2 is defined as the Borel parameter [36]. This transformation effectively removes the polynomials and makes(9)1s-p2→e-s/M2.Then,(10)∑h0jhp2e-mh2/M2+higherstates=∫0∞ρQCDse-s/M2,which resembles QCD parameters and hadronic properties. This equation still shows presence of unknown parameters. The e-mh2/M2 factor makes the contribution of small masses dominant. To parameterize contributions of higher states, quark-hadron duality approximation is used. According to quark-hadron duality, for s>s0, ρphen(s)≃ρQCD(s). ρphen(s) has contribution of higher states and heavier hadrons when s>s0. s0 is called as continuum threshold and is related mass of the hadron that is studied in sum rules. So, we can write (10) as follows:(11)0jmhp2e-mh2/M2=∫0s0ρQCDse-s/M2.In this equation mh is the hadron of the smallest mass which can be created by j.Physical properties extracted from the sum rules must be independent of Borel parameter, (M2). Here we assume that there exist a range of M2, called Borel window, in which two sides have a good overlap and information on the lowest state can be extracted. Minimum and maximum values of Borel window can be extracted in a way that QCD side convergence gives the minimum value, and the condition that pole contribution should be bigger than the continuum contribution gives the maximum value of Borel window [38].


### 2.1. Mass Sum Rule

The mass sum rule can be obtained by matching QCD and phenomenological sides of correlation function [36, 38, 39]. Here we will give the following formula:(12)m2=∫smins0dse-s/M2sρQCDs∫smins0dse-s/M2ρQCDs.


### 2.2. Decay Constant

The decay constant can be obtained from the formula [40] as follows:(13)fmh2=emh/M21mh2∫smins0dse-s/M2ρQCDs,where mh is the hadron mass extracted from sum rules.


## 2.1. Mass Sum Rule

The mass sum rule can be obtained by matching QCD and phenomenological sides of correlation function [36, 38, 39]. Here we will give the following formula:(12)m2=∫smins0dse-s/M2sρQCDs∫smins0dse-s/M2ρQCDs.


## 2.2. Decay Constant

The decay constant can be obtained from the formula [40] as follows:(13)fmh2=emh/M21mh2∫smins0dse-s/M2ρQCDs,where mh is the hadron mass extracted from sum rules.


## 3. Quark Model

Also known as potential model or quark potential model, quark model considers one or more interacting particles under a given potential. In the early 60s quarks were modelled and experimental evidences were found subsequently. This approach provided a reliable basis to study and investigate particle physics and gave compatible results with the experiments.The most important part of the quark model is the potential. After the November revolution of particle physics in 1974, the year in which charmonium states were observed, new models were proposed to calculate spectrum and radiative transitions [41–43]. The so-called Cornell potential, proposed in [42], reads as(14)Vr=-ar+br+c,wherea, b, and c are some parameters to extract from fit to the experimental data. This potential is still used with some modifications to account, for example, for hyperfine splittings in the energy levels. The other potentials such as power law potential [44], logarithmic potential [45], Richardson potential [46], Buchmüller-Tye potential [47], and Song-Lin potential [48] were used to fit quarkonium spectra and gave good results in agreement with experiments. These were phenomenological spin-independent potentials and not directly QCD motivated. The interquark potential was not derived from first principles of QCD in the early quarkonium phenomenology. This means, in terms of QCD, that potential is universal (flavour independent) and since quarks are colorless particles, it was reasonable to assume the universality as valid, despite the fact that gluons couple to color charge. These spin-independent potential models performed good but not complete explanation of the energy level splittings. If we want to accommodate these splittings in the theory, we have to take care, i.e., ofspin-spin andspin-orbit interactions in the model. Reference [49] reports an example of a QCD-motivated, spin- and velocity-dependent potential. These potentials deliver reliable results.


## 4. Elaboration of the Problem

### 4.1. QCD Sum Rules

In QCD sum rules, the choice of thej(x) current is important, since it creates hadrons from vacuum. We used the following current:(15)jx=iQ¯axγ5qax,whereQ is heavy quark, q is light quark, a is the color index, and γ5 is the Dirac matrix. We take care of mq→0 limit. In the limit of mq→0, there appears a flavor symmetry between b and c quarks. By this symmetry it is possible to extract information about c and b sector with the same current. b and c quarks are heavy quarks so that it cannot be expected to be in the vacuum by themselves. So it is possible to ignore such condensate terms like 〈αs(GG/π)〉 and 〈q¯gsσGq〉. By introducing the current term into (2), one can obtain the following spectral density:(16)Fs0,M2=-qq¯e-mQ2/M2mQ+6e-s0/M2e-su/M2u1-u2×e-s0/M2M2mQ2+su+M2-e-su/M2mQ4+M2s0+M2where〈qq¯〉 is the condensate, and u1 and u2 are solutions of s(u)=mQ2/(1-u)+mq2/u=s0.The mass sum rule can be obtained by taking derivative with respect to1/M2 and dividing the result by (16):(17)mh2=M41Fs0,M2dFs0,M2dM2.The decay constant sum rule can be obtained as follows:(18)fmh2=emh2/M21mh2Fs0,M2.The mass values and decays constants for heavy-light mesons are presented in Tables1 and 2 and Figures 1–8.Table 1
Mass spectra of heavy-light mesons in MeV. QM denotes quark model and SR denotes sum rule calculations. The parameters areκ=0.471, a=0.192GeV2, mc=1.320GeV, mb=4.740GeV [52], 〈q-q〉=0.241GeV3, mu=md=0.340GeV, and ms=0.600GeV.
Meson Exp. [2] QM SR [9] [14] D 0 / D + 1869.3± 0.4 1859 1972± 94 1870.82 1854.7 Ds+ 1968.2± 0.4 2056 2118± 75 1966.62 1974.5 B + / B 0 5279.0± 0.5 5260 5259± 109 5273.50 5277.2 B s 0 5367.7± 1.8 5442 5488± 76 5365.99 5384.8Table 2
Pseudoscalar and vector decay constants of heavy-light mesons in MeV. QM denotes quark model and SR denotes sum rule calculations.
Meson Exp. QM SR [9] [16] [53] D 0 / D + 206± 8.9 199 210.25± 11.60 205.14 206.2± 7.3 ± 5.1 207.53 Ds+ 249 253 245.70± 7.46 241.84 245.3± 15.7 ± 4.5 262.56 B + / B 0 204± 31 209 223.45± 12.4 201.09 193.4± 12.3 ± 4.3 208.13 B s 0  275 277.22± 11 292.04 232.5± 18.6 ± 2.4 262.39Figure 1
Borel parameter dependence of theD0/D+ masses.Figure 2
Borel parameter dependence of theDs+ mass.Figure 3
Borel parameter dependence of theB+/B0 masses.Figure 4
Borel parameter dependence of theBs0 mass.Figure 5
Borel parameter dependence of theD0/D+ decay constants.Figure 6
Borel parameter dependence of theDs+ decay constant.Figure 7
Borel parameter dependence of theB+/B0 decay constants.Figure 8
Borel parameter dependence of theBs0 decay constant.


### 4.2. Quark Model

Energy eigenvalues can be obtained by solving the Schrödinger equation in the quark model. The Schrödinger equation reads as follows:(19)HΨn=EnΨn,wheren denotes the principal quantum number. We can separate the wave function into radial Rnl and angular parts Ylm(θ,ϕ) as follows:(20)Ψnlmr,θ,ϕ=RnlrYlmθ,ϕ.R n l is the radial wave function given as follows:(21)Rnl=Nnlrle-νr2Ln-l/2l+1/22νr2,with the associated Laguerre polynomials L(n-l)/2l+1/2 and the normalization constant:(22)Nnl=2ν3π2n-l/2!νln+l/2+1!!.With the wave function in hand one can obtain masses as well as decay constants for heavy and light mesons. The mass spectra can be obtained by solving (19). For the decay constants we employ the following formulas:The results are as follows:(23)fp=3mp×∫d3k2π31+mqEk×1+mq¯Ek¯×1-k2Ek+mqEk¯+mq¯ϕk→,for pseudoscalar mesons and(24)fv=3mv×∫d3k2π31+mqEk×1+mq¯Ek¯×1+k23Ek+mqEk¯+mq¯ϕk→,for the vector mesons [50].In the nonrelativistic limit, these two equations take a simple form, which is known to be Van Royen and Weisskopf relation [51]. For the meson decay constants,(25)fp/v2=12Ψp/v02mp/v.Here mp/v denotes the pseudoscalar and vector mass of the related meson.The results are shown in Tables1 and 2.


## 4.1. QCD Sum Rules

In QCD sum rules, the choice of thej(x) current is important, since it creates hadrons from vacuum. We used the following current:(15)jx=iQ¯axγ5qax,whereQ is heavy quark, q is light quark, a is the color index, and γ5 is the Dirac matrix. We take care of mq→0 limit. In the limit of mq→0, there appears a flavor symmetry between b and c quarks. By this symmetry it is possible to extract information about c and b sector with the same current. b and c quarks are heavy quarks so that it cannot be expected to be in the vacuum by themselves. So it is possible to ignore such condensate terms like 〈αs(GG/π)〉 and 〈q¯gsσGq〉. By introducing the current term into (2), one can obtain the following spectral density:(16)Fs0,M2=-qq¯e-mQ2/M2mQ+6e-s0/M2e-su/M2u1-u2×e-s0/M2M2mQ2+su+M2-e-su/M2mQ4+M2s0+M2where〈qq¯〉 is the condensate, and u1 and u2 are solutions of s(u)=mQ2/(1-u)+mq2/u=s0.The mass sum rule can be obtained by taking derivative with respect to1/M2 and dividing the result by (16):(17)mh2=M41Fs0,M2dFs0,M2dM2.The decay constant sum rule can be obtained as follows:(18)fmh2=emh2/M21mh2Fs0,M2.The mass values and decays constants for heavy-light mesons are presented in Tables1 and 2 and Figures 1–8.Table 1
Mass spectra of heavy-light mesons in MeV. QM denotes quark model and SR denotes sum rule calculations. The parameters areκ=0.471, a=0.192GeV2, mc=1.320GeV, mb=4.740GeV [52], 〈q-q〉=0.241GeV3, mu=md=0.340GeV, and ms=0.600GeV.
Meson Exp. [2] QM SR [9] [14] D 0 / D + 1869.3± 0.4 1859 1972± 94 1870.82 1854.7 Ds+ 1968.2± 0.4 2056 2118± 75 1966.62 1974.5 B + / B 0 5279.0± 0.5 5260 5259± 109 5273.50 5277.2 B s 0 5367.7± 1.8 5442 5488± 76 5365.99 5384.8Table 2
Pseudoscalar and vector decay constants of heavy-light mesons in MeV. QM denotes quark model and SR denotes sum rule calculations.
Meson Exp. QM SR [9] [16] [53] D 0 / D + 206± 8.9 199 210.25± 11.60 205.14 206.2± 7.3 ± 5.1 207.53 Ds+ 249 253 245.70± 7.46 241.84 245.3± 15.7 ± 4.5 262.56 B + / B 0 204± 31 209 223.45± 12.4 201.09 193.4± 12.3 ± 4.3 208.13 B s 0  275 277.22± 11 292.04 232.5± 18.6 ± 2.4 262.39Figure 1
Borel parameter dependence of theD0/D+ masses.Figure 2
Borel parameter dependence of theDs+ mass.Figure 3
Borel parameter dependence of theB+/B0 masses.Figure 4
Borel parameter dependence of theBs0 mass.Figure 5
Borel parameter dependence of theD0/D+ decay constants.Figure 6
Borel parameter dependence of theDs+ decay constant.Figure 7
Borel parameter dependence of theB+/B0 decay constants.Figure 8
Borel parameter dependence of theBs0 decay constant.


## 4.2. Quark Model

Energy eigenvalues can be obtained by solving the Schrödinger equation in the quark model. The Schrödinger equation reads as follows:(19)HΨn=EnΨn,wheren denotes the principal quantum number. We can separate the wave function into radial Rnl and angular parts Ylm(θ,ϕ) as follows:(20)Ψnlmr,θ,ϕ=RnlrYlmθ,ϕ.R n l is the radial wave function given as follows:(21)Rnl=Nnlrle-νr2Ln-l/2l+1/22νr2,with the associated Laguerre polynomials L(n-l)/2l+1/2 and the normalization constant:(22)Nnl=2ν3π2n-l/2!νln+l/2+1!!.With the wave function in hand one can obtain masses as well as decay constants for heavy and light mesons. The mass spectra can be obtained by solving (19). For the decay constants we employ the following formulas:The results are as follows:(23)fp=3mp×∫d3k2π31+mqEk×1+mq¯Ek¯×1-k2Ek+mqEk¯+mq¯ϕk→,for pseudoscalar mesons and(24)fv=3mv×∫d3k2π31+mqEk×1+mq¯Ek¯×1+k23Ek+mqEk¯+mq¯ϕk→,for the vector mesons [50].In the nonrelativistic limit, these two equations take a simple form, which is known to be Van Royen and Weisskopf relation [51]. For the meson decay constants,(25)fp/v2=12Ψp/v02mp/v.Here mp/v denotes the pseudoscalar and vector mass of the related meson.The results are shown in Tables1 and 2.


## 5. Summary and Conclusions

In this paper, we calculated mass spectra and decay constants of pseudoscalar and vector heavy-light mesons (B, Bs, D, and DS) in the framework of QCD sum rule and quark model. Obtained results for masses of B and D mesons are in good agreement with the available experimental data. In the mass spectra, the extrapolation via quark model gave close results to experimental data compared to the QCD sum rule consideration. The QCD sum rules approach gives reasonable but not very good-matching results compared to the experimental values, because of the adopted approximation when evaluating the current, whereas the higher dimension of that operator could improve the estimates. Other potentials and further studies should be taken into consideration for a better understanding.The heavy-light mesons under study in this paper are well established indeed, and any prediction or reproduction of mass spectrum does not directly guarantee the validity of the model, but shows a possible path to follow for a further investigation. Therefore other physical observables such as decay constants should be experimentally investigated to give more inputs to the theory. For example the only precise value of decay constant is known forD mesons, and systematics are evaluated. The other mesons in this study need more experimental data. For Bs there is no available experimental data. We predicted for the first time decay constant value for Bs in this manner.Decay constants give information about short distance structure of hadrons. The obtained results for decay constants are in agreement with the other studies and available data. We did not consider in this work relativistic corrections.In QCD Sum Rule calculations, physical observables must be independent of the Borel parameter. In Figures1–8 the smoothness of the graphs is compatible with existing data. It is worthy to note that in Figures 1 and 5 the ’slope’ of the two curves of D0/D+ is not in the same range. The reason for that could be the smallness of the Borel parameter and continuum threshold energy, since correlation function receives main contribution at s≠M2. On the other hand, the smallness of Borel parameter can blow up the corrections to the perturbative part of the correlation function.In summary, we obtained good results in accordance with the available data and theoretical studies. As mentioned before, other potential models and interpolating currents can be used to study mass spectra and decay constants. Heavy-light systems in view of the quark model are important to study hadronic interactions. In particular, Heavy Quark Spin Symmetry can play an essential role in heavy-light systems. The higher dimensions of the operators in interpolating currents would deliver more accurate results.



---

*Source: 8095653-2018-09-25.xml*