On the (1+3) threading of spacetime with respect to an arbitrary timelike vector field

We develop a new approach on the (1+3) threading of spacetime $(M, g)$ with respect to a congruence of curves defined by an arbitrary timelike vector field. The study is based on spatial tensor fields and on the Riemannian spatial connection $\nabla^{\star}$, which behave as $3D$ geometric objects. We obtain new formulas for local components of the Ricci tensor field of $(M, g)$ with respect to the threading frame field, in terms of the Ricci tensor field of $\nabla^{\star}$ and of kinematic quantities. Also, new expressions for time covariant derivatives of kinematic quantities are stated. In particular, a new form of Raychaudhuri's equation enables us to prove Lemma 6.2, which completes a well known lemma used in the proof of Penrose-Hawking singularity theorems.Finally, we apply the new $(1+3)$ formalism to the study of the dynamics of a Kerr-Newman black hole.


Introduction
The (1+3) threading of spacetime by a congruence of curves determined by a unit timelike vector field ξ (4-velocity) is by now a well established theory which studies the geometry, dynamics, and observational properties of some well known cosmological models. Most of the important results on this theory, and an exhaustive list of references wherein these results have been published, can be found in the excellent monograph of Ellis, Maartens and MacCallum [1].
Our work on this matter is motivated by the simple remark that it is difficult to apply the above theory to the metrics of general form presented in (2.8). This is because, in this case, ξ = ∂/∂x 0 is not a unit vector field and thus it should be normalized. But this process leads to complicated formulas for kinematic quantities and Ricci tensor field, which of course makes difficult their study. The question is: Are there important cosmological models whose metrics have the general form (2.8). The answer is in the affirmative and it is based on the following two examples. First, the study of the cosmological perturbations of the FLRW universes is developed with respect to the metric (cf. (10.12) of [1]) (1.1) Then, we recall that the metric of a Kerr-Newman black hole is given by (cf.(12.3.1 in [2]) where we put ∆ = (x 1 ) 2 + a 2 + e 2 − 2mx 1 , Σ = (x 1 ) 2 + a 2 (cos x 2 ) 2 .
The metric (1.1) was intensively studied with respect to the (1+3) threading of almost FLRW universes. From Chapters 10 and 11 of [1] we can see that the study is not an easy one in literature. Also, as far we know, very little has been done with respect to the (1+3) threading theory for the metric (1.2) (cf. [2,3]).
In this paper we present a new approach on the (1+3) threading of spacetime with respect to a congruence of curves defined by an arbitrary timelike vector field ξ. We develop a method that is based on the follwing concepts: (i) Threading frame and coframe fields.
(iii) Riemannian spatial connection. The threading frame and coframe fields are naturally constructed from the coordinate fields (cf. (2.3) and (2.4)), and have a great role throughout the paper. The spatial tensor fields have been used in earlier literature, but in here we work only with their 3D local components with respect to the above special frames (cf. (3.1)). This brings a substantial simplification into the study of such general metrics. Finally, the Riemannian spatial connection ∇ ⋆ (cf. (3.12)) is a metric linear connection on the spatial distribution, which introduces both the spatial and time covariant derivatives. It is important to note that throughout the paper, all geometric objects and equations involved into the study, are expressed in terms of spatial tensor fields and their spatial or time covariant derivatives. As the metrics (1.1) and (1.2) fall into the class of the general metrics given by (2.8), the (1+3) threading theory developed here can be easily applied to their study. Now, we outline the content of the paper. In Section 2 we consider the orthogonal decomposition (2.1) of the tangent bundle of the spacetime (M, g), and construct the threading frame and coframe fields {∂/∂x 0 , δ/δx i } and {δx 0 , dx i }, respectively (cf. (2.3) and (2.4)). Also, we consider the Riemannian metric h on the spatial distribution SM given by its 3D local components (cf. (2.15)). Then, in Section 3 we introduce the notion of spatial tensor field via 3D local components (cf. (3.1)), and show that the vorticity, expansion and shear tensor fields ω ij , θ ij and σ ij given by (3.5a) and (3.8), are indeed spatial tensor fields. Also, we define the Riemannian spatial connection ∇ ⋆ on SM (cf. (3.11)) and express the Levi-Civita connection ∇ in terms of the local coefficients of ∇ ⋆ and the above kinematic quantities (cf. (3.17)). A comparison between the concepts defined in this paper and the corresponding ones from earlier literature is done in Section 4. In particular, for a unit timelike vector field we obtain (4.14) for kinematic quantities, and deduce that they do not depend on the Levi-Civita connection of the spacetime. In Section 5 we express both the curvature tensor field and the Ricci tensor field of (M, g) by spatial tensor fields and their spatial and time covariant derivatives (cf. (5.3), (5.5a), (5.11), (5.12a)). Next, in Section 6 we obtain the Raychaudhuri's equation (6.1) with respect to an arbitrary timelike vector field, which for a congruence of timelike geodesics takes the forms (6.5) or (6.28). It is important to note that (6.28) is the main ingredient used in the proof of Lemma 6.2, which should be considered as a completion of Lemma 6.1 that has been the key in the proof of Penrose-Hawking singularity theorems. Also, we express the non-zero local components of the electric Weyl curvature tensor field in terms of spatial tensor fields (cf. (6.20)), and deduce new formulas for time covariant derivatives of the kinematic quantities (cf.(6.9), (6.12), (6.24), (6.25)). Finally, the last three sections are devoted to the study of a Kerr-Newman black hole via the new approach on the (1 + 3) threading of spacetime developed in the paper. In particular, we characterize spatial geodesics and obtain the 3D force identity (cf.(9.10)).

Threading Frame and Coframe Fields
Let (M, g) be a 4D spacetime, and ξ be a timelike vector field that is globally defined on M . Note that ξ is not necessarily a unit timelike vector field, as it was considered in early literature. The timelike congruence determined by ξ is tangent to the fibres of the line bundle V M , that we call the time distribution. Also, we consider the spatial distribution SM , which is complementary orthogonal to V M in T M , that is, we have (2.1) Throughout the paper we use the ranges of indices: i, j, k, ... ∈ {1, 2, 3} and a, b, c, ... ∈ {0, 1, 2, 3}. Also, for any vector bundle E over M denote by Γ(E) the F(M )-module of smooth sections of E, where F(M ) is the algebra of smooth functions on M .
The foliation by curves that is tangent to V M , induces a special coordinate system (x a ) such that ξ = ∂/∂x 0 . If (x a ) is another coordinate system, then we have since ∂/∂x 0 and ∂/∂x 0 represent the same vector field ξ, and hence ∂x 0 /∂x 0 = 1. Then, from (2.1) we deduce that for each ∂/∂x i there exist a unique δ/δx i ∈ Γ(SM ) and a unique function A i , such that This enables us to consider the threading frame field {∂/∂x 0 , δ/δx i }, and the threading coframe field {δx 0 , dx i }, where we put Note that {δ/δx i } are transformed exactly as {∂/∂x i } on a 3D manifold, while {A i }, in general, do not satisfy some 3D tensorial transformations.
Next, we consider the 1-form ξ ⋆ given by The local components of ξ ⋆ with respect to the natural frame field {∂/∂x i } are given by where Φ is a non-zero function on M which is independent of x 0 . The above condition on Φ is not restrictive for our theory, because most of the important cosmological models satisfy it. According to (2.7), the line element of g is expressed as follows: where we put Taking into account that and therefore Multiply (2.5c) by Φ 2 , and using (2.11), we obtain Hence ξ i do not define a 3D 1-form on M . Now, denote by h the Riemannian metric induced by g on SM , and put (2.14) Then by using (2.14), (2.3), (2.11), (2.9) and (2.7), we infer that Thus ds 2 from (2.8) is expressed in terms of the threading coframe field {δx 0 , dx i } as follows, Note that h ij and the entries h ij of the inverse of the matrix [h ij ] are transformed exactly like 3D tensor fields, that is, we have Kinematic Quantities as Spatial Tensor Fields on (M, g) The purpose of this section is to define the vorticity tensor field, expansion tensor field, expansion scalar and shear tensor field, as spatial tensor fields on the spacetime (M, g). First, we give the following definition. A spatial tensor field T of type (p, q) on M , is locally given by 3 p+q locally defined functions T with respect to the transformations (2.2). In other words, the local components of a spatial tensor field on M should satisfy the same transformations as the local components of a tensor field on a 3-dimensional manifold. From (2.17) we see that h ij (resp. h ij ) define a spatial tensor field of type (0,2) (resp.(2,0)) on M . By using (2.5a) and taking into account that we deduce that define a spatial tensor field of type (0,1). Next, by direct calculations using (2.3), (2.11) and (3.3), we deduce that where we put (3.5) Now, apply δ/δx j and ∂/∂x 0 to (2.5c) and by using (2.5a) and (2.5b), we infer that Then, by using (3.6) into (3.5), we obtain Hence ω ij and a i define spatial tensor fields of type (0.2) and (0,1), respectively. We call ω = (ω ij ) the vorticity tensor field for the timelike congruence defined by ξ on M . From (3.4a) we see that the spatial distribution SM is integrable if and only if ω ij vanish identically on M . Next, we define Then, take derivatives with respect to x 0 in (2.17a), and obtain is, Θ ij define a spatial tensor field of type (0,2). We call it the expansion tensor field. Clearly Θ from (3.8b) is a function, and σ ij from (3.8c) define a trace-free spatial tensor field of type (0,2). We call Θ the expansion scalar and σ ij the shear tensor filed for the congruence. According to the terminology from earlier literature, we call {ω ij , Θ ij , Θ, σ ij } given by (3.5a) and (3.8), the kinematic quantities with respect to the congruence of curves defined by the timelike vector field ξ = ∂/∂x 0 .
Raising and lowering latin indices is done by using h ij and h ij , as follows In order to define covariant derivatives of the above kinematic quantities, we consider the Levi-Civita connection ∇ on (M, g) given by (cf. [4], p.61) (3.10) for all X, Y, Z ∈ Γ(T M ). Then define the linear connection ∇ ⋆ on the spatial distribution as the spatial projection of ∇ on SM , that is, we have where s is the projection morphism of T M on SM with respect to (2.1). Note that ∇ ⋆ is a metric linear connection on SM . We call it the Riemannian spatial connection. Remark 3.1 The Riemannian spatial connection ∇ ⋆ is different from the three-dimensional operator∇ that has been used in earlier literature (cf. (4.19) of [1]). ∇ ⋆ is a linear connection on SM and therefore defines covariant derivatives of any spatial tensor field with respect to vector fields on M . On the contrary,∇ is an operator which acts on tensor fields on M , but in general, does not define a linear connection on M . ✷ Locally, we put Then, take X = δ/δx j , Y = δ/δx i and Z = δ/δx h in (3.10) and by using (3.11), (3.12a), (2.14) and (3.4a), we obtain Similarly, we deduce that Now, consider a spatial tensor field T of type (p, q).
T are spatial tensor of type (p, q + 1) and (p, q), respectively. As an example, we consider T = (T i j ), and obtain We call (3.15a) (resp. (3.15b)) the spatial (resp. time) covariant derivative of T . As ∇ ⋆ is a metric connection on SM , we have Finally, by using (3.10), the above spatial tensor fields and the local coefficients of ∇ ⋆ , we express the Levi-Civita connection ∇ on (M, g) as follows (3.17)

Comparison with Concepts from Earlier Literature
In the previous section we introduced the kinematic quantities on a spacetime (M, g) with respect to the congruence that is tangent to an arbitrary timelike vector field ξ. If in particular, ξ is a unit timelike vector field, the configuration of the spacetime with respect to the congruence of timelike curves determined by ξ is known in literature as (1+3) threading of spacetime (cf. [1,5]). In this section we show that for Φ 2 = 1 in (3.5) and (3.8) (that is, ξ is a unit vector field), we obtain the well known kinematic quantities from earlier literature. First, by using (2.6) and taking into account that ∇ is a metric connection, we obtain Then, consider the threading frame {∂/∂x 0 , δ/δx i } and using (3.17), we infer that Next, taking into account (2.12) and (2.7b), we express the natural frame field as follows Then, consider the covariant acceleration vector field and using (4.2b) and (4.2d), we obtain where we put Thus, we deduce that the congruence defined by ξ is a congruence of timelike geodesics, if and only if, we have For this reason we call b i the geodesic spatial covector field of the congruence. Now, by direct calculations, using (4.3) and (4.2), we infer that Taking into account (4.3), (4.8) and (4.5), we find Now, we suppose that ξ is a unit vector field. Then, according to the formula (4.38) from [1], p.85, we have where σ ab , θ, h ab and ω ab are quantities defined in earlier literature. On the other hand, in this case we have Φ 2 = 1, and from (4.9) we obtain since c j = 0, for all j ∈ {1, 2, 3}.
Comparing the symmetric and skew-symmetric parts in (4.10) and (4.11) we deduce that (4.12) Finally, comparing (4.12a) with (4.31) from [1], p.81, we obtain According to (3.8), (3.5a), (4.12) and (4.13), we conclude that in case ξ is a unit vector field, the only possible non-zero local components of expansion, shear and vorticity tensor fields from earlier literature are the following (4.14) As far as we know, the formulas (4.14) do not appear in earlier literature. Due to them we can state that the expansion, shear and vorticity tensor fields do not depend on the Levi-Civita connection of the spacetime (M, g). Of course, due to (3.5a) and (3.8), this conclusion is still valid for the general case of a congruence defined by an arbitrary timelike vector field ξ.
We close this section with an interesting property of vorticity tensor field. Suppose that we have a congruence of geodesics defined by ξ. Then, according to (4.7), (4.6), (3.3) and (3.5b), we have (4.15) By using (4.15), (3.2) and (3.3), we infer that (a) Then, take the derivative with respect to x 0 in (3.5a) and by (3.2), (3.3) and (4.16), we obtain Thus, we can state that the vorticity tensor field for a congruence of timelike geodesics of a spacetime, is independent of time. In particular, this is true for a congruence of geodesics with respect to a unit vector field ξ. However, we did not see this result in earlier literature. This is because the formulas (3.5a), (3.8) and (4.14) we deduced for the kinematic quantities are much simpler than the ones by means of Levi-Civita connection.

Curvature and Ricci Tensor Fields of a Spacetime via Spatial Tensor Fields
In this section we show that the curvature tensor field of (M, g) is completely determined by three spatial tensor fields R ijkh , R i0kh and R i0k0 (cf. (5.3), (5.5a)). A similar result we obtain for the Ricci tensor of (M, g) (cf. (5.11), (5.12a)). Note that all these spatial tensor fields are expressed in terms of the curvature and Ricci tensor fields of the of the Riemannian spatial connection, and of all kinematic quantities introduced in Section 3.
In what it follows, R denotes both the curvature tensor field of (M, g) of type (0,4) and (1,3), given by for all X, Y, Z, U ∈ Γ(T M ). Then the curvature tensor field of (M, g) is completely determined by its local components By direct calculations, using (5.2), (5.1), (3.17), (3.4), (3.3) and (4.6), we obtain where R ⋆ ijkh are the local components of the curvature tensor field of the Riemannian spatial connection defined as in (5.2a), and given by Taking the symmetric and skew-symmetric parts, in (5.3c) we deduce that Then, we deduce that According to [4], p.87, the Ricci tensor of (M, g) is given by Then, by using (5.8), (5.6), (5.7) and (5.2), we obtain where we put (5.10) By using (5.3a), (5.3b) and (5.5a) into (5.9), we deduce that where Θ k i| k , ω k i | k and b k | k are spatial divergences given by formulas deduced from (3.15a). Now, take symmetric and skew-symmetric parts in (5.11a) and obtain where we put We call R ⋆ ik the spatial Ricci tensor of the spacetime (M, g). From (5.12b) we see that if the spatial distribution is integrable, then we have (5.14) In this case, we have Also, note that because in this particular case the vorticity vanishes identically, from (5.3a), (5.3b), (5.5a), (5.11a), (5.11b) and (5.12a), we deduce that the curvature and Ricci tensors on (M, g) are expressed as follows and where R ⋆ ik is given by (5.15).

Raychaudhuri's Equation and Time Covariant Derivatives of Kinematic Quantities
First, by using (3.8b) and (3.8c) into (5.11c), we infer where we put In particular, if ξ is a unit timelike vector field, that is which is Raychaudhuri's equation expressed in terms of local components of spatial tensor fields introduced in the present paper. Thus, we are entitled to call (6.1) the generalized Raychaudhuri's equation with respect to a congruence defined by an arbitrary timelike vector field ξ.
According to (4.7), in case of a timelike congruence of geodesics, (6.1) and (6.4) become and respectively. Remark 6.1. Formally, (6.6) looks like (9.2.11) in [2], but we should note that ω, σ and Θ from (6.6) are calculated via their 3D spatial components (see (3.5a), (3.8) and (6.2)), while in [2] they are calculated in terms of the 4D local components with respect to the natural frame field {∂/∂x a }. ✷ Next, observe that (5.5b) gives a formula for the time covariant derivative of vorticity tensor field. By using (3.8b) and (3.8c) into (5.5b) we find Now, from (5.5a) we deduce that the time covariant derivative of the expansion tensor field is given by Another formula in terms of Ricci tensors is deduced from (5.12a): Taking time covariant derivative in (3.8c), and using (3.16c), (6.8) and (6.1), we infer that where R i0k0 is the trace-free part of the spatial tensor field R i0k0 , given by In a similar way, but using (6.9) instead of (6.8), we obtain (6.12) Now consider the Weyl tensor field in (M, g) given by , (6.14) and R is the scalar curvature of (M, g). Then, we consider the electric Weyl curvature tensor field E = (E ac ) given by (6.15) and taking into account that ξ = ∂/∂x 0 , we obtain (6.16) Then by direct calculations, using (6.16), (6.13), (2.7), (2.9) and (2.15), we deduce that the only possible non-zero local components of E with respect to the natural frame field are Note that due to (2.3) and (6.16), we have By using (2.5a) , from (6.18) we infer that E ik define a spatial tensor field of type (0,2). Using (2.12) in (6.14b) and (6.14c), we obtain (6.19) where R i0k0 and {R ik , R i0 , R 00 } are given by (5.2c) and (5.10), respectively. Taking into account of (6.19) into (6.17) and using (2.15), we express E ik in terms of spatial tensor fields, as follows The scalar curvature R of (M, g) is given by Then, replace R from (6.21) into (6.20) and taking into account (6.11), we deduce that (6.22) where R ik is the trace-free part of R ik , that is, we have Finally, by using (6.11) and (6.22) into (6.8) and (6.10), we obtain (6.24) and (6.25) respectively. It is interesting to note that the generalized Raychaudhuri's equation (6.1) can be expressed by using the scalar curvature R of (M, g) and the spatial scalar curvature R ⋆ of ∇ ⋆ given by Indeed, contracting (6.9) by h ik and using (3.16d), (3.8b), (6.21) and (6.26) we deduce that (6.27) In particular, if ξ is a unit vector field that defines a timelike congruence of geodesics (see (6.3) and (4.7)), then (6.27) becomes This is a new form of the Raychaudhuri's equation (6.6) for a congruence of timelike geodesics. It is well known that (6.6) is the key equation used in the proof of Penrose-Hawking singularity theorems. More precisely, it was proved the following Lemma. Lemma 6.1. (see Lemma 9.2.1 in [2]) Let ξ be the tangent field of a hypersurface orthogonal timelike geodesic congruence. Suppose the following conditions are satisfied: (i) Ric(ξ, ξ) ≥ 0, which is the case if Einstein's equations hold in the spacetime and the strong energy condition is satisfied by the matter.
(ii) The expansion Θ takes the negative value Θ 0 at a point on a geodesic in the congruence corresponding to the proper time τ = 0.
Then Θ goes to −∞ along that geodesic within the proper time τ ≤ 3 |Θ 0 | . Note that in the above lemma, ξ is tangent field of a hypersurface orthogonal timelike geodesic congruence means that SM is an integrable distribution. Now, by using the new form (6.28) of Raychaudhuri's equation we can complete Lemma 6.1 with the following Lemma.

Kinematic Quantities for Kerr-Newman Black Holes
The new point of view developed here on the (1 + 3) threading of spacetime, is applied in this section to the charged Kerr black hole (also called Kerr-Newman black hole). We show that the curvature and Ricci tensor fields of (M, g) are simply expressed in terms of curvature and Ricci tensor fields of the Riemannian spatial connection, via the kinematic quantities. Now, according to the notations used in Sections 2 and 3, for the metric of a Kerr-Newman black hole given by (1.2), we have The spatial distribution SM of (M, g) is locally spanned by and it is the kernel of the 1-form By using (2.15) (1.2), (7.1a) and (7.1b), we deduce that the only non-zero local components of the Riemannian metric h on SM with respect to the threading frame field from (7.2), are the following Hence the line element from (1.2) becomes with respect to the threading coframe field {δx 0 , dx i }.
Next, by using (3.3), (4.6), (7.1a), and (7.1c), we deduce that the geodesic spatial tensor field b = (b i ) is given by Due to (4.7) and (7.6) we conclude that the curves from the congruence defined by ξ = ∂/∂x 0 which sit in the surface given by the equations for a = 0, or in the hypersurface for a = 0, are the only geodesics of (M, g) that are tangent to ξ. Moreover, we see that such geodesics have the equations for a = 0, and for a = 0, where k and c are constants. In particular, the integral curves of ξ can not be geodesics in the Schwarzschild spacetime. Now, taking into account (3.8) and (7.4), we obtain Also, by using (3.5a), (7.1b), (7.2) and (7.6), we deduce that the only nonzero local components of the vorticity tensor field (ω ij ), are given by (7.8) From (7.8) we see that the spatial distribution SM is not integrable for both Kerr-Newman and Kerr black holes. On the contrary, for the Reissner-Nordstrom and Schwarzschild solutions, the timelike vector field ξ = ∂/∂x 0 is hypersurface orthogonal.
Finally, we note that the local components of the curvature and Ricci tensor fields with respect to the threading frame field have very simple expressions. Indeed, by using (7.6) and (7.7) into (5.3a), (5.3b) and (5.5a), we obtain (7.9) where R ⋆ ijkh is the curvature tensor field of the Riemannian spatial connection. In a similar way, from (5.12a), (5.11b) and (5.11c), we infer that the local components of the Ricci tensor of (M, g) with respect to the threading frame field are given by where R ⋆ ik is the Ricci tensor of the Riemannian spatial connection. As far as we know, (7.9) and (7.10) have been not stated in earlier literature. They can bring more information and ideas in the study of the geometry and physics of the black holes. In particular, for Reissner-Nordstrom and Schwarzschild solutions, (7.9) and (7.10) become and (7.12) respectively.

Equations of Motion in a Kerr Black Hole
In this section and in the next one, we take e = 0 in (1.2), that is, we consider the Kerr black hole (M, g). As it is well known, the geodesic equations in a Kerr black hole have been explicitly integrated for first time by Carter [6].
In the present section we will state a new form of the equations of motion in a Kerr black hole, and obtain informations about the position of geodesics in M with respect to the spatial distribution. In particular, we show that the geodesics of (M, g) which are tangent to SM coincide with the autoparallel curves of the Riemannian spatial connection. Let C be a smooth curve in M given by parametric equations where λ does not necessarily represents the time in (M, g). The velocity vector field d/dλ for C is expressed in terms of the threading frame {∂/∂x 0 , δ/δx i } as follows: where we put Now, by using (7.7) and (7.1c) in (3.17), we express the Levi-Civita connection ∇ on (M, g), as follows: Then, by direct calculations using (8.4) and (8.2), we obtain which leads to the following equations of motion Note that λ from (8.5) is an affine parameter for geodesics in (M, g). A geodesic of (M, g) which is tangent at any of its points to the spatial distribution SM , is called a spatial geodesic. Then by using (8.2), (8.3) and (8.5) we deduce that a curve C given by (8.1) is a spatial geodesic, if and only if, it is a solution of the system Now, we remark that (7.7a) implies that the Kerr spacetime has bundle-like metric with respect to the foliation determined by ξ (cf. [7], p.112). Thus, we have the following interesting property: If a geodesic of a Kerr black hole is tangent to the spatial distribution at one point, then it remains tangent to it at all later times. Also, due to (8.6) we may state the following: The spatial geodesics in (M, g) coincide with autoparallel curves for the Riemannian spatial connection.
Next, we suppose that C is a geodesic in (M, g) which is not spatial, that is, we have δx 0 /δλ = 0. Without loss of generality we suppose δx 0 /δλ > 0. Then by using (3.3) in (8.5b), we obtain which is equivalent to where K is a positive constant. By using (8.7) into (8.5a), we deduce that a geodesic of a Kerr black hole (which is not a spatial geodesic), must be a solution of the system formed by (8.7) and the equations Finally, taking into account (8.6b) and (8.7), we conclude that the system of differential equations for the geodesics in a Kerr black hole can be arranged in such a way that one of the equations is of first order.

A 3D Identity Along a Geodesic in a Kerr Black Hole
The Riemannian spatial connection given by (3.11), enables us to define a 3D force in a Kerr black hole, and to deduce what we call the 3D force identity (cf.(9.10)). Note that this 3D force is a direct consequence of the existence of the fourth dimension (time), and emphasizes an important difference between Newtonian gravity and Einstein's general relativity. Let C be a geodesic in (M, g), and U (λ) be the projection of velocity d/dλ on SM . Then by (8.2) we obtain which we call the 3D velocity along C. Now, consider the 3D arc-length parameter s ⋆ on C given by and obtain is a unit spatial vector field, that is, we have Since ds ⋆ /dλ is positive, we can take s ⋆ as a new parameter on C. Then, define the 3D force along C, as the spatial vector field F (s ⋆ ) given by (9.5) Here, ∇ ⋆ is the Riemannian spatial connection, and d/ds ⋆ is given by (9.6) Taking into account that ∇ ⋆ is a metric connection on SM , and using (9.4) and (9.5), we deduce that F (s ⋆ ) is orthogonal to both U (s ⋆ ) and U (λ). By (9.5), (9.6) and (8.6), we see that the 3D force vanishes along a spatial geodesic C, if and only if, s ⋆ is an affine parameter on C.
Next, we put F (s ⋆ ) = F k (s ⋆ ) δ δx k , and by using (9.5), (9.6), (3.12), (3.14) and (7.7a), obtain provided C is not a spatial geodesic. Now, by using (8.7) and taking into account that , from (9.7) we deduce that the local components of the 3D force with respect to the affine parameter λ , are given by (9.8) Finally, by using (8.8) into (9.8), we infer that Taking into account that h kh F k (λ) dx h dλ = 0, and using (9.9), (3.3) and (9.2), we obtain Note that the identity (9.10) is a direct consequence of the existence of the 3D force F given by (9.5). For this reason we call it the 3D force identity.
Such an identity could be useful in a study of motions in (M, g), and even for solving the equations of motion. For example, from (9.10) we deduce that: The 3D arc length parameter s ⋆ is an affine parameter on the geodesic C, if and only if, Φ is constant along C.

Conclusions
In the present paper we develop a theory for a (1+3) threading of spacetime (M, g) with respect to a congruence of curves determined by an arbitrary timelike vector field ξ = ∂/∂x 0 . The generality of the study is not the only difference between our approach and what is known in literature for the case of the unit vector field ξ. The main differences consist in the following: (i) Our approach is entirely developed with geometric objects expressed by their local components with respect to the threading frames {∂/∂x 0 , δ/δx i } and threading coframes {δx 0 , dx i }.
(ii) The spatial distribution is not supposed to be necessarily integrable, and therefore this theory can be easily applied to the study of any cosmological model with non-zero vorticity.
(iii) All the equations and results we state, are expressed in terms of the spatial tensor fields (cf. (3.1)) and their spatial and time covariant derivatives (cf. (3.15)) induced by the Riemannian spatial connection given by (3.12).
(v) The proof of Lemma 6.2 which completes the well known Lemma 6.1, is entirely based on a new form of Raychaudhuri's equation for a congruence of timelike geodesics (cf.(6.28)).
(vi) It is the first time in literature when the spatial geodesics of a Kerr black hole are investigated (see (8.6) and the assertions which follow it).
(vii) The 3D force (9.5) and the 3D force identity (9.10) are new objects in the general theory of Kerr back holes, and illustrate the differences between Newtonian gravity and Einstein's general relativity.