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

Authors: Lorenzo Sebastiani; Sunny Vagnozzi; Ratbay Myrzakulov
Journal: Advances in High Energy Physics
 (2017)

---

## 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 mimeticf(R) 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.


---

## Body

## 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 Section2, 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 L=L(gμν), where gμν is the metric tensor, is given by the following:(1)Eμν=β-gRμν-12gμνR+κ-ggμν,where β and κ are constants and Rμν and R 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 invarianceTherefore, 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 ofscalar 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 f(T) 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 addvector 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 wheretensor 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 byf(R) 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 f(R) gravity in [303, 305, 311–313]. Belonging to the same family are also variations of the former such as f(RμνRμν), f(□R), fR,T, f(R,T,RμνTμν) 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 asf(R/□) and f(Rμν□-1Rμν) 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 satisfiescs=0, 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.


## 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 L=L(gμν), where gμν is the metric tensor, is given by the following:(1)Eμν=β-gRμν-12gμνR+κ-ggμν,where β and κ are constants and Rμν and R 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 invarianceTherefore, 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 ofscalar 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 f(T) 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 addvector 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 wheretensor 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 byf(R) 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 f(R) gravity in [303, 305, 311–313]. Belonging to the same family are also variations of the former such as f(RμνRμν), f(□R), fR,T, f(R,T,RμνTμν) 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 asf(R/□) and f(Rμν□-1Rμν) 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 satisfiescs=0, 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.


## 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 ofscalar 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 f(T) 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 addvector 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 wheretensor 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 byf(R) 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 f(R) gravity in [303, 305, 311–313]. Belonging to the same family are also variations of the former such as f(RμνRμν), f(□R), fR,T, f(R,T,RμνTμν) 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 asf(R/□) and f(Rμν□-1Rμν) 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 satisfiescs=0, 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 gμν in terms of an auxiliary metric g~μν and a scalar field ϕ, dubbed mimetic field, as follows:(2)gμν=-g~μνg~αβ∂αϕ∂βϕ.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 g~μν→Ω(t,x)2g~μν, Ω(t,x) 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:(3)gμν∂μϕ∂νϕ=-1.Thus, the gravitational action, taking into account the reparametrization given by (2) now takes the form:(4)I=12∫Md4x-gg~μν,ϕRg~μν,ϕ+2Lm,where M is the space-time manifold, R≡R(gμν,ϕ) is the Ricci scalar, Lm is the matter Lagrangian, and g≡g(g~μν,ϕ) 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](5)Gμν-Tμν+G-T∂μϕ∂νϕ=0,where Gμν=Rμν-gμνR/2 is the Einstein tensor, with Rμν being the Ricci tensor, while G(=-R) and T 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:(6)T~μν=-G-T∂μϕ∂νϕ.We note that both energy-momentum tensors, Tμν and T~μν, are covariantly conserved, that is, ∇μTμν=∇μT~μν=0 (with ∇μ the covariant derivative), whereas the continuity equation for T~μν with the mimetic constraint (3) leads to(7)∇κG-T∂κϕ≡1-g∂κ-gG-Tgκσ∂σϕ=0.Finally, the trace of (5) is found to be(8)G-T1+gμν∂μϕ∂νϕ=0.It is clear that the above is automatically satisfied if one takes into account (3) even if G≠T. Thus, the theory admits nontrivial solutions and the conformal degree of freedom becomes dynamical even in the absence of matter (T=0 but G≠0) [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 p is given by(9)Tμν=ρ+puμuν+pgμν,uμuμ=-1.Notice that the mimetic stress-energy tensor in (6) assumes the same form of the one of a perfect fluid with pressure p=0 and energy density ρ=-(G-T), 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:(10)ds2=-dt2+at2dx2,where a≡a(t) 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 (G-T), which corresponds to the energy density of the mimetic stress-energy tensor, decays with the scale factor of the FLRW universe as (G-T)∝1/a3. Recall that the energy density of a component with equation of constant state parameter w evolves as a-3(w+1) in an FLRW universe, and hence the evolution in the energy density of the mimetic field corresponds to w=0, 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 g~ 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 as(11)I=12∫Md4x-gR+λgμν∂μϕ∂νϕ+1+2Lm.Variation of the Lagrangian with respect to the Lagrange multiplier field λ leads to (3), while variation with respect to the physical metric gμν yields(12)Gμν-Tμν+λ∂μϕ∂νϕ=0,whose trace, when one takes into account (3), is given by(13)λ=G-T.Thus, 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 that(14)∫tintfindtΩx=tfin-tin,Ωx≡g~αβ∂αϕ∂βϕ.Thus, for a spatially homogeneous mimetic field ϕ=t, ϕ˙=Ω. 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 x≡f(y˙) into an action I(x,t): the class of functions over which the variation is admissible does not only comprise those for which the variation of x is vanishing at the boundary but also those where the integral of the variation of x 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 gμν=gμν(g~μν,ϕ) is noninvertible even for fixed ϕ, owing to the fact that the map g→g~,ϕ 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 I) of the form:(15)∂αg~μνx=αxgμνx,∂αIgμνg~μν,ϕ=0,where α(x) 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 Igμν(g~μν,ϕ). Mimetic gravity can thus be seen as a conformal extension of GR which is Weyl-invariant in terms of the auxiliary metric g~μν. 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 gμν in terms of a fiducial metric g~μν and a scalar field ϕ. The map between the two is defined as a “disformal transformation,” or a “disformation,” and takes the following form:(16)gμν=Aϕ,Xg~μν+Bϕ,X∂μϕ∂νϕ,where X≡g~μν∂μϕ∂νϕ. The functions A and B are referred to as conformal factor and disformal factor, respectively. In general the functions A(ϕ,X), B(ϕ,X) are arbitrary, with A≠0. 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 g~μν and ϕ) reduce to those of obtained by varying with respect to the metric gμν [79].To make progress, it is useful to contract the two equations of motion withg~μν and ∂μϕ∂νϕ. Although we will not perform the steps explicitly, it is easy to show that this leads to the following two equations of motion:(17)ΩA-X∂A∂X-ωX∂B∂X=0,ΩX2∂A∂X-ωA-X2∂B∂X=0,where the two quantities Ω and ω are defined as(18)Ω≡Gμν-Tμνg~μν,ω≡Gμν-Tμν∂μϕ∂νϕ.The determinant of system (17) is given by(19)det=X2A∂∂XB+AX.If the above is nonzero, it is trivial to obtain that the resulting set of equations consists of Einstein’s equation (Gμν=Tμν) 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 A≠0, this immediately determines B, which is of the form:(20)BX,ϕ=-AX,ϕX+Eϕ,where E(ϕ)≠0 is an arbitrary function. We will not show the steps explicitly, which are instead discussed in detail in Section IV 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 A=X and B=0 in (16), and correspondingly E=1 in (20). In general, when the relation defined by (20) exists between the conformal factor A and the disformal factor B, 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]:(21)Ig~,ϕ=∫d4x-g~XRg~+32g~μν∂μX∂νXX,where we have defined X≡g~μν∂μϕ∂νϕ/2. One immediately sees that (21) corresponds to the singular/conformal Brans-Dicke action [270] with density parameter ω=-3/2. 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 g~μν but to 2Xg~μν [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 metricgμν, 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 ϵ=T-G=T+R 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, pϕ, satisfies pϕ=0.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:(22)gμν=-g~μνg~αβuαuβ.The Proca field is made dynamical by adding a Maxwell kinetic term F2 to the action, where F 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) to(23)gμν=-fϕg~μνg~αβuαuβ,where 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 by(24)δϕ=0,δgμν=ϵ∂A∂Xgμν+∂B∂X∂μϕ∂νϕ,where as usual A and B correspond to the conformal and disformal factors. Being A=X and B=0 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.


## 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 that(14)∫tintfindtΩx=tfin-tin,Ωx≡g~αβ∂αϕ∂βϕ.Thus, for a spatially homogeneous mimetic field ϕ=t, ϕ˙=Ω. 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 x≡f(y˙) into an action I(x,t): the class of functions over which the variation is admissible does not only comprise those for which the variation of x is vanishing at the boundary but also those where the integral of the variation of x 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 gμν=gμν(g~μν,ϕ) is noninvertible even for fixed ϕ, owing to the fact that the map g→g~,ϕ 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 I) of the form:(15)∂αg~μνx=αxgμνx,∂αIgμνg~μν,ϕ=0,where α(x) 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 Igμν(g~μν,ϕ). Mimetic gravity can thus be seen as a conformal extension of GR which is Weyl-invariant in terms of the auxiliary metric g~μν. 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 gμν in terms of a fiducial metric g~μν and a scalar field ϕ. The map between the two is defined as a “disformal transformation,” or a “disformation,” and takes the following form:(16)gμν=Aϕ,Xg~μν+Bϕ,X∂μϕ∂νϕ,where X≡g~μν∂μϕ∂νϕ. The functions A and B are referred to as conformal factor and disformal factor, respectively. In general the functions A(ϕ,X), B(ϕ,X) are arbitrary, with A≠0. 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 g~μν and ϕ) reduce to those of obtained by varying with respect to the metric gμν [79].To make progress, it is useful to contract the two equations of motion withg~μν and ∂μϕ∂νϕ. Although we will not perform the steps explicitly, it is easy to show that this leads to the following two equations of motion:(17)ΩA-X∂A∂X-ωX∂B∂X=0,ΩX2∂A∂X-ωA-X2∂B∂X=0,where the two quantities Ω and ω are defined as(18)Ω≡Gμν-Tμνg~μν,ω≡Gμν-Tμν∂μϕ∂νϕ.The determinant of system (17) is given by(19)det=X2A∂∂XB+AX.If the above is nonzero, it is trivial to obtain that the resulting set of equations consists of Einstein’s equation (Gμν=Tμν) 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 A≠0, this immediately determines B, which is of the form:(20)BX,ϕ=-AX,ϕX+Eϕ,where E(ϕ)≠0 is an arbitrary function. We will not show the steps explicitly, which are instead discussed in detail in Section IV 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 A=X and B=0 in (16), and correspondingly E=1 in (20). In general, when the relation defined by (20) exists between the conformal factor A and the disformal factor B, 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]:(21)Ig~,ϕ=∫d4x-g~XRg~+32g~μν∂μX∂νXX,where we have defined X≡g~μν∂μϕ∂νϕ/2. One immediately sees that (21) corresponds to the singular/conformal Brans-Dicke action [270] with density parameter ω=-3/2. 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 g~μν but to 2Xg~μν [129].


## 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 gμν in terms of a fiducial metric g~μν and a scalar field ϕ. The map between the two is defined as a “disformal transformation,” or a “disformation,” and takes the following form:(16)gμν=Aϕ,Xg~μν+Bϕ,X∂μϕ∂νϕ,where X≡g~μν∂μϕ∂νϕ. The functions A and B are referred to as conformal factor and disformal factor, respectively. In general the functions A(ϕ,X), B(ϕ,X) are arbitrary, with A≠0. 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 g~μν and ϕ) reduce to those of obtained by varying with respect to the metric gμν [79].To make progress, it is useful to contract the two equations of motion withg~μν and ∂μϕ∂νϕ. Although we will not perform the steps explicitly, it is easy to show that this leads to the following two equations of motion:(17)ΩA-X∂A∂X-ωX∂B∂X=0,ΩX2∂A∂X-ωA-X2∂B∂X=0,where the two quantities Ω and ω are defined as(18)Ω≡Gμν-Tμνg~μν,ω≡Gμν-Tμν∂μϕ∂νϕ.The determinant of system (17) is given by(19)det=X2A∂∂XB+AX.If the above is nonzero, it is trivial to obtain that the resulting set of equations consists of Einstein’s equation (Gμν=Tμν) 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 A≠0, this immediately determines B, which is of the form:(20)BX,ϕ=-AX,ϕX+Eϕ,where E(ϕ)≠0 is an arbitrary function. We will not show the steps explicitly, which are instead discussed in detail in Section IV 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 A=X and B=0 in (16), and correspondingly E=1 in (20). In general, when the relation defined by (20) exists between the conformal factor A and the disformal factor B, 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]:(21)Ig~,ϕ=∫d4x-g~XRg~+32g~μν∂μX∂νXX,where we have defined X≡g~μν∂μϕ∂νϕ/2. One immediately sees that (21) corresponds to the singular/conformal Brans-Dicke action [270] with density parameter ω=-3/2. 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 g~μν but to 2Xg~μν [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 metricgμν, 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 ϵ=T-G=T+R 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, pϕ, satisfies pϕ=0.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:(22)gμν=-g~μνg~αβuαuβ.The Proca field is made dynamical by adding a Maxwell kinetic term F2 to the action, where F 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) to(23)gμν=-fϕg~μνg~αβuαuβ,where 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 by(24)δϕ=0,δgμν=ϵ∂A∂Xgμν+∂B∂X∂μϕ∂νϕ,where as usual A and B correspond to the conformal and disformal factors. Being A=X and B=0 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:(25)I=12∫Md4x-gR+λgμν∂μϕ∂νϕ+1-Vϕ+2Lm.The equations of motion of the theory are then given by(26)Gμν-2λ∂μϕ∂νϕ-gμνVϕ=Tμν,which, by taking the trace, can be used to determine the Lagrange multiplier:(27)λ=12G-T-4V.When plugged into (26), (27) yields(28)Gμν=G-T-4V∂μϕ∂νϕ+gμνVϕ+Tμν.Of 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 p=-V and ρ=G-T-3V, respectively [73]. One further equation of motion can be derived by varying the action with respect to the mimetic field, which gives(29)∇νG-T-4V∂νϕ=-∂V∂ϕ,when 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](30)1a3ddta3ρ-V=-dVdt,which can be integrated to give(31)ρ=3a3∫daa2V,whereas the pressure remains p=-V. The Friedmann equation can instead be manipulated to the form:(32)2H˙+3H2=Vt,where as usual the Hubble parameter is defined as H≡a˙/a. Further progress can be made by performing the substitution y≡a3/2, which yields the following equation [73]:(33)y¨-34Vty=0.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 y satisfies a linear differential equation. We now proceed to study a few interesting potentials and the corresponding solutions.


#### 3.1.1.V∝1/ϕ2

Let us consider the following potential [73]:(34)Vϕ=αϕ2=αt2.Solving the corresponding (33) and substituting for the scale factor, a≡y2/3, yield the following solution for α≥-1/3:(35)at=t1/31+1+3α1+βt-1+3α2/3,where β is an integration constant. For α<-1/3 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) w≡p/ρ if we recall that the energy density is given by (31) whereas the pressure reads p=-V. Explicit calculation gives the following:(36)w=-3α1+1+3α1-βt-1+3α1+βt-1+3α-2.It is interesting to note that, for α≫1, the EoS approaches w=-1, 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 EoSw~. The scale factor in this situation evolves as(37)a∝t2/31+w~,and hence (31) can be used to deduce the evolution of the energy density of mimetic matter, which decays as(38)ρ=-αw~t2.Given that the pressure of mimetic matter reads p=-V, we immediately see that the EoS for mimetic matter is w=w~, 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 α/w≪1. 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:(39)Vϕ=αϕn=αtn,for which the solution of (33) can be written in terms of the Bessel functions [73]:(40)y=t1/2Z1/n+2-3αn+2tn+2/2.For n<-2 (with n=-2 corresponding to the case we have studied previously) the limiting behaviour of the scale factor is that of a dust-dominated universe, with EoS w=0. For n>-2 and α<0 (which corresponds to a positive pressure), the corresponding solution is a singular oscillating universe. For n>-2 and α>0 instead, the pressure is negative and hence we expect accelerating solutions [73]. In fact, n=0 corresponds to a cosmological constant as expected (the potential is simply a constant), whereas n=2 gives an inflationary expansion solution with scale factor:(41)a∝t-1/3eα/12t2,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 parameterH or equivalently in the scale factor a(t)], find the corresponding y-parameter through a=y2/3, and then invert (33) to find the corresponding potential which can provide the desired expansion:(42)Vϕ=Vt=4y¨3y.As an example, we can consider the following potential [73]:(43)Vϕ=αϕ2eϕ+1,whose corresponding solution for the scale factor is exponential (a∝e-t2) for large negative times and corresponding to EoS w=0 (i.e., dust, a∝t2/3) 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]:(44)Vϕ=αe-κϕ=αe-κt.In this case, the scale factor can be expressed in terms of Bessel functions of the first and second kinds [83]:(45)at=βJ0-3ακe-kt/2+γY0-3ακe-kt/22/3,where β and γ are integration constants and J0 and Y0 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 w=0 (a∝t2/3). At early times the behaviour of the scale factor depends on the sign of α, in particular providing us with an inflationary solution for α>0, whereas the solution is a bouncing nonsingular universe for α<0. A similar behaviour can be obtained if one chooses the potential:(46)Vϕ=αϕ2neκϕ+1=αt2neκt+1,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:(47)Vϕ=4311+ϕ22=4311+t22.As usual, the scale factor can be determined by solving (33), which yields [73](48)at=t2+11+βarctant2/3.If we set the integration constant β to 0, the corresponding energy density and pressure (one again, refer to (31) and recall that p=-V) are given by(49)ρ=43t21+t22,p=-4311+t22.At very early times (large negative t) the EoS approaches w→0; the universe is dominated by dust and contracts. At a certain time corresponding to |t|~1, 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 |β|<2/π [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](50)Vϕ=43αϕ02+ϕ22=43αt02+t22.Although we will not show the solution explicitly, in this case the scale of the bounce is reduced to α/t02 and the duration of the bounce is now t0 [73].


### 3.2. MimeticF(R) Gravity

The next step which was performed by Nojiri and Odintsov is to generalize mimetic gravity to mimeticF(R) 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 F(R) term. The action of the theory is given by [81](51)I=∫d4x-gg~μν,ϕFRg~μν,ϕ+Lm,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 F(R) gravity can be equally written employing a Lagrange multiplier field analogously to (11) [81]:(52)I=∫d4x-gFRgμν-Vϕ+λgμν∂μϕ∂νϕ+1+Lm,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]:(53)0=12gμνFR-RμνF′R+∇μ∇νF′R-gμν□F′R+12gμν-Vϕ+λgαβ∂αϕ∂βϕ-λ∂μϕ∂νϕ+12Tμν.Variation with respect to the mimetic field instead yields the following equation [81]:(54)2∇μλ∂μϕ+dVdϕ=0.As usual, by construction, variation with respect to the Lagrange multiplier gives the mimetic constraint:(55)gμν∂μϕ∂νϕ=-1.As we have mentioned previously, in mimeticF(R) 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 F(R) 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 mimeticF(R) 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 p. Assuming that the mimetic field depends only on time, (53), (54), and (55) can be expressed as follows [81]:(56)0=-FR+6H˙+H2F′R-6HdF′Rdt-λϕ˙2+1+Vϕ+ρ,(57)0=FR-2H˙+3H+2d2F′Rdt2+4HdF′Rdt-λϕ˙2-1-Vϕ+p,(58)0=2ddtλϕ˙+6Hλϕ˙-dVdϕ,(59)0=ϕ˙2-1.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:(60)0=FR-2H˙+3H+2d2F′Rdt2+4HdF′Rdt-Vt+p,which, if we assume that the contribution of ordinary matter is negligible (ρ=p=0), reduces to(61)Vt=FR-2H˙+3H+2d2F′Rdt2+4HdF′Rdt.On the other hand, (57) can be solved for λ as follows:(62)λt=-12FR+3H˙+H2F′R-3HdF′Rdt,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 F(R) which gives the wanted evolution [encoded in H(t) or, equivalently, a(t)]. Alternatively, one can start from a given form of F(R) 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 F(R) gravity is also a solution in mimetic F(R) 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 F(R) gravity. Of course, the exact forms of the mimetic potential or the F(R) 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 F(R) 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, ns and r, respectively. Both the reconstruction (determination of the potential once the form of F(R) and the evolutionary history of the universe are given) and the inverse reconstruction (determination of the form of F(R) 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 F(R) gravity in this respect: namely, the forms of both the mimetic potential and the F(R) function can become extremely complicated. The forms of both the mimetic potential and the Lagrange multiplier increase in complexity as the complexity of the F(R) form increases. Therefore, these and other studies on mimetic F(R) 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 F(R) 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 F(R) 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 F(R), 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 F(R) gravity with the desired inflationary properties would be to choose the simplest possible functional form of F(R) 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 F(R,ϕ) 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 mimeticF(R) 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 F(R) gravity as well. In other words, stable solutions in mimetic F(R) 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 F(R) gravity. An immediate implication of this finding is that, although mimetic F(R) gravity could drive inflation differently from F(R) gravity, the late-time acceleration of the universe in these theories has to coincide with the usual F(R) 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 F(R) gravity compared to conventional F(R) 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 F(R) gravity.As we mentioned above, the conclusions reached about the late-time evolution in mimeticF(R) gravity hold only at the background level. However, in conventional F(R) 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, dF(R)/dR) 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 F(R) gravity the problem is solved by adding power-law modifications by hand.In [153], it was instead argued that in mimetic F(R) 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 F(R) model whose corresponding conventional F(R) model suffered from the oscillations problems. The oscillations die out for redshifts z≤3, 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 F(R) 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 F(R) gravity as a cosmological framework.


#### 3.2.2. MimeticF(R,ϕ) Gravity

The question of constructing a theoretically motivated but at the same time simple model for mimeticF(R) 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 F(R,ϕ) 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 F(R) 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 F(R).


#### 3.2.3. Nonlocal MimeticF(R) Gravity

A further extension of mimeticF(R) 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 F(R) gravity is given by the action [133]:(63)I=∫d4x-gR1+f□-1R+λgμν∂μϕ∂νϕ+1-Vϕ,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]:(64)I=∫d4x-gR1+fψ+ξ□ψ-R+λgμν∂μϕ∂νϕ+1-Vϕ,where ξ is an additional Lagrange multiplier which enforces the constraint on the scalar field ψ:(65)□ψ=R.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]:(66)Rμν1+fψ-ξ-12gμνR1+fψ-ξ-∂ρξ∂ρψ=12∂μξ∂νψ+∂μψ∂νξ-gμν□-∇μ∇νfψ-ξ-λ∂μϕ∂νϕ-gμνVϕ.Instead, variation with respect to the two scalar fields leads to the following equations of motion [133]:(67)□ξ+f′ψR=0,-1-g∂ν-g∂νϕ=12dVdϕ.In [133], this model was studied in detail making use of the reconstruction technique. In particular, two forms for the f(ψ) 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:(68)-g=1.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](69)I=∫d4x-gR-Vϕ-ηgμν∂μϕ∂νϕ+1-λ+λ,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]:(70)0=12gμνR-Vϕ-ηgαβ∂αϕ∂βϕ+1-λ-Rμν+η∂μϕ∂νϕ+12Tμν,whereas variation with respect to the mimetic field yields the usual equation of motion:(71)0=2∇μλ∂μϕ-dVdϕ.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 V(ϕ)=V(t) (as usual, the mimetic field can be identified with time):(72)Vϕ=Vt=at23∫0ϕdtat-2-6Htpt-2dptdt+2-18Ht3-6HtdHtdt+4d2Htdt2.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 H(t) or a(t). 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]:(73)Vϕ=12e2H0ϕH03ϕ.It can be easily shown that it leads to the following solution for the scale factor a(t) and the Hubble parameter H(t):(74)at=eH0t,Ht=H0,that is, a de Sitter cosmology. The functional forms of the Lagrange multiplier are also quite simple and are given by the following:(75)λt=6H021+2e2H0tH0t,ηt=3H02.Another choice for the potential which leads to a physically interesting solution is the following [136]:(76)Vϕ=-8ϕ-2+4/31+w1+5w+2w291+w3,for which the scale factor and the Hubble parameter read(77)at=t2/31+w,Ht=23t1+w.The above solution is that corresponding to an universe dominated by a fluid with EoS w. The two Lagrange multipliers are given by(78)λt=8-3w1+w+t4/31+w1+5w+2w29t21+w3,ηt=42+w3t21+w2.Thus, 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 p as follows:(79)ρ=G-T-4V~,p=-V~,where V~ is defined as(80)V~=-λt-Vt.Furthermore, 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 -g). It is nonetheless possible to present a covariant formulation of unimodular mimetic gravity via the following action:(81)I=∫d4x-gR-ηgμν∂μϕ∂νϕ+1-λ+λϵμνρσ∂μaνρσ,where aνρσ is a three-form field, variation of which gives the constraint ∂μλ=0, 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 with(82)-g=ϵμνρσ∂μaνρσ.Further 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 MimeticF(R) Gravity

A minimal extension of the unimodular mimetic gravity framework we have discussed so far is to consider unimodular mimeticF(R) gravity, which is described by the action [138]:(83)I=∫d4x-gFR-Vϕ-ηgμν∂μϕ∂νϕ+1-λ+λ.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]:(84)gμν2FR-Vϕ+ηgαβ∂αϕ∂βϕ+1-λ-RμνF′R-η∂μϕ∂νϕ+∇μ∇νF′R-gμν□F′R=0.Finally, variation with respect to the mimetic field yields the following equation:(85)0=2∇μλ∂μϕ-dVdϕ.The more complex structure of the equations of motion complicates the reconstruction procedure, that is, the equivalent of (72), which now reads [138](86)Vt=at3/22∫dtat-3/2ft,where the function f(t) is given by(87)ft=-18HtH˙+H2F′R-6H2dF′Rdt-18H˙H3+6Hd3F′Rdt3+6H¨+2H˙HF′R+HH˙dF′Rdt+6H2dF′Rdt-2Hd2F′Rdt2-2H¨+6H˙H+2d2F′Rdt2.Despite 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 F(R). 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]:(88)I=∫d4x-gLgμν,∂λ1gμν,…,∂λ1⋯∂λpgμν,ϕ,∂λ1ϕ,…,∂λ1⋯∂λqϕ+λgμν∂μϕ∂νϕ+1,where p,q≥2 are integers and L 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 b(ϕ)gμν∂μϕ∂νϕ=-1, but for the sake of simplicity we will set b(ϕ)=1 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 b(ϕ)≠1 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:(89)I=∫d4x-gLH=∫d4x-g∑nLn,where Lns read(90)L0=KX,ϕ,L1=-G3X,ϕ□ϕ,L2=G4,XX,ϕ□ϕ2-∇μ∇νϕ2+RG4X,ϕ,L3=-16G5,XX,ϕ□ϕ3-3□ϕ∇μ∇νϕ2+2∇μ∇νϕ3+Gμν∇μ∇νϕG5X,ϕ,where X≡-gμν∇μϕ∇νϕ/2, (∇μ∇νϕ)≡∇μ∇νϕ∇μ∇νϕ, (∇μ∇νϕ)3≡∇μ∇νϕ∇μ∇ρϕ∇ν∇ρϕ, the functions K(X,ϕ), G3(X,ϕ), G4(X,ϕ), G5(X,ϕ) are free, and X denotes differentiation with respect to X.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 K, G3, G4, and G5, as well as in the function b(ϕ) (which, when ≠1, 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:(91)KX,ϕ=c2X,G3X,ϕ=0,G4X,ϕ=12,G5X,ϕ=0.In this case, on a flat FLRW background the solution is given by(92)at=t2/31+w,ϕt=±-αc2lntt0,bϕ=-1ϕ˙2,with t0 and integration constant and α=-8w/3(1+w)2. Thus, we see that the scenario under consideration has reproduced the expansion history of a universe filled with a perfect fluid with EoS w [104]. Another case considered in [104] is the mimetic cubic Galileon model, where the functions take the following form:(93)KX,ϕ=c2X,G3X,ϕ=2c3Λ~3X,G4X,ϕ=12,G5X,ϕ=0,where the cut-off scale is subsequently set to Λ~=1. 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, g~μν 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 by(94)IH=∫d4x-gαXR+□ϕ2-∇μ∇νϕ∇μ∇νϕ+γϕGμν∇μ∇νϕ-βϕ□ϕ,which corresponds to the following choice for the functions discussed previously:(95)KX,ϕ=0,G3X,ϕ=βϕ,G4X,ϕ=αX,G5X,ϕ=0.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],(96)Vϕ=-β+2V0+3V02c1ϕ-ϕ02,gives rise to the following solution:(97)Ht=V0c1t-t0,which represents a regular bounce solution. Another bounce solution can be obtained by considering the following potential [114]:(98)Vϕ=b2c1sinh2bϕcosh2bϕ,for which the scale factor and correspondingly the Hubble parameter read:(99)at=a0coshbt,Ht=bsinhbtcoshbt.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 uμ (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, N=N(t), the theory takes the name of projectable Hořava-Lifshitz gravity. Recall that in the Arnowitt-Deser-Misner decomposition of space-time [445]:(100)ds2=-N2c2dt2+gijdxi+Nidtdxj+Njdt,the function N 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,uμ=N∇μT, where T is a scalar field and N is chosen in such a way to ensure uμ has unit norm. In this case T is responsible for the preferred foliation in Hořava gravity, and N is the lapse function. If T is set equal to coordinate time, the resulting action is that of nonprojectable Hořava gravity [357, 443].Further requiring thatNdT=S, where S is a scalar field, reduces the vector Einstein-aether theory to the scalar Einstein-aether theory [446]. This condition implies that N=N(T), which upon identification of T 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 N but has to be imposed at the level of the action, for instance, via a Lagrange multiplier term ∝λ(∇μS∇μS-1). Therefore, the condition under which dS=NdT is invariant is that S be invariant under a shift symmetry: S→S+dS, 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 as(101)S=SEH+∫d4x-gΣ2gμν∂μϕ∂νϕ+1+γ2□ϕ2,which 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](102)I=∫d4x-gR2-α∂μϕ∂νϕ∇μ∇ν+2U0∇ρ∇ρn∂μϕ∂νϕRμν+U0R2-λ12∂μϕ∂μϕ+U0,where U0 and α are constants and U0=1/2 to recover the mimetic constraint. In particular, the above case corresponds to the z=2n+2 version of CRG, where z 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]:(103)I=12∫d4x-gR-2Λ-αRμν-R2gμν∇μϕ∇νϕn-λ2gμν∇μϕ∇νϕ+1-Vϕ.In the above action, the Horndeski-like coupling of the scalar field to curvature (for the case n=1) 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.


## 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:(25)I=12∫Md4x-gR+λgμν∂μϕ∂νϕ+1-Vϕ+2Lm.The equations of motion of the theory are then given by(26)Gμν-2λ∂μϕ∂νϕ-gμνVϕ=Tμν,which, by taking the trace, can be used to determine the Lagrange multiplier:(27)λ=12G-T-4V.When plugged into (26), (27) yields(28)Gμν=G-T-4V∂μϕ∂νϕ+gμνVϕ+Tμν.Of 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 p=-V and ρ=G-T-3V, respectively [73]. One further equation of motion can be derived by varying the action with respect to the mimetic field, which gives(29)∇νG-T-4V∂νϕ=-∂V∂ϕ,when 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](30)1a3ddta3ρ-V=-dVdt,which can be integrated to give(31)ρ=3a3∫daa2V,whereas the pressure remains p=-V. The Friedmann equation can instead be manipulated to the form:(32)2H˙+3H2=Vt,where as usual the Hubble parameter is defined as H≡a˙/a. Further progress can be made by performing the substitution y≡a3/2, which yields the following equation [73]:(33)y¨-34Vty=0.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 y satisfies a linear differential equation. We now proceed to study a few interesting potentials and the corresponding solutions.


### 3.1.1.V∝1/ϕ2

Let us consider the following potential [73]:(34)Vϕ=αϕ2=αt2.Solving the corresponding (33) and substituting for the scale factor, a≡y2/3, yield the following solution for α≥-1/3:(35)at=t1/31+1+3α1+βt-1+3α2/3,where β is an integration constant. For α<-1/3 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) w≡p/ρ if we recall that the energy density is given by (31) whereas the pressure reads p=-V. Explicit calculation gives the following:(36)w=-3α1+1+3α1-βt-1+3α1+βt-1+3α-2.It is interesting to note that, for α≫1, the EoS approaches w=-1, 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 EoSw~. The scale factor in this situation evolves as(37)a∝t2/31+w~,and hence (31) can be used to deduce the evolution of the energy density of mimetic matter, which decays as(38)ρ=-αw~t2.Given that the pressure of mimetic matter reads p=-V, we immediately see that the EoS for mimetic matter is w=w~, 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 α/w≪1. 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:(39)Vϕ=αϕn=αtn,for which the solution of (33) can be written in terms of the Bessel functions [73]:(40)y=t1/2Z1/n+2-3αn+2tn+2/2.For n<-2 (with n=-2 corresponding to the case we have studied previously) the limiting behaviour of the scale factor is that of a dust-dominated universe, with EoS w=0. For n>-2 and α<0 (which corresponds to a positive pressure), the corresponding solution is a singular oscillating universe. For n>-2 and α>0 instead, the pressure is negative and hence we expect accelerating solutions [73]. In fact, n=0 corresponds to a cosmological constant as expected (the potential is simply a constant), whereas n=2 gives an inflationary expansion solution with scale factor:(41)a∝t-1/3eα/12t2,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 parameterH or equivalently in the scale factor a(t)], find the corresponding y-parameter through a=y2/3, and then invert (33) to find the corresponding potential which can provide the desired expansion:(42)Vϕ=Vt=4y¨3y.As an example, we can consider the following potential [73]:(43)Vϕ=αϕ2eϕ+1,whose corresponding solution for the scale factor is exponential (a∝e-t2) for large negative times and corresponding to EoS w=0 (i.e., dust, a∝t2/3) 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]:(44)Vϕ=αe-κϕ=αe-κt.In this case, the scale factor can be expressed in terms of Bessel functions of the first and second kinds [83]:(45)at=βJ0-3ακe-kt/2+γY0-3ακe-kt/22/3,where β and γ are integration constants and J0 and Y0 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 w=0 (a∝t2/3). At early times the behaviour of the scale factor depends on the sign of α, in particular providing us with an inflationary solution for α>0, whereas the solution is a bouncing nonsingular universe for α<0. A similar behaviour can be obtained if one chooses the potential:(46)Vϕ=αϕ2neκϕ+1=αt2neκt+1,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:(47)Vϕ=4311+ϕ22=4311+t22.As usual, the scale factor can be determined by solving (33), which yields [73](48)at=t2+11+βarctant2/3.If we set the integration constant β to 0, the corresponding energy density and pressure (one again, refer to (31) and recall that p=-V) are given by(49)ρ=43t21+t22,p=-4311+t22.At very early times (large negative t) the EoS approaches w→0; the universe is dominated by dust and contracts. At a certain time corresponding to |t|~1, 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 |β|<2/π [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](50)Vϕ=43αϕ02+ϕ22=43αt02+t22.Although we will not show the solution explicitly, in this case the scale of the bounce is reduced to α/t02 and the duration of the bounce is now t0 [73].


## 3.1.1.V∝1/ϕ2

Let us consider the following potential [73]:(34)Vϕ=αϕ2=αt2.Solving the corresponding (33) and substituting for the scale factor, a≡y2/3, yield the following solution for α≥-1/3:(35)at=t1/31+1+3α1+βt-1+3α2/3,where β is an integration constant. For α<-1/3 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) w≡p/ρ if we recall that the energy density is given by (31) whereas the pressure reads p=-V. Explicit calculation gives the following:(36)w=-3α1+1+3α1-βt-1+3α1+βt-1+3α-2.It is interesting to note that, for α≫1, the EoS approaches w=-1, 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 EoSw~. The scale factor in this situation evolves as(37)a∝t2/31+w~,and hence (31) can be used to deduce the evolution of the energy density of mimetic matter, which decays as(38)ρ=-αw~t2.Given that the pressure of mimetic matter reads p=-V, we immediately see that the EoS for mimetic matter is w=w~, 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 α/w≪1. 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:(39)Vϕ=αϕn=αtn,for which the solution of (33) can be written in terms of the Bessel functions [73]:(40)y=t1/2Z1/n+2-3αn+2tn+2/2.For n<-2 (with n=-2 corresponding to the case we have studied previously) the limiting behaviour of the scale factor is that of a dust-dominated universe, with EoS w=0. For n>-2 and α<0 (which corresponds to a positive pressure), the corresponding solution is a singular oscillating universe. For n>-2 and α>0 instead, the pressure is negative and hence we expect accelerating solutions [73]. In fact, n=0 corresponds to a cosmological constant as expected (the potential is simply a constant), whereas n=2 gives an inflationary expansion solution with scale factor:(41)a∝t-1/3eα/12t2,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 parameterH or equivalently in the scale factor a(t)], find the corresponding y-parameter through a=y2/3, and then invert (33) to find the corresponding potential which can provide the desired expansion:(42)Vϕ=Vt=4y¨3y.As an example, we can consider the following potential [73]:(43)Vϕ=αϕ2eϕ+1,whose corresponding solution for the scale factor is exponential (a∝e-t2) for large negative times and corresponding to EoS w=0 (i.e., dust, a∝t2/3) 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]:(44)Vϕ=αe-κϕ=αe-κt.In this case, the scale factor can be expressed in terms of Bessel functions of the first and second kinds [83]:(45)at=βJ0-3ακe-kt/2+γY0-3ακe-kt/22/3,where β and γ are integration constants and J0 and Y0 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 w=0 (a∝t2/3). At early times the behaviour of the scale factor depends on the sign of α, in particular providing us with an inflationary solution for α>0, whereas the solution is a bouncing nonsingular universe for α<0. A similar behaviour can be obtained if one chooses the potential:(46)Vϕ=αϕ2neκϕ+1=αt2neκt+1,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:(47)Vϕ=4311+ϕ22=4311+t22.As usual, the scale factor can be determined by solving (33), which yields [73](48)at=t2+11+βarctant2/3.If we set the integration constant β to 0, the corresponding energy density and pressure (one again, refer to (31) and recall that p=-V) are given by(49)ρ=43t21+t22,p=-4311+t22.At very early times (large negative t) the EoS approaches w→0; the universe is dominated by dust and contracts. At a certain time corresponding to |t|~1, 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 |β|<2/π [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](50)Vϕ=43αϕ02+ϕ22=43αt02+t22.Although we will not show the solution explicitly, in this case the scale of the bounce is reduced to α/t02 and the duration of the bounce is now t0 [73].


## 3.2. MimeticF(R) Gravity

The next step which was performed by Nojiri and Odintsov is to generalize mimetic gravity to mimeticF(R) 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 F(R) term. The action of the theory is given by [81](51)I=∫d4x-gg~μν,ϕFRg~μν,ϕ+Lm,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 F(R) gravity can be equally written employing a Lagrange multiplier field analogously to (11) [81]:(52)I=∫d4x-gFRgμν-Vϕ+λgμν∂μϕ∂νϕ+1+Lm,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]:(53)0=12gμνFR-RμνF′R+∇μ∇νF′R-gμν□F′R+12gμν-Vϕ+λgαβ∂αϕ∂βϕ-λ∂μϕ∂νϕ+12Tμν.Variation with respect to the mimetic field instead yields the following equation [81]:(54)2∇μλ∂μϕ+dVdϕ=0.As usual, by construction, variation with respect to the Lagrange multiplier gives the mimetic constraint:(55)gμν∂μϕ∂νϕ=-1.As we have mentioned previously, in mimeticF(R) 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 F(R) 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 mimeticF(R) 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 p. Assuming that the mimetic field depends only on time, (53), (54), and (55) can be expressed as follows [81]:(56)0=-FR+6H˙+H2F′R-6HdF′Rdt-λϕ˙2+1+Vϕ+ρ,(57)0=FR-2H˙+3H+2d2F′Rdt2+4HdF′Rdt-λϕ˙2-1-Vϕ+p,(58)0=2ddtλϕ˙+6Hλϕ˙-dVdϕ,(59)0=ϕ˙2-1.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:(60)0=FR-2H˙+3H+2d2F′Rdt2+4HdF′Rdt-Vt+p,which, if we assume that the contribution of ordinary matter is negligible (ρ=p=0), reduces to(61)Vt=FR-2H˙+3H+2d2F′Rdt2+4HdF′Rdt.On the other hand, (57) can be solved for λ as follows:(62)λt=-12FR+3H˙+H2F′R-3HdF′Rdt,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 F(R) which gives the wanted evolution [encoded in H(t) or, equivalently, a(t)]. Alternatively, one can start from a given form of F(R) 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 F(R) gravity is also a solution in mimetic F(R) 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 F(R) gravity. Of course, the exact forms of the mimetic potential or the F(R) 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 F(R) 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, ns and r, respectively. Both the reconstruction (determination of the potential once the form of F(R) and the evolutionary history of the universe are given) and the inverse reconstruction (determination of the form of F(R) 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 F(R) gravity in this respect: namely, the forms of both the mimetic potential and the F(R) function can become extremely complicated. The forms of both the mimetic potential and the Lagrange multiplier increase in complexity as the complexity of the F(R) form increases. Therefore, these and other studies on mimetic F(R) 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 F(R) 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 F(R) 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 F(R), 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 F(R) gravity with the desired inflationary properties would be to choose the simplest possible functional form of F(R) 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 F(R,ϕ) 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 mimeticF(R) 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 F(R) gravity as well. In other words, stable solutions in mimetic F(R) 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 F(R) gravity. An immediate implication of this finding is that, although mimetic F(R) gravity could drive inflation differently from F(R) gravity, the late-time acceleration of the universe in these theories has to coincide with the usual F(R) 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 F(R) gravity compared to conventional F(R) 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 F(R) gravity.As we mentioned above, the conclusions reached about the late-time evolution in mimeticF(R) gravity hold only at the background level. However, in conventional F(R) 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, dF(R)/dR) 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 F(R) gravity the problem is solved by adding power-law modifications by hand.In [153], it was instead argued that in mimetic F(R) 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 F(R) model whose corresponding conventional F(R) model suffered from the oscillations problems. The oscillations die out for redshifts z≤3, 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 F(R) 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 F(R) gravity as a cosmological framework.


### 3.2.2. MimeticF(R,ϕ) Gravity

The question of constructing a theoretically motivated but at the same time simple model for mimeticF(R) 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 F(R,ϕ) 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 F(R) 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 F(R).


### 3.2.3. Nonlocal MimeticF(R) Gravity

A further extension of mimeticF(R) 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 F(R) gravity is given by the action [133]:(63)I=∫d4x-gR1+f□-1R+λgμν∂μϕ∂νϕ+1-Vϕ,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]:(64)I=∫d4x-gR1+fψ+ξ□ψ-R+λgμν∂μϕ∂νϕ+1-Vϕ,where ξ is an additional Lagrange multiplier which enforces the constraint on the scalar field ψ:(65)□ψ=R.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]:(66)Rμν1+fψ-ξ-12gμνR1+fψ-ξ-∂ρξ∂ρψ=12∂μξ∂νψ+∂μψ∂νξ-gμν□-∇μ∇νfψ-ξ-λ∂μϕ∂νϕ-gμνVϕ.Instead, variation with respect to the two scalar fields leads to the following equations of motion [133]:(67)□ξ+f′ψR=0,-1-g∂ν-g∂νϕ=12dVdϕ.In [133], this model was studied in detail making use of the reconstruction technique. In particular, two forms for the f(ψ) 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.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 mimeticF(R) 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 F(R) gravity as well. In other words, stable solutions in mimetic F(R) 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 F(R) gravity. An immediate implication of this finding is that, although mimetic F(R) gravity could drive inflation differently from F(R) gravity, the late-time acceleration of the universe in these theories has to coincide with the usual F(R) 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 F(R) gravity compared to conventional F(R) 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 F(R) gravity.As we mentioned above, the conclusions reached about the late-time evolution in mimeticF(R) gravity hold only at the background level. However, in conventional F(R) 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, dF(R)/dR) 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 F(R) gravity the problem is solved by adding power-law modifications by hand.In [153], it was instead argued that in mimetic F(R) 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 F(R) model whose corresponding conventional F(R) model suffered from the oscillations problems. The oscillations die out for redshifts z≤3, 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 F(R) 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 F(R) gravity as a cosmological framework.


## 3.2.2. MimeticF(R,ϕ) Gravity

The question of constructing a theoretically motivated but at the same time simple model for mimeticF(R) 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 F(R,ϕ) 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 F(R) 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 F(R).


## 3.2.3. Nonlocal MimeticF(R) Gravity

A further extension of mimeticF(R) 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 F(R) gravity is given by the action [133]:(63)I=∫d4x-gR1+f□-1R+λgμν∂μϕ∂νϕ+1-Vϕ,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]:(64)I=∫d4x-gR1+fψ+ξ□ψ-R+λgμν∂μϕ∂νϕ+1-Vϕ,where ξ is an additional Lagrange multiplier which enforces the constraint on the scalar field ψ:(65)□ψ=R.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]:(66)Rμν1+fψ-ξ-12gμνR1+fψ-ξ-∂ρξ∂ρψ=12∂μξ∂νψ+∂μψ∂νξ-gμν□-∇μ∇νfψ-ξ-λ∂μϕ∂νϕ-gμνVϕ.Instead, variation with respect to the two scalar fields leads to the following equations of motion [133]:(67)□ξ+f′ψR=0,-1-g∂ν-g∂νϕ=12dVdϕ.In [133], this model was studied in detail making use of the reconstruction technique. In particular, two forms for the f(ψ) 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:(68)-g=1.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](69)I=∫d4x-gR-Vϕ-ηgμν∂μϕ∂νϕ+1-λ+λ,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]:(70)0=12gμνR-Vϕ-ηgαβ∂αϕ∂βϕ+1-λ-Rμν+η∂μϕ∂νϕ+12Tμν,whereas variation with respect to the mimetic field yields the usual equation of motion:(71)0=2∇μλ∂μϕ-dVdϕ.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 V(ϕ)=V(t) (as usual, the mimetic field can be identified with time):(72)Vϕ=Vt=at23∫0ϕdtat-2-6Htpt-2dptdt+2-18Ht3-6HtdHtdt+4d2Htdt2.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 H(t) or a(t). 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]:(73)Vϕ=12e2H0ϕH03ϕ.It can be easily shown that it leads to the following solution for the scale factor a(t) and the Hubble parameter H(t):(74)at=eH0t,Ht=H0,that is, a de Sitter cosmology. The functional forms of the Lagrange multiplier are also quite simple and are given by the following:(75)λt=6H021+2e2H0tH0t,ηt=3H02.Another choice for the potential which leads to a physically interesting solution is the following [136]:(76)Vϕ=-8ϕ-2+4/31+w1+5w+2w291+w3,for which the scale factor and the Hubble parameter read(77)at=t2/31+w,Ht=23t1+w.The above solution is that corresponding to an universe dominated by a fluid with EoS w. The two Lagrange multipliers are given by(78)λt=8-3w1+w+t4/31+w1+5w+2w29t21+w3,ηt=42+w3t21+w2.Thus, 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 p as follows:(79)ρ=G-T-4V~,p=-V~,where V~ is defined as(80)V~=-λt-Vt.Furthermore, 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 -g). It is nonetheless possible to present a covariant formulation of unimodular mimetic gravity via the following action:(81)I=∫d4x-gR-ηgμν∂μϕ∂νϕ+1-λ+λϵμνρσ∂μaνρσ,where aνρσ is a three-form field, variation of which gives the constraint ∂μλ=0, 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 with(82)-g=ϵμνρσ∂μaνρσ.Further 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 MimeticF(R) Gravity

A minimal extension of the unimodular mimetic gravity framework we have discussed so far is to consider unimodular mimeticF(R) gravity, which is described by the action [138]:(83)I=∫d4x-gFR-Vϕ-ηgμν∂μϕ∂νϕ+1-λ+λ.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]:(84)gμν2FR-Vϕ+ηgαβ∂αϕ∂βϕ+1-λ-RμνF′R-η∂μϕ∂νϕ+∇μ∇νF′R-gμν□F′R=0.Finally, variation with respect to the mimetic field yields the following equation:(85)0=2∇μλ∂μϕ-dVdϕ.The more complex structure of the equations of motion complicates the reconstruction procedure, that is, the equivalent of (72), which now reads [138](86)Vt=at3/22∫dtat-3/2ft,where the function f(t) is given by(87)ft=-18HtH˙+H2F′R-6H2dF′Rdt-18H˙H3+6Hd3F′Rdt3+6H¨+2H˙HF′R+HH˙dF′Rdt+6H2dF′Rdt-2Hd2F′Rdt2-2H¨+6H˙H+2d2F′Rdt2.Despite 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 F(R). 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.3.1. Unimodular MimeticF(R) Gravity

