Abstract

An analytic solution of the N-dimensional radial Schrödinger equation with the combination of vector and scalar potentials via the Laplace transformation method (LTM) is derived. The current potential is extended to encompass the spin hyperfine, spin-orbit, and tensor interactions. The energy eigenvalues and the corresponding eigenfunctions have been obtained in the N-dimensional space. The present results are employed to study the different properties of the heavy-light mesons (HLM). The masses of the scalar, vector, pseudoscalar, and pseudovector for B, Bs, D, and Ds mesons have been calculated in the three-dimensional space. The effect of the dimensional number space is discussed on the masses of the HLM. We observed that the meson mass increases with increasing dimensional space. The decay constants of the pseudoscalar and vector mesons have been computed. In addition, the leptonic decay widths and branching ratio for the B+, D+, and mesons have been studied. Therefore, the used method with the current potential gives good results which are in good agreement with experimental data and are improved in comparison with recent theoretical studies.

1. Introduction

One of the most important tasks in nonrelativistic quantum mechanics is to get the solution of the Schrödinger equation. The solution of the Schrödinger equation with spherically symmetric potentials plays a significant role in many fields of physics such as hadronic spectroscopy for understanding the quantum chromodynamics theory. Numerous works have been introduced to get the solution of Schrödinger equation using different methods like the operator algebraic method [1], path integral method [2], the conventional series solution method [3, 4], Fourier transform [5, 6], shifted () expansion [7, 8], point canonical transformation [9], quasi-linearization method [10], supersymmetric quantum mechanics (SUSQM) [11], Hill determinant method (HDM) [12], and other numerical methods [1315].

Recently, the study of the different topics has received a great attention from theoretical physicists in the higher dimensional space. In addition, the study is more general and one can obtain the required results in the lower dimensions directly, such as the hydrogen atom [1618], harmonic oscillator [19, 20], random walks [21], Casimir effects [22], and the quantization of angular momentum [2327]. The -dimensional Schrödinger equation has been studied with different forms of spherically symmetric potentials [2833]. The -dimensional Schrödinger equation has been investigated with the Cornell potential and extended Cornell potential [3438] using different methods such as the Nikiforov-Uvarov (NU) method [32, 36, 39, 40], power series technique (PST) [41], the asymptotic iteration method (AIM) [34], Pekeris type approximation (PTA) [41, 42], and the analytical exact iteration method (AEIM) [43, 44].

The LTM is one of the useful methods that contributed to finding the exact solution of Schrödinger equation in one-dimensional space for Morse potential [45, 46], the harmonic oscillator [47], and three-dimensional space with pseudoharmonic and Mie-type potentials [48] and with noncentral potential [49]. The N-dimensional Schrödinger equation has been solved via the LTM in many studies for Coulomb potential [28], harmonic oscillator [50], Morse potential [51], pseudoharmonic potential [52], Mie-type potential [53], anharmonic oscillator [54], and generalized Cornell potential [38].

