Space-time singularities in Weyl manifolds

We extend one of the Hawking–Penrose singularity theorems in general relativity to the case of some scalar-tensor gravity theories in which the scalar field has a geometrical character and space-time has the mathematical structure of a Weyl integrable space-time. We adopt an invariant formalism, so that the extended version of the theorem does not depend on a particular frame.


Introduction
Until the mid-1960s it was argued by some cosmologists that the presence of space-time singularities in general relativistic cosmological models was not an essential property of the model, being, in fact, a consequence of the high degree of symmetries of the distribution of matter assumed in these models [1]. Accordingly, it was believed that in more realistic situations these singularities would disappear. However, this scenario was to change drastically after a series of general mathematical results, namely, the so-called Hawking-Penrose singularity theorems, were proved [2][3][4][5]. As is widely known, these theorems use methods of global analysis to show that, under the assumption of the validity of general relativity and a reasonable physical behavior of matter, space-time singularities are general phenomena which occur in gravitational collapse and cosmology (such as the big bang) irrespective of the symmetry of the models.
The investigation of space-time singularities has touched on some deep philosophical issues and it is the view taken by many scientists that these issues seem to call for a conceptual revision of general relativity that at least take quantum mechanics into account. As far as cosmological singularities are concerned it seems that the current view held by most cosmologists is that general relativity must break down at a e-mail: iarley.pereiralobo@icranet.org b e-mail: adrianobraga@fisica.ufpb.br c e-mail: cromero@fisica.ufpb.br times less than the Planck time. In fact, there has been a fair deal of work on classical and quantum new approaches to gravity in which mechanisms that naturally prevent the formation of a cosmological singularity are present. In fact, from the standpoint of quantum mechanics the mere existence of a cosmological singularity would be in contradiction with Heisenberg's uncertainty principle. In connection with this point, many efforts have been made to understand what could be called the quantum structure of space-time. Among the numerous attempts to make progress in this rather difficult issue, we would like to call attention for a recent proposal which argues that the quantum structure of space-time may even be related to the possibility of the space-time geometry possessing a non-Riemannian character [6]. On the other hand, non-Riemannian geometries have appeared in physics mainly as a way to modify Einstein's gravity while keeping the idea that the gravitational field is a manifestation of the space-time geometry [7,8]. As is well known, one of the earlier attempts in this direction was the Weyl unified field theory [9][10][11][12][13][14][15][16].
Concerning space-time singularities, we know that modified theories of gravity call for a new mathematical treatment of the problem, as the approach provided by Hawking and Penrose is suitable only for general relativity. It is our purpose here to pursue this question further in the light of some modified gravity theories formulated in a particular kind of space-time geometry, namely, the so-called Weyl integrable space-time (WIST) [17][18][19][20][21][22][23][24][25][26][27][28][29].
The present article is organized as follows. In Sect. 2, we outline the fundamental ideas of the geometry developed by Weyl, which underlies his unified field theory. We then proceed to briefly review the main features of Weyl's attempt to unify gravity and electrodynamics in a single geometric framework. In Sect. 3, we consider the extension of Raychaudhuri equation from a Riemannian setting to the case of a Weyl integrable space-time. This extension is invariant under Weyl transformations. In Sect. 4, we work out a gen-eralization of one of the Hawking-Penrose singularity theorems. Some applications of the formalism are given in Sect. 5. These correspond to cosmological solutions of scalar-tensor theories framed in a Weylian space-time. We conclude with some remarks in Sect. 6.

Weyl geometry
Weyl geometry arises from the weakening of one of the postulates of the Riemannian geometry, similarly as non-Euclidean geometry was conveived after the fifth postulate of Euclidean geometry was relaxed. The postulate we are referring to is the so-called the Riemannian compatibility condition, which states the following.
As is well known, this condition is equivalent to the requirement that the covariant derivative of the metric vanish, thereby implying that the length of a vector remain unaltered by parallel transport.
On the other hand, in Riemannian geometry, it will also be assumed that the connection ∇ be torsionless (or symmetric), i.e., that for any U, V ∈ T (M) the following condition holds: From these two postulates we are led to the important Levi-Civita theorem, which states that the affine connection is entirely determined from the metric [30]. In a certain sense, this theorem characterizes a Riemannian manifold. In 1918 the mathematician Hermann Weyl generalized the geometry of Riemann by introducing the possibility of change in the length of a vector through parallel transport. To implement this idea Weyl conceived the following new compatibility condition.

