Mimetic Gravity: A Review of Recent Developments and Applications to Cosmology and Astrophysics

Abstract

Mimetic gravity is a Weyl-symmetric extension of General Relativity, related to the latter by a singular disformal transformation, wherein the appearance of a dust-like perfect fluid can mimic cold dark matter at a cosmological level. Within this framework, it is possible to provide a unified geometrical explanation for dark matter, the late-time acceleration, and inflation, making it a very attractive theory. In this review, we summarize the main aspects of mimetic gravity, as well as extensions of the minimal formulation of the model. We devote particular focus to the reconstruction technique, which allows the realization of any desired expansionary history of the universe by an accurate choice of potential or other functions defined within the theory (as in the case of mimetic gravity). We briefly discuss cosmological perturbation theory within mimetic gravity. As a case study within which we apply the concepts previously discussed, we study a mimetic Hořava-like theory, of which we explore solutions and cosmological perturbations in detail. Finally, we conclude the review by discussing static spherically symmetric solutions within mimetic gravity and apply our findings to the problem of galactic rotation curves. Our review provides an introduction to mimetic gravity, as well as a concise but self-contained summary of recent findings, progress, open questions, and outlooks on future research directions.

1. Introduction

The past decade has seen the astounding confirmation of the “dark universe” picture, wherein the energy budget of our universe is dominated by two dark components: dark matter and dark energy [1–50]. The race to determine the nature and origin of these components is in progress on both the observational and theoretical fronts. Theories of modified gravity appear quite promising in this respect, particularly given that gravity remains the least understood of the four fundamental forces. For an incomplete list of comprehensive reviews, as well as seminal works on the subject, we refer the reader to [51–68] and references therein.

A particularly interesting theory of modified gravity which has emerged in the past few years is mimetic gravity [69]. In mimetic gravity, as well as minimal modifications thereof, it is possible to describe the dark components of the universe as a purely geometrical effect, without the need of introducing additional matter fields. In the past three years, interest in this theory has grown rapidly, with over 90 papers following up on the original idea or at least touching on it in some way [70–161]. For this reason, we believe it is timely to present a review on the progress achieved thus far in the field of mimetic gravity. This review is not intended to provide a detailed pedagogical introduction to mimetic gravity but rather summarize the main findings and directions in current research on the theory, as well as providing useful directions to the reader, should she/he wish to deepen a given topic in mimetic gravity. For this reason, this review should not be seen as a complete introduction to mimetic gravity, nor should it substitute consultation of the original papers. No prior knowledge on the subject is assumed.

This review is structured as follows: in this section we will provide a historical and technical introduction to modified gravity, which shall justify our subsequent endeavour in mimetic gravity. In Section 2, we will introduce mimetic gravity. Given that understanding the reason behind the equations of motion of mimetic gravity differing from those of General Relativity is not obvious, a major goal of the section will be to provide a clear and concise explanation for this fact. Subsequently, in Section 3 we will explore a few solutions of mimetic gravity and correspondingly some extensions of the theory, as well as correspondences with related theories of modified gravity. Section 4 will provide a brief interlude focusing on perturbations in mimetic gravity. In Section 5 we will present a case study of a mimetic-like model, namely, mimetic covariant Hořava-like gravity, with a focus on its solutions and cosmological perturbations, and the need to extend the model beyond its basic formulation. Section 6 will be devoted to studying spherically symmetric solutions in mimetic gravity. In Section 7 we will touch upon the issue of rotation curves in mimetic gravity and how this issue is addressed. Finally, we will conclude in Section 8.

1.1. Why Modify General Relativity?

General Relativity (GR henceforth), first formulated by Einstein in 1915 [162–166] (for a pedagogical review, see, e.g., [167]), is an extremely successful and predictive theory, and together with Quantum Field Theory forms one of the pillars of modern physics. The traditional picture of GR is a geometrical one, with the theory being one of space-time and its metric. A more modern view is free of geometrical concepts and sees GR as the unique theory of massless spin-2 particles.

Confirmations of GR abound (see, e.g., [168] for a complete review), ranging from gravitational lensing [169] to the precession of Mercury’s orbit [165]. Shortly after its centennial in 2015 one of the pillars of GR, the existence of gravitational waves, was grandiosely confirmed by the detection of GW150914 [170] and GW151226 [171] by LIGO (see also [172]). Before we even embark on a review of mimetic gravity, then, it is worth reminding the reader why one should even consider questioning a theory as successful as GR. Aside from the philosophical perspective that questioning theories and exploring other approaches are a sensible route in science, provided of course that there is agreement with observations, hints persist in the literature that complicating the gravitational action may indeed have its merits. In fact, the reader should be reminded that as early as 1919 (four years after GR had been formulated), proposals started to be put forward as to how to extend this theory, notably, in the form of Weyl’s scale independent theory [173] and Eddington’s theory of connections [174]. These early attempts to modify GR were driven solely by scientific curiosity with no formal theoretical, or let alone experimental, motivation.

Nonetheless, theoretical motivation for modifying the gravitational action came quite soon. The underlying reason is that GR is nonrenormalizable and thus not quantizable in the way conventional Quantum Field Theories are quantized. However, it was proven that 1-loop renormalization requires the addition of higher order curvature terms to the Einstein-Hilbert action. In fact, it was later demonstrated that, while actions constructed from invariants quadratic in curvature are renormalizable [175], the addition of higher order time derivatives which follows from the addition of terms higher order in the curvature leads to the appearance of ghost degrees of freedom, which entail a loss of unitarity. More recent results show that when quantum corrections or string theory are taken into account, the effective low-energy gravitational action admits higher order curvature invariants (see, e.g., [176] for a general review). However, all these early attempts to modify GR had a common denominator in the fact that terms which modified the gravitational action were only considered to be relevant in proximity of the Planck scale, thus not affecting the late universe.

With the emergence of the “dark universe” picture in recent years, the limits of GR have been fully exposed, and further motivations to modify this theory have emerged. A series of experiments and surveys, including but not limited to CMB experiments, galaxy redshift surveys, cluster surveys, supernovae surveys, lensing experiments, and quasar surveys, have depicted a peculiar picture of our universe [1–50]. This scenario suggests that our näive picture of the world we live in being described by the Standard Model of Particle Physics (SM) supplemented by General Relativity is, at best, incomplete. The concordance cosmology model suggests a scenario where only ~4% of the energy budget of the universe consists of baryonic matter, whereas ~24% consists of nonbaryonic dark matter (DM), and the remaining ~76% consists of dark energy (DE). Of the last two extra components, dark matter is (presumably) the one with properties most similar to ordinary matter. It shares the clustering properties of ordinary matter, but not its couplings to SM gauge bosons (e.g., electromagnetic ones), and is believed to be responsible for structure formation during the matter domination era of the cosmological history. As ordinary matter, DM satisfies the strong energy condition. Dark energy, instead, is more peculiar still, given that it does not share the clustering properties of ordinary matter or DM, as it violates the strong energy condition. It is believed to be responsible for the late-time speed-up of the universe, which has been inferred from a variety of cosmological and astrophysical observations, ranging from type Ia-supernovae to the CMB. Whereas evidence for DE is somewhat indirect and exclusively of cosmological origin, clues as to the existence of DM are instead present on a wide variety of scales, from cosmological to astrophysical (galactic and subgalactic) ones. For technical reviews on DM, see, for instance, [177–181]; for similar reviews on DE, see, for instance, [182–186].

The late-time acceleration of our universe, however, is most likely not the only period of accelerated expansion that our universe has experienced. A period of accelerated (exponential) expansion during the very early universe, prior to the conventional radiation and matter domination epochs, is required to solve the horizon, flatness, and monopole problems. This period of accelerated expansion is known as inflation (see, e.g., [187–196] for pioneering work), and a vast class of models attempting to reproduce such period exists in the literature (for an incomplete list, see, e.g., [197–206] and references therein, see also, e.g., [207–231] for more recent inflationary model-building which is extremely closely relevant to mimetic gravity and variations of it). For reviews on inflation, see, for instance, [232–236]. Inflation also purports to be the mechanism generating primordial inhomogeneities which are quantum in origin [237, 238], providing the seeds which grow under gravitational instability to form the large-scale structure of the universe. The fact that the universe presumably undergoes acceleration at both early and late times or, equivalently, at high and low curvature is very puzzling and might be hinting to a more profound structure.

It thus appears that concordance cosmology requires at least three extra (possibly dark) cosmological components: one or more dark matter components, some form of dark energy, and one or more inflaton fields. There is no shortage of ideas as to what might be the nature of each of these components. Nonetheless, adding these three or more components opens another set of questions, which include but are not limited to the compatibility with the current SM and the consistency of formulation. On the other hand, gravity is the least understood of the four fundamental interactions and the most relevant one on cosmological and astrophysical scales. If so, it could be that our understanding of gravity on these scales is inadequate or incomplete, and modifying our theory of gravitation could indeed be the answer to the dark components of the universe. One could argue that this solution is indeed more economical and possibly the one to pursue in the spirit of Occam’s razor. In other words, modifications to Einstein’s theory of General Relativity might provide a consistent description of early and late-time acceleration and of the dark matter which appears to pervade the universe. Modified theories of gravity not only can provide a solution to the “dark universe riddle” but also possess a number of alluring features such as unification of the various epochs of acceleration and deceleration (matter domination) of the universe’s evolution, transition from nonphantom to phantom phase being transient (and thus without Big Rip), solution to the coincidence problem, and also interesting connections to string theory.

Having presented some motivation to modify our theory of gravitation, we now proceed to briefly discuss systematic ways by means of which this purpose can be achieved.

1.2. How to Modify General Relativity?

Essentially all attempts to modify General Relativity are guided by Lovelock’s theorem [239]. Lovelock’s theorem states that the only possible second-order Euler-Lagrange expression obtainable in 4D space from a scalar Lagrangian density of the form , where is the metric tensor, is given by the following:where and are constants and and are the Ricci tensor and scalar, respectively. It follows that constructing metric theories of gravity whose equations differ from those of GR requires at least one of the following to be satisfied:(i)Presence of other fields apart from or in lieu of the metric tensor(ii)Work in a number of dimensions different from 4(iii)Accept metric derivatives of degree higher than 2 in the field equations(iv)Giving up locality or Lorentz invariance

Therefore, we can imagine broadly classifying the plethora of modified gravity theories according to which of the above assumptions is broken.

1.2.1. New Degrees of Freedom

Relaxing the first assumption leads to what is probably the largest class of modified gravity theories. Theories corresponding to the addition of scalar degrees of freedom include quintessence (e.g., [240–242]) and coupled quintessence (e.g., [243]) theories, the Chern-Simons theory (e.g., [244]), Cuscuton cosmology (e.g., [245–247]), Chaplygin gases (e.g., [248, 249]), torsion theories such as theories (e.g., [250–253], see also, e.g., [254–264] for recent work) or the Einstein-Cartan-Sciama-Kibble theory (e.g., [265–268]), scalar-tensor theories (e.g., [269]) such as the Brans-Dicke theory [270], ghost condensates (e.g., [271]), galileons (e.g., [272, 273]), KGB [274], Horndeski’s theory [275], and many others.

One can instead choose to add vector degrees of freedom, as in the case of the Einstein-aether theory (e.g., [276–280]). The addition of a vector field leads to the introduction of a preferred direction in space-time, which entails breaking Lorentz invariance.

