This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Spinning test particle motion around a rotating wormhole

Farrux Abdulxamidov farrukhabd63@gmail.com National University of Uzbekistan, Tashkent 100174, Uzbekistan Ulugh Beg Astronomical Institute, Astronomy St. 33, Tashkent 100052, Uzbekistan    Carlos A. Benavides-Gallego cabenavidesg20@shao.ac.cn Shanghai Astronomical Observatory, 80 Nandan Road, Shanghai 200030, P. R. China    Wen-Biao Han wbhan@shao.ac.cn Shanghai Astronomical Observatory, 80 Nandan Road, Shanghai 200030, P. R. China Hangzhou Institute for Advanced Study, University of Chinese Academy of Sciences, Hangzhou 310124, China School of Astronomy and Space Science, University of Chinese Academy of Sciences, Beijing 100049, China Shanghai Frontiers Science Center for Gravitational Wave Detection, 800 Dongchuan Road, Shanghai 200240, China    Javlon Rayimbaev javlon@astrin.uz Ulugh Beg Astronomical Institute, Astronomy St. 33, Tashkent 100052, Uzbekistan Akfa University, Kichik Halqa Yuli Street 17, Tashkent 100095, Uzbekistan National University of Uzbekistan, Tashkent 100174, Uzbekistan Tashkent State Technical University, Tashkent 100095, Uzbekistan    Ahmadjon Abdujabbarov ahmadjon@astrin.uz Shanghai Astronomical Observatory, 80 Nandan Road, Shanghai 200030, P. R. China National University of Uzbekistan, Tashkent 100174, Uzbekistan Ulugh Beg Astronomical Institute, Astronomy St. 33, Tashkent 100052, Uzbekistan Tashkent Institute of Irrigation and Agricultural Mechanization Engineers, Kori Niyoziy, 39, Tashkent 100000, Uzbekistan Institute of Nuclear Physics, Tashkent 100214, Uzbekistan
Abstract

In this work, we investigated the motion of spinning test particles around a rotating wormhole, extending, in this way, the previous work of Benavides-Gallego et al. in [Phys. Rev. D 101, no.12, 124024] to the general case. Using the Mathisson-Papapetrous-Dixon equations, we study the effective potential, circular orbits, and the innermost stable circular orbit (ISCO) of spinning test particles. We found that both the particle and wormhole spins affect the location of the ISCO significantly. On the other hand, Similar to the non-rotating case, we also found two possible configurations in the effective potential: plus and minus. Furthermore, the minimum value of the effective potential is not at the throat due to its spin aa, in contrast to the motion of the non-spinning test particles in a non-rotating wormhole, where the effective potential is symmetric, and its minimum value is at the throat of the wormhole. In the case of the ISCO, we found that it increases as the spin of the wormhole aa increases, in contrast to black holes where the presence of spin decreases the value of the ISCO. Finally, since the dynamical four-momentum and kinematical four-velocity of the spinning particle are not always parallel, we consider the superluminal bound, finding that the allowed values of ss change as the wormhole’s spin aa increases.

I Introduction

In 1935, Einstein and Rosen explored the possibility of an atomistic theory of matter and electricity without singularities Einstein:1935tc . The main idea behind their work was to consider the physical space as two identical sheets connected by a “bridge.” In this geometrical representation, particles were the very bridges connecting these sheets. Hence, using the metric tensor gμνg_{\mu\nu} of general relativity and the fields φμ\varphi_{\mu} of electromagnetism, they modified the gravitational equations, demonstrating that it is possible to obtain regular spherically symmetric solutions. These solutions are currently known as Einstein-Rosen bridges or “wormholes”.

The representation of particles as “bridges” could have been discovered back in 1916. According to G. W. Gibbons G.W.Gibbons2015 , previous to Einstein-Rosen’s paper, and a few months after the Schwarzschild solutions Schwarzschild:1916uq ; Schwarzschild:1916ae , Ludwig Flamm submitted a manuscript in which he explores some geometrical aspects of both exterior and interior solutions of Schwarzschild space-time L.Flamm1916 . In the former case (the exterior solution), Flamm was able to show that “the planar section is isometric to a surface of revolution, where the meridional curve is a parabola.” However, Flamm never contemplates the possibility of interpreting this result as a “bridge” connecting to regions of space-time. For this reason, we assume that the work of Einstein and Rosen gave birth to the modern study of “wormholes”.

Initially, the idea of “bridges” was considered more attractive than that of black holes until Wheeler and Fuller showed that the Einstein-Rosen bridge111Also known as the Schwarzschild wormhole. is unstable. Using a proper analysis, they discovered that the bridge would pinch-off in a finite time Fuller:1962zza . Therefore, Schwarzschild wormholes are not traversable. Nevertheless, the possibility of traversable wormholes was considered years latter in several works Ellis:1973yv ; Bronnikov:1973fh ; Morris:1988tu ; Morris:1988cz ; Visser:1989kh ; Visser:1989kg ; Visser:1989vq ; AzregAinou:1989wr ; Visser:1989am ; Poisson:1989zz ; Visser:1990wj ; Frolov:1990si ; Visser:1990wi . In particular, in Ref Morris:1988cz , Morris and Thorne used a different approach which allowed them to propose a set of “basic wormhole criteria.” Thus, by assuming the wormhole geometry, they used the Einstein field equations to compute physical quantities, such as the total energy density, the tension per unit area, and the pressure, in terms of what they call the “redshift” and “shape” functions. Their analysis establishes three properties to describe a wormhole. Firstly, the space-time is assumed to be static and spherically symmetric. Secondly, any wormhole solution should contain a throat that connects two asymptotically flat regions of space-time. This property is deeply related to the “shape function”. Finally, the solution must be horizonless. Since wormholes causally connect two different portions of the space-time by the throat, the presence of a horizon would prevent the two universes from being connected causally. This condition is satisfied by demanding the “redshift function” to be finite everywhere.

On the other hand, the most relevant conclusion in Morris and Thorne’s analysis is the necessity of “exotic matter” to generate a traversable wormhole. According to the authors, an observer passing through the throat with a radial velocity close to the speed of light would perceive negative energy. Therefore, from the classical point of view, traversable wormhole solutions violate the well-known energy conditions, established precisely to avoid negative energy densities. However, from the quantum point of view, there are some situations in which such violations may be physically valid. The quantum mechanical creation of particles, for example Martin-Moruno:2013sfa . In this sense, one can not entirely rule out the possibility of the existence of the exotic material required for the throat of a traversable wormhole to hold.

The idea of wormholes has been considered in different scenarios, such as Einstein’s gravity Echeverria:1991nk ; Deser:1992ts ; Deser:1993jr ; Hochberg:1997wp ; Abdujabbarov:2016efm ; Abdujabbarov:2009ad ; Pugliese:2011xn ; Willenborg:2018zsv ; Arganaraz:2019fup ; Stuchlik:2021tcn ; Benavides-Gallego:2021lqn and alternatives theories of gravity Moffat:1991xp ; Carlini:1992jda ; Bhawal:1992sz ; Letelier:1993cj ; Vollick:1998qf . In Ref. Benavides-Gallego:2021lqn , Benavides-Gallego et al. investigated the motion of spinning test particles around traversable wormholes. Using the Mathisson-Papapetrous-Dixon (MPD) equations, the authors computed the effective potential and showed that it is affected by the adimensional spin ss of the test particle. When s=0s=0, the particle follows the geodesic equation, and its effective potential is symmetric: its behavior is the same in both the lower and upper universes. Therefore, the innermost stable circular orbit (ISCO) is at the same distance from the wormhole’s throat.

On the other hand, when the spin s0s\neq 0, there is a different kind of symmetry which depends on the sing of the particle’s spin ss and \mathcal{L}, the particle’s dimensionless angular momentum. In that case, there are two possible configurations which define two different symmetries: “plus” and “minus”. In this sense, the effective potential has the same behavior if it has the same configuration. For example, a particle with negative ss and negative \mathcal{L} will have the same behavior as a particle with both ss and \mathcal{L} positive. In other words, the effective potential of a particle with ss parallel to the symmetry axis (s>0s>0) and moving counterclockwise (>0\mathcal{L}>0) is the same as that of a particle with ss antiparallel to the symmetry axis (s<0s<0) and moving clockwise around the wormhole (<0\mathcal{L}<0), see Fig. 3 of Ref. Benavides-Gallego:2021lqn . The authors also found a “mirror behavior” when the “plus” and “minus” configurations are considered at the same time. To explain this behavior, the authors consider two spinning particles moving counterclockwise with the same angular momentum \mathcal{L}, one in the lower universe (with s>0s>0) and the other in the upper universe (with s<0s<0). In this case, the effective potential has the same behavior for each particle, resembling the case of a non-spinning test particle, see the first panel in figure 3 of Ref. Benavides-Gallego:2021lqn .

The existence of two configurations affects the location of the ISCO. As mentioned before, in the case of non-spinning test particles, one finds a single value for lISCOl_{ISCO} located at the same distance from the wormhole’s throat in both universes222ll is the radial coordinate in Ref. Benavides-Gallego:2021lqn . The value l=0l=0 represents the wormhole’s throat.. However, the behavior is different when we consider spinning test particles. If |s|1|s|\geq 1, one finds only one possible value for the ISCO. On the other hand, if ss belongs to the interval 1<s<1-1<s<1, one encounters two possible values for the ISCO. One of these values is always closer to the wormhole’s throat.

One significant conclusion from Ref. Benavides-Gallego:2021lqn is the constrain obtained for the particle’s spin ss. It is well-known that the dynamical four-momentum pαp^{\alpha} and the kinematical four-velocity uαu^{\alpha} of a spinning test particle are not always parallel. In this sense, although pαpα=m2p_{\alpha}p^{\alpha}=-m^{2} holds, the normalization uαuα=1u_{\alpha}u^{\alpha}=-1 does not. Therefore, while the spinning particle moves closer to the center of symmetry, uαu^{\alpha} increases, and eventually, for certain values of the spin ss and radius ll, some components of the four-velocity may diverge. This means that for a certain value of ss, the particle’s trajectory changes from time-like to space-like, becoming, in this way, superluminal. From the physical point of view, the space-like motion does not have any meaning because the transition to uαuα>0u_{\alpha}u^{\alpha}>0 is not allowed for real particles. Therefore, one must impose an additional constrain defined by the relation uαuα=0u_{\alpha}u^{\alpha}=0, the superluminal bound. In the case of a non-rotating wormhole, spinning test particles will always move in the time-like region as long as |s|<1.5|s|<1.5 Benavides-Gallego:2021lqn .

In this paper, we consider the motion of spinning test particles around a rotating wormhole, generalizing in this way the results obtained by Benavides-Gallego et al. in Ref. Benavides-Gallego:2021lqn . We organize our work as follows. In Sec. II, we review the properties of rotating wormholes discussed in Ref. Teo:1998dp by E. Teo and its connection with the ideas of traversable wormholes described in Ref. Morris:1988cz . Next, in Sec. III, we discuss the motion of spinning particles in a static and axially symmetric space-time, where we use the MPD equations to obtain the effective potential and the superluminal condition. Then, in Sec. IV, we apply the results of Sec. III to Teo’s wormhole solution. We compute the effective potential, circular orbits, and the innermost stable circular orbit (ISCO). We also use the superluminal bound to find a constraint for the allowed spin of the particle. Finally, in Sec. V, we summarize our work and discuss the implications of the results. Throughout the manuscript, we use geometrized units setting G=c=1G=c=1 and b0=M=1b_{0}=M=1.

II Stationary, axisymmetric space-times

In Ref. Teo:1998dp , Teo follows the same paradigm as Morris-Thorne in Ref. Morris:1988cz , i.e. he first assumes the space-time geometry and then uses the Einstein field equations to deduce the form of the matter required to maintain the wormhole. Hence, the author begins by considering a stationary and axially symmetric space-time.

It is well-known that space-times are said to be stationary if it possesses a time-like Killing vector field ξα(/t)α\xi^{\alpha}\equiv(\partial/\partial t)^{\alpha}, which generates invariant time translations. On the other hand, a space-time is axisymmetric if it possesses a space-like killing vector field ψα(/φ)α\psi^{\alpha}\equiv(\partial/\partial\varphi)^{\alpha}, related to invariant rotations with respect to φ\varphi. Therefore, a space-time is stationary and axisymmetric if it possesses both ξα\xi^{\alpha} and ψα\psi^{\alpha} killing vector fields satisfying the following commutation relation Wald:1984rg ,

[ξ,ψ]=0.[\xi,\psi]=0. (1)

The commutativity of ξα\xi^{\alpha} and ψα\psi^{\alpha} in Eq. (1) allow us to choose a coordinate system in such a way that the Killing vectors represent the directions in which the space-time has symmetries Wald:1984rg . Hence, one can set x0=tx^{0}=t, x1=φx^{1}=\varphi, x2x^{2}, and x3x^{3} as a coordinate system. The stationary and axisymmetric character of the space-time requires the metric components to be independent of tt and φ\varphi. Therefore, the metric takes the form,

ds2=gμν(x2,x3)dxμdxν.ds^{2}=g_{\mu\nu}(x^{2},x^{3})dx^{\mu}dx^{\nu}. (2)

From the physical point of view, stationary and axisymmetric space-times have been of considerable interest in the study of black holes and stars since this kind of geometry describes the exterior gravitational field of rotating bodies, see Refs Hartle:1967he ; Hartle:1968si ; Thorne:1971R and references therein.

Thorne discusses the properties of static and axisymmetric space-times in Ref. Thorne:1971R . There, he stars by pointing out that (t,φ,x2,x3)(t,\varphi,x^{2},x^{3}) and (t,φ+2π,x2,x3)(t,\varphi+2\pi,x^{2},x^{3}) represent the same point. This behavior is due to the fact that φ\varphi is an angular coordinate about the rotation axis. As a consequence, φ\varphi belongs to the interval [0,2π)[0,2\pi).

On the other hand, because the space-time is axially symmetric, it must be invariant under a simultaneous inversion of tt and φ\varphi, i.e. the space-time does not change if ttt\rightarrow-t and φφ\varphi\rightarrow-\varphi. As a consequence, the metric coefficients gt2g_{t2}, gt3g_{t3}, gφ2g_{\varphi 2} and gφ3g_{\varphi 3} must vanish because they change the sign under simultaneous inversion of φ\varphi and tt. Therefore, the line element of Eq. (2) simplifies even more and reduces to Papapetrou:1966zz ; Carter:1969zz

ds2=g00dt2+2g01dtdφ+g11dφ2+gijdxidxj,ds^{2}=g_{00}dt^{2}+2g_{01}dtd\varphi+g_{11}d\varphi^{2}+g_{ij}dx^{i}dx^{j}, (3)

with i,j=2,3i,j=2,3. Here, the presence of the term g01g_{01} is related to the well-know dragging effect, see Appendix A.

The coordinates in Eq. (3) are uniquely determine up to a coordinate transformation of the form Thorne:1971R

x¯2=x¯2(x2,x3)andx¯3=x¯3(x2,x3).\begin{array}[]{ccc}\overline{x}^{2}=\overline{x}^{2}(x^{2},x^{3})&\text{and}&\overline{x}^{3}=\overline{x}^{3}(x^{2},x^{3}).\end{array} (4)

