Abstract
We study the heavy-light mesons within basis light-front quantization. The resulting mass spectra of D, \(D_s\), B, and \(B_s\) agree reasonably well with experiments. We also predict states which could be measured in the near future. In the light-front formalism, we calculate the light-front wave functions and additional experimental observables, such as parton distribution functions, distribution amplitudes, and decay constants by means of integrations over light-front wave functions. We also provide ratios of decay constants for selected pseudoscalar meson decays (\(D_s\) to D and \(B_s\) to B) as they may prove to be theoretically more robust and more reliably determined in experiments. We find that our ratios are systematically smaller than existing experiments and other approaches by 5–18%.
Similar content being viewed by others
Explore related subjects
Discover the latest articles and news from researchers in related subjects, suggested using machine learning.Avoid common mistakes on your manuscript.
1 Introduction
The motivation for studying heavy-light mesons (D, \(D_s\), B, \(B_s\)) with basis light-front quantization (BLFQ) [1] is twofold. First, the combination of a heavy and a light quark is the closest QCD analogue of the hydrogen atom in QED, so that similarities and differences in spectroscopic features could inform discussions of the relative roles of gauge-boson exchange and confinement. Second, the successful applications of BLFQ to heavy meson systems [2, 3] provides a foundation for understanding the roles of key elements of the quark–antiquark effective Hamiltonian adopted for the heavy-light system. We employ the same form of Hamiltonian with only two fit parameters in order to test the validity of our model as well as some novel behavior for the heavy-light mesons.
Heavy-light systems have been studied widely by various methods for a long time, and different approaches emphasize certain features of the system and have their own advantages. Examples include the heavy-quark effective theory [4]; models that incorporate chiral dynamics owing to the light constituent quark [5, 6]; QCD sum rules [7]; Dyson-Schwinger equations approach [8]; etc. Nevertheless, there are detailed properties of the excitation spectrum of heavy-light mesons that are yet to be fully understood. In addition, recent discoveries of new excited states, including the first heavy flavored spin-3 resonance \(D^*_{s3}(2860)^\pm \) state, and \(B_J(5840)^{0,+},\ B_J(5960)^{0,+}\) [9, 10] have stimulated renewed quests for theoretical understandings.
Additional motivation stems from the fact that the leptonic decay constant of heavy-light mesons plays a significant role in determining the Cabibbo–Kobayashi–Maskawa (CKM) quark-mixing matrix elements within Standard Model [11, 12]. Moreover, with continually improving experimental precision, these observables constitute powerful tools to probe physics beyond Standard Model [13].
In this paper, we discuss the heavy-light mesons within BLFQ, a non-perturbative Hamiltonian framework. The paper is organized as follows. Section 2 outlines the elements of BLFQ, the theoretical approach to the heavy-light systems. Then in Sect. 3, we exhibit the numerical results of the mass spectrum, light-front wave functions, parton distribution functions, distribution amplitudes, and the decay constants. Section 4 provides our discussions and outlook.
2 Basis light-front quantization for heavy-light systems
In this work, we adopt the effective Hamiltonian within the \(|q{\bar{q}}\rangle \) Fock sector in the form introduced for the heavy mesons [2, 3]. It comprises two parts \(H_\mathrm{tot} = H_0 + V^\mathrm{eff}_g\), where
is the holographic QCD Hamiltonian [14, 15] augmented by massive quark kinematics and the longitudinal confinement [2, 3]. In addition, the spin structure of the hadrons is generated by the effective one-gluon-exchange potential \(V^\mathrm{eff}_g\). Along with a running coupling \(\alpha _s(Q^2)\), \(V^\mathrm{eff}_g\) governs the short-range physics:
With the Hamiltonian, one can solve the light-front eigenvalue problem \(H_\mathrm{tot} |\varPsi _h(P,J,m_J)\rangle = M^2 |\varPsi _h(P,J,m_J)\rangle \) for the bound state eigenmasses and light-front wave functions (LFWFs). In the BLFQ framework, we expand the LFWFs in terms of basis functions for the transverse and longitudinal directions, respectively:
where \(k_\perp \equiv |{\mathbf {k}_\perp }|\). The transverse basis functions,
are the 2D harmonic oscillator functions where n and m are the principal and orbital quantum numbers, respectively; \(\mathbf {q}_\perp = \mathbf {k}_\perp /\sqrt{x(1-x)}\), with x being the longitudinal momentum fraction carried by the quark \(x=p^+_q/(p^+_q+ p^+_{{\bar{q}}})\); \(\theta _q = \arg \mathbf {q}_\perp \), b sets the scale of the basis, and \(L_n^{|{m}|}\) is the associated Laguerre polynomial. The longitudinal basis functions are given by
where \(P_l^{(\alpha ,\beta )}(2x-1)\) is the Jacobi polynomial with quantum number l; \(\alpha \) and \(\beta \) are two dimensionless parameters associated with the constituent (anti-) quark masses [15]:
In practical calculations, we truncate the infinite basis by restricting the quantum numbers:
Due to the transverse cutoff \(N_\mathrm{max}\), the ultraviolet (UV) and infrared (IR) regulators are implicit through \(\varLambda _\mathrm{UV} \approx b\sqrt{N_\mathrm{max}}\), and \(\varLambda _\mathrm{IR} \approx b/\sqrt{N_\mathrm{max}}\), respectively. \(L_\mathrm{max}\) represents the basis resolution in the longitudinal direction. In this work, we choose \(N_\mathrm{max} = L_\mathrm{max} \) for simplicity. Notice that the total angular momentum projection is defined by \(m_J = m + s +{\bar{s}}\), which is conserved in the constructed basis. Solving the eigenequation by diagonalizing the \(H_\mathrm{tot}\) matrix, one obtains the eigenvalues which indicate the spectra as squared masses, and the eigenvectors which are the coefficients \( \psi _h(n,m,l,s,{\bar{s}}) \) in Eq. (3).
3 Numerical results
In the framework of BLFQ, the confining strength and constituent (anti-) quark masses are typically taken as free parameters that are fitted to experiments. Following the practice adopted for the heavy quarkonia systems [3, 15], we take the confining strength of longitudinal and transverse directions to be the same as reflected in Eq. (1), and we set the 2D harmonic oscillator scale parameter by \( b = \kappa \). We adopt the running coupling as regularized in Ref. [15] now with the flavor number (\(N_f\)) dependence listed in Table 1. Our convention is to set the heavier quarks, charm (\({\mathsf {c}}\)) and bottom (\({\mathsf {b}}\)) to be the quarks, while setting the lighter quarks, up (\({\mathsf {u}}\)), down (\({\mathsf {d}}\)), and strange (\({\mathsf {s}}\)) to be the antiquarks in the heavy-light systems. We take the masses of \({\mathsf {u}}\) and \({\mathsf {d}}\) antiquarks to be the same. Unless otherwise stated, we use B to denote \(B^0\) and \(B^-\), and D for \(D^0\) and \(D^+\).
Following the scheme adopted in Ref. [3], the confining strength of the flavored meson is implemented as \(\kappa _{{\mathsf {x}}\bar{{\mathsf {y}}}} = \sqrt{(\kappa ^2_{{\mathsf {x}}\bar{{\mathsf {x}}}} +\kappa ^2_{{\mathsf {y}}\bar{{\mathsf {y}}}})/2}\). We adopt \(\kappa _{{\mathsf {u}}\bar{{\mathsf {u}}}/{\mathsf {d}}\bar{{\mathsf {d}}}} = \kappa _{{\mathsf {s}}\bar{{\mathsf {s}}}} = 0.59\text { GeV}\)Footnote 1 as the confining strength for the light mesons [16], while the values for heavy quarkonia, \(\kappa _{{\mathsf {b}}\bar{{\mathsf {b}}}} \) and \(\kappa _{{\mathsf {c}}\bar{{\mathsf {c}}}} \), are taken from Ref. [2]. Following our previous works on heavy systems [2, 3], we keep the same heavy quark masses (\(m_{\mathsf {b}}\), \(m_{\mathsf {c}}\)), while we fit the light-quark masses (\(m_{\mathsf {u}}/m_{\mathsf {d}}\), \(m_{\mathsf {s}}\)) to reproduce the heavy-light meson mass spectra from experiments.
More specifically, we tune only the \(m_{\mathsf {u}}/m_{\mathsf {d}}\ (m_{\mathsf {s}})\) to minimize the r.m.s. mass deviation between experiments and theory for the lowest pseudoscalar and vector states, \(D^0\), \(D^*(2007)\), \(B^\pm \), \(B^*\) (\(D_s^\pm ,\ D_s^{*\pm }\), \(B_s^0,\ B_s^{*0}\)). The \(m_{\mathsf {u}}/m_{\mathsf {d}}\) quark mass we obtain agrees with the value fitted directly to the light mesons excluding the pion [17].
The model parameters of heavy-light mesons are listed in Table 1. We will elaborate on the details of our heavy-light mass spectra in the following section. While heavy (e.g. \(\eta _c\)) and light mesons (e.g. \(\rho , \pi \)) are not the subject of interests in this work, we nevertheless include the obtained values for comparison purpose only. Note that due to our omission of chiral dynamics, these obtained light meson masses differ substantially from the experimental measurement, as expected. Within the same framework, Ref. [18] shows that incorporation a NJL type chiral symmetry breaking interaction dramatically improves the predictions for \(\rho \), \(\pi \) masses and additional observables.
3.1 Mass spectrum
Here we present mass spectra of the four heavy-light meson systems, with all 12 low-lying states in Fig. 1. For each of the BLFQ calculated states, we show the mass as a black box with a dashed line, where the dashed line represents the mean value of the invariant mass,
and the box indicates the spread of eigenmasses obtained from different \(m_J\)’s: \(\delta _J M \equiv \max (M_{m_J}) - \min (M_{m_J})=\) box height. In many cases, the box height is small and not visible in the figure. Since rotational symmetry would imply degeneracy (zero box height), we introduce an overall mean spread for hadrons (h) within a fixed flavor to characterize the violation of rotational symmetry,
listed in the last column of Table 1. For each system, \(\overline{\delta _J M}\) is obtained with 9 states whose \(J > 0\), which includes three \( J^P = 1^-\) states, two \(1^+\), one \(2^+\), two \(2^-\), and one \(3^-\). We observe that the mean spread has the tendency of decreasing with increasing meson mass, which is in agreement with our expectation. Viewed as a percentage deviation (see the parentheses in the last column of Table 1), the decrease with increasing meson mass is in a trend consistent with the trend found for all light and heavy mesons [2, 3].
The hyperfine splitting of ground state pseudoscalar and vector is a particularly sensitive test of the spin-sensitive component of the Hamiltonian – the effective one-gluon exchange interaction. It is well-known that for the light mesons, this mass splitting is driven by the chiral symmetry breaking [19]. However, in our approach, the splitting is brought in by the one-gluon exchange interaction and its interplay with the confining strength \(\kappa \) when the basis size is fixed by \(N_\mathrm{max}\) and \(L_\mathrm{max}\). Specifically, smaller \(\kappa \) tends to lead a larger mass splitting between the two states. The size of the hyperfine splittings are reasonable for heavy mesons [2, 3, 15] and heavy-light mesons within this model. However, we observe that for BLFQ to reproduce the properties of the light mesons, such as the Goldstone nature of the pions in the chiral limit, the Hamiltonian needs to include chiral dynamics [18].
We compare our mass spectra with experimental values summarized by the Particle Data Group (PDG) [20], and Lattice QCD [21,22,23,24] in Fig. 1. Most of our results are within the quoted uncertainties of either experiments or Lattice. Some states, \(J^P = 0^+\) or \(1^+\) for instance, show somewhat larger differences among the states compared. In BLFQ, the accuracy may be enhanced by introducing higher Fock sectors, such as including a dynamical gluon in the Fock sectors. Nevertheless, our work provides the LFWFs that can be used to calculate some hadron observables, which may be tested by ongoing and forthcoming experiments, such as the Belle II experiment at the SuperKEKB accelerator [25].
Mass spectra of the four heavy-light meson systems, each with 12 low-lying states, using \(N_\mathrm{max}=L_\mathrm{max}=32\). Horizontal axis refers to the \(J^P\) values of the states. The green dashed lines indicate the corresponding open flavor threshold. Black boxes are the results of this work, indicating the highest and lowest mass from different \(m_J\)’s, while the dashed lines in between are the averages \({\overline{M}}\). Red bars are the experimental masses compiled by PDG [20]. Specifically, we use the values of \(D^0 \ ({\mathsf {c}}\bar{{\mathsf {u}}})\) and \(B^\pm \ ({\mathsf {b}}{\bar{{\mathsf {u}}}}/{\mathsf {b}}\bar{{\mathsf {d}}})\) from experiments for comparison. Shaded blue bars are the Lattice QCD results from Refs. [21,22,23,24]. For both PDG and Lattice, we use the solid lines for the center values, while shaded boxes to indicate the uncertainties
3.2 Light-front wave function
The LFWFs obtained from solving the light-front Hamiltonian matrix eigenvalue problem provide direct access to hadron structure observables. The orthonormality condition of the LFWFs reads,
In fact, the LFWFs themselves, especially the unequal mass systems, illustrate interesting features. The equal-mass quarkonium LFWFs exhibit symmetry with respect to the transverse momentum and longitudinal momentum fraction. We anticipate and observe that asymmetry emerges when the constituent masses differ. In Fig. 2, we present the ground state spin-singlet LFWFs \(\psi ^{m_J=0}_{\uparrow \downarrow -\downarrow \uparrow }(\mathbf {k}_\perp ,x) \) of three meson states, \(\eta _b({\mathsf {b}}\bar{{\mathsf {b}}})\), \(B_c({\mathsf {b}}\bar{{\mathsf {c}}})\), and \(B({\mathsf {b}}\bar{{\mathsf {u}}})\), where we define \(\psi ^{m_J=0}_{\uparrow \downarrow \pm \downarrow \uparrow }(x, \vec {k}_\perp ) \equiv 1/\sqrt{2}[ \psi ^{m_J=0}_{\uparrow \downarrow }(x,\vec {k}_\perp ) \pm \psi ^{m_J=0}_{\downarrow \uparrow }(x,\vec {k}_\perp ) ] \). For the sake of convenience, we visualize the LFWFs at \(k_y=0\), i.e. \(\theta _q=0\) or \(\pi \) in the phase \(\exp ({\mathrm {i}}m \theta _q)\) (cf. Eq. 4).
The LFWFs of three listed mesons are symmetric in the transverse direction. The equal-mass (quarkonium) case, Fig. 2a, also shows the expected symmetry in the longitudinal direction. However, the unequal-mass cases, Fig. 2b, c, show the anticipated asymmetry in the longitudinal direction. Specifically, the peak of LFWFs for the quarkonia is located at \(x=1/2\) on the longitudinal direction, which is the same as the quark mass fraction \(m_q/(m_q+m_{{\bar{q}}})\) where \(m_q = m_{{\bar{q}}}\). However, for the unequal-mass system, the peak location is nontrivial as shown in Fig. 3, where we present the density plot of the ground state spin-singlet LFWFs of D and B. At \(k_\perp = 0\) the asymmetry is maximal, and the peak is located at \(x > m_\mathsf {c(b)} / (m_\mathsf {c(b)} + m_{\mathsf {u}})\), but as \(k_\perp \) increases the peak location (along fixed \(k_\perp \)) shifts to smaller values of x. In the limit \(k_\perp \rightarrow \infty \) the peak location approaches \(x = 1/2\). This can easily be understood because for \(k_\perp \gg m_\mathsf {c(b)}+m_{\mathsf {u}}\), the quark masses become irrelevant, and the LFWF approaches that of equal-mass constituents. These features confirm the prediction from parton gas model [26]. The nontrivial dependence of the LFWF on the quark masses can also be found in \(B_c\) meson [27], and it affects light-front observables such as the parton distribution function (PDF) and distribution amplitude (DA) as we will see in the following sections.
Density plot of the ground state (\(1^1S_0\)) LFWFs of D and B mesons. The dashed black line represents the quark mass fraction, which is \(x=m_\mathsf {c(b)}/(m_\mathsf {c(b)}+m_{\mathsf {u}})\) for D (B) meson. Red dots indicate the local maxima along x for fixed \(k_\perp \), and the orange curve is fitted from the dots based on the light-front parton gas model [26, 27], to illustrate the trend of peaks: for the heavy-light mesons, the peak is located at \(x>m_q/(m_q+m_{{\bar{q}}})\) when \(k_\perp =0\), and it shifts to smaller x as \({k_\perp }\) increases and finally approaches to \(x=1/2\) at the limit \({k_\perp } \rightarrow \infty \)
3.3 Parton distribution function
The structure of hadrons, as probed in inclusive deep inelastic scattering (DIS) at large momentum transfer, can be described by the PDF \(f(x; \mu )\). The PDF within collinear factorization defines the probability for finding a quark carrying a longitudinal momentum fraction x at resolution scale \(\mu \). In the LFWF representation, the PDF is calculated simply by integrating out the transverse momentum of the square of the wave function modulus,
and the truncation parameter \(N_\mathrm{max}\) of the LFWFs provides a UV scale for model \(\mu \approx \kappa \sqrt{N_\mathrm{max}}\), which is denoted schematically by the upper limit in the integration. That is, the integration is performed without any additional cutoff. Within the two-body Fock sector truncation, the PDF is normalized to unity \(\int _0^1 f(x; \mu ) {\mathrm {d}}x =1\), which relates to the orthonormal condition of the LFWFs (cf. Eq. 10).
We study the PDFs of the heavy-light systems, and present the results of pseudoscalar and vector states in Fig. 4. The location of peak reflects the asymmetry of the quark mass distribution in the heavy-light systems. Note that the peak of the PDFs is not at the quark mass fraction, due to the nontrivial behavior of the maxima in x along fixed \(k_\perp \) in the LFWFs. However, we find that the peak in the PDF for the radial excited states \(2^1S_0\) and \(2^3S_1\) is at a significantly smaller value of x than that of \(1^1S_0\) and \(1^3S_1\), and in fact it is very close to \(x = m_q / (m_q+m_{{\bar{q}}})\). The latter may be coincidental. The width of PDFs decreases as the mass of the system increases, and approaches a \(\delta \) function in the nonrelativistic limit. The \(2^1S_0\) and \(2^3S_1\) states contain bumps on both sides of the peaks incorporating features arising from radial excitations. The difference between solid and dashed curves reveals the spin excitation. Those differences between \(1^1S_0\) and \(1^3S_1\) states are larger than between \(2^1S_0\) and \(2^3S_1\), that is due to the significant hyperfine splittings in lower excited states.
Note that we use the basis size as \(N_\mathrm{max} = L_\mathrm{max} = 32\), which corresponds to different resolution scales: \(\mu _{D/D_s} \approx 4.5\) GeV and \(\mu _{B/B_s}\approx 6.0\) GeV. The PDFs obtained here are only for valence partons with a low resolution where the valence Fock sector approximation is reasonable. Ref. [28] shows that combining the low resolution PDFs with the Dokshitzer–Gribov–Lipatov–Altarelli–Parisi (DGLAP) evolution provides access to experiment-relevant PDFs (valence, sea and gluon) at higher scales. We anticipate that DGLAP could also be applied to our results but that is beyond the scope of the present effort.
PDFs of pseudoscalar (solid curves) and vector (dashed curves) states of the heavy-light systems at \(N_\mathrm{max} = L_\mathrm{max} = 32\), which is equivalent to the UV regulators \(\mu _{D/D_s} \approx 4.5\) GeV and \(\mu _{B/B_s}\approx 6.0\) GeV. The vertical dashed gray lines indicate the mass fraction of the quark, \(x=m_q/(m_q+m_{{\bar{q}}})\)
3.4 Distribution amplitude
Another light-cone distribution we investigate is the DA. DA controls the exclusive process at large momentum transfer that is defined from the lightlike vacuum-to-meson matrix elements. In the LFWF representation, it can be written as [29]
where \(f_{P(V)}\) are the decay constantsFootnote 2 of the pseudoscalar (vector) state, and is associated with the minus (plus) sign in the subscript of the wave function. Like in the PDF, we have the UV cutoff \(\mu \) taken as \(\mu \approx \kappa \sqrt{N_\mathrm{max}}\) in the basis representation. In these definitions, DAs are normalized to unity when one uses the same basis cutoff for the LFWFs and decay constants. We compare the DAs of four heavy-light systems in Fig. 5.
Noticeably, DAs share some features with the PDFs: the DA spreads wider along x in the lighter system as it is more relativistic; dips in \(2^1S_0\) and \(2^3S_1\) states reflect their character as radial excitations; the discrepancy between pseudoscalar and vector which is caused by the different spin configuration mixing appears more significant in lower exited states; and lastly, the asymmetries shown in the heavy-light meson DAs. However, different from PDFs, the DAs of the ground state (\(1^1S_0\)) peaked approximately at the quark mass fraction \(x=m_q/(m_q+m_{{\bar{q}}})\), which is in accordance with the other unequal mass system \(B_c\) [3]. Noticeable wiggles that appear in DAs especially of B and \(B_s\) are due to the limited range of basis spaces employed. Similar patterns are also found in the light system [18], and are understood to be resolved by increasing the basis size.
We find interesting similarities as well as differences of our DAs with those from the Dyson-Schwinger equations approach [30]. For example, the DAs for our D and \(D_s\) ground states are more widely spread in x. On the other hand, the peak heights of our DAs for B and \(B_s\) are very similar to the corresponding DAs in Ref. [30].
3.5 Decay constants
The decay constants provide important information of the internal structure of the mesons. \(f_P\) and \(f_V\) are the decay constants of the pseudoscalar and vector mesons, which characterize the strength of the matrix elements of the electroweak current between the corresponding meson and the vacuum. By definition, they are,
where \(p^\mu \) is the four-momentum of the meson, \(\epsilon ^\mu _{m_J}\) and \(M_V\) are the polarization vector and mass of the vector meson, respectively. These decay constants correspond to the LFWFs at the origin in coordinate space, and can be calculated by integrating the LFWFs as follows [2, 3],
We calculate the decay constants of the ground state pseudoscalar and vector of the heavy-light mesons, and compare with the experiments [31, 32] and other published approaches [33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56]. As mentioned earlier, the basis cutoff \(N_\mathrm{max}\) is associated with the UV regulator by \(\varLambda _\mathrm{UV} \approx \kappa \sqrt{N_\mathrm{max}}\). For calculating the decay constants, we employ different \(N_\mathrm{max}\) values for different systems so that \(\varLambda _\mathrm{UV} \approx 0.85 (m_q + m_{{\bar{q}}})\) in accordance with our previous work of heavy quarkonia [2]: we apply \(N_\mathrm{max} = 8\) for D and \(D_s\), while \(N_\mathrm{max} =16\) for B and \(B_s\). Uncertainties quoted for BLFQ are used to show the basis cutoff sensitivity, where we take \(\varDelta f =|{ f(N_\mathrm{max}=8) - f(N_\mathrm{max}=16)}|\) for all heavy-light systems. The choice of \(N_\mathrm{max}\) is a result of compromise between the needs for better basis resolution and a lower UV scale owing to the omitted radiative corrections.
The difference of the decay constants of pseudoscalar and vector heavy-light mesons is due to differences in the internal spin configurations. There is a significant tension among results obtained from different methods. Within the same approach, Lattice QCD for example, the tension still exists among the results reported by different groups [37,38,39]. In this work, \(f_P/f_V\) is greater than 1 for all heavy-light mesons, whereas the QCD sum rules (QCDSR) and relativistic quark model (RQM) quoted here favor a value less than 1. Details have been discussed in Refs. [38, 57,58,59]. In addition, we check the ratio for the light mesons \(f_{K}/f_{K^*}\) based on the same model, and obtain a result slightly smaller than 1, which agrees with experiments and also suggests this ratio is not systematically larger than 1 as seen in heavy systems. However, more work is needed to bring in additional physics, including the role of chiral symmetry [8, 19], before one can extend this approach to light systems [17, 18].
The ratios of the decay constants for the mesons which consist of only one different flavored quark is another interesting quantity. These quantities can reduce some of the systematic bias in theoretical and experimental approaches. The ratios \({\mathcal {R}}(f_{B_s}/f_{B})\) and \({\mathcal {R}}(f_{D_s}/f_{D})\) are presented in Fig. 6, with a comparison between results from other methods and from experiments. The vertical black lines indicate the central values of the ratio, and rectangles indicate the corresponding uncertainties. For the BLFQ results, the central lines are the ratios calculated with \(N_\mathrm{max} = 8\) for D and \(D_s\), \(N_\mathrm{max} =16\) for B and \(B_s\), respectively. Meanwhile the uncertainty is quoted as \(\varDelta _{\mathcal {R}} = |{{\mathcal {R}}(N_\mathrm{max}=8)-{\mathcal {R}}(N_\mathrm{max}=16)}|\), because \(\varDelta f\) is strongly correlated. Unlike the decay constant itself, the ratio does not have a monotonic dependence on the basis size. Among the several methods, ratios from BLFQ are smaller than those from other approaches by 5–18%. We mention that the decay constant of a pseudoscalar meson, taking the flavor symmetry breaking into account, is discussed in Ref. [60].
Comparison of the ratios \({\mathcal {R}}(f_{D_s}/f_{D})\) and \({\mathcal {R}}(f_{B_s}/f_{B})\) with other approaches [35, 47, 56] and experimental average [20]. See Table 2 for the specific values and their sources for each result plotted here. The vertical black lines indicate the central value of the ratios, while rectangles are the corresponding uncertainties
4 Summary and outlook
In this work, a model consisting of a holographic confining interaction, a longitudinal confinement and a one-gluon exchange is applied to mesons containing a light quark and is solved within the BLFQ method. The mass spectroscopy shows reasonable agreement with the experimental measurements and Lattice QCD calculations. This extends the previously successful application of the same model to heavy mesons. In this work, we fitted two parameters \(m_\mathsf {u/d}\) and \(m_{\mathsf {s}}\), and all those parameters such as \(m_{\mathsf {c}}\) and \(\kappa _{\mathsf {c}\bar{\mathsf {c}}}\) are universal among these systems. Without incorporating a chiral effective term in the Hamiltonian, we still obtain reasonable hyperfine splitting between the ground state pseudoscalar and vector, which suggests that the chiral dynamics is not dominant in the spectra of heavy-light mesons as much as in the pure light systems.
We investigate the LFWFs of the heavy-light systems and found non-trivial distribution of peaks along x. Due to this attribute of the LFWFs, the light-cone distributions such as PDF and DA, calculated by integrating the wave function, will be effected. Therefore we studied the PDFs and DAs of the heavy-light systems, and noticed that the peak of these functions, were not always locate at the quark mass fraction as one might naively expect. In particular, the peak of ground state DA is located close to the \(x=m_q/(m_q+m_{{\bar{q}}})\), but that of PDF is located well above \(x=m_q/(m_q+m_{{\bar{q}}})\).
We calculated the decay constants and the ratios which are also of great experimental and theoretical interest. The values of decay constant vary with the basis cutoffs due to the lack of radiative corrections in our model. Nevertheless, the ratio of heavier to lighter meson decay constants remains relatively stable. Even though the absence of chiral dynamics seems not to spoil the mass spectrum of heavy-light mesons, one could add an effective chiral potential in the Hamiltonian [18] in a future work, in order to address the decay constant ratio of vector to pseudoscalar more rigorously and investigate the role of chiral effects. In the mean time, incorporating the self-energy correction by including higher Fock sectors is another direction of improvement, which could reduce the sensitive dependence on basis size of the decay constant.
The fitted masses of light quarks of this work agree well with the fits from light mesons [17], both of which provide \(m_{\mathsf {u}}/m_{\mathsf {d}} \sim 550\) MeV in the BLFQ approach. The obtained LFWFs for a wide range of meson systems provide opportunities to explore a broader class of observables. Notable examples include the dilepton decays, diphoton decays, radiative transitions, and weak decays.
Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors’ comment: The authors plan to publish the light-front wave functions later in a data repository.]
Notes
We have tested our model with another confining strength \(\kappa = 0.54\) GeV provided by Ref. [16], and found the overall spectra were not significantly affected: the deviation is no more than \(20\%\) in r.m.s or \(6\%\) in \(\overline{\delta _J M} \) through the four heavy-light systems presented in Table 1.
Details of decay constant will be discussed in Sect. 3.5.
References
J.P. Vary, H. Honkanen, J. Li, P. Maris, S.J. Brodsky, A. Harindranath, G.F. de Téramond, P. Sternberg, E.G. Ng, C. Yang, Hamiltonian light-front field theory in a basis function approach. Phys. Rev. C 81, 035205 (2010)
Y. Li, P. Maris, J.P. Vary, Quarkonium as a relativistic bound state on the light front. Phys. Rev. D 96, 016022 (2017)
S. Tang, Y. Li, P. Maris, J.P. Vary, \({B}_{c}\) mesons and their properties on the light front. Phys. Rev. D 98, 114038 (2018)
N. Isgur, M.B. Wise, Weak decays of heavy mesons in the static quark approximation. Phys. Lett. B 232(1), 113–117 (1989)
M.A. Nowak, M. Rho, I. Zahed, Chiral effective action with heavy-quark symmetry. Phys. Rev. D 48, 4370–4374 (1993)
W.A. Bardeen, C.T. Hill, Chiral dynamics and heavy quark symmetry in a solvable toy field-theoretic model. Phys. Rev. D 49, 409–425 (1994)
M. Nielsen, F.S. Navarra, S.H. Lee, New charmonium states in QCD sum rules: a concise review. Phys. Rep. 497(2), 41–83 (2010)
P. Maris, P.C. Tandy, QCD modeling of hadron physics. Nucl. Phys. B Proc. Suppl. 161, 136–152 (2006). (Proceedings of the Cairns Topical Workshop on Light-Cone QCD and Nonperturbative Hadron Physics)
R. Aaij, B. Adeva, M. Adinolfi, A. Affolder, Z. Ajaltouni, S. Akar, J. Albrecht, F. Alessio, M. Alexander, S. Ali et al., Observation of overlapping Spin-1 and Spin-3 \(\overline{D}^{0}{K}^{-}\) resonances at mass \(2.86\,\,{\rm GeV} /{c}^{2}\). Phys. Rev. Lett. 113, 162001 (2014)
The LHCb collaboration, R. Aaij, B. Adeva, M. Adinolfi, A. Affolder, Z. Ajaltouni, S. Akar, J. Albrecht, F. Alessio, M. Alexander, S. Ali, and et. al. Precise measurements of the properties of the \({B}_1(5721)^{0,+}\) and \({B}_2^*(5747)^{0,+}\) states and observation of \({B}^{+,0}\), \(\pi ^{-,+}\) mass structures. J. High Energy Phys., 2015(4):24 (2015)
N. Cabibbo, Unitary symmetry and leptonic decays. Phys. Rev. Lett. 10, 531–533 (1963)
M. Kobayashi, T. Maskawa, CP-violation in the renormalizable theory of weak interaction. Progress Theoret. Phys. 49(2), 652–657 (1973)
S. Narison, \(|{V}_{cd}|\), \(|{V}_{cs}|\) and \(f_{D_{(s)}}\) from (semi-) leptonic \({D}_{(s)}\)-decays: Signals of new physics? Phys. Lett. B 668(4), 308–311 (2008)
G.F. de Téramond, S.J. Brodsky, Light-front holography: a first approximation to QCD. Phys. Rev. Lett. 102, 081601 (2009)
Y. Li, P. Maris, X. Zhao, J.P. Vary, Heavy quarkonium in a holographic basis. Phys. Lett. B 758, 118–124 (2016)
S.J. Brodsky, G.F. de Téramond, H.G. Dosch, J. Erlich, Light-front holographic QCD and emerging confinement. Phys. Rep. 584, 1–105 (2015)
W. Qian, S. Jia, Y. Li, J.P. Vary. Light mesons within the basis light-front quantization framework (In preparation)
S. Jia, J.P. Vary, Basis light front quantization for the charged light mesons with color singlet Nambu–Jona–Lasinio interactions. Phys. Rev. C 99, 035206 (2019)
P. Maris, C.D. Roberts, \(\pi \)- and \({K}\)-meson Bethe–Salpeter amplitudes. Phys. Rev. C 56, 3369–3383 (1997)
M. Tanabashi, K. Hagiwara, K. Hikasa, K. Nakamura, Y. Sumino, F. Takahashi, J. Tanaka, K. Agashe, G. Aielli, C. Amsler et al., Review of particle physics. Phys. Rev. D 98, 030001 (2018)
E.B. Gregory, C.T.H. Davies, I.D. Kendall, J. Koponen, K. Wong, E. Follana, E. Gámiz, G. Peter Lepage, E.H. Müller, H. Na, J. Shigemitsu, Precise \({B}\), \({B}_{s}\), and \({B}_{c}\) meson spectroscopy from full lattice QCD. Phys. Rev. D 83, 014506 (2011)
R.J. Dowdall, C.T.H. Davies, T.C. Hammant, R.R. Horgan, Precise heavy-light meson masses and hyperfine splittings from lattice QCD including charm quarks in the sea. Phys. Rev. D 86, 094510 (2012)
C.B. Lang, D. Mohler, S. Prelovsek, R.M. Woloshyn, Predicting positive parity \({B}_s\) mesons from lattice QCD. Phys. Lett. B 750, 17–21 (2015)
M. Kalinowski, M. Wagner, Masses of \({D}\) mesons, \({D}_{s}\) mesons, and charmonium states from twisted-mass lattice QCD. Phys. Rev. D 92, 094508 (2015)
E. Kou, P. Urquijo, W. Altmannshofer, F. Beaujean, G. Bell, M. Beneke, I.I. Bigi, F. Bishara M. Blanke, C. Bobeth, M. Bona, N. Brambilla, et al. The Belle II physics book (2018)
S. Jia, J.P. Vary. Parton distribution functions from scalar light front parton gas model (2018)
P. Maris, S. Jia, M. Li, Y. Li, S. Tang, J.P. Vary. On the light-front wave functions of quarkonia. In: Light Cone 2019 (LC2019) Palaiseau, France, September 16–20, 2019 (2020)
J. Lan, C. Mondal, S. Jia, X. Zhao, J.P. Vary, Parton distribution functions from a light front hamiltonian and QCD evolution for light mesons. Phys. Rev. Lett. 122, 172001 (2019)
G. Peter Lepage, S.J. Brodsky, Exclusive processes in perturbative quantum chromodynamics. Phys. Rev. D 22, 2157–2198 (1980)
D. Binosi, L. Chang, M. Ding, F. Gao, J. Papavassiliou, C.D. Roberts, Distribution amplitudes of heavy-light mesons. Phys. Lett. B 790, 257–262 (2019)
B.I. Eisenstein, I. Karliner, S. Mehrabyan, N. Lowrey, M. Selen, E.J. White, J. Wiss, R.E. Mitchell, M.R. Shepherd, D. Besson et al., Precision measurement of \(\cal{B}({D}^{+}\rightarrow {\mu }^{+}\nu )\) and the pseudoscalar decay constant \({f}_{{D}^{+}}\). Phys. Rev. D 78, 052003 (2008)
A. Zupanc, I. Adachi, H. Aihara, K. Arinstein, D.M. Asner, T. Aushev, A.M. Bakich, A. Bala, B. Bhuyan, G. Bonvicini et al., Measurements of branching fractions of leptonic and hadronic \({D}_s^+\) meson decays and extraction of the \({D}_s^+\) meson decay constant. J. High Energy Phys. 2013(9), 139 (2013)
A. Bazavov, C. Bernard, N. Brown, C. DeTar, A.X. El-Khadra, E. Gámiz, S. Gottlieb, U.M. Heller, J. Komijani, A.S. Kronfeld, J. Laiho, P.B. Mackenzie, E.T. Neil, J.N. Simone, R.L. Sugar, D. Toussaint, R.S. Van de Water, \({B}\)- and \({D}\)-meson leptonic decay constants from four-flavor lattice QCD. Phys. Rev. D 98, 074512 (2018)
H. Na, C.J. Monahan, C.T.H. Davies, R. Horgan, G. Peter Lepage, J. Shigemitsu, \({B}\) and \({B}_{s}\) meson decay constants from lattice QCD. Phys. Rev. D 86, 034506 (2012)
S. Aoki, Y. Aoki, D. Becirevic, T. Blum, G. Colangelo, S. Collins, M. Della Morte, P. Dimopoulos, S. Dürr, H. Fukaya, et al. Flag review 2019 (2019)
V. Lubicz, A. Melis, S. Simula, Masses and decay constants of \({D}_{(s)}^{*}\) and \({B}_{(s)}^{*}\) mesons with \({N}_{f}=2+1+1\) twisted mass fermions. Phys. Rev. D 96, 034524 (2017)
R.J. Dowdall, C.T.H. Davies, R.R. Horgan, C.J. Monahan, J. Shigemitsu, \({B}\)-meson decay constants from improved lattice nonrelativistic QCD with physical \(u\), \(d\), \(s\), and \(c\) quarks. Phys. Rev. Lett. 110, 222003 (2013)
B. Colquhoun, C.T.H. Davies, J. Kettle, J. Koponen, A.T. Lytle, R.J. Dowdall, G.P. Lepage, \({B}\)-meson decay constants: a more complete picture from full lattice QCD. Phys. Rev. D 91, 114509 (2015)
D. Becirevic, P. Boucaud, J.P. Leroy, V. Lubicz, G. Martinelli, F. Mescia, F. Rapuano, Nonperturbatively improved heavy-light mesons: masses and decay constants. Phys. Rev. D 60, 074501 (1999)
D. Bečirević, V. Lubicz, F. Sanfilippo, S. Simula, C. Tarantino, \({D}\)-meson decay constants and a check of factorization in non-leptonic \({B}\)-decays. J. High Energy Phys. 2012(2), 42 (2012)
P. Dimopoulos, R. Frezzotti, G. Herdoiza, V. Lubicz, C. Michael, D. Palao, G.C. Rossi, F. Sanfilippo, A. Shindler, S. Simula, C. Tarantino, M. Wagner, Lattice QCD determination of \(m_b\), \(f_{B}\) and \(f_{B_s}\) with twisted mass wilson fermions. J. High Energy Phys. 2012(1), 46 (2012)
N.H. Christ, J.M. Flynn, T. Izubuchi, T. Kawanai, C. Lehner, A. Soni, R.S. Van de Water, O. Witzel, \({B}\)-meson decay constants from \(2+1\)-flavor lattice QCD with domain-wall light quarks and relativistic heavy quarks. Phys. Rev. D 91, 054502 (2015)
A. Bazavov, C. Bernard, C.M. Bouchard, C. DeTar, M. Di Pierro, A.X. El-Khadra, R.T. Evans, E.D. Freeland, E. Gámiz, S. Gottlieb et al., \({B}\)- and \({D}\)-meson decay constants from three-flavor lattice QCD. Phys. Rev. D 85, 114506 (2012)
N. Carrasco, P. Dimopoulos, R. Frezzotti, P. Lami, V. Lubicz, F. Nazzaro, E. Picca, L. Riggio, G.C. Rossi, F. Sanfilippo, S. Simula, C. Tarantino, Leptonic decay constants \({f}_{K}\), \({f}_{D}\), and \({f}_{{D}_{s}}\) with \({N}_{f}=2+1+1\) twisted-mass lattice QCD. Phys. Rev. D 91, 054507 (2015)
W.-P. Chen, Y.-C. Chen, T.-W. Chiu, H.-Y. Chou, T.-S. Guu, T.-H. Hsieh, Decay constants of pseudoscalar \({D}\)-mesons in lattice QCD with domain-wall fermion. Phys. Lett. B 736, 231–236 (2014)
H. Na, C.T.H. Davies, E. Follana, G. Peter Lepage, J. Shigemitsu, \(|{V}_{cd}|\) from \({D}\) meson leptonic decays. Phys. Rev. D 86, 054510 (2012)
P. Gelhausen, A. Khodjamirian, A.A. Pivovarov, D. Rosenthal, Decay constants of heavy-light vector mesons from QCD sum rules. Phys. Rev. D 88, 014015 (2013)
W. Lucha, D. Melikhov, S. Simula, OPE, charm-quark mass, and decay constants of \({D}\) and \({D}_s\) mesons from QCD sum rules. Phys. Lett. B 701(1), 82–88 (2011)
W. Lucha, D. Melikhov, S. Simula, Decay constants of the charmed vector mesons \({D}^*\) and \({D}_s^*\) from QCD sum rules. Phys. Lett. B 735, 12–18 (2014)
S. Narison, A fresh look into \(\overline{m}_{c, b} (\overline{m}_{c, b}) \) and precise \(f_{D_{(s)}, B_{(s)}}\) from heavy-light QCD spectral sum rules. Phys. Lett. B 718(4), 1321–1333 (2013)
Z.-G. Wang, Analysis of the masses and decay constants of the heavy-light mesons with QCD sum rules. Eur. Phys. J. C 75(9), 427 (2015)
M.-Z. Yang, Wave functions and decay constants of \({B}\) and \({D}\) mesons in the relativistic potential model. Eur. Phys. J. C 72(2), 1880 (2012)
H.-K. Sun, M.-Z. Yang, Wave functions and leptonic decays of bottom mesons in the relativistic potential model. Phys. Rev. D 99, 093002 (2019)
D. Ebert, R.N. Faustov, V.O. Galkin, Relativistic treatment of the decay constants of light and heavy mesons. Phys. Lett. B 635(2), 93–99 (2006)
G. Cvetič, C.S. Kim, G.-L. Wang, W. Namgung, Decay constants of heavy meson of \(0^-\) state in relativistic Salpeter method. Phys. Lett. B 596(1), 84–89 (2004)
C.-W. Hwang, \({SU}(3)\) symmetry breaking in decay constants and electromagnetic properties of pseudoscalar heavy mesons. Phys. Rev. D 81, 054022 (2010)
D. Becirevic, A. Le Yaouanc, A. Oyanguren, P. Roudeau, F. Sanfilippo, Insight into \(D/B\rightarrow \pi \ell \nu _\ell \) decay using the pole models (2014)
W. Lucha, D. Melikhov, S. Simula, QCD sum-rule results for heavy-light meson decay constants and comparison with lattice QCD. In: 8th International Workshop on the CKM Unitarity Triangle (CKM 2014) Vienna, Austria, September 8–12, 2014 (2014)
S. Narison, Decay constants of heavy-light mesons from QCD. Nucl. Particle Phys. Proc. 270–272, 143–153 (2016). [18th Montpellier International Conference on Quantum Chromodynamics (QCD 15)]
S.S. Gershtein, M.Y. Khlopov., Leptonic decays of heavy pseudoscalar mesons. JETP Lett. 23(6), 338 (1976)
Acknowledgements
We wish to thank S. Jia, M. Li, W. Qian, and A. Yu for fruitful discussions. S.T. and J.P.V. wish to thank Dr. Soeren Prell for valuable advice. This work was supported in part by the U.S. Department of Energy under Grants Nos. DE-FG02-87ER40371 and DE-SC0018223 (SciDAC-4/NUCLEI). Computational resources were provided by the National Energy Research Supercomputer Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Funded by SCOAP3
About this article
Cite this article
Tang, S., Li, Y., Maris, P. et al. Heavy-light mesons on the light front. Eur. Phys. J. C 80, 522 (2020). https://doi.org/10.1140/epjc/s10052-020-8081-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-020-8081-9