Fab Four: When John and George play gravitation and cosmology

Scalar-tensor theories of gravitation have recently regained a great interest after the discovery of the Chameleon mechanism and of the Galileon models. The former allows, in principle, to reconcile the presence of cosmological scalar fields with the constraints from experiments at the Solar System scale. The latter open up the possibility of building inflationary models that, among other things, do not need ad hoc potentials. Further generalizations have finally led to the most general tensor-scalar theory, recently dubbed the"Fab Four", with only first and second order derivatives of the fields in the equations of motion and that self-tune to a vanishing cosmological constant. This model has a very rich phenomenology that needs to be explored and confronted with experimental data in order to constrain a very large parameter space. In this paper, we present some results regarding a subset of the theory named"John", which corresponds to a non-minimal 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 a non-minimal, but non-derivative, coupling between scalar field and Ricci scalar, a term named"George"in the Fab Four terminology. In this way, we find a more sensible inflationary model, and, by performing a post-newtonian expansion of spherically symmetric solutions, we derive the set of equations that constrain the parameter space with data from experiments in the solar system.


I. INTRODUCTION
The Galileon theory has recently emerged as an effective theoretical realization of the Dvali-Gabadaze-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 f (R) theories (as in, e.g. [3,9]), also non-minimal derivative couplings to the curvature [11,12]. 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 non-linearities of the scalar sector is that they trigger the Vainshtein mechanism that makes 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 [18,19].
Generally speaking, the Galileon model has opened the way to new models for cosmology, including inflationary or late-accelerated ones [21][22][23]. In fact, almost all sort of cosmological scenarios are possible depending on which Galileon model is chosen, see e.g. [25].
In this paper, we wish to focus on a subclass of models dubbed "purely kinetic gravity" [25][26][27][28], or also"John" in the "Fab Four's" terminology [22]. This also appears as a special case of the very general Galileon class considered in [24].
We begin with a review of the equations of motion in a flat and empty Universe. We then complete the analysis made in [26] by computing the number of e-folds in Sec. II A. We show that a kinetically driven inflationary phase requires highly transplanckian values for the initial field velocity, which rule out the model. Moreover, in Sec. II B we provide for a detailed analysis of the noghost and causality conditions during cosmic evolution, and we show that the theory becomes acausal for such arXiv:1203.4446v3 [gr-qc] 6 Nov 2012 transplanckian values unless the coupling constant is vanishingly small or if the initial Hubble constant is transplanckian. 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. Indeed, although ϕ = 0 is only one possible solution of the equations of motion in the vacuum, the theory is forced to reduce to pure GR when the same equations are solved inside a compact body.
In Sec. III, we extend our considerations to a more general model, in which we include a coupling of the scalar field to the Ricci scalar (and thus, eventually, to matter). This corresponds to the "John plus George" combination in the Fab Four terminology. The equations of motion, together with no-ghosts and causality conditions, are derived and studied numerically for a cosmological background. We find that inflation is possible provided both coupling constants are positive. In Sec. III B we derive the field equations in spherical and static symmetry, and we put some constraints on the parameter space using Solar System tests. We conclude in Section IV with some remarks and open problems.

II. JOHN LAGRANGIAN
We begin our analysis by considering the action where R is the curvature scalar, G µν the Einstein tensor, ψ collectively denotes the matter degrees of freedom coupled to the metric g µν , and κ = 8πG N . 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 can one can see by setting K(X) = X, G 3 = 0, G 4 = 1/(2κ), G 5 = κγϕ/2.

