Advances in High Energy Physics
Volume 2018 (2018), Article ID 4549705, 7 pages
https://doi.org/10.1155/2018/4549705

Research Article

Nonrelativistic Arbitrary -States of Quarkonium through Asymptotic Iteration Method

Hakan Ciftci¹ and Hasan Fatih Kisoglu²

¹Gazi Üniversitesi, Fen-Edebiyat Fakültesi, Fizik Bölümü, 06500 Teknikokullar-Ankara, Turkey
²Department of Basic Sciences, Faculty of Maritime, Mersin University, Mersin, Turkey

Correspondence should be addressed to Hakan Ciftci; hciftci@gazi.edu.tr and Hasan Fatih Kisoglu; hasanfatihk@mersin.edu.tr

Received 15 February 2018; Accepted 27 March 2018; Published 29 May 2018

Academic Editor: Chun-Sheng Jia

Copyright © 2018 Hakan Ciftci and Hasan Fatih Kisoglu. 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 SCOAP³.

Abstract.

The energy eigenvalues with any states and mass of heavy quark-antiquark system (quarkonium) are obtained by using Asymptotic Iteration Method in the view of nonrelativistic quantum chromodynamics, in which the quarks are considered as spinless for easiness and are bounded by Cornell potential. A semianalytical formula for energy eigenvalues and mass is achieved via the method in scope of the perturbation theory. The accuracy of this formula is checked by comparing the eigenvalues with the ones numerically obtained in this study and with exact ones in literature. Furthermore, semianalytical formula is applied to c, b, and c meson systems for comparing the masses with the experimental data.

1. Introduction

Investigation of an atomic or subatomic system is done by achieving an energy spectrum of the system. This is carried out, generally, for the events in which the system is bounded by a potential function. Besides, the scattering states or resonance states can also be observed in the investigation of the system. The eigenvalues (or eigenenergies) of Hamiltonian of this system is obtained for a given potential function. In order to do this, various mathematical methods are used in quantum mechanics. One of these, named Asymptotic Iteration Method (AIM), has been commonly used since 2003 [1]. AIM can be used for analytically as well as numerically (or approximately) solvable problems [2–4]. Moreover, it can be used for obtaining the perturbative energy eigenvalues of the system without any need of the unperturbative eigenstate [5, 6].

As a subatomic system, a quarkonium that is composed of a heavy quark-antiquark (q) pair has attracted attention of particle physicists since the first half of 1970, and [7–11] are just a few studies of them. In most of these studies, for easiness, the system is examined via Schrödinger equation in nonrelativistic quantum chromodynamics (NRQCD), assuming that the quarks are spinless [12–15]. Cornell potential is one of the potential functions that represent interactions between the quarks in such a q system. It is used for obtaining the mass and energy spectrum of the quarkonium and obtaining the hadron decay widths [7–9, 16]. Cornell potential is given as

where and are positive constants ( is in energy dimension). As is seen in (1), Cornell potential has two parts: one is the Coulombic term and the other is the linear part. For obtaining the energy levels and mass of the quarkonium, and may be fitted to the first-few states. Therefore, the full spectrum of the quarkonium can be constructed through these potential parameters.

In literature, it is possible to find many studies in which the solutions of Schrödinger equation for Cornell potential have been obtained. For example, in [17], Hall has found an approximate energy formula to construct an energy spectrum of Schrödinger equation for Cornell potential, under some conditions. Jacobs et al. [13] have compared the eigenvalues of Schrödinger and spinless Salpeter equations in the cases of Cornell potential and Wisconsin potential [18]. Vega and friends have obtained, for l=0 states, the energy spectrum, mass, and wavefunctions at the origin for c, b, and b mesons by using the usual variation method in the scope of supersymmetric quantum mechanics (SUSYQM) [19, 20], in [12]. They have also compared their results with the exact ones in literature and with the experimental data.

