Fermat Metrics

: In this paper we present a survey of Fermat metrics and their applications to stationary spacetimes. A Fermat principle for light rays is stated in this class of spacetimes and we present a variational theory for the light rays and a description of the multiple image effect. Some results on variational methods, as Ljusternik-Schnirelmann and Morse Theory are recalled, to give a description of the variational methods used. Other applications of the Fermat metrics concern the global hyperbolicity and the geodesic connectedeness and a characterization of the Sagnac effect in a stationary spacetime. Finally some possible applications to other class of spacetimes are considered.


Introduction
In this paper we present a review of the Fermat metrics associated with standard stationary spacetime and some applications for the studying of the geometrical and physical properties of such spacetimes. The Fermat metrics are two Finsler metrics of Randers type and are so called since they allow to extend the Fermat principle of classical optics. In an isotropic medium, light rays can be characterised as critical points of the optical length functional, or equivalently as geodesics of the conformal metric obtained dividing the ambient metric by the refraction index, see for instance [1,2], where the geometric aspects of the classical Fermat principle are emphasised.
The problem to extend the Fermat principle to relativistic optics was settled soon after the formulation of General Relativity by Albert Einstein in 1915, and the first formulations of the Fermat principle were obtained by Hermann Weyl for static spacetimes, see [3] and also [4], and by Tullio Levi-Civita for stationary spacetimes, see [5]. The basic role played by light rays, modelled by light-like geodesics in the four dimensional spacetime, in the physical properties of the spacetime was confirmed by observation the deflection of light shown by Eddington during the 1919 solar eclipse. It was the first experimental validation of general relativity.
Around the same time, a new branch of the calculus of variations appeared because of the work of Morse [6] and Ljusternik-Schnirelmann [7], starting with previous results from Poincaré and Birkhoff, the so called calculus of variations in the large. The aim of the calculus of variations in the large was to relate the set of critical points of a smooth function on a finite dimensional manifold to the topological properties of the manifold itself, see [8]. The objective of these researches was to relate the number and the type of geodesics joining two points on a Riemannian manifold, or the existence of a closed geodesic on a compact Riemannian manifold, to the topology of the underlying manifold. The variational methods introduced by Morse, Ljusternik, and Schnirelmann fit well with the geodesics, because they are the critical points of the action integral on the space of curves with fixed boundary conditions on the curves. Since variational problems involving functions and curves are infinite dimensional in nature, Morse developed a shortening procedure to reduce the infinite dimensional variational problem to a finite dimensional one.
In the 1960s, with the fundamental paper by Richard Palais and Steve Smale [9,10], the bases of an infinite dimensional critical point theory were settled, where the relations between the critical points of a smooth function defined on a Hilbert (or more generally a Banach manifold), and the topological properties of the manifold were exploited. Palais applied its theory to develop an infinite dimensional theory for the geodesics on a Riemannian manifold. The critical point theory has been applied also to obtain deep results in the study of non-linear ordinary and partial differential equations, with applications to physics, geometry, and engineering, see [11].
On the other hand, in the 1970s, the first examples of a gravitational lens were detected. The gravitational lensing is a relativistic optical phenomena, due to the deflection of light, for which an observer in spacetime can detect multiple images of the same light source, as a star, at the same time instant. From the point of view of the geometry of spacetime, an observer at a point p of the spacetime is joined by multiple light-like geodesics starting by the world line γ of the source. The gravitational lens effect is still an active research field in general relativity, for instance in the study of detection of dark matter, see for instance [12] for a wide geometric and physical description of the gravitational lensing effect and to [13] for a more recent review with many historic recalls.
The multiple image effect was predicted already in the 1930s by the Swiss astronomer Zuicky, and it can be caused both by physical properties of spacetime (as the presence of a massive object), encoded in the geometric properties of the spacetime metric, but also by the topological properties of the spacetime, for instance if the spacetime has a compact Cauchy surface. If a Fermat principle holds, we can characterise light-like geodesics as critical points of a suitable functional, and apply the methods of the calculus of variations in the large to relate the multiple image effect to the topological and geometric properties of spacetime.
In this paper, we focus our attention on standard stationary spacetimes (or, in general, Lorentzian manifolds without any restriction on the dimension of the manifold), which are defined at Section 4. In this class of spacetimes, a Fermat principle can be stated in terms of the Fermat metrics associated to such manifolds. They are two Finsler metrics of Randers type F + and F − , defined not on the whole spacetime, but only on its spatial component. We recall that Randers metrics wee introduced by Randers in the paper [14] in order to define a Lagrangian for a charged particle in relativistic spacetime. We will state the following version of the Fermat principle, see Section 4. It states that the light-like geodesics on standard stationary spacetimes can be put into bijective correspondence with the (pre)geodesics of one of the Fermat metrics, whose geodesics are the critical points of the action integral for the Fermat metrics.
A light-like curve is a future (respectively past) pointing light-like geodesic if and only its spatial projection is a pregeodesic for the Fermat metric F + (respectively, F − ).
In the literature, it is usually written lightlike A first version of the Fermat principle in stationary spacetimes was proven by Tullio Levi-Civita, see [5]. Other versions were proved in [15,16], while in [17] the Fermat metrics are explicitly introduced for the first time. We point out that, since the notions of light-like curve and light-like geodesics are conformally invariant, the Fermat principle above holds for conformally standard stationary spacetimes.
In this paper we present some applications of the Fermat metrics to the physical and geometric properties of a standard stationary spacetime. The paper is organised in the following way. In Section 2 we recall the main properties of Finsler metrics.
In Section 3 we recall the main properties of Lorentzian manifolds, because spacetimes of general relativity are modelled on four dimensional Lorentzian manifolds.
In Section 4, standard stationary Lorentzian manifold and the Fermat metrics associated are introduced. Moreover the Fermat principle for light rays involving the Fermat metrics is stated.
In Section 5 we present some physically relevant examples of standard stationary spacetimes and we evaluate their Fermat metrics.
In Section 6 we recall the main results on critical point theory, in particular both the Ljusternik-Schnirelmann and the Morse theories.
In Section 7 the main results on critical point theory for the geodesics on Finsler manifolds, and also Riemannian manifolds, are presented, using the infinite dimensional approach by Palais. In Section 8 an alternative variational principle for geodesics on Randers metrics are presented. Sometimes this principle could be useful, because it involves a functional of class C 2 , while in Finsler geometry the action functional is only of class C 1,1 .
In Section 9 we present the variational theory for light rays joining a point (the observer) with a vertical line (the world line of the light source) in a standard stationary spacetime. By the Fermat principle, this variational theory is equivalent to the variational theory for the geodesics of one of the Fermat metrics, according the light rays are future or past pointing. We produce conditions which guarantee the multiple image effect and a proof of the odd image theorem in the gravitational lens effect. The case of manifolds with boundary is treated because many physically relevant spacetimes of general relativity have boundaries, similarly to Schwarzschid, Reissner-Nordstrom, and Kerr spacetimes.
In Section 10 we relate the global hyperbolicity of a stationary spacetime to the completeness of the Fermat metrics. This is a very interesting topic because global hyperbolicity is a crucial assumption in the study of evolution equations on a spacetime, as Einstein equations and wave maps. The global hyperbolicity depends on the global behaviour of the Fermat metrics, and this gives further merit to these Finsler metrics.
In Section 11, we study a more geometric problem, the geodesic connectedness of a standard stationary spacetime, that is the existence (and, eventually, also the multiplicity) of geodesics joining two arbitrary points. We shall present the main results on this problem and we shall also see that the Fermat metrics now play a role.
Finally, in Section 12 we give a geometric description of the Sagnac effect in terms of the Fermat metrics. The Sagnac effect is an interference effect caused by an asymmetric behaviour of light rays moving on a rotating platform, dependent on if the light rays move clockwise or counterclockwise. Geometrically it is a manifestation of an asymmetry, and this asymmetry is measured by the Fermat metrics. When the platform does not rotate, the light rays have the same behaviour, the Sagnac effect is null and there is no asymmetry because in this case the metric is static and the two Fermat metrics are equal to a Riemnanian conformal metric. So in this case the Sagnac effect appears when the optical geometry is of Finsler type, while when the optical geometry is Riemannian, the Sagnac effect does not appear. Therefore, the description of the Sagnac effect by Finsler geometry is not a mere generalisation of a Riemannian one, because it appears only in the presence of a Finslerian optical geometry. In Section 13 we present some new research direction on Fermat metrics.

