Motion of charged particle in Reissner-Nordstrom spacetime: A Jacobi metric approach

The present work discusses motion of neutral and charged particles in Reissner - Nordstrom spacetime. The constant energy paths are derived in a variational principle framework using the Jacobi metric which is parameterized by conserved particle energy. Of particular interest is the case of particle charge and Reissner-Nordstrom black hole charge being of same sign since this leads to a clash of opposing forces - gravitational (attractive) and Coulomb (repulsive). Our paper aims to compliment the recent works of Pugliese, Quevedo and Ruffini [1,2]. The energy dependent Gaussian curvature (induced by Jacobi metric), plays an important role in classifying the trajectories.


Introduction
Point particle dynamics serves the important purpose of initiating the study more complex motion of extended bodies. A major area of interest, especially in recent times, is the motion of particles in the presence of black holes and naked singularities. The exhaustive topical works by Pugliese, Quevedo and Ruffini [1,2] (as well as their series of earlier papers [3]) have revealed clearly the possibility of distinguishing between black holes and naked singularities in the case of charged particle dynamics in Reissner -Nordström background. The systematic analysis is based on the nature of particle motion for different charge to mass ratio in the Reissner -Nordström metric, (which is expected), but the non-trivial effect of the charge to mass ratio of the probe particle is indeed a surprising element. In this perspective, the aim of the present paper, dealing once again on the motion of point charge in the presence of a charged black hole, is not to add any more detail on the results of [1,2,4], but to compliment their observations and bring in new insights from an entirely different, (and hitherto less explored), point of view, that of Jacobi metric approach. This framework has generated a lot of interest after the series of works by Gibbons [5] and by Chanda, Gibbons and Guha [6].
The primary major difference is the following. On the one hand, Pugliese et.al. [1,2,3] relies on the covariant framework where geodesic motion takes place in Reissner-Nordström spacetime. The Lagrangian equation of motion is derived from the unrestricted variational principle with the charge and mass parameters of the source and probe dictating the particle motion for arbitrary particle energy, which is a derived quantity. On the other hand, in Jacobi metric framework, the equation of motion is obtained from a restricted variational principle with fixed energy where the Lagrangian depends on the Jacobi metric that 2 Maupertuis principle and Jacobi metric The Maupertuis transformation and the associated Jacobi metric approach has generated a lot of recent interest [5,6]. To understand the context let us start by considering an obvious fact that a curve can be parameterized in infinite number of ways (see for example [8] for a detailed discussion). For example, an ellipse x 2 a 2 + y 2 b 2 = 1 (2) can be parameterized by x = a sin(t), y = b cos(t), and also by a different parametert related to t bỹ t = (a 2 + b 2 )t + (a 2 − b 2 )sin(t).
Clearly the former parameterization (3) is connected to a two dimensional oscillator whereas the latter (4) is related to the Kepler problem. Generically the same trajectory can be attributed to distinct integrable systems where the parameter t plays the role of time variable, conjugate to the respective Hamiltonian. Hence the question is how to derive the particular parameterization that matches with a known or interesting problem. The Maupertuis principle recasts this problem in a dynamical setup. The Maupertuis variational principle states that, given an n-dimensional configuration space Lagrangian L(q i ,q i ), extremals of the action coincide with extremals of the reduced action in a 2n − 1-dimensional phase space that is a level set of the Hamiltonian function H(q i , p i ) = E. Here p i constitute the conjugate momenta to q i and the variation takes place on a fixed energy (E) hypersurface. In a more explicit way let us suppose a natural Hamiltonian in a Riemanniam manifold with metric g ij , The trajectories for H on a fixed energy smooth submanifold will coincide with trajectories of a new HamiltonianH given byH along with a transformation of the parameters, Together (7,8) constitute the Jacobi transformation. A scaling of the effective metric provides the Jacobi metric. The two major novelties and advantages of the Maupertuis-Jacobi framework are: (i) The HamiltonianH yields geodesic motion. This allows one to treat a dynamical problem as geodesic motion and well known machinery of geodesic motion can be directly exploited.
(ii) The Maupertuis principle and subsequent Jacobi transformation preserve integrability that is under this map an integrable system with a natural Hamiltonian goes over to another integrable system in the same phase space. Thus this scheme provides a method to search for new integrable systems.
(iii) The Gaussian curvature induced by Jacobi metric allows one to classify the trajectories based on their energy since it appears explicitly in Jacobi metric. For a generic Lagrangian of the form, (following the notation of [5]) it was shown by Jacobi that the constrained motion of a particle with energy E is provided by geodesics of the rescaled metric It is interesting to observe that particle interactions can induce a curvature in the Jacobi metric through the potential function in an otherwise flat Newtonian space. One of the early workers in this topic was Pin [9] who considered many body systems and in particular showed that the Gaussian curvature of the Jacobi metric has opposite sign to the particle energy E (non-relativistic, without the rest energy). In later times Gibbons and coworkers [10,11] have considered the optical metric in various physical situations which is a closely related concept for massless particles. Discussions on Jacobi metric approach in modern perspective can be found in [12]. Again very recently the Jacobi metric formalism in relativistic scenario has been applied by Gibbons [5] in an elegant study of massive particle motion in Schwarzschild spacetime. In the present paper we closely follow and extend this work to massive particle motion in Reissner -Nordström background. Indeed, Jacobi metric for a more general metric, that is the Kerr-Newman, has been derived in [6] although the probe particle dynamics was not considered there. We consider both cases of the probe particle being neutral and charged. The results show a qualitative difference between the two cases since in the latter one needs to consider the additional Coulomb interaction term between the source and probe particle. For neutral particle the correction terms depend on Q 2 , Q being the charge of the black hole but interestingly for the charged probe correction terms involve qQ terms as well, q being charge of the particle, showing that the relative sign between the particle and black hole charge becomes important.
For the generic metric the action for a massive particle in this background can be written as, The canonical momentum leads to the Hamiltonian This provides the Hamilton-Jacobi equation for the geodesics, parameterized by the energy E, where p i = ∂ i S . Finally, the Hamiltonian-Jacobi equation for geodesics of Jacobi-metric j ij is given by, where j ij is defined as Infact f ij = V −2 g ij turns out to be the optical or Fermat metric. For massless particles (m = 0) the Jacobi metric becomes equal to the Fermat metric modulo a factor of E 2 and subsequently the geodesics do not depend upon energy E. However in the massive case, m = 0, the geodesics are E-dependent.
Let us put our work in its proper perspective. The explicit results and observations of the present work are not entirely new. Some of these are discussed in the book by Chandrasekhar [13]. More recent and detailed studies are provided in [1,2,3,14]. However we have revisited the system from the Jacobi metric point of view. The Jacobi metric construction for charged massive particle and the subsequent analysis is completely new. Furthermore in the present work the Gaussian curvature induced by the Jacobi metric plays an essential role since the motion is geodesic in nature. We have compared results computed from Gaussian curvature consideration (in Jacobi metric approach) with similar results obtained via conventional scheme. For the charged probe, the interplay between the gravitational and Coulomb forces proves to be interesting.
3 Jacobi metric for neutral particle in Reissner -Nordström Geometry Reissner-Nordström metric is a spherically symmetric solution of the coupled Maxwell Einstein gravity. It represents a black hole with a mass M and a charge Q and is given by For Q = 0 the Schwarzschild metric is recovered. As is well known it has two horizons at r = r ± = M ± M 2 − Q 2 . The nature of the horizon singularities are different in Reissner -Nordström and Schwarzschild geometries. The latter is spacelike whereas the former is timelike and thus yielding richer possibilities regarding the nature of trajectories. There are timelike worldlines for particles that can cross r + horizon and skirting the singularity can move out to another spacetime region after crossing r − . On the contrary for Schwarzschild geometry, after crossing the event horizon at r = 2M the particle has no option but to fall towards the singularity. The r + -horizon acts like the r = 2M event horizon of Schwarzschild whereas r − -horizon is termed as the Cauchy horizon. Now, generalizing the result of Gibbons [5], the Jacobi metric corresponding to Reissner -Nordström solution is The first part is the conformal factor whereas the second factor is the optical metric. Due to spherical symmetry, we are allowed to study the system in the equatorial plane θ = π 2 without any loss of generality. This reduces the Jacobi metric to the form, Because of axial symmetry φ is a cyclic coordinate so that the angular momentum l is conserved, Now, together with (20), (21) yields that can be rewritten as This satisfies the standard result, Here τ is the proper time along geodesic and the angular momentum is Incidentally the above relation connects the proper time τ to the Jacobi path length s. Conventionally the trajectory is expressed in terms of a new variable u = 1 r : where, for the Reissner-Nordström case, we find, In the above h = l/m is the conserved angular momentum per unit mass. Thus (27,28) constitute the particle worldline or trajectory equation. It is straightforward to solve (27) and generate the following first order differential equation, where, C is a constant related to the energy per unit mass ∈= E m by This constitutes the first integral of motion and is one of our major results. Incidentally h ∈ is the impact parameter. We immediately notice a qualitative change the black hole charge has brought about in the trajectory. In comparison to the Schwarzshild case [5] the present result involves a quartic term in u. Writing the quartic polynomial in terms of roots, the following identities are recovered, Clearly these are extensions of analogous relations given in [5] for the Schwarzshild metric.
Our aim is to propose explicit solutions for the trajectories in the same manner as those derived in [5]. To that end let us quickly recapitulate the orbit equation for Schwarzschild black hole [15,5], with the explicit solution, The constant parameters A G , B g , C G , ω G are given in appendix, with the subscript G standing for the Gibbons solution [5].
Let us now come to our work. An explicit solution for the trajectory equation for Reissner -Nordström case (29) studied here is given by where A, B, k all are constants. One can solve for the constants perturbatively for small Q to first nontrivial order in Q 2 . The result is given in Appendix A. This constitutes one of our new results. This is a generic form of closed orbit. As we have pointed out in the Introduction, this type of analytic result for a generic orbit, indeed in an approximate sense of of small charge Q on the black hole) can compliment the graphical analysis of [1,2] without any approximations. It is straightforward to reduce our analysis to circular orbit orbits for which u = u c , a constant. We rewrite (29) as, Presence of a biquadratic term of u in f (u), compared to a cubic one in Schwarzschild geometry, is the qualitative change that leads to significant difference only for the orbits that cross the event horizon at r + and can skirt the singularity and come out of r = r + , (as discussed earlier), instead of terminating at the singularity at r = 0 as in the case of Schwarzschild black hole (see for example [13] for details).
For the occurrence of circular orbits the conditions are as follows, The relations change for the null geodesics which we are not considering at present (see for example [13]). We can easily obtain the expressions for energy and angular momentum of a circular orbit of radius r c = 1 uc from the above two equations. The expressions are, and, Incidentally, these results agree with [1,2]. The minimum radius for a stable circular orbit will occur at the point of inflection of the function f (u), i.e., Eliminating h 2 from the above equation using (41) we obtain, or, in terms of r c , From (44) , for Q 2 = 0 , we recover the well known result for Schwarzschild geometry, r c = 6M . However, neglecting Q 4 , a leading order correction to the radius is easily obtained, Using (38), (39) and (42), (37) takes the form, and the solution is given by, For Reissner-Nordstrom case, to O( Q 2 M ) we obtain This constitutes the first part of our work.

