Timelike bound orbits and pericenter precession around black hole with conformally coupled scalar hair

We investigate the geodesic motions of timelike particles around a static hairy black hole with conformally coupled scalar field. We mainly focus on the effects of the scalar charge and electric charge on the marginally bound orbits (MBO), innermost stable circular orbits (ISCO) and on the precessing orbits around this black hole. Our results show that both the scalar and electric charges suppress the energy as well as the angular momentum of the particles in the bound orbits. Then, we study the relativistic periastron precessions of the particles and constrain the charge parameters by employing the observational result of the S2 star's precession in SgrA*. It is found that the constraints on the charge parameters from S2 star's motion are tighter than those from black hole shadow. Finally, we analyze the periodic motions of the particles and figure out samples of periodic orbits' configurations around the hairy black hole.

Black hole provides natural platform to test gravity in the strong field regime.Recent breakthroughs in astrophysical observation including gravitational waves (GW) [1][2][3], black hole shadows [4][5][6][7][8] and the supermassive object in the center of our Milky Way [9] give strong evidences of the existence of black holes in our Universe.Those observations further demonstrate the success of Einstein's general relativity (GR), but the uncertainties in the data also leave some space for alternative theories of gravity.Therefore, plenty of research on the properties in the vicinity of black holes have been carried on since those remarkable progresses provide potential tools to test or distinguish GR and modified gravity theories.Among them, the structures of geodesics for massive particles or photons around the black holes have been attracting extensive interest, since the orbits of photons are the first step to study the black hole shadow related physics, while the bound timelike orbits are key ingredients to study the motions of stars around SgrA* as well as GW radiation during the black hole merger.Those studies could also inversely help us further understand the essential structure of gravity from the central black holes.From the timelike geodesic, one important type of bound orbit is the precessing orbit.It is what the Mercury undergoes, which gives one of the early evidences of GR [10].Later, the precessing of the stars orbiting around the supermassive black hole in the central of SgrA* has been disclosed in [11][12][13][14][15].The measurement of the precessing orbit of S2 star [16] has then been fitting in different theories of gravity to constrain the model parameters, see for examples [17][18][19][20].Thus, these precessing orbits may serve as potential tools to test the alternative theories of gravity.Another important type of bound orbit is the periodic orbit.It is known that the generic orbits can be treated as small deviations from periodic orbits which encode fundamental information about orbits around a central black hole [21].Especially, a zoom-whirl classification for the periodic orbits of the massive particles was proposed in [21], which is characterized by three topological integers (z, w, v) representing zoom, whirl and vertex behaviors of the orbit, respectively.It was previously found in [22] that the periodic orbits with zoom-whirl pattern may be helpful for faster computation of adiabatic extreme mass ratio inspiral process, which more recently was further testified in two models [23,24].It is believed that periodic orbits could give more significant information about strong-field properties of the spacetime than the precessing orbits, so the taxonomy of the periodic orbits for timelike particles has inspired considerable extensive studies  and references therein.
On the other hand, it is commonly accepted that a more general theory of gravity beyond GR is eagerly required, and so lots of modified gravitational theories have been proposed [46][47][48].One of the remarkable ways is to introduce an additional scalar field into the action of GR.The underlying motivations mainly stem from three aspects: (i) scalar field is a potential candidate for dark matter, dark energy and inflation which were proposed to explain some observations of our Universe, but their essences are still unclear.(ii) scalar field may be a kind of ubiquitous composite in nature, for example, the ultralight axions are indispensable in string theory [49].(iii) the hypothesis of no hair theorem in classical GR states that black holes have only three characterized quantities, i.e., the black hole mass, electric charge and angular momentum [50,51].However, the no hair theorem of black hole is still not be verified and from the current observation there is no hint for the absence of other fundamental quantity describing black holes.So the introducing of scalar field into the action could be a direct way to testify it.A minimally coupled scalar field with gravity usually does not obey the Gauss-law, so in this case a black hole cannot have a non-trivial regular scalar hair [52] such that the no hair theorem holds.However, the output is different when the non-minimal couplings between the gravity and scalar field are introduced.As the first counterexample to the no hair theorem using scalar field, the authors of [53,54] found that an Einstein-conformally coupled scalar theory could lead to a secondary scalar hair around the Bocharova-Bronnikov-Melnikov-Bekenstein (BBMB) black hole which was numerically reproduced in [55].But the scalar field in this sector diverges at the horizon which is difficult to be physically acceptable.This disadvantage was then erased by introducing a cosmological constant into the solution, named as the "Martinez-Troncoso-Zanelli (MTZ)" black hole, in which the scalar field singularity was pushed behind the event horizon [56,57].It is noted that the MTZ black hole only has a spherical or hyperbolic horizon depending on the sign of the cosmological constant, but the planar solution is not allowed, which inspired some extensive studies [58][59][60].In particular, it is natural to add the Maxwell field into Einstein-conformally coupled scalar theory and extend the theory into Einstein-Maxwell-conformally coupled scalar theory of which the action is [56] where κ = 8πG, R is the Ricci scalar, F µν is the strength of Maxwell field A µ and φ is the scalar field.The author of [61] then constructed a static spherically symmetric black hole with the conformally coupled scalar field which we will review in next section.Many physical phenomena on this black hole have been extensively studied, for instance, the stability of the black hole against perturbations [53,62], the Hawking radiation of charged particles [63], the black hole shadow [64], the gravitational lensing effects [65] and photon ring bound [66], in which the conformally coupled scalar field reflects significant influences.Thus, the aim of this paper is to study the bound timelike geodesics, precessing and periodic orbits around the static spherically symmetric black hole with the conformally coupled scalar field to further understand the gravitational aspects of the corresponding theory of gravity.Starting for the timelike geodesic equation in this black hole, we will first analyze the effects of the scalar charge and electric charge on the characterized quantities of the marginally bound orbits and innermost stable circular orbits.We will see that the energy and the angular momentum of the particles in the bound orbits are both suppressed by the larger charge parameters of the hairy black hole.Then, we shall calculate the relativistic periastron precessions of the particles, and connect the obtained precession angle with S2 star's precession in SgrA* to constrain the charge parameters, which we will find is tighter than those from black hole shadows addressed in [64].Finally, we will move on to the periodic motions of the particle and analyze periodic orbits' configurations around the black hole.Our current studies could be helpful for us to better understand the gravitational structure of black hole with scalar and electric hairs.
The remaining of this paper is organized as follows.In section II, starting from the timelike geodesic, we will analyze the properties of bound orbits around the black hole with conformally coupled scalar hair.In section III, we will focus on the precessing orbits to give preliminary constraints on charge parameters with the use of the orbits of S2 star in SgrA*.We figure out samples of periodic orbit and discuss the effects charge parameters in section IV.The final section contributes to our conclusions and discussions.In this paper, we shall use the units G = c = 1.Additionally, in our theoretical evaluations, all the physical quantities are rescaled by the black mass parameter M to be dimensionless, therefore, we set M = 1 unless we especially restate.

II. BOUND TIMELIKE ORBITS AROUND BLACK HOLE WITH CONFORMALLY COUPLED SCALAR HAIR
The static and spherically symmetric metric for the spacetime with conformally coupled scalar hair derived from the action (1) is given by [61] with And the profiles of the matter fields are In the above solution, M , q and s denote the mass parameter, electric charge parameter and conformally-coupled scalar charge parameter, respectively.According to these parameters' relation, the metric (2) and (3) can describe different spacetimes: (i) When s > M 2 − q 2 , f (r) = 0 has no real root, so the metric describes a naked singularity without horizon.
(ii) When M 2 − q 2 ≥ s ≥ 0, f (r) = 0 has two positive real roots r ± = M ± √ M 2 − q 2 − s.So the metric describes a black hole with Cauchy horizon r − and event horizon r + .In particular, when s = M 2 − q 2 the black hole becomes extreme with r − = r + , and the geometry is then a Reissner-Nordström (RN) black hole with an electric charge Q 2 = q 2 + s.
(iii) When −q 2 < s < 0, the profile of scalar field in (4) is imaginary, which is not considered in this article.
(iv) When s < −q 2 , the metric also describes a black hole with two horizons r ± , in which the coefficient of 1/r 2 -term is negative differentiating from the RN black hole.Therefore, the black hole in this case is also dubbed a mutated RN black hole.
Thus, in this paper, we will work in the spherical black hole with comformally coupled scalar field of which the charge parameters should satisfy M 2 − q 2 ≥ s ≥ 0 or s < −q 2 .It is noted that the black hole spacetime ( 2) is very close to the "BBMB black hole" constructed in [53,54], except that the scalar field here is regular instead of non-regular.The Ricci scalar, the square of the Ricci tensor and the Kretschmann scalar for the black hole are then calculated as which denote that the charge parameters, q and s, have influence on R 2 and K 2 , but the Ricci scalar always vanishes.Then, we shall start from the timelike geodesic around the spherically symmetric black hole (2).To this end, we consider one massive particle moving on the equatorial plane (θ = π/2) , and the corresponding Lagrangian is with the dot represents the derivatives with respect to the affine parameter λ.Then the Euler-Lagrangian equation gives us where E and L denote the conserved energy and orbital angular momentum per unit mass of the particle.Then considering that g µν ẋµ ẋν = −1 for timelike geodesic, we can derive the radial equation of motion as from which we can get the effective potential of the radial motion of the particle It is obvious that the effective potential V ef f depends explicitly on the angular momentum L, the radius r, the charge parameters q and s.In particular, V eff | r→∞ = 1 such that E = 1 is the upper limit for the bound orbits because in this case we have ṙ2 = E 2 − V eff | r→∞ = 0.The particle's motion is determined by the balance of the effective potential and the energy E. Therefore, according to [67], we can then qualitatively discuss the properties of various possible motions of the particles around the black hole by analyzing the effective potential's profile, including the scattering orbits (E > 1) for which particle can escape into infinity, and the bound orbits when E ≤ 1.Here we are interested in bound orbits for which the energy and angular momentum of the particle should be in certain region, as we will see later.Before analyzing the properties of general bound orbits, we will first study two particularly noteworthy bound orbits, i.e, the marginally bound orbits and innermost stable circular orbits [68].
• Marginally bound orbits (MBO) : MBO is characterized by a balance between the particle's kinetic energy and the gravitational potential energy, resulting in a zero radial velocity.Marginally bound orbits are of particular interest in black hole physics as they correspond to the innermost orbits a particle can have without falling into the black hole.It is a specific unstable circular orbit around the black hole, of which the maximum effective potential is 1, so it yields to Plugging Eq.( 12) into Eq.(13), and adopting the convention M = 1, we directly obtain two equations Then, by simultaneously solving the above equations, we can derive the radial distance r M BO and the angular momentum L M BO of the particle around the black hole with conformally coupled scalar hair with respect to the charge parameters.The results are depicted in FIG. 1.For positive scalar charge parameters, an increase in |s| leads to a gradual decrease in both the radius r M BO and angular momentum L M BO .Conversely, for negative scalar charge parameters, the trends in their variations are reversed.It is noteworthy that the blue, green, and red curves in the figure represent q = 0, q = 1/4, and q = 1/3, respectively, which show that for the black hole with larger q, the MBO always has smaller radius and angular momentum.It is noticed that the gap in the figure raises because the black hole exists only for M 2 − q 2 ≥ s ≥ 0 or s < −q 2 as we addressed in the previous section.FIG.1: The radius and angular momentum as functions of scalar charge parameter for the particles in the marginally bound orbit around the black hole with conformally coupled scalar hair.Here we choose samples of electric charge parameters.
• Innermost stable circular orbit (ISCO): ISCO is the closest stable circular orbits that a particle can orbit the black hole.So a small perturbation could not cause the particle to either fall inward or escape outward.ISCO encodes key information of the dynamics of matter and energy in extreme gravitational fields, so it is of great importance in black hole physics such that studying accretion disks around black holes, as well as for modeling and predicting gravitational wave signatures from merging black holes.ISCO should satisfy the following conditions, which are deduced as By solving the above equation group, we can extract the key information of the ISCO.The radius r ISCO , angular momentum L ISCO and energy E ISCO with respect to the scalar charge s for different electric charges q are plotted in FIG. 2, which shows that r ISCO , L ISCO and E ISCO all decrease as the increasing of both charge parameters in the allowed parameter space.Then, we shall analyze the properties of general bound orbits.Since the bound orbit only exists when the particle's energy is no more than 1 as we previously mentioned, so here we shall focus on the cases between MBO and ISCO to discuss the properties of bound orbits.The analysis mainly involves studying the radial motion of timelike particles and profiles of its effective potential.For the black hole with fixed parameters, the angular momenta for the general bound orbits should be between the corresponding L ISCO and L M BO , and the particle in a bound orbit with selected angular momentum can only have certain region of energy.To explicitly show these properties, we show the sample of effective potential for various bound orbits with q = 1/4 and s = −1/4 in FIG. 3, which has two extremal points moving closer to each other as the angular momentum decreases from L M BO = 4.09067 and finally merging together when it reaches L ISCO .The orbits with angular momentum between L ISCO and L M BO could be bound by the corresponding potential well, depending on the particle's energy.As shown in FIG. 4, we figure out the allowed (L − E) space for the bound orbits.It shows that the two charge parameters will both alter the allowed (L − E) region of bound orbits.Then, we shall analyze the radial motions in the bound orbits for the particles.The sample results with fixed L = (L BM O + L ISCO )/2 for different energies are shown in FIG. 5.In this case, only the particles with the energy satisfying 0.956375 ≤ E ≤ 0.969664 can have the bound orbits, where ṙ2 = 0 has at least two roots.In details, for the purple curve (the first solid line from top) with E = 0.969664, ṙ2 = 0 has only two roots and the bound orbits exist in the regime between the two roots.Then for the particle with smaller energy, ṙ2 = 0 has three roots and the bound orbits exist in the regime with ṙ2 > 0 between the last two roots, and the bound orbits converge to the ISCO with E = 0.956375.Additionally, when the particle's energy is beyond the aforementioned regime, ṙ2 = 0 has no or single root so no bound orbits exist.
Next, we shall study the special subclasses of bound orbits, the precessing and periodic orbits, around the hairy black hole with conformally coupled scalar field.