Postulate 2.2 (Weyl)
Let M be a differentiable manifold endowed with an affine connection ∇, a metric tensor g and a one-form field σ , called a Weyl field. It is said that ∇ is compatible (W-compatible) with g if for any vector fields U, V, Clearly, this represents a generalization of the Riemannian compatibility condition. Naturally, if the one-form σ vanishes, we reobtain (1). In this way, we have a generalized version of the Levi-Civita theorem given by the following proposition (see Ref. [31] for a proof): Theorem 2.1 (Levi-Civita extended) Let M be a differentiable manifold endowed with a metric g and a differentiable one-form field σ defined on M, then there exists one and only one affine connection ∇ such that: It follows that in a coordinate basis {x a } one can express the components of the affine connection completely in terms of the components of g and σ : where { a bc } = 1 2 g ad [g db,c + g dc,b − g bc,d ] denotes the Christoffel symbols of second kind. The next proposition gives a helpful insight on the geometrical meaning of the Weyl parallel transport.

Corollary 2.1
Let M be a differentiable manifold with an affine connection ∇, a metric g and a Weyl field of one-forms σ . If ∇ is W-compatible, then for any smooth curve α = α(λ) and any pair of two parallel vector fields V and U along α, we have where d dλ denotes the vector tangent to α.
By integrating the above equation along the curve α from a point p 0 = α(λ 0 ), we get Let us set U = V and denote by L(λ) the length of the vector V (λ) at an arbitrary point p = α(λ) of the curve. It is easy to check that in a local coordinate system {x a } Eq. (5) reduces to Consider now the set of all closed curves α : [a, b] ∈ R → M, i.e., with α(a) = α(b). The equation defines a holonomy group. If we want the elements of this group to correspond to an isometry, then we must require that for any loop. From Stokes' theorem it follows that σ must be an exact form, that is, there exists a scalar function φ, such that σ = dφ. In this case we have a Weyl integrable manifold. Weyl manifolds are completely characterized by the set (M, g, σ ), which will be called a Weyl frame. Let us remark that the Weyl compatibility condition (5) remains unchanged when we go to another Weyl frame (M, g, σ ) by carrying out the following simultaneous transformations in g and σ : where f is an arbitrary scalar function defined on M. The conformal map ( 10) and the gauge transformation (11) define classes of equivalences in the set of Weyl frames. It is important to mention that the discovery that the compatibility condition (5) is invariant under this group of transformations was essential to Weyl's attempt at unifying gravity and electromagnetism, extending the concept of space-time to that of a collection of manifolds equipped with a conformal structure, leading to the notion that space-time might be viewed as a class [g] of conformally equivalent Lorentzian metrics [9][10][11][12][13]15,16].