In this study, we attempted to get the energy eigenvalues (for any states) and masses of heavy mesons by using Asymptotic Iteration Method in the view of NRQCD, in which the quarks are considered as spinless for easiness and are bounded by Cornell potential. We achieved a semianalytical formula for constructing the energy spectrum and obtaining the masses of the mesons, using the method in scope of the perturbation theory. The accuracy of this formula was cross-checked by comparing the eigenvalues with the ones numerically obtained in this study and with the exact ones in literature. Furthermore, semianalytical formula was applied to c, b, and c heavy mesons for comparing the masses with the experimental data.

AIM has been firstly applied to Schrödinger equation for Cornell potential by Hall and Saad in [21]. They have used Airy function as an asymptotic form of the wavefunction and have got highly-accurate numerical results in their study. Alternatively, we obtained a semianalytical mass-energy formula for quarkonium by having differential equation which gives polynomial solutions for asymptotic forms of the wavefunction of the system.

This paper is organized as follows: we give a short summary of AIM in Section 2, while Section 3 includes the main problem. In Section 4, we give numerical results for the eigenenergies and obtain semianalytical energy formula by applying perturbation theory to our problem in the view of AIM. Furthermore, in Section 4, we compare our energy spectrum and masses with the exact ones in literature and with the experimental data. Finally, Section 5 includes some comments about our results.

2. The Asymptotic Iteration Method (AIM)

According to the organization of the paper, we summed up AIM in this section, while it is comprehensively introduced in [1]. The AIM is used to solve second-order homogeneous linear differential equations in the following form:

where and have continuous derivatives in the defined interval of the independent variable. If there is an asymptotic condition such as

for , where is large enough, the general solution of (2) is obtained as

with the functions

As a field of application, AIM can be used to deal with Schrödinger equation (or energy eigenvalue problem) in mathematical physics. The eigenvalues can be obtained through the following quantization condition:

If the energy eigenvalues () can be obtained from (6), independently from the variable, the problem is exactly solvable. In this case, the eigenvalue and eigenfunction of th energy level can be derived in explicit algebraic form via iterations. However, there are limited numbers of suitable potentials for this case.

As for the approximately (or numerically) solvable problems, depends on both and . In this case, an appropriate value, , should be determined to solve with respect to [2, 6]. The energy eigenvalue of an th level is obtained through iterations where .

3. Formulation of the Problem

Consider the following Cornell potential:

where , are real and positive constants ( is in energy dimension) and . If we substitute into Schrödinger equation in three dimensions, we have

where , , and . and are energy eigenvalue of th level and reduced mass of the q system, respectively ( and are quark masses). Besides, we study in natural units (i.e., , ) for the system. After changing the variable, in (8), as , then substituting , we get

If one puts into (9), in accordance with the domain of the problem, we have

where , , , and . The final equation is suitable for applying AIM. After this point, we can apply AIM to the problem in two different ways: one is direct application (i.e., approximate solution) to get the numerical results and the other is usage of the method in scope of perturbation theory to obtain perturbative energies through a perturbation expansion as follows:

where , , ,… are perturbation expansion coefficients. These can be obtained independently from the potential parameters. Thus, we can get a semianalytical formula for the energy eigenvalues. One can also achieve the mass-energy of the system by using this formula, as given in Section 4.

3.1. Numerical Results

In this section, we directly apply AIM to (10) to get the energy eigenvalues for different potential parameters, and we compare our results with the perturbative energies, for which (28) in the next section has been used.

From (12), it is easily seen that and according to (2). We tabulate the results of direct application of AIM in Tables 1, 2, and 3. For simplicity, in the calculations, the reduced mass has been considered GeV. In Table 1 the potential parameters have been chosen as and GeV while , GeV in Table 2, and , GeV in Table 3. , seen in the tables, is for the comparison and has been obtained by using (28).

Table 1: Comparisons of the perturbative energy eigenvalues () (in GeV) with those obtained by direct application of AIM () for the potential parameters and GeV, and for the reduced mass GeV.

En0
Epert

E0l
Epert

0
1.39788
1.41015
0
1.39788
1.41015