III. ORBITAL PERICENTER PRECESSION AND CONSTRAINTS ON CHARGE PARAMETERS FROM S2 STAR'S MOTION
Due to the spherically symmetric nature of the background black hole, the motion of particles is solely determined by the radial coordinate r and the azimuthal angle ϕ.For a precessing orbit with two turning points r 1 and r 2 , the timelike particle will be reflected between the turning points.So, the apsidal angle ∆ϕ passed by the particle in this process is Especially, for the periastron precession, the particle's trajectory is described by [67] r of which the turning points are the periastron, r p = a(1 − e), and apastron, r a = a(1 + e), respectively.Here a and e are the semi-major axis and the eccentricity of the orbit, while Ψ denotes the intersection angle between the semi-major axis and radial of the orbit.Subsequently, the above angle ∆ϕ in this case is evaluated as with where the expression of dϕ/dr is derived from ( 10) and ( 11), and dΨ/dr is obtained from (21).Recalling that the turning points satisfy ṙ| r=rp = ṙ| r=ra = 0, we can obtain Then, to perform the integral in ( 22), we shall do some approximation on the integral function (23).We consider the weak gravitational field case, i.e., M /r is treated as small quantity.Then we rescale the charge parameters as in which the first line represents the result for Schwarzschild black hole in GR.Then, substituting (25) into ( 22), we will obtain the precessing angle as which will recover the result for Schwarzschild black hole, △ω GR = 6πM a−ae 2 , if the charge parameters vanish.Now, we are ready to use the observations of Schwarzschild precession of the S2 star around Sgr A* to test the current theory.The Schwarzschild precession of the S2 star around Sgr A* was measured by the GRAVITY project with the technics of spectroscopic and astrometric measurements, which is the first measurement of periastron precession for a star around supermassive black hole [16].The best-fit ratio between the measured result to that predicted by GR is given by within 1σ uncertainty which leaves some space for alternative gravity beyond GR.So, we shall presuppose that the central supermassive black hole in SgrA* is the hairy black hole with conformally coupled scalar field, and the precessing motion of s2 star follows the trajectory we studied above.Subsequently, the shift f sp is then evaluated by ) . ( Then we can calculate f sp by considering the realistic values of various physical quantities in the S2 star's motion, such as the mass of SgrA*, M SgrA = 4.3×10 6 M ⊙ , the luminosity distance between the S2 star and SgrA*, D = 8.35 kpc, and the eccentricity and the semi-major axis of the orbit e = 0.884649 and a = 125.058mas [16].The results in (S, Q) parameter space are depicted in FIG.6, the left of which we focus on hairy black hole with M 2 − q 2 ≥ s ≥ 0 while the right panel corresponds to s < −q 2 .It is obvious that the S2 star's motion can further constrain the charge parameters.In details, for s ≥ 0 in the left plot, only the charge parameters in the lower left region of the green solid curve, which corresponds to lower bound f sp = 1.1 − 0.19, match the S2 star's observation, while for s < −q 2 in the right plot, the charge parameters are constrained into the upper right region of the blue solid curve which gives the upper bound f sp = 1.1 + 0.19.By carefully comparing the constrained parameter region from black hole shadow addressed in [64], we find that our constraints on the charge parameters from S2 star's motion are tighter.

