Abstract

The relativistic wave equations determine the dynamics of quantum fields in the context of quantum field theory. One of the conventional tools for dealing with the relativistic bound state problem is the Klein-Fock-Gordon equation. In this work, using a developed scheme, we present how to surmount the centrifugal part and solve the modified Klein-Fock-Gordon equation for the linear combination of Hulthén and Yukawa potentials. In particular, we show that the relativistic energy eigenvalues and corresponding radial wave functions are obtained from supersymmetric quantum mechanics by applying the shape invariance concept. Here, both scalar potential conditions, which are whether equal and nonequal to vector potential, are considered in the calculation. The energy levels and corresponding normalized eigenfunctions are represented as a recursion relation regarding the Jacobi polynomials for arbitrary states. Beyond that, a closed form of the normalization constant of the wave functions is found. Furthermore, we state that the energy eigenvalues are quite sensitive with potential parameters for the quantum states. The nonrelativistic and relativistic results obtained within SUSY QM overlap entirely with the results obtained by ordinary quantum mechanics, and it displays that the mathematical implementation of SUSY quantum mechanics is quite perfect.

1. Introduction

The exactly solvable problems for quantum systems have long been a subject of intense study in many branches of quantum physics. The main aim of an analytical solution of wave equations for this attention is that the wave function contains all the requisite information for the full description of a quantum system [16]. In physics, especially the relativistic quantum mechanical applications to particle and nuclear physics, the relativistic wave equations predict particles’ reaction at high energies [57]. The analytical solution of the Klein-Fock-Gordon (KFG) equation with physical potentials plays a central role in relativistic quantum mechanics since this wave equation perfectly defines the spinless pseudoscalar pions and Higgs boson.

In principle, numerous methods were developed, and they are still successfully implemented in solving the nonrelativistic and relativistic wave equations with some familiar potentials. The Nikiforov-Uvarov method [8], factorization method [9], Laplace transform approach [10], and the path integral method [11] and shifted 1/N expansion approach [12, 13] are used for solving radial and azimuthal parts of the wave equations exactly or quasiexactly in for various potentials. Additionally, there are numerous interesting research works about the KFG equation with physical potentials by using different methods in the literature [1426]. Among them, as an example, in Ref. [22], the -wave KFG equation with the vector Hulthén type potential was treated by the standard method. As reported by Talukdar et al., the scattering state solutions of the -wave KFG equation with the vector and scalar Hulthén potentials were obtained for the irregular and regular boundary conditions [25]. Besides, the supersymmetry method (SUSY) was also proposed for solving the wave equations analytically [2730]. Nonetheless, Okon et al. reported analytical solutions of the Schrödinger equation for the Hulthén-Yukawa plus inversely quadratic potential. [31] In Ref. [3239], the scalar potential, which is nonequal and equal to the vector potential, was supposed to get the bound states of the KFG equation for some typical potential from the ordinary quantum mechanics. Furthermore, KFG equation with the ring-shaped potential was investigated by Dong et al. [36]. If the condition where the interaction potential is insufficient to form antiparticle-particle pairs is considered, the KFG and Dirac equations can be utilized for the investigation of zero- and 1/2-spin particles, respectively.

When a particle is in a strong field, the relativistic wave equations should be considered in the quantum system. In any case, it can be corrected quickly for nonrelativistic quantum mechanics. The Hulthén potential is one of the essential short-range potentials in physics, extensively used to describe the continuum and bound states of the interaction systems. It has been applied to several research areas such as nuclear and particle, atomic, chemical, and condensed matter physics, so analyzing relativistic effects for a particle under this potential could become significant, especially for strong coupling. The Hulthén potential is defined as where and are the atomic number and the screening parameter, respectively. They determine the range for the Hulthén potential [40]. The Yukawa potential was proposed in 1935 as an operative potential to describe the strong interactions between nucleons [41]. It takes the following form: where describes the strength of the interaction and , its range. Unfortunately, for arbitrary states (), the KFG equation cannot get an exact solution with these potentials due to the centrifugal term of potentials. The numerous research works reveal the SUSY QM method’s power and simplicity in solving wave equations of the central and noncentral potentials for arbitrary states [4249].

