Bounds on spinning particles in their innermost stable circular orbits around rotating braneworld black hole

We study the innermost stable circular orbit (ISCO) of a spinning test particle moving in the vicinity of an axially symmetric rotating braneworld black hole (BH). We start with the description of the event horizon, static limit surface and ergosphere region of such BH and bring out the effect of tidal charge parameter on ergosphere. It is found that the ISCO of rotating braneworld BH is very sensitive to braneworld BH parameter C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ C $$\end{document} (also known as tidal charge parameter) in addition to its rotation parameter. We further discovered that the orbital radius of the spinning test particles changes non-monotonously with the braneworld BH tidal charge parameter. It is found that for rotating brane-world BH the allowed range of the particle spin grows as the tidal charge parameter C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ C $$\end{document} decreases, in contrast with the Kerr–Newman BH. We also found the similar behavior of the particle’s spin for the braneworld Reissner–Nordstro¨\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\ddot{o}$$\end{document}m (C<0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ C <0$$\end{document}) BH in contrast with its counterpart having (C>0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ C >0$$\end{document}).


Introduction
In this era of advanced laser interferometers (i.e., Laser Interferometer Gravitational Wave Observatory (LIGO) [1], Virgo [2], and KAGRA [3]) and high resolution telescopes (i.e. EventHorizon Telescope (EHT)) [4][5][6][7][8][9], the dynamics of bodies (point-like or extended) around a central black hole (BH) is of indispensable interest because a compact body such as a BH or neutron star of few solar mass orbiting around a central massive BH (of mass ∼ 10 4 to 10 6 M ) is usually referred to an extreme mass ratio inspiral and is a propitious source of gravitational waves for the space-based interferometer eLISA [10] and DECIGO [11] in the near future. In a e-mail: unucamendi@gmail.com b e-mail: ricardo.becerril@umich.mx c e-mail: pankaj.sheoran@umich.mx these scenarios of extreme mass ratios, the numerical relativity approach is not very efficient [12] and the approximation methods like the effective one body formalism [13,14] become handy.
Previous to these theoretical advancements, the dynamics of a non-spinning particle moving on a geodesic around a Schwarzschild BH was first studied by Kaplan in [15]. Since then the dynamics of non-spinning particles has been studied vastly by researchers [16][17][18][19][20][21][22][23][24][25][26][27][28][29][30][31][32][33][34]. The study of non-spinning particles in the vicinity of different BHs showed that the motion of the massive or massless test particles gets affected by the BHs parameters like mass, charge and rotation. In addition to these parameters, it is also found that if the central object (i.e. BH) is inspired by alternative theories of gravity then the motion of the non-spinning particles in the vicinity of BH gets influenced by the extra parameter known as the deviation parameter [29,[31][32][33][34]. On the hand, the dynamics of spinning particles around the BHs has not received as much attention as it needs to be.
The study of spinning particles started with the pioneering works of Mathisson [35,36], Papapetrou [37,38] and Dixon [39] (MPD), where they developed the dynamical equations of the spinning particles moving in curved backgrounds by considering the "pole-dipole" approximation only. The theory of the spinning particles was further developed by several authors [40][41][42][43][44][45][46], where it was showed that the four-velocity vector and the corresponding conjugate momentum vector are not parallel to each other which is in contrast with the case of the nonspinning particle where these two corresponding quantities are parallel. In [41,46] it was also showed that, for the case of the spinning particle, the four velocity might change from timelike to spacelike if the spin of the particle is greater than a certain critical limit whereas the conjugate momentum always remains timelike along its trajectory and satisfying the conservation of mass relation (P μ P μ = −m 2 ). It is relevant to point out here that the study of the dynamics of the spinning particle in a curved spacetime by considering the "pole-dipole" approximation is valid only for homogeneous fields and fails when nonhomogeneous fields are taken into consideration. However, the study of spinning particles in curved spacetimes is extended to higher order-corrections in [47,48]. It was showed in [49][50][51][52] that the four velocity of a spinning particle will always remain timelike and never transforms to spacelike (i.e. it avoids the Superluminal problem) if one considers the coupling between spin and gravity via the gravitomagnetic moment, that is to say, that the superluminal regime can be avoided if multipole effects are taken into account. Taking into consideration the effect of multipole moments becomes important when one wishes to describe the self-gravitating compact bodies and gravitational radiation emitted by them. Recently, in [53][54][55] the authors have studied the gravitational radiation emitted by the spinning particles moving in curved backgrounds as well as the associated chaos effect.
From the literature [41][42][43][44][45][46][49][50][51][52][53][54][55][56][57][58][59][60][61][62][63][64][65][66][67][68], it is well understood that a nonspinning particle will follow a geodesic in the vicinity of a BH. On the contrary, if the interactions (i.e. gravitational self-force) of a test particle are taken into account then its trajectory is no longer a geodesic [69][70][71]. Also, in the case of a spinning test particle, the motion is non geodesic due to an additional force known as the spin-curvature force [40,72] which comes into play due to the interaction of spin of the particle and the curvature of spacetime around a massive central object. In this work we only consider the spin-curvature coupling and discard all other reaction of the particle with the background of a BH and study the ISCO of a spinning test particle in its vicinity. The ISCO of a spinning particle around a Schwarzchild and Kerr BH was first studied by Suzuki and Maeda in their pioneering work [56]. Later on, the ISCO parameters of the spinning particles were obtained for a slowly rotating Kerr BH up to quadratic order, in terms of rotation parameter a and particle's spin by Jefremov et al [64].
The study of the inner most stable circular orbit (ISCO) is very important from the point of view of the gravitational wave astronomy because the circular orbits which are located inside the ISCO are unstable under a perturbation, and then they can be taken as an initial point for the final merger of any binary system [73][74][75][76][77][78][79]. Also, the study of ISCO for a given BH tells about the properties of its background geometry because the motion of a particle depends on its mass, charge, rotation and extra deviation parameters (coming from alternative theories of gravity).
Since the study of spinning particle dynamics has been done mostly in the context of Einstein's theory but scarcely for alternative theories of gravity, in this work we study the ISCO of a spinning massive test particle with arbitrary spin s in the vicinity of a rotating braneworld BH [80][81][82][83][84] which has an additional parameter C known as the braneworld tidal charge parameter in addition to the usual mass M and rotation a parameters.
The braneworld models are an effective four-dimensional version of higher-dimensional string theory [85,86]. According to these models, our physical universe nests on a 3brane of a higher dimensional spacetime while gravity enters as an extra spatial dimension [87]. Therefore, studying the behavior of gravity in braneworld models shades light on the physical signature of higher dimensions on our fourdimensional physical world. In this context, the braneworld BHs are interesting to study. Additionally, the study of braneworld BHs is fascinating in many other ways, one that the non-rotating braneworld BHs can be characterized with a Reissner-Nordström type geometry [80], and the axially symmetric rotating braneworld BHs as a Kerr-Newman type geometry [81]. Second, the tidal charge parameter can have both positive and negative values [80,81] unlike in Einstein-Maxwell theory where the square of the electric charge is always positive. These exciting properties of braneworld BHs gained a lot of attention from researchers in the past decade and were studied in works related to accretion phenomenon [83], solar system tests [88], quasiperiodic oscillations (QPOs) [89], shadows of BH [90], gravitational lensing [91] and many more topics [82,84,92,93]. More recently, in [94] the shadow of M87* was used to constrain the curvature radius of rotating braneworld BH and in particular about the tidal charge parameter C.
In this context and for simplicity, we consider only the "pole-dipole" estimation and numerically investigate the physical behavior for the ISCO parameters (radius r , energy E and orbital angular momentum component in the z-direction L z ) using the superluminal constraint condition and imposing the Tulczyjew spin-supplementary condition (TSSC) for the spinning massive test particle in the rotating braneworld BH background. To be consistent, we showed that our results for the ISCO parameters of non-extremal cases of Kerr and Kerr-Newman like BHs (C > 0) match exactly with the results obtained in [64,67] for Kerr and Kerr-Newman BHs.
We organized our paper as follows, in Sect. 2 we start with a brief review of the MPD equations (equations of motion (EOM) of a spinning particle) in a curved background. Next, we divide the Sect. 3 into two subsections: In the first subsection, we study the behavior of the event horizon, static limit surface and ergosphere ( Fig. 1) for different values of the tidal charge parameter C and numerically present the bound on the rotation a and tidal charge parameters C (Fig. 2) of rotating braneworld BH and in the second subsection, we find its conserved quantities and present the EOM for the spinning particles moving around it. In Sect. 4, we find the expression for the effective potential and show its behavior in Fig. 3 for different values of the particle spin S and the parameter C. In Sect. 5, we study numerically the behavior of ISCO parameters for spinning particles moving in the vicinity of rotating braneworld BH. It comes out that for rotating braneworld BH the range of the particle spin, for which the behavior of ISCO parameters are physical, increases as the parameter C decreases. In this section we also present a summary of the results that we have obtained from Figs. 4, 5, 6 and 7. We summarize and give conclusions about our work in Sect. 6. Finally, in the A we showed the explicit form of the equations that are used to study the behavior of ISCO parameters.
Throughout the paper, we have worked with (−, +, +, +) signature and fixed the fundamental constants (c and G) to unity. Additionally, the transformation we have used for projecting any four-vector in the tetrad frame is as follows: ν x ν where the indices in curve brackets indicate tetrad components while greek indices mean spacetime components.

The general equations of a spinning test particle
In this section, we briefly review the general EOM of an extended body for a spinning test particle in a curved spacetime in the pole-dipole approximation. These EOMs were first obtained by Mathisson [35,36] and Papapetrou [37,38] and later fine tuned by Dixon [39]. The final form of the set of EOMs reads where, τ , U μ ≡ dx μ /dτ , u μν , P μ , S μν and R μ νρσ represent the affine parameter of the orbit, four-velocity, angular velocity tensor, conjugate momentum vector, spin angular momentum tensor and Riemann tensor of a curved spacetime, respectively. The EOMs mentioned above were obtained by considering the curved spacetime lagrangian (L = L (c 1 , c 2 , c 3 , c 4 ) with the "pole-dipole" approximation. The quantities c 1 , c 2 , c 3 and c 4 are four independent invariants [41,43] defined as: Particles in curved spacetime do not follow geodesic motion due to the coupling between the spin tensor and the curvature of the background geometry and hence present non zero acceleration as shown in Eq. (1). The case of spinning particles is considered as an open system because in this case the P μ and U μ are not proportional to each other which gives rise to 14 unknown variables (i.e. four for P μ and U μ each and six for the spin tensor S μν ), while we have only ten equations at hand. Hence, to close this open system we need extra constraint conditions (also known as the spin supplementary condition (SSC)). In the literature available, there are three widely studied SSCs that can be used to close the open system of equations of motion, namely: (i) the Papapetrou and Corinaldesi SSC S 0i = 0 leading to the no dipolar mass moment condition, (ii) The Mathisson-Pirani SSC S μν U ν = 0 implying that orbits are helical in nature, (iii) the Tulczyjew-Dixon SSC S μν P ν = 0 which leads to an exact solution of Eqs.
(1) and (2) which in turn gives conservation of the mass M and the spin S of the spinning particle as shown in [46] and defined by the relation: In this work, we use the Tulczyjew-Dixon SSC. For convenience we also normalize the affine parameter τ with the help of a normalized momentum (V μ ≡ P μ /m) as Now, using the Eqs. (5) and (6) together with the Tulczyjew-Dixon SSC, one can obtain the following relation between It is easy to see from the above equation that the four-velocity and four-momentum are not parallel anymore due to existence of the particle spin. Thus, to study the dynamics of the spinning particle in the background of BH inspired by the modified theories of gravity we need to solve the Eqs. (1), (2), (7) and the Tulczyjew-Dixon SSC.

The Metric
The static BH solution for the braneworld scenario was first found by Dadhich et al. in [80]. The metric is obtained by solving the 5D Einstein field equations constrained to the 3Dbrane. We start this section by introducing the metric for a rotating BH in the Randall-Sundrum braneworld scenario and localized to a 3D-brane which was found in [81], the metric can be written in terms of Boyer-Lindquist coordinates as where Here a and C are the rotation and tidal charge parameters respectively. The parameter C gives rise to three class of BH metrics: (i) C > 0 corresponding to the Kerr-Newman like metric of Einstein-Maxwell theory, (ii) C = 0 corresponding to the standard Kerr metric of GR and (iii) C < 0 corresponding to the Kerr-Newman like braneworld metric with negative tidal effects.
In complete analogy to the Kerr-Newman BH, the metric (8) possesses two major surfaces, namely the event horizon (EH) and the stationary limit surface (SLS). The EH is a null surface characterized by the radius r E H or r + determined by equating the contravariant component g rr to zero (i.e. g rr = 0), The event horizon does exist if M 2 ≥ a 2 + C, where the equality leads to the extremal BH case (i.e., r + = r − ). Clearly, for the rotating braneworld BH the parameter C can be negative [81].
The SLS (characterized by the radius r SL S ) is obtained by setting the prefactor of the term dt 2 in the metric equal to zero (g tt = 0). The largest root of this equation gives the location of SLS around the BH It is clear from Eq. (11) that the SLS lies outside the EH except for the values θ = 0 and π where it touches it. The region between the SLS and EH is known as the ergosphere. The negative value of tidal charge parameter C for the case of braneworld-Kerr BH leads to the possibility of a greater EH and SLS in comparison with the corresponding Kerr-Newman like BH where charge parameter Q 2 = C > 0. It can be seen from the Eq. (10) that when the parameter C is negative, the EH do exists for extremal Kerr-Newman like braneworld BH if and only if the rotation parameter a 2 = M 2 − C > M. Hence, there exists a super-spinning case for the rotating braneworld BH and makes it more interesting in comparison with its Einstein-Maxwell theory counterpart (i.e., KN BH) [81] where this situation is impossible. Using Eqs. (10) and (11), we found that the ergoregion lies within the limit M < r < M +sin θ √ M 2 − C for the extremal case and from this, it is easy to conclude that the negative tidal charge parameter C (i.e., C < 0) of the rotating braneworld BH makes it more energetic than the usual Kerr-Newman BH which enables us to extract more rotational energy via the Penrose Process. Figure 1 shows the effect of rotation and tidal charge parameters on the ergosphere. In each one of the first three rows of the figure, we fix the tidal charge parameter and vary the rotation parameter. The ergoregion (shown in green) in the first three rows increases with the increment in the rotation parameter. However, in the last row we fixed the rotation parameter a > M and vary the tidal charge parameter. Here also, the ergoregion (shown in blue) increases as we increase the value of tidal charge for a fixed rotation parameter. Figure 2 illustrates the bounds on the rotation and tidal charge parameters for the rotating braneworld BH. The gray region gives the values of parameters a and C for which we get an event horizon whereas the yellow region corresponds to the values of the same parameters for which there exists a naked singularity. The boundary of these two regions corresponds to the extremal rotating braneworld BH case.

Conserved quantities
Conserved quantities play a very important role in the study of a test body around a rotating BH as these help in simplifying the equations of motion. We have a conserved quantity K ψ , if there exists a Killing vector ψ which satisfies the Killing equation ψ μ;ν + ψ ν;μ = 0 (the semicolon means covariant derivative), the conserved quantity along the trajectory of a spinning particle is obtained with the help of master equation: Since we are dealing with stationary and axially symmetric spacetimes, there are two Killing vectors ξ μ and Φ μ , one time-like and the other space-like respectively (thus in addition to the mass m and spin S of a spinning particle, there exists two more conserved quantities). These Killing vectors in a covariant tetrad base are given by where the quantities within the bracket denote their tetrad components.
For convenience hereafter, we restrict ourselves to the equatorial plane (θ = π/2) and use spin vector S μ instead of spin tensor S μν , satisfying the following relation in tetrad frame where, ε (b)(c)(d)(e) is the completely antisymmetric tensor which satisfies the relation ε (t)(r )(θ)(φ) = 1. As the spinning particle is restricted to the equatorial plane, the spin direction is always perpendicular to the plane θ = π/2. Hence, S (θ) is the only nonvanishing component of the spin vector S (b) . For simplicity, we set S (θ) ≡ −S. Here, the bold face S represents both the magnitude and the direction of a spin for the spinning particle. S > 0 represents that the direction of spin of a spinning particle is parallel to the rotation axis of the rotating braneworld BH, whereas S < 0 represents the antiparallel spin with respect to this axis. Furthermore, by using Eq. (16), the nonzero tetrad components of the spin tensor S (b)(c) in terms of spin vector S (b) are The conserved energy and the z component of the total angular momentum J z 1 of the spinning particle defined by using Eqs. (12), (13), (14) and (17) read Using the Eqs. (7) and (18), we obtained the following nonvanishing components of four velocity in terms of fourmomentum where Next, with the help of tetrad relations (15), the components of velocity fields for a spinning test particle in the background of a rotating braneworld BH can be written as Now, by using the Eqs. (18), (19) and (20), the explicit form of Eqs. (21)-(23) reads where σ s = r 2

Effective potential
Analyzing the effective potential is important for studying the dynamics of the spinning particle moving in the vicinity of a BH. Hence, in this section we wish to bring out the effect of rotation and tidal charge parameters on the effective potential of the rotating braneworld BH. We consider equatorial motion solely, taking into account only the equation for the radial velocity. We rewrite the Eq. (26) in quadratic form for the parameter E as where, Equation (28) can also be written in the form where Hereafter, we work only with V e f f (+) because this is the physical effective potential corresponding to the spinning particles with future pointing four momentum [95]. For the above stated reason in Fig. 3, we only plot V e f f (+) for some values of parameters (a, L z , C and S). The effective potential V e f f (+) of the spinning particle can have one or two extreme points as shown in Fig. 3, where it has circular orbits. The extreme point corresponding to maximum value of V e f f (+) leads to unstable circular orbits while the extreme point corresponding to its minimum value leads to stable circular orbits. From Fig. 3, it is evident that V e f f (+) is very sensitive to the test particle's spin S because the maximum and the minimum values of V e f f (+) increases with the rise of S. In Fig. 3, we fixed the tidal charge C parameter and varied the rotation parameter a along each row except the last row where parameter a is fixed and parameter C varies. It is clear from each of the first three rows of Fig. 3 that the unstable circular orbits shift closer to the EH of rotating braneworld BH as the rotation parameter a goes higher for the respective value of S. Additionally, we notice a similar kind of behavior in the last row, when the a is fixed and the parameter C grows from left to right along the row.

Inner most stable circular orbit (ISCO) and the superluminal constraint
In this section, we numerically study the behavior of ISCO parameters (i.e., r, E and L z ) for different values of rotation a and tidal charge C parameters as a function of the test particle spin S. It is known from literature [56,64,66,67] that to investigate the ISCO of the spinning test particle in a curved background we need to solve a system which is comprising of the following three conditions: Here, the Eq. (32) signifies that the radial velocity of spinning particle is zero and Eq. (33) means that the radial velocity for the spinning particle is constant, which means, in other words, that the radial acceleration should vanish. If both Eqs. (32) and (33) hold true simultaneously, the spinning particle will move in a circular orbit around the rotating braneworld BH. Now, in order to find the location of the ISCO, the point   Fig. 4 Behavior of ISCO parameters (r, E and L z ) and superluminal constraint (U 2 ) as a function of spin (S) for the spinning test particle in the rotating braneworld BH for corotating case (J z > 0). Here also, we keep the mass parameter M equals to unity where the maximum and minimum of the V e f f (+) meets, we need to introduce one more condition (i.e., Eq. (34)) in addition to the conditions described by Eqs. (32) and (33). The explicit form of these three equations are shown in Appendix is used to find the behavior of the ISCO parameters as a function of the particle spin S. It is worth mentioning here that for the case of spinning particle, the four-velocity and the four-momentum are not parallel and hence the four velocity can be timelike (U 2 < 0) or spacelike (U 2 > 0), where U 2 ≡ U μ U μ . Therefore, to study the behavior of the ISCO parameters as the function of the parameter S we need to take into consideration the superluminal constraint as well, which gives information of the region where the circular motion of the spinning particle will be superluminal (unphysical behavior) or subluminal (physical behavior). In Figs. 4 and 5, we numerically explore the behavior of ISCO parameters (r, L z and E) together with the superluminal constraint as a function of the particle spin S for both corotating (J z > 0) and counter-rotating (J z < 0) cases in the rotating braneworld BH. We divide Figs. 4 and 5 into two parts to bring out the effect of tidal charge parameter C on the ISCO parameters. In both Figs. 4 and 5, the parts (a) and (b) correspond to the parameter C = −2.0 and −8.0, respectively for fixed value of the parameter a = 1.2. With the help of Figs. 4, 5, 6 and 7, a brief summary of the results about the behavior of ISCO parameters(r, L z and E), superluminal constraint (35) and the particle spin S is presented as follows: • For the Kerr-Newman like braneworld BH case (i.e. a > 0 and C < 0), the physical values of ISCO parameters decrease as the spin S of the particle increases for the corotating cases, whereas for the counter-rotating cases, the ISCO parameter r increases and the orbital angular momentum parameter L z decreases with the increases in the spin S of the particle. It is found from Fig. 4, that the ISCO radius r and the corresponding orbital angular momentum L z increase as the parameter C decreases for corotating cases, wehereas for the counter-rotating cases as presented in Fig. 5, the ISCO parameter r increases, while the ISCO parameter L z decreases with decrease in parameter C. • Corresponding to the same value of parameter C, the ISCO radius for counter-rotating case is more than the corotating case (compare the corresponding part of both the Figs. 4 and 5 for reference). • The allowed range of the particle spin S, where the ISCOs exist for the Reissner-Nordström like braneworld BH (a = 0 and C < 0) increases monotonously with decrease in parameter C, whereas for the Reissner-Nordström like BH (a = 0 and C > 0), the allowed range of the particle spin S, where the ISCOs exist first decreases and then increases as the parameter C decreases. In Fig. 6, the light yellow and light green regions represented the Reissner-Nordström like braneworld BH and the dark yellow and dark green regions represented the Reissner-Nordström like BH. • Figure 7, showed the allowed range of the particle spin S for the rotating braneworld BH. The light blue and the light grey regions of rotating braneworld BH correspond to Kerr-Newman like braneworld BH (a > 0 and C < 0) whereas dark blue and dark grey regions correspond to Kerr-Newman like BH (a > 0 and C > 0). Similar to braneworld Reissner-Nordström BH case, the range of the particle spin for rotating braneworld BH increases with decrease in parameter C. However, for the fixed value of parameter C the allowed range of the parameter S decreases as parameter a of BH increases.

Summary and conclusion
In this paper, we have studied in detail the properties of the ergosphere for a rotating braneworld BH and highlight the effect of tidal charge parameter C on it. In Fig. 1, we showed that its ergoregion (blue) increases with the rise of the tidal charge parameter C for a fixed value of rotation parameter a. We also showed that when C < 0, it is possible to have a rotation parameter greater than unity unlike Kerr and Kerr-Newman BH where the BH rotation parameter a is always less than unity. Hence, for the rotating braneworld BH it is possible to have a super spinning case (ie., a > 1) as it is also mentioned in [82][83][84]. We also gave the numerical bounds on the rotation parameter a and tidal charge parameter C of a rotating braneworld BH as shown in Fig. 2. It is clear from the figure that the upper bound on the parameter a can be greater than the corresponding Kerr and Kerr-Newman BHs parameter. In Fig. 3, we showed the behavior of effective potential V e f f (+) as a function of radial parameter r of the orbit. Fig. 6 Phase plot between the parameters S and C, shows the bound on the particle's spin S for different set of values of parameter C in the case a = 0. The permitted region is between the curves formed by the dots.
Here, light and dark yellow regions represent corotating (J z > 0) cases of the spinning particle around the Reissner-Nordström like braneworld BH (C < 0) and the Reissner-Nordström like BH (C > 0), respectively. The same applies for the green region representing the counter-rotating (J z < 0) cases. The red-dot in the yellow and green regions shows the limiting value of parameter S for corotating and counter-rotating (J z < 0) ISCOs around Schwarzschild BH. The mass parameter M = 1 Fig. 7 Phase plot between the parameters S and C, shows the bound on the particle's spin S for different set of values of parameters a and C. The permitted region is between the curves formed by the dots.
Here, the blue and the grey regions represent corotating (J z > 0) and counter-rotating (J z < 0) cases of the spinning particle around rotating braneworld BH, respectively. The mass parameter M sets to unity Most importantly, we numerically probe the motion of spinning test particle in the rotating braneworld BH background, by solving the system of Eqs. (A.1)-(A.3) in general ((32)-(34)), where we study the behavior of ISCO parameters for different combinations of values of rotation a and tidal charge C parameters of both corotating and counter-rotating orbit cases around rotating braneworld BH as shown in Figs. 4 and 5. While considering the ISCO parameter r for the rotating braneworld BH, we showed that this ISCO parameter was always greater for counter-rotating orbits than that obtained for corotating orbits for the same value of rotation parameter a = 1.2.
Our analysis showed some interesting results about the Reissner-Nordström and Kerr-Newman BHs, which were not reported in earlier studies [64,67]. Such as (i) it is shown in Fig. 6 that the behavior of the particle spin parameter S as a function of parameter C is totally symmetric for corotating and counter-rotating orbits around Reissner-Nordström like BH (i.e., Q 2 = C > 0) under the change L z to −L z and S to −S. (ii) Also, the parameter S first decreases and then increases as the parameter C becomes smaller for Reissner-Nordström like BH. (iii) Similar to Reissner-Nordström like BH case the parameter S first decreases and then increases with decrease in parameter C, for the small values of rotation parameter a (say when a ≤ 0.6) as shown in Fig. 7, while the allowed range of parameter S always increases with decrease in parameter C when rotation parameter a ≥ 0.8 for the rotating braneworld BH. Here, it is worth reminding once again that rotating braneworld BH becomes Kerr-Newman like BH when parameter C > 0.
Acknowledgements P.S. would like to thanks Olivier Sarbach for useful discussions at some stage of this work. U. N. acknowledges support from PRODEP-SEP and the CONACYT thematic network project 280908 'Agujeros Negros y Ondas Gravitatorias' for financial support. R. B. and U. N. acknowledge support from CIC-UMNSH and SNI-CONACYT. P. S. would like to thank Programa de Desarrollo Profesional Docente (PRODEP) of the Secretaría de Educacíon Pública (SEP) of the Mexican government, for providing the financial support through post-doctoral program.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: This is a theoretical study and no experimental data has been listed.] 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 .

Appendix A: The explicit form of ISCO equations
Here, we present the explicit form of the ISCO Eqs. (32)- (34) in terms of x, a, C, J z and E, after doing the transformation x = 1/r (for the sake of simplicity): +E −24a 2 C J z Sx 5 + 20a 2 J z Sx 4 + 16aC 2 J z S 2 x 7 −28aC J z S 2 x 6 − 12aC J z S 2 x 5 −8aC J z x 3 + 12a J z S 2 x 5 + 10a J z S 2 x 4 +12a J z x 2 − 16C J z Sx 3 + 18J z Sx 2 − 4J z Sx +10a 2 C 2 S 4 x 9 − 18a 2 CS 4 x 8 + 12a 2 CS 2 x 5 It is worth mentioning here that working with parameter x instead of parameter r does not change the form of Eqs.