IV. EFFECTS OF THE CHARGE PARAMETERS ON PERIODIC MOTIONS
In this section, we will move on to study the periodic orbits, which repeat themselves over time.They describe the regular and predictable motion of objects in space, ranging from planets orbiting the Sun to stars moving within galaxies.The study of periodic orbits is crucial for understanding the stability and dynamics of celestial systems, as well as for making predictions about their long-term behavior.As suggested in [21], the bound periodic orbit can be described by a unique rational number b as where the angle ∆ϕ is defined in (20), and the integers z, w and v are the zoom number, whirl number and vertex number in the orbit, respectively.Therefore, the values of b can be further calculated via We show the results of b for periodic orbits as functions of energy E as well as angular momentum L for selected charge parameters in FIG.7 In FIG. 7, we choose ϵ = 0.3 and ϵ = 0.5 to describe b as a function of energy.For all curves, b first increases slowly with the increase of energy E and then experiences a sudden surge when E reaches a maximal value, but larger ϵ suppress the value of b.Moreover, we find that increasing both charge parameters will shift the curve into smaller energy.In FIG. 8, we fix E = 0.955 and E = 0.962, respectively, to describe b as function of angular momentum.For each curve with chosen charge parameters, as the angular momentum decreases, b increases smoothly and then explodes as L approaches its minimal values.Similarly, for black hole with both larger charge parameters, the curve is shifted into left corresponding smaller angular momentum of the orbits.Then in order to figure out the trajectories of the periodic orbit (z, w, v) around the black hole with conformally coupled scalar field, we shall first determine the corresponding energy E (z,w,v) and angular momentum L (z,w,v) .To this end, we firstly list the corresponding energy E (z,w,v) with fixed L = L M BO +L ISCO 2 (i.e., ϵ = 0.5) in Table I, and the corresponding angular momentum L (z,w,v) with fixed E = 0.96 in Table II, respectively.The effects of charge parameters indeed follow the findings in the previous figures of b.In details, as the electric charge q increases, both E (z,w,v) and L (z,w,v) gradually decrease.While when s > 0, both E (z,w,v) and L (z,w,v) gradually increase with the increasing amount of scalar hair, but when s < 0, the situation is exactly the opposite.This suggests that, under otherwise identical conditions, a massive particle orbiting around a hairy black hole in current model will have higher energy E and angular momentum L due to the presence of the electric charge q compared to that in Schwarzschild case.While the effect of the scalar hair s on the energy E and angular momentum L of the massive particle compared to the Schwarzschild black hole depends on its sign.With the data in hands, we can figure out the periodic orbits with (z, w, v) in the polar coordinates (r, ϕ).By selecting several charge parameters, we show the exemplified results with fixed E = 0.96 in FIG. 9 and FIG. 10 of which the corresponding L is present in Table II.It is obvious that the zoom number z describes the number of the leaf pattern of the orbit, the whirl number w describes the number of nearly circular whirls close to periastron per leaf, and the vertex number v gives the order of which the leaves are traced out.It is found by comparing the orbits with different (z, w, v) in each row that for orbits with larger z, the leaf pattern grows such that the orbit becomes more complicate.Moreover, by comparing the orbits with the same (z, w, v) in each column, we see that both larger charge parameters will correspond to larger outermost trajectory in the orbit.We consider two values of energy, E = 0.955 (left column) and E = 0.962 (right column).The energy E for the samples of periodic orbits (z,w,v) with selected charge parameters.Here we fix the particle's energy as L = L M BO +L ISCO 2 or ϵ = 0.5.

V. CLOSING REMARKS
Searching for alternative theories beyond GR still attracts the attention of physicists, especially by introducing additional scalar fields into the action of GR.An intuitive approach for investigating partial properties of such gravity theories is to study the appearance of the black hole images in specified spacetime backgrounds.In this article, we  The angular momentum L for the samples of periodic orbits (z,w,v) with selected charge parameters.Here We fix the energy of the particle as E = 0.96.FIG.9: The trajectories of the periodic orbits (z,w,v) with selected scalar charges by the electric charge q = 1/4.We fix E = 0.96 for all orbits and the corresponding L were listed in Table .II.
( We fix E = 0.96 for all orbits and the corresponding L were listed in Table .II.
have considered a static spherically symmetric black hole with the conformally coupled scalar field, and study the properties of the bound timelike geodesics, precessing and the periodic orbits.Our results suggest that the scalar and electric charges have significant influences on the motions of the timelike particles.We also provide tighter constraints on the charge parameters from the observational result of the S2 star's precession in SgrA* than those from the black hole shadow.
We firstly investigated the general properties of bound orbits of the timelike particles, which exist between the MBOs and the ISCOs, in the considered hairy black hole.We demonstrate that under the specified parameter condition, the increasing of scalar and charge parameters will both gradually decrease the radius of MBOs and ISCOs and angular momentum of the particles on these orbits.The allowed (L − E) regions for the bound orbits are then altered by the increasing parameters.
Then we studied the orbital pericenter precession and constraint on charge parameters by employing the observational results from S2 star's motion.The error margin of the ratio between the measured precessing angle to that predicted by GR leaves some space for alternative gravity beyond GR.By considering the realistic values of S2 star's motion, we have found that in the weak field case the constraints on the charge parameters are tighter than the constraints from the black hole shadow.
Finally, we investigated the effects of the scalar and charge parameters on the periodic orbits.The results indicate that under the identical conditions a massive particle orbiting around the considered hairy black hole will have smaller energy and angular momentum with a larger charge parameter.These conclusions are also consistent with the previous sections.We also illustrated the trajectories of the periodic orbits and the larger parameters correspond to larger outermost trajectory in the orbits.In general, during the periodic motions, the particle's energy and angular momentum will decrease, which is accompanied by the emission of gravitational waves.Thus, the periodic orbits play important role in the study of gravitational wave emitted in the process of the two initial black holes of extreme mass ratio approaching to each other (EMRI), which was firstly proposed in [69] and further emphasized in [21].More recently, by considering the central black hole as a polymer black hole in loop quantum gravity [23] and Einsteinian cubic black hole [24], the connections between the periodic orbits of timelike particles and the gravitational wave radiated from EMRI were further preliminarily studied.So along this line, based on the current work, it is interesting to study the effects of charge parameters on the gravitational radiation in the EMRI system via the periodic orbits around the hairy black hole.Moreover, it is also of interest to explore the possible effect of this gravitational radiation on the evolution of periodic orbits [70].
Using black hole images to study the characteristics of different gravitational theories has always been one of the most concise and intuitive methods to distinguish between GR and alternative theories.With the development of observation technology, more and more details of black hole images will bring us increasingly strict parameter constraints.This will provide evidence support for further excluding certain theories of gravity.Although in this work we have investigated the timelike bound orbits and pericenter precession around black hole with conformally coupled scalar hair, there are still other directions that could provide more insights about the current gravitational theory, for example, black hole illuminated by thin accretion disks or immersed in dark matter halo, which we will leave to future work.
timelike orbits around black hole with conformally coupled scalar hair III.Orbital pericenter precession and constraints on charge parameters from S2 star's motion IV.Effects of the charge parameters on periodic motions V. Closing remarks I. INTRODUCTION

3 FIG. 2 :FIG. 3 :
FIG.2:The radius, angular momentum and energy as functions of scalar charge parameter for the particles in the ISCO around the black hole with conformally coupled scalar hair.

3 FIG. 6 :
FIG. 6: The density plot of fsp of the S2 star's motion in the parameter (q/M, s/M 2) space by treating the current hairy black hole as the central supermassive black hole in SgrA*.The left plot shows the case with S > 0, while the right one illustrates the case with S < 0.

3 FIG. 7 :
FIG. 7:The number b of periodic bound orbit as a function of energy E for selected black hole parameters.We choose two angular momenta by setting ϵ = 0.3 (left column) and ϵ = 0.5 (right column).

3 FIG. 8 :
FIG.8:The number b of periodic bound orbit as a function of the angular momentum L for selected black hole parameters.We consider two values of energy, E = 0.955 (left column) and E = 0.962 (right column).

3 FIG. 10 :
FIG.10:The trajectories of the periodic orbits (z,w,v) with selected electric charges by the scalar charge s = −1/4.We fix E = 0.96 for all orbits and the corresponding L were listed in Table.II.