Finsler Metrics
In this section we introduce the main notions and properties of Finsler geometry. The roots of Finsler geometry are already as a generalised metric in the Habilitationsvortrag of Bernhard Riemann, where a local arc length was defined as a homogeneous function of degree 1 in the fibre coordinates [18]. The main properties of Finsler geometry was developed by Finsler in his 1918 thesis [19]. Finsler geometry had a parallel development, but it was essentially considered as a generalisation of Riemannian geometry. Perhaps one of the main arguments against Finsler was the difficulty in representing tensor fields or differential operators in a free-coordinate form. The things began to change when Katok discovered a family of Finsler metrics of Randers type on the two sphere S 2 admitting exactly two closed geodesics [20], while it was conjectured that any Riemannian metric on S 2 should admit infinitely many geometrically distinct closed geodesics, as successively proved by Bangert and Franks. Therefore, Finsler geometry has its own peculiarity with respect to Riemannian geometry. Moreover, many applications of Finsler geometry to many fields of applied sciences have appeared, creating further interest in the study of Finsler geometry, see for instance [21][22][23][24][25]. A big step in the development of Finsler geometry was obtained by Chern and their school in the search for global operators, as the Chern connection [26,27]. One of the motivations of this paper is to contribute to present geometrical and physical models, in which Finsler geometry has a valuable role.
In this section, any smooth manifold has to be considered as class C ∞ , even if it suffices to be of class C 2 to introduce the main notions. We refer to the book of Bao-Chern-Shen [28] for the details.
Let M 0 be a smooth manifold and let TM 0 be the tangent bundle of M 0 . Moreover, we shall denote by 0 the zero-section of TM 0 and by TM 0 \ 0, the so called slit tangent bundle. Definition 1. A Finsler metric on a smooth manifold M 0 is a function F : TM 0 → R + , such that for any (x, v) ∈ TM, v = 0, the quadratic form is positive definite on TM 0 \ 0. It is sometimes called the fundamental tensor associated to the Finsler metric F. The couple (M 0 , F) is called Finsler manifold.
Then a generic Finsler metric is a positively homogeneous function of degree 1 in the tangent variable, while a reversible Finsler metric is absolutely homogeneous.
On a Finsler manifold, metric properties can be introduced as for Riemannian manifolds. Let (M 0 , F) be a Finsler manifold, let γ : [0, 1] → M 0 NO, never also in the rest of the paper, where enlighted be a (piecewise) smooth curve, the length L(γ) of γ with respect to the Finsler metric F is defined setting We remark that the length L(γ) of the curve γ depends on the orientation of γ, unless F is reversible.
We can define the distance d(p, q) between two points p, M 0 setting where C(p, q) is the set of the piecewise smooth curves γ, such that γ(0) = p and γ(1) = q.
The function d F satisfies all the properties of a distance, except for the symmetry d F (p, q) = d F (q, p) for any couple of points p and q. The distance d F is symmetric if the metric F is reversible. Moreover, the metric d F induces a topology on M 0 which coincides with the initial topology of the manifold.
The asymmetry of a Finsler metric F permits to introduce two different notions of completeness on a Finsler manifold, the forward and backward completeness, corresponding to forward and backward Cauchy sequences. Obviously, the two notions of completeness coincide if the metric F is reversible. Let (M 0 , F) be a Finsler manifold and let p ∈ M 0 , for any r > 0, the forward metric ball B + (p, r) and the inward metric ball B − (p, r) are defined setting Obviously, the balls B + (p, r) and B − (p, r) are equal if the metric is reversible. A sequence (x k ) k∈N is a forward Cauchy sequence (respectively, backward Cauchy sequence) if for any > 0, there exists a number ν ∈ N, such that for .The Finsler manifold (M 0 , F) is said to be forward complete (respectively, backward complete) if any forward (respectively, backward) Cauchy sequence is convergent. The two notions are not equivalent and there examples of Finsler manifolds which are complete only on a direction, see [28]. Of course, the two notions of completeness agree if the Finsler metric is reversible.
In this paper, we shall consider two fundamental classes of Finsler metrics, the square root of a Riemannian metric and Randers metrics. For other many interesting examples of Finsler metrics we refer to [28].
Example 1 (Riemannian metrics). Any Riemannian metric h on a manifold M 0 induces a reversible Finsler metric F setting F = √ h.
In this paper, the most important example of Finsler metrics is given by Randers metrics Example 2 (Randers metrics). A Randers metric on a manifold M 0 is a Finsler metric F of the form where h is a Riemannian metric and W(x) is a smooth vector field on M 0 , such that for any Assumption (8) guarantees that F is a Finsler metric, moreover it is not reversible, unless the vector field W(x) is null. Randers metrics were introduced as Gunnar Randers in the paper [14] and represent the Lagrangian of a relativistic charged particle under the action of an electromagnetic field, when h is the Minkowski metric. Another equivalent way to define a Randers metric is using differential forms. Indeed, let ω the differential form metrically equivalent to the vector field W, that is As in Riemannian geometry, we can define geodesic curves and they play a basic role in the study of the geometric properties of a Finsler manifold. Here we introduce geodesics using a "variational approach", for a more geometric approach, see [28]. For a smooth (or at least piecewise smooth) curve x(s) : [0, 1] → M 0 , the length functional L(x) is defined at (6). Moreover, the action functional (sometimes called energy functional) E(x) is defined setting The geodesics for the Finsler metric F are the curves which are stationary points of the action integral E(x), with suitable boundary conditions, for instance curves joining two fixed points of M 0 , closed curves, or curves joining two fixed submanifolds of M 0 , see [17,28].
The geodesic equations for the the Finsler metric F represent the Euler-Lagrange equations for the action functional (10) and can be written as where Dẋẋ is the Chern connection for the Finsler metric F with referenceẋ, see [28] for the definition of the Chern connection. In local coordinates, a geodesic satisfies the system of non-linear differential equationsẍ Here the functions Γ k ij (x,ẋ) are the Christoffel symbols of the Chern connection. We point out that the Christoffel symbols for the Chern connection are not defined on the manifold M 0 , as for the Christoffel symbols of the Levi-Civita associated to a Riemannian metric, but rather on the slit tangent bundle TM 0 \ 0.
As we shall see in Section 7, the action functional E(x) is only of class C 1,1 on the Sobolev manifold of H 1 curves on the manifold M 0 , and not of class C 2 as for the action integral of a Riemannian metric.
We point out that, as in Riemannian geometry, the property of a curve to be a geodesic depends on the parameterisation of the curve itself. A definition independent on parameterisation of the curve is that of pregeodesic curve. A curve is a pregeodesic if it admits a reparemetrisation which is a geodesic.
Moreover, we also point out that if the Finsler metric is not reversible, the notion of geodesic depends on the orientation of the curve. In other words if a curve is a geodesic, the reverse oriented curve could be not a geodesic, too. Clearly this property depends on the asymmetry of the Finsler metric.

Lorentzian Manifolds
In this section we introduce some basic notions in Lorentzian geometry, we refer to [29][30][31][32] for the details about Lorentzian geometry and its applications to general relativity.

DONE
A Lorentzian manifold is a couple (M, g), where M is a connected smooth manifold and g is a Lorentzian metric on M, that is g is a smooth symmetric tensor field such that for any z ∈ M, the bilinear form g(z) : T z M × T z M → R is non-degenerate and has index equal to 1. In general relativity, a spacetime is described by a four dimensional Lorentzian manifold. The geometric properties of the Lorentzian manifold are related to the physical properties of the relativistic spacetime by the Einstein equations.
where Ric g and S g are, respectively, the Ricci tensor and the scalar curvature of the metric g, while the double symmetric tensor T is called the energy-impulse tensor and it encodes all the physical properties of the spacetime. The bilinear form g(z) is not positive definite so this property allows to introduce the so called causal character of vectors. Let z ∈ M, a vector v ∈ T z M is said to be: This tripartition of the tangent vectors is called causal character of tangent vectors. On a Lorentzian manifold a time orientation, introducing the past and the future, can be defined using a continuous time-like vector field Y(p) on M, that is Y(p) ∈ T p M is a time-like vector tangent to M at p. It is well known that any non-compact Lorentzian manifold M admits a continuous time-like vector field, while some topological restriction is needed, namely the Euler-Poincaré characteristic must be 0, to have the existence of a continuous time-like vector field on a compact Lorentzian manifold.
We point out that the scalar product g(p)[v, w] of two causal vectors, one of which is time-like, is never null, see for instance [30].
The causal character of vectors can be extended also to define the causal character of a curve on a Lorentzian manifold. Moreover, the time orientation of vectors can be extended to causal curves.  As in Riemannian or Finslerian geometry, an important class of curves in a Lorentzian manifolds is given by geodesics. where D sż (s) is the covariant derivative of the vector fieldż along z induced by the Levi-Civita connection of the Lorentzian metric g.
In local coordinates of the manifold M, a curve z is a geodesic if, and only if, it satisfies the system of ordinary differential equations . . , dim M}, and the functions Γ k ij are the Christoffel symbols for the metric g. We remark that the property of a curve to be a geodesic depends on its parameterisation. A notion independent on the parameterisation is the notion of pregeodesic curve. It is not difficult to show that the notion of pregeodesic in independent on its parameterisation. Moreover, it is not difficult to show that a non-constant curve is a pregeodesic if, and only if, all it admits is a reparemeterisation, which is a geodesic.
The causal character of a geodesic can be defined. Indeed, by the main properties of the Levi-Civita connection associated to a Lorentzian metric, it can be proved that if z : [0, 1] → M is a geodesic, then there exists a constant E z , such that for any s ∈ [0, 1],

Definition 7.
A non-constant geodesic z : [0, 1] → M on the Lorentzian manifold (M, g) is said to be: Notice that we have not defined the causal character of non-constant geodesics, which are usually considered as space-like ones. The causal character of geodesics plays a fundamental role in general relativity. Indeed time-like geodesics represent the world lines of free falling particles for which only gravity, which is encoded in the metric g, acts.
Light-like geodesics represent the world lines of light rays, while space-like geodesics have no physical meaning because traveling on a space-like geodesic, a particle should be faster than light.
One of the most important properties of geodesics is that they satisfy a variational principle. Let (M, g) be a semi-Riemannian manifold, we consider the action integral defined on the set of C 1 (or piecewise C 1 ) curves z : [0, 1] → M. The action integral is smooth and its first variation is given by where ζ : [0, 1] → TM is a smooth vector field long z and D s ζ is the covariant derivative of ζ along z. Now, let p and q two points of M and let Ω(p, q; M) be the set of the smooth curves z : [0, 1] → M, such that z(0) = p and z(1) = q. Moreover, let T z Ω(p, q; M) be the set of the smooth vector fields ζ along z, such that ζ(0) = 0 and ζ(1) = 0. Then the following variational principle gives a characterisation of the geodesics joining p and q, see for instance [33].