Weyl integrable space-time
We now consider the particular case of a WIST, where σ = dφ. As already mentioned, the set (M, g, φ) consisting of a differentiable manifold M endowed with a metric g and a Weyl scalar field φ will be referred to as a Weyl frame. In this case (11) becomes If we set f = −φ in the above equation, we get φ = 0. We refer to the set (M, g = e − f g, φ = 0) as the Riemann frame, because in this frame the manifold becomes Riemannian. On the other hand, it can easily be verified that Eq. (4) follows directly from ∇ α g μν = 0. This result has interesting and useful consequences. In fact, the metric γ = e −φ g, defined for any frame (M, g, φ), is invariant under the Weyl transformations (10) and (11) any geometric quantity built exclusively with γ is invariant. More generally, geometric objects such as the components of the curvature tensor R α βμν , the components of the Ricci tensor R μν , the scalar e φ R are invariant under the Weyl transformations (10) and (11).
It is important to note that because the Weyl transformations (10) and (11) define an equivalence relation between frames (M, g, φ) it seems more appropriate to look into the equivalence class of such frames rather than on a particular frame. In other words, a Weyl manifold should be regarded as a frame (M, g, φ) that is only defined "up to a Weyl transformation". In this way, when dealing with a certain Weyl manifold we choose a particular frame in the equivalence class and consider that only geometric entities defined in that frame which are invariant are of interest, since they can be regarded as representative of the whole class. From this point of view, it is more natural to redefine some Riemannian concepts to meet the requirements of invariance. This procedure is analogous to the one adopted in conformal geometry, a branch of geometry, in which the geometric objects of interest are those invariant under conformal transformation, such as, say, the angle between two directions [32]. Accordingly, one should naturally generalize the definition of all invariant integrals when dealing with the integration of exterior forms. For example, the Riemannian q-dimensional volume form defined as = √ −gdx 1 ∧ · · · ∧ dx q is not invariant under Weyl transformations, hence it should be replaced by = √ −γ e − q 2 φ dx 1 ∧ · · · ∧ dx q , and so on. 1 Likewise, in a Weyl integrable manifold it is more natural to define the concept of "length of a curve" in an invariant way. It follows that our familiar notion of proper time as the arc length of worldlines in four-dimensional Lorentzian space-time should be modified. Because of this, we shall redefine the proper time τ measured by a clock moving along a parametrized timelike curve x μ = x μ (λ) between x μ (a) and x μ (b), in such a way, that τ is the same in all frames. This suggests the following definition: It must be noted that the above expression can also be obtained from the special relativistic definition of proper time if we adopt the prescription η μν → e −φ g μν . It is clear that the right-hand side of this equation is invariant under Weyl transformations and that, in the Riemann frame, it reduces to the definition of proper time in general relativity. We, therefore, take τ , as given above, as the extension to an arbitrary Weyl frame of general relativistic clock hypothesis, i.e., the assumption that τ measures the proper time measured by a clock attached to the particle. It is now easy to see that the extremization of the functional (13) leads to the equations where μ αβ designates the Christoffel symbols calculated with g μν . Recall that in the derivation of the above equations the parameter λ must be chosen such that along the curve, which, up to an affine transformation, permits us to identify λ with the proper time τ . Surely, these equations are exactly those that yield the affine geodesics in a Weyl integrable space-time, since they can be rewritten as (4), are identified with the components of the Weyl connection. Thus the extension of the geodesic postulate by requiring that the functional (13) be an extremum is equivalent to assuming that the particle motion must follow affine geodesics defined by the Weyl connection μ αβ . Let us note that, as a consequence of the Weyl compatibility condition (3) between the connection and the metric, Eq. (15) holds automatically along any affine geodesic determined by (16). Since both the connection components As is well known, the geodesic postulate is not only concerned with the motion of particles, but it also determines the propagation of light rays in space-time. On the other hand, since the paths of light rays are null curves, one cannot use the proper time as a parameter to describe these curves. Thus light rays are supposed to follow null affine geodesics, which cannot be defined in terms of the functional (13), but, instead, they must be characterized by their behavior with respect to parallel transport. We naturally extend this postulate by simply assuming that light rays follow Weyl null affine geodesics.

