Heavy-light mesons on the light front

We study the heavy-light mesons within basis light-front quantization. The resulting mass spectra of D, Ds\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D_s$$\end{document}, B, and Bs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_s$$\end{document} 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 (Ds\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D_s$$\end{document} to D and Bs\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_s$$\end{document} 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%.


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 a e-mail: leeyoung@iastate.edu (corresponding author) 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 heavylight 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) ± 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.

Basis light-front quantization for heavy-light systems
In this work, we adopt the effective Hamiltonian within the |qq Fock sector in the form introduced for the heavy mesons [2,3]. It comprises two parts 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 gen-erated by the effective one-gluon-exchange potential V eff g . Along with a running coupling α s (Q 2 ), V eff g governs the short-range physics: With the Hamiltonian, one can solve the light-front eigenvalue problem H tot |Ψ h (P, J, m J ) = M 2 |Ψ h (P, J, m J ) 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 ⊥ ≡ |k ⊥ |. The transverse basis functions, are the 2D harmonic oscillator functions where n and m are the principal and orbital quantum numbers, respectively; , with x being the longitudinal momentum fraction carried by the quark x = p + q /( p + q + p + q ); θ q = arg q ⊥ , b sets the scale of the basis, and L |m| n is the associated Laguerre polynomial. The longitudinal basis functions are given by where P (α,β) l (2x − 1) is the Jacobi polynomial with quantum number l; α and β 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 max , the ultraviolet (UV) and infrared (IR) regulators are implicit through Λ UV ≈ b √ N max , and Λ IR ≈ b/ √ N max , respectively. L max represents the basis resolution in the longitudinal direction. In this work, we choose N max = L max for simplicity. Notice that the total angular momentum projection is defined by m J = m + s +s, which is conserved in the constructed basis. Solving the eigenequation by diagonalizing the H tot matrix, one obtains the eigenvalues which indicate the spectra as squared masses, and the eigenvectors which are the coefficients ψ h (n, m, l, s,s) in Eq. (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 = κ. 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 (c) and bottom (b) to be the quarks, while setting the lighter quarks, up (u), down (d), and strange (s) to be the antiquarks in the heavylight systems. We take the masses of u and 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 κ xȳ = (κ 2 xx + κ 2 yȳ )/2. We adopt κ uū/dd = κ ss = 0.59 GeV 1 as the confining strength for the light mesons [16], while the values for heavy quarkonia, κ bb and κ cc , are taken from Ref. [2]. Following our previous works on heavy systems [2,3], we keep the same heavy quark masses (m b , m c ), while we fit the light-quark masses (m u /m d , m s ) to reproduce the heavylight meson mass spectra from experiments.
More specifically, we tune only the m u /m d (m s ) to minimize the r.m.s. mass deviation between experiments and theory for the lowest pseudoscalar and vector states, D 0 , The m u /m 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. η c ) and light mesons (e.g. ρ, π) 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 interac-tion dramatically improves the predictions for ρ, π masses and additional observables.

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: 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, δ 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 κ when the basis size is fixed by N max and L max . Specifically, smaller κ 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 heavylight 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].

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 . For the sake of convenience, we visualize the LFWFs at k y = 0, i.e. θ q = 0 or π in the phase exp(imθ 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 + mq ) where m q = mq . 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 ⊥ = 0 the asymmetry is maximal, and the peak is located at x > m c(b) /(m c(b) + m u ), but as k ⊥ increases the peak location (along fixed k ⊥ ) shifts to smaller values of x. In the limit k ⊥ → ∞ the peak location approaches x = 1/2. This can easily be understood because for k ⊥ m c(b) + m 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 Table 1 Summary of the model parameters with the basis truncation N max = L max = 32. Among them m q (mq ) = 0.553 GeV and 0.647 GeV are the two fitted parameters. The calculated meson masses of the two low-lying states, known as pseudoscalar (PS) and vector (V), are listed in the table. The r.m.s. are the root-mean-square differences of our results from the experimental measured masses; the number of compared states N exp are shown in the square brackets. The mean spread δ J M is the spread in the masses over allowed M J values and the devia-tion from zero reflects the violation of rotational symmetry. Following the absolute mean spread value, we provide the relative spread with respect to the total mass of constituent (anti-) quarks in the parenthesis. Note that the masses of light quarks listed in the table are fitted from heavy-light mesons, and applied to light mesons for supplementary calculations. Hence we do not provide the r.m.s. of the light mesons due to the absence of chiral dynamics in this model, which is outside the scope of this work as the parton distribution function (PDF) and distribution amplitude (DA) as we will see in the following sections.

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; μ). The PDF within collinear factorization defines the probability for finding a quark carrying a longitudinal momentum fraction x at resolution scale μ. In the LFWF representation, the PDF is calculated simply by integrating out the transverse momentum of the square of the wave function modulus, (11) and the truncation parameter N max of the LFWFs provides a UV scale for model μ ≈ κ √ N 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 1 0 f (x; μ)dx = 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 ⊥ in the LFWFs. However, we find that the peak in the PDF for the radial excited states 2 1 S 0 and 2 3 S 1 is at a significantly smaller value of x than that of 1 1 S 0 and 1 3 S 1 , and in fact it is very close to x = m q /(m q + mq ). The latter may be coincidental. The width of PDFs decreases as the mass of the system increases, and approaches a δ function in the nonrelativistic limit. The 2 1 S 0 and 2 3 S 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 1 S 0 and 1 3 S 1 states are larger than between 2 1 S 0 and 2 3 S 1 , that is due to the significant hyperfine splittings in lower excited states.
Note that we use the basis size as N max = L max = 32, which corresponds to different resolution scales: μ D/D s ≈ 4.5 GeV and μ B/B s ≈ 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.

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  [20]. Specifically, we use the values of D 0 (cū) and B ± (bū/bd) 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 elements. In the LFWF representation, it can be written as [29] f P, where f P(V ) are the decay constants 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 μ taken as μ ≈ κ √ N max in the basis representation. In these definitions, DAs are normalized to unity when one uses the same basis cutoff for the LFWFs 2 Details of decay constant will be discussed in Sect. 3.5. and decay constants. We compare the DAs of four heavylight 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 1 S 0 and 2 3 S 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 1 S 0 ) peaked approximately at the quark mass fraction x = m q /(m q +mq ), 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].

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 μ is the four-momentum of the meson, μ 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 . As mentioned earlier, the basis cutoff N max is associated with the UV regulator by Λ UV ≈ κ √ N max . For calculating the decay constants, we employ different N max values for different systems so that Λ UV ≈ 0.85(m q + mq ) in accordance with our previous work of heavy quarkonia  Red dots indicate the local maxima along x for fixed k ⊥ , 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 + mq ) when k ⊥ = 0, and it shifts to smaller x as k ⊥ increases and finally approaches to x = 1/2 at the limit k ⊥ → ∞ 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 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 max = 8 for D and D s , N max = 16 for B and B s , respectively. Meanwhile the uncertainty is quoted as Δ R = |R(N max = 8) − R(N max = 16)|, because Δ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].

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 u/d and m s , and all those parameters such as m c and κ cc 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 Table 2 The decay constant of ground state pseudoscalar and vector states of four heavy-light systems. We list our results (BLFQ) and compare to existing experimental results, and other theoretical  (20) 313 (67) 306 (39) 233 (50) 202 (32) 259 (54) 230 (36) (10) 263 (21) 240 (10) 308 (21) 194 (15) 213 (18) 231 (16) 255 (19) [51] RQM … … … … 210 (10) 223 (16) 229 (11) 242 (17 [35,47,56] and experimental average [20]. See Table 2 for the specific values and their sources for each result plotted here. The ver-tical black lines indicate the central value of the ratios, while rectangles are the corresponding uncertainties 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 + mq ), but that of PDF is located well above x = m q /(m q + mq ).
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 selfenergy 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 u /m d ∼ 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. were provided by the National Energy Research Supercomputer Cen-ter, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231.

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.] 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://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .