Advances in High Energy Physics
Volume 2016 (2016), Article ID 7910341, 12 pages
http://dx.doi.org/10.1155/2016/7910341

Research Article

Relativistic Energy Analysis of Five-Dimensional q-Deformed Radial Rosen-Morse Potential Combined with q-Deformed Trigonometric Scarf Noncentral Potential Using Asymptotic Iteration Method

Subur Pramono,¹ A. Suparmi,² and Cari Cari²

¹Physics Department, Graduate Program, Sebelas Maret University, Jl. Ir. Sutami 36A Kentingan, Surakarta 57126, Indonesia
²Physics Department, Faculty of Mathematics and Fundamental Science, Sebelas Maret University, Jl. Ir. Sutami 36A Kentingan, Surakarta 57126, Indonesia

Received 28 March 2016; Revised 20 July 2016; Accepted 31 July 2016

Academic Editor: Edward Sarkisyan-Grinbaum

Copyright © 2016 Subur Pramono et al. 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.

We study the exact solution of Dirac equation in the hyperspherical coordinate under influence of separable -deformed quantum potentials. The -deformed hyperbolic Rosen-Morse potential is perturbed by -deformed noncentral trigonometric Scarf potentials, where all of them can be solved by using Asymptotic Iteration Method (AIM). This work is limited to spin symmetry case. The relativistic energy equation and orbital quantum number equation have been obtained using Asymptotic Iteration Method. The upper radial wave function equations and angular wave function equations are also obtained by using this method. The relativistic energy levels are numerically calculated using Matlab, and the increase of radial quantum number causes the increase of bound state relativistic energy level in both dimensions and . The bound state relativistic energy level decreases with increasing of both deformation parameter and orbital quantum number .

1. Introduction

Dirac equation as relativistic wave equation was formulated by P. A. M Dirac in 1928; the exact solution of Dirac equation for some quantum potentials plays a fundamental role in relativistic quantum mechanics [1]. In order to investigate nuclear shell model, spin symmetry and pseudospin symmetry solutions of Dirac equations have been an important field of study in nuclear physics. The concept of spin symmetry and pseudospin symmetry limit with nuclear shell model has been used widely in explaining a number of phenomena in nuclear physics and related field [2]. In nuclear physics, spin symmetry and pseudospin symmetry concepts have been used to study the aspect of deformed and super deformation nuclei. The concept of spin symmetry has been applied to the level of meson and antinucleon [3]. Pseudospin symmetry has been observed in deformed nuclei and can be enhanced in heavy proton-rich nuclei [4].

Solutions of Dirac equation for some potentials under limit case of spin symmetry and pseudospin symmetry have been investigated intensively whether in three- [5, 6], two-, or one- [7–13] dimensional space some -dimensional spherical symmetric spacetimes [14, 15]. However, The -dimensional Dirac equation with ()-dimensional separable noncentral potential has not been investigated yet; therefore, it may be worthy to investigate Dirac equation in 5 dimensions with separable 4-dimensional noncentral potential in this study.

In recent years, some researchers have studied solution of Dirac equation with quantum potentials with different application and methods. These investigations include Eckart potential and trigonometric Manning-Rosen potential using Asymptotic Iteration Method (AIM) [5], -deformed hyperbolic Pöschl-Teller potential and trigonometric Scarf II noncentral potential using Nikiforov-Uvarov method [6], -deformed trigonometric Scarf potential with -deformed Trigonometric Tensor Coupling Potential for Spin and Pseudospin Symmetries Using Romanovski Polynomial [7], generalized nuclear Wood-Saxon potential under relativistic spin symmetry limit [8], relativistic bound states of particle in Yukawa field with Coulomb tensor interaction [9], Rosen-Morse potential including the spin-orbit centrifugal term using Nikiforov-Uvarov (NU) method [3], pseudospin symmetric solution of the Morse potential for any state using AIM [11], Scalar, Vector, and Tensor Cornell Interaction using Ansatz method [12], Scalar and Vector Generalized Isotonic Oscillators and Cornell Tensor Interaction using Ansatz method [13], Mie-type potentials for energy dependent pseudoharmonic potential via SUSYQM [16], trigonometric Scarf potential in -dimension for spin and pseudospin symmetry using Nikiforov-Uvarov (NU) method [15], Coulombic potential and its thermodynamics properties in -dimensional space using NU method [17], and hyperbolic tangent potential and its application in material properties in -dimensional space [18].

Asymptotic Iteration Methods (AIM) have small deviation for determination of eigenenergies and eigenfunctions of Dirac equation. The separable -dimensional quantum potentials are not studied yet by some researchers. In this paper, we use Asymptotic Iteration Method (AIM) to solve the Dirac equation under influence of separable -dimensional quantum potentials. The relativistic energy levels can be obtained from calculation of relativistic energy equation using Matlab R2013a. In Section 2, we present basic theory of Dirac equation in hyperspherical coordinate with -dimensional separable quantum potential. In this section, deformed quantum potential is also included which is proposed by Dutra in 2005 [19]. In Section 3, we present Asymptotic Iteration Method. Result and discussion are included in Section 4, and in Section 5 we present the special case in 3-dimensional space. In the last section, we present conclusion.

2. Dirac Equation with Separable q-Deformed Quantum Potential in the Hyperspherical Coordinates

For single particle, Dirac equation with vector potential and scalar potential in the hyperspherical coordinate can be expressed as follows (in the unit ) [13, 14]:

where , , and are -dimensional momentum operator, total relativistic energy, and relativistic mass of the particle, respectively:

where are Pauli’s matrices and 1 is the unit matrix. Here we use relations between Pauli’s matrices as

The wave function of Dirac spinor can be classified in two forms, upper spinor and lower spinor as follows [15, 16]:

By substituting (2), (3), and (5) into (1), we get

Exact spin symmetry limit is characterized with , but for spin symmetry limit which is constant, and for pseudospin symmetry limit which is also constant. For spin symmetry limit we have since ; therefore, (6) can be rewritten as

If (8) is substituted into (7) we get

where

Using momentum operator definition in quantum mechanics, where , the hyperspherical Laplacian is given as [15]

The eigenvalue of is and the angular momentum operator is expressed as [15]

By inserting (11) into (9), we get

The separable variable potential used in this study is -deformed hyperbolic Rosen-Morse potential plus -deformed noncentral Scarf trigonometric potential in hyperspherical coordinate space. The effective potential can be written as

with angular potentials taken as

By substituting (14) into (13) and using variable separation method, we get the radial part and the angular part of Dirac equations in hyperspherical coordinate with .

2.1. The Radial Part

The -dimensional Dirac equation with -deformed hyperbolic Rosen-Morse potential plus -deformed trigonometric Scarf noncentral potentials can be resolved into the form of radial part and angular part equations. The radial part of -dimensional Dirac equation in this case can be expressed as

with .

2.2. The Angular Part

The angular part of 5-dimensional Dirac equation obtained from (13)-(14) can be resolved into four parts, and for , we get

Equations (20)–(23) are the angular part of Dirac equation for until , respectively. The -dimensional relativistic wave functions and orbital quantum numbers are obtained from those equations.

3. Review of Asymptotic Iteration Method (AIM)

Asymptotic Iteration Method (AIM) is an alternative method which has accuracy and high efficiency to determine eigenenergies and eigenfunctions for analytically solvable hyperbolic-like potential. Asymptotic Iteration Method is also giving solution for exactly solvable problem [17].

AIM is used to solve the second-order homogeneous linear equation as follows [5, 20–23]:

where and prime symbol refers to derivation along . The other parameter is interpreted as radial quantum number. Variables and are variables that can be differentiated along . To get the solution of (24), we have to differentiate (24) along , and then we get

where

Asymptotic Iteration Method (AIM) can be applied exactly in the different problem if the wave function has been known and fulfills boundary conditions zero () and infinity ().

Equation (2) can be simply iterated until () and (), and then we get

which is called as recurrence relation. Eigenvalue can be found using equation given as

where is the iteration number and is representation of radial quantum number.

Equation (2) is the second-order homogeneous linear equation which can be solved by comparing it with the second-order linear equation as follows [20]:

where

The wave function of (2) is the solution of (32) which is given as [20]

The -deformed hyperbolic and trigonometric functions are used as one of the parameters in the modified Rosen-Morse potential and noncentral Scarf trigonometric potentials were defined by Arai [24] some years ago as follows:

Deformation with -parameter in the hyperbolic function can be extended into trigonometric function. Definition of trigonometric function can be arranged by the same way as in the hyperbolic function introduced by Suparmi et al. [7] as follows:

By a convenient translation of spatial variable, one can transform the deformed potentials into the form of nondeformed potentials or vice-versa. In analogy to the translation of spatial variable for hyperbolic function introduced by Dutra [19], we propose the translation of spatial variable for hyperbolic and trigonometric function as follows:

And then by inserting (21) into (18) and (19), we have

The translation of spatial variable in (38) can be used to map the energy and wave function of nondeformed potential toward deformed potential of Scarf potential.

4. Result and Discussion

4.1. Radial Part

The radial part of Dirac equation in hyperspherical space when can be expressed as

Equation (39) cannot be solved directly when ; in this condition, we use approximation to solve centrifugal term with Pekeris approximation. Because we have used -deformed quantum potential, the Pekeris approximation in this condition can be expressed as [3]

with

By inserting (38) into (37a) and (37b) and changing from exponential form to hyperbolic form, we get

with

where is equilibrium distance that can be derived by potential parameter.

Let us assume

Equations (44) are inserted into (40) so we obtain

Let

By using (45)-(46), we get

By substituting (46) and into (47), we get

Equation (48) has two regular singular points for and so the general solution from (48) is .

Let and and then substitute them into (48) to get

Equation (49) is the second differential equation that can be manipulated into the form as in (22):

Let and using (29)-(30), together with (51) and (52), we get

The eigenvalue of (50) can be found by using (31):

where is th eigenvalue when and if is radial quantum number. By using (29)-(30) and the last equation in (55), we have

The last equation in (55) can be rewritten as

From (57), we get that the relativistic energy equation of this system is

By comparing (32) and (50) and using (34), we obtain as

From (59), we determine unnormalized radial wave function with for various as shown in Table 6, here is normalization constant.

4.2. Solution of Angular Part

In this study, the four angular parts of Dirac equations are presented in (20)–(23), so we have to solve each equation of angular Dirac equation using AIM.

4.2.1. Equation of Angular Part for

We can solve (20) by changing it into the form of the second-order hypergeometric-type differential equation that is similar to (22) after we insert (15) in (20),

and substitute into (60) so we get

Equation (62) has two regular singular points for and , and then the solution of is set as

If we replace in (61) with (62) and simplify it by using appropriate variable substitution as follows:

then (61) reduces to

From (64), we get

By using (65)–(68) and using relation of (31), we get

where is orbital quantum number. The angular wave function of (64) is determined by using (34) and then we get

where is normalization constant and is hypergeometry function. From (70) we solve unnormalized angular wave function as function of completely for variation of as shown in Table 7.

4.2.2. Equation of Angular Part for , , and

Solutions for , , and are determined by the same way as the solution for , by setting subscript with 2, 3, and 4. So we get the solution for orbital quantum number for , , and , respectively, as follows:

where

The unnormalized angular wave functions for , , and are in the same pattern as the angular wave function for so we obtain the last three angular wave functions by simple change of parameters in (70) with parameters in (74)–(79).

The relativistic energy levels are calculated numerically using Matlab program R2013a. Table 1 shows the relativistic energy levels as a function of deformation parameter ; the relativistic energy decreases when deformation parameter increases. Here we apply the value of from 0.2 until 1.2 with step 0.2. The relativistic energy levels as a function of orbital quantum numbers are also shown in Table 1. The magnitude of energies decreases when the orbital quantum numbers increase. The relativistic energy levels as a function of radial quantum number are shown in Table 2. The relativistic energies increase when radial quantum number increases. The relativistic energies numerically also change as a function of potential parameters and where are the components of th noncentral potential (Table 3). The relativistic energies increase when both potential parameters in each potential component increase. This suggests that the bounded energies become less bounded with increasing of potential parameters. The unnormalized angular wave functions are listed in Table 4. The unnormalized radial wave functions are plotted by using (46) and (59) as shown in Figure 1. From Figures 1(a)–1(c), it is seen that the amplitude of the wave function increases when the orbital quantum number increases. This suggests that the probability of finding particles is larger for higher orbital quantum number .