Theories where tensor degrees of freedom are added include Eddington-Born-Infeld gravity (e.g., [281–284]) and bimetric MOND (e.g., [285, 286]) among others. TeVeS (tensor vector scalar gravity, [287]) is instead a theory which features the addition of all three types of degrees of freedom together.

Broadly speaking, mimetic gravity belongs to the class of theories of modified gravity where an additional scalar degree of freedom is added. Caution is needed with this identification though because, as we shall see later, mimetic gravity does not possess a proper scalar degree of freedom, but rather a constrained one.

1.2.2. Extra Dimensions

Relaxing the second assumption instead brings us to consider models with extra dimensions, the prototype of which is constituted by Kaluza-Klein models (e.g., [288–290]). Models of modified gravity with extra dimensions abound when considering string theory, including Randall-Sundrum I [291] and II [292] models, Einstein-Dilaton-Gauss-Bonnet gravity (e.g., [293, 294]), cascading gravity (e.g., [295–298]), and Dvali-Gabadadze-Porrati gravity (e.g., [299, 300]). Another interesting theory that falls within the extra dimension category is represented by 2T gravity [301].

1.2.3. Higher Order

The most famous and studied example of a theory falling within this category is undoubtedly represented by gravity ([302], see also, e.g., [60, 65–67, 303–309] or [310] for black holes phenomenology). In fact, unification of inflation and late-time acceleration was proposed in the context of gravity in [303, 305, 311–313]. Belonging to the same family are also variations of the former such as , , , gravity (see, e.g., [132, 314–323]), but also Gauss-Bonnet (see, e.g., [324–334]) and conformal gravity (see, e.g., [335–337]). Another well-known theory which lies within this family is represented by Hořava-Lifshitz gravity ([338], see also [339–363], which also violates Lorentz invariance explicitly), and correspondingly the vast category of Hořava-like theories, including those which break Lorentz invariance dynamically (e.g., [356, 361, 362, 364–372]).

1.2.4. Nonlocal

If we choose to relax the assumption of locality (we have already seen cases where the assumption of Lorentz invariance is relaxed above), we can consider nonlocal gravity models whose action contains the inverse of differential operators of curvature invariants, such as and gravity (e.g., [373–377]). Some degravitation scenarios belong to this family as well (e.g., [378, 379]).

Broadly speaking, mimetic gravity belongs to the class of theories of modified gravity where an additional scalar degree of freedom is added. Caution is needed with this identification though because mimetic gravity does not possess a proper scalar degree of freedom. Instead, the would-be scalar degree of freedom is constrained by a Lagrange multiplier, which kills all higher derivatives. As such, the mimetic field cannot have oscillating solutions and the sound speed satisfies , confirming that there is no propagation of scalar degrees of freedom (however, this is true only in the original mimetic model but does not necessarily hold in extensions thereof). Furthermore, the same Lagrange multiplier term introduced a preferred foliation of space-time, which breaks Lorentz invariance (although this is preserved at the level of the action). These aspects of mimetic gravity will be discussed in more detail in the following section.

2. Mimetic Gravity

The expression “mimetic dark matter” was first coined in a 2013 paper by Chamseddine and Mukhanov [69] although, as we shall see later, the foundation for mimetic theories had actually been developed a few years earlier in three independent papers [380–382]. In [69], the proposed idea is to isolate the conformal degree of freedom of gravity by introducing a parametrization of the physical metric in terms of an auxiliary metric and a scalar field , dubbed mimetic field, as follows:From (2) it is clear that, in such a way, the physical metric is invariant under conformal transformations of the auxiliary metric of the type , being a function of the space-time coordinates. It is also clear that, as a consistency condition, the mimetic field must satisfy the following constraint:Thus, the gravitational action, taking into account the reparametrization given by (2) now takes the form:where is the space-time manifold, is the Ricci scalar, is the matter Lagrangian, and is the determinant of the physical metric.

By varying the action with respect to the physical metric one obtains the equations for the gravitational field. However, this process must be done with care, for the (variation of the) physical metric can be written in terms of the (variation of the) auxiliary metric and the (variation of the) mimetic field. Taking this dependency into account, variation of the action with respect to the physical metric yields [69]where is the Einstein tensor, with being the Ricci tensor, while and are the trace of Einstein’s tensor and the stress-energy tensor of matter, respectively. Notice that the auxiliary metric does not enter the gravitational field equation explicitly, but only through the physical metric, whereas the mimetic field enters explicitly. In fact, the mimetic field contributes to the right-hand side of Einstein’s equation through an additional stress-energy tensor component:We note that both energy-momentum tensors, and , are covariantly conserved, that is, (with the covariant derivative), whereas the continuity equation for with the mimetic constraint (3) leads toFinally, the trace of (5) is found to beIt is clear that the above is automatically satisfied if one takes into account (3) even if . Thus, the theory admits nontrivial solutions and the conformal degree of freedom becomes dynamical even in the absence of matter ( but ) [69].

Let us examine the structure of the mimetic stress-energy tensor. Recall that the stress-energy tensor of a perfect fluid whose energy density is and pressure is given byNotice that the mimetic stress-energy tensor in (6) assumes the same form of the one of a perfect fluid with pressure and energy density , while the gradient of the mimetic field, , plays the role of 4-velocity. Thus, the mimetic fluid is pressureless, suggesting it could play the role of dust in a cosmological setting. To confirm whether the mimetic fluid can indeed play the role of dust, it is necessary to investigate cosmological solutions. In fact, this is easy to do on a Friedmann-Lemaître-Robertson-Walker (FLRW) setting, with a metric of the form:where is the scale factor. If we take the hypersurfaces of constant time to be equal to those of constant , we immediately see that the constraint equation (3) is automatically satisfied if the mimetic field is identified with time up to an integration constant (which we arbitrarily set to 0). Thus, the mimetic field plays the role of “clock” on an FLRW background. It is then easy to show that (7) implies that , which corresponds to the energy density of the mimetic stress-energy tensor, decays with the scale factor of the FLRW universe as . Recall that the energy density of a component with equation of constant state parameter evolves as in an FLRW universe, and hence the evolution in the energy density of the mimetic field corresponds to , namely, the equation of state for dust. In other words, the conformal degree of freedom of gravity can mimic the behaviour of dark matter at a cosmological level, hence the name “mimetic dark matter” [69].

Lagrange Multiplier Formulation. Before further discussing some fundamental aspects of mimetic dark matter (or mimetic gravity henceforth), such as the reason behind the different solutions from GR despite the seemingly innocuous parametrization given by (2), let us discuss an alternative but equivalent formulation of mimetic gravity. The equations of motion obtained from the action written in terms of the auxiliary metric are equivalent to those one would conventionally obtain from the action expressed in terms of the physical metric with the imposition of an additional constraint on the mimetic field. This suggests that the mimetic constraint given by (3) can actually be implemented at the level of the action by using a Lagrange multiplier. That is, the action for mimetic gravity (4) can be written asVariation of the Lagrangian with respect to the Lagrange multiplier field leads to (3), while variation with respect to the physical metric yieldswhose trace, when one takes into account (3), is given byThus, one recovers (5) again. In this review, we will always make use of the action given by (11) to introduce the mimetic field.

A remark is in order here. Actions such as (11) had actually been introduced three years before the term “mimetic dark matter” was first coined. Three independent papers in 2010, by Lim et al. [380], Gao et al. [381], and Capozziello et al. [382], respectively, presented models with two additional scalar fields, one of which playing the role of Lagrange multiplier enforcing a constraint on the derivative of the other (see also [383] for recent work on the role of Lagrange multiplier constrained terms in cosmology). In fact, it was shown that these types of models can produce a unified theory describing dark energy and dark matter, because the term inside the Lagrange multiplier can always be arranged in such a way to reproduce the conventional expansion history of CDM. Thus, it is fair to state that mimetic theories had actually seen birth prior to the 2013 paper by Chamseddine and Mukhanov.

2.1. Understanding Mimetic Gravity

Before we can make further progress in exploring solutions in mimetic gravity, generalizing the theory, or studying connections to other theories, we need to touch on two very important points: first, why the seemingly innocuous parametrization given by (5) has led to a completely new class of solutions not contemplated by GR, and second, whether the theory is stable or not. As we shall see, the first point can ultimately be explained in terms of singular disformal transformations.

It might appear puzzling at first that, only by rearranging parts of the metric, one is faced with a different model altogether. A first explanation appeared in [70], which explained this property in terms of variation of the action taking place over a restricted class of functions. This is the case in mimetic gravity, precisely because the consistency equation (3) enforces an additional condition for any admissible variation of the action, in particular demanding thatThus, for a spatially homogeneous mimetic field , . Varying over a restricted class of functions now provides less conditions for the stationarity of the action and hence more freedom in the dynamics. This is a well-known property relevant when one makes derivative substitutions into an action : the class of functions over which the variation is admissible does not only comprise those for which the variation of is vanishing at the boundary but also those where the integral of the variation of is zero. This extra restriction leads to less conditions for stationarity of the action and the appearance of additional dynamics with respect to the original case.

Another explanation was presented in [71], which identified mimetic gravity as a conformal (Weyl-symmetric) extension of GR. The first important point to notice is that the parametrization is noninvertible even for fixed , owing to the fact that the map is a map from ten variables to eleven. With this parametrization, the theory is manifestly conformally invariant, that is, invariant with respect to transformations of the auxiliary metric (and correspondingly the action ) of the form:where is a function of the space-time coordinates. Two immediate corollaries of the theory’s conformal invariance are its yielding of identically traceless equations of motion for the gravitational field and requiring conformal gauge fixing. In fact, recall that the equation of motion for the gravitational field (5) is automatically traceless if the consistency condition given by (3) is satisfied. Therefore, we can identify (3) with the conformal gauge condition in the locally gauge-invariant theory with action . Mimetic gravity can thus be seen as a conformal extension of GR which is Weyl-invariant in terms of the auxiliary metric . However, this theory is quite different from the off-shell conformal extensions of GR proposed in [384] (which preserves GR on-shell but modifies its effective action off-shell); here the gravitational action is already modified at a classical level, by adding an extra degree of freedom provided by a collisionless perfect fluid. This additional degree of freedom, which can mimic collisionless cold dark matter for cosmological purposes, arises from gauging out local Weyl invariance.

2.1.1. Singular Disformal Transformations

As we anticipated above, mimetic gravity and the appearance of the extra degree of freedom which can mimic cosmological dark matter are rooted into the role played by singular disformal transformations. As was shown by Bekenstein [385], because GR enjoys invariance under diffeomorphisms, one is free to parametrize a metric in terms of a fiducial metric and a scalar field . The map between the two is defined as a “disformal transformation,” or a “disformation,” and takes the following form:where . The functions and are referred to as conformal factor and disformal factor, respectively. In general the functions , are arbitrary, with . It is easy to show that, provided the transformation is invertible, the equations of motion for the theory (obtained by variation of the action with respect to and ) reduce to those of obtained by varying with respect to the metric [79].

To make progress, it is useful to contract the two equations of motion with and . Although we will not perform the steps explicitly, it is easy to show that this leads to the following two equations of motion:where the two quantities and are defined asThe determinant of system (17) is given byIf the above is nonzero, it is trivial to obtain that the resulting set of equations consists of Einstein’s equation () and a second empty equation. Therefore, the theory does not feature new solutions with respect to GR [79].

The situation is quite different if the determinant in (19) is zero, which corresponds to the physical case when the disformal transformation given by (16) is noninvertible or singular. In this case, being the function , this immediately determines , which is of the form:where is an arbitrary function. We will not show the steps explicitly, which are instead discussed in detail in Section of [79], but it is not hard to obtain the equations of motion and notice that they differ from those of GR, due to the presence of an extra term on the right-hand side of Einstein’s equation (i.e., an additional contribution to the stress-energy tensor). Therefore, when the disformal transformation is singular, one is faced with the appearance of extra degrees of freedom which result in equations of motion differing from those of GR [79].

The parametrization (5) defining mimetic gravity can be identified with a singular disformal transformation, with and in (16), and correspondingly in (20). In general, when the relation defined by (20) exists between the conformal factor and the disformal factor , the resulting disformal transformation is singular and, as a result, the system possesses additional degrees of freedom, explaining the origin of the extra degree of freedom in mimetic gravity which mimics a dust component. This aspect has been at the center of a number of studies recently; see, for instance, [79, 90, 96, 102, 104, 107, 128, 140, 142, 386–389]. Moreover, [104] has shown that the two approaches towards mimetic gravity (and further extensions to be discussed later, such as mimetic Horndeski theories), namely, Lagrange multiplier (11) and singular disformal transformation (16), are in fact equivalent.

2.1.2. Mimetic Gravity from the Brans-Dicke Theory

There actually exists a third route to mimetic gravity, apart from disformal transformations and Lagrange multiplier, whose starting point is the singular Brans-Dicke theory. Namely, by starting from the action (4) (neglecting matter terms) and by performing the conformal transformation given by (5), we end up with the action [129]:where we have defined . One immediately sees that (21) corresponds to the singular/conformal Brans-Dicke action [270] with density parameter . Thus, we conclude that a third way of obtaining mimetic gravity is by substituting the kinetic term in lieu of the scalar field in the singular Brans-Dicke action [129]. In case matter fields are included, the substitution has to be performed on the matter part of the Lagrangian as well, which means that matter will not be coupled to but to [129].

2.2. Stability

Is mimetic gravity stable? In other words, does its spectrum contemplate the presence of states with negative norm, or fields whose kinetic term has the wrong sign (corresponding to negative energy states), which could possibly destabilize the theory? This is an important question which has yet to find a definitive answer. Recall that the original mimetic theories formulated in 2010 were found to suffer from a tachyonic instability [390]; therefore, the question of whether mimetic gravity is stable is totally pertinent.

If we formulate the theory of mimetic gravity using the physical metric , we inevitably incur into the risk of the appearance of higher derivatives of the mimetic field, which could entail the emergence of ghosts. Addressing this question requires performing a Hamiltonian analysis of the theory, identifying all constraints and counting the local degrees of freedom. A preliminary analysis of this problem was conducted in [71], which concluded that the theory is stable if the energy density of the mimetic field is positive. This condition is, of course, easy to understand physically. Moreover, it indicates a preference for de Sitter-type backgrounds with a positive cosmological constant, since in that case both contributions to the energy density, given by curvature and trace of the matter stress-energy tensor, are positive. Therefore, it is presumed that mimetic is gravity is stable provided the time evolution of the system preserves the positivity of the energy density stored in the mimetic degree of freedom. The work of [71] also identified another possible instability issue, namely, caustic singularities (which are not dangerous at the quantum level, unlike ghost instabilities). These are presumably due to the pressureless nature of the mimetic field and can possibly be circumvented if one modifies the theory with higher derivative terms, as we will discuss later.

The analysis of [71] imposed the conformal gauge condition (3) prior to proceeding to the canonical formulation. A proper analysis should instead take place in full generality and has been conducted in [75]. This analysis finds that the Hamiltonian constraint depends linearly on the momentum, which in most cases signals that the Hamiltonian density of the theory is unbounded from below, a classical sign of instability. As anticipated, this occurs frequently for higher derivative theories, which are prone to the Ostrogradski instability. The work in [75] concluded, as [71], that mimetic gravity is stable as long as the energy density in the dust degree of freedom in the theory remains positive. However, this is not always consistent with the dynamics of the theory, given that for some initial configurations, the energy density could start its evolution with a positive value but then end up with a negative value, which would cause instability. In fact, [75] finds that the requirement that the theory be stable correspond to the requirement that initial configurations do not cross the surface for which the momentum conjugate to the mimetic field, , satisfies .

A possible solution to these instability issues was presented in [71] and studied in [75]. The idea is to modify the parametrization (5) by making use, instead of the gradient of a scalar field, of a dynamical vector (Proca) field:The Proca field is made dynamical by adding a Maxwell kinetic term to the action, where is the field-strength tensor of the vector field. It is beyond the scope of our review to provide details of the analysis, conducted in [75] which finds that the Hamiltonian in the Proca mimetic model shows no sign of instability. Furthermore, [75] proposes an interesting extension of mimetic gravity to a mimetic tensor vector scalar model, by generalizing (2) towhere now both the scalar and vector degrees of freedom contribute to mimetic matter. It is furthermore demonstrated that the theory is free of ghosts [75].

Another recent work confirmed in all generality that the original mimetic gravity theory suffers from ghost instability [140], in the following way. It is immediate to show that mimetic gravity is invariant under the local symmetry defined bywhere as usual and correspond to the conformal and disformal factors. Being and for mimetic gravity, (24) corresponds as expected to invariance of the physical metric under conformal transformations of the auxiliary metric. In the Hamiltonian description, this symmetry is associated with a first class constraint. In fact, one can show that the primary constraint, which corresponds to the generator of infinitesimal conformal transformations, is first class, with its Poisson commuting with the Hamiltonian and momentum constraints. This leaves no place for a secondary constraint which could eliminate the Ostrogradski ghost. Thus, this result confirms indeed that the original mimetic gravity proposal suffers from a ghost instability.

3. Solutions and Extensions of Mimetic Gravity

Having discussed the underlying physical foundation of mimetic gravity, and its stability, we can now proceed to study solutions and extensions of this theory.

3.1. Potential for Mimetic Gravity

Recall that, in a cosmological setting, the mimetic field plays the role of “clock.” Therefore, one can imagine making the mimetic field dynamical by adding a potential for such field to the action. A field-dependent potential corresponds to a time-dependent potential which, by virtue of the Friedmann equation, corresponds to a time-varying Hubble parameter (and correspondingly scale factor). Therefore, by adding an appropriate potential for the mimetic field, one can in principle reconstruct any desired expansion history of the universe. This is the idea behind the minimal extension of mimetic gravity first proposed in [73]. The action of mimetic gravity (in the Lagrange multiplier formulation) is thus extended to include a potential for the mimetic field:The equations of motion of the theory are then given bywhich, by taking the trace, can be used to determine the Lagrange multiplier:When plugged into (26), (27) yieldsOf course, variation with respect to the Lagrange multiplier as usual yields the constraint (3). Thus, when a potential is added to the action, the mimetic field contributes a pressure and energy density of and , respectively [73]. One further equation of motion can be derived by varying the action with respect to the mimetic field, which giveswhen taking into account the expression for the Lagrange multiplier (27).

To study cosmological solutions, it is useful to consider a flat FLRW background (10), since therein the mimetic field can be identified with time. In this case, it is not hard to show that (29) reduces to [73]which can be integrated to givewhereas the pressure remains . The Friedmann equation can instead be manipulated to the form:where as usual the Hubble parameter is defined as . Further progress can be made by performing the substitution , which yields the following equation [73]:It should be noticed that the equations of motion simplify greatly because of the identification of the mimetic field with time on an FLRW background. In this way, the pressure becomes a known function of time and satisfies a linear differential equation. We now proceed to study a few interesting potentials and the corresponding solutions.

3.1.1.

Let us consider the following potential [73]:Solving the corresponding (33) and substituting for the scale factor, , yield the following solution for :where is an integration constant. For the solution describes an oscillating flat universe with amplitude of oscillations which grows with time; however, the solution presents singularities and for this reason we will not write it down explicitly [73]. We can furthermore determine the equation of state parameter (EoS) if we recall that the energy density is given by (31) whereas the pressure reads . Explicit calculation gives the following:It is interesting to note that, for , the EoS approaches , that is, a cosmological constant, at late times.

We can consider the case where mimetic matter is a subdominant energy component in the universe, which is instead dominated by another form of matter with EoS  . The scale factor in this situation evolves asand hence (31) can be used to deduce the evolution of the energy density of mimetic matter, which decays asGiven that the pressure of mimetic matter reads , we immediately see that the EoS for mimetic matter is , demonstrating that mimetic matter, when subdominant, can imitate the EoS of the dominant energy component [73]. A comment is in order here. Mimetic matter can only be subdominant if . If this condition is not satisfied, mimetic matter will only start imitating the dominant matter component at late times, while acting as a cosmological constant at earlier times.

3.1.2. Power-Law Potential

We can consider an arbitrary power-law potential:for which the solution of (33) can be written in terms of the Bessel functions [73]:For (with corresponding to the case we have studied previously) the limiting behaviour of the scale factor is that of a dust-dominated universe, with EoS . For and (which corresponds to a positive pressure), the corresponding solution is a singular oscillating universe. For and instead, the pressure is negative and hence we expect accelerating solutions [73]. In fact, corresponds to a cosmological constant as expected (the potential is simply a constant), whereas gives an inflationary expansion solution with scale factor:which resembles that of chaotic inflation sourced by a quadratic potential [73].

3.1.3. Inflation in Mimetic Gravity

One can always reconstruct the appropriate potential for the mimetic field which can provide an inflationary solution. The method is very simple: choose a desired expansion history of the universe [encoded in the Hubble parameter or equivalently in the scale factor ], find the corresponding -parameter through , and then invert (33) to find the corresponding potential which can provide the desired expansion:As an example, we can consider the following potential [73]:whose corresponding solution for the scale factor is exponential () for large negative times and corresponding to EoS (i.e., dust, ) at late times. Thus, we see that with this potential mimetic gravity can provide us with an inflationary solution with graceful exit to the matter dominated era [73].

Another interesting possibility is given by an exponential potential [83]:In this case, the scale factor can be expressed in terms of Bessel functions of the first and second kinds [83]:where and are integration constants and and are the modified Bessel functions of order zero, of the first and second kinds, respectively. At late times, the behaviour of the scale factor is that of a matter dominated universe, that is, with EoS (). At early times the behaviour of the scale factor depends on the sign of , in particular providing us with an inflationary solution for , whereas the solution is a bouncing nonsingular universe for . A similar behaviour can be obtained if one chooses the potential:that is, a solution with inflation at early times and matter domination at late times. More complicated potentials which can reproduce qualitatively similar behaviours were studied in [103].

3.1.4. Bouncing Universes in Mimetic Gravity

As we have already seen in previous cases, one can easily construct bouncing solutions in mimetic gravity. Let us work through one further example here. Consider a potential of the form:As usual, the scale factor can be determined by solving (33), which yields [73]If we set the integration constant to 0, the corresponding energy density and pressure (one again, refer to (31) and recall that ) are given byAt very early times (large negative ) the EoS approaches ; the universe is dominated by dust and contracts. At a certain time corresponding to , the energy density drops suddenly to zero, after which the universe begins expanding. During the first instants of the expansion (within one Planckian time), the energy density of the universe is Planckian but subsequently drops as the expansion proceeds as a conventional expansion in a dust-dominated universe. The interesting feature of this potential is that the EoS crosses the phantom divide without singularity. This remains true even in the general case where the integration constant is nonzero, provided [73].

