Relativistic Two-Dimensional Harmonic Oscillator Plus Cornell Potentials in External Magnetic and AB Fields

The Klein-Gordon (KG) equation for the two-dimensional scalar-vector harmonic oscillator plus Cornell potentials in the presence of external magnetic and Aharonov-Bohm (AB) flux fields is solved using the wave function ansatz method. The exact energy eigenvalues and the wave functions are obtained in terms of potential parameters, magnetic field strength, AB flux field, and magnetic quantum number. The results obtained by using different Larmor frequencies are compared with the results in the absence of both magnetic field ( 𝜔 𝐿 = 0) and AB flux field ( 𝜉 = 0 ) cases. Effect of external fields on the nonrelativistic energy eigenvalues and wave function solutions is also precisely presented. Some special cases like harmonic oscillator and Coulombic fields are also studied.


Introduction
The exact solution of Schrödinger equation (SE) and the relativistic wave equations for some physical potentials are very important in many fields of physics and chemistry since they contain all the necessary information for the quantum system under investigation. The hydrogen atom and the harmonic oscillator are usually given in textbooks as two of several exactly solvable problems in both classical and quantum physics [1]. The exact -state solutions of the SE are possible only for a few potentials, and hence approximation methods are used to obtain their solutions. According to the Schrödinger formulation of quantum mechanics, a total wave function provides implicitly all relevant information about the behaviour of a physical system. Hence, if it is exactly solvable for a given potential, the wave function can describe such a system completely. Until now, many efforts have been made to solve the stationary SE with anharmonic potentials in three dimensions (3D) and two dimensions (2D) [2][3][4][5][6][7] with many applications to molecular and chemical physics. However, the study of SE with some of these potentials in arbitrary dimension D is also solved in the following (cf. [8] and the references therein). The study of bound states is fundamental in the understanding of molecular spectrum of a diatomic molecule and provides us with insight into the physical problem under consideration in quantum mechanics [9].
It is well-known that the non-relativistic quantum mechanics is an approximate theory of the relativistic one. However, when a particle moves in a strong potential field, the relativistic effects must be considered which give the corrections for non-relativistic quantum mechanics [18]. Therefore, the motion of spin-0 and spin-1/2 particles satisfies the KG and the Dirac equations, respectively. In particular, solutions to relativistic equations play an important role in many aspects of modern physics. For instance, the Dirac equation has been used to explain the antinucleon bound in a nucleus [19], deformed nuclei [20], and super deformation [21] and to establish an effective nuclear shell model scheme [22][23][24], while the KG equation has been used in describing a wide variety of phenomena, which include classical wave systems, such as the displacement of a string attached to an elastic bed [25] and quantum system based on scalar field theories [26].
Recently, the Schrödinger equation is solved exactly for its bound states (energy spectrum and wave functions) [27][28][29] to study the spectral properties in a 2D charged particle (electron or hole) confined by a harmonic oscillator in the presence of external strong uniform magnetic field → along the direction and Aharonov-Bohm (AB) flux field created by a solenoid. The energy levels and the wave functions of an electron confined by 2D harmonic and pseudo-harmonic oscillators have been studied in presence of external fields [30,31] using the Nikiforov-Uvarov (NU) method. Overmore, it is natural to study the relativistic effects of the external magnetic and AB flux fields on the KG equation for pseudo-harmonic oscillator potential, especially for a strong coupling by means of the NU method [32]. The 2D solution of Schrödinger equation for the Kratzer potential with and without the presence of a constant magnetic field is studied within the framework of the asymptotic iteration method [33]. The energy eigenvalues are obtained analytically (numerically) for the absence (presence) of uniform magnetic field case. These results have been obtained by using different Larmor frequencies ( ̸ = 0) and potential parameters are compared with the results in the absence of magnetic field case ( = 0). Overmore, the spectral properties of the 2D Schrödinger equation for the pseudo-harmonic-Coulomb-linear potential are studied using the analytical iteration method [34]. Under spin symmetry, the energy states and wave functions of the Dirac equation for the Killingbeck (harmonic oscillator plus Cornell) potential have been carried out using the wave function ansatz method [35]. The spin and pseudospin symmetry in Dirac equation have been studied in the Killingbeck potential within the context of quasi-exact solution [36]. Further, the energy eigenvalues and normalized eigenfunctions of the radial Schrödinger equation in -dimensional Hilbert space for the quark-antiquark interaction Killingbeck potential have been obtained using the power series technique via a suitable choice of ansatz to the wave function [37,38].
Very recently, we studied the scalar charged particle exposed to relativistic scalar-vector Killingbeck potentials in presence of magnetic and Aharonov-Bohm flux fields and obtained its energy eigenvalues and wave functions using the analytical exact iteration method [39]. Therefore, the behavior of a spinless relativistic particle moving under the Killingbeck potential in a static magnetic and AB flux fields has not been investigated yet, and in this work, we aim to solve the KG equation in 2D for equal mixture of scalar and vector Killingbeck potentials with and without constant magnetic and AB flux fields for the first time. We present the exact energy eigenvalues and wave functions of the Killingbeck potential for any and quantum numbers in a constant magnetic field with the different Larmor frequencies using the wave function ansatz method [36]. The non-relativistic energy eigenvalues and wave functions of our solution are presented by making an appropriate mapping of parameters. Overmore, special cases of KG for equal scalar-vector Killingbeck potentials are also presented in the presence ( ̸ = 0, ̸ = 0) and absence ( = 0, = 0) of uniform fields.
The structure of this paper is as follows. We study effect of external uniform magnetic and AB flux fields on a relativistic spinless particle (antiparticle) under equal mixture of scalar and vector Killingbeck potentials in Section 2. We discuss some special cases in Section 3. Finally, we give our concluding remarks in Section 4.

