Probing gravitational non-minimal coupling with dark energy surveys

We investigate observational constraints on a specific one-parameter extension to the minimal quintessence model, where the quintessence field acquires a quadratic coupling to the scalar curvature through a coupling constant ξ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi $$\end{document}. The value of ξ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi $$\end{document} is highly suppressed in typical tracker models if the late-time cosmic acceleration is driven at some field values near the Planck scale. We test ξ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi $$\end{document} in a second class of models in which the field value today becomes a free model parameter. We use the combined data from type-Ia supernovae, cosmic microwave background, baryon acoustic oscillations and matter power spectrum, to weak lensing measurements and find a best-fit value ξ>0.289\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi {>}0.289$$\end{document} where ξ=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi = 0$$\end{document} is excluded outside the 95% confidence region. The effective gravitational constant Geff\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$G_\mathrm{eff}$$\end{document} subject to the hint of a non-zero ξ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi $$\end{document} is constrained to -0.003<1-Geff/G<0.033\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-0.003< 1- G_\mathrm{eff}/G < 0.033$$\end{document} at the same confidence level on cosmological scales, and it can be narrowed down to 1-Geff/G<2.2×10-5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1- G_\mathrm{eff}/G < 2.2 \times 10^{-5}$$\end{document} when combining with Solar System tests.


Introduction
Cosmological observations in support of the late-time cosmic acceleration, such as measurements for the Type-Ia supernovae [1,2], the cosmic microwave background (CMB) radiation [3,4], the large-scale structure and baryon acoustic oscillations [5,6], strongly suggest that the current universe density must consist of some unfamiliar "negative pressure matter," namely, dark energy. Dark energy can be explained by a cosmological constant ( ) with an extremely tiny density ρ ≈ 10 −47 GeV 4 , compared with the typical energy a e-mail: geng@phys.nthu.edu.tw b e-mail: g9522545@oz.nthu.edu.tw c e-mail: ypwu@phys.sinica.edu.tw scales in particle physics. Another possible candidate for such a source, as legitimate as imposing , is referred to a quintessence component whose dynamics is described by a scalar field φ [7].
Quintessence with sufficiently flat potentials exhibits attractor solutions such that wide ranges of initial conditions approach the scalar field dominated universe ( φ = 1). One of the common forms is the exponential potential V (φ) ∼ e −λκφ with κ 2 = 8π G, where the cosmic acceleration is realized if λ 2 < 2 [8][9][10][11]. These late-time attractor solutions, however, require nearly fixed values of the quintessence density ρ φ all along the evolution history of the universe. The fact that ρ φ coincides with the (dust-like) matter density today implies that the initial densities of the two species have to be in some huge yet precise hierarchy, leading to the cosmic coincidence problem [12].
On the other hand, tracker fields are introduced as a specific class of quintessence aiming to solve the cosmic coincidence problem [12,13]. The main advantage of a tracker field is that the initial value of φ can be close to unity to within few orders of magnitude, and thus it accepts the assumption of equipartition after inflation [12,13]. Typical tracker solutions have been found in potentials of the inverse power-law V (φ) ∼ φ −α with α > 0 [8,14,15]. Unfortunately, recent observational constraints on the tracker fields require α 1 [16][17][18][19], which significantly narrows down the tracking basin of attraction [20]. To mimic a cosmological constant by taking α → 0, the present energy density must be realized from the fine-tuned initial condition ρ φ 10 −47 GeV 4 , only to duplicate the cosmic coincidence problem.
Given these considerations, many efforts have been made to address the extended scenario in which the quintessence field is non-minimally coupled to the scalar curvature R ( [21,22]; see also [23][24][25][26][27][28][29][30][31]) through a non-minimal coupling ξ Rφ 2 , where ξ is the coupling constant. There exist track-ing solutions purely induced by the non-minimal coupling such that a very wide range of initial values of ρ φ can evolve toward the late-time attractor of the cosmic acceleration, even if the potential is extremely flat [32]. This gravitational effect is manifest under conditions of the R-boost [22] or the slowroll [33]. In either case, the gravitational tracker solutions are obtained without assuming particular forms of the potentials.
A fine-tuned potential, however, is needed to achieve a preferable cosmic acceleration, despite that the gravitational tracker solutions help to relax the cosmic coincidence problem. In particular, when the potential is exactly flat the extended quintessence model [21] can be cast into the prototype Brans-Dicke (PBD) theory [34,35] but with a timevarying Brans-Dicke parameter ω depending on the value of φ. Constraints on the value of ω have been widely studied from the CMB anisotropy and structure formations [36][37][38][39][40], the parametrized post-Newtonian parameters [41][42][43][44], and the big-bang nucleosynthesis [45][46][47]. These constraints in terms put a tight bound on the non-minimal coupling as |ξ | < 10 −2 in the inverse power-law model of the [14] type, given that the scalar field must reach a Planck scale by now (φ 0 ∼ 10 19 GeV) after rolling on the track from arbitrary initial conditions [22].
In this work, we show that the restriction on ξ shall be significantly reduced in other types of potentials where the present day value of field can be freely determined by its initial conditions. As a concrete example, we put the quintessence field in the slow-roll exponential potential type [8][9][10][11] where the field value is nearly frozen all along the expansion history when λ → 0. We investigate the effects of ξ from an analytic approach to the background equation of state w and to the effective gravitational constants on cosmological scales or inside the Solar System. These modified gravitational constants are sketched by the phenomenological functions μ and , which can be useful for testing the deviation from the CDM model [48][49][50]. We then confront the exponential model with the recent observational data to fit the most preferable value of ξ on cosmological scales. We pick up the absolutely flat model (λ = 0) as our paradigm for which its background dynamics is identical to that of CDM where ξ only modifies the growth dynamics through μ or . It is also interesting that the flat model with a constant potential has only the same free parameters as that of the PBD theory. Finally, we study the possible near future cosmic evolution based on a non-zero ξ .
This article is organized as follows: In Sect. 2, we review basic equations of a non-minimally coupled scalar field. In Sect. 3, we provide the analytic and numerical results of the field evolution in the exponential potential model, including the constant potential model as a special case. We use observational data to obtain the best-fit or constraints on the model parameters, as shown in Sect. 4. We present a summary and discussions in Sect. 5.