In the case we have just examined, the bounce occurs at the Planck scale, and hence the classical analysis we provided might not be valid as quantum gravity effects would be playing an important role. However, a minimal modification allows lowering the scale of the bounce and correspondingly increases the duration of the bounce (which now lasts more than a Planckian time). The corresponding potential which can provide this behaviour is given by [73]Although we will not show the solution explicitly, in this case the scale of the bounce is reduced to and the duration of the bounce is now [73].

3.2. Mimetic Gravity

The next step which was performed by Nojiri and Odintsov is to generalize mimetic gravity to mimetic gravity [81]. In this theory we expect to have two additional degrees of freedom compared to GR: the constrained (nonpropagating) scalar degree of freedom of mimetic gravity and the additional scalar degree of freedom arising from the term. The action of the theory is given by [81]where as usual the relation between the physical and auxiliary metric and the mimetic field is given by (2), and the constraint equation (3) has to be satisfied. Because of this, the action of mimetic gravity can be equally written employing a Lagrange multiplier field analogously to (11) [81]:where the Lagrange multiplier enforces the constraint on the gradient of the mimetic field and in addition we have added a potential for the mimetic field.

The equations of motion of the theory are slightly more complicated than that of conventional mimetic gravity. Varying with respect to the metric gives the gravitational field equations [81]:Variation with respect to the mimetic field instead yields the following equation [81]:As usual, by construction, variation with respect to the Lagrange multiplier gives the mimetic constraint:

As we have mentioned previously, in mimetic gravity one has two additional degrees of freedom. Therefore, by appropriately tuning either or both the potential for the mimetic field, or the form of the function, one can reconstruct basically any desired expansion history of the universe. On the other hand, the interpretation of the cosmological role played by mimetic dark matter remains the same as in mimetic gravity [81].

Let us proceed to study some of the properties of mimetic gravity in a cosmological setting. As usual, we consider a flat FLRW universe, and we model the matter contribution as that of a perfect fluid with energy density and pressure . Assuming that the mimetic field depends only on time, (53), (54), and (55) can be expressed as follows [81]:The last equation shows that up to an integration constant, which we can set to zero, the mimetic field can be identified with time just as in ordinary mimetic gravity. Thus, (58) can be expressed as follows:which, if we assume that the contribution of ordinary matter is negligible (), reduces toOn the other hand, (57) can be solved for as follows:which shows that (59) is automatically satisfied.

The above equations put on a quantitative footing the statement we previously made: namely, that by tuning the behaviour of either or both the two additional scalar degrees of freedom, we can reconstruct any possible expansion history of the universe [81]. For instance, one can imagine fixing the form of the scalar potential and then reconstruct the form of which gives the wanted evolution [encoded in or, equivalently, ]. Alternatively, one can start from a given form of which might not admit the wanted evolution (e.g., matter dominated-like expansion followed by accelerated expansion) and reconstruct the form of the scalar potential which can allow for such expansion. It should also be remarked that any solution in conventional gravity is also a solution in mimetic gravity, but the converse is not true. In [81] specific solutions which allow unification of early-time inflation and late-time acceleration with intermediate matter domination era, as well as bouncing Universes, are studied and it is shown that they can be implemented in mimetic gravity. Of course, the exact forms of the mimetic potential or the function in these cases are quite complicated; nonetheless, the study serves as a proof of principle that, in such theories, one can realize any given expansion history of the universe without the need for dark components, which remains the main goal of modified theories of gravity.

Three further recent studies by Odintsov and Oikonomou [112, 120, 130] have demonstrated how one can, in mimetic gravity, realize inflationary cosmologies which are compatible with Planck and BICEP2/Keck Array constraints on the scalar spectral index and on the tensor-to-scalar ratio, and , respectively. Both the reconstruction (determination of the potential once the form of and the evolutionary history of the universe are given) and the inverse reconstruction (determination of the form of once the potential and the evolutionary history of the universe are given) are studied in detail and it is demonstrated that several viable options for realizing inflation compatibly with observational constraints are possible. However, the studies also point out a possible weakness of mimetic gravity in this respect: namely, the forms of both the mimetic potential and the function can become extremely complicated. The forms of both the mimetic potential and the Lagrange multiplier increase in complexity as the complexity of the form increases. Therefore, these and other studies on mimetic gravity and extensions thereof should not be viewed as the ultimate cosmological theory of everything, but rather as a proof of principle that within these theories one can reproduce basically any cosmological scenario and thus solve the “dark universe” problems, although this might come at the cost of sacrificing simplicity.

In addition, [112, 120, 130] also remarked that, although in principle the forms of the mimetic potential and the function can be arbitrary, it must be kept in mind that, in order to realize viable inflation, a mechanism for graceful exit to the conventional radiation dominated era must be achieved. This entails ensuring that the theory contains unstable de Sitter vacua, which eventually becomes the cosmological attractor of the dynamical system. It is precisely the functional form of which has to ensure that graceful exit takes place. Therefore, in mimetic gravity, although in principle the form of the potential is arbitrary, the same cannot be said about the functional form of , which has to be such as to ensure graceful exit from inflation. Therefore, in the interest of simplicity, a practical approach to constructing a minimal model of mimetic gravity with the desired inflationary properties would be to choose the simplest possible functional form of which ensures graceful exit from inflation, then performing the reconstruction technique to determine the form of the mimetic potential which allows the desired expansion history following inflation to be realized. Another possible solution, which we will discuss shortly, is to consider inflation [95, 116], where a dynamical scalar field is coupled to gravity.

3.2.1. Late-Time Evolution in Mimetic Gravity

So far we have discussed mimetic gravity and variants thereof at early times, that is, at the epoch when primordial curvature perturbations were generated. However, it is also interesting to consider late-time evolution in mimetic gravity. The equations of motion are incredibly complex and in principle do not allow for analytical solutions. However, this complexity can be bypassed by means of the method of dynamical analysis (see, e.g., [391–394]), which gives information about the global behaviour of solutions. In particular, one proceeds by transforming the equations of motion into their autonomous form and extract the critical points. Subsequently perturbations are linearized around these critical points and expressed in terms of the perturbation matrix, the eigenvalues of which determine the type and stability of the critical points.

A detailed dynamical analysis of mimetic gravity was presented in [86]. This type of analysis allows us to bypass the complexity of the equations of motion by extracting critical points and studying the corresponding observables, such as the energy densities of the various energy components, the corresponding EoS, and the deceleration parameter. In particular, the analysis finds that the only stable critical points, that is, those that can play the role of attractors at late times, are those that exist in gravity as well. In other words, stable solutions in mimetic gravity can only affect the expansion history of the universe at early and intermediate times, whereas at late times the expansion history has to coincide with that driven by conventional gravity. An immediate implication of this finding is that, although mimetic gravity could drive inflation differently from gravity, the late-time acceleration of the universe in these theories has to coincide with the usual gravity one [86]. However, these conclusions have been reached only by studying the theory at the level of the background. It is expected that different conclusions would be reached if the same study would be performed at the level of perturbations. This is true because the new terms present in the equations of motion of mimetic gravity compared to conventional gravity can contribute to the perturbation equations, although they do not contribute at the background level [86]. Finally, the energy conditions required to avoid the Dolgov-Kawasaki instability in mimetic gravity were studied in [108], which found that these are the same as in conventional gravity.

As we mentioned above, the conclusions reached about the late-time evolution in mimetic gravity hold only at the background level. However, in conventional gravity, there exists a serious problem during the late-time evolution at the perturbation level, namely, that of dark energy oscillations [395] (see also, e.g., [396–399]). The degree of freedom associated with the modification of GR (that is, ) leads to high frequency oscillations of the dark energy around the line of the phantom divide during matter era. As a consequence, some derivatives of the Hubble parameter may diverge and become singular, and the solution is unphysical. Usually in conventional gravity the problem is solved by adding power-law modifications by hand.

In [153], it was instead argued that in mimetic gravity it is possible to overcome the problem by a suitable choice of the potential. By appropriately choosing the potential and the Lagrange multiplier, it is possible to damp the oscillations within a mimetic model whose corresponding conventional model suffered from the oscillations problems. The oscillations die out for redshifts , so there is no issue with dark energy oscillations at our current epoch. Moreover, the values of the dark energy equation of state and the total equation of state are very close to the observed values. The model can in principle be discerned from gravity in that the predicted growth factor is lower in magnitude, a very testable prediction in view of future experiments, further supporting the viability of mimetic gravity as a cosmological framework.

3.2.2. Mimetic Gravity

The question of constructing a theoretically motivated but at the same time simple model for mimetic gravity which provides an early-time inflation epoch, but at the same time graceful exit from the latter, was addressed in [95]. Here, a model of mimetic gravity was considered, where is a scalar field coupled to gravity. We will not go into the details of the work, for which we refer the reader to the original paper [95]. The basic idea is to use the sector to reproduce a variety of cosmological scenarios: among the ones considered in the paper were accelerated cosmologies at high and low curvatures (thus unifying inflation and late-time acceleration), with Einstein gravity at intermediate curvatures. In particular, the accelerated cosmologies are realized by making use of a “switching-on” cosmological constant. The dynamical field evolves in such a way to allow for graceful exit from the inflationary period, thus making the vacua of the first de Sitter period (corresponding to inflation) unstable. Entry into the late-time accelerated epoch, represented by a stable de Sitter attractor, is also made possible by the dynamical field, which thus links all epochs of the expansion history of the universe in a unified way. In the minimal case studied in [95], the mimetic component ensures the presence of cosmological nonbaryonic dark matter although, as we have extensively discussed, it is possible by adding a suitable potential for the mimetic field to obtain similar solutions, but with a different form of .

3.2.3. Nonlocal Mimetic Gravity

A further extension of mimetic gravity was presented in [133], which embeds the theory into the framework of nonlocal theories of gravity. Recall that these theories were first presented in [373], inspired by quantum loop corrections. These theories feature nonlocal operators (i.e., inverse of differential operators) of the curvature invariants. The prototype of nonlocal mimetic gravity is given by the action [133]:where as usual is the mimetic field and the Lagrange multiplier term enforces the constraint on its gradient. It is actually more useful to introduce an additional scalar field , which allows us to translate the action given by (63) to a local scalar-tensor form, as follows [133]:where is an additional Lagrange multiplier which enforces the constraint on the scalar field :Aside from the two constraint equations obtained by varying the action with respect to the Lagrange multipliers, variation of the action with respect to the yields the gravitational field equations [133]:Instead, variation with respect to the two scalar fields leads to the following equations of motion [133]:In [133], this model was studied in detail making use of the reconstruction technique. In particular, two forms for the function have been studied: exponential and power-law. It was shown that appropriate choices for the mimetic potential, as expected, can give the desired expansion history of the universe, which in the cases studied included viable inflation unified with late-time acceleration with intermediate epoch compatible with Einstein’s gravity and cosmological dark matter provided by the mimetic field, as well as solutions with cosmological bounces.