A minimal extension of the unimodular mimetic gravity framework we have discussed so far is to consider unimodular mimeticF(R) gravity, which is described by the action [138]:(83)I=∫d4x-gFR-Vϕ-ηgμν∂μϕ∂νϕ+1-λ+λ.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]:(84)gμν2FR-Vϕ+ηgαβ∂αϕ∂βϕ+1-λ-RμνF′R-η∂μϕ∂νϕ+∇μ∇νF′R-gμν□F′R=0.Finally, variation with respect to the mimetic field yields the following equation:(85)0=2∇μλ∂μϕ-dVdϕ.The more complex structure of the equations of motion complicates the reconstruction procedure, that is, the equivalent of (72), which now reads [138](86)Vt=at3/22∫dtat-3/2ft,where the function f(t) is given by(87)ft=-18HtH˙+H2F′R-6H2dF′Rdt-18H˙H3+6Hd3F′Rdt3+6H¨+2H˙HF′R+HH˙dF′Rdt+6H2dF′Rdt-2Hd2F′Rdt2-2H¨+6H˙H+2d2F′Rdt2.Despite 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 F(R). 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]:(88)I=∫d4x-gLgμν,∂λ1gμν,…,∂λ1⋯∂λpgμν,ϕ,∂λ1ϕ,…,∂λ1⋯∂λqϕ+λgμν∂μϕ∂νϕ+1,where p,q≥2 are integers and L 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 b(ϕ)gμν∂μϕ∂νϕ=-1, but for the sake of simplicity we will set b(ϕ)=1 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 b(ϕ)≠1 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:(89)I=∫d4x-gLH=∫d4x-g∑nLn,where Lns read(90)L0=KX,ϕ,L1=-G3X,ϕ□ϕ,L2=G4,XX,ϕ□ϕ2-∇μ∇νϕ2+RG4X,ϕ,L3=-16G5,XX,ϕ□ϕ3-3□ϕ∇μ∇νϕ2+2∇μ∇νϕ3+Gμν∇μ∇νϕG5X,ϕ,where X≡-gμν∇μϕ∇νϕ/2, (∇μ∇νϕ)≡∇μ∇νϕ∇μ∇νϕ, (∇μ∇νϕ)3≡∇μ∇νϕ∇μ∇ρϕ∇ν∇ρϕ, the functions K(X,ϕ), G3(X,ϕ), G4(X,ϕ), G5(X,ϕ) are free, and X denotes differentiation with respect to X.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 K, G3, G4, and G5, as well as in the function b(ϕ) (which, when ≠1, 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:(91)KX,ϕ=c2X,G3X,ϕ=0,G4X,ϕ=12,G5X,ϕ=0.In this case, on a flat FLRW background the solution is given by(92)at=t2/31+w,ϕt=±-αc2lntt0,bϕ=-1ϕ˙2,with t0 and integration constant and α=-8w/3(1+w)2. Thus, we see that the scenario under consideration has reproduced the expansion history of a universe filled with a perfect fluid with EoS w [104]. Another case considered in [104] is the mimetic cubic Galileon model, where the functions take the following form:(93)KX,ϕ=c2X,G3X,ϕ=2c3Λ~3X,G4X,ϕ=12,G5X,ϕ=0,where the cut-off scale is subsequently set to Λ~=1. 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, g~μν 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 by(94)IH=∫d4x-gαXR+□ϕ2-∇μ∇νϕ∇μ∇νϕ+γϕGμν∇μ∇νϕ-βϕ□ϕ,which corresponds to the following choice for the functions discussed previously:(95)KX,ϕ=0,G3X,ϕ=βϕ,G4X,ϕ=αX,G5X,ϕ=0.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],(96)Vϕ=-β+2V0+3V02c1ϕ-ϕ02,gives rise to the following solution:(97)Ht=V0c1t-t0,which represents a regular bounce solution. Another bounce solution can be obtained by considering the following potential [114]:(98)Vϕ=b2c1sinh2bϕcosh2bϕ,for which the scale factor and correspondingly the Hubble parameter read:(99)at=a0coshbt,Ht=bsinhbtcoshbt.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 uμ (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, N=N(t), the theory takes the name of projectable Hořava-Lifshitz gravity. Recall that in the Arnowitt-Deser-Misner decomposition of space-time [445]:(100)ds2=-N2c2dt2+gijdxi+Nidtdxj+Njdt,the function N 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,uμ=N∇μT, where T is a scalar field and N is chosen in such a way to ensure uμ has unit norm. In this case T is responsible for the preferred foliation in Hořava gravity, and N is the lapse function. If T is set equal to coordinate time, the resulting action is that of nonprojectable Hořava gravity [357, 443].Further requiring thatNdT=S, where S is a scalar field, reduces the vector Einstein-aether theory to the scalar Einstein-aether theory [446]. This condition implies that N=N(T), which upon identification of T 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 N but has to be imposed at the level of the action, for instance, via a Lagrange multiplier term ∝λ(∇μS∇μS-1). Therefore, the condition under which dS=NdT is invariant is that S be invariant under a shift symmetry: S→S+dS, 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 as(101)S=SEH+∫d4x-gΣ2gμν∂μϕ∂νϕ+1+γ2□ϕ2,which 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](102)I=∫d4x-gR2-α∂μϕ∂νϕ∇μ∇ν+2U0∇ρ∇ρn∂μϕ∂νϕRμν+U0R2-λ12∂μϕ∂μϕ+U0,where U0 and α are constants and U0=1/2 to recover the mimetic constraint. In particular, the above case corresponds to the z=2n+2 version of CRG, where z 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]:(103)I=12∫d4x-gR-2Λ-αRμν-R2gμν∇μϕ∇νϕn-λ2gμν∇μϕ∇νϕ+1-Vϕ.In the above action, the Horndeski-like coupling of the scalar field to curvature (for the case n=1) 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.


## 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 uμ (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, N=N(t), the theory takes the name of projectable Hořava-Lifshitz gravity. Recall that in the Arnowitt-Deser-Misner decomposition of space-time [445]:(100)ds2=-N2c2dt2+gijdxi+Nidtdxj+Njdt,the function N 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,uμ=N∇μT, where T is a scalar field and N is chosen in such a way to ensure uμ has unit norm. In this case T is responsible for the preferred foliation in Hořava gravity, and N is the lapse function. If T is set equal to coordinate time, the resulting action is that of nonprojectable Hořava gravity [357, 443].Further requiring thatNdT=S, where S is a scalar field, reduces the vector Einstein-aether theory to the scalar Einstein-aether theory [446]. This condition implies that N=N(T), which upon identification of T 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 N but has to be imposed at the level of the action, for instance, via a Lagrange multiplier term ∝λ(∇μS∇μS-1). Therefore, the condition under which dS=NdT is invariant is that S be invariant under a shift symmetry: S→S+dS, 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 as(101)S=SEH+∫d4x-gΣ2gμν∂μϕ∂νϕ+1+γ2□ϕ2,which 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](102)I=∫d4x-gR2-α∂μϕ∂νϕ∇μ∇ν+2U0∇ρ∇ρn∂μϕ∂νϕRμν+U0R2-λ12∂μϕ∂μϕ+U0,where U0 and α are constants and U0=1/2 to recover the mimetic constraint. In particular, the above case corresponds to the z=2n+2 version of CRG, where z 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]:(103)I=12∫d4x-gR-2Λ-αRμν-R2gμν∇μϕ∇νϕn-λ2gμν∇μϕ∇νϕ+1-Vϕ.In the above action, the Horndeski-like coupling of the scalar field to curvature (for the case n=1) 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 Section2 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 satisfycs=0. 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(104)ds2=-1+2Φt,xdt2+a2t1-2Ψt,xδijdxidxj,i,j=1,2,3.Φ≡Φ(t,x) and Ψ≡Ψ(t,x) are functions of the space-time coordinates such that |Φ(t,x)|,|Ψ(t,x)|≪1, and g00(t,x)≃-1+2Φ(t,x), g11(t,x)≃a(t)-2(1+2Ψ(t,x)). Here, we used the conformal Newtonian gauge. Moreover, from field equations we also have(105)Φt,x=Ψt,x.Correspondingly, given that around a flat FLRW background the mimetic field plays the role of “clock,” we perturb it as(106)ϕ=t+δϕt,x,where |δϕ|≡|δϕ(t,x)|≪1, which together with (3) implies that(107)Φ=δϕ˙.Then, by perturbing the 0i 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:(108)δϕ¨+Hδϕ˙+H˙δϕ=0.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 cs2Δδϕ 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]:(109)I=d4x-gR+λgμν∂μϕ∂νϕ+1-Vϕ+12γ□ϕ2.The corresponding equations of motion read(110)Gνμ=V+γ∂αϕ∂αχδνμ+2λ∂μϕ∂νϕ-γ∂μχ∂νϕ+∂μϕ∂νχ,where χ≡□ϕ. Thus, identification of the mimetic field with time on a flat FLRW background implies that χ=3H. 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](111)2H˙+3H2=22-3γV,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]:(112)δϕ¨+Hδϕ˙-cs2a2Δδϕ+H˙δϕ=0,where(113)cs2=γ2-3γ.Therefore, 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., ns<1) scalar perturbations which are enhanced over gravity waves (implying a small value of r), 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 (□ϕ)2 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 betweenN-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:(114)I=12∫d4x-gR-2Λ-αRμν-β2Rgμν∇μϕ∇νϕn-λ12gμν∂μϕ∂νϕ+U0,where α, β are arbitrary constants, Λ is the cosmological constant, n is a natural number, λ is a Lagrange multiplier, and U0 determines the constraint imposed on the gradient of the cosmological field ϕ. Recall that the Horndeski-like coupling in the action above, for n=1, has been considered in several works (e.g., [420, 433, 447–455]). If we set U0=1/2, 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 reads(115)I=12∫d4x-gR-2Λ-αRμν-R2gμν∇μϕ∇νϕn-λ2gμν∂μϕ∂νϕ+1-Vϕ,where we have set the dimensionless parameter β=1 and added a potential for the mimetic field, V(ϕ). Variation with respect to λ immediately leads to (3), whereas variation with respect to the field modifies (7) to(116)dVϕdϕ=∇μ2nαFn-1Gμν+λgμν∂νϕ=1-g∂μ-g2nαFn-1Gμν+λgμν∂νϕ,where we define the following quantities:(117)F≡TμνRμν-RT2,Tμν≡∇μϕ∇νϕ,T≡gμνTμν=-1.Finally, variation of the action with respect to the metric leads to the gravitational field equations, which for this theory read(118)Gμν+Λgμν+α2Fngμν=nαFn-1RμρTρν+RνρTρμ-12TRμν+RTμν+λ2Tμν+nαDαβμνTαβFn-1-12DμνTFn-1+ΩαβδTαβδgμν-gμνVϕ2,where we have defined the differential operators:(119)Dαβμν≡14gμαgνβ+gναgμβ□+gμν∇α∇β+∇β∇α-gμα∇β∇ν+gνα∇β∇μ+gμβ∇α∇ν+gνβ∇α∇μ,Dμν≡gμν□-12∇μ∇ν+∇ν∇μ,where □≡∇i∇i is the d’Alembertian operator. Note that in the above (118), Ωμν is a tensor that will not play any role if Tμν 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):(120)-R+4Λ-λ2T=2αFnn-1+nα2gμν□+∇μ∇ν+∇ν∇μTμνFn-1-3nα2□TFn-1-2Vϕ.We will consider the case with n=1, 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:(121)ϕ=t.As in the original mimetic gravity, the scalar field can induce an effective cold dark matter component given that, from (7), one has -(G-T)∝a3, and recall that -(G-T) 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 tensorTμν in (117), which reads(122)T00=ϕ˙2=1,T0i=Ti0=Tij=0,i,j=1,2,3.The, from (116) we derive(123)1a3∂0a32nα3H2n-λ=dVϕdϕ.Given the mimetic constraint (3), this equation is a consequence of the two EOMs inferred from (118):(124)0=Λ-3H2+α21-4n3H2n+λ2+Vϕ2,(125)0=Λ-3H2-2H˙+α21-2n3H2n+3n-1αn1-2nH˙H2n-2+Vϕ2.In our analysis, as the second independent equation, together, with (123), we choose the trace equation (120): namely,(126)λ2=6H˙+12H2-4Λ+α5n-23H2n+3nnα2n-1H2n-2H˙-2Vϕ.For the simplest case V(ϕ)=0, from (123) we get(127)2nα3H2n-λ=C0a3,where C0 is an integration constant. Thus, in the limit α=0, we recover the original mimetic gravity model of [69], while when α≠0 we may interpret this equation as a generalized Friedmann equation, with C0 determining the amount of dark matter in the universe (given that its contribution scales with scale factor a as a-3, as is expected for dust). Moreover, in the limits given by(128)1≪α,n=1,1αn-1≪R,1<n,if we neglect the cosmological constant Λ, as well as the integration constant C0, one gets from (126)(129)-H˙H2=2n.Thus, if 1/αn-1 is set at the scale of the early-time acceleration, models with 1≪n lead to cosmological solutions for inflation.In the case where we setn=1, 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) that(130)dVdt=1a3∂0a36αH2-λ,(131)λ=62+αH˙+64+3αH2-4Vϕ,where we have taken into account that V≡V(ϕ)=V(t) 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 gives(132)2H˙+3H2=V2+α,showing 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:(133)u¨-342+αVu=0,by introducing the Sturm-Liouville canonical substitution(134)H=23u˙u,a=u2/3.Let us discuss some examples by starting from (133). If V=V0 is constant, we recover the de Sitter solution:(135)u~exp3H0t2,a~expH0t,H0=233V042+α.Another well-motivated choice for the potential is a quadratic one:(136)Vϕ=32+αH02+β22ϕ-ϕ0-H0+β242ϕ-ϕ0-23β2,where H0 is a constant Hubble parameter and β, ϕ0 are dimensional constants (in the specific case [β]=[ϕ0-1]=[H]). After the identification ϕ=t the explicit solution for u is found to be(137)ut=u0e3/2H0t-3β2/4tt-2t0,t0≡ϕ02,with u0 constant and t0 a fixed time. The Hubble parameter is given by(138)H≡23u˙u=H0-β2t-t0,and we see that, for t close to t0, one has a quasi-de Sitter expansion, while for large t0≪t, 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 by(139)Vϕ=4A22+α3coshAϕ1+coshAϕ,where 0<A is a constant with mass-dimension 1. The corresponding solution is given by(140)ut=1+coshAt,a=1+coshAt2/3,H=2A3sinhAt1+coshAt.This solution represents a cosmological bounce with -∞<t,ϕ<+∞ 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 n=1, 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 that(141)T00=1+2δϕ˙,T0i=∂iδϕ,T=-1+OΦ2.From the (1,2)-component of (118), we obtain(142)G121-α2=αDαβ12Tαβ,with(143)G12=-∂x∂yΦ-Ψ,Dαβ12Tαβ=H∂x∂yδϕ+∂x∂yδϕ˙.Therefore, we obtain(144)Ψ=Φ+2α2-αHδϕ+δϕ˙.From (0,1)-component of (118) we derive(145)G011+α2=αH˙∂xδϕ+λ2∂xδϕ+αDαβ01Tαβ,with(146)G01=2∂xΨ˙+HΦ,Dαβ01Tαβ=-H2+H˙∂xδϕ,λ=6αH2-4H˙-2αH˙,where the last equality is a consequence of (131), (132).Finally, we can obtain a closed equation forδϕ, which reads(147)δϕ¨+Hδϕ˙+H˙δϕ=0.From 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 n≠1 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 n=1 as(148)I=12∫d4x-gR1+agμν∇μϕ∇νϕ-c2□ϕ2+b2∇μ∇νϕ2-λ2gμν∇μϕ∇νϕ+1-Vϕ.The original model given by (103) for n=1 is recovered for a=α/2 and b=c=2α by using the following identity [508]:(149)-12gμν∇μϕ∇νϕR+□ϕ2-∇μ∇νϕ2=Gμν∇μϕ∇νϕ+totalderivative. On the other hand, for generic values of a, b, c, the action (148) describes a higher order derivative model in the scalar sector with field equations at the fourth order; namely,(150)1-aGμν=12gμνb2ϕαβϕαβ-c2□ϕ2-Vϕ+λ∇μϕ∇νϕ-bϕμρϕνρ+b2gαβ∇αϕμν∇βϕ-∇αϕμβ∇νϕ-∇αϕνβ∇μϕ+cϕμν+gμνgαβ∇α□ϕ∇βϕ-∇μ□ϕ∇νϕ-∇ν□ϕ∇μϕ,where we adopted the notation(151)ϕαβ≡∇α∇βϕ.Let us explore some cosmological applications. On a flat FRW metric (10), (132) is modified as follows:(152)2H˙+chH2=cvVϕ,where(153)ch≡12a+3b-9c-124a+b-3c-4,cv≡24-4a-b+3c.It is understood that by setting a=α/2 and b=c=2α, we recover (132). Similarly to how we proceeded in the previous section, we can rewrite (152) as(154)u¨-chcvVϕ4u=0,where we introduced the auxiliary function u:(155)H=2chu˙u,at=u2/ch.By choosing the following potential:(156)Vϕ=9chcvH02+β22ϕ-ϕ0-H0+β242ϕ-ϕ0-23β2, with H0, β, ϕ0 dimensional constants, we recover the Starobinsky-like solution (137) and (138). If we choose the following potential:(157)Vϕ=4A2cvchcoshAϕ1+coshAϕ,with 0<A 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 (i,j)-components, i,j=1,2,3 of (150), we obtain now(158)Ψ=Φ+b2-2aδϕ˙+Hδϕ.Moreover, from the components (0,i) or (i,0), i=1,2,3 of (150), we derive(159)δϕ¨+Hδϕ˙-cs2a2∇2δϕ+H˙δϕ=0,where the nonvanishing squared sound speed cs2 reads(160)cs2≡b-c2c2,c2≡2+b-2a4+3c-4a-b4a-1.We immediately note that cs2=0 when b=c and we recover (147). Since only for b=c does the model fall within the Horndeski class, we can consider (b-c) 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, b≠c. The result is analogous to that obtained in [73, 84].


