1. Introduction
One hundred years ago, Albert Einstein completed the mathematical formulation of his revolutionary view of the gravitational interactions in terms of curved space-time. The spirit, elegance, and experimental successes of the original theory has captivated the international scientific community and the theory has been accepted as the
standard model for gravity [
1]. At its centennial, pushed by new observational evidences, the theory is at a dramatic crossroad. Its continuation as the reference gravitational framework will imply that the universe is mainly composed by exotic sources of matter and energy whose existence is purely inferred from their gravitational effects at the largest astrophysical and cosmological scales [
2,
3,
4,
5,
6,
7,
8]. However, if such sources are not detected in any direct way, then we might be facing a failure of one of the most original theories of the twentieth century. The important implications of the two opposed alternatives,
i.e., the search for unknown dark side constituents or the revision of gravitational theory, demand a careful scrutiny of the different possible scenarios. In this work, we consider the second case, namely, the situation in which the gravitational dynamics may depart from that predicted by Einstein’s theory of General Relativity (GR) at ultraviolet and infrared scales. We mainly focus on the latter regime.
Given the success of GR at relatively short scales (such as the Solar System, stellar models, or compact binary systems), the idea that modified dynamics could arise at larger scales has been investigated in much detail over the last years. Theories in which the gravitational action consists of more general combinations of curvature invariants than the pure Einstein-Hilbert term have been investigated with special emphasis [
9,
10,
11,
12,
13,
14,
15,
16,
17]. From these investigations it was soon noticed that the usual metric formulation of alternative theories of gravity is generically different from its Palatini (or metric-affine) counterpart (see [
18] for a recent review on the Palatini approach). Whereas the metric approach typically leads to higher-order derivative equations, in the Palatini formulation the resulting field equations are always second-order. The appealing character of the second-order equations of the Palatini formalism, however, is accompanied by certain algebraic relations between the matter fields and the affine connection, which is now determined by a set of equations coupled to the matter fields and the metric. The case of
theories is particularly useful to illustrate the differences between these two approaches. In the metric formulation, the object
behaves as a dynamical scalar field, which satisfies a second-order equation with self-interactions that depend on the form of the Lagrangian
. In order to have an impact at large astrophysical and cosmological scales, the scalar field
ϕ should have a very low mass, implying a long interaction range. It is well known, however, that light scalars do have an impact at shorter scales, where their presence is strongly constrained by laboratory and Solar System observations unless some kind of screening mechanism is invoked [
19,
20,
21,
22,
23]. In the Palatini case, a scalar-tensor representation is also possible, but with the scalar field satisfying an algebraic rather than a differential equation. One then finds that the scalar
ϕ turns out to be an algebraic function of the trace of the stress-energy tensor of the matter,
, which may lead, in models of late-time cosmic speed-up, to undesired gradient instabilities at various contexts, as has been shown by studies of cosmological perturbations [
24,
25] and atomic physics [
26,
27].
In this article we will review the
hybrid variation of these theories, in which the (purely metric) Einstein-Hilbert action is supplemented with (metric-affine) correction terms constructed à la Palatini [
28,
29]. Given that metric and Palatini
theories allow the construction of simple extensions of GR with interesting properties and, at the same time, suffer from different types of drawbacks, we initiated a program to establish bridges between these two seemingly disparate approaches hoping to find ways to cure or improve their individual deficiencies. For that purpose, in a number of works we considered a hybrid combination of metric and Palatini elements to construct the gravity Lagrangian and found that viable models sharing properties of both formalisms are possible. An interesting aspect of these theories is the possibility to generate long-range forces without entering into conflict with local tests of gravity and without invoking any kind of screening mechanism (which would however require that at the present time the cosmological evolution reduces to GR). The possibility of expressing these hybrid
metric-Palatini theories using a scalar-tensor representation simplifies the analysis of the field equations and the construction of solutions. In some sense, considering a theory like
means that we retain all the positive results of GR, represented by the Einstein-Hilbert part of the action
R, while the further “gravitational budget” is endowed in the metric-affine
component. In fact it is well known that metric-affine and purely metric formalisms coincide in GR,
i.e., considering the action
R. On the contrary, the two formalisms lead to different results considering more generic functions
[
18].
A related approach to study
theories that interpolate between the metric and Palatini families is that of the so-called C-theories proposed in [
30,
31]. There the spacetime connection is associated to the metric
that is conformally related to the spacetime metric
, but the relation may depend upon the a scalar curvature
. This framework contains the metric,
, and the Palatini,
, formalisms as special limits, and one also finds that even when
physically distinct theories are possible. For further studies on variations of variational principles see e.g., [
32,
33,
34,
35,
36,
37,
38].
Other extensions of the
framework modify the coupling of matter to gravity by defining an action which depends linearly [
39] or nonlinearly upon the matter Lagrangian [
40,
41,
42,
43,
44,
45,
46,
47], or its trace [
48,
49,
50]. The new couplings generally induce non-geodesic motion mediated by an extra force orthogonal to the four-velocity [
51], which may have nontrivial effects already in flat Minkowski space. Instabilities due to new nonlinear interactions within the matter sector are thus common in these theories [
52,
53]. We note, from this perspective, that in the hybrid metric-Palatini approach, considered in this article, it can be expected that such instabilities in the matter sector are absent, because the usual conservation laws are satisfied.
In this paper, we review the formulation and the main applications of hybrid gravity models in late-time cosmological and astrophysical scenarios. The article is organized in three main parts considering the general formalism, the cosmology and the astrophysical applications. In
Section 2, we start the discussion considering the action and the field equations of the hybrid metric-Palatini formalism. In particular, we discuss the scalar-tensor representation, the Cauchy problem, and more general hybrid theories than
.
Section 3 is devoted to hybrid-gravity cosmology. We derive the Friedmann equations, construct the related dynamical system, and briefly consider some solutions. Furthermore, we analyse the cosmological perturbations in order to understand structure formation in these theories. We focus on the evolution of perturbations in the matter dominated era and vacuum fluctuations relevant to inflation. The weak field behaviour that is crucial for the Solar system precision tests of gravity is considered in
Section 4, where we also discuss the galactic phenomenology and the astrophysical applications of hybrid gravity as an alternative to dark matter. In particular, we study the stellar dynamics and the theory of orbits, the generalization of virial theorem, the flat rotation curves of spiral galaxies, and the galactic clusters starting from the relativistic Boltzmann equation for collisionless systems of particles. The conclusions are drawn in
Section 5.
4. Astrophysical Applications
In this Section we develop methods to study dark matter phenomenology in the hybrid models; for studies of dark matter generated by metric
modification, see e.g., [
88,
89].
Hybrid gravity allows to address several issues related to dark matter dynamics ranging from galaxies to galaxy clusters. These self-gravitating structures can be probed by studying the motion of test particles (stars for galaxies and galaxies themselves for galaxy clusters) moving into a gravitational potential. The behaviour of rotational and dispersion velocities of such test particles can be explained within the framework of the gravitational potential derived from the theory. For example, the tangential velocity can be explicitly obtained as a function of the scalar field of the equivalent scalar-tensor description. The model predictions can be compared with samples of rotation curves of spiral galaxies and galaxy clusters, respectively. The possibility of constraining the form of the scalar field and the parameters of the model by using the stellar velocity dispersions is also analysed. Furthermore, the Doppler velocity shifts are also obtained in terms of the scalar field. Finally suitable generalizations of the virial theorem and the relativistic Boltzmann equation allow to construct a self-consistent theory for galaxy clusters. In conclusion, all the physical and geometrical quantities and the numerical parameters in the hybrid metric-Palatini model can be expressed in terms of observable/measurable quantities, such as the tangential velocity, the baryonic mass of the galaxy, the Doppler frequency shifts, the dispersion velocity, the geometrical quantities characterizing the clusters of galaxies respectively. These results open the possibility of testing the hybrid metric-Palatini gravitational models at the galactic or extra-galactic scale by using direct astronomical and astrophysical observations.
Let us start our considerations by dealing with the weak field limit of the theory.
4.1. The Weak Field Limit
It is of paramount importance to determine the post-Newtonian parameters of the theory as they determine the compatibility of the theory with the local precision gravity tests. For post-Newtonian analysis of the metric and Palatini
theories, see e.g., References [
55,
56], and for a unified analysis Reference [
90]. Here, in particular, we are interested in the parameter
γ that is basically the fractional difference of the Newtonian potentials in Equation (
81) in the limit where the cosmological expansion can be neglected,
.
To this end, we need to consider the perturbations of Equations (
19) and (
23) in a Minkowskian background. The usual procedure is to assume
, where
is the asymptotic value of the field far away from the local system (and should be given by the cosmological background solution), and to take a quasi-Minkowskian coordinate system in which
, with
. Then it is easy to see that the quadratic terms
and
in Equations (
19) and (
23) do not contribute to the linear order. The potential terms in Equation (
23) can be linearized as follows (in this subsection we denote
):
The zeroth-order term in this equation is due to the background, and can be absorbed into a coordinate redefinition. The coefficient of the first-order term can be interpreted as a mass squared. The linearized scalar field equation is thus given by
where, as usual, in this order of approximation we have neglected the time derivatives and the pressure terms.
The linearization of the metric field equations is a bit more complicated because we need to establish suitable gauge conditions. Since the background is Minkowskian, the perturbed Ricci tensor is given by
where
. The term
on the righthand side of Equation (
19) can be combined with the terms
to give the following gauge conditions
With this choice, the linearized equations for the metric become
For consistency, to this order
,
,
, which leads to
Before solving Equations (
98), (
102), and (
103), it is worth noting that while the connection Equation (
14) is invariant under constant rescalings of the field
, the other field equations do not share this invariance. This is manifest in the combinations
in the above perturbation equations.
Using the generic solution
we find that
In spherical symmetry and far from the sources, the above equations become
where we have denoted
These results represent the standard post-Newtonian metric up to second order for this class of theories.
We emphasize a striking feature of
gravity. Note that in
gravity, to obtain
from Equation (
112), there is only one possibility [
55,
56], namely,
from millimetres to astronomical scales,
i.e., the range of the scalar interaction,
, should be smaller than a few millimetres. In the current case, however, there are two possibilities to obtain
. The first one is the same as in
theories and involves a very massive scalar field. The second possibility implies a small value
. If
, then the Yukawa-type corrections are very small regardless of the magnitude of
. This could allow for the existence of a long-range scalar field able to modify the cosmological and galactic dynamics, but leaving unaffected the Solar System. Subtle modifications could in the most optimistic case be detected as anomalies in the local gravitational field [
91].
4.2. Galactic Phenomenology: Stable Circular Orbits of Test Particles around Galaxies
The most direct method for studying the gravitational field inside a spiral galaxy is provided by the galactic rotation curves. They are obtained by measuring the frequency shifts
z of the 21-cm radiation emission from the neutral hydrogen gas clouds. The 21-cm radiation also originates from stars. The 21-cm background from the epoch of reionization is a promising cosmological probe: line-of-sight velocity fluctuations distort redshift, so brightness fluctuations in Fourier space depend upon angle, which linear theory shows can separate cosmological from astrophysical information (for a recent review see [
92]). Instead of using
z the resulting redshift is presented by astronomers in terms of a velocity field
[
93,
94].
In the following, we will assume that the gas clouds behave like test particles, moving in the static and spherically symmetric geometry around the galaxy. Without a significant loss of generality, we assume that the gas clouds move in the galactic plane , so that their four-velocity is given by , where the overdot stands for derivation with respect to the affine parameter s. In this subsection we find it illustrative to restore the units of c.
The static spherically symmetric metric outside the galactic baryonic mass distribution is given by the following line element
where the metric coefficients
and
are functions of the radial coordinate
r only. The motion of a test particle in the gravitational field with the metric given by Equation (
114), is described by the Lagrangian [
95,
96]
where
, which simplifies to
along the galactic plane
. From the Lagrange equations it follows that we have two constants of motion, namely, the energy
E per unit mass, and the angular momentum
l per unit mass, given by
and
, respectively. The normalization condition for the four-velocity
gives
, from which, with the use of the constants of motion, we obtain the energy of the particle as
From Equation (
116) it follows that the radial motion of the test particles is analogous to that of particles in Newtonian mechanics, having a velocity
, a position dependent effective mass
, and an energy
E. In addition to this, the test particles move in an effective potential provided by the following relationship
The conditions for circular orbits, namely,
and
lead to
respectively. Note that the spatial three-dimensional velocity is given by
For a stable circular orbit
, and the tangential velocity of the test particle can be expressed as
In terms of the conserved quantities, and along the galactic plane
, the angular velocity is given by
and taking into account Equation (
118), we finally obtain the following relationship [
95,
96]
Therefore, once the tangential velocity of test particles is known, the metric function
outside the galaxy can be obtained as
The tangential velocity
of gas clouds moving like test particles around the center of a galaxy is not directly measurable, but can be inferred from the redshift
observed at spatial infinity, for which
[
95,
96]. Due to the non-relativistic velocities of the gas clouds, with
, we observe that
, as the first part of a geometric series. The observations show that at distances large enough from the galactic center the tangential velocities assume a constant value,
i.e.,
constant [
93,
94]. In the latter regions of the constant tangential velocities, Equation (
123) can be readily integrated to provide the following metric tensor component
where
is an arbitrary constant of integration. If we match the metric given by Equation (
124) with the Schwarzschild metric on the surface of the galactic baryonic matter distribution, having a radius
,
, we obtain the following relationship
An important physical requirement for the circular orbits of the test particle around galaxies is that they must be stable. Let
be the radius of a circular orbit and consider a perturbation of it of the form
, where
[
97] . Taking expansions of
and
about
, it follows from Equation (
116) that
The condition for stability of the simple circular orbits requires
[
97] . Hence, with the use of the condition
, we obtain the condition of the stability of the orbits as
. By taking into account Equation (
122), it immediately follows stable circular orbits always exist for massive test particles.
4.2.1. Galactic Geometry and Tangential Velocity Curves in Hybrid Metric-Palatini Gravity
The rotation curves only determine one, namely
, of the two unknown metric functions,
and
, which are required to describe the gravitational field of the galaxy [
98]. Hence, in order to determine
we proceed to solve the gravitational field equations for the hybrid metric-Palatini gravitational theory outside the baryonic matter distribution. This allows us to take all the components of the ordinary matter stress-energy tensor as being zero. Furthermore, our task becomes easier when we restrict to perturbative weak-field treatment.
The weak field limit of the gravitational theories at the Solar System level is usually obtained by using isotropic coordinates as above in
Section 4.1. However, it is useful to apply Schwarzschild coordinates in studying exact solutions and in the context of galactic dynamics, and this is what we do here. We assume that the gravitational field inside the halo is weak, so that
, which allows us to linearise the gravitational field equations retaining only terms linear in
and again consider the scalar field as
, where
is a small perturbation around the background value
. The Klein-Gordon equation was already solved in
Section 4.1 at this limit and the result was that the interaction range is given by
, where
is given in Equation (
113). In a spherically symmetric configuration, the general solution then has the form
where
is an integration constant. Comparing this expression with the results obtained in [
28] for the weak-field limit (taking into account the transformation from isotropic to Schwarzschild coordinates), we find that8
where
and
are the mass and the radius of the galactic baryonic distribution, respectively.
Within this linear approximation the stress-energy tensor of the scalar field is given by
where
α and
β are defined by
Therefore the linearized gravitational field equations take the form
Equation (
131) can be immediately integrated to provide
where
is an integration constant. Comparing again with the results obtained in [
28] for the weak-field limit, we find that
. The tangential velocity of the test particles in stable circular orbits moving in the galactic halo can be derived immediately from Equation (
132), and is given by
which in terms of the solutions found above becomes
where
. The term proportional to
corresponds to the cosmological background, namely the de Sitter geometry, and we assume that it has a negligible contribution on the tangential velocity of the test particles at the galactic level.
On the surface of the baryonic matter distribution the tangential velocity must satisfy the boundary condition
which, with the use of Equation (
128), gives the following constraint on the parameters of the model,
In order to satisfy the above condition would require that
, or, equivalently,
In the regions near the galactic baryonic matter distribution, where
, we have
, to a very good approximation. Hence in this region the tangential velocity can be approximated as
If the parameters of the model satisfy the condition
the term proportional to
becomes negligible, while for small values of
, and
, the term proportional to
r can also be neglected. Therefore for the tangential velocity of test particles rotating in the galactic halo we obtain
Since according to our assumptions,
, the coefficient
α can be approximated as
, which provides for the rotation curve, in the constant velocity region, the following expression
Since
, the scalar field potential must satisfy the condition
. In the first order of approximation, with
, for the tangential velocity we obtain the expression
Alternatively, in general we can write the tangential velocity as follows,
As compared to our previous results, in this representation we have instead of . Since we are working in a regime in which , the choice of the constants or does not seem very relevant, since it just amounts to a rescaling of . From now on we will also assume that .
From the above equation we want to find the constraints on the model parameters that arise from the expected behavior at different scales. For that purpose, it is convenient to write the equation, equivalently, as follows:
At intermediate scales, the asymptotic tangential velocity tends to a constant. If we expand the exponential as
, hen we obtain the following three constraints on the free parameters of the model,
With increasing r, and by assuming that the condition still holds, the rotation curves will decay, at very large distances from the galactic center, to the zero value.
4.2.2. On Astrophysical Tests of Hybrid Metric-Palatini Gravity at the Galactic Level
In [
98], some observational astrophysical tests of hybrid metric-Palatini gravity at the galactic level were discussed. More specifically, a comparison of the theoretical predictions of the model with a sample of rotation curves of low surface brightness galaxies was analysed. Indeed, the predictions of the theoretical model with the observational results show that the contribution of the scalar field energy density to the tangential velocity of the test particles can explain the existence of a constant rotational velocity region around the baryonic matter, without requiring the presence of the dark matter. The possibility of observationally determining the functional form of the scalar field
φ by using the velocity dispersion of stars in galaxies, and the red and blue shifts of gas clouds moving in the galactic halo could also be worth considering. It can be shown that one can constrain the explicit functional form of the scalar field, and the free parameters of the model, in order to adequately fit the observational data. We refer the reader to [
98] for more details. Of course, was particle matter detected directly, we could exclude the gravitational
alternative to dark matter.
4.3. Galactic Clusters: The Generalized Virial Theorem in Hybrid Metric-Palatini Gravity
A large number of astronomical and astrophysical observations confirm that galaxies form complex hierarchical structures, in which galaxies concentrate in large objects, called clusters of galaxies, bounded by the gravitational interaction. The total mass of the galaxy cluster ranges from
for the so-called groups (formed by a few hundred of galaxies) up to a few
for very large clusters, containing thousands of galaxies. From a morphological point of view galactic clusters are usually by a main component, which is regular and centrally peaked [
99,
100]. For fundamental physics research the main importance of the galactic clusters consists in the fact that they are considered to be “dark matter” dominated astrophysical objects. Moreover, their formation and evolution is almost entirely controlled by the gravitational force, a property which allows the testing of different dark matter models. On the other hand the mass distribution inside the clusters is fully determined by the initial conditions of the mass distribution that originate in the early universe [
101,
102,
103], thus allowing the use of galaxy cluster properties to also test cosmological models. In this subsection we discuss the possibility of observationally testing the metric-Palatini gravity theory by using observational evidence from galaxy clusters [
104].
4.3.1. Galaxy Cluster As a System of Identical and Collisionless Point Particles
As a first step necessary to obtain our main result, consisting in the generalization of the relativistic virial theorem for galaxy clusters in the hybrid metric-Palatini gravitational we write down the gravitational field equations and the Boltzmann for a static and spherically symmetric distribution of matter. We adopt a simplified physical model, in which the galactic cluster consists of a self-gravitating system of identical, collisionless point particles (the galaxies) in random dynamical motion. The metric will thus be described by Equation (
114), and the fluid by a distribution function
that obeys the general relativistic Boltzmann equation [
105].
The energy-momentum tensor of the matter in the cluster, is thus determined by the distribution function
, and its components are given by the equation [
106,
107]
where
m is the mass of the particle (galaxy),
is the corresponding galactic four-velocity, with
denoting the temporal component. Finally, by
we denote the invariant volume element in the velocity space. Alternatively, the energy-momentum tensor
describing the matter distribution in a cluster of galaxies can be represented in terms of an effective energy density
, and of two effective anisotropic thermodynamic pressures, the radial
and the tangential
pressure, respectively. These thermodynamic parameters of the cluster are obtained by averaging over the matter and velocity distributions, and are given by
where
ρ is the mass density of the ordinary baryonic matter, and
,
denotes the average value of
,
, representing the square of the components of the four-velocities of the galaxies in the cluster [
108]. The full field equations for the metric Equation (
114) with the source Equation (
151) were listed in Reference [
104], but here we shall only need their trace:
where we have denoted
.
Since in the following we are interested only in astrophysical applications at the extra-galactic cluster scale, we will adopt a Newtonian type approximation, which consists in assuming that the deviations from standard general relativity (corresponding to the background value
) are small for the systems we are considering. This approximation implies that
. Thus, Equation (
152) can be approximated as
where
corresponds to an effective, geometric type “energy” of the scalar field in the hybrid metric-Palatini gravitational model.
4.3.2. The Relativistic Boltzmann Equation
Next we proceed to the second step in the derivation of the virial theorem for galaxy clusters and determine the general relativistic Boltzmann equation that governs the evolution of the galactic distribution function .
A basic result in statistical physics is the transport equation for the distribution function for a system of particles in a curved arbitrary Riemannian space-time. This transport equation is given by the Boltzmann equation without collision term, and which can be formulated as [
106,
107]
where
is the four-momentum of the galaxy (particle), and
are the Christoffel symbols associated to the metric. An important consequence of the collissionless Boltzmann equation is that the local phase space density, as measured by an observer in a frame co-moving with a galaxy, is conserved.
An important simplification of the mathematical formulation of the Boltzmann equation can be achieved by introducing an appropriately chosen orthonormal frame, or tetrad
,
. The tetrad fields vary smoothly over some coordinates neighborhood
U, and they satisfy the general condition
for all
, where
denotes the Minkowski metric tensor [
106,
107,
108]. A basic property of the tetrad fields is that any tangent vector
defined at an arbitrary point
x can be represented as
, a relation which defines the tetrad components
.
For the spherically symmetric line element given by Equation (
114), the frame of orthonormal vectors can be chosen in an appropriate way as [
106,
107,
108]:
The tetrad components of the four velocity are
. In the tetrad components corresponding to our present choices the relativistic Boltzmann equation Equation (
155) is given by
where the distribution function
is a function of
and
, respectively, and
are the Ricci rotation coefficients [
106,
107,
108]. Due to the spherical symmetry of the problem the distribution function depends only on the radial coordinate
r, and hence Equation (
157) becomes [
108]
Due to the spherical symmetry of our astrophysical system the coefficient of
in Equation (
158) must be zero. From a mathematical point of view this implies that the distribution function
is only a function of
r,
and
. As a next steps in our analysis we multiply Equation (
158) by
, and we integrate over the velocity space. Then, by taking into account that the distribution function
vanishes sufficiently rapidly as the velocities tend to
, we find the equation
Now we multiply Equation (
159) by
, and by integrating over the entire cluster volume, corresponding to a radius
R, gives [
108]
4.3.3. Geometrical Quantities Characterizing Galactic Clusters
In order to obtain some analytical estimations of the main geometrical and physical quantities characterizing the galactic clusters we introduce some approximations to the motion of both test particles in stable circular orbits around galaxies, as well as to the motion of galaxies in galactic clusters. As a first approximation we assume that
ν and
λ are slowly varying functions of the radial coordinate
r. Then in Equation (
153) we can neglect all the quadratic terms as being negligibly small as compared to the first order terms. Secondly, astronomical observations show that the motion of the galaxies in clusters is non-relativistic. Hence the galactic velocities are much smaller than the velocity of the light, that is,
. Thus, Equations (
153) and (
160) can be written as
and
respectively, where
is the total kinetic energy of the galaxies in the cluster. We define the total baryonic mass
of the galactic cluster as
We assume that the main contribution to the baryonic mass is due to the presence of the intra-cluster gas and of the stars. On the other hand we also include in the mass contribution of other particles, like, for example, massive neutrinos, which may also give a significant contribution to .
By multiplying Equation (
161) by
, and integrating from 0 to
r we obtain
where we have introduced the notation
It is interesting to note at this point that in hybrid metric-Palatini gravity, the quantity
has essentially a geometric origin. Hence it is natural to call it as the
geometric mass of the galactic cluster. In the following we introduce the gravitational potential energies of the cluster by means of the definitions
where
R is the cluster radius. By multiplying Equation (
165) with
, and by integrating from 0 to the cluster radius
R, we obtain the important relation
4.3.4. The Generalized Virial Theorem in Hybrid Metric-Palatini Gravity
As a last step in our analysis, with the help of Equation (
162), we obtain the generalization of the virial theorem in hybrid metric-Palatini gravity, which can be formulated in the familiar form
In the above equation the total gravitational potential energy of the system, Ω, defined as
contains a supplementary term
, which has a purely geometric origin.
It is useful to represent the generalized virial theorem, given by Equation (
170), in a more transparent physical form, which can be obtained by introducing the radii
and
, defined by
and
respectively. The quantity
, having a geometric origin similarly to the geometric mass considered above, may be called as the
geometric radius of the cluster of galaxies in the hybrid metric-Palatini gravity theory. Hence, finally we obtain the baryonic potential energy
and the effective scalar field potential energy
as given by
respectively. Another important observational quantity, the virial mass
of the cluster of galaxies is defined as follows
Eventually, the fundamental relation between the virial and the baryonic mass of the galaxy cluster can be obtained after substitution of Equation (
175) into the virial theorem as
If
, a condition which is satisfied by the astrophysical parameters of most of the observed galactic clusters, then from Equation (
176) we obtain the virial mass of galactic clusters in hybrid metric-Palatini gravity, which can be approximated as
The virial mass is determined observationally from the study of the velocity dispersion of the stars and of the galaxies in the galactic clusters. An important consequence of the virial theorem in hybrid metric-Palatini gravity is that in a cluster with mass most of the mass is in the form of the geometric mass . Hence we can use the approximation in order to study the cluster dynamics. A fundamental question related to the possibility of the observational testing of the hybrid metric-Palatini gravity is to find out what astrophysical processes can detect the presence of the geometric mass. Such an observational possibility may be provided by gravitational lensing. Through the study of the lensing properties of the galactic clusters one can obtain direct evidence of the existence of the geometric mass, of its distribution properties, as well as of the gravitational effects associated to the presence of the scalar field. It is interesting to note that gravitational lensing can give us theoretical information even at cosmical scales extending far beyond of the virial radius of the matter distribution of the galactic cluster.
4.3.5. On Astrophysical Tests of Hybrid Metric-Palatini Gravity at the Cluster Level
In concluding, in the framework of hybrid
gravity theory we have established the the existence of a strict proportionality between the virial mass of the cluster and its baryonic mass, a relation which can also be tested observationally. One of the important, and observationally testable, predictions of the hybrid metric-Palatini gravitational “dark matter” model is that the geometric masses associated to the clusters, as well as its gravitational effects, extend beyond the virial radii of the clusters. Observationally, the virial mass
is obtained from the study of the velocity dispersions of the stars in the cluster. Due to the observational uncertainties, this method cannot give a reliable estimation of the numerical value of the total mass
in the cluster. However, a much more powerful method for the determination of the total mass distribution in clusters is the gravitational lensing of light, which may provide direct evidence for the gravitational effects at large distances from the cluster, and for the existence of the geometric mass. The presence of hybrid metric-Palatini modified gravity effects at large distances from the cluster, and especially the large extension of the geometric mass, may lead to significantly different lensing observational signatures, as compared to the standard relativistic/dark matter model case. The bending angle in the hybrid metric-Palatini gravity models could be larger than the one predicted by the standard dark matter models. Therefore, the detailed observational study of the gravitational lensing could discriminate between the different theoretical models introduced to explain the motion of galaxies (“particles”) in the clusters of galaxies, and the standard dark matter models. We refer the reader to Reference [
104] for more details.
Finally, it is worth pointing out that hybrid gravity can be precisely tested also at smaller scales like those around the Galactic Centre. As reported in [
109], the observed peculiar orbit of S2 star, moving around the centre of our Galaxy, is theoretically reconstructed if one adopts the gravitational potential coming from hybrid gravity. This result opens new perspectives in achieving precision tests for the theory.
5. Conclusions
In this work we have presented a hybrid metric-Palatini framework for theories of gravity, and have tested the new theories it entails using a number of theoretical consistency checks and observational constraints. From the field theory perspective, we found that the
class of theories, where
, enjoys a similar special status amongst the more general hybrid metric-Palatini theories as the
theories within the narrower framework of purely metric gravity [
70]. This is so because when one excludes theories inhabited by ghost-like, superluminally propagating and otherwise pathological degrees of freedom, there is evidence, as shown in
Section 2.4, that the
family is singled out as the only viable form of an action one can construct using the metric (and thus the metric Levi-Civita connection) and an independent “Palatini connection”. The underlying reason is that in the special case of
actions the higher derivatives in the gravity sector can be separated into a scalar mode, thus avoiding an Ostrogradskian instability. Indeed the
gravity represents a generic case within the one-parameter family of the Algebraic Scalar-Tensor theories (recall Equation (
17)), at one end of which lies the pure Palatini
(wherein the scalar field is a function of the stress-energy trace
T) and at the other end the pure metric
(where the field is a function of the metric curvature
R). Furthermore, the propagating degrees of freedom have proven to be healthy also on curved backgrounds as confirmed also by our cosmological perturbation analysis in
Section 3.2. Concerning the Cauchy problem, it was shown that in this class of theories the initial value problem can always be well-formulated and well-posed depending on the adopted matter sources.
Having established the theoretical consistency and interest on the hybrid metric-Palatini
family of theories, we considered applications in which these theories provide gravitational alternatives to dark energy. As shown by our post-Newtonian analysis in
Section 4.1, the hybrid theories are promising in this respect as they can avoid the local gravity constraints but modify the cosmological dynamics at large scales. This is simply because as a scalar-tensor theory, the hybrid
gravity is characterised by an evolving Brans-Dicke coupling, which allows to introduce potentially large deviations from GR in the past (and future) as long as the coupling at the present epoch is strong enough to hide the field from the local gravity experiments. In contrast, in the metric
models the Brans-Dicke coupling is a finite constant and one needs to invoke some of the various “screening mechanisms” (workings of which remain to be studied in the hybrid theories) in order to reconcile the Solar system experiments with cosmology.
Cosmological perturbations have been also analysed in these models up to the linear order [
28,
54,
78], and the results imply that the formation of large-scale structure in the aforementioned accelerating cosmologies is viable though exhibits subtle features that might be detectable in future experiments. In
Section 3.2, we derived the full perturbations equations and extracting their Newtonian limit, describing the observable scales of the matter power spectrum, the growth of matter overdensities was shown to be modified by a time-dependent effective fifth force that is expected to modify the redshift evolution of the growth rate of perturbations. We also note that numerical studies of the perturbations imply that the difference of the gravitational potentials can exhibit oscillations at higher redshifts even when the background expansion and the full lensing potential are indistinguishable from the standard ΛCDM predictions [
78]. Such features could potentially be observed in cross-correlations of the matter and lensing power spectra, but detailed calculations of the cosmic microwave background anisotropies and other perturbation observables remain to be carried out. This is especially worthwhile in view of the potential of the forthcoming Euclid mission to experimentally test different accelerating cosmological models by their impact to the large-scale structure [
110].
At an effective level, the
modifications involve both (the trace of) the matter stress energy and (the Ricci scalar of) the metric curvature, and from this point of view it appears appealing to speculate on the possible relevance of these theories to both the problems of dark energy and dark matter, in a unified theoretical framework and without distinguishing
a priori matter and geometric sources. Various aspects of dark matter phenomenology from astronomical to galactic and extragalatic scales were discussed in
Section 4. The generalised virial theoreom can acquire, in addition to the contribution from the baryonic masses, effective contributions of geometrical origin to the total gravitational potential energy, which may account for the well-known virial theorem mass discrepancy in clusters of galaxies. In the context of galactic rotation curves, the scalar-field modified relations between the various physical quantities such as tangential velocities of test particles around galaxies, Doppler frequency shifts and stellar dispersion velocities were derived. More recently, observational data of stellar motion near the Galactic centre was compared with simulations of the hybrid gravity theory, which turned out particularly suitable to model star dynamics. Yet, to promote the
theory into a convincing alternative to particle dark matter, one should produce also the cosmological successes of the ΛCDM model without the CDM component.
Though wormhole solutions have already been found in these theories [
111], the nature of possible black hole solutions remains an outstanding open question. Though no-hair theorems for scalar-tensor theories exclude the simplest nontrivial solutions, their assumptions are restrictive. Another interesting though yet unstudied issue is the strong field regime of hybrid gravity and the constraints that one can put on the theories from astrophysical data such as measurements of binary pulsars. Finally, the possible relevance of the hybrid gravity framework to the early universe cosmology has not been explored in any detail. The Einstein frame formulation of the scalar-tensor description Equation (
25) could provide a useful starting point to investigate how the inflaton potential changes due to finite
, in order to understand how the “hybrid” nature of spacetime modifies the standard predictions of e.g., the
model of inflation.
To conclude, whilst the physics of the metric and the Palatini versions of
gravity have been uncovered in exquisite detail in a great variety of different contexts [
9,
10,
11,
12,
13,
14,
15,
16,
17,
18], those studies largely wait to be extended for the hybrid
version of the theory. We believe the results this far, as reported in this review, provide compelling motivation for the further exploration of these particular theories.