#### Abstract

Scalar-tensor theories of gravitation attract again a great interest since the discovery of the Chameleon mechanism and of the Galileon models. The former allows reconciling the presence of a scalar field with the constraints from Solar System experiments. The latter leads to inflationary models that do not need ad hoc potentials. Further generalizations lead to a tensor-scalar theory, dubbed the “Fab Four,” with only first and second order derivatives of the fields in the equations of motion that self-tune to a vanishing cosmological constant. This model needs to be confronted with experimental data in order to constrain its large parameter space. We present some results regarding a subset of this theory named “John,” which corresponds to a nonminimal derivative coupling between the scalar field and the Einstein tensor in the action. We show that this coupling gives rise to an inflationary model with very unnatural initial conditions. Thus, we include the term named “George,” namely, a nonminimal, but nonderivative, coupling between the scalar field and Ricci scalar. We find a more natural inflationary model, and, by performing a post-Newtonian analysis, we derive the set of equations that constrain the parameter space with data from experiments in the Solar System.

#### 1. Introduction

The Galileon theory has recently emerged as an effective theoretical realization of the Dvali-Gabadadze-Porrati model (DGP) [1]. The subsequent developments eventually led to the definition of the most general second order scalar-tensor theory that includes, besides the usual terms of scalar-tensor or theories (as in, e.g., [2–10]), also nonminimal derivative couplings to the curvature [11–13]. This general theory, explored for the first time by the pioneering work of Horndeski many years ago [14], provides a wide framework that virtually encompasses all scalar models analyzed so far in the literature.

One interesting aspect of the nonlinearities of the scalar sector is that they trigger the Vainshtein mechanism that makes it possible to build viable cosmological models with sufficiently small effects at local scales to evade Solar System constraints. Therefore, we have at hand an interesting alternative to chameleons [15], in which the parameter space for local gravity does not overlap with the one allowed by cosmic acceleration, as recently shown in [16]. In addition, inflationary phases are permitted without the introduction of ad hoc scalar potentials, making these models more “natural.” Nonlinearities are also responsible for new phenomena in the dark sector including, for instance, sub/superluminality and/or effective violation of the Null Energy Condition, thus allowing for a stable and well-defined phantom-like phase [17–20]. Generally speaking, the Galileon model has opened the way to new models for cosmology, including inflationary or late-accelerated ones [21–24]. In fact, almost all sorts of cosmological scenarios are possible depending on which Galileon model is chosen; see, for example, [25].

In this paper, we first focus on a subclass of models dubbed “purely kinetic gravity” [25, 27–29] or also similar to “John Lagrangian” in the “Fab Four’s” terminology. (Note that in the Fab Four Lagrangian, the scalar field has no standard kinetic term. It propagates however, via its derivative coupling to the curvature tensors.) [22, 30]. More generally speaking, it is a special case of the very general Galileon class considered in [23]. We mostly investigate the inflationary phase and its naturalness, taking into account no-ghost constraints and causality conditions. In a next step, a generalization of the theory including a coupling of the scalar field to the matter—via, in the Jordan frame, a Brans-Dicke-like coupling to the scalar curvature—is introduced. The analysis of the inflationary phase is then reconducted in such a more realistic theory, while a preliminary analysis of solar system constraints is discussed.

We begin in Section 2 with a review of John’s equations of motion in a flat and empty Universe. We then complete the analysis made in [27] by computing the number of e-folds in Section 2.1. We show that a kinetically driven inflationary phase requires highly trans-Planckian values for the initial field velocity, which rule out the model. Moreover, in Section 2.2 we provide a detailed analysis of the no-ghost and causality conditions during cosmic evolution, and we show that the theory becomes acausal for such transplanckian values unless the coupling constant is vanishingly small or if the initial Hubble constant is trans-Planckian. In passing, we prove that some claims made in the literature regarding the sign of the coupling constant are wrong. We argue that this model is to be discarded also because there are no reasons why the scalar field should be generated at all in the first place. This is particularly clear in an astrophysical context where the theory is forced to reduce to pure GR (), when the same equations are solved inside a compact body.

