Advances in High Energy Physics

Volume 2012 (2012), Article ID 873619, 15 pages

doi:10.1155/2012/873619

Research Article

Exact Solutions of a Spatially Dependent Mass Dirac Equation for Coulomb Field plus Tensor Interaction via Laplace Transformation Method

M. Eshghi,¹ M. Hamzavi,² and S. M. Ikhdair3,4

¹Physics Department, Imam Hossein Comprehensive University, Tehran, Iran

²Department of Basic Sciences, Shahrood Branch, Islamic Azad University, Shahrood, Iran

³Physics Department, Near East University, 922022 Nicosia, Northern Cyprus, Turkey

⁴Physics Department, Faculty of Science, An-Najah National University, Nablus, West Bank, Palestine

Received 9 October 2012; Accepted 4 November 2012

Academic Editor: Kingman Cheung

Copyright © 2012 M. Eshghi 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.

Abstract.

The spatially dependent mass Dirac equation is solved exactly for attractive scalar and repulsive vector Coulomb potentials including a tensor interaction potential under the spin and pseudospin (p-spin) symmetric limits by using the Laplace transformation method (LTM). Closed forms of the energy eigenvalue equation and wave functions are obtained for arbitrary spin-orbit quantum number . Some numerical results are given too. The effect of the tensor interaction on the bound states is presented. It is shown that the tensor interaction removes the degeneracy between two states in the spin doublets. We also investigate the effects of the spatially-dependent mass on the bound states under the conditions of the spin symmetric limit and in the absence of tensor interaction .

1. Introduction

Dirac equation has become one of the most appealing relativistic wave equations for spin-1/2 particles. However, solving such a wave equation is still a very challenging problem even if it has been derived more than 80 years ago and has been utilized profusely. It is always useful to investigate the relativistic effects [1–4]. For example, in the relativistic treatment of nuclear phenomena the Dirac equation is used to describe the behavior of the nuclei in nucleus and also in solving many problems of high-energy physics and chemistry. For this reason, it has been used extensively to study the relativistic heavy ion collisions, heavy ion spectroscopy, and more recently in laser-matter interaction (for a review, see [5] and references therein) and condensed matter physics [6].

On the other hand, systems with position-dependent mass (PDM) have been found to be very useful in studying the physical properties of various microstructures [7–14]. Recently, there has been increased interest in searching for analytical solutions of the Dirac equation with PDM and with constant mass under the spin and p-spin symmetries [15–34].

Here, we shall attempt to solve the Dirac equation by using the Laplace transform method (LTM). The LTM is an integral transform and is comprehensively useful in physics and engineering [35] and recently used by many authors to solve the Schrödinger equation for different potential forms [36–40]. This method could be a nearly new formalism in the literature and serve as a powerful algebraic treatment for solving the second-order differential equations. As a result, the LTM describes a simple way for solving of radial and one-dimensional differential equations. The other advantage of this method is that a second-order equation can be converted into more simpler form whose solutions may be obtained easily [36]. In this paper, we obtain solution of the Dirac equation both PDM and tensor interaction for attractive scalar and repulsive vector Coulomb potential under the spin and p-spin symmetry limits. We give some numerical results.

2. Review to Dirac Equation including Tensor Coupling

The Dirac equation which describes a nucleon in repulsive vector and attractive scalar and a tensor potentials is written as

where is the effective mass of the fermionic particle, is the relativistic energy of the system, is the three-dimensional momentum operator. and are the Dirac matrices give as

where is unitary matrix and are three-vector spin matrices

The total angular momentum operator and spin-orbit , where is orbital angular momentum, of the spherical nucleons commute with the Dirac Hamiltonian. The eigenvalues of spin-orbit coupling operator are and for unaligned spin and the aligned spin , respectively. () can be taken as the complete set of the conservative quantities. Thus, the spinor wave functions can be classified according to their angular momentum , spin-orbit quantum number and the radial quantum number can be written as follows:

where is the upper (large) component and is the lower (small) component of the Dirac spinors. and are spin and p-spin spherical harmonics, respectively, and is the projection of the angular momentum on the -axis. Substituting (2.4) into (2.1) and using the following relations [41] as

and with the following properties:

one obtains two coupled differential equations for upper and lower radial wave functions and as

where

Eliminating and from (2.5e) and (2.5f), we finally obtain the following two Schrödinger-like differential equations for the upper and lower radial spinor components, respectively:

where and . These radial wave functions are required to satisfy the necessary boundary conditions. The spin-orbit quantum number is related to the quantum numbers for spin symmetry and p-spin symmetry as

and the quasi-degenerate doublet structure can be expressed in terms of a p-spin angular momentum and pseudo-orbital angular momentum , which is defined as

where . For example, and can be considered as p-spin doublets. In the next section, we will consider the spin and p-spin symmetry cases.

3. Relativistic Bound State Solutions

3.1. Spin Symmetry Case

Equation (2.5i) cannot be solved analytically because of the last term . Therefore, in solving the mathematical relation [42, 43], we can then exactly solve (2.5i). In this stage, we take the vector potential in the form of an attractive Coulomb-like field [18] as

where is being a vector dimensionless real parameter coupling constant and is being a constant with dimension. Also, it is convenient to take the mass function [18] as

where and stand for the rest mass of the fermionic particle and the perturbed mass, respectively. Further, is the dimensionless real constant to be set to zero for the constant mass case and is the Compton-like wavelength in fm units. Further, the tensor interaction takes the simple form:

where is the coulomb radius, and stand for the charges of the projectile particle and the target nucleus , respectively.

Substituting (3.1)–(3.3) into (2.5i) considering the spin symmetry case where ., that is, [44, 45]. Thus, the equation obtained for the upper component of the Dirac spinor becomes

Further, defining the new parameters

and introducing , then (3.4) turns into the form

Setting with is a constant and then inserting into (3.6a), we have

Now, to obtain a finite wave function when , if we take in (3.7) then it becomes

The LTM [46, 47]

leads to an equation

Equation (3.10) is a first-order differential equation and therefore we may directly make use of the integral to get the expression

where is a constant. Noting that is a multi-valued function and the wave functions are required to be single-valued, we must take

which gives single-values wave functions. Using this requirement and further expanding (3.11) into series, we obtain

where is a constant. In terms of a simple extension of the inverse Laplace transformation [46, 47], we can immediately obtain

and from , we then obtain

where is a constant. On the other hand, the confluent hyper-geometric functions is defined as a series expansion [48]

So, we find the upper-spinor component of wave function as

where is normalization constant. By using the normalization condition given as , and the relation between the Laguerre polynomials and confluent hyper-geometric functions as [47], the normalization constant in (3.17) is written as

where we have used [47]

Inserting the parameters in (3.5) and (3.6b) into (3.12), one obtain the transcendental energy eigenvalue equation as follows

which is identical to [18] when and the tensor interaction is removed, that is, . For the constant mass case, we have

In case when ,  ,  ,  ,   and , the above equation reduces to the non-relativistic energy limit (in units of ):

The wave function (3.17) is finite in the entire range and at the boundaries, that is, and .

The lower spinor wave function can be obtained via

where we have used

In Tables 1, 2, and 3 with , we give some numerical results for the energy eigenvalues from energy formula (3.20).

Table 1: The bound state energy eigenvalues of the Coulomb potential under the pseudo-spin symmetry limit for several values of and for , , and .

1
1, −2

0.0
−4.86154
−4.75676
−4.47059
−4.86154

0.2
−4.80508
−4.66644
−4.36211
−4.80508

0.4
−4.74207
−4.57172
−4.26086
−4.74207

0.5
−4.70859
−4.52376
−4.21286
−4.70859

2
1, −3