Table 1: Relativistic energies with variation of and .

, , ,

(fm)

(fm⁻¹)

(fm)
(fm⁻¹)

0.2
0.0333
1
1
−5.2311
0
1
1
0.1671
−5.0977

0.4
0.0667
1
1
−5.4766
1
1
1
0.1671
−6.4721

0.6
0.1001
1
1
−5.7371
2
1
1
0.1671
−9.1544

0.8
0.1335
1
1
−6.0240
3
1
1
0.1671
−13.1643

1
0.1671
1
1
−6.4721
4
1
1
0.1671
−18.5283

1.2
0.2007
1
1
−6.9767
5
1
1
0.1671
−25.2837

Table 2: Relativistic energies with variation of radial quantum number .

, , , (fm⁻¹), and (fm)

(fm⁻¹)

0
1
1
−6.4725

1
1
1
−6.4721

2
1
1
−6.4717

3
1
1
−6.4713

4
1
1
−6.4709

Table 3: The relativistic energies with variation of potential parameters and with ; ; ; ; ; ; ; ; and .

(fm)

2
2
2
2
2
2
2
2
−5.0980

4
4
2
2
2
2
2
2
−5.0821

6
6
2
2
2
2
2
2
−5.0561

8
8
2
2
2
2
2
2
−5.0382

10
10
2
2
2
2
2
2
−5.0269

2
2
4
4
2
2
2
2
−5.0528

2
2
6
6
2
2
2
2
−5.0319

2
2
8
8
2
2
2
2
−5.0211

2
2
10
10
2
2
2
2
−5.0171

2
2
2
2
4
4
2
2
−5.0383

2
2
2
2
6
6
2
2
−5.0136

2
2
2
2
8
8
2
2
−5.0084

2
2
2
2
10
10
2
2
−5.0055

2
2
2
2
2
2
4
4
−5.0453

2
2
2
2
2
2
6
6
−5.0266

2
2
2
2
2
2
8
8
−5.0200

2
2
2
2
2
2
10
10
−5.0156

Table 4: The unnormalized upper angular wave function for higher dimension .

Unnormalized angular wave function

0

1

(a)

(b)

(c)

Figure 1: Unnormalized radial wave functions versus hyperspherical radius with variation of radial quantum number, where , , , , , and , . (a) For , (b) for , and (c) for .

5. Special Case for Three-Dimensional System

