Probing a black-bounce-Reissner–Nordström spacetime with precessing and periodic motion

We investigate the precessing and periodic orbits of a test timelike particle around the black-bounce-Reissner–Nordström spacetime which is characterized by its charge and bounce parameter. Its marginally bound orbit and innermost stable circular orbit are obtained in the exact forms. We pay closely attention to its precessing orbits and find the resulting relativistic periastron advance. We also study its periodic orbits and demonstrate that small variations of the charge and bounce parameter can make the motion jump among the periodic and precessing orbits. In these two kinds of orbits, we find a distinct degeneracy that some specific combinations of the charge and bounce parameters can generate exactly the same orbital motion in the black-bounce-Reissner–Nordström spacetime, which can also mimic those of the Schwarzschild black hole. In order to break such a degeneracy, we make use of the precession of S2 star around Sgr A* detected by GRAVITY together with the shadow diameter of Sgr A* measured by Event Horizon Telescope and find preliminary bounds on the charge and bounce parameter.


Introduction
The detection of gravitational waves [1][2][3][4][5][6] and the direct images of the supermassive black holes in the centers of galaxy M87 [7][8][9][10][11][12] and of the Galaxy [13][14][15][16][17][18] reveal that black holes, predicted by the Einstein's general relativity more than one hundred years ago, are very common in the Universe. As an ideal laboratory for examining and testing the theories of gravity in the strong fields, a black hole might be poisoned by its event horizon that would cause the information-loss problem, and by its central singularity that would break the general relativity down. In order to erase the a e-mail: yixie@pmo.ac.cn (corresponding author) singularity, a number of ways have been proposed, such as bouncing by quantum pressure [19][20][21], building a regular core [22][23][24][25], and forming a quasi-black hole [26][27][28][29] (see Ref. [30] for a review).
In the past few years, the family of black-bounce spacetimes has been paid much attention (see Ref. [31] for a brief review). They are built under the general relativity, globally free from curvature singularities, and smoothly interpolate between regular black holes and traversable wormholes. The simplest black-bounce spacetime is the black-bounce-Schwarzschild spacetime [32], which can reduce to the Schwarzschild spacetime and the Ellis wormhole [33] as its bounce parameter and mass vanish, respectively. Its gravitational properties and astrophysical features have been intensively investigated [34][35][36][37][38][39][40][41], and it has been extended to several new classes [42][43][44].
In this work, we focus on the black-bounce-Reissner-Nordström spacetime [45], which is a charged extension of black-bounce-Schwarzschild spacetime and can be interpreted as standard Maxwell electromagnetism with an anisotropic fluid. It is well accepted that astrophysical objects must be without electric charges because of quickly neutralizing by environmental plasma. But the sounds that a black hole might be charged by some mechanisms, such as accumulating charged matter, induction through the rotation in the external magnetic field [46] and inheriting from its charged collapsed progenitor [47], can also be heard. It is likely that the supermassive black hole Sgr A* in the Galactic center might have transient and small positive charge [48], implying studies on the charged spacetimes may go beyond the territory of theoretical interests. Although we have investigated the gravitational-lensing properties of the black-bounce-Reissner-Nordström spacetime, it turns out that current observations might constrain its charge only but leave its bounce parameter untouched [49]. Therefore, in the present work, we will consider the bound orbits around the black-bounce-Reissner-Nordström spacetime, especially the precessing and periodic motion, hoping to provide more clues for testing such a spacetime.
As another subclass of the bound orbits, the periodic orbits will emerge when a timelike particle is in the vicinity of a black hole, showing zoom-whirl patterns [85][86][87][88]. This strong-field feature can be indicated by the ratio of the average angular frequency to the radial frequency per radial cycle [89]. When the ratio is irrational, the particle will depict a precessing orbit; otherwise, it will go back to the initial location exactly in the finite time and trace out a perfectly periodic orbit. Since the rational numbers are dense on the real number domain, a generic precessing orbit with an irrational ratio can be sufficiently approximated by a close-by periodic one. Though the zoom-whirl pattern of a periodic orbits has yet been observed, the periodic orbits might be helpful for faster computation of adiabatic extreme mass-ratio inspirals [90] and for providing information about strong-field properties of the spacetime that is unavailable from the precessing orbits. Therefore, the periodic orbits around a Schwarzschild black hole [89], a Kerr black hole [90][91][92][93][94] and other black holes [95][96][97][98][99][100][101][102][103][104][105][106], as well as those in binary black holes [107,108] have been intensively examined.
Triggered by these theoretical and observational progresses and in order to further explore the black-bounce-Reissner-Nordström spacetime, we will intensively investigate the precessing and periodic orbits around it in this work, providing more information about its signatures in the geodesic motion. The outline of this paper is as follows. In Sect. 2, we briefly review the black-bounce-Reissner-Nordström spacetime and investigate the bound orbits, including the marginally bound orbit and innermost stable circular orbit of a timelike particle around it. We focus on its precessing motion and estimate the bounds on the charge and bounce parameter based on the observations of GRAVITY and Event Horizon Telescope (EHT) in Sect. 3. In Sect. 4, we examine the periodic orbits around it and pay close attention to their transition with respect to small variations of the charge and the bounce parameter. Conclusions and discussion are presented in Sect. 5.