0.0
−4.91089
−4.86154
−4.75676
−4.75676

0.2
−4.87335
−4.80508
−4.66644
−4.66644

0.4
−4.83043
−4.74207
−4.57172
−4.57172

0.5
−4.80716
−4.70859
−4.52376
−4.52376

3
1, −4

0.0
−4.93793
−4.91089
−4.86154
−4.47059

0.2
−4.91138
−4.87335
−4.80508
−4.36211

0.4
−4.88067
−4.83043
−4.74207
−4.26086

0.5
−4.86385
−4.80716
−4.70859
−4.21286

4
1, −5

0.0
−4.95431
−4.93793
−4.91089
−4.47059

0.2
−4.93461
−4.91138
−4.87335
−4.36211

0.4
−4.91166
−4.88067
−4.83043
−4.26086

0.5
−4.89903
−4.86385
−4.80716
−4.21286

1
2, −2

−4.91089
−4.86154
−4.75676
−4.91089

−4.87402
−4.80794
−4.69083
−4.87402

−4.83243
−4.75015
−4.62526
−4.83243

−4.81014
−4.72024
−4.59287
−4.81014

2
2, −3

−4.93793
−4.91089
−4.86154
−4.86154

−4.9116
−4.87402
−4.80794
−4.80794

−4.88134
−4.83243
−4.75015
−4.75015

−4.86487
−4.81014
−4.72024
−4.72024

3
2, −4

−4.95431
−4.93793
−4.91089
−4.75676

−4.93469
−4.9116
−4.87402
−4.69083

−4.91194
−4.88134
−4.83243
−4.62526

−4.89945
−4.86487
−4.81014
−4.59287

4
2, −5

−4.96498
−4.95431
−4.93793
−4.75676

−4.94984
−4.93469
−4.91160
−4.69083

−4.93218
−4.91194
−4.88134
−4.62526

−4.92244
−4.89945
−4.86487
−4.59287

Table 2: The bound state energy eigenvalues of the Coulomb potential under the pseudo-spin symmetry limit for several values of and for , , and .

2
0, 2

0.0
−4.75676
−4.86154
−4.91089
−4.96498

0.2
−4.66644
−4.80508
−4.87335
−4.94979

0.4
−4.57172
−4.74207
−4.83043
−4.93205

0.5
−4.52376
−4.70859
−4.80716
−4.92225

3
0, 3

0.0
−4.86154
−4.91089
−4.93793
−4.97231

0.2
−4.80508
−4.87335
−4.91138
−4.96026

0.4
−4.74207
−4.83043
−4.88067
−4.94614

0.5
−4.70859
−4.80716
−4.86385
−4.93833

4
0, 4

0.0
−4.91089
−4.93793
−4.95431
−4.97756

0.2
−4.87335
−4.91138
−4.93461
−4.96777

0.4
−4.83043
−4.88067
−4.91166
−4.95628

0.5
−4.80716
−4.86385
−4.89903
−4.94991

5
0, 5

0.0
−4.93793
−4.95431
−4.96498
−4.98144

0.2
−4.91138
−4.93461
−4.94979
−4.97334

0.4
−4.88067
−4.91166
−4.93205
−4.96381

0.5
−4.86385
−4.89903
−4.92225
−4.95852

2
1, 2

0.0
−4.86154
−4.91089
−4.93793
−4.97231

0.2
−4.80794
−4.87402
−4.91160
−4.96028

0.4
−4.75015
−4.83243
−4.88134
−4.94621

0.5
−4.72024
−4.81014
−4.86487
−4.93843

3
1, 3

0.0
−4.91089
−4.93793
−4.95431
−4.97756

0.2
−4.87402
−4.91160
−4.93469
−4.96778

0.4
−4.83243
−4.88134
−4.91194
−4.95632

0.5
−4.81014
−4.86487
−4.89945
−4.94997

4
1, 4

0.0
−4.93793
−4.95431
−4.96498
−4.98144