Theorem 1.
A curve z(s) ∈ Ω(p, q; M) is a geodesic for the metric g joining the points p and q, if, and only if, z(s) is a critical point of the action integral E on Ω(p, q; M), that is for any ζ ∈ T z Ω(p, q; M), In the language of the calculus of variations, the geodesics equation D sż = 0 is the Euler-Lagrange equation of the action integral. Sometimes the vector fields ζ ∈ T z Ω(p, q; M), that are with null boundary conditions, are called variational vector fields.
A similar variational principle holds for geodesics satisfying more general boundary conditions, as for geodesics with endpoints on two submanifolds of (M). In particular, a variational principle holds for closed geodesics. Let Λ(M) be the set of the smooth (or at least piecewise smooth) closed curves z : [0, 1] → M, so that z(0) = z(1). Moreover let T z Λ(M) be the set of the vector fields ζ along z, such that ζ(0) = ζ (1). Then the following variational principle holds:

Theorem 2. A closed curve z(s) ∈ Λ(M) is a closed geodesic for the metric g, if, and only if, z(s) is a critical point of the action integral E on
In particular the constant curves are closed geodesics for the metric g, so the search of closed geodesics needs to avoid the trivial constant curves.
We now analyze the second variation of the action integral. Let z(s) ∈ Ω(p, q; M) be a geodesic for the metric g joining the points p and q, then the second variation E (z) of the action integral E at the geodesic z is the bilinear form is well defined and for any ζ, ζ ∈ T z Ω(p, q; M), we have where R(·, ·)· denotes the Riemann tensor for the metric g.
We define now the notion of Jacobi vector field along a geodesic and of non-conjugate points for the metric g.

