{'fraction_non_alphanumeric': 0.03369179163378058, 'fraction_numerical': 0.011000394632991317, 'mean_word_length': 4.404692082111437, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 6, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The Ricci Rotation Coefficients in the description of trajectories of spinning test particles off-equatorial planes in a rotational gravitational field

We describe the trajectories of circular orbits of spinless and spinning test particles around of rotating bodies in equatorial and non-equatorial planes via the Mathisson-Papapetrou-Dixon equations which include the Ricci rotation coefficients with the purpose of describing not only the curvature of space time, but also the rotation of the spinning test particles that orbit around the rotating masive bodies.


Introduction
For to study the physics of a rotational gravitational field, one of the ways is describe the motion of test particles in this field. Tipically in the literature, the description of the motion of spinning test particles in a rotating gravitational field in non-equatorial plane we finds two main approaches. The first one is given by Mathisson-Papapetrou-Dixon Equations (MPD) [1], [2], [3] which yields the equations of motion of a spinning test particle in a given gravitational field. The second takes the first integrals of motion given by Carter [4] and derives the orbits of spinless test particles around of rotating massive bodies. The term of particle is because the size of the test body is small compared to the scale of the curvature. The motion of particles on the metric type Kerr have been discussed by many authors [5], [6], [7], [8]. Many of these works study the motion in the equatorial plane for spinless and spinning test particles via MPD Equations [9], [10]. Using this last formalism, we calculate the trajectories of spinless and spinning test particles in non-equatorial planes. On the other hand, in the literature, there are studies that use the Carter s equations for both spinless and spinning test particles in equatorial planes [11], [12], but there are few works that describe the motion of spinless test particles around on rotating bodies in non-equatorial planes [13], [14]. This is the first novelty in relation to previous works. On the other hand, the study of trajectories of spinning test particles in non-equatorial planes for a Kerr spacetime under the Carter s equations needs further investigation. We focus our work in the formalism given by the MPD equations [2]. In this via, it is given a distribution of mass (m) with spin tensor (S ρσ ) around on central source (M ) which has a metric tensor g µν . These equations of motion for a spinning test particle are obtained in terms of an expansion that depends on the derivatives of the metric and the multipole moments of the energy-momentum tensor (T µν ) [2], and are given by where D/dτ means the covariant derivative, the vector p µ and the antisymmetric tensor S µβ are the linear and spin angular momenta respectively. R µ νρσ is the curvature tensor. Since p µ p µ = constant and S ρσ S ρσ = constant along the particle trajectory [15]; then, we may set and with these expressions, we obtain the center of mass condition and the relation between the spin tensor and the vector spin which is called spin supplementary condition (SSC). For this case, we used the Tulczyjew supplementary condition [16]: In the present paper, we calculate numerically, via MPD equations, the trajectories for spinless and spinning test particles in circular orbits when the test particles are in non-equatorial planes. We are interested in to study the effects from spin of the test particle orbiting in a rotating gravitational field which is generated by a massive rotating body. In the majority of cases, the description of trajectories is given by the study of the Christoffel symbols, the novelty of our work is to take the Ricci Rotation Coefficients (RRC) in orden to study the trajectories of spinning test particles in non-equatorial planes. On the other hand, when we study these trajectories, we can describe the space time around a rotating massive body and its phenomena.
The Carter´s equations are given by the first integrals of the equations of motion of a spinless test particle around on a rotating body for a Kerr metric both in an equatorial plane and in non-equatorial planes [4]. These equations use the symmetries of the Kerr geometry and the conserved quantities of energy (E), angular momentum (L), rest mass (M ) and a fourth integral of motion called Carter´s constant (Q). For the case of spinning test particles in the equatorial plane, when the space-time possess a Killing field, there exists a linear combination of the components of momenta (p µ ) and the spin tensor (S µν ) [17]. The study for orbits of spinning test particles off to equatorial plane by the Carter´s equations needs further investigation.
We shall compare the results by numerical integration in three cases: First, we compare between the trajectories of spinless test particles in non equatorial planes given by the Carter´s equations and the MPD equations. Second, the trajectories both of spinless and spinning test particles in non-equatorial planes given by the MPD equations. Third, we compare the behavior of two Boyer Linquist coordinates ( θ, ϕ) in regard to the proper time (τ ) both to spinless and to spinning test particles.
The paper is organized as follows. In Section 2, the equations of motion MPD are reduced to expressions that include the Ricci Rotation Coefficients (RRC). These equations describe the motion of spinless and spinning test particles around on a rotating body. We study the cases when the spinless and spinning test particles are in the equatorial plane and give the set of equations when the spinless and spinning test particles are out of the equatorial plane. In Section 3, there is a brief description of the Carter´s equations for spherical orbits. These equations take three constants of motion given by the symmetries of a rotating body. A fourth constant is obtained from the separability of the Hamilton-Jacobi equation. These orbits are calculated numerically. Then, in Section 4, we present a comparison numerical between the MPD equations and Carter´s equations for the spinless test particles both in the equatorial plane and the non-equatorial plane. We take the initial values from the Carter´s equations and are replaced in the MPD equations. In the last section, the conclusions and some futures works are formulated in order to describe spinning test particles in a Kerr type metric.

MPD equations and the Ricci Rotation Coefficients
The equations MPD had been traditionally used in the equatorial plane. We focus our study not only in the equatorial plane, but also in the non equatorial plane. In this paper, we take the work from Tanaka et al. [9] who reduce the equations of motion using the RRC. They solve the equations of motion for a circular orbit in the equatorial plane while we, with the same via, solve the case for trajectories non-equatorial plane. Given this, the equations of motion (1) and (2) for the tetrad components are reduced to where ω βγ α are the RRC and ζ α is the unit spin vector, which is defined by ζ α := S a /S . Here R α and S µν are defined by the last equation yields the relation between tensor spin (S γλ ) and the vector spin (S κ ). For our calculation, we have an orthonormal frame with the RRC and we can extract information about the curvature of the spacetime in question. It will be important when we describe the trajectories of test particles off-equatorial plane.
On the other hand, the RRC are related with the the coordinates of the derivative vectors that measure not only the deviation, but also the rotation of all the frame vector when moved in various directions. Furthermore, we can observe that the RRC carry information about how the frame rotates and not just information about the curvature of space time.

Calculation of the Ricci Rotation Coefficients
Next, it is to define the Ricci rotation coefficients for a Kerr spacetime [18] as ω ab c = e a µ e b ν e c ν; µ where the semicolon indicates the covariant derivative (Appendix A).
In order to find a solution, it is convenient to introduce the tetrad frame [9] defined by e i µ = e i t , e i r , e i θ , e i where the Greek letters distinguish the tensor indices (µ = t, r, θ, ϕ) from the tetrad indices which are in the Latin alphabet (i = 0, 1, 2, 3). Here it is used the metric signature (−, +, +, +). The RRC are given by [19], [20] where λ (a)(b)(c) are the rotation pre-coefficients and are defined by

Spinless test particles in a Kerr metric
In this section, the equations of motion for spinless test particles in a Kerr metric with radius constant (r = r 0 ) are solved. The aim is to write down the equations with the help of the Eq. (7) and the RRC (Appendix A). In this case, the magnitude of the spin is equal to zero (S = 0 ) and r is constant, i.e. u 1 = 0. Since the test particles is spinless and according to Eq. (9), it can be identified the velocity v µ with the velocity u µ in Eq. (7). Therefore the system of equations is given by Let be replaced the values of the RRC and the tetrads in Eqs. (15) - (18). These latter equations for spinless test particles will be solved by numerical integration with Mathematica.
For orbits in the equatorial plane, we assume that corresponds to: r = constant (u 1 = 0), and close to the equatorial plane (u 2 = O( ∼ θ)), that is, Under this assumption, the nontrivial equation is [9] d dτ u 1 = ω 00 In this case, the orbital angular velocity is given by where M is the mass of central body and a is the angular momentum of the source.

Spinning test particles in a Kerr metric
First of all, for the case of spinning test particles in circular orbits and in the equatorial plane, we assume that corresponds to: r = constant (u 1 = 0), close to the equatorial plane (u 2 = O( ∼ θ)), and with the magnitude of spin (S 2 = S µ S µ ). In this case, according to Eq. (7), the nontrivial equations of the orbital motion are [9] d dτ u 1 = ω 00 1 u 0 u 0 + ω 30 1 u 3 u 0 + ω 03 where the components of R a are given by In the equatorial plane, the orbital angular velocity for spinning test particles is given by this new expression corresponds to the spin contribution (S) at O( ∼ θ). Now for the orbits non-equatorial, first, we take the work by Tanaka et al. [9] which studies the motion of the spinning test particle in the equatorial plane, and then, we go one step further and numerically calculate the trajectories off-equatorial orbits. For this treatment, we replace the "dynamical velocity v" Eq (9) in Eqs. ( 7) and (8), and obtain the set of equations: For the numerical integration, we need the initial values of the four velocity which will be obtained from the Carter´s equations.

Equations of Motion and Carter´s constant
In a manifold type Kerr, the symmetries provide three constants of motion: the energy (E), the angular momentum (L), and the rest central mass (M ). In addition, there is another, one constant which is due by the separability of the Hamilton -Jacobi equation and is called Carter´s constant (Q). The equation of Lagrange for a Kerr metric leads immediately to the first integrals of the t and ϕ equations and for the other two integrals of (r) and (θ) are obtained by a separable solution of the Hamilton -Jacobi equation. The set of equations is given by [21] where L, E, and Q are constants and Σ : = r 2 + a 2 cos 2 θ, ∆ : = r 2 + a 2 − 2M r, M , and a = J/M c denote the central mass and specific angular momentum of the central source which gives rise to the gravitational field represented by the Kerr space-time.
In the study of orbits of test particles around on Kerr metric, there exists a kind of orbits called spherical, i.e., constant radius. This last orbit intersects the equatorial plane in a point called node. Since the metric has angular momentum, the nodes of spherical orbits are dragged in the sense of the spin of the rotating body. When there is a particle orbiting in a nonequatorial orbit, this traces a kind of helix until a maximum of latitude and when reachs the maximum this begins to descend until a minimum latitude which is symmetric to the maximum [22].
When the space-time admits a Killing vector ξ υ , there exists a property that includes the covariant derivative and the spin tensor, which gives a constant and is given by the expression [17] where p ν is the linear momentum, ξ ν,µ is the covariant derivative of Killing vector, and S νµ is the spin tensor of the particle. In the case of Kerr metric, there are two Killing vectors, owing to its stationary and axisymmetric nature. As consequence, Eq. (30) yields two constants of motion: E is the total energy and L is the component of its angular momentum along the axis of symmetry. Unfortunately, there is not a constant of motion associated to the coordinate r; therefore, for this way, there does not exist a equation of motion that connects the constant of motion linked to coordinate r and the spin tensor. This mathematical expression needs further investigation. On the other hand, people avoid the error that shows the numerical integration due to the explicit square roots in the r and θ in Eqs.(27) and (28). In this formulation, they take the Hamilton´s canonical equations which denote differentiation of the Hamiltonian with respect to r and θ. This demostration is based on a taxonomy of all periodic orbits around black holes [23]. The set of equations constitutes a smooth system of ordinary differential equations and can be integrated without change of variables. These equations were integrated numerically to generate the periodic orbits. These orbits describe the motion of test particles out of equatorial plane, but this approach does not take into account the spinning test particles. So far, there is not a work that connect successfully the equations of motion given by Carter and spin of the test particle. Within Carter´s equations, one finds in the literature a number of works, but none considers the spinning test particle out of the equatorial plane [24].

Numerical integration of the Equations of motion
In this section of this paper, we worked the MPD Equation not so much with the Christoffel symbols, but with the RRC, which describe the trajectories of spinning test particles around on rotating massive bodies. Now, we shall compare numerically the trajectories of the test particles both the spinless and spinning. For this comparison, first we need to compare between the trajectories of the spinless test particles given by the Carter´s equations with the trajectories given by MPD equations. In this case, we take as parameters: radius constant r = 10, E = 0.9525, a = 1, Q = 4.224806, L = 2.810974, maximum latitude = 53.92928 • and magnitude of spin equal to zero. We replace these values in the Equations of Carter (26) -(29) and find the initial values of the four velocity: dt 0 /dτ = 1.00748, dθ 0 /dτ = 2.854 × 10 −3 , and dϕ 0 /dτ = 8.255 × 10 −1 . Then, these initial values of the four velocity are replaced in the set of MPD equations for spinless test particles (15), (17) and (18). We plot in 3D the trajectory for a spinless test particle in Mathematica ( Figure 1). Finally, we compare the data table of the cartesian coordinates both the Carter´s equations and MPD equations, and find that the difference of values is in the order of 10 −7 . Next step, we compare between the trajectory of a spinless test particle ( 15,17,18), with the trajectory of a spinning test particle (24,25) in nonequatorial planes via MPD equations. We replace the values of the RRC and the components nonvanishing of curvature tensor given by Mino [18]. Then, we plot in 3D the trajectories both for a spinless test particle and a spinning test particle in Mathematica (Figure 2). And finally, we compare the data table of the cartesian coordinates both for a spinless test particle and a spinning test particle, and find that the difference of values is in the order of 10 −4 . This small difference in the trajectories is given by the effect of spin in a test particle [25]. We make another numerical comparison of our results given by the MPD equations in regard to the trajectories of spinless and spinning test particles. In this case, we take the same initial values for the four velocity given above which were obtained by the Carter´s equations (26) -(29), and replace these initial values in the MPD equations (24,25). Then, we take two Boyer Lindquist coordinates (θ, ϕ) and obtain a graph for each coordinate in regard to the proper time (τ ) which show the difference in the trajectories between the spinless and spinning test particles. The first one is the figure of the polar angle versus proper time (θ vs τ ) (Figure 3). In this graph, we find the period of the spinning test particle is longer than the period of the spinless test particle. Moreover, in each orbit, the difference between the trajectories of each test particle increases, that is, it is an accumulative process. This is due to the contribution of the value of spin in the trajectory. The test particles starts in the polar initial angle (θ 0 = 53.92928 • ) and reaches the maximum value in the opposite side of the equatorial plane (θ op = 143.2928 • ). One can see in this graph that both the spinless and spinning test particle reaches the same maximum o minimum of latitude, but with a time difference for the two paths. This is due to the precession of the spinning test particle. The second figure shows the relation between the azimuthal angle and the proper time ( Figure  4). We find the slope of the graph for the spinning test particle is winding with relation to the graph of the spinless particle. This phenomenum is another effect of the spin of the particle [26] and like the previous graph this winding is due to the nutation of the spinning particle. The MPD equations described with the help of the RRC allow characterizing these phenomena.

Conclusions
We found that the equations of motion for spinning test particles in the MPD equations (7 -9) describe the relation between the spin of the particle and the angular momentum of the central mass. For this case, we take the variation of spin vector in the time and the first term includes as much from the density of angular momentum (a) as from the spin of the test particle (S). For the case of spin -orbit perturbation, the equations yields the relation between the spin of the particle given by the magnitude of spin (S) and the orbit given to both the curvature tensor (R α βγλ ) and the velocity of center of mass (v µ ). For circular orbits close to equatorial plane, the difference of trajectories between spinless test particles and spininng test particles is very small; but when the particles have spin the difference in the time is bigger that when the particles does not have it.
The majority of works that include the Christoffel symbols in the MPD equations does a description of the coupling between spin and field while in this paper, since that we include the RRC in the MPD equations, we study not only this coupling, but also the rotation of the frame that travels with the test particle. These phenomena are studied with graphs both of the polar angle and the azimuth angle.
Additionally the results obtained in this paper may be important guideline to quantitative estimates of gravitational waves in numerical relativistic simulations. When there is a spinning test particle in a gravitational field the description of its geodesics is given, in the same time, the characteristics of the space-time for where is traveling the particle. Now, if the gravitational waves are detected for instruments as the Laser Interferometric Gravitational Wave Observatory (LIGO) or the Laser Interferometer Space Antenna (LISA), it is possible comparating the description of the orbits of the spinning test particles with the signal of these instruments and giving some conclusions. In the majority of these researches, they use techniques in interferometry. Our results are good material when it comes to face the construction of a gravitational wave detector that takes account the effect of spin in the template of interferometry. The description of this phenomenom can be an extension of the work done in this paper with help of the MPD equations.