Background equations
In this work we shall focus on the one-parameter extension of the minimal dark energy scenario [7] that the quintessence field φ has an explicit coupling to the scalar curvature R given by the action where S M is the action for both the relativistic and the nonrelativistic matter, and with the constant ξ representing the non-minimal coupling parameter which exhibits the minimal value at ξ = 0 and the conformal one at ξ = −1/6. If the potential is simply a constant, say V (φ) = , the model (2) can be cast into the PBD theory [34,35] through the field redefinition = (1 + ξκ 2 φ 2 )/2, leading to with the time evolving Brans-Dicke parameter Therefore a massless extended quintessence is a special scalar-tensor theory that has only the same number of model parameters as that of the PBD theory. Note that the positive energy condition 1 + ξκ 2 φ 2 > 0 guarantees ω > 0 so that the model (2) can satisfy the no-ghost condition ω > −3/2 in the Brans-Dicke theory [35]. The Solar System measurements report a very strong constraint ω > 40,000 on the Brans-Dicke parameter ( [41][42][43]), but this constraint is significantly reduced to ω > 692 on cosmological scales [40]. The variation of the action (2) with respect to the metric g μν results in the modified Einstein field equation where μν = ∇ μ φ∇ ν φ − g μν (∇φ) 2 Here G μν is the Einstein tensor and T μν is the energymomentum tensor of matter. Taking the spatially flat Friedmann-Robertson-Walker (FRW) background, with the scale factor a(t), the background field equations are given by where ρ M and P M are the energy density and pressure of matter, while ρ φ and P φ are the effective energy density and pressure of the scalar field, respectively. The equation of motion for the scalar field is governed by the Klein-Gordon equation, where ≡ ∇ μ ∇ μ and V φ ≡ dV /dφ. In the homogeneous FRW background, the Klein-Gordon equation takes the form of where we have used R = 6Ḣ + 12H 2 . The background equation (13) can be regained from the continuity equation of ρ φ and P φ : so that the non-minimally coupled scalar field is effectively a perfect fluid with an equation of state w φ = P φ /ρ φ . In terms of the e-folding N = ln a, the Klein-Gordon Eq. (13) becomes where the primes are e-folding derivatives, and the Planck unit has been used for φ. Taking H = −3/2(1 + w)H , Eq. (15) can be rewritten as Similarly, the first Friedmann equation (8) can be given in terms of the density fraction as where r ( m ) is the density fraction of radiation (dust-like matter) with today's value m0 = 0.308 [4], and φ is the density of the scalar field, which is further divided into where Note that we have included the kinetic termφ 2 /2 into nc φ since it depends nontrivially on the parameter ξ .