A. Inflation with John
As it was realized in [26], 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 t = −∞ (see [26]), 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 complete the analysis found in [26] by providing the number of e-foldings as a function of the free parameters of the theory. The cosmological equations in vacuum derived from Eq. (1) read where a is the scale factor and α = ln a. The system can be partially decoupled to allow for a numerical integration whose results are shown in Sec. II C. Isolating the second order derivatives, we find which 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 Fig. 1 for both positive and negative γ. For both signs of γ, the EoS tends to −1 in the high energy limit (κ|φ| 1), 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, Eq. (2a) can be inverted to 3α 2 = κφ 2 /(2 + 3κ 2 γφ 2 ), which needs to be positive since the Hubble constant is a real number. Thus, γ < 0 implies that |φ| < −2/3γ, which, in turn, means that we always have w ϕ > 0. Therefore, the scalar field cannot even start in the w ϕ < 0 region if γ < 0. More generally speaking, accelerated phases driven by a scalar field in this model require γ > 0.
In view of these considerations, in the following we shall assume that γ is positive. Eq. (3b), together with the condition that 1 − 3κγα 2 > 0 (see Sec. II B), shows thatφ < 0 for an initially expanding Universe (α > 0). Hence, the velocity of the field decreases with time and ω ϕ is driven towards ω = +1. We characterize the end of inflation by the instant t end , at which ω ϕ = −1/3 1 . Under the assumption that κγφ 2 is initially large, one may derive an analytical (approximate) solution for the scale factor a(t) = e α(t) and the scalar field at early times. First, we write Eq. (2a) as 3α 2 = κφ 2 /(2 + 3κ 2 γφ 2 ) so that where H is the Hubble constant. Integration yields the approximate scale factor Using Eq. (5) in Eq. (3b), and expanding according to Now recall that inflation ends at t end such that ω ϕ (t end ) = −1/3. This corresponds to κ 2 γφ 2 = ζ end , where as can be shown by solving Eq. (4). Using Eq. (7), the condition κ 2 γφ 2 = ζ end reduces to from which one finds Replacing this in the expression (6) for the scale factor leads to a end a i ∼ κ 2 γφ 2 i ζ end 1 6 .
We finally impose that inflation lasts for a number of e-folds N = ln(a end /a i ) greater than 60. This gives a relation between the initial velocity of the fieldφ i 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, Eq. (5) is also of interest, as it fixes the Hubble constant at the beginning of the inflationary phase H i ∼ 1/3 √ κγ. Therefore, the last equation might also be written as It follows that the "natural" initial conditions in Planckian units are not allowed. On the contrary, a natural value for the initial expansion H i = 1 (and thus γ ≈ 0.11) requires an extremely high transplanckian value for the initial velocity of the fielḋ ϕ i ∼ 10 78 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, namelẏ ϕ i 10 51 in natural units.

B. 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 (non hyperbolic, ie. non causal, or carrying negative energy, ie. 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 welldefined, 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 c 2 ϕ ≥ 0 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 superluminal [17,30]. The two conditions above are indeed equivalent to the requirement that the effective metric along which the scalar field propagates, is hyperbolic with the same (mostly +) signature than g µν . The conditions (16) and (17) are best analyzed in terms of the reduced dimensionless variables x(t) = κφ and y(t) = √ κα. The first one requires γy 2 < 1/3 and it is automatically satisfied if γ < 0, but also if γ > 0 provided that there is a maximal value for the Hubble constant H =α < 1/ √ 3κγ. The second one reduces to an algebraic condition with the help of Eqs. (2) and (3), namely (20) In summary, Eqs. (16)(17) require γ > 0 and |y| < 1/ √ 3γ or γ < 0 and the two possibilities |x| > 1/ √ −γ or |x| < −2/3γ. These conditions are therefore less restrictive than the one implied by the Friedmann equation (see previous section).
The curvature background implies non-standard 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 [31], whose conventions we follow. These conditions, namely Eqs. (23), (25)(26)(27) of [31] reduce to rather simple algebraic constraints in our case, after the necessary manipulation using the equations of motion (2) and (3) where Q T and c T are defined as their scalar counterpart Q ϕ and c 2 ϕ above, 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 x defined above, and one typically finds that both positive and negative values for γ are allowed on a given range |x| < ξ γ , where typically ξ γ behaves as O 1/ |γ| , see e.g. Figs. 2 and Fig. 3 below. Hence, large (transplanckian) values for |x| are only allowed for small |γ| 1. This means that the space for possible velocities of the field x = κφ needs to be typically subplanckian, unless γ is vanishingly small. This will be linked to the results found earlier, where transplanckian initial velocity were required for a successful inflation, leading to negative squared speeds c 2 S and c 2 T in that epoch. This is shown in Fig. 5.  In passing, we note that the claim made in the literature (see e.g. [24,28,32]) according to which only the subclass γ < 0 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 γ < 0), 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.

C. 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 non-causal behavior for the scalar and tensor perturbations of the metric. These plots (Figs. 4-7) have been obtained by numerical integration for an initial condition ofφ i = 10 in natural units in the case γ = 1, anḋ ϕ i = 0.1 for γ = −1. Scale Factor a Figure 4. Evolution of the acceleration parameter q = +äa/ȧ 2 vs. the scale factor, for γ = 1. Also shown in the graph: the evolution of the scalar field EoS wϕ and Qϕ and c 2 ϕ . The scalar field becomes superluminal during the transition from a de Sitter Universe to a to stiff matter-dominated one. However, as discussed in the text, it remains hyperbolic and, thus, causal. Moreover, it carries positive energy.

D. 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  Figure 5. Evolution of the parameters for metric perturbations (scalar and tensorial parts), respectively QS, c 2 S , QT , c 2 T , as defined in the text vs the scale factor, for γ = 1. As found theoretically, the initial high velocity of the field drives both the speed of (scalar and tensor) metric perturbations to imaginary values, thus signaling a breakdown of hyperbolicity for metric perturbations. (Here it happens when the corresponding curves terminate, since the y axis is in logarithmic scale).  transplanckian values for the initial velocity of the field, which, in turn, implies various instabilities, as discussed above.
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, ϕ = 0 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, i.e. ϕ = 0 for all r 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 non-vanishing 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: R(1 + ϕ), namely the "George" term in the Fab Four terminology.

III. 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 G eff = G/(1 + √ κϕ) might easily become negative in this model, meaning that the action chosen here shall trivially lead to dynamical pathologies for ϕ sufficiently large and negative 2 .
Such an argument would call in favor of defining a better coupling function F in the Georges term F (ϕ)R. 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 F 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 F .
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 with 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 .

A. Cosmological behavior
The equations of motion in a flat, empty Universe, derived from Eq. (23) are given in Appendix A, see Eqs. (A1-A3). 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 Appendix A and B.
The numerical results are the following. The case = 1 and γ = 1 is pretty similar to the case John alone, see Figs. (8) and (9) below. Inflation thus occurs in the case > 0 and γ > 0, 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 te scope of this paper would determine whether the addition of the George term help in solving the naturalness problem encountered with John alone in Sec. II.
The case = −1, γ = 1 is clearly pathological for the various no-ghosts and no-acausal conditions, as seen on conformal transformation, this is equivalent to the usual Brans-Dicke condition 2w + 3 > 0, where w = 2 (1 + √ κϕ) here. Then, our model with γ = 0 would indeed be pathological if ϕ ≤ −(3/(2 2 ) + 1)/( √ κ) . However the γ term introduces new terms in the equation of motion for the scalar field which invalidate such a conclusion in the general case γ = 0.  Figure 8. Evolution of q, wϕ, Qϕ and c 2 ϕ with the scale factor. We observe the same transition from inflation to stiff matter for the scalar field. Initial conditions areφi = 100 and ϕi = 1.  Figure 9. Evolution of QT , Qϕ and c 2 T , c 2 ϕ with the scale factor. These two speeds are negative in the early universe. Fig. (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 fail to exhibit any acceleration at all, see Figs. (10).

B. 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 A 2 = 1 + i a i /r i , B 2 = 1 + i b i /r i and ϕ = p 0 + i p i /r i and of the field equations in powers of 1/r. By equating the coefficients of equal powers of r we find p zr 5 , is the Schwarzschild radius of the central body, z = 1 + p0 Mp , andγ = γp 2 1 . In the expansion above, we neglected higher order terms in r s /r, in , in p 1 /(rM p ) and inγ/(rM p ) (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 p 0 = ϕ(r → ∞) (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 p 1 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 [33] with suitable asymptotic and regularity conditions, together with, e. g. 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 systeme constraints can be derived.
The final step consists then in computing observables effects from the metric Eqs. (24)(25), 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 3 . 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  Figure 10.
Cosmological evolution of q, Qϕ, c 2 ϕ , QS, c 2 S , QT , c 2 T with the scale factor. (Top): the case = −1 and γ = 1 corresponds to a "double inflation". The scenario is pathological in many respects: singularity in the scalar field velocity (see top left) while the scalar field EoS is imaginary. Also (see top right, the y axis is logarithmic): c 2 T < 0 and there are periods for which c 2 s < 0 and QT < 0. (Bottom): the case = 1 and γ = −1. The model is well behaved but does not accelerate the expansion, ie. as in the John alone model. The universe is actually in a super-stiff regime, and hence in a highly decelerating phase.
the change of the argument of the perihelion is exhibited by a time average (over an orbital period T ) The results of this whole procedure is given by where a is the semi-major axis of the orbit, e its eccentricity, and n its mean motion n = 2π T . The advance of perihelion of Solar System planets is very tightly constrained by planetary ephemerides. In particular, IN-POP10a gives constraints on supplementary advances of perihelia (see Table 5 of [34]). Table I gives the value of the α i coefficients appearing in the expression of the advance of perihelion (27) 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 (24)(25). One way to obtain such constraint is to follow the strategy and to use the software of [35][36][37]. The main idea presented in these papers is to simulate radioscience observables directly from the spacetime 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 (24)(25) 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 (10 −14 ) allows us to give order of magnitude of constraints on the theory. Fig. 11 represents the incompressible signatures produced by parameters entering the metric (24-25) 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 (24)(25). Fig. 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 It should be noted that radioscience constraints are significantly better forγ M 4 p z and p 2 1 M 2 p z as compared to perihelia advances. On the other hand, the constraint from the ephemerides onγ 2 M 4 p z 2 is significantly better while the constraint on p1 zMp 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 of the Georges coupling) is necessary.

IV. CONCLUSIONS
In this paper, we have explored some phenomenology associated to 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 α1 (mas/cy/m 4 ) α2 (mas/cy/m 4 ) α3 (mas/cy/m 2 ) α4 (mas/cy/m) α5 (mas/cy/m 2 ) dω/dt (mas/cy    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 non-minimal 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 infla-tion. 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 ϕ = 0 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, i.e. 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 are non-trivial 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 above, this paper is the first of a series that aim at a systematic study of the "Fab Four" phenomenology, in order to isolated 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 adressed in the litterature 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.