In Section 3, we thus extend our considerations to a more general model, in which we also include a coupling of the scalar field to the Ricci scalar. This corresponds to the “John plus George” combination in the Fab Four’s terminology. We derive and study numerically the equations of motion for the matter-free cosmological background as well as the no-ghosts and causality conditions. We find that inflation is possible provided both coupling constants are positive. In Section 3.2 we derive the field equations in spherical and static symmetry, and we briefly sketch how to put some constraints on the parameter space using Solar System tests. We conclude in Section 4 with some remarks and perspectives.

#### 2. John Lagrangian

We begin our analysis by considering the action where is the curvature scalar, the Einstein tensor, collectively denotes the matter degrees of freedom coupled to the metric , and . In this paper we use the mostly plus signature; is a dimensionless parameter whereas has the dimension of an inverse length. This action is a special case of the generalized Galileon one presented in [23], as one can see by setting , , , .

##### 2.1. Inflation with John

As it was realized in [27], this model allows for a quasi de Sitter inflation with a graceful exit without the need for any specific scalar potential. Inflation is essentially driven by kinematics, and it crucially depends on the initial high velocity of the field, as we will shortly see. Although, in principle, the inflationary solutions begin at (see [27]), we will consider the action as an effective model only valid from few Planck times after an unknown transplanckian phase. Our first concern is to establish whether the model accommodates an inflationary phase together with reasonable assumptions for the initial conditions at that time. This section thus completes the analysis found in [27] by providing the number of e-foldings as a function of the free parameters of the theory. The cosmological equations in vacuum derived from (1) read where is the scale factor and . The system can be partially decoupled to allow for a numerical integration whose results are shown in Section 2.3. Isolating the second order derivatives, we findwhich can be integrated numerically in a straightforward way. The effective equation of state (EoS) for the scalar field can be obtained from its stress-energy tensor or, more simply, by comparing our equations of motion directly to the standard Friedmann equations. After some algebra we find which is plotted in Figure 1 for both positive and negative . For both signs of , the EoS tends to in the high energy limit (), so that a large initial velocity for the scalar field will result in a quasi de Sitter phase. However, only the case of positive can lead to inflation. Indeed, (2a) can be inverted to , which needs to be positive since the Hubble constant is a real number. Thus, implies that , which, in turn, means that we always have . Therefore, the scalar field cannot even start in the region if . More generally speaking, accelerated phases driven by a scalar field in this model require .

In view of these considerations, in the following we shall assume that is positive. Equation (3b), together with the condition that (see Section 2.2), shows that for an initially expanding Universe (). Hence, the velocity of the field decreases with time and is driven towards . We characterize the end of inflation by the instant , at which (throughout this paper we assume a vanishing cosmological constant). Under the assumption that is initially large, one may derive an analytical (approximate) solution for the scale factor and the scalar field at early times. First, we write (2a) as so that where is the Hubble constant. Integration yields the approximate scale factor Using (5) in (3b), and expanding according to , gives and Now recall that inflation ends at such that . This corresponds to , where as can be shown by solving (4). Using (7), the condition reduces to from which one finds Replacing this in the expression (6) for the scale factor leads to We finally impose that inflation lasts for a number of e-folds greater than . This gives a relation between the initial velocity of the field and , namely: which is the crucial condition for a successful (purely kinetic-driven) inflationary phase. We see that it involves a rather unusual very large pure number. In order to discuss naturalness, (5) is also of interest, as it fixes the Hubble constant at the beginning of the inflationary phase . Therefore, the last equation might also be written as It follows that the “natural” initial conditions and in Planckian units are not allowed. On the contrary, a natural value for the initial expansion (and thus ) requires an extremely high transplanckian value for the initial velocity of the field in natural units.

It is not even possible to obtain a Planckian value for the initial velocity in this model, since, in any event, the initial Hubble constant will be greater than the one today. This implies that in such an inflationary scenario, must be less than the Hubble radius today, still implying a very unnatural bound for the initial velocity, namely, in natural units.

##### 2.2. Theoretical Constraints