The Klein-Gordon Atom
The Klein-Gordon atom for the spinless particle with mass and charge − moving in external electromagnetic field and AB flux field given by potentials ( ), ( ) and → reads [40,41] The scalar and vector potentials are chosen in the following form [27][28][29][30]: where ( ) is the Killingbeck potential, that is, harmonic oscillator potential plus Cornell potential [35][36][37][38], which is extensively used in particle physics [42,43]. Moreover, the vector potential in the symmetric gauge is defined by where the applied magnetic field → = (0, 0, ) is perpendicular to the plane of transversal motion of the particle and → 2 describes the additional AB flux field Φ AB created by a solenoid in cylindrical coordinates [32]. The wave function in (1) is defined by where is the eigenvalue of angular momentum. The relationship between the attractive scalar and repulsive vector potentials is given by ( ) = ( ), where −1 ≤ ≤ 1 is arbitrary constant and hence the KG equation could be reduced to a Schrödinger-type second-order differential equation as follows: Now, we will treat the bound-state solutions of the two cases in (4) as follows.
Advances in High Energy Physics 3 2.1. The Positive Energy Solution. The positive energy states require that ( ) = ( ) (i.e., = 1 case) which, in the non-relativistic limit, corresponds to the solution of the wave equation: where ( , ) stands for non-relativistic wave function. Thus, the choice ( ) = ( ) produces a non-relativistic limit with a potential function 2 ( ) and not ( ). Accordingly, it would be natural to scale the potential term in (4) and (5) so that in the non-relativistic limit the interaction potential becomes ( ) not 2 ( ). Thus, we need to recast (4) for the ( ) = ( ) as [32,39,40,44] where and in order to simplify (6), we introduce new parameters 1 = + 2 and 2 = − 2 so that it can be reduced to the following form: By inserting (2) and (3) into (8), we obtain where the effective potentials depending on the magnitudes of two fields strength with and are the Larmor frequency and a new eigenvalue of angular momentum (magnetic quantum number), respectively. It is worthy to mention that the frequency Ω is called the cyclotron frequency [45]. This is the frequency of rotation corresponding to the classical motion of a charged particle in a uniform magnetic field and Ω/2 is the Larmor frequency in units of Hz [45]. Moreover, we take as integer with the flux quantum Φ 0 = ℎ / . Here, ( ) is a pure Killingbeck potential (2), the second term is the harmonic oscillator-type potential, and the other terms are the rotational potential creating the rotational energy levels. Equations (9a), (9b), and (9c) can be alternatively expressed as with the asymptotic behaviors (0) = 0 and (∞) → 0. It is interesting to look at the results obtained from (10a) for a special case 2 = 4 = 0; that is, we have which corresponds to the differential equation of a harmonic oscillator with a centrifugal term. It is a form similar to the one given by (7) in [46] with the equivalence 1 → − ( +1), 5 → 2 and 3 → − ( , , and are the parameters used in [46]). So by putting the parameter values into (14) of [45], we obtain the energy equation as Now, we find a solution to (10a) by making the following choice of the wave function [37,38,47,48]: where and are parameters whose values are to be determined in terms of the potential parameters , , and . Substituting (14) into (10a), we obtain the series [36] 4 Advances in High Energy Physics with = ( + ) ( + − 1) + 1 , where the two parameters and are determined in order for the radial wave functions ( ) to be finite everywhere and vanish at = 0 and as → ∞. To obtain the recurrence relation which can connect various expansion coefficients , we make identical powers of in (15), that is, equate the coefficients of + −2 to zero. Thus, the relations become where = 0, 1, 2, . . ., with 0 ̸ = 0. Here, the positive sign of parameter = + 1/2 has been selected. The power series for large values of or is convergent. For convenience, we take the ratio of two successive terms, that is, +1 / , which becomes It is apparent from the above relations that the power series converges to zero when → ∞. Hence, the series must be truncated (bounded) for = max . At this value of , we obtain the following equations: where the negative sign of the coefficient has been chosen in (19a). It is worth noting that the potential parameter must be positive when = 0. Hence, (19c) gives a restriction on the potential parameters , , and as follows: and it follows that (19d) gives the energy formula as where = 0, 1, 2, . . .. We may find a solution to the above transcendental equation as = (+) KG = ( , ). Overmore, the wave function (14) with the help of (17) where , is the normalization constant. In the non-relativistic limit, when + 2 → 2 , + 2 → and = 1, (21) and (22a) give the energy formula and wave function respectively, for the Schrödinger-Killingbeck system with the following restriction on the potential parameters , , and given by Notice that the present model has been solved in 2D space with an external uniform magnetic field since it is perpendicular to the plane, where the vector and scalare Cornell potentials have the dimensions of = and ; that is, = = ( , ) [32]. However, without the magnetic field, the model can be solved in any desired dimensional space by considering the change → + ( − 2)/2 [49,50], where refers to spatial dimension and also the existence of interdimensional degeneracy.
As shown in Figure 1(a) and (9b), the effective potential function changes in shape when the magnetic field strength increases, say = 8 and in absence of AB flux field. The energy levels are raised when the strength of the magnetic field increases and in absence of AB flux field = 0. We see that the effective potential changes gradually from the pure pseudo-harmonic oscillator potential, when = 0, to a pure harmonic oscillator type behavior in short potential range when = 8. In Figure 1(b), the effective potential (9b) which is pseudo-harmonic oscillator when = 0 becomes sensitive to the increasing AB flux field = 8 in the short range region; that is, 0 < < 4 a.u.