0.2
−4.91160
−4.93469
−4.94984
−4.97757

0.4
−4.88134
−4.91194
−4.93218
−4.96383

0.5
−4.86487
−4.89945
−4.92244
−4.95856

5
1, 5

0.0
−4.95431
−4.96498
−4.97231
−4.98440

0.2
−4.93469
−4.94984
−4.96028
−4.97758

0.4
−4.91194
−4.93218
−4.94621
−4.96957

0.5
−4.89945
−4.92244
−4.93843
−4.96512

Table 3: The bound state energy eigenvalues of the Coulomb potential for several states under the pseudo-spin symmetry limit with values of with , , and .

1

−1.00
−4.75676
−4.52376
−4.96498
−4.92225

−1.25
−4.76351
−4.53699
−4.96595
−4.92441

−1.50
−4.77027
−4.55022
−4.96693
−4.92657

−1.75
−4.77703
−4.56345
−4.96790
−4.92873

−2.00
−4.78378
−4.57668
−4.96887
−4.93089

−2.25
−4.79054
−4.58991
−4.96984
−4.93305

−2.50
−4.79730
−4.60313
−4.97082
−4.93521

−2.75
−4.80405
−4.61636
−4.97179
−4.93736

3

−1.00
−4.91089
−4.80543
−3.20000
−2.90681

−1.25
−4.91337
−4.81083
−3.25000
−2.96496

−1.50
−4.91584
−4.81624
−3.30000
−3.0231

−1.75
−4.91832
−4.82164
−3.35000
−3.08124

−2.00
−4.92079
−4.82705
−3.40000
−3.13939

−2.25
−4.92327
−4.83245
−3.45000
−3.19753

−2.50
−4.92574
−4.83786
−3.50000
−3.25568

−2.75
−4.92822
−4.84326
−3.55000
−3.31382

1

−1.00
−4.86154
−4.70859
−4.86154
−4.70859

−1.25
−4.86538
−4.71669
−4.86538
−4.71669

−1.50
−4.86923
−4.72478
−4.86923
−4.72478

−1.75
−4.87308
−4.73288
−4.87308
−4.73288

−2.00
−4.87692
−4.74097
−4.87692
−4.74097

−2.25
−4.88077
−4.74907
−4.88077
−4.74907

−2.50
−4.88462
−4.75716
−4.88462
−4.75716

−2.75
−4.88846
−4.76526
−4.88846
−4.76526

3

−1.00
−4.93793
−4.86385
−4.47059
−4.21286

−1.25
−4.93966
−4.86764
−4.48529
−4.23472

−1.50
−4.94138
−4.87142
−4.50000
−4.25659

−1.75
−4.9431
−4.8752
−4.51471
−4.27845

−2.00
−4.94483
−4.87898
−4.52941
−4.30032

−2.25
−4.94655
−4.88276
−4.54412
−4.32218

−2.50
−4.94828
−4.88654
−4.55882
−4.34405

−2.75
−4.95000
−4.89033
−4.57353
−4.36591

3.2. p-Spin Symmetry Case

To avoid repetition in the solution of (2.5j), the negative energy solution for p-spin symmetry can be obtained directly from those of the above positive energy solution for spin symmetry by using the parameter mapping [18]:

Following the previous results with the above transformations, we finally arrive at the transcendental energy equation as

which is identical to [18] when and the tensor interaction is removed, that is, . For the constant mass case, we have

the lower-spinor component of wave function as

where , and is normalization constant

The upper-spinor wave function can be obtained via

4. Numerical Results

In Tables 1 to 3, we see that energies of bound states such as: , , , (where each pair is considered as a spin doublet) in the absence of the tensor potential are degenerate but in the presence of the tensor potential, the degeneracies are removed. Also, we investigate the effects of the and parameters on the bound states under the conditions of the spin symmetry limits for . The results are given in Tables 1 to 3. It is readily seen that if and parameters increases, the value of the bound state energy eigenvalues of this potential increases for several states. It is shown that the energy eigenvalues decrease with decreasing and when increase with increasing the tensor strength for and . In Table 3, for constant values of and , the energy increases when increasing. The decrease in the energy values is large without tensor interaction while small in presence of tensor interaction and then being large. Further, when increases, the energy increasing.

