Kerr–Sen–Taub–NUT spacetime and circular geodesics

We present a 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, and the separability of Hamilton–Jacobi equation for a test body in the spacetime is also investigated. Some numerical results related to the circular motions on equatorial plane are presented, but there is none that supports the existence of such geodesics.


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 Q-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 Eina e-mail: haryanto.siahaan@unpar.ac.id (corresponding author) stein 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 C-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 l 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][21][22][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 Sect. 2, we provide a quick review on the low energy heterotic string theory action and obtain the KSTN solution. In Sect. 3, we obtain the electric charge, mass, and angular momentum in KSTN spacetime. Then, in Sect. 4 we study the timelike and null circular geodesics in KSTN spacetime, where some numerical plots related to these type geodesics are presented. Separability of Hamilton-Jacobi equation in the spacetime is investigated in Sect. 5, and effective potentials related to the radial motion from far away approaching the horizon are given in Sect. 6. Finally, we present conclusions and discussions.

Symmetry of fields in low energy effective action
In this paper we consider the stationary and axially symmetric spacetimes, g μν = g μν (r, x), with the coordinates 1 x μ = [t, r, x, φ]. A set of fields G μν , , A μ , B μν which solves the equations of motion derived from the four dimensional low energy effective action for heterotic string can been expressed in terms of the seed solution [7] ds 2 =g tt dt 2 + 2g tφ dtdφ +g φφ dφ 2 +g rr dr 2 +g x x dx 2 . (2. 2) The seed metricg μν 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 g μν to a set of fields G μν , , A μ , B μν can be simply considered as a mapping between solutions in the theory described by action (2.1). Explicitly, the fields solutions G μν , , A μ , B μν can be expressed as the followings. In the string frame, the spacetime metric reads 3) 1 The coordinate x = cos θ in Boyer-Lindquist type. where = 1 + (1 +g tt ) sinh 2 (α/2). The components of U (1) gauge field are which give the field strength tensor F μν = ∂ μ A ν − ∂ ν A μ in action (2.1) above. The dilaton field is = − ln , (2.5) and the non vanishing components of the second-rank antisymmetric tensor field are The third-rank tensor appearing in action (2.1) consists of the gauge field A μ and tensor field B μν , i.e.
where H αβμ is just a Chern-Simons term in the absence of B μν . All non-gravitational fields , A μ , B μν disappear as sinh (α/2) → 0, 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, ∂ t and ∂ ϕ respectively, as in the case of seed metric. The conserved quantities related to these two Killing symmetries are given in Sect. 3. In the following sections, the geodesic calculations will be performed in Einstein frame where the corresponding action reads [5] The relation between string and Einstein frames metrics is given by G μν = e g μν .

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 metricg μν , where r = r 2 − 2mr + a 2 − l 2 , x = 1 − x 2 , and ρ 2 = r 2 +(l + ax) 2 . Taking l → 0 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 According to (2.5) and (2.4), the corresponding dilaton and vector fields are and respectively. Finally, the non-zero components of B μν are given by (2.6), i.e. (2.14) The static limit a → 0 limit of the solutions (2.10) and (2.14) above were discovered in [5]. In the next section, we obtain the conserved charges with respect to ∂ t and ∂ φ symmetries, i.e. mass and angular momentum measured at infinity, respectively. For KSTN spacetime, the mass is found to be just like in the case of Kerr-Sen black holes [1]. The electric charge is also the same as that of Kerr-Sen, and we will represent the charge in b instead of Q, The Einstein frame version of the metric (2.10) is the one that we use in the following section, which can be written as Note that we have expressed the equations above in terms of the mass M and charge Q of Kerr-Sen black holes. 2 Alternatively, the non-gravitational fields can be expressed as and The set of fields solutions (2.16)-(2.22) reduce to that of Kerr-Sen [1] after setting l = 0, and becomes the standard Kerr-NUT (2.9) in the absence of electric charge Q. 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 spacetime 3 [17]. The seed metric (2.9) has the semi-infinite line singularity, where it is regular on the half axis x = 1, i.e.
but conical defect occurs at x = −1. 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˜ r = 0. They are located at (2.23) If we allow a black hole can exist in this Kerr-Sen-Taub-NUT spacetime, 4 then the radius of that black hole horizon would be r + . From Eq. (2.23), we can see the maximal rotation or extremality is achieved at a 2 = (M − b) 2 + l 2 , i.e. to maintain r ± 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 which is the outer radial solution to g tt = 0 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 ζ (t) ≡ ∂ t and ζ (φ) ≡ ∂ φ , which hint the existence of two associated conserved quantities, namely the energy E and angular momentum L. 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 r and x is given by , (2.25) which at the horizon reduces to . (2.26) In the vanishing of the NUT parameter l in the last equation, we have the angular velocity of the Kerr-Sen black hole horizon [1] KS The area covered by the sphere with radius r + can be computed as which at l → 0 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 r → 0. This property is due to the presence of NUT parameter l. In a spacetime where black holes can live, the singular Kretschmann scalar at r → 0 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 r = 0 and on equatorial plane. The reading is where In expressions above, m = M − b, and we prefer to use m instead of M − b for economical reason. Note that the squared of Riemann tensor R κλμν R κλμν evaluated at r = 0 computed in (2.29) is finite for a non-vanishing l.

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 ∂ t 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 cen-tral charge. Clearly, the Barnich-Brandt method can also be applied for the exact diffeomorphisms ∂ t 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 ∂ t and ∂ φ Killing vectors of KSTN spacetime using the Barnich-Brandt method. As the starting point, let us reparameterize the mass m as one parameter family m = sm, where s ∈ [0, 1]. In such consideration, one can define the variation of each fields as follows 5 Accordingly, we have and The charge formula as a result of diffeomorphisms of each fields read In equation above, k (ζ ; ) = k μν (ζ ; ) d S μν where is the fields {g, A, }, and we have used 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. k (ζ ;g) integration.
On the other hand, employing the charge formula for the axial Killing vector ζ (φ) yields As one would expect, the last formula matches the angular momentum of Kerr-Sen black hole if one sets l → 0. 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 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).