3.3. Unimodular Mimetic Gravity

As we have seen, mimetic gravity provides a geometric explanation for dark matter in the universe, with dark matter emerging as an integration constant as a result of gauging local Weyl invariance, without the need for additional fluids. An older theory, known as unimodular gravity [400] (see also [98, 401–419]), had instead been proposed much earlier to solve, in a geometrical fashion as well, one more of the conundrums of modern cosmology: the dark energy problem. In this framework, dark energy emerges in the form of a cosmological constant from the trace-free part of Einstein’s field equations, with the trace-free part which results in turn by enforcing the condition that the square root of (minus) the determinant of the metric is equal to 1, or in general a constant. It would therefore be interesting to combine the two different approaches of mimetic gravity and unimodular gravity into a single framework which could geometrically explain both dark matter and dark energy by a vacuum theory, without need for additional fluids. This is the proposal of Nojiri et al. in [136].

In order to combine mimetic gravity and unimodular gravity it is necessary to enforce two constraints. The first is the constraint on the gradient of the mimetic field (3), whereas the second is the unimodular constraint:In order to enforce these two constraints, it is conceptually simple to make use of two Lagrange multipliers. This approach has two advantages. First, it keeps the two concepts of mimetic and unimodular gravity separate and facilitates the extraction of physical information. Second, it is as usual more convenient to have the two constraints emerge from the equations of motion. The action for unimodular mimetic gravity thus reads [136]where variation with respect to the Lagrange multiplier enforces the mimetic constraint (3), whereas variation with respect to the Lagrange multiplier enforces the unimodular constraint (68).

The equations of motion for the gravitational field are obtained by varying the action with respect to the metric and are given by [136]:whereas variation with respect to the mimetic field yields the usual equation of motion:Although we do not show the steps explicitly, for which we refer the reader to the original paper [136], in the usual FLRW setting it is possible to manipulate the Einstein equations in order to get the following reconstruction equation for the mimetic potential (as usual, the mimetic field can be identified with time):The content of the above equation is clear: as in all extensions so far discussed of mimetic gravity, one can always reconstruct the potential for the mimetic field which can provide the desired expansion encoded in or . The reconstruction technique is very powerful although, as we have seen, the corresponding potentials are complicated and somewhat hard to justify from first principles, although the reconstruction technique serves in this case as a proof of principle tool.

Having made this consideration, let us consider a few examples where the reconstruction technique is applied. Let us consider the following simple potential [136]:It can be easily shown that it leads to the following solution for the scale factor and the Hubble parameter :that is, a de Sitter cosmology. The functional forms of the Lagrange multiplier are also quite simple and are given by the following:Another choice for the potential which leads to a physically interesting solution is the following [136]:for which the scale factor and the Hubble parameter readThe above solution is that corresponding to an universe dominated by a fluid with EoS . The two Lagrange multipliers are given byThus, we see that with two relatively simple choices of potential it is possible to reconstruct two important expansion histories of the universe: the late-time de Sitter phase and the expansion dominated by a perfect fluid with arbitrary EoS. Although we will not discuss this case explicitly here, for which instead we redirect the reader to [136], it is possible by a choice of a more complicated potential, to realize a viable inflationary model within unimodular mimetic gravity, which is compatible with bounds from Planck and BICEP2/Keck Array. Moreover, it has been shown that graceful exit from these types of inflationary periods can be achieved by ensuring that the corresponding de Sitter vacua which drives the period of accelerated expansion is unstable.

Two further comments are in order here. First, it is possible to provide an effective fluid description of unimodular mimetic gravity [136]. Namely, manipulation of the Einstein equations shows that the contribution of the unimodular and mimetic parts of the action can be interpreted as that of a perfect fluid carrying energy density and pressure as follows:where is defined asFurthermore, the effective energy density and pressure defined as per above satisfy the continuity equation.

The second comment is related to the fact that the unimodular constraint is enforced in a noncovariant way (cf. the action given by (69), where one of the terms in is not multiplied by ). It is nonetheless possible to present a covariant formulation of unimodular mimetic gravity via the following action:where is a three-form field, variation of which gives the constraint , implying that the Lagrange multiplier is constant. On the other hand, the covariant version of the unimodular constraint is obtained by variation with respect to , from which one is left withFurther manipulation, for which we refer the reader to [136], shows that the Friedmann equations one obtains from the covariant version of unimodular mimetic gravity are equivalent to those of the noncovariant version, and thus one can reproduce precisely the same cosmological scenarios in both theories.

3.3.1. Unimodular Mimetic Gravity

A minimal extension of the unimodular mimetic gravity framework we have discussed so far is to consider unimodular mimetic gravity, which is described by the action [138]:As expected, the equations of motion are slightly more complicated than in the unimodular mimetic case, but no conceptual difficulty is added. Specifically, variation with respect to the two Lagrange multipliers enforces the usual unimodular and mimetic constraints, whereas variation with respect to the metric gives rise to the equations for the gravitational field [138]:Finally, variation with respect to the mimetic field yields the following equation:The more complex structure of the equations of motion complicates the reconstruction procedure, that is, the equivalent of (72), which now reads [138]where the function is given byDespite the increased complexity of the equations of motion, the same considerations apply as for unimodular mimetic gravity, as well as all extensions of mimetic gravity hereto considered. Namely, it is always possible to reconstruct any viable cosmological expansion scenario, including unification of inflation and late-time acceleration with intermediate radiation and matter domination, with graceful exit from inflation triggered by unstable de Sitter vacua. This can be achieved by appropriately choosing either or both the form of the mimetic potential or the function . The price to pay would eventually be a considerable complexity in the functional form of both, which of course does not represent a first principle obstacle [138].

3.4. Mimetic Horndeski Gravity

One can further consider more general scalar-tensor theories, which can be “mimetized” according to the procedures we have described so far, namely, through a singular disformal transformation or through a Lagrange multiplier term in the action enforcing the mimetic constraint. In fact, analogously to GR, one can show that the most general scalar-tensor model is invariant under disformal transformations, provided the latter is invertible. This has been shown in all generality in [104]. One can then “mimetize” such theories by considering the following action [104]:where are integers and is the Lagrangian density which is a function of the metric and the mimetic field. In general, the constraint enforced by the Lagrange multiplier can be generalized to , but for the sake of simplicity we will set here and redirect the reader to the work of [104] for more general discussions, and for the explicit form of the equations of motion. In fact, setting is basically equivalent to assigning a potential to the mimetic field. In [104] it was shown that the two approaches to mimetic Horndeski gravity, namely, singular disformal transformation and Lagrange multiplier, are equivalent.

Of course, the considerations made above can be applied in the case of a specific scalar-tensor model, namely, Horndeski gravity [275] (see, e.g., [420–438] for further recent work on the topic of Horndeski gravity and theories beyond Horndeski). Recall that Horndeski gravity is the most general 4D local scalar-tensor theory with equations of motion no higher than second order. The Horndeski action can be written as a sum of four terms:where s readwhere , , , the functions , , , are free, and denotes differentiation with respect to .

The mimetic version of the above Horndeski model has been studied in a variety of papers recently (e.g., [104, 114, 128]). We report some particular cases taken into consideration. We remark that the freedom in the free functions , , , and , as well as in the function (which, when , is equivalent to providing a potential for the mimetic field), results in the possibility of reproducing basically any given expansion scenario of the universe. A specific model studied in [104] is one where the functions take the form:In this case, on a flat FLRW background the solution is given bywith and integration constant and . Thus, we see that the scenario under consideration has reproduced the expansion history of a universe filled with a perfect fluid with EoS [104]. Another case considered in [104] is the mimetic cubic Galileon model, where the functions take the following form:where the cut-off scale is subsequently set to . It is then found that the model can reproduce the expansion history of a universe filled with nonrelativistic matter, followed by a cosmological constant dominated expansion analogous to the late-time acceleration we are experiencing [104]. The case of a nonminimal coupling to the auxiliary metric, was examined as well, which we will not report on here and for which we will redirect the reader to the original paper [104].

To conclude, we report on the following specific case of mimetic Horndeski model which was studied in [114]. The Horndeski part of the action of the theory is given bywhich corresponds to the following choice for the functions discussed previously:A number of solutions, including cosmological bounces, inflation unified with late-time acceleration, and future singularities, have been discussed. As a specific example, on a flat FLRW background, the following choice of potential for the mimetic field [114],gives rise to the following solution:which represents a regular bounce solution. Another bounce solution can be obtained by considering the following potential [114]:for which the scale factor and correspondingly the Hubble parameter read:In Section 5 we will present a detailed discussion on a specific mimetic Horndeski model, constructed as an extension of a covariant Hořava-like theory of gravity.

3.5. Einstein-Aether Theories, Hořava-Lifshitz Gravity, and Covariant Renormalizable Gravity

In closing, we comment on some connections between mimetic gravity and other theories of modified gravity, such connections having been identified recently: namely, the scalar Einstein-aether theory, Hořava-Lifshitz gravity, and a covariant realization of the latter, that is, covariant renormalizable gravity.

3.5.1. Einstein-Aether Theories

An interesting connection which can be identified is that between mimetic gravity and Einstein-aether theories [276, 277] (see also [278–280, 439–444]). These are a class of Lorentz-violating generally covariant extensions of GR. By Lorentz-violating generally covariant we mean that Lorentz invariance is preserved at the level of the action, only to be broken dynamically. In fact, the theory contains a unit time-like vector (which is called the aether) whose norm is fixed by a Lagrange multiplier term in the action. This entails the fixing of a preferred rest frame at each space-time point. In fact, mimetic gravity itself dynamically violates Lorentz symmetry because the gradient of the mimetic field fixes a preferred direction in space-time.

To be precise, mimetic gravity is in correspondence with a particular version of the Einstein-aether theory, namely, the scalar Einstein-aether theory [77, 145]. In this theory, the role of aether is played by the gradient of a scalar field, precisely as occurs in mimetic gravity: the role of aether is played by the four-velocity vector which in turn is the gradient of the mimetic field. It should be noted that the scalar Einstein-aether theory is quite different from the original vector theory. Moreover, in the case where the potential for the potential for the scalar field in such theories (corresponding to the potential for the mimetic field in mimetic gravity) is constant, the model corresponds to the IR limit of projectable Hořava-Lifshitz gravity, which we will comment on further in Section 3.5.2.

3.5.2. Hořava-Lifshitz Gravity

Recall that Hořava-Lifshitz gravity [338] (HLG hereafter) is a framework and candidate theory of quantum gravity, wherein gravity is made power-counting renormalizable by altering the graviton propagator in the UV. This is achieved by abandoning Lorentz symmetry as a fundamental symmetry of nature, in favour of a Lifshitz anisotropic scaling in the UV. For an incomplete list of references concerning further work in Hořava-Lifshitz gravity, see, for instance, [339–363] and references therein. If the lapse function in HLG is only a function of time, that is, , the theory takes the name of projectable Hořava-Lifshitz gravity. Recall that in the Arnowitt-Deser-Misner decomposition of space-time [445]:the function takes the name of lapse.

Previous work has shown that the IR limit of the nonprojectable version of Hořava-Lifshitz gravity can be obtained from the Einstein-aether theory (the vector version) by requiring that the aether be hypersurface orthogonal: that is, , where is a scalar field and is chosen in such a way to ensure has unit norm. In this case is responsible for the preferred foliation in Hořava gravity, and is the lapse function. If is set equal to coordinate time, the resulting action is that of nonprojectable Hořava gravity [357, 443].