1
3.47509
3.47509
1
2.82565
2.82690

2
5.03291
5.03224
2
3.85058
3.85089

3
6.37015
6.36948
3
4.72675
4.72687

4
7.57493
-
4
5.51698
-

5
8.68791
-
5
6.24840
-

Table 2: Comparisons of the perturbative energy eigenvalues () (in GeV) with those obtained by direct application of AIM () for the potential parameters, and GeV, and for the reduced mass GeV.

En0
Epert
l
E0l
Epert

0
-0.221031
-0.164433
0
-0.221031
-0.164433

1
0.0347222
0.033627
1
0.0174006
0.023501

2
0.141913
0.138477
2
0.102472
0.104008

3
0.220287
0.217229
3
0.159831
0.160406

4
0.286111
-
4
0.206238
-

5
0.344602
-
5
0.246681
-

Table 3: Comparisons of the perturbative energy eigenvalues () (in GeV) with those obtained by direct application of AIM () for the potential parameters, and GeV, and for the reduced mass GeV.

n
En0
Epert
l
E0l
Epert

0
46.4022
46.4047
0
46.4022
46.4047

1
85.3393
85.3394
1
70.0161
70.0165

2
116.729
116.729
2
89.7154
89.7154

3
144.315
144.315
3
107.334
107.334

4
169.461
-
4
123.562
-

5
192.851
-
5
138.761
-

As can be seen from Tables 1–3, the perturbative energy eigenvalues are in very good agreement with the numerically obtained ones, even for small values of the parameter . Furthermore, they are in accordance with each other for GeV, while (see in Tables 1 and 3). Additionally, this agreement is much better for higher quantum states. The perturbative eigenvalues are a little bit different from that obtained as numerically, for GeV, , and the lower quantum states (see in Table 2). However, they are in agreement for the higher levels.

4. Perturbation Theory

Although the usage of perturbation method in the frame of AIM is comprehensively introduced in [5], we give a summary about the methodology in this section, assuming that the potential of a system is written as

where is solvable (unperturbed Hamiltonian) potential. and are potential of the perturbed Hamiltonian and perturbation expansion parameter, respectively. The Schrödinger equation then reads

where eigenvalues are written as a series expansion of th-order correction as follows:

After substituting in (14), one can obtain the following equation for :

and the termination condition in this case can be written as

Once is expanded about , we obtain

where .

According to perturbation method in the framework of AIM, solving the equation with respect to (unknown) gives (eigenvalues of unperturbed Hamiltonian), equation gives (first-order correction to ), gives (second-order correction to ), and so on. Besides, the perturbative eigenfunctions can be achieved in the same vein with the eigenvalues. This is an alluring feature of the AIM usage in the perturbation theory for obtaining the eigenfunctions given as follows:

where . is expanded about in a similar manner, done for obtaining the eigenvalues. So,

where . Thus, perturbation expansion of is written as follows:

where th-order correction to is

4.1. Perturbation Theory for the Cornell Potential

For our problem, we may apply the perturbation expansion which has been elucidated in previous section to the following differential equation:

Suppose that is written as follows:

where is the perturbation expansion parameter. So, the energy eigenvalue is yielded as

and more clearly

In the above expansion, the general form of the zeroth-order correction is obtained via

The first-order correction, , is obtained by using the equation in the same manner with , while is used for . Numerical results of , , and coefficients, obtained by AIM, are reported in Table 4 for some energy levels. Besides, for GeV, comparisons of the perturbative energy eigenvalues with the ones obtained by direct application of AIM have been given in Tables 1, 2, and 3, in previous section. We emphasize, in Table 4, that corrections to the perturbation expansion do not depend on the potential parameters.

Table 4: Perturbation coefficients of the expansion given as (24) and (26). Notice that corrections to the perturbation expansion do not depend on the potential parameters.

l
n
σ₀(n,l)
σ₁(n,l)
σ₂(n,l)

0
0
3.71151
-0.525933
-0.0232729

1
6.48922
-0.366743
-0.00767365