Generalities
In this section, we will show that innermost stable circular orbit (ISCO) does not exist on equatorial plane of Kerr-Sen-Taub-NUT spacetime, due to the more stringent constraint on that type of geodesics in the presence of NUT parameter.
As pointed out in [40], equatorial geodesics studies are not too straightforward for the spacetime with NUT parameter. In addition to the requirementṙ = 0 andr = 0 for the circular motion, we also need to verify that there is no velocity and acceleration in x direction. 6 In [27], the authors employ the separable Hamilton-Jacobi equation in Kerr-Newman-NUT spacetime to studyθ = 0 andθ = 0. However, the separability of Hamilton-Jacobi equation does not occur in Kerr-Sen-Taub-NUT background, as we show in the next section.
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 r . As a start, we consider a general Lagrangian of a neutral test particle living in a curved background where the spacetime metric is a function of r and x. Constants of motion related to this Lagrangian are − E = g ttṫ + g tφφ and L = g φφφ + g tφṫ .
In the last equations, E and L are interpreted as energy and angular momentum of a test particle, respectively. Furthermore, a little algebra on the last equations yieldṡ t x r = g tφ L + g φφ E and −φ x r = g tt L + g tφ E. (4.3) Note that, in the last equation we have used a general relation that applies to Kerr-Sen-Taub-NUT spacetime, i.e. g tt g φφ − g 2 tφ = x r . If one considers an equatorial motion, then this relation reduces to the useful identity that normally employed in equatorial circular geodesics, namely g tt g φφ − g 2 tφ = r . 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 where ≡ g rrṫ 2 + 2g tφṫφ + g φφφ 2 . Alternatively, one can show = g φφ E 2 + 2g tφ L E + g tt L 2 x r (4.5) after making use of (4.3). The metric signature that we are using yields δ = 0 for null case and δ = −1 for the timelike. 6 Or in θ direction for the standard Boyer-Lindquist type coordinate. Accordingly, we can introduce an effective potential V r,eff = −ṙ 2 associated to (4.6), namely A circular motion denoted byṙ = 0 requires this potential to be vanished. One small subtlety that needs attention here is to make sure the non-vanishing of 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 b < M, we can safely claim the non-vanishing of (4.8).
Similarly, we can also define an effective potential V x,eff = −ẋ 2 that coming from a fixed radius consideration in (4.4). Instead of having (4.6) as a start, here we have that leads us to (4.10) Provided that r | r =r 0 = 0, which is guaranteed by our interest in studying circular geodesics with radii r 0 outside horizon, the vanishing V x,eff r =r 0 on equatorial plane can be achieved for the particular r = r 0 (E, L , M, b, l) that solves V r,eff x=0 = 0 in Eq. (4.7). Furthermore, to ensure that there is no acceleration in r or x direction, the following constraints and dV x,eff dx = 0, (4.12) 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 V r .

