Abstract
In order to classify and understand structure of the spacetime, investigation of the geodesic motions of massive and massless particles is a key tool. So the geodesic equation is a central equation of gravitating systems and the subject of geodesics in the black hole dictionary attracted much attention. In this paper, we give a full description of geodesic motions in three-dimensional spacetime. We investigate the geodesics near charged BTZ black holes and then generalize our prescriptions to the case of massive gravity. We show that electric charge is a critical parameter for categorizing the geodesic motions of both lightlike and timelike particles. In addition, we classify the type of geodesics based on the particle properties and geometry of spacetime.
Similar content being viewed by others
1 Introduction
Regarding the low energy limit of gravitational interactions, general relativity (GR) is a successful theory describing various phenomena. Unlike the successes of the mentioned theory, some issues such as accelerated expansion of the universe, the existence of dark matter [1], massive gravitons and several other subjects show the necessity of modifying this theory. There are various attempts for modifying GR, such as F(R) gravity [2,3,4,5,6,7], Lovelock gravity [8,9,10,11,12], Horava-Lifshitz gravity [13, 14], brane world cosmology [15,16,17], scalar-tensor theories [18,19,20,21,22,23,24,25,26], massive gravity [27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44], rainbow gravity [45,46,47,48,49,50,51,52], and massive gravity’s rainbow [53,54,55]. The mentioned generalized theories lead to interesting results in various aspects of black objects, especially their geometrical and physical properties. In this regard, we are going to study the effects of massive modification of Einstein gravity on the geometric structure of low dimensional black holes, especially, the geodesic motions of a test particle. The primary motivation of studying the geodesic motions of test particles around the massive objects come from some interesting astrophysical phenomena, such as the precession of the perihelion of Mercury and gravitational lensing. Due to the strong curvature effects, black holes have considerable influence on the geodesic motions.
One of the predictions of GR is the existence of black holes. These mysterious singular solutions can be described by the first static spherically symmetric solution of Einstein’s field equations found by Karl Schwarzschild in 1916. Since the interpretation of the Schwarzschild horizon as a one way membrane by Finkelstein in 1958 [56], considerable efforts have been made to study exact solutions of Einstein gravity with a singularity and their interesting properties. One major achievement in the theoretical analysis of these objects was finding the first three dimensional solution of Einstein field equations by Banados-Teitelboim-Zanelli (BTZ) in 1992 [57]. Later on, many people studied various aspects of BTZ black hole as a typical laboratory of high energy physics [27, 58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76] . BTZ black holes provide good tools for handling some conceptual questions in the context of AdS/CFT correspondence, quantum gravity, string theory models and gauge field theory [77, 78]. In the context of classical gravity, BTZ black holes have been studied in the presence of other gauge fields such as Maxwell and power Maxwell fields [79,80,81], Born–Infeld theory [82,83,84,85,86] and other nonlinear electrodynamics [87,88,89], and also massive gravity [90], gravity’s rainbow [91], massive gravity’s rainbow [92] and dilaton field [81, 93, 94]. Also higher dimensional black hole solutions with BTZ analogy (BTZ like black holes) have been investigated in Refs. [95, 96].
We can investigate the black hole properties in various ways by considering the motion of a test particle in the spacetime with a black hole, described by the geodesic equation. Exact analytical solutions to the geodesic equation provide us with the best basic understanding of geodesic motion, but may not always be possible. We can then employ analytical approximation schemes or numerical solutions. A first seminal paper for studying the geodesic equations and its analytical solutions was by Hagihara in 1931, when he solved the geodesic equations of Schwarzschild black holes in terms of Weierstrass elliptic functions [97]. Afterwards, Darwin solved the geodesic equations by using the Jacobian elliptic functions [98, 99]. Also, the analytical solution of geodesic equations have been studied in four-dimensional Schwarzschild–de Sitter [100], Kerr [101, 102], Kerr–Newman [103], and Kerr–de Sitter [104] spacetime. Moreover, the higher dimensional solution of Schwarzschild, Schwarzschild-(anti-) de Sitter, Reissner–Nordström, Reissner Nordström–(anti-) de Sitter [105] and Myers-Perry spacetimes [106] have been investigated. Studying the geodesic equations have been extended to F(R) gravity and conformal gravity in BTZ and GMGHS (Gibbons-Maeda-Garfinkle-Horowitz-Strominger) black holes [107,108,109,110,111], regular and modified Hayward black holes [112, 113], charged dilatonic black holes [114], the singly spinning and (charged) doubly spinning black ring [115, 116], brane-world black hole [117], (rotating) black string [118], Schwarzschild, and Kerr pierced by black string spacetimes [119, 120]. Moreover, regarding three dimensional spacetime, the stability and existence of circular geodesics in a family of asymptotically AdS black holes in new massive gravity theory [121] and the null geodesics in a static circularly symmetric black hole spacetime [122] have been investigated before. Furthermore, in Ref. [123], the exact solutions of null and timelike geodesics are found and it is shown that the rotating (static) black hole is geodesically complete (incomplete).
In this paper, we consider (charged) BTZ black holes and its generalization to massive gravity. We provide a complete classification of the geodesic motion, and solve the geodesic equation. As far as possible, we follow the method which is introduced in Ref. [100]. For all of the previous works on the solutions of geodesic equations, the metric function was a polynomial function. In this paper, we investigate both polynomial and non-polynomial metric functions. A non-polynomial metric function appears in the charged cases and since there is no known exact analytical solution for this type, we resort to numerical solutions for this case.
The outline of this paper is as follow. First, we present the geodesic equation and effective potential in BTZ and its extension to massive gravity black holes in Sect. 2, and define acceptable regions of motion. After that, we investigate possible regions of motion and orbits for charged and uncharged black holes in Sect. 3. Then, we completely classify geodesic motions and solve the geodesic equation around the (charged) BTZ black holes in Sect. 4. The neutral black holes in massive gravity are treated in Sect. 5 and the charged ones in Sect. 6. Finally, we end the paper with some concluding remarks.
2 Three dimensional line element
In \((2+1)-\)dimensional spacetimes, the line element takes the following form
where \(\psi (r)\) is an arbitrary function of coordinate r which should be obtained based on the field equations. In this paper, we will consider four different cases for the metric function, \(\psi (r)\). The first one describes static BTZ black holes [57]
where \(\Lambda \) denotes the cosmological constant and \(m_{0}>0\) is an integration constant proportional to the total mass of the black hole. Note that here \(\Lambda <0\) is necessary to have a black hole solution (for a discussion of the singularity at \(r=0\) see [58]). Geodesic motions of massive/massless particles around the rotating uncharged BTZ black holes were investigated in Refs. [124, 125]. In the mentioned papers, it was shown that for the static case, there is no bound orbit for the positive geometrical mass of the black holes (\(m_{0}\)).
The line element (1) was generalized to the linearly charged solution of BTZ black holes with the metric function [95]
Both \(m_{0}>0\) and q are integration constants which are, respectively, related to mass and electric charge of the black hole solutions. In order to have an event horizon, here, \(\Lambda \) should be negative. It is also worth mentioning that \(r_{0}\) is an arbitrary constant with length dimension which is necessary to obtain dimensionless argument for the logarithmic function. In general, \(r_{0}\) is different from the length scale related to the cosmological constant (\(\Lambda \propto -l^{-2}\)), however, in [126], it was shown that the equality \(r_{0}=l\) is necessary to avoid an ensemble dependency.
The geodesic equation of charged BTZ black holes have been studied in [107]. Since the charge term of the metric function in Ref. [107] is positive, its related black hole has a Schwarzschild like horizon (spacelike singularity). Here, we consider real valued charged BTZ black hole with two horizon (timelike singularity). The roots of the metric function (3) have been reported in Ref. [95]
where \(\mathrm {W}(x)\) denotes the principal branch of the Lambert W function and \(\mathrm {W}(-1,x)\) the branch with \(\mathrm {W}(-1,x)\le -1\). The only physically acceptable solution of Eqs. (4) and ( 5) is
The metric function of BTZ black holes in massive gravity is [90]
where \(m_{0}>0\) as before and m, c and \(c_{1}\) are three constants related to the massive gravity (see Ref. [90] for more details). Here, we define a new massive parameter \(m^{\prime }\) to combine all massive parameters, as follow
Note that in the massive case \(\Lambda <0\) is not necessary for black hole solutions and we may as well consider \(\Lambda >0\).
Following Ref. [90], we find that by considering specific values for different parameters, the metric function (7) may have two roots or one extreme root (we are not interested in naked singularity). In addition, since the constant c is positive, the sign of \( m^{\prime }\) depends on the positive or negative sign of \(c_{1}\) (see [30] for more details). The uncharged BTZ black holes in massive gravity have a curvature singularity at \(r=0\).
In the charged black holes in massive gravity, the metric function \(\psi (r)\) is obtained as [90]
Following Refs. [79,80,81, 90], one finds the mentioned metric functions describe a three dimensional spacetime with a singularity located at \(r=0\). It is also notable that the singularity of the charged solutions is timelike, while it is spacelike for uncharged ones (for more details regarding the horizon and geometry of the mentioned spacetimes, we refer the reader to [79,80,81, 90]).
We will analyze all four cases introduced above in terms of the motion of test particles, both for massless and massive particles. Since we are working in static and spherical symmetric spacetimes, we can immediately obtain two conserved quantities, energy and angular momentum, as
The Lagrangian \({\mathcal {L}}\) of a test particle is given by
where \(\epsilon \) takes the values 1 and 0 for massive and massless particles, respectively, and \(\lambda \) is the affine parameter for massless particles and the proper time for massive particles.
From this, we then find the following geodesic equations
In addition, considering Eq. (12), an effective potential can be introduced as
Before solving the equations of motion we analyze the structure of possible types of geodesic motion for the test particles. The major point in these analyzes is that Eq. (13) implies \(P(r)\ge 0\) as a necessary condition for the existence of a geodesic. The real roots of P(r) are related to intersection points of \(E^2\) and \(V_{\mathrm {eff}}\). Equivalently, from (14) the acceptable region of motion is \( E^2\ge V_{\mathrm {eff}}\). The number of real roots of P(r) characterize the shape of the orbit [100]. For a given set of parameters, P(r) may have a certain number of roots. By varying E and L, and fixing other parameters the number of real roots and, therefore, the possible types of orbits, will change [101].
3 Classification of geodesic motion
Changes in the possible orbit configurations happen if two real zeros of P(r) merge to a double zero. The corresponding constants of motion can be obtained by solving \(P(r)=0\) and \(\frac{dP(r)}{dr}=0\) for \(E^{2}\) and \(L^{2}\) . Obtaining \(E^{2}\) and \(L^{2}\) for these limit cases is therefore crucial for studying the motion of particles. In all four cases introduced in Sect. 2, the function P(r) has a double root \(r=0\) for all parameters. Therefore, we can introduce a new function \({\widetilde{P}}(r)\) by
Instead of searching for roots of P(r), we may solve \({\widetilde{P}}(r)=0\) and \(\frac{d{\widetilde{P}}(r)}{dr}=0\) for \(E^{2}\) and \(L^{2}\), for both massive and massless particles.
The asymptotic behavior of P(r) is very important since it determines whether particles have flyby or bound orbits. A particle reaches infinity if P(r) is positive for \(r\rightarrow \infty \). In the limit of large r we find
For negative \(\Lambda \), massive particles may therefore never reach infinity.
3.1 Uncharged cases
In the uncharged case, the metric function is a polynomial function. By considering P(r), we can discuss about possible motions. In the uncharged cases that we study in this paper, P(r) is a polynomial in r of order not more than 6. Therefore, \({\widetilde{P}}(r)\) is of order not more than 4.
For all uncharged cases the following kinds of orbits can be identified.
Flyby orbits: r starts from \(\infty \), then approaches a periapsis \( r=r_{\min }\) and goes back to \(\infty \).
Bound orbits: r changes between two boundary values \(r_{1} \le r \le r_{2}\) with \(0< r_{1} \le r_{2} < \infty \).
Terminating bound orbits: r starts in \(\left( 0,r_{m}\right] \) for \(0<r_{m} < \infty \) and it falls into the singularity at \(r=0\).
Terminating escape orbits: r comes from \(\infty \) and falls into the singularity at \(r=0\).
The polynomial \({\widetilde{P}}(r)\) can have up to four real zeros, which together with the asymptotic behaviour could give rise to ten different orbit configurations. It will however turn out that \(P(r) \rightarrow +\infty \) for \( r\rightarrow \infty \) always corresponds to an even number of roots whereas \(P(r) \rightarrow -\infty \) corresponds to an odd number. Therefore, only the following five different orbit configurations are possible:
In region (0), \({\widetilde{P}}(r)\) has no positive real roots and \({\widetilde{P}}(r)>0\) for \(r\ge 0\). The possible orbit types are terminating escape orbits.
In region (1), \({\widetilde{P}}(r)\) has one positive real root \( r_{1}\) with \({\widetilde{P}}(r)\ge 0\) for \(0\le r\le r_{1}\) with possible orbit types being terminating bound orbits.
In region (2), \({\widetilde{P}}(r)\) has two positive real roots \(r_{1}\le r_{2}\) with \({\widetilde{P}}(r)\ge 0\) for \(0\le r\le r_{1}\) and \( r_{2}\le r\) with possible orbit types being flyby and terminating bound orbits.
In region (3), \({\widetilde{P}}(r)\) has three positive real zeros \(r_{1}\le r_{2}\le r_{3}\) with \({\widetilde{P}}(r)\ge 0\) for \(0\le r\le r_{1}\) and \(r_{2}\le r\le r_{3}\) with possible orbit types being bound and terminating bound orbits.
In region (4), \({\widetilde{P}}(r)\) has four positive real zeros \(r_{1}\le r_{2}\le r_{3}\le r_{4}\) with \({\widetilde{P}}(r)\ge 0\) for \(0\le r\le r_{1}\) and \(r_{2}\le r\le r_{3}\) and \(r_{4}\le r\) with possible orbit types being bound, terminating bound and flyby orbits .
Note that the regions (1) and (3) only appear for \(\Lambda <0\), whereas the regions (0), (2), and (4) are only possible if \(\Lambda >0\) in the case of massive particles (\(\epsilon =1\)), and if \((1/b^2)+\Lambda >0\) in the case of massless particles (\(\epsilon =0\)). Here \(b=L/E\) is the impact parameter.
3.2 Charged cases
For the charged case, the presence of the logarithmic function spoils the polynomial form of P(r). Note that although P(r) always has a double zero at \(r=0\), for the charged case, the function \({\widetilde{P}}(r)\) approaches \(-\infty \) for \(r\rightarrow 0\). Therefore, based on Bolzano’s theorem [127], since \({\widetilde{P}}(r)\) is continuous in \( (0,+\infty )\), if \(\Lambda >0\) there is always a real and positive root.
Since the logarithmic charge term changes the asymptotic behaviour of \( {\widetilde{P}}(r)\) for \(r\rightarrow 0\) from a positive finite value to negative infinity, this suggests the presence of a small additional root \(r_{ \mathrm {extra}}\) as compared to the uncharged case. If \(r_{extra}<r_{-}\) (where \(r_{-}\) is inner horizon), then a flyby or bound orbit will be reflected by the singularity and cross the horizons multiple times and entering in a new copy of the universe. These orbits are called two-world escape orbit and many-world bound orbit, respectively (see [105, 118, 128, 129] for more detail).
Numerical calculations show that the following orbits can be identified for the charged cases
Flyby orbits: r starts from \(\infty \), then approaches a periapsis \( r=r_{\min }\) for \(r_{\min }>r_{+}\), and goes back to \(\infty \).
Bound orbits: r changes between two boundary values \(r_{1}\) \(\le r\) \(\le \) \(r_{2}\), with \(r_{1},r_{2}>r_{+}\) or \(0<r_{1},r_{2}<r_{-}\) where \( r_{+}\) is event horizon and \(r_-\) the inner horizon.
Many-world bound orbits: r changes between to boundary values \(r_{1}\) \(\le r\) \(\le \) \(r_{2}\) , with \(0<r_{1}\) \(<r_{-}\) and \(r_{2}>r_{+}\).
Two-world escape orbits: with \(r>r_{1}\) and \(0<r_{1}\) \(<r_{-}\).
For the charged case, we can find six regions with different number of roots which leads to following orbits:
In region (0), \({\widetilde{P}}(r)\) has no real positive root with \({\widetilde{P}}(r)\) \(<0\) and the motion is impossible .
In region (1), \({\widetilde{P}}(r)\) has one real positive root \( 0<r_{1}<r_{-}\) with \({\widetilde{P}}(r)\ge 0\) for \(\,r_{1}\le r\) and possible orbit types are two-world escape orbits.
In region (2), \({\widetilde{P}}(r)\) has two real positive roots \(0<r_{1}<r_{-}<r_{+}<r_{2}\) with \({\widetilde{P}}(r)\ge 0\) for \(r_{1}\le r\le r_{2}\) and possible orbit types are many-world bound orbits.
In region (3), \({\widetilde{P}}(r)\) has three real positive roots \(0<r_{1}<r_{-}<r_{+}<r_{2}\le r_{3}\) with \({\widetilde{P}}(r)\ge 0\) for \(r_{1}\le r\le r_{2}\) and \(r_{3}\le r\) with possible orbit types being many-world bound and flyby orbits.
In region (4), \({\widetilde{P}}(r)\) has four real positive roots \(0<r_{1}<r_{-}<r_{+}<r_{2}\le r_{3}\le r_{4}\) with \({\widetilde{P}} (r)\ge 0\) for \(r_{1}\le r\le r_{2}\) and \(r_{3}\le r\le r_{4}\) with possible orbit types being bound and many-world bound orbits.
In region (5), \({\widetilde{P}}(r)\) has five real positive roots \(0<r_{1}<r_{-}<r_{+}<r_{2}\le r_{3}\le r_{4}\le r_{5}\) with \( {\widetilde{P}}(r)\ge 0\) for \(r_{1}\le r\le r_{2}\) and \(r_{3}\le r\le r_{4}\) and \(r_{3}\le r\), with possible orbit types being many-world bound, bound and flyby orbits.
Note that region (2) is the only possible region for the charged BTZ black holes. Region (4) only appears for particle motion in the charged massive BTZ black hole with \(\Lambda <0\), \(m^{\prime }>0\).
4 BTZ black holes
In this section, we study the test particle’s motion in charged BTZ black holes. Although the geodesic motion of a particle around the uncharged black hole has been studied before [124, 125], for the sake of comparison, we will first derive the equations of motion for the uncharged BTZ black holes before proceeding to the charged case.
4.1 Uncharged BTZ black holes
4.1.1 General classification of motion
Substituting Eq. (2) into Eq. (13), we obtain
Let us first consider the case of massive test-particles with \(\epsilon =1\). Upon inspection of the effective potential (18) and keeping in mind that \(\Lambda <0\), we see that it diverges to infinity for \( r\rightarrow \infty \) and to minus infinity for \(r\rightarrow 0\). Its derivation with respect to r is
which is always positive for \(r>0\). Referring to (12), this implies that all massive particle trajectories have some outer turning point \(r_{0}>0\) and eventually have to cross the black hole horizon at \( r_{+}=m_{0}/\sqrt{-\Lambda }\).
The same result can be inferred from (19). By Descartes’ rule of signs, P(r) posses at most one positive real zero \( r_{0} \), and \(P(r)\rightarrow -\infty \) for \(r\rightarrow \infty \). Massive test particles are therefore bound to the region \(0\le r\le r_{0}\).
Now let us turn to massless particles, \(\epsilon =0\). Here, it is convenient to rescale the affine parameter \(\lambda \) such that the equation of motion simplifies to
where \(b=L/E\) is the impact parameter. Then, the effective potential \({\hat{V}}_{\mathrm {eff}}\) approaches \(-\Lambda b^{2}\) from below for \(r\rightarrow \infty \) and diverges to minus infinity for \(r\rightarrow 0\). As for derivation with respect to r of \({\hat{V}}_{\mathrm {eff}}\),
is positive for \(r>0\), this implies that a photon may reach infinity if \( -\Lambda b^{2}<1\), and otherwise it is bounded by an outer turning point \( r_{0}\). From this analysis, it is also clear that circular orbits cannot exist.
4.1.2 Analytic solution of geodesic equations
We use the substitution \(r^{2}=\frac{1}{u}\) and slightly rewrite Eq. (19) to obtain
where \(c_{1}=\frac{E^{2}}{L^{2}m_{0}}+\frac{\Lambda }{m_{0}}+\frac{\epsilon }{L^{2}}\), \(c_{2}=\frac{\epsilon \Lambda }{m_{0}L^{2}}\). This integral has the solution
Equivalently, we find u as a function of \(\varphi \),
where \(\xi \) is related to \(u_{{0}}\) as
Therefore, \(r(\varphi )\) is
which is valid for both timelike and lightlike geodesics.
All the possible types of null geodesics in BTZ black holes are plotted in Fig. 1. Types of orbits mentioned in these plots have been discussed in Sect. 3.
4.2 Charged BTZ black holes
For the line element (1), the linearly charged solution of BTZ black holes has been obtained as [95]
so from Eq. (13), P(r) and \(V_{\mathrm {eff}}\) are
In contrast to the uncharged case, here we can find circular orbits, which are given by \(\frac{dr}{d\lambda }=0\) and \(\frac{d^{2}r}{d\lambda ^{2}}=0\). These two conditions are equivalent to \(P(r)=0\) and \(\frac{dP}{dr}=0\). We can solve these two equations for the squared energy and angular momentum as
for massive particles \((\epsilon =1)\). Let us discuss these two equations. From Eq. (28), it is evident that \(E^{2}>0\) is valid only for a negative denominator and thus
In addition, the numerator of Eq. (29) has to be negative, and therefore, keeping in mind that \(\Lambda <0\)
Considering the constraints (30) and (31), simultaneously, we find that circular orbits exist if the following relation is satisfied,
It is notable that such inequality implies the absence of horizon in the domain \(r>0\). This confirms that for the charged BTZ spacetime circular orbits may only exist around naked singularities, but not for the case of a black hole.
Therefore, similar to the previous section, P(r) has the same number of roots in each combination of E, L and constant parameters. For \(r\ne r_{0}\), we can adjust arbitrary combination of E, L and other constant parameters, in which we obtain two real positive roots for P(r). However, for \(r=r_{0}\), the uncharged case is recovered in which P(r) has always one root. All the possible types of timelike orbits have been plotted in the Fig. 2. Properties of the orbits mentioned in these figures have been presented in Sect. 3.
Let us now turn to massless particles. In this case (\(\epsilon =0\)) circular orbits exist for
where \(b=L/E\) is the impact parameter. Due to the fact that \(\Gamma \) is real, in order to have physically acceptable motion for the massless particles, the denominator of Eq. (33) must be negative. Therefore, the following constraint must hold, keeping again in mind that \( \Lambda <0\),
This results in the same inequality (32) as for the case of massive particles which indicates absence of an event horizon.
The many-world bound orbit of null geodesics are shown in Fig. 3.
5 Uncharged Massive BTZ black hole
5.1 General classification of motion
The metric function for the BTZ black holes in massive gravity is presented in Eqs. (7) and (8) [90]. By substituting Eqs. (7) and (8) into Eqs. (13) and (14), we have
Let us first discuss massive particles (\(\epsilon =1\)). From the form of P(r) , according to the Descartes rule of signs, it is clear that for \( \Lambda <0\), there may be one or three positive real zeros, whereas for \( \Lambda >0\) there are four, two, or no positive real zeros, depending on the values of E and L as well as the sign of \(m^{\prime }\). This points to the existence of circular orbits, which are given by \(\frac{dr}{d\lambda }=0\) and \(\frac{d^{2}r}{d\lambda ^{2}}=0\). These two conditions are equivalent to \(P(r)=0\) and \(\frac{dP}{dr}=0\). We may solve these two equations for \(E^{2}\) and \(L^{2}\), which gives for massive particles \((\epsilon =1)\)
It is clear from the expression for \(E^{2}\) that circular orbits can only exist if \(2m_{0}-m^{\prime }r<0\), which is for \(r>0\) only possible if \( m^{\prime }>0\), which then gives us \(r>2m_{0}/m^{\prime }\). The equation (39) for L then implies that \(m^{\prime }-2\Lambda r >0\) has to hold, which is automatically fulfilled for \(\Lambda <0\). For \( \Lambda >0 \) we find \(r<m^{\prime }/(2\Lambda )\).
The circular orbit is stable if its radius corresponds to a maximum of P(r), i.e. if the second derivative of P(r) is negative. The second derivative of P(r) together with Eqs. (38) and (39) reads
This expression can be solved for the radius of the circular orbit,
For \(\Lambda <0\), we find one positive zero (for the negative sign before the root in (40)) with stable orbits for larger radii, which implies that an innermost stable circular orbit (ISCO) exists. The radius of this ISCO approaches \(3m_0/m^{\prime }\) for small \(\Lambda \) and \( 8m_0/3m^{\prime }\) for large negative \(\Lambda \). We plotted the ISCO for fixed \(m_0\) in Fig. 4.
Interestingly, for \(\Lambda >0\) stable orbits only exist if \(r_c\) from (40) is real, that is, if and only if \((m^{\prime 2}>16\Lambda m_0\). Note that \((m^{\prime 2}<4 \Lambda m_0\) can be rewritten as \(m^{\prime }/(2\Lambda )<2m_0/m^{\prime }\), which implies that no circular orbits exist according to our discussion of (38) and (39) above. Summarized, for \(m^{\prime }>0\) and \(\Lambda >0\) circular orbits exist for \(r \in [2m_0/m^{\prime },m^{\prime }/(2\Lambda )]\), and may be stable if \(m^{\prime }/(2\Lambda )>8m_0/m^{\prime }\). Then the circular orbits are stable in between the two radii given by (40) and we have an innermost and an outermost stable circular orbit. We show the important radii for the case \(m^{\prime }>0\), \(\Lambda >0\) in Fig. 5.
Now let us turn to massless particles \((\epsilon =0)\), where we introduce \(b=L/E\) as the impact parameter. We then see that the leading coefficient of P(r) may change its sign for \(1/b^{2}=-\Lambda \). If \( 1/b^{2}<-\Lambda \) the polynomial P(r) diverges to \(-\infty \) for \( r\rightarrow \infty \) which implies that the photon may not reach infinity. Moreover, P(r) may have at most one positive real zero (\(r_{0}\)), and therefore, it is bound in a region \(0\le r\le r_{0}\). On the other hand, if \(1/b^{2}>-\Lambda \), then P(r) is positive for \(r\rightarrow \infty \), and one finds the photon may reach infinity, and P(r) has two or no positive real zeros. This again points to the existence of a circular orbit for this case. If we solve the two conditions \(P(r)=0\) and \(\frac{dP}{dr}=0\) for the impact parameter b and the radius of the circular orbit \(r_{c}\), we find
From Eq. (41) for b we infer that \(4\Lambda m_{0}-m^{\prime 2}<0\), or equivalently \(\Lambda <m^{\prime 2}/(4m_{0})\) is necessary for the existence of circular orbits. This circular orbit is stable if it corresponds to a maximum of the polynomial P(r). For its second derivative with respect to r, we find
which implies that the circular photon orbit is unstable. Interestingly, the radius of the photon orbit only depends on the ratio \(m_{0}/m^{\prime }\), but not on \(\Lambda \) or the individual parameters \(m_{0}\) and \(m^{\prime }\) . For \(\Lambda \) this is to be expected, as it does not enter as an independent parameter in the equation of motion just as in four dimensional Schwarzschild-de Sitter spacetime, but this seems surprising for the independent parameters \(m_{0}\) and \(m^{\prime }\). However, we correctly recover that the circular orbit vanishes (its radius shifts to infinity) for \(m^{\prime }\rightarrow 0\).
The results of Eqs. (38), (39) and ( 41) for both massive and massless particles are given in Figs. 6, 7 and 8. Additional properties of these figures have been presented in Sect. 3.
5.2 Analytic solution of geodesic equations
Here, we discuss the analytical solution of Eq. (13). The function P(r) in given in Eq. (37) is a polynomial of degree 6, which can be written in the following form
where the coefficients \(a_{i}\) are
By substitution \(r=u^{-1}+r_{M}\) into Eq. (13) , where \(r_{M}\) is a root of P(r), for instance \(r_{M}=0\), we find
and the coefficients \(b_{j}\) can be calculated as
Now, we are in a position to obtain analytic solutions of Eq. (46 ) for massless and massive particles.
5.2.1 Null geodesics \((\epsilon =0)\)
By substituting \((\epsilon =0)\) into Eq. (46), it transforms to the following
where \(\alpha _{j}=\left. b_{j}\right| _{\epsilon =0}\).
The roots of Eq. (48) lead to an elliptic type differential equation. Another substitution \(u=\frac{1}{\alpha _{3}}(4y- \frac{\alpha _{2}}{3})\) transforms \(P_{3}(u)\) into the following form
in which the coefficients \(g_{2}\) and \(g_{3}\) are
It is known that the solution of Eq. (49) is the Weierstrass function \(\wp (\varphi )\) in the following form
where
in which \(r_{0}\) and \(\varphi _{0}\) are the initial values of the differential equation. Therefore, the general solution of Eq. (48) is
which is in agreement with the result of Ref. [107].
5.2.2 Timelike geodesics \((\epsilon =1)\)
By substituting \((\epsilon =1)\) into Eq. (46), it can be transformed to
where \(\beta _{j}=\left. b_{j+2}\right| _{\epsilon =1}\). Equation (54) is hyperelliptic type and its analytical solution is given by the Kleinian sigma function [101, 106, 130]
in which \(\sigma _{i}\left( z\right) \) denotes the derivative of the Kleinian sigma function with respect to the \(i^{th}\) component of z
where \(\left( 2\omega ,2\omega ^{^{\prime }}\right) \) is the period matrix, \( \left( 2\eta ,2\eta ^{^{\prime }}\right) \) is the period matrix of the second kind, C is a constant that can be given explicitly but does not enter (55) and \(\tau =\omega ^{-1}\omega ^{^{\prime }}\). The theta function is defined as follow
in which g, h are two dimensional vectors related to the vector of Riemann constants. They are defined as \(g=(1/2,1/2)\), \(h=(0,1/2)\). The argument \( \varphi _{\sigma }\) in (55) is defined as
where f is given by the condition \(\sigma \left( \varphi _{\sigma }\right) =0\). The constant \(\varphi _{in}\) reads
and only depends on the initial values \(\varphi _{0}\) and \(u_{0}\). The general solution of the radial coordinate r is given by
Equations (53) and (60) completely describe the motion of massive and massless particles in this spacetime. All the possible type of orbits in BTZ black holes of massive gravity have been plotted in Figs. 9, 10 and 11. Types of orbits mentioned here have been introduced in the total classification Sect. 3.
6 Charged black holes in massive gravity
6.1 General classification of motion
In the charged black holes in massive gravity, the metric function \(\psi (r)\) is obtained as [90]
From Eq. (13), we then find P(r) as
Again, we can look for circular orbits by solving \(P(r)=0\) and \(\frac{dP(r)}{ dr}=0\) for the squared energy \(E^{2}\) and angular momentum \(L^{2}\). For massive particles \((\epsilon =1)\) we find
As expected, for \(q=0\) Eqs. (63) and (64) reduce to Eqs. (38) and (39), respectively, and for \(m^{\prime }=0\) they reduce to Eqs. (28) and (29). From equation (63) we infer the inequality
For \(m^{\prime }>0\) the function R diverges to infinity at \(r=0\) and \( r=\infty \), and has only a single extrema, a minimum at \(r=4q^2/m^{\prime }\). It may have two zeros \(\Gamma _{\mathrm {m},0} < \Gamma _{\mathrm {m},-1}\) given by
If \(m^{\prime }\le 0\) the denominator R is monotonically decreasing on \( (0,\infty )\) and introduces therefore an upper bound on r given by \( r<\Gamma _{\mathrm {m},0}\). We can determine the relative position of the zeros of R and the horizons: let \(r_H\) be one of the horizons, i.e. \( \psi _4(r_H)=0\). Then we find from this
Therefore, R is positive at the inner horizon and negative at the event horizon. For \(\Lambda >0\) a cosmological horizon also exists, where R is positive again. We conclude that \(r_{H,\mathrm {inner}}<\Gamma _{\mathrm {m},0}<r_{H,\mathrm {event}}<\Gamma _{\mathrm {m},-1}<r_{H,\mathrm {cosmo}}\), if the respective horizons and zeros of R exist.
We can directly infer that for \(\Lambda <0\) and \(m^{\prime }\le 0\) circular orbits outside of a black hole cannot exist, as for charged BTZ black holes discussed in Sect. 4.2 and uncharged BTZ black holes in massive gravity discussed in Sect. 5.
In addition, the nominator in (64) needs to be negative. For \(\Lambda < 0\) this implies a lower bound on r given by
Note that this bound is the minimum of \(\psi _4\) and, therefore, is smaller than the event horizon. For \(\Lambda <0\) and \(m^{\prime }>0\) circular orbits around a black hole exist at \(r>\Gamma _{\mathrm {m},-1}\), as in the case of uncharged BTZ black holes in massive gravity discussed in Sect. 5.
From the condition that the nominator in (64) has to be negative, in the case \(\Lambda >0\) circular orbits are not allowed for \(m^{\prime }\le 0\), and have only a limited range for \( m^{\prime }>0\) given by
Note that the two bounds correspond to the minimum and the maximum of the metric function \(\psi _{4}\). Therefore, it only remains to show that \(\Gamma _{\mathrm {m},-1}\) is smaller than the upper bound in (70): as \(r_{*}:=\Gamma _{\mathrm {m},-1}\) is a zero of R we find
It is clear that \(\psi _4\) is positive at \(r_*=\Gamma _{\mathrm {m},-1}\) which implies that \(r_*\) is smaller than the upper bound in (70). We therefore find circular orbits for \( \Lambda >0\) and \(m^{\prime }>0\) in the range
Let us now turn to massless particles \((\epsilon =0)\). We find for the circular orbits
where \(r_c\) is a solution of \(R=0\), see Eq. (65). Therefore, again for \(m^{\prime }<0\) there is no circular photon orbit outside of a black hole. If \(m^{\prime }>0\) a circular photon orbit exists at \(r=\Gamma _{\mathrm {m},-1}\). In the limit \(m^{\prime }=0\) this reduces to the results of Sect.4.2, and for \(q=0\), the zero of R is \(r=2m_0/m^{\prime }\) as in Sect. 5. The circular photon orbit is always unstable: if we insert \(b^2\) and r into the second derivative of P, we find stability for
which is however incompatible with a real \(\Gamma _{\mathrm {m},-1}\).
The results of Eqs. (63), (64) and (74) for both massive and massless particles are given in Figs. 12 and 13. (For more details we refer the reader to Sect. 3). We observe that a variation of q as well as the value and sign of cosmological constant have considerable effect on the geodesic motion of massless particles. More clearly, we find that increasing the cosmological constant or the charge leads to increasing the possibility of two world escape obits rather than flyby and many world bound orbits. According to these figures, for some values of the charge or the cosmological constant, there is no physical motion for large values of b. For massive particles, we see that the possible types of orbits are rather different for \(\Lambda <0\) and \(\Lambda >0\). In the case \(\Lambda <0\), it is impossible to reach radial infinity and bound orbits outside the horizons do not exist. On the other hand, for \( \Lambda >0\) all parameter combinations allow orbits reaching infinity and bound orbits outside the horizons.
6.2 Numerical solution of geodesic equations
Now, we are in a position to discuss both null and timelike geodesics of charged BTZ solutions in massive gravity. Since the metric function is no longer a polynomial function, to our knowledge, there is no analytical solution for the charged cases. Therefore, we use the Runge–Kutta–Fehlberg numerical method. It is worthwhile to mention that we have checked this method for uncharged cases which we have analytical solutions, and the results were the same with high accuracy.
All the possible type of orbits in charged BTZ black holes in massive gravity with a positive cosmological constant are plotted by using the numerical analysis in Figs. 14 and 15 (see also Sect. 3 for more details). We restricted to positive \(\Lambda \) due to the interesting possibility of bound orbital motion.
As already familiar from the case of the charged BTZ solution, for both timelike and null geodesics, we find many-world bound orbits as well as the two-world escape orbits, that cross both the inner and the event horizon. For a positive cosmological constant, it is well known from four-dimensional black holes like the Schwarzschild-de Sitter solution, that flyby orbits can exist that are deflected from the cosmological barrier. This can be seen in Fig. 14b, e, where a massive test particle approaches the black hole rather straight and is reflected back without revolving around the black hole.
7 Summary and conclusions
In this paper, we analyzed the geodesic motion in (un)charged BTZ black hole and its generalization to massive gravity. We provided a complete classification of the possible types of geodesic motion for the four different spacetimes considered here, thereby investigating the effects of the various parameters of the metric on the geodesics. In particular, to our knowledge for the first time, we investigated the geodesic motion in a spacetime whose equation of motion for test particles has a non-polynomial structure. Moreover, we presented analytical solutions for the equations of motion for the uncharged cases, and used a Runge–Kutta–Fehlberg method for the spacetimes of the charged black holes, where an analytical solution could not be found.
While the uncharged BTZ black hole has rather a poor variety of orbital motions, in particular lacking bound orbital motion outside the horizon, the inclusion of a charge introduces a potential barrier within the inner horizon. This barrier reflects particles and light such that they have to cross the inner horizon for a second time, thereby entering another copy of the universe. This behaviour is known from some four-dimensional spacetimes, and the corresponding orbits are called two-world flyby or many-world bound orbits, respectively.
In three-dimensional massive gravity, the black hole solutions show an even richer structure of geodesic motion. First, we are no longer restricted to a negative cosmological constant, which gives rise to new types of geodesics. Second, the massive gravity parameter \(m^{\prime }\) introduced in Eq. (8) further enriches the structure of geodesic motion. We found that \(m^{\prime }>0\) is a necessary condition for the existence of circular orbits for both uncharged and charged black holes.
For the uncharged BTZ black hole in massive gravity, in addition to the general classification of geodesic motion, we derived the radius of the innermost stable and (for \(\Lambda >0\)) the outermost stable circular orbit. The results are plotted in Figs. 4 and 5. In particular, it was pointed out that stable bound orbital motion outside the event horizon is possible, and we showed an example in Fig. 9. Moreover, we calculated the radius of the unstable photon sphere, that is an important radius for the calculation of the shadow of the black hole.
Finally, for the most general case discussed in this paper, the charged BTZ black hole in massive gravity, we derived the range of radii where circular orbits may exist. Based on this, we presented a complete classification of geodesic motion. However, due to the complicated structure of the equations, we did not determine the innermost stable circular orbit, and leave this point for future research.
Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors’ comment: This is a theoretical study and no experimental data has been listed.]
References
N. Jarosik et al., Astrophys. J. Supp. 192, 14 (2011)
M. Akbar, R.G. Cai, Phys. Lett. B 648, 243 (2007)
J.C.C. de Souza, V. Faraoni, Class. Quantum Gravity 24, 3637 (2007)
S. Capozziello, R. Cianci, C. Stornaiolo, S. Vignolo, Class. Quantum Gravity 24, 6417 (2007)
G. Cognola, E. Elizalde, S. Nojiri, S.D. Odintsov, L. Sebastiani, S. Zerbini, Phys. Rev. D 77, 046009 (2008)
T.P. Sotiriou, V. Faraoni, Rev. Mod. Phys. 82, 451 (2010)
S.H. Hendi, B. Eslam Panah, C. Corda, Can. J. Phys. 92, 76 (2013)
D. Lovelock, J. Math. Phys. 12, 498 (1971)
D. Lovelock, J. Math. Phys. 13, 874 (1972)
N. Deruelle, L. Farina-Busto, Phys. Rev. D 41, 3696 (1990)
M.H. Dehghani, N. Alinejadi, S.H. Hendi, Phys. Rev. D 77, 104025 (2008)
M.H. Dehghani, N. Bostani, S.H. Hendi, Phys. Rev. D 78, 064031 (2008)
P. Horava, Phys. Rev. D 79, 084008 (2009)
P. Horava, Phys. Rev. Lett. 102, 161301 (2009)
P. Brax, C. van de Bruck, Class. Quantum Gravity 20, R201 (2003)
L.A. Gergely, Phys. Rev. D 74, 024002 (2006)
M. Demetrian, Gen. Relativ. Gravit. 38, 953 (2006)
C. Brans, R.H. Dicke, Phys. Rev. 124, 925 (1961)
T.P. Sotiriou, Class. Quantum Gravity 23, 5117 (2006)
K.I. Maeda, Y. Fujii, Phys. Rev. D 79, 084026 (2009)
R.G. Cai, Y.S. Myung, Phys. Rev. D 56, 3466 (1997)
M.H. Dehghani, J. Pakravan, S.H. Hendi, Phys. Rev. D 74, 104014 (2006)
S.H. Hendi, J. Math. Phys. 49, 082501 (2008)
J.A. Agudelo, J.R. Nascimento, AYu. Petrov, P.J. Porfirio, A.F. Santos, Phys. Lett. B 762, 96 (2016)
S. Bhattacharya, K.F. Dialektopoulos, A. Enea Romano, T.N. Tomaras, Phys. Rev. Lett. 115, 181104 (2015)
S.H. Hendi, S. Panahiyan, B. Eslam Panah, Z. Armanfard, Eur. Phys. J. C 76, 396 (2016)
D. Birmingham, S. Mokhtari, I. Sachs, Phys. Rev. D 82, 124059 (2010)
L. Alberte, A.H. Chamseddine, V. Mukhanov, JHEP 04, 004 (2011)
C. de Rham, G. Gabadadze, A.J. Tolley, Phys. Rev. Lett. 106, 231101 (2011)
D. Vegh, arXiv:1301.0537
F. Capela, P.G. Tinyakov, JHEP 04, 042 (2011)
R. Banerjee, S. Gangopadhyay, D. Roy, JHEP 10, 121 (2011)
F. Sbisa, G. Niz, K. Koyama, G. Tasinato, Phys. Rev. D 86, 024033 (2012)
Q.G. Huang, Y.S. Piao, S.Y. Zhou, Phys. Rev. D 86, 124014 (2012)
C. Deffayet, J. Mourad, G. Zahariade, JCAP 01, 032 (2013)
P. Martin-Moruno, M. Visser, Class. Quantum Gravity 30, 155021 (2013)
E. Babichev, A. Fabbri, Class. Quantum Gravity 30, 152001 (2013)
Y.F. Cai, E.N. Saridakis, Phys. Rev. D 90, 063528 (2014)
Yl Zhang, R. Saito, Dh Yeom, M. Sasaki, JCAP 02, 022 (2014)
R.G. Cai, Y.P. Hu, Q.Y. Pan, Y.L. Zhang, Phys. Rev. D 91, 024032 (2015)
S.H. Hendi, B. Eslam Panah, S. Panahiyan, JHEP 11, 157 (2015)
R. Dong, D. Stojkovic, Phys. Rev. D 92, 084045 (2015)
T.Q. Do, Phys. Rev. D 93, 104003 (2016)
S.H. Hendi, G.H. Bordbar, B. Eslam Panah, S. Panahiyan, JCAP 07, 004 (2017)
J. Magueijo, L. Smolin, Class. Quantum Gravity 21, 1725 (2004)
G. Amelino-Camelia, M. Arzano, G. Gubitosi, J. Magueijo, Phys. Rev. D 88, 041303 (2013)
R. Garattini, M. Sakellariadou, Phys. Rev. D 90, 043521 (2014)
A.F. Ali, M. Faizal, M.M. Khalil, JHEP 12, 159 (2014)
S.H. Hendi, G.H. Bordbar, B. Eslam Panah, S. Panahiyan, JCAP 09, 013 (2016)
A. Chatrabhuti, V. Yingcharoenrat, P. Channuie, Phys. Rev. D 93, 043515 (2016)
B. Eslam Panah, G.H. Bordbar, S.H. Hendi, R. Ruffini, Z. Rezaei, R. Moradi, Astrophys. J. 848, 24 (2017)
Z.-W. Feng, S.Z. Yang, Phys. Lett. B 772, 737 (2017)
S.H. Hendi, B. Eslam Panah, S. Panahiyan, Phys. Lett. B 769, 191 (2017)
S.H. Hendi, M. Momennia, B. Eslam Panah, S. Panahiyan, Phys. Dark Universe 16, 26 (2017)
Y. Heydarzade, P. Rudra, F. Darabi, A. Farag Ali, M. Faizal, Phys. Lett. B 774, 46 (2017)
D. Finkelstein, Phys. Rev. 110, 965 (1958)
M. Banados, C. Teitelboim, J. Zanelli, Phys. Rev. Lett. 69, 1849 (1992)
M. Banados, M. Henneaux, C. Teitelboim, J. Zanelli, Phys. Rev. D 48, 1506 (1993)
S. Nojiri, S.D. Odintsov, Mod. Phys. Lett. A 13, 2695 (1998)
R. Emparan, G.T. Horowitz, R.C. Myers, JHEP 01, 021 (2000)
S. Hemming, E. Keski-Vakkuri, P. Kraus, JHEP 10, 006 (2002)
B. Sahoo, A. Sen, JHEP 07, 008 (2006)
M. Cadoni, M.R. Setare, JHEP 07, 131 (2008)
M. Park, Phys. Rev. D 77, 026011 (2008)
J. Parsons, S.F. Ross, JHEP 04, 134 (2009)
M. Akbar, H. Quevedo, K. Saifullah, A. Sanchez, S. Taj, Phys. Rev. D 83, 084031 (2011)
Y.S. Myung, Y.W. Kim, T. Moon, Y.J. Park, Phys. Rev. D 84, 024044 (2011)
L. Hodgkinson, J. Louko, Phys. Rev. D 86, 064031 (2012)
T. Moon, Y.S. Myung, Phys. Rev. D 86, 124042 (2012)
E. Frodden, M. Geiller, K. Noui, A. Perez, JHEP 05, 139 (2013)
M. Eune, W. Kim, S.H. Yi, JHEP 03, 020 (2013)
J.P.S. Lemos, G.M. Quinta, Phys. Rev. D 89, 084051 (2014)
M. Bravo-Gaete, M. Hassaine, Phys. Rev. D 90, 024008 (2014)
M.R. Setare, H. Adami, Phys. Rev. D 91, 104039 (2015)
S.A. Hosseini Mansoori, B. Mirza, M. Davoudi Darareh, S. Janbaz, Int. J. Mod. Phys. A 31, 1650067 (2016)
X.H. Wu, R. Li, J.K. Zhao, Phys. Rev. D 93, 064008 (2016)
E. Witten, arXiv:0706.3359
S. Carlip, Class. Quantum Gravity 22, R 85 (2005)
O. Gurtug, S. Habib Mazharimousavi, M. Halilsoy, Phys. Rev. D 85, 104004 (2012)
S.H. Hendi, B. Eslam Panah, R. Saffari, Int. J. Mod. Phys. D 23, 1450088 (2014)
S.H. Hendi, B. Eslam Panah, S. Panahiyan, A. Sheykhi, Phys. Lett. B 767, 214 (2017)
M. Cataldo, A. Garcia, Phys. Lett. B 456, 28 (1999)
Y.S. Myung, Y.W. Kim, Y.J. Park, Phys. Rev. D 78, 044020 (2008)
S.H. Mazharimousavi, M. Halilsoy, Mod. Phys. Lett. A 30, 1550177 (2015)
S.H. Hendi, S. Panahiyan, B. Eslam Panah, Adv. High Energy Phys. 2015, 743086 (2015)
Z.Y. Tang, C.Y. Zhang, M.K. Zangeneh, B. Wang, J. Saavedra, Eur. Phys. J. C 77, 390 (2017)
M. Cataldo, A. Garcia, Phys. Lett. B 456, 28 (1999)
R. Yamazaki, D. Ida, Phys. Rev. D 64, 024009 (2001)
S.H. Hendi, JHEP 03, 065 (2012)
S.H. Hendi, B. Eslam Panah, S. Panahiyan, JHEP 05, 029 (2016)
S.H. Hendi, B. Eslam Panah, S. Panahiyan, Prog. Theor. Exp. Phys. 2016, 103A02 (2016)
S.H. Hendi, S. Panahiyan, S. Upadhyay, B. Eslam Panah, Phys. Rev. D 95, 084036 (2017)
K.C.K. Chan, R.B. Mann, Phys. Rev. D 50, 6385 (1994)
K.C.K. Chan, Phys. Rev. D 55, 3564 (1997)
S.H. Hendi, Eur. Phys. J. C 71, 1551 (2011)
S.H. Hendi, Prog. Theor. Phys. 127, 907 (2012)
Y. Hagihara, Jpn. J. Astron. Geophys. 8, 67 (1931)
C. Darwin, Proc. R. Soc. Lond. A 249, 180 (1959)
C. Darwin, Proc. R. Soc. Lond. A 263, 39 (1961)
E. Hackmann, C. Lammerzahl, Phys. Rev. D 78, 024035 (2008)
E. Hackmann, PhD thesis, Universitat Bremen (2010)
R. Fujita, W. Hikida, Class. Quantum Gravity 26, 135002 (2009)
E. Hackmann, H. Xu, Phys. Rev. D 87, 124030 (2013)
E. Hackmann, C. Lammerzahl, V. Kagramanova, J. Kunz, Phys. Rev. D 81, 044020 (2010)
E. Hackmann, V. Kagramanova, J. Kunz, C. Lammerzahl, Phys. Rev. D 78, 124018 (2008) [Erratum-ibid. 79, 029901 (2009)]
V.Z. Enolski, E. Hackmann, V. Kagramanova, J. Kunz, C. Lammerzahl, J. Geom. Phys. 61, 899 (2011)
S. Soroushfar, R. Saffari, A. Jafari, Phys. Rev. D 93, 104037 (2016)
S. Soroushfar, R. Saffari, E. Sahami, Phys. Rev. D 94, 024010 (2016)
B. Hoseini, R. Saffari, S. Soroushfar, J. Kunz, S. Granau, Phys. Rev. D 94, 044021 (2016)
S. Soroushfar, R. Saffari, J. Kunz, C. Lmmerzahl, Phys. Rev. D 92, 044010 (2015)
S. Soroushfar, R. Saffari, S. Kazempour, S. Grunau, J. Kunz, Phys. Rev. D 94, 024052 (2016)
G. Abbas, U. Sabiullah, Astrophys. Space Sci. 352, 769 (2014)
J.P. Hu, Y. Zhang, L.L. Shi, G.H. Guo, P.F. Duan, Mod. Phys. Lett. A 34, 1950182 (2019)
T. Maki, K. Shiraishi, Class. Quantum Gravity 11, 227 (1994)
S. Grunau, V. Kagramanova, J. Kunz, C. Lammerzahl, Phys. Rev. D 86, 104002 (2012)
S. Grunau, V. Kagramanova, J. Kunz, Phys. Rev. D 87, 044054 (2013)
G. Abbas, N. Yousaf, M. Zubair, R. Saleem, Int. J. Mod. Phys. A 34, 1950208 (2019)
S. Grunau, B. Khamesra, Phys. Rev. D 87, 124019 (2013)
E. Hackmann, B. Hartmann, C. Laemmerzahl, P. Sirimachan, Phys. Rev. D 81, 064016 (2010)
E. Hackmann, B. Hartmann, C. Lammerzahl, P. Sirimachan, Phys. Rev. D 82, 044024 (2010)
A. Acena, E. Lopez, F. Aldas, Circular geodesics stability in a static black hole in new massive gravity. Galaxies 8, 14 (2020)
K. Nakashi, S. Kobayashi, S. Ueda, H. Saida, Null geodesics and repulsive behavior of gravity in (2+ 1) D massive gravity. Prog. Theor. Exp. Phys. 2019, 073E02 (2019)
N. Cruz, C. Martinez, L. Pena, Geodesic structure of the (2+ 1)-dimensional BTZ black hole. Class. Quantum Gravity 11, 2731 (1994)
C. Farina, J. Gamboa, A.J. Segui-Santonja, Class. Quantum Gravity 10, L193 (1993)
N. Cruz, C. Martinez, L. Pena, Class. Quantum Gravity 11, 2731 (1994)
S.H. Hendi, S. Panahiyan, R. Mamasani, Gen. Relativ. Gravit. 47, 91 (2015)
S.B. Russ, Hist. Math. 7, 156 (1980)
S. Zhou, J. Chen, Y. Wang, Int. J. Mod. Phys. D 21, 1250077 (2012)
S. Grunau, H. Neumann, S. Reimers, Phys. Rev. D 97, 044011 (2018)
V.M. Buchstaber, V.Z. Enolskii, D.V. Leykin, Reviews in Mathematics and Mathematical Physics, vol. 10 (Gordon and Breach, London, 1997)
Acknowledgements
We would like to thank the anonymous referee for constructive comments. We also thank Shiraz University Research Council. BEP thanks National Elites Foundation of Iran. E.H. is grateful for support from the Research Training Group RTG-1620 Models of Gravity, and the Cluster of Excellence EXC-2123 QuantumFrontiers—390837967, both funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation).
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Funded by SCOAP3
About this article
Cite this article
Hendi, S.H., Tavakkoli, A.M., Panahiyan, S. et al. Simulation of geodesic trajectory of charged BTZ black holes in massive gravity. Eur. Phys. J. C 80, 524 (2020). https://doi.org/10.1140/epjc/s10052-020-8065-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-020-8065-9