Perturbation functions
Let us characterize the effects of the non-minimal coupling to the gravitational constants in the linear density perturbations. The full perturbed Einstein equations are provided in Sect. Appendix A: (see also [51,52]). For the study of non-linear perturbations one may refer to [53,54] In a non-minimally coupled theory, the scalar field fluctuations give rise to anisotropy between the curvature perturbation and the Newtonian potential (as defined in Eq. (A.1)). On subhorizon scales where k a H and in the Newtonian limit where time derivatives are negligible with respect to spatial derivatives, the anisotropy parameter γ (φ, k) = / is determined by where The modified Poisson equation for the dust-like matter and the relativistic matter are given by [49] respectively, where is the comoving matter density perturbation. In the extended quintessence model (2), we have a scale-independent function for the lensing effect (φ in Planck unit) and a scale-dependent function for the matter growth It is possible to choose a present value φ 0 = 0 such that γ 0 = μ 0 = 0 = 1, which coincide with the predictions in general relativity. Constraints on the deviation from general relativity are 1−μ 0 = −0.05±0.25 and 1− 0 = 0.00±0.14 for the fiducial CDM background expansion at the 68% confidence level [55]. Some tension with the CDM prediction is reported by the combined Planck CMB polarization and low multipole data with the BAO and weak lensing measurement [56]. We can define κ 2 eff = 8π G eff = 8π G to rewrite the Einstein equation (5) as so that the Friedmann equations (8) and (9) take the same forms as those in the minimal quintessence model up to the effective constant G eff whenφ is negligible. Let us consider the usual slow-roll potential of the exponential form where where ξ = 0 leads to BD = μ BD = 1 with ω → ∞. In the PBD case (λ = 0) [40] report 0.981 ≤ 0 ≤ 1.285 at the 99% confidence level from cosmological tests. The Solar System bound on the time variation of the Newton gravi- In general, the local value μ local is determined by the present day value of φ in our Solar System. Yet, if no screening mechanism is assumed, μ local may coincide with the background value μ 0 = μ(φ = φ 0 ) under the weak-field condition with the quasi-static approximation (see for example [57,58]).

Model parameters
Here we outline the generic feature of the parameter dependence in the exponential potential model (26), and we refer the interested reader to subsections for the analytic solutions up to the asymptotic future. In a constant potential case (λ = 0), the background equation of state w is basically indistinguishable from that of the CDM model, but w may notably deviate from −1 in the near future if ξ > 0; see the upper panel of Fig. 1. In cases where λ > 0, the actual value of w given by the numerical result in the lower panel of Fig.  1 can differ from that of the CDM model, and the difference is enhanced if ξ > 0. Viable cosmic expansions can be realized even if ξ > 1.5, but the tracking basin of attraction is tightly restricted. For λ ≥ 0, the function tends to deviate from (converge to) unity since matter became to dominate the universe, given that φ is governed by an increasing (decreasing) mode when ξ > 0 (ξ < 0), as shown in Fig. 2.