The Raychaudhuri equation
The Raychaudhuri equation played a fundamental role in the derivation of the Hawking-Penrose singularity theorems. In the derivation, however, it is assumed right from the beginning that the geometry underlying the space-time is Riemannian. In this section, we investigate the extension of this equation to the case of a Weyl integrable space-time.
Let us first remark that the extension of the geodesic postulate to WIST assumes that the particle motion must follow Weyl time-like geodesics. In the following we shall consider a smooth congruence of time-like geodesics corresponding to the worldlines of a class of observers, parametrized by the invariant proper time τ defined in (13). Hence, the tangents to the congruence generate a tangent vector field V normalized to unit length. In order to keep the formalism invariant under Weyl transformations we shall choose the affine parameter of the congruence as the Weyl invariant arc length, i.e., we normalize V with respect to the invariant metric γ μν = e −φ g μν , Therefore, in a local coordinate system, a geodesic curve described by and we have the affine geodesics in a Weyl integrable spacetime shown in (16). Furthermore, once we have normalized the tangent vector field with the invariant metric γ μν , we obtain since ∇ α γ μν = 0.
Let us now consider, at some point p of M, the hypersurface orthogonal to the vector field V . As in the standard procedure, we define the operator of projection onto the hypersurface as As is well known, μν represents to the first fundamental form of the hypersurface induced by the metric γ and its role is to project any vector of T p M at p onto T p , the tangent space to the submanifold . 2 We proceed with the derivation of the Raychaudhuri equation in this new setting. We first need to consider a smooth one-parameter subfamily α s (τ ) of geodesics in the congruence of V and then define a deviation vector η that represents an infinitesimal spatial displacement from a given geodesic α o (τ ) to a neighboring geodesic in this subfamily. Once we have been given η in , we define the vector field η along this subfamily by Lie dragging it along V , that is, by requiring that From the definition of Lie derivative and the fact that the connection ∇ is assumed to be torsionless we are led to the following equation: 2 Note that we are using the invariant metric γ in order to ensure a invariant projection tensor. Accordingly, the operation of raising and lowering tensorial indices must be always carried out with γ . This guarantees that the duality between covariant and contravariant vectors is not modified by Weyl transformations.
which shows how the deviation vector changes along the congruence. An important role played in the investigation of the behavior of neighboring geodesics as we go along the congruence is played by the so-called deformation tensor Q μν , defined as In terms of Q μν we can rewrite (21) as which clearly means that the deformation tensor measures the failure of η μ to be parallel transported [33]. Furthermore, it is easy to see that this tensor is purely spatial, since Q μν V μ = 0 = Q μν V ν . Thus, Q μν is a tensor defined in the subspace of the tangent space perpendicular to V . Finally, in order to have a physical interpretation of some kinematical aspects of the congruence , we can decompose Q μν in its irreducible parts: where the parameters , σ μν , and ω μν are known, respectively, as the expansion, shear and vorticity of the congruence (for instance, Ref. [33]). From the above equation we have It is easy to see that in the case where the congruence is locally orthogonal to the hypersurface , we have ω μν = 0.
We are now particularly interested in the behavior of , which measures the expansion of the congruence and can tell us about the existence of conjugated points. Thus, in order to obtain the Raychaudhuri equation, we need to know the rate of change of along the congruence. Therefore, since, by definition, On the other hand, from the definition of the curvature tensor, we have one can readily obtain where we denote σ 2 = σ μν σ μν , ω 2 = ω μν ω μν , and the term R μν V μ V ν is usually referred to as the Raychaudhuri scalar. Setting˙ = V μ ∇ μ , we can write (30) in the following form, known as the Raychaudhuri equation: It is worth noting that (31) has the same form as in the case of a Riemannian space-time, although it must be recalled that the Ricci tensor is built with the Weylian connection. In fact, this is not surprising as we have redefined the proper time in an invariant way, using the invariant metric γ μν = e −φ g μν .
In the next section, we shall analyze the conditions that lead to singularities in the space-time.

Extending the singularity theorem
Because of the form of the Raychaudhuri equation takes in a Weyl integrable space-time, the description of conjugate points is the same as in Riemannian geometry. Thus, we have the following statement: if θ 0 = θ(τ 0 ) < 0 for some τ = τ 0 , R μν V μ V ν ≤ 0, then in a finite invariant proper time τ ≤ 3/|θ 0 |, the congruence will develop a conjugate point θ → −∞. As a matter of fact, the extension of any theorem from Riemannian geometry to a Weyl integrable can be trivially carried out by simpling considering the invariant metric γ = e −φ g. For instance, as we have already mentioned the extremization of the functional (13) leads directly to the Weylian geodesics. On the other hand, results coming from differential topology and concerning the causal structure of space-time that are valid in a Riemannian space-time are also valid in a Weyl integrable space-time since conformal transformations do not affect the light-cone structure nor the manifold orientability. Indeed, the proof of one of the most important results we are going to enunciate now, namely, the generalized Jacobi theorem, proceeds along the same lines of reasoning employed in the Riemannian case, where we merely replace g by the invariant metric γ = e −φ g [33].

