Journal of Nuclear and Particle Physics

p-ISSN: 2167-6895    e-ISSN: 2167-6909

2020;  10(1): 13-22

doi:10.5923/j.jnpp.20201001.03

 

Spectroscopy of the Quarkonium Systems for Heavy Quarks

Hesham Mansour, Ahmed Gamal

Physics Department, Faculty of Science, Cairo University, Giza, Egypt

Correspondence to: Hesham Mansour, Physics Department, Faculty of Science, Cairo University, Giza, Egypt.

Email:

Copyright © 2020 The Author(s). Published by Scientific & Academic Publishing.

This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/

Abstract

The study of the spectroscopy of bound states of quarkonium systems like and meson in the quark model framework with phenomenological potentials is motivated by quantum chromodynamics (QCD). It is found that analysis of the mass spectra of these systems is effectively given by the nonrelativistic Schrödinger's equation. There are several methods which are used to solve Schrödinger's equation with a general polynomial potential one of them is the Nikiforov-Uvarov (NU) method. It’s one of the effective methods which gives the energy eigenvalues and eigenstates for our potential. The results obtained are in good agreement with the experimental data and are better than previous theoretical studies.

Keywords: Quarkonium spectroscopy, Meson Mass Spectrum

Cite this paper: Hesham Mansour, Ahmed Gamal, Spectroscopy of the Quarkonium Systems for Heavy Quarks, Journal of Nuclear and Particle Physics, Vol. 10 No. 1, 2020, pp. 13-22. doi: 10.5923/j.jnpp.20201001.03.

1. Introduction

In quark and anti-quark system, the quantitively description is given by (quantum chromodynamics (QCD) spectroscopy and the standard model theory) [1-4] is important to specify the mechanism of that system and its nature of being bound systems.
The Schrödinger's equation describes quarkonium systems [5-9] with a heavy quark and anti-quark interaction (two-body problem). We solve the Schrödinger equation in a spherical-symmetric coordinate and using radial potentials which can be described by the asymptotic limits of QCD which has been qualitatively verified by Lattice QCD calculations [2-4].
The main purpose of this paper is that an interaction potential in the quark-antiquark bound system is taken as a general polynomial to get the general eigenvalue and eigenfunction solution then choosing a specific potential according to the description of the physical mechanism of the system.
The Nikiforov-Uvarov (NU) method [10-15], gives asymptotic expressions for the eigenfunctions and eigenvalues of the Schrödinger's equation. Hence one can calculate the energy eigenvalues and eigenstates for the spectrum of the quarkonium systems [17-26].

2. The Schrödinger Equation with Polynomial Potential

The Schrödinger equation reads
(1)
We use the generalized potential
(2)
(3)
By substituting in Schrödinger equation (1), we get
(4)
(5)
Where
(6)
and
(7)
Let and substituting in equation (5), we get
(8)
Because of singularity, we expand the polynomial function by Taylor’s series around x=0
Let
(9)
(10)
Neglecting higher order terms, we get
(11)
By substituing , we get
(12)
By substituting in equation (9), we obtain
(13)
(14)
By substituting in equation (8), we get
(15)
We rearrange equation (15), and divide by where to obtain,
(16)
(17)
Equation (16) becomes
(18)
We use the Nikiforov-Uvarov (NU) method [12,14] as mentioned before,
are symbols used in the Nikiforov-Uvarov (NU) method.
(19)
(20)
(21)
We chose the value of which make the square root in equation (21) as a quadratic term
(22)
(23)
Substituting in equation (21), we get
(24)
(25)
We take the negative value
(26)
(27)
(28)
(29)
(30)
(31)
By equating equations (29) and (31), we get the eigen value equation
(32)
(33)
(34)
By substituting equation (17) in equation (34), we get
(35)
By substituting equation (7) in equation (35), we get
(36)
Putting and substituting in equation (36), we obtain
(37)
So, the total energy eigen value is
(38)
To find the eigenfunctions for the general potential form
(39)
First, we calculate
(40)
(41)
(42)
Second, we calculate
(43)
(44)
(45)
(46)
(47)
By substituting in equation (39), we obtain
(48)
Where is a normalization constant.
So, the radial wavefunction of Schrödinger equation (1) for a polynomial potential is
(49)
To find the normalization constant
(50)
The angular part of the spherical symmetric potential is
(51)
So, the total wavefunction in spherical symmetric potential is
(52)
In our interquark potential, and using the natural units
(53)
Putting
(54)
The energy eigen values equation according to such potentials become
(55)