The freedom in such transformations can be used to simply the mathematics in the Einstein’s field equations or adapt the geometry to specific problems. Under transformation of the form given in Eq. (4), the components g00g_{00}, g01g_{01}, and g11g_{11} are invariant. Hence, g00=ξαξαg_{00}=\xi_{\alpha}\xi^{\alpha}, g01=ψαξαg_{01}=\psi_{\alpha}\xi^{\alpha}, and g11=ψαψαg_{11}=\psi_{\alpha}\psi^{\alpha} Wald:1984rg ; Thorne:1971R .

Finally, the space-time described in Eq. (3) must be asymptotically flat. This means that g001g_{00}\rightarrow 1, g011/rg_{01}\rightarrow 1/r and g11r2sin2θg_{11}\rightarrow r^{2}\sin^{2}\theta as rr\rightarrow\infty. The asymptotically-flatness of the line element (3) is important to define the mass, and the angular momentum333Here θ\theta and rr are the usual spherical coordinates, but not necessarily the same as x2x^{2} and x3x^{3}Thorne:1971R .

II.1 Canonical form of the rotating wormhole space-time

In the particular case of the rotating wormhole solution, by considering g22=g33=g11/sin2x2g_{22}=g_{33}=g_{11}/\sin^{2}x^{2} and g23=0g_{23}=0, the line element of Eq. (3) can be expressed as Teo:1998dp

ds2=N2dt2+eμdr2+r2K2[dθ2+sin2θ(dφωdt)2],ds^{2}=-N^{2}dt^{2}+e^{\mu}dr^{2}+r^{2}K^{2}\left[d\theta^{2}+\sin^{2}\theta(d\varphi-\omega dt)^{2}\right], (5)

where the functions444Also known as the four gravitational potentials. NN, μ\mu, KK and ω\omega only depend on the (x2,x3)(θ,r)(x^{2},x^{3})\equiv(\theta,r). Following Ref. Chandrasekhar:1985kt , it is possible to show the dragging effect in the space-time described by the line element in Eq. (5), see Appendix A.

Once the line element of a stationary and axisymmetric space-time is defined, Teo discusses some important features related to the metric. For example, he points out that the function K(r,θ)K(r,\theta) is a positive, nondecreasing function of rr, which he uses to define the “proper distanceRrK(r,θ)R\equiv rK(r,\theta) (with R/r>0\partial R/\partial r>0) measure at (r,θ)(r,\theta) from the origin. In this sense, one can interpret the value 2πRsinθ2\pi R\sin\theta as the proper circumference of the circle located at the point with coordinates (r,θ)(r,\theta).

On the other hand, the metric defined in Eq. (5) has the discriminant555From now on, we associate the coordinates tt, ϕ\phi, rr, θ\theta with the numbers, 0, 1, 2 and 3, respectively. Teo:1998dp ,

D2=gtφ2gttgφφ=(N(r,θ)K(r,θ)sinθ)2.D^{2}=g^{2}_{t\varphi}-g_{tt}g_{\varphi\varphi}=(N(r,\theta)K(r,\theta)\sin\theta)^{2}. (6)

According to Teo, the existence of horizons is determined by the function N(r,θ)N(r,\theta), which plays the role of the “redshift function”. Whenever, N=0N=0, D2=0D^{2}=0, implying the presence of an event horizon. Therefore, to avoid a singular behavior of the metric (D20D^{2}\neq 0) on the rotation axis θ=0\theta=0 and θ=π\theta=\pi, Teo imposes the regularity conditions on the gravitational potential. These conditions state that the partial derivatives with respect to θ\theta of N(r,θ)N(r,\theta), μ(r,θ)\mu(r,\theta), and K(r,θ)K(r,\theta) must vanish on the rotation axis. The regularity conditions are, in this sense, crucial to establish the wormhole geometry in the line element of Eq. (5) since wormholes, according to Ref. Morris:1988cz , must be horizonless.

In analogy to Morris-Thorne function μ\mu, Teo defines a similar function by including the dependence on θ\theta. Therefore, μ(r,θ)\mu(r,\theta) is given by Teo:1998dp

μ(r,θ)=ln(1b(r,θ)r);\mu(r,\theta)=-\ln\left(1-\frac{b(r,\theta)}{r}\right); (7)

with b(r,θ)b(r,\theta) playing the role of the “shape function”. In this sense, the radial coordinate must be constrained to rbr\geq b, where the throat is located at r=br=b. In this way, Eq. (5) reduces to the Morris-Thorne case when there is no rotation, e. i. N(r,θ)eΦ(r)N(r,\theta)\rightarrow e^{\Phi}(r), b(r,θ)b(r)b(r,\theta)\rightarrow b(r), K(r,θ)1K(r,\theta)\rightarrow 1 and ω(r,θ)0\omega(r,\theta)\rightarrow 0. Moreover, Teo assumes the gravitational potentials to be well-behaved at the wormhole’s throat. The reason for such an assumption has to do with the singularity-free behavior at the throat. If one computes the curvature scalar of the metric in Eq. (5) (evaluated at the throat), it is possible to see that it has the form Harko:2009xf

\displaystyle\mathcal{R} =1r2K2(μθθ+12μθ2)μθNr2K2(Nsinθ)θsinθ\displaystyle=-\frac{1}{r^{2}K^{2}}\left(\mu_{\theta\theta}+\frac{1}{2}\mu^{2}_{\theta}\right)-\frac{\mu_{\theta}}{Nr^{2}K^{2}}\frac{(N\sin\theta)_{\theta}}{\sin\theta} (8)
2Nr2K2(Nθsinθ)θsinθ2r2K3(Kθsinθ)θsinθ\displaystyle-\frac{2}{Nr^{2}K^{2}}\frac{(N_{\theta}\sin\theta)_{\theta}}{\sin\theta}-\frac{2}{r^{2}K^{3}}\frac{(K_{\theta}\sin\theta)_{\theta}}{\sin\theta}
eμμr[ln(Nr2K2)]r+sin2θωθ22N2+2r2K4(K2+Kθ2),\displaystyle e^{-\mu}\mu_{r}[\ln(Nr^{2}K^{2})]_{r}+\frac{\sin^{2}\theta\omega^{2}_{\theta}}{2N^{2}}+\frac{2}{r^{2}K^{4}}(K^{2}+K^{2}_{\theta}),

where the subscripts denote partial derivatives with respect to θ\theta and rr. From the last equation, one can see that \mathcal{R} could have a singular behavior due to terms Teo:1998dp ; Harko:2009xf

μθθ+12μθ2\displaystyle\mu_{\theta\theta}+\frac{1}{2}\mu^{2}_{\theta} =bθθrb+32bθ2(rb)2,\displaystyle=\frac{b_{\theta\theta}}{r-b}+\frac{3}{2}\frac{b^{2}_{\theta}}{(r-b)^{2}}, (9)
μθ\displaystyle\mu_{\theta} =bθrb.\displaystyle=\frac{b_{\theta}}{r-b}.

Therefore, to avoid singularities in the curvature scalar at the throat, it is necessary that bθ=bθθ=0b_{\theta}=b_{\theta\theta}=0. Ergo, the throat is at a constant value of rr.

As mentioned before, if one wants the line element in Eq. (5) to describe a wormhole, it is crucial to satisfy the so-called “flare-out” condition at the throat. Following the same process described in Ref. Morris:1988cz , Teo embeds the space-time in a higher-dimensional space by considering a constant value of θ\theta in a slice of constant tt, what Morris-Thorne think of as a picture of the whole space-time at a fixed moment tt.

After embedding the metric, Teo found the following “flare-out” condition at the throat Teo:1998dp

d2rdz2=bbrr2b2>0,\frac{d^{2}r}{dz^{2}}=\frac{b-b_{r}r}{2b^{2}}>0, (10)

which is the same condition as in the Morris-Throne wormhole Morris:1988cz . Since bθ=0b_{\theta}=0, it is possible to define a new radial coordinate l2=r2+b2l^{2}=r^{2}+b^{2} in the vicinity of the throat, satisfying the relation

dldr±(1br)1/2.\frac{dl}{dr}\equiv\pm\left(1-\frac{b}{r}\right)^{-1/2}. (11)

Hence, in the immediate vicinity of the throat666To first order in rr0r-r_{0}, with r0r_{0} the location of the throat. the line element in Eq. (1), reduces to Teo:1998dp

ds2\displaystyle ds^{2} =N2(l,θ)dt2+dl2+r2(l)K2(l,θ)\displaystyle=-N^{2}(l,\theta)dt^{2}+dl^{2}+r^{2}(l)K^{2}(l,\theta) (12)
×[dθ2+sin2θ(dφω(l,θ)dt)2].\displaystyle\times[d\theta^{2}+\sin^{2}\theta(d\varphi-\omega(l,\theta)dt)^{2}].

The metric expressed in this way smoothly connects two asymptotic regions of the space-time across the throat, in contrast to Eq. (1), where the radial coordinate rr is singular. If the shape function does not depend on θ\theta, Eq. (11) is valid everywhere and the coordinate ll takes the range (,)(-\infty,\infty). As a consequence, the metric in Eq. (11) covers the whole space-time, and we can assume the wormhole throat is at l=0l=0 and define the upper universe when l>0l>0, and the lower universe when l<0l<0.

According to Morris-Thorne, one can use the space-time in Eq. (1) to compute the non-vanishing components of the stress-energy tensor. To do so, one needs to consider a “local Lorentz frame”, where physical “observations” are performed by a local observer, who remains at rest with respect to the coordinate system (t,θ,r,φ)(t,\theta,r,\varphi). In this frame, the components are 777Here we use the same notation as in Ref. Chandrasekhar:1985kt . T(t)(t)T_{(t)(t)}, T(t)(φ)T_{(t)(\varphi)}, T(φ)(φ)T_{(\varphi)(\varphi)}, and T(i)(j)T_{(i)(j)}. These components have the usual physical meaning. For example, T(t)(t)T_{(t)(t)} is the mass-energy density, while T(t)(φ)T_{(t)(\varphi)} represents the rotation of the matter distribution. In Ref. Harko:2009xf , Harko et al. obtained the expressions evaluated at the throat.

Finally, using the null energy condition Wald:1984rg

Rαβκακβ0,R_{\alpha\beta}\kappa^{\alpha}\kappa^{\beta}\geq 0, (13)

where RαβR_{\alpha\beta} is the Ricci tensor and κα\kappa^{\alpha} a null vector given by Teo:1998dp

κα=(1N,eμ/2,0,ωN),\kappa^{\alpha}=\left(\frac{1}{N},-e^{-\mu/2},0,\frac{\omega}{N}\right), (14)

Teo have found that Teo:1998dp

Rαβκακβ\displaystyle R_{\alpha\beta}\kappa^{\alpha}\kappa^{\beta} =eμμr(rK)rrKωθ2sin2θ2N214μθ2(rK)2\displaystyle=e^{-\mu}\mu_{r}\frac{(rK)_{r}}{rK}-\frac{\omega^{2}_{\theta}\sin^{2}\theta}{2N^{2}}-\frac{1}{4}\frac{\mu^{2}_{\theta}}{(rK)^{2}} (15)
12(μθsinθ)θ(rK)2sinθ+(Nθsinθ)θ(rK)2Nsinθ<0.\displaystyle-\frac{1}{2}\frac{(\mu_{\theta}\sin\theta)_{\theta}}{(rK)^{2}\sin\theta}+\frac{(N_{\theta}\sin\theta)_{\theta}}{(rK)^{2}N\sin\theta}<0.

Nevertheless, He remarks that by choosing NN, and μ\mu appropriately, RαβR_{\alpha\beta} could be positive at some point in the interval (0,π)(0,\pi). Consequently, it is possible to move the exotic matter that supports the wormhole around the throat so that an infalling observer would not encounter it. One example of such space-time is given by considering Teo:1998dp

N\displaystyle N =K=1+(4Jcosθ)2r,\displaystyle=K=1+\frac{(4J\cos\theta)^{2}}{r}, (16)
b\displaystyle b =b02r,\displaystyle=\frac{b^{2}_{0}}{r},
μ\displaystyle\mu =ln(1br)\displaystyle=-\ln\left(1-\frac{b}{r}\right)
ω\displaystyle\omega =2Jr3.\displaystyle=\frac{2J}{r^{3}}.

III Equations of motion

In this section, we review the theoretical background necessary to investigate the motion of spinning test particles. It is well-known that this problem was considered for the first time by Mathisson in 1937 when he studied the problem of extended bodies in general relativity (GR). According to Mathisson, the motion of spinning test particles does not follow the usual geodesic equation of GR due to the coupling between the Riemann curvature tensor and the spin of the moving particle Mathisson:1937zz . Papapetrou considers the same problem in Refs. Papapetrou:1951pa ; Corinaldesi:1951pb , where he established a similar approach. Later, Tulczyjew improved on the methods of Mathisson tulczyjew1959motion ; BWTulzcyjew1962 while Moller and others made improvements in the definition of center-of-mass moller1949definition ; beiglbock1967center ; dixon1964covariant ; Dixon:1970zza ; Dixon:1970zz ; ehlers1977dynamics . Today, the equations that describe the motion of extended bodies with spin and mass are known as the Mathisson-Papapetrous-Dixon (MPD) equations. Recently, some authors modified the MPD equations, see Refs. Deriglazov:2017jub ; Deriglazov:2018vwa .

Mathematically, the MPD equations are given by

Dpαdλ\displaystyle\frac{Dp^{\alpha}}{d\lambda} =12RβδσαuβSδσ,\displaystyle=-\frac{1}{2}R^{\alpha}_{\;\;\beta\delta\sigma}u^{\beta}S^{\delta\sigma}, (17)
DSαβdλ\displaystyle\frac{DS^{\alpha\beta}}{d\lambda} =pαuβpβuα,\displaystyle=p^{\alpha}u^{\beta}-p^{\beta}u^{\alpha},

where D/dλuααD/d\lambda\equiv u^{\alpha}\nabla_{\alpha} is the projection of the covariant derivative along the particle’s trajectory, uμ=dxμ/dλu^{\mu}=dx^{\mu}/d\lambda is the 4-velocity of the test particle, pαp^{\alpha} is the canonical 4-momentum, RβδσαR^{\alpha}_{\;\;\beta\delta\sigma} is the Riemann curvature tensor, and λ\lambda is an affine parameter. The second rank tensor SαβS^{\alpha\beta} is antisymmetric, Sαβ=SβαS^{\alpha\beta}=-S^{\beta\alpha}.

Expressed in that form, Eq. (17) shows us the coupling between the Riemann tensor RβδσαR^{\alpha}_{\;\;\beta\delta\sigma} and the spin of the moving particle. To see this coupling clearly, let’s consider the well-known geodesic equation of GR

uββuα+Γσβαuσuβ=0.u^{\beta}\partial_{\beta}u^{\alpha}+\Gamma^{\alpha}_{\;\;\;\sigma\beta}u^{\sigma}u^{\beta}=0. (18)

In terms of the particle 4-momentum and the projection of the covariant derivative along the particle’s trajectory, Eq. (18) reduces to

Dpαdλ=0.\frac{Dp^{\alpha}}{d\lambda}=0. (19)

Therefore, comparing Eqs. (17) and (19), one can see how the interaction between the Riemann curvature tensor and the antisymmetric tensor SαβS^{\alpha\beta} does affect the motion of spinning test particles in curved space-times.

A crucial aspect in the MPD equations is related to the center of mass of the spinning test particle. In this sense, to solve the system in Eq. (17), one needs to fix its center of mass. This is done by including the condition tulczyjew1959motion ; Saijo:1998mn

Sαβpα=0.S^{\alpha\beta}p_{\alpha}=0. (20)

This condition is known as the Tulczyjew Spin Supplementary Condition (SSC) Saijo:1998mn . From Eq. (20), the canonical momentum and the spin of the particle provide two independent conserved quantities given by the relations

pαpα\displaystyle p^{\alpha}p_{\alpha} =m2,\displaystyle=-m^{2}, (21)
SαβSαβ\displaystyle S^{\alpha\beta}S_{\alpha\beta} =2S2\displaystyle=2S^{2}

Nevertheless, in contrast to the spinning test particle’s canonical momentum conservation, it is important to point out that the squared velocity does not necessarily satisfy the condition

uαuα=1,u_{\alpha}u^{\alpha}=-1, (22)

because the 4-vectors pαp^{\alpha} and uαu^{\alpha} are not always parallel. In this sense, to ensure that the particle’s 4-velocity is always smaller than the speed of light, one needs to impose an additional condition: the superluminal bound, see Sec. III.2.

Additionally, to the conserved quantities resulting from the Tulczyjew-SSC condition, one also has the conserved quantities associated to the space-time symmetries. As remarked in Sec.II, the line element in Eq. (5) has associated two Killing vector fields. One generating invariant time translations (ξα\xi^{\alpha}) and the other related to invariant rotations with respect to φ\varphi (ψα\psi^{\alpha}). Hence, the conserved quantities associated to them can be obtained from the following relation,

pακα12Sαββκα=pακα12Sαββκα=constant,p^{\alpha}\kappa_{\alpha}-\frac{1}{2}S^{\alpha\beta}\nabla_{\beta}\kappa_{\alpha}=p^{\alpha}\kappa_{\alpha}-\frac{1}{2}S^{\alpha\beta}\partial_{\beta}\kappa_{\alpha}=\text{constant}, (23)

where κα\kappa^{\alpha} is the Killing vector field. In the last expression, we used the fact that the term SαβΓβαγS^{\alpha\beta}\Gamma^{\gamma}_{\;\;\beta\alpha} in the covariant derivative SαββκαS^{\alpha\beta}\nabla_{\beta}\kappa_{\alpha} vanishes because SαβS^{\alpha\beta} is antisymmetric while Γβαγ\Gamma^{\gamma}_{\;\;\beta\alpha} is symmetric.

III.1 The effective potential

In the case of static axially symmetric space-times, Toshmatov and Malafarina obtained the most general expression for VeffV_{\text{eff}} in terms of the metric components gμνg_{\mu\nu} and its derivatives888Where denotes the derivative with respect to the radial coordinate gμνg^{\prime}_{\mu\nu} Ref. Toshmatov:2019bda . Here, we review and discuss the most important steps in the calculation.

Let’s start considering the line element of a stationary and axisymmetric space-time. As mentioned Sec. II, Eq. (3) give us its general form

ds2=gttdt2+grrdr2+2gtφdtdφ+gθθdθ2+gφφdφ2.ds^{2}=g_{tt}dt^{2}+g_{rr}dr^{2}+2g_{t\varphi}dtd\varphi+g_{\theta\theta}d\theta^{2}+g_{\varphi\varphi}d\varphi^{2}. (24)

As we pointed out above, due to the Killing vectors ξα=δtα\xi^{\alpha}=\delta^{\alpha}_{t} and ψα=δφα\psi^{\alpha}=\delta^{\alpha}_{\varphi}, this space-time has associated two constants of motion: the energy EE and the angular momentum LL. Hence, after using Eq. (23), we obtain999There is a typo in the second relation of Eq. (11) of Ref. Toshmatov:2019bda . It should be SαβS^{\alpha\beta}.

E\displaystyle-E =\displaystyle= pt12gtα,βSαβ\displaystyle p_{t}-\frac{1}{2}g_{t\alpha,\beta}S^{\alpha\beta}
=\displaystyle= pt12(gtt,rStr+gφt,rSφr),\displaystyle p_{t}-\frac{1}{2}\left(g_{tt,r}S^{tr}+g_{\varphi t,r}S^{\varphi r}\right),
L\displaystyle L =\displaystyle= pφ12gtα,βSαβ\displaystyle p_{\varphi}-\frac{1}{2}g_{t\alpha,\beta}S^{\alpha\beta}
=\displaystyle= pφ12(gtφ,rStr+gφφ,rSϕr).\displaystyle p_{\varphi}-\frac{1}{2}\left(g_{t\varphi,r}S^{tr}+g_{\varphi\varphi,r}S^{\phi r}\right).

In the last system of equations, we assume that the particle’s motion is constrained to the equatorial plane (θ=π/2\theta=\pi/2). Due to this assumption, the metric functions depend only on the radial coordinate and pθ=0p_{\theta}=0. Furthermore, since Sθα=0S^{\theta\alpha}=0, the number of independent components of the spin tensor is reduced to three, i.e. StrS^{tr}, StφS^{t\varphi}, and SrφS^{r\varphi}.

To solve the system in Eq. (III.1), it is necessary to express the components StφS^{t\varphi} and SrφS^{r\varphi} in terms of StrS^{tr}. To do so, one uses the Tulczyjew-SSC contidion (20), from which one obtains

Stφ\displaystyle S^{t\varphi} =prptSφr=prpφStr\displaystyle=\frac{p_{r}}{p_{t}}S^{\varphi r}=-\frac{p_{r}}{p_{\varphi}}S^{tr} (27)
Srφ\displaystyle S^{r\varphi} =ptpφSrt=ptpφStr.\displaystyle=-\frac{p_{t}}{p_{\varphi}}S^{rt}=\frac{p_{t}}{p_{\varphi}}S^{tr}.

Using these relations, Eq. (III.1) takes the form

E\displaystyle-E =pt12(gttpφgφtpt)Strpφ\displaystyle=p_{t}-\frac{1}{2}\left(g^{\prime}_{tt}p_{\varphi}-g^{\prime}_{\varphi t}p_{t}\right)\frac{S^{tr}}{p_{\varphi}} (28)
L\displaystyle L =pφ12(gtφpφgφφpt)Strpφ.\displaystyle=p_{\varphi}-\frac{1}{2}\left(g^{\prime}_{t\varphi}p_{\varphi}-g^{\prime}_{\varphi\varphi}p_{t}\right)\frac{S^{tr}}{p_{\varphi}}.

Now, from the spin coservation and the normalization conditions in Eq. (21), we have that

Str=pφsgrr(gtφ2gφφgtt).S^{tr}=\frac{p_{\varphi}s}{\sqrt{g_{rr}(g^{2}_{t\varphi}-g_{\varphi\varphi}g_{tt})}}. (29)

Here, s=S/ms=S/m represents the specific angular momentum of the particle, which can be positive or negative with respect to the direction of pφp_{\varphi}. After replacing Eq. (29) into Eq. (28) and solving the system for ptp_{t} and pφp_{\varphi}, one obtains Toshmatov:2019bda

pt\displaystyle p_{t} =E+s(AL+BE)1Ds2\displaystyle=\frac{-E+s(AL+BE)}{1-Ds^{2}} (30)
pφ\displaystyle p_{\varphi} =L+s(BL+CE)1Ds2.\displaystyle=\frac{L+s(BL+CE)}{1-Ds^{2}}.

with Toshmatov:2019bda

A\displaystyle A =gtt2grr(gtφ2gφφgtt),\displaystyle=\frac{g^{\prime}_{tt}}{2\sqrt{g_{rr}(g^{2}_{t\varphi}-g_{\varphi\varphi}g_{tt})}}, (31)
B\displaystyle B =gtφ2grr(gtφ2gφφgtt),\displaystyle=\frac{g^{\prime}_{t\varphi}}{2\sqrt{g_{rr}(g^{2}_{t\varphi}-g_{\varphi\varphi}g_{tt})}},
C\displaystyle C =gφφ2grr(gtφ2gφφgtt),\displaystyle=\frac{g^{\prime}_{\varphi\varphi}}{2\sqrt{g_{rr}(g^{2}_{t\varphi}-g_{\varphi\varphi}g_{tt})}},
D\displaystyle D =(gtφ)2gttgφφ4grr(gtφ2gφφgtt).\displaystyle=\frac{(g^{\prime}_{t\varphi})^{2}-g^{\prime}_{tt}g^{\prime}_{\varphi\varphi}}{4g_{rr}(g^{2}_{t\varphi}-g_{\varphi\varphi}g_{tt})}.

Now, from the normalization condition in Eq. (21), we obtain Toshmatov:2019bda

pr2=grr(gttpt2gφφpφ22gtφptpφm2).p^{2}_{r}=g_{rr}(-g^{tt}p^{2}_{t}-g^{\varphi\varphi}p^{2}_{\varphi}-2g^{t\varphi}p_{t}p_{\varphi}-m^{2}). (32)

After substituting Eqs. (30) and (21) into Eq. (32), ones obtains a second-degree polynomial equation for the energy EE

pr2=βα(E2+δLβE+σL2βρβ),p^{2}_{r}=\frac{\beta}{\alpha}\left(E^{2}+\frac{\delta L}{\beta}E+\frac{\sigma L^{2}}{\beta}-\frac{\rho}{\beta}\right), (33)

which can be expressed as

pr2=βα(EVeff+)(EVeff),p^{2}_{r}=\frac{\beta}{\alpha}(E-V_{\rm eff}^{+})(E-V_{\rm eff}^{-}), (34)

where

Veff±=δL2β±(δL2β)2+(ρβσL2β).V_{\rm eff}^{\pm}=-\frac{\delta L}{2\beta}\pm\sqrt{\left(\frac{\delta L}{2\beta}\right)^{2}+\left(\frac{\rho}{\beta}-\frac{\sigma L^{2}}{\beta}\right)}. (35)

The expressions for α\alpha, β\beta, δ/2\delta/2, σ\sigma and ρ\rho are given in Ref. Toshmatov:2019bda . Nevertheless, it is important to note that we found some typos in the expressions defining α\alpha, δ\delta and σ\sigma. In the definition of α\alpha, it should be grrg^{rr} instead grrg_{rr}. In the expression for δ\delta, the second term in the numerator where s2s^{2} is the factor, one should have gtφgttg^{t\varphi}g^{\prime}_{tt} instead of gtφgttg^{t\varphi}g_{tt}. Finally, in the expression defining σ\sigma, the minus sign in the second term of the numerator, where ss is the factor, one should has +gφφgtφ+g^{\varphi\varphi}g^{\prime}_{t\varphi} instead of gφφgtφ-g^{\varphi\varphi}g^{\prime}_{t\varphi}. Hence, we obtain the following expressions

α\displaystyle\alpha =grr(1Ds2)2,\displaystyle=g^{rr}(1-Ds^{2})^{2}, (36)
β\displaystyle\beta =gtt+s(gtφgφφ+gttgtφ)grr(gtφ2gφφgtt)s2[gtt(gtφ)2+gφφ(gφφgφφ+2gtφgtφ)]4grr(gtφ2gφφgtt),\displaystyle=-g^{tt}+\frac{s(g^{t\varphi}g^{\prime}_{\varphi\varphi}+g^{tt}g^{\prime}_{t\varphi})}{\sqrt{g_{rr}(g^{2}_{t\varphi}-g_{\varphi\varphi}g_{tt})}}-\frac{s^{2}\left[g^{tt}(g^{\prime}_{t\varphi})^{2}+g^{\prime}_{\varphi\varphi}(g^{\varphi\varphi}g^{\prime}_{\varphi\varphi}+2g^{t\varphi}g^{\prime}_{t\varphi})\right]}{4g_{rr}(g^{2}_{t\varphi}-g_{\varphi\varphi}g_{tt})},
δ2\displaystyle\frac{\delta}{2} =gtφ+s(gttgttgφφgφφ)2grr(gtφ2gφφgtt)s2[gtφ(gttgtt+gtφgtφ)+gφφ(gtφgtt+gφφgtφ)]4grr(gtφ2gφφgtt),\displaystyle=g^{t\varphi}+\frac{s(g^{tt}g^{\prime}_{tt}-g^{\varphi\varphi}g^{\prime}_{\varphi\varphi})}{2\sqrt{g_{rr}(g^{2}_{t\varphi}-g_{\varphi\varphi}g_{tt})}}-\frac{s^{2}\left[g^{\prime}_{t\varphi}(g^{tt}g^{\prime}_{tt}+g^{t\varphi}g^{\prime}_{t\varphi})+g^{\prime}_{\varphi\varphi}(g^{t\varphi}g^{\prime}_{tt}+g^{\varphi\varphi}g^{\prime}_{t\varphi})\right]}{4g_{rr}(g^{2}_{t\varphi}-g_{\varphi\varphi}g_{tt})},
σ\displaystyle\sigma =gφφs(gtφgtt+gφφgtφ)grr(gtφ2gφφgtt)s2[gtt(gtt)2+gtφ(gφφgtφ+2gtφgtt)]grr(gtφ2gφφgtt),\displaystyle=-g^{\varphi\varphi}-\frac{s(g^{t\varphi}g^{\prime}_{tt}+g^{\varphi\varphi}g^{\prime}_{t\varphi})}{\sqrt{g_{rr}(g^{2}_{t\varphi}-g_{\varphi\varphi}g_{tt})}}-\frac{s^{2}\left[g^{tt}(g^{\prime}_{tt})^{2}+g^{\prime}_{t\varphi}(g^{\varphi\varphi}g^{\prime}_{t\varphi}+2g^{t\varphi}g^{\prime}_{tt})\right]}{\sqrt{g_{rr}(g^{2}_{t\varphi}-g_{\varphi\varphi}g_{tt})}},
ρ\displaystyle\rho =m2(1Ds2)2.\displaystyle=m^{2}(1-Ds^{2})^{2}.

III.2 Superluminal bound

In this section, we focus our attention on the well-known superluminal bound. This constraint will be crucial to find the allowed values of ss to keep a spinning test particle moving in a trajectory with physical meaning. In this sense, we will obtain analytic expressions in the case of static and axially symmetric space-times. Then, in Sec. IV, we will use these formulas in the particular case of a traversable rotating wormhole.

As we pointed out above, although pαpα=m2p_{\alpha}p^{\alpha}=-m^{2} is satisfied, the normalization uαuα=1u_{\alpha}u^{\alpha}=-1 does not necessarily hold because the four-momentum and the four-velocity are not always parallel to each other. Therefore, as the spinning test particle moves closer to the center of symmetry, the four-velocity uαu^{\alpha} increases, and for specific values of spin ss and radius, some components of uαu^{\alpha} may diverge. In this sense, the motion of the spinning particle crosses the boundary between time-like and space-like trajectories, becoming in this way superluminal.

Particles moving in a space-like trajectory (superluminal motion) do not have a physical meaning. Therefore, the transition to uαuα>0u_{\alpha}u^{\alpha}>0 is not allowed for real particles. As a consequence, one must impose a further constraint: the superluminal bound, defined by the condition uαuα=0u_{\alpha}u^{\alpha}=0. Hence, to keep spinning test particles moving in space-like trajectories, it is necessary to impose the following constrain Toshmatov:2019bda ; Conde:2019juj

uαuα(ut)2=gtt+grr(drdt)2+2gtφdφdt+gφφ(dφdt)20.\frac{u_{\alpha}u^{\alpha}}{(u^{t})^{2}}=g_{tt}+g_{rr}\left(\frac{dr}{dt}\right)^{2}+2g_{t\varphi}\frac{d\varphi}{dt}+g_{\varphi\varphi}\left(\frac{d\varphi}{dt}\right)^{2}\leq 0. (37)

To consider the superluminal bound in our investigation, we need to obtain analytical expressions for dr/dtdr/dt and dφ/dtd\varphi/dt. In order to do so, we follow a method proposed by Hojman and Asenjo in Ref. Hojman:2012me .

In Sec.III.1, we constrain our calculations to the equatorial plane (θ=π/2\theta=\pi/2). Ergo, the non-vanishing components of SαβS^{\alpha\beta} are StrS^{tr}, StφS^{t\varphi} and SrφS^{r\varphi}. Hence, from the second MPD equation in Eq. (17), we obtain the following system of equations

DStrdλ\displaystyle\frac{DS^{tr}}{d\lambda} =pturutpr,\displaystyle=p^{t}u^{r}-u^{t}p^{r}, (38)
DStφdλ\displaystyle\frac{DS^{t\varphi}}{d\lambda} =ptuφutpφ,\displaystyle=p^{t}u^{\varphi}-u^{t}p^{\varphi},
DSφrdλ\displaystyle\frac{DS^{\varphi r}}{d\lambda} =pφuruφpr.\displaystyle=p^{\varphi}u^{r}-u^{\varphi}p^{r}.

In order to solve the last system, we note that it is possible to reduce it into a system of two equation by expressing the first and second equations in terms of SφrS^{\varphi r}. To do so, one needs to consider the Tulczyjew-SSC condition Sαβpα=0S^{\alpha\beta}p_{\alpha}=0. In the case of StrS^{tr}, for example, Eq. (20) reduces to

Strpt+Sφrpφ=0.S^{tr}p_{t}+S^{\varphi r}p_{\varphi}=0. (39)

Then, after applying the operator D/dλD/d\lambda, using the second relation in Eq. (27), solving for DStr/dλDS^{tr}/d\lambda, and replacing in the first equation of Eq. (38), we obtain Benavides-Gallego:2021lqn

DStrdλ\displaystyle\frac{DS^{tr}}{d\lambda} =\displaystyle= Sφrpt(pφptDptdλDpφdλ)pφptDSφrdλ\displaystyle\frac{S^{\varphi r}}{p_{t}}\left(\frac{p_{\varphi}}{pt}\frac{Dp_{t}}{d\lambda}-\frac{Dp_{\varphi}}{d\lambda}\right)-\frac{p_{\varphi}}{p_{t}}\frac{DS^{\varphi r}}{d\lambda} (40)
=\displaystyle= pturutpr.\displaystyle p^{t}u^{r}-u^{t}p^{r}.

We proceed similarly with DStφ/dλDS^{t\varphi}/d\lambda in Eq. (38), obtaining Benavides-Gallego:2021lqn

DStφdλ\displaystyle\frac{DS^{t\varphi}}{d\lambda} =\displaystyle= prptDSφrdλ+Sφrpt(DprdλprptDptdλ)\displaystyle\frac{p_{r}}{p_{t}}\frac{DS^{\varphi r}}{d\lambda}+\frac{S^{\varphi r}}{p_{t}}\left(\frac{Dp_{r}}{d\lambda}-\frac{p_{r}}{p_{t}}\frac{Dp_{t}}{d\lambda}\right) (41)
=\displaystyle= ptuφutpφ.\displaystyle p^{t}u^{\varphi}-u^{t}p^{\varphi}.

With the MPD equation for StrS^{tr} and StφS^{t\varphi} expressed in terms of DSφr/dλDS^{\varphi r}/d\lambda, we can now use the MPD for SφrS^{\varphi r} (third equation in Eq. (38)) to reduce the system from three to only two equations. This new system is given by Benavides-Gallego:2021lqn

Sφrpt(pφDptdλptDpφdλ)\displaystyle\frac{S^{\varphi r}}{p_{t}}\left(p_{\varphi}\frac{Dp_{t}}{d\lambda}-p_{t}\frac{Dp_{\varphi}}{d\lambda}\right) =ur(pφpφ+ptpt)\displaystyle=u^{r}\left(p_{\varphi}p^{\varphi}+p_{t}p^{t}\right) (42)
pr(pφuφ+ptut),\displaystyle-p^{r}(p_{\varphi}u^{\varphi}+p_{t}u^{t}),
Sφrpt(ptDprdλprDptdλ)\displaystyle\frac{S^{\varphi r}}{p_{t}}\left(p_{t}\frac{Dp_{r}}{d\lambda}-p_{r}\frac{Dp_{t}}{d\lambda}\right) =uφ(prpr+ptpt)\displaystyle=u^{\varphi}\left(p_{r}p^{r}+p_{t}p^{t}\right)
pφ(prur+ptut).\displaystyle-p^{\varphi}\left(p_{r}u^{r}+p_{t}u^{t}\right).

Now, using the first MPD equation in Eq. (38), we obtain the following relations (see Appendix B)

Dptdλ\displaystyle\frac{Dp_{t}}{d\lambda} =Sφrpt[(pφRtrtrptRtrφr)urprRtφtφuφ],\displaystyle=\frac{S^{\varphi r}}{p_{t}}\left[\left(p_{\varphi}R_{trtr}-p_{t}R_{tr\varphi r}\right)u^{r}-p_{r}R_{t\varphi t\varphi}u^{\varphi}\right], (43)
Dprdλ\displaystyle\frac{Dp_{r}}{d\lambda} =[(pφRrttr+ptRrtrφ)ut+(pφRrφtr+ptRrφrφ)uφ]\displaystyle=\left[(p_{\varphi}R_{rttr}+p_{t}R_{rtr\varphi})u^{t}+(p_{\varphi}R_{r\varphi tr}+p_{t}R_{r\varphi r\varphi})u^{\varphi}\right]
×Sφrpt\displaystyle\times\frac{S^{\varphi r}}{p_{t}}
Dpφdλ\displaystyle\frac{Dp_{\varphi}}{d\lambda} =Sφrpt[ur(pφRφrtr+ptRφrrφ)prRφttφut].\displaystyle=\frac{S^{\varphi r}}{p_{t}}\left[u^{r}(p_{\varphi}R_{\varphi rtr}+p_{t}R_{\varphi rr\varphi})-p_{r}R_{\varphi tt\varphi}u^{t}\right].

After replacing Eq. (43) into Eq. (42), we obtain the following system

ur[pφ2𝒜+2pφpt𝒟+pt2]=\displaystyle u^{r}\left[p^{2}_{\varphi}\mathcal{A}+2p_{\varphi}p_{t}\mathcal{D}+p^{2}_{t}\mathcal{B}\right]= uφ𝒞pφpr+ut𝒞ptpr,\displaystyle u^{\varphi}\mathcal{C}p_{\varphi}p_{r}+u^{t}\mathcal{C}p_{t}p_{r}, (44)
uφ[pt2+ptpφ𝒟+pr2𝒞]=\displaystyle u^{\varphi}\left[p^{2}_{t}\mathcal{B}+p_{t}p_{\varphi}\mathcal{D}+p^{2}_{r}\mathcal{C}\right]= ur[𝒟prpt+𝒜prpφ]+\displaystyle u^{r}\left[\mathcal{D}p_{r}p_{t}+\mathcal{A}p_{r}p_{\varphi}\right]+
ut[𝒜ptpφ+pt2𝒟],\displaystyle u^{t}\left[\mathcal{A}p_{t}p_{\varphi}+p^{2}_{t}\mathcal{D}\right],

where

𝒜^\displaystyle\mathcal{\hat{A}} =gφφ+(Sφrpt)2Rtrrt,\displaystyle=g^{\varphi\varphi}+\left(\frac{S^{\varphi r}}{p_{t}}\right)^{2}R_{trrt}, (45)
^\displaystyle\mathcal{\hat{B}} =gtt+(Sφrpt)2Rφrrφ,\displaystyle=g^{tt}+\left(\frac{S^{\varphi r}}{p_{t}}\right)^{2}R_{\varphi rr\varphi},
𝒞^\displaystyle\mathcal{\hat{C}} =grr+(Sφrpt)2Rφttφ.\displaystyle=g^{rr}+\left(\frac{S^{\varphi r}}{p_{t}}\right)^{2}R_{\varphi tt\varphi}.
𝒟^\displaystyle\mathcal{\hat{D}} =gtφ+(Sφrpt)2Rtrφr.\displaystyle=g^{t\varphi}+\left(\frac{S^{\varphi r}}{p_{t}}\right)^{2}R_{tr\varphi r}.

Hence, after following the gauge choices and invariant relations in Ref Hojman:2012me , we can solve the above system of equations to obtain

drdt=urut\displaystyle\frac{dr}{dt}=\frac{u^{r}}{u^{t}} =𝒞prpt+𝒟pφ,\displaystyle=\frac{\mathcal{C}p_{r}}{\mathcal{B}p_{t}+\mathcal{D}p_{\varphi}}, (46)
dφdt=uφut\displaystyle\frac{d\varphi}{dt}=\frac{u^{\varphi}}{u^{t}} =𝒟pt+𝒜pφpt+𝒟pφ.\displaystyle=\frac{\mathcal{D}p_{t}+\mathcal{A}p_{\varphi}}{\mathcal{B}p_{t}+\mathcal{D}p_{\varphi}}.

Note that Eqs. (45) and (46) reduces to Eqs. (42) and (41) of Ref. Benavides-Gallego:2021lqn when Rtrφr=0R_{tr\varphi r}=0 and gtφ=0g^{t\varphi}=0, e.i. when 𝒟=0\mathcal{D}=0.

IV Dynamics of spinning test particles around a rotating wormhole

The results of Sec.II and Sec.III describe the dynamics of a spinning test particle in a static and axially symmetric space-time. In this section, we apply these results to the geometry of a rotating wormhole.

IV.1 Effective potential

We begin by considering the canonical form of Eq. (5) with the following functions

N=eΦ,Φ=b0r,K=1,ω=2Jr3,b=b02r,μ=ln(1br).\begin{array}[]{cccc}N&=e^{\Phi},&\Phi=-\frac{b_{0}}{r},&\\ \\ K&=1,&\omega=\frac{2J}{r^{3}},&\\ \\ b&=\frac{b^{2}_{0}}{r},&\mu=-\ln\left(1-\frac{b}{r}\right).&\end{array} (47)

Hence, the space-time takes the form

ds2=\displaystyle ds^{2}= (e2Φω2r2sin2θ)dt22ωr2sin2θdtdφ\displaystyle-\left(e^{2\Phi}-\omega^{2}r^{2}\sin^{2}\theta\right)dt^{2}-2\omega r^{2}\sin^{2}\theta dtd\varphi (48)
+(1b02r2)1dr2+r2(dθ2+sin2θdφ2).\displaystyle+\left(1-\frac{b^{2}_{0}}{r^{2}}\right)^{-1}dr^{2}+r^{2}(d\theta^{2}+\sin^{2}\theta d\varphi^{2}).

Here, b0b_{0} is the wormhole throat, which is interpreted as the wormhole mass. When ω=0\omega=0, note that the space-time reduces to the Morris-Thorne space-time Morris:1988cz . Using the coordinate transformation r2=l2+b02r^{2}=l^{2}+b^{2}_{0}, Eq. (48) reduces to

ds2=\displaystyle ds^{2}= (e2Φω2(l2+b02)sin2θ)dt2\displaystyle-\left(e^{2\Phi}-\omega^{2}(l^{2}+b^{2}_{0})\sin^{2}\theta\right)dt^{2} (49)
2ω(l2+b02)sin2θdtdφ\displaystyle-2\omega(l^{2}+b^{2}_{0})\sin^{2}\theta dtd\varphi
+dl2+(l2+b02)(dθ2+sin2θdφ2),\displaystyle+dl^{2}+(l^{2}+b^{2}_{0})(d\theta^{2}+\sin^{2}\theta d\varphi^{2}),

with Φ\Phi and ω\omega now functions of ll. In this new radial coordinate, the throat of the wormhole is at l=0l=0 (r0=b0r_{0}=b_{0}).

Before computing the effective potential, we want to express the results in terms of dimensionless variables, where

llb0,ssb0=Smb0,b0=Lmb0,aJb02.\begin{array}[]{cccc}l\rightarrow\frac{l}{b_{0}},&s\rightarrow\frac{s}{b_{0}}=\frac{S}{mb_{0}},&\mathcal{L}\rightarrow\frac{\mathcal{L}}{b_{0}}=\frac{L}{mb_{0}},&a\rightarrow\frac{J}{b^{2}_{0}}.\end{array} (50)

Therefore, the effective potential has the following form

𝒱eff=δ2β±(δ2β)2+(ρβσ2β).\mathcal{V}_{\rm eff}=-\frac{\delta\mathcal{L}}{2\beta}\pm\sqrt{\left(\frac{\delta\mathcal{L}}{2\beta}\right)^{2}+\left(\frac{\rho}{\beta}-\frac{\sigma\mathcal{L}^{2}}{\beta}\right)}. (51)

with 𝒱eff=Veff/m\mathcal{V}_{\rm eff}=V_{\rm eff}/m101010From now on, we use 𝒱eff\mathcal{V}_{\text{eff}} with ++ in Eq. (51).. Then, after using the definitions in Eq. (36), we obtain

β\displaystyle\beta =e2Φ6ale3Φ(l2+1)2s+l2e4Φ(9a2(l2+1)2e2Φ)(l2+1)4s2,\displaystyle=e^{-2\Phi}-\frac{6ale^{-3\Phi}}{\left(l^{2}+1\right)^{2}}s+\frac{l^{2}e^{-4\Phi}\left(9a^{2}-\left(l^{2}+1\right)^{2}e^{2\Phi}\right)}{\left(l^{2}+1\right)^{4}}s^{2}, (52)
δb0δ\displaystyle\delta\rightarrow b_{0}\delta =4ae2Φ(l2+1)3/2+2e3Φ(12a2l+(l2+1)2e2Φ((l2+1)Φl))(l2+1)7/2s\displaystyle=-\frac{4ae^{-2\Phi}}{\left(l^{2}+1\right)^{3/2}}+\frac{2e^{3\Phi}\left(12a^{2}l+\left(l^{2}+1\right)^{2}e^{2\Phi}\left(\left(l^{2}+1\right)\Phi^{\prime}-l\right)\right)}{(l^{2}+1)^{7/2}}s
+2ale4Φ(l)(18a2l(l2+1)2e2Φ(3(l2+1)Φ+l))(l2+1)11/2s2,\displaystyle+\frac{2ale^{-4\Phi(l)}\left(-18a^{2}l-\left(l^{2}+1\right)^{2}e^{2\Phi}\left(3\left(l^{2}+1\right)\Phi^{\prime}+l\right)\right)}{\left(l^{2}+1\right)^{11/2}}s^{2},
σb02σ\displaystyle\sigma\rightarrow b^{2}_{0}\sigma =4a2e2Φ(l2+1)31l2+12ae3Φ(12a2l+(l2+1)2e2Φ(2(l2+1)Φ+l))(l2+1)5s\displaystyle=\frac{4a^{2}e^{-2\Phi}}{\left(l^{2}+1\right)^{3}}-\frac{1}{l^{2}+1}-\frac{2ae^{-3\Phi}\left(12a^{2}l+\left(l^{2}+1\right)^{2}e^{2\Phi}\left(2\left(l^{2}+1\right)\Phi^{\prime}+l\right)\right)}{(l^{2}+1)^{5}}s
+a2le4Φ(36a2l+(l2+1)2e2Φ(12(l2+1)Φl))+(l2+1)6Φ2(l2+1)7s2,\displaystyle+\frac{a^{2}le^{-4\Phi}\left(36a^{2}l+\left(l^{2}+1\right)^{2}e^{2\Phi}\left(12\left(l^{2}+1\right)\Phi^{\prime}-l\right)\right)+\left(l^{2}+1\right)^{6}\Phi^{\prime 2}}{(l^{2}+1)^{7}}s^{2},
ρρm2\displaystyle\rho\rightarrow\frac{\rho}{m^{2}} =s2(18a2l2e2Φ(l2+1)42lΦl2+1)+s4(81a4l4e4Φ(l2+1)8+18a2l3e2ΦΦ(l2+1)5+l2Φ2(l2+1)2)+1.\displaystyle=s^{2}\left(-\frac{18a^{2}l^{2}e^{-2\Phi}}{\left(l^{2}+1\right)^{4}}-\frac{2l\Phi^{\prime}}{l^{2}+1}\right)+s^{4}\left(\frac{81a^{4}l^{4}e^{-4\Phi}}{\left(l^{2}+1\right)^{8}}+\frac{18a^{2}l^{3}e^{-2\Phi}\Phi^{\prime}}{\left(l^{2}+1\right)^{5}}+\frac{l^{2}\Phi^{\prime 2}}{\left(l^{2}+1\right)^{2}}\right)+1.

It is straightforward to check that the last expressions reduce to those in Ref. Benavides-Gallego:2021lqn when a=0a=0.

From Eqs. (50) and (52), it is possible to see the symmetries in 𝒱eff\mathcal{V}_{\text{eff}} depending on the signs of ss, aa and \mathcal{L}. We call 𝒱effP\mathcal{V}^{P}_{\text{eff}} the “plus” configuration111111Since 𝒱eff\mathcal{V}_{\text{eff}} depends on the wormhole’s spin aa, we also found the following relations: 𝒱effM(l,s,a,)\displaystyle\mathcal{V}^{M}_{\text{eff}}(l,-s,a,\mathcal{L}) =𝒱eff(l,s,a,)=𝒱eff(l,s,a,)\displaystyle=\mathcal{V}_{\text{eff}}(l,s,-a,-\mathcal{L})=\mathcal{V}_{\text{eff}}(-l,s,a,\mathcal{L}) =𝒱eff(l,s,a,),\displaystyle=\mathcal{V}_{\text{eff}}(-l,-s,-a,-\mathcal{L}), and 𝒱effP(l,s,a,)\displaystyle\mathcal{V}^{P}_{\text{eff}}(l,s,a,\mathcal{L}) =𝒱eff(l,s,a,)=𝒱eff(l,s,a,)\displaystyle=\mathcal{V}_{\text{eff}}(l,-s,-a,-\mathcal{L})=\mathcal{V}_{\text{eff}}(-l,-s,a,\mathcal{L}) =𝒱eff(l,s,a,).\displaystyle=\mathcal{V}_{\text{eff}}(-l,s,-a,-\mathcal{L}).

𝒱effP(l,s,a,)\displaystyle\mathcal{V}^{P}_{\text{eff}}(l,s,a,\mathcal{L}) =𝒱effP(l,s,a,),\displaystyle=\mathcal{V}^{P}_{\text{eff}}(l,-s,-a,-\mathcal{L}), (53)

and 𝒱effM\mathcal{V}^{M}_{\text{eff}} the “minus” configuration

𝒱effM(l,s,a,)\displaystyle\mathcal{V}^{M}_{\text{eff}}(l,-s,a,\mathcal{L}) =𝒱effM(l,s,a,).\displaystyle=\mathcal{V}^{M}_{\text{eff}}(l,s,-a,-\mathcal{L}). (54)

The behavior of each configuration is shown in Fig. 1, where we plot together 𝒱effP\mathcal{V}^{P}_{\text{eff}} (black) and 𝒱effM\mathcal{V}^{M}_{\text{eff}} (red), as functions of ll with |a|=0.1|a|=0.1 and ||=2.24359426)|\mathcal{L}|=2.24359426). Note that the change from “plus” to “minus” configurations also changes the location of the innermost stable circular orbit (ISCO), shown in the figure with black and red dots. According to Eqs. (53) and (54), we obtain the following relations

𝒱effP(l,+0.3,0.1,2.24359)\displaystyle\mathcal{V}^{P}_{\text{eff}}(l,+0.3,0.1,-2.24359) =𝒱effP(l,0.3,0.1,2.24359),\displaystyle=\mathcal{V}^{P}_{\text{eff}}(l,-0.3,-0.1,2.24359), (55)
𝒱effM(l,0.3,0.1,2.24359)\displaystyle\mathcal{V}^{M}_{\text{eff}}(l,-0.3,0.1,-2.24359) =𝒱effM(l,+0.3,0.1,2.24359).\displaystyle=\mathcal{V}^{M}_{\text{eff}}(l,+0.3,-0.1,2.24359).
Refer to caption
Figure 1: The “plus” (black) and “minus” (red) configurations for the effective potential. The ISCO for each configuration is shown using dots. We assume b0=M=1b_{0}=M=1.
Refer to caption
Figure 2: The effective potential in different situations. The “plus” and “minus” configurations are shown in black and red colors, respectively. The values for the ISCO in each configuration are shown using dots. We assume b0=M=1b_{0}=M=1 and =2.0\mathcal{L}=2.0, 1.897882451.89788245, 2.40282807-2.40282807 and 2.46920218-2.46920218 for the first, second, third, and fourth panels respectively.

Hence, from the physical point of view, in the first configuration (“plus”), we have a system in which the particle’s spin ss and the wormhole’s spin aa align with the symmetry axis of the space-time, while the particle’s angular momentum \mathcal{L} is antiparallel. The “minus” configuration, on the other hand, corresponds to the case where the particle’s spin ss and its angular momentum \mathcal{L} are parallel to the symmetry axis of the space-time, while antiparallel to the wormhole’s angular momentum \mathcal{L}. In other words, in the “plus” configuration, the particle’s spin and the wormhole’s spin are always antiparallel to the particle’s angular momentum \mathcal{L}, in contrast to the “minus” configuration, where the particle’s spin ss and its angular moment \mathcal{L} are always antiparallel to the wormhole’s spin aa.

In Fig.2, we show the behavior of the effective potential in different situations as a function of ll. In the first panel, we show the behavior in the simplest case, i.e. when s=a=0s=a=0. In this situation, the effective potential is symmetric, and the ISCO radius locates at the same distance from the throat for both the lower and upper universes (lISCO=±1.73205081l_{ISCO}=\pm 1.73205081). In the second panel, we show the shape of 𝒱eff\mathcal{V}_{\text{eff}} when s=0.3s=0.3 and a=0.0a=0.0. In this case, as showed by Benavides-Gallego et al. in Ref. Benavides-Gallego:2021lqn , the effective potential is not symmetric and has two possible configurations: “plus” and “minus”, depending on the sign of the particle’s spin ss. In the “plus” configuration, the ISCO situates in the upper universe. When we change to the “minus” configuration, the ISCO changes from lISCO=+2.25990751l_{ISCO}=+2.25990751 to lISCO=2.25990751l_{ISCO}=-2.25990751.

In the third panel, we consider the case in which s=0.0s=0.0 and a=0.1a=0.1. From the figure, it is possible to see that 𝒱eff\mathcal{V}_{\text{eff}} is once again symmetric, with the ISCO located at the same distance from the wormhole’s throat in both universes (lISCO=±2.9936153l_{ISCO}=\pm 2.9936153). Also, note that the effective potential has a minimum value at the throat. In this particular case, changing the wormhole spin from a=0.1a=0.1 to a=0.1a=-0.1, the effective potential changes its shape drastically, see Fig. 3. There is no ISCO in this situation.

Refer to caption
Figure 3: The “minus” configuration for the third panel in Fig. 2. We assume b0=M=1b_{0}=M=1 and =2.40282807\mathcal{L}=-2.40282807.
Refer to caption
Figure 4: The effective potential in different situations. Left panel: 𝒱eff\mathcal{V}_{eff} vs. ll for different values of aa. Center: 𝒱eff\mathcal{V}_{eff} vs. ll for different values of ss. Right panel: 𝒱eff\mathcal{V}_{eff} vs. ll for different values of \mathcal{L}. In black color we show the case s=0.3s=0.3, a=0.1a=0.1 and =2.46920218\mathcal{L}=-2.46920218. We assume b0=M=1b_{0}=M=1.

In the fourth panel, we show the behavior of 𝒱eff\mathcal{V}_{\text{eff}} in the most general case, i.e. a spinning test particle moving around a rotating wormhole. The figure shows a similar behavior as the case investigated in Ref. Benavides-Gallego:2021lqn : the effective potential also has two configurations (“plus” and “minus”), which change the position of the ISCO. Nevertheless, it is possible to see a difference in the shape of 𝒱eff\mathcal{V}_{\text{eff}}. According to the figure, the minimum value of the effective potential shifts to the lower or upper universe depending on the configuration. In contrast to the case shown in the second panel of Fig. 2, where the effective potential has a minimum value at the throat. For example, in the case shown in the fourth panel of Fig. 2, 𝒱min\mathcal{V}_{\text{min}} is shifted from the lower universe to the upper universe when we change the configuration from “plus” to “minus”. This effect is mainly due to the wormhole’s spin aa.

Refer to captionRefer to captionRefer to caption\begin{array}[]{ccc}\includegraphics[scale={0.28}]{figure6aC.png}&\includegraphics[scale={0.28}]{figure6bC.png}&\includegraphics[scale={0.28}]{figure6cC.png}\end{array}

Figure 5: The first (continuous black line) and second derivatives (dashed black line) of 𝒱eff\mathcal{V}_{\text{eff}} as a function of ll. The intersection between the plots corresponds to lISCOl_{ISCO}, see the black dot in each panel. Left panel: lISCOl_{ISCO} in the lower universe. Center: lISCOl_{ISCO} in the upper universe. Right panel: Contour plot of the first (black) and second (green) derivatives of 𝒱eff\mathcal{V}_{\text{eff}}. lISCOl_{ISCO} is shown using black dots. The shape of 𝒱eff\mathcal{V}_{\text{eff}} vs. ll is shown in the small panel for the first and second panels in the figure. We assume b0=M=1b_{0}=M=1.

Refer to captionRefer to captionRefer to captionRefer to caption\begin{array}[]{cc}\includegraphics[scale={0.55}]{figure7a.pdf}&\includegraphics[scale={0.55}]{figure7b.pdf}\\ \includegraphics[scale={0.55}]{figure7c.pdf}&\includegraphics[scale={0.55}]{figure7d.pdf}\end{array}

Figure 6: lISCOl_{ISCO}, ISCO\mathcal{L}_{ISCO}, ISCO\mathcal{E}_{ISCO} and the function \mathcal{F} vs. ss for different values of a=0.1a=0.1 (blue), a=0.2a=0.2 (red), a=0.3a=0.3 (black), a=0.6a=0.6 (green) and a=0.9a=0.9 (orange) in the upper (continuous line) and lower (dashed line) universes. We assume b0=M=1b_{0}=M=1.

Refer to captionRefer to captionRefer to captionRefer to caption\begin{array}[]{cc}\includegraphics[scale={0.55}]{figure7an.pdf}&\includegraphics[scale={0.55}]{figure7bn.pdf}\\ \includegraphics[scale={0.55}]{figure7cn.pdf}&\includegraphics[scale={0.55}]{figure7dn.pdf}\end{array}

Figure 7: lISCOl_{ISCO}, ISCO\mathcal{L}_{ISCO}, ISCO\mathcal{E}_{ISCO} and the function \mathcal{F} vs. ss for different values of a=0.1a=-0.1 (blue), a=0.2a=-0.2 (red), a=0.3a=-0.3 (black), a=0.6a=-0.6 (green) and a=0.9a=-0.9 (orange) in the upper (continuous line) and lower (dashed line) universes. We assume b0=M=1b_{0}=M=1.

In Fig.4, we plot 𝒱eff\mathcal{V}_{\text{eff}} vs. ll for different values of aa, ss, and \mathcal{L}. In the left panel, we show the behavior of the effective potential for different values of aa. The figure shows how the effective potential decreases as the wormhole’s spin aa increases. In this case, the minimum value for the effective potential decreases from 𝒱min0.5\mathcal{V}_{\text{min}}\approx 0.5 to 𝒱min1\mathcal{V}_{\text{min}}\approx-1 as the wormhole’s spin aa changes from 0.10.1 to 0.40.4. In the central panel, we show the behavior of 𝒱eff\mathcal{V}_{\text{eff}} for different values of ss. From the figure, it is possible to see that 𝒱eff\mathcal{V}_{\text{eff}} does not change significantly as ss increases from 0.10.1 to 0.40.4. The figure also shows that 𝒱eff\mathcal{V}_{\text{eff}} in the lower universe is larger than those in the upper universe for all values of ss. Furthermore, the value of the 𝒱eff\mathcal{V}_{\text{eff}} at the throat is the same for all values of ss; in contrast to the first and third panels, where its value decreases as aa increases and decreases as \mathcal{L} increases, respectively.

In the right panel of Fig. 4, we plot the effective potential as a function of ll for different values of the angular momentum \mathcal{L} while keeping constant aa and ss. From the figure, we can see how 𝒱eff\mathcal{V}_{\text{eff}} decreases as the particle’s angular momentum \mathcal{L} increases. The figure also shows that the effective potential is always larger in the lower universe than in the upper one.

IV.2 Innermost stable circular orbits

Now, we focus our attention on circular orbits of spinning test particles in the geometry of a rotating wormhole given by Eq. (49). It is well-known that circular orbits occur when the radius of an orbit of particles is constant. Therefore, the radial velocity of a spinning test particle vanishes dl/dλ=0dl/d\lambda=0. In this sense, according to Eq. (34), =𝒱eff\mathcal{E}=\mathcal{V}_{\text{eff}}, where we have defined the energy of the test particle per unit mass as =E/m\mathcal{E}=E/m. On the other hand, the radial acceleration of the particle also vanishes, i.e. d2l/dλ2=0d^{2}l/d\lambda^{2}=0, which implies that d𝒱eff/dl=0d\mathcal{V}_{\text{eff}}/dl=0. Nevertheless, these conditions do not guarantee the stability of circular orbits. To ensure stability, the second radial derivative of the effective potential must be positive, namely,

d2𝒱effdl20.\frac{d^{2}\mathcal{V}_{\text{eff}}}{dl^{2}}\geq 0. (56)
Table 1: Values for lISCOl_{\text{ISCO}}, ISCO\mathcal{L}_{\text{ISCO}} and ISCO\mathcal{E}_{\text{ISCO}} for different values of aa and ss. Here, UU and LL subscripts denote upper and lower universes, respectively.
aa ss lISCOl_{\text{ISCO}} ISCO\mathcal{L}_{\text{ISCO}} ISCO\mathcal{E}_{\text{ISCO}} lISCOl_{\text{ISCO}} ISCO\mathcal{L}_{\text{ISCO}} ISCO\mathcal{E}_{\text{ISCO}}
0.1 0.0 \mp 2.994 -2.403 0.900
±0.1\pm 0.1 \mp3.155 -2.356 0.904 ±\pm2.803 -2.440 0.895
±\pm 0.4 \mp 3.519 -2.180 0.911 ±\pm 1.901 -2.436 0.865
±\pm 0.8 \mp 3.816 -1.893 0.917 ±\pm 0.895 -1.649 0.682
±\pm 1.31 \mp 3.962 -1.472 0.919 ±\pm 0.601 -0.874 0.502
±\pm 1.6 \mp 3.925 -1.210 0.919
±\pm 2.0 \mp 3.663 -0.814 0.915
0.2 0.0 \mp 3.714 -2.635 0.916
±\pm 0.1 \mp 3.870 -2.581 0.919 ±\pm 3.534 -2.682 0.913
±\pm 0.55 \mp 4.376 -2.289 0.926 ±\pm 1.209 -2.496 0.840
±\pm 1.0 \mp 4.675 -1.946 0.930
±\pm 1.5 \mp 4.829 -1.526 0.932
±\pm 2.0 \mp 4.796 -1.073 0.932
±\pm 2.5 \mp 4.482 -0.581 0.929
0.3 0.0 \mp 4.285 -2.813 0.926
±\pm 0.1 \mp 4.439 -2.755 0.928 ±\pm 4.109 -2.865 0.923
±\pm 0.52 \mp 4.928 -2.476 0.933 ±\pm 2.639 -2.948 0.900
±\pm 1.2 \mp 5.388 -1.948 0.938
±\pm 1.8 \mp 5.555 -1.435 0.940
±\pm 2.4 \mp 5.497 -0.884 0.940
±\pm 3.0 \mp 5.076 -0.290 0.937
0.6 0.0 \mp 5.619 -3.204 0.942
±\pm 0.1 \mp 5.758 -3.147 0.943 ±\pm 5.448 -3.264 0.940
±\pm 0.61 \mp 6.377 -2.782 0.947 ±\pm 3.826 -3.456 0.925
±\pm 1.6 \mp 7.055 -1.987 0.952
±\pm 2.4 \mp 7.290 -1.289 0.953
±\pm 3.2 \mp 7.259 -0.554 0.954
±\pm 4.0 \mp 6.818 0.224 0.952
0.9 0.0 \mp 6.684 -3.494 0.951
±\pm 0.1 \mp 7.598 -2.961 0.955 ±\pm 6.515 -3.558 0.950
±\pm 0.74 \mp 6.839 -3.428 0.951 ±\pm 4.051 -3.815 0.935
±\pm 1.8 \mp 8.351 -2.098 0.959
±\pm 2.7 \mp 8.677 -1.312 0.961
±\pm 3.6 \mp 8.749 -0.492 0.961
±\pm 4.5 \mp 8.498 0.362 0.961

Hence, the marginally stable circular orbit, corresponding to the smallest allowed value for stable circular orbits, also known as the innermost stable circular orbit or ISCO, can be obtained when d2𝒱eff/dl2=0d^{2}\mathcal{V}_{\text{eff}}/dl^{2}=0.

Using the conditions =𝒱eff\mathcal{E}=\mathcal{V}_{\text{eff}} and d𝒱eff/dl=0d\mathcal{V}_{\text{eff}}/dl=0 we obtain the values of \mathcal{E} and \mathcal{L} in terms of the circular orbit radius ll and then from d2𝒱eff/dl2=0d^{2}\mathcal{V}_{\text{eff}}/dl^{2}=0 the value of the radius of the ISCO for a spinning test particle. Since the process involves the solution of a non-linear system of equations for ll and \mathcal{L}, we must solve it numerically, see table 1.

In Fig. (5), we show the first and second derivatives of 𝒱eff\mathcal{V}_{\text{eff}} as a function of ll. The intersection between the curves correspond to the innermost stable circular orbit lISCOl_{ISCO}. According to the figure, when a=0.1a=0.1 and s=1.3s=1.3, there are two possible values for lISCOl_{ISCO}. One in the lower universe lISCOL=3.961l^{L}_{\text{ISCO}}=-3.961, and one in the upper universe lISCOU=0.607l^{U}_{\text{ISCO}}=0.607121212We use the upscript LL and UU for lower and upper universes, respectively..

In Figs. 6 and 7, we plot the behavior of the radius (top panel left), angular momentum (top panel right), energy (lower panel left), and \mathcal{F} (lower panel right) at the ISCO for different values of the wormhole spin aa, co-rotating (a>0a>0) and counter-rotating (a<0a<0), respectively131313The behavior between co-rotating and counter-rotating is similar. For this reason, we focus our discussion to the former case.. In Fig. 6, we can see how lISCOl_{ISCO} increases as aa increases. For example, if we consider the upper universe, the values of the ISCO for a=0.9a=0.9 (orange) are larger than those with a=0.1a=0.1. In the lower universe, on the other hand, lISCOl_{ISCO} becomes smaller as the spin parameter aa increases. See the dashed lines.

For constant values of aa, the figure shows how lISCOUl^{U}_{ISCO} increases as ss increases, reaching a maximum value. Then, the ISCO radius decreases as ss increases; see the continuous lines in the top-left panel of Fig. 6. In the lower universe, the situation is the opposite: lISCOl_{ISCO} decreases until it reaches a minimum value, and then it increases as ss increases. For example, when a=0.2a=0.2, the maximum/minimum value for lISCOl_{ISCO} is around 55 and 5-5, respectively. While for a=0.9a=0.9, the maximum/minimum value for lISCOl_{ISCO} is around 66 and 6-6, respectively. The figure also shows that there is an interval in which there is only one value of lISCOl_{ISCO}, and one interval in which there are two values. For example, in the upper universe, when a=0.1a=0.1, the figure shows only one value of lISCOl_{ISCO} if 2.5s<1.25-2.5\geq s<-1.25. When |s|1.25|s|\leq 1.25, the figure shows two values: one in the upper universe and the other in the lower universe. Then, when s>1.25s>1.25, the figure shows only one value of lISCOl_{ISCO} located in the lower universe. One can see a similar behavior for other values of aa with different ranges for ss.

In the top-right panel of Fig. 6, we show the behavior of ISCO\mathcal{L}_{ISCO} as a function of ss for different values of aa. According to the figure, the value of ISCO\mathcal{L}_{ISCO} in the upper universe decreases as ss increases. Nevertheless, for values of a<0.3a<0.3, ISCO\mathcal{L}_{ISCO} reaches a minimum value and starts to increase again, see the blue and red (continuous) lines. On the other hand, when we consider the lower universe (dashed lines), the behavior is the contrary. For values of a0.3a\geq 0.3, ISCO\mathcal{L}_{ISCO} decreases as ss increases. While for values of a<0.3a<0.3, ISCO\mathcal{L}_{ISCO} decreases, reaches a minimum, and starts to increase again, see the dashed blue and red lines in Fig. 8.

Refer to caption
Figure 8: ISCO\mathcal{L}_{ISCO} vs. ss for a=0.1a=0.1 (blue) and 0.20.2 (red). We assume b0=M=1b_{0}=M=1.

In the bottom-left panel of Fig. 6, we plot the behavior of the energy at the ISCO for different values of the wormhole’s spin aa. From the figure, for constant values of the spin ss, we see that the energy ISCO\mathcal{E}_{ISCO} increases as aa increases. Moreover, the figure also shows that ISCO\mathcal{E}_{ISCO} is always lower than =1\mathcal{E}=1. In the upper universe, while ss increases, ISCO\mathcal{E}_{ISCO} increases. Then, it reaches a maximum value and decreases until some value of ss, the lowest limit for the particle’s spin. On the other hand, when we consider the lower universe, the behavior is the opposite (as before). ISCO\mathcal{E}_{ISCO} increases as ss increases. Then, it reaches a maximum value and decreases until some value of ss.

IV.3 Superluminal bound

As mentioned before, we can obtain the superluminal bound using Eq. (37). In this sense, after replacing dr/dtdr/dt and dφ/dtd\varphi/dt (see Eq. 46), we obtain

\displaystyle\mathcal{F} =gtt(^pt+𝒟^pφ)2+grr(𝒞^pr)2+2gtφ(𝒟^pt+𝒜^pφ)\displaystyle=g_{tt}(\hat{\mathcal{B}}p_{t}+\hat{\mathcal{D}}p_{\varphi})^{2}+g_{rr}(\hat{\mathcal{C}}p_{r})^{2}+2g_{t\varphi}(\hat{\mathcal{D}}p_{t}+\hat{\mathcal{A}}p_{\varphi}) (57)
×(^pt+𝒟^pφ)+gφφ(𝒟^pt+𝒜^pφ)20.\displaystyle\times(\hat{\mathcal{B}}p_{t}+\hat{\mathcal{D}}p_{\varphi})+g_{\varphi\varphi}(\hat{\mathcal{D}}p_{t}+\hat{\mathcal{A}}p_{\varphi})^{2}\leq 0.

Therefore, the limit value of ss for which the spinning test particle’s motion is time-like can be obtained with the condition =0\mathcal{F}=0. Then, using Eqs. (32) for prp_{r} and replacing into Eq. (57), we get

=𝒳1(ptm)2+2𝒳2(ptm)(pφm)+𝒳3(pφm)2𝒳40,\mathcal{F}=\mathcal{X}_{1}\left(\frac{p_{t}}{m}\right)^{2}+2\mathcal{X}_{2}\left(\frac{p_{t}}{m}\right)\left(\frac{p_{\varphi}}{m}\right)+\mathcal{X}_{3}\left(\frac{p_{\varphi}}{m}\right)^{2}-\mathcal{X}_{4}\leq 0, (58)

with141414Here, we keep the notation of Eq. (24), where rr is the radial coordinate. Nevertheless, in the rotating wormhole, the radial coordinate changes to ll, so rlr\rightarrow l in all the equations.

𝒳1\displaystyle\mathcal{X}_{1} =gtt^2+2gtφ^𝒟^+gtφ𝒟^2gttgrr2𝒞^2\displaystyle=g_{tt}\mathcal{\hat{B}}^{2}+2g_{t\varphi}\mathcal{\hat{B}}\mathcal{\hat{D}}+g_{t\varphi}\mathcal{\hat{D}}^{2}-g^{tt}g_{rr}^{2}\mathcal{\hat{C}}^{2} (59)
𝒳2\displaystyle\mathcal{X}_{2} =gtt^𝒟^+gtφ(𝒜^^+𝒟^2)+gφφ𝒜^𝒟^gtφgrr2𝒞^2\displaystyle=g_{tt}\mathcal{\hat{B}}\mathcal{\hat{D}}+g_{t\varphi}(\mathcal{\hat{A}}\mathcal{\hat{B}}+\mathcal{\hat{D}}^{2})+g_{\varphi\varphi}\mathcal{\hat{A}}\mathcal{\hat{D}}-g^{t\varphi}g_{rr}^{2}\mathcal{\hat{C}}^{2}
𝒳3\displaystyle\mathcal{X}_{3} =gtt𝒟^2+2gtφ𝒜^𝒟^+gφφ𝒜^2gφφgrr2𝒞^2\displaystyle=g_{tt}\mathcal{\hat{D}}^{2}+2g_{t\varphi}\mathcal{\hat{A}}\mathcal{\hat{D}}+g_{\varphi\varphi}\mathcal{\hat{A}}^{2}-g^{\varphi\varphi}g_{rr}^{2}\mathcal{\hat{C}}^{2}
𝒳4\displaystyle\mathcal{X}_{4} =grr2𝒞^2.\displaystyle=g_{rr}^{2}\mathcal{\hat{C}}^{2}.

Note that Eq. (58) reduces to Eqs. (45) of Ref. Benavides-Gallego:2021lqn when we consider the space-time of a non-rotating wormhole.

Table 2: Superluminal bound for upper and lower universes in the co-rotating case (a>0a>0). In the table s+s_{+} and ss_{-} are the limits for positive and negative values of ss, respectively.
aa Universe s+s_{+} ss_{-} +\mathcal{F}_{+} \mathcal{F}_{-}
0.1 Upper 1.315 -2.30305 4.8×106-4.8\times 10^{-6} 4.1×106-4.1\times 10^{-6}
Lower 2.30305 -1.315 4.1×106-4.1\times 10^{-6} 4.8×106-4.8\times 10^{-6}
0.2 Upper 0.558933 -2.89066 4.2×107-4.2\times 10^{-7} 6.3×106-6.3\times 10^{-6}
Lower 2.89066 -0.558933 6.3×106-6.3\times 10^{-6} 4.2×107-4.2\times 10^{-7}
0.3 Upper 0.52 -3.38559 -0.931192 4.5×106-4.5\times 10^{-6}
Lower 3.38559 -0.52 4.5×106-4.5\times 10^{-6} -0.931192
0.6 Upper 0.61 -4.60353 -0.9693 3.0×106-3.0\times 10^{-6}
Lower 4.60353 -0.61 3.0×106-3.0\times 10^{-6} -0.9693
0.9 Upper 0.74 -5.4 -0.9495 -0.7493
Lower 5.4 -0.74 -0.7493 -0.9495

In the bottom-left panel of Fig. (6), we plot the behavior of \mathcal{F} as a function of ss for different values of aa. From the figure, it is possible to identify an interval for the particle’s spin ss in both the upper and lower universes. There are always two limits, one for negative values of the spin ss_{-} and the other for positive values, s+s_{+}. For example, in the upper universe, when the value of a=0.1a=0.1, the function \mathcal{F} reaches the superluminal bond at s2.3s_{-}\approx-2.3 and s+1.3s_{+}\approx 1.3, see the continuous blue line. One can see similar behavior when a=0.2a=0.2. In that case, \mathcal{F} reaches the superluminal bound when s2.8s_{-}\approx-2.8 and s+0.5s_{+}\approx 0.5. Finally, when a=0.3a=0.3, the spinning test particle reaches the superluminal bound when s3.3s_{-}\approx-3.3. Nevertheless, when s>0.52s>0.52, the system of non-linear equations

d𝒱eff/dl=0andd2𝒱eff/dl2=0.\begin{array}[]{ccc}d\mathcal{V}_{\text{eff}}/dl=0&\text{and}&d^{2}\mathcal{V}_{\text{eff}}/dl^{2}=0.\end{array} (60)

does not have a solution. As a consequence, \mathcal{F} does not reach the superluminal bound =0\mathcal{F}=0. However, although <0\mathcal{F}<0, we can consider s0.52s\approx 0.52 as the positive limit value for a spinning test particle in the upper universe, s+s_{+}.

On the other hand, for large values of aa (e. i. 0.9), the figure shows that \mathcal{F} does not reach the superluminal limit neither for negative nor for positives values of ss. Once again, although \mathcal{F} is always negative in these cases, we can use these values to set an interval of motion for a spinning test particle. Hence, when a=0.9a=0.9 the limit values are s5.4s_{-}\approx-5.4 and s+0.74s_{+}\approx 0.74.

When we consider the lower universe (dashed lines in Fig. 6), the figure shows the same behavior as in the upper universe. Once again, \mathcal{F} gives two limit values for s<0s<0 and s>0s>0. Nevertheless, due to symmetries in 𝒱eff\mathcal{V}_{\text{eff}}, these values are the same and only change in the sign. For example, in the upper universe, we found that s+1.3s_{+}\approx 1.3 and s2.3s_{-}\approx-2.3 when a=0.1a=0.1. Therefore, due to the symmetries, in the lower universe s1.3s_{-}\approx-1.3 while s+2.3s_{+}\approx 2.3. In table 2, we show some of the values.

V Conclusions

In this work, we investigated the motion of spinning test particles around a rotating wormhole, an extension of the previous work presented by Benavides-Gallego et al. in Ref. Benavides-Gallego:2021lqn , where the authors take into account a spinning test particle moving around a non-rotating traversable wormhole, the well-known Morris-Thorne wormhole Morris:1988cz . There, the authors showed how the effective potential depends on the wormhole’s spin aa, the particle’s spin ss, the proper distance ll, the angular momentum \mathcal{L}, and the wormhole’s throat b0b_{0}. According to the authors, the effective potential 𝒱eff\mathcal{V}_{\text{eff}} shows some symmetries in its behavior, represented by two configurations: plus and minus. In the plus configuration, 𝒱effP(l,s,)=𝒱effP(l,s,)\mathcal{V}^{P}_{\text{eff}}(l,s,\mathcal{L})=\mathcal{V}^{P}_{\text{eff}}(l,-s,-\mathcal{L}). In the minus configuration, on the other hand, we have that 𝒱effM(l,s,)=𝒱effM(l,s,)\mathcal{V}^{M}_{\text{eff}}(l,-s,\mathcal{L})=\mathcal{V}^{M}_{\text{eff}}(l,s,-\mathcal{L}). Therefore, a spinning test particle moving with clockwise spin and angular momentum is equivalent to a particle moving with counterclockwise spin and angular momentum (plus configuration); while a spinning test particle moving with clockwise spin and counterclockwise angular momentum is equivalent to a particle moving with counterclockwise spin and clockwise angular momentum (minus configuration).

In the rotating case, we also found that 𝒱eff\mathcal{V}_{\text{eff}} has “plus” and “minus” configurations defined as 𝒱effP(l,s,a,)\mathcal{V}^{P}_{\text{eff}}(l,s,a,\mathcal{L}) and 𝒱effP(l,s,a,)\mathcal{V}^{P}_{\text{eff}}(l,-s,a,\mathcal{L}), respectively. Moreover, since 𝒱eff\mathcal{V}_{\text{eff}} depends on the wormhole’s spin aa, there are more relations of symmetry, see footnote 1010. In each configuration, the ISCO has the same value. However, if we change the configuration, the ISCO also changes. Hence, if lISCOl_{ISCO} is in the upper universe for the “plus” configuration, it will change to the lower universe once we chance to the “minus” configuration, see Fig. 1. Furthermore, from the dynamical point of view, it is important to point out that these symmetries allow some equivalences in the motion of spinning test particles regarding the sign of ss, \mathcal{L} and aa.

On the other hand, the shape of 𝒱eff\mathcal{V}_{\text{eff}} is symmetric only in two particular cases: s=a=0s=a=0 or s=0s=0 and a0a\neq 0. In these situations, considering the upper and lower universes, the ISCO is at the same distance from the throat. Moreover, 𝒱eff\mathcal{V}_{\text{eff}} reaches a minimum value at l=0l=0. On the other hand, when aa and ss are different from zero, 𝒱eff\mathcal{V}_{\text{eff}} is non-symmetric and looks like the effective potential of a spinning test particle in a non-rotating wormhole. Nevertheless, in contrast with the non-rotating case, where the effective potential reaches a minimum value at the throat (l=0l=0), 𝒱eff\mathcal{V}_{\text{eff}} has a minimum value shifted to the right or left depending on the configuration, see Fig. 1. It is important to point out that 𝒱min\mathcal{V}_{\text{min}} decreases in three different situations: when the wormhole’s spin, the particle’s spin or angular momentum increases while keeping constant ss and \mathcal{L}, aa and \mathcal{L}, and aa and ss, respectively. See Fig. 4.

In the non-rotating case, Ref. Benavides-Gallego:2021lqn shows that there is only one value for lISCOl_{ISCO} (in the upper or lower universes) when |s|1|s|\geq 1 and two values when 1<s<1-1<s<1 (one in the upper universe and the other in the lower universe), see Fig. 8 in Ref. Benavides-Gallego:2021lqn . We found a similar behavior for the rotating case. However, the interval for ss in which there are two values of lISCOl_{ISCO} changes as the wormhole’s spin aa increases.

If we consider the lower universe, the analysis performed in Ref. Benavides-Gallego:2021lqn shows that lISCOl_{ISCO} decreases as ss increases. Then, it reaches a minimum value at s0.8s\approx-0.8 and increases again up to the throat (when s=1s=1), see Fig. 8 in Ref. Benavides-Gallego:2021lqn . The opposite behavior occurs when we consider the upper universe. In the rotating case, on the other hand, our analysis shows that lISCOl_{ISCO} in the lower universe decreases as ss decreases. Then, it reaches a minimum value and increases again. Nevertheless, in contrast to the non-rotating case, lISCOl_{ISCO} does not arrive at the wormhole’s throat. See the upper-left panel of Fig. 6, This behavior may be a consequence of the shifting effect on 𝒱eff\mathcal{V}_{\text{eff}} at the minimum.

In Ref. Zhang:2017nhl , Y. P. Zhang et al. investigated the ISCO orbit for a classical spinning test particle in the background of a Kerr-Newman black hole. There, the authors show that the motion of the spinning test particle is related to its spin, and it will be superluminal if its spin is large. According to their analysis, the authors found that lISCOl_{ISCO} decreases as the spin of the particle increases. We also found the same behavior in the rotating wormhole when the upper universe is considered (the behavior is the opposite in the lower universe). Nevertheless, in contrast with the results of Ref. Zhang:2017nhl , where a spinning test particle can orbit in a smaller circular orbit than a non-spinning test particle, we found that |lISCO||l_{ISCO}| increases as the wormhole’s spin aa increases (while keeping constant the particles’ spin ss). In this sense, a spinning test particle can move in larger circular orbits than the non-rotating case.

In this work, we found that the energy of the innermost stable circular orbit ISCO\mathcal{E}_{ISCO} is always below unity, i.e. ISCO<1\mathcal{E}_{ISCO}<1. This result agrees with the non-rotating case, where ISCO<0.9\mathcal{E}_{ISCO}<0.9 Benavides-Gallego:2021lqn . Furthermore, it is important to remark that the energy increases as the wormhole’s spin aa increases. Hence, ISCO\mathcal{E}_{ISCO} in the non-rotating case is smaller than the energy in the rotating case, where 0.8<ISCO<0.960.8<\mathcal{E}_{ISCO}<0.96. In Ref. Benavides-Gallego:2021lqn , the authors found that ISCO\mathcal{E}_{ISCO} has the same value for both the lower and upper universes when s=0s=0. In our analysis, we found the same behavior, which is a consequence of the mirror-like symmetric of 𝒱eff\mathcal{V}_{\text{eff}}151515Recall that 𝒱eff\mathcal{V}_{\text{eff}} is symmetric when s=0s=0 in both the rotating and non-rotating wormhole, see the first and third panels in Fig 2..

In this work, we also found that the angular momentum at the ISCO behaves differently from the non-rotating case. In Ref. Benavides-Gallego:2021lqn , the authors showed that the angular momentum reaches the maximum value ISCO2.002\mathcal{L}_{ISCO}\approx 2.002 when |s|0.25|s|\approx 0.25 Benavides-Gallego:2021lqn and it is always positive, ISCO>0\mathcal{L}_{ISCO}>0. Moreover, in the lower universe, the angular momentum increases as ss increases, reaching its maximum value at s0.25s\approx-0.25. Then, it decreases for 0.25<s<2-0.25<s<2. See the left-bottom panel of Fig. 8 in Ref. Benavides-Gallego:2021lqn . The opposite behavior occurs when we consider the upper universe. In the rotating case, on the other hand, we found that 4<ISCO<1.4-4<\mathcal{L}_{ISCO}<1.4. Furthermore, for small values of ss (0.10.1, 0.20.2 or 0.30.3), we show ISCO\mathcal{L}_{ISCO} in the upper universe decreases as ss increases. Then, it reaches a minimum value and starts to increase again to a certain value of ss. The opposite behavior occurs in the lower universe. Furthermore, for larger values of ss (0.60.6 or 0.90.9), the ISCO\mathcal{L}_{ISCO} in the lower universe decreases as ss decreases, reaching a minimum value.

Finally, we consider the superluminal bound to investigate the constraints for the particle’s spin ss. In the non-rotating case, Benavides-Gallego et al. showed that the motion of a spinning test particle has physical meaning (the trajectory is time-like) if 1.5<s<1.5-1.5<s<1.5. Nevertheless, for values of |s|>1.5|s|>1.5, the particle’s trajectory is superluminal (space-like) and it does not have a physical meaning. In the rotating case, the shape of \mathcal{F} is similar to the non-rotating case; but we found some differences. For example, when one considers small values of aa (0.10.1, or 0.20.2), =0\mathcal{F}=0 for two values of the particle’s spin: ss_{-} and s+s_{+}; see the right-bottom panel of Fig. 6. Moreover, for |s|<1|s|<1 and a=0.3a=0.3, the function \mathcal{F} stops because the non-linear system in Eq. (60) does not have a solution, setting a constraint value for ss which is positive/negative for the upper/lower universe, see table 2. For larger values of aa (0.90.9), the function \mathcal{F} never reaches the superluminal bound (=0\mathcal{F}=0) because the non-linear system Eq. (60) does not have solution. In this sense, the spin is constrained to those values for which the non-linear system has a solution, see table 2.

It is important to remark that our analysis uses the MPD equations. Therefore, we considered the approximation in which the mass and size of the spinning test particle are negligible in relation to the mass of the central object and must not affect the geometry background. Nevertheless, from the astrophysical point of view, the motion of spinning test particles may still determine some features that enable us to distinguish black holes from wormholes. As we have shown in this paper, the spin does affect the motion of test particles around a rotating wormhole. On the other hand, observationally speaking, the spinning test particles may form the accretion disk of black holes. These particles could be larger objects, such as asteroids, planets/exoplanets orbiting stellar-mass objects, or rapidly rotating black holes and neutron stars orbiting supermassive candidates. Hence, the spin could be a crucial parameter to consider when describing the motion of such objects and both the electromagnetic and gravitational wave observations that would allow us to conclude if these objects are black holes or wormholes.

Acknowledgements.
The work of F. A and C. A. B. G is supported, respectively, by Grant F-FA-2021-432 of the Uzbekistan Ministry for Innovative Development and the PIFI program of the Chinese Academy of Sciences. W. H. is supported by CAS Project for Young Scientists in Basic Research YSBR-006, and acknowledges the support of The National Key R&D Program of China (Grant No. 2021YFC2203002), NSFC (National Natural Science Foundation of China) No. 12173071 and No. 11773059. J.R, F.A. and A.A. acknowledges the support of the grants F-FA-2021-432, F-FA-2021-510, and MRB-2021-527 of the Uzbekistan Ministry for Innovative Development. A.A. also acknowledges the support of the PIFI program and J.R. thanks to the ERASMUS+ project 608715-EPP-1-2019-1-UZ-EPPKA2-JP (SPACECOM).

Appendix A The dragging of the inertial frame

In this appendix, we follow Chandrasekhar to show the existence of a dragging effect in the line element given by Eq. (5). Hence, the contravariant components of the metric tensor gμνg^{\mu\nu} can be expressed in the matrix form as161616Recall that we set t0t\rightarrow 0, θ2\theta\rightarrow 2 and r3r\rightarrow 3.

(gμν)=(1N2ωN200ωN2AK2N2r200001K2r20000eμ),(g^{\mu\nu})=\left(\begin{array}[]{cccc}-\frac{1}{N^{2}}&-\frac{\omega}{N^{2}}&0&0\\ -\frac{\omega}{N^{2}}&-\frac{A}{K^{2}N^{2}r^{2}}&0&0\\ 0&0&\frac{1}{K^{2}r^{2}}&0\\ 0&0&0&e^{-\mu}\\ \end{array}\right), (61)

where is defined as

A=(KrωNcsc(θ))(Krω+Ncsc(θ)).A=(Kr\omega-N\csc(\theta))(Kr\omega+N\csc(\theta)). (62)

This space-time has associated the following tetrad171717We use the same notation as in Ref. Chandrasekhar:1985kt

e(0)μ\displaystyle e_{(0)\mu} =(N,0,0,0),\displaystyle=(-N,0,0,0), (63)
e(1)μ\displaystyle e_{(1)\mu} =(rKωsinθ,rKsinθ,0,0),\displaystyle=(-rK\omega\sin\theta,rK\sin\theta,0,0),
e(2)μ\displaystyle e_{(2)\mu} =(0,0,rK,0),\displaystyle=(0,0,rK,0),
e(3)μ\displaystyle e_{(3)\mu} =(0,0,0,eμ2).\displaystyle=(0,0,0,e^{\frac{\mu}{2}}).

using the relation e(a)μ=gμνe(a)νe^{\;\;\;\;\mu}_{(a)}=g^{\mu\nu}e_{(a)\nu}, the contravariant vectors are given by

e(0)μ\displaystyle e^{\;\;\;\;\mu}_{(0)} =(1N,ωN,0,0)\displaystyle=\left(\frac{1}{N},\frac{\omega}{N},0,0\right) (64)
e(1)μ\displaystyle e^{\;\;\;\;\mu}_{(1)} =(0,1rKsinθ,0,0),\displaystyle=\left(0,\frac{1}{rK\sin\theta},0,0\right),
e(2)μ\displaystyle e^{\;\;\;\;\mu}_{(2)} =(0,0,1rK,0),\displaystyle=\left(0,0,\frac{1}{rK},0\right),
e(3)μ\displaystyle e^{\;\;\;\;\mu}_{(3)} =(0,0,0,eμ2).\displaystyle=\left(0,0,0,e^{-\frac{\mu}{2}}\right).

Therefore, for the tetrad so defined, we have

e(a)μe(b)μ=η(a)(b)=(1000010000100001).e^{\;\;\;\;\mu}_{(a)}e_{(b)\mu}=\eta_{(a)(b)}=\left(\begin{array}[]{cccc}-1&0&0&0\\ 0&1&0&0\\ 0&0&1&0\\ 0&0&0&1\end{array}\right). (65)

This means that the chosen frame is Minkowskian, i.e. it represents locally an inertial frame. It is straightforward to show that the components of the line element (5) can be obtained using the relation

gμν=η(a)(b)e(a)μe(b)ν.g_{\mu\nu}=\eta^{(a)(b)}e_{(a)\mu}e_{(b)\nu}. (66)

The components of the four-velocity are given by

u0=dtdλ,u1=Ωu0,uα=vαu0,\begin{array}[]{ccc}u^{0}=\frac{dt}{d\lambda},&u^{1}=\Omega u^{0},&u^{\alpha}=v^{\alpha}u^{0}\ ,\end{array} (67)

where α=2\alpha=2, 33, vα=dxα/dtv^{\alpha}=dx^{\alpha}/dt, and Ωdφ/dt\Omega\equiv d\varphi/dt. In the inertial frame, the components of the four-velocity are computed using the following relation

u(a)=η(a)(b)e(b)μuμ,u^{(a)}=\eta^{(a)(b)}e_{(b)\mu}u^{\mu}, (68)

from which

u(0)\displaystyle u^{(0)} =Nu0,\displaystyle=Nu^{0}, (69)
u(1)\displaystyle u^{(1)} =(Ωω)rKsinθu0,\displaystyle=(\Omega-\omega)rK\sin\theta u^{0},
u(2)\displaystyle u^{(2)} =rKv2u0,\displaystyle=rKv^{2}u^{0},
u(3)\displaystyle u^{(3)} =eμ2v3u0.\displaystyle=e^{\frac{\mu}{2}}v^{3}u^{0}.

From the second relation in Eq. (69), we can conclude that a point moving in circular motion with angular velocity Ω\Omega in the coordinate system (t,φ,θ,r)(t,\varphi,\theta,r), will move with an angular velocity (Ωω)rKsinθu0(\Omega-\omega)rK\sin\theta u^{0} in the inertial frame. Similarly, a point which is considered at rest in the local inertial frame (i.e. u(1)=u(2)=u(3)=0u^{(1)}=u^{(2)}=u^{(3)}=0), will have an angular velocity ω\omega in the coordinate frame. Therefore, the non-vanishing of ω\omega is said to describe a dragging of the inertial frame. Since the space-time is asymptotically flat, then ω=2J/r3\omega=2J/r^{3}.

Appendix B Calculation of Dpα/dλDp^{\alpha}/d\lambda for tt, rr and φ\varphi

From the first MPD equation in Eq. (17), we obtain

Dpνdλ=12RνβδσuβSδσ.\frac{Dp_{\nu}}{d\lambda}=-\frac{1}{2}R_{\nu\beta\delta\sigma}u^{\beta}S^{\delta\sigma}. (70)

from which one can compute the Dpt/dλDp_{t}/d\lambda, Dpr/dλDp_{r}/d\lambda, and Dpφ/dλDp_{\varphi}/d\lambda in terms of the components of the Riemann tensor.

In the case of Dpt/dλDp_{t}/d\lambda, one obtains the following expression,

Dptdλ=12[2RttδσutSδσ+2RtrδσurSδσ+2RtφδσuφSδσ].\frac{Dp_{t}}{d\lambda}=-\frac{1}{2}\left[2R_{tt\delta\sigma}u^{t}S^{\delta\sigma}+2R_{tr\delta\sigma}u^{r}S^{\delta\sigma}+2R_{t\varphi\delta\sigma}u^{\varphi}S^{\delta\sigma}\right]. (71)

The factor 22 comes because RνβδσR_{\nu\beta\delta\sigma} and SδσS^{\delta\sigma} are skew-symmetric tensors. Therefore, we have to count twice in the sum because RabδσSδσ=RabσδSσδR_{ab\delta\sigma}S^{\delta\sigma}=R_{ab\sigma\delta}S^{\sigma\delta}. Hence, from Eq. (71), we obtain

Dptdλ=\displaystyle\frac{Dp_{t}}{d\lambda}= (RtttrStr+RtttφStφ+RttrφSrφ)ut\displaystyle-\left(R_{tttr}S^{tr}+R_{ttt\varphi}S^{t\varphi}+R_{ttr\varphi}S^{r\varphi}\right)u^{t} (72)
(RtrtrStr+RtrtφStφ+RtrrφSrφ)ur\displaystyle-\left(R_{trtr}S^{tr}+R_{trt\varphi}S^{t\varphi}+R_{trr\varphi}S^{r\varphi}\right)u^{r}
(RtφtrStr+RtφtφStφ+RtφrφSrφ)uφ.\displaystyle-\left(R_{t\varphi tr}S^{tr}+R_{t\varphi t\varphi}S^{t\varphi}+R_{t\varphi r\varphi}S^{r\varphi}\right)u^{\varphi}.

Then, after using Eq. (27) and considering the non-vanishing components of the Riemann tensor, the last expression reduces to

Dptdλ=Sφrpt[(pφRtrtrptRtrφr)urprRtφtφuφ].\frac{Dp_{t}}{d\lambda}=\frac{S^{\varphi r}}{p_{t}}\left[\left(p_{\varphi}R_{trtr}-p_{t}R_{tr\varphi r}\right)u^{r}-p_{r}R_{t\varphi t\varphi}u^{\varphi}\right]. (73)

In a similar way, we obtain

Dprdλ\displaystyle\frac{Dp_{r}}{d\lambda} =[RrtδσutSδσ+RrrδσurSδσ+RrφδσuφSδσ],\displaystyle=-\left[R_{rt\delta\sigma}u^{t}S^{\delta\sigma}+R_{rr\delta\sigma}u^{r}S^{\delta\sigma}+R_{r\varphi\delta\sigma}u^{\varphi}S^{\delta\sigma}\right], (74)
Dpφdλ\displaystyle\frac{Dp_{\varphi}}{d\lambda} =[RφtδσutSδσ+RφrδσurSδσ+RφφδσuφSδσ].\displaystyle=-\left[R_{\varphi t\delta\sigma}u^{t}S^{\delta\sigma}+R_{\varphi r\delta\sigma}u^{r}S^{\delta\sigma}+R_{\varphi\varphi\delta\sigma}u^{\varphi}S^{\delta\sigma}\right].

From which, after using Eq. (27) and considering the non-vanishing components of the Riemann tensor, we obtain

Dprdλ\displaystyle\frac{Dp_{r}}{d\lambda} =Sφrpt[(pφRrttr+ptRrtrφ)ut+(pφRrφtr+ptRrφrφ)uφ],\displaystyle=\frac{S^{\varphi r}}{p_{t}}\left[(p_{\varphi}R_{rttr}+p_{t}R_{rtr\varphi})u^{t}+(p_{\varphi}R_{r\varphi tr}+p_{t}R_{r\varphi r\varphi})u^{\varphi}\right], (75)
Dpφdλ\displaystyle\frac{Dp_{\varphi}}{d\lambda} =Sφrpt[ur(pφRφrtr+ptRφrrφ)prRφttφut].\displaystyle=\frac{S^{\varphi r}}{p_{t}}\left[u^{r}(p_{\varphi}R_{\varphi rtr}+p_{t}R_{\varphi rr\varphi})-p_{r}R_{\varphi tt\varphi}u^{t}\right].

Note that in the case of a non-rotating wormhole, Rtrφr=0R_{tr\varphi r}=0. Therefore, Eqs. (73), (75) reduces to the Eqs. (A11) and (A12) of Ref. Benavides-Gallego:2021lqn , respectively.

References

  • (1) A. Einstein and N. Rosen, Phys. Rev. 48, 73-77 (1935)
  • (2) G. W. Gibbons, Gen. Rel. Grav. 47 (2015) 71.
  • (3) K. Schwarzschild, Sitzungsber. Preuss. Akad. Wiss. Berlin (Math. Phys. ) 1916, 189-196 (1916) [arXiv:physics/9905030 [physics]].
  • (4) K. Schwarzschild, Sitzungsber. Preuss. Akad. Wiss. Berlin (Math. Phys. ) 1916, 424-434 (1916) [arXiv:physics/9912033 [physics.hist-ph]].
  • (5) L. Flamm, Gen. Rel. Grav. 47 (1916) 72.
  • (6) R. W. Fuller and J. A. Wheeler, Phys. Rev. 128, 919-929 (1962) doi:10.1103/PhysRev.128.919
  • (7) H. G. Ellis, J. Math. Phys. 14, 104-118 (1973) doi:10.1063/1.1666161
  • (8) K. A. Bronnikov, Acta Phys. Polon. B 4, 251-266 (1973)
  • (9) M. S. Morris and K. S. Thorne, Am. J. Phys. 56, 395-412 (1988) doi:10.1119/1.15620
  • (10) M. S. Morris, K. S. Thorne and U. Yurtsever, Phys. Rev. Lett. 61, 1446-1449 (1988) doi:10.1103/PhysRevLett.61.1446
  • (11) M. Visser, LA-UR-89-1008.
  • (12) M. Visser, Phys. Rev. D 39, 3182-3184 (1989) doi:10.1103/PhysRevD.39.3182 [arXiv:0809.0907 [gr-qc]].
  • (13) M. Visser, Nucl. Phys. B 328, 203-212 (1989) doi:10.1016/0550-3213(89)90100-4 [arXiv:0809.0927 [gr-qc]].
  • (14) M. Azreg-Ainou and G. Clement, NTH-89-9.
  • (15) M. Visser, Phys. Lett. B 242, 24-28 (1990) doi:10.1016/0370-2693(90)91588-3
  • (16) E. Poisson and W. Israel, Phys. Rev. Lett. 63, 1663-1666 (1989) doi:10.1103/PhysRevLett.63.1663
  • (17) M. Visser, PRINT-90-0513 (WASH.U.,ST.LOUIS).
  • (18) V. P. Frolov and I. D. Novikov, Phys. Rev. D 42, 1057-1065 (1990) doi:10.1103/PhysRevD.42.1057
  • (19) M. Visser, Phys. Rev. D 43, 402-409 (1991) doi:10.1103/PhysRevD.43.402
  • (20) P. Martín-Moruno and M. Visser, Phys. Rev. D 88, no.6, 061701 (2013) doi:10.1103/PhysRevD.88.061701 [arXiv:1305.1993 [gr-qc]].
  • (21) F. Echeverria, G. Klinkhammer and K. S. Thorne, Phys. Rev. D 44 (1991), 1077-1099 doi:10.1103/PhysRevD.44.1077
  • (22) S. Deser and R. Jackiw, Comments Nucl. Part. Phys. 20 (1992) no.6, 337-354 [arXiv:hep-th/9206094 [hep-th]].
  • (23) S. Deser, Class. Quant. Grav. 10 (1993), S67-S73 doi:10.1088/0264-9381/10/S/006
  • (24) D. Hochberg and M. Visser, Phys. Rev. D 56, 4745-4755 (1997) doi:10.1103/PhysRevD.56.4745 [arXiv:gr-qc/9704082 [gr-qc]].
  • (25) A. Abdujabbarov, B. Juraev, B. Ahmedov and Z. Stuchlík, Astrophys. Space Sci. 361, no.7, 226 (2016) doi:10.1007/s10509-016-2818-9
  • (26) A. A. Abdujabbarov and B. J. Ahmedov, Astrophys. Space Sci. 321, 225-232 (2009) doi:10.1007/s10509-009-0023-9 [arXiv:0903.0446 [gr-qc]].
  • (27) D. Pugliese, H. Quevedo and R. Ruffini, Phys. Rev. D 84, 044030 (2011) doi:10.1103/PhysRevD.84.044030 [arXiv:1105.2959 [gr-qc]].
  • (28) F. Willenborg, S. Grunau, B. Kleihaus and J. Kunz, Phys. Rev. D 97, no.12, 124002 (2018) doi:10.1103/PhysRevD.97.124002 [arXiv:1801.09769 [gr-qc]].
  • (29) M. Argañaraz and O. Lasso Andino, Class. Quant. Grav. 38, no.4, 045004 (2021) doi:10.1088/1361-6382/abcf86 [arXiv:1906.11779 [gr-qc]].
  • (30) Z. Stuchlík and J. Vrba, Universe 7, no.8, 279 (2021) doi:10.3390/universe7080279 [arXiv:2108.09562 [gr-qc]].
  • (31) C. A. Benavides-Gallego, W. B. Han, D. Malafarina, B. Ahmedov and A. Abdujabbarov, Phys. Rev. D 104, no.8, 084024 (2021) doi:10.1103/PhysRevD.104.084024 [arXiv:2107.07998 [gr-qc]].
  • (32) J. W. Moffat and T. Svoboda, Phys. Rev. D 44 (1991), 429-432 doi:10.1103/PhysRevD.44.429
  • (33) A. Carlini, “Spacetime Wormholes, Stringy Black Holes and 5-th Time Gravity,”
  • (34) B. Bhawal and S. Kar, Phys. Rev. D 46 (1992), 2464-2468 doi:10.1103/PhysRevD.46.2464
  • (35) P. S. Letelier and A. Wang, Phys. Rev. D 48 (1993), 631-646 doi:10.1103/PhysRevD.48.631
  • (36) D. N. Vollick, Class. Quant. Grav. 16 (1999), 1599-1604 doi:10.1088/0264-9381/16/5/309 [arXiv:gr-qc/9806096 [gr-qc]].
  • (37) E. Teo, Phys. Rev. D 58, 024014 (1998)
  • (38) R. M. Wald, “General Relativity”, (University of Chicago Press, Chicago, 1984)
  • (39) J. B. Hartle, Astrophys. J. 150, 1005-1029 (1967)
  • (40) J. B. Hartle and K. S. Thorne, Astrophys. J. 153, 807 (1968)
  • (41) K. S. Thorne, “Relativistic stars, black holes and gravitational waves (including an in-depth review of the theory of rotating, relativistic stars)”, in General Relativity and Cosmology, ed. R.K. Sachs (Academic Press, New York, 1971), pp. 237-283.
  • (42) A. Papapetrou, Ann. Inst. H. Poincare Phys. Theor. 4, 83-105 (1966)
  • (43) B. Carter, J. Math. Phys. 10, 70-81 (1969).
  • (44) S. Chandrasekhar, “The mathematical theory of black holes”, (Clarendon Press, Oxford, 1992)
  • (45) T. Harko, Z. Kovacs and F. S. N. Lobo, Phys. Rev. D 79, 064001 (2009)
  • (46) R. M. Wald, “General Relativity,” doi:10.7208/chicago/9780226870373.001.0001
  • (47) M. Mathisson, Acta Phys. Polon. 6, 163-2900 (1937)
  • (48) A. Papapetrou, Proc. Roy. Soc. Lond. A 209, 248-258 (1951) doi:10.1098/rspa.1951.0200
  • (49) E. Corinaldesi and A. Papapetrou, Proc. Roy. Soc. Lond. A 209, 259-268 (1951) doi:10.1098/rspa.1951.0201
  • (50) W. Tulczyjew, Acta Phys. Pol18, 393 (1959).
  • (51) B. Tulzcyjew and W. Tulzcyjew, “Recent Developments in General Relativity.” (1962)
  • (52) Møller, Christian, “On the definition of the centre of gravity of an arbitrary closed system in the theory of relativity,” (1949)
  • (53) W. Beiglböck, Communications in Mathematical Physics 5, 106-130 (1967)
  • (54) W. G. Dixon, Il Nuovo Cimento 34, 317-339 (1964)
  • (55) W. G. Dixon, Proc. Roy. Soc. Lond. A 314, 499-527 (1970) doi:10.1098/rspa.1970.0020
  • (56) W. G. Dixon, Proc. Roy. Soc. Lond. A 319, 509-547 (1970) doi:10.1098/rspa.1970.0191
  • (57) J. Ehlers , E. Rudolph, Gen. Relat. Gravit. 8, 197-217 (1977) https://doi.org/10.1007/BF00763547
  • (58) A. A. Deriglazov and W. Guzmán Ramírez, Adv. Math. Phys. 2017, 7397159 (2017)
  • (59) A. A. Deriglazov and W. Guzmán Ramírez, Phys. Lett. B 779, 210-213 (2018)
  • (60) M. Saijo, K. i. Maeda, M. Shibata and Y. Mino, Phys. Rev. D 58, 064005 (1998) doi:10.1103/PhysRevD.58.064005
  • (61) B. Toshmatov and D. Malafarina, Phys. Rev. D 100, no.10, 104052 (2019)
  • (62) C. Conde, C. Galvis and E. Larrañaga, Phys. Rev. D 99 (2019) no.10, 104059 doi:10.1103/PhysRevD.99.104059 [arXiv:1905.01323 [gr-qc]].
  • (63) S. A. Hojman and F. A. Asenjo, Class. Quant. Grav. 30, 025008 (2013) doi:10.1088/0264-9381/30/2/025008 [arXiv:1203.5008 [physics.gen-ph]].
  • (64) Y. P. Zhang, S. W. Wei, W. D. Guo, T. T. Sui and Y. X. Liu, Phys. Rev. D 97, no.8, 084056 (2018) doi:10.1103/PhysRevD.97.084056 [arXiv:1711.09361 [gr-qc]].