Finally, we plot the relativistic energy eigenvalues under spin and p-spin symmetry limitations in Figures 1 to 4. In Figure 2, we plot the energy eigenvalues of spin symmetry limit versus the perturbated mass . It is seen that when increases, the energy increases too. In Figure 2, we have shown the variation of the energy as a function of . We can see the degeneracy removes between spin doublets and also they become far from each other, when the parameter increases. In Figures 3 to 4, we plot the energy states of the pseudospin symmetry limit for different levels as functions of parameters and , respectively. The variation of energy can also be seen from these figures.

Figure 1: The variation of the energy levels as a function in the presence of spin symmetry.

Figure 2: The variation of the energy levels as a function in the presence of spin symmetry.

Figure 3: The variation of the energy levels as a function in the presence of pseudospin symmetry.

Figure 4: The variation of the energy levels as a function in the presence of pseudospin symmetry.

5. Conclusion

In this paper, the relativistic equation for particles with spin 1/2 was solved exactly with both spatially-dependent mass and tensor interaction for attractive scalar and repulsive vector Coulomb potentials under the spin symmetry limit via the Laplace transformation method. Some numerical results are given for specific values of the model parameters. Effects of the tensor interaction on the bound states were presented that tensor interaction removes degeneracy between two states in spin doublets. We also investigated the effects of the spatially-dependent mass on the bound states under the conditions of the spin symmetry limits for .

References

I. C. Wang and C. Y. Wong, “Finite-size effect in the Schwinger particle-production mechanism,” Physical Review D, vol. 38, pp. 348–359, 1988.

G. Mao, “Effect of tensor couplings in a relativistic Hartree approach for finite nuclei,” Physical Review C, vol. 67, Article ID 044318, 12 pages, 2003.

R. J. Furnstahl, J. J. Rusnak, and B. D. Serot, “The nuclear spin-orbit force in chiral effective field theories,” Nuclear Physics A, vol. 632, no. 4, pp. 607–623, 1998.

P. Alberto, R. Lisboa, M. Malheiro, and A. S. de Castro, “Tensor coupling and pseudospin symmetry in nuclei,” Physical Review C, vol. 71, Article ID 034313, 7 pages, 2005.

Y. I. Salamin, S. X. Hu, K. Z. Hatsagortsyan, and C. H. Keitel, “Relativistic high-power laser-matter interactions,” Physics Reports, vol. 427, no. 2-3, pp. 41–155, 2006.

M. I. Katsnelson, K. S. Novoselov, and A. K. Geim, “Chiral tunnelling and the Klein paradox in graphene,” Nature Physics, vol. 2, no. 9, pp. 620–625, 2006.

M. R. Galler and W. Kohn, “Quantum mechanics of electrons in crystals with graded composition,” Physical Review Letters, vol. 70, pp. 3103–3106, 1993.

A. Puente and M. Gasas, “Non-local energy density functional for atoms and metal clusters,” Computational Materials Science, vol. 2, pp. 441–449, 1994.

F. Arias De Saavedra, J. Boronat, A. Polls, and A. Fabrocini, “Effective mass of one ⁴He atom in liquid ³He,” Physical Review B, vol. 50, no. 6, pp. 4248–4251, 1994.

M. Barranco, M. Pi, S. M. Gatica, E. S. Hernandez, and J. Navarro, “Structure and energetics of mixed ⁴He-³He drops,” Physical Review B, vol. 56, pp. 8997–9003, 1997.