Metric
The black-bounce-Reissner-Nordström spacetime can be regarded as the solution to the electrovac Einstein field equations of the general relativity in the curvature coordinates (t, r, θ, ϕ) with an interpretation as standard Maxwell electromagnetism with an anisotropic fluid [45]. It can be obtained with the Reissner-Nordström spacetime by keeping dr in the metric tensor unchanged and replacing the radial coordinate r with r 2 + α 2 • , where α • is the bounce parameter with the dimension of length. Thus, the metric of the black-bounce-Reissner-Nordström spacetime with mass m • and the charge Q • can be written as (G = c = 1) [45] It is asymptotically flat as |r | → ±∞ and is globally regular due to the existence of the bounce parameter. The general relativity is thought to be valid at least sufficiently far away from its core region [45]. When α • = 0, the black-bounce-Reissner-Nordström spacetime goes back to the Reissner-Nordström spacetime. As the charge vanishes, i.e., Q • = 0, one recovers the black-bounce-Schwarzschild spacetime [32]. When m • = 0 and Q • = 0, the spacetime (1) arrives at the Ellis wormhole [33]. As Q • = 0 and α • = 0, one returns to the Schwarzschild black hole. The radius of the outer and inner event horizons for the spacetime (1) are [45] respectively, which requires for the existence of real r ± H . For later convenience, we define the following dimensionless quantities as and the rescaled bounce parameter as The condition for the existence of the outer event horizon (5) can be rewritten as and the dimensionless radius of the outer event horizon can be defined as As our first step in the study of geodesic motions around the black-bounce-Reissner-Nordström spacetime, we focus on such a spacetime with its event horizon(s) in this work and adopt the domain (8) in the following sections. The geodesic motions around the horizonless one will be left in our next moves.