Constant potential models
Let us study the case with an absolutely flat potential λ = 0, where V (φ) = V 0 and V φ = 0. In this case, the Klein-Gordon equation (16) is reduced to 1 If w is a constant, this equation can be exactly solved as 1 We assume the background equation of state is a constant with w = 1/3 during the epoch of radiation domination (RD) and w = 0 during matter domination (MD).
where C ± are constants to be determined by the initial conditions and For ξ > 0, one finds L + > 0 and thus e L + N = a L + corresponds to an increasing mode, while e L − N = a L − is a decreasing mode as L − < 0. For ξ = 0, L + = 0 holds, showing the constant mode of a massless scalar field. For ξ < 0, no increasing mode exists since L ± ≤ 0. In what follows we examine the evolution of φ with (a) ξ ≤ 0, (b) 0 < ξ ≤ 3/2 and (c) ξ > 3/2, respectively.
(a) ξ ≤ 0 If ξ = 0, φ is governed by a constant mode and the kinetic energy decays rapidly. Eventually, we find ρ φ = V 0 and w φ = −1 for arbitrary initial conditions. In this limit the quintessence field reproduces the result of a cosmological constant. For If φ is released from rest, the decreasing mode (∝ e L − N ) becomes negligible with the expansion of the universe and the solution (30) can be approximated to be 2 Taking this solution into Eq. (10), the energy density reads Assuming that the initial condition satisfies ρ φ V 0 , the energy density is then reduced to Here, indicates that the equation of state takes the form Given that w φ > 0 in MD when ξ < 0, ρ φ decays faster than the matter density as φ approaches zero. Eventually, the value of φ is frozen around zero and one reproduces the behavior of a cosmological constant with w φ → −1 as In this case, the initial density of φ may be negative. Given that L + > 0 holds during both RD and MD, e L + N is always increasing with time. Assuming that φ is released from rest, we can adopt the solution (32), and the evolution essentially has three stages: (i) During RD or MD where ρ φ is subdominant, the energy density is increasing with time. To see this, we may denote the solution at this stage as where the decreasing mode has been neglected and The energy density with φ = φ 1 takes the form where it can be checked that ρ φ < 0 in both RD and MD. One may deduce w φ = w 1 from Eq. (37) as w 1 (w, ξ ) = w − 2 3 L 1 (w, ξ ), which gives where −1 ≤ w 1 < 0 in MD for 0 < ξ ≤ 3/2. (ii) Once the density of V 0 is dominant, the universe will undergo a de Sitter expansion. Consequently, we can assume that w = −1 with H = 0 so that H = H ds is a constant. Accordingly, Eq. (15) is simplified as The corresponding solution (keeping only the increasing mode) reads where and the energy density (10) is of the form Note that L 2 > 0 since ξ > 0 and therefore, φ remains increasing with time during the V 0 domination epoch. (iii) Since φ 2 2 ∝ a 2L 2 is increasing with time, the subdominant energy density of ρ φ (the second term in the righthand side of Eq. (42)) eventually becomes comparable with V 0 , where the de Sitter expansion is interrupted. To calculate the cosmological evolution after the V 0domination, we study the asymptotic universe, which has a constant equation of state w = w 3 , satisfying Again, by neglecting the decreasing mode, the solution is given by where Provided that the scalar field is the dominant species of the universe, the Friedmann equation yields where the energy density For a sufficiently large N such that φ 3 6/(L 2 3 − 6ξ − 12ξ L 3 ), the left-hand side of Eq. (46) is negligible and the lowest order of the Friedmann equation gives ρ φ = V 0 + ρ nc φ = 0. The asymptotic solution shows a fine cancelation between V 0 and ρ nc φ (the second term in the right-hand side of Eq. (47)). This cancelation implies φ 2 3 ∼ H −2 ∝ a 3(1+w 3 ) , which leads to Given that −1 < w 3 ≤ 0 with 0 < ξ ≤ 3/2, the scalar field will keep overtaking the matter density after the epoch of the de Sitter expansion (the phase of the V 0 -domination), as seen from Fig. 1.
If the initial value of ρ φ is too small such that it remains negative during MD, Eq. (38) shows that w 1 < −1 for ξ > 3/2 and therefore, the energy density ρ φ will always be dominated by the non-minimal coupling term (ρ nc φ ) as given by Eq. (37). In this case, ρ φ is always negative, and the constant potential V 0 is never important so that there exhibits no epoch of the cosmic acceleration.
On the other hand, if ρ φ can be positive before MD, V 0 will dominate the universe and drive a de Sitter expansion. Similarly, this de Sitter expansion shall be interrupted with the growths of φ and |ρ nc φ |, where the resulting solution has w = w φ = (−3 + 2ξ)/[3(1 + 2ξ)]. However, since w φ > 0 for ξ > 3/2, the matter density ρ M may eventually catch up with ρ φ and become one of the dominant density components.
Let us assume that the asymptotic solution for ξ > 3/2 has a constant equation of state w 3 with φ 3 taken the form of Eq. (44). For arbitrary initial conditions, the final attractor is a fine cancelation between ρ M and ρ nc φ , leading to the lowest order Friedmann equation ρ M +ρ nc The condition of the late-time cosmic acceleration restricts the value of φ at the initial time. Since φ is a constant during RD, a non-negative ρ φ shall satisfy 3ξ H 2 eq φ 2 i ≤ κ 2 V 0 3H 2 0 , where H eq is the Hubble parameter at the matter-radiation equality. This condition is approximately φ 2 i ≤ ξ −1 × 10 −17 for ξ > 3/2, assuming that the equality temperature T eq ≈ 5.5 eV.
Since V 0 corresponds to the energy scale of the dark energy domination, one finds that the potential term is negligibly small in the early time (N 0), where the equation of motion (50) coincides with that of the constant potential, as given by Eq. (29). For simplicity, we assume that φ is released from rest so that the solution (32) can be applied with the evolutions discussed as follows.
(a) ξ < 0 For ξ < 0, the non-minimal coupling term in the Klein-Gordon equation tends to pull back the scalar field from rolling down the potential, and the solution of φ is a constant mode in RD and a decreasing mode in MD (see Sect. III A). Even if φ is initially at rest, ρ φ can be much greater than the potential energy, and it takes the form of Eq. (34) with a nearly constant equation of state as given by Eq. (35). In this case ρ φ is decaying with time for N 0, and eventually approaches V (φ) as φ → 0.
Whenever ρ φ = V (φ) is reached, the φ-field will be temporarily frozen by the friction of the Hubble expansion with w φ −1 until ρ φ starts to dominate the universe.
The asymptotic attractor solution is the stable point balanced between V (φ) and the non-minimal coupling term. Taking φ to be asymptotically stabilized at some value where the final energy density reads It can be checked that φ f < 1/ √ −ξ is satisfied for an arbitrary value of λ, and thus ρ φ will not reach the singularity when 1 + ξφ 2 = 0.
The non-minimal coupling term with ξ > 0 acts as an additional force in the Klein-Gordon equation that pushes φ to run downhill, and thus φ is governed by an increasing mode during RD and MD (see Eq. (36)). Let us schematically divide the evolutions into five stages: (i) For ξ > 0, a wide range of initial conditions can lead to a negative energy density during RD and MD with |ρ φ | > V (φ), where ρ φ is given by Eq. (37) and increasing toward V (φ). Once ρ φ = V (φ) is reached, the energy density becomes positive and the value of φ remains increasing due to the effect of the non-minimal coupling. (ii) Once the φ-field dominates the total energy density, the universe experiences a (quasi) de Sitter expansion for a sufficiently small λ. The solution and the energy density of φ at this stage are given by Eqs. (40) and (42), respectively. (iii) Similar to the case of a constant potential, this de Sitter expansion is going to be interrupted with the increase of φ, once the contribution of the non-minimal coupling term in ρ φ is comparable with V (φ). For 0 < ξ < 3/2, the solution that follows the de Sitter expansion is a power-law expansion with w = w φ = (−3+2ξ)/ [3(1+2ξ)] from a temporary fine cancelation of the energy density (that is, ρ φ = 0; see Eq. (48)). (iv) As φ keeps growing, V (φ) is exponentially decaying with time so that at some epoch the potential energy suddenly transits to the kinetic energy, leading to a phase of the super kinetic-energy domination with w φ > 1.
Given that V (φ) becomes negligible, the lowest order of the Friedmann equation gives where φ 4 = C 4 e L 4 N is the solution of the Klein-Gordon equation with the potential term neglected, and L 4 = −3/4(1 − w 4 ) + 9/16(1 − w 4 ) 2 + 3ξ(1 − 3w 4 ). Since Eq. (53) implies L 4 = 6ξ + √ 6ξ(1 + 6ξ), one can find that (v) The density ρ φ once again turns into negative when V (φ) is subsidiary with the increase of φ. We may denote the solution at the final stage as The asymptotic attractor is a scaling solution in which the density of the φ-field closely tracks that of the dust-like matter, that is, ρ φ = −ρ M ∝ a −3 , while the total energy density of the universe ρ tot ρ M is always positive. Since ρ φ ∝ H 2 φ 2 5 ∝ a −3 , we can obtain , and L 5 = 2ξ 1 + 4ξ , (ξ > 0).

