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 f(R) 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 (φ=0), 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
(1)S=∫-gd4x[R2κ-12(gμν+κγGμν)∂μφ∂νφ]+Smat[ψ,gμν],
where R is the curvature scalar, Gμν the Einstein tensor, ψ collectively denotes the matter degrees of freedom coupled to the metric gμν, and κ=8πGN. 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 K(X)=X, G3=0, G4=1/(2κ), G5=κγφ/2.
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 t=-∞ (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 (2a)3α˙2=κφ˙22(1-9κγα˙2),(2b)2α¨+3α˙2=-κφ˙22(1+κγ[3α˙2+2α¨+4α˙φ¨φ˙-1]),(2c)1a3ddt(a3φ˙(1-3κγα˙2))=0,where a is the scale factor and α=lna. 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 find(3a)α¨=(3κγα˙2-1)2×3α˙2+(κφ˙2/2)(1-9κγα˙2)1-3γκα˙2+(κ2γφ˙2/2)(1+9κγα˙2),(3b)φ¨=-3α˙φ˙(1+κ2γφ˙2)1-3γκα˙2+(κ2γφ˙2/2)(1+9κγα˙2),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
(4)ωφ=(2+3κ2γφ˙2)(1-κ2γφ˙2)2+3κ2γφ˙2+3κ4γ2φ˙4,
which is plotted in Figure 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, (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.
log linear plot of ωφ as a function of the dimensionless variable κφ˙ for γ=1 and γ=-1 (dashed).
In view of these considerations, in the following we shall assume that γ is positive. Equation (3b), together with the condition that 1-3κγα˙2>0 (see Section 2.2), 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 tend, at which ωφ=-1/3 (throughout this paper we assume a vanishing cosmological constant). 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 (2a) as 3α˙2=κφ˙2/(2+3κ2γφ˙2) so that
(5)H=α˙≃13κγ,
where H is the Hubble constant. Integration yields the approximate scale factor
(6)a~aiexp(t-ti3κγ).
Using (5) in (3b), and expanding according to κ2φ˙2≫1, gives φ¨/φ˙≃-1/κγ and
(7)φ˙~φi˙exp(-t-tiκγ).
Now recall that inflation ends at tend such that ωφ(tend)=-1/3. This corresponds to κ2γφ˙2=ζend, where
(8)ζend=16(3+57)≈1.76,
as can be shown by solving (4). Using (7), the condition κ2γφ˙2=ζend reduces to
(9)κ2γφ˙i2exp(-2(tend-ti)κγ)~ζend,
from which one finds
(10)tend-ti=κγ2ln(κ2γφ˙i2ζend).
Replacing this in the expression (6) for the scale factor leads to
(11)aendai~(κ2γφ˙i2ζend)1/6.
We finally impose that inflation lasts for a number of e-folds N=ln(aend/ai) greater than 60. This gives a relation between the initial velocity of the field φ˙i and γ, namely:
(12)φ˙i2≳ζendκ2γexp(360),
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 Hi~1/3κγ. Therefore, the last equation might also be written as
(13)κφ˙i2Hi2≳9ζendexp(360)~10157.
It follows that the “natural” initial conditions Hi=𝒪(1) and φ˙i=𝒪(1) in Planckian units are not allowed. On the contrary, a natural value for the initial expansion Hi=1 (and thus γ≈0.11) requires an extremely high transplanckian value for the initial velocity of the field φ˙i~1078 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, φ˙i≳1051 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
(14)ds2=-dt2+e2α(t)dx2,
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
(15)Sφ=∫a3dtd3xQφ(φ˙2-cφ2a2∇φ2),
where
(16)Qφ=12(1-3κγα˙2)>0
as it needs to be positive for the scalar field to carry positive energy. The field propagates at a squared speed given by
(17)cφ2=1-κγ(2α¨+3α˙2)1-3κγα˙2≥0.
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 super-luminal [17, 31]. The two conditions previously mentioned are indeed equivalent to the requirement that the effective metric
(18)g~μν=gμν+κγGμν,g~μν∇μνφ=0,
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 γy2<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 (2a)–(3b), namely:
(19)1-3γy21-9γy2+54γ2y4≥0⟺(1+γx2)(2+3γx2)2+3γx2+3γ2x4≥0.
In summary, (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 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, (23), (25)–(27) of [32], reduce to rather simple algebraic constraints in our case, after the necessary manipulation using the equations of motion (2a)–(3b):(20a)QT>0⟺1+γx22>0,(20b)cT2≥0⟺1-γx22≥0,where QT and cT are defined as their scalar counterparts Qφ and cφ2, but stand for the tensor perturbations of the metric field. Similar conditions need to hold for the scalar part of the metric perturbations, namely:(21a)QS>0⟺4+6γx2+6γ2x42+3γx2>0,(21b)cS2≥0⟺12+36γx2+19γ2x4-12γ3x6-3γ4x82+3γx2+3γ2x4≥0.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 previously, and one typically finds that both positive and negative values for γ are allowed on a given range |x|<ξγ, where typically ξγ behaves as 𝒪(1/|γ|); see, for example, Figures 2 and 3. 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 was required for a successful inflation, leading to negative squared speeds cS2 and cT2 in that epoch. This is shown in Figure 5.
Plot of the six conditions Qi and ci derived previously, as a function x. Here γ=1. Allowed values for the field velocity are typically |x|<𝒪(γ-1/2), as made clear by drawing similar plots while varying γ.
Plot of the six conditions Qi and ci derived previously, as a function x. Here γ=-1. Allowed values for the field velocity are typically |x|<𝒪(|γ|-1/2).
In passing, we note that the claim made in the literature (see, e.g., [24, 29, 33]) 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.
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 φ˙i=10 in natural units in the case γ=1, and φ˙i=0.1 for γ=-1.
Evolution of the acceleration parameter q=+a¨a/a˙2 versus 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.
Evolution of the parameters for metric perturbations (scalar and tensorial parts), respectively, QS, cS2, QT, cT2, as defined in the text versus 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).
Evolution of the acceleration parameter q versus the scale factor, for γ=-1 with initial condition satisfying |xi|=0.1<-2/3γ. The field starts with an EoS wφ~0.5 and the Universe only decelerates.
Evolution of the parameters for metric perturbations (scalar and tensorial parts), respectively, QS, cS2, QT, cT2, as defined in the text versus the scale factor, for γ=-1. This model with γ<0 is well behaved but does not accommodate inflation.
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, φ=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; that is, φ=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 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: F(φ)R, namely, the “George” term in the Fab Four terminology, and where only a particular function F will be analyzed here.
3. George and John
We now consider the extended model given by
(22)S=∫-gd4x×[R2κ(1+ϵκφ)-12(gμν+κγGμν)∂μφ∂νφ]+Smat[ψ,gμν],
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 Geff=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. (In fact, what matters in the case γ=0 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 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.).
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 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 (A.1). 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 ϵ=1 and γ=1 is pretty similar to the case John alone; see Figures 8 and 9. 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 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.
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.
Evolution of QT, Qφ, and cT2, cφ2 with the scale factor. These two speeds are negative in the early universe.
The case ϵ=-1, γ=1 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.
Cosmological evolution of q, Qφ, cφ2, QS, cS2, QT, cT2 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 (a)) while the scalar field EoS is imaginary. Also (see (b), the y axis is logarithmic): cT2<0 and there are periods for which cs2<0 and QT<0. Bottom: the case ϵ=1 and γ=-1. The model is well behaved but does not accelerate the expansion, that is, as in the John alone model. The universe is actually in a super stiff regime and hence in a highly decelerating phase.
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
(23)ds2=-A(r)2dt2+B(r)2(dr2+r2d2Ω),
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 A2=1+∑iai/ri, B2=1+∑ibi/ri and φ=p0+∑ipi/ri and of the field equations in powers of 1/r. By equating the coefficients of equal powers of r we find (24a)A2=1-rsr+rs22r2+ϵ2p124Mp2z2r2+p1ϵrs224Mpzr3-p12rs24Mp2zr3+34γ-ϵ2Mp4z2r4-rsγ-8Mp4zr5,(24b)B2=1+rsr-2ϵp1Mpzr-p124Mp2zr2-γ-4Mp4zr4,where rs=2GM/c2 is the Schwarzschild radius of the central body, z=1+(ϵp0/Mp), and γ-=γp12. In the expansion previously mentioned, we neglected higher order terms in rs/r, in ϵ, in p1/(rMp) and in γ-/(rMp) (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 p0=φ(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 p1 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 T)
(25)〈dωdt〉=1T∫0Tdωdtdt.
The results of this whole procedure are given by
(26)〈dωdt〉=-α1γ-Mp4z-α2γ-ϵ2Mp4z2-α3p12Mp2z-α4p1ϵzMp-α5p12ϵ2z2Mp2,
with (27a)α1=n16a48+24e2+3e4(1-e2)4,(27b)α2=9n8a3rs4-3e2-e4(1-e2)4,(27c)α3=n8a24+e2(1-e2)2,(27d)α4=2na(1-e2),(27e)α5=n4ars(1-e2),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, INPOP10a gives constraints on supplementary advances of perihelia (see Table 5 of [26]). Table 1 gives the value of the αi coefficients appearing in the expression of the advance of perihelion (26) and the constraints coming from INPOP10a.
Values of the different αi coeffcients (27a)–(27e) for planets and constraints given by INPOP10a for supplementary advances of perihelia (values given in [26]).
α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) INPOP10a
Mercury
3.2×10-32
4.85×10-24
8.83×10-11
19.4
8.2×10-4
0.4±0.6
Venus
7.69×10-34
2.53×10-25
9×10-12
3.89
1.65×10-4
0.2±1.5
Earth
1.3×10-34
5.9×10-26
2.9×10-12
1.73
7.33×10-5
-0.2±0.9
Mars
1.36×10-35
9.13×10-27
6.77×10-13
6.1×10-1
2.58×10-5
-0.04±0.15
Jupiter
1.51×10-38
3.56×10-29
9.06×10-15
2.81×10-2
1.19×10-6
-41±42
Saturn
5.33×10-40
2.3×10-30
1.08×10-15
6.16×10-3
2.61×10-7
0.15±0.65
The most stringent constraints in the case considered here are obtained by data from Mercury and read (28a)-3.12×1031m4<γ-Mp4z<6.25×1030m4,(28b)-2.06×1023m4<γ-ϵ2Mp4z2<4.12×1022m4,(28c)-1.13×1010m2<p12Mp2z<2.26×109m2,(28d)-5.16×10-2m<p1ϵzMp<1.03×10-2m.
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 (10-14) 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.
Incompressible anomalous signal due to the different terms in the metric (23)–(24b) on Cassini Doppler signal with indicated values. The produced residuals are larger than Cassini accuracy (which is 10-14); therefore these values for the parameters are ruled out.
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 (10-14) gives the boundary values (29a)|γ-Mp4z|=|γp12Mp4(1+(ϵp0/Mp))|<3.65×1026m4,(29b)|γ-ϵ2Mp4z2|=|γp12ϵ2Mp4(1+(ϵp0/Mp))2|<1.15×1026m4,(29c)|p12Mp2z|=|p12Mp2(1+(ϵp0/Mp))|<3.53×108m2,(29d)|p1ϵzMp|=|p1ϵMp(1+(ϵp0/Mp))|<5.56×10-2m.
Representation of the maximal Doppler residuals signal due to John + George modification of gravity parametrized by the four parameters in the metric (23)–(24b) for the Cassini mission between Jupiter and Saturn. The red (dashed) lines represent the assumed Cassini accuracy.
It should be noted that radioscience constraints are significantly better for γ-/Mp4z and p12/Mp2z as compared to perihelia advances. On the other hand, the constraint from the ephemerides on γ-ϵ2/Mp4z2 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 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 φ=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, 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.
AppendicesA. Cosmological Equations
In terms of the reduced variables x(t)=κφ˙, y(t)=κα˙ and z(t)=1+ϵκφ(t), the equations of motion for a flat and empty universe derived from action Equation (22) are
(A.1)6ϵxy+x2(-1+9γy2)+6y2z=0,4x(ϵ+γκx˙)y+x2(1+3γy2+2γκy˙)+6y2z+2κ(ϵx˙+2y˙z)=0,3y(x-2ϵy-3γxy2)+κ(x˙-3γx˙y2-3(ϵ+2γxy)y˙)=0,
which can be decoupled in the following way:
(A.2)x˙=-3x[ϵx+4(ϵ2+γx2)y+7ϵγxy2]+6y(-2x+ϵy)z2κ1/2(3ϵ2+12ϵγxy+γx2(1+9γy2)+(2-6γy2)z)y˙=(2ϵxy(1-15γy2)+x2[-1+3γy2(4-9γy2)]-6y2(2ϵ2+z-3γy2z))×(2κ1/2[3ϵ2+12ϵγxy+γx2(1+9γy2)+(2-6γy2)z]2κ1/2)-1.
The scalar field EoS is given by
(A.3)ωφ=-x(3γx2+2z)ND
with
(A.4)N=6x(γx2-z)(3γx2+2z)2+63ϵ3×x2(3ϵ2+3γx2+2z)(7γx2+2z)-18ϵ4(7γx3+2xz)-23ϵx2(3ϵ2+3γx2+2z)×(33γ2x4+16γx2z-4z2)+9ϵ2(15γ2x5+4γx3z-4xz2)D=[-3ϵx+3x2(3ϵ2+3γx2+2z)]2×[18γ3x6+30γ2x4z+24γx2z2+8z3+63ϵγx(γx2+2z)x2(3ϵ2+3γx2+2z)+3ϵ2(3γ2x4+4z2)].
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:
(B.1)Qφ=12(1-3γy2)>0cφ2=1-γ(3y2+2κy˙)1-3γy2≥0.
For the metric perturbations, we derive, based on equations (23), (25), (26), and (27) of [32]
(B.2)QT>0⇒z+γx22>0,cT2≥0⇒z-γx22≥0
for the tensorial part, and also
(B.3)QS>0⇒3ϵ2+12ϵγxy+9γ2x2y2+2z+γ(x2-6y2z)>0
for the scalar part of the metric perturbations, while their squared speed cs2≥0 leads to
(B.4)2y(γx22+z)2(ϵx+3γx2y+2yz)+2x(ϵ+γκx˙)(γx2+2z)(ϵx+3γx2y+2yz)+12(γx2-2z)(ϵx+3γx2y+2yz)2-2κ(γx22+z)2×[ϵ(x˙+2xyκ)+3γx(2x˙y+xy˙)+2y˙z]≥0.
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 A, B, R, and φ. Finally, we impose the gauge R=r, and we find three equations plus a Hamiltonian constraint that read
(C.1)0=(2φ′2Br2B′′+4φ′B2rφ′′+4φ′Br2φ′′B′-5φ′2B′2r2+2φ′2B2)γκ2+(B′2-4B3φB′′r2-2B3B′r2φ′-4B4φ′r-8B3φB′r-2B4r2φ′′+2B2φB′2r2)ϵκ1/2-8B3rB′-4B3r2B′′+2B2r2B′2-φ′2κ2B4r2,0=(4φ′2B′ArB+2φ′2B2A-3φ′2B′2Ar2+4φ′Br2φ′′AB′+8φ′2B2rA′+2φ′2B′A′Br2+2φ′2A′′r2B2+4φ′B2rφ′′A+4φ′A′r2φ′′B2+2φ′2ABr2B′′)γκ2+(-4B3AB′r2φ′+2B2AφB′2r2-8B4Aφ′r-8φrA′B4-8B3AφB′r-4B3AφB′′r2-4φB′A′r2B3-6A′r2φ′B4-4φA′′r2B4-4B4Ar2φ′′B′2)ϵκ1/2-4B3Ar2B′′-4r2A′′B4-8rA′B4-8B3ArB′-4A′r2B′B3+2B2Ar2B′2-φ′2κ2B4Ar2,0=(B′24φ′A′B3+4φ′AB′B2+4φ′B′Ar2B′′B+4φ′B′′A′r2B2+4φ′′rA′B3+4φ′rA′′B3+4φ′′B′A′r2B2+8φ′B′A′rB2+4φ′B′′ArB2+4φ′′B′ArB2-4φ′B′2ArB-6φ′B′2A′r2B-6φ′B′3Ar2+4φ′B′A′′r2B2+2φ′′B′2Ar2B)γκ+(B′2-8B′ArB4-4B5rA′-4r2AB′′B4-2B′A′r2B4-2B5r2A′′+2B′2Ar2B3)ϵκ-1/2+2A′r2φ′B5+2B5r2φ′′A+4φ′B5Ar+2B′r2φ′AB4,0=(2φ′2A′B2+2rφ′2A′′B2-6φ′2B′2Ar+4φ′Aφ′′B2-2φ′2B′AB-4φ′2A′B′rB+4φ′rB′Aφ′′B+4rφ′A′φ′′B2+2φ′2rB′′AB)γκ2+(-4φB′AB3+4rφB′2AB2-4B4rA′φ′-4B4φA′′r-4B4φ′′Ar-4B4φA′-4B4φ′A-4φB′′ArB3B′2)ϵκ1/2+4rB′2AB2-4rAB′′B3-4B′AB3-4B4A′-4B4A′′r-2B4φ′2Arκ2.
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.
NicolisA.RattazziR.TrincheriniE.Galileon as a local modification of gravityWagonerR. V.Scalar-tensor theory and gravitational wavesDamourT.Esposito-FareseG.Tensor multiscalar theories of gravitationDamourT.NordtvedtK.Tensor-scalar cosmological models and their relaxation toward general relativityPerrottaF.BaccigalupiC.MatarreseS.Extended quintessenceNojiriS.OdintsovS. D.Modified gravity with negative and positive powers of curvature: unification of inflation and cosmic accelerationChibaT.1/R gravity and scalar-tensor gravityCapozzielloS.FrancavigliaM.Extended theories of gravity and their cosmological and astrophysical applicationsSotiriouT. P.FaraoniV.f(R) theories of gravityCliftonT.FerreiraP. G.PadillaA.SkordisC.Modified gravity and cosmologyDeffayetC.Esposito-FareseG.VikmanA.Covariant GalileonDeffayetC.GaoX.SteerD.ZahariadeG.From k-essence to generalized GalileonsDeffayetC.DeserS.Esposito-FareseG.Generalized Galileons: all scalar models whose curved background extensions maintain second-order field equations and stress tensorsHorndeskiG. W.Second-order scalar-tensor field equations in a four-dimensional spaceKhouryJ.WeltmanA.Chameleon cosmologyHeesA.FuzfaA.Combined cosmological and solar system constraints on chameleon mechanismBabichevE.MukhanovV.VikmanA.k-Essence, superluminal propagation, causality and emergent geometryCreminelliP.NicolisA.TrincheriniE.Galilean genesis: an alternative to inflationBabichevE.Galileon accretionDeffayetC.PujolasO.SawickiI.VikmanA.Imperfect dark energy from kinetic gravity braidingChowN.KhouryJ.Galileon cosmologyCharmousisC.CopelandE. J.PadillaA.SaffinP. M.General second order scalar-tensor theory, self tuning, and the Fab FourKobayashiT.YamaguchiM.YokoyamaJ.Generalized g-inflationTsujikawaS.Observational tests of inflation with a field derivative coupling to gravitySaridakisE. N.SushkovS. V.Quintessence and phantom cosmology with nonminimal derivative couplingFiengaA.LaskarJ.KuchynkaP.MancheH.DesvignesG.GastineauM.CognardI.TheureauG.MechanicsC.AstronomyD.The INPOP10a planetary ephemeris and its applications in fundamental physicsSushkovS. V.Exact cosmological solutions with nonminimal derivative couplingGubitosiG.LinderE. V.Purely kinetic coupled gravityGermaniC.KehagiasA.Ultraviolet-protected inflationCharmousisC.CopelandE. J.PadillaA.SaffinP. M.Self-tuning and the derivation of the Fab FourBrunetonJ. P.Causality and superluminal behavior in classical field theories: applications to k-essence theories and modified-Newtonian-dynamics-like theories of gravityDe FeliceA.TsujikawaS.Conditions for the cosmological viability of the most general scalar-tensor theories and their applications to extended Galileon dark energy modelsGermaniC.KehagiasA.New model of inflation with nonminimal derivative coupling of standard model higgs boson to gravityDamourT.Esposito-FareseG.Nonperturbative strong-field effects in tensor-scalar theories of gravitationHeesA.LamineB.ReynaudS.JaekelM. T.Le Poncin-LatteC.LaineyV.FuzfaA.CourtyJ. M.DehantV.WolfP.Radioscience simulations in general relativity and in alternative theories of gravityHeesA.WolfP.LamineB.ReynaudS.JaekelM. T.Le Poncin-LatteC.LaineyV.FuzfaA.DehantV.Radioscience simulations in General Relativity and in alternative theories of gravityProceedings of Les Rencontres de Moriond 20112011HeesA.WolfP.LamineB.ReynaudS.JaekelM. T.Le Poncin-LatteC.LaineyV.DehantV.AlecianG.BelkacemK.SamadiR.Valls-GabaudD.Testing Gravitation in the Solar System with Radio Science experimentsProceedings of the Annual Meeting of the French Society of Astronomy and Astrophysics (SF2A '11)2011653658CopelandE. J.PadillaA.SaffinP. M.The cosmology of the Fab-Fourhttp://arxiv.org/abs/1208.3373