Bound orbits
For a timelike test particle moving freely around the blackbounce-Reissner-Nordström spacetime, its Lagrangian can be expressed as where the particle is assumed to be confined in the equatorial plane θ = π/2 and a dot represents the differentiation with respect to the proper time. We can obtain two motion constants as and the equation of radial motion aṡ with the effective potential defined as where the relation A(r )B(r ) = 1 is used. Based on Eq. (6), the effective potential can be rewritten as If a timelike particle falls to the black-bounce-Reissner-Nordström spacetime freely from r = ∞ andṙ = 0, it would be finally captured into an unstable circular orbit around the black hole withṙ =r = 0 at r = r mb which is called marginally bound orbit, also known as zero binding energy zoom-whirl orbit. Such a marginally bound orbit satisfies that where denotes the derivative against r . From Eq. (17), we can obtain the dimensionless angular momentum and the dimensionless radius of the marginally bound orbit as (see Appendix A for details) where with Due to the sole dependence of y mb on the charge q, l mb is unaffected by the dimensionless bounce parameter α and it decreases with the increment of q, as seen in Fig. 1a. Based on Fig. 1b, we find that x mb decreases with the growth of q and λ. Thus, for a given (q, λ) ∈ D H , the ranges of the marginally bound orbit are It also suggests that the marginally bound orbit in the blackbounce-Reissner-Nordström spacetime could be closer to the central object than the one of the Schwarzschild spacetime and than the one of the Reissner-Nordström black hole for a given q. As the charge of the spacetime (1) vanishes, i.e., q = 0, l mb and x mb go back to those of the black-bounce-Schwarzschild spacetime [37]. When both q and α disappear, Another special case for the geodesic motion of a timelike particle is the stable circular motion with a minimal radius around the black-bounce-Reissner-Nordström spacetime which is known as innermost stable circular orbit. It is important in astrophysics because it is usually regarded as the inner edge of an accretion disk around the black hole and plays as a transition from inspiral to plunge for a binary merger. The innermost stable circular orbit emerges when the maximum and minimum points of the effective potential coincide by demanding that After some calculations (see Appendix B for details), we have where with Since y isco only depends on the charge q, the bounce parameter α does not affect the E isco and l isco . As seen in Fig. 2a and b, both of them decrease with the increment of q, while x isco shrinks with the growth of q and λ, see Fig. 2c. Thus, the ranges of the innermost stable circular orbit on the domain D H can be found as It shows that the innermost stable circular orbit in the blackbounce-Reissner-Nordström spacetime could be closer to the central object than the one of the Schwarzschild spacetime and than the one of the Reissner-Nordström black hole for a given q. The innermost stable circular orbit returns to the one of the black-bounce-Schwarzschild spacetime as q = 0 [37]. When q = λ = 0, it goes back to the one of the Schwarzschild spacetime with E isco,Sch = 2 √ 2/3, l isco,Sch = 2 √ 3 and x isco,Sch = 6. For a timelike test particle with given E and l in the blackbounce-Reissner-Nordström spacetime, its bound orbit can be described by a unique number p as [89] where ϕ is the accumulated azimuth between successive periastron during a radial period and can be expressed as where Eqs. (12) and (13) have been used and r ± are two turning points of the bound orbit, i.e., the roots ofṙ 2 = 0. When p is a rational number, the test particle will move in Fig. 2 The energy E isco , the dimensionless angular momentum l isco and the dimensionless radius x isco of the innermost stable circular orbit for the black-bounce-Reissner-Nordström spacetime are shown on the domain D H where E isco and l isco are independent on λ a periodic orbit and go back to its initial location exactly after a finite time. Otherwise, it runs a precessing orbit with an irrational number p. The precession per revolution of a timelike particle can be obtained as After having these essentials about the bound motion of a timelike particle around the black-bounce-Reissner-Nordström spacetime, we will focus on the precessing and periodic orbits in the following sections.

Precessing orbits
When p is an irrational number, the timelike particle will run a precessing orbit. As a particular subclass of the bound orbits, precessing orbits are usually used to test the theories of gravitation. In order to evaluate the periastron advance of the particle around the black-bounce-Reissner-Nordström spacetime and given the relation between the standard coordinate r s and the radial coordinate r that r 2 s = r 2 + α 2 • , we parametrize its orbit as 1 where the standard coordinate r s has a well-known parametrization as with a being the semi-major axis, e being the eccentricity and χ being the relativistic true anomaly. Then, the periastron r − and apoastron r + can be obtained for χ = 0 and π , respectively, which are With the fact that the radial velocity at the turning points r ± is zero, i.e.ṙ 2 ± = 0, the energy and angular momentum can be determined by Eq. (13) as Thus, Eq. (31) can be rewritten as where with Although the precession of a timelike particle around the black-bounce-Reissner-Nordström spacetime can be calculated numerically based on Eq. (38) with the given E and L (or r ± ), it would be helpful to obtain it analytic approximation, which might provide some intuitive information about how the charge and bounce parameter affect the orbital precession.
Treating the effects of the charge and the bounce parameter as perturbation and keeping their leading terms, we can obtain the precession of a timelike test particle around the blackbounce-Reissner-Nordström spacetime as Here, the first term is the well-known post-Newtonian precession in the Schwarzschild spacetime, and the second and third ones are caused by the charge and the bounce parameter, respectively. When α • = 0, the precession ω returns to the one of the Reissner-Nordström spacetime [109][110][111][112]. As Q • = 0, it goes back to the one of the black-bounce-Schwarzschild spacetime [37].
We can clearly see that the contributions of the charge Q • and the bounce parameter α • in the precession ω have opposite signs, leading to opposite precessing directions. Therefore, even though both of them have some quite big values, they would cancel out in the precession as long as  Fig. 3 where and For S2 star, we adopt its eccentricity e = 0.885 and semimajor axis a = 1031 au and mark the inferred values of f SP provided by GRAVITY's observations on it as [84] f SP = 1.10 ± 0.19.
The black, red and blue dashed lines represent the best-fit estimation of, the lower and the upper limits of f SP , respectively, while the dash-dotted line indicates the line of f SP = 1. From Fig. 3, the measured S2's precession might rule out two regions of the parameter space for the black-bounce-Reissner-Nordström spacetime: one with α • 100 m • and the other with Q • ≈ m • . However, it is unable to separate Q • and α • by the single measurement of S2 star. It demands measurements on the precession of more S stars or other kinds of observables to break the degeneracy of Q • and α • . Recently, directly imaging Sgr A* by EHT [13][14][15][16][17][18] might provide such an independent observable. We showed in Ref. [49] that the size of the photon sphere for the blackbounce-Reissner-Nordström spacetime is determined by the charge Q • and it is immune to the bounce parameter α • . The size of the shadow cast by such a photon sphere also relies on the mass and distance of Sgr A*, which was adopted from GRAVITY's measurement [113] for self-consistently joint analysis with the precession of S2 measured by GRAVITY [84]. Following the approach of the Ref. [18], we find the bound on Q • according to the measured diameter of Sgr A*'s shadow and then obtain the range of α • based on the precession of S2. These results are which are consistent with the one of the Reissner-Nordström spacetime (α • = 0) [112] and the one of the black-bounce-Schwarzschild spacetime (Q • = 0) [37]. Such a big value of α • does not suggest that Sgr A* is a traversable wormhole, but indicates the insufficient precision of current measurements. In estimation of the range of Q • , we did not take the spin of the black-bounce-Reissner-Nordström spacetime into account, since it is well known that the shape and size of the shadow have a very weak dependence on the spin and its inclination [114]. We expect it is valid at least at the leading order and we will leave this subtle issue as our next move. In estimation of α • , we used the best-fit values of a and e of S2 provided by GRAVITY [84], which were determined with assuming Sgr A* as a Schwarzschild black hole. It would make the range of α • be overestimated due to the lack of the correlations among the orbital parameters and Fig. 3 The color-indexed f SP for the black-bounce-Reissner-Nordström spacetime against the charge Q • and bounce parameter α • is shown. The black, red and blue dashed lines represent the best-fit estimation of, the lower and the upper limit of S2 star's f SP , respectively, while the dash-dotted line marks the line of f SP = 1 for the Schwarzschild precession the model parameters in our approach. Given these facts, we suggest that one should take our results (46) as preliminary bounds on Q • and α • and they mainly manifest the potential of multiple kinds of observations to break the degeneracy of parameters in constraining the properties of a spacetime.