Observational constraints
In this section we confront the extended quintessence model with cosmological observations. We define ε ≡ V / φ − 1 to parametrize the fraction of the dark energy density in addition to the potential. The density fraction induced by the non-minimal coupling is then given by | nc φ | ≡ |ε φ |, where ε is positive (negative) for ξ > 0 (ξ < 0). Given that in the potential (26) a constant shift of the field value φ can be absorbed into a redefined V 0 , we will fix V 0 and treat φ 0 as a free parameter to be determined by the present day value ε = ε 0 . To do so we have to further fix λ so that the exponential model has in total 2 + 6 parameters, which are ξ and ε 0 plus the six parameters of the standard spatially flat CDM model [3]. Our optimal choice for the potential is λ = 0 where the background expansion closely reproduces that of the CDM model. The prior ξ > 0 is chosen due to the enhancement on the value of φ at low redshifts, which results in interesting modifications of w and to the values in the CDM model (see Sect. 3). On the other hand each sample with ξ < 0 eventually converges to the result of CDM.
We fit the model parameters by virtue of the Markov-chain Monte Carlo approach through the cosmomc program [59]. The background and linear perturbation equations given in Sect. 2 and Appendix Appendix A: are numerically solved by the mgcamb package [60,61], where formalisms in both the conformal Newtonian gauge and the synchronous gauge are used. We have included both the scalar and the tensor mode perturbations in our code, as considered in [62,63].
We first fix the value of ξ to test the allowed region of 0 by using the combined data of PLANCK + BAO + MPK with supernova (SNIa) and weak lensing (WL). We find 0 < {0, 016, 0.051, 0.085} and 0 < {0, 035, 0.114, 0.201} at the 68 and 95% confidence levels for ξ = {0.1, 0.3, 0.5}, as shown in Table 1. Results of the main model parameters with fixed ξ values are given in Fig. 3 and Table 2. We note that the increase of the ξ value tends to suppress the active neutrino mass sum, as seen by Table 1. We then treat ξ as a free parameter and perform the fitting with respect to various combinations of data sets. The bestfit value for 0 is less than 0.2 in the first two sets of results based on PLANCK + BAO and their combination with MPK. The best-fit value becomes less than 0.1 when the SNIa and WL data are included in sequence, as shown in Fig. 4. In Table 1 Upper bounds of the model parameters ε 0 and m ν in the constant potential model with the fixed non-minimal coupling parameter ξ by using the combined PLANCK, BAO, MPK, SNIa, and WL data particular, we find 0 < {0.223, 0.456} at the 68 and 95% confidence levels from the PLANCK + BAO + MPK data, and the constraints are narrow down to 0 < {0.077, 0.180} by adding SNIa. Given that 0 = 0 is located inside the 1σ limit in all four sets of results, we conclude that there is no hint for the existence of the additional density nc φ , or namely φ = V today. Fig. 3 Contours of the model parameters in the constant potential model with the fixed non-minimal coupling parameter ξ by using combined PLANCK, BAO, MPK, SNIa, and WL data, where the prior of the model parameters is given in Table 1  Table 2 Constraints of the model parameters at the 95% confidence level in the constant potential model with the fixed non-minimal coupling parameter ξ by using the combined PLANCK, BAO, MPK, SNIa, and WL data We now discuss the constraint on the non-minimal coupling parameter ξ and the today's value of , where the full data combinations are listed in Table 3 For example, taking 0 = 0.989 we find the upper bounds φ 0 < {0.19, 0.11} for ξ > {0.3, 1.0}. Let us consider more specifically that the extended quintessence model remains a valid approximation on very small scales, such as inside our Solar system. In this case we can use ω 0 > 40,000 as an independent constraint, where Eq. (4) leads to φ * = φ 0 e L 2 N * , where L 2 is given by Eq (41). Given that the background curvature is decreasing with the increase of , we may take 3H 2 ≈ |L 2 2 /2 − 3ξ − 6ξ L 2 |H 2 φ 2 * to solve N * , which reads Setting ξ = 0.3, we find N * > 2.0 for φ 0 < 0.19 from the cosmological tests and N * > 5.4 for φ 0 < 8.3 × 10 −3 from the local constraint. Therefore the onset of the ξ -induced power-law expansion could occur as early as some 27.6 (74.5) billion years later from the present for N * > 2.0 (N * > 5.4). For the time unit, we have used H 0 dt = d ln a = d N with H −1 0 ≈ 13.8 × 10 9 years.