Further requiring that , where is a scalar field, reduces the vector Einstein-aether theory to the scalar Einstein-aether theory [446]. This condition implies that , which upon identification of with coordinate time corresponds exactly to the defining condition for projectable Hořava gravity. In this case, the unit norm constraint cannot be solved for a generic but has to be imposed at the level of the action, for instance, via a Lagrange multiplier term . Therefore, the condition under which is invariant is that be invariant under a shift symmetry: , which shows why the equivalence between scalar Einstein-aether theory and mimetic gravity fails if a nonzero potential for the mimetic field is included [446]. Notice also that, as is known, dark matter emerges as an integration constant in the IR limit of projectable Hořava-Lifshitz gravity [354]. Given the correspondence between this theory and mimetic gravity, then, it should not come as a surprise that dark matter emerges in a similar fashion within the framework of mimetic gravity, as a purely geometrical effect.

A complete proof of the equivalence between mimetic gravity and the IR limit of projectable Hořava-Lifshitz gravity was presented in [135]. In particular, it was shown that the action for the IR limit of projectable Hořava-Lifshitz gravity can be written aswhich we immediately recognize as the action for mimetic gravity without potential, with the addition of a higher derivative term which, as we shall see shortly, was actually added shortly after the original formulation of the theory on purely phenomenological grounds, to cure some undesirable properties at the perturbation level. In Hořava-Lifshitz gravity, is a Lagrange multiplier which had been introduced to enforce the projectability condition. This equivalence was demonstrated rigorously in Appendix A of [135]. Therefore, mimetic gravity appears in the IR limit of a candidate theory of quantum gravity.

3.5.3. Covariant Renormalizable Gravity

Recall that HLG achieves power-counting renormalizability by breaking diffeomorphism invariance. However, this breaking appears explicitly at the level of the action. This has been at the center of criticism, which has related this explicit breaking to the appearance of unphysical modes in the theory which is coupled strongly in the IR [350, 351, 357–359]. Therefore, in order to circumvent these issues it would be desirable to have a theory which preserves the renormalizability properties of Hořava gravity in the UV but retains diffeomorphism invariance at the level of the action. Therefore, diffeomorphism invariance should be broken dynamically in the UV.

An example of such theory has been presented by Nojiri and Odintsov [356, 361, 362] (see also [364–372]) and comes under the name of covariant renormalizable gravity (CRG henceforth). The theory features a nonstandard coupling of a perfect fluid to gravity. When considering perturbations around a flat background, this nonstandard coupling dynamically breaks diffeomorphism invariance. The price to pay is the presence of this exotic fluid, which could have a stringy origin. In particular, one formulation of CRG introduces the fluid via a Lagrange multiplier term, precisely as done in mimetic gravity. For this reason, in the limit where the nonstandard coupling of the fluid to gravity disappears, one recovers the action of mimetic gravity. In one of its equivalent formulations, the action of CRG is given by [356, 361, 362]where and are constants and to recover the mimetic constraint. In particular, the above case corresponds to the version of CRG, where is the dynamical critical exponent which quantifies the degree of anisotropy between space and time in the UV regime of the theory.

Following the initial proposal by Nojiri and Odintsov, other CRG-like models were studied in recent years. For instance, one particular CRG-like model was studied by Cognola et al. in [367], where black hole and de Sitter solutions were also studied. The action of such CRG-like theory takes the following form [367]:In the above action, the Horndeski-like coupling of the scalar field to curvature (for the case ) has been considered copiously in the literature; for an incomplete list see, for instance, [420, 426, 433, 447–455] and references therein.

In the continuation of our review, after a brief interlude on perturbations in mimetic gravity, we shall consider a case study of a mimetic-like theory. Our choice of case study will fall upon the CRG-like model of Cognola et al. as defined by the action in (103). We shall study in detail its cosmological solutions and perturbations around a flat background. A pathological behaviour of the model when considering perturbations around a flat FLRW background will be cured by appropriately modifying the model.

4. A Brief Interlude: Cosmological Perturbations in Mimetic Gravity

Before proceeding to the case study of a specific mimetic-like model, it is mandatory to discuss the issue of perturbations in mimetic gravity. Recall in Section 2 that we remarked that mimetic gravity does not properly belong to the class of modified theories of gravity with a full extra scalar degree of freedom. Instead, the additional scalar is constrained by the Lagrange multiplier, a situation quite different from, for instance, that of Horndeski models which have a proper scalar degree of freedom.

However, it is clear that the Lagrange multiplier kills the wave-like parts of the scalar degree of freedom: in other words, given that the constraint takes out any higher derivative, it is not possible to have oscillating (wave-like solutions). As a consequence, we can already envisage that the sound speed in the minimal mimetic gravity model will satisfy . This implies, as we have anticipated, that there are no propagating scalar degrees of freedom in the theory.

The property of vanishing sound speed in mimetic gravity can be rigorously demonstrated, by considering small longitudinal perturbations around a flat background [73]. The line element is then given by and are functions of the space-time coordinates such that , and , . Here, we used the conformal Newtonian gauge. Moreover, from field equations we also haveCorrespondingly, given that around a flat FLRW background the mimetic field plays the role of “clock,” we perturb it aswhere , which together with (3) implies thatThen, by perturbing the components of the Einstein equations, a relatively straightforward calculation [73] shows that the evolution equation for the perturbation to the mimetic field, , satisfies the following equation:The aforementioned pathological property of perturbations in mimetic gravity can be noticed from (108). The evolution equation for perturbations to the mimetic field does not depend on the Laplacian (or in general on spatial derivatives) of the latter. In other words, there is no term of the form in (108), which implies that the sound speed is identically zero. This means that, even when pressure is nonvanishing, the dust degree of freedom induced in mimetic gravity behaves as dust with zero sound speed, and as such quantum perturbations to the mimetic field cannot be defined in the usual fashion. Else, they would fail in providing the seeds for the observed large-scale structure of the universe which grow via gravitational instability. We remark once more that this behaviour is not unexpected, given that the condition enforced by the Lagrange multiplier eliminates wave degrees of freedom. Moreover, this fact has been shown in all generality for mimetic Horndeski models in [128].

In order to have a theory whose quantum perturbations can be defined in a sensible way, the minimal action for mimetic gravity has to be modified, for instance, by introducing higher derivative (HD) terms. As an example, consider the following action [73]:The corresponding equations of motion readwhere . Thus, identification of the mimetic field with time on a flat FLRW background implies that . There are two main effects of the introduction of the specific HD term on the theory: one at the level of the background and one at the level of perturbations. At the level of the background, it can easily be shown that the Friedmann equation is modified from (32) to [73]which corresponds to a renormalization of the amplitude of the potential. At the level of perturbations, it can be shown that (108) is modified to the following [73]:whereTherefore, the addition of the higher derivative term results in a small but nonvanishing sound speed, which implies that the behaviour of mimetic matter deviates from the usual perfect fluid dust. The nonvanishing sound speed also results in the possibility of defining quantum perturbations in a sensible way. It is beyond the scope of our review to provide further technical details on the matter, but it can be shown that the simple model we have just described is capable of producing red-tilted (i.e., ) scalar perturbations which are enhanced over gravity waves (implying a small value of ), which is consistent with observations from Planck and BICEP2/Keck Array [73].

In concluding this brief interlude, let us also spare a few words on further modifications of mimetic gravity involving higher derivative terms. These have been studied in [84, 85], in particular by adding the previously discussed term, as well as a term proportional to . It was noticed that these terms affect the growth of perturbations below the sound horizon, in particular suppressing the growth of those with large momenta. The result is the presence of a cut-off in the matter power spectrum for perturbations below a certain wavelength. On larger scales, instead, the predictions for the matter power spectrum match those of collisionless cold dark matter.

The suppression of power on small scales is particularly intriguing in the light of the observation that the collisionless cold dark matter paradigm appears to suffer from a number of shortcomings on subgalactic scales. The core-cusp problem refers to the discrepancy between -body simulations of collisionless cold dark matter, which predict a cuspy profile for the dark matter halo in galaxies, and observations which instead suggest a cored profile towards the center. The discrepancy is particularly large for dwarf galaxies, but indications that a cored profile is favoured for larger galaxies as well persist (see, e.g., [456–460]). Moreover, these same simulations predict an abundance of substructure which is approximately 10 times larger than what we actually observe, an issue which is referred to as the “missing satellites problem” (see, e.g., [461–463]). Although this problem is most acute for satellite galaxies, it exists for field galaxies as well (e.g., [464]). To make matters worse, the most massive subhalos in these simulations find no observed counterpart, despite one would expect star formation to be more efficient within them: this is known as the “too big to fail problem” (see, e.g., [465–468]). For an incomplete list of comprehensive reviews on these issues, refer, for instance, to [469, 470] and references therein.

Several approaches to solving these problems exist in the literature. If one insists that dark matter is cold and collisionless, then an important role must be played by the baryonic content of the universe. In fact, it has been argued in several works that baryonic feedback processes (see, e.g., [471–480]) due to supernovae or to the stellar and gaseous content of galaxies, or dynamical friction between dark matter and baryon clumps (see, e.g., [481, 482]), can in principle solve or at least alleviate these problems. Alternatively, the small-scale discrepancies might be taken as an indication that something is lacking in the collisionless cold dark matter picture, with DM possibly having sizeable self-interactions. This approach has been undertaken in a number of works; see, for instance, [483–492] and references therein. In particular, a possibility which has received a lot of attention recently is that where the paucity of structure on small scales is explained by modifying the properties of dark matter in such a way that the resulting matter power spectrum is suppressed at large wavenumbers, by coupling the dark matter to a bath of dark radiation (e.g., a massless or light dark photon) or to a light scalar, which delays kinetic decoupling; refer, for instance, to [493–507] and references therein.

A different mechanism, but with similar outcomes, occurs in mimetic gravity. Namely, the suppression of small-scale power, operated by the higher derivative terms, has the potential to solve the missing satellite problem and the too big to fail problem, as shown in [84]. Moreover, the same higher derivative terms could potentially also cure the caustic singularities from which mimetic gravity suffers. The reason for this is that the higher derivative terms effectively correspond to terms parametrizing dissipation, or viscosity, which emerges due to the fact that the velocity dispersion of the dust is now nonzero. Actually, in [84] it was also argued that these same dissipative terms could alleviate the core-cusp problem. Although these discussions remain preliminary, it is very interesting to note that modifying the action for mimetic gravity by the addition of higher derivative terms could provide a solution to the small-scale structure puzzles of collisionless cold dark matter.

5. A Case Study: Mimetic Horndeski Covariant Hořava-Like Gravity

Having discussed in detail the physics behind mimetic gravity, and many of its extensions, we now provide a detailed case study of a specific mimetic model. Our choice falls on the covariant Hořava-like theory of gravity first discussed by Cognola et al. [367], with action defined by (103). We will argue that this model can be viewed as a mimetic Horndeski model. We will study the model, its background solutions, and scalar perturbations in detail. Not unexpectedly, we will find that the sound speed of the minimal model is vanishing, and thus quantum perturbations cannot be defined in the usual way. To circumvent this problem, we modify the model by the addition of higher order terms, repeating the analysis of background solutions and scalar perturbations, and show that it will be necessary to go beyond the Horndeski framework in order to have a nonvanishing sound speed. The discussion in this section is largely based on the work of [131].