## 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:(114)I=12∫d4x-gR-2Λ-αRμν-β2Rgμν∇μϕ∇νϕn-λ12gμν∂μϕ∂νϕ+U0,where α, β are arbitrary constants, Λ is the cosmological constant, n is a natural number, λ is a Lagrange multiplier, and U0 determines the constraint imposed on the gradient of the cosmological field ϕ. Recall that the Horndeski-like coupling in the action above, for n=1, has been considered in several works (e.g., [420, 433, 447–455]). If we set U0=1/2, 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 reads(115)I=12∫d4x-gR-2Λ-αRμν-R2gμν∇μϕ∇νϕn-λ2gμν∂μϕ∂νϕ+1-Vϕ,where we have set the dimensionless parameter β=1 and added a potential for the mimetic field, V(ϕ). Variation with respect to λ immediately leads to (3), whereas variation with respect to the field modifies (7) to(116)dVϕdϕ=∇μ2nαFn-1Gμν+λgμν∂νϕ=1-g∂μ-g2nαFn-1Gμν+λgμν∂νϕ,where we define the following quantities:(117)F≡TμνRμν-RT2,Tμν≡∇μϕ∇νϕ,T≡gμνTμν=-1.Finally, variation of the action with respect to the metric leads to the gravitational field equations, which for this theory read(118)Gμν+Λgμν+α2Fngμν=nαFn-1RμρTρν+RνρTρμ-12TRμν+RTμν+λ2Tμν+nαDαβμνTαβFn-1-12DμνTFn-1+ΩαβδTαβδgμν-gμνVϕ2,where we have defined the differential operators:(119)Dαβμν≡14gμαgνβ+gναgμβ□+gμν∇α∇β+∇β∇α-gμα∇β∇ν+gνα∇β∇μ+gμβ∇α∇ν+gνβ∇α∇μ,Dμν≡gμν□-12∇μ∇ν+∇ν∇μ,where □≡∇i∇i is the d’Alembertian operator. Note that in the above (118), Ωμν is a tensor that will not play any role if Tμν 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):(120)-R+4Λ-λ2T=2αFnn-1+nα2gμν□+∇μ∇ν+∇ν∇μTμνFn-1-3nα2□TFn-1-2Vϕ.We will consider the case with n=1, 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:(121)ϕ=t.As in the original mimetic gravity, the scalar field can induce an effective cold dark matter component given that, from (7), one has -(G-T)∝a3, and recall that -(G-T) 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 tensorTμν in (117), which reads(122)T00=ϕ˙2=1,T0i=Ti0=Tij=0,i,j=1,2,3.The, from (116) we derive(123)1a3∂0a32nα3H2n-λ=dVϕdϕ.Given the mimetic constraint (3), this equation is a consequence of the two EOMs inferred from (118):(124)0=Λ-3H2+α21-4n3H2n+λ2+Vϕ2,(125)0=Λ-3H2-2H˙+α21-2n3H2n+3n-1αn1-2nH˙H2n-2+Vϕ2.In our analysis, as the second independent equation, together, with (123), we choose the trace equation (120): namely,(126)λ2=6H˙+12H2-4Λ+α5n-23H2n+3nnα2n-1H2n-2H˙-2Vϕ.For the simplest case V(ϕ)=0, from (123) we get(127)2nα3H2n-λ=C0a3,where C0 is an integration constant. Thus, in the limit α=0, we recover the original mimetic gravity model of [69], while when α≠0 we may interpret this equation as a generalized Friedmann equation, with C0 determining the amount of dark matter in the universe (given that its contribution scales with scale factor a as a-3, as is expected for dust). Moreover, in the limits given by(128)1≪α,n=1,1αn-1≪R,1<n,if we neglect the cosmological constant Λ, as well as the integration constant C0, one gets from (126)(129)-H˙H2=2n.Thus, if 1/αn-1 is set at the scale of the early-time acceleration, models with 1≪n lead to cosmological solutions for inflation.In the case where we setn=1, 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) that(130)dVdt=1a3∂0a36αH2-λ,(131)λ=62+αH˙+64+3αH2-4Vϕ,where we have taken into account that V≡V(ϕ)=V(t) 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 gives(132)2H˙+3H2=V2+α,showing 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:(133)u¨-342+αVu=0,by introducing the Sturm-Liouville canonical substitution(134)H=23u˙u,a=u2/3.Let us discuss some examples by starting from (133). If V=V0 is constant, we recover the de Sitter solution:(135)u~exp3H0t2,a~expH0t,H0=233V042+α.Another well-motivated choice for the potential is a quadratic one:(136)Vϕ=32+αH02+β22ϕ-ϕ0-H0+β242ϕ-ϕ0-23β2,where H0 is a constant Hubble parameter and β, ϕ0 are dimensional constants (in the specific case [β]=[ϕ0-1]=[H]). After the identification ϕ=t the explicit solution for u is found to be(137)ut=u0e3/2H0t-3β2/4tt-2t0,t0≡ϕ02,with u0 constant and t0 a fixed time. The Hubble parameter is given by(138)H≡23u˙u=H0-β2t-t0,and we see that, for t close to t0, one has a quasi-de Sitter expansion, while for large t0≪t, 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 by(139)Vϕ=4A22+α3coshAϕ1+coshAϕ,where 0<A is a constant with mass-dimension 1. The corresponding solution is given by(140)ut=1+coshAt,a=1+coshAt2/3,H=2A3sinhAt1+coshAt.This solution represents a cosmological bounce with -∞<t,ϕ<+∞ 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 n=1, 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 that(141)T00=1+2δϕ˙,T0i=∂iδϕ,T=-1+OΦ2.From the (1,2)-component of (118), we obtain(142)G121-α2=αDαβ12Tαβ,with(143)G12=-∂x∂yΦ-Ψ,Dαβ12Tαβ=H∂x∂yδϕ+∂x∂yδϕ˙.Therefore, we obtain(144)Ψ=Φ+2α2-αHδϕ+δϕ˙.From (0,1)-component of (118) we derive(145)G011+α2=αH˙∂xδϕ+λ2∂xδϕ+αDαβ01Tαβ,with(146)G01=2∂xΨ˙+HΦ,Dαβ01Tαβ=-H2+H˙∂xδϕ,λ=6αH2-4H˙-2αH˙,where the last equality is a consequence of (131), (132).Finally, we can obtain a closed equation forδϕ, which reads(147)δϕ¨+Hδϕ˙+H˙δϕ=0.From 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 n≠1 in (103) or we must modify the original action along the lines of [73, 84].


## 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:(121)ϕ=t.As in the original mimetic gravity, the scalar field can induce an effective cold dark matter component given that, from (7), one has -(G-T)∝a3, and recall that -(G-T) 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 tensorTμν in (117), which reads(122)T00=ϕ˙2=1,T0i=Ti0=Tij=0,i,j=1,2,3.The, from (116) we derive(123)1a3∂0a32nα3H2n-λ=dVϕdϕ.Given the mimetic constraint (3), this equation is a consequence of the two EOMs inferred from (118):(124)0=Λ-3H2+α21-4n3H2n+λ2+Vϕ2,(125)0=Λ-3H2-2H˙+α21-2n3H2n+3n-1αn1-2nH˙H2n-2+Vϕ2.In our analysis, as the second independent equation, together, with (123), we choose the trace equation (120): namely,(126)λ2=6H˙+12H2-4Λ+α5n-23H2n+3nnα2n-1H2n-2H˙-2Vϕ.For the simplest case V(ϕ)=0, from (123) we get(127)2nα3H2n-λ=C0a3,where C0 is an integration constant. Thus, in the limit α=0, we recover the original mimetic gravity model of [69], while when α≠0 we may interpret this equation as a generalized Friedmann equation, with C0 determining the amount of dark matter in the universe (given that its contribution scales with scale factor a as a-3, as is expected for dust). Moreover, in the limits given by(128)1≪α,n=1,1αn-1≪R,1<n,if we neglect the cosmological constant Λ, as well as the integration constant C0, one gets from (126)(129)-H˙H2=2n.Thus, if 1/αn-1 is set at the scale of the early-time acceleration, models with 1≪n lead to cosmological solutions for inflation.In the case where we setn=1, 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) that(130)dVdt=1a3∂0a36αH2-λ,(131)λ=62+αH˙+64+3αH2-4Vϕ,where we have taken into account that V≡V(ϕ)=V(t) 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 gives(132)2H˙+3H2=V2+α,showing 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:(133)u¨-342+αVu=0,by introducing the Sturm-Liouville canonical substitution(134)H=23u˙u,a=u2/3.Let us discuss some examples by starting from (133). If V=V0 is constant, we recover the de Sitter solution:(135)u~exp3H0t2,a~expH0t,H0=233V042+α.Another well-motivated choice for the potential is a quadratic one:(136)Vϕ=32+αH02+β22ϕ-ϕ0-H0+β242ϕ-ϕ0-23β2,where H0 is a constant Hubble parameter and β, ϕ0 are dimensional constants (in the specific case [β]=[ϕ0-1]=[H]). After the identification ϕ=t the explicit solution for u is found to be(137)ut=u0e3/2H0t-3β2/4tt-2t0,t0≡ϕ02,with u0 constant and t0 a fixed time. The Hubble parameter is given by(138)H≡23u˙u=H0-β2t-t0,and we see that, for t close to t0, one has a quasi-de Sitter expansion, while for large t0≪t, 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 by(139)Vϕ=4A22+α3coshAϕ1+coshAϕ,where 0<A is a constant with mass-dimension 1. The corresponding solution is given by(140)ut=1+coshAt,a=1+coshAt2/3,H=2A3sinhAt1+coshAt.This solution represents a cosmological bounce with -∞<t,ϕ<+∞ 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 n=1, 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 that(141)T00=1+2δϕ˙,T0i=∂iδϕ,T=-1+OΦ2.From the (1,2)-component of (118), we obtain(142)G121-α2=αDαβ12Tαβ,with(143)G12=-∂x∂yΦ-Ψ,Dαβ12Tαβ=H∂x∂yδϕ+∂x∂yδϕ˙.Therefore, we obtain(144)Ψ=Φ+2α2-αHδϕ+δϕ˙.From (0,1)-component of (118) we derive(145)G011+α2=αH˙∂xδϕ+λ2∂xδϕ+αDαβ01Tαβ,with(146)G01=2∂xΨ˙+HΦ,Dαβ01Tαβ=-H2+H˙∂xδϕ,λ=6αH2-4H˙-2αH˙,where the last equality is a consequence of (131), (132).Finally, we can obtain a closed equation forδϕ, which reads(147)δϕ¨+Hδϕ˙+H˙δϕ=0.From 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 n≠1 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 n=1 as(148)I=12∫d4x-gR1+agμν∇μϕ∇νϕ-c2□ϕ2+b2∇μ∇νϕ2-λ2gμν∇μϕ∇νϕ+1-Vϕ.The original model given by (103) for n=1 is recovered for a=α/2 and b=c=2α by using the following identity [508]:(149)-12gμν∇μϕ∇νϕR+□ϕ2-∇μ∇νϕ2=Gμν∇μϕ∇νϕ+totalderivative. On the other hand, for generic values of a, b, c, the action (148) describes a higher order derivative model in the scalar sector with field equations at the fourth order; namely,(150)1-aGμν=12gμνb2ϕαβϕαβ-c2□ϕ2-Vϕ+λ∇μϕ∇νϕ-bϕμρϕνρ+b2gαβ∇αϕμν∇βϕ-∇αϕμβ∇νϕ-∇αϕνβ∇μϕ+cϕμν+gμνgαβ∇α□ϕ∇βϕ-∇μ□ϕ∇νϕ-∇ν□ϕ∇μϕ,where we adopted the notation(151)ϕαβ≡∇α∇βϕ.Let us explore some cosmological applications. On a flat FRW metric (10), (132) is modified as follows:(152)2H˙+chH2=cvVϕ,where(153)ch≡12a+3b-9c-124a+b-3c-4,cv≡24-4a-b+3c.It is understood that by setting a=α/2 and b=c=2α, we recover (132). Similarly to how we proceeded in the previous section, we can rewrite (152) as(154)u¨-chcvVϕ4u=0,where we introduced the auxiliary function u:(155)H=2chu˙u,at=u2/ch.By choosing the following potential:(156)Vϕ=9chcvH02+β22ϕ-ϕ0-H0+β242ϕ-ϕ0-23β2, with H0, β, ϕ0 dimensional constants, we recover the Starobinsky-like solution (137) and (138). If we choose the following potential:(157)Vϕ=4A2cvchcoshAϕ1+coshAϕ,with 0<A 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 (i,j)-components, i,j=1,2,3 of (150), we obtain now(158)Ψ=Φ+b2-2aδϕ˙+Hδϕ.Moreover, from the components (0,i) or (i,0), i=1,2,3 of (150), we derive(159)δϕ¨+Hδϕ˙-cs2a2∇2δϕ+H˙δϕ=0,where the nonvanishing squared sound speed cs2 reads(160)cs2≡b-c2c2,c2≡2+b-2a4+3c-4a-b4a-1.We immediately note that cs2=0 when b=c and we recover (147). Since only for b=c does the model fall within the Horndeski class, we can consider (b-c) 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, b≠c. The result is analogous to that obtained in [73, 84].