In principle, the radial function nature at the origin was investigated particularly for singular potentials by Khelashvili et al. [50, 51]. While the Laplace operator is portrayed in spherical coordinates, the radial wave equation’s exact derivation demonstrates the perspective of a delta function term. Thus, the delta function term of the Laplace operator yields an essential contribution to the energy level. Although the various research attempts have provided satisfactory bound state energies using Hulthén and Yukawa potentials separately [1, 5257], we first considered these potentials under the linear combination form. [58] It is also worth mentioning that this potential can be used in nuclear physics to investigate the interaction between the deformed pair of the nucleus and spin-orbit coupling for the particle motion in the potential fields. Another fascinating perspective of this potential can be used as a mathematical model in the description of vibrations on the hadronic system’s side, and it can constitute a convenient model for other physical situations. The investigation of the relativistic bound states in the arbitrary -wave KFG equation with the linear combination of Hulthén and Yukawa potentials is quite interesting, and it can provide the deeper and accurate appreciations of the physical properties of the wave functions and energies in the continuum and bound states of the interacting systems. Inspired by all developments and works, in this paper, we present the solution of the relativistic radial KFG equation for the linear combination of Hulthén and Yukawa potentials, defined as where and is the screening parameter.

To study the system, we use an improved scheme to overcome the centrifugal term and the SUSY quantum mechanics [59, 60]. Despite our previous research effort on this potential [58], the investigation of this potential still needs to be clarified in detail. Accordingly, the main goal is to solve the KFG equation for the linear combination of Hulthén and Yukawa potentials by considering two cases, i.e., the scalar potential which is equal and unequal to vector potential by using SUSY QM. Thereby, the energy eigenvalues and corresponding radial wave functions are found for any orbital angular momentum case. Then, we compare the obtained results with the results obtained by the NU method in ordinary quantum mechanics to present the legitimacy and feasibility of this SUSY QM method. The remainder of the paper is structured as follows. In Section 2, we introduce the analytical solution of the radial KFG equation for the linear combination of Hulthén and Yukawa potentials from SUSY quantum mechanics. Next, the analysis of the results is presented in Section 3. Finally, Section 4 contains the conclusions.

2. Bound State Solution of the Radial Klein-Fock-Gordon Equation

2.1. Implementation SUSY Quantum Mechanics

Two different types of potential can be introduced into the KFG equation, which contains two objects: (i) the four-vector linear momentum operator and (ii) the scalar rest mass. Hence, the first one is a vector potential , which is introduced via minimal coupling, and the second one is a scalar potential , which is introduced via scalar coupling [5]. At this moment, they allow one to introduce two types of potential coupling: the vector potential and the space-time scalar potential . The natural units () are set throughout this study. In the spherical coordinates systems, the KFG equation with vector potential and scalar potential has the form where is the relativistic energy and denotes the rest mass of the system’s scalar particle. For the separation of the angular and radial parts of the wave function, in the stationary KFG equation with the linear combination of Hulthén and Yukawa potentials, the wave function should utilize the following wave function: and substituting this into Equation (4), the radial KFG equation is defined in the following form:

As it is known that the KFG equation with this potential can be solved exactly using a suitable approximation scheme to surmount the centrifugal term. To solve Equation (6) for , we ought to approximate the centrifugal term of the Yukawa potential in this system. As a result of this, while , the improved approximation scheme [6165] must be used as

Next, the vector and scalar potential forms for the general Hulthén and Yukawa potentials can be considered in the following forms:

Then, Equation (6) becomes

Thereby, the effective potential of the Hulthén and Yukawa potential linear combination has the following form: where

For investigation in detail, the nonrelativistic limit of the formula must be studied for the energy level. When , Equation (4) reduces to a Schrödinger equation for the potential . Based on supersymmetric quantum mechanics, the eigenfunction of ground state in Equation (6) should be in the following form: where and are normalised constant and superpotential, respectively. The connection between the supersymmetric partner potentials and of the superpotential is as follows [27, 28]:

The particular solution of the Riccati equation (Equation (13)) must be in the following form: where and are unknown constants. Having inserted Equation (14) into Equation (13) and taking into account that , we obtain

After small simplification, it can be rewritten as

From comparison of compatible quantities in the left and right sides of the equation (Equation (16)), we find the following relations for and constants:

Considering extremity conditions for wave functions, we obtain and . Solving Equation (19) yields and considering from Equations (18) and (19), we find that

From Equation (17) and Equation (21), we find that

After inserting (22) into (11) for the definitions of the energy eigenvalue of ground state for general case , we obtain the following energy level equation:

When , the chosen superpotential -. Inserting Equation (14) into Equation (13), the supersymmetric partner potentials and can be found in the following forms:

By using the superpotential from Equation (14), we can find radial eigenfunction of ground state in the following form: where ; , , and ; . Two partner potentials and which differ from each other with additive constants and have the same functional form are called the invariant potentials [59, 60]. Hence, for the partner potentials and given with Equations (13) and (14), the invariant forms are defined as where the reminder is independent of . If we keep going this procedure and make the following substitution , the whole discrete level of Hamiltonian can be written as and we obtain the following form:

In the following, we obtain the energy level equation in ordinary quantum mechanics

As seen from Equation (30), it is in a perfect agreement with the result obtained in Equation (27) of Ref. [49]. If we consider Equation (11) into Equation (30) and do some simple algebraic derivation, we can obtain the energy eigenvalue equation in the simplest form with .

Based on the SUSY QM method and knowing the ground state eigenvalues and eigenfunctions , all energy eigenvalues and eigenfunctions can be easily obtained. Briefly, using the following equation:

can be easily obtained in terms of the ground state wave functions. The superpotential depends on two parameters and the first partner potential has like that parameter . Hence, Equation (32) will be in the following form:

We define a new variable and factoring out the ground state wavefunction

Substituting into Equation (33) and using the ground state wavefunction (Equation (25)), we get

Based on comparison, it with the recursion relation in Ref [66]: it is seen that is proportional to the Jacobi polynomial . Thus, the normalized eigenfunction for this potential is taken in the following form: where . The normalization constant can be found by using the normalization condition by utilizing the following integral formula in Ref. [66]: where and . After making simple calculations, we arrive at the following expression for the normalization constant:

3. Results and Discussion

In this section, we present the numerical evaluation for the bound state solutions of the -wave KFG equation with the vector and scalar form of the linear combination of Hulthén and Yukawa potentials. To study the property of the energy levels regarding potential parameters in some quantum states, (see Figure 1) we take , , , , and and , , , , and . The little difference (0.01) in model potential parameters , , , and is quite sufficient, in order to see the energy level of quantum states displaye completely different behavior. In Figures 1(a) and 1(b), the energy levels of quantum states first are decreasing until some of small values (0.075-0.1), then the energy levels increase in the . In Figures 1(c) and 1(d), the energy levels of quantum states have very little variation for an interval of , and it causes the degenerate of all quantum states, then the energy levels of quantum states continue to gradually increase with increments of . These behaviors are better recognized in higher quantum states.

