Distinguishing rotating naked singularities from Kerr-like wormholes by their deflection angles of massive particles

We study the gravitational deflection of relativistic massive particles by Janis–Newman–Winicour (JNW) spacetimes (also known as a rotating source with a surface-like naked singularity), and a rotating Kerr-like wormholes. Based on the recent article (Jusufi in Phys Rev D 98:064017, 2018), we extend some of these results by exploring the effects of naked singularity and Kerr-like objects on the deflection of particles. We start by introducing coordinate transformation leading to an isotropic line element which gives the refraction index of light for the corresponding optical medias. On the other hand, the refraction index for massive particles is found by considering those particles as a de Broglie wave packets. To this end, we apply the Gauss–Bonnet theorem to the isotropic optical metrics to find the deflection angles. Our analysis shows that, in the case of the JNW spacetime the deflection angle is affected by the parameter 0<γ<1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0<\gamma <1$$\end{document}, similarly, we find that the deformation parameter λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lambda $$\end{document} affects the deflection angle in the case of Kerr-like wormholes. In addition to that, we presented a detailed analysis of the deflection angle by means of the Hamilton–Jacobi equation that lead to the same results. As a special case of our results the deflection angle of light is recovered. Finally, we point out that the deflection of particles by Kerr-like wormholes is stronger compared to JNW spacetime, in particular this difference can be used to shed some light from observational point of view in order to distinguish the two spacetimes.


Introduction
Deflection of light ray in weak gravitational field was the first experimental proof of the theory of general relativity for the Schwarzschild spacetime since 1919 [1]. Now it is serving as one of the important tool of probing a number of interesting phenomena in astronomy and cosmology. It is widely used for investigating the astrophysical objects, e.g., black holes, wormholes, super-dense neutron stars, naked singularity etc. This phenomenon has also been considered to study exoplanet detection from the dark energy, dark matter, and Hubble parameter measurements. Moreover, the gravitational lensing (GL) provides an alternative approach towards mass measurement that leads to measure the mass of the white dwarf Stein 2051 B by using the Hubble space telescope [2]. Earlier works focused on the weak field theory of gravitational lensing which is based on the first order expansion of the small deflection angle [3]. Although the main disadvan-tage is when the lens is a very compact object (e.g. black hole), in this case the weak field approximation is no longer valid. The weak field approximation can not able to distinguish between various different solutions that are asymptotically flat. Therefore, after the pioneering work of Darwin [4], the strong gravitational lensing has received much attentions. Note that in the case of strong lensing, the effect of lensing is strong, which means multiple imaging or strong deformation of the image. On the other hand, in case of weak lensing, the effect is weak, and the lensing images are characterized by weak distortions and small magnifications. In this case the fact of lensing cannot be identified by using an individual source, but can be revealed by statistical analysis of many objects. In addition, we can distinguish situations of lensing on the basis of the magnitude of the deflection angle: weak deflection of photons (with small deflection angles, sometimes referred as weak field situation) and strong deflection of photons (deflection angles are not small, sometimes referred as strong field situation), correspondingly. However, the crucial thing is that all observational areas of GL (strong, weak and micro) are based on using only weak deflection approximation meaning that, the deflection angle remains to be small. Strong deflection case on the other hand, is used for another problems: for example, for investigations of socalled relativistic images. To this end, we wish to point out the strong gravitational lensing for the Schwarzschild black hole [5,6], then extended to the Reissner-Nordström black hole [7,8], regular black holes [9], and wormholes [10][11][12][13], while the investigations of so-called relativistic images see also Refs. [14][15][16][17][18][19][20].
An interesting method has immersed for computing the angle of light deflection by Gibbons and Werner [21,22]. They considered the source and the receiver are located at an asymptotic region and then they applied the Gauss-Bonnet theorem (GBT) to the optical metric of a lens which allows us to describes a light ray as a spatial curve. Later on Werner [23] extended this approach to the stationary spacetime and computed the deflection angle for the Kerr black hole. This method has also been applied in the weak limit approximation to compute the deflection of light for global monopole and rotating cosmic string spacetime [24,25]. Note that the deflection angle of conical singularities has been pointed out in the old paper [26] . Nowadays the scientific community is discussing this phenomenon with great interest. Along this line of thinking, Ishihara et al. [27] obtained the bending angle of light by using the GBT for finite distance from a lens object to a light source and a receiver. Further, this approach was extended by considering the stationary, axisymmetric and asymptotically flat spacetimes [28]. Jusufi and Övgün [29] studied the effect of the cosmological constant on the deflection angle by a rotating cosmic string. For the influence of cosmological constant on the light deflection see the recent review by Lebedev and Lake [30]. On the other hand the GBT method has been considered for the several wormholes [31][32][33] and black holes [34][35][36] solutions.
It is worth noticing that in the gravitational field the bending trajectories not only found for light or massless particles but also for massive particles. For a review on massive particles whose trajectories can go through a bending process, we refer the reader to [37][38][39], in particular the deflection angle for massive particles in weak deflection case is calculated in Ref. [40], and for strong deflection case see [41,42]. Thus the important fact is trajectory bending due to the strong gravity near the black holes and the large velocity of the hypervelocity stars needs a full relativistic treatment, which will differ from the calculations made using only Newtonian gravity. In recent years, the existence of strong gravitationally lensed supernovae has long been predicted [43][44][45][46]. Thus it is natural to expect that the neutrinos and cosmic rays emitted by these supernova and AGN are also lensed. Neutrino gravitational lensing by astrophysical objects has already been done in [47]. Besides their positions, they also estimated the depletion of the neutrino flux after crossing a massive object [47]. At the same time, Eiroa and Romero [48] have studied gravitational lensing of neutrinos by Schwarzschild black holes.
Recently, Crisnejo and Gallo [49] showed that how one can successfully study the GBT theorem in a plasma medium for a static and spherically symmetric gravitational field. The gravitational deflection in plasma medium was considered in details before, see [50,51]. Even in the GBT method, they computed the deflection of massive particles with the same symmetries. On the other hand, Jusufi [52] introduced a new method based on GBT to compute the deflection angle for massive particles in a rotating spacetime geometry such as the Kerr black hole [53] and Teo wormhole [54]. In doing so, he used an isotropic type metrics for a linearized rotating gravitational field where the refractive index of the corresponding optical media can be obtained. Motivated by the works we study the gravitational deflection of massive particles by Janis-Newman-Winicour (JNW) spacetimes [55] (also known as a rotating source with a surface-like naked singularity), and a rotating Kerr-like wormholes [56]. Along this line of thought, recently different aspects of naked singularities and Kerr-like wormholes have been studied, including the iron line spectroscopy [57], accretion disc properties [58], spin precession of a test gyroscope [59,60], echoes of Kerrlike wormholes [56], shadow images [61] the strong and the weak deflection of light by Kerr-like wormholes [62,63]. In the present article, we study the deflection of massive particles by assuming the propagating particles as a de Broglie wave packets introduced in [64].
The paper is organized as follows. In Sect. 2, we use a coordinate transformation to recover an isotropic metric form for the rotating naked singularities, then derive the refractive index of the corresponding optical media. In Sect. 3, we consider the problem of the gravitational deflection of massive particles by Kerr-like wormholes. We present a detailed analyses of the deflection of massive particles by using the Hamilton-Jacobi approach in Sect. 4. We conclude our results in Sect. 5. We use the geometrized unit G = c =h = 1 throughout this paper.

Deflection of massive particles by rotating naked singularities
In this section, we calculate the deflection angle of rotating naked singularities thereby we consider Janis-Newman-Winicour (JNW) spacetime. It is a solution of the Einstein equation described with the assumption of spherical symmetry, asymptotic flatness with the static massless scalar field . The action associated with this spacetime [55,57] is the variation of action turns out the following field equations where the scalar field is given by The corresponding rotating spacetime metric [57] has the form where m = M/γ , M is the Arnowitt-Deser-Misner (ADM) mass, and = r 2 + a 2 cos 2 θ, ω = a sin 2 θ, = r 2 + a 2 − 2mr, Note that a = J/M, where J is the spin angular momentum of the source. For an instant if we set γ = 1 and q = 0, one can recover the Kerr black hole solution and the larger root of = 0 provides the radius of the event horizon, It can be easily checked that for 0 < γ < 1, and q = 0, the metric (4) describes rotating naked singularities with a scalar curvature [57] R = 2(γ 2 − 1)m 2 5 1 − 2mr γ −3 × (r 2 − a 2 cos 2 θ) 2 + (ra 2 sin 2 (2θ)) 2 , and diverges on the surface r is the surface-like singularity, where curvature invariants diverge and the spacetime is geodetically incomplete as well.
In addition, the scalar curvature diverges where g tt = 0 and represent the presence of a curvature singularity at as the domain of variation of θ depends on the ratio γ 2 a 2 /m 2 . Also, it follows that in the transition from a black hole (γ = 1) to a naked singularity (γ < 1), the singularity is indeed naked, not the event horizon.

Deflection angle
Since, we are mainly interested in the weak limit, we shall simplify the problem by considering the deflection of particles in the equatorial plane (θ = π/2) with a linearized metric in a, which results in Now introducing a coordinate transformation via we deduce an isotropic metric with the line element given by where The velocity of light in this particular optical media can be found by setting ds 2 = 0, and using the definition v = |dρ|/dt, where |dρ| 2 = dρ 2 + ρ 2 dϕ 2 . With this information in hand, we can find the effective refractive index of the optical media by using n(ρ) = c/v(ρ) and c = 1, yielding We can, therefore, write down the optical metric in terms of the refractive index as follows By inverting the coordinate transformation (10), we obtain which shows leading terms in ρ r − m. The above relations holds for the deflection of light. However, we are interested in studying the deflection of massive particles. One way to incorporate this is to consider the propagating massive particles as the de Broglie wave packets. In particular, following [52,64], we can argue that by using p =hk and H = E =hω, (note that we have temporarily introducedh) results where μ is the rest mass of the particle. We rearrange this equation to obtain a constant quantity of the given optical media In other words, this equation is a generalization of wellknown result in wave-optics. For massive particles, in a given isotropic metric, the expression should be constant everywhere in the optical medium, i.e., λN = const. Thus the refractive index relation for massive particles is given by modifying the optical metric Furthermore, the energy of the particles measured at infinity is given by where w represents the relativistic velocity. The angular momentum can be written as where b is the impact parameter. Following the definition of the impact parameter, we can write which reduces to b in the case of light w = c = 1. Therefore, the quantity dϕ/dt, in case of massive particles is modified as follow yielding This result clearly indicates that the refractive index is modified due to the angular momentum parameter a.
Let us continue to calculate the Gaussian optical curvature using metric (19). Rewriting the optical metric (19), in terms of a new coordinates, thus the optical metric reads In terms of these coordinates, the Gaussian optical curvature K can be expressed as: We can, therefore, express the Gaussian optical curvature for massive particles in terms of the modified refractive index as In particular, by using (15) we can rewrite the last equation, in terms of the old coordinate r , which in leading order terms results with With these results in hand, we can chose a non-singular region D R with boundary ∂D R = γg ∪ C R . The GBT provides a connection between geometry (in a sense of optical curvature) and topology (in a sense of Euler characteristic number) stated as follow with κ is the geodesic curvature, and K is the Gaussian optical curvature. Now we can choose a non-singular domain with Euler characteristic number χ(D R ) = 1. The geodesic curvature κ, in the case of massive particles is defined by the relation but also keeping in mind an additional unit speed conditioñ g(γ ,γ ) = 1, withγ being the unit acceleration vector. The two corresponding jump angles in the limit R → ∞, reads θ O +θ S → π . Therefore the GBT now is simplified as follow By construction, there is a zero contribution from the geodesics i.e. κ(γg) = 0, therefore, we end up with a contribution only for the the curve C R , thus we can write To evaluate the above expression one can express the curve in terms of the radial coordinate C R := r (ϕ) = R = const, resulting with a nonzero contribution for the radial component The first term vanishes, however, using the unit speed condition and the optical metric, we find the final result We see that the final result of the geodesic curvature is slightly modified, letting w = 1, we recover κ(C R ) → R −1 , as expected. For an observer located at the coordinate R, we find lim Putting together these results, from the GBT, we find the following expression for the deflection anglê where the surface element is given by d S = g dr dϕ N 2 (r )rdrdϕ. In addition to that, from Eq. (24) the following approximation holds N (r ) 2 w 2 . This integral can easily be evaluated, where the plus and minus sign refers to the retrograde and prograde light rays, respectively. That being said, it is clear that the deflection angle is bigger in the case of retrograde orbit, and smaller in the case of prograde orbit. Now we can consider three special cases in the last equation: (i) If γ = 1, we find the deflection angle of massive particles in Kerr black hole geometry. (ii) If w = 1, then the deflection angle of light by a JNW spacetime is recovered. (iii) Finally, letting γ = w = 1, we recover the deflection of light in a Kerr geometry. Our result for the deflection angle is consistent with previous studies [65][66][67][68][69][70][71][72][73]. It is important to realize that our method applies only for relativistic particles, i.e. 0 < w ≤ 1, in other words, there is an apparent singularity in the deflection angle in the limit w → 0. Of course, this apparent singularity can be removed but a more general setup is needed (Fig. 1).

Deflection of massive particles by Kerr-like wormholes
We consider a Kerr-like wormhole spacetime, recently proposed in [56]. The spacetime metric can be written, in Boyer-Lindquist coordinates (t, r, θ, φ) as where andˆ are expressed by = r 2 + a 2 cos 2 θ,ˆ = r 2 + a 2 − 2M(1 + λ 2 )r. (40) It contains a family of parameters where a and M corresponds to the spin and the mass of wormhole, and λ 2 is the which represents a special region that connects two different asymptotically flat regions. To simplify the problem, we can consider the deflection in the equatorial plane in a linearized Kerr-like wormhole metric in spin a. Such a metric can also describe the gravitational field around the rotating star or planet given by

Deflection angle
We follow the same GBT approach to calculate the deflection angle of Kerr-like wormhole as discussed in previous section.
Here we use the following coordinate transformation after substitution the metric takes the following form The last equation represents an isotropic metric, so that the line element can be written in the form The isotropic coordinate speed of light v(ρ) can be found from the relation Using n(ρ) = c/v(ρ) and c = 1, the last equation yields the effective refractive index for light in the Kerr gravitational field In this way the optical metric reads By inverting the coordinate relation (43), we find which suggests that in leading order terms ρ r −M(1+λ 2 ). The refractive index relation for massive particles in that case reads We can apply this expression to our Kerr optical media which results with in which Eq. (43) has been used. Therefore, the quantity dϕ/dt, in the case of massive particles is modified as follow yielding This result clearly indicates that the refractive index is modified due to the angular momentum parameter a. In particular, in leading order terms, we find Going through the same procedure we find From the optical metric is quite easy to see that at a very large distance Now from the GBT we find the following expression for the deflection anglê This integral can easily be evaluated, yieldinĝ Note that the plus and minus signs refers to the retrograde and prograde light ray, respectively. Setting w = 1, in the last equation we recover the deflection angle of light (Fig.  2).

Geodesics approach
In this section, we will derive the equations of motion of an uncharged test particle of mass μ and then we will calculate the deflection angle of the particle moving in the rotating naked singularity or the Kerr-like wormhole spacetime, given by (4) and (39), respectively. The geodesic equations may be derived using Hamilton-Jacobi approach, as we look for a function of the coordinates and the geodesic affine parameter σ , which is a solution of the Hamilton-Jacobi equation In general, the solution S of the Hamiltonian-Jacobi equation depends on four integration constants and satisfy the relations for the conjugate momenta which can be defined by the Lagrangian of the system L = 1 2 g μνẋ μẋ ν through the equations p μ − ∂L/∂ẋ μ = 0. Here, the dot represents differentiation with respect to the geodesic affine parameter σ . In order to involve the mass of the test particle, σ must be related to the proper time by τ = μσ , which turns the conjugate momenta into where u μ ≡ dx μ /dσ stands for the four-velocity of the particle with normalization condition g σρ u σ u ρ = 1. Therefore, one can obtain the Hamiltonian which does not depend explicitly on geodesic affine parameter σ and is a constant of motion. By taking zero or positive values of μ 2 , the Hamiltonian (63) can be used to give null or timelike geodesics, respectively. The corresponding Hamilton-Jacobi equation is where S is the Jacobi function (59). From the symmetries one can obtain two constants of motion corresponding to conservation of energy, E, and angular momentum with respect to the rotational symmetry axis of the naked singularity or Kerr-like wormhole, J . Therefore, we have It follows that for the time and azimuthal coordinates we obtain the equationṡ In addition, we have the constant of motion given by Eq. (63), corresponding to the conservation of the rest mass of the particle. The energy E and the angular momentum of the particle J can be expressed by relativistic velocity w and impact parameter b at infinity, where the spacetime is asymptotically flat, as it is shown by Eqs. (20) and (21).
Hereafter, we will consider only a light ray laying on θ = π/2 hypersurface of the considered rotating naked singularity or wormhole geometry expecting the largest ray deviation. Therefore, we shall seek a solution of the Hamilton-Jacobi equation of the form Taking into account the relationship between a metric tensor and its inverse, thus the nonzero contravariant metric coeffi-cients are The Hamilton-Jacobi equation than takes the form Using the relation p r = g rrṙ = ∂ S/∂r = d S r /dr, it may be related directly to the radial equation of motion of a particle with rest mass μ dr dσ 2 = g φφ E 2 + 2g tφ E J + g tt J 2 g rr (g 2 tφ − g tt g φφ ) where the solutions of the equation in E represent the effective potentials Evaluating the Hamiltonian at the minimum distance of approach of the particle, whereṙ = 0, we find the precise relation between the impact parameter ξ = J/E = wb and the closest approach distance r 0 which reduces to expected photon impact parameter b for stationary axially symmetric metrics, when w = c = 1 [74]. The sign in front of the square root is chosen to be positive when the orbital motion of the particle and rotation of the compact object are in the same direction. The orbital equation of the particle in terms of azimuthal φ and radial r coordinates can be obtained by using the fact that the partial derivative of the Jacobi function with respect to the constant of motion J is a constant Solving Eq. (69) or taking into account Eq. (70), we obtain the azimuthal shift function of the particle by Eq. (74) where the sign in front of the integral is changing when the particle passes thought a turning point r 0 of the orbital motion. V ± (r, ξ) ≡ V ± (r, J )/E stand for the effective potentials per unit of the energy and depend directly on the impact parameter ξ . It can be verified that an Eq. (75) is consistent with Eqs. (66) and (70).
Having the azimuthal shift of the photon as a function of the coordinate distance r , we can calculate the deflection angle of the test-particle as we assume that the source r S and observer r O are placed in the asymptotically flat region of the spacetime, so that r S → ∞ and r O → ∞. Then, after calculating the total azimuthal angle according to the orbital equation (75), the deflection angle iŝ

Deflection angle by rotating naked singularities
In order to calculate the deflection angle of a test-particle in the spacetime of rotating naked singularities, firstly we will calculate the impact parameter of the test-particle ξ = wb as a function of the distance of closest approach r 0 . According to Eq. (73) the equatorial impact parameter of a particle with relativistic velocity w is given by The impact parameter (78) reduces to that for photons, obtained in [75], letting w = 1. For a > 0 the naked singularity rotates counterclockwise, while for a < 0 the naked singularity and the particles rotate in converse direction. Expecting small deviations in the orbital motion of the test-particles, we will assume that the distance of closest approach r 0 is of the same order as the impact parameter b due to the validity of inequalities, b m and b a. Therefore, considering Eq. (70) in the framework of the assumptions adopted, we suggest the following solution where we are introducing two independent expansion parameters, as the former is in term of metric parameter m related to the ADM mass M and scalar charge q, while the latter is in term of specific angular momentum a of the rotating naked singularity Here the coefficients c m,a are real numbers, and the summation is over all possible combinations of epsilon powers i, j up to and including second order terms. Applying the procedure to the equationṙ (r 0 ) = 0 as discussed above, we acquire the following result The calculation of the deflection angle of massive particles in the spacetime of rotating naked singularities requires integration of the first derivative of the azimuthal shift function dφ/dr, as follows from Eq. (75). Since, we are expecting small deviation in the orbital motion of the particle, we will consider only terms with the most significant contribution, for simplicity. Introducing new small enough parameters ε m = m/r 0 and ε a = a/r 0 , we will perform power series expansion of dφ/dr, in the spacetime under consideration (3), up to and including second order terms of ε m and ε a . After the convenient change of variables r = r 0 /u, we obtain dφ du where we are introducing the quartic polynomial in powers of u Completing the integration of (82) from the minimal distance of closest approach r 0 to observer and source positioned at infinite distance with respect to the center of the compact object, we obtain the following approximate expression for the deflection angle of a massive particle, according to Eq.
The most significant contribution to the deflection angle of massive particles is given by massive term proportional to O(ε m ), while the order proportional to angular momentum O(ε a ) is vanishing. The next second order terms proportional to the mass O(ε 2 m ) and to the mass and angular momentum O(ε m ε a ) give also a contribution, as opposed to the second order term proportional to the square of the angular momentum O(ε 2 a ), which is missing. Moreover, taking into account the power series (81), we can represent the deflection angle of the relativistic particles in terms of the impact parameter b up to and including the second powers of the small parameters m and a . After performing the calculations, we get the following invariant representation of the deflection for a given relativistic velocity w of the massive particlê By comparison of the calculated approximate expression for the deflection angle of massive particles (84) via geodesic approach with the result (38) found by GBT it appears that up to the second formal order in the expansion parameters the both angles coincide. Letting w = 1, we recover the light deflection anglê presented in the static case in [5,75,76]. Similarly, the plus (minus) sign represents the retrograde (prograde) light ray, respectively.

Deflection angle by Kerr-like wormholes
In the case of the Kerr-like wormhole (39) the impact parameter of the massive particles ξ = wb with relativistic velocity w and the minimal distance of closest approach r 0 obey the following relation Since, we are considering propagation of massive particles in the asymptotically flat spacetime regions with geometrical curvature, causing small deflection angles of trajectories, then r 0 is expected to be of the same order as the impact parameter b. Therefore, under assumptions b M and b a, we are taking advantage of the two small enough expansion parameters where M and a are the mass and spin of the Kerr-like wormhole, respectively. Hence, we could find an approximate solution of the radial Eq. (70), as for itself we suggest the ansatz Here the coefficients c M,a are real numbers. Solving the equationṙ (r 0 ) = 0, for the closest approach distance we obtain the following power series expansion Now, we can proceed ahead toward finding of the perturbative expansion of the deflection angle of a massive particle moving in the spacetime of the Kerr-Like wormhole. Since, we are interested in the weak deflection limit, we will perform power series expansion of the derivative dφ/dr, adopting the new small parameters ε M = M/r 0 and ε a = a/r 0 up to and including second order terms. In order to perform the integration afterwards more easily, we introduce a new variable u = r 0 /r in terms of which we obtain dφ du where it is introduced the fourth degree polynomial in variable u Therefore, the bending angle of the light ray, defined by Eq. (76), up to and including the second orders of the expansion parameters ε M and ε a is given bŷ Taking advantage of the power series (89) we can represent the deflection angle of the relativistic particles in terms of the impact parameter b, up to and including second order terms of the small parameters M and a . Therefore, the deviation of trajectories of the massive particles by the Kerr-like wormhole leads to the following components of the deflection anglê As in the previous case of the rotating naked singularities, here as well the most significant contribution to the deflection angle is given by the massive term proportional to O( M ).
The second order terms proportional to the mass O( 2 M ) and to the mass and angular momentum of the wormhole O( M a ) give a contribution to the deflection angle of the relativistic massive particles. The first and second order terms based on the wormhole's angular momentum, O( a ) and O( 2 a ) respectively, do not contribute in the perturbative analysis. In the special case w = c = 1, from expression (92) the deflection angle of light is recovered [63] Comparing between the obtained results for the deflection angles (58) and (92) via the GBT and geodesic approach respectively, shows coincidence of the perturbative expressions up to the leading terms proportional to O( 2 M ) and O( M a ).

Conclusions
We studied the gravitational deflection of massive/massless particles in a rotating naked singularities and Kerr-Like wormholes backgrounds. In the first part of this paper, we used a recent geometric approach to compute the gravitational deflection angle using the refractive index of the optical media and by utilizing the GBT applied to the optical metric. In order to incorporate the effect of relativistic speed on the deflection angle, we have assumed the propagating massive particles as a de Broglie wave packets yielding a modification of the refractive index N , and a constant quantity of a given optical media, λN = const. In particular, we have shown that the refractive index and the Gaussian optical curvature are affected by the mass M, angular momentum parameter a, naked singularity parameter γ , the deformation parameter λ, and the relativistic velocity of the particle w. Although the geodesic curvature κ(C R ) is slightly modified due the velocity of the particle, i.e. κ(C R ) → w R −1 , we have shown that the expression for the deflection angle readŝ In other words we have found exactly the same expression as in the case of deflection of light, which suggests that the above result can be viewed as a form-invariant quantity for asymptotically flat spacetimes. Again, the total deflection angle is found by integrating over a domain outside the light ray described by r γ , hence the global effects are important. Note that in our setup we have simplified the problem by considering a linearized spacetime metrics and we have also assumed that the source and the observer are located at infinity. Applying then the GBT to the isotropic metrics we have found two interesting results for the deflection angles, namelŷ α Naked sing. 2mγ in the case of rotating naked singularity, and in the case of Kerr-like wormholes. From the last two equations we see that an apparent singularity appears when w → 0, therefore an additional constraint should be imposed, namely 0 < w ≤ 1. Now we specialize further the case of λ = 0 and γ = 1, in that case both results coincide with the Kerr deflection angle for massive particles reported recently in [52], but also in [77,78]. Note that our approach works perfectly well for relativistic particles, in this sense the apparent singularity is not problematic. Furthermore, we did a careful analyses by studying the problem of deflection of particles using the Hamilton-Jacobi equation and we find agreement of results in leading order terms. Our results show that, while the deflection angle decreases with γ in the case of JNW spacetime, contrary to this, the deflection angle increases with the deformation parameter λ 2 in the case of Kerr-like wormholes as can be seen from Fig. 3. Our results are not important only from conceptual point of view, namely the difference between deflection angels may shed light from observational point of view in order to distinguish the two spacetimes. Finally, we plan in the near future to extend these results by including finite distance corrections on the deflection angle, in other words by assuming that the distance between the source and the observer is finite.