A. Puente, L. I. Serra, and M. Casas, “Dipole excitation of Na clusters with a non-local energy density functional,” Zeitschrift Für Physik D, vol. 31, pp. 283–286, 1994.

L. Serra and E. Lipparini, “Spin response of unpolarized quantum dots,” Europhysics Letters, vol. 40, no. 6, pp. 667–672, 1997.

G. Bastard, Wave Mechanics Applied to Semiconductor Heterostructures, Les Editions de physique, Les Ulis, France, 1988.

A. G. M. Schmidt, “Wave-packet revival for the Schrödinger equation with position-dependent mass,” Physics Letters A, vol. 353, pp. 459–462, 2006.

M. Eshghi and H. Mehraban, “Solution of the dirac equation with position-dependent mass for q-parameter modified pöschl-teller and coulomb-like tensor potential,” Few-Body Systems, vol. 52, pp. 41–47, 2012.

E. Maghsoodi, H. Hassanabadi, and S. Zarrinkamar, “Spectrum of dirac equation under deng-fan scalar and vector potentials and a coulomb tensor interaction by SUSYQM,” Few-Body Systems, vol. 53, no. 3-4, pp. 525–538, 2012.

C. S. Jia and A. de Souza Dutra, “Extension of PT-symmetric quantum mechanics to the Dirac theory with position-dependent mass,” Annals of Physics, vol. 323, no. 3, pp. 566–579, 2008.

S. M. Ikhdair and R. Sever, “Solutions of the spatially-dependent mass Dirac equation with the spin and pseudospin symmetry for the Coulomb-like potential,” Applied Mathematics and Computation, vol. 216, no. 2, pp. 545–555, 2010.

Y. Xu, S. He, and C.-S. Jia, “Approximate analytical solutions of the Dirac equation with the Pöschl-Teller potential including the spin-orbit coupling term,” Journal of Physics A, vol. 41, no. 25, p. 255302, 8, 2008.

H. Akcay, “Dirac equation with scalar and vector quadratic potentials and Coulomb-like tensor potential,” Physics Letters A, vol. 373, no. 6, pp. 616–620, 2009.

M. Eshghi and M. Hamzavi, “Spin symmetry in dirac-attractive radial problem and tensor potential,” Communications in Theoretical Physics, vol. 57, article 355, 2012.

M. Eshghi, “Pseudo-harmonic oscillatory ring-shaped potential in a relativistic equation,” Chinese Physics Letters, vol. 29, no. 11, Article ID 110304, 2012.

O. Panella, S. Biondini, and A. Arda, “New exact solution of the one-dimensional Dirac equation for the Woods-Saxon potential within the effective mass case,” Journal of Physics A, vol. 43, Article ID 325302, 2010.

A. D. Alhaidari, “Solution of the Dirac equation with position-dependent mass in the Coulomb field,” Physics Letters A, vol. 322, no. 1-2, pp. 72–77, 2004.

X.-L. Peng, J.-Y. Liu, and C.-S. Jia, “Approximation solution of the Dirac equation with position-dependent mass for the generalized Hulthén potential,” Physics Letters A, vol. 352, no. 6, pp. 478–483, 2006.

M. Eshghi and H. Mehraban, “Eigen spectra for Manning-Rosen potential including a Coulomb-like tensor interaction,” International Journal of Physical Sciences, vol. 6, no. 29, pp. 6643–6652, 2011.

M. Hamzavi, M. Eshghi, and S. M. Ikhdair, “Effect of tensor interaction in the Dirac-attractive radial problem under pseudospin symmetry limit,” Journal of Mathematical Physics, vol. 53, Article ID 082101, 10 pages, 2012.

M.-C. Zhang and G.-Q. Huang-Fu, “Solution of the Dirac equation in the tridiagonal representation with pseudospin symmetry for an anharmonic oscillator and electric dipole ring-shaped potential,” Annals of Physics, vol. 327, no. 3, pp. 841–850, 2012.