The study of different properties of HLM is very vital for understanding the structure of hadrons and dynamics of heavy quarks. Thus, many theoretical and experimental efforts have been done for understanding distinct characteristics of HLM. In [4, 34, 55], the authors calculated the mass spectra of quarkonium systems as charmonium and bottomonium mesons with the quark-antiquark interaction potential using various methods in many works. Al-Jamel and Widyan [56] studied the spin-averaged mass spectra of heavy quarkonia with Coulomb plus quadratic potential using (NU) method. Abou-Salem [57] has computed the masses and leptonic decay widths of , , , , , and numerically using Jacobi method. The strong decays, spectroscopy, and radiative transition of heavy-light hadrons have been computed using the quark model predictions [58]. The decay constant of HLM has been calculated using the field correlation method [59]. Moreover, the spectroscopy of HLM has been investigated in the framework of the QCD relativistic quark model [60]. The spectroscopy and Regge trajectories of HLM have been obtained using quasi-potential approach [61]. The decay constants of heavy-light vector mesons [62] and heavy-light pseudoscalar mesons [63] have been calculated with QCD sum rules. A comparative study has been introduced for the mass spectrum and decay properties for the D meson with the quark-antiquark potential using hydrogeometric and Gaussian wave function [64]. In framework of Dirac formalism the mass spectra of Ds [65] and D [66] mesons have been obtained using Martin-light potential in which the hadronic and leptonic decays of D and Ds mesons have been evaluated [67]; besides the rare decays of and mesons into dimuon (μ+μ-) [68] and the decay constants of B and Bs have been calculated [69]. The mass spectra and decay constants for ground state of pseudoscalar and vector mesons have been obtained using the variational analysis in the light quark model [70]. The spectroscopy of bottomonium and B meson has been studied using the free-form smearing in [71]. The variational method has been employed to compute the masses and decay constants of HLM in [72]. In addition, the decay properties of D and Ds mesons have been investigated using the quark-antiquark potential in [73]. The B and mesons spectra and their decays have been studied with a Coulomb plus exponential type potential in [74]. The leptonic and semileptonic decays of B meson into τ have been studied [75]. The degeneracy of HLM with the same orbital angular momentum has been broken with the spin-orbit interactions [76]. The relativistic quark model has been investigated to study the properties of B and mesons [77] and the excited charm and charm-strange mesons [78]. The perturbation method has been employed to determine the mass spectrum and decay properties of HLM with the mixture of harmonic and Yukawa-type potentials [79]. In [80], the authors have investigated the leptonic decays of seven types of heavy vector and pseudoscalar mesons. The spectra and wave functions of HLM have been calculated within a relativistic quark model by using the Foldy-Wouthuysen transformation [81]. The isospin breaking of heavy meson decay constants had been compared with lattice QCD from QCD sum rules [82]. The decay constants of pseudoscalar and vector B and D mesons have been studied in the light-cone quark model with the variational method [83]. In [84], the authors have calculated the strong decays of newly observed (3000) and (3040) with two 2P (1+) quantum number assignments. The leptonic and semileptonic decays have been analyzed using the covariant quark model with infrared confinement within the standard model framework [85]. The weak decays of B, , and into -wave heavy-light mesons have been studied using Bethe-Salpeter equation [86]. In [87], the decay constant and distribution amplitude for the heavy-light pseudoscalar mesons have been evaluated using the light-front holographic wavefunction. By using the Gaussian wave function with quark-antiquark potential model, the Regge trajectories, spectroscopy, and decay properties have been studied for and mesons [88], D and Ds mesons [89], and also the radiative transitions and the mixing parameters of the -meson have been obtained [90]. The dimensional space dependence of the masses of heavy-light mesons has been investigated using the string inspired potential model [91].

The goal of this work is to get the analytic solution of the N-dimensional Schrödinger equation for the mixture of vector and scalar potentials including the spin-spin, spin-orbit, and tensor interactions using LTM in order to obtain the energy eigenvalues in the N-dimensional space and the corresponding eigenfunctions. So far no attempt has been made to solve the N-dimensional Schrödinger equation using the LTM when the spin hyperfine, spin-orbit, and tensor interactions are included. To show the importance of present results, the present results are employed to calculate the mass spectra of the HLM in three-dimensional space and in the higher dimensional space. In addition, the decay constants, leptonic decay widths, and branching fractions of the HLM are calculated.

The paper is systemized as follows: the contributions of previous works are displayed in Section 1. In Section 2, a brief summary of Laplace transformation method is introduced. In Section 3, an analytic solution of the N-dimensional Schrödinger equation is derived. In Section 4, the obtained results are discussed. In Section 5, summary and conclusion are presented.

2. Overview of Laplace Transformation Method

The Laplace transform or of a function is defined by [92]If there is some constant such that for sufficiently large , the integral in (1) exists for Re for . The Laplace transform may fail to exist because of a sufficiently strong singularity in the function as . In particularwhere is the gamma function. The Laplace transform has the derivative propertieswhere the superscript stands for the -th derivative with respect to for and with respect to for . If is the singular point, the Laplace transform behaves near as and then for On the other hand, if near origin behaves like with , then behaves near as

