On the trajectories of null and timelike geodesics in different wormhole geometries

The paper deals with an extensive study of null and timelike geodesics in the background of wormhole geometries. Starting with a spherically symmetric spacetime, null geodesics are analyzed for the Morris-Thorne wormhole(WH) and photon spheres are examined in WH geometries. Both bounded and unbounded orbits are discussed for timelike geodesics. A similar analysis has been done for trajectories in a dynamic spherically symmetric WH and for a rotating WH. Finally, the invariant angle method of Rindler and Ishak has been used to calculate the angle between radial and tangential vectors at any point on the photon's trajectory.


Introduction
In general relativity, a wormhole (WH) is considered to be a tunnel through which two distant regions of spacetime can be connected [1]. Long back in 1916, Flamm [2] introduced the idea of wormhole, analyzing at that time the recently discovered Schwarzschild solutions. In 1935, Einstein and Rosen [3] constructed WH type solution considering an elementary particle model as a bridge connecting two identical sheets. This mathematical representation of space being connected by a WH type solution is known as "Einstein-Rosen bridge". Wheeler [4,5] in the 1950s considered WHs as objects of quantum foam connecting different regions of spacetime and operating at the Planck scale. Subsequently, using this idea, Hawking [6] and collaborators introduced the idea of Euclidean wormholes. But these types of WHs are not traversable and, in principle, would develop some type of singularity [7]. However, these hypothetical shortcut paths, i.e., traversable WHs, have been rekindled by the pioneering work of Morris and Thorne [8] which is considered as the modern renaissance of WH physics. Subsequently, it was claimed that there is no strong ground [9,10] for the energy conditions and hence one considered WH, with two mouths and a throat, to be an object of nature, i.e., an astrophysical object.
On the other hand, in general relativity, WH physics is a specific example where the matter stress-energy tensor components are evaluated from the spacetime geometry by solving Einstein's field equations. But for a traversable WH, the stress-energy tensor components so obtained always violate the null energy condition [1,8]. As the null energy condition (NEC) is the weakest of all the classical energy conditions, its violation signals that the other energy conditions are also violated. In fact, they violate all the known pointwise energy conditions and averaged energy conditions, which are fundamental to the singularity theorems and theorems of classical black hole thermodynamics. Generally, it is believed that a classical matter obeys energy conditions [11] but, in fact, it is known that they also get violated by some quantum fields (namely as regards the Casimir effect and Hawking evaporation [12]). Further, for a quantum system in classical gravity, it is found that the averaged weak or null energy condition(ANEC), which states that the integral of the energy density as measured by a geodesic observer is non-negative, could also be violated by a small amount [13,14].
Finally, it is worth to mention a few important dynamical WH solutions. Hochberg and Visser [15] and Hayward [16] independently formulated the dynamical WH solutions, choosing a quasi local definition of the WH throat in a dynamical spacetime. Accordingly, WH throat is a trapping horizon [17] of different kind but again matter in both of them violates the NEC. On the other hand, Maeda,et al. [18] have developed another class of dynamical WHs (cosmological WHs) which are asymptotically FRW spacetimes with big bang singularity at the beginning. This class of WHs contain matter which not only obey NEC but also the dominant energy condition everywhere. These two classes of dynamical WHs are distinct from the geometrical point of view. For the former one, the WH throat is a 2D surface of non-vanishing minimal area of a null hypersurface, while for the later one, there is no past null infinity due to the initial singularity. Hence, the WH throat is defined only on a space-like hypersurface and the spacetime is trapped everywhere without any trapping horizon [19]. Recently, Lobo et al. [20][21][22] formulated wormhole solutions which are dynamically generated using a single charged fluid. Also, dynamical WHs are considered with a two-fluid system [23,24], for a matter distribution relevant to present day observations [25] and using the mechanism of particle creation [26]. Then for evolving WH 1 , one may refer to Refs. [27][28][29][30][31].
The particle motion in wormhole spacetimes is an important issue related to traversable WHs. It is interesting to examine whether a timelike or null geodesic can tunnel through the throat of the WH. Cataldo et al. [32] studied motion of test particles in the background of zero tidal force Schwarzschild-like WH spacetime. They showed that particles moving along the radial geodesics reach the throat with zero tidal velocity in finite time while the particle velocity reaches maximum at infinity if it travels along a radially outward geodesic. For non-radial geodesics on the other hand, the particles may cross the throat with some restrictions. Olmo et al. [33] carried out a detailed investigation of the geodesic structure for three possible WH configurations, namely: Reissner-Nordström-like WH, Schwarzschild-like WH and Minkowski-like WH. They have shown that it is possible to have geodesically complete paths for all these WH spacetimes. Culetu [34] examined both timelike and null geodesics for a WH belonging to the Planck world ( WHs whose throat size is of the order of the Planck length l P ) where quantum fluctuations are supposed to exist and the spacetime smoothness seems to break down. Muller [35] also studied null and timelike geodesics in WH configuration using elliptic and Jacobian integral functions. He showed that it is possible to connect two distant events geodesically. Regarding geodesic study in non-static WHs, recently Chakraborty and Pradhan [36] have studied the geodesic structure of the rotating traversable Teo WH. Also, Nedkova et al. [37] discussed the shadow of a class of rotating traversable WH in the framework of general relativity. They showed that the images depend on the angular momentum of the WH and the inclination angle of the observer. Finally, it is worthy to mention the work of Ellis [38]. He constructed a static, spherically symmetric, geodesically complete, horizonless spacetime manifold with a topological hole (drainhole) at its center by coupling the geometry of Schwarzschild spacetime to a scalar field. It is found that on one side of the drainhole the manifold is asymptotic to a Schwarzschild manifold with positive mass parameter 'm', and on the other to a Schwarzschild manifold with negative mass parameter 'm', with the condition −m > m. As a consequence, there is attraction of particles on one side while there is repulsion on the other side (with higher strength).
The present work presents a detailed investigation of both timelike and null geodesics both for static and dynamical WHs. The paper is organized as follows: Sect. 2 deals with static spherical WHs in which null and timelike geodesics are studied in great detail. A similar geodesic analysis is presented for dynamical WH in Sect. 3 and rotating WH in Sect. 4. Sect. 5 uses the invariant angle method of Rindler and Ishak to calculate the angle between radial and tangential vectors at a point on the photon's trajectory. Finally, the paper ends with a short discussion and concluding remarks in Sect. 6. Throughout our analysis, we have chosen to work with wormholes whose material extends all the way from the throat out to infinity.

Trajectories in a spherically symmetric and static geometry
The metric for a general spherically symmetric and static metric can be written as (Ref. [39,40]), where, lim An important relation between momenta one-forms of a freely falling body and the background geometry is given by the geodesic equation [39], dp β dλ where λ is some affine parameter. This relation tells us immediately that if all the components of g αν are independent of x β for some fixed index β, then p β is a constant along any particle's trajectory, i.e., a constant of motion. Now, if we work in the equatorial plane by setting θ = π/2, then, in Eq.(1), all the g αβ become independent of t, θ, φ (cyclic coordinates). That means that we can find the respective Killing vector fields δ ν α ∂ ν with α as cyclic coordinates. Now, since p t and p φ are constants of the motion, we will set them as where E is the energy and L is the angular momentum of the photon or a particle as measured by observers at asymptotically flat regions far from the source. Thus, we get , and let, p r = dr dλ =ṙ (4) where the dot represents the derivative w.r.t. some affine parameter λ.

Null geodesics
Now, for null-geodesics, we have p α p α = 0. Thus, Using Eq.(4) and (5), we can write the equation of the photon trajectory in terms of the impact parameter, µ = L/E, as: If we assume that the geometry is caused by a source of radius r s , then the photon coming from infinity will not hit the surface if there exists a solution r o > r s for whichṙ 2 = 0. We then call r o as the distance of closest approach or the turning point. In that case, The impact parameter then becomes, Using Eq. (6), we can write, Now, if a photon coming from the polar coordinate lim r→∞ (r, −π/2 − α/2) passes through a turning point at (r o , 0) before approaching the point lim r→∞ (r, π/2 + α/2), then this α, which is a function of r o , is what we refer to as the deviation/deflection angle, given by (Ref. [41] α(r o ) = −π + 2 ∞ ro B(r)dr However, it is possible that a photon might get trapped in a sphere of constant r and thus may not approach lim r→∞ (r, π/2 + α/2). In that case, the integral will diverge. Such spheres are called photon spheres; they are discussed in sec.2.3.

Morris-Thorne wormhole
The Morris-Thorne wormhole metric(Ref. [8]) is given by, where Φ(r) is the redshift function and b(r) is the shape function of the wormhole for which b(r) ≤ r and equality holds only at the throat. Both the functions are such that they also satisfy asymptotic flat conditions. Thus, the equation of the trajectory for null geodesics, Eq.(6), becomes However, note that the coordinate r cannot be used for describing the whole spacetime since it accounts for a coordinate singularity at the throat and is therefore valid for describing geometry only at one side of the throat. Thus, for geodesics that actually reach and pass through the throat, one should not use this formula for the trajectory equation. Instead, one can always work with the proper distance(l) which must be valid everywhere and throughout the wormhole. As an example, for the metric given in Eq.  where, by definition, this proper radial distance is positive for the upper universe, negative for the lower universe and is zero at the throat. Using this, Eq.(6) can be generalized for wormholes as: where we have substituted r in terms of l which, in principle, could be obtained by inverting Eq.(13) to get r ≡ r(l). In this paper, however, we will mostly be interested in the behavior of trajectories on one side of a throat and so we will mostly work with r for our convenience. Now, using Eq.(10), the null-geodesics coming from infinity and not reaching the throat gets deflected by an angle, It turns out that, for stationary observers in the r, θ, φ system, the radial tidal forces can be made to vanish if we have Φ (r) = 0, which we can do by simply choosing Φ(r) ≡ 0, say. This condition gives us a simple class of solutions and corresponds to precisely zero tidal forces. Using Eq. (8), it can also be deduced that for these wormholes, light can reach the throat only if |µ| < b o , where µ is the impact parameter and b o is the radius of throat. Thus, for these ultra-static wormholes, the light deflection angle becomes, Now, for an asymptotically flat geometry, a good choice for the function b(r) is, where b o = b(r t ) = r t corresponds to the throat radius and n=2 gives us the famous Ellis-wormhole [38]. We will call this parameter 'n' the shape exponent. Now, the deflection angle for this choice of b(r) in terms of r o and n becomes, We can see how the deflection angle depends upon the value of shape exponent and the distance of closest approach as given in Fig. 1. For the Schwarzschild Metric, the deviation angle becomes, It turns out that, for Ellis wormhole, we can write the exact expression for α(Ref. [42]), as where we have written |µ| = r o . In the weak-field regime where |µ| << b o , the deflection angle becomes, Now, we wish visualize the geometry of such a wormhole, i.e. Eq.(11), (Ref. [8]). Since the geometry is spherically symmetric and static, we can confine our attention to an equatorial slice through our wormhole at any instant in time. Then the metric becomes Now to visualize this 2D geometry, we can embed it into a higher dimensional space of IR 3 ,i.e., ordinary 3D Euclidean space whose metric can be written in the form: Then we can show that, For a wormhole of type Eq. (17), it becomes Thus, after integrating, the embedding function becomes(Ref. [43]), where 2 F 1 represents the hypergeometric function. Note that due to spherical symmetry, there is no dependence on the coordinate φ. Thus, we can easily visualize the geometry through a 2D plot as in Fig. 2 which represents the embedding function for various values of the shape exponent and throat radius. We can also study how proper length l depends upon the radial coordinate r and the shape exponent n. If we keep b o = 1, then we have A comparision of proper lengths for different geometries is shown in Fig.3.  (b) The plot shows how dr/dl varies with r. Figure 3: In both plots, it can be seen that with increasing r the slope tends asymptotically to 1, thus corresponding to an asymptotically flat geometry. Also, the slope approaches 1 even faster when the value of n increases. Note that the throat radius has been set to b o = 1.

Photon spheres
A photon sphere is a location where the curvature of spacetime is such that even null geodesics can travel in circles. In other words, a photon sphere is a region where bothl andl vanish for a photon. We will calculate the possibility for such a region considering a general static metric as given in Eq.(1). Rewriting Eq.(5) in terms of l gives,l Thus, a photon sphere will exist at a location where, Another way of looking at it is that, it is the region for which the deviation angle of a photon diverges. Thus, for Morris-Thorne wormhole, Eq. 15 implies that the throat itself is a photon sphere which is evident from Fig.  1a. Also, by Eq. 30, any other place which satisfies the following condition contains a photon sphere, rΦ (r) = 1 {condition for photon sphere to exist at any r > b o }

Timelike geodesics
For timelike-geodesics, we have p µ p µ = −m 2 where m is the mass of the particle. If we define the quantities E and L as the energy per unit mass (E/m) and the angular momentum per unit mass (L/m) respectively, theṅ Notice that a timelike particle always reaches the throat with zero radial velocity, independent of the value of the impact parameter µ (∵ B(r) −1 = 0 at throat). Now, the general equation of the trajectory becomes where µ = L/ E and l is the proper length. If we differentiate Eq.(32) with respect to the affine parameter, we obtain for the second derivative of the radial coordinatë The dependence on L 2 is a consequence of spherical symmetry as it tells that the orientation of the angular momentum does not affect the radial acceleration. For a Morris-Thorne wormhole, it becomes For a particle with zero initial velocity (ṙ = L = 0),r ∝ −Φ (r). Thus in ultra-static wormholes, a particle stays at the same position if not given any initial velocity. Also, at the throat,ṙ = 0 also impliesr = 0. Thus, a particle reaching throat not only attains a zero radial velocity but also has vanishing radial acceleration. The expression for the radial acceleration in an ultra-static wormhole with shape exponent reduces to: The case n=1 is studied in detail in Ref. [32]. However, note that if L = 0, thenr ≡ 0 for E = 1, while on the other hand,r > 0 for E > 1. Thus, this family of geometries correspond to repulsive gravity.

Unbounded orbits
If the particle falling from infinity does not hit the throat, it will get deflected after approaching a closest distance of r o , where r o is then the real solution of the equation, Using Eq. (33) and Eq. (37), we can then write Now, if the particle does not fall into the throat, the total deflection angle (α) for a particle falling from infinity will be, For Morris-Thorne wormhole, it becomes,

Bounded orbits
For a Morris-Thorne wormhole, Eq.(32) becomes, For ultra-static wormholes, we can simply write where, dl is the differential proper length and V 2 (l) can be thought of as the effective potential. This case is studied in detail in Ref. [43]. However, we will choose a different form of e 2Φ(r) and will try to study the trajectories it allows. Let us define where (r) is a continuous function which is significant only near the throat and is vanishingly small otherwise. Now for this choice, we can write a simplified form of Eq.(41), for distances far from the throat, thus: Note that we should not choose (r) ≡ 0, because then the throat of the wormhole will be a horizon which will make the wormhole non-traversable. Now, we can define an effective potential, Therefore, which tells us immediately that the allowed region for a particle with energy E (as measured at infinity) can be determined from the inequality : In other words, the radial range of a particle, depending upon its conserved energy E, is bounded within those radii for which V is smaller than E. Also note that, since r > b o , we must have lim r→ bo It is important to note that any bound orbit that exists around a spherically symmetric source can be of only two types. It can be either a circular orbit (stable or unstable) or an orbit that oscillates around the radius of a stable circular orbit (Ref. [44]). So, let us study the possibility of circular orbits in our geometry.

Circular orbits
Now for circular orbits, we require that bothṙ andr vanish for at least some r. Therefore, It means, for circular orbits, that the energy of a particle should be an extremum of the effective potential. Precisely, if the conserved energy corresponds to a maximum or a saddle point of the potential, then it will be an unstable orbit, while if it corresponds to a minimum of the potential, it will be a stable orbit. Now, if we choose b(r) = b n o r 1−n as described in Eq. (17), we can write After simplification, we get where we have defined, Also note that, Now if r c is some real root of the Eq. (49), then a circular orbit is possible only when Let r p be the point where the first derivative vanishes. Then, Now, we will consider three cases: Case I : n > 2 For this case, it is easy to see that Thus the behavior of f(r) is such that it will start from a negative value at r = 0 and will grow further negative with the increase in r until it hits a turning point at r = r p , after which it increases monotonically. Thus, it can be inferred that f (r) will have only one positive real root r c say. Then it is clear that f (r) > 0 ∀ r > r c . Thus, from Eq.(51), we can say that this root must also satisfy condition III. Since there is only one turning point of f (r), it is obvious from Eq.(48) that this corresponds to the maximum of the potential, in which case it will always lead to an unstable orbit. Hence, for n > 2, there will be only one unstable circular orbit for a particular value of E and L.
Case II : n = 2 For n=2, f (r) becomes As we can see, if the conserved angular momentum L is larger than the throat radius, then we definitely have one root, r c , which satisfies condition III. By the same argument as above, it is clear that it must correspond to a maximum of the potential which can only lead to an unstable circular orbit. Hence, for n=2, there is a possibility of only one unstable circular orbit depending upon L and b o and E.
Case III : 0 < n < 2 For this case, we have lim  where equality holds for f (r p ) = 0 and inequality for f (r p ) > 0. Note that when f (r p ) = 0, then, by using Eq.(50) & (54), we also have d 2 V 2 (r p )/dr 2 = 0. Thus, it will correspond to a saddle point of the potential at r = r p . For the condition given in Eq.(59), we always have r p > b o which means condition III is also satisfied. So, we will have the possibility of one unstable circular orbit for this case.
Meanwhile if f (r p ) > 0, we will have two real roots. Now, from Eq. (48), it is clear that the smaller of these roots will correspond to a local maximum and the larger root will correspond to a local minimum. And, using Eq.(51) and the condition in Eq.(59), it can be inferred that both of these roots will also satisfy cond. III.
Hence, for 0 < n < 2, there is a possibility for one unstable circular orbit or a combination of one unstable circular orbit and a stable circular orbit. Note that this is the only case where we have the possibility of stable circular orbits. It is also interesting to note that the Schwarzschild geometry, for which n=1, lies in this case.
As mentioned before, any bound orbit, which is not circular, is possible only when it oscillates around the radius of a stable circular orbit. Thus, we can say that we surely have no non-circular bound orbits when n ≥ 2 for any E and L of the particle.
For the Schwarzschild case (n=1), we can substitute b o = 2M (Schwarzschild radius) with r satisfying r > b o . Then, the condition for circular orbit,Eq.(59), becomes which we know is the correct limit for the Schwarzschild case. Thus, what we have done in this section is a general treatment for any shape exponent. But physically, we can say that (Ref. [45]), where, V g (r) is the gravitational potential for a Newtonian like gravitational force given by, So, all our conclusions are valid for this interesting analogy as well. Hence, we have proved, using GR, that in a universe where "Newtonian like gravity" dies out as r −3 or faster, no stable orbits are possible. In other words, the existence of planets will itself be almost impossible.

Time period of circular orbits
We have seen that there is at least one unstable circular orbit possible for any value of the shape exponent. So now, we will try to calculate the time period of a circular orbit of a particle at a distance r c in terms of n and b o . Using Eq.
It means that the velocity required for a satellite to be set in an orbit of radius r c around a wormhole is given by As we can see, the time period is always proportional to the radius of the orbit but is inversely proportional to the throat radius. The latter condition signifies that increasing the throat radius can be thought of as keeping the throat radius fixed but decreasing the radius of the circular orbit itself; in which case it is logical that its time period will decrease. Again, we can recognize Eq.(65) as the generalization of Kepler's third law for an attractive force law given by Eq. (61). It can be proved immediately by scaling arguments if we put, say, r = λr and t = µt in Eq.(61), to get µ 2 ∝ λ n+2 (Ref. [46]). So, we can retrieve the Kepler's third law in its original form by putting n=1, so that For n=2, the time period becomes: Choice of (r) As we have mentioned, Fig.4 is not valid near the throat as it does not consider the significance of (r) in that region. Now, we shall try to guess a physically reasonable form for (r). First, let us consider the problem of tidal forces. For a spaceship whose one end is at r = a and the other end is at r = b, the magnitude of the tidal force experienced by the ship would just be the difference between the forces at r=a and r=b. If the spaceship is far away from the throat where (r) << 1, then the tidal force can be written in terms of the effective potential as: where τ represents the tidal force. The above expression for the tidal force is obtained by considering the gradient of the potential throughout the spaceship and using the fact that the gradients, except at the two ends of the spaceship, should almost vanish for it to not deviate much from its rigid structure. Pictorially, it is just the absolute difference between the slopes at the two points on V 2 (r) vs. r curve in Fig.4. It can also be noted that without considering the significance of (r), we get a discontinuity in the slopes of the potential at the throat. That would correspond to an impulse of force which will be experienced by a particle at the throat while traversing through the wormhole. It would be like hitting a thin membrane of a tough material. However, it should be noted that even in the presence of infinite tidal forces, causal contact is never lost among the elements making up the observer; this suggests that curvature divergences may not be as pathological as traditionally thought (Ref. [47]). We would definitely want to remove this problem by choosing a reasonable (r). If we want a particle to smoothly traverse through the wormhole without any impulse of force, then we can do so by demanding that slope of V 2 (r) goes to zero as it reaches the throat. It would imply that lim r→ bo And since we want it to vanish for large r, a simple choice of (r) might be, The larger the κ, the faster it will die out. We can sketch the plot forṙ 2 vs. r as shown in Fig.5.
It is clearly visible that we have removed the discontinuity in the slope at the throat. Also note that there exists a local minimum ofṙ 2 at the throat. Such a local minimum will correspond to an unstable bound orbit. Thus, for our choice of (r), the throat will correspond to a region of unstable circular orbits.

Trajectories in a dynamic spherically symmetric wormhole
The metric for a spherically symmetric and dynamic Morris-Thorne wormhole can be written as [Ref. [26]], which corresponds to a 3-geometry with a time dependent scale factor a(t).

Null geodesics
Due to spherical symmetry, we should expect to find the same answer for the deflection angle as that of the static case. For geodesics in equitorial plane, θ = π/2 & p θ = p θ = 0. And since φ is a cyclic coordinate, p φ is a constant of motion. So let, For null geodesics, ds 2 = 0, where, E is a positive constant of integration. Substituting this into Eq.(73), we getṙ And sinceφ = L/a 2 r 2 , the equation of trajectory becomes, where, µ = L/E. Note that the time-independence of this equation is just an artifact of our poorly chosen coordinate system. This is because r is itself a comoving coordinate. We should define a new coordinate, r (r, t) = a(t).r, so that any surface r = const., t = const. is a two-sphere of area 4πr 2 and circumference 2πr . This coordinate r can then be called as the 'curvature coordinate'. In this coordinate, equation of trajectory becomes, However, the total deflection angle can be calculated using the coordinate r by the following equation, which is same as Eq. (16). Also, light will reach the throat only if |µ| < b o as in the static case. We can make above conclusions due the fact that the geometry, inspite of a time-dependent scale factor, is always spherically symmetric.

Timelike geodesics
For simplicity, we will work for the ultra-static case, i.e., in which Φ(r) = 0. Then, for timelike geodesics, we have And the time-component of geodesic equation becomes, Now, from the above two equations, we getṫẗ which upon integration gives,ṫ where E is a constant of integration. Using it, we geṫ Thus, we will get the same equation of motion as Eq.(80). And the deflection angle will be same as Eq.(81).

Trajectories in a rotating wormhole
The metric for a rotating wormhole can be written as (Ref. [28]), where N, K, ω and µ are functions of r and θ, and ω(r, θ) may be interpreted as the angular velocity dφ/dt of a particle that falls freely from infinity to a point (r, θ). Assume that K(r, θ) is a positive, non-decreasing function of r that determines the proper radial distance R, i.e., R ≡ rK. We also require this metric to be asymptotically flat, which implies However, the metric (86) was initially derived for slowly rotating stars [48] and hence it implicitly assumes the absence of effects due to centrifugal forces [49]. Now, at the equatorial plane, the metric becomes, We can compare it with the metric far from a rotating source of mass M and angular momentum S as given by (Ref. [40]), In cylindrical coordinates, x 1 = r cos φ, x 2 = r sin φ, x 3 = z. Assuming axial symmetry, only S 3 term survives. Let's call it J. Then, Thus, in asymptotically flat limit, comparing it with the g tφ metric term of the (88), we get Now since the metric terms in (88) are independent of t and φ, the corresponding momenta one forms are conserved. Thus, we can write (Ref. [50]) Thus, the expressions for E and L in terms ofṫ andφ becomeṡ

Null geodesics
For null geodesics, using Eq.(88) and Eq.(93), we get ds dλ It can be rewritten asṙ where V ± are the roots of the equation Now, since r, K(r), N (r) are all non-negative functions, we get Now, a photon will make its closest transit from the wormhole at a distance r o if at that point the condition, E = V ± (r o ), is satisfied. If there's no such point, then the photon will definitely fall into the throat. Without loss of generality, we can assume J > 0, where J is the angular momentum of the rotating wormhole. This assumption also implies that ω(r) > 0. Now, we have two possibilities for the conserved angular momentum (L) of the photon: it can be either L > 0 or L < 0. This will determine whether the light ray is traversing along the direction of frame dragging or opposite to it. Also, since we have assumed that there are no horizons, the g tt term of the metric can never change sign. Thus, So, considering the above inequality and using the fact that E > 0, the condition, E = V ± (r o ) becomes, where, ω o = ω(r o ), N o = N (r o ) and so on. If we denote r o the distance of closest approach when L < 0 and by r o when L > 0, then from the above equation, we can write If N (r) is a smooth decreasing function, then this equation proves that the distance of closest approach is greater when the light ray is moving in the direction of frame dragging than that of light moving opposite to it. Now, from Eq.(93) and Eq.(94), we can write the equation of motion of photon trajectory as, As we can see, the equation of motion of a photon along the direction of frame dragging is different from that of the opposite direction. Now, if a photon does not fall into the throat, it will get deflected according to where, According to the above equation, the deflection angle for a photon moving along the direction of the frame dragging will be larger than that of a photon coming the other way [51].

Timelike geodesics
For timelike geodesics, we haveṙ where for simplicity, we have denoted E as simply E and L as L. Now, we can rewrite the above equation aṡ where, Now, the equation of a particle's trajectory can be written So, for timelike geodesics, if it does not reach the throat, it will be follow a trajectory given by the equation, , if L < 0 (108)

Invariant angle method of Rindler and Ishak
In this section, we will calculate the angle between radial and tangential vectors at a point on the photon's trajectory by Invariant angle Method which was proposed by Rindler and Ishak (Ref. [52]). Let δ represent the radial direction and d represent the tangential direction at any point on the photon's trajectory. Let ψ be the angle between them. Then the invariant formula for cos ψ becomes For a photon coming from far left of a source and heading toward the far right while being deflected, the directions d and δ in the (r, φ) basis can be written as d ≡ (±dr, −dφ) = (∓dr/dφ, 1)dφ = (∓A, 1)dφ, where A = dr/dφ δ ≡ (dr, 0) = (1, 0)dr ⇒ cos ψ = g rr d r δ r + g φφ d φ δ φ (g rr δ r δ r ) 1/2 (g rr d r d r + g φφ d φ d φ ) 1/2 = |A| √ g rr Rewriting this in the form of tan ψ (Ref. [53]), we get For a general Morris-Thorne wormhole, it becomes In the ultra-static limit, it simplifies to It is interesting to note that the expression for ψ is independent of the shape function b(r). This independence is true in the case of dynamical and rotating wormhole geometries as well.
The expression for ψ in the Schwarzschild geometry takes the form, tan ψ = 1 6 Short discussion and concluding remarks A detailed study of particle and photon trajectories has been conducted in the background of wormhole geometry. Starting with the Morris-Thorne wormhole, null geodesics and photon spheres have been analyzed, while for particle trajectories both bounded and unbounded orbits are considered. Subsequently, both null and timelike geodesics are analyzed in the geometry of dynamic spherically symmetric WH and rotating WH. Finally, using the invariant angle method of Rindler and Ishak, the angle between radial and tangential vectors on the photon's trajectory has been evaluated. Based on the above study, we have found that in a Morris-Thorne wormhole and its dynamic and rotating counterparts, the throat itself is a photon sphere. We have also seen that in such geometries, the angle between tangential and radial vectors at any point on a photon's trajectory is independent of the shape function b(r). The geodesics in ultra-static wormholes with shape exponents have already been studied in great detail in Ref. [43]. Also, in Ref. [32] Cataldo et al. studied a Schwarzschild-like traversable WH which is obtained by putting n=1 with some slight modification. For geodesics, they showed that a test particle which is radially moving towards the throat always reaches it with zero velocity and at a finite time, while for radially outward geodesics the particle velocity tends to a maximum value, reaching infinity. However, in this paper we have shown that it is true for all possible n. Also, general conditions for non-radial geodesics were derived which are required to be satisfied in order for it to cross the throat. These results are in agreement with our study and can, roughly, be obtained by putting n = 1 in our general equations for arbitrary n. Similarly, in Ref. [38], the Ellis wormhole (n = 2) is studied in great detail including the behavior of geodesics in such geometry. For Ellis wormhole, the particles are attracted on one side and are repelled on the other and so the throat is of saddle nature. In our paper, we have mainly stressed on the geodesics that remain on one side of the wormhole, unlike the above mentioned references where geodesics through the throat are studied in detail. Furthermore, we have analyzed the possibility of bounded timelike orbits for different shape exponents in a different WH geometry, which can be regarded as the generalization of the Schwarzschild geometry far from the throat. There, we used the fact that any bounded timelike orbit in a spherically symmetric geometry is either a circular orbit or an orbit that oscillates around the radius of a stable circular orbit. For this geometry, we found that, for a wormhole with shape exponent n > 2, there always exists the possibility of one unstable circular orbit while for n = 2, there exists one unstable circular orbit only when L > b o and no bound orbits otherwise. That means, no non-circular bound orbits exist when shape exponent n ≥ 2 for any value of the impact parameter. For 0 < n < 2, we found that depending upon the value of L it can either have the possibility of one unstable circular orbit or a combination of one unstable circular orbit and a stable circular orbit. While studying trajectories in a rotating wormhole geometry, we have seen that the equations of motion of both photon and particle depend upon whether it is traveling in the direction of frame dragging or opposite to it.