Behind that of these results, we can investigate some special cases. (i)In case , and , namely, , we obtainwhere (ii)If , , we obtain the energy level equation for Hulthén potential casewhere and . This result is in good agreement with the expression obtained in Equation (A.1) of Ref [58] (iii)In case , , but , we obtain the energy level equation for Yukawa potential, which is defined as the following form:where , and . This result is in good agreement with the expression obtained in Equation (A.3) of Ref [58]. Furthermore, this result is also the same with the expression for the constant mass case obtained in Ref [67]. One can easily see this by setting =1 and in Equation (39) of Ref [67] (iv)Also, if , and , namely, where (v)If and or in Equation (31), the potential reduces to Coulomb potential, , and the corresponding energy spectrum is obtained asand this result is the same with Equation (A.2) of Ref [67] (vi)If we take (the -wave case), the centrifugal term in Equation (9) disappears because and the equation turns to the -wave KFG equation. By setting in Equation (31), its energy spectrum equation is the following form:(vii)If we take = and =, and based on the following transformations: , , , and , we obtain the energy level equation of Equation (31) for the nonelativistic case. Briefly, because of the following relation:the energy level equation of Equation (31) for the nonrelativistic case, can be written the following form: which is good agreement with the result in Equation (28) (If B and C are considered zero as a special case) of Ref [31]. Generally, it is obviously seen from Equation. (31) that the bound states show more stability in the case of the linear combination of Hulthén and Yukawa potentials than Hulthén and Yukawa potential cases. Furthermore, the energy eigenvalues of the quantum states are considerably sensitive with depending potential parameters

4. Conclusion

To conclude, we admit that the SUSY QM method was presented to solve the KFG equation for the linear combination of Hulthén and Yukawa potentials. Hence, the energy eigenvalues and corresponding eigenfunctions of a mentioned quantum system were analytically obtained for arbitrary angular momentum and radial quantum numbers. Next, a closed form of the normalization constant of the wave functions was also found. Beyond that, it was also shown that the energy eigenvalues are considerably sensitive respecting quantum states. Finally, the results obtained within SUSY QM are in excellent agreement with the results obtained by ordinary quantum mechanics, and it confirms that the mathematical application of SUSY quantum mechanics is ideal for similar systems.

It is worth mentioning that the main results of this paper are the explicit and closed form expressions for the energy eigenvalues and the normalized wave functions. The method presented in this study is systematic, and in many cases, one of the most definite works in this field. In particular, the linear combination of Hulthén and Yukawa potentials can be one of the essential exponential potentials, and it probably provides a promising avenue in many branches of physics, especially in hadronic and nuclear physics.

Appendix

A. SUSYQM Method

For in SUSYQM, it is possible to define two nilpotent operators, and . They satisfy the following anticommutation relations:

Here, is the supersymmetric Hamiltonian operator and conventionally and . The and are also known as the supercharges operators. Here, is bosonic operator and is its adjoint. In terms of these operators, the Hamiltonian can be defined as [27, 28]: where the is named as the Hamiltonian of supersymmetric partner. Note also that and operators commute with . If we have zero ground state energy for (i.e., ), we can always represent the Hamiltonian as a product of a linear differential operators pairs in a factorable form. Therefore, the ground state obeys the Schrödinger equation as follows: hence,

This result makes us possible to globally reconstruct the above potential from the information of its ground state wave function that contain zero nodes. Hence, factorizing of is quite easy by using the following ansatz [27, 28]: where

After that, the Riccati equation for can be written as

Solving for from this equation, we can express it in terms of by

We obtain this solution by noticing that when is satisfied, we have We then introduce the operator which is written by reversing the order of the components. After a bit simplification, we find that is nothing but the Hamiltonian for new potential .

We call as supersymmetric partner potentials. For example, when the ground state energy of is with eigenfunction , from Equation (A.5), we can always write where

The SUSY partner Hamiltonian is defined by [27, 28]. where

Using Equation (A.12), for and , the energy eigenvalues and eigenfunctions are obtained as

Here represents the energy eigenvalue, where and denote the energy level and the th Hamiltonian , respectively. Hence, it is clear that if has bound states with corresponding eigenvalues , as well as eigenfunctions defined in , then we can always generate a hierarchy of Hamiltonians, i.e., such that the has the same spectrum of eigenvalue as , apart from the fact that the first eigenvalues of are absent in [27, 28]: where

We also have such that, by knowing all the eigenfunctions and eigenvalues of we also obtain the corresponding eigenfunctions and energy eigenvalues of the Hamiltonians .

Data Availability

Data availability statement does not exist.

Conflicts of Interest

All authors have no conflicts of interest.