3. Results and Discussion

In this section, we will calculate the spectra for the bound states of heavy quarks such as charmonium, bottomonium and meson. The mass spectra equation is
(56)
By substituting equation (55) in equation (56), we get
(57)
Equation (57) depends on the potential parameters which will be obtained from the experimental data.
In the case of charmonium the rest mass equation is
(58)
And we get the following masses in GeV
The charmonium system
     
     
In the following we draw the radial wave functions as a function of the radius,
Figure 1. Graphs represent the relation between the radial wave function and the radius in different n, L states according to experimental energy states in the charmonium system
In the case of bottomonium , the rest mass equation is
(59)
And we get the following masses in GeV
The bottomonium system
     
     
In the following we draw the radial wave functions as a function of the radius,
Figure 2. Graphs represent the relation between the radial wave function and the radius in different n, L according to experimental energy states in the bottomonium system
Figure 3. Graphs represent the relation between the radial wave function and the radius in different n, L according to experimental energy states in the bottomonium system
In the case of the meson the rest mass equation is
(60)
And we get the following masses in GeV
The meson
     
     
In the following we draw the radial wave functions as a function of the radius,
Figure 4. Graphs represent the relation between the radial wave function and the radius in different n, L according to experimental energy states in the meson system
In conclusion, form the tables, we found that our theoretical work is comparable with the experimental data and explains the behavior of the quarkonium systems. The difference between the experimental data and theoretical work may be because we neglect the spin terms, so, the spin can also be considered if one uses relativistic corrections and the appropriate relativistic Schrödinger's equation. From the figures which represent the quarkonium radial state wave functions, one can calculate physical parameters like the decay parameter.

References