The characteristic feature of Galileon models is the derivative coupling of the scalar field to the metric. This implies a direct coupling between scalar field and metric degrees of freedom or, in other words, that the scalar field propagation explicitly depends on the metric background and vice versa. Therefore, there might exist backgrounds for which the propagation becomes pathological (nonhyperbolic, i.e., noncausal, or carrying negative energy, i.e., ghosts). In the following, we restrict ourselves to the cosmological background in a flat universe with line element and we explore the conditions for the theory to be well defined, for both scalar field and metric perturbations. We start with the scalar field, whose action can be written as where as it needs to be positive for the scalar field to carry positive energy. The field propagates at a squared speed given by The condition is necessary to ensure that the scalar field equation of motion remains hyperbolic, ultimately expressing its causal behavior regardless of whether the scalar field perturbations are sub- or super-luminal [17, 31]. The two conditions previously mentioned are indeed equivalent to the requirement that the effective metric along which the scalar field propagates, is hyperbolic with the same (mostly +) signature than .

The conditions (16) and (17) are best analyzed in terms of the reduced dimensionless variables and . The first one requires , and it is automatically satisfied if , but also if provided that there is a maximal value for the Hubble constant . The second one reduces to an algebraic condition with the help of (2a)–(3b), namely: In summary, (16)-(17) require and or and the two possibilities or . These conditions are therefore less restrictive than the one implied by the Friedmann equation (see previous section).

The curvature background implies nonstandard propagation for the scalar degree of freedom. In a quite similar fashion, the scalar field background modifies the standard spectrum of metric perturbations. Similar conditions for the avoidance of ghosts and euclidean metrics also exist. These have been derived in full generality in a very wide class of Galileon models in [32], whose conventions we follow. These conditions, namely, , of [32], reduce to rather simple algebraic constraints in our case, after the necessary manipulation using the equations of motion (2a)–(3b):where and are defined as their scalar counterparts and , but stand for the tensor perturbations of the metric field. Similar conditions need to hold for the scalar part of the metric perturbations, namely:This whole set of equations is difficult to reduce algebraically because of the last one. However, one might easily plot the six functions of defined previously, and one typically finds that both positive and negative values for are allowed on a given range , where typically behaves as ; see, for example, Figures 2 and 3. Hence, large (transplanckian) values for are only allowed for small . This means that the space for possible velocities of the field needs to be typically subplanckian, unless is vanishingly small. This will be linked to the results found earlier, where transplanckian initial velocity was required for a successful inflation, leading to negative squared speeds and in that epoch. This is shown in Figure 5.

In passing, we note that the claim made in the literature (see, e.g., [24, 29, 33]) according to which only the subclass is a ghost-free theory is wrong (at least in the background considered here). Notice that the scalar field is well defined although being a phantom in certain regime (in the case ), a situation reminiscent of the one discussed in [18]. However, as shown previously, the Friedmann equation actually prevents the scalar field to enter this regime.

##### 2.3. Numerical Results

In this section, we quickly show the cosmological behavior in the John model, for both positive and negative . As discussed before, the negative case leads only to a decelerating Universe: the phantom regime is not an acceptable initial condition (as it entails an imaginary Hubble constant), and neither can be reached. Only positive leads to acceleration and to an inflationary phase in the early Universe, a drawback being the presence of noncausal behavior for the scalar and tensor perturbations of the metric. These plots (Figures 4, 5, 6, and 7) have been obtained by numerical integration for an initial condition of in natural units in the case , and for .

##### 2.4. Discussion

We have established that kinetically driven inflation in the Galileon theory involving the simplest coupling to Einstein tensor is not viable. It requires unnatural transplanckian values for the initial velocity of the field, which, in turn, implies various instabilities.

This model has anyway another serious drawback. In the absence of any direct coupling to the Ricci scalar, there is no reason why the scalar field should be generated at all (even in presence of cosmological matter fluid). In other words, is always a solution in this class of models, whatever the matter content is. At local scales, this problem appears in the following way. We checked numerically that a relativistic star (with flat asymptotic conditions) must be described by pure GR; that is, for all is the only solution that is regular at the center of the star. To conclude, the model considered so far is trivial in the sense that it cannot be different than GR, except if one imposes nonvanishing initial conditions for the scalar field at early times.

