Table of Contents Author Guidelines Submit a Manuscript
Advances in High Energy Physics
Volume 2018, Article ID 5961031, 12 pages
https://doi.org/10.1155/2018/5961031
Research Article

S-Wave Heavy Quarkonium Spectra: Mass, Decays, and Transitions

Physics Department, Faculty of Arts and Sciences, Ondokuz Mayis University, 55139 Samsun, Turkey

Correspondence should be addressed to Halil Mutuk; moc.liamg@kutumlilah

Received 20 August 2018; Revised 18 October 2018; Accepted 30 October 2018; Published 25 November 2018

Guest Editor: Xian-Wei Kang

Copyright © 2018 Halil Mutuk. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP3.

Abstract

In this paper we revisited phenomenological potentials. We studied S-wave heavy quarkonium spectra by two potential models. The first one is power potential and the second one is logarithmic potential. We calculated spin averaged masses, hyperfine splittings, Regge trajectories of pseudoscalar and vector mesons, decay constants, leptonic decay widths, two-photon and two-gluon decay widths, and some allowed M1 transitions. We studied ground and 4 radially excited S-wave charmonium and bottomonium states via solving nonrelativistic Schrödinger equation. Although the potentials which were studied in this paper are not directly QCD motivated potential, obtained results agree well with experimental data and other theoretical studies.

1. Introduction

Heavy quarkonium is the bound state of and and one of the most important playgrounds for our understanding of the strong interactions of quarks and gluons. Quantum chromodynamics (QCD) is thought to be the true theory of these strong interactions. QCD is a nonabelian local gauge field theory with the symmetry group . In principle, one should be able to calculate hadronic properties such as mass spectrum and transitions by using QCD principles. But QCD does not readily supply us these hadronic properties. This challenge can be attributed to the several features that are not present in other local gauge field theories.

Foremost, being a nonabelian gauge theory, gluons which are gauge bosons, have color charge and interact among themselves. Unlike from quantum electrodynamics (QED), where a photon does not interact with other photon, in QCD one must consider interactions among gluons. This nonabelian nature of the theory makes some calculations complicated, for example, loops in propagators.

There are three other important features of QCD: asymptotic freedom, confinement, and dynamical breaking of chiral symmetry. Asymptotic freedom says that strong interaction coupling constant, , is a function of momentum transfer. When the momentum transfer in a quark-quark collision increases (at short distances), the coupling constant becomes weaker whereas it becomes larger when momentum transfer decreases (at large distances). The idea behind confinement is that, there are no free quarks outside of a hadron; i.e., color charged particles (quarks and gluons) cannot be isolated out of hadrons. Flux tube model gives a reasonable explanation of confinement. When the distance between quark-antiquark (or quarks) pair increases, the gluon field between a pair of color charges forms a flux tube (or string) between them resulting a potential energy which depends linearly on the distance, where is the string constant. As distance increases between quarks, the potential energy can create new quark-antiquark pairs in colorless forms instead of a free quark. Up to now, nobody has been able to prove that confinement from QCD. Lattice QCD calculations simulate this confinement well and give a value for the string tension [1]. The last feature of QCD is the dynamical breaking of chiral symmetry. The QCD Lagrangian with quark flavor has an exact chiral symmetry but breaks down to symmetry because of the nonvanishing expectation value of QCD vacuum [2, 3]. The Goldstone bosons corresponding to this symmetry breaking are the pseudoscalar mesons.

The present aspects of the QCD caused other approaches to deal with these challenges. QCD sum rules, Lattice QCD, and potential models (quark models) are examples of these approaches. These approaches are nonperturbative since the strong interaction coupling constant, which should be the perturbation parameter of QCD is of the order one in low energies, hence the truncation of the perturbative expansion cannot be carried out. Since perturbation theory is not applicable, a nonperturbative approach has to be used to study systems that involve strong interactions. QCD sum rules and lattice QCD are based on QCD itself whereas in potential models, one assumes an interquark potential and solves a Schrödinger-like equation. The advantage of potential model is that, excited states can be studied in the framework of potential models whereas in QCD sum rules and lattice QCD, only the ground state or in some exceptional cases excited states can be studied.

After the discovery of charmonium () states, potential models have played a key role in understanding of heavy quarkonium spectroscopy [4, 5]. These potentials were in type of Coulomb plus linear confining potential with spin dependent interactions. The discovery of bottomonium () states were well described by the potential model picture which was used in the charmonium case. Heavy quarkonium spectroscopy was studied since that era with fruitful works [618]. A general review about potential models can be found in [19, 20] and references therein.

In the potential models, many features such as mass spectra and decay properties of heavy quarkonium could be described by an interquark potential in two-body Schrödinger equation. Interquark potentials are obtained both from phenomenology and theory. In the phenomenological method, it is assumed that a potential exist with some parameters to be determined by fits to the data. In the theory side, one can use perturbative QCD to determine the potential form at short distances and use lattice QCD at long distances [19]. These potentials can be classified as QCD motivated potentials [2125] and phenomenological potentials [2631]. The most commonly used phenomenological potentials are power-law potentials, for example [26] and logarithmic potentials, for example [30]. The detailed properties of these type potentials are studied extensively in [29]. All the potentials which are mentioned here have almost similar behaviour in the range of which is characteristic region of charmonium and bottomonium systems [32, 33]. Outside the range, the behaviour of potentials differ. Up to now, no one was able to obtain a potential which is compatible at the whole range of distances by using QCD principles.

The potential model calculations have been quite successful in describing the hadron spectrum. Most of the phenomenological potentials must satisfy the following conditions:It means that static potential is a monotone nondecreasing and concave function of which is a general property of gauge theories [34].

The great success of quarkonium phenomenology was somehow cracked at 2003 after the observation of [35]. The properties of this exotic particle are not compatible with the conventional quark model, the reason why it is named exotic. For example in [36], the authors studied near threshold zero in the S-wave. There are other exotic states, , and the exotic particle zoo is growing. In this paper we will present some exotic states in the framework of quark model.

Energy spectra of heavy quarkonium are a rich source of the information on the nature of interquark forces and decay mechanisms. The prediction of mass spectrum in accordance with the experimental data does not verify the validity of a model for explaining hadronic interactions. Different potentials can produce reliable spectra with the experimental data. Thus other physical properties such as decay constants, leptonic decay widths, radiative decay widths, etc. need to be calculated.

A specific form of the QCD potential in the whole range of distances is not known. Therefore one needs to use potential models. In this work we revisited a power-law potential [26] and a logarithmic potential [30] to study S-wave heavy quarkonium. These potentials satisfy Eqn. (1), i.e. having nonsingular behaviour for . For our purposes, it must be mentioned that power-law and logarithmic potentials have nice scaling properties when used with a nonrelativistic Schrödinger equation [19]. We generated S-wave charmonium and bottomonium mass spectrum with the decays and M1 transitions. At Section 2 we give out theoretical model. In Sections 3 and 4, we generate S-wave heavy quarkonium spectrum, decays and transitions. In Section 5 we discuss our results and in Section 6 we conclude our results.

2. Formulation of the Model

When quark model was proposed, many authors treated baryons in detail with the harmonic oscillator quark model by using harmonic oscillator wave functions [3739]. Mesons comparing to baryons are simpler objects since they are composites of two quarks. The reason for using harmonic oscillator wave function is that they form a complete set for a confining potential [40].

In order to obtain mass spectra, we solved Schrödinger equation by variational method. The variational method by using harmonic oscillator wave function gave successful results for heavy and light meson spectrum [15, 41, 42]. The procedure for this method is calculating expectation value of the Hamiltonian via the trial wave function:The mass of the meson is found by adding two times the mass of quark to the eigenenergyThe Hamiltonian we consider iswhere , is the relative momentum, is the reduced mass, and is the potential between quarks. The spectrum can be obtained via solving Schrödinger equationwith the harmonic oscillator wave function defined asHere is the radial wave function given aswith the associated Laguerre polynomials and the normalization constant is the well-known spherical harmonics.

Armed with these, the expectation value of the given Hamiltonian can be calculated. In the variational method, one chooses a trial wave function depending on one or more parameters and then finds the values of these parameters by minimizing the expectation value of the Hamiltonian. It is a good tool for finding ground state energies but as well as energies of excited states. The condition for obtaining excited states energies is that the trial wave function should be orthogonal to all the energy eigenfunctions corresponding to states having a lower energy than the energy level considered. In (7), is treated as a variational parameter and it is determined for each state by minimizing the expectation value of the Hamiltonian.

In the following sections we study power-law and logarithmic potentials in order to obtain full spectrum.

3. Mass Spectra of Power-Law and Logarithmic Potentials

Power-law potential is given by [26]They showed that upsilon and charmonium spectra can be fitted with that potential. The small power of refers to a situation in which the spacing of energy levels is independent of the quark masses. This situation is also valid for the purely logarithmic potential [30]

At first step we obtained spin averaged mass spectrum for and systems, respectively. The constituent quark masses are and for power-law potential and and for logarithmic potential. Table 1 shows the charmonium spectrum and Table 2 shows the bottomonium spectrum.

Table 1: Spin-averaged mass spectrum of charmonium (in MeV).
Table 2: Spin-averaged mass spectrum of Bottomonium (in MeV).

Since the interquark potential does not contain the spin dependent part, (2) gives the spin averaged mass for the corresponding states. The calculated masses agree well with the available experimental data and with the values obtained from other theoretical studies. A general potential usually includes spin-spin interaction, spin-orbit interaction, and tensor force terms. To obtain whole picture, it is necessary to consider spin dependent terms within the potential. For , there are spin-orbit and tensor force terms which contribute to the fine structure. For equal mass , the spin-orbit interaction is given by and is responsible for the wave splittings. Again for equal mass , the tensor potential is given byFor , there is spin-spin term which we will consider in the present work. In the model of the spin averaged mass spectra discussion, all the spin dependent effects are ignored and hence it fails to take into account the splittings due to spin. For example, such splitting exist between the (1S) and mesons by . These mesons occupy the level. The in the (1S) have , while in the , . As a result of this, the mass difference should be related to spin dependent interaction.

3.1. Spin-Spin Interaction

Mass splitting is closely connected with the Lorentz-structure of the quark potential [45]. The origin of the spin-spin interaction term lies in the one-gluon exchange term which is related to . Spin is proportional of the magnetic moment of a particle. Magnetic moments generate short range fields . In the case of heavy quarkonium systems which are nonrelativistic, wave functions of two particles overlap in a significant amount. This means that particles are very close to each other. So spin-spin interactions play a significant role in the dynamics. The spin-spin interaction term of two particles can be written asThis term can explain wave splittings and has no contribution to states. Putting this term into Schrödinger equation we getImplementing Dirac-delta function property we obtainThe matrix element of spin products can be obtained via so thatTherefore we obtain hyperfine splittings energy asHere is the wave function at the origin and can be obtained by the following relation:Expectation value is obtained by the wave function given in (6). S-wave charmonium and bottomonium masses can be seen in Tables 3 and 4. In this calculation, is taken to be 0.37 for charmonium and 0.26 for bottomonium [15].

Table 3: Charmonium mass spectrum (in MeV). In [18] LP denotes linear potential and SP denotes screened potential.
Table 4: Bottomonium mass spectrum (in MeV).

The mass differences are shown in Tables 5 and 6 for charmonium and bottomonium, respectively.

Table 5: Mass differences of S-wave charmonium states (in MeV).
Table 6: Mass differences of S-wave bottomonium states (in MeV).

As can be seen from Tables 3, 4, 5, and 6 our results are compatible with both experimental and theoretical results.

The Regge trajectories for pseudoscalar and vector mesons are shown in Figures 1 and 2 for charmonium and in Figures 3 and 4 for bottomonium.

Figure 1: Regge trajectories of pseudoscalar charmonium in plane. The polynomial fit is for power potential and for logarithmic potential.
Figure 2: Regge trajectories of vector charmonium in plane. The polynomial fit is for power potential and for logarithmic potential.
Figure 3: Regge trajectories of pseudoscalar bottomonium in plane. The polynomial fit is for power potential and for logarithmic potential.
Figure 4: Regge trajectories of vector bottomonium in plane. The polynomial fit is for power potential and for logarithmic potential.

As can be seen from figures, Regge trajectories show nonlinear behaviour.

4. Dynamical Properties

4.1. Decay Constants

Leptonic decay constants give information about short distance structure of hadrons. In the experiments this regime is testable since the momentum transfer is very large. The pseudoscalar () and the vector () decay constants are defined, respectively, through the matrix elements [12]andIn the first relation, is meson momentum and is pseudoscalar meson state. In the second relation, is mass, is the polarization vector, and is the state vector of meson.

The matrix elements can be calculated by quark model wave function in the momentum space. The result isfor pseudoscalar meson andfor the vector meson [12].

In the nonrelativistic limit, these two equations take a simple form which is known to be Van Royen and Weisskopf relation [46] for the meson decay constants

The first-order correction which is also known as QCD correction factor is given bywhere is given by [47]Here for pseudoscalar mesons and for vector mesons. Decay constants are given in Tables 7 and 8 for pseudoscalar and vector mesons, respectively.

Table 7: Pseudoscalar decay constants (in MeV).
Table 8: Vector decay constants (in MeV).
4.2. Leptonic Decay Widths

Leptonic decay of a vector meson with quantum numbers can be pictured by the following annihilation via a virtual photonThis state is neutral and in principle can decay into a different lepton pair rather than electron-positron pair. The above amplitude can be calculated by the Van Royen and Weisskopf relation [46]where is the fine structure constant, is the quark charge, is the mass of state, and is the wave function at the origin of initial state. The term in the parenthesis is the first-order QCD correction factor while represents higher corrections. The obtained values for leptonic decay widths can be found in Tables 9 and 10 for charmonium and bottomonium, respectively.

Table 9: Charmonium leptonic decay widths (in keV). The widths calculated with and without QCD corrections are denoted by and .
Table 10: Bottomonium leptonic decay widths (in keV). The widths calculated with and without QCD corrections are denoted by and .
4.3. Two-Photon Decay Width

states with quantum number of charmonium and bottomonium can decay into two photons. In the nonrelativistic limit, the decay width for state can be written as [48]The term in the parenthesis is the first-order QCD radiative correction. The results are listed in Table 11.

Table 11: Two-photon decay widths (in keV). The widths calculated with and without QCD corrections are denoted by and .
4.4. Two-Gluon Decay Width

Two-gluon decay width is given by [48]

The terms in the parenthesis refer to QCD corrections. The obtained results are given in Table 12.

Table 12: Two-gluon decay widths (in MeV). The widths calculated with and without QCD corrections are denoted by and .
4.5. M1 Transitions

M1 (magnetic dipole transition) decay widths can give more information about spin-singlet states. Moreover M1 transition rates show the validity of theory against experiment [11]. Magnetic transitions conserve both parity and orbital angular momentum of the initial and final states but in the M1 transitions the spin of the state changes. M1 width between two S-wave states is given by [51]where is the photon energy and is the zeroth-order spherical Bessel function. In the case of small , spherical Bessel function tends to 1, . Thus transitions between the same principal quantum numbers, , are favored and usually known to be allowed. In the other case, when the overlap integral between initial and final state is 0 and generally designated as forbidden transitions. The obtained transition rates for the allowed ones of S-wave charmonium and bottomonium states are given in Tables 13 and 14, respectively.

Table 13: Radiative M1 decay widths of charmonium. In [18] LP stands for linear potential and SP stands for screened potential.
Table 14: Radiative M1 decay widths of bottomonium.

5. Results and Discussion

In the present paper we studied S-wave heavy quarkonium spectra with two phenomenological potentials. We have computed spin averaged masses, hyperfine splittings, Regge trajectories for pseudoscalar and vector mesons, decay constants, leptonic decay widths, two-photon and gluon decay widths, and allowed M1 partial widths of S-wave heavy quarkonium states.

In general, most of the quark model potentials tend to be similar, having a Coulomb term and a linear term. For example, in [11] they used standard color Coulomb plus linear scalar form, and also included a Gaussian smeared contact hyperfine interaction in the zeroth-order potential. In [13], the authors used a nonrelativistic potential model with screening effect. In [18] nonrelativistic linear potential and screened potential, in [14, 16, 44] a modified of nonrelativistic potential models and in [15] Hulthen potential are used. Potential models give reliable results with the appropriate parameters in the model. Therefore, the shape of the potential at the limits and have similar behaviours.

Spin averaged mass spectra give idea about the formulation of model since the results are close to experimental values due to contributions from spin dependent interactions are small compared to contribution from potential part. If one ignores all spin dependent interactions, obtained results under this assumption are thought to be averages over related spin states for principal quantum number. Including hyperfine interaction, we obtained the mass spectra for pseudoscalar and vector mesons. The obtained spectra for both charmonium and bottomonium are in good agreement with the experimentally observed spectra and other theoretical studies.

Both power and logarithmic potentials produced approximately same mass differences and are in agreement with experiment for the lowest state in charmonium sector. But for the highest states, the shift is not compatible with the references. The reason for this should be the behaviour of linear part of the potential. In the case of bottomonium sector, mass differences of both power and logarithmic potentials are in accord with the given studies except the lowest state.

The fundamental point in the Regge trajectories is that they can predict masses of unobserved states. For the hadrons constituting of light quarks, the Regge trajectories are approximately linear but for the heavy quarkonium case Regge trajectories can be nonlinear. In the present work, we found that all Regge trajectories show nonlinear properties.

The decay constants are calculated for pseudoscalar and vector mesons by equating their field theoretical definition with the analogous quark model potential definition. This is valid in the nonrelativistic and weak binding limits where quark model state vectors form good representations of the Lorentz group [52, 53]. For pseudoscalar mesons, the corrected value of power and logarithmic potentials are a few MeV above than the available experimental data. For the rest of the pseudoscalar mesons, obtained results are compatible with other studies. In the case of vector mesons, logarithmic potential gave higher values than power potential. In the meson, when the radially states are excited, both two potential gave similar results within the error of experimental value. Computations of the vector decay constant beyond the weak binding limit can be important in the quark potential model frame and need more elaboration [12].

Obtained leptonic decay widths are comparable with the experimental values and other theoretical studies. The QCD corrected factors are more close to experimental values for power and logarithmic potential and this can be referred as the importance of the QCD correction factor in calculating the decay constants and other short range phenomena using potential models.

levels of charmonium and bottomonium states can decay into two photons or gluons. Especially two- photon decays of these levels are important for understanding the accuracy of theoretical models. Obtained results are smaller than the nonrelativistic widths including the one-loop QCD correction factor. For example, results of power and logarithmic potentials in (1S) are not in accord with experimental data. The reason of these differences can be due to the static potential between quarks that we used in the solution of two-body Schrödinger equation. For higher states, power and logarithmic potentials results are comparable with other studies. Two-photon decays are complicated processes such as pseudoscalar meson decay to two photons is governed by an intermediate vector meson followed by a meson dominance transition to a photon [12]. These schematic diagrams must be added to calculations to obtain a whole picture. For two-gluon decay widths, two phenomenological potentials gave comparable results with the available experimental data. Notice that in some cases QCD corrected factor is in accord with the experimental data whereas in some cases it is not. The reason for this can be that, there are significant radiative corrections from three-gluon decays so computing only two-gluon decay width could not explain the mechanism in all details.

Finally M1 transitions are calculated. The M1 radiative decay rates are very sensitive to relativistic effects. Even for allowed transitions relativistic and nonrelativistic results differ significantly. An important example is the decay of . The nonrelativistic predictions for its rate are more than two times larger than the experimental data [10]. In the charmonium sector, the available experimental data for is comparable with the power potential result, while logarithmic potential result is 1 eV higher. In the bottomonium sector, there is no experimental data available on M1 transitions. Since photon energies and transition rates are very small, the detection of these transitions is an objection. And this can be a reason why no spin-singlet S-wave levels have been observed yet [10]. The obtained values for M1 transitions are comparable with the references.

Some states in the charmonium and bottomonium sector show properties different from the conventional quarkonium state. Some examples are , , and . For , there is not much available experimental data and more is needed. Wang et al. studied two-body open charm OZI-allowed strong decays of and considered as and , respectively, by the improved Bethe-Salpeter method combined with the [54]. They calculated strong decay width of as and as where the experimental values are for and for [43]. They concluded that is a good candidate of and is a not good candidate of due to larger decay width of comparing to experimental data. We give our results comparing to these exotic states in Table 15.

Table 15: Exotic states. Experimental data are taken from [43] unless stated. The units for mass and strong decays are in MeV and two-photon decay is in keV.

Looking at Table 15, we can deduce that, according to our model and results, we can assign as , as , and as . To be more accurate, more data is needed to corroborate whether these states are conventional quarkonium or not.

6. Conclusions

Quark potential models have been very successful to study on various properties of mesons. The short distance behaviour of interquark potential appears to be similar where QCD perturbation theory can be applied where at large distance the potential is linear in where nonperturbative methods are need to be used. The improvements on the potentials can be made and spin-spin, spin-orbit type interactions can be added to model to arrive high accuracy. The potential model approach is a valuable task, which has given to us many insights into the nature of both heavy and light quarkonium physics. Using a relativistic approach together with a model in which and thresholds are taken into account, detailed analysis can be made on various aspects of heavy quarkonium.

Data Availability

No data were used to support this study.

Conflicts of Interest

The author declares that they have no conflicts of interest.

References

  1. T. Kawanai and S. Sasaki, “Charmonium potential from full lattice QCD,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 85, no. 9, 2012. View at Publisher · View at Google Scholar
  2. A. Le Yaouanc, L. Oliver, O. Pène, and J. C. Raynal, “Spontaneous breaking of chiral symmetry for confining potentials,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 29, no. 6, pp. 1233–1257, 1984. View at Publisher · View at Google Scholar
  3. S. Bolognesi, K. Konishi, M. Shifman, and D. Phys. Rev, “Patterns of symmetry breaking in chiral QCD,” Physical Review D, Covering Particles, Fields, Gravitation, and Cosmology, vol. 97, no. 9, Article ID 094007, 2018. View at Google Scholar
  4. T. Appelquist and H. D. Politzer, “Heavy Quarks and,” Physical Review Letters, vol. 34, no. 1, pp. 43–45, 1975. View at Publisher · View at Google Scholar
  5. A. De Rújula and S. L. Glashow, “Is bound charm found?” Physical Review Letters, vol. 34, no. 1, pp. 46–49, 1975. View at Publisher · View at Google Scholar · View at Scopus
  6. E. Eichten, K. Gottfried, T. Kinoshita, K. D. Lane, and T. M. Yan, “Interplay of confinement and decay in the spectrum of charmonium,” Physical Review Letters, vol. 36, p. 500, 1976. View at Publisher · View at Google Scholar
  7. D. P. Stanley and D. Robson, “Nonperturbative potential model for light and heavy quark-antiquark systems,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 21, no. 11, pp. 3180–3196, 1980. View at Publisher · View at Google Scholar · View at Scopus
  8. S. Godfrey and N. Isgur, “Mesons in a relativized quark model with chromodynamics,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 32, no. 1, pp. 189–231, 1985. View at Publisher · View at Google Scholar
  9. L. P. Fulcher, “Matrix representation of the nonlocal kinetic energy operator, the spinless Salpeter equation and the Cornell potential,” Physical Review D, vol. 50, p. 447, 1994. View at Google Scholar
  10. D. Ebert, R. N. Faustov, and V. O. Galkin, “Properties of heavy quarkonia and,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 67, no. 1, 2003. View at Publisher · View at Google Scholar
  11. T. Barnes, S. Godfrey, and E. S. Swanson, “Higher charmonia,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 72, article 054026, 2005. View at Publisher · View at Google Scholar
  12. O. Lakhina and E. S. Swanson, “Dynamic properties of charmonium,” Physical Review D, vol. 74, 2006. View at Google Scholar
  13. B. Q. Li and K. T. Chao, “Higher charmonia and X, Y, Z states with screened potential,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 79, Article ID 013011, 2009. View at Publisher · View at Google Scholar
  14. B. Q. Li and K. T. Chao, “Bottomonium Spectrum with Screened Potential,” Communications in Theoretical Physics, vol. 52, pp. 653–661, 2009. View at Google Scholar
  15. K. B. V. K. Bhaghyesh and A. P. Monteiro, “Heavy quarkonium spectra and its decays in a nonrelativistic model with Hulthen potential,” Journal of Physics G: Nuclear and Particle Physics, vol. 38, Article ID 085001, 2011. View at Google Scholar
  16. J. Segovia, P. G. Ortega, D. R. Entem, and F. Fernandez, “Bottomonium spectrum revisited,” Physical Review D, vol. 93, 2016. View at Publisher · View at Google Scholar
  17. S. Godfrey, K. Moats, and E. S. Swanson, “B and Bs meson spectroscopy,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 94, no. 5, 2016. View at Publisher · View at Google Scholar
  18. W.-J. Deng, H. Liu, L.-C. Gui, and X.-H. Zhong, “Spectrum and electromagnetic transitions of bottomonium,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 95, Article ID 074002, 2017. View at Google Scholar
  19. D. B. Lichtenberg, “Excited quark production at Hadron colliders,” International Journal of Modern Physics A, vol. 2, no. 6, pp. 1669–1705, 1987. View at Publisher · View at Google Scholar
  20. N. Brambilla, Quarkonium Working Group, CERN Yellow Report, 2005, arXiv:hep-ph/0412158.
  21. E. Eichten, K. Gottfried, T. Kinoshita, J. Kogut, K. D. Lane, and T.-M. Yan, “Spectrum of charmed quark-antiquark bound states,” Physical Review Letters, vol. 34, no. 6, pp. 369–372, 1975. View at Publisher · View at Google Scholar
  22. J. Richardson, “The heavy quark potential and the ϒ, J/ψ systems,” Physics Letters B, vol. 82B, p. 272, 1979. View at Google Scholar
  23. Y.-Q. Chen and Y.-P. Kuang, “Improved QCD-motivated heavy-quark potentials with explicit ΛMS dependence,” Physical Review D, vol. 46, p. 1165, 1992. View at Google Scholar
  24. Y.-Q. Chen and Y.-P. Kuang, “Erratum: “Improved QCD-motivated heavy-quark potentials with explicit ΛMS dependence”,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 47, p. 35, 1993. View at Publisher · View at Google Scholar · View at MathSciNet
  25. V. Khruschev, V. Savrin, and S. Semenov, “On the parameters of the QCD-motivated potential in the relativistic independent quark model,” Physics Letters B, vol. 525, no. 3-4, pp. 283–288, 2002. View at Publisher · View at Google Scholar
  26. A. Martin, “A simultaneous fit of bb, cc, ss (bcs Pairs) and cs spectra,” Physics Letters B, vol. 100, p. 511, 1981. View at Publisher · View at Google Scholar
  27. M. Machacek and Y. Tomozawa, “ψ Phenomenology and the nature of quark confinement,” Annals of Physics, vol. 110, p. 407, 1978. View at Google Scholar
  28. G. Fogleman, D. B. Lichtenberg, and J. G. Wills, “Heavy-meson spectra calculated with a one-parameter potential,” Lettere al Nuovo Cimento, vol. 26, p. 369, 1979. View at Google Scholar
  29. C. Quigg and J. L. Rosner, “Quantum mechanics with applications to quarkonium,” Physics Reports, vol. 56, no. 4, pp. 167–235, 1979. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  30. C. Quigg and J. L. Rosner, “Quarkonium level spacings,” Physics Letters B, vol. 71B, pp. 153–157, 1977. View at Google Scholar
  31. S. Xiaotong and L. Hefen, “A new phenomenological potential for heavy quarkonium,” Zeitschrift für Physik C Particles and Fields, vol. 34, no. 2, pp. 223–231, 1987. View at Publisher · View at Google Scholar · View at Scopus
  32. C. Quigg, H. B. Thacker, and J. L. Rosner, “Constructive evidence for flavor independence of the quark-antiquark potential,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 21, Article ID 3393, p. 234, 1980. View at Google Scholar
  33. W. Buchmüller, Y. Jack Ng, and S.-H. H. Tye, “Hyperfine splittings in heavy-quark systems,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 24, no. 12, pp. 3312–3314, 1981. View at Publisher · View at Google Scholar
  34. C. Bachas, “Concavity of the quarkonium potential,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 33, no. 9, pp. 2723–2725, 1986. View at Publisher · View at Google Scholar
  35. S. K. Choi, “Belle collaboration,” Physical Review Letters, vol. 91, Article ID 262001, 2003. View at Google Scholar
  36. X.-W. Kang and J. A. Oller, “Different pole structures in line shapes of the X(3872),” The European Physical Journal C, vol. 77, no. 6, 2017. View at Google Scholar
  37. R. R. Horgan, “The construction and classification of wavefunctions for the harmonic oscillator model of three quarks,” Journal of Physics G: Nuclear Physics, vol. 2, p. 625, 1976. View at Publisher · View at Google Scholar
  38. N. Isgur and G. Karl, “P-wave baryons in the quark model,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 18, no. 11, pp. 4187–4205, 1978. View at Publisher · View at Google Scholar
  39. A. W. Hendry, “Decays of high spin Δ* and N* resonances in the quark model,” Annals of Physics, vol. 140, no. 65, 1982. View at Publisher · View at Google Scholar · View at MathSciNet
  40. A. W. Hendry and D. B. Lichtenberg, “Properties of hadrons in the quark model,” Fortschritte der Physik/Progress of Physics banner, vol. 33, no. 3, pp. 139–231, 1985. View at Google Scholar
  41. K. B. V. Kumar, B. Hanumaiah, and S. Pepin, “Meson spectrum in a relativistic harmonic model with instanton-induced interaction,” The European Physical Journal A - Hadrons and Nuclei, vol. 19, no. 9, p. 247, 2004. View at Publisher · View at Google Scholar
  42. K. B. V. Kumar, M. Y.-L. Bhavyashri, and A. P. Monteiro, “P-wave meson spectrum in a relativistic model with instanton induced interaction,” International Journal of Modern Physics A, vol. 24, 2009. View at Google Scholar
  43. M. Tanabashi, “Particle Data Group,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 98, Article ID 030001, 2018. View at Google Scholar
  44. S. Godfrey and K. Moats, “Erratum:,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 92, no. 11, 2015. View at Publisher · View at Google Scholar
  45. V. Lengyel, Y. Fekete, I. Haysak, and A. Shpenik, “Calculation of hyperfine splitting in mesons using configuration interaction approach,” The European Physical Journal C, vol. 21, no. 2, pp. 355–359, 2001. View at Publisher · View at Google Scholar · View at Scopus
  46. R. Van Royen and V. F. Weisskopf, “Protsyessy raspada adronov i modyel cyrillic small soft sign kvarkov,” Il Nuovo Cimento A, vol. 50, no. 3, pp. 617–645, 1967. View at Publisher · View at Google Scholar · View at Scopus
  47. E. Braaten and S. Fleming, “QCD radiative corrections to the leptonic decay rate of the Bc meson,” Physical Review D, vol. 52, p. 181, 1995. View at Google Scholar
  48. W. Kwong, P. B. Mackenzie, R. Rosenfeld, and J. L. Rosner, “Quarkonium annihilation rates,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 37, no. 11, pp. 3210–3215, 1988. View at Publisher · View at Google Scholar
  49. J. T. Laverty, S. F. Radford, and W. W. Repko, “∖ga∖ga and g g decay rates for equal mass heavy quarkonia,” High Energy Physics - Phenomenology, 2009, arXiv:hep-ph/0901.3917v3. View at Google Scholar
  50. M. Ablikim, J. Z. Bai, Y. Ban, X. Cai, H. F. Chen, and H. S. Chen, “Determination of the ψ(3770), ψ(4040), ψ(4160) and ψ(4415) resonance parameters,” Physics Letters B, vol. 660, no. 3, pp. 315–319, 2008. View at Publisher · View at Google Scholar
  51. E. Eichten, K. Gottfried, T. Kinoshita, K. D. Lane, and T. -. Yan, “Charmonium: The model,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 17, no. 11, pp. 3090–3117, 1978. View at Publisher · View at Google Scholar
  52. C. Hayne and N. Isgur, “Beyond the wave function at the origin: Some momentum-dependent effects in the nonrelativistic quark model,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 25, no. 7, pp. 1944–1950, 1982. View at Publisher · View at Google Scholar
  53. N. Isgur, D. Scora, B. Grinstein, and M. B. Wise, “Semileptonic B and D decays in the quark model,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 39, no. 3, pp. 799–818, 1989. View at Publisher · View at Google Scholar · View at Scopus
  54. Z.-H. Wang, Y. Zhang, L.-B. Jiang, T.-H. Wang, Y. Jiang, and G.-L. Wang, “The strong decays of X(3940) and X(4160),” European Physical Journal C, vol. 77, no. 1, article 43, 2017. View at Publisher · View at Google Scholar