Definition 8.
A vector field ζ ∈ T z Ω(p, q; M along a geodesic z joining two points p and q is called Jacobi field I use the bold to distinguish from the italic if, and only if, for any ζ ∈ T z Ω(p, q; M), By an integration by parts and the fundamental lemma of calculus of variations, we have that ζ ∈ T z Ω(p, q; M) is a Jacobi field along z if, and only if, it satisfies the Jacobi equation In other words, the Jacobi fields along z are the variational vector fields along which satisfies the Jacobi equations. Equivalently the set of Jacobi fields along z is the kernel of E (z) on T z Ω(p, q; M), so it has the structure of a vectorial subspace of T z Ω(p, q; M). Since the Jacobi equations consist of a linear system of differential equations of second order, we find that the dimension of such vector space is finite, and it can be proved that it is at most n = dim M, see [30].
We give now the important definition of conjugate points along a geodesic, which holds both for Riemannian and Lorentzian manifolds, and more in general for any semi-Riemannian manifold, see [30]. Clearly the multiplicity of a conjugate point is finite and it is equal at most to n = dimM. We give now the definition of geometric index of a geodesic.

Definition 10.
As in the definition above, the geometric index µ(z) of the geodesic z is the number of conjugate points z(s), s ∈]0, 1[, counted with their multiplicity.
One of the most important results on geodesics on Riemannian manifold is the Morse index theorem, which states that the geometric index µ(z) of a Riemannian manifold is finite and it is equal to the Morse index of the geodesic as a critical point of the action functional, see for instance [34]. The Morse index theorem still holds for time-like and light-like geodesics on Lorentzian manifolds, see [31], while it does not hold for space-like, for which infinitely many conjugate points may exist.
We give now the definition of non-degenerate geodesic and non-conjugate points for a semi-Riemannian metric g.

Definition 11.
A geodesic z ∈ Ω(p, q; M) joining two points p and q in the semi-Riemannian manifold (M, g) is said to be non-degenerate if the point q is non-conjugate to p along z, that is the null field ζ = 0 is the unique Jacobi field along z with ζ(0) = 0 and ζ(1) = 0, or equivalently KerE (z) = {0}. Definition 12. Two points p and q of a semi-Riemannian manifold (M, g) are said to be nonconjugate if any geodesic z ∈ Ω(p, q; M) joining p and q is non-degenerate, or equivalently p and q are non-conjugate along any geodesic joining them.
From the variational point of view, the fact that p and q are non-degenerate is equivalent to say that the critical points of the action integral E(z) on Ω(p, q; M) are nondegenerate, and this is equivalent to say that E(z) is a Morse function.
It can be proved, using the Sard theorem, that any couple of points p and q, except for a residual set of M × M, are non-degenerate, see [34].

Standard Stationary Lorentzian Manifolds and the Fermat Metrics
We introduce the class of standard stationary spacetimes and the Fermat metrics in these spacetimes. We fix a Riemannian manifold (M 0 , ·, · ) and we consider M 0 on a smooth vector field δ(x) and a smooth strictly positive function β(x) : M 0 → R. Then the standard stationary Lorentian manifold defined by the data (M 0 , ·, · , δ(x), β(x)) is the Lorentzian manifold (M, g), where M is the product manifold M = M 0 × R is the product manifold, and the metric g is defined in the following way: for any = (x, t) ∈ M = M 0 × R, and for any ζ = (ξ, τ), It can be proven that the bilinear form g(z) is non-degenerate and its index is equal to 1, so it defines a Lorentzian manifold, which we shall call the standard stationary Lorentzian manifold with initial data (M 0 , ·, ·, , δ(x), β(x)). The vector field is sometimes called the shift, while the positive function β(x) is called the lapse of the standard stationary metric g, see, for instance, [33] for other details.

Definition 13.
A standard stationary Lorentzian manifold is said to be static if the shift vector δ(x) is null, so the metric tensor of a standard static Lorentzian manifold has the form: Many physically relevant spacetimes of general relativity are standard static or stationary, as the Minkowski spacetime of special relativity, the Schwarzschild spacetime outside the event horizon and the Reissner-Nordstrom spacetime outside the external event horizon are static, while the Kerr spacetime outside the stationary surface is stationary. In Section 5, we shall define these spacetimes and we evaluate their Fermat metrics.

Remark 2.
The notion of standard stationary is a particular case of the more general notion of stationary Lorentzian manifold. A Lorentzian manifold (M, g) is said to be stationary if it admits a time-like and killing vector field Y(z). Therefore Y(z) satisfies g(z)[Y(z), Y(z)] < 0 and the flow generated by the vector field Y(z) is given by isometries at any stage. A stationary Lorentzian manifold is said to be static if the orthogonal distribution Y ⊥ on the tangent bundle TM is integrable, see for instance [30,31] for the details. In a standard stationary Lorentian manifold the time-like vector field Y(z) = Y(x, t) = (0, 1) is killing, essentially because the shift δ(x) and the lapse β(x) does not depend on the variable t (they are invariant by the flow of Y). Moreover, a standard static Lorentzian manifold is static. An example of a non-standard stationary Lorentzian manifold is the Godel spacetime, see [29].
We now define the Fermat metrics associated to a standard stationary Lorentzian manifolds. They are two Finsler metrics of Randers type, which reduce to a Riemannian metric when the metric is standard static. Such metrics are related to light-like vectors, light-like curves and light-like geodesics by a version of the Fermat principle, which is an extension of the classical Fermat principle to light rays in standard stationary spacetimes.
Let (M, g) be a stationary spacetime, with M = M 0 × R and the metric g is given by (18). We shall assume that all curves are defined in the integral [0, 1].
Let z(s) = (x(s), t(s)) : [0, 1] → M be a smooth curve and assume that z(s) is a light-like curve, then z(s) satisfies We solve Equation (20) with respect toṫ and we obtain two solutions: Consider the Riemannian conformal metric g 0 = ·, · /β(x), then we have alsȯ Then, in a fixed spatial curve x(s) : [0, 1] → M 0 , there are exactly two light-like curves z + (s) = (x(s), t + (s)) and z − (s) = (x(s), t − (s)) on M, having the curve x(s) as a spatial component. We notice that by (21) and (22), the light-like curve z + is future pointing, while the light-like z − is past pointing.
In other words, we have shown that any spatial curve x(s) : [0, 1] → M 0 can be lifted up to exactly two light-like curves z + and z − , being z + pointing to the future and z − pointing to the past. Therefore, we have found two Finsler metrics of Randers type F + , F − : TM 0 → R + , defined for any (x, v) ∈ TM 0 , x ∈ M 0 and v ∈ T x M 0 : and We call them the Fermat metrics associated to the standard stationary Lorentzian manifold (M, g) and were explicitly introduced in [17]. Now, let z + = (x(s), t + (s)) : [0, 1] → M be a future pointing light-like curve, then the travel time T + (z) of z + is a defined setting.
so we have T(z + ) = 1 0 F + (x,ẋ)ds, that is the travel time T(z + ) is equal to the length of the spatial projection x with respect to the Finsler metric F + . Analogously so that T(z − ) = − 1 0 F(x,ẋ)ds, that is the travel time T(z − ) is equal to minus the length of the spatial projection x with respect to the Finsler metric F − .
So we have characterised the travel time of a light-like curve z(s) = (x(s), t(s)) of a standard stationary Lorentzian metric as the length of its spatial part x with respect to a Finsler metric of Randers type, the metric This characterisation recalls the Fermat principle in classical geometric optics. In an isotropic medium, modelled with a smooth manifold M 0 with Riemannian metric g 0 , the light speed depends only on the position and not on the direction, and is proportional to the refraction index n(x) > 0 of the medium. The travel time of a curve x(s) : [0, 1] → M 0 in the medium is given by so the the travel time T(x) is equal to the length of the curve x with respect to the conformal metric g 0 /n, which is often called the optical metric.
Theorem 3 (Fermat principle in classical optics ). The light rays joining two fixed points in the isotropic medium M 0 with refraction index n(x) are the stationary points of the travel time T(x) and they are the geodesics of the optical metric g 0 /n.
Inspired by the Fermat principle in classical optics, we give the following definition Definition 14. Let (M, g) be a standard stationary spacetime, the Fermat metrics associated to (M, g) are the Finsler-Randers metrics defined at (23) and (24) and are often called optical metrics associated to (M, g).

Remark 3.
If the standard stationary spacetime (M, g) is also static, the vector field δ is null and the two Fermat metrics F + and F − are both equal to the reversible Finsler metric √ g 0 , induced by the Riemannian metric g 0 , that is Therefore, in a static spacetime we have only one optical metric and it is a Riemannian metric.
We state an extension of the Fermat principle of classical optics to standard stationary Lorentzian manifolds. Let (M, g) be a stationary Lorentzian manifold, with M = M 0 × R and the standard stationary metric g given by (18). Let p = (x 0 , 0) ∈ M and let γ(s) = (x 1 , s) : R → M be a vertical line which represents the world line of a stationary observer, with x 0 , x 1 ∈ M 0 . We give a characterisation of light-like geodesics (light rays) The following theorem holds.
Theorem 4 (Fermat principle for future pointing light-like geodesics). A light-like curve z + (s) = (x + (s), t + (s)) : [0.1] → M is a future pointing light-like geodesic for g joining p and the curve γ if, and only if, the temporal component t + (s) satisfies (21) with x = x + , and the spatial component x + is a pregeodesic for the Fermat metric F + . Moreover, the travel time T(z + ) of z + is equal to the length F + (x + ) of x + with respect to F + . Finally, the spatial part x + is a pregeodesic for the Fermat metric F + parameterised in order that g 0 (ẋ + ,ẋ + ) + g 0 (δ(x + ),ẋ + ) 2 is a constant.
An analogous result holds for past pointing light-like geodesics.  (22), with x = x − and the spatial part x − is a pregeodesic for the Fermat metric F − . Moreover, the travel time The previous theorems was first proved in [17], exploiting the role of the Fermat metrics in the statement. Previous result on variational characterisations of light rays for stationary spacetimes were obtained by H. Weyl in [3], who stated a Fermat principle in static spacetimes, see also [4]. A Fermat principle for light rays on stationary spacetimes was stated by Tullio Levi Civita, see [5] and the paper [15].

Remark 4.
We point out that all the previous results for the Fermat metrics and the Fermat principles hold for the class of conformally standard stationary Lorentzian manifoldḡ = α(x, t)g, where g is a standard stationary metric and α : M → R is a smooth and strictly positive function. This is due to the fact that the light-like curves and light-like pregeodesics are invariant by a conformal change of the metric. This fact permits to apply the results also to wider classes of Lorentzian metrics.

Remark 5.
A Fermat principle holds for arbitrary spacetimes, and it was first stated by Kovner [35], see also [36].

Some Example of Fermat Metrics
In this section, we present some examples of Fermat metrics, associated to physically relevant spacetimes of general relativity, for more details on the physical properties of these spacetimes see [29].

Example 3.
[Minkowski spacetime] The (n + 1)-Minkowski spacetime (R n+1 , g M ) of special relativity is a static Lorentzian manifold, with where ·, · is the standard, positive definite Euclidean scalar product in R n . Clearly, the Minkowski spacetime (R n+1 , g M ) is a static Lorentzian manifold, and the two Fermat metrics F + and F − are equal to the Euclidean metric ·, · . Example 4 (Swarzschild spacetime). Let m be a positive constant, let (r, θ, φ) the polar coor- The Scharzschild metric ds 2 on the manifold M is defined setting in the coordinates (r, θ, φ, t), where The couple (M, ds 2 ) is called the Schwarzschild spacetime outside the event horizon. From the physically point of view, the Schwarzschild spacetime represents the empty spacetime outside a spherically symmetric massive body of mass m, as a star. The number r = 2m is called the Schwarzschild radius. The Swarzschild spacetime is static, so the two Fermat metrics F + and F − are equal to the Riemannian conformal metric which is a complete Riemannian metric on M 0 .
Example 5 (Reissner-Nördstrom spacetime). Another physical relevant example of static spacetime is the Reissner-Nördstorm spacetime, which represents the spacetime outside a spherically symmetric charged body of mass m and electric charge e. The Reissner-Nördstorm metric is defined as in (27), but now with β(r) defined as Assuming that m > e, then the Reissner-Nördstrom spacetime is defined and is static on the Again, the two Fermat metrics are equal.
We present now a stationary and non-static spacetime, the Kerr spacetime, for which the Fermat are defined and distinct.
Example 6 (Kerr spacetime). The Kerr spacetime outside the stationary limit surface represents the stationary, axis-symmetric, and asymptotically flat spacetime outside a rotating massive object. Let m and a be two constant, with m > 0, which represent, respectively, the mass and the angular momentum of the massive body as measured from infinity. In the Boyer-Lindquist coordinates (r, θ, φ, t) coordinates, the Kerr metric is given by where λ(r, θ) = r 2 + a 2 cos 2 θ, and ∆(r) = r 2 − 2mr + a 2 .
Notice that, when a = 0, the Kerr metric reduces to the Schwarzschild metric. Assume that m 2 > a 2 , and consider the stationary limit surface then the metric (28) is stationary on the Kerr spacetime outside the stationary limit surface and, equipped with the metric (28), is called the Kerr spacetime outside the stationary limit surface. Now we evaluate the Fermat metrics in the spacetime outside the limit surface. Let z(s) = (r(s), θ(s), φ(s), t(s)) : [0, 1] → M a be a light-like curve with image in the Kerr spacetime outside the stationary limit surface, so it satisfies: Setting and solving with respect toṫ we obtaiṅ Therefore, we obtain the two Fermat metrics for any x = (r, θ, φ) and for any tangent vector and where A(r, θ) = λ(r, θ) β(r, θ)∆(r) ,

A Review of Critical Point Theory
In this section, we review the main results of critical point theory that we shall apply to Finsler manifolds and light rays. The critical point theory on infinite dimensional manifolds was developed by Richard Palais and Stephen Smale in a series of papers in the 1960s, see [9,10,37,38]. Then the theory has been applied to the study of the solutions of non-linear ordinary or partial differential equations, and geometric analysis. Since, from [9], the theory has been applied to the study of geodesics on Riemannian manifolds. Successively geodesics on Finsler manifolds and Lorentzian manifolds have been studied. We refer to [11,33,39,40] for the proof of the main results on critical point theory and several applications to non-linear differential equations.
Let (X, h) be a smooth (C 3 ), possibly infinite dimensional Riemannian manifold, we consider a functional f : We introduce now the Palais-Smale (PS) compactness condition, which plays a basic role in infinite dimensional variational problems.
Definition 15. Let f : X → R be a functional of class C 1 on a Riemannian manifold (X, h), we say that f satisfies the Palais-Smale (PS) compactness condition, if any sequence (x k ) k∈N of points of X, such that: admits a converging subsequence. Here · is the norm induced on the tangent bundle TX by the Riemannian metric h on X We state now that the theorem on the existence of a minimum point of a smooth functional satisfying the (PS) condition, see [39,40].
Theorem 6. Let f : (X, h) → R be a C 1 functional, defined on a complete Riemannian manifold (X, h), and assume that f is bounded from below and satisfies the Palais-Smale condition on X. Then the infimum is attained by f and there exists a point x 0 ∈ X, such that f (x 0 ) = inf X f .
We now state a multiplicity result of critical points of the functional f , for which the topological properties of the manifold X plays a basic role. We now introduce a topological invariant, the Ljusternik-Schnirelmann category of a topological space, which gives a measure of how much a topological space has a rich topology. The Ljusternik-Schnirelmann category gives a lower bound on the number of critical points of a C 1 functional. Definition 16. Let X be a topological space and let A be a subset of X, the Ljusternik-Schnirelmann category of the set A in X, denoted by cat X (A) is equal to the minimal number of closed and contractible subsets of X which cover A. If such a minimal number does not exist, we set cat X (A) = +∞. Moreover we set cat(X) = cat X (X).
Theorem 7. Let f : (X, h) → R be a C 1 functional defined on a complete Riemannian manifold (X, h), bounded from below and satisfying the Palais-Smale condition. Then the functional f has at least cat(X) critical points. Moreover, if cat(X) = +∞, then there exists a sequence x n of critical points of f , such that f (x n ) → +∞.
The proof of this theorem was obtained by Ljusternik and Schnirelmann at the end of the 1920s. If the manifold X is not contractible, for instance X is a Hilbert space, the previous theorem reduces to existence of a minimum point for f .
We now present the main results on Morse theory for a functional f defined on a complete Riemannian manifold, bounded from below and satisfying the Palais-Smale condition. Morse theory gives more precise estimates on the critical points of a functional, with respect to the Ljusternik-Schnirelmann theory, in particular on the number of critical points. However, in order to develop the Morse theory, we need a major regularity assumption on the functional f , namely it has to be of class C 2 . Moreover, we assume that the critical points of f have finite Morse index and are non-degenerate, and sometimes we shall say that f is a Morse function.
Therefore, assume that the functional f is of class C 2 on the Riemannian manifold (X, h) and let x be a critical point of f , we can define the Hessian H f (x) of the functional f at the point x in the following way: denote by T x X the tangent space at x to X, the Hessian where γ :] − ε, ε[→ X is a smooth curve, such that γ(0) = x,γ(0) = ξ, and then we extend H f (x) by polarisation to any couple of tangent vectors.
Clearly the Morse index and the augmented Morse index may be equal to +∞ if X is an infinite dimensional Riemannian manifold. The critical point x is said non-degenerate if the linear operator f (x) : T x X → T x X induced by H f (x) on T x X, equipped of the Hilbert space structure h(x), is an isomorphism. The functional f is said to be a Morse function if it is of class C 2 on the manifold (X, h) and all its critical points are non-degenerate.
We give now some recalls on algebraic topology that we need to state the main results of Morse theory, see for instance [41] for the details.
We consider a topological pair (A, B), where A is a topological space and B is a subspace of A, and we consider an algebraic field K a field. For any k ∈ N, H k (A, B; K) the k-th relative homology group, with coefficients in K, of the topological pair (A, B) are well defined. Moreover, since K is a field, the homology group H k (A, B; K) is also a vector field and its dimension β k (A, B; K) ∈ N ∪ {+∞} is called the k-th Betti number of (A, B) with respect to K.
The Poincarè polynomial P (A, B; K)(r) of the pair (A, B), with respect to the field K is a formal series whose coefficients are positive cardinal numbers belonging to N ∪ {+∞}, and it defined setting We state now the Morse relations for a Morse functional of class C 2 , and having critical points of finite index. The Morse relations link the set of the critical points of the functional to the Poincaré polynomial of the manifold with respect to the field K. Again, the Morse relations allow to relate the set of critical points of the functional f to the topological properties of the manifold X measured by the Poincaré polynomial. As corollaries of the Morse relations we obtain the classical Morse inequalities and the total Betti number formula for the critical points of f . For the proofs see [39,40].
Theorem 8. Let f : X → R be a functional of class C 2 , defined on a complete Riemannian manifold (X, h) of class C 3 . Assume that f is bounded from below, satisfies the Palais-Smale (PS) condition on X. Moreover, assume that all the critical points of f are non-degenerate and have finite Morse index m(x, f ), and denote by K( f ) the set of the critical points of f .
Then for any field K there exists a formal series Q(r) with coefficients in N ∪ {+∞}, such that Moreover, let for any k ∈ N, β k (X; K) the k-th Betti number of the manifold X with respect to the field K and denote by M( f , k) the number of critical points x of f , such that the Morse index m(x, f ) is equal to k. Then the Morse inequalities hold: Finally, denote by B(X, K) the total Betti number of f with respect to the field K, that is Then the total Betti number formula on the number K( f ) of the critical points of f holds: We point out that, under the assumption of the previous theorem, it can be proved the set of the critical points of the functional f is countable, because nondegenerate critical points are isolated and (PS) condition holds, so that Equation (32) is well defined.

Remark 6.
We point out that in order to state the Morse relations (32), we need that the functional f is of class C 2 , for instance in order to define the Morse index. In the next section we shall see that the action integral of a Finsler is only of class C 1,1 , that is of class C 1 with locally Lipschitz gradient. In spite of this, we shall see that Morse relations hold also for Finsler metrics.

Variational Theory for Geodesics on a Finsler Manifold
In this section, we state some results on the existence and multiplicity of geodesics joining two points on a Finsler manifold. We refer to [17] for the details of the proof using the infinite dimensional variational approach and general boundary conditions. Moreover, we state the Morse relations for geodesics joining two non-conjugate points, see [42,43]. The variational approach was already introduced by Franco Mercuri in the study of closed geodesics on a compact Finsler manifold, see [44]. An approach using the geodesic shortening flow was developed by Matthias in [45], and finally a variational approach for reversible Finsler metrics was previously developed by Kozma, Kristaly, and Varga, see [46].
We fix the interval I = [0, 1]. Moreover, let n ∈ N, n ≥ 1, we shall consider the Sobolev space H 1,2 (I, R n ) of the absolutely continuous curves x(s) : I → R n , such that the derivativeẋ(s) ∈ L 2 (I, R n ), see for instance [47]. It is well known that the space H 1,2 (I, R n ) has a structure of Hilbert space with norm We consider also the subspace H 1,2 0 (I, R n ) of H 1,2 (I, R n ) of the curves x(s), such that x(0) = x(1) = 0. On the subspace H 1,2 0 (I, R n ) an equivalent norm is given by The Sobolev manifolds were first introduced in the fundamental paper in [9], and it was also proven that H 1,2 (I, M) is equipped of a structure of infinite dimensional Hilbert manifold modelled on the Hilbert space H 1,2 (I, R n ). For any x ∈ H 1,2 (I, M), the tangent space T x H 1,2 (I, M) to H 1,2 (I, M) at the curve x is given by where TM is the tangent bundle of M and π = TM → M is its bundle projection. In other words, an element ξ(s) ∈ T x H 1,2 (I, M) is a continuous vector field along the curve x(s), such that, locally, ξ(s) is of class H 1,2 .
The Sobolev manifold H 1,2 (I, M) admits two important submanifolds for the study of variational problems of geodesics. Both submanifolds are still infinite dimensional. The first consists of curves joining two fixed points. Let p and q on M and set It can be proved, see for instance [48] that Ω 1,2 (p, q; M) is a smooth infinite dimensional submanifold of H 1,2 (I, M) and the tangent space at a point x ∈ Ω 1,2 (p, q; M) is given by Additionally, Λ 1,2 (M) is an infinite submanifold of H 1,2 (I, M) and the tangent space at a point x ∈ Λ 1,2 (p, q; M) is given by We have introduced the Sobolev manifolds where the variational problems are well defined for geodesics on a Riemannian manifold, or more generally on a Finsler manifold and a Lorentzian manifold. We state now some results on the existence and multiplicity for geodesics on a Finsler manifold. Moreover we state the Morse relations for the nonconjugate points. We fix a Finsler manifold (M, F) as in Definition 1, we consider the action integral E(x) : H 1,2 (I, M) → R, defined on the Hilbert-Sobolev manifold of the H 1 -curves on the manifold M. The action functional E is of class C 1,1 , that is, it is of class C 1 and its differential is locally Lipschitz. Moreover, if the Finsler metric F is the square root of a Riemannian metric, then E(x) is of class C 2 , see for instance [49].
The following variational principles hold which characterise the geodesics joining two points or the closed geodesics joining two points, see [17]. Theorem 9. Let (M, F) be a Finsler manifold and let p and q be two points of M, then geodesics joining p and q for the Finsler metric F are the critical points of the action integral E(x) restricted to the submanifold Ω 1,2 (p, q; M) of H 1,2 (I, M). Moreover the closed geodesics for F are the critical points of the restriction of E(x) to the sumbmanifold Λ 1,2 (M) of H 1,2 (I, M). In both cases the geodesics satisfy a variational principle, indeed if x(s) : [0, 1] → M is a geodesic joining two fixed points p and q, or it is a closed geodesic, then there exists a non-negative constant E x , such that for any s ∈ [0, 1], The previous variational principle holds for more general boundary conditions, see [17]. We point out that the constant E x is null if, and only if, x is a trivial geodesic, that is, it is reduced to a point.
The previous theorem and the abstract critical theorems stated in Section 6 allow to produce some results on the existence and multiplicity for geodesics joining two points on a Finsler manifold, satisfying some completeness assumptions.
We first state an existence result for geodesics joining two points. A proof of the previous can be found in [17], where it is shown that the action integral E(x), which is obviously bounded from below, satisfies the (PS) condition, under the assumption that F is forward or backward complete.
In order to state a multiplicity result for the geodesics joining two points, the topological richness of the infinite dimensional manifold Ω 1,2 (p, q; M) has to be studied. This was performed by Fadell and Husseini in the two papers [50,51], which study the topological properties of the based loop space of a manifold. The main results of Fadell and Husseini are stated in the following theorem: Theorem 11 (Fadell-Husseini). Let M be a smooth and non-contractible manifold and let p and q two points of M, then the Ljusternik-Scnirelmann category of the Sobolev manifold Ω 1,2 (p, q; M) is equal to +∞. Moreover, there exists a sequence (K m ) m∈N of compact subsets of Ω 1,2 (p, q; M), such that lim m→∞ cat Ω 1,2 (p,q;M) (K M ) = +∞.
By applying the previous theorem to the restriction of the action integral on Ω 1,2 (p, q; M) we attain the following multiplicity result on the geodesics joining two points on a Finsler manifold, see [17]. We examine now the Morse theory for geodesics on a Finsler manifold. Now there are remarkable differences between the Riemannian geometry and Finsler geometry. Indeed, in Riemannian geometry the action functional (10) is of class C 2 and this fact permits the finding of the classical Morse Relations developed by Palais in [9], see also [34] for a proof using the classical curve shortening procedure of Marston Morse.
Instead, for a generic Finsler metric the action integral is only of class C 1,1 and so we cannot apply Theorem 8 to obtain the Morse relations for geodesics. In spite of this, the Morse relations can be obtained also in the case of Finsler manifolds.
We now start the Riemannian case. Let (M, g) be a Riemannian manifold and let p and q two points of M. The action integral E(x) : Ω 1,2 (p, q; M) → R is of class C 2 (see for instance [48]). Let x be a critical point of E on Ω 1,2 (p, q; M) (that is x is a geodesic joining p and q), then the Hessian E (x) : T x Ω 1,2 (p, q; M) × T x Ω 1,2 (p, q; M) → R ig given for any ξ, ξ ∈ T x Ω 1,2 (p, q; M), by where ∇ is the Levi-Civita connection associated to the Riemannian metric g, R(·, ·)· denotes the Riemann tensor for the metric g, see also (15). Now assume that p and q are non-conjugate, that is there are no non-trivial ξ(s) ∈ T x Ω 1,2 (p, q; M) solutions of the Jacobi equations D 2 ξ + R(ẋ, ξ)ξ = 0, then the action integral is a Morse function, that is all its critical points are non-degenerate. Moreover, it can be proved that any geodesic joining p and q has finite Morse index m(x, E) and co-index M * (X, E), because the Hessian E (x) is a perturbation of the positive definite bilinear form Q(x)[ξ, ξ ] = Putting together these results and applying Theorem 8, we obtain the Morse relations for the geodesics joining two non-conjugate points of a Riemannian manifold. Theorem 14. Let (M, g) be a complete Riemannian manifold and let p, q ∈ M be two nonconjugate points for the metric g.
Then, for any field K there exists a formal series Q(r) whose coefficient belong to N ∪ {+∞} such that ∑ x∈K(E) r m(x,E) = P (Ω 1,2 (p, q; M)) + (1 + r)Q(r), where we have denoted by K(p, q) the set of the geodesics joining p and q.
We consider the case now of a Finsler manifold (M, F), the action integral is only of class C 1,1 on the Sobolev manifold H 1,2 (I, M), see, for instance, [17,44,49], so we cannot obtain the Morse relations for geodesics directly from the abstract Theorem 8. However the Morse relations for the geodesics joining two non-conjuagte points still hold. Fix two points p and q in the manifold M and consider the restriction of the action integral E to the manifold Ω 1,2 (p, q; M), so that the critical points of E on Ω 1,2 (p, q; M) are the geodesics joining p and q. Fix a geodesic x ∈ Ω 1,2 (p, q; M), we can define, at least in a formal way, the Hessian E (x) : T x Ω 1,2 (p, q; M) × T x Ω 1,2 (p, q; M) → R, setting, for any ξ, ξ ∈ T x Ω 1,2 (p, q; M), see [28] where g(ẋ,ẋ) is the fundamental tensor associated to the Finsler metric F, defined at (4), D s denotes the Chern connection with referenceẋ of ther Finsler metric F and R(x,ẋ) denotes the flag curvature associated to the Finsler metric F. The Hessian can be defined first for smooth vector fields ξ along the geodesic x and then can be extended by density on T x Ω 1,2 (p, q; M). Now, as for Riemannian metrics, the Hessian E (x) is a bilinear form on the tangent space T x Ω 1,2 (p, q; M), and the notions of non-degenerate geodesics, Morse index and co-index, Jacobi equations and conjugate points can be extended to Finsler metrics. In particular any geodesic has finite index m(x; E) and coindex m * (x, E), because E (x) is a bilinear form which is a compact perturbation of a positive definite bilinear form. Moroever, as for the Riemannain case, the Morse index m(x, E) of a geodesic is equal to the geometric index µ(x), the number of conjugate points counted with their multiplicity.
Using a fundamental result by Richard Palais, see [37], on the homotopical equivalence between the Sobolev manifold H 1,2 (I, M) and the Banach manifold C 0 (I, M) of the continuous curves on M, the following result has been proven in [42,43].
Theorem 15. Let (M, F) be a Finsler manifold, and assume that F is forward (or backward) complete. Let p and q on M and assume that they are non-conjugate, that is any geodesic x joining p and q is a non-degenerate critical point of the action integral E.
The Morse relations for the geodesics joining two non-conjugate points of a forward (or backward) Finsler manifolds hold as in Riemannian geometry. A consequence of all the previous results and which has an interesting applications in general relativity is the following theorem on the number of geodesics joining two non-conjugate points on a Finsler manifold (obviously it holds also for Rimennian manifolds). Theorem 16. Let (M, F) be a forward (or backward) complete Finsler manifold and let p and q be two non-conjugate points for F. Then the number of geodesics joining p and q is infinite or it is an odd number.
Indeed, if the manifold M is non-contractible, the Ljusternik-Schnirelmann category of Ω 1,2 (p, q; M) is equal to +∞ and so that there exists infinitely many geodesics joining p and q. Assume that M is contractible into itself, then also the Hilbert manifold Ω 1,2 (p, q; M) is contractible into itself, so its Poincaré polynomial is equal to 1. Setting r = 1 in (40), we obtain: If Q(1) ∈ N ∪ {+∞} is equal to +∞, then the set of the geodesics joining p and q is still equal to +∞. However, if Q(1) is finite, then the number of geodesics joining p and q is equal to 1 + 2Q(1), so it is an odd number. A condition which guarantees that the number of geodesics is finite is the existence of a function L : M → R, which is strictly convex with respect to the Finsler structure F, and with only one non-degenerate at the absolute minimum, see [52,53].
We conclude this section recalling some results about the existence of closed geodesics. Additionally, for this problem, there are deep differences between the Riemannian and the Finslerian case. A long standing conjecture about Riemannian geometry was the existence of infinitely many, geometrically closed geodesics on a compact Riemannian manifold, see for instance [48]. Two closed geodesics are geometrically distinct if they have different image on the manifold. This conjecture is still open for a generic Riemannian manifold, but it was solved, for instance, on the two dimensional sphere S 2 , where Bangert and Franks, see [54,55], proved that any Riemannian metric admits infinitely geometrically distinct closed geodesics. However, this result is false for Finsler manifolds, by the very famous Katok counterexamples [20,56], which are perturbations of the standard metric on S 2 admitting only 2 geometrically closed geodesics. The problem on the multiplicity of closed geodesics on Riemannian and Finsler manifolds is still a very active research field, see for instance [57], which gives also a survey on the closed geodesic problem on Finsler manifolds.

An Alternative Variational Principle for Geodesics for a Randers Metric
We have seen in the last section that one of the main differences between Riemannian geometry and Finsler geometry is that the regularity of the action integral on the Hilbert-Sobolev manifold H 1,2 (I, M) and its submanifolds. It is of class at least C 2 in Riemannian geometry and it is of class C 1,1 in Finsler geometry. However, for the geodesics of Randers metrics, it can be proved an alternative variational involving a functional of class C 2 . Let (M 0 , F), be a Randers metric, with (7). The action integral for a Randers metric has the form and it clearly of class C 1,2 on the H 1,2 (I, M 0 ). Now consider the functional The functional I(x) is of class C 1,1 on H 1,2 (I, M 0 ) and it is of class C 2 on the open submanifold of non-constant curves. In particular I(x) is of class C 2 on the submanifold Ω 1,2 (p, q; M 0 ), if p = q, while it is of class C 1,1 on the submanifold Lambda 1,2 (I, M 0 ). The following alternative variational principle has been proven in [58]. We restrict to the case of two points p and q. Theorem 17. (M 0 , F) be a Finsler manifold, let F be a Randers metric and let p and q two points of M 0 . Then a curve x(s) ∈ Ω 1,2 (p, q; M 0 ) is a critical point of the functional I on Ω 1,2 (p, q; M 0 ) if, and only if, x is a pregeodesic for the Randers metric F and the parameterisation of x(s) satisfies the property that h(x)[ẋ,ẋ] is constant.
An analogous result holds for closed geodesics for the Randers metric F. By the previous theorem, the results on the existence and multiplicity for geodesics on a Randers can be obtained using the C 2 functional I(x), rather than the action integral E.

Multiplicity of Light Rays and Applications to the Gravitational Lensing
In this section, we shall apply the results of the previous sections to obtain some results on the existence and multiplicity of light rays on a standard stationary Lorentzian manifolds. In particular the multiplicity results allow to prove under which geometric conditions the physical phenomena of multiple images to the gravitational lens effect can occur.
We fix a standard stationary Lorentzian manifold (M, g), with M = M 0 × R and g is given by (18). Let p = (x 0 , 0) be a point of M (the choice of the time coordinate equal to 0 is not restrictive) and a vertical line γ(s) = (x 1 , s) : R → M on the manifold M, with x 0 = x 1 . The following results on the existence and multiplicity of light rays joining p and γ hold, see [17]. Theorem 18. Assume that the Fermat metric F + , see (23) associated to (M, g) is forward complete, then there exists a future pointing light ray z(s) = (x(s), t(s)) joining p and γ, such that the travel time T(z) = t(1) = F + (x).
Moreover assume that the manifold M is non-contractible into itself, then there exists a sequence z m = (x m , t m ) of future pointing light rays joining p and γ with arrival times T(z m ), such that lim m→∞ T(z m ) = F + (x m ) = +∞.
An analogous result holds for past pointing light rays. Theorems 18 and 19 can be proven by applying the critical point theory on geodesics on Finsler metrics to the Fermat metrics F + and F − . Indeed, by the Fermat principle, the existence and multiplicity of the future pointing (respectively, past pointing) light rays joining p and γ is equivalent to the same problem for the geodesics of the Fermat metrics F + (respectively, F − ), so we can apply the critical point theory to the action functionals defined on the Sobolev manifold Ω 1,2 (M 0 , x 0 , x 1 ). The functionals are of class C 1,1 on Ω 1,2 (M 0 , x 0 , x 1 ), are bounded from below and satisfies the Palais-Smale condition, assuming the forward completeness of the Fermat metrics F + (respectively, F − ).
The multiplicity result follows by the fact that if the manifold M, or equivalently the manifold M 0 is non-contractible, then the Ljusternik-Schnirelmann category oh the Sobolev manifold Ω 1,2 (M 0 , x 0 , x 1 ) is equal to +∞, and so we can apply Theorem 7. Notice that these results can be obtained also using the alternative variational principle for geodesics for a Randers metric, see [58,59].
This multiplicity result allows to obtain a mathematical description of the so called gravitational lensing effect in a standard stationary spacetime, that is multiple images of a light source in the spacetime can be detected by an observer, if the spacetime is non-contractible.
We state now the Morse relations for the light rays joining p and γ, assuming that p and γ are not light-conjugate. We still consider a standard stationary Lorentzian manifold (M, g), with M = M 0 × R and g is given by (18) There is a strict relation between the light-non-degeneracy of the point p = (x 0 , 0) and the vertical line γ(s) = (x 1 , s) on the manifold M, with the non-degeneracy of the points x 0 and x 1 with respect to the Fermat metrics F + and F − on F − . Indeed, the following theorem has been proved in [42], see also [59] for a proof involving the alternative variational principle. Theorem 20. The point p = (x 0 , 0) and the vertical line γ(s) = (x 1 , s) are light-non-conjugate if, and only if, x 0 and x 1 are not conjugate with respect to the Fermat metric F + . Moreover, let z(s) = (x(s), t(s)) be a light-like geodesic joining p and γ, then the Morse index of the reparameterisation of x(s) which is a geodesic for the Fermat metric F + , is equal to the geometric index µ(z) of the light-like geodesic z(s), that is the number of conjugate points of z, counted with their multiplicity.
Thanks to the previous theorem, we can obtain the Morse relations for the light rays joining p and γ, assuming that p and γ are light-conjugate, by applying the Morse relations to the action integral of the Fermat metric F + for the non-conjugate points x 0 and x 1 . We shall denote by L 1,2 + (p, γ; M) the Hilbert manifold of the future pointing light-like curves joining p and γ. Clearly L 1,2 + (p, γ; M) is homotopically equivalent to Ω 1,2 (x 0 , x 1 ; M 0 ).

Theorem 21.
Assume that the Fermat metric F + is forward complete and assume that the point p and the vertical γ are light-non-conjugate and let Z be the set of future pointing light-like geodesics joining p and γ.
Then, for any field K, there exists a formal series Q(r) with coefficients in N ∪ {+∞}, such that ∑ z∈Z r µ(z) = P (L 1,2 An analogous result hold for past pointing light-like geodesics joining p and γ, assuming that F − is forward complete. The proof follows from Theorems 15 and 20. Indeed we have the Morse relations for the action integral for the Fermat metric F + and the points x 0 and x 1 and then we can replace the set of the geodesics joining x 0 and x 1 with the set Z of the future pointing light-like geodesics joining p and γ. Moreover, the Hibert manifolds L 1,2 + (p, γ; M) and Ω 1,2 (x 0 , x 1 ; M 0 ) have the same Poincaré series, because they are homotopically equivalent.
In particular, applying Theorem 16 to the Fermat metric, we produce the following on the number of light rays joining p and γ, assuming that they are light-non-conjugate Theorem 22. Assume that p and γ are light-non-conjugate, then the number of future pointing light rays joining p and γ are infinite or it is an odd number.
The previous theorem gives a theoretical proof of the odd image theorem, which claims that in the gravitational lens effect the number of images is odd. We see that this is true if the number of images is finite. We point out that the assumption that p and γ are light-non-conjugate is true except for a residual set. This follows by the Sard theorem.
Geometric conditions which gurantee that the number of images is finite has been studied in [52,53].
The results obtained on the existence and multiplicity of light rays joining the point p and the curve γ can be extended to the case the manifold M has a boundary. This case is physically very relevant, since many spacetimes in general relativity, as the Scwarzschild, the Reissner-Nordstrom, Kerr, and Kerr-Newman spacetimes have a boundary, or are a stationary structure only on some open subset of the spacetime. The boundary has to satisfy the following condition of light-convexity introduced in [59]. If the boundary ∂N of N is a smooth submanifold of M, then the light-convexity can be expressed by the Hessian, evaluated along light-like vectors, of a smooth function Φ : M → R, such that N = Φ −1 (]0, +∞[) and ∂N = Φ −1 (0). Light-convexity can be described also using the second fundamental form of the submanifold ∂N , see [60] for the details and other applications to general relativity. In [59] it has been proved that some open subsets of the Kerr spacetimes outside the stationary limit surface have lightconvex boundary. The light-convexity of open stationary subsets in a stationary spacetime can be related to the geometric properties of the Fermat metrics, see [61]. The same results hold for past pointing light rays, assuming that the Fermat metric F − is forward complete on N 0 ∪ ∂N 0 .
Applications of this theorem are presented in the papers [28,59,62,63], where it is shown that some open subsets of physically relevant spacetimes of general relativity as the Schwarzschild, the Reissner-Nordstrom, the Kerr, and the Kerr-Newman spacetimes contain non-contractible and open subsets with light-convex boundary, so that multiplicity of light rays hold and the gravitational lensing effect occur. In all these cases, in particular the Schwarzschild spactime, the manifold are non-contractible, and infinitely many images occur. Conditions which guarantee the finiteness of images can be found in the papers [53,64].
We conclude this section presenting another interesting class of light rays in a standard stationary Lorentzian manifold, the spatially closed light rays, which are light-like geodesics z(s) = (x(s), t(s)) : [0, 1] → M, such that the spatial component x(s) is a closed curve, that is x(0) = x(1). If we assume without any restriction that t(0) = 0, then z is called a T-spatially periodic light rays if t(1) = T. These classes of light rays are very interesting, because they are a generalisation to an arbitrary stationary spacetime of circular geodesics in static or stationary axis symmetric spacetimes, where circular geodesics define the light sphere and the shadow of a black hole, see for instance [32]. Additionally, spatially closed light rays satisfy a Fermat principle, because their spatial parts have to be closed pregeodesics for one of the Fermat metrics F + or F − , so the existence of spatially closed light rays can be studied searching for non-constant critical points of the action integral for the Fermat metrics on the manifold Λ 1,2 (M 0 ) of the closed curves. Some results on the existence of spatially closed light rays have been obtained in [60,65,66].

Causal Properties of a Stationary Spacetime and Fermat Metrics
In this section we shall see how the Fermat metrics influence the causal structure of standard stationary spacetime. We first recall some notions from causality theory of Lorentzian manifolds, we refer to [29][30][31] for details.
Let (M, g) be a Lorentzian manifold, two points p and q on (mathcalM are said causally related, there exists a curve γ : [0, 1] → M, such that γ(0) = p, γ(1) = q and γ is causal, that isγ(s) is a causal vector for any s ∈ [0, 1], section 3 for the definition of causal curve. Let p ∈ M, the causal future I + (p) is the subset of M given by the points q such that there exists a causal curve joining p and q. The causal past I − (p) is the subset of M given by the points q, such that there exists a causal curve joining q and p. We introduce now the notion globally hyperbolicity of a Lorentzian manifolds. This notion was introduced by Jean Leray in the 1950s in order to study the well posedness of strictly hyperbolic systems of partial differential equations. The notion of global hyperbolicity is a fundamental notion in general relativity in order to study evolution equations. Definition 20. A Lorentzian manifold (M, g) is said globally hyperbolic if (M, g) is causal, that is, it does not contain causal closed curves, and for any couple of causally related points ts p and q, the intersection of the future I + (p) of p with the past I − (q) of q is a compact subset of M.
It is possible to prove that a Lorentzian manifold (M, g) is globally hyperbolic if, and only if, it contains a Cauchy hypersurface M 0 , that is a subset of M such that any inextendible causal curve z of the manifold intersects M 0 exactly one, and only one, time.
Global hyperbolicity is an important assumption in the study of evolution equations in general relativity and, in particular, for Einstein equations.
The global hyperbolicity of a standard stationary are deeply related to the completeness of the Fermat metrics. Indeed the following result was proved in [17]. Conversely, if one of F + or F − is complete in a direction, then (M, g) is globally hyperbolic.
A necessary and sufficient condition for the global hyperbolicity of (M, g) has been successively proved in [67].

Theorem 25.
A standard stationary spacetime (M, g) is globally hyperbolic if, and only if, for any p ∈ M and for any r > 0, the intersection B + (p, r) ∩ B − (p, r) is a compact subset of M, where B + (p, r) and B − (p, r) are, respectively, the forward balls centred at p of radius r with respect to the Fermat metrics F + and F − .
In particular it follows that if the standard spacetime is globally hyperbolic, all the results on the existence, multiplicity and Morse theory joining a point p and a vertical line γ(s) = (x 1 , s) hold.

Geodesic Connectedness of a Lorentzian Manifold
A Lorentzian manifold (M, g) is said geodesically connected if for any couple of points p and q there is a geodesic for g joining p and q. Geodesics joining p and q are the critical points of the action functional E(z) = 1 0 g(z(s))[ż,ż]ds defined on the Hilbert manifold Ω 1,2 (M; p, q). Since the metric g is not positive definite, the action integral E is not bounded not from below nor from above, so a geodesic joining p and q cannot be found by minimisation of the functional E as for complete Riemannian manifolds. Moreover, it can be proven that all the stationary of the action functional E are saddle points of infinite Morse index and co-index, and this fact creates many difficulties to apply the critical point techniques developed in the last years to find a critical point of E.
The problem has been widely studied in a standard stationary spacetime, using critical point theorems for strongly indefinite functionals, or reducing the problem to the search of the critical points of E on a natural constraint, on which the action functional is bounded from below, which exists thanks to the stationarity of the metric. Let M, g be a standard stationary Lorentzian manifold, with M = M 0 × R and g given by (18). The action integral has the form defined on the Sobolev manifold H 1,2 (I, M) = H 1,2 (I, M 0 ) × H 1,2 (I, R). Let z 0 = (x 0 , t 0 ) and z 1 = (x 1 , t 1 ) be two points of M, then the geodesics joining z 0 and z 1 for the metric g are the critical points of the action integral on the Sobolev submanifold Ω 1,2 (z 0 , z 1 ; M) = Ω 1,2 (x 0 , x 1 ; M 0 ) × Ω 1,2 (t 0 , t 1 ; R). From (43), it immediately appears that the action integral E is unbounded both from below and from above, so a critical point cannot be found by a minimum argument. In spite of this the problems of the existence and multiplicity of the geodesics joining z 0 and z 1 have been intensively studied, and a quite definitive answer can be given to this problem. Here we present the main results obtained about these problems

•
The first result about the geodesic connectedness of a standard stationary Lorentzian manifolds was obtained by Vieri Benci and Dino Fortunato in [68], assuming that the spatial hypersurface M 0 is compact and using techniques of strongly indefinite functionals; • Successively the geodesic connectedness was proved in [69], assuming that the Riemannian manifold (M 0 , ·, · ) is complete and some assumptions on the growth of the shift vector field δ(x) and the lapse function β(x) are satisfied. In [69] a natural constraint for the action integral E is proven and it allows to prove that the existence of the geodesics joining z 0 = (x 0 , t 0 ) and z 1 = (x 1 , t 1 ) is equivalent to the search of the critical points of a new functional J defined only on the submanifold Ω 1,2 (x 0 , x 1 ; M 0 ). The functional J : Ω 1,2 (x 0 , x 1 ; M 0 ) → R is defined in the following way: Under the growth assumptions on δ(x) and β(x), it can be proved that J is bounded from below. Moreover, the completeness of the Riemannian manifold (M 0 , ·, · ) guarantees that J satisfies the Palais-Smale condition, so a minimum of the functional J exists and it corresponds, by the variational principle, to a geodesic joining z 0 and z 1 . Finally, if the manifold is non-contractible, Theorem 11 guarantees the existence of infinitely many critical points of the functional J and also the existence of infinitely many geodesics joining z 0 and z 1 . Moreover in [69], the case of manifold with convex boundary is also considered; • In [70], it is studied the geodesic connectedness of a non-necessarily standard stationary spacetime; • Finally in [71] it is proved the geodesic connectedness of globally hyperbolic standard stationary Lorentzian manifold, assuming that the Riemannian manifold (M 0 , ·, · ) is complete, see (18). Now, by the characterisation of the global hyperbolicity of a standard stationary Lorentzian manifold in terms of the Fermat metrics, we state the following proposition, where the assumptions are given only on the initial conditions, the spatial component (M 0 , ·, · ), the shift δ(x) and the lapse β(x) and not on the whole stationary metric, i.e., the global hyperbolicity. Proposition 1. Let (M, g) be a standard stationary Lorentzian manifold as in (18) and assume that (M 0 , ·, · ) is a complete Riemannian manifold and one of the Fermat metrics F + or F − is forward complete. Then (M, g) is geodesically connected.

Finsler Geometry and the Sagnac Effect
In this section we present another application of Fermat metrics, a geometric description of the celebrated Sagnac effect using the Fermat metrics. We shall see that the Sagnac effect is a measure of the asymmetry of the Fermat metrics F + and F − and, indeed, if the metric is static, the Sagnac effect does not appear. This observation is not new and has already been exploited in many papers. Here we emphasise the role of the Fermat metrics in the Sagnac and we think that it strongly put in evidence the value of Finsler geometries in the study of physical phenomena. In this case, the application of Finsler geomtery to the Sagnac effect is not a mere generalisation of a phenomena geometrically described using Riemannian geometry, because in that case the Sagnac effect does not appear.
We present, now, the geometric description of the Sagnac effect. We refer to [72] for the original paper by George Sagnac, to [73] for the Paul Langevin interpretation of the Sagnac effect in special relativity, to [74] for a geometric approach to the Sagnac effect, where the role of the Fermat metrics is not exploited, to [75] for a more intrinsic approach to the Sagnac effect, and to [76] for technological applications to the the Sagnac effect as in the global positioning system (GPS). Let (M, g) be a standard stationary spacetime with M = M 0 × R and the metric is given by (18).
Let x(s) : [0, 1] → M 0 be a spatial closed curve, so x(0) = x(1) and take the lift z(s) = (x(s), t(s)) of x to a future pointing light-like curve, so g(z)[ż,ż] = 0 and t(s) solves (21), that iṡ Now let y(s) = x(1 − s) : [0, 1] → M 0 the closed curve obtained reversing the orientation of x(s) and let w(s) = (y(s), k(s)) the lift of y(s) to a future pointing light-like curve, so thatk (s) = g 0 (y)[ẏ,ẏ] + g 0 (y)[δ(y),ẏ] 2 + g 0 (y)[δ(y),ẏ] > 0. Now we evaluate the travel times of the two future pointing light-like curves z(s) and w(s) and we obtain On the other hand, since y(s) has been defined reversing the orientation of x(s), we have T(w) = F + (y) = The difference of travel times causes an optical diffraction, which is called the Sagnac effect along the closed spatial curve x(s).

Remark 7.
We point out that both the Fermat metrics, F + and F − , contribute to the evaluation of the Sagnac effect. Moreover if the metric is static, and there is only one Fermat metric g 0 which is Riemannian, then ∆T = 0 and the Sagnac effect does not appear.
In particular, in the Boyer-Lindquist coordinates the Sagnac effect does not appear on the Schwarzschild and Reissner-Nordstrom. On the Kerr spacetime outside the stationary limit surface, the Sagnac effect appear.
Let x(s) = (r(s), θ(s), φ(s)) : [0, 1] → M a be a closed curve whose support is contained on the spatial part of Kerr spacetime outside the stationary limit surface, then by (44) and (30)

Conclusions and New Directions
In this paper we have introduced the Fermat metrics associated to a standard stationary spacetime. They are two Finsler metrics of Randers type and their definition is obtained evaluating the travel time of light-like curve of the spacetime. We have recalled the Fermat principle which relates the light rays of a standard stationary spacetimes to the geodesics of the Fermat metrics and we have presented some applications to the study of the geometrical and physical properties of this class of spacetimes. In particular we have studied the multiplicity of light rays joining an event and a time-like curve, giving a mathematical description of the multiple image effect. We have studied some causal property of a standard stationary spacetime, in particular, we related the global hyperbolicity to the completeness of the Fermat metrics. A result on the geodesic connectedness has been obtained involving again the completeness of the Fermat metrics. Finally, a geometric formulation of the Sagnac effect with the Fermat metrics has been presented.
Those are only some of the possible applications of the Fermat metrics. Some other application of the Fermat metrics concerns the strict relations between the Randers-Finsler geometry and the Zermelo navigation problem, see [77][78][79] for the formulation of the Zermelo navigation problem and its link with Randers geometry.
Another very interesting application of Fermat metrics is the description of the gravitational lensing in Kerr-Randers geometry using the Gauss-Bonnet theorem, see [80] and also the recent paper [81]. Finally, in the paper [82], the case in which the Fermat metrics have no sign restriction is considered. Now we present other research fields in which some optical metrics, that we always define as Fermat metrics, appear. These research fields are actually strongly active and have some interesting physical applications.

Spacetimes with a Light-like Killing Field
We have introduced the Fermat metrics in the class of standard stationary Lorentzian manifolds, where the vector field Y(z) = (0, 1) is killing and time-like, see Remark 1. Recently it has been considered the case in which Y(z) = (0, 1) is a light-like killing vector field, see [83]. As in Section 4, let (M 0 , ·, · ) be a Riemannian manifold and let δ(x) be a smooth surface, such that δ(x) = 0, for all x ∈ M 0 . It can be proved that the following covariant symmetric tensor field g(z)[ζ, ζ ] = g(x, t)[(ξ, τ), (ξ , τ )] = ξ, ξ + δ(x), ξ τ + δ(x), ξ τ is non-degenerate and it defines a Lorentzian metric on the product M = M 0 × R. These class of metrics are often used to model gravitational waves. The differences between the standard stationary Lorentzian manifolds are that the lapse function β(x) is null and the shift δ(x) cannot be equal to 0 to any point. This fact imposes some topological restrictions on the manifold M 0 . Indeed, if M 0 is a compact manifold, then its Euler-Poincaré characteristic χ(M 0 ) must be different from 0. The vector field Y(z) = (0, 1) is still a killing vector field, because the coefficients of the metric g do not depend on the variable t, but now it is light-like, as it can be easily proved showing that g(x, t)[(0, 1), (0, 1)] = 0. The Lorentzian (M, g) is still time-oriented and a global time-like vector field is the gradient ∇T of the temporal function T(x, t) = t : M → R. The fact that the shift δ(x) is never null needs to prove that ∇T is time-like everywhere. Now, if we consider a future light-like curve z(s) = (x(s), t(s)) : [0, 1] → M, for the metric g, arguing as in Section 4, we obtain: g(z)[ż,ż] = ẋ,ẋ + 2 δ(x),ẋ ṫ = 0 and solving with respect toṫ, we obtain: The travel time of the light-like curve z(s) is given by So in this case the Fermat metric is given by the function which a positively homogeneous function of degree 1 as a Finsler metric, but it is not defined on the whole tangent bundle, because it is not defined on the set The metrics (46) are known in literature as Kropina metrics and are a very important class of the so called conical Finsler metrics, see [84,85], where Kropina metrics were introduced first.
A Fermat principle relating the light-like geodesics for the metric g defined at (45) and the geodesics for the Kropina metric (46) has been proved in [83]. Some results on the existence of geodesics connecting two points or closed geodesics for a Kropina metric have been recently proved in [86]. We also mention the paper [87], where the existence of geodesics joining two points for the Lorentzian metric (45) is studied. We point out that, since a Kropina metric is not defined on the whole tangent bundle, there always exists couple of point which are not joined by any smooth curve and by no geodesic, see [86]. Instead, the problem on the existence of a closed geodesic for an arbitrary Kropina metric on a compact manifold is still open.

Finsler Spacetimes
In the recent years, some alternative models of the gravitational field have been introduced, in order to justify some observational results in astronomy. The first alternative model consists in the modification of the Einstein-Hilbert functional, whose Euler-Lagrange equations are the Einstein equations, with functionals whose Lagrangian density is a function of the scalar curvature. Another model, involved especially in observations of an asymmetric behaviour of matter, requires the substitution of the classical metric tensor given by a Lorentzian metric, with a Lorentz-Finsler metric, that is an homogeneous function defined on the tangent bundle of a manifold, such that the fundamental tensor defined as in (4) is a non-degenerate bilinear form of index 1, that is a Minkowski scalar product. This research field is actually very hot and it has many open problems, as the same definition of a Lorentz-Finsler metric F. Indeed some authors consider 1-homogeneous functions F, as in classical Finsler geometry, while authors prefer to define 2-homogeneous functions F, viewing F a sort of Lagrangian density and avoiding deep regularity problems given by the choice of 1-homogeneous functions. We refer to the review papers [88][89][90] and to the recent papers [91][92][93][94][95] for the main definitions and recent developments of Lorentz-Finsler geometry and its applications to gravitational physics.
A Fermat principle for light rays on Lorentz-Finsler spacetime was stated in [96], while in the papers [97,98], the notions of static and stationary spacetimes are extended to Lorentz-Finsler geometry, see also [99,100] for further results in this direction. For these classes of Lorentz-Finsler a detailed discussion of optical geometry and an analysis on the geometric and physical assumptions which guarantee the multiple image effect is still an open problem.