Summary and discussions
We have considered one of the simplest extensions to the CDM model based on a quintessence field with a nonminimal coupling ξ to gravity. Our principal goal has been to identify the compatibility of such gravitational effect in current observations. We have in particular focused on the class of models with very weak restrictions on the today's field value φ 0 so that both ξ and φ 0 can be treated as free model parameters. This condition is suitable for the slow-roll type model V (φ) = V 0 e −λκφ , where the potential energy is insensitive to the field value φ with a sufficiently small λ. In the absolutely flat limit where λ = 0, the theory coincides with the α = 0 limit of the tracker model V (φ) = V 0 (κφ) −α , which has a completely vanishing tracking basin of attraction. However, in general cases the tracker models with α > 0 only realize the cosmic acceleration when φ 0 ∼ 1/κ, provided that V 0 is responsible for the energy scale of the dark energy domination [12,13,22].
We have derived analytic solutions for the background evolution in exponential models, including the possible future time dynamics beyond the de Sitter expansion. For ξ < 0, we have found that the asymptotic attractor is always the de Sitter like solution Eq. (51) with the constant Hubble parameter H f , the coupling constant ξ , and the final field value φ f . In the special case where λ = 0, Eq. (51) leads to φ f = 0 and H f = κ √ V 0 /3. For ξ > 0, we have shown that the current potential-driven expansion becomes a ξ -induced expansion with w = (−3 + 2ξ)/[3(1 + 2ξ)]. This transition is due to the fine cancelation between the potential energy and the negative density component corresponding to the non-minimal coupling (namely ρ nc φ ). In the λ = 0 case, the ξ -induced power-law expansion can occur even if ξ > 3/2, provided that ρ φ is positive before the matter domination. However, in this case the solution w = (−3 + 2ξ)/[3(1 + 2ξ)] is a saddle point and the final  Table 1 attractor reads w = 4ξ/[3(1 + 4ξ)]. In general cases where λ > 0, the final attractor can be a scaling solution in which |ρ φ | tracks closely with ρ M . In addition, we have numerically checked that the total density of the universe is always positive even if |ξ | is very large.
To probe the existence of ξ , we have taken the combined data from type-Ia supernovae (Union 2.1), cosmic microwave background (Planck 2015), baryon acoustic oscillations and matter power spectrum (6dF Galaxy Survey, BOSS and Wig-gleZ), to the weak lensing (CFHTLenS) measurements. With the optimal choice λ = 0 of the model, the background evolution is identical to the fiducial CDM scenario and the number of the model parameters is the same as that of the PBD research [40]. We have confimed that the most preferable value of 0 is zero so that there is no hint for the dark energy density additional to that of a pure cosmological con- stant. This result is true whether we fix a constant value to ξ or treat it as a free parameter. Nevertheless, we have demonstrated that the current observations prefer a non-zero value ξ > 0.289 for a variable ξ , where ξ = 0 is excluded outside the 95% confidence region. This is in contrast to the extended tracker model [22] where a strong constraint |ξ | 1 was found. Meanwhile, the central value of the effective gravitational constant = G eff /G today is found to be slightly less than unity as 0 = 0.989 +0.014 −0.023 at 95% level, consistent with the hint of a positive ξ . It is noteworthy that the best-fit sample with a variable ξ parameter brings a slightly better χ 2 results than those with fixed values of ξ . However, a ξ > 0 prior always leads to marginally larger χ 2 than the CDM (ξ = 0) case.
It is noteworthy that the current observations [4] (see also [73]) on the dark energy equation of state indicate a nearly 2σ preference in the phantom domain, i.e. w DE < −1. This result implies the importance to probe the constraint on ξ with a slightly phantom-type background expansion, as a straightforward extension to our current study. Note that the non-minimal theory given by (2) cannot cross the phantom divide due to its conformal equivalence to the coupled quintessence model [52]. In addition, it has been found that the gravitational non-minimal coupling ξ can act as a cosmological attractor for not only the early-time [74][75][76][77][78][79] but also for the late-time [80] cosmic accelerations. In this class of models one may go beyond the slow-roll conditions for V (φ), and one may obtain a significant density fraction | nc φ / V | ∼ O(1) today. We also emphasize that the working frame used in this study is based on the matter point of view [81,82] (see also [83,84]) in which both matter and radiation are minimally coupled with gravity and the dark energy field. This physical frame appears to be the most natural choice when performing the fitting with the cosmological priors given in Table 1, whereas the possible distinction on observables in different conformal frames is an issue goes beyond the scope of the current paper.
Finally, we remark that if a future measurement would confirm 1 − 0 > 0, it can support the existence of a positive coupling ξ but will falsify at the same time the contribution of many types of screening mechanisms near the Solar System [72]. On the other hand, a future measurement with 1− 0 < 0 (as already hinted in [56]) will not be in favor of a negative ξ but instead can rule out the whole class of massless extended quintessence models, given that ξ < 0 only leads to the decay of the field value so that always converges to the CDM limit where 0 = 1.
We decompose the scalar field into a homogeneous part and a perturbation one as φ(x) = φ(τ ) + δφ(x). Let us put the perturbed metric, the perturbed energy-momentum tensor and the perturbed φ in Eq. (5) to obtain the perturbed field equations as (κ 2 eff = κ 2 ) where in this section primes are derivatives with respect to τ , and Similarly, the perturbed Klein-Gordon equation is derived as δφ + 2H δφ − 6ξ(H + H 2 )δφ − ∂ 2 δφ + a 2 V φφ δφ + 2ξφ[6(H + H 2 ) + 9H + 3H + 3 The tensor mode perturbation given by the spatial part of the modified Einstein equation (5) is simply D i j + 2 H + ξ φφ D i j − ∂ 2 D i j = κ 2 a 2 π i j , (A. 11) where π i j is assumed to be vanishing in this work.