The Structure Functions in the Weak and Semileptonic Decay of Meson B to the Mesons D and D*: A Computational Study ()
1. Introduction
In this paper, we study the Beth-Salpeter equation for the bound state of two-piont particles with a confining kernel. Accordingly, we achieve the wave function of heavy meson in the framework of heavy quark effective theory assumptions. We compute the semileptonic weak decay of B to D and D* using the BS wave function to give q2- dependence of the form factors. The obtained results for the form factors, their slopes and curvature are consistent with the empirical data. The QCD as the strong-interactions theory has succeeded in describing the high-energy physics. The discovery of asymptotic freedom has made it possible to perturbationally calculate a substantial number of physical quantities. Nevertheless, the QCD behaves inappropriately at the infrared region. For instance, there is not a tool to describe the enclosure of quarks so that either the lattice gauge theory should be utilized to describe it or the QCD should be utilized to solve it through creating effective theories.
Nowadays, the weak decays of hadrons including heavy quarks are used for experimenting on the standard model and measuring its parameters. Especially, these decays are an appropriate tool to obtain the CKM matrix’s entries. Considering the QCD’s problems, the employment of completely independent ways of investigating such decays is of great significance [1-3].
To acquire the zero-order wave function, we conform to the references [1-3] that are based on the zero-order HQET. More, a confining kernel [4-6] is introduced for the Bethe-Salpeter equation; through utilizing this kernel, the first-order Bethe-Salpeter wave functions are achieved. Then, the structure functions of and transitions are obtained by means of the zeroand first-order Bethe-Salpeter wave functions.
2. The Bethe-Salpeter Equation and Confining Kernel
In the momentum space, the Bethe-Salpeter equation for the two-fermion bound state (quarks) is as follows:
(1)
where is the amplitude of B-S for the bound state with the momentum P, is a kernel of B-S equation, and the s are the free particle propagators. In general, the kernel K is a sum of all irreducible renormalized diagrams representing the Green’s function with two incident fields and two outgoing fields.
Therefore, if the amplitude of the zero-order wave function (and the B-S equation’s kernel are available, the bound-state wave function can be achieved at every desirable order.
The references [1-3], that are based on the zero-order HQET, are conformed to obtain the zero-order wave function. The heavy quark is free and moves as fast as hadrons; consequently, it can be written as the form below:
(2)
where P and M are respectively the meson’s momentum and mass, and are the heavy meson’s momentum and mass, respectively. When the heavy quark is free, the quark spin can be generally indicated in a fourdimensional space of Dirac indexes. For example, a heavy meson is demonstrated by a two-index wave function (). However, when the heavy quark moves as fast as meson, is true for the Bergman-Wigner equation.
(3)
The above equation brings about the following quasiscalar wave function:
(4)
where is a matrix indicating the small degrees of freedom (soft gluons and light quark) that is a Lorentz scalar matrix function of light quark momentum. Herein, we choose A(k) as follows:
(5)
where a is a normalizing constant. k and m are the light quark’s momentum and mass, respectively. Hence, the zero-order wave function is written as the form below:
(6.1)
(6.2)
Now, to obtain the first-order wave function, it is sufficient to have the kernel of. In the reference [4-6], through taking a confining kernel, and after the Wick rotation, they achieved for the B-S equations the scalar kernel superposition in the momentum space:
(7)
and applying the below conditions:
(8)
A kernel is obtained as the following form:
(9)
into which a new effective coupling constant (U), that has mass, is inserted to be simplified. The Fourier transform of the kernel of Equation (9) is a constant in the fourdimensional space (). This means that it does not behave same as a non-relativistic constant potential for which the like behavior is expected. Therefore, the kernel of Equation (9) is not exactly similar to the nonrelativistic constant potential, and the effective constant () does not relate to a constant of such potential.
By the impact of the kernel of Equation (9) on the B-S equation of Equation (1), and under the Wick rotation, it can be written as:
(10)
3. Results and Discussion
3.1. The Computation of Structure Functions Using the Zero-Order Bethe-Salpeter
In this paper, for the first time, we proposed a novel analytical approach toward the structure functions in the weak and semileptonic decay of meson B to the mesons D and D*. The matrix of the flow producing the transition of mesons is as follows:
(11)
where is the ith meson’s wave function that can be calculated at every desirable order by virtue of Equations (1) or (10) and (6); expresses the weak-vertex flow, and is proportional to the CKM matrix, and equals in which is the corresponding entry to this transition. Accordingly, the matrix of flow for transitions can be written as follows:
(12.1)
(12.2)
The integration of k is logarithmically divergent; therefore, it is regulated through the dimensional regularizetion. Performing this, and calculating the Trs, these equations can be easily written as the following forms:
(13.1)
(13.2)
In the normalization, the term of is removed by means of the Luke’s theorem, and it is not necessary to compute it. Hence, after the normalization, the structure functions at this state equal:
(14)
This conclusion is predictable since the wave functions we chose were the static wave functions of extremely heavy quark; and as the HQET predicts for this state, the structure functions fulfill Equation (14) at the maximum recoil [7-14].
3.2. The Computation of Structure Functions Using the First-Order Bethe-Salpeter Wave Function
The zero-order wave functions (Equations (6.1) and (6.2)) are put into the first-order B-S equation with a confining kernel (Equation (10)) to obtain the first-order wave function () and actually to gain better results. Consequently, the first-order wave functions for mesons are as follows:
(15.1)
(15.2)
Through replacing the above equations with the matrixes of flow producing the transitions and with regards to Equation (11), these matrixes are as follows:
(16.1)
(16.2)
assuming that the heavy quarks are mass-shell:
(17)
and through algebra the γs, they are simplified as the forms below:
(18.1)
(18.2)
If the Feynman parameterization is applied:
(19)
where the s are the first and second (or third) fractions’ dominator. The s and s are defined as follows:
(20)
First, each of rejections in the integrand’s numerators (Equations (18.1) and (18.2)) is calculated. The numerator of the first fraction of Equation (18.1) equals:
(21.1)
and the numerator of the first fraction of Equation (18.2) equals:
(21.2)
which through inputting to the above equations, and removing the odd terms in proportion to, they are respectively simplified as follows:
(22.1)
(22.2)
The second fractions’ numerators are easily achieved as the following forms:
(23.1)
(23.2)
and the third fractions’ numerators are calculated as follows:
(24.1)
(24.2)
At the third step, the change of variable is applied to both above equations. Also, writing the odd terms with respect to is neglected since the integration of them becomes zero. Therefore, the structure functions are transformed into the following forms:
(25.1)
(25.2)
where C and are the normalizing constants that are computed using the Luke’s theorem [15].
For the next computations, it is essential to have the structure functions as the functions of parameters without, while the structure functions achieved in Equations (25.1) and (25.2) are not only the functions of but also the functions of x and y which should be first integrated. However, since the Mathematica software is unable to parametrically solve these integrals, and also regarding the fact that the structure functions are normal merely up to order, before the integration of x and y, they are expanded according to, and then x and y are integrated.
Generally, the product of is determined as a function of; and is obtained through extrapolating the data at the recoil of zero () (and are a combination of and functionsrespectively). For this purpose, the values of for various s are empirically achieved, a linear or square curve is fitted to them, and thereby obtaining at; consequently,. Therefore, the structure functions’ slope and curve achieved theoretically are of importance. Accordingly, the main reason of inaccuracy in the theoretical computations for determining pertains to the value of and the shape of curve used for fitting the laboratory data. In other words, various values of slope and curve result in different values of. Hence, the structure functions are also written as an expansion around. Since the ω’s range of values available is narrow (), the expansion of structure functions around is logical [16]:
(26)
The slope and the curve are considered as parameter.
4. Conclusion
The results obtained using the zero-order Bethe-Salpeter wave function are similar to the results the HQET predicts for extremely heavy quark at the maximum recoil. Table 1 provides the results achieved using the first-order B-S wave function, some of empirical values, and the results obtained through other theoretical methods. As can be observed, the results acquired by means of the first-order B-S wave function for heavy-to-heavy transitions are appropriately consistent with the empirical values and values obtained in other theoretical methods.
5. Acknowledgement
The work described in this paper was fully supported by grants from the Institute for Advanced Studies of Iran. The authors would like to express genuinely and sincerely thanks and appreciated and their gratitude to Institute for Advanced Studies of Iran.
NOTES