The Singularity Theorems of General Relativity and Their Low Regularity Extensions

On the occasion of Sir Roger Penrose’s 2020 Nobel Prize in Physics, we review the singularity theorems of General Relativity, as well as their recent extension to Lorentzian metrics of low regularity. The latter is motivated by the quest to explore the nature of the singularities predicted by the classical theorems. Aiming at the more mathematically minded reader, we give a pedagogical introduction to the classical theorems with an emphasis on the analytical side of the arguments. We especially concentrate on focusing results for causal geodesics under appropriate geometric and initial conditions, in the smooth and in the low regularity case. The latter comprise the main technical advance that leads to the proofs of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$C^{1}$\end{document}C1-singularity theorems via a regularisation approach that allows to deal with the distributional curvature. We close with an overview on related lines of research and a future outlook.


Introduction
The singularity theorems of General Relativity (GR) are commonly counted among the 20th century milestones in mathematical physics.They comprise a body of rigorous results in Lorentzian differential geometry that, under physically reasonable conditions, imply the occurrence of a "singularity" in the sense of causal geodesic incompleteness of the spacetime manifold.At the time of its appearance, the first singularity theorem proved by Roger Penrose in [Pen65] came as a surprise to the community, debunking the widely held belief that the singularities encountered in many of the exact solutions of GR were a mere artefact of their high degree of symmetry.It has even been argued e.g. in [Sen12,Sec. 15.1] that Penrose's theorem, that appeared exactly 50 years after the birth of GR, was the first true post-Einsteinian contribution in the sense that it was not foreseen by its founding father.And indeed, it is this paper that won Roger Penrose (one half of) the 2020 Nobel Prize in Physics.
Despite its brevity of less than three pages, it not only shaped the way we now primarily think of singularities-via geodesic incompleteness-but also put forward the fundamental new idea of a trapped surface, which has stimulated mathematical GR up to the present day.The theorem itself roughly says that, in a situation that physically amounts to the formation of a black hole-intuitively a region in spacetime of so strong gravity that not even light can escape-there necessarily exists an incomplete null geodesic, i.e., a light ray that ends suddenly1 .
A short time later, Stephen Hawking realised that one could apply a similar reasoning to an everywhere expanding universe and showed that in such a spacetime there has to be a timelike geodesic that is past incomplete.This theorem, known as the Hawking singularity theorem [Haw67], is generally thought of as mathematical evidence for the occurrence of a big bang.Then, during the next couple of years, Robert Geroch, George Ellis and others helped to shape a body of results which we now call the classical singularity theorems-for some historical details see e.g.[SG15,Sec. 5].Notably, in 1970 Hawking and Penrose, in their only joint paper [HP70], proved the most refined of these results which is known as the Hawking-Penrose singularity theorem.
Over the decades the singularity theorems have not only become an integral part of GR, but are still an area of active research, cf.e.g.[SG15,Sec. 8].There is, of course, an extensive textbook coverage available and we only mention the early classics [Pen72,HE73], and the more mathematically oriented standard accounts [BEE96, O'N83, Cla93, Kri99], as well as the review articles [Sen98, SG15, MAS15, Sen21, Daf21, Lan21], the latest three of which have appeared on the occasion of the 2020 Nobel Prize.
In this work we offer a guided tour for the mathematically oriented reader to the classical singularity theorems and to their recent extensions for Lorentzian metrics of low regularity.The study of the latter is motivated by the discussion of the "character of the singularities" (cf.[HE73,Sec. 8.4]) predicted by the theorems, which we will recall after presenting the classical results.
In some more detail we briefly collect the necessary preliminaries on GR and Lorentzian geometry in Section 2, in a way especially suited to the thread of this account.In Section 3 we discuss the main classical theorems, namely the one of Penrose, the one of Hawking, and the one of Hawking-Penrose, together with the main arguments that enter their proofs.In particular, we will discuss the focusing of geodesics under certain curvature and other geometric conditions (Subsection 3.2) as well as the relevant parts of Lorentzian causality theory (Subsection 3.4).Then, in a brief analysis of the conditions and statements of the theorems in Subsection 4.1, we motivate the quest for their low regularity extensions which are the topic of Section 4.Here we present an overview of these results which were obtained in the last couple of years, and, again provide some key insights on the techniques and arguments employed in their proofs.The main thread is laid out in Subsections 4.2-4.6 and 4.10, while the remaining rather technical subsections are directed towards the more initiated reader.Finally, in Section 5 we provide a summary, some conclusions, and an outlook to current and future lines of research in the direction at hand.

Preliminaries: General Relativity, geometry, singularities
Here we first collect in a nutshell the basics of GR, Albert Einsteins theory of space, time, gravitation, and matter.Then we turn to its geometric foundations, that is Lorentzian differential geometry and review the necessary preliminaries for our approach to the singularity theorems.Finally, we briefly discuss the very notion of a singular spacetime in GR.

Some basics of General Relativity
The stage of the theory is an n-dimensional2 Lorentzian manifold (M, g) where the metric g is a symmetric, non-degenerate (0, 2)-tensor field, assigning to any point p in the smooth manifold3 M a scalar product g p on the tangent space T p M at p.This scalar product, however, is not positive definite, but has signature (−1, +1, . . ., +1), which gives rise to the following distinction of tangent vectors: We call v ∈ T p M spacelike (timelike) if v, v ≡ g p (v, v) > 0 (< 0).We call it null if v, v = 0 but v = 04 and causal if it is timelike or null.These notions extend naturally to vector fields and, via their tangent vectors, to sufficiently smooth curves.
Given the existence of a timelike vector field on M we can use it to define a time orientation, i.e. a smooth choice of one of the two connected components of the set of causal vectors at any p ∈ M , called the future cone.A spacetime then is a smooth Lorentzian manifold (M, g) together with a time orientation.
In this note we will also be interested in non-smooth spacetimes, that is smooth Lorentzian manifolds with a time oriented metric of regularity below C 2 , i.e., twice continuously differentiable.In fact, it makes sense to regard such metrics as being of low regularity since the bulk of classical smooth Lorentzian geometry extends verbatim to C 2 -metrics.We will in particular be interested in metrics of regularity C 1,1 , i.e., metrics that are continuously differentiable (C 1 ) and, in addition, have first derivatives that are locally Lipschitz continuous5 , as well as in merely C 1 -or C 0 -metrics.We will always assume the time orientation to be induced by a smooth timelike vector field and we do not reduce the regularity of the differential structure itself, since this is no loss of generality: Any C k -manifold (k ≥ 1) possesses a unique smooth structure compatible with the given C k -structure.
At the heart of GR are the field equations which connect the curvature of spacetime to the matter and energy it contains.To make this statement precise, recall the Riemannian curvature tensor of (M, g) where X, Y, Z ∈ X(M ) = Γ(M, T M ) are smooth vector fields and ∇ denotes the Levi-Civita connection of g.Then we have the Ricci tensor where here and in the following (E i ) n i=1 denotes a (local) orthonormal frame field.Finally, writing S for the curvature scalar, i.e., the contraction of the Ricci tensor, we arrive at the Einstein equations Ric − 1 2 S g + Λ g = 8 π T. (2.3) Here T is the energy momentum tensor which encodes information on the matter and energy content of spacetime.We avoid a discussion of T since we are mostly interested in the vacuum case, i.e., when there is no matter or energy and we have T = 0. Also we will generally assume the cosmological constant to vanish and set Λ = 0.The simplest vacuum solution of (2.3), of course, is flat Minkowski space, i.e., R n with g = diag(−1, 1, . . ., 1), which is the stage of the special theory of relativity.Much more interesting is already the famous Schwarzschild metric which gives the entire one-parameter family of spherically symmetric vacuum solutions of (2.3).It models spacetime outside a spherically symmetric body of mass m but the Kruskal extension also is the basic model of a static black hole.
Given these foundations, GR over its first decades was developed mainly as a geometric theory and a large zoo of solutions to the field equations were explicitly found and studied.This branch of GR known as exact solutions (see e.g.[SKM + 03, GP09]) was complemented by more analytic approaches that came into gear with the celebrated local existence result for the Einstein equations by Yvonne Choquet-Bruhat [FB52].During the past decades the PDE aspect of the theory has become more and more prevailing, producing many remarkable existence and stability results, see e.g.[Rin15] for an overview.The past decades have also seen a growing impact of numerical methods in GR (see e.g.[BS10]) and it finally was a combination of analytical and numerical methods that provided the theoretical background for the celebrated first direct observation of gravitational waves in 2015 [Aea16], see e.g.[BGY17].
Returning to the foundations of GR we proceed with a famous quote by John Wheeler cf.[MTW73, p. 5], in which he concisely describes its core: Space tells matter how to move.Matter tells space how to curve.
While the latter statement refers to the field equations (2.3) introduced above, we now turn to the first one.Indeed matter, more precisely test particles of negligible (rest) mass and inner structure, often called (freely falling) observers, move along timelike geodesics, while light rays trace out null geodesics in spacetime.Collecting the basic facts on geodesics in Lorentzian manifolds we now begin our brief primer on the geometrical background of GR.

Some basics of Lorentzian geometry
Geodesics are curves γ : I → M defined on an interval I that are self parallel, i.e., their tangents γ satisfy ∇ γ γ = 0.The geodesic equation written in some local chart (U, x i = (x 0 , . . ., x n−1 )) takes the form ẍi + Γ i jk ẋj ẋk = 0, (2.4) with the Einstein summation convention in effect.Here the Christoffel symbols are where g jk and g jk denote the coefficients of the spacetime metric and its inverse, respectively.Also we have followed the common habit to denote the components x i • γ of γ simply by x i .Equation (2.4) is a secondorder, non-linear ODE, hence given any point p ∈ M and any tangent vector v ∈ T p M there is a unique geodesic γ v defined on a maximal half-open interval I v := [0, β) assuming this data, i.e., γ v (0) = p, and γv (0) = v.We say that γ v is complete if β = ∞.More generally, we call any geodesic γ : I → M extendible if there exists a geodesic γ : J → M with J I and γ | I = γ.Otherwise γ is called inextendible.Finally, γ is called complete if it may be extended to all values of its (or any other affine) parameter, that is γ : R → M , and we call M itself geodesically complete it this is true for all geodesics.
Next we discuss the exponential map at p defined via where the domain is It is a basic fact that exp p for any p is a diffeomorphism from some open zero neighbourhood U ⊆ T p M onto U = exp p (U ) and we call U a normal neighbourhood of p if U is star shaped.Any point q in a normal neighbourhood U of p is connected to p by a unique radial geodesic γ : [0, 1] → U for which we have γ(0) = exp −1 p (q).A neighbourhood U is called (geodesically) convex if it is normal for all its points, which implies the existence of a unique geodesic in U between any pair of its points.Finally, every p ∈ M possesses a base of convex neighbourhoods.
Up to this point there do not occur any differences between Riemannian and Lorentzian manifolds.However, the Gauss Lemma, which states that the exponential map is a radial isometry, has quite different consequences in the two cases.While in the Riemannian case it implies that radial geodesics minimise the Riemannian distance in convex neighbourhoods, it here leads to radial causal geodesics maximising the Lorentzian distance (sometimes called the time separation function).
To explain this statement in some detail, we need to introduce some basic notions from causality theory, i.e., the theory of futures and pasts of points in a spacetime.For more details see [MS08] and the authoritative source [Min19b].For p, q ∈ M we write p ≪ q if there is a future directed timelike curve 6 from p to q, and q ≤ q if there is a future directed causal curve from p to q or if p = q.We then define the chronological and causal future of p via respectively.Analogously one defines the chronological and causal past I − (p) and J − (p) of a point, and for a set A ⊆ M one defines I ± (A) = ∪ p∈A I ± (p) and analogously for J ± (A).It is a fundamental fact that I ± (p) is open but J ± (p) need neither be closed nor open.Also the so-called push-up principle holds, which can bee seen as an improved transitivity of the causality relations: if p ≪ q and q ≤ r (or if p ≤ q and q ≪ r), then p ≪ r, and the name refers to the fact that a curve that has a timelike and a null part can be "pushed up" to give an overall timelike curve.If γ : I → M is a (sufficiently smooth) curve, its length is defined as Given two timelike related points p ≪ q there are always arbitrarily short future directed timelike curves connecting them-just take curves arbitrarily close to a piecewise null zig-zag curve, which always has vanishing length.On the other hand, if p and q lie in a convex neighbourhood U then the unique future directed timelike radial geodesic between them is the longest curve in U from p to q.So it makes sense to define the Lorentzian distance on M × M by where the sup runs over all future directed causal curves γ from p to q. Contrary to its Riemannian sister, the Lorentzian distance is not symmetric and it satisfies the reverse triangle inequality, d(p, q) ≥ d(p, r) + d(r, q) for all p ≤ r ≤ q, (2.8) i.e., detours make curves shorter rather than longer.Also the sup is not finite if the causality of the spacetime behaves badly, e.g. if there are closed timelike curves.
There are even more striking new phenomena in Lorentzian geometry when contrasted with the Riemannian case.Recall that the Riemannian distance d : M ×M → R given by d(p, q) = inf{L(γ)} (where the inf runs over all curves γ connecting p and q) is a continuous metric that induces the manifold topology.Moreover, the Hopf-Rinow theorem asserts that a Riemannian manifold M is complete as a metric space iff it is geodesically complete and iff M is geodesically connected, i.e., every pair of points can be joined by a minimising geodesic γ, that is In the Lorentzian case, as discussed above, d is not a metric and it fails to be upper semicontinuous in general, while it is lower semicontinuous where it is finite.There is, however, a class of spacetimes where d is finite and continuous, namely the globally 6 While most textbooks base causality theory on piecewise smooth curves it has turned out to be more economically [Chr11,Min19b] to use locally Lipschitz curves.Hence for us a curve is timelike (causal, null, future or past directed) if it is locally Lipschitz and γ(t), which exists almost everywhere by Rademacher's theorem, is timelike (causal, null, future or past directed) almost everywhere.hyperbolic ones.These are defined to be causal7 (i.e., there are no closed causal curves) with the so-called causal diamonds J(p, q) := J + (p) ∩ J − (q) being compact for all p, q.Also, by the Avez-Seifert theorem, globally hyperbolic spacetimes are causally geodesically connected, that is, any pair of points p ≤ q can be joined by a causal geodesic γ that is maximising, i.e., L(γ) = d(p, q).Conversely, and just as for minimising curves in the Riemannian case, maximising curves may be reparametrised as causal geodesics.
The notion of geodesic completeness in the Lorentzian case can be defined individually for spacelike, timelike, and null geodesics, giving rise to three independent notions.Also, causal geodesic completenss does not imply causal geodesic connectedness.So, in some sense, global hyperbolicity is the Lorentzian counterpart of Riemannian completeness, while incompleteness is a central notion in the singularity theorems, as we shall discuss next.

Singularites and incompleteness
Despite vivid examples such as the Schwarzschild metric and cosmological models of expanding universes, it is a somewhat subtle matter to define the general notion of a "singularity" in GR.Intuitively it should be a point where the curvature becomes unbounded but such a scenario is literally incompatible with the dynamical picture of spacetime in GR: We can only ever speak of a point in spacetime, if we have solved the field equations around it and hence found the corresponding manifold and metric structure.So, rather than being a point in spacetime a singularity should be seen as some kind of singular boundary point of spacetime8 .The quest then is, how to detect the occurrence of a singularity from within spacetime, that is by properties of the spacetime itself-for a detailed discussion see e.g.[HE73,Sec. 8.1] and the classic [G68].
To cut a long story (see e.g.[Lan21, Sec.2.3]) short, after Penrose's seminal paper [Pen65] it has-mainly under the influence of Hawking-become standard to define singularities via causal geodesic incompleteness.More precisely, we call a spacetime singular, if it contains an incomplete causal geodesic.
In addition to giving a clear geometric condition, this definition is also physically reasonable: A future incomplete causal geodesic corresponds to a freely falling observer or to a light ray that suddenly ends its existence.In the past case it corresponds to an observer or light ray that suddenly pops into existence from nowhere.Both situations have to be considered as being even more objectionable than a blow up of curvature, and the general point of view has become to regard causal geodesic completeness as a minimal condition for a spacetime to be "free of singularities".
The above definition, however, also has severe drawbacks, see e.g.[Sen98, Sec.3] for a detailed discussion.In particular, the link between singularities and divergence of the curvature is lost.In fact, incompleteness could occur for trivial reasons, e.g.cutting out one point from an otherwise perfectly fine spacetime.This immediately leads us to the following notion: We call a spacetime extendible if it can be isometrically embedded as an open submanifold into a larger spacetime9 .It is a fact, [BEE96, Prop.6.16] that timelike (or null) geodesic completeness implies inextendability.So we see that also in general, extendability can be a source of incompleteness and it becomes an important issue in the wake of the singularity theorems to exclude it.We shall return to this discussion in Section 4.1, but we now head on to discuss the classical theorems.

The classical singularity theorems
In this section we review the classical theorems.We start with the one of Hawking, which corresponds to the cosmological situation, since it is technically somewhat easier to formulate and to prove.We then proceed to the Penrose theorem which covers the gravitational collapse scenario, and finally, cover the most refined of the classical results, the Hawking-Penrose theorem.We start by revealing the general structure behind essentially all singularity theorems.

The basic structure of the singularity theorems
It has long been observed that all classical singularity theorems share the same structure.This point has been made most clearly by José Senovilla, who formulated a "pattern singularity theorem" in [Sen98, Thm.6.1] to analyse the various statements, their conditions and their conclusions in a well organised way.
Theorem 3.1 (Pattern singularity theorem).Let (M, g) be a spacetime such that the following hold: (E) a condition on the curvatue, also called energy condition, (C) a condition on the causality, and (I) an initial or boundary condition.
Then M is causal geodesically incomplete.
Without specifying these conditions in detail we can nevertheless give a first idea of how the corresponding proofs work, that is, we explain how the conditions (E), (C), and (I) conspire to contradict geodesic completeness.
The initial condition (I) serves the purpose that some causal geodesics start to focus towards each other.The energy condition (E) then implies that this focusing goes on until eventually a causal geodesic develops a focal or conjugate point and consequently stops maximising the Lorentzian distance.On the other hand the causality condition (C) implies the existence of maximising geodesics, at least in some region of spacetime.Now to dissolve this contradiction we have to accept that some geodesics stop existing before they reach a conjugate or focal point, that is, they are incomplete.
Let us now examine each of the three conditions in some more detail and introduce the corresponding notions in a meaningful and precise manner.

Energy conditions and the focusing of geodesics
We start with condition (E) which actually is a condition on the curvature of spacetime.It is, however, called energy condition since via the field equations (2.3) it corresponds to a condition on the energy-momentum tensor T .Physically these conditions roughly amount to the fact that gravity is attractive, at least on average, and most "reasonable" matter models will satisfy them-for an extensive discussion see e.g.[HE73, Sec.4.3] and, in the context of the singularity theorems, [Sen98, Sec.6.2].
Here we only introduce the most important of these conditions.
We say that a spacetime statisfies the strong energy condition (SEC), if Ric(X, X) ≥ 0 for all timelike vectors X, (SEC) and we say it satisfies the null energy condition (NEC) if Ric(X, X) ≥ 0 for all null vectors X. (NEC) Note that by continuity (SEC) implies (NEC) and hence nonnegativity of Ric(X, X) for all causal X.So this terminology is consistent with the habit to formulate the (SEC) for all causal vectors10 .It is worthwhile to note that the energy conditions are the only place where the field equations come into play, so that the validity of the singularity theorems widely exceeds GR and continue to hold for alternative theories, as long as they imply the corresponding curvature conditions.
In the context of the singularity theorems the curvature condition is needed to exert a focusing effect on geodesics and we shall explain this in some detail, for more information see e.g.[Kri99, Sec.4.6], [BEE96, Ch. 9].To begin with, we introduce the central notion of Jacobi fields, i.e., vector fields J along a geodesic γ that satisfy the Jacobi equation J + R(J, γ) γ = 0 , (3.1) where J denotes the (iterated induced) covariant derivative along γ.Jacobi fields are in a one-to-one correspondence with geodesic variations of γ, which picture the situation more vividly.Let γ : [a, b] → M be any curve, then a variation of γ is a two-parameter map (with δ some small positive number) such that x(t, 0) = γ(t) for all a ≤ t ≤ b, which we also call the base curve of x.For fixed s we call the t-parameter curves t → x(t, s) longitudinal, and the s-parameter curves for fixed t, s → x(t, s) transverse.The tangents of the transverse curves at the base curve give the variational vector field V (t) = d ds | s=0 x(t, s) along γ.Now given a geodesic γ then any Jacobi field is the variational vector field of a geodesic variation of γ, that is, a variation where all longitudinal curves are geodesics as well.
Essential for the focusing of geodesics is the notion of conjugate points: two points γ(a), γ(b) are called conjugate along γ if there exists a nontrivial Jacobi field vanishing at a and b.This condition is equivalent to the existence of a geodesic variation x with x(a, s) = γ(a) for all s and V (b) = 0, which is further equivalent to the exponential map exp γ(a) being singular at b γ(a).
A geodesic variation as above with V (b) = 0 can be pictured as an "almost meeting point" of nearby geodesics, all starting out at γ(0).Suppose for the moment that a meeting point actually occurs along a maximising causal geodesic γ, so that there is another maximising causal geodesic σ with the same endpoints (i.e., σ(a) = γ(a), σ(b) = γ(b)).Then γ fails to maximise the Lorentzian distance beyond γ(b) since the concatenation of σ | [a,b] with γ | [b,b+η) is a geodesic from γ(0) to some γ(b + η) with a break point at γ(b) and the same length as γ.However, being broken it cannot be maximising and so γ can't be either.Making this argument precise one arrives at the following central statement: Proposition 3.2.A causal geodesic fails to maximise the Lorentzian distance after its first conjugate point.
From here it is evident that the curvature can be linked to the global structure of a spacetime by means of conjugate points.In fact, in the Riemannian case the analogous observation leads e.g. to the theorems of Hadamard and Myers, while in the Lorentzian case it is intimately related to the singularity theorems, cf.e.g.[MAS15].
At this point it should perhaps be pointed out that a causal geodesic from a point stops being maximising if and only if either (a) there exists a distinct causal geodesic between the same endpoints of the same length (as argued above) or (b) the geodesic encounters a conjugate point.But the cut-locus of a point p in a Lorentzian manifold (i.e., the set of points where the geodesics emanating from p stop maximising) is a closed set of measure zero with the set of points which can be reached from p by two distinct maximising geodesics being dense.Therefore, almost all geodesics that stop maximising do so due to (a) and hence do so even before they encounter their first conjugate point.However, since (a) is related to the global geometry of the manifold, there is no way to detect such points via estimates on the curvature.So the power of conjugate points lies the fact that they lead to geodesics no longer being maximising and that we can estimate when they occur.
Motivated by this observation we now introduce the analytical tools to find conjugate points, for a detailed and systematic analysis based upon Lorentzian Morse index theory see [BEE96, Ch. 10], [Oha22].
To begin with, we observe that the relevant information on conjugate points is contained in the (n − 1)-dimensional subspace of Jacobi fields vanishing at a given point and taking values in γ(t) ⊥ := {v ∈ T γ(t) M : v, γ(t) = 0}.Now, to also properly deal with the null case (where γ(t) ∈ γ(t) ⊥ ) it is useful to work on the quotients respectively.Note that in the timelike case [ γ(t)] ⊥ coincides with γ(t) ⊥ hence is of dimension n − 1, while in the null case its dimension is n − 2. We will henceforth use the letter d to denote the dimension of [ γ] ⊥ simultaneously in both cases.Observe that the restriction of the metric g | [ γ] ⊥ is positive definite in both cases.
Next, to represent the information contained in Jacobi fields most economically we collect them into Jacobi tensor classes.That is, classes of (1, 1)-tensor fields [A] : [ γ] ⊥ → [ γ] ⊥ for which the tensor Jacobi equation where the shear σ is defined as The key observation now is the following: The second and the third term on the right hand side of (3.6) are non-positive.If we assume the first term to be nonpositive as well by demanding that Ric( γ, γ) ≥ 0, e.g.via the energy conditions (SEC) or (NEC), we can easily "generate" conjugate points.Indeed, if θ(a) is negative at some parameter value a, then it will diverge to −∞ in finite parameter time.More precisely, (SEC)/(NEC) and (3.6) imply θ ≤ −θ 2 /d which upon integration from a to some t > a gives To sum up, we have seen that if along a causal geodesic a negative expansion occurs, then given the energy conditions, eventually a conjugate point will arise implying the failure of the geodesic to be maximising beyond it.Now, to "generate" a point with negative expansion is the purpose of the initial conditions to which we turn next.

Initial conditions
We now proceed to condition (I), which is the ingredient that gets the focusing effect of geodesics started, i.e., which produces a negative θ(a) at some point γ(a) along a causal geodesic γ.
Beginning with the simplest case, consider a spacelike hypersurface S in spacetime (i.e., g | T S Riemannian) and its future directed unit normal ν, i.e., ν, ν = −1.Now at a given point p in S we consider the geodesic γ ν starting in p with future normal unit velocity, i.e., γ ν (0) = p, γν (0) = ν(p).Further consider a Jacobi tensor11 A along γ ν with A(0) = id and Ȧ(0) = −S ν (p), where S ν is the shape operator of S given by S ν (X) := −∇ X ν.A standard calculation now gives where H is the mean curvature vector of S and we call k the future convergence of S. Now, if we suppose that H(p) is past pointing timelike or, equivalently, that the convergence k(p) is positive, we have achieved θ(0) < 0 and we have kicked off focusing.
To formulate a corresponding statement we still have to introduce the analogue of a conjugate point along a geodesic in the current "endmanifold" case, i.e., the concept of a focal point.More precisely, a point q is called focal to a hypersurface S along a timelike geodesic γ that starts at some p ∈ S and normal to it, if there is a nontrivial normal Jacobi field J along γ with J(0) = 0 and J ′ (0) = ∇ J(0) ν that vanishes at q = γ(b).Analogously to the case of a conjugate point, such a Jacobi field corresponds to the variational vector field V of a variation x of γ through geodesics that all start at S and normal to it with V (b) = 0. Equivalently the normal exponential map exp ⊥ p : , is singular at b γ(0).Now, just as in the conjugate-point case, a geodesic γ ν starting normal to a hypersurface S stops maximising the Lorentzian distance to S, denoted by d(S, .), after its first focal point.Also observe that if a geodesic does not start out normally from S, it is not maximising at all.A typical focusing result then reads as follows: Proposition 3.4 (Focal points from negative expansion).Let γ : [0, β) → M be an inextendible timelike geodesic starting at p normally from a spacelike hypersurface S and suppose Ric( γ, γ) ≥ 0 for all t.If θ(0) is negative (equivalently if the convergence k(p) is positive) then there is a focal point γ(b) to S for some b ≤ −(n − 1)/θ(0), provided b < β.Consequently, γ stops maximising the Lorentzian distance to S (the latest) at b if it exists that long.
The above focusing result will turn out to be at the analytical core of Hawking's singularity theorem, and we now turn to the initial condition for Penrose's theorem.In fact, the notion of a trapped surface-to be introduced now-is one of the great innovations of the celebrated paper [Pen65] and its influence on the development of GR cannot be overestimated, see e.g.[SG15, Sec.7.2] for a brief description of its impact.It gives a precise mathematical formulation to the idea of a spacetime region of so strong gravity that not even light can escape the gravitational pull.Its key model are the "trapped spheres" inside the horizon of a Schwarzschild black hole.
To begin with, we consider a spacelike submanifold P of codimension 2. Since its normal bundle in M is Lorentzian, there are at any point of P two linearly independent future pointing null vectors normal to P .We call them ν 1 , and ν 2 and normalise them via the condition ν 1 , ν 2 = −1.The corresponding null vector fields on P give rise to two distinct families of null geodesics which, physically speaking, represent the two families of light rays emanating orthogonally from S into the future.Interpreting P for the moment as the surface of a star at a given instant of time, our experience leads us to think that the congruence of light rays sent towards the center of the star should converge, while the second "outgoing" congruence should diverge.However, if the effect of gravity is strong-which in the present picture manifests itself in the extrinsic curvature of the spacelike hypersurface representing the instant of time-both congruences will converge.In fact, this is the condition of P being future trapped.Mathematically we can expressed this by demanding the expansions of both the families of normal future pointing null geodesics θ 1 and θ 2 to be negative on all of P .Equivalently, we can demand the mean curvature vector of S to be past pointing timelike.
This alone is not a useful property as is demonstrated by the trivial example of the intersection of two past lightcones in Minkowski space.However if we ask the hypersurface to be closed, i.e., compact without boundary (as suggested in the above example) it becomes a mighty tool as we shall see below.
To get a vivid picture of the contraction property, we express it via the variation of the area A of P along the flow of a vector field ξ, which by a standard formula is given by (3.9) This equation, in the Riemannian case, of course expresses the fact that minimal surfaces (H = 0) are critical points of the area functional.In the Lorentzian situation, however, H being past pointing and timelike implies that the variation of the area along any future directed null vector field ξ is negative.
We finally give the official definition of a closed trapped submanifold of arbitrary codimension 1 < m < n, and reserve the term closed trapped surface for the case m = 2.
Definition 3.5 (Closed trapped submanifold).A compact without boundary, spacelike submanifold P is called closed future trapped submanifold if its mean curvature vector field is past pointing timelike on all of P .
The idea of lightrays being trapped can also be expressed in purely causal terms via the notion of a future trapped set, i.e., a closed, achronal set A such that its future horismos Here achronal means that no two points in A are timelike related.By push up E + (A) is achronal, and if it is compact, then A has to be compact itself since , although equality holds locally and also in globally hyperbolic spacetimes.Now, using a focusing argument much like the one that lead to Proposition 3.4 one may show that the analytic concept of a closed trapped surface P under appropriate conditions implies the causal concept of trappedness.Indeed, by a variational argument E + (P ) is generated by conjugate-free null geodesics: to any point q ∈ E + (P ) there runs a null geodesic from P in E + (P ) to q that has no conjugate points before q.Now by focusing all these null geodesics eventually do develop conjugate points which implies that E + (P ) is contained in the normal exponential image of a compact set and by closedness is compact itself.We have hence argued for the following result to hold, which can be generalised to trapped submanifolds of codimension m > 1 under an additional curvature condition, see [GS10, Prop.3]: Proposition 3.6 (Trapped set from trapped surface).Let (M, g) be a future null complete spacetime where (NEC) holds.Then any achronal closed future trapped surface is also a future trapped set.
The occurrence of geodesic focusing as such, however, does not lead to a singularity, and that is where some more causality theory comes into play.

Causality theory
In the discussion leading to Proposition 3.6 we have used the inexistence of conjugate points in E + (P ).More generally, we need to exclude the existence of focal or conjugate points in certain subdomains of spacetime to infer the existence of an incomplete geodesic from geodesic focusing.
To give a proper account we first introduce the notion of a Cauchy surface, which informally can be interpreted as an "instance of time" serving as "initial surface" when formulating the Einstein equations as an evolutionary system.Formally it is an achronal (i.e., no two points are chronologically related) closed (Lipschitz) hypersurface of spacetime that is hit by every inextendible causal curve.A central fact is that a spacetime admits a Cauchy surface iff it is globally hyperbolic, a statement which connects causal properties of spacetime with PDE-theory.Finally, we need to consider the boundary of the Cauchy development, called the Cauchy horizon.More precisely, this notion is defined in causal terms as follows.The future Cauchy horizon of an achronal set A is given by (3.11) It is now a basic fact that H + (S) = I + (S)∩ ∂D + (S), and with these preparations we are now ready to proceed from the pattern theorem to the "real" theorems.

The three classical theorems
The first one of the classical statements we wish to discuss is the one by Hawking, which technically is the easiest, and we will be able to provide a sketch of its proof along the above discussion.In the literature one actually finds several (versions of) singularity theorems associated with the name of Hawking and we wish to discuss the one that supposes the existence of a compact Cauchy surface.It is, in particular, applicable to a spatially closed universe and hence gives evidence for a big bang in such models.It has, however, become a custom12 to formulate it in a time-reversed manner, predicting a future singularity.
Theorem 3.8 (Hawking).Let (M, g) be a spacetime such that (E) Ric(X, X) ≥ 0 for all timelike vectors X, i.e., (SEC) holds, (C) there is a compact spacelike Cauchy surface S in M , with (I) everywhere positive future convergence k.Then M is future timelike geodesically incomplete.Sketch of proof.By compactness there is a positive minimum k 0 of k on S. So by Proposition 3.4 every timelike geodesic starting normally from S encounters a focal point at t = 1/k 0 the latest.On the other hand by Proposition 3.7 every p ∈ D + (S) \ S is reached by a normal timelike geodesic without focal point and therefore If a future directed timelike curve starting in S were to leave D + (S) it had to pass through ∂D + (S) and hence the Cauchy horizon H + (S) would be non-empty.However, being a Cauchy surface, S has empty Cauchy horizon and so I + (S) ⊆ D + (S).This, in particular, forces the above geodesics to have finite length, which makes them incomplete.
Next we turn to the theorem of Penrose.As already indicated above, its initial condition is the existence of a trapped surface, but the theorem has been generalised to trapped submanifolds of arbitrary codimension in [GS10, Thm.1].We will, however, include this advancement explicitly only in in our statement of the Hawking-Penrose theorem.
At the heart of the proof of the Penrose theorem lies the following fact: The topological condition (existence of a trapped set) that is derived from the existence of a trapped surface, contradicts the existence of a non-compact Cauchy surface.The latter condition physically amounts to the fact that we are considering an isolated system and, in particular, the gravitational collapse of an isolated body.
Theorem 3.9 (Penrose).Let (M, g) be a spacetime such that (E) Ric(X, X) ≥ 0 for all null vectors X, i.e., (NEC) holds, (C) there is a non-compact Cauchy surface S, and (I) there is an achronal closed future trapped surface P .Then M is future null geodesically incomplete.Sketch of proof.Indirectly assuming completeness, we first establish that E + (P ) is a non-empty, compact, and achronal topological hypersurface: By Proposition 3.6 E + (P ) is compact and global hyperbolicity implies that E + (P ) = ∂J + (P ).The latter is nonempty (by achronality of P ) and being the boundary of a so called future-set (i.e., a set containing its own chronological future) it is a C 0 -hypersurface (see also below).Now we take any timelike vector field Y on M (whose existence is guaranteed by time-orientability) and define the map ρ : M → S, taking each p to the unique intersection point of the maximal integral curve of Y though p with S. Obviously ρ leaves S invariant, and one may show that it is continuous and open.The restriction of ρ to ∂J + (P ) is injective by achronality and open between topological hypersurfaces and so by invariance of domain ρ(∂J + (P )) is open.On the other hand by compactness it is also closed and so ρ(∂J + (P )) = S.But since the latter set is non-compact, we have reached a contradiction.
For the final part of this section we turn to the most sophisticated of the classical singularity theorems, namely to the one by Hawking and Penrose [HP70].It collects the various developments and variants of singularity theorems that have appeared in the years before and recovers most of them under much weaker assumptions.In particular, the assumption of global hyperbolicity is avoided throughout.Consequently the necessary focusing results require a deeper analysis of the influences of curvature on causal geodesics and we will discuss this issue in some detail.In particular, a new assumption called genericity condition is introduced, which guarantees that all causal geodesics γ, at least at one point γ(t) "feel" the effect of curvature in the sense that at γ(t) the tidal force operator is nontrivial, for more details see [BEE96, Sec.

2.5].
Also the Hawking-Penrose theorem features a third possible initial condition: a (future) trapped point13 , i.e., a point p such that the expansion becomes negative for any (future directed) null geodesic starting in p.Also the initial condition of the Hawking theorem is generalised from the existence of a compact Cauchy surface to merely a compact partial Cauchy surface, i.e., an achronal set without edge.Here the edge of an achronal set A consists of all points p ∈ A possessing arbitrarily close pairs of points x ∈ I − (p), y ∈ I + (p) that can be connected by a timelike curve not intersecting A. An achronal set can be shown to be a closed topological hypersurface iff its edge is empty.This also implies that the boundary of a future set is a closed achronal topological hypersurface-a result already used in the proof of the Penrose theorem.
Later Dennis Gannon [Gan75] and Charles Walter14 Lee [Lee76] independently established results under yet another kind of initial condition: the occurrence of some nontrivial topology in a compact region of a (partial) Cauchy surface.Hence this body of results-often termed the Gannon-Lee theorems-establish a firm link between the topological and the singularity structure of spacetime.
Another more technical innovation of the Hawking-Penrose "singularity theorem par excellence"15 is that the causal part of the argument was outsourced to the following separate statement16 : Lemma 3.10 (Hawking & Penrose).In any spacetime (M, g) the following three statements cannot simultaneously hold: (C1) M is chronological, i.e., it contains no closed timelike curve.
(C2) Every inextendible causal geodesic in M contains a pair of conjugate points.
(C3) There is trapped set A.
This Lemma goes well beyond the causality arguments presented in this review so far and we will not attempt to sketch a proof here.Note that the causality condition (C1) is way more general than global hyperbolicity which was used in the above theorems.Nevertheless, it can be generalised somewhat further, see the discussion in [Sen98,p. 793].
We now give the analytic result in an extended form which, as already announced, includes yet another class of initial conditions, namely trapped submanifolds of codimension 2 < m < n, which is due to [GS10, Thm.3]17 .Theorem 3.11 (Hawking & Penrose).Let (M, g) be a spacetime such that (E) the strong energy condition (SEC) holds as well as the genericity condition along any causal geodesic γ, i.e., there is a point γ(t 0 ) such that (C) it is chronological.Moreover, assume it contains at least one of the following: (I1) a compact achronal set without edge, (I2) a closed future trapped surface P , (I3) a closed future trapped submanifold P of co-dimension 2 < m < n such that additionally for any future directed null geodesic with γ(0) orthogonal to P , or (I4) a future trapped point, i.e., p ∈ M such that on every future directed null geodesic from p the expansion θ becomes negative.
Then M is causal geodesically incomplete.
Observe that we have not stated future causal incompleteness, which, however, holds in the respective future cases of the assumptions, see [Sen98,p. 792,bottom] for details.
The line of arguments proving Theorem 3.11 from Lemma 3.10 then is as follows: Obviously the causality condition (C) and (C1) agree.The genericity condition (GC) and (SEC) are used to show that: (1) every inextendible causal geodesic in M contains a pair of conjugate points, i.e., (C2) holds, and (2) any of the initial conditions (I1)-(I4) imply that there is a trapped set, i.e., that (C3) holds.
Let us elaborate somewhat on these two items.Starting with (1), the trick is done via the following advanced focusing result: Proposition 3.12 (Focal points from genericity).Let γ be a complete causal geodesic and assume (GC) and Ric( γ(t), γ(t)) ≥ 0 for all t.Then γ has a pair of conjugate points.
Observe that this result is much stronger than Proposition 3.3 since we do not assume the existence of a point with negative expansion.This is actually why a more detailed analysis of the influence of the curvature on causal geodesics is needed and the full matrix Riccati equation (3.5) has to be used rather than just its trace, the (vorticity free) Raychaudhuri equation (3.6).The key step in the proof, which is long and technical (cf.e.g.[BEE96, p. 436-443]), is to establish that along any geodesic as in the statement with the tidal force operator nontrivial at some γ(t 0 ) (which exists due to (GC)) we have
Proof of (3.14).First note that (3.6) and (SEC) give θ(t) ≤ 0 for all t and so θ(t) ≤ 0 for all t ≥ t 0 .If θ(t) < 0 for some t ≥ t 0 then Proposition 3.3 gives the claim.So assume that θ(t) = 0 for all t ≥ t 0 implying that also θ(t) vanishes there.Inserting again into (3.6) and using (SEC) we obtain 0 ≤ −tr(σ 2 (t)) for all t ≥ t 0 .But tr(σ 2 (t)) is non-negative, and so it has to vanish for t ≥ t 0 .Consequently for t ≥ t 0 the shear σ = 1/2([B] + [B † ]) being self-adjoint, has to vanish itself and since [B] is self-adjoint as well, it too has to vanish.But this implies by the Riccati equation (3.5) that the tidal force operator is trivial at t 0 , which contradicts our initial assumption.
We now turn to the discussion of (2) and begin with the following extension of Proposition 3.6: Proposition 3.13 (Comapact horismos).Suppose (NEC) and the existence of either (i) a closed future trapped surface P , or (ii) a closed future trapped submanifold P as in (I3), or (iii) a future trapped point P .Then E + (P ) is compact, or the spacetime is null incomplete.
The proof uses again the focusing argument Proposition 3.3 to establish the main point, i.e., that E + (P ) is contained in the compact subset exp([0, T ]K) for some compact K. Now given (I4) condition (C3) follows immediately.In the cases (I2), (I3) one is only almost there, since in general P need not be achronal.The trick here is to establish that S = E + (P )∩P is achronal with compact horismos, hence a trapped set.This can be done under the additional assumption of strong causality, which, however in chronological spacetimes is a consequence of (SEC), (GC), and null completeness.
Finally in case (I1) one shows that a compact achronal set P without edge is a topological hypersurface with E + (P ) = P which, again gives (C3).This finishes our sketch of the arguments that establish the Hawking-Penrose theorem and also our journey into the classical singularity theorems.We next turn to their low-regularity extensions.

Low regularity singularity theorems
In this section we wish to provide an overview of the extension of the classical singularity theorems to metrics of low regularity that have emerged over the last couple of years.Indeed the three key theorems discussed above have first been generalised to Lorentzian metrics of regularity C 1,1 [KSSV15, KSV15, GGKS18] and then in a further effort to C 1 -metrics [Gra20, KOSS22] with a Gannon-Lee theorem proved in [SS21].We shall review the main mathematical advances that were developed to arrive at these results, mainly concentrating on the analytic side of the arguments, but we shall also comment on the recent advancements of the causality parts of the theorems [Min19a].The corresponding extensions of the singularity theorems in a purely causal setting, as well as in a synthetic setting [GKS19] and those using methods from optimal transport [CM20] will be briefly described in the final section 5.
To begin with, we discuss the motivation behind this endeavor.

Why low regularity
Taking up the discussion from Section 2.3 we now take a closer look at the conclusions of the singularity theorems, which-despite all their power and their glory-are generally considered to be a weak spot, cf.e.g.[SG15, Sec.5.1.5].Indeed, they assert merely the existence of incomplete causal geodesics and, in general, there is no way to link such a singularity to curvature blow-up in a suitable sense.Also there is the issue of extensions of spacetime: As we have seen, it is essential for the physical interpretation of the theorems that they are applied to (maximally) extended spacetimes 18 .
Taking a step back we also see that there is a regularity assumption which is implicit in the classical theorems.In fact, the theorems assert causal geodesic incompleteness of the spacetime provided the metric is smooth, and, since the bulk of Lorentzian geometry remains valid there, if it is C 2 .In particular, they do not exclude the possibility that the spacetime is complete but of lower regularity.
If the regularity of the metric was just C 1,1 the curvature would become discontinuous rather than unbounded hence would hardly be considered 'singular' on physical grounds.Indeed, via the field equations, this just corresponds to a finite jump in the matter variables.And there are many interesting systems of that type, such as the Oppenheimer-Snyder model of a collapsing star [OS39], to give a classical example, and general matched spacetimes, see e.g.[Lic55,Isr66,MS93].
Moreover, if the regularity was even lower, one could be inclined to accept such a scenario as long as there is an analytic way to define the curvature and to make sense of the field equations (2.3).This is indeed possible (in a stable way) for the Geroch-Traschen (or GT-)class of metrics [GT87], i.e., metrics of regularity H 2 loc ∩ L ∞ loc that are uniformly nondegenerate in a suitable sense [LM07,SV09].There the Riemann tensor is a tensor distribution and if a sequence of metrics converges to a GT-regular one in H 2 loc then the respective curvatures converge in distributions (for some more details see Section 4.3 below).While current techniques certainly do not carry that far, it would still be interesting to prove singularity theorems in the regularity classes used in the classical existence results for the Einstein equations which is H 5/2+ε or those used in current formulations of the strong cosmic censorship conjecture [Chr09], which demand a locally square integrable connection.However, a first substantial conceptual problem arises with the very notion of geodesic incompleteness: Below C 1,1 the initial value problem for the geodesic equation is no longer uniquely solvable and below C 1 not even classically meaningful.Thus one would have to resort to some non-classical solution concept, like e.g.Filippov solutions [Fil88], which have been used in this wider context e.g. in [Ste14,LAC22].
A first step in this direction is to lower the differentiability of the metric to C 1 for which the curvature is a distribution of order one (again see Section 4.3 below for details) and the i.v.p. for the geodesic equation is at least classically solvable if not uniquely so.In fact, we will discuss results in C 1,1 and C 1 -regularity below.
Of course, the regularity issue connected with the singularity theorems was already noted early on and extensively discussed in [HE73, Sec.8.4].There the authors argue that at least the Hawking theorem should continue to hold for C 1,1 -metrics and express their expectation that this should also be true for C 0,1 (i.e.locally Lipschitz continuous) metrics and also for all the other classical theorems.In fact, a natural next step seems to be to extend the recent C 1 -results to this class and current research is directed at this goal, see also Section 5.
Let us finally come back to the issue of extensions of spacetime.Certainly the classical theorems assert that the incomplete spacetime cannot be extended to a complete one keeping the assumptions and the C 2 -regularity of the metric.Likewise the results in C 1,1 and C 1 can be read as obstructions to such extensions keeping the respective regularity of the spacetime.This point of view nicely complements recent work by Jan Sbierski [Sbi18] who showed that the Schwarzschild solution cannot be extended as a continuous19 spacetime.In a similar vein, it has been established in [GLS18] that timelike geodesic completeness remains an obstruction to extendability also in the class of C 0 -spacetimes.

Low regularity: issues and strategies
To begin our technical account, note that the bulk of Lorentzian geometry remains valid for C 2 -metrics, since the main tools such as normal and convex neighbourhoods as well as normal coordinates are still available.However, slightly below, things begin to worsen gradually.First, the exponential map retains maximal regularity also for C 1,1 -metrics, being a bi-Lipschitz homeomorphism [Min15,KSS14].While this secures the existence of normal and convex neighbourhoods, normal coordinates are of limited use since there the metric is only continuous.Moreover, for metrics of Hölder regularity C 1,α for any α < 1, convexity breaks down completely, since the exponential map needs not be injective on any zero-neighbourhood of the origin as is demonstrated by a Riemannian example in [HW51], which easily gives rise to a static Lorentzian example, see [SS18].There it is also shown that, in general, causal geodesics fail to locally maximise the Lorentzian distance.
Returning to the context of the singularity theorems, let us briefly collect the issues at hand and discuss the strategies we will employ in the sections to come.For a longer and more technically detailed list of places where the classical proofs rest on the C 2 -differentiability see [Sen98, Ch. 6.1].
(A) The curvature is only locally bounded for g ∈ C 1,1 and merely a distribution of order one for g ∈ C 1 .
(B) Normal neighbourhoods are not useful for g ∈ C 1,1 and the exponential map is no longer even defined for g ∈ C 1 .
(C) The geodesic equation fails to be uniquely solvable for g ∈ C 1 .
Item (A) and, more precisely, the failure of the curvature to be defined pointwise causes several difficulties.First the energy conditions (E) have to be adapted, where especially (NEC) and (GC) turn out to be a delicate matter as we will discuss in detail below.Then, Jacobi fields cannot be reasonably defined which, of course, means that one has to do without using the central concepts of conjugate and focal points.A clear strategy to address this issue-as was already pointed out in [HE73, Sec.8.3]-is regularisation of the metric by smooth approximations for which the classical tools are still available.We will generally pursue this path and explain it in some detail in the following Section 4.3.Then in Section 4.4 and in Section 4.7 we will see how appropriate distributional versions of the energy conditions (E) lead to the focusing of causal geodesics for the approximating metrics.
Next, item (B) makes it necessary to revisit the whole machinery of causality theory and to extended it to the regularity at hand.We will refrain from going into any technical details here and just summarise the necessary background in Section 4.5 Now, (C) first of all forces us to make a choice concerning the conclusion of the theorems, namely on the notion of incompleteness.It has turned out to be favourable to use the more stringent alternative20 put forward in [Gra20]: Definition 4.1 (C 1 -completeness).A C 1 -spacetime is called timelike (respectively null or causal) geodesically complete if all inextendible timelike (respectively null or causal) solutions of the geodesic equation are defined on all of R.However, despite the loss of uniqueness many properties of geodesics extend from the smooth to the C 1 -setting.In particular, geodesics do have a fixed causal character [Gra20, Cor.2.3] and so the above definition makes sense.
One main issue will be to adjust the following vital aspect of the regularisation approach: We need to approximate maximising causal geodesics of the low regularity metric by maximising causal geodesics of the approximating metrics, which becomes a delicate matter in the absence of an exponential map (B) and unique solvability of the geodesic equation (C).This issue was addressed for globally hyperbolic C 1spacetimes in [Gra20, Sec.2] by establishing that between any pair of points p < q there is (at least) one maximising causal geodesics that can be suitably approximated by a sequence of maximising causal geodesics of approximating metrics.This and corresponding results for points in the chronological future of a Cauchy surface and in the horismos of a closed spacelike codimension 2 surface, allow for corresponding proofs of the Hawking and the Penrose theorem, which we will sketch in Section 4.6.However, in the context of the Hawking-Penrose theorem more general results are needed, which will ultimately force us to introduce a new condition, namely a nonbranching assumption for maximising causal geodesics, to be detailed in Section 4.8.This will finally allow us to discuss the recent C 1 -version of the Hawking-Penrose theorem [KOSS22] in Sections 4.9 and 4.10.

Distributional curvature & regularisation
Here we briefly discuss the general distributional framework in which to understand the curvature of metrics 21 of regularity below C 2 , with a special emphasis on the C 1 -case.For more information consult [Mar68, GT87, LM07, GKOS01, Ste08].
The space of (scalar) distributions of order k on M is the topological dual of the space of compactly supported C k -one-densities Γ k c (M, Vol(M )) (here Vol(M ) denotes the volume bundle), i.e., Similarly, the space of distributional (r, s)-tensor fields of order k is defined as Here T r s (M ) denotes the space of smooth tensor fields of rank (r, s), i.e., T r s (M ) = Γ ∞ (M, T r s ) = Γ(M, T r s ), since we generally omit k if it is infinite.For the spaces of vector fields, one forms and tensor fields of finite differentiability we will write X C k , Ω 1 C k , and (T r s ) C k , respectively.With this notation we have 22 The first line says that distributional tensor fields are sections of the corresponding tensor bundle with distributional coefficients 23 , and the second line reveals them as C k multilinear maps on one forms and vector fields of regularity C k that give scalar distributions of order k.This fact for k = 1 will turn out to be essential in formulating the genericity condition for C 1 -metrics below.
A distributional connection is a map ∇ : X(M ) × X(M ) → D ′ T 1 0 (M ) satisfying the usual computational rules: All that is said here also applies to metrics of arbitrary signature.The isomorphisms in (4.2) are algebraic and bornological but not topological, cf.[Nig13].Here ⊗ C k (M) denotes the balanced tensor product over the module Denoting by G any of the spaces C k (0 ≤ k) or L p loc (1 ≤ p), we call a distributional connection a G-connection, if ∇ X Y is a G-vector field for any X, Y ∈ X(M ).A particularly important case are L 2 loc -connections since they allow to define the curvature tensor R : Moreover, if E i is a local frame in X(U ) and E j ∈ Ω 1 (U ) is its dual frame, then the Ricci tensor is given by (4.4) The significance of the GT-regular metrics is then rephrased by saying that their Levi-Civita connections actually are L 2 loc -connections.More specifically, metrics g ∈ C 1 have C 0 -Levi-Civita connections which implies that their Riemann tensor R ∈ D ′(1) T 1 3 (M ) as well as their Ricci tensor and scalar curvature are of order 1.Finally, the standard local formulae hold in D ′(1) and in (4.4) we even can use E i of regularity C 1 , which is of great technical importance since it makes it possible to use frames derived via parallel transport w.r.t. a C 1 -metric.
As already indicated above our method of choice to deal with the analytic arguments in the proofs of the singularity theorems is regularisation and we will detail our convolution based approach below.Prior, we outline our overall strategy to implement an analytic machinery that forces causal geodesics of the rough metric to stop maximising: We will formulate suitable energy conditions for the low regularity metric g ∈ C 1 and derive from it surrogate energy conditions for a sequence of smooth approximating metrics g ε .These will be weaker than the classical conditions (SEC) and (NEC) in so far as the corresponding expressions24 Ric[g ε ](X, X) will be shown to be (only) mildly negative.We will then extend the arguments explained in Sections 3.2, 3.5 in order to still prove the occurrence of focal/conjugate points along causal g ε -geodesics.(Note, that due to (A) we cannot resort to standard results using averaged energy conditions as put forward e.g. in [FG11].)This will eventually force the geodesics of g itself to stop maximising.To achieve these goals we have to take the following into account: (R1) When deriving the surrogate energy conditions for g ε from conditions on g we face the problem that while Ric[g ε ] → Ric[g] distributionally, we cannot achieve local uniform convergence even for g ∈ C 1,1 .Therefore we will have to compare Ric[g ε ] to a regularisation of Ric[g] instead.
(R2) Since we want to use the classical arguments on smooth approximations g ε to derive focusing results, we have to control their causality in terms of the causality of g.This can be done thanks to an adapted regularisation procedure put forward by Chrusciel and Grant [CG12], that provides us with approximations ǧε and ĝε with narrower and wider lightcones25 than g, respectively.
(R3) We will have to relate the geodesics of the approximating metrics to the geodesics of g.More precisely we will have to show that maximising causal g-geodesics are C 1 -limits of g ε -causal g ε -maximising g ε -geodesics.
We proceed by introducing the regularisations to be used in detail.The basic ingredient is chartwise convolution and we begin by choosing a mollifier, i.e., a smooth, nonnegative function ρ on R n , supported in the unit ball, and with unit integral.Then we cover M by a countable and locally finite family of relatively compact chart neighbourhoods (U i , ψ i ) (i ∈ N) and let (ζ i ) i be a subordinate partition of unity with supp(ζ i ) ⊆ U i .Then we choose a family of cut-off functions χ i ∈ C ∞ c (U i ) with χ i ≡ 1 on a neighbourhood of supp(ζ i ).Finally, for ε ∈ (0, 1] we set ρ ε (x) := ε −n ρ x ε .Then, denoting the push-forward and pull-back of distributions with upper and lower stars, consider for any T ∈ D ′ T r s (M ) the expression Here, ψ i * (ζ i • T ) is viewed as a compactly supported distributional tensor field on R n , so componentwise convolution with ρ ε yields a smooth field on R n .The cut-off functions χ i secure that (ε, x) → T ⋆ M ρ ε (x) is a smooth map on (0, 1] × M .For any compact set and small ε, equation (4.5) reduces to a finite sum with all χ i ≡ 1, hence to be omitted from the formula.Just as is the case for smoothing via convolution in the local setting we obtain optimal convergence, that is 26 .In particular, for g ∈ C 1 we now set to obtain a sequence (actually a net) of smooth Lorentzian metrics that converges 26 in C 1 to g.It is now essential for our purposes to tweak this construction to obtain a regularisation adapted to the causality as indicated above in (R2).We will use the version of [Gra20, Lem.4.2, Cor.4.3]: Lemma 4.2 (Regularisations and convergence).Let (M, g) be a C 1 -spacetime.Then for any ε > 0 there exist smooth Lorentzian metrics ǧε , ĝε on M satisfying ǧε ≺ g ≺ ĝε for all ε and ǧε , ĝε → g in C 1 (ε → 0).(4.7) Moreover we have control on the speed of convergence of g ε and compatibility between ǧε and g ε as follows: For any compact K there is c K > 0 such that for small enough An analogous statement holds for ĝε as well as for the inverse metrics g −1 , g −1 ε , (ǧ ε ) −1 , and (ĝ ε ) −1 .
With these preparations we may now have a look at the distributional energy conditions for metrics g ∈ C 1 and how they imply focusing of causal geodesics for approximating metrics.

Distributional energy conditions & focusing
In this section we introduce the energy conditions akin to (SEC) and (NEC) for C 1metrics and show that they imply suitable surrogate energy conditions on the approximations, cf.(R1) above.Here "suitable" means that we can still use the surrogate conditions to show focusing results (cf.(R2)) for the geodesics of the approximating metrics-despite the fact that they are weaker than the classical conditions, which are manifestly violated, but only by a controlled margin.
We will be more detailed in case of the (SEC) where we point out the main ideas in the proofs, and will be more sketchy in case of the technically more demanding (NEC).
To begin with, recall that a scalar distribution u ∈ D ′ (M ) is nonnegative, u ≥ 0, if u(ω) ≥ 0 for all nonnegative test densities ω ∈ Γ c (M, Vol(M )).A nonnegative distribution is always a measure and hence a distribution of order 0.Moreover, nonnegativity is stable with respect to regularisation,27 for details see [H 03, Thms.2.1.7 and 2.1.9].Finally, for u, v We say that a C 1 -spacetime satisfies the distributional strong energy condition The new condition (DSEC) is compatible with the usual classical conditions as well with their obvious reformulation in L ∞ for g ∈ C 1,1 .The surrogate energy condition for the approximation now is, cf.[Gra20, Lem.4.6]: Lemma 4.3 (Surrogate (SEC)).Let (M, g) be a C 1 -spacetime satisfying (DSEC) and let K ⊆ M be compact.Then for all δ > 0 and all smooth X ∈ T M | K with ǧε (X, X) = −1 we have for small enough ε.
Sektch of proof.As already remarked above in (R1) the main problem is that the convergence of Ric[ǧ ε ] to Ric[g] is not good enough to directly carry the positivity of Ric[g](X, X) through the argument.Rather we proceed as follows: By standard properties of the convolution (DSEC) implies Ric[g]⋆ M ρ ε (X, X) > 0 and the result will follow from the compatibility of the distinct regularisations we have used.Indeed we are done, if we can show that Ric To establish these statements, we have to estimate the convolution of a product a f (basically the components of the inverse metric times a derivative of the components of the metric as occurring in the Christoffel symbols) to a corresponding product of convolutions.After reducing everything to the local situation, this is done by the following Friedrichs-type lemma which takes as an essential input the final estimate in Lemma 4.2, i.e., (4.8) for the inverse metrics.Indeed, the components of the inverse of the regularised metric precisely possess the convergence properties assumed for a ε in the Lemma below, cf.[Gra20, Lem.4.9].
This focusing result will play a mayor role in the proof of the C 1 -Hawking theorem.Let us now turn to the corresponding focusing result for the C 1 -Penrose theorem.First a proper distributional formulation of (NEC) needs some more care due to the fact that vectors that are ǧε -null are only almost g-null, cf.[Gra20, Sec.5].
We say that a C 1 -spacetime satisfies the distributional null energy condition (DNEC), if for any compact set K and any δ > 0 there exists for any local smooth vector field X ∈ X(U ) (U ⊆ K open) with X h = 1 and | g(X, X) |< ε on U .Here h denotes the norm with respect to some complete Riemannian background metric h.Note that this and all future such conditions and estimates will be local and hence in fact be independent of the choice of h.
Again the new condition (DNEC) is compatible with both (NEC) in the smooth case, and the almost everywhere condition used in the C 1,1 -case.The following analogue of Lemma 4.3 shows that the above definition of the null energy condition is the correct one in the sense that it produces the following surrogate energy condition on the level of approximations, see [Gra20, Lem.5.5]: Lemma 4.6 (Surrogate (NEC)).Let (M, g) be a C 1 -spacetime satisfying (DNEC).Let K ⊆ M be compact and let c 1 , c 2 > 0. Then for all δ > 0 there is Now the following result is both the generalisation of the arguments that lead to smooth focusing in Propositions 3.3 and 3.6 to the C 1 -setting, and the null version of Proposition 4.5.It says that negative expansion still leads to focusing even under (SNEC).We formulate it in a slightly different and quantified manner which assumes the geodesic to be maximising and then restricts its length, for details see [Gra20, Lem.5.6].
Proposition 4.7 (Smooth focusing under (SNEC)).Let g be smooth and let γ : [0, β] → M be a future directed null geodesic starting at p from a spacelike submanifold P of codimension 2. If γ is maximising the distance to P (and hence starts normal to P ) we have: We will return to the issue of focusing in Section 4.7 where we will also introduce an appropriate distributional version of the genericity condition.First we discuss the C 1 -versions of the first two of the classical theorems in Section 4.6.But to do so, we have to briefly turn to causality theory in low regularity in the next section.

A brief word on causality theory
The regularisation of Lorentzian metrics with controlled causality (cf.Lemma 4.2) put forward in the seminal paper [CG12] was in turn used to study causality theory of continuous metrics.In fact, this paper together with [FS12] initiated the recent systematic study of causality theory in low regularity, see [Min19a, Sec.1] for a brief overview.
In this way the bulk of Lorentzian causality theory has been transferred to C 1,1spacetimes.While convexity fails below that regularity (cf.Section 4.2) nevertheless most aspects of causality theory can be maintained even under Lipschitz regularity of the metric.Further below some significant changes occur [CG12, GKSS20], while some robust features continue to hold even in more general settings [Min19a, KS18, BS18, GKS19], see also Section 5 below.
In particular, for C 1 -spacetimes we may still build the causality relations on local Lipschitz curves28 .Also the push-up principle is still valid and I + (A) is open for any set A ⊆ M .Moreover, even for continuous metrics one may consistently define global hyperbolicity via causality and compactness of causal diamonds J(p, q) = J + (p) ∩ J − (q), [Min19a, Prop.2.20], which is still equivalent to the existence of a Cauchy surface, see [Säm16, Thms.5.7 and 5.9].Also it is then clear that if g is globally hyperbolic, so is ǧε .
From there all further ingredients needed in our approach to the singularity theorems in C 1 can be derived in C 1 -regularity, see [KOSS22, Appendix A].These results, however, also follow by compatibility from the more general approaches of [CG12,Säm16] and, in particular [Min19a], where cone structures on manifolds are studied and many results are derived under minimal regularity assumptions.We will briefly return to this issue in the final Section 5.

The Hawking and the Penrose theorems in C 1
We are now finally ready to present the extensions of the first two of the three classical singularity theorems we have been dealing with throughout.We start with the Hawking theorem which was given for C 1,1 -metrics in [KSSV15, Thm.1.1] and generalised to C 1 in [Gra20, Thm.4.11] as follows: Theorem 4.8 (C 1 -Hawking).Let (M, g) be a C 1 -spacetime such that (E) the timelike Ricci curvature is nonnegative in D ′ , i.e., (DSEC) holds, (C) there is a compact spacelike Cauchy surface S in M , with (I) everywhere positive future convergence k.Then M is future timelike geodesically incomplete.Sektch of proof.First we assume by contradiction that (M, g) is timelike geodesically complete.Then, by compactness of S and continuity of the mean curvature k := (n − 1) min S H, ν (with H = H[g] and ν = ν[g] the future directed g-unit normal) exists and is positive.
Next we want to apply Proposition 4.5 to these geodesics in a uniform way.In order to do so we need to apply (SSEC) with a uniform constant δ and hence we have to make sure that (the images of) all the ǧε -geodesics from above stay in one compact set (at least for small ε), cf.Lemma 4.3.To explain in some detail how this is done we introduce the following notation: For any compact K ⊆ T S and any parameter value t we write where the union runs over all g-geodesics with γ(0) ∈ K. Similarly we write F ε,K,t if the γ-geodesics are replaced by ǧε -geodesics with data in K. Now we consider and set b = 4(n − 1)/k 0 .Then (cf.[Gra20, Prop.2.9]) there is ε 0 (K, b) such that is a relatively compact subset of T M .Denote by L a compact neighbourhood of F ≤ε 0 ,K,b .Then setting δ 0 = k 2 0 /(2(n − 1)), Lemma 4.3 gives us ε 1 (L, δ 0 ) such that (SSEC) holds on L with δ 0 and for all ε ≤ ε 1 .Therefore, for all ε ≤ min(ε 0 , ε 1 ) we may apply Proposition 4.5 to conclude that the respective ǧε -geodesics stop maximising the ǧε -distance to S the latest at b 0 = 2(n − 1)/k 0 < b.Therefore, as in the proof of the classical Hawking theorem 3.8, we have that for ε small D + [ǧ ε ](S) is contained in a compact set.Recall that S is a Cauchy surface also for ǧε .But the Cauchy horizon Note that we have neither shown nor needed that also D + (S) ⊆ {p ∈ M : d(S, p) ≤ 1/k}.This fact, however, can be established (cf. the discussion preceding [Gra20, Thm.4.13]) using that also for globally hyperbolic g ∈ C 1 maximising geodesics exist.More precisely, the Avez-Seifert theorem was established even for g ∈ C 0 in [Säm16], establishing the existence of maximising causal curves between any pair of points p < q.In C 1,1 -and in C 1 -spacetimes these maximisers are also geodesics by [Min15,Thm. 6] and [SS21, Thm.3.3], respectively, but see [LAC22] for more general results.Hence we have existence of maximising geodesics but as already indicated at the end of Section 4.2, we also need to approximate these by a sequence of ǧε -maximising ǧε -geodesics, cf.(R3).In the present globally hyperbolic C 1 -setting this was established for any pair of points p < q and between a Cauchy surface S and any p ∈ I + (S) in [Gra20, Prop.2.12, Cor.2.14] (even before [SS21, Thm.3.3] was available).
For the proof of the Penrose theorem we will need the following version for null maximisers from (n − 2)-surfaces given in [Gra20, Cor.2.15]: Proposition 4.9.Let (M, g) be a globally hyperbolic C 1 -spacetime and let P ⊆ M be a closed, spacelike (n − 2)-dimensional submanifold.Then for any q ∈ E + (P ) there exists at least one null geodesic from P to q maximising the distance to P .Further, such a geodesic can be obtained as a C 1 -limit of a sequence of ǧεn -null ǧεn -geodesics γ εn maximising the ǧεn -distance to P .
We now proceed to the Penrose theorem which was generalised to C 1,1 -metrics in [KSV15, Thm.1.1] and further to the C 1 -setting in [Gra20, Thm.5.7] as follows: Theorem 4.10 (C 1 -Penrose).Let (M, g) be a C 1 -spacetime such that (E) the null Ricci curvature is nonegative in D ′ , i.e., (DNEC) holds, (C) there is a non-compact Cauchy surface S, and (I) there is an achronal future trapped surface P .Then M is future null geodesically incomplete.Sketch of proof.We assume null completeness and first establish that the future horismos E + (P ) of P is compact.Since P is trapped, H, ν > 0 for all future pointing null normals ν.Defining the compact set K = {X ∈ T P ⊥ : 0 ≤ H, X ≤ 2} we show that E + (P ) ⊆ F K,1 which is compact.This gives the claim since by global hyperbolicity E + (P ) is also closed.So let's assume by contradiction that there is p ∈ E + (P )\F K,1 .Since p ∈ E + (P ), by Proposition 4.9 it is reached by a normal maximising null geodesic γ : [0, 1] → M which is the uniform C 1 -limit of a sequence of ǧε k -null geodesics γ ε k : [0, 1] → M maximising the ǧε k -distance to P .Since p ∈ F K,1 we have H, γ(0) > 2 and so for k large −(n − 2)ǧ ε k (H ε k , γε k (0)) < −2(n − 2) =: θ 0 .By convergence we may choose a compact set and constants c 1 , c 2 to apply Lemma 4.6 to obtain (SNEC) for all large enough k with an appropriate choice of δ < −3θ 0 /4.Then Proposition 4.7 gives 1 = β ≤ − 4(n−2) 3θ 0 = 2 3 , a contradiction.So, ∂I + (P ), which by global hyperbolicity equals E + (P ), is compact and we aim at topological argument as in the classical proof.To this end we choose a smooth metric g ′ ≺ g.Then S is a non-compact Cauchy surface also for (M, g ′ ).Further, I + (P ) is a g ′ -future set by I + [g ′ ](I + (P )) ⊆ I + (I + (P )) = I + (P ).So its boundary ∂I + (P ) is a compact g ′ -achronal topological hypersurface.Now we proceed exactly as in the proof of 3.9 to obtain a homeomorphism between S and ∂I + (P ).The latter set is compact while the former is not, again a contradiction.
In the following we begin to turn towards the low regularity extension of the third of the classical theorems, namely the one of Hawking-Penrose.We start with the causal result, i.e., Lemma 3.10.Here we have, of course, to remove the explicit occurrence of conjugate points in condition (C2).However, as we have seen in our previous discussions, what really matters is that inextendible causal geodesics stop maximising the Lorentzian distance.Moreover, a closer inspection of the classical proofs reveals that in the timelike case this property is only needed on open globally hyperbolic subsets of spacetime (cf.also the discussion at the end of Section 4.8).Using the extensions of causality theory briefly discussed in Section 4.5 one may establish the following result, see [KOSS22, Thm.6.2]: (C3 ) There is trapped set A.
While in establishing the C 1 -Hawking-Penrose theorem we will follow the general layout of the classical arguments used in Section 3.5, the steps from Lemma 4.11 to the analytic result is now considerably more involved.Considering the sophisticated analysis of the influence of curvature on causal geodesics needed in the classical proof (cf.Proposition 3.12) this is not surprising at all: The technical issues arising form the lack of a suitable concept of conjugate points become more pronounced here.In fact, it will take us three more technical sections to prepare the statement and a sketch of proof which we will finally provide in Section 4.10 30 .We start our account by discussing the distributional genericity condition in the next Section 4.7.We will show that it allows to still derive a focusing result for the causal geodesics of approximating smooth metrics.Then, following our general plan, we want to establish that this leads to focusing of the causal geodesics of the C 1 -metric g.In particular, we have to establish that, under appropriate assumptions, conditions (C2') and (C2") of Lemma 4.11 hold.While this can be achieved using the methods discussed so far in the case of a C 1,1 -metric (cf.[GGKS18, Thms.5.1, 5.3]), it needs a more careful analysis of how in the C 1 -case maximising causal geodesics can be approximated by respective geodesics of the approximating metrics (R3).We will discuss this issue in Section 4.8 where we introduce a new condition that prohibits the branching of maximal causal geodesics for g.Finally, in Section 4.9 we will deal with appropriately generalised initial conditions and discuss how they each-once again using focusing-lead to the formation of a trapped set.

Distributional genericity & advanced focusing
We now wish to introduce an appropriate distributional version of the genericity condition which will allow us to prove the advanced focusing results needed to establish the C 1 -Hawking-Penrose theorem.The general strategy is the same as employed in Section 4.4: From the distributional genericity condition we will establish estimates on the tidal force operator for the approximations, which are still strong enough to produce pairs of conjugate points along causal geodesics of the approximating metrics that are "long enough".
As in the case of the (DNEC) we will have to use extensions of vector fields to small neighbourhoods.However, since in the course of our later arguments we are bound to use g-frames which possess mere C 1 -regularity we need to formulate the condition for C 1 -fields.It is here that we crucially rely on the fact that the Ricci tensor is a first order distribution and hence allows to insert C 1 -vector fields, cf.(4.2).
We say that along a causal geodesic γ the distributional genericity condition (DGC) holds at γ(t 0 ), if it possesses a neighbourhood U with C 1 -vector fields X restricting to γ, and V restricting to a vector field normal to γ at all γ(t) ∈ U and there exist c > 0 and δ > 0 such that for all C 1 -vector fields X, Ṽ on U with X − X h < δ and V − Ṽ h < δ we have This condition again is consistent with its smooth and C 1,1 -counterparts.Also it implies useful estimates on the tidal force operators of the approximations.The precise results are [KOSS22, Lem.2.16, 2.17], which we summarise next.At the heart of the technical proofs is again positivity of the direct convolution of the positive term in (DGC) and the compatibility of the various regularisations, established by the Friedrichs lemma 4.4.Lemma 4.12 (Bounds on the tidal force operator from (DGC)).Let γ be a causal geodesic in a C 1 -spacetime and assume (DGC) at γ(0).Let γε be ǧε -geodesics, whose ǧε -causal character is the same as the g-causal character of γ and which converge in C 1 to γ.Then there are vector fields E ε i on some neighbourhood U of γ(0) such that (i) E ε i • γ ε is a g ε -orthonormal frame along γ ε , and (ii) , where E i • γ is an orthonormal frame along γ.Finally, there are c > 0, r > 0, and C = C(ε) > 0 such that along γ ε the tidal force operator for small enough ε satisfies (4.11) The announced focusing result asserting the existence of a pair of conjugate points along causal geodesics of smooth metrics given an estimate of the form (4.11) is [KOSS22,Lem. 4.1]31 : Proposition 4.13 (Advanced smooth focusing).Let g be smooth.Given some c > 0 and 0 < r < π Finally, since both [ R] and [ B](t 1 ) are diagonal, the Riccati equation for [ B] decouples and can be explicitly solved by where and H −C,f is of the same form with the hyperbolic tangent replacing the cotangent.Now again using the Raychaudhuri equation for θ and explicitly analysing the function H c,f it turns out that under the sole condition 4r √ c < π, the constants δ and T (depending only on c) can be chosen small, respectively large enough to arrive at a contradiction to [B] ≤ [ B].Now, as discussed at the end of Section 4.6, we want to establish via the above Proposition that also geodesics for the C 1 -metric stop maximising if they become too long.We do so in the next section introducing some new ideas.

Geodesic branching
To motivate the introduction of the new non-branching condition we wish to sketch the essential argument that allows us to pass from Prop.4.13 to the fact that no causal geodesic of g ∈ C 1 is globally maximising under (DGC) and (DSEC) resp.(DNEC).(4.12) In fact, we will only do so for timelike geodesics and g ∈ C 1 globally hyperbolic.Assume indirectly that we are given a complete timelike geodesic γ that is maximising between any of its points.The idea is to first construct a sequence of maximising timelike ǧε -geodesic γ ε from γ(−T ) to γ(T ) for some suitably large T which converges to γ in C 1 .Then assuming (DGC) at γ(0) and (DSEC) we want to employ Proposition 4.13 to show that actually γ ε cannot be maximising.This allows to reach a contradiction, since the C 1 -limit γ of γ ε was assumed to be a maximiser.
More precisely the C 1 -convergence of the approximating geodesics will allow us to apply Lemma 4.12 to conclude from (DGC) that (4.11) holds, which actually is assumption (ii) in Proposition 4.13 for ǧε .Furthermore C 1 -convergence of γ ε will also allow us to choose an appropriate compact set in T M to turn (DSEC) into assumption (i) of Proposition 4.13 for ǧε via Lemma 4.3.Then for an appropriate choice of δ (small) and T (large) Proposition 4.13 asserts that γ ε has conjugate points and hence cannot be maximising, which gives the desired contradiction.
So everything boils down to construct a sequence of maximising timelike ǧεgeodesic γ ε from γ(−T ) to γ(T ) that converges to γ in C 1 .In case g ∈ C 1,1 , this can actually be achieved in the following way: By global hyperbolicity of ǧε there is a maximising timelike ǧε -geodesic γ ε from γ(−T ) to γ(T ).Suitably reparametrising γ and γ ε we may achieve that γ(−T ) = γ ε (−T ) and that γ(−T ) and v ε := γε (−T ) have the same h-norm bound.Consequently v ε has a subsequence v ε k converging to some w with w h = γ(−T ) h .Therefore γ ε k converges in C 1 to a g-geodesic γ w which is the unique solution to the geodesic equation with initial data γ w (−T ) = γ(−T ) and γw (−T ) = w which, as a limit of maximisers, is maximising and can be shown to reach γ(T ).Now if γ w = γ then both are maximising from γ(−T ) to γ(T ) and hence γ is not maximising beyond T , which contradicts our assumption.Hence we may assume that γ w = γ and we have constructed the desired sequence of approximating geodesics, which allows us to follow the arguments laid out above.
However, if g is merely a C 1 -metric uniqueness of solutions of the geodesic equation fails and the above construction does not produce a sequence of approximating geodesics.So we have to look for an alternative.
First, note that in the present situation, i.e., given an ODE with merely continuous r.h.s. and possibly different solutions to the same initial data, the following issue arises in any approximation approach: Regularising the coefficients of the equation results in a smooth situation with unique solvability of the initial value problem.Now to approximate a given (non-unique) solution of the original i.v.p. by (unique) solutions of the regularised problem seems unfeasible unless a specific regularisation is constructed to force exactly the desired convergence.Such an approach is of course completely unsuited to the situation at hand, where we want to approximate a given but arbitrary maximising geodesic.
Instead we will employ a non-branching condition for maximising geodesics, which is well motivated by similar conditions used in (Riemanian) metric geometry.There, in absence of a differentiable structure, geodesics are defined as (local) minimisers of the length functional and local uniqueness of geodesics is expressed in the form of a non-branching condition.Here, a branch point is defined as an element of a minimiser at which the curve splits into two minimisers that on some positive parameter interval do not have another point in common, cf., e.g., [Shi93,Vil09].Similarly, in the synthetic Lorentzian setting [KS18], the role of causal geodesics is taken on by maximising causal curves, and non-branching is formulated analogously.In both cases, lower synthetic sectional curvature bounds (formulated via triangle comparison in constant curvature model spaces) imply non-branching [Shi93,KS18].
While in the present C 1 -setting, the coincidence between causal local maximisers and geodesics, that is familiar from smooth Lorentzian geometry ceases to hold ([KOSS22, Ex. 3.2], [SS18]), it still seems reasonable to assign a privileged role to causal geodesics that are locally maximising and to preclude them from branching.We now explicitly introduce our (rather weak) non-branching conditions, where we do not require the second branch to be maximising as well: (4.13) A C 1 -spacetime is called maximally causally (resp.timelike, resp.null) non-branching (MCNB, MTNB, MNNB), if no maximal causal (resp.timelike, resp.null) geodesic branches in the above sense.
behaved with respect to approximation, cf.proof of Thm.5.3 in [GGKS18].So these have to be excluded in the statement by assuming the spacetime to be causal instead of merely chronological, as was sufficient in the classical results.

Initial conditions
In this last of the preparatory sections we generalise the four distinct initial conditions of the Hawking-Penrose Theorem to the present C 1 -setting.More precisely, we provide suitable extensions of the four initial conditions (I1)-(I4) of the classical theorem 3.11 and discuss how they each lead to the formation of a trapped set, i.e., condition (C3) of the C 1 -Hawking-Penrose Lemma 4.11.
The easiest case to deal with is the first, i.e., (I1), since it needs no special attention: Given the extension of causality theory discussed in Section 4.5 it follows just as in the smooth case (cf. the end of Section 3.5) that a compact achronal set P without edge is a topological hypersurface and E + (P ) = P , which gives (C3).
Next we deal with the trapped submanifold-cases, i.e., (I2) and (I3).In the present regularity class it is natural to extend the corresponding initial conditions to trapped C 0 -submanifolds defined in the support sense.Definition 4.17.(Closed trapped C 0 -submanifolds) A compact without boundary, spacelike C 0 -submanifold P ⊆ M of codimension 1 < m < n is called a future trapped submanifold if (1) any point p ∈ P possesses a neighbourhood U with U ∩ P achronal in U , and (2) the mean curvature vector field is past pointing timelike in the sense of support submanifolds on all of P .
Condition (2) means that for any q ∈ P there exists a future C 2 -support submanifold P for P at q whose mean curvature vector at q is past-pointing timelike.Here a future support submanifold P for P at q ∈ P is a submanifold of the same dimension containing q and such that there is a neighbourhood U of q in M such that P ∩ U ⊆ J + U (P ), the causal future of P within U .Now, to show that, given a closed trapped C 0 -submanifold P , a trapped set forms, i.e. that (C3) holds, we have to again establish that lightrays from P stop maximising, now under (DNEC) (in case m = 2) and a suitable distributional generalisation of condition (3.13) (in case of general m).This condition is very much in the spirit of (DGC), in the sense that it also asks for a C 1 -stability of the condition on the curvature: We consider a future directed normal null geodesic γ starting at a point p ∈ P .Choose a frame e 1 (p), . . ., e n−m (p) of T p P and denote by E 1 , . . .E n−m ist (C 1 -)parallel transport along γ.Now we assume that for any b in the domain of γ there is a neighbourhood U of γ | [0,b] and C 1 -extensions E i of E i and N of γ to U such that for each δ > 0 there exists η > 0 such that for all collections of C 1 -vector fields { Ẽ1 , . . ., Ẽn−m , Ñ } on U with Ẽi − E i h < η for all i and Ñ − N h < η, we have Finally we turn to (I4) and introduce trapped points in the current setting.In [GGKS18, Sec.6.3], a faithful generalisation of the classical condition to the C 1,1setting is given, motivated by Jacobi tensor classes and the mean curvature of spacelike 2-surfaces, given as the level sets of the exponential map that generate the light cone.Although these tools are no longer at our disposal we can use the very formulation which, once more, is given in the support sense: Definition 4.19 (Trapped points in C 1 ).A point p ∈ M is called future trapped, if for any future pointing null vector v ∈ T p M and for any null geodesic γ with γ(0) = p, γ(0) = v, there exists a parameter t and a spacelike C 2 -submanifold S of codimension m = 2 with S ⊆ J + (p), γ(t) ∈ S and convergence k S ( γ(t)) > 0. Now the techniques established above allow one to also prove ([KOSS22, Prop.5.13]) that in any strongly causal, null geodesically complete, MNNB C 1 -spacetime satisfying (DNEC) the horismos E + (p) is compact for any future trapped point p, and we have established (C3) once more.Now we are finally in a position to put everything together and establish the main result in our account.

The Hawking-Penrose theorem in C 1
Finally, we are in the position to formulate the long seeked result, namely the analytical Hawking-Penrose theorem for maximally causally non-branching C 1 -spacetimes.We first give the precise statement and comment on some of its specifics, and then put together the arguments that allow to derive it from the causal result given in Lemma 4.11.
Theorem 4.20 (C 1 -Hawking-Penrose).Let (M, g) be a C 1 -spacetime such that (E) the distributional strong energy condition (DSEC) holds as well as the distributional genericity condition (DGC) along any causal geodesic, (C) it is causal, and (B) it is maximally causally non-branching (MCNB).Moreover, assume it contains at least one of the following: Then M is causally geodesically incomplete.
We now briefly compare the present result with its classical counterpart.The distributional energy conditions used in (E) are faithful generalisations of the classical conditions.The causality condition (C) is slightly stronger than the condition of chronology used in the classical result and we have commented on our use of causality at the end of Section 4.8.At the moment it is not clear whether this is a mere technical point or whether causality is strictly necessary for the result to hold.The initial conditions (I1)-(I4) are again faithful generalisations of their classical counterparts and their formulation in the support sense seems only natural in the C 1 -regularity class.
Finally, we come to condition (B), which is entirely new.As discussed in Section 4.8 it is, on the one hand, necessary to secure a main argument in the proof of the the advanced focusing result (4.12), which is the fundamental ingredient to establish the inexistence of lines under (DSEC)/(DNEC) and (DGC).On the other hand, it is well motivated by similar conditions used in metric geometry, see also Section 5, below.Moreover, it adds a novel aspect to the interpretation of the C 1 -theorem: Under the given conditions the result predicts either geodesic incompleteness or branching of maximising causal geodesics.The latter alternative physically signifies an event equally catastrophic for the corresponding observer or light ray: instead of suddenly beginning or ending its existence it splits in two or, in the past case, two observers/light rays are merged into one.
Of course, there is still the alternative that the regularity of the metric drops even further, i.e., below C 1 , which renders the curvature a distribution of higher order.Again, the result also forbids the extension of the spacetime to a complete one of regularity C 1 without (maximal causal) geodesic branching and we shall briefly return to this discussion in the concluding Section 5.

Sektch of proof.
Assuming causal geodesic completeness, we once again deduce a contradiction, now using Lemma 4.11.
First note that (C) implies (C1) there.Then Theorem 4.16 implies the assumptions (C2') and (C2") of Theorem 4.11 and it remains to establish the existence of a trapped set, i.e. (C3) for which we have to distinguish the four different cases.
The case for (I1) has been made at the beginning of Section 4.8.For the remaining cases we need-just as in the classical proof of Theorem 3.11-an "upgrade" of (C) to strong causality.In the present setting this can be achieved from the fact that there are no null lines, cf.[KOSS22, Lem.A.31].Then Proposition 4.18 covers the cases (I2) and (I3), while we have argued for case (I4) already at the end of Section 4.8.
With these arguments we have finished our account on the extension of all three of the major singularity theorems of GR to regularity C 1 and we proceed to some conclusions and a general discussion.

Conclusions, alternatives & outlook
We begin this final section with a summary and some conclusions drawn from the results presented here.Then we will discuss some further perspectives of the approach at hand as well as alternative approaches to singularity theorems beyond the smooth setting.
In this review we have discussed the extension to C 1 -spacetimes of the classical singularity theorems of GR, which under physically reasonable condition assert causal geodesic incompleteness.These results, as well as the C 1 -extension of the Gannon-Lee theorem in [SS21], which we have avoided to discuss here, are entirely in the spirit of the classical theorems.That is, they extend the original conditions in a natural way, use essentially the same line of arguments and come to the same conclusions as the classical results, however, with one noteworthy extension: The Hawking-Penrose and the Gannon-Lee theorem add a further alternative to causal geodesic incompleteness, namely the branching of maximising causal geodesics.
The proofs are, at the one hand, based on the recent extensions of causality theory to low regularity Lorentzian metrics and, on the other hand, rely on an extension of the focusing results for causal geodesics.Indeed the latter, which is achieved via a regularisation approach, is the main technical advance presented here.
Next we briefly discuss the further prospects of this approach.It is generally expected that the causality parts of the results extend to C 0,1 -metrics,33 since such metrics still belong to the so-called causally plain ones, cf.[CG12, Def.1.16].The latter allow for essentially the same causality theory as smooth spacetimes, see [CG12,Thm. 1.25].Indeed, it is only below Lipschitz regularity that such core features of causality theory as the push-up principle and the openness of I + become an issue [GKS19], and the lightcones may form subsets of full measure [CG12, Ex. 1.11].
On the analytic side, it seems also feasible to extend the recent techniques to locally Lipschitz metrics.Here, one primary task is to extend the Friedrichs Lemma 4.4, which lies at the analytical core of the regularisation techniques, as it allows to derive from the distributional energy conditions of the singular metric useful surrogate energy conditions for the smooth approximations.Next, Lemma 4.2 which gives the precise speed of convergence of the regularised metrics to the singular one would have to be revisited.Also, the formulation of the energy conditions becomes more subtle since the curvature ceases to be an order-one distribution which forecloses the insertion of C 1 -vector fields which was essential at least for the genericity condition.Finally, one faces the problem that the right hand side of the geodesic equation now is merely locally bounded and it seems unavoidable to resort to non-classical solution concepts such as Filippov solutions [Fil88], see also the discussion in Section 4.1.
Still this is actually a long way from the largest possible class where the curvature can be (stably) defined in an analytical (distributional) way, i.e. the GT-class H 2 loc ∩ L ∞ loc , see Section 4.2.However, as remarked there, the quest is to at least go into the direction of regularity classes more closely linked to the PDE-approach to GR, e.g.g ∈ H 5/2+ε or ∇ ∈ L 2 .
A somewhat related topic is singularities in semi-classical and quantum theories of gravity.There the energy conditions are expected to be violated and to hold only in some averaged sense, see e.g.[Vis95].The question then arises whether under such assumptions the focusing effect persists, or whether the singularity theorems vanish altogether in the quantum regime, see [SG15, Sec.8.2], [Sen98, Sec.6.2].Indeed, focusing can be maintained under energy conditions averaged along causal geodesics, see e.g.[FG11] and the references therein.Similar results using index form techniques have appeared in [FK20] and recent work is concerned with energy conditions directly related to quantum energy inequalities [FFK21] and worldvolume energy inequalities [GKOS22].There is certainly a technical proximity of the methods used there and the ones described in this review, and future research will investigate their interrelations more closely.Now shifting away from the more analytic parts of the singularity theorems we first turn to causality theory.It has long been clear that causality theory is very robust in general.Most arguments are rather topological in nature and can actually be seen as belonging to frameworks more abstract than Lorentzian differential geometry.Indeed, Ettore Minguzzi in [Min19a] has recently put forward a very general theory of causal cone structures, that is a theory of upper semi-continuous distributions of cones over manifolds (which generalise the lightcone of Lorentzian metrics).In this setting it is indeed possible to establish the causal core of some of the singularity theorems: One may see the analytic concepts like the energy conditions, focusing results, etc. used throughout this note merely as tools that produce subsets in M that possess specific causality properties.Completely removing them from the arguments one arrives at purely causal results.To give some flavour of these, we quote a version of the Penrose theorem which appeared as [Min19a, Thm.2.67], for more details see Sec. 2.15 there.Theorem 5.1 (Causal Penrose Theorem).Let (M, C) be a globally hyperbolic closed cone structure admitting a non-compact stable Cauchy hypersurface.Then there are no compact future trapped sets and if S is non-empty and compact there is an inextendible future null geodesic entirely contained in E + (S).
Another and quite different approach to singularity theorems has been opened up in the context of the recently developed synthetic approach to Lorentzian ge-ometry put forward by Michael Kunzinger and Clemens Sämann in [KS18].The Lorentzian length spaces introduced there are the analogue of metric length spaces, which have long been used as an essential tool to extract the metric core of many notions and results in Riemannian geometry, see e.g.[BBI01].Lorentzian pre-length spaces (X, d, ≤, ≪, τ ) are metric spaces (X, d) together with with a preorder ≤ and a transitive relation ≪ contained in ≤ (which model the causal and timelike relations of Lorentzian geometry) and a lower semi-continuous map τ : X × X → [0, ∞] that satisfies the reverse triangle inequality (and models the Lorentzian distance function).Such a space is called a Lorentzian length space, if, in addition to some technical conditions, τ is intrinsic in the sense that the distance between points defined via the sup of the τ -length of connecting causal curves coincides with their τ -distance.Causality theory in Lorentzian length spaces [KS18, GKS19, AHCPS20, BORS22] extends standard causality theory beyond the spacetime setting to which it reduces for continuous spacetimes with strongly causal and causally plain metric.
In this setting it becomes possible to formulate synthetic versions of the singularity theorems, in the sense that the energy conditions are implemented as synthetic curvature bounds.A first theorem in Lorentzian length spaces that are warped products and which uses suitable sectional curvature bounds (implying Ricci curvature bounds in such geometries), based on triangle comparison is the following version of the Hawking theorem, cf.[AGKC22, Cor 6.2(ii)]: Theorem 5.2 (Synthetic Hawking Theorem).Let X be a geodesic length space, 34 and let Y = I × f X be a warped procduct with I = (a, b) and f : I → (0, ∞) smooth. 35.Assume that Y has timelike sectional curvature bounded below by 0 and that f is nonconstant.Then a > −∞ or b < ∞ and hence Y is past or future timelike geodesically incomplete.
Of course, implementing the classical energy conditions (SEC) or (NEC) rather amounts to Ricci curvature bounds than to sectional ones.Indeed, synthetic Ricci curvature bounds have been intensively studied in Riemannian geometry using optimal transport, see e.g.[Vil09].These techniques have recently been transferred to the smooth Lorentzian setting setting in [McC20,MS22] and further extended to the synthetic setting of Lorentzian length spaces by Fabio Cavaletti and Andrea Mondino in [CM20].The basic idea is that timelike lower Ricci bounds can be characterised in terms of the convexity of an entropy functional along l p -geodesics in the space of probability measures, where l p is the Lorentz-Wasserstein distance.The corresponding timelike curvature-dimension conditions TCD(K,N) and its weaker variant the timelike measure contraction property TMCP(K,N) then allow to formulate a version of the Hawking singularity theorem which we here quote in a loose way omitting technicalities, cf.[CM20, Thm.5.2] for the precise version: 34 That is a metric length space where each pair of points can be joined by a minimising curve. 35In the space Y , τ is defined via the sup of the length of future directed causal curves γ = (α, β), given by L(γ) = α2 − (f • α) 2 v 2 β with v β the metric derivative.Observe that this precisely models the smooth situation.
Theorem 5.3 (TMCP-Hawking Theorem).Let X be a timelike non-branching, globally hyperbolic Lorentzian pre-length space satisfying a TMCP-property.Let V be a Borel achronal future timelike complete subset with mean curvature bounded above.Then every future timelike geodesic starting in V has a bounded maximal domain of existence.
To sum up, in this review we have discussed the classical singularity theorems of GR and sketched the main arguments leading to their proofs.One may actually identify two main lines in these arguments, the analytical and the causal one.The former is concerned with providing focusing results for causal geodesics using the energy conditions which lead to the occurrence of conjugate or focal points and hence provide estimates on when causal geodesics stop maximising the Lorentzian length.The second, causal line of arguments gives criteria for maximising causal geodesic to exist.Confronting these two threads leads to a contradiction unless some of the causal geodesics become incomplete.
In the main part of this work we have extended the causal and the analytic line of arguments to Lorentzian metrics of regularity C 1 and have presented corresponding extensions of the classical theorems.The main achievements presented are on the analytic side using a regularisation approach that allows to deal with the distributional curvature associated with a C 1 -metric.
In the final discussion we have complemented these results with an overview of recent versions of singularity theorems in a purely causal setting as well as in a synthetic setting.It is clear that these results and techniques are still fresh, and that many interesting lines of research emerge from here.Also, it is unclear to date how the synthetic results precisely relate to the analytical approach put forward here and its possible extension to even lower regularity as discussed above, but see [KOV22] for a first work addressing this issue.
In any case, it can firmly be stated that the singularity theorems are not only an integral part of GR and Lorentzian geometry but, even more than half a century after they first emerged, they are still an interesting field of research holding many quests to be resolved in the future-both from the physical side, see e.g.[SG15, Sec.8], as well as from a mathematical perspective as laid out here.
Now global hyperbolicity is a strong condition and in many situations one wants certainly to do without it.Therefore the following generalisation comes in handy: We define the future Cauchy development D + (A) of any achronal set A as D + (A) := {p ∈ M : Every past inextendible causal curve through p meets A}.We define the past Cauchy development D − (A) analogously, set D(A) := D + (A) ∪ D − (A) and call it the Cauchy development of A. It is then a result that the interior D(A) • of the Cauchy development is globally hyperbolic, generalising the situation of a Cauchy surface, whose Cauchy development is the entire, globally hyperbolic spacetime.More specifically, for any acausal (i.e., no two points are causally related) topological hypersurface S, its Cauchy development is open and globally hyperbolic.Now, compactness of the causal diamonds allows one to show: Proposition 3.7 (Existence of maximiser).Let S be a closed, achronal, spacelike hypersurface.Then to any point p in the future Cauchy development D + (S) there runs a future directed maximising geodesic γ from S.Moreover, γ starts normal to S, has no focal point before p, and it is timelike unless p ∈ S.

(
I1) a compact achronal set without edge, (I2) a closed future trapped C 0 -surface P , (I3) a closed future trapped C 0 -submanifold P of co-dimension 2 < m < n such that (4.14) holds for any support submanifold and any future pointing normal null geodesic staring from it, (I4) a future trapped point in the sense of Definition 4.19.
This condition leads to a surrogate version for smooth approximating metrics which, just as (SNEC), secures focusing of approximating geodesics, cf.[KOSS22, Props.5.5, 5.7].From there we apply the machinery of Section 4.8 to "lift" the result to MNNB C 1 -spacetimes.Finally, we obtain [KOSS22, Prop.5.10, 5.11]: Proposition 4.18 (Trapped set from trapped C 0 -submanifold).Let (M, g) be a strongly causal, MNNB C 1 -spacetime and let P ⊆ M be a trapped C 0 -submanifold of codimension 1 < m < n.In case (i) m = 2 suppose (DNEC), and in case ii) 2 < m < n suppose (4.14) for any support submanifold P of P and any future pointing normal null geodesic staring from P Then E + (P )∩P is achronal, and E + (E + (P )∩P ) is compact or M is null geodesically incomplete.