Decay Constants of S Wave Heavy Quarkonia

We derive expressions for the decay constants of S wave pseudoscalar and vector quarkonium (QQ̄\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$Q\bar {Q}$\end{document}) states through approximating the quark model definition by including terms upto k2/c2. Using the derived formula, the decays constants for ground and excited states are calculated. Our predictions are found to be in good agreement with the available experimental and lattice results.

potential models have been largely successful in explaining hadronic properties. Importance of potential models stems from the fact that masses (upto rather high excitations) and decay properties can be investigated in an unified way.
One of the important quantities that can be obtained through potential models are the vector and pseudoscalar decay constants of heavy mesons. Decay constants are important because it can be used to estimate the hadronic matrix elements and provide information about the short distance structure of hadrons and non-perturbative QCD dynamics [19,27]. Precise knowledge of their values play very important roles, for example in determination of CKM matrix elements, study of weak decays, meson mixing, etc. [44,50]. They are input parameters which constrains various decays and processes and their precise knowledge is important in various experimental measurements [51,52]. Vector (f V ) and pseudoscalar (f P ) decay constants are defined through the matrix elements [27]: The quark model definition of decay constants are through [25,27]: The above expressions are mostly used in potential models which use the momentum space representation. In the non-relativistic limit (k → 0), one obtains the well-known Van Royen-Weisskopf formula which relates the decay constants to the value of the meson wave-function at the origin [53]: In (3), ψ P /V (0) is the wave function of the pseudoscalar/vector meson calculated at the origin (r = 0) and M P /V is the mass of the pseudoscalar/vector meson. This formula is widely used in potential model calculations [22-25, 27, 29, 32-35]. One reason that it is widely used is that the coordinate space wave function is sufficient to calculate the decay constants. In potential models, in which the vector and pseudoscalar states differ only by the spindependent corrections, one obtains f P ≈ f V [22-25, 32, 35, 54]. The relativistic effects which enters as the v 2 /c 2 corrections in hyperfine spltting cannot be neglected in quarkonia. They produce significant difference in the wavefunction at origin for pseudoscalar and vector mesons [55], and hence f P becomes very different from f V . Ebert et al. [25] have shown that inclusion of relativistic effects produces considerable (∼70 MeV) difference between vector and pseudoscalar decay constants. Also the present experimental estimates (f J/ψ = 416 ± 6 MeV [3] and f η c = 335 ± 75 MeV [56]), lattice [19,52] and other theoretical studies [19,44,45,48,50,57,58] predict that f V is greater than f P . Hence the non-relativistic Van Royen-Weisskopf formula fails to provide exact results in accordance to experimental and lattice results. We suggest that the higher order terms (O(k 2 /c 2 )) in (1) and (2) are important. In this paper, we derive expressions for the decay constants of S-wave quarkonia by approximating the quark model definition by including terms upto k 2 /c 2 .
This paper is organised as follows: in Section 2 we derive expressions for pseudoscalar and vector decay constants for S-wave quarkonia including the k 2 /c 2 corrections. In Section 3 we discuss the theoretical model used in the present analysis. Results and discussions of our present analysis are given in Section 4.

Decay Constants
For quarkonia m Q = mQ = m and hence E k = Ek = E. Also using the relation, Now, Substituting (5) in (4), we get Now using the three dimensional Fourier transform of φ(k) the second integral in (6) can be written as: For S-wave states, there is no angular dependence: Therefore, (7) becomes: Similarly the integral Substituting (9) and (10) in (6), we get the final expression for vector decay constant of S-wave quarkonia as: Similarly setting Substituting (5) in (12), we get: Fourier tranforming φ(k) and performing the integrals as done earlier, we obtain the final expression for pseudoscalar decay constant of S-wave quarkonia as: Taking into account the QCD corrections [59][60][61], the decay constants of vector and pseudoscalar states equal tof where δ V = 8/3 and δ P = 2.

Theoretical Model
In order to obtain the bound state wavefunctions of heavy mesons, we have used a screened potential model suggested by Li and Chao [30,31]. The quark-antiquark potential in this model is given by The hyperfine interaction term which provide the spin singlet-triplet splitting is where δ σ (r) = (σ/ √ π) 3 e −σ 2 r 2 . The hyperfine interaction is treated perturbatively. The parameters used in (15) and (16) are same as in References [30,31] and are listed in Table 1 for completeness.
Using the potential (15), we solve the corresponding non-relativistic Schrodinger equation numerically [62] and obtain the bound state wavefunctions for heavy mesons. Substituting the wavefunctions in (14) we obtain respectively the vector and pseudoscalar decay constants.
To check the validity of our numerical computation method, we calculate the mass spectra and the radial wavefunction at the origin (|R(0)| 2 ) and compare the results with that of References [30,31]. The results are presented in Tables 2 and 3 along with the spinaveraged masses (M SA ) and hyperfine contributions. H V hf and H P hf are respectively the hyperfine contributions for spin-triplet (vector) and spin-singlet (pseudoscalar) states. The |R(0)| 2 of References [30,31] are extracted from the corresponding leptonic decay widths. The mass spectra and |R(0)| 2 obtained in the present analysis are in good agreement with the corresponding results of References [30,31], confirming the validity of our numerical computation method.

Results and Discussions
The computed values of pseudoscalar and vector decay constants are presented in Tables 4  and 5. f NR andf NR are respectively the decay constants calculated using the nonrelativistic Van Royen-Weisskopf formula (3) without and with the QCD correction factor. f andf are respectively the decay constants calculated from our present analysis ((11) and (13)) without and with the QCD correction factor (14).
Vector decay constants are experimentally extracted from the leptonic decay widths through the relation M V Substituting the PDG averages for Γ (V → e + e − ) in the above expression, one can estimate the vector decay constants. As far as the pseudoscalar decay constants are concerned, only the f n c is experimentally measured through the B → n c K decay [56]. f n b is not yet experimentally determined.   As seen from Tables 6-9, the predictions of decay constants from various models show a wide range of variation. For most of the states, the experimental results are yet to be available. Results from our analysis are found to be in good agreement with the available experimental results. For both vector and pseudoscalar states, the decay constants keep decreasing with the radial excitation number which is also in accordance with existing experimental/lattice results and other theoretical studies.
The ratios of pseudoscalar to vector decay constant (f P /f V )) for the 1S states are shown in Table 10. For both cc and bb mesons, we obtainf P /f V < 1, which is consistent with the available experimental data (f η c /f J/ψ = 0.81 ± 0.18 [3,56]) and other theoretical predictions. This result is in sharp contrast with predictions of some of the nonrelativistic potential models [22-25, 32, 35] which predict f P /f V ≈ 1. This indicate that for quarkonia the relativistic effects are important in calculating the decay constants.
In summary, we have derived expressions for pseudoscalar and vector decay constants for S-wave quarkonia including the k 2 /c 2 corrections. Our predictions are consistent with the experiment and predictions from other theoretical studies.