Periodic orbits
Another specific subclass of the bound orbits is the periodic orbit which possesses a rational number p. Since the distribution of the rational numbers are dense on the real number domain, a generic orbit with an irrational p can be well approximated by a nearby periodic one. A periodic orbit can be specified by three integers (z, w, v) as [89] where z is the zoom number representing the number of closed leaf in the orbit, w is the whirl number representing the number of nearly circular whirls close to periastron per leaf, and v is the vertex number indicating the order in which the z leaves are traced out. Figure 4 shows an example of the variation of p with respect of q and λ under the energy E = 0.966 and the angular momentum l = 4. We can see that p increases slowly with the growth of λ and decreases very fast with q. Moreover, Fig. 4 clearly demonstrates the degeneracy of p between the charge q and the bounce parameter λ. It means that given E and l, two different sets of (q, λ) can generate the same orbit with the same p. Similar behavior in the orbital precession was also discussed in Sect. 3. In order to reveal the degeneracy for the periodic orbits, Fig. 5 shows some examples of the bound orbits around the black-bounce-Reissner-Nordström spacetime with various E but a fixed l = 3.98 where (X, Y ) are defined as 2 and it is easy to compare our results with those for other metrics in the standard coordinates. In each row, all orbits share the same E, and both q and λ for each panel increase from 0 to some specific values. The periodic orbits in the first, third and fifth columns have the same shape with the same p, although they have very different (q, λ). It means that the periodic orbits in the third and fifth columns can mimic those around the Schwarzschild black hole in the first column (q = λ = 0). For (q, λ) between those of the columns with the periodic orbits, the orbits might either become precessing orbits in most cases, or transition to another periodic orbit with different p, such as the cases in panel Fig. 5(4a) and (4b) where p jumps from 5/6 to 4/5 after small changes of q and λ. Figure 6 also shows some examples of the bound orbits around the black-bounce-Reissner-Nordström spacetime. It is similar to Fig. 5 except that its each row shares the same l under the fixed energy E = 0.986. Likewise, the periodic orbits in the first, third and fifth columns have the same p and shape but with different (q, λ). It also demonstrates that, small variations of q and λ might not only make bound motions change from the periodic orbits to the precessing ones, but also cause them to transition between the periodic ones with different p, such as those in the panels of (4a)→(4b)→(4c) and of (5a)→(5b)→(5c) in Fig. 6. In summary, there is a degeneracy of bound motion between q and λ, whose small changes can also make the motion jump among the periodic orbits and precessing ones.