. Then a necessary and sufficient condition for γ to locally maximize the invariant arc length between p and q is that γ be a Weyl geodesic without conjugate points between p and q.
Note the relation between the Raychaudhuri scalar and the extrinsic curvature of the submanifold orthogonal to geodesic congruence, which is represented by the mixed tensor Q For the case of a congruence of geodesics orthogonal to hypersurface and parameterized by the invariant arc length we have Q μ ν = ∇ ν V μ . Thus, the trace of Q μ ν is equal to , i.e., Q = Q α α = . Now we are ready to prove the following proposition. (M, g, φ) be a globally hyperbolic Weylian space-time with R μν V μ V ν ≤ 0 for all time-like vectors V . Suppose there exists a space-like Cauchy surface such that the trace of its extrinsic curvature (for the orthonormal congruence of past directed geodesics) satisfies B ≤ C < 0 over the whole surface, B and C being constants. Then no past directed time-like curve coming from can have an invariant arc length greater than 3/|C|. In particular, all past directed time-like geodesics are incomplete.

Theorem 5.2 Let
Proof Let us prove this theorem by contradiction. Suppose there exists a time-like curve λ = λ(τ ) coming from whose value of the invariant arc length τ at some point p ∈ λ is greater than 3/|C|. It is well known that as the set of curves joining two points in a globally hyperbolic manifold is compact, the invariant arc length function must have a maximum value for a given curve. Then this curve must be a Weyl geodesic. Therefore, there must be a geodesic γ (with invariant arc length greater than 3/|C|) joining p to . This means there are no conjugated points between p and . But, from Raychaudhuri's inequality, we know that γ must have conjugated points between p and , which is a contradiction. We then conclude that the original curve λ cannot exist.
In the case of general relativity where the space-time mathematical structure is that of Riemannian geometry, which is a special case of Weyl geometry when φ is a constant of motion, the geometric condition is equivalent to the so-called strong energy condition and requiring that Q < 0 is equivalent to assuming that in the course of the cosmic evolution the Universe underwent an expansion period, which seems to be a rather reasonable assumption. In view of the above, this leads to the conclusion that the Universe, as modeled by GR, must necessarily have had a beginning starting from a singular state. Let us now consider some different scenarios offered by two alternative gravitational theories, namely, the Novello-Oliveira-Salim-Elbaz's theory (N) [17,34] and a recent geometrical approach to scalar-tensor theory (GST), both inspired by the idea that space-time can be described by Weyl integrable geometry [35].

Novello's theory
Novello's theory starts with the action where R denotes the Weylian scalar curvature, φ represents the Weyl scalar field, ξ is a dimensionless parameter, while L m stands for the Lagrangian of the matter fields. The form of L m is determined from the corresponding Lagrangian in special relativity by replacing the ordinary derivatives by covariant derivatives with respect to the Weyl connection. The field equations are given by where we have set λ = 3 2 − 2ξ.
From the above equations it is not difficult to verify that Clearly, in this theory (32) does not imply the violation of the strong energy condition (33), but rather we have the following situation. Consider the conditions below: If any one of these conditions is violated, then the solution may correspond to a non-singular space-time. Consider, as an example, the vacuum solution of the field equations (35) and (36) obtained in Ref. [34]. By assuming homogeneity and isotropy we can write This leads to the equationṡ 2aä (a 3φ ) ,0 = 0.
By integrating the last of these equations we getφ = ζ a −3 , where ζ is a constant. Now, if 4ξ − 3 > 0, it follows thaṫ with (a 0 ) 4 = ( 4ξ −3 12 )ζ 2 , thus implying that a(t) ≥ a 0 . The model describes a non-singular bouncing universe, a scenario in which after undergoing a period of contraction, dominated by the scalar field, the universe reaches a minimum value, and then starts to expand, at an inflationary rate, until the radiation dominates over the scalar field and the scale factor begins to follow the standard cosmic evolution [34]. Clearly, the non-singular behavior of the model is a consequence of the fact that the condition ξ ≤ 1/2 is violated. In order to display the bouncing character of this model, it is convenient to solve (46) in terms of the conformal time η, defined as dη = dt/a(t). It can be shown that the scale factor and the scalar field are given, respectively, by