## 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 (i,j)-components, i,j=1,2,3 of (150), we obtain now(158)Ψ=Φ+b2-2aδϕ˙+Hδϕ.Moreover, from the components (0,i) or (i,0), i=1,2,3 of (150), we derive(159)δϕ¨+Hδϕ˙-cs2a2∇2δϕ+H˙δϕ=0,where the nonvanishing squared sound speed cs2 reads(160)cs2≡b-c2c2,c2≡2+b-2a4+3c-4a-b4a-1.We immediately note that cs2=0 when b=c and we recover (147). Since only for b=c does the model fall within the Horndeski class, we can consider (b-c) 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, b≠c. 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 to(161)Gμν=λ2gμν∂μϕ∂νϕ-gμνVϕ2,namely, (118) with α=0 and Λ=0. The trace of this equation reads(162)R=λ2+2Vϕ.Thus, from (161), we get(163)Gμν=R-2Vϕ∂μϕ∂νϕ-12gμνVϕ.The continuity equation of the mimetic field is derived as(164)1-g∂ν-gλ∂νϕ=dVϕdϕ,which is automatically satisfied when (163) holds true.In this chapter we will consider pseudo-SSS space-times, whose general topological formulation is given by(165)ds2=-ar2brdt2+dr2br+r2dρ21-kρ2+ρ2dϕ2,where a(r), b(r) are functions of the radial coordinate r and the manifold is a sphere when k=1, a torus when k=0, or a compact hyperbolic manifold when k=-1. The Ricci scalar in this case reads(166)R=-1r23r2b′ra′rar+r2b′′r+2r2bra′′rar+4rb′r+4rbra′rar+2br-2k,where the prime denotes a derivative with respect to r. The symmetries of the EOMs require the mimetic field to be a function of r only: namely, ϕ≡ϕ(r), and from (3) one has(167)ϕ′r=-1br,leading 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 uμ in (9) cannot be physical. This has been explored in [119] and will be discussed in Section 7.The(0,0)- and (1,1)-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], L=2a(1-rb′(r)-b(r))+a(r)r2λbrϕ′(r)2+1-a(r)r2V(ϕ). Thus, the derivatives with respect to a(r) and b(r) with (3) lead to (168) and (170), while the derivation with respect to ϕ leads to (171)):(168)k-b′rr-br=Vϕr22,(169)b′rr+2ra′rarbr+br-k=λ2brr2ϕ′r2-Vϕr22,where λ is given by (162).We can also rewrite (169) with (167)-(168) as(170)4a′rbr=-λarr.Finally, from (164), one finds the following:(171)ddrarbrλr2ϕ′=arr2dVϕdϕ.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 setV(ϕ)=0. From (168) we immediately see that(172)br=k-rsr,with rs being a mass scale, either positive or negative, while the second field equation (169) with (162) and (166) leads to(173)ar=a1+a21-krs/r1-krsrlogrr01+1-krsr-1,k=±1,(174)ar=a1+a22rr03/2+3,with a1, a2 being constants and r0 a length scale. If a2=0, namely, R=0 and λ=0, we recover the topological Schwarzschild solution of General Relativity. When a2≠0, we can pose a1=0. Thus, the solution for flat topology (k=0) is as follows:(175)ds2=-a222rr03/2+32r~srdt2+dr2r~s/r+r2dρ2+ρ2dϕ2,0<r~s,where r~s=-rs has to be positive to preserve the metric signature. The Ricci scalar is nonzero and reads(176)R=-6r~2r3+3r/r03/2r03.This solution presents a naked singularity at r=0, as for the corresponding topological case k=0 of the Schwarzschild metric.The spherical case (k=1) is the more interesting and is given by(177)ds2=-a221-rsrlogrr01+1-rsr-12dt2+dr21-rs/r+r2dθ2+sin2θdϕ2,where we have introduced the polar coordinates θ, ϕ in the angular part (see also [79]). In this case the Ricci scalar reads(178)R=-1r21-rs/rlogr/r01+1-rs/r-1.For rs<0 the metric is regular everywhere, but when 0<rs, the solution is regular and preserves the signature; only in the region rs<r the solution is regular and preserves the signature; since for r<rs the metric coefficient g00(r)≡-a2(r)b(r) acquires an imaginary part. In this respect, we note that the special choice rs=r0 when r<rs leads to(179)g00r≡ar2br=-a22rsr-1arctanrsr-1+12,which is real and negative.In general, we observe that when0<rs, the point r=rs cannot represent a horizon like for the Schwarzschild space-time, due to the fact that g11(r)≡1/b(r) diverges but g00(r)=-a22 and the Ricci scalar are regular (we also cannot associate with r=rs any thermodynamical quantity such as temperature T=(a(r)/2)db(r)/dr which diverges at r=rs); moreover, for r<rs the metric becomes imaginary (or, in the special case rs=r0, it acquires the signature (−+++)) and does not allow a wormhole description.Finally, the topological casek=-1 reads(180)ds2=-a221+rsrlogrr01+1+rsr-12dt2+dr2-1-rs/r+r2dρ21+ρ2+ρ2dϕ2,and the Ricci scalar is given by(181)R=-1r21+rs/rlogr/r01+1+rs/r-1.In this case, for 0<rs, the metric coefficient g11(r)≡1/b(r) is negative and the solution is unphysical. On the other hand, for rs<0, we obtain 0<g11(r) when r<-rs, but g00(r)≡-a2(r)b(r) becomes imaginary (except for the special choice r0=-rs).


### 6.2. Nonvacuum Solutions

In this subsection, we will consider the caseV(ϕ)≠0. We immediately see that, for V(ϕ)=2Λ with Λ a cosmological constant, one solution of (168) and (169) with Λ=0 is the topological Schwarzschild de Sitter solution:(182)br=k-rsr-Λr23,ar=a1,with rs, a1 being constants.In considering other solutions, we will take the spherical casek=1 in (165). By using (167) one readily obtains from (170) and (171) that(183)ϕr=±i∫drbr,4ddra′rbr3/2r=arr2brdVrdr,where we treat the potential as a function of r. These two equations with (168) can be used to reconstruct the potential when a choice for b(r) is made. We will now provide some examples of the reconstruction technique.Let us consider a linear modification to the Schwarzschild metric:(184)br=1-rsr+γr,with rs, γ being constants whose mass-dimension is positive, such that from (168) we obtain(185)Vr=-4γr.The corresponding solution for the mimetic field is found to be an elliptic function, and explicit expressions can be given only in limiting cases. When r≈rs, one can neglect the linear correction in (184) to recover the solution in (172) and (173) with k=1. Thus, the mimetic field reads(186)ϕr≪rsγ≃±ir1-rsr+rs2log2r1+1-rsr-rs,or, by expanding the result around r=rs,(187)ϕr≃rs≃ϕs±2irsr-rs,r≃rs-ϕs-ϕ24rs,with ϕs=±(irs/2)log(rs). The explicit form of the potential at r≃rs, ϕ≃ϕs is found to be(188)Vϕ≃ϕs≃-4γrs-γϕs-ϕ2rs3.On the other hand, for large distances, we may ignore the Newtonian term in (184) and from the second equation in (183) we derive(189)arsγ≪r≃c14+6γr+3c21+γr-c22+3γrarctan1+γr1+γr,with c1,2 dimensional constants. We can choose c1=1/4 and c2=0, and in the given limit the metric simply reads(190)ds2rsγ≪r≃-1+3γr22dt2+dr21+γr+r2dθ2+sin2θdϕ2.The corresponding expressions for the field and the potential are given by(191)ϕr≃±2i1+rγγ,r≃-4+γ2ϕ24γ,Vϕ≃16γ24+γ2ϕr2.Here, we must require 4/γ2<|ϕ|2 in order to guarantee the positivity of r. 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 as(192)g00rsγ≪r≃-1+3γr22.Therefore, the corresponding Newtonian potential Φ(r)=-(g00(r)+1)/2 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 γ2r2 is sufficiently small. On the other hand, the linear term (for γ>0) 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:(193)br=1-rs2r2,with rs being once more a positive dimensional constant. From (168) we get(194)Vr=-2rs2r4,and the metric is fixed by making use of the second equation in (183) which leads to(195)ar=c1+c2arctanr/rs2-r21-rs2/r2.The metric signature for rs<r is preserved when c2=0, and by choosing c1=1 we obtain(196)ds2=-dt2+dr21-rs2/r2+r2dθ2+sin2θdϕ2.The field and the potential are found to be(197)ϕ=±ir1-rs2r2,r=rs2-ϕ2,Vϕ=-2rs2rs2-ϕ22.The 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 r=rs. Moreover, the metric satisfies the following requirements:(1)
g00(r) and g11-1(r) are well defined for all rs≤r.(2)
g00(r) is regular on the throat with g00+(rs)=g00-(rs) and g00+′(rs)=g00-′(rs).(3)
g11-1(rs)=0 and 0<g11-1(r) for all rs<r.(4)
Giveng11(r)-1=[1-b~(r)/r], we have b~+′(rs)=b~-′(rs)<1.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 fixa(t)2b(r), which encodes the physical Newtonian potential, in order to apply mimetic gravity to reproduce the phenomenology of galactic rotation curves.


## 6.1. Vacuum Solutions