5.1. Mimetic Covariant Hořava-Like Gravity

Let us start from the action of the CRG-like model first discussed by Cognola et al. [367], which is given by the following:where , are arbitrary constants, is the cosmological constant, is a natural number, is a Lagrange multiplier, and determines the constraint imposed on the gradient of the cosmological field . Recall that the Horndeski-like coupling in the action above, for , has been considered in several works (e.g., [420, 433, 447–455]). If we set , we see that the constraint on the gradient of the scalar field corresponds precisely to that of mimetic gravity (3). Thus, we can extend the model to include the mimetic field by adding a potential for the scalar field, which is now made dynamical. The action now readswhere we have set the dimensionless parameter and added a potential for the mimetic field, . Variation with respect to immediately leads to (3), whereas variation with respect to the field modifies (7) towhere we define the following quantities:Finally, variation of the action with respect to the metric leads to the gravitational field equations, which for this theory readwhere we have defined the differential operators:where is the d’Alembertian operator. Note that in the above (118), is a tensor that will not play any role if does not depend on the metric, which we assume is our case; thus, we can drop it from the gravitational field equations. Finally, the form of the Lagrange multiplier can be determined from the trace of (118):We will consider the case with , which is particularly interesting given that it is equivalent to a specific instance of a mimetic Horndeski model.

5.1.1. Cosmological Solutions

Let us now consider cosmological solutions, in particular considering a flat FLRW metric (10). If we take the hypersurfaces of constant time to be equal to those of constant , by making use of the mimetic constraint (3), we see that the field can be identified (up to an integration constant) with time:As in the original mimetic gravity, the scalar field can induce an effective cold dark matter component given that, from (7), one has , and recall that gives the energy density of the perfect dust-like fluid induced by mimetic field (or, more precisely, by its gradient which plays the role of 4-velocity field). However, the introduction of additional terms depending on the field in the action (115) changes these features. Here, we would like to study the behaviour of mimetic field and the cosmological solutions of the more involved theory (115).

Let us begin by considering the tensor in (117), which readsThe, from (116) we deriveGiven the mimetic constraint (3), this equation is a consequence of the two EOMs inferred from (118):In our analysis, as the second independent equation, together, with (123), we choose the trace equation (120): namely,For the simplest case , from (123) we getwhere is an integration constant. Thus, in the limit , we recover the original mimetic gravity model of [69], while when we may interpret this equation as a generalized Friedmann equation, with determining the amount of dark matter in the universe (given that its contribution scales with scale factor as , as is expected for dust). Moreover, in the limits given byif we neglect the cosmological constant , as well as the integration constant , one gets from (126)Thus, if is set at the scale of the early-time acceleration, models with lead to cosmological solutions for inflation.

In the case where we set , the parameter is dimensionless and the EOMs are second order, which is not surprising given that our model is a special case of Horndeski’s theory. Therefore, we can incorporate the cosmological constant in the potential and obtain from (123), (126) thatwhere we have taken into account that on the FLRW space-time, so that we can replace the potential derivative of the field with its time derivative.

Let us manipulate (125) further, which givesshowing that the model we are considering is essentially equivalent to the model proposed in [73]. Since (130) is a consequence of (131) and (132), we may choose to infer the cosmological solutions from (132) only. This equation is a nonlinear Riccati type equation and can be transformed in the linear second-order differential equation:by introducing the Sturm-Liouville canonical substitutionLet us discuss some examples by starting from (133). If is constant, we recover the de Sitter solution:Another well-motivated choice for the potential is a quadratic one:where is a constant Hubble parameter and , are dimensional constants (in the specific case ). After the identification the explicit solution for is found to bewith constant and a fixed time. The Hubble parameter is given byand we see that, for close to , one has a quasi-de Sitter expansion, while for large , the Hubble parameter tends to vanish. This solution corresponds to a Starobinsky-like accelerated expansion [188] in the Jordan frame and gives an interesting inflationary solution.

Let us provide one final example with a potential given bywhere is a constant with mass-dimension 1. The corresponding solution is given byThis solution represents a cosmological bounce with and shows that the mimetic field may act as a phantom fluid.

5.1.2. Cosmological Scalar Perturbations

In this section, we will consider the scalar perturbations around the FLRW metric (10) in the model defined by (115) with , analyzed in the previous section. The perturbed metric in conformal Newtonian gauge is given by (104). Concerning perturbations to the mimetic field, we obtain once more (106) which leads to (107) once more. Note that in this case the identity (105), which is valid for the original mimetic dark matter model, is no longer true. Furthermore, we notice thatFrom the -component of (118), we obtainwithTherefore, we obtainFrom -component of (118) we derivewithwhere the last equality is a consequence of (131), (132).

Finally, we can obtain a closed equation for , which readsFrom the above we immediately read that the sound speed is vanishing, given that there is no dependence on the Laplacian of . The implications are that in these kinds of models, scalar perturbations do not propagate (as in the original mimetic model of [73]), rendering the usual definition of quantum perturbations quite problematic. Given that this feature is rooted in the model being of the mimetic Horndeski form [128], to address the problem, either we must take in (103) or we must modify the original action along the lines of [73, 84].

5.2. Modified Higher Order Mimetic Horndeski Model

With our goal being that of addressing the problem of scalar perturbations, we modify the model given by (115) for asThe original model given by (103) for is recovered for and by using the following identity [508]: On the other hand, for generic values of , , , the action (148) describes a higher order derivative model in the scalar sector with field equations at the fourth order; namely,where we adopted the notationLet us explore some cosmological applications. On a flat FRW metric (10), (132) is modified as follows:whereIt is understood that by setting and , we recover (132). Similarly to how we proceeded in the previous section, we can rewrite (152) aswhere we introduced the auxiliary function :By choosing the following potential: with , , dimensional constants, we recover the Starobinsky-like solution (137) and (138). If we choose the following potential:with being a constant, we recover the bounce solution (140). Therefore, we can recover all the solutions of the Horndeski-like model previously analyzed. We will now study the scalar perturbations around FRW space-time.

5.2.1. Cosmological Scalar Perturbations

If we consider the perturbed metric (104) and the perturbed field (106), we still recover (107) and (141). Therefore, from the -components, of (150), we obtain nowMoreover, from the components or , of (150), we derivewhere the nonvanishing squared sound speed readsWe immediately note that when and we recover (147). Since only for does the model fall within the Horndeski class, we can consider as being the Horndeski breaking parameter. This also confirms that, in order to obtain a nonvanishing sound speed, it is necessary to go beyond the mimetic Horndeski framework: that is, . The result is analogous to that obtained in [73, 84].

6. Spherically Symmetric Solutions in Mimetic Gravity

In this section, we will explore static spherically symmetric solutions (SSS) in mimetic gravity. To do so, let us return to the general formulation of mimetic gravity with action given by (11). Variation of the action with respect to the metric with the mimetic constraint (3) simply leads tonamely, (118) with and . The trace of this equation readsThus, from (161), we getThe continuity equation of the mimetic field is derived aswhich is automatically satisfied when (163) holds true.

In this chapter we will consider pseudo-SSS space-times, whose general topological formulation is given bywhere , are functions of the radial coordinate and the manifold is a sphere when , a torus when , or a compact hyperbolic manifold when . The Ricci scalar in this case readswhere the prime denotes a derivative with respect to . The symmetries of the EOMs require the mimetic field to be a function of only: namely, , and from (3) one hasleading to a pure imaginary expression for the field, which is to be expected from a time-like vector with temporal component equal to zero. Therefore, it is clear that the correspondence with the dark matter is only formal, which justifies the introduction of the potential to make the mimetic field dynamical and possibly reproduce the dark matter phenomenology, since the four-velocity vector in (9) cannot be physical. This has been explored in [119] and will be discussed in Section 7.

The - and -components of the field equations (164) lead to (to derive the EOMs, one may plug the expression for the Ricci scalar directly into the action to obtain, after integration by parts [509], . Thus, the derivatives with respect to and with (3) lead to (168) and (170), while the derivation with respect to leads to (171)):where is given by (162).

We can also rewrite (169) with (167)-(168) asFinally, from (164), one finds the following:As a first investigation, we will consider the case where potential is equal to zero, that is, vacuum solutions.

6.1. Vacuum Solutions

In this subsection we set . From (168) we immediately see thatwith being a mass scale, either positive or negative, while the second field equation (169) with (162) and (166) leads towith , being constants and a length scale. If , namely, and , we recover the topological Schwarzschild solution of General Relativity. When , we can pose . Thus, the solution for flat topology () is as follows:where has to be positive to preserve the metric signature. The Ricci scalar is nonzero and readsThis solution presents a naked singularity at , as for the corresponding topological case of the Schwarzschild metric.

The spherical case () is the more interesting and is given bywhere we have introduced the polar coordinates , in the angular part (see also [79]). In this case the Ricci scalar readsFor the metric is regular everywhere, but when , the solution is regular and preserves the signature; only in the region the solution is regular and preserves the signature; since for the metric coefficient acquires an imaginary part. In this respect, we note that the special choice when leads towhich is real and negative.

In general, we observe that when , the point cannot represent a horizon like for the Schwarzschild space-time, due to the fact that diverges but and the Ricci scalar are regular (we also cannot associate with any thermodynamical quantity such as temperature which diverges at ); moreover, for the metric becomes imaginary (or, in the special case , it acquires the signature (−+++)) and does not allow a wormhole description.

Finally, the topological case readsand the Ricci scalar is given byIn this case, for , the metric coefficient is negative and the solution is unphysical. On the other hand, for , we obtain when , but becomes imaginary (except for the special choice ).

6.2. Nonvacuum Solutions

In this subsection, we will consider the case . We immediately see that, for with a cosmological constant, one solution of (168) and (169) with is the topological Schwarzschild de Sitter solution:with , being constants.

In considering other solutions, we will take the spherical case in (165). By using (167) one readily obtains from (170) and (171) thatwhere we treat the potential as a function of . These two equations with (168) can be used to reconstruct the potential when a choice for is made. We will now provide some examples of the reconstruction technique.

Let us consider a linear modification to the Schwarzschild metric:with , being constants whose mass-dimension is positive, such that from (168) we obtainThe corresponding solution for the mimetic field is found to be an elliptic function, and explicit expressions can be given only in limiting cases. When , one can neglect the linear correction in (184) to recover the solution in (172) and (173) with . Thus, the mimetic field readsor, by expanding the result around ,with . The explicit form of the potential at , is found to beOn the other hand, for large distances, we may ignore the Newtonian term in (184) and from the second equation in (183) we derivewith dimensional constants. We can choose and , and in the given limit the metric simply readsThe corresponding expressions for the field and the potential are given byHere, we must require in order to guarantee the positivity of . The behaviour of the potential can be derived by interpolating the expressions in (188) and (191).

The metric under investigation reduces to the usual Schwarzschild space-time for short distances, while at large distances its 00-component behaves asTherefore, the corresponding Newtonian potential acquires linear and quadratic contributions with respect to the Newtonian solution. The quadratic correction can be viewed as a negative cosmological constant in the background and can be ignored if is sufficiently small. On the other hand, the linear term (for ) could help in explaining the inferred flatness of galactic rotation curves, which has been interpreted as one of the key evidences for the presence of dark matter. We will return to this issue in the next section and analyze the problem more closely there, where we will provide a metric whose cosmological constant-like contribution is independent from the linear one.