3. Analytic Solution of the -Dimensional Radial Schrödinger Equation

The N-dimensional radial Schrödinger equation that describes the interaction between quark-antiquark systems takes the form [41]where represent the angular quantum number and the dimensional number greater than one, respectively, and is the reduced mass of the quark-antiquark system.

In the nonrelativistic quark model, the quark-antiquark potential consists of the spin independent potential and the spin dependent potential , respectively:The spin independent potential is taken as a combination of vector and scalar parts [93]: where and are the vector and scalar parts, respectively, and stands for the mixing coefficient. ɑ, b, and c are arbitrary parameters where a, b, and c > 0 which are fitted with experimental data. The harmonic and linear terms represent the confining part at long distance and the Coulomb term stands for the quark-antiquark interactions through one gluon exchange at short distances which gives better description of quark-antiquark interaction.

The spin dependent potential is extended to three types of interaction terms as [94]while the spin-orbit and tensor terms give the fine structure of the states, the spin-spin interaction term describes the hyperfine splitting of the state, and is an angular quantum operator, and is a spin operator (for detail, see [94]).where is radial Laplace operator. The diagonal elements of the are defined. Substituting (11)-(16) into (9) then the nonrelativistic quark-antiquark potential takes the formwhereSubstituting (21) into (8), then whereThe complete solution of (25) takes the formwhere the term confirms that the solution is bounded at . The function is yet to be determined. From (27) we getSubstituting (27), (28), and (29) into (25), then, whereIn order to apply the Laplace transform of the above differential equation, the parametric condition is taken as in [52, 54]. Thus, (32) has a solution We take the physical solution of (32) as in [52, 54].

Substituting (34) into (30) yieldsBy expanding the term around , where and is a parameter as in [36, 56], we get Substituting (37) into (36) yieldswhereThe Laplace transform is defined as and taking boundary condition yields HereThe singular point of (40) is By using the condition of (5), the solution of (40) takes the form

From (42), Substituting (42)-(44) into (40), we obtain the following relations:Using (26), (39), and (41) and the set of (45)-(47), then, the energy eigenvalue of (8) in the N-dimensional space is given by the relationTake the inverse Laplace transform such that The function takes the following form:Using (11), (13), and (23), the eigenfunctions of (9) take the following form:From the condition , the normalization constant can be computed. In addition, the wave equation satisfies the boundary condition .

4. Discussion of Results

In Figure 1, the current potential has been plotted in comparison to other potential models; we see that the present potential is in a qualitative agreement with other potential models [72, 74, 79], in which the confining part is clearly obtained in comparison to Cornell and Coulomb plus exponential potentials. The different states of and mesons have been shown in Figures 2 and 3, respectively, in which the principal number of states plays an important role in confining part of potential.

In the following subsections, we employ the obtained results in the previous section to determine the mass spectra of scalar, vector, pseudoscalar, and pseudovector of , , , and mesons in the -dimensional space in comparison with the experimental data (PDG 2016) [95] and with other recent studies. In addition, the decay properties such as decay constants, leptonic decay width, and the branching ratio of HLM are calculated.

4.1. Mass Spectra of Heavy-Light Mesons

