Analytical solutions for the Klein–Gordon equation with combined exponential type and ring-shaped potentials

In this study, we have successfully obtained the analytical solutions for the Klein–Gordon equation with new proposed a non-central exponential potential \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$V\left( r \right) = D\left[ {1 - \sigma_{0} \coth (\alpha r)} \right]^{2} + (\eta_{1} + \eta_{2} \cos \theta )/r^{2} \sin^{2} \theta$$\end{document}Vr=D1-σ0coth(αr)2+(η1+η2cosθ)/r2sin2θ. Our approach involves a proper approximation of the centrifugal term, with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${l}{\prime}$$\end{document}l′ representing the generalized orbital angular momentum quantum number, and the utilization of the Nikiforov–Uvarov method. The resulting radial and angular wave functions are expressed in terms of Jacobi polynomials, and the corresponding energy equation is also derived. Our calculations of the eigenvalues for arbitrary quantum numbers demonstrated significant sensitivity to potential parameters and quantum numbers. Additionally, we evaluate the dependence of energy eigenvalues on screening parameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}α for arbitrary quantum numbers \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n_{r}$$\end{document}nr and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N$$\end{document}N to establish the accuracy of our findings. Furthermore, we determine the non-relativistic limits of the radial wave function and energy equation, which align with corresponding previous results in the case where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta_{1} = \eta_{2} = 0$$\end{document}η1=η2=0.

In this study, we have successfully obtained the analytical solutions for the Klein-Gordon equation with new proposed a non-central exponential potential V(r) = D 1 − σ 0 coth(αr) 2 + (η 1 + η 2 cos θ)/r 2 sin 2 θ .Our approach involves a proper approximation of the centrifugal term, with l′ representing the generalized orbital angular momentum quantum number, and the utilization of the Nikiforov-Uvarov method.The resulting radial and angular wave functions are expressed in terms of Jacobi polynomials, and the corresponding energy equation is also derived.Our calculations of the eigenvalues for arbitrary quantum numbers demonstrated significant sensitivity to potential parameters and quantum numbers.Additionally, we evaluate the dependence of energy eigenvalues on screening parameter α for arbitrary quantum numbers n r and N to establish the accuracy of our findings.Furthermore, we determine the non-relativistic limits of the radial wave function and energy equation, which align with corresponding previous results in the case where η 1 = η 2 = 0.
One of the main tasks of relativistic and nonrelativistic quantum mechanics is the description of the motion of microparticles (molecules, atoms, atomic nuclei, etc.) in an external potential field [1][2][3][4][5] .In this case, as is known, depending on the type of potential energy of particle interaction with the field V (r) , the equation of motion can have solutions that belong to both discrete and continuous energy spectra.The solutions corresponding to the discrete spectrum describe the bound states of the particle, the states of the continuous spectrum describe the scattering of the particle in the potential field.Bound states correspond to the finite motion of a particle, when the particle is held by the forces of attraction of the field in a limited region of space.In other words, bound states arise when the interaction potential has the form of a well.If the height of the potential well is V 0 , then at particle energies E < V 0 the energy spectrum is discrete, and at energies E > V 0 it is continuous.Continuous states correspond to the infinite motion of a particle, when the particle comes from infinity and, scattering, goes to infinity.We emphasize that when a particle is in a potential field, in addition to bound states and scattering states, there can also be unstable states that decay with time.Indeed, if the energy of the particle is greater than the value of the potential at infinity V ∞ , i.e. belongs to the region V ∞ < E < V 0 , then sooner or later, due to the tunnel effect, the particle will fly out of the region of attraction of potential forces.
The basic equation of motion of nonrelativistic quantum mechanics is the Schrödinger equation.Based on this equation, nonrelativistic quantum mechanics explains the main properties of molecular, atomic and nuclear physicochemical phenomena 1,2,5 .Nonrelativistic quantum mechanics underlies the theory of solids.In all these cases, different phenomenological interaction potentials are used.The synthesis of non-relativistic quantum mechanics with the theory of relativity led to the creation of relativistic quantum mechanics and quantum field theory.If the energy of a particle is comparable to its rest energy, but the interaction potential is still insufficient for the formation of particle-antiparticle pairs, then the motion of such a particle must be described by relativistic quantum mechanics.The basic equations of motion of relativistic quantum mechanics are the Klein-Gordon (KG) and Dirac equations [3][4][5] .In this regard, we also note that there is another version, namely the finite-difference version of relativistic quantum mechanics, which was developed and applied to solve a number of problems in particle physics in papers (see 6 and the references therein).An analogue of the Schrödinger equation in finite-difference relativistic quantum mechanics is a Schrödinger-type finite-difference equation in configuration space.Finitedifference relativistic quantum mechanics has many important features of non-relativistic quantum mechanics.In particular, the wave functions here admit a probabilistic interpretation.
In the framework of relativistic and nonrelativistic quantum mechanics, potential models have always played an important role and will play, perhaps, even the main role in the physics of the nucleus and elementary particles, in the physics of atoms and molecules.The importance of a potential model is primarily determined by how well it describes certain properties of the physical system under study.The physical properties of micro-objects in potential models are described and interpreted using various wave equations, such as the Schrödinger equation, the Dirac equation, the KG equation and the relativistic finite difference equation.In this case, exact or phenomenologically introduced central and non-central interaction potentials are widely used.There are many such potentials.The most well-known potentials widely used in the relativistic and nonrelativistic energy regions are the Coulomb potential and the harmonic oscillator potential, as well as their various combinations.][9][10][11][12][13][14][15][16][17] .
Non-central potentials are good models in quantum chemistry, in nuclear and atomic physics, and in molecular physics.Ring potentials can be used in quantum chemistry to describe organic ring-shaped molecules such as benzene, and in nuclear physics to study the interactions of a deformed pair of nuclei and spin-orbit coupling when particles move in potential fields.This potential is also used as a mathematical model for describing the vibrations of diatomic molecules and is a convenient model in other physical situations.For example, the Hartmann potential 15 is one of the useful combinations in the form of Eq. ( 1), which is used to describe organic molecules in the framework of nonrelativistic quantum mechanics.Other important off-center potentials were proposed by Hautot 16 to describe the motion of a charged non-relativistic particle in an electric field.The Hartmann potential is a special case of one of the Hautot potentials.In 17 , Quesne investigated a new non-central potential, obtained by replacing the Coulomb part of the Hartmann potential by a harmonic oscillator term.In 51 , the relativistic problem of the motion of a scalar particle in a non-central field, which is equal to the sum of the Coulomb and one of the Hautot potentials, was solved.Here also a dynamic symmetry group for the radial wave functions is constructed.
On the other hand, within the framework of nonrelativistic quantum mechanics, the authors of 52 found the energy eigenvalues for a hyperbolic potential of the form 53 (also called the exponential type potential) where D , σ 0 andα are positive constant parameters.It was shown in [52][53][54] that this potential is closely related to the Morse, Kratzer, Coulomb potentials, the harmonic oscillator, and other interaction potentials.The hyperbolic potential is an important exponential type potential 53 .For more information about its properties and possible applications, the reader can refer to [52][53][54] .
As well known that an analytical and exact solution, for example, of the KG equation is possible only in the s-wave case with an orbital quantum number l = 0 for some potentials.In this regard, we emphasize that, as is known, one of the frequently used approximate methods for solving equations for any l = 0 value of orbital quantum number of motion with a central or noncentral potential are the Greene and Aldrich and Pekeris approximations 55,56 .These approximations widely applied to the solution of the KG equation for different potentials, such as Hulthén potential 57,58 , generalized Hulthén potential [59][60][61] , Manning-Rosen potential 62,63 , Wood-Saxon potential 64 , generalized hyperbolic potential 65 , Deng-Fan molecular potentials 66 , inversely quadratic Hellman potential 67 and Kratzer Potential 68 and similarly for the case of combined potentials like Hulthén plus Yukawa potential [69][70][71] , Manning-Rosen plus a class of Yukawa 72 , Hellmann plus modified Kratzer potential 73 , Mobius squared plus Eckart potential 74 , and Eckart plus a ring-shaped like potential studied in Refs. 75.
The purpose of this work is to analytically solve the KG equation for the new proposed as sum of hyperbolic potential and a ring-shaped potential of the form where η 1 and η 2 are positive constant parameters that characterize the properties of interaction potentials.The combined potential considered in this paper is obtained as a linear combination of an exponential potential with a ring-shaped potential in the form: (1) Here we generalize the results of 52 : firstly, to the relativistic case, and secondly, to the case of noncentral potential (4).
Let us note that some authors [41][42][43][45][46][47][48] have assumed that in the KG equation the scalar potential is equal to the vector potential and have solved exactly or approximately the KG equation. In thi paper, we also assume that the scalar and vector potentials are equal to each other, and we solve the KG equation, and present the approximately analytical expressions for the wave functions and for the discrete energy spectrum.
Thus, new potential (4) as self-behavior is a non-central potential.
In this study, we try to study analytical solutions of the KG equation for any value of the orbital quantum numbers l = 0 in the framework of the Greene and Aldrich approximation 55 for the new proposed combined potential (4) using the Nikiforov-Uvarov (NU) method 76 .This method is based on solving a second-order linear differential equation by reducing it to a generalized hypergeometric type equation, which is a homogeneous second-order differential equation with polynomial coefficients of degree not exceeding the corresponding order of differentiation.Many studies show the strength and simplicity of the NU method in solving central and noncentral potentials for arbitrary l = 0 states.
We emphasize that in 77 the relativistic problem of bound states for the hyperbolic potential ( 2) was also considered.However, there are two differences from our work: firstly, we consider the case of a non-central potential, and secondly, we use a different approximation for the centrifugal term.
This paper is organized as follows.In Section "Brief summary of the Nikiforov-Uvarov method" given a brief summary about the Nikiforov-Uvarov method.In Section "Klein-Gordon equation for non-central potential" separates the variables in the KG equation.In Section "Radial Klein-Gordon equation" we find analytically the bound-states solutions of the radial KG equation with the non-central potential (4).In Section "Angular part of the Klein-Gordon equation", we present the solution of the angular part of the KG equation.We discuss the nonrelativistic limits of the radial wave function and the energy equation and some of their particular cases in Section "Limiting case".We present numerical estimates for the energy eigenvalues and their corresponding normalized wave functions in Section "Discussion of the numerical results".Finally, we conclude our work in Section "Conclusion".

Brief summary of the Nikiforov-Uvarov method
The NU method solves a second-order differential equation of the following form 76 : Here σ (g), and σ (g) are polynomials of degree at most 2, a τ (g) is a is a polynomial of degree at most 1.Then, using the following factorization for the function χ(s): we reduce Eq. ( 5) to a hypergeometric type equation of the form Function ϕ(g) must satisfy the logarithmic derivative with where the primes denote the derivative with respect to g , and it can be at most first order.The term inside the square root is transposed as the zero discriminant of the second-order polynomial.Therefore, the equation reduces to an equation of hypergeometric type, where one of its solutions is y g .Therefore, the expression for k is found after solving such an equation by the NU method.
Therefore, the equation reduces to an equation of hypergeometric type, where one of its solutions is y g .Therefore, the polynomial expression σ (g) = σ (g) + π 2 (g) + π(g) τ (g) − σ ′ (g) + π ′ (g)σ (g) can be divided by a factor σ (g) , such that σ (g)/σ (g) = .Here we use the following relations where τ (g) satisfy the condition τ ′ g < 0.Here k is parameter, the definition of which is essential in the calcu- lation π g .This parameter is simply determined from expression (9) by setting the square root discriminant (4) to zero.Hence the general quadratic equation for k can be obtained.Then k can be used to calculate the energy eigenvalues using the formula Polynomial solutions y g are determined by the Rodrigues formula: where C n normalization constant, and ρ(g) denotes a weight function that obeys [81], Usually this is called Pearson's differential equation.