Jacobi metric for charged particle in Reissner -Nordström Geometry
The next level of generalization is to consider the trajectory of a probe with charge q in the presence of a charged black hole. Indeed this is a non-trivial extension to the previous case since an additional Coulomb interaction term of the form ∼ (qQ)/r is involved. In Reissner-Nordstrom geometry where a test particle has a charge per unit mass q, the only non vanishing component of the vector potential is A 0 and its motion is determined by the Lagrangian as of the form, However, we need not attempt to construct a generalized Reissner -Nordström solution starting from the Einstein-Maxwell point charge action. The Jacobi metric formalism provides a quick answer. Thus for the Reissner -Nordström case with a charged probe, the Jacobi metric is given by The only distinct feature which arises due to the probe charge q in this case is that the energy for a particle (having a turning point) that arrives at the event horizon will be , E = mqQ r + and this can be negative if qQ < 0 which gives rise to the theoretical speculation of generating energy from a black hole. Once again a restriction of the motion to the equatorial plane, i.e. θ = π 2 , reduces the Jacobi metric to In an identical fashion, as done in previous cases, we derive the trajectory equation for u = 1 r , where h = l/m is the conserved angular momentum per unit mass, C is a constant related to the energy per unit mass ∈= E m and the angular momentum per unit mass h is by the relation C = ∈ 2 −1 h 2 (as defined earlier in (30)).
Again we discuss our work of constructing the orbit. As we are considering upto the first order correction terms of Q 2 so we can drop the Q 2 q 2 term from (52) in our approximation. The solution of (52) can be written as, where A q , B q , k q , ω q all are constant and their approximate expressions are once again provided in the Appendix B.
The trajectory equation for the charged probe is, This is the other principal result of our paper. We now concentrate on circular trajectories. For the occurrence of circular orbits the conditions are as follows, and We can easily obtain the expressions for energy and angular momentum of a circular orbit of radius r c = 1 uc from the above two equations. The expressions for energy is, Now, upto order qQ (ignoring the term q 2 Q 2 ), the expression can be written as, Apart from some minor mismatch in numerical factors this result essentially agrees with [1,2]. Similarly for angular momentum, we have Thus, upto order of qQ , There is an interesting observation regarding a possible scaling of the charges following [13] where variations of 2 = (E/m) 2 and h = l/m against M u c are discussed with the scaling Q 2 = pM 2 , p being a numerical constant. The resulting relation for the latter for neutral probe is [13] However if we consider an identical scaling in our present case of with a charged probe, the relation turns out to be, Appearance of the parameter q is indicative of the fact that the Coulomb force is essentially non-geometric and hence the trajectories are not pure geodesic in nature. The effect of the probe charge, especially whether it is of same or opposite sign as the black hole charge, is quite striking. Intuitively we can argue that for the opposite sign case the probe charge effect will not be very significant because both the gravitational force and Coulomb force will be attractive and so qualitatively similar behavior to the neutral case will be observed. This is shown in Figure 1.
On the other hand, if the probe and black hole charges are of same sign the Coulomb force will be repulsive whereas the gravitational force is attractive as before. Interplay between these two forces produces an upper bound of the q parameter above which the results become unphysical. This is demonstrated in Figure 2.
Similar behavior for 2 = (E/m) 2 vs. r c /M is observed in Figure 3 and in Figure 4 where, negative and positive values of q are considered respectively. The minimum radius for a stable circular orbit will occur at the point of inflection of the function f (u), i.e., Eliminating h 2 from the above equation using (59) , we can write, or, in terms of r c , Let us define a parameter Λ = (qQE)/(mM ), in terms of which the above equation is rewritten as  Ignoring O(Q 4 ) terms we obtain Note that for q = 0 that is neutral probe the earlier result (48) is recovered. Considering q and Λ to be small (which is quite natural) we find This shows that the effective charge of the black hole increases for positive Λ that is when the probe charge q and black hole charge Q are of same sign but it decreases when they are of opposite sign. This concludes our discussion on charged particle trajectory in the presence of a charged black hole.

