Abstract
We derive expressions for the decay constants of S wave pseudoscalar and vector quarkonium (\(Q\bar {Q}\)) 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.
Similar content being viewed by others
1 Introduction
Ever since the discovery of first charmonium state (J/Ψ) at SLAC and BNL, quarkonia spectroscopy has always been a subject of continuous research, both experimentally and theoretically [1, 2]. Thanks to the experimental facilities at various colliders, new states have been constantly discovered [2, 3]. Recently discovered conventional quarkonium states include ηb(1S), ηb(2S), hb(1P), hb(2P), χb1(3P) and χb2(3P) [3]. There are also lot of states that are discovered (named X, Y, Z) that do not fit into the conventional quarkonia spectra and are assumed to be strong candidates for non-conventional hadronic states [2,3,4,5,6]. Comprehensive reviews on the status of heavy quarkonium physics can be found in References [5,6,7,8,9,10,11,12,13].
Even though quantum chromodynamics (QCD) is the theory of strong interactions, there is no direct method to extract the properties of observed hadrons directly from QCD Lagrangian. Lattice QCD provides great hope in this direction and significant results are being reported [1, 2]. Various other approaches widely used in the literature to study hadronic systems include QCD sum rules [14,15,16,17,18,19], nonrelativistic and relativistic potential models [20,21,22,23,24,25,26,27,28,29,30, 32,33,34,35,36], effective field theory [37,38,39] and formalisms based on Bethe-Salpeter equation and Dyson-Schwinger equation [40,41,42,43,44,45,46,47,48,49]. Among theses QCD inspired 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 (fV) and pseudoscalar (fP) 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\rightarrow 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 MP/V is the mass of the pseudoscalar/vector meson. This formula is widely used in potential model calculations [22,23,24,25, 27, 29, 32,33,34,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 spin-dependent corrections, one obtains fP ≈ fV [22,23,24,25, 32, 35, 54]. The relativistic effects which enters as the v2/c2 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 fP becomes very different from fV. Ebert et al. [25] have shown that inclusion of relativistic effects produces considerable (\(\sim \)70 MeV) difference between vector and pseudoscalar decay constants. Also the present experimental estimates (fJ/ψ = 416 ± 6 MeV [3] and \(f_{\eta _{c}}=335\pm 75\) MeV [56]), lattice [19, 52] and other theoretical studies [19, 44, 45, 48, 50, 57, 58] predict that fV is greater than fP. 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 (\(\mathcal {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 k2/c2.
This paper is organised as follows: in Section 2 we derive expressions for pseudoscalar and vector decay constants for S-wave quarkonia including the k2/c2 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.
2 Decay Constants
For quarkonia \(m_{Q}=m_{\bar {Q}}=m\) and hence \(E_{k}=E_{\bar {k}}=E\). Also using the relation, k2 = (E + m)(E − m), (1) becomes
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 \(m_{Q}=m_{\bar {Q}}=m\), \(E_{k}=E_{\bar {k}}=E\) and k2 = (E + m)(E − m), (2) becomes
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 to
where δV = 8/3 and δP = 2.
3 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 \(\delta _{\sigma } (r)=(\sigma /\sqrt {\pi })^{3} e^{-\sigma ^{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 spin-averaged masses (MSA) and hyperfine contributions. \(\langle H_{hf}^{V} \rangle \) and \(\langle H_{hf}^{P} \rangle \) 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.
4 Results and Discussions
The computed values of pseudoscalar and vector decay constants are presented in Tables 4 and 5. fNR and \(\bar {f}^{NR}\) are respectively the decay constants calculated using the non-relativistic Van Royen-Weisskopf formula (3) without and with the QCD correction factor. f and \(\bar {f}\) are respectively the decay constants calculated from our present analysis ((11) and (13)) without and with the QCD correction factor (14).
In Tables 6, 7, 8 and 9, pseudoscalar (\(\bar {f}_{P}\)) and vector (\(\bar {f}_{V}\)) decay constants for charmonium and bottomonium states obtained from our analysis are compared with the available experimental results [3, 56] and predictions from lattice [19, 52, 63,64,65,66,67,68], nonrelativistic potential models [27, 69], QCD sum rules (QCDSR) [19, 57, 58], light front quark model (LFQM) [50] and formalisms based on Bethe-Salpeter and Dyson-Schwinger equations (BS/DSE) [44,45,46,47,48,49].
Vector decay constants are experimentally extracted from the leptonic decay widths through the relation
Substituting the PDG averages for \({\varGamma } (V\rightarrow 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\rightarrow 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 (\(\bar {f}_{P}/\bar {f}_{V}\))) for the 1S states are shown in Table 10. For both \(c\overline c\) and \( b\overline b\) mesons, we obtain \(\bar {f}_{P}/\bar {f}_{V} < 1\), which is consistent with the available experimental data (\(f_{\eta _{c}}/f_{J/\psi }=0.81\pm 0.18\) [3, 56]) and other theoretical predictions. This result is in sharp contrast with predictions of some of the nonrelativistic potential models [22,23,24,25, 32, 35] which predict fP/fV ≈ 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 k2/c2 corrections. Our predictions are consistent with the experiment and predictions from other theoretical studies.
References
Bodwin, G.T., et al.: A Snowmass White Paper. arXiv:1307.7425v3
Brambilla, N, et al., (Quarkonium Working Group): . Eur. Phys. J. C 71, 1534 (2011). arXiv:1010.5827v3
Patrignani, C, et al., (Particle Data Group): . Chin. Phys. C 40, 100001 (2016)
Godfrey, S., Olsen, S.L.: Ann. Rev. Nucl. Part. Sci. 58, 51 (2008). arXiv:0801.3867v1
Olsen, S.L.: Front. Phys. (Beijing) 10, 121 (2015). arXiv:1411.7738
Lebed, R.F., Mitchell, R.E., Swanson, E.S.: Prog. Part. Nucl. Phys. 93, 143 (2017). arXiv:1610.04528
Esposito, A, et al.: Int. J. Mod. Phys. A 30, 1530002 (2015). arXiv:1411.5997
Chen, H.X., Chen, W., Liu, X., Zhu, S.L.: Phys. Rept. 639, 1 (2016). arXiv:1601.02092
Ali, A., Lange, J.S., Stone, S.: Prog. Part. Nucl. Phys. 97, 123 (2017). arXiv:1706.00610
Guo, F.K, et al.: Rev. Mod. Phys. 90, 015004 (2018). arXiv:1705.00141
Olsen, S.L., Skwarnicki, T., Zieminska, D.: Rev. Mod. Phys. 90, 015003 (2018). arXiv:1708.04012
Brambilla, N, et al.: Eur. Phys. J. C 71, 1534 (2011). arXiv:1010.5827
Brambilla, N., et al.: Eur. Phys. J. C 74, 2981 (2014). arXiv:1404.3723
Shifman, M.A., Vainshtein, A.I., Zakharov, V.I.: Nucl. Phys. B 147, 385 (1979). Nucl. Phys. B 14, 448 (1979)
Reinders, L.J., Rubinstein, H.R., Yazaki, S.: Nucl. Phys. B 186, 109 (1981)
Novikov, V.A., Okun, L.B., Shifman, M.A., Vainshtein, A.I., Voloshin, M.B., Zakharov, V.I.: Phys. Rept. 41, 1 (1978)
Beilin, V.A., Radyushkin, A.V.: Nucl. Phys. B 260, 61 (1985)
Nielsen, M., Navarra, F.S., Lee, S.H.: Phys. Rept. 497, 41 (2010). arXiv:0911.1958v2
Becirevic, D., Duplancic, G., Klajn, B., Melic, B., Sanfilippo, F.: Nucl. Phys. B 883, 306 (2014). arXiv:1312.2858v3
Eichten, E., Gottfried, K., Kinoshita, T., Kogut, J., Lane, K.D., Yan, T.M.: Phys. Rev. Lett. 34, 369 (1975)
Eichten, E., Gottfried, K., Kinoshita, T., Lane, K.D., Yan, T. M.: Phys. Rev. D 17, 3090 (1978). Erratum Phys. Rev. D 21, 313 (1980)
Eichten, E.J., Quigg, C.: Phys. Rev. D 49, 5845 (1994). arXiv:hep-ph/9402210v1
Gershtein, S.S., Kiselev, V.V., Likhoded, A.K., Tkabladze, A.V.: Phys. Rev. D 51, 3613 (1995). arXiv:hep-ph/9406339v1
Fulcher, L.P.: Phys. Rev. D 60, 074006 (1999). arXiv:hep-ph/9806444v2
Ebert, D., Faustov, R.N., Galkin, V.O.: Phys. Rev. D 67, 014027 (2003). arXiv:hep-ph/0210381v2
Barnes, T., Godfrey, S., Swanson, E.S.: Phys. Rev. D 72, 054026 (2005). arXiv:hep-ph/0505002v3
Lakhina, O., Swanson, E.S.: Phys. Rev. D 74, 014012 (2006). arXiv:hep-ph/0603164v2
Radford, S.F., Repko, W.W.: Phys. Rev. D 074031, 75 (2007). arXiv:hep-ph/0701117v3
Ikhdair, S.M., Sever, R.: Int. J. Mod. Phys. A 21, 3989 (2006). arXiv:hep-ph/0508144v1
Li, B.-Q., Chao, K.-T.: Phys. Rev. D 79, 094004 (2009). arXiv:0903.5506v2
Li, B.-Q., Chao, K.-T.: Commun. Theor. Phys. 52, 653 (2009). arXiv:0909.1369v1
Rai, A.K., Patel, B., Vinodkumar, P.C.: Phys. Rev. C 78, 055202 (2008). arXiv:0810.1832v1
Patel, S., Vinodkumar, P.C., Bhatnagar, S.: nopunct. Chin. Phys. C 40, 053102 (2016). arXiv:1504.01103v3
Monteiro, A.P., Bhat, M., Vijaya Kumar, K.B.: Phys. Rev. D 95, 054016 (2017). arXiv:1608.05782v3
Kher, V., Rai, A.K.: Chin. Phys. C 42, 083101 (2018). arXiv:1805.02534v1
D’Souza, P.P., Monteiro, A.P., Vijaya Kumar, K.B.: Commun. Theor. Phys. 71, 192 (2019). arXiv:1703.10413
Brambilla, N., Pineda, A., Soto, J., Vairo, A: arXiv:hep-ph/0410047v2. Rev. Mod. Phys. 77, 1423 (2005)
Bodwin, G.T., Braaten, E., Lepage, G.P.: Phys Rev. D 51, 1125 (1995). Erratum:, Phys. Rev. D 55, 5853 (1997)
Guo, F.-K., Hanhart, C., Li, G., Meißner, U.G., Zhao, Q.: Phys. Rev. D 83, 034013 (2011). arXiv:1008.3632v2
Smith, C.H.L.: Ann. Phys. 53, 521 (1969)
Alkofer, R., von Smekel, L.: Phys. Rep. 353, 281 (2001). arXiv:hep-ph/0007355v2
Maris, P., Roberts, D.: Int. J. Mod. Phys. E 12, 297 (2003). arXiv:nucl-th/0301049v1
Bhatnagar, S., Alemu, L.: Phys. Rev. D 97, 034021 (2018). arXiv:1610.03234v5
Wang, G.-L.: Phys. Lett. B 633, 492 (2006). arXiv:math-ph/0512009v1
Cvetic, G., Kim, C.S., Wang, G.-Li., Namgung, W.: Phys. Lett. B 596, 84 (2004)
Wang, Z.-G., Yang, W.-M., Wan, S.-L.: Phys. Lett. B 615, 79 (2005). arXiv:hep-ph/0411142v2
Negash, H., Bhatnagar, S.: Intl. J. Mod. Phys. E 25, 1650059 (2016). arXiv:1508.06131v3
Blank, M., Krassnigg, A.: Phys. Rev. D 84, 096014 (2011). arXiv:1109.6509v1
Bhagwat, M.S., Maris, P.: Phys. Rev. C 77, 025203 (2008). arXiv:nucl-th/0612069v1
Choi, H.-M.: Phys. Rev. D 75, 073016 (2007). arXiv:hep-ph/0701263v2
Mojica, F.F., Vera, C.E., Rojas, E., El-Bennich, B.: Phys. Rev. D 96, 014012 (2017). arXiv:1704.08593v1
Bailas, G., Blossier, B., Morenas, V.: Eur. Phys. J. C 78, 1018 (2018). arXiv:1803.09673v1
Van Royen, R., Weisskopf, V.F.: Nuovo Cim. 50, 617 (1967). ibid. 51, 583 (1967)(E)
Hwang, D.S., Kim, G.-H.: Z. Phys. C 76, 107 (1997). arXiv:hep-ph/9703364v1
Ahmady, M.R., Mendel, R.R.: Phys. Rev. D 51, 141 (1995). arXiv:hep-ph/9401315v2
Edward, K.W., et al., CLEO Collaboration: . Phys. Rev. Lett. 86, 30 (2001). arXiv:1010.3110v2
Veli Veliev, E., et al.: J. Phys. G Nucl. Part. Phys. 39, 015002 (2012). arXiv:1010.3110v2
Veli Veliev, E., et al.: Eur. Phys. J. A 47, 110 (2011). arXiv:1103.4330v1
Braaten, E., Fleming, S.: Phys. Rev. D 52, 181 (1995). arXiv:hep-ph/9501296v2
Kiselev, V.V.: Int. J. Mod. Phys. A 11, 3689 (1996)
Gershtein, S.S., et al.: arXiv:hep-ph/9803433v1 (1998)
Lucha, W., Schoberl, F.F.: Int. J. Mod. Phys. C 10, 607 (1999). arXiv:hep-ph/9811453v2
Dudek, J.J., Edwards, R.G., Richards, D.G.: Phys. Rev. D 73, 074507 (2006). arXiv:hep-ph/0601137v2
Davies, C.T.H., et al., HPQCD collaboration: . Phys. Rev. D 82, 114504 (2010). arXiv::1008.4018v2
Donald, G.C., et al., HPQCD collaboration: . Phys. Rev. D 86, 094501 (2012). arXiv:1208.2855v2
McNeile, C., et al., HPQCD collaboration: . Phys. Rev. D 86, 074503 (2012). arXiv:1207.0994v1
Colquhoun, B., et al., HPQCD collaboration: . Phys. Rev. D 91, 074514 (2015). arXiv:1408.5768v1
Chiu, T.W., et al., TWQCD Collaboration: . Phys. Lett. B 651, 171 (2007). arXiv:0705.2797v1
Soni, N.R., Joshi, B.R., Shah, R.P., Chauhan, H.R., Pandya, J.N.: Eur. Phys. J. C 78, 592 (2018). arXiv:1707.07144v2
Acknowledgments
One of the authors (Bhaghyesh) would like to thank Manipal Academy of Higher Education (MAHE), Manipal for covering the article processing charges for open access publication.
Funding
Open access funding provided by Manipal Academy of Higher Education, Manipal.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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/.
About this article
Cite this article
Azhothkaran, B., K., N.V. Decay Constants of S Wave Heavy Quarkonia. Int J Theor Phys 59, 2016–2028 (2020). https://doi.org/10.1007/s10773-020-04474-5
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10773-020-04474-5