[1]  F.Halzem, A.Martin, Quarks and Leptons an Introductory Course in Modern Particle Physics, Wiley publication, 1984.
[2]  David Griffiths, introduction to elementary particles, Wiley publication, 1987.
[3]  Donald H. Perkins, Introduction to High Energy Physics, Cambridge university press, 2000.
[4]  B. R. Martin, Graham G Shaw, 2nd Edition, Particle Physics, the Manchester Physics Series, 1997.
[5]  B. H. Bransden and C. J. Joachain, Quantum Mechanics, Prentice Hall, 2000.
[6]  L. D. Landau and E. M. LIFSHITZ, Quantum mechanics, Butterworth Heinemann, 1981.
[7]  David Griffiths, introduction to elementary particle, Wiley -VCH, 2008.
[8]  Nouredine Zettili, Quantum Mechanics Concepts and Applications, John Wiley, 2009.
[9]  Cüneyt Berkdemir, theoretical concept of quantum mechanics, In Tech, 2012.
[10]  F Yaşuk, C Berkdemir and A Berkdemir, Exact solutions of the Schrodinger equation with non-central potential by the Nikiforov–Uvarov method, 38, 6579 (2005).
[11]  Yan-Fu Cheng and Tong-Qing Dai, Exact solution of the Schrödinger equation for the modified Kratzer potential plus a ring-shaped potential by the Nikiforov–Uvarov method, Phys. Scr. 75, 274 (2007).
[12]  [Özlem Ye, silta, PT/non-PT symmetric and non-Hermitian Pöschl–Teller-like solvable potentials via Nikiforov–Uvarov method, Phys. Scr. 75, 41 (2007).
[13]  H. Hassanabadi1, S. Zarrinkamar2 and A.A. Rajabi, Exact Solutions of D-Dimensional Schrödinger Equation for an Energy-Dependent Potential by NU Method, Communications in Theoretical Physics journal, 55, 541 (2011).
[14]  Mohammad R Setare and S Haidari, Bound states of the Dirac equation with some physical potentials by the Nikiforov–Uvarov method, Phys. Scr. 81, 015201 (5pp) (2010).
[15]  Birkhauser-Boston; H. Karayer, D. Demirhan, and F. Büyükkılıç, math.phys.j. 56, 063504 (2015).
[16]  A. Kumar Ray and P. C. Vinodkumar, J. Phys. 66, 953 (2006),; Said Laachir and Aziz Laaribi, Inter. J. of. Math. Compu. Phys. Electr. and Compu. Eng. 7(1) (2013); Hakan Ciftci, Richard L. Hall, and Nasser Saad, Phys. Rev. A 72, 022101 (2005).
[17]  E. J. Eichten and C Quigg, Phys. Rev. D 49, 5845 (1994); Mary Alberg, L.Wilets, Phys Lett. A 286, 7 (2001).
[18]  D. Ebert, R. N. Faustov and V. O. Galkin, Phys. Rev. 67, 014027 (2003); B. Ita, P. Tchoua, E. Siryabe, G. E. Ntamack, Inter. J. of Theor. And Math. Phys. 4(5), 173 (2014); H. Goudarzi 1 and V. Vahidi, Adv. Studies Theor. Phys. 5(10), 469 (2011).
[19]  R. Kumar and F. Chand, Commun. Theor. Phys. 59, 528 (2013); M. Abu-Shady, Inter. J. of Appl. Math. and Theor. Phys. 2(2), 16 (2016).
[20]  A. Al-Jamel and H. Widyan, Appl. Phys. Rese. 4, 94 (2013); CHANG ChaoHsi and WANG GuoLi, SCI. chinese 53(11), 2025 (2010).
[21]  N. V. Masksimenko and S. M. Kuchin, Russ. Phy. J. 54, 57 (2011); E. Eichten, T. Kinoshita, K. Gottfried, et al., Phys. Rev. Lett. 34, 369 (1975).
[22]  Z. Ghalenovi, A. A. Rajabi, A. Tavakolinezhad, Mod. Phys. inter. J. 5(6), 1250057 (2012); Juan-Juan Niu, a Lei Guo, b, Hong-Hao Ma,c, Shao-Ming Wang, Eur. Phys. J. C 78, 657 (2018).
[23]  S. M. Kuchin and N. V. Maksimenko, Univ. J. Phys. Appl. 7, 295 (2013); R. Faccini, BABAR-CONF-07/035 SLAC-PUB-13080.
[24]  R. Kumar and F. Chand, Phys. Scr. 85, 055008 (2012); Yang Li, Pieter Maris, and James P.Vary, Phys. Rev. D 96, 016022 (2017).
[25]  J. Beringer et al. [Particle Data Group], Phys. Rev. D 86, 1 (2012); Byungsik Hong, Overview of quarkonium production in heavy-ion collisions at LHC, published by EDP Sciences (2015).
[26]  Yu-Qi Chen 1'3, Yu-Ping Kuang 1'2'3, General relations of heavy quark-antiquark potentials induced by reparameterization invariance, Z. Phys. C 67, 627 (1995).
[27]  Q ¯ Q (Q ∈ {b, c}) spectroscopy using Cornell potential, N. R. Soni,∗ B. R. Joshi,† R. P. Shah,‡ H. R. Chauhan,§ and J. N. Pandya, European Physical Journal C 78(7), (2017).
[28]  D. Ebert, R. N. Faustov, and V. O. Galkin, Spectra of heavy mesons in the Bethe-Salpeter approach,Eur. Phys. J. C 71, 1825 (2011).
[29]  W.-J. Deng, H. Liu, L.-C. Gui, and X.-H. Zhong, Properties of Low-Lying Charmonium States in a Phenomenological Approach, Phys. Rev. D 95, 034026, 1608.00287 (2017).
[30]  C. S. Fischer, S. Kubrak, and R. Williams, Spectra of heavy mesons in the Bethe-Salpeter approach,Eur. Phys. J. A 51, 1409.5076, 10 (2015).
[31]  B.-Q. Li and K.-T. Chao, Higher Charmonia and X,Y,Z states with Screened Potential, Phys. Rev. D 79, 094004 (2009).
[32]  B.-Q. Li and K.-T. Chao, Bottomonium Spectrum with Screened Potential, Commun. Theor. Phys. 52, 0909. 1369, 653 (2009).
[33]  N. Devlani, V. Kher, and A. K. Rai, Masses and electromagnetic transitions of the Bc mesons, Eur. Phys. J. A 50, 154 (2014).
[34]  S. Godfrey, Phys. Rev. D 70, Spectroscopy of Bc mesons in the relativized quark model, 054017, hep-ph/0406228, (2004).
[35]  A. P. Monteiro, M. Bhat, and K. B. Vijaya Kumar, cb¯ spectrum and decay properties with coupled channel effects, Phys. Rev. D 95, 1608.05782, 054016 (2017).