Gaussian curvature for Jacobi metric
Study of the Gaussian curvature for the Jacobi metric for Reissner-Nordstrom black hole, as experienced by a particle of fixed energy E is the other major contribution of our work. The Jacobi metric corresponding to Reissner-Nordstrom geometry (20) (considering the system in the equatorial plane θ = π 2 ) vanishes when The equation saturates for a critical value r c : This leads to an inequality of the form m 2 m 2 −E 2 ≥ Q 2 4M 2 and since m 2 m 2 −E 2 < 1, one has Q < 2M . It is interesting to recall that Q = M is the extremality condition and moreover Q < M is generally assumed to avoid the presence of a naked singularity (or violation of Cosmic Sensorship hypothesis by Penrose). Hence Q < 2M derived above is a weaker condition and does not add any further restrictions on the charge mass ratio.
Interpreting the following expression as Jacobi circumference [5] 2π the boundary r c actually reduces to a point since the circumference vanishes there. Condition for circular geodesics is derived from the extrema of the Jacobi circumference, Circular orbits exist provided (1 − 3M r + 2Q 2 r 2 ) ≥ 0. This means that to O(Q 2 /M ) the roots are Several examples of radii (to O(Q 2 )) for some specific values of energy are provided below: for m 2 = 0, the term (1 − 3M r + 2Q 2 r 2 ) must be equal to zero leading to r ≈ 3M − 2Q 2 3M . The outermost and innermost circular orbit radii are computed by extremizing r.h.s of (72). These are respectively Incidentally the last example agrees with our previously derived result (48). For Q = 0 the above reduce to the Schwarschild geometry results [5].
Let us now come to the explicit structure of Gaussian curvature K. As a warmup exercise let us compute K S for Schwarschild case. Expressing the Jacobi metric pertaining to Schwarschild geometry as, where, r ) 2 the Gaussian curvature K S is given by, For massless probe m = 0, the curvature simplifies to, For r ≥ 2M , K S is always negative [10].
Uncharged probe: Let us start with massless neutral probe. In a similar way as described above, for the Reissner-Nordstrom case, the Gaussian curvature to O(Q 2 ) is given by, where For m = 0 massless probe, the Gaussian curvature reduces to the form, Comparing with Q = 0 case (76) , the situation becomes more complicated and indeed, it is possible that the Q-contribution might reverse the sign of K RN . To verify this we define a convenient scaling E 2 = c 1 M 2 , r = aM, Q = bM and rewrite K RN in (79) as, To see the effect of Q we plot K RN vs.  The graphs indicate that, to the order of approximation we are considering, for massless probe, K RN stays negative, as in Schwarschild case [5] and as the energy of the particle increases, the magnitude of the Gaussian curvature decreases and asymptotically tends to become zero. Thus for massless probe, the presence of black hole charge has not too much effect except it increases the negativity of the Gaussian curvature.
Let us now consider a probe with mass m where the dynamics changes qualitatively from the massless probe discussed above. We introduce the parameterization E 2 = cm 2 , r = aM, Q = bM and plot K RN vs. c = E 2 /m 2 (Fig.s 6,7,8). The Gaussian curvature, (to order Q 2 ), near outer boundary a = 6 − 3 2 b 2 , is From Fig. 6 we notice that, for the energy range 0.7 < E 2 /m 2 < 0.8, for small b 2 (small Q), the Gaussian curvature K RN is positive but it changes sign and becomes negative for larger Q ∼ 0.5 and from Fig. 7 we see that the transition occurs around b 2 = 0.3. Interestingly in this case the K RN starts with a negative value and briefly reaches positive values before vanishing asymptotically. From Fig.8, we find that for c = E 2 /m 2 ≤ 0.2, K RN remains positive both for Q = 0 and for non-zero Q.
For inner boundary, we put a = 4 − b 2 and K RN is We have chosen the same value b 2 = 0.3 for which the sign change occurred at outer boundary. The plot in Fig. 9 indicates that, near inner boundary, for b 2 = 0.3, K RN tends to become more negative but the graphs shows similar qualitative nature for both zero and non zero Q. Charged probe: We now discuss the case of charged probe where an additional Coulomb interaction needs to be taken into account. The Jacobi metric corresponding to Reissner Nordstrom case with a charged probe (q) (in the equatorial plane) is The metric vanishes when where the q 2 Q 2 -term has been dropped for small q. This leads to, Again this can be thought of as a point where the Jacobi circumference vanishes. The circumference can be written from the metric as, Circular geodesics correspond to the extrema of the Jacobi circumference for which, For the Reissner Nordstrom geometry where the probe charge q is present, the Gaussian curvature can be written as, K RN P = K S + (C/D)qQ (88) where It is worthwhile to point out that the leading correction term depends on qQ that both linearly on q and Q which leads to interesting consequences. Firstly, the relative sign of black hole charge Q and probe charge q becomes important as it dictates the nature of Coulomb interaction, that is whether it is repulsive or attractive. Secondly, unlike the Reissner Nordstrom geometry with neutral probe, here we can drop Q 2 -terms due to the presence of qQ-terms. Now where the probe charge q is present, we parametrize Q = M/2, r = 15 8 pM and E 2 /m 2 = 1 and find the Gaussian Curvature as, From the above expression we can easily observe that for positive values of probe charge q or in the absence of the probe charge, the graphs will show similar qualitative nature (shown in fig.10) where as for negative values of probe charge i.e. for (2 + q) < 0 or, q < −2 , it will reverse the sign of Gaussian curvature K and the curvature starts with negative value and tends to become zero. Similarly, if we plot K vs. c = E 2 /m 2 for negative and positive values of probe charge q , the graphs shows similar qualitative nature for positive values of probe charge whereas it shows reverse nature for negative values of probe charge(shown in fig.11).