In three-dimensional case, , the bound state relativistic energy levels are calculated numerically from (58), (43), (44), (56), (69), (74), and (77) using Matlab software R2013a and are presented in Tables 5, 6, and 7. By using (44) and (62), the unnormalized radial wave functions are listed in Table 8. In Table 5, the relativistic energy levels decrease with the increase of both deformation parameter and orbital quantum number , where is orbital quantum number for and is orbital quantum number for . From Table 7, we can conclude that the value of relativistic energies increases with the increase of the radial quantum number and with either the presence or absence of angular potential parameter. The angular potential parameters and influence the relativistic energy level, where the relativistic energy levels increase by increasing the angular potential parameter. The unnormalized upper radial wave functions are listed in Table 8 when the radial quantum numbers are and and is the normalization constant. The unnormalized angular wave functions as functions of and are listed in Table 9. Without the presence of the noncentral potential, degeneracy energy spectra occur in the spin doublets, with quantum numbers () and (), where , , and are the radial, the orbital, and the total angular momentum quantum numbers, respectively, for example, (, ) for but for difference values of and , (, ) for with ( and , (, ) for with ( and , and (, ) for with ( and . The degeneracy energies can be removed with the presence of noncentral potential, by changing , where and for and .

Table 5: Relativistic energy with variation of and .

, , (fm⁻¹)  

(fm)

(fm⁻¹)

(fm)
(fm⁻¹)

0.2
0.0333
1
1
−5.1360
0
1
1
0.1671
−5.1394

0.4
0.0667
1
1
−5.2867
1
1
1
0.1671
−6.1034

0.6
0.1001
1
1
−5.4552
2
1
1
0.1671
−7.9374

0.8
0.1335
1
1
−5.6616
3
1
1
0.1671
−10.6614

1
0.1671
1
1
−6.1034
4
1
1
0.1671
−14.3053

1.2
0.2007
1
1
−6.5579

Table 6: Relativistic energies with variation of and potential parameter.

, , , , (fm⁻¹), and = 0.1671 (fm)

(fm⁻¹)

(fm⁻¹)

0
1
2
−6.1042
0
1
0
−0.0616

1
1
2
−6.1034
1
1
0
−0.0564

2
1
2
−6.1024
2
1
0
−0.0509

3
1
2
−6.1014
3
1
0
−0.0454

4
1
2
−6.1004
4
1
0
−0.0398

Table 7: Relativistic energies for  fm⁻¹, , , and without presence of -deformed noncentral potentials.

(fm⁻¹), ,
(fm⁻¹), ,
(fm⁻¹), ,

1
0
−2
0p3/2
−0.5028
−0.5014
−0.5006

2
0
−3
0d5/2
−0.5008
−0.5002
−0.4999

3
0
−4
0f7/2
−0.5002
−0.4999
−0.4996

4
0
−5
0g9/2
−0.5000
−0.4998
−0.4995

1
1
−2
1p3/2
−0.5042
−0.5028
−0.5020

2
1
−3
1d5/2
−0.5013
−0.5007
−0.5003

3
1
−4
1f7/2
−0.5005
−0.5002
−0.4999

4
1
−5
1g9/2
−0.5002
−0.4999
−0.4997

1
0
1
0p1/2
−0.5028
−0.5014
−0.5006

2
0
2
0d3/2
−0.5008
−0.5002
−0.4999

3
0
3
0f5/2
−0.5002
−0.4999
−0.4996

4
0
4
0g7/2
−0.5000
−0.4998
−0.4995

1
1
1
1p1/2
−0.5042
−0.5028
−0.5020

2
1
2
1d3/2
−0.5013
−0.5007
−0.5003

3
1
3
1f5/2
−0.5005
−0.5002
−0.4999

4
1
4
1g7/2
−0.5002
−0.4999
−0.4997

Table 8: Unnormalized upper radial wave function for dimension .

Radial wave function

0

1

Table 9: The unnormalized upper angular wave function for dimension .

Unnormalized angular wave function

0

1

6. Conclusion

In this study, we have obtained the bound state solution of the 5-dimensional Dirac equation with separable -deformed quantum potential under condition of spin symmetry. The upper component of Dirac spinors and relativistic energy have been obtained using Asymptotic Iteration Method (AIM). The numerical result shows that the bound state relativistic energy level for spin symmetry case with both dimensions and increases with increasing of radial quantum number and decreases with increasing of both deformation parameter and orbital quantum number . The degeneracy energy states occur in the absence of -deformed noncentral potential.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This Research is partly supported by Higher Education Project Grant with contract no. 632/UN27/27.21/LT/2016.

References

W. Greiner, Relativistic Quantum Mechanics, Springer, New York, NY, USA, 3rd edition, 2000.

B. J. Falaye and K. J. Oyewumi, “Solutions of the Dirac equation with spin and pseudospin symmetry for the trigonometric Scarf potential in D-dimensions,” The African Review of Physics, vol. 6, article 0025, 2011.

S. M. Ikhdair, “Approximate solutions of the Dirac equation for the Rosen-Morse potential including the spin-orbit centrifugal term,” Journal of Mathematical Physics, vol. 51, no. 2, Article ID 023525, 2010.

J. N. Ginocchio, “Pseudospin as a relativistic symmetry,” Physical Review Letters, vol. 78, no. 3, pp. 436–439, 1997.

R. A. Sari, A. Suparmi, and C. Cari, “Solution of Dirac equation for Eckart potential and trigonometric Manning Rosen potential using asymptotic iteration method,” Chinese Physics B, vol. 25, no. 1, Article ID 010301, 2015.

A. Kurniawan, A. Suparmi, and C. Cari, “Approximate analytical solution of the Dirac equation with q-deformed hyperbolic Pöschl-Teller potential and trigonometric Scarf II non-central potential,” Chinese Physics B, vol. 24, no. 3, Article ID 030302, 2015.

A. Suparmi, C. Cari, U. A. Deta, A. S. Husein, and H. Yuliani, “Exact solution of Dirac equation for q-deformed trigonometric Scarf potential with q-deformed trigonometric tensor coupling potential for spin and pseudospin symmetries using Romanovski polynomial,” Journal of Physics: Conference Series, vol. 539, no. 1, Article ID 012004, 2014.

M. Hamzavi and A. A. Rajabi, “Generalized nuclear Woods-Saxon potential under relativistic spin symmetry limit,” ISRN High Energy Physics, vol. 2013, Article ID 987632, 7 pages, 2013.

M. Sameer, “Approximate relativistic bound states of a particle in Yukawa field with Coulomb tensor interaction,” Physica Scripta, vol. 87, no. 3, Article ID 035002, 2013.

Suparmi and Cari, “Bound state solution of Dirac equation for generalized Pöschl-Teller plus trigomometric Pöschl-Teller non-central potential using SUSY quantum mechanics,” Journal of Mathematical and Fundamental Sciences, vol. 46, no. 3, pp. 205–223, 2014.

O. Bayrak and I. Boztosun, “The pseudospin symmetric solution of the Morse potential for any κ state,” Journal of Physics A, vol. 40, no. 36, pp. 11119–11127, 2007.

H. Hassanabadi, E. Maghsoodi, S. Zarrinkamar, and H. Rahimov, “Dirac equation under scalar, vector, and tensor cornell interactions,” Advances in High Energy Physics, vol. 2012, Article ID 707041, 17 pages, 2012.

H. Hassanabadi, E. Maghsoodi, A. N. Ikot, and S. Zarrinkamar, “Dirac equation under scalar and vector generalized isotonic oscillators and cornell tensor interaction,” Advances in High Energy Physics, vol. 2014, Article ID 831938, 7 pages, 2014.

A. Lopez-Ortega, “The Dirac equation in D-dimensional spherically symmetric spacetimes,” https://arxiv.org/abs/0906.2754

R. L. Hall and P. Zorin, “Refined comparison theorems for the Dirac equation in d dimensions,” Annalen der Physik, vol. 527, no. 5-6, pp. 408–422, 2015.

A. N. Ikot, H. Hassanabadi, E. Maghsoodi, and S. Zarrinkamar, “D-dimensional Dirac equation for energy-dependent pseudoharmonic and Mie-type potentials via SUSYQM,” Communications in Theoretical Physics, vol. 61, no. 4, pp. 436–446, 2014.

C. Cari, A. Suparmi, M. Yunianto, and A. S. Husein, “Solution of D-dimensional Dirac equation for Coulombic potential using NU method and its thermodynamics properties,” AIP Conference Proceeding, vol. 1710, Article ID 030009, 2016.

A. Suparmi, C. Cari, B. N. Pratiwi, and U. A. Deta, “Solution of D dimensional Dirac equation for hyperbolic tangent potential using NU method and its application in material properties,” AIP Conference Proceedings, vol. 1710, Article ID 030010, 2016.

A. S. Dutra, “Mapping deformed hyperbolic potentials into nondeformed ones,” Physics Letters A, vol. 339, no. 3–5, pp. 252–254, 2005.

T. Das, “Exact solutions of the Klein-Gordon equation for q-deformed manning-rosen potential via asymptotic iteration method,” https://arxiv.org/abs/1409.1457v1

A. Soylu, O. Bayrak, and I. Boztosun, “An approximate solution of Dirac-Hulthén problem with pseudospin and spin symmetry for any κ state,” Journal of Mathematical Physics, vol. 48, no. 8, Article ID 082302, 2007.

A. Soylu, O. Bayrak, and I. Boztosun, “K state solutions of the Dirac equation for the Eckart potential with pseudo-spin and spin symmetry,” Journal of Physics A: Mathematical and Theoretical, vol. 41, Article ID 065308, 2008.

G. Kocak, O. Bayrak, and I. Boztosun, “Supersymmetric solution of Schrödinger equation by using the asymptotic iteration method,” Annalen der Physik, vol. 524, no. 6-7, pp. 353–359, 2012.

B.-H. Nam, J.-Y. Moon, E.-H. Park et al., “Conserved gene structure and function of interleukin-10 in teleost fish,” Journal of Animal and Veterinary Advances, vol. 13, no. 12, pp. 774–782, 2014.