Klein-Gordon equation for non-central potential
We know well that in relativistic quantum mechanics the dynamics of spin-zero particles is described by the KG equation.For scalar V 1 (r) and vector V 2 (r) potentials, the KG equation in atomic units ( = c = 1) in a spherical coordinate system has the form 3 where M and E are the mass and energy of the relativistic particle.For simplicity, we assume that V 1 (r) and V 2 (r) are equal: Replacement (16) in Eq. ( 15) gives where, as we know, Taking into account (3) and (4), Eq. ( 17) can be written as This equation in spherical coordinates admits separation of variables.Therefore, the wave function ψ(r) can be written in the factorized form: It should be immediately noted that, since the operator L = −i∂ φ commutes with the Hamiltonian of Eq. ( 19), the function �(φ) has the standard form After substituting function (20) into Eq.( 19), the following system of second-order differential equations is obtained: where µ and m are separation constants.Let's define the generalized orbital quantum number l′ , which in the general case is not an integer: µ = l′(l′ + 1) .It dependents on the energy, т.e.l′ = l′(E) .In the nonrelativistic limit, this dependence is vanish.Obviously, in the absence of a ring-shaped potential, the number l′ will coincide with the orbital quantum number l =0, 1, 2, …, i.e. at η 1 = η 2 = 0 we have where P m l (cosθ ) is Legendre function.From Eq. ( 22) we can extract the effective potential in the following form: www.nature.com/scientificreports/In Fig. 1 shows the dependences of potential (4) V (r, θ) both on the separation distance r and on the angle θ .Also in Fig. 2, we visualized the dependence of the potential as a function of the separation distance for some fixed values of θ.

Radial Klein-Gordon equation
In this section, applying the NU method, we find an analytical solution of the radial KG Eq. ( 12).The NU method is one of the effective methods for solving a certain type of second-order differential equations and is widely used to explicitly obtain the polynomial parts of the wave functions and discrete bound-state energy levels for quantum systems for some solvable interaction potentials.This has been shown by many studies, as noted above.
Let us rewrite the differential Eq. ( 22) in expanded form: As can be seen, the effective potential contains the centrifugal term.Therefore, the KG equation with this potential cannot be solved analytically exactly for µ = 0 .To solve analytically Eq. ( 25) for the case µ = 0 , we apply a special approximation to the centrifugal potential.If αr ≪ 1 we use the following scheme for approach- ing the centrifugal potential 56 :  www.nature.com/scientificreports/After applying approximation (27) in Eq. ( 26), we get We will find an analytical solution of the radial Eq. ( 28) using the NU method.
To apply the NU method to the solution of Eq. ( 28), we introduce a new variable in the form g = e −2αr ∈ [0, 1] , and r ∈ [0, ∞) .Then we can represent Eq. ( 28) in the following form: or We are interested in bound states: E > M , ε > 0 .As can be seen, Eq. ( 30) is a suitable type for the NU method.Now we can successfully apply this method to determine the energy spectrum.Comparing Eq. ( 30) with Eq. ( 5) for polynomials σ (g) , τ (g) and σ (g) , we obtain: where By inserting the factorization (6) in Eq. (30), we obtain a hypergeometric type equation similar to Eq. (7).Given condition (8) for ϕ(g) , the function π(g) can be determined from formula (9): To find the parameter k , we use the condition that the discriminant of the radical square trinomial in (32) is equal to zero.From here we get After substituting expression (34) in Eq. (33), we obtain the following expressions for π(g) Of the four possible expressions for π(g) , we choose the case when the function τ (g) has a negative deriva- tive, since other forms are physically unacceptable.Thus, we choose the appropriate function for π(g) and τ (g) as follows: 10), the parameter can be written as where Vol.:(0123456789) www.nature.com/scientificreports/Note that according to the NU method for a non-negative integer n r = 0, 1, 2, 3, . . ., the hypergeometric type Eq. ( 30) has a unique polynomial y = y n r (g) solution of degree n r .On the other hand, according to Eq. ( 12), we have In expressions (38) and (39), the left parts are equal, therefore, by equating the right parts of these expressions and solving the equation for √ a 3 , we find an analytical expression for the energy eigenvalues in the form: where p = γ D α 2 , q = 1 4 + µ + pσ 2 0 , the parameter µ is determined by the equation of the angular part of the KG equation.Now we are trying to obtain the corresponding radial wave function χ(g) for the considered potential.σ (g) and π(g) from expressions (31) and ( 36) substituting into Eq.( 8) and then solving the first order differential equation, we find the function ϕ(g) .She becomes equal The polynomial part of the wave function χ(g) is a polynomial y n r (g) of hypergeometric type, which can be obtained from the Rodrigues relation (13).Solving Eq. ( 14) for the weight function ρ(g) , we find Now substituting ρ(g) and σ (g) into Eq.( 13), we get Let us now compare formula (43) with the definition of the Jacobi polynomials 78 It follows from the results of this comparison that Substitution ϕ n r (g) from Eq. ( 41) and y n r (g) from Eq. (45) into Eq.( 6) for the radial wave function gives the expression The Jacobi polynomials are expressed in terms of the hypergeometric function as follows 78 If to substitute expression (47) in (46) it is possible to express χ n r (g) radial wave function through hypergeo- metrical polynoms as: where K = 1 2 + q = 1 2 + 1 4 + µ + pσ 2 0 .The normalization constant C n r in Eq. ( 48) is found from the condition of orthonormality of the radial wave function χ n r (g) , as To calculate the integral, we use the following master formula 78 : Vol:.( 1234567890) Here we take the potential parameters as σ 0 = 0.20 , D = 2 , α = 0.30 , η 1 = 2 and η 2 = 1.