Conclusion
In the present work we have considered particle trajectories that are parameterized by constant energy value. This feature helps to visualize quickly the bounded and unbounded nature of particle orbits related  to particle energy. This characteristics is very succinctly incorporated in the Jacobi extension of least action principle. The formalism starts with the construction of the Jacobi metric where the (conserved) particle energy appears explicitly in the metric. As has been proved by Pin [9] a restricted variational principle a la Maupertuis with conventional metric and constant particle energy is equivalent to an unrestricted variational principle with Jacobi metric, which explicitly involves the particle energy. Hence the constant energy paths are still geodesics but of the Jacobi metric.
Exploiting this formalism we have studied the worldlines of both uncharged and charged probes in Reissner Nordstrom background. The former is a straightforward generalization of the Schwarzschild black hole as given by Gibbons [5] whereas the latter is a non-trivial extension since it has an additional Coulomb interaction. In both cases we have derived the circular orbit condition. The relative sign of the the probe and black hole (whether both have same sign or opposite sign) plays a significant role in determining nature of particle trajectory. Indeed, it should be stressed that there already exist series of recent [1] and earlier works [2,3] by Pugliese, Quevedo and Ruffini where this problem has been treated in detail in an exact way from graphical perspective. We, on the other hand, have attempted to compliment their works by providing analytic form of the orbits in an approximate scenario of small charge for the black hole and probe. We have also given expressions for the generic closed orbit and its reduction to circular orbits.
The Gauss curvature pertaining to Jacobi metric comprises another new and interesting aspect of our work. In earlier works [9,5], Gaussian curvature of the Jacobi metric played an important role in characterizing the nature of the particle worldlines in terms of open or closed orbits related to the particle energy. We have studied the properties of Gaussian curvature related to Jacobi metric for uncharged probe (massless and massive case) and also for the charged probe.
For massless uncharged probe, the presence of black hole charge Q only increases the negativity of Gaussian curvaturen whereas for the massive case it plays a significant role i.e. the black hole charge can reverse the sign of Gaussian curvature for holes with charge-mass ratio above a certain value.
When the probe is charged, an additional Coulomb interaction comes into play and the sign of the probe charge, whether it is of same sign or opposite to the black hole charge, becomes crucial. It is expected that if the charges are of the same sign, the effect becomes pronounced due to the competition between the repulsive Coulomb force and attractive gravitational force. On the other hand, the effect will be weaker if the charges are of opposite sign since both the Coulomb and gravitational forces will be attractive in nature. This is manifested in the angular momentum vs. circular radius and energy vs. circular radius graphs. However, curiously enough, the Gaussian curvature behaves in a different way. For a positively charged black hole, if the probe charge is also positive, the Gaussian curvature starts with a negative value and asymptotically tends to zero. This is similar to the case of a neutral probe. But if the probe has a negative charge (q < −2), (opposite to the positive black hole charge), the curvature stays positive throughout and asymptotes to zero.