2
8.76334
-0.297538
-0.00400191

3
10.7732
-0.256486
-0.00251618

1
0
5.33566
-0.322683
-0.00554189

1
7.75358
-0.258925
-0.00282569

2
9.85399
-0.222298
-0.00176295

3
11.7558
-0.197751
-0.00122526

2
0
6.74357
-0.244191
-0.00241586

1
8.93661
-0.208300
-0.00148846

2
10.9037
-0.184664
-0.00102765

3
12.7146
-0.167585
-0.000761053

3
0
8.01784
-0.200753
-0.00134507

1
10.0516
-0.177251
-0.000921458

2
11.9129
-0.160449
-0.000679139

3
13.6471
-0.147666
-0.000525832

As a practice, we have applied our perturbation expansion formula (up to second-order correction) to get the ground-state energies of quarkonium in Table 5, for various values of the parameter , while GeV and GeV. In Table 5, we also report comparisons of the perturbative energy eigenvalues with the ones of s-wave heavy quarkonium from [10, 21].

Table 5: Comparisons of energy eigenvalues (in GeV) obtained by using the perturbation expansion formula in (26) () with the ones of s-wave heavy quarkonium from [10, 21]. The potential parameter is taken as GeV, while the reduced mass is GeV in this case. The eigenvalues of [10, 21] are exact results.

A
E₀₀ (Ref. [10])
E₀₀ (Ref. [21])
Epert
A
E₀₀ (Ref. [10])
E₀₀ (Ref. [21])
Epert

0.2
2.16732
2.16732
2.16741
0.1
2.25368
2.25368
2.25369

0.4
1.98850
1.98850
1.98923
0.3
2.07895
2.07895
2.07927

0.6
1.80107
1.80107
1.80367
0.5
1.89590
1.89590
1.89740

0.8
1.60441
1.60441
1.61063
0.7
1.70394
1.70393
1.70808

1
1.39788
1.39788
1.41015
0.9
1.50242
1.50242
1.51132

1.2
1.18084
1.18083
1.20221
1.1
1.29071
1.29071
1.30711

1.4
0.95264
0.95264
0.98683
1.3
1.06817
1.06817
1.09545

1.6
0.71266
0.71266
0.76400
1.5
0.83416
0.83416
0.87635

1.8
0.46027
0.46026
0.53373
1.7
0.58805
0.58805
0.64980

As is seen from Table 5, the results for which our perturbation expansion (up to second-order correction) has been used are in very good agreement with [10, 21] for small values of . However, our analytical results are a little bit different from the exact ones as gets larger values. It seems that the perturbation expansion, which includes third-order correction, may give more accurate results. The more correction term we add to the perturbative expansion, the more compatible results we get. Nevertheless, we can say that (26) can be used as an eigenvalue formula of the Schrödinger equation in case of Cornell potential, for practical purposes. So, one can use the following formula:

for obtaining the eigenvalues and mass of the quarkonium for Cornell potential. Besides, it can be fit to mass formula of experimental values for determining the potential parameters and . The advantage of (28) is that the coefficients , , and are independent of the potential parameters.

4.2. Energy Eigenvalues and Mass Spectrum for Heavy Quarkonium

In this section, we tested our formula through cross-checking with the exact results in literature and with the experimental data. For comparing our energy eigenvalues with the exact ones, the parameters of Cornell potential have been considered and GeV. Besides, we have chosen the quark masses as GeV and GeV, in this case [12].

Also, we tested our formula by comparing our results, for the masses of heavy mesons, with the experimental data. For doing this, we have taken the quark masses as GeV and GeV and the potential parameters as and GeV. All these values have been obtained by fitting our formula to the experimental data in [22].

In Table 6, we compared our energy eigenvalues calculated by using (28) with the ones of [12]. Furthermore, in Table 7, we gave our results for the masses of the mesons obtained by the same equation. Table 7 also includes the experimental data got from [22].