The masses of HLM in the N-dimensional space are defined [44]: Substituting (48) into (51), then the mass spectra of HLM in the N- dimensional space can be found from the relationIn Tables 2–6, we have calculated the masses of the HLM in the three-dimensional space in comparison with the experimental data and other recent studies [7274, 81, 88, 89, 96]. The parameters used in the present calculations are shown in Table 1. In addition, the masses at N = 4 and N = 5 are calculated. In Tables 2 and 3, we observe that and meson masses close to experimental data and other meson masses are in good agreement with experimental data and become better in comparison to the results in recent studies [7274, 81, 88, 89, 96]. In comparison with [72], they used the variational method for the Cornell potential to study the HLM with including the spin-spin and spin-orbit interactions. They ignored the tensor interactions in their calculations. The present results are good in comparison to the results in [72]. In addition, we used the LTM in the present calculations. Yazarloo and Mehiraban used the variational method to study D and mesons for the Cornell potential [73] and used the Nikiforov-Uvarov (NU) method to study B and mesons for the Coulomb plus exponential type potential [74]. The present results are in good agreement with the results of [73, 74]. Kher et al. [89] used a Gaussian wave function to calculate the mass spectra of D and in addition to B and mesons [88] for the Cornell potential. Jing-Bin [81, 96] obtained the spectra of the HLM in the relativistic model from the Bethe-Salpeter equation using the Foldy-Wouthuysen transformation in his works.

We note that the present results for D and meson masses become better in comparison to the results of [81, 88, 89, 96], where the values of pseudoscalar D and mesons are close to the experimental data in Table 2. The values of vector and mesons close to the experimental data and the values of vector Ds and mesons are good in comparison to the experimental results in Table 3.

The masses of the scalar mesons are presented in Table 4; the value of D meson is close to the experimental value. The values of Ds and B are in agreement with the experimental values and the value of meson is in good agreement with the theoretical studies [7274, 81, 88, 89, 96]. In Table 5, we observe that all the values of pseudovector mesons are close to the experimental results except the value of B meson which is in good agreement with the experimental value. The values of vector Ds and B mesons are in good agreement with the experimental results. In Table 6, the results of the p-wave state for the HLM are reported.

The present predictions of D, , B, and Bs mesons are in agreement in comparison to the experimental data and the theoretical studies [73, 74, 81, 88, 89, 96].

In addition, we have investigated the masses of the HLM in the higher dimensions at N=4 and N=5. In Tables 2–6, the effect of the dimensional number is investigated on the masses of the HLM. One can see that the masses increase with increasing dimensional number. The influence of the dimensional number is not considered on the masses of the HLM in the works [7274, 81, 88, 89, 96]. Roy and Choudhury [91] have studied the masses of heavy flavor mesons in the higher dimensional space using string inspired potential. They found that an increase of the dimensional number leads to increase the meson masses. Therefore, the present results of the mass spectra of HLM are in good agreement in comparison with the results of [91].

4.2. Decay Constants

The study of the decay constants is one of the very significant characteristics of the HLM, as it provides a direct source of information on the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. Many theoretical studies have been done for determining the decay constants with different models as relativistic quark model [9799], lattice QCD [100102], QCD sum rules [62, 97, 103], and nonrelativistic model [7274, 79, 97].

The Van Royen-Weisskopf formula [104] can be used to calculate the decay constants of the pseudoscalar and vector mesons, respectively, in the nonrelativistic limit which is defined asThe Van Royen-Weisskopf formula with the QCD radiative corrections taken into account can be written as [105]whereand and , for pseudoscalar and vector mesons, respectively.

In Tables 7 and 8, we have determined the decay constants of the pseudoscalar and vector B and D mesons obtained from (53) and (54) in comparison with the results of other recent works. In [87], the authors evaluated the decay constant for the heavy-light pseudoscalar mesons using the helicity-improved light-front holographic wavefunction. In [83], the authors applied the variational method to study the decay constants of the pseudoscalar and vector B and D mesons in the light-cone quark model for the relativistic Hamiltonian with the Gaussian-type function.

In [72], the authors used the variational method to compute the decay constants of HLM from the radial Schrödinger equation with the Cornell potential. Zhi-Gang Wang [97] introduced an analysis of the decay constants of HLM with QCD sum rules. Yazarloo and Mehiraban [79] used the perturbation method to study the decay constants of , , B, and mesons with the combination of harmonic and Yukawa-type potentials.