This acts as a leitmotiv for a more realistic model, studied in the next section, where we add to the Lagrangian a direct coupling to the Ricci scalar: , namely, the “George” term in the Fab Four terminology, and where only a particular function will be analyzed here.

#### 3. George and John

We now consider the extended model given by where is a dimensionless, free parameter. Of course this is not the most general coupling one might consider but it is anyway reminiscent of Brans-Dicke coupling. Notice that one might worry that the effective gravitational constant might easily become negative in this model, meaning that the action chosen here shall trivially lead to dynamical pathologies for sufficiently large and negative. (In fact, what matters in the case is that the scalar field propagates positive energy in the Einstein frame. Performing a conformal transformation, this is equivalent to the usual Brans-Dicke condition , where here. Then, our model with would indeed be pathological if . However the term introduces new terms in the equation of motion for the scalar field which invalidate such a conclusion in the general case .).

Such an argument would call in favor of defining a better coupling function in the Georges term . However, this would be a misleading conclusion here, since the John term couples the derivatives of the metric and of the scalar field, thus impacting their propagation. Therefore, only the entire set of no-ghost conditions together with causal propagation conditions (positivity of the squared velocities) for both the scalar and the metric perturbations can decide which regions of the configuration space are well behaved. This is done in the following sections (on a cosmological background), based on the conditions derived in Appendix B. In this light, the function chosen above is just the simplest one could have chosen and might furthermore be understood as retaining only the first term in a series expansion of a more general function .

The cosmological evolution in this theory is typically a function of four parameters, the initial value of the field, its velocity, and of the two dimensionless parameters and . It goes beyond the scope of the present paper to provide a comprehensive study of this parameter space. However we highlight some essential features of the model thanks to numerical results displayed in the next section. We provide both the cosmological evolution and the analysis of causality and positivity of energy within subclasses of the model, depending on the signs of and .

##### 3.1. Cosmological Behavior

The equations of motion in a flat, empty Universe, derived from (22), are given in Appendix A; see . We extended the analysis of the no-ghost and causality conditions to this more general framework, and we also provide the scalar field EoS; see Appendices A and B.

The numerical results are the following. The case and is pretty similar to the case John alone; see Figures 8 and 9. Inflation thus occurs in the case and , but the acausal behavior still shows up in the very early Universe. The number of e-folds is a function of the two initial conditions for the field and its velocity and the dimensionless parameters and . A further analysis that goes beyond the scope of this paper would determine whether the addition of the George term helps in solving the naturalness problem encountered with John alone in Section 2.

The case , is clearly pathological for the various no-ghosts and no-acausal conditions, as seen in Figure 10. Actually this theory leads to a double inflation scenario (see the acceleration parameter): the Universe transits from one de Sitter phase to another one and experiences in between a super acceleration phase. Finally, the case with negative is similar to what we found for John alone: the theory is well defined, ghost-free, and causal, but fails to exhibit any acceleration at all; see Figure 10.

**(a)**

**(b)**

**(c)**

**(d)**

##### 3.2. Solar System

In this section we show how to derive Solar System constraints on the free parameters of the model “John + George.”

In a first step we solve asymptotically the field equations in spherical symmetric and isotropic coordinates. The line element reads and the equations of motion for the components and for are given in Appendix C. The post-Newtonian analysis begins with the expansion of the metric components as , and and of the field equations in powers of . By equating the coefficients of equal powers of we find where is the Schwarzschild radius of the central body, , and . In the expansion previously mentioned, we neglected higher order terms in , in , in and in (which means we suppose these terms to be smaller than 1). We recall that, in our conventions, and are dimensionless parameters. The asymptotic scalar field value (in GeV) is a free parameter that can eventually be connected to the cosmological evolution of the scalar field.