Discussions
In this section, we discuss some special cases of interest from our general solution.
(i) If we set = = 0, the effective Killingbeck potential turns to effective harmonic oscillator potential in relativistic case as eff ( , , ) = 2 + 2 2 2 The bound state solutions (with the change = 2( + 1)) are [32,36] 2 − 2 4 = 2 2 ℎ + 4ℎ 2 ( + + 1 2 ) Advances in High Energy Physics 7 where , is the normalization constant. On the other hand, the nonrelativistic effective harmonic oscillator potential eff ( , , ) = 2 + has bound states as [32] ( , ) = ℎ (21), the killingbeck potential turns to the Coulomb potential case (i.e., ( ) = − / and = + 1). Thus, we obtain the relativistic energy formula and the wave function Overmore, in the nonrelativistic limit, we have = ℎ + ℎ ( + ) − 2 4 2ℎ 2 ( + + 1/2) 2 , (39) and the wave function becomes Notice that in the absence of external fields = = 0, the problem can be solved in any dimension. Thus, applying the transformation | | → + 1/2, our previous results are reduced to the well-known three-dimensional bound state solutions for the harmonic oscillator and the hydrogen-like atoms in the Coulomb fields [40]: respectively. According to (41), the energy spectrum can be found for the scalar particle and antiparticle. In Figure 2(a), we plot the effective potential for the case of low vibrational ( = 0, 1, 2, 3) and rotational ( = 1) levels for various Larmor frequencies = 0, 1, 5, 8 and = 0 case. As shown in Figure 2(a) and (31), the effective potential function changes in shape as well as the bound state energy eigenvalues increase when = 8. It is shown that the energy levels are raised when the strength of the magnetic field increases and in the absence of AB flux field. It is also obvious that the effective potential changes gradually from the pure pseudo-harmonic oscillator potential, which is a nonmagnetic ( = 0) and AB flux ( = 0) fields case, to a pure harmonic oscillator type behavior in short potential range when the strength of the applied magnetic field is increased to = 8. If we consider a strong magnetic field case = 8 which has the shape of pure harmonic oscillator potential function, the energy difference between adjacent energy levels is nearly equal which is a known characteristic of the pure harmonic oscillator potential. In Figure 2(b), we plot the effective potential (34) in the absence of the magnetic field and in the presence of AB flux field in the short range region.
In Tables 1 and 2, we show the effect of magnetic field and AB flux field, respectively, on the low vibrational 8 Advances in High Energy Physics and rotational relativistic energy states of the harmonic oscillator potential. As shown in Table 1, when the magnetic field is not applied and without AB flux field ( = 0, = 0), the spacing between the energy levels of the effective potential is narrow and decreases with increasing . But when the magnetic field strength increases, the energy levels of the effective potential increase and the spacings between states also increase. In Table 2, when the AB flux field is applied and without magnetic field, the energy states become degenerate for various values of and and for various AB flux field strength values. In Tables 3 and 4, we show the effect of magnetic field and AB flux field, respectively, on the low vibrational and rotational nonrelativistic energy states of the harmonic oscillator potential. As shown in Table 3, when the magnetic field is not applied and without AB flux field ( = 0, = 0), the energy states are equally spaced (the pure harmonic oscillator case). But when the magnetic field strength is raised, the energy levels of the effective potential increase and the spacings between states also increase. In Table 4, when the AB flux field is applied and without magnetic field, the energy states become degenerate and equally spaced for various values of and and for various AB flux field strength values.