Angular part of the Klein-Gordon equation
The eigenvalues and eigenvectors of the azimuthal part of the KG Eq. ( 23) can also be obtained by the NU method.To do this, we introduce a new variable z = cos θ .Its substitution into Eq. ( 23 gives This equation is an equation of type (5) in which Therefore, using the NU method, for the function π(z) depending on the angle θ , we find the following expression: The constant parameter k takes the form: From the condition of applicability of the NU method, it is known that the parameter τ (z) must satisfy the condition τ ′ (z) < 0 .Therefore, in this case τ (z) takes the form By introducing the notation ξ = m 2 +γ η 1 +ϑ 2 , rewrite ( 56) and (57)
Similarly to the radial wave function, we now find the explicit form of the wave function depending on the polar angle.First factorize the function �(z): Using Eq. ( 8), for ϕ(z) obtain where On the other hand, to find an explicit expression for polynomials y N (z), N = 0, 1, 2, . . .we first need to find the weight function ρ(z) .From the Pearson Eq. ( 14) for the weight function, we get: Accounting for (64) in formula ( 13) makes it possible to obtain the Rodrigues formula for the polynomial y N (z) 76,78 : Based on the definition (44) of Jacobi polynomials, we can write: Substituting expressions ( 64), ( 67) and ( 68) into Eq.( 63), we obtain for �(z) the following expression: where CN -normalization constant determined from the normalization condition.Using the orthogonality rela- tion for Jacobi polynomials 78 : after simple calculations of the normalization constant, we obtain: Thus, after substituting Eqs. ( 21), ( 44) and (67) into Eq.( 20), we directly obtain the total wave function for the considered quantum system with a combined potential defined as the sum of potentials of exponential and ring types Wave function (72) satisfies the normalization condition with C n r N = C n r • CN .

Limiting case
In this section, we find the nonrelativistic limits of the radial wave function (46) and energy Eq.(62) of the system under consideration and compare them with the corresponding results in 52 .In doing so, we take into account that the Eq. ( 17) reduces to the Schrödinger equation for the potential 2 V (r) .First, we restore the constants and c in formulas ( 46) and (62), i.e. we will make replacements in them ( 63) www.nature.com/scientificreports/and give the nonrelativistic limits of the parameters a 1 , a 2 , a 3 and the energy E: Here we use the following notation (compare with 52 ): where µ N obtained from µ by replacing γ with γ N = 4M/ 2 .Now it is easy to calculate the nonrelativistic limit of function (46), which has the form lim c→∞ = χ n r = χ n r N , where The calculation of the nonrelativistic limit of the energy equation is also not difficult, and as a result we get Thus, expressions ( 77) and ( 78) are nonrelativistic wave functions and energy spectrum of the non-central potential V (r) = D[1 − σ 0 coth(αr)] 2 + (η 1 + η 2 cos θ)/r 2 sin 2 θ for the arbitrary l′-wave, where l ′ = N + ξ.
It is easy to check that for η 1 = η 2 = 0 formulas ( 77) and ( 78) coincide with the corresponding formulas in 53 with normalization constant:

Discussion of the numerical results
In this paper, energy eigenvalues and their corresponding eigenfunctions are obtained for arbitrary states l = 0 by analytically solving the modified Klein-Gordon equation for the sum of the exponential and ring-shaped potential applying the Nikiforov-Uvarov method.
In this section, we carry out numerical evaluations of our analytical results.We carry out the numerical calculations using the MATLAB R2023a and MATHEMATICA 8 package programs.The figures in this work have been drawn with help of the program ORIGIN 8.5.1.
In Tables 1 and 2, we present several energy levels E n r N the dependence of the screening parameter for few values of the radial n r and N quantum numbers.In numerical estimation the energy eigenvalues of the screening parameter α we used the parameters of the potential as η 1 = η 2 = 1 , D = 4 and σ 0 = 0.2 , M = 1 , m = 0 .These parameters used in the numerical calculations are taken arbitrarily for illustrative purposes.In Figs. 3 and 4 are visualized the dependence energy spectrum of the screening parameter α for different radial n r and N = 0, 1, 2, 3 quantum numbers.Here, n r is the number of nodes of radial wave functions, and the usual principal quantum number is given by n = n r + l ′ + 1 .The screening parameter is changed in the interval α ∈ [0.01 ÷ 0.30] .It is also from the Figs. 3 and 4 seen that the energy spectrum are monotonically increasing with an increase in the screening parameter.Such an increase is almost logarithmic with increasing screening parameter α .Moreover, the energy eigenvalues increase with the increment as n r and N .Should be noted that states with the same total value n r + N have close energy values to each other.For example, the following states ( n r ,N) = (0, 1), (1, 0), ( n r ,N) = (0, 2), (1, 1), (2, 0).In addition, it was shown that the energy eigenvalues and the corresponding wave eigenfunctions are sensitive to the choice of radial n r , orbital l and N quantum numbers.Therefore, the study of the analytical solution of the modified KG equation for a linear combination of exponential potential and ring-shaped potential in the framework of quantum mechanics can contribute to obtaining valuable information about the dynamics in nuclear, atomic and molecular physics and open a new window for deeper research.Must be noted that central potential describe quantum system in the radial direction (region), but noncentral potential is describe quantum system in the azimuthal region.Thus, the combined (noncentral) potential keep a more information regarding the central potential, and give us for getting the better results about quantum mechanical system.

Conclusion
In this study, we presented solutions of the KG equation in the bound state with a new combined exponential type potential plus ring-shaped potential, using an approximation scheme to work with the centrifugal term.For an analytical solution in the framework of ordinary quantum mechanics, the well-established method is used.Equations for the numerical calculation of the energy eigenvalues and the corresponding normalized wave functions of the specified quantum system for any n r and N are obtained analytically, and they are presented in a closed and compact form.The normalized wave functions are expressed in terms of the hypergeometric function and of the Jacobi polynomial.It is obvious that the solutions in the bound state are more stable for the potential we are considering than individual cases.We have estimated the change in energy eigenvalues for various parameters and quantum numbers.The results obtained for the energy spectrum show that the energy eigenvalues are sensitive to the screening parameter α as well as the quantum numbers n r and N.
Therefore, obtaining results can provide important information about the dynamics in atomic and molecular physics and makes it possible to study this problem in more depth.Thus, we can conclude that the results obtained by us will be of interest not only to a theoretical physicist, but also to an experimental physicist, due to accurate and more general results.Table 1.Bound state energy eigenvalues E n r N as a function screening parameter α for various values of n r and N.For potential parameters used η 1 = η 2 = 1 , D = 4 and σ 0 = 0.2 , M = 1 , m = 0 , ( n r ; N) = (0, 1; 0, 1, 2, 3).
As a result, we obtain the following expression for C n r Reports | (2024) 14:5527 | https://doi.org/10.1038/s41598-024-53650-8www.nature.com/scientificreports/ 1) .Finally, the results obtained are of interest not only to a theoretical physicist, but also to an experimental physicist, since they are accurate and more general.