Timelike case
In timelike consideration, from Eq. (4.12) we can have 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 L, 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 L yields the reading of V r,eff = 0, V r,eff = 0, and V r,eff = 0 equations can be expressed as 14) It is not easy to extract some qualitative results from these there equations, (4.14)-(4.16), to claim that the timelike equatorial circular motion can occur or the opposite. However, some numerical plots can be evaluated based on We have used a * to denote the ratio a/M. The peculiar behavior of larger radius as a * increases is understood since the events depicted in the plots is not physical these three equations, which we present in Fig. 1. From these numerical results, we learn that not only the expected circular radii are smaller than the corresponding horizons, the associated NUT parameters l are also purely imaginary. The case of timelike circular geodesics in Kerr-Taub-NUT spacetime is represented by the plots for b = 0, which cannot take place on equatorial plane x = 0. This is in agreement to a finding in [40], which states that the circular motion in spacetime with NUT parameter is not in equatorial plane. Furthermore, from the numerical evaluations of b = 0.05 M, and b = 0.1 M, the same conclusions can also be drawn, namely the timelike circular geodesics cannot take place on the plane x = 0 if the NUT parameter is non-zero. More numerical evaluations of for particular b 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.

Null case
Now let us turn our discussion to the null case. Similarly, Eq. (M{5r + 6b} − 3{r + b}{r + 2b})} = 0, (4.19) respectively. In equations above, 0 ≡ (x = 0) and 0 ≡ (x = 0) appearing in the metric (2.16). Here we consider non-zero 0 to maintain the effective potentials V x,eff and V r,eff to be non-singular. Nevertheless, one cannot get a real valued non-zero solution for the NUT parameter l from the last two equations. Therefore we can conclude that the null circular orbit which leads to the so called light ring does not exist on equatorial plane in Kerr-Sen-Taub-NUT spacetime.

Hamilton-Jacobi equation and separability
As an alternative to obtain the geodesic equations discussed in Sect. 4, one can employ the Hamilton-Jacobi equation (HJE) [43] g μν ∂ μ S∂ ν S = −m 2 . (5.1) For a spacetime with stationary and axial Killing vectors, the usual action ansatz is If Eq. (5.1) is separable after employing the ansatz (5.2), this equation splits into two independent equations, namely the radial and angular parts. These two parts which come from (5.1) can later give the effective potentials V r,eff and V x,eff studied in Sect. 4. However, not all stationary and axial symmetric spacetime enjoy the separability for HJE, in massive or massless cases. The work presented in [44] studied the separability of black holes in string theory, and one of the subject is separability of HJE. Spacetimes which possess the Conformal Killing Stackel Tensor (CKST) is guaranteed to have separability in the corresponding massless HJE, and those which exhibit Exact Killing Stackel Tensor (EKST) can be shown to have separability in massive HJE.
To see whether the spacetime (2.16) has CKST, let us rewrite the metric in the general form considered in [44], namely As shown in [44], if the metric can be expressed in terms of CKSTs X μν (x) and R μν (r ), namely where and It is obvious some terms in r.h.s. of Eq. (5.6) are already in the separable form between the functions that depend on r and x, except the last one that couples to ∂ 2 t . The existence of as a denominator in that term causes the problem, exactly similar to the issue for HJE separability of black hole spacetime in Kaluza-Klein theory [45]. Consequently, (5.6) cannot take the separable form as aimed in (5.5), which then allow us to conclude that massless HJE in Kerr-Sen-Taub-NUT background (2.16) is not separable. Furthermore, the failure of constructing CKSTs X μν (x) and R μν (r ) for the background (2.16) forbids the possibility to build the EKSTs This implies that the massive HJE in for the spacetime (2.16) cannot be separable as well. Nevertheless, the subjects of discussions in [44] are asymptotically flat black holes in string theory, which is definitely not the case of spacetime solution studied in this paper. However, the analysis should apply to non-asymptotic spacetime as well, provided that the spacetime possesses the stationary and axial Killing symmetries. Note that the existence of CKSTs and EKSTs is more obvious if the function A in metric (5.3) is r dependent only. This condition is satisfied in the metric (5.3) if NUT parameter l is null, and the term that couples to ∂ 2 t in (5.6) becomes separable as shown in [46]. In fact, it is the existence of NUT parameter which costs the spacetime (2.16) to lose the asymptotic flatness.
Obviously, the existence of NUT parameter in a spacetime does not imply the lost of separability of HJE that belongs to that spacetime. A familiar example is the separable massive HJE in Kerr-Newman-Taub-NUT spacetime [26]. An alternative check whether HJE is separable or not in a particular stationary and axial symmetric background is whether the corresponding metric can be expressed in canonical form or not [43]. The canonical form of the metric is where C x , D x , and Z x are x dependent functions, and C r , D r , and Z r are r dependent. Incompatibility of solution (2.16) to be expressed in the canonical form (5.10) could also be considered as the origin in separability problem of HJE in the background (2.16). Hence, studies presented in this section lead us to the conclusion that separable HJE in both massless and massive case is not the feature of Kerr-Sen-Taub-NUT spacetime.
It is interesting to find that the separability of massive HJE in the rotating and charged black hole with NUT parameter in Einstein-Maxwell theory, i.e. Kerr-Newman-Taub-NUT spacetime [26], does not repeat here for the analogous case in the low energy limit of heterotic string theory. Despite the nonexistence of separability in HJE does not hinder the study of test body geodesics in the spacetime, the absence of this separability would mean that the Carter constant does not exist. Note that the seed solution (2.9) and also the null NUT case of (2.16) known as the Kerr-Sen spacetime both possess the separability of massive HJE in the corresponding background. In other words, the limit Q → 0 or l → 0 case of (2.16) enjoys the separability of HJE. This separability disappears when both parameters, namely the charge Q and NUT parameter l, present at the same time. We may argue that this lost of separability occurs due to the twisting procedure that we employ in transforming the seed solution (2.9) to get (2.16), which acts on the (t, t), (t, φ), and (φ, φ) sectors of the metric while g rr and g x x remain invariant. This is obvious in the string frame (2.10) where the conformal factor e − has not been applied to get the Einstein frame version of the metric. The NUT parameter gives rise to some extra terms in g t,t , g t,φ , and g φ,φ in the seed solution which after Hassan-Sen twisting procedure in Sect. 2 yields some non-separable terms in (5.1).

