Kerr-Sen-Taub-NUT spacetime and circular geodesics
Haryanto M. Siahaan

TL;DR
This paper introduces a new solution in heterotic string theory describing a rotating, charged spacetime with gravitomagnetic monopole moment, and analyzes circular geodesics within this geometry.
Contribution
It presents a novel Kerr-Sen-Taub-NUT solution with detailed conserved charges and investigates test particle circular motions in this spacetime.
Findings
New exact solution in heterotic string theory
Conserved charges characterized for the spacetime
Analysis of circular geodesics and particle motion
Abstract
We present a novel solution obeying classical equation of motion in the low energy limit of heterotic string theory. The solution represents a rotating mass with electric charge and gravitomagnetic monopole moment. The corresponding conserved charges are discussed. We also study the test body circular equatorial motions in the spacetime.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
**Kerr-Sen-Taub-NUT spacetime and circular geodesics
**
Haryanto M. Siahaan***[email protected]
Center for Theoretical Physics, Department of Physics,
Parahyangan Catholic University,
Jalan Ciumbuleuit 94, Bandung 40141, Indonesia
We present a novel solution obeying classical equation of motion in the low energy limit of heterotic string theory. The solution represents a rotating mass with electric charge and gravitomagnetic monopole moment. The corresponding conserved charges are discussed. We also study the test body circular equatorial motions in the spacetime.
1 Introduction
String theory offers a consistent quantum description of gravity and unification of all forces in nature. However, the lacking for experimental verifications yield the final answer for quantum gravity explanation is still opened. While waiting for proofs or insights from experimental aspects, theoretical works related to string theory constantly appear, including ones which make predictions from the low energy limit effective field theory.
In particular of our interest is the low energy limit effective field description of heterotic string theory. In this framework, Ashoke Sen has constructed a solution describing rotating and charged black holes, known as Kerr-Sen black holes, analogous to the Kerr-Newman black holes of Einstein-Maxwell theory. Despite the similarities between Kerr-Sen and Kerr-Newman solutions, there are several features which distinguish between the two solutions. For example, entropy of an extremal Kerr-Sen black hole can be expressed solely on its angular momentum [1], the lacking of -picture hidden conformal symmetry for the generic Kerr-Sen spacetime [2, 3], and its unique behavior of light deflection [4].
Subsequently Johnson and Myers discovered a static dyonic Taub-NUT solution of the low energy heterotic string [5], starting with the Taub-NUT solution to the vacuum Einstein equation as the seed metric. In their works, the authors of [1, 5] employed the Hassan-Sen [6] transformation that can map any stationary and axial symmetric spacetime which solves the vacuum Einstein equation to a new solution in the low energy limit of heterotic string theory. Recently, the author of [7] employ the same method to obtain the acccelerating spacetime in the low energy limit of heterotic string using the -metric as the seed solution. Interestingly, aspects of black holes in the low energy limit of heterotic string theory have been studied quite comprehensively in literature [8, 9, 10, 11, 12, 13, 14, 15, 16].
Taub-NUT solution in vacuum Einstein theory can be considered as an extension of Schwarzschild solution, where in addition to the mass there exist a continuous NUT parameter which nowadays quite common to be referred as the magnetic mass or the gravitomagnetic monopole moment. The Taub-NUT solution can also be generalized to have an electric charge and rotation, known as the Kerr-Newman-Taub-NUT (KNTN) solution that solves the Einstein-Maxwell equations [17]. The KNTN spacetime possesses the timelike and axial Killing symmetries, just like the Kerr or Kerr-Newman solutions. Recently, there appear a proposal where the gravitomagnetic monopole could be non-zero [18], hence the Kerr-Taub-NUT (KTN) or KTNT solutions may have some physical relevances.
Since KTN spacetime solves the vacuum Einstein equations and has the timelike and axial symmetries, one may wonder the outcome of using this solution as a seed in the Hassan-Sen transformation. We could expect the result is something similar to the KNTN in Einstein-Maxwell theory, which is true and then can be called the Kerr-Sen-Taub-NUT (KSTN) solution. Some features of KNTN spacetime have been studied quite extensively in literature[18, 19, 20, 21, 22, 21, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32], then investigating some features of KSTN spacetime can be potentially appealing. Some of our particular interests are the conserved charges associated to the solution and circular geodesics in KSTN spacetime.
The organization of this papers is as the followings. In section 2, we provide a quick review on the low energy heterotic string theory action and obtain the KSTN solution. In section 3, we obtain the electric charge, mass, and angular momentum in KSTN spacetime. Then, in section 4 we study the timelike and null circular geodesics in KSTN spacetime, followed by the accompanying effective potential plots in the next section. Finally, we present conclusions and discussions.
2 Twisting the Kerr-Taub-NUT spacetime
2.1 Symmetry of fields in low energy effective action
In this paper we consider the stationary and axially symmetric spacetimes, , with the coordinates222The coordinate in Boyer-Lindquist type. . A set of fields which solves the equations of motion derived from the four dimensional low energy effective action for heterotic string
[TABLE]
can been expressed in terms of the seed solution [7]
[TABLE]
The seed metric is stationary and axially symmetric solution to the vacuum Einstein equation. Indeed, the vacuum Einstein system is a special case of action (2.1), i.e. when all fields but the spacetime metric vanish. Therefore, the Hassan-Sen transformation which maps the vacuum Einstein metric to a set of fields can be simply considered as a mapping between solutions in the theory described by action (2.1).
Explicitly, the fields solutions can be expressed as the followings. In the string frame, the spacetime metric reads
[TABLE]
where . The components of gauge field are
[TABLE]
which give the field strength tensor in action (2.1) above. The dilaton field is
[TABLE]
and the non vanishing components of the second-rank antisymmetric tensor field are
[TABLE]
The third-rank tensor appearing in action (2.1) consists of the gauge field and tensor field , i.e.
[TABLE]
where is just a Chern-Simons term in the absence of .
All non-gravitational fields disappear as , and the metric (2.3) reduces to the original seed solution (2.2). The resulting metric (2.3) has the timelike and axial Killing vectors, and respectively, as in the case of seed metric. The conserved quantities related to these two Killing symmetries are given in section 3. In the following sections, the geodesic calculations will be performed in Einstein frame where the corresponding action reads [5]
[TABLE]
The relation between string and Einstein frames metrics is given by .
2.2 Kerr-Sen-Taub-NUT solution
In constructing a set of fields representing a massive rotating charged object equipped with NUT parameter in the low energy limit of heterotic string theory, namely the Kerr-Sen-Taub-NUT (KSTN) solution, we take the Kerr-NUT metric as the seed metric ,
[TABLE]
[TABLE]
where , , and . Taking in the last equation yields the Kerr metric. Spacetime (2.9) is stationary and axial symmetric, and it solves the vacuum Einstein equation. Using the prescription (2.3), we have
[TABLE]
[TABLE]
where
[TABLE]
According to (2.5) and (2.4), the corresponding dilaton and vector fields are
[TABLE]
and
[TABLE]
respectively. Finally, the non-zero components of are given by (2.6), i.e.
[TABLE]
In the next section, we obtain the conserved charges with respect to and symmetries, i.e. mass and angular momentum measured at infinity, respectively. For KSTN spacetime, the mass is found to be
[TABLE]
just like in the case of Kerr-Sen black holes [1]. The charge is also the same as that of Kerr-Sen,
The Einstein frame version of the metric (2.10) is the one that we use in the following section, which can be written as
[TABLE]
[TABLE]
where
[TABLE]
[TABLE]
and
[TABLE]
Note that we have expressed the equations above in terms of the mass and charge of Kerr-Sen black holes333In section 3 we obtain the conserved charges, including mass and electric charge, in Kerr-Sen-Taub-NUT geometry.. Alternatively, the non-gravitational fields can be expressed as
[TABLE]
[TABLE]
and
[TABLE]
The set of fields solutions (2.16) - (2.22) reduce to that of Kerr-Sen [1] after setting , and becomes the standard Kerr-NUT (2.9) in the absence of electric charge .
It has been well known that Taub-NUT or Kerr-Taub-NUT solutions suffer conical singularities, which according Misner can be removed by introducing a periodic timelike coordinate in the stationary spacetime444Clearly the conical singularities or periodic timelike coordinate are not the ones that we may expect to experience in reality. [17]. The seed metric (2.9) has the semi-infinite line singularity, where it is regular on the half axis , i.e.
[TABLE]
but conical defect occurs at . Therefore, it is not surprising to find that the metric (2.16) also suffers the same problem. Consequently, the conical defect yields a black hole is not a quite well defined object in the KSTN spacetime. Nevertheless, there exist quite a number of works in literature which keep using the black hole terminology in a spacetime with conical defect coming from the presence of NUT parameter [26, 33].
Regardless that we have a well defined description of black holes in KSTN spacetime or not, the celestial kinematics of objects in this geometry are still fascinating to be investigated since they might have some relevances to astronomical observations. In order to do this, we need to explore more the properties of this spacetime. Despite we cannot have an event horizon in normal sense living in this spacetime, one can still locate the radius which yields the metric (2.16) to be singular, that is given by . They are located at
[TABLE]
If we allow a black hole can exist in this Kerr-Sen-Taub-NUT spacetime555This situation is similar to that in Kerr-Newman-Taub-NUT spacetime, where there is a radius which yields the metric to be singular, and one could have a horizon which give exactly Kerr-Newman black hole’s horizon in the absence of NUT parameter., then the radius of that black hole horizon would be . From eq. (2.23), we can see the maximal rotation or extremality is achieved at , i.e. to maintain to be real valued. Violation of this bound leads to a production of the naked ring-singularity. Furthermore, for the KSTN spacetime, there also exist the static limit surface
[TABLE]
which is the outer radial solution to and widely known as the outer region of ergosphere.
The KSTN spacetime (2.16) is stationary and axial symmetric, symmetries inherited from the seed solution. Therefore, the corresponding Killing vectors are and , which hint the existence of two associated conserved quantities, namely the energy and angular momentum . As it is expected in a rotating spacetime, one can also show the dragging effect in KSTN spacetime. The angular velocity of a stationary observer with constant and is given by
[TABLE]
which at the horizon reduces to
[TABLE]
In the vanishing of the NUT parameter in the last equation, we have the angular velocity of the Kerr-Sen black hole horizon [1]
[TABLE]
The area covered by the sphere with radius can be computed as
[TABLE]
which at is just the area of Kerr-Sen black hole.
Just like the analogous solution in Einstein-Maxwell theory, i.e. KNTN spacetime, the squared of Riemann tensor or Kretschmann scalar for KSTN spacetime is not singular at . This property is due to the presence of NUT parameter . In a spacetime where black holes can live, the singular Kretschmann scalar at is interpreted as the physical singularity that exist in the center of a black hole. Although a mathematical symbolic manipulation programs like MAPLE or Mathematica can do the computation to get the Kretschmann scalar for a spacetime, its full expression can be lengthy. This is exactly the case for KSTN spacetime, which is the geometry under investigation in this paper. Nevertheless, to support our future claim related to the curve of effective potential, let us just show the Kretschmann scalar for KSTN spacetime evaluated at and on equatorial plane. The reading is
[TABLE]
where
[TABLE]
In expressions above, , and we prefer to use instead of for economical reason. Note that the squared of Riemann tensor evaluated at computed in (2.29) is finite for a non-vanishing .
3 Mass and angular momentum
Defining mass and angular momentum in a spacetime with NUT parameter is not quite simple. In an asymptotically flat spacetime, for example Kerr-Newman family without NUT parameter in Einstein-Maxwell theory, one can employ the Komar integral to get the mass and angular momentum. However, when the NUT parameter is present, one needs to carefully treat the surface integral at infinity to get the desired conserved charges associate to the and Killing vectors[23].
Under some general diffeomorphisms, we can employ the Barnich-Brandt method [34] to compute the conserved charges in any field theory, including those containing gravity. In Kerr/CFT correspondence, this method is used to obtain charges in the theory under consideration [35, 36, 37]. Virasoro algebra between these charges contains the central charge which is crucial in reproducing the Bekenstein-Hawking entropy for black holes in the spacetime whose asymptotic symmetries are associated to the obtained central charge. Clearly, the Barnich-Brandt method can also be applied for the exact diffeomorphisms and .
The low energy limit of heterotic string discussed in this paper contains four fields, where there exist diffeomorphisms associated to each of these fields. Consequently, each of these diffeomorphisms may contribute to the central charge according to Barnich-Brandt method. In this section, we compute the conserved charges associated to and Killing vectors of KSTN spacetime using the Barnich-Brandt method. As the starting point, let us reparameterize the mass as one parameter family , where . In such consideration, one can define the variation of each fields as follows666At the end of calculation, we can perform the integration for the obtained conserved charges.
[TABLE]
Accordingly, we have
[TABLE]
and
[TABLE]
The charge formula as a result of diffeomorphisms of each fields read
[TABLE]
where
[TABLE]
[TABLE]
[TABLE]
and
[TABLE]
[TABLE]
In equation above, where is the fields , and we have used
[TABLE]
as the two-surface orthonormal to the timelike and radial directions with . In writing equations (3.5) - (3.8), we have omitted the contributions from the gauge freedoms of and fields. The mass can be computed by taking in (3.4), which gives
[TABLE]
This is exactly the mass of Kerr-Sen black holes [1, 37], which can be computed using the standard Komar integral. In fact, the same mass can also be obtained using Abbott-Deser-Tekin method [38]. Furthermore, in calculation above, the conserved mass is given solely from the gravity contribution, i.e. integration.
On the other hand, employing the charge formula for the axial Killing vector yields
[TABLE]
As one would expect, the last formula matches the angular momentum of Kerr-Sen black hole if one sets . In other words, the angular momentum measured at infinity in KSTN spacetime gets contribution from the NUT parameter, which resembles the similar situation in KNTN case [39], and in the boosted Kerr-Taub-NUT study [23]. To compute the electric charge, one can use the formula
[TABLE]
which gives the same electric charge to that of Kerr-Sen solution [1]. Note that the electric charge above is obtained in the Einstein frame described by the action (2.8).
4 Equatorial circular geodesics
4.1 Generalities
In this section, we present some investigations on the equatorial and circular geodesics in the Kerr-Sen-Taub-NUT spacetime. As pointed out in [40], equatorial geodesics studies are not too straightforward for the spacetime with NUT parameter. In addition to the requirement and for the circular motion, we also need to verify that there is no velocity and acceleration in direction777Or in direction for the standard Boyer-Lindquist type coordinate. . In [27], the authors employ the Hamilton-Jacobi equation in Kerr-Newman-NUT spacetime to study and .
In this paper, we prefer to compute and by using the effective potential approach from the test body geodesics, just as in the same fashion we deal with the quantities and . As a start, we consider a general Lagrangian of a neutral test particle living in a curved background
[TABLE]
where the spacetime metric is a function of and . Constants of motion related to this Lagrangian are
[TABLE]
In the last equations, and are interpreted as energy and angular momentum of a test particle, respectively. Furthermore, a little algebra on the last equations yields
[TABLE]
Note that, in the last equation we have used a general relation that applies to Kerr-Sen-Taub-NUT spacetime, i.e. . If one considers an equatorial motion, then this relation reduces to the useful identity that normally employed in equatorial circular geodesics, namely .
In the followings, we look for some setups which allow circular geodesics on equatorial plane. From the metric that associate to the Lagrangian (4.1), we can have
[TABLE]
where . Alternatively, one can show
[TABLE]
after making use of (4.3). The metric signature that we are using yields for null case and for the timelike. Now let us first consider test body motions on equatorial plane, i.e. , hence eq. (4.4) now reads888Here means a function evaluated at some fixed .
[TABLE]
Accordingly, we can introduce an effective potential associated to (4.6), namely
[TABLE]
A circular motion denoted by requires this potential to be vanished. One small subtlety that needs attention here is to make sure the non-vanishing of
[TABLE]
hence the effective potential (4.7) cannot be singular. However, guided by the physical insight that it is very unlikely for a massive body to maintain a huge amount of electric charge, namely , we can safely claim the non-vanishing of (4.8).
Similarly, we can also define an effective potential that coming from a fixed radius consideration in (4.4). Instead of having (4.6) as a start, here we have
[TABLE]
that leads us to
[TABLE]
Provided that , which is guaranteed by our interest in studying circular geodesics with radii outside horizon, the vanishing on equatorial plane can be achieved for the particular that solves in eq. (4.7). Furthermore, to ensure that there is no acceleration in or direction, the following constraints
[TABLE]
and
[TABLE]
must be fulfilled as well. The last equation is normally satisfied in most spacetimes without NUT parameter, for example those in Einstein-Maxwell-dilaton [41] and low energy limit of heterotic string [42] theories. However, when NUT parameter exists, eq. (4.12) could add another constrain to the familiar ones obtained from . Particularly for timelike case, the stability of circular geodesics is dictated by the behavior of second derivative with respect to radius, and to find the so called innermost stable circular orbit (ISCO) we make use the condition
[TABLE]
In the followings, we employ the generalities that have been developed here to verify the existence of equatorial circular orbits in both timelike and null cases.
4.2 Circular timelike orbits
In timelike consideration, from eq. (4.12) we can have
[TABLE]
This equation is satisfied for zero NUT parameter, or the term inside the curly bracket vanish. The first case is the reason why equatorial geodesics in Kerr-Sen spacetime is guaranteed [14]. If we are interested in the non-vanishing of NUT parameter , then the energy and angular momentum of test particles get an additional constraint from the term in curly bracket above. Using this condition to constrain the angular momentum yields the reading of , , and equations can be expressed as
[TABLE]
[TABLE]
[TABLE]
where
[TABLE]
[TABLE]
[TABLE]
[TABLE]
[TABLE]
[TABLE]
[TABLE]
[TABLE]
[TABLE]
[TABLE]
[TABLE]
[TABLE]
[TABLE]
It is not easy to extract some qualitative results from these there equations, (4.15), (4.16), and (4.17). However, we can obtain some numerical plots from these three equations, which are presented in fig. 4.1. From these numerical results, we learn that not only the outcome for ISCO radius is smaller than the corresponding horizons, the associated NUT parameter is also purely imaginary, i.e. . The case of timelike circular geodesics in Kerr-Taub-NUT spacetime is represented by the plots for , which cannot take place on equatorial plane . This is in agreement to a result in [40], which states that the circular motion in spacetime with NUT parameter is not in equatorial plane. Furthermore, from the numerical evaluations of , and , the same conclusions can also be drawn, namely the timelike circular geodesics cannot take place on the plane if the NUT parameter is non-zero. More numerical evaluations of for particular would likely to possess the same features, leading us to the same conclusion where circular timelike geodesics cannot exist on equatorial plane in Kerr-Sen-Taub-NUT spacetime.
4.3 Light ring
Now let us turn our discussion to the null case. Similarly, eq. (4.12) in null consideration
[TABLE]
demands a constraint for if we let the NUT parameter to be non-zero. Solving an equation obtained from the term in curly bracket for and substituting the result to and give us
[TABLE]
and
[TABLE]
respectively. In equations above, and appearing in the metric (2.16). Here we consider non-zero to maintain the effective potentials and to be non-singular. Nevertheless, one cannot get a real valued non-zero solution for the NUT parameter from the last two equations. Therefore we can conclude that light ring does not exist on equatorial plane in Kerr-Sen-Taub-NUT spacetime.
5 Effective potentials
Also, quite frequently the complexity of full expression for effective potential does not allow us to easily extract the desired informations from the system. Hence, we could rely on the numerical plots for to help in giving insights for some particular questions, despite the answer cannot be general. For example, by looking at the one can tell whether a test particle can travel from infinity to the vicinity of black holes. In this section, we do not aim to extract all the information we can get from the plots, but rather to verify the non-singular value of Kretschmann scalar in KSTN spacetime at . We also like to learn the effect of non-vanishing NUT parameter to the profile of . Therefore, limiting ourself to the case of timelike test particle only would be adequate, and the corresponding effective potential is (4.7) with . In fig. 5.1, we can observe that the behavior of effective potential in the presence of NUT parameter is distinguishable from the vanishing counterparts, i.e. Kerr and Kerr-Sen. The effective potentials for Kerr and Kerr-Sen are singular at physical singularity , but it is not for KSTN case. This resembles the analogous property of Kerr-Newman-Taub-NUT spacetime, which also differs from the case where Taub-NUT parameter is absent, i.e. the family of Kerr-Newman black holes. The finiteness of effective potential for the case represented by fig. 5.1 of KSTN spacetime can be understood from the non-singular Kretschmann scalar (2.29) evaluated at . Nevertheless, the plots for each cases in fig. 5.1 start to coincide at the larger radius.
Now let us see what happens when the electric charge parameter increases, as shown in fig. 5.2. Here we observe that near the plots is steeper as decreases999For economical reason we use the starred notation for all quantities in consideration to denote their ratio to black hole mass, namely , , , , and as in fig. 4.1. . Alternatively, one can state that the gravitational attraction is stronger when vanishes. Note that our test timelike particle is electrically neutral, hence Coulomb interaction does not exist between the object that curves the spacetime and the test particle. This finding is similar to that of Kerr-Sen case, where the effective potential plots for a neutral test particle raise as one increases the electric charge of the black hole.
On the other hand, as shown in fig. 5.3, we find the plots is deeper for larger values of Taub-NUT parameter near the center of radius. In this region, the gravitational attraction is stronger as increases. Similar situations exist in the case of Einstein-Maxwell theory, for example in the plots presented in [20], but the lines start to coincide at larger radius. Another similarity between KSTN and KNTN cases for timelike particle geodesics is the plots of for different angular momentum . Near the center of coordinate , is steeper for lower angular momentum, but the curves start to intersect at some distance. Moreover, from fig. 5.4 we can learn that the curves for different angular momentum intersect at a point which is the positive root of quadratic equation101010One can also find an equation dictating the intersection points between curves in fig. 5.3, which turn out to be some quintic polynomials in . Hence, we do not think it is necessary to express the equation here.
[TABLE]
[TABLE]
In equation above, and are the values of under consideration.
6 Conclusion
In this paper, we obtain a solution describing rotating charged mass with NUT parameter in the low energy limit of heterotic string theory, namely the Kerr-Sen-Taub-NUT spacetime. This solution resembles the Kerr-Newman-Taub-NUT spacetime in Einstein-Maxwell theory. As we expect, there exist several similarities between the two solutions, for example the non-singular squared Kretschmann scalar at and the profile of plots as NUT parameter increases.
In section 3, we also study the conserved charges associated to the new solution. The mass and angular momentum are calculated using Barnich-Brandt method, where the obtained mass is in agreement to that computed using standard Komar integral. However, we find a contribution from the NUT parameter in the angular momentum, which resembles the cases of Kerr-Taub-NUT or Kerr-Newman-Taub-NUT spacetimes. The electric charge is the same to that of Kerr-Sen black holes.
Related to the equatorial circular geodesics in both timelike and null considerations given in section 4, we learn that such motions cannot exist in Kerr-Sen-Taub-NUT. This conclusion is obvious in the null geodesic case, but is taken based on several numerical evaluations in the timelike consideration. However, we are convinced that any numerical results similar to the plots in fig. 4.1 lead to the same conclusion for any possible numerical values of . Our results presented in section 4 agree to the prediction made in [40], namely circular geodesics in spacetime with NUT parameter do not exist in equatorial plane.
There are future potential projects that can be done related to Kerr-Sen-Taub-NUT solution presented in this paper. First is the non-equatorial circular motions in Kerr-Sen-Taub-NUT spacetime. This is interesting since it may have some astronomical interests related to the observation of black hole image [43]. Since in this work we point out that circular null geodesics cannot take place on equatorial plane in Kerr-Sen-Taub-NUT spacetime, one might wonder where they can exist. Also, here we have not obtained the separable Hamilton-Jacobi equation for test particles in Kerr-Sen-Taub-NUT spacetime. The constraints that yields and for the test particle geodesics can be derived from this equation, and in general are independent each other. We aim to include this project in our future work. Another interesting problem that can be worked out is the Banados-Silk-West effect [44, 26] in Kerr-Sen-Taub-NUT spacetime. Indeed, the latter project is related to the previous one, namely obtaining the separable Hamilton-Jacobi equation for test particles in the spacetime.
Acknowledgement
This research is supported by LPPM-UNPAR under the grant no. PL72018027. I thank my colleagues from physics dept. UNPAR for their encouragements.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] A. Sen, Phys. Rev. Lett. 69 , 1006 (1992)
- 2[2] A. M. Ghezelbash and H. M. Siahaan, Class. Quant. Grav. 30 (2013) 135005
- 3[3] D. Chen, P. Wang and H. Wu, Gen. Rel. Grav. 43 (2011) 181
- 4[4] G. N. Gyulchev and S. S. Yazadjiev, Phys. Rev. D 75 , 023006 (2007)
- 5[5] C. V. Johnson and R. C. Myers, Phys. Rev. D 50 , 6512 (1994)
- 6[6] S. F. Hassan and A. Sen, Nucl. Phys. B 375 , 103 (1992)
- 7[7] H. M. Siahaan, Phys. Lett. B 782 , 594 (2018)
- 8[8] K. Düztaş, Int. J. Mod. Phys. D 28 (2018) no.02, 1950044