The next step then amounts to relate the dimensionless parameter to the scalar charge of the central body and thus eventually to the couplings and . This could be determined numerically by solving the equations in the interior of the body [34] with suitable asymptotic and regularity conditions, together with, for example, a polytropic equation of state for the star’s interior. Such a discussion about possible spontaneous scalarization goes however beyond the scope of the present paper and is left for future studies. In this section we only sketch how solar system constraints can be derived.

The final step consists then in computing observable effects from the metric Equations (23)–(24b), which show a deviation from standard General Relativity in the Solar System. A first constraint comes from the anomalous perihelion shift. First of all, we determine the geodesic equation (to first order in the metric deviation) to find the planetary equations of motion (when computing the geodesics, we suppose that matter is minimally coupled to the metric). The anomalous term is treated perturbatively to find the perihelion shift with the Gauss equations, which determines the evolution of the different orbital parameters due to perturbative forces. In particular, we derive the rate of change of the argument of the perihelion . Finally, the secular term of the change of the argument of the perihelion is exhibited by a time average (over an orbital period ) The results of this whole procedure are given by with where is the semi major axis of the orbit, its eccentricity, and its mean motion . The advance of perihelion of Solar System planets is very tightly constrained by planetary ephemerides. In particular, INPOP10a gives constraints on supplementary advances of perihelia (see Table 5 of [26]). Table 1 gives the value of the coefficients appearing in the expression of the advance of perihelion (26) and the constraints coming from INPOP10a.

The most stringent constraints in the case considered here are obtained by data from Mercury and read

These constraints are obtained by considering only deviations on the dynamics (i.e., on the equations of motion). Other constraints can be derived using propagation of light rays in the solar system. For example, radioscience experiments include light propagation through a Shapiro-like term that can be derived from the expression of the metric (23)–(24b). One way to obtain such constraint is to follow the strategy and to use the software of [35–37]. The main idea presented in these papers is to simulate radioscience observables directly from the space-time metric so that the software includes both deviations on the dynamics and deviations on light propagation. In particular, we use this software to simulate a two-way Doppler link between Earth and Cassini spacecraft from May 2001 on the metric (23)–(24b) and to analyze them in GR by fitting the initial conditions of the spacecraft. The residuals that emerge are the incompressible signature produced by the alternative theory considered on Doppler signal for Cassini. Comparing this signal to the Doppler accuracy of the mission () allows us to give order of magnitude of constraints on the theory. Figure 11 represents the incompressible signatures produced by parameters entering the metric (23)–(24b) on Cassini Doppler. The three sharp peaks occurred at solar conjunctions. The order of magnitude of the residuals observed in this figure is larger than Cassini accuracy, which means that the values of the parameters should be smaller than the indicated values.

We run a set of simulations with different values for the parameters appearing in the metric (23)–(24b). Figure 12 represents the evolution of the maximal Doppler residuals obtained in Cassini signal as function of metric parameters. Requiring the residuals to be lower than Cassini accuracy () gives the boundary values

It should be noted that radioscience constraints are significantly better for and as compared to perihelia advances. On the other hand, the constraint from the ephemerides on is significantly better while the constraint on is of the same order of magnitude.

We thus see that solar systems tests of GR will put severe constraints on the free parameters of the model. As already stressed; however, the full numerical study of Solar’s interior is required, while a thorough analysis of possible spontaneous scalarization (depending also on the Georges coupling) is necessary.

#### 4. Conclusions

In this paper, we have explored some phenomenology associated with a subset of the “Fab Four” scalar-tensor theory. The philosophy behind this preliminary work is that we cannot forget about solar system constraints on the parameter space, even when we deal with inflationary solutions. Traditional models of inflation rely upon the fact that the inflaton field decays, at some stage, into ordinary matter through some reheating mechanism. Therefore, the scalar-tensor nature of inflationary gravity is lost very soon in the evolution of the Universe. On the opposite, in the models studied in this paper the scalar field should live and show its effects until nowadays. Therefore, the parameter space determined by constraints from cosmological observations must overlap with the one determined by solar system tests. The “Fab Four” theory has many parameters with a very rich phenomenology, and with this paper we begin an ambitious plan for its systematic study.