Conclusions and discussion
In this work, we investigate the bound motion of a timelike test particle around the black-bounce-Reissner-Nordström spacetime. We obtain its marginally bound orbit and innermost stable circular orbit in the exact forms, whose energy and angular momentum depend on the charge Q • only and whose radii are affected by both Q • and the bounce parameter α • . Then, we pay closely attention to its precessing orbits and find the resulting relativistic periastron advance. We also study its periodic orbits and demonstrate that small variations of Q • and α • can make the motion jump among the periodic orbits and precessing ones. In both precessing and periodic orbits, we find a distinct degeneracy so that orbital motion with non-vanishing Q • and α • can mimic those of the Schwarzschild black hole. In order to test the black-bounce-Reissner-Nordström spacetime and constrain its parameters, it requires other independent measurements to break such a degeneracy. Therefore, based on the precession of S2 star around Sgr A* detected by GRAVITY [84] together with the shadow diameter of Sgr A* measured by EHT [13][14][15][16][17][18], we find preliminary bounds on Q • and α • , see Eq. (46).
The black-bounce-Reissner-Nordström spacetime we considered is without rotation, while celestial objects are rotating in the real Universe. Although the nonrotational approximation might be adequate for bound motion sufficiently far from it or for very slowly rotating cases, one has to take its spin into account for those much close to its center, bringing more complicated and new features. As the spinning version of a charged black bounce spacetime, the black-bounce-Kerr-Newman spacetime is known, but such a solution requires non-trivial matter content [45]. Other effects, such as the gravitational radiation and its backreaction on the motion of the particle, will be interesting topics for testing this spacetime with gravitational waves in the future. These sophisticated studies might demand dedicated numerical methods [115][116][117][118][119][120] to investigate the dynamics of particles which would be necessary for extracting meaningful physical information from their motion.
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 . SCOAP 3 supports the goals of the International Year of Basic Sciences for Sustainable Development.

Appendix A: Determination of marginally bound orbit
With the dimensionless quantities x = m −1 • r and l = m −1 • L, the conditions for the marginally bound orbit (17) can be rewritten as From Eq. (A.1), we can obtain l 2 as Considering the timelike geodesics must stay outside the event horizon, i.e., x > x H and A(x) > 0, we can find where Eqs. (9) and (7)  Due to exclusive dependence of y mb on q, Eq. (A.12) makes l mb immune to the dimensionless bounce parameter α. As q = 1 and q = 0 respectively, l mb arrive at its smallest and biggest values that where Eq. (A.11) has been used. Using the Eq. (A.6) and (A.9), we can find the dimensionless radius of the marginally bound orbit as where α is replaced with the rescaled bounce parameter λ by Eq. (7). On D H , x mb reaches its the biggest value at q = 0 and λ = 0 and its smallest one at q = 1 and λ = 1, respectively, that As q = 0, l mb and x mb return to those of the black-bounce-Schwarzschild spacetime [37]. When both q and λ disappear, they go back to those of the Schwarzschild spacetime with l mb,Sch = 4 and x mb,Sch = 4.

Appendix B: Determination of innermost stable circular orbit
With the dimensionless quantities x = m −1 • r and l = m −1 • L, the conditions for the innermost stable circular orbit (23) can be rewritten as We can solve l 2 from Eq. (B.16) as and substitute it into Eq. (B.17) to obtain Since A(x) is positive for any timelike particles outside the event horizon, we can find E 2 from Eq. (B.20) as It clearly shows that y isco is merely determined by the charge q. As q = 0 and q = 1, y isco arrives at its biggest and smallest values, respectively, that y isco ∈ [4,6] . Since y isco only depends on q, Eq. (B.28) ensures that E isco is unaffected by the dimensionless bounce parameter α. E isco reaches its biggest and smallest values as q = 0 and q = 1, respectively, that . where Eq. (7) has been used to replace α with the rescaled bounce parameter λ. On D H , x isco reaches its biggest value at q = 0 and λ = 0, and its smallest one at q = 1 and λ = 1, respectively, that x isco ∈ √ 15, 6 . (B.33) When q = 0, E isco , l isco and x isco return to those of the black-bounce-Schwarzschild spacetime [37]. When both q and λ disappear, they reduce to those of the Schwarzschild spacetime with E isco,Sch = 2 √ 2/3, l isco,Sch = 2 √ 3 and x isco,Sch = 6.