Table 6: Comparisons of the energy eigenvalues (in GeV) of the mesons c, b, and b calculated by using (28) with the exact ones of [12]. The parameters of Cornell potential are and GeV, while the quark masses are   =1.84 GeV and   =5.18 GeV.

c
b
b

En
Exact [12]
Ref. [12]
AIM
Exact [12]
Ref. [12]
AIM
Exact [12]
Ref. [12]
AIM

1s
0.2575
0.2578
0.2660
-0.1704
-0.1702
-0.1216
0.1110
0.1113
0.1269

2s
0.8482
0.8096
0.8481
0.4214
0.3579
0.4203
0.6813
0.6324
0.6803

3s
1.2720
1.1427
1.2715
0.7665
0.5612
0.7635
1.0686
0.9065
1.0668

Table 7: Comparisons of the masses (in GeV), obtained via AIM, of the heavy mesons c, b, and c with the ones of [12], and with the experimental data from [22]. In this case, we have taken the quark masses as   =1.44 GeV and   =4.87 GeV, and the potential parameters as and GeV, for our calculations. All these parameters have been obtained by fitting our formula, given in (28), to the experimental data.

c
b
c

Mn
Exp.
Ref. [12]
AIM
Exp.
Ref. [12]
AIM
Exp.
Ref. [12]
AIM

1s
3.097
3.097
3.096
9.460
9.350
9.462
6.275
6.291
6.362

2s
3.686
3.649
3.672
10.023
9.878
10.027
6.842
6.812
6.911

3s
4.039
3.963
4.085
10.355
10.081
10.361
-
7.087
7.284

4s
-
-
4.433
10.579
-
10.624
-
-
7.593

1p
3.511
-
3.521
9.899
-
9.963
-
-
6.792

2p
3.927
-
3.951
10.260
-
10.299
-
-
7.178

3p
-
-
4.310
10.512
-
10.564
-
-
7.494

1d
-
-
3.800
10.164
-
10.209
-
-
7.051

It can be seen from Table 6 that the energy eigenvalues of the mesons c, b, and b, obtained by (28), are more compatible with the exact ones, than those of [12]. The difference between AIM and [12] becomes clearer as the energy level increases. Similar things can be said for the masses in Table 7: the results obtained via AIM are closer to the experimental data than those of [12].

5. Conclusion

We have used AIM to obtain both, the eigenvalues of Schrödinger equation and mass of q system for Cornell potential, in three dimensions. AIM has some advantages such as being used for either exactly or numerically (or approximately) solvable problems. Furthermore, one can use AIM in the frame of perturbation theory. Once it is performed to obtain perturbative solutions, the wavefunction of unperturbed Hamiltonian is not needed to get the corrections to the perturbation expansion.

In the present study, the energy eigenvalues in the case of Cornell potential have been achieved by direct application of the method. Besides, we have performed perturbation theory in the view of AIM for the problem and found a semianalytical formula for energy eigenvalues. Numerical results obtained by using this formula, for the reduced mass GeV, conform with the exact results of [10, 21], in a wide spectrum of the potential parameters and (especially for ). Furthermore, the results are compatible with the ones obtained directly, in Section 3. It is also possible to see from the results that the perturbative eigenvalues fit in with the exact ones for higher quantum states, even for the large values of . For any values of and , the higher quantum states are more consonant with the exact ones than the lower states. The perturbation expansion, which includes third-order correction, may give more accurate results. The more correction terms we add to the perturbative expansion, the more compatible results we may get.

We have also tested our semianalytical formula, by cross-checking it with the exact results in literature and with the experimental data. It can be seen, from Table 6, that our energy eigenvalues calculated by using (28) are more compatible with the exact ones than those of [12]. Furthermore, the difference between our results and [12] becomes clearer as the energy level increases. By using AIM, we have also obtained mass results which are closer to the experimental data than [12].

As a consequence, semianalytical formula achieved for energy eigenvalues and mass of quarkonium can be used for practical purposes in the case of Cornell potential. If our formula is fitted to the experimental data, the potential parameters (and masses of the quarks, if it is needed) can also be obtained.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This academic work was supported by the Mersin Technology Transfer Office Academic Writing Center of Mersin University.

References

H. Çiftçi, R. L. Hall, and N. Saad, “Asymptotic iteration method for eigenvalue problems,” Journal of Physics A: Mathematical and General, vol. 36, no. 47, pp. 11807–11816, 2003.

N. Saad, R. L. Hall, and N. Saad, “Exact and approximate solutions of Schrödinger's equation for a class of trigonometric potentials,” Central European Journal of Physics, vol. 11, no. 1, pp. 37–48, 2013.

E. Olğar, R. Koç, and H. Tütüncüler, “The exact solution of the s-wave Klein-Gordon equation for the generalized Hulthén potential by the asymptotic iteration method,” Physica Scripta, vol. 78, no. 1, Article ID 015011, 2008.

F. M. Fernández, “On an iteration method for eigenvalue problems,” Journal of Physics A: Mathematical and General, vol. 37, no. 23, pp. 6173–6180, 2004.

H. Ciftci, R. L. Hall, and N. Saad, “Perturbation theory in a framework of iteration methods,” Physics Letters A, vol. 340, no. 5-6, pp. 388–396, 2005.

H. Ciftci and H. F. Kisoglu, “Accidental Degeneracies in N dimensions for Potential Class αr2d−2−βrd−2 via Asymptotic Iteration Method (AIM),” Communications in Theoretical Physics, vol. 67, p. 350, 2017.

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.

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.

E. Eichten, K. Gottfried, T. Kinoshita, K. D. Lane, and T. M. Yan, “Charmonium: comparison with experiment,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 21, no. 1, pp. 203–233, 1980.

H. Chung S and J. Lee, “Cornell potential parameters for S-wave heavy quarkonia,” Journal of the Korean Physical Society, vol. 52, p. 1151, 2008.

M. Hamzavi and A. A. Rajabi, “Scalar–vector–pseudoscalar Cornell potential for a spin-1/2 particle under spin and pseudospin symmetries: 1 + 1 dimensions,” Annals of Physics, vol. 334, pp. 316–320, 2013.

A. Vega and J. Flores, “Heavy quarkonium properties from Cornell potential using variational method and supersymmetric quantum mechanics,” Pramana – Journal of Physics, vol. 87, Article ID 73, 2016.

S. Jacobs, M. G. Olsson, and C. Suchyta, “Comparing the Schrödinger and spinless Salpeter equations for heavy-quark bound states,” Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 33, no. 11, pp. 3338–3348, 1986.

B. Grinstein, “A modern introduction to quarkonium theory,” International Journal of Modern Physics A, vol. 15, p. 461, 2000.

W. Lucha, F. Schoberl, and D. Gromes, “Bound states of quarks,” Physics Reports, vol. 200, no. 4, pp. 127–240, 1991.

P. Evans, C. Allton, and J. Skullerud, Physical Review D: Particles, Fields, Gravitation and Cosmology, vol. 89, no. 7, 2014.

R. L. Hall, “Simple eigenvalue formula for the Coulomb-plus-linear potential,” Physical Review D, vol. 30, p. 433, 1984.

K. Hagiwara J and S. Jacobs, “Probing QCD with heavy quark bound states,” Physics Letters B, vol. 131, pp. 455–460, 1983.

E. Gozzi, M. Reuter, and W. Thacker D, “Variational methods via supersymmetric techniques,” Physics Letters A, vol. 183, p. 29, 1993.

F. Cooper, A. Khare, and U. Sukhatme, “Supersymmetry and quantum mechanics,” Physics Reports, vol. 251, no. 5-6, pp. 267–385, 1995.

L. R. Hall and S. Nasser, “Schrödinger spectrum generated by the Cornell potential,” Open Physics, vol. 13, p. 83, 2015.

C. Patrignani and etal. (Particle Data Group), “Review of Particle Physics,” Chinese Physics C, vol. 40, 2016, 100001.