Effective potentials
Also, quite frequently the complexity of full expression for effective potential V r,eff does not allow us to easily extract the desired informations from the system. Hence, we could rely on the numerical plots for V r,eff to help in giving insights for some particular questions, despite the answer cannot be general. For example, by looking at the V r,eff 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 V r,eff plots, but rather to verify the non-singular value of Kretschmann scalar in KSTN spacetime at r = 0. We also like to learn the effect of non-vanishing NUT parameter to the profile of V r,eff . Therefore, limiting ourself to the case of timelike test particle only would be adequate, and the corresponding effective potential is (4.7) with δ = −1. In Fig. 2, we can observe that the behavior of effective potential in the presence of NUT parameter l is distinguishable from the vanishing l counterparts, i.e. Kerr and Kerr-Sen. The effective potentials for Kerr and Kerr-Sen are singular at physical singularity r = 0, 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. 2 of KSTN spacetime can be understood from the non-singular Kretschmann scalar (2.29) evaluated at r = 0. Nevertheless, the plots V r,eff for each cases in Fig. 2 start to coincide at the larger radius. Now let us see what happens when the electric charge parameter b increases, as shown in Fig. 3. Here we observe that near r = 0 the plots is steeper as b * decreases. 8 Alternatively, one can state that the gravitational attraction is stronger when b 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. 4, we find the plots is deeper for larger values of Taub-NUT parameter l near the center of radius. In this region, the gravitational attraction is stronger as l 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 V r,eff for different angular momentum L. Near the center of coordinate r = 0, V r,eff is steeper for lower angular momentum, but the curves start to intersect at some distance. Moreover, from Fig. 5 we can learn that the V r,eff curves for different angular momentum intersect at a point which is the positive root of quadratic r * equation 9 In equation above, p and q are the values of L * under consideration. 9 One can also find an equation dictating the intersection points between curves in Fig. 4, which turn out to be some quintic polynomials in r * . Hence, we do not think it is necessary to express the equation here.

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 r = 0 and the profile of V eff plots as NUT parameter l increases. In Sect. 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 l 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 Sect. 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. 1 lead to the same conclusion for any possible numerical values of b. Our results presented in Sect. 4 agree to the prediction made in [40], namely circular geodesics in spacetime with NUT parameter are not in equatorial plane. Geodesics of test body either massive or massless are discussed without using some Carter constants, since the absence of Hamilton-Jacobi equation separability in Kerr-Sen-Taub-NUT spacetime as showed in Sect. 5.
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, which can be some particular interests due to to the observation of black hole image [47]. Another interesting problem that can be worked out is the Banados-Silk-West effect [26,48] in Kerr-Sen-Taub-NUT spacetime.
Acknowledgements This research is supported by LPPM-UNPAR under the Grant no. PL72018027. I thank my colleagues from physics dept. UNPAR for their encouragements.

Data Availability Statement
This manuscript has associated data in a data repository. [Authors' comment: This is a theoretical study and there is no associated experimental data in the article.] 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://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .