Journal of Modern Physics Vol.06 No.03(2015), Article ID:54289,23
pages
10.4236/jmp.2015.63034
A Modified Method for Deriving Self-Conjugate Dirac Hamiltonians in Arbitrary Gravitational Fields and Its Application to Centrally and Axially Symmetric Gravitational Fields
M. V. Gorbatenko, V. P. Neznamov*
Russian Federal Nuclear Center, All-Russian Research Institute of Experimental Physics, Sarov, Russia
Email: *neznamov@vniief.ru
Copyright © 2015 by authors and Scientific Research Publishing Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
Received 6 February 2015; accepted 24 February 2015; published 27 February 2015
ABSTRACT
We have proposed previously a method for constructing self-conjugate Hamiltonians Hh in the h- representation with a flat scalar product to describe the dynamics of Dirac particles in arbitrary gravitational fields. In this paper, we prove that, for block-diagonal metrics, the Hamiltonians Hh can be obtained, in particular, using “reduced” parts of Dirac Hamiltonians, i.e. expressions for Dirac Hamiltonians derived using tetrad vectors in the Schwinger gauge without or with a few summands with bispinor connectivities. Based on these results, we propose a modified method for constructing Hamiltonians in the h-representation with a significantly smaller amount of required calculations. Using this method, here we for the first time find self-conjugate Hamiltonians for a number of metrics, including the Kerr metric in the Boyer-Lindquist coordinates, the Eddington- Finkelstein, Finkelstein-Lemaitre, Kruskal, Clifford torus metrics and for non-stationary metrics of open and spatially flat Friedmann models.
Keywords:
Self-Conjugate Hamiltonian, Dirac Particle, Arbitrary Gravitational Field, Schwinger Gauge, Kerr Metric
1. Introduction
In [1] , we proposed a method for constructing self-conjugate Hamiltonians in the -representation with a flat scalar product to describe the dynamics of Dirac particles in arbitrary gravitational fields.
Using the algorithm proposed in [1] , we calculated Hamiltonians in the -representation for the Schwarzschild and Friedmann-Robertson-Walker cosmological model metrics. However, application of the algorithm to the Kerr metric necessitated a large amount of calculations to find Christoffel symbols, bispinor connectivities etc., and cumbersome algebraic transformations of arising expressions.
We made attempts to simplify the algorithm [1] . First, we proved the theorem, according to which a Hamiltonian in the -representation for an arbitrary gravitational field, including a time-dependent one, is a Hermitian part of the initial Dirac Hamiltonian derived using tetrad vectors in the Schwinger gauge1.
(1)
Then, for block-diagonal metrics, using Equation (1), we proved the second theorem, according to which the Hamiltonians and in Equation (1) can be replaced by their “reduced” parts without or with a few summands with bispinor connectivities:
(2)
Block-diagonal metrics are understood to be metric tensors of the form
(3)
Apparently, the cases belong to the same kind as (3) when, and also when or are used instead of.
In Equation (2), is part of the initial Dirac Hamiltonian, which contains only the mass term and terms with momentum operator components (i.e. with coordinate derivatives).
The summand in (2) equals
(4)
One can see that is a fairly simple expression, which in some cases differs from zero for the block-dia- gonal metrics with. For example, the Kerr metric in Boyer-Lindquist coordinates [2] , [3] belongs to such case. Of course, application of Equation (2) makes the procedure of deriving self-conjugate Hamiltonians in the -representation much less complicated.
Equations (1) and (2) are proven in Sections 3, 4 of this paper.
In the second part of the paper, we use (1) and (2) to find self-conjugate Hamiltonians for the Kerr [2] , [3] , Eddington-Finkelstein [4] [5] , Painlevè-Gullstrand [6] [7] , Finkelstein-Lemaitre [5] , Kruskal [8] [9] , Clifford torus [10] metrics, and for non-stationary metrics of open and spatially flat Friedmann models. For all these metrics, except for the Painlevè-Gullstrand one, self-conjugate Hamiltonians are derived for the first time2.
At the end of the paper, we prove that self-conjugate Dirac Hamiltonians in a weak Kerr field are physically equivalent in both harmonic Cartesian and Boyer-Lindquist coordinates.
In the Conclusions, we discuss the outcome of this study and the results of applying the developed algorithm to the evolution of bound atomic and quark states in the expanding universe.
2. Reducing the Dirac Equation to the Schrödinger Form. An Algorithm for Finding a Self-Conjugate Hamiltonian in the -Representation
Let us recall the line of corresponding reasoning and introduce the notation. Tetrad vectors are defined by the relation
(5)
where
(6)
In addition to the system of tetrad vectors, one can introduce three other systems of tetrad vectors, , , , which differ from in the location of the global and local (underlined) indices. The global indices are raised up and lowered by means of the metric tensor and inverse tensor, and the local indices, by means of the tensors,.
We assume that the quantum mechanical motion of particles is described by the Dirac equation, which is written in the units of as
(7)
Here, is the particle mass, is a four-component “column” bispinor, and are Dirac ma- trices, satisfying the relation
(8)
in (8) means a unity matrix. The round parentheses in (7) contain a covariant bispinor derivative,:
(9)
Equation (9) for contains the bispinor connectivity, for finding which one should fix some system of tetrad vectors defined as (5). Upon that, the quantity can be expressed through “Christoffel” vector derivatives in the following way (the “Christoffel” derivatives are denoted by a semicolon):
(10)
The expression for in (10) is defined below, see (14). The bispinor connectivity given by (10) provides invariance of the covariant derivative with respect to the transition from one system of tetrad vectors to another.
In what follows, along with Dirac matrices with global indices, we will use Dirac matrices with local indices. The relation between and is given by the expression
(11)
It follows from (11), (8), (5) that
(12)
In terms of the matrices, the Dirac Equation (7) can be written as follows:
(13)
It is convenient (but not necessary) to choose the quantities so that they have the same form for all local frames of reference. Both systems and can be used to construct a full system of matrices. The full system is, for example, the system
(14)
Any system of Dirac matrices provides for several discrete automorphisms. We restrict ourselves to the automorphism
(15)
The matrix will be called anti-Hermitizing.
It follows from (7) that the initial Hamiltonian is given by the following expression:
(16)
The operator (16) has a meaning of the evolution operator of the wave function of a Dirac particle in the chosen global frame of reference.
Ref. [1] formulates the rules of finding a Hamiltonian in the -representation for a Dirac particle in an arbitrary gravitational field. A-priori information, which is assumed to be known, is information about the metric
tensor, Christoffel symbols, local metric tensor and local Dirac matrices. These
rules are the following:
1) For a gravitational field described by the metric, we find a system of tetrads, satisfying the Schwinger gauge. Recall that components of the tetrads, in this gauge correlate with components of the tensor as follows:
(17)
components are identically zero
(18)
In order to find, we introduce a tensor, , with components
(19)
The tensor satisfies the condition
(20)
As we can use any triplet of three-dimensional vectors, satisfying the relation
(21)
In what follows, the quantities, dependent on the choice of tetrad vectors, are denoted by a tilde, if they are calculated in the system of tetrad vectors in the Schwinger gauge.
2) In accordance with (16), we write a general expression for the Hamiltonian.
(22)
Here
(23)
(24)
3) The expression for the Hamiltonian equals
(25)
where the operator is defined by the relation
(26)
As distinct from [1] , in Equation (26) we use only the gravitational part of the determinant, which is there due to the presence of an external gravitational field. An additional multiplier arises, if we use curvilinear coordinates in accordance with the equality of scalar products for wave functions in the initial and -representations [1] :
(27)
Hence,
(28)
where is the determinant, which arises when the volume element is written in curvilinear coordinates. Given that the conditions of coordinate harmonicity [12] [13] are satisfied, for Cartesian coordinates, for cylindrical coordinates, is for spherical coordinates etc.
Equations (25), (26) define the operator, which is the sought Hermitian Hamiltonian in the representation.
Thus,
(29)
Note that procedures for constructing self-conjugate Hamiltonians with a flat scalar product that could be used for studying the dynamics of spin particles in gravitational fields of particular form have been proposed in literature more than once [14] -[17] . These attempts are not general, but they produce correct results as applied to the choice of particular metrics and tetrad vectors.
3. Proving the Equality for an Arbitrary, including Time-Dependent, Gravitational Field
We start the proof from transforming the right side in (25).
(30)
After inserting (22), (26) into (25) we have:
(31)
The next step in the proof is to find an expression for. For this purpose, we employ relations (77) from [1] :
(32)
(33)
(34)
We insert (33), (34) into (32):
(35)
Using (22) and (35), we calculate the quantity:
(36)
By comparing (36) with (31) we conclude that Equation (1) is valid.
(37)
4. Proving the Equality for Gravitational Fields with the Block-Diagonal Metrics
An expression for the “reduced” Hamiltonian is derived from (22) by deleting the terms with bispinor connectivities. Thus,
(38)
This expression can also be written as
(39)
Taking the Hermitian conjugation from (39), we obtain
(40)
It follows from (39), (40) that
(41)
Considering Equation (37), we obtain
(42)
Let us introduce the following notation:
(43)
(44)
Some transformations give the following expressions for
(45)
and for:
(46)
where
(47)
(48)
(49)
(50)
In (50), is a totally antisymmetric third-rank tensor.
Then, we calculate and using the relations (3), (5), (17), (18), and diagonal representation of. Direct calculations show that
(51)
(52)
Thus, it turns out that for the block-diagonal metrics of the form (3) we can find the Hamiltonian using the fairly simple formula (2).
5. Centrally Symmetric Gravitational Field
This section presents Hamiltonians in the -representation for Dirac particles in centrally symmetric gravitational fields, when the metrics are written in various coordinates.
5.1. The Schwarzschild Metric
Writing the Schwarzschild solution in the coordinates
(53)
gives:
(54)
In Equation (54), is the gravitational radius.
The resulting expression for derived in [1] and revised to include (26) is
(55)
In Equation (55),.
It is easy to verify that Equation (55) can be found in a comparatively straightforward manner using formula (2), if we take into account that
(56)
In Refs. [1] [14] , the authors also derived a Hamiltonian for the Schwarzschield metric in isotropic coor- dinates
(57)
(58)
The expression for can be easily derived from (2) using
(59)
5.2. Eddington-Finkelstein Metric
The Eddington-Finkelstein solution ([4] [5] ) in the coordinates
(60)
is given by
(61)
(62)
The inverse tensor has the following form:
(63)
Calculations of a “reduced” Hamiltonian using (38) gives
(64)
The Hamiltonian in the -representation is calculated using (2) given that for the metric of interest. in Table 1. We obtain:
Table 1.―like tetrad vectors.
(65)
5.3. Painlevé-Gullstrand Metric
In this section, we find a self-conjugate Hamiltonian for a Dirac particle in a spherically symmetric gravitational field described by the Painlevé-Gullstrand metric. The Hamiltonian for this metric is calculated first using the algorithm of [1] and then using (1) and (2).
In the coordinates, the Painlevé-Gullstrand metric [6] has the following form:
(66)
The determinants equal
(67)
Tetrad vectors in the Schwinger gauge:
(68)
(69)
Christoffel symbols:
(70)
(71)
(72)
(73)
Bispinor connectivities are calculated by the formula (24) using (66)-(73). We obtain:
(74)
In order to find the Hamiltonian in the -representation, , expressions for and expressions (74) for are put into the primary Dirac Hamiltonian,
(75)
These transformations give
(76)
The operator for the Painlevè-Gullstrand metric equals
(77)
so the -representation coincides with the representation of the Hamiltonian.
The Hamiltonian (76) is self-conjugate. It is evident that the formula (1) is valid in this case.
It is easy to obtain (76) from (2) given that, and is written as
(78)
Thus, as applied to the Painlevé-Gullstrand metric, the same Hamiltonian was obtained both by the standard algorithm and in a simpler manner using (2).
In [11] , a self-conjugate Hamiltonian was obtained for the Painlevé-Gullstrand metric using tetrad vectors in the Schwinger gauge with a set of local Dirac matrices written in spherical coordinates.
(79)
The set is related to the set through a unitary matrix,
(80)
The Hamiltonian from [11] can be written as
(81)
The Hamiltonians (76) and (81) are physically equivalent, because they are related through a unitary transformation,
(82)
Generally speaking, all Hamiltonians in the Schwinger gauge are connected with each other by physically equivalent matrices of spatial rotation. This is what we meant [1] speaking about the uniqueness of Hamiltonians in the -representation (see the comments by M. Arminjon in [18] ).
5.4. Finkelstein-Lemaitre Metric
It is of independent interest to study the motion of a Dirac particle in the nonstationary Finkelstein-Lemaitre metric [5] , because the time coordinate in this metric coincides with the proper time.
(83)
The determinants equal
(84)
Non-zero components of tetrad vectors in the Schwinger gauge:
(85)
For this metric, in (2),.
“Reduced” Hamiltonian:
(86)
We insert (86) into (2) and obtain
(87)
The Hamiltonian (87) is self-conjugate with a fairly complex time dependence.
5.5. Hamiltonian in the -Representation for Dirac Particles in the Kruskal Gravitational Field
The Kruskal metric [8] is a further development of the Lemaitre-Finkelstein metric to build the most complete frame of reference for a point-mass field. The formula below, in which the frame of reference is synchronous, has been developed by I.D. Novikov [9] . In the coordinates,
(88)
(89)
The determinants equal
(90)
Equations (88), (89) show that the metric (89) is related to the radial coordinate and proper time through the parameter.
The non-zero components of the tetrad vectors in the Schwinger gauge equal
(91)
“Reduced” Hamiltonian:
(92)
According to (2), with, we have
(93)
The derivative in the last summand of (93) should allow for the dependence (see (89)).
6. Axially Symmetric Gravitational Field
6.1. Kerr Metric in the Boyer-Lindquist Coordinates
The Kerr solution in the Boyer-Lindquist coordinates [3]
(94)
is given by
(95)
(96)
The inverse tensor has the following form:
(97)
Here,
(98)
6.2. Tetrad Vectors in the Schwinger Gauge
We will need expressions for tetrad vectors in the Schwinger gauge. The results of calculating the components of tetrad vectors are presented in Table 2. Table 3 shows the components of vectors.
Table 2.―like tetrad vectors.
Table 3.―tetrad vectors.
6.3. Hamiltonian
First, from formula (38), we obtain
(99)
For the metric under consideration, in (2) differs from zero:
(100)
The Hamiltonian is calculated using (2).
(101)
We put the tetrad vector components
and the metric components
(102)
into (101). Finally,
(103)
The quantities are defined by (102), (98).
In order to turn to the Schwarzschild Hamiltonian, one should assume that
(104)
After such a replacement, from (103), we obtain the Hamiltonian (55) for the Schwarzschild field:
(105)
If in the expression for (103) we restrict ourselves to the terms not higher than of the first order of smallness in the parameters, we will obtain a self-conjugate Hamiltonian for the weak Kerr field.
(106)
In Section 7, is found by the general algorithm [1] , and using (1).
Previously, a self-conjugate Hamiltonian for a weak Kerr field has been obtained in Refs. [16] , [19] for the metric written in isotropic coordinates.
7. Weak Axially Symmetric Gravitational Field
7.1. Kerr Metric in the Boyer-Lindquist Coordinates
For our purposes, we write Equations (95)-(98), leaving the terms not exceeding the first order of smallness in
the quantities and. In this approximation,
(107)
(108)
(109)
(110)
7.2. Tetrad Vectors in the Schwinger Gauge
We will need expressions for tetrad vectors in the Schwinger gauge. The results of calculating the components of the tetrad vectors are presented in Table 4. Table 5 shows the components of the vectors.
Table 4. Tetrad vectors.
Table 5. Tetrad vectors.
7.3. Christoffel Symbols
Christoffel symbols:
(111)
7.4. Bispinor Connectivities
Bispinor connectivities are calculated by the formula
(112)
We obtain:
(113)
7.5. Hamiltonian
Taking into account SubSections 7.1 - 7.4, we derive an expression for using (22):
(114)
Since the Kerr solution is stationary, the general formula for,
(115)
in our case will be written as
(116)
where is defined by (110).
As a result, the Hamiltonian can be written as
(117)
Equation (117) coincides with Equation (106), derived by expanding the general expression for (103). In turn, the general expression (103) is obtained using (2).
Analogously, using (114), we can easily check if the formula (1) is valid for the metric under consideration (108).
Thus, the same expression for for a weak Kerr field is in fact derived in three different ways.
For the block-diagonal metrics like (3), as exemplified by the Kerr metric with the formula (2), we can see that the algorithm for finding the Dirac self-conjugate Hamiltonians with a flat scalar product becomes significantly simpler.
8. Open Friedmann Model
Consider the case of the open Friedmann model in the coordinates
For this model, the non-stationary metric takes the following form:
(118)
(119)
The non-zero components of the tetrad vectors in the Schwinger gauge equal
(120)
Calculations of the Hamiltonian give
(121)
The Hamiltonian is defined by (2). For this metric,
(122)
In the quasi-stationary approximation, for the cosmological time, the energy operator for a particle moving in the -direction equals
(123)
Here,
Let us denote
(124)
where; zero subscripts correspond to the current time.
If the radius of the spatial curvature of the universe currently goes to infinity, , then
(125)
In this case, for the spatially flat Friedmann model, the Hamiltonian (122) becomes equal to
(126)
In Cartesian coordinates, the expression for is
(127)
9. Clifford Torus Metric
The metric proposed in [10] in the coordinates is given by
(128)
In (128), the prime denotes the derivative with respect to the coordinate.
Tetrad vectors in the Schwinger gauge:
(129)
In (129),. In accordance with (38), the “reduced” Hamiltonian equals
(130)
For this metric, in (2),. Then, the Hamiltonian in accordance with (2) equals
(131)
10. Equivalence of Hamiltonians with Harmonic Cartesian or Boyer-Lindquist Coordinates in a Weak Kerr Field
As we know, harmonic coordinates satisfy the condition formulated by Th. De-Donder and V. A. Fock [12] [13] .
In Refs. [16] [19] , the following form of self-conjugate Dirac Hamiltonian was derived using harmonic Cartesian coordinates for a weak Kerr field:
(132)
Similar to Subsects 6, 7,; in (132) is angular momentum of a rotating source of the Kerr field.
When using Boyer-Lindquist coordinates [3] , the self-conjugate Hamiltonian for a weak Kerr field is defined by (106), (117):
(133)
In (132), (133), the summands without correspond to the Schwarzschild metric. These parts of (132), (133) are physically equivalent to each other.
The last but one summands in (132), (133) are also physically equivalent to each other. Indeed, in a weak Kerr field, Boyer-Lindquist coordinates are reduced to spherical coordinates.
Hence,
(134)
(135)
Given (135), we can see the desired equivalence.
As for the last summands in (132), (133), they do not seem to be equivalent at first. Suffice it to note that the last summand in (132) contains three spin matrices, while the last summand in (133) has only one spin matrix. This may put the correctness of (133) into question, although it was derived in Subsection 7 in three different ways.
In order to resolve this, let us write the Hamiltonian (133) using the representation of Dirac matrices in spherical coordinates (79). The matrices (79) are related to the matrices by the unitary transformation (80).
The Hamiltonian (133) with the local Dirac matrices takes the following form:
(136)
Given (134), the last summands in (136) and (132) coincide with each other. The Hamiltonian (136) is physically equivalent to the Hamiltonian (133), since it is obtained using the unitary transformwation (80):
(137)
The analysis indirectly proves that Equation (103) is valid for the general Hamiltonian in Boyer-Lindquist coordinates. Equation (103) can be used for Kerr gravitational field of arbitrary strength and angular momentum of the field source of arbitrary magnitude.
The results obtained above demonstrate that for clear physical interpretation of individual summands of Dirac Hamiltonians one should use harmonic Cartesian coordinates. Classical interpretation of individual Hamiltonian terms requires transition to the Foldy-Wouthuysen representation [20] [21] .
11. Conclusions
This study develops the algorithm proposed in [1] for constructing self-conjugate Hamiltonians in the - representation with a flat scalar product to describe the dynamics of Dirac particles in arbitrary gravitational fields. We prove that a Hamiltonian in the h-representation for any gravitational field, including a time-dependent field, is a Hermitian part of the initial Dirac Hamiltonian derived using tetrad vectors in the Schwinger gauge. We also prove that for the block-diagonal matrices like (3), the Hamiltonian can be calculated by the formula (2) using “reduced” parts of the Hamiltonians and without or with a small number of su- mmands with bispinor connectivities. Using this method, we for the first time find self-conjugate Hamiltonians for the Kerr metric in the Boyer-Lindquist form and for the Eddington-Finkelstein, Finkelstein-Lemaitre, Kruskal, Clifford torus metrics and also for non-stationary metrics of open and spatially flat Friedmann models.
In this paper, we also prove physical equivalence of Dirac Hamiltonians in a weak Kerr field in harmonic Cartesian and Boyer-Lindquist coordinates. We point at the necessity of using harmonic Cartesian coordinates for clear physical interpretation of individual terms in the Hamiltonians.
In [22] , the algorithm for deriving self-conjugate Dirac Hamiltonians in the -representation is extended to the electromagnetic case. The Hamiltonian derived is applied to the case when the nonstationary gravitational field describes the spatially flat Friedmann model, and the electromagnetic field is an extension of the Coulomb potential to the case of this model.
Following other authors [17] , we demonstrate that energy levels in atomic systems are invariable in cosmological time.
Spectral lines of atoms in the spatially flat Friedmann model are identical at different points of cosmological time, and redshift is attributed completely to the growth of the wavelength of photons in the expanding universe.
At the same time, we observed that interaction forces and physical dimensions of atomic and quark bound systems vary with universe expansion.
The expressions for Hamiltonians, derived in this paper, can also be employed to study the behavior of Dirac particles in the vicinity of black holes, and scattering and absorption of such particles by black holes.
Acknowledgements
The authors would like to thank Prof. P. Fiziev for the useful discussions, advices and criticism.
References
- Gorbatenko, M.V. and Neznamov, V.P. (2011) Physical Review D, 83, Article ID: 105002. http://dx.doi.org/10.1103/PhysRevD.83.105002
- Kerr, R.P. (1963) Physical Review Letters, 11, 237. http://dx.doi.org/10.1103/PhysRevLett.11.237
- Boyer, R.H. and Lindquist, R.W. (1967) Journal of Mathematical Physics, 8, 265-281. http://dx.doi.org/10.1063/1.1705193
- Eddington, A.S. (1924) Nature, 113, 192. http://dx.doi.org/10.1038/113192a0
- Finkelstein, D. (1958) Physical Review, 110, 965-967. http://dx.doi.org/10.1103/PhysRev.110.965
- Painlevè, P. (1921) C. R. Acad. Sci. (Paris), 173, 677-680.
- Gullstrand, A. (1922) Arkiv för Matematik, Astronomi och Fysik, 16, 1-15.
- Kruskal, M. (1960) Physical Review, 119, 1743-1745. http://dx.doi.org/10.1103/PhysRev.119.1743
- Novikov, I.D. (1963) AJ, 40, 772.
- Fiziev, P.P. arxiv:1012.3520.
- Lasenby, A., Doran, C., Pritchard, J., Caceres, A. and Dolan, S. (2005) Physical Review D, 72, Article ID: 105014. http://dx.doi.org/10.1103/PhysRevD.72.105014
- De-Donder, Th. (1921) La gravifique Einshtenienne. Gauthier-Villars, Paris.
- Fock, V.A. (1939) Journal of Physics, 1, 81-116.
- Obukhov, Y.N. (2001) Physical Review Letters, 86, 192. Obukhov, Y.N. (2002) Fortschritte der Physik, 50, 711-716. http://dx.doi.org/10.1002/1521-3978(200205)50:5/7<711::AID-PROP711>3.0.CO;2-Z
- Leclerc, M. (2006) Classical and Quantum Gravity, 23, 4013-4019. http://dx.doi.org/10.1088/0264-9381/23/12/001
- Obukhov, Y.N., Silenko, A.J. and Teryaev, O.V. (2009) Physical Review D, 80, Article ID: 064044. http://dx.doi.org/10.1103/PhysRevD.80.064044
- Andretsch, J. and Schäfer, G. (1978) General Relativity and Gravitation, 9, 243, 489.
- Arminjon, M. arxiv:1107.4556v2 [gr-qc].
- Gorbatenko, M.V. and Neznamov, V.P. (2010) Physical Review D, 82, Article ID: 104056. http://dx.doi.org/10.1103/PhysRevD.82.104056
- Foldy, L.L. and Wouthuysen, S.A. (1950) Physical Review, 78, 29-36. http://dx.doi.org/10.1103/PhysRev.78.29
- Neznamov, V.P. and Silenko, A.J. (2009) Journal of Mathematical Physics, 50, Article ID: 122302. http://dx.doi.org/10.1063/1.3268592
- Gorbatenko, M.V. and Neznamov, V.P. arxiv: 1105.4709 [gr-qc].
NOTES
*Corresponding author.
1We use the same notations as in [1] .
2For the Painlevè-Gullstrand metric, a physically equivalent self-conjugate Hamiltonian has been derived and studied earlier in [11] .