To obtain closed-form solutions for the radial Schrödinger wave equation with non-solvable potential models, we use a simple, easy, and fast perturbation technique within the framework of the asymptotic iteration method (PAIM). We will show how the PAIM can be applied directly to find the analytical coefficients in the perturbation series, without using the base eigenfunctions of the unperturbed problem. As an example, the vector Coulomb and the harmonic oscillator plus linear scalar potential parts implemented with their counterpart spin-dependent terms are chosen to investigate the meson sectors including charm and beauty quarks. This approach is applicable in the same form to both the ground state and the excited bound states and can be easily applied to other strongly non-solvable potential problems. The procedure of this method and its results will provide a valuable hint for investigating tetraquark configuration.
One of the challenging problems in non-relativistic quantum mechanics is to find closed-form solutions to the radial Schrödinger wave equation for non-solvable potentials that are used in different fields of physics. In particular, the perturbed Coulomb potentials are representing the most important potentials in many scientific models, especially in particle physics .
As an example, the experimental observations of the exotic states in the charm sector with content at Belle  and BESIII , as well as the structures of the form in ref. [4,5], have opened new interest to re-study the hadronic structures and their spectroscopy with different phenomenological potential models.
However, the lack of closed-form solutions to the radial Schrödinger wave equation with many potential models makes the study of such potentials one of the most popular theoretical laboratories for examining the validity of the various approximation techniques based on perturbative and non-perturbative approaches [6,7,8, 9,10,11].
Many of the existing theoretical work in the literature can obtain the hadron masses by solving numerically the two-body radial Schrödinger wave equation with Cornell-inspired potentials, and then the respective spin-dependent interaction terms are added and treated by the standard perturbation technique . Nevertheless, the computational schedule in those works is very complicated and can be hardly extended to study other potential models.
On this basis, the present work is aimed to solve the radial Schrödinger wave equation for heavy quark–antiquark structures with the vector Coulomb potential plus linear and harmonic scalar potentials parts, implemented in the same form with their counterpart spin-dependent terms.
We will show how the PAIM can be used directly to find analytically the coefficients in the energy perturbation series as a solution to the radial Schrödinger wave equation, without the need of using the base eigenfunctions of the unperturbed problem. Furthermore, PIAM is unlike the other perturbation methods, with no constraints on the coupling constants or the quantum numbers involved in the phenomenological potential, and even it avoids sticking to numerical computations at the early stages of its procedure. Besides, PIAM can yields very accurate and rapidly converging eigenvalues, and it is applicable in the same form to both the ground and excited energy states.
With this in mind, this article is organized as follows. Section 2 is devoted to the theoretical framework of the phenomenological potential model. We presented the static potential, which is composed of scalar and vector parts, implemented with their counterpart spin-dependent correction terms. In Section 3, the formulation of the PAIM for finding the energy eigenvalues is given. The analytical expressions for PAIM are cast in such a way that allows the reader to use them without proceeding into their derivation. In Section 4, as an example, we tested numerically the method to find some of the P-wave mass spectrums of the charmonium , and bottonium mesons, and therein for comparison we presented the available experimental and theoretical results in the literature. Finally, we give a summary and concluding remarks for possible future works on tetraquark masses.
It is generally believed that the phenomenological potential models used in the hadronic spectroscopy with the Coulomb potential perturbed by a term or terms involving various powers of seem to provide the clearest and most accurate picture for meson states. The Coulomb potential arises from the one-gluon exchange associated with a Lorentz vector structure, and the various powers of parts in the potential are responsible for the confinement. These parts are known as the static spin-independent potential parts. For our further discussion, we used a fairly flexible mixing type of such potentials. We took the potential interaction parts as a sum of vector and scalar parts, plus the spin-dependent terms:
To implement the spin-dependent terms within the pre-assumed phenomenological potential parts, we based on the Breit-Fermi one-gluon exchange Hamiltonian , where the spin-dependent part for quark masses and have three distinctive types of radial-interaction parts known as the spin–spin , the spin–orbit , and the tensor part , respectively:
The expectation values of the spin–spin interaction operator in equation (4) are calculated in terms of the spin quantum numbers: , where and are the spins of the two quarks and , respectively, and stands for the total spin of the structure under consideration. Since each quark (anti-quark) has spin , thus they can couple into singlets with or triplet with .
On the other hand, the expectation values of the spin–orbit interaction operator are calculated in terms of the total angular momentum quantum number , defined by the vector sum . For S-wave states ( ) and/or for the singlet state case , the spin–orbit term is always zero. However, for the triplet state cases with , we have , which, in general, can yield a non-zero value. Similarly, the expectation values of the tensor operator are calculated in terms of the total spin , and the orbital angular momentum quantum number , and it is given by ref. :
The radial Schrödinger wave equation that describes the internal structure between any quark–antiquark with the phenomenological potential model in equation (13) for (in atomic units, ) is as follows:
The solution of equation (16) within the framework of PAIM is started with the change of variables by letting , and further we removed the first derivative by employing the ansatz , provided that , which in turn implies a Schrödinger-like radial equation:
It should be noted here that the choice of the ansatz in equation (19) is motivated, that is, when we switch off, equation (19) is reduced to an exactly solvable eigenvalue problem within the framework of AIM [14,15]. On the other hand, if we switch on, a straightforward application of AIM to equation (19) does not give any reasonable results, where it is observed that the sequence appears not to converge, a result that violates the principle behind the AIM.
Accordingly, to have a closed-form solution for equation (19) within the framework of PAIM, we proposed be expressed as a sum of two parts:
Following the systematic procedure of PAIM, the energy eigenvalues are solved by imposing the asymptotic aspect of the method for sufficiently large iteration number :
To obtain the analytical leading energy term , one should simply switch off in equation (19), and for further simplifications, we recall the mathematical definitions of equations (21)–(23). Then the roots of leads to an exactly solvable eigenvalue in terms of :
To find the other analytical perturbed expansion terms, one should simply go to equation (19) switch on, therein, replace with , and terminate the iterations by imposing the condition = 0. The first root of the resulting equation gives . Proceeding exactly in the same way, one can obtain the second-order closed-form perturbed expansion terms and so on. The analytical expressions of this procedure are displayed in Table 1. Of course, in Table 1, and are changeable according to the total spin and total angular momentum of the structure state.
Consequently, the final closed-form solution for the energy eigenvalues in conjunction with equation (19) within PAIM framework is as follows:
A further interesting quantity in this respect is the transformed eigenfunctions . Knowledge of will help to construct the radial eigenfunctions . Inspired by equation (18), given that , and following again exactly the same systematic procedure of PAIM explained above, by expanding this time :
Having solved the radial Schrödinger wave equation with the proposed phenomenological potential model, the corresponding mass of a particular quark–antiquark structure within this model can be written as:
Few P-wave mass spectrums of the charmonium are calculated and listed in Table 2. The calculated masses are found to be in good agreement with the recently observed ones. Moreover, with the same input parameters, in a unified way, we extend this study and we predicted the P-wave bottonium mass spectrums. In Table 3, we present our final findings on the P-wave mass spectrums of bottonium . In the tables, we reported also on the predictions of the other methods together with the experimental results. A direct comparison of our results with the experimental results reveals that they are in good agreement.
Evidently, our numerical analysis shows that the contribution of the non-perturbative terms is about 90%, while the maximum contribution of the perturbative terms is about 10% of the total sum.
Finally, although in this work we have tested with PIAM the P-wave mass spectrum of the potential model of equation (1), one can easily find the analytical and numerical solutions for other potential models using the same procedure. It is also important to emphasize that the PAIM described in this article and its analytical expressions are much more useful than pure numerical calculations.
The present work is focused on the theoretical implementation of the PAIM and to show how it is simple and easy to solve the radial Schrödinger wave equation with non-solvable potential models. Although we aimed not to make the fittings with the available experimental we only emphasized on the method and how it works, and we believe that the interested readers if necessary can do better fitting. Also, we aimed to show how it is an easy task to use the PAIM without having to worry about the ranges of the couplings in the introduced potential model. Moreover, the degree of precision of the results can also be drastically improved by raising the perturbative order in the energy expansion series one step more when it is necessary without any technical difficulties. We believe that this method and its results will provide a valuable future hint for investigating tetraquark configuration.
H. I. Alrebdi extends her appreciation to the Deanship of Scientific Research at Princess Nourah bint Abdulrahman University, where this research was funded by the Deanship of Scientific Research at Princess Nourah bint Abdulrahman University through the Fast-track Research Funding Program.
Conflict of interest: Authors state no conflict of interest.
 Roychoudhury RK, Varshni YP. Shifted 1/N expansion and exact solutions for the potential V(r)=−Z∕r+gr+λr2. J Phys A. 1988;21:3025. Search in Google Scholar
 Liu ZQ, Shen CP, Yuan CZ, Adachi I, Aihara H, Asner Dm, Liu ZQ, et al. Study of e+e−→π+π−J/ψ and observation of a charged charmonium-like state at Belle. Phys Rev Lett. 2013;110:252002 Search in Google Scholar
 Ablikim M, Achasov MN, Ai XC, Albayrak O, Ambrose DJ, An FF, et al. Observation of a charged charmonium like structure in e+e−→π+π−J/ψ at s=4.26GeV. Phys Rev Lett 2013;110:252001. Search in Google Scholar
 Brambilla N, Eidelman S, Heltsley BK, Vogt R, Bodwin GT, Eichten E, et al. Heavy quarkonium: progress, puzzles, and opportunities. Eur Phys J C. 2011;71:1534. Search in Google Scholar
 Bevan AJ, Golob B, Mannel T, Prell S, Yabsley BD, Yabsley H, et al. The Physics of the B Factories. Phys J C. 2014;74:3026. Search in Google Scholar
 Bessis D, Vrscay ER, Handy CR. Hydrogenic atoms in the external potential V(r)=gr+λr2: exact solutions and ground-state eigenvalue bounds using moment methods. J Phys A. 1987;20:419. Search in Google Scholar
 Datta DP, Mukherjee S. Analytic continued fraction technique for bound and confined states for a class of confinement potentials. J Phys A. 1980;13:3161. Search in Google Scholar
 Chaudhuri RN. False eigenvalues of the Hill determinant method. J Phys A. 1988;21:567. Search in Google Scholar
 Adhikari R, Dutt R, Varshni YP. Exact solutions for polynomial potentials using super symmetry inspired factorization method. Phys Lett A. 1989;141:1. Search in Google Scholar
 Roychoudhury RK, Varshni YP. Family of exact solutions for the Coulomb potential perturbed by a polynomial in r. Phys Rev A. 1990;42:184. Search in Google Scholar
 Ciftci H, Hall RL, Saad N. Construction of exact solutions to eigenvalue problems by the asymptotic iteration method. Phys Lett A. 2005;340:388. Search in Google Scholar
 Saxena RP, Varma VS. Some exact solutions for the potential −r−1+2λr+2λ2r2. J Phys A. 1982;15:L149. Search in Google Scholar
 Lucha W, Schöberl FF, Gromes D. Bound states of quarks. Phys Rep. 1991;200:127–240. Search in Google Scholar
 Ciftci H, Hall RL, Saad N. Asymptotic iteration method for eigenvalue problems. J Phys A 2003;36:11807. Search in Google Scholar
 Barakat T. The asymptotic iteration method for the eigenenergies of the Schrödinger equation with the potential V(r)=−Z∕r+gr+λr2. J Phys A. 2006;39:823. Search in Google Scholar
 Patrignani C, Agashe K, Aielli G, Amsler C, Antonelli M, Asner DM, et al. Review of particle physics. Chinese Phys C. 2016;40:100001. Search in Google Scholar
 Bhavsar T, Shah M, Vinodkumarc PC. Status of quarkonia-like negative and positive parity states in a relativistic confinement scheme. Eur Phys J C. 2018;78:227. Search in Google Scholar
 Shah M, Parmar A, Vinodkumar PC. Leptonic and Digamma decay Properties of S-wave quarkonia states. Phys Rev D. 2012;86:034015. Search in Google Scholar
 Li B-Q, Chao K-T. Higher charmonia and X, Y, Z states with screened potential. Phys Rev D. 2009;79:094004. Search in Google Scholar
 Barnes T, Godfrey S, Swanson ES. Higher charmonia. Phys Rev D. 2005;72:054026. Search in Google Scholar
© 2021 Haifa Ibrahim Alrebdi and Thabit Barakat, published by De Gruyter
This work is licensed under the Creative Commons Attribution 4.0 International License.