O. Aydogdu and R. Sever, “Exact solution of the Dirac equation with the Mie-type potential under the pseudospin and spin symmetry limit,” Annals of Physics, vol. 325, pp. 373–383, 2010.

M. Eshghi and H. Mehraban, “Eigen spectra in the Dirac-hyperbolic problem with tensor coupling,” Chinese Journal of Physics, vol. 50, no. 4, p. 533, 2012.

M. Eshghi, “Dirac-hyperbolic scarf problem including a coulomb-like tensor potential,” Acta Scientiarum, vol. 34, p. 207, 2012.

G. F. Wei and S. H. Dong, “A novel algebraic approach to spin symmetry for Dirac equation with scalar and vector second Pöschl-Teller potentials,” The European Physical Journal A, vol. 43, pp. 185–190, 2010.

V. G. C. S. dos Santos, A. de Souza Dutra, and M. B. Hott, “Real spectra for the non-Hermitian Dirac equation in dimensions with the most general coupling,” Physics Letters A, vol. 373, no. 38, pp. 3401–3406, 2009.

M. Hamzavi, A. A. Rajabi, and H. Hassanabadi, “Exact spin and pseudospin symmetry solutions of the Dirac equation for mie-type potential including a coulomb-like tensor potential,” Few-Body Systems, vol. 48, no. 2, pp. 171–182, 2010.

G. Doetsch, Guide to the Applications of Laplace Transforms, Princeton University press, Princeton, NJ, USA, 1961.

G. Chen, “The exact solutions of the Schrödinger equation with the Morse potential via Laplace transforms,” Physics Letters A, vol. 326, no. 1-2, pp. 55–57, 2004.

A. Arda and R. Sever, “Exact solutions of the Schrödinger equation via Laplace transform approach: pseudoharmonic potential and Mie-type potentials,” Journal of Mathematical Chemistry, vol. 50, no. 4, pp. 971–980, 2012.

E. Schrodinger, “Quantisierung als Eigenwertproblem. I and II,” Annalen Der Physik, vol. 79, pp. 361–376, 1926.

R. A. Swainson and G. W. F. Drake, “A unified treatment of the nonrelativistic and relativistic hydrogen atom. I. The wavefunctions,” Journal of Physics A, vol. 24, no. 1, pp. 79–94, 1991.

G. Chen, “Exact solutions of N-dimensional harmonic oscillator via Laplace transformation,” Chinese Physics, vol. 14, no. 6, p. 1075, 2005.

J. D. Bjorken and S. D. Drell, Relativistic Quantum Mechanics, McGraw-Hill, New York, NY, USA, 1964.

D. Agboola, “Dirac-Hulthen Problem with Position-dependent Mass in D-dimensions,” http://arxiv.org/abs/1011.2368

A. Arda, R. Sever, and C. Tezcan, “Approximate analytical solutions of the effective mass Dirac equation for the generalized Hulthén potential with any κ-value,” Central European Journal of Physics, vol. 8, no. 5, pp. 843–849, 2010.

J. Meng, K. Sugawara-Tanabe, S. Yamaji, P. Ring, and A. Arima, “Pseudospin symmetry in relativistic mean field theory,” Physical Review C, vol. 58, no. 2, pp. R628–R631, 1998.

J. Meng, K. Sugawara-Tanaha, S. Yamaji, and A. Arima, “Pseudospin symmetry in Zr and Sn isotopes from the proton drip line to the neutron drip line,” Physical Review C, vol. 59, pp. 154–163, 1999.

M. R. Spiegel, Schaum's Outline of Theory and Problems of Laplace Transforms, Schaum Publishing, New York, NY, USA, 1965.

F. Doetsch, Guides to the Application of Laplace Transforms, Princeton University press, Princeton, NJ, USA, 1961.

M. Abramowitz and I. A. Stegun, Handbook Mathematical Functions With Formulas, Graphs, and Mathematical Tables, Dover, New York, NY, USA, 1970.