As a second example of reconstruction procedure, we consider the following ansatz:with being once more a positive dimensional constant. From (168) we getand the metric is fixed by making use of the second equation in (183) which leads toThe metric signature for is preserved when , and by choosing we obtainThe field and the potential are found to beThe radial coordinate is always real and positive.

The above metric (196) is very interesting, given that it may be used to describe a traversable wormhole (see, e.g., [510–534] and references therein; see also [535–549] for recent work on the subject), where the space is divided into two spherical holes connected by a “throat” located at . Moreover, the metric satisfies the following requirements:(1) and are well defined for all .(2) is regular on the throat with and .(3) and for all .(4)Given , we have .

The above correspond to the traversability conditions [512, 515], and thus we conclude that the wormhole described by our solution is traversable. We also note that our space-time is asymptotically flat.

In the next section we will use a different approach to fix , which encodes the physical Newtonian potential, in order to apply mimetic gravity to reproduce the phenomenology of galactic rotation curves.

7. Rotation Curves of Galaxies in Mimetic Gravity

So far we have seen that, within mimetic gravity, dark matter emerges as a geometrical effect at a cosmological level, that is, in the form of a perfect fluid whose energy density decays as . However, one of the first (and probably most renown) clues as to the presence of dark matter came not from cosmological scales but from astrophysical ones. In fact, it was the observations of Vera Rubin, along with her colleagues Kent Ford, David Burstein, and Norbert Thonnard, that galactic rotation curves were asymptotically flat or even slightly growing with radius , far beyond the region where luminous matter is present, which signalled the presence of something unexpected [4–6]. In fact, if luminous matter was the only thing responsible for the shape of the inferred rotation curves, simple classical mechanics calculations dictate that such curve should fall as , which clearly did not match what was observed. Therefore, it would be interesting to see whether it is possible, within mimetic gravity, to explain the shape of galactic rotation curves, thus fulfilling the requirements of a successful theory of dark matter on both cosmological and astrophysical scales.

The first solution to this problem was found in [119]. The key idea is closely related to that of MOND [550]: to introduce a new scale in the theory, relating it to either the scale where predictions of Newtonian gravity fail or to a scale intrinsically present in the rotation curves data. The work of [119] has drawn from two classes of examples in the literature: the static solution of conformal Weyl gravity [551], used in such context in [552–555], and solutions in gravity [556], used in [557–559]. Let us begin by considering the first cases, which feature linear and quadratic corrections to the Schwarzschild metric.

For our purpose, it is convenient to redefinewithin the metric (165) with : namely,Sinceby using (168) the second equation in (183) readsWe choose the following ansatz for :with , , and being positive dimensional constants. Given the metric element , the Newtonian potential readsThus, the Newtonian potential associated with (202) can be described in the following way: (1)At small distances, the metric leads to a classical Newtonian term .(2)At very large distances, the “cosmological constant” term emerges, reflecting the fact that the metric is immersed in a cosmological de Sitter background.(3)At intermediate distance the linear term appears, revealing a new feature with respect to the Schwarzschild de Sitter metric at galactic scales.

At intermediate galactic scales, we can safely assume that the term is negligible; hence, the Newtonian potential (203) readsCorrespondingly, the rotational velocity profile readswhere we reintroduced the speed of light and is the contribution expected from the luminous matter component. Therefore, on sufficiently large scales, does not fall-off as per the Keplerian result but increases slightly as . This is particularly true for galaxies where the falling Newtonian contribution cannot compete with the rising determined by the term (depending of course on the size of ), which occurs for small and medium sized low surface brightness (LSB) galaxies. In fact rotation curve results for such galaxies exhibit precisely this behaviour (see, e.g., discussion in [554] in the context of conformal Weyl gravity): that is, the rotation curves of these galaxies start rising immediately.

The situation is different for sufficiently extended galaxies, for instance, large high surface brightness (HSB) galaxies. For these galaxies the Newtonian contribution might be sufficient to complete with the rising linear term, . This leads to a region of approximate flatness before any rise starts and is consistent with the data for such galaxies (see, e.g., [554]). Moreover, for these galaxies might be sufficiently large that the de Sitter term should be taken into account. Because of the negative sign, the effect of this term is to reduce the velocity; thus, the rotational velocity profile is given byClearly, sufficiently far from the center of such galaxies, the quadratic term takes over and arrests the rising behaviour driven by the linear term. This is in perfect agreement with data from HSB galaxies, which are large enough to feel the effect of the de Sitter term. Moreover, the negative sign in front of this term has another important implication: given that cannot go negative, bound orbits are no longer possible on scales greater than . This could provide a dynamical explanation for the maximum size of galaxies, determined by the interplay between the linear () and the quadratic () terms.

Let us turn to the question of reproducing such behaviour in mimetic gravity. In order to reconstruct the complete form of the metric (199), we must use (201) to derivewith being a constant. The (on-shell) form of the potential inferred from (168) is quite involved and results inMoreover, one uses the first equation in (183) to recover the behaviour of the field which can be found only in the limiting cases [119].

In the limit (corresponding to small distances), one hasThus, if we setwe havenamely, we recover the vacuum Schwarzschild solution of General Relativity. Analogously to (186) and (187) we derivewith . In this case the potential behaves as given in the following:In the limit (cosmological scales) we obtainwhich corresponds to the static patch of the de Sitter solution. Then, the field assumes the formWe note that as long as , where is the cosmological horizon of the de Sitter solution with positive cosmological constant. The corresponding behaviour of the potential isThis result with is consistent with [73].

Finally, in the limit of galactic scales, the metric readsIn this case, the field is given bywhile the potential in (208) behaves asThe explicit reconstruction of the potential at the galactic scale must take into account that only the solutions with the positive sign inside yield a positive radius for . Thus, we obtainIn conclusion, we are able to describe, in limiting cases, the behaviour of the potential which leads to the solution in (202) and (207) with the corresponding Newtonian potential in (204). The considered solution turns out to be the Schwarzschild solution at small distances, the static patch of the de Sitter space-time at cosmological distance, and, most intriguingly, presents a linear term at the galactic scales which can address the problem of galactic rotation curves in mimetic gravity. We note that, from the dependence of the potential on the field, the potential is always real, as it is required for consistency.

To fix the values of , in the metric, one needs the data from rotation curves of galaxies. Since at the galactic scale our metric reproduces the Newtonian potential (204), while asymptotically it turns out to coincide with the de Sitter metric, we can follow the analyses of [554, 555], where the same potential has been found as a result of the Riegert solution [551] of conformal Weyl gravity. Thus, we adopt the results in [554, 555], where the same parameters were fitted to rotation curves. The total sample fitted consists of 138 galaxies, 25 of them being dwarf galaxies. Some of the galaxies in the sample are sufficiently extended as to be sensitive to the de Sitter quadratic term. We refer the reader to [554, 555] and references therein for data on the galaxies included.

As per the analysis of [554, 555] the potential we reconstructed yields an excellent fit to the rotation curves, with a reduced of ≈1. In fact, the linear and quadratic corrections to the Newtonian potential capture the falling and rising features in the rotation curves quantitatively rather than barely qualitatively. The best-fit values to our and parameters in mimetic gravity (202) turn out to be [554, 555]:The parameter is best expressed as ~(100 Mpc)⁻², suggesting that it is most important on scales of large galaxies or clusters.

We previously mentioned that the idea adopted to fit rotation curves resembles that of MOND, that is, to introduce a new scale in the theory, which could be a scale intrinsically present in the data. Let us elaborate on this point more quantitatively. Considering the measured distance and rotational velocity of the outermost data in the rotation curves, it can be shown that the combination for each of the galaxies in the sample is, within better than an order of magnitude, very close to the best-fit value for . In other words, the rotation curve data contain a preferred scale, which we introduced through the non-Newtonian correction.

We conclude this chapter by considering, for completeness, the case of a general power-law correction:with being a positive real parameter, with the full metric given by (201); namely,In deriving this last expression we have set the integration constant in such a way as to recover the Schwarzschild solution in the limit , namely, at short distances. On the other hand, at large distances, in the limit , we once more find the static patch of the de Sitter solution. Instead, at galactic scales, in the limit , we obtainThe potential can be found only in an implicit way and its on-shell expression results inThe trivial case corresponds to a de Sitter-like solution with negative cosmological constant (for ):In this case the field is given by and the potential reads .

The non-Newtonian correction we have considered has recently been studied in the context of gravity in [556]. In the low-energy limit of power-law gravity, a correction to the Newtonian potential emerges, with related to the power asStability of the potential at large distances and constraints from Solar System tests require . In [556] the above model was fitted to 15 LSB galaxies, and an excellent fit is achieved, with best-fit value (corresponding to ). Interestingly, the study [557] performs a similar fit to two objects in [557]. The two objects, whose rotation curves cannot be explained successfully by particle dark matter, are the following: the dwarf galaxy Orion and the low luminosity spiral NGC 3198 (which had not been analyzed in [556]). In this case, an excellent fit is obtained for .

Let us conclude by making an important remark. Although the potential leading to the chosen non-Newtonian corrections has only been given implicitly, the form of has been provided explicitly, demonstrating how the desired behaviour for galaxy rotation curves, consistent with data, can be obtained in mimetic gravity. This is a notable result, as it implies that it is possible to reproduce the behaviour of dark matter not only on cosmological scales in mimetic gravity, but on astrophysical scales as well, solving the galaxy rotation curves problem, which is one of the pillars of evidence for dark matter.

8. Conclusion

Mimetic gravity has emerged as an interesting and viable alternative to General Relativity, wherein the dark components of the universe (underlying dark matter, the late-time acceleration, and inflation) can find unified geometrical explanation and interpretation. The theory is related to General Relativity by a singular disformal transformation, which is the reason behind its exhibiting a wider class of solutions. Here, we have reviewed the main aspects of mimetic gravity, beginning by placing it in the wider context of theories of modified gravity. After having reviewed the underlying theory behind mimetic gravity, we have studied some of its solutions and extensions, such as mimetic gravity, mimetic unimodular gravity, and others, focusing on the reconstruction technique which allows the realization of numerous wishful cosmological expansion scenarios. After having discussed the issue of perturbations within the theory, we have considered a specific mimetic-like model, namely, mimetic covariant Hořava-like gravity, wherein we applied the concepts discussed in the first part of the review. The final part of the review has been devoted to the study of static spherically symmetric solutions within mimetic gravity, and the application of these to the study of rotation curves within such theory.

The dark components of our universe remain as mysterious as ever. It is possible that we might shed light on the nature of dark matter, dark energy, and inflation, as more data from experiments and surveys pours in the coming years. Thus far, theories of modified gravity, despite their prima facie complications, appear as viable and theoretically well-motivated routes to pursue. This is true for mimetic gravity as well, and thus we expect it to remain an active arena of research in the field of modified gravity in the coming years.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

Sunny Vagnozzi thanks the Niels Bohr Institute, where the majority of this work was completed, for hospitality. The authors have benefited from many discussions with Sergei Odintsov and Sergio Zerbini, whom we wish to thank.