In this work, we contained ourselves to the cases John and John plus George. The John case represents a theory with a nonminimal derivative coupling between scalar field and Einstein tensor. It was already known that this model admits inflationary solutions with a graceful exit. Here, we constrained the coupling constant by showing that it must be positive in order to have successful inflation. If this is the case, however, very unnatural initial conditions are required. In particular the field velocity, which is related to the energy density, must be huge compared to the Planck scale, rendering the theory no longer trustworthy. Negative values for are permitted but do not allow for inflation while ghost states might appear. The most serious problem however comes from the fact that the model turns out to be trivial when one tries to solve the equations of motion inside a compact object. Indeed, we found that the only solution with a finite field at the centre is everywhere.

These facts have convinced us to extend the theory to include the term named “George,” which is nothing but a coupling between the scalar field and the Ricci scalar. The parameter space is now two-dimensional, and we solved numerically the equations of motion for a cosmological background. The main result is that the sign of the two coupling constants must be positive in order to have an inflationary phase with graceful exit. We reserve the analysis of the naturalness problem, that is, the need for extreme initial conditions, for future work. We also performed a post-Newtonian analysis of the theory by solving the equations of motion by imposing static and spherical symmetry and expanding the fields for large radial distance. The aim was to put some constraint on the free parameters, which now include also the asymptotic value of the scalar field and the scalar charge, which now makes sense as there is nontrivial solution also for the interior of a compact object. The results still allow for a large parameter space; therefore future work is necessary in order to improve the constraints.

As mentioned previously, this paper is the first of a series that aim at a systematic study of the “Fab Four” phenomenology, in order to isolate the experimentally viable sectors of the theory. Besides these aspects, there are several issues that deserve further investigations. For example, if these models truly lead to inflationary cosmology and graceful exit, we will need an alternative reheating mechanism. Maybe, the complexity of the theory acts as an effective potential for the scalar field, which resembles the usual power-law forms of usual inflation. But if this is not the case, then we need to find alternative explanations. Another aspect that needs to be studied is the relevance of these modifications of gravity in terms of late-time cosmology, as we expect modifications driven by the scalar field. This would include, at the background level, the study of tracking solutions and of the convergence mechanism towards GR, if any. This has not yet been addressed in the literature for the John Lagrangian (but see [38]). The study of cosmological perturbations, in particular CMB spectra and large-scale structures, might then further reduce the parameter space. Finally, gravitational effects might be relevant at galactic scales and give rise to alternative explanations to the anomalous galactic rotation curves.

In conclusion, we believe that the recent developments in scalar-tensor theories of gravity have opened the door to new and intriguing research directions, and we are confident that many interesting results will be obtained in the near future.

#### Appendices

#### A. Cosmological Equations

In terms of the reduced variables , and , the equations of motion for a flat and empty universe derived from action Equation (22) are which can be decoupled in the following way: The scalar field EoS is given by with

#### B. Ghost Conditions

The coupling to the Ricci scalar, “George,” does not change the analysis made for the scalar field sector of the theory. Thus the two following conditions still hold: For the metric perturbations, we derive, based on equations , , , and of [32] for the tensorial part, and also for the scalar part of the metric perturbations, while their squared speed leads to

#### C. Spherically Symmetric Equations of Motion

We derive the equations of motion for the action (22) with a spherically symmetric and static field configuration. We consider the metric (23), and we replace its components in the Lagrangian. With the Noether theorem, we find the equations of motion for the fields , , , and . Finally, we impose the gauge , and we find three equations plus a Hamiltonian constraint that read

#### Acknowledgments

J.-P. Bruneton is FSR/COFUND postdoctoral researcher at naXys. M. Rinaldi is supported by a grant of ARC 11/15-040 convention. A. Kanfon thanks ARC 11/15-040 for travel support. A. Hees is FNRS Research Fellow. A. Hees thanks B. Lamine, P. Wolf, C. Le Poncin-Lafitte, V. Lainey, S. Reynaud, and M. T. Jaekel for useful discussions about radioscience simulations within Solar System.