Concluding Remarks
To sum up, in this paper, we have studied the solution of two-dimensional KG and Schrödinger equations with the Killingbeck potential for low vibrational and rotational energy levels without and with a constant magnetic field having arbitrary Larmor frequeny and AB flux field. We have applied the wave function ansatz method for ̸ = 0 (with magnetic field) and ̸ = 0 (with AB flux field) to obtain analytical expressions, in closed form, for bound state energies and wave functions of the spinless relativistic particle subject to a Killingbeck interaction expressed in terms of external uniform magnetic and AB flux fields in any vibrational and rotational states. The above results show that the problems of relativistic quantum mechanics can be also solved exactly as in the non-relativistic ones. We considered the solution of both positive (particle) and negative (antiparticle) KG energy states. It is noticed that the solution with equal mixture of scalar-vector potentials can be easily reduced into the well-known Schrödinger solution for a particle with an interaction potential field and a free field, respectively. We have also studied the bound-state solutions for some special cases including the non-relativistic limits (Schrödinger equation for harmonic oscillator and Coulomb potentials under external magnetic and AB flux fields) and also the KG equation for harmonic oscillator and Coulomb interactions. The results show that the splitting is not constant and is mainly dependent on the strength of the external magnetic field and AB flux field. In order to show the effect of constant magnetic and AB flux fields on the vibrational and rotational energy levels of the harmonic oscillator, we plot the effective potential and corresponding energy levels with increasing Larmor frequency and flux field for special potential parameters. We have seen that the effective potential function and corresponding energy levels are raised in energy when magnetic and AB flux field strengths increase. The effective potential function behavior gradually changes from the pure pseudo-harmonic oscillator to a pure harmonic oscillator shape in short potential range as the magnetic and AB flux fields strengths increase.