In this subsection we setV(ϕ)=0. From (168) we immediately see that(172)br=k-rsr,with rs being a mass scale, either positive or negative, while the second field equation (169) with (162) and (166) leads to(173)ar=a1+a21-krs/r1-krsrlogrr01+1-krsr-1,k=±1,(174)ar=a1+a22rr03/2+3,with a1, a2 being constants and r0 a length scale. If a2=0, namely, R=0 and λ=0, we recover the topological Schwarzschild solution of General Relativity. When a2≠0, we can pose a1=0. Thus, the solution for flat topology (k=0) is as follows:(175)ds2=-a222rr03/2+32r~srdt2+dr2r~s/r+r2dρ2+ρ2dϕ2,0<r~s,where r~s=-rs has to be positive to preserve the metric signature. The Ricci scalar is nonzero and reads(176)R=-6r~2r3+3r/r03/2r03.This solution presents a naked singularity at r=0, as for the corresponding topological case k=0 of the Schwarzschild metric.The spherical case (k=1) is the more interesting and is given by(177)ds2=-a221-rsrlogrr01+1-rsr-12dt2+dr21-rs/r+r2dθ2+sin2θdϕ2,where we have introduced the polar coordinates θ, ϕ in the angular part (see also [79]). In this case the Ricci scalar reads(178)R=-1r21-rs/rlogr/r01+1-rs/r-1.For rs<0 the metric is regular everywhere, but when 0<rs, the solution is regular and preserves the signature; only in the region rs<r the solution is regular and preserves the signature; since for r<rs the metric coefficient g00(r)≡-a2(r)b(r) acquires an imaginary part. In this respect, we note that the special choice rs=r0 when r<rs leads to(179)g00r≡ar2br=-a22rsr-1arctanrsr-1+12,which is real and negative.In general, we observe that when0<rs, the point r=rs cannot represent a horizon like for the Schwarzschild space-time, due to the fact that g11(r)≡1/b(r) diverges but g00(r)=-a22 and the Ricci scalar are regular (we also cannot associate with r=rs any thermodynamical quantity such as temperature T=(a(r)/2)db(r)/dr which diverges at r=rs); moreover, for r<rs the metric becomes imaginary (or, in the special case rs=r0, it acquires the signature (−+++)) and does not allow a wormhole description.Finally, the topological casek=-1 reads(180)ds2=-a221+rsrlogrr01+1+rsr-12dt2+dr2-1-rs/r+r2dρ21+ρ2+ρ2dϕ2,and the Ricci scalar is given by(181)R=-1r21+rs/rlogr/r01+1+rs/r-1.In this case, for 0<rs, the metric coefficient g11(r)≡1/b(r) is negative and the solution is unphysical. On the other hand, for rs<0, we obtain 0<g11(r) when r<-rs, but g00(r)≡-a2(r)b(r) becomes imaginary (except for the special choice r0=-rs).


## 6.2. Nonvacuum Solutions

In this subsection, we will consider the caseV(ϕ)≠0. We immediately see that, for V(ϕ)=2Λ with Λ a cosmological constant, one solution of (168) and (169) with Λ=0 is the topological Schwarzschild de Sitter solution:(182)br=k-rsr-Λr23,ar=a1,with rs, a1 being constants.In considering other solutions, we will take the spherical casek=1 in (165). By using (167) one readily obtains from (170) and (171) that(183)ϕr=±i∫drbr,4ddra′rbr3/2r=arr2brdVrdr,where we treat the potential as a function of r. These two equations with (168) can be used to reconstruct the potential when a choice for b(r) is made. We will now provide some examples of the reconstruction technique.Let us consider a linear modification to the Schwarzschild metric:(184)br=1-rsr+γr,with rs, γ being constants whose mass-dimension is positive, such that from (168) we obtain(185)Vr=-4γr.The corresponding solution for the mimetic field is found to be an elliptic function, and explicit expressions can be given only in limiting cases. When r≈rs, one can neglect the linear correction in (184) to recover the solution in (172) and (173) with k=1. Thus, the mimetic field reads(186)ϕr≪rsγ≃±ir1-rsr+rs2log2r1+1-rsr-rs,or, by expanding the result around r=rs,(187)ϕr≃rs≃ϕs±2irsr-rs,r≃rs-ϕs-ϕ24rs,with ϕs=±(irs/2)log(rs). The explicit form of the potential at r≃rs, ϕ≃ϕs is found to be(188)Vϕ≃ϕs≃-4γrs-γϕs-ϕ2rs3.On the other hand, for large distances, we may ignore the Newtonian term in (184) and from the second equation in (183) we derive(189)arsγ≪r≃c14+6γr+3c21+γr-c22+3γrarctan1+γr1+γr,with c1,2 dimensional constants. We can choose c1=1/4 and c2=0, and in the given limit the metric simply reads(190)ds2rsγ≪r≃-1+3γr22dt2+dr21+γr+r2dθ2+sin2θdϕ2.The corresponding expressions for the field and the potential are given by(191)ϕr≃±2i1+rγγ,r≃-4+γ2ϕ24γ,Vϕ≃16γ24+γ2ϕr2.Here, we must require 4/γ2<|ϕ|2 in order to guarantee the positivity of r. 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 as(192)g00rsγ≪r≃-1+3γr22.Therefore, the corresponding Newtonian potential Φ(r)=-(g00(r)+1)/2 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 γ2r2 is sufficiently small. On the other hand, the linear term (for γ>0) 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:(193)br=1-rs2r2,with rs being once more a positive dimensional constant. From (168) we get(194)Vr=-2rs2r4,and the metric is fixed by making use of the second equation in (183) which leads to(195)ar=c1+c2arctanr/rs2-r21-rs2/r2.The metric signature for rs<r is preserved when c2=0, and by choosing c1=1 we obtain(196)ds2=-dt2+dr21-rs2/r2+r2dθ2+sin2θdϕ2.The field and the potential are found to be(197)ϕ=±ir1-rs2r2,r=rs2-ϕ2,Vϕ=-2rs2rs2-ϕ22.The 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 r=rs. Moreover, the metric satisfies the following requirements:(1)
g00(r) and g11-1(r) are well defined for all rs≤r.(2)
g00(r) is regular on the throat with g00+(rs)=g00-(rs) and g00+′(rs)=g00-′(rs).(3)
g11-1(rs)=0 and 0<g11-1(r) for all rs<r.(4)
Giveng11(r)-1=[1-b~(r)/r], we have b~+′(rs)=b~-′(rs)<1.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 fixa(t)2b(r), 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 asa-3. 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 r, 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 vrot(r)∝1/r, 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 f(R)=Rn 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 redefine(198)a~r=ar2br,within the metric (165) with k=1: namely,(199)ds2=-a~tdt2+dr2br+r2r2dθ2+sin2θdϕ2.Since(200)ar=a~rbr,a′r=12a~rbra~′r-a~rb′rbr,by using (168) the second equation in (183) reads(201)ddra~′rbr-a~rb′rra~r=a~r-b′′rr-2r1-br.We choose the following ansatz for a~(r):(202)a~r=1-rsr+γr-λ0r2,with rs, λ0, and γ being positive dimensional constants. Given the metric element g00(r)=-a~(r), the Newtonian potential reads(203)Φr=-g00r+12.Thus, the Newtonian potential associated with (202) can be described in the following way: (1)
At small distances, the metric leads to a classical Newtonian termrs/r.(2)
At very large distances, the “cosmological constant” termλ0r2 emerges, reflecting the fact that the metric is immersed in a cosmological de Sitter background.(3)
At intermediate distance the linear termγr 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λ0r2 term is negligible; hence, the Newtonian potential (203) reads(204)Φr≃-rs2r1-γr2rs.Correspondingly, the rotational velocity profile reads(205)vrot2≃vNewt2+γc2r2,where we reintroduced the speed of light c and vNewt is the contribution expected from the luminous matter component. Therefore, on sufficiently large scales, vrot does not fall-off as per the Keplerian result 1/r but increases slightly as r. 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 γ0), 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,∝γr. 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 r might be sufficiently large that the de Sitter term ∝r2 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 by(206)v2≃vNewt2+γc2r2-λ0c2r2.Clearly, 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 v2 cannot go negative, bound orbits are no longer possible on scales greater than R~γ0/2λ0. This could provide a dynamical explanation for the maximum size of galaxies, determined by the interplay between the linear (γ) and the quadratic (λ0) 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 derive(207)br=1-rs/r+γr-λ0r21-3rs/r+γr/3+c0/r21-3rs/2r+γr/22,with c0 being a constant. The (on-shell) form of the potential inferred from (168) is quite involved and results in(208)Vr=-23r22r-3rs+γr2354rs2r-27rs3+171γrs2r2-8γ2λ0r7+r416γ+7rsγ2+324rsλ0+4rsr3-17γ-108rsλ0+r6γ3-44γλ0+6r5γ2-12λ0+12rsγλ0-12c0-rs+2r1+rsγ+2r3γ2+λ0-γλ0r4+3r2γ-3rsλ0.Moreover, 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γ=λ0=0 (corresponding to small distances), one has(209)a~r≃1-rsr,br≃4c0+rr-3rsr-rsr2r-3rs2.Thus, if we set(210)c0=9rs24,we have(211)a~r≃1-rsr,br≃1-rsr;namely, we recover the vacuum Schwarzschild solution of General Relativity. Analogously to (186) and (187) we derive(212)ϕr≃rs≃ϕs±2irsr-rs,r≃rs-ϕs-ϕ24rs,with ϕs=±(irs/2)log[rs]. In this case the potential behaves as given in the following:(213)Vϕ≃ϕs≃-32γ3rs+13γϕ-ϕs2rs3.In the limit rs=γ=c0=0 (cosmological scales) we obtain(214)a~r=br≃1-λ0r2,which corresponds to the static patch of the de Sitter solution. Then, the field assumes the form(215)ϕ≃±iarcsinλ0rλ0,r≃±sinλ0ϕλ0.We note that 0<r as long as 0<r<H0-1, where H0-1=1/λ0 is the cosmological horizon of the de Sitter solution with positive cosmological constant. The corresponding behaviour of the potential is(216)Vϕ≃6λ0∓4γ3λ0sinλ0ϕ+4λ0sinλ0ϕ.This result with γ=0 is consistent with [73].Finally, in the limitrs=λ0=c0=0 of galactic scales, the metric reads(217)a~r≃1+γr,br≃41+γr3+γr32+γr2.In this case, the field is given by(218)ϕ≃±i2γ33+4γr+γ2r2,r≃-6∓9-12γ2ϕ23γ,while the potential in (208) behaves as(219)Vr≃-2γ16+6γr+γ2r23r2+γr3.The explicit reconstruction of the potential at the galactic scale must take into account that only the solutions with the positive sign inside r yield a positive radius for 0<γ. Thus, we obtain(220)Vϕ=-23γ227-4γ2ϕ2+29-12γ2ϕ23-4γ2ϕ23/2-6+9-12γ2ϕ2.In conclusion, we are able to describe, in limiting cases, the behaviour of the potential V(ϕ) 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 χ2 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 λ0 parameters in mimetic gravity (202) turn out to be [554, 555]:(221)γ≃3.06×10-30cm-1,λ0≃9.54×10-54cm-2.The parameter λ0 is best expressed as ~(100 Mpc)−2, 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 distanceRlast and rotational velocity vlast of the outermost data in the rotation curves, it can be shown that the combination γlast≡vlast/c2Rlast 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:(222)a~r=1-rsr+γrm-λ0r2,with m being a positive real parameter, with the full metric given by (201); namely,(223)br=42+mr2-122+mrsr+92+mrs2-4m-2γr2+mr-rs+γr1+m-λ0r32+mr3rs-2r+m-2γr1+m2.In deriving this last expression we have set the integration constant in such a way as to recover the Schwarzschild solution in the limit γ=λ0=0, namely, at short distances. On the other hand, at large distances, in the limit rs=γ=0, we once more find the static patch of the de Sitter solution. Instead, at galactic scales, in the limit rs=λ0=0, we obtain(224)a~r=1+γmr,br≃41+γrm2+m+2γrm-mγrm2+mm-2γrm-22.The potential can be found only in an implicit way and its on-shell expression results in(225)Vr=2m2γrm-2γ2m2r2m+m8-6γrm-4γ2r2m+42+3γrm+γ2r2m2+mm-2γrm-23.The trivial case m=2 corresponds to a de Sitter-like solution with negative cosmological constant (for 0<γ):(226)a~r=br=1+γr2,m=2.In this case the field is given by ϕ=±iarcsinh[γr]/γ and the potential reads V(ϕ)≃-6γ.The non-Newtonian correction we have considered has recently been studied in the context off(R) gravity in [556]. In the low-energy limit of power-law F(R)∝Rn gravity, a rm correction to the Newtonian potential emerges, with m related to the power n as(227)m=12n2-7n-1-36n4+12n3-83n2+50n+16n4-4n+2.Stability of the potential at large distances and constraints from Solar System tests require 0<m<1. In [556] the above model was fitted to 15 LSB galaxies, and an excellent fit is achieved, with best-fit value m=0.817 (corresponding to n=3.5). 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 m=0.7.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 ofb(r) 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 mimeticf(R) 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 theirprima 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.



---

*Source: 3156915-2017-03-02.xml*