Geometrical scalar-tensor theory
The geometrical approach to scalar-tensor theory starts with the Brans-Dicke action where φ is assumed a priori to be a geometrical field, i.e., an intrinsic part of the space-time geometry, ω is a dimensionless parameter and k * = 8π c 4 . By applying the Palatini variational method, one obtains the Weyl integral compatibility condition [35], Naturally, the action (49) can easily be extended to accommodate a scalar potential V(φ) and to allow for a functional dependence of ω on φ, thus leading to [36] The field equations obtained from (49) are given by where denotes the d'Alembertian operator calculated with the Weyl connection.
By taking into account (52) one easily obtains As we can see, here again the (32) does not require the violation of the strong energy condition (33). Thus, if we assume that (33) holds, then any space-time described by these equations will satisfy R μν V μ V ν ≤ 0, as long as We conclude that in this case the singularity theorem applies.
On the other hand, if (55) is violated, then the solution may correspond to a non-singular space-time.
In the following let us consider some solutions to (52, 53), for some choices of the potential V(φ) in the case ω(φ) = constant and T μν = 0. These solutions correspond to homogeneous and isotropic models in GST theory and are obtained in the Riemann frame (M, g = e − f g, φ = 0), in which case |V | 2 = e φ . If we take the line element written as in (42), then the field equations (52) reduce to while (53) gives where = 0, ±1, according to the curvature of the spatial section. 3 Solutions to Eqs. (56)-(58) with flat spatial section = 0 were obtained for the following choices of the scalar potential: e −2φ , e −(2+λ)φ V o , e −2φ (mφ 2 + ), and 2λe −2φ (φ 2 + 2ω/3) 2 , where , V o , λ, and m are constants. These, in the Riemann frame, correspond to the cosmological constant, the dilaton field, the massive scalar field, and a field with a quartic interaction. These three cases are displayed in the tables below, where the presence of a singular behavior is determined according to whether or not (55) holds (Table 1). 3 It is interesting to note here that from the above we can obtain the following equation: This equation might be useful to set the possible values of ω.  Table 2 The parameters are the same as in Table 1 and a o is also a constant

Final remarks
The Hawking-Penrose singularity theorems are a direct consequence of Einstein's theory of gravity. Given the important role they have played in our understanding of the Universe, as modeled by general relativity, it is also of interest to find the analogous of these theorems in alternative gravity theories. Singularity theorems and energy conditions have been studied in connection with Brans-Dicke theory, perhaps the most popular scalar-tensor gravity theories. However, their meaning remains still controversial due to the question of whether or not the Einstein frame and the Jordan frame are physically equivalent [37,38]. In the present geometrical approach, this controversy does not arise, as the physical entities defined in the theory are naturally invariant under frame transformations [39][40][41][42] (Table 2). Finally, it has been shown that Weyl integral space-time seems to be the natural geometrical scenario that arises in the context of scalar-tensor theories, at least, whenever one has a theory that assumes a non-minimal coupling of the scalar field to the Ricci curvature [35]. The basic principle which underlies the determination of the space-time geometry would be, in this case, the adoption of the Palatini variational principle, which has played an important role in recently proposed modified theories of gravity [43]. The investigation of the behavior of space-time singularities in this new framework was one of our main motivations to write the present paper.