Int J Theor Phys (2011) 50:2811–2818 DOI 10.1007/s10773-011-0780-4
Dirac Equation and Ground State of Solvable Potentials: Supersymmetry Method H. Panahi · Z. Bakhshi
Received: 22 December 2010 / Accepted: 30 March 2011 / Published online: 15 April 2011 © Springer Science+Business Media, LLC 2011
Abstract The supersymmetry in non-relativistic quantum mechanics is applied as an algebraic method to obtain the solutions of the Dirac equation with spherical symmetry electromagnetic potentials. We show that some of the superpotentials related to ground state of the solvable potentials in non-relativistic quantum mechanics can be used for studying of the Dirac equation. Keywords Dirac equation · Solvable potentials · Supersymmetric quantum mechanics 1 Introduction Solutions of the Dirac equation with physical potential are very useful to investigate the relativistic effects [1]. The Dirac-Coulomb problem is an exactly solvable problem in relativistic quantum mechanics and its solution can be found in text books on quantum mechanics [2, 3]. The Dirac-Coulomb problem has also been studied via supersymmetry quantum mechanics (SUSY QM) [4–11]. In 1989, the relativistic Dirac-Oscillator potential was introduced by adding an off-diagonal linear radial term to the Dirac operator. The relativistic bound states spectrum and its eigen states were also obtained explicitly [12–16]. The DiracMörse problem was solved successfully by Alhaidari [17, 18]. After that, a lot of shapeinvariance potentials such as Dirac-Rosen-Mörse, Dirac-Eckart, Dirac-Pöschel-Teller and Dirac-Scarf potentials have been solved exactly [19, 20]. An effective approach for solving the Dirac equation for a charged spinor in spherically symmetric four component electromagnetic potential was introduced by Alhaidari [17–20]. In this method the gauge invariance and spherically symmetry of the electrostatic potential are used to arrive at the radial Dirac equation. In fact, an unitary transformation can be applied to Dirac equation such that the resulting second order differential equation becomes H. Panahi () · Z. Bakhshi Department of Physics, University of Guilan, Rasht 51335-1914, Iran e-mail:
[email protected] Z. Bakhshi e-mail:
[email protected]
2812
Int J Theor Phys (2011) 50:2811–2818
Schrödinger-like equation. There are a lot of methods for solving the Schrödinger-like equation, such as SUSY QM [11], shape invariance [21], point canonical transformation (PCT) [22] and etc. In this work, in Sect. 2, we give a brief introduction for the approach which transform the Dirac equation with spherically symmetric potentials to two Schrödinger equations for upper and lower spinor fields. We factorize these obtained Hamiltonians in terms of two incretion and annihilation differential operators and then try to solve them by SUSY QM method. According to concepts of SUSY in non-relativistic quantum mechanics, the solvable potentials can be stated in terms of superpotentials by a Riccati equation and so we try to relate the well-known superpotentials to our obtained potentials in Dirac equation. It is seen that by inducing one assumption over electrostatic potential and gauge field of Dirac equation, a large class of solvable potentials can be used for solving Dirac equation. We show that the relativistic energy of all of these potentials are constant and the spinor fields of them are obtained by ground state of the solvable potentials. We complete our calculation for two examples in Sect. 3 and we give a table for other potentials. In Sect. 4, the paper ends with a brief conclusion.
2 Dirac Equation and Supersymmetry Method In atomic units (m = e = = 1) and taking the speed of light c = α −1 , the Hamiltonian for a can be written as [2]: Dirac spinor in the four-component electromagnetic potential (A0 , A) 1 + αA0 −iα σ .∇ + α σ .A H= , (2.1) −iα σ .∇ + α σ .A −1 + αA0 where α is fine structure constant and σ are the three 2 × 2 Pauli matrices [17]. Taking into consideration gauge invariance, the form of electromagnetic potential for static charge distribution with spherical symmetry is: = (αυ(r), rˆ ω(r)), (A0 , A)
(2.2)
where rˆ is radial unit vector, υ(r) is electrostatic potential and ω(r) is a gauge field does not contribute to the magnetic field. By substituting the two off-diagonal term α σ .A by ±iα σ .A in (2.1), the Hamiltonian leads to the following two component radial Dirac equation [2]: α( kr + ω(r)) − drd 1 + α 2 υ(r) G(r) G(r) =ε , (2.3) d F (r) F (r) α( kr + ω(r)) + dr −1 + α 2 υ(r) where ε is the relativistic energy and k is the spin-orbit coupling parameter defined as k = ±(j + 12 ) for l = j ± 12 . Equation (2.3) gives two coupled first order differential equations for the radial spinor components. Separating differential term in (2.3); we will have: k d 1−ε 0 0 −( r + ω(r)) −αυ(r) F (r) F (r) F (r) dr α + = . d 1+ε G(r) G(r) G(r) αυ(r) ( kr + ω(r)) 0 dr 0 α (2.4) Now, let D be the operator which diagonalizes the matrix that is appeared in the interaction term as [17–20]: cos ρ(r) sin ρ(r) D= (2.5) − sin ρ(r) cos ρ(r)
Int J Theor Phys (2011) 50:2811–2818
2813
and multiplying it in the Dirac equation (2.4), with assumption ρ(r) = constant, we can obtain: k S ˜ C k ε d + + ω(r) C + − αυ(r) + − + ω(r) S F˜ (r), (2.6) G(r) = dr r α α α r k S ˜ C k d ε ˜ + ω(r) C + − αυ(r) − − + ω(r) S G(r), − + F (r) = dr r α α α r (2.7) where S = sin(2ρ) and C = cos(2ρ) and: F˜ (r) F (r) =D . (2.8) ˜ G(r) G(r) By defining the intertwining operators as: d + U (r), dr d A+ = − + U (r), dr A=
where
U (r) =
S k + ω(r) C + , r α
(2.9) (2.10)
(2.11)
then we can write (2.6) and (2.7) in the following simple form: ˜ AG(r) = (f1 + f2 )F˜ (r),
(2.12)
˜ A F˜ (r) = (f1 − f2 )G(r),
(2.13)
+
where f1 = αε − αυ(r) and f2 = Cα − ( kr + ω(r))S. Multiplying A+ and A in the left of (2.12) and (2.13) respectively, we have: ˜ (2.14) = A+ (f1 + f2 ) F˜ (r) + (f1 + f2 )A+ F˜ (r), A+ AG(r) ˜ ˜ + (f1 − f2 )AG(r). AA+ F˜ (r) = A(f1 − f2 ) G(r)
(2.15)
A+ (f1 + f2 ) = 0,
(2.16)
A(f1 − f2 ) = 0,
(2.17)
Now by assuming
then after some calculation, we obtain f12 − f22 = constant ≡ η,
(2.18)
that is, the functions f1 and f2 are related to each other by a constant parameter. Hence, as f1 and f2 are written in terms of υ(r), ω(r) and ε, then, we can obtain a relation between the electrostatic potential υ(r), the gauge field ω(r) and the relativistic energy ε for the radial Dirac equation (2.1).
2814
Int J Theor Phys (2011) 50:2811–2818
Also, it is easy to see that, the following eigenvalue equations can be obtained from the above equations ˜ ˜ ˜ A+ AG(r) = f12 − f22 G(r) = ηG(r), (2.19) AA+ F˜ (r) = f12 − f22 F˜ (r) = ηF˜ (r),
(2.20)
that is, by inserting of the constant parameter η, we have obtained two isospectral eigenvalue equations for the upper and lower spinor fields. Also, according to (2.9) and (2.10), we see that the Hamiltonians H1 = A+ A and H2 = AA+ have the following potentials V1 (r) = U 2 (r) −
dU (r) , dr
(2.21)
dU . dr
(2.22)
V2 (r) = U 2 (r) +
For instance, for upper spinor field, (2.19) leads to following Schrodinger equation −
d 2 G(r) k 2 2 2CS k S2 + ω(r) + C + ω(r) + + dr 2 r α r α2 k dω(r)
= ηG(r).
− 2 C G(r) − dr r
(2.23)
On the other hand, from factorizing of H1 and H2 in terms of intertwining operators A and A+ , it is seem that, we can obtain the solution of the radial Dirac equation from the concept of supersymmetry in quantum mechanics. In fact, according to Ref. [11], the Hamiltonians H1 and H2 , except for the ground state, are the partner supersymmetry each other and so the spectrum of non-zero energies are the same and except to the ground state, the corresponding eigenfunctions can obtain form each other by acting the ladder or rasing operator on each other. In other word, except for the ground state, the isospectral states of the supersymmetry partner Hamiltonians can be obtained from each other by an algebraic method [11]. Hence, if we could find the electrostatic potential υ(r) and the gauge field ω(r) for radial Dirac equation (2.1), such that η in (2.18) became a constant, then the upper and lower spinor fields of (2.1) can be solved by supersymmetry method in non-relativistic quantum mechanics. Unfortunately, we could not find such functions and it is might only for η = 0, that is, only for ground state of the Hamiltonians (2.19) and (2.20). With this condition, these Hamiltonians are not the supersymmetry partner of each other and so, there eigenfunctions are not related to each other in concept of supersymmetric quantum mechanics. On the other hand, it is easy to see that for η = 0, we will have k S, υ(r) = ± ω(r) + r
(2.24)
ε = ±C,
(2.25)
that is, the electrostatic potential υ(r) and the gauge field ω(r) are related to each other by (2.24) and also the relativistic energy ε becomes a constant, the constant C which obtained from the unitary transformation (2.5). Therefore the eigenfunctions of system or the upper
Int J Theor Phys (2011) 50:2811–2818
2815
and lower spinor fields can obtain from the ground state of (2.19) and (2.20). In other words, they should be solved from (2.19) and (2.20) by inserting η = 0: −
˜ d 2 G(r) ˜ = 0, + V1 (r)G(r) dr 2
(2.26)
˜ d 2 F (r) ˜ = 0. + V2 (r)F (r) (2.27) dr 2 Since each of the potentials V1 (r) and V2 (r), associated to upper and lower spinor fields, are related to superpotential U (r), by Riccati equation (2.11), then the ground state with zero energy all of the solvable potentials in non-relativistic quantum mechanics can be used in our work. This means that, the Dirac equation for those spherically symmetry electrostatic potentials and gauge fields, which are related to each other by (2.24), have constant relativistic energy and the spinor fields are obtained from the ground state of potentials V1 (r) and V2 (r). In this work, first we consider (2.26) for upper spinor field and try to compare all of the well-known solvable potentials, given in Ref. [22], to potential V1 (r) which according to (2.23) and (2.24) could be written in terms of υ(r) or ω(r). Then, for each superpotential given in Ref. [22], we try to obtain υ(r) and ω(r) for radial Dirac equation (2.1) and so the ground state wave function of solvable potential will be as the upper spinor field of the Dirac equation with relativistic energy ε = ±C. In other word, the ground spinor field of the radial Dirac equation (2.1) can be obtained from the ground state of the solvable potentials, given in Ref. [22], which the relativistic energy all of them are constant. The lower spinor field is also obtained from the upper by (2.12). We can also consider the eigenvalue equation (2.27) for lower spinor and do similar formalism for obtaining υ(r) and ω(r) with constant relativistic energy for radial Dirac equation (2.1). In this case, the upper spinor field is obtained from the lower by (2.13). In the next section, we calculate the above formalism for two examples, The DiracShifted Oscillator and The Dirac-Coulomb, and end we give a table for a large class of electrostatic potentials and gauge fields which can be obtained from solvable potentials. −
3 Examples 3.1 The Dirac-Shifted Oscillator First, we consider the following superpotential according to Ref. [22] 1 U (r) = ωr − b, 2
(3.1)
where ω and b are the constant value. Corresponding to (3.1) and by using of (2.11) and (2.24), the electrostatic potential υ(r) and the gauge field ω(r) for radial Dirac equation (2.1) are calculated as: k ω (3.2) r− , ω(r) = ± √ 2 2 r 2 1−α b αωb υ(r) = ∓ √ r, 2 1 − α 2 b2
(3.3)
2816
Int J Theor Phys (2011) 50:2811–2818
where S b=− . α
(3.4)
Also, according to (2.25), the relativistic energy for the above electromagnetic potentials is obtained as: 1
ε = ±(1 − α 2 b2 ) 2 .
(3.5)
The upper spinor field is also obtained from the ground state of the given superpotential (3.1), which according to Ref. [22], is given as:
= exp − 1 g 2 , G(r) 2
(3.6)
12 2b ω r− g(r) = . 2 a
(3.7)
where
By using the upper field, we can calculate the lower spinor field from (2.12) too. 3.2 The Dirac-Coulomb If we consider the superpotential as [22]: U (r) =
(l + 1) e2 − , 2(l + 1) r
(3.8)
then, after some calculation we will have the following gauge field ω(r) and the electrostatic potential υ(r) for radial Dirac equation from (2.11) and (2.24) ω(r) = ∓
2(l + 1)2 1 2 4 2
[4(l + 1)2 − α e ]
υ(r) = ±
1 k − , r r
e2 1 2 4 2
[4(l + 1)2 − α e ]
l+1 , r
(3.9)
(3.10)
where S e2 = . α 2(l + 1)
(3.11)
We can also obtain the relativistic energy as: ε=± 1−
α 2 e4 4(l + 1)2
12 ,
(3.12)
and for the upper spinor, we have:
= g l+1 (r) exp − g(r) , G(r) 2
(3.13)
Int J Theor Phys (2011) 50:2811–2818
2817
Table 1 Upper spinor of the ground state and relativistic energy of the radial Dirac equation for electrostatic potential and gauge field, associated with superpotential and internal function. The lower spinor can be calculated from (2.12) Type of Dirac-potentials
Shifted oscillator
Superpotential U (r)
Gauge field ω(r)
Ground state G(r)
Internal function g(r)
Electrostatic potential υ(r)
Relativistic energy ε
1 ωr − b 2
±
ω k 1 r− r 2(1−α 2 b2 ) 2 bω ∓[ 1 ]r 2(1−α 2 b2 ) 2
exp(− 21 g 2 (r))
1 2b 2 g(r) = ( ω 2 ) (r − a ) ω = constant Three- dimensional oscillator
Coulomb
1 ωr − l+1 r 2 g(r) = 12 ωr 2 ω = constant e2 − l+1 r 2(l+1) e2 )r g(r) = ( n+l+1
±( ω 2 )r k = −(l + 1) 0
1 ε = ±(1 − α 2 b2 ) 2 S = −αb (l+1) g(r) g 2 (r) exp(− 2 ) ε = ±1 S=0
2(l+1)2 k 1 1 )r − r [4(l+1)2 −α 2 e4 ] 2 e2 l+1 ±( 1 ) r [4(l+1)2 −α 2 e4 ] 2 ∓(
g(r) g l+1 (r) exp(− 2 ) 2 4 1 ε = ±(1 − α e 2 ) 2 4(l+1) αe2 2(l+1) g(r) s−n g (r) exp(− 2 ) S=
Mörse
B
A − B exp(−ar)
k ∓ 1 exp(−ar) − r (1−α 2 A2 ) 2
g(r) = 2B a exp(−ar)
±
a = constant Scarf-II (hyperbolic)
Rosen-Mörse-II(hyperbolic)
B cosh(ar) g(r) = i sinh(ar) a = constant
0
A tanh(ar) + B A
±
A tanh(ar) +
±
a = constant Eckart
B − A coth(ar) A B > A2 g(r) = coth(ar)
Generalized Pöchl-Teller
Rosen-Mörse-I (trigonometric)
A cot(ar) − B A
±
−A cot(ar) + B csc(ar) π −π 2 ≤ ar ≤ 2 g(r) = cos(ar) a = constant
AB 1 tanh(ar) (A2 −α 2 B 2 ) 2
AB 1 coth(ar) (A2 −α 2 B 2 ) 2
1 ε = ±(1 − α 2 A2 ) 2 S = αA −1 (1 − g 2 (r)) 2 × exp(−λ tan−1 (−ig(r))) ε = ±1 S=0 (1 − g(r))
S = αB A (g(r) − 1)
S = αB A
0
ε = ±1 S=0
±
A2 k 1 cot(ar) − r 2 (A −α 2 B 2 ) 2 AB 1 cot(ar) (A2 −α 2 B 2 ) 2
−(s+n+a) −(s+n−a) 2 2 × (g(r) + 1)
2 2 1 )2 ε = ±(1 − α B A2
(g(r) − 1)
±
(s−n−a) (s−n+a) 2 2 × (1 + g(r))
2 2 1 )2 ε = ±(1 − α B A2
B ±(A coth(ar) − ) − kr sinh(ar)
a = constant Scarf-I (trigonometric)
k A2 1 tanh(ar) − r 2 (A −α 2 B 2 ) 2
k ∓ 1 coth(ar) − r (A2 −α 2 B 2 ) 2
0 ≤ ar ≤ π g(r) = −i cot(ar)
B ) − kr cosh(ar)
A2
a = constant B A coth(ar) − sinh(ar) A