In Table 7, the obtained results are in good agreement in comparison to the results of [72, 83, 87, 97]. In Table 8, the present results are compatible with the results of [73, 74, 79, 83]. In addition, the ratio of decay constants for D mesons is . This value is in good agreement with the experimental value [95]. The present result is in agreement with the obtained values in [72] and in [83]. Also, we have which is in agreement with the calculated values in [87] and in [97]. The calculated ratio of decay constants for B mesons and are in good agreement in comparison with and in [83].

4.3. Leptonic Decay Widths and Branching Ratio

The charged HLM can decay to a charged lepton pair via a virtual W± boson. The leptonic decay widths of the HLM can be obtained from the relation [106]where is the Fermi constant and the relevant CKM elements are taken from the PDG [95] as , ,  and   The leptonic masses are taken as , , and We obtain the decay constants of the HLM from Tables 7 and 8 into (56) to compute leptonic decay widths of the HLM. The obtained results of the leptonic decay width of B+, D+, and mesons are shown in Tables 9, 10, and 11, respectively. Vinodkumar et al. [107] calculated the leptonic decay widths of B, Bs mesons besides, D and Ds mesons [66, 67, 108] for the Martin-like potential with Dirac formalism. We have determined the leptonic decay widths of B+ meson in Table 9 in comparison with the results of the [74, 79, 107], as well as the leptonic decay widths of D+ meson in Table 10 in comparison with the results of [66, 79, 108] and the leptonic decay widths of meson in Table 11 compared with the results of [66, 79, 108]. We note that the present results are in good agreement with the results of [66, 67, 74, 107, 108].

The branching ratio of the HLM is defined aswhere the lifetime of B+, D+, and mesons is taken as , , and [95]. We have determined the branching ratio for the , , and mesons compared with the experimental data and with the results of other recent studies [7274, 88, 89].

In Table 12, we note that the present values of the branching ratio for the B+ meson are close to experimental data and are in agreement in comparison with the theoretical results [72, 74, 79, 88, 107]. In addition, in Tables 13 and 14, we note that the evaluated results of branching ratio for the D+ and mesons are close to the experimental data and become better in comparison with works [72, 73, 79, 89, 108].

5. Summary and Conclusion

In this work, we have presented an approximate-analytic solution of the N-dimensional radial Schrödinger equation for the mixture of vector and scalar potentials via the LTM. The spin-spin, spin-orbit, and tensor interactions have been included in the extended Cornell potential model. The energy eigenvalues and the corresponding eigenfunctions have been determined in the N-dimensional space. In three-dimensional space, we have employed the obtained results to study the different properties of the HLM that are not considered in many recent studies. The masses of the scalar, vector, pseudoscalar, and pseudovector for B, Bs, D, and Ds mesons have been calculated in the three-dimensional space and in the higher dimensional space in Tables 2–6. Most of the present calculations are close to the experimental data and are improved in comparison with the recent calculations [7274, 81, 88, 89, 96]. As well, we have computed the masses of the HLM in the higher dimensional space at N=4 and N=5. The dependence of the masses of HLM on the dimensional number is discussed. We found that the masses increase with increasing dimensional number. This result is obtained in [91]. In Tables 7 and 8, the decay constants of the pseudoscalar and vector mesons have been determined in comparison with the results of [7274, 79, 83, 87, 97]. The calculated ratios of the decay constants of D mesons and are close to the experimental ratio

The present results of the decay ratio of B mesons are in good agreement with the results of [72, 83]. The leptonic decay widths of B+ meson have been studied in comparison with the results of [74, 79, 107] and the leptonic decay widths of D+ meson in comparison with the results of [66, 79, 108]. In addition, the leptonic decay widths of meson have been studied in comparison with the results of [66, 79, 108].

The obtained results of the leptonic decay widths are compared with the results of [66, 67, 74, 107, 108]. We have determined the branching ratio for the B+, D+, and mesons that are in good agreement with the experimental data and with the recent studies [7274, 88, 89]. Therefore, the current potential with used method gives very good predictions for the heavy-light meson properties. We hope to extend this work to include external force as a future work.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.