Constraining the attractive fifth force in the general free scalar-tensor gravity with solar system experiments

In this paper, we focus on the general free scalar-tensor gravity with three free coupling functions, which in the near-field region looks like general relativity (GR) plus a fifth force of Yukawa-type induced by the scalar field. We show that the fifth force is always attractive in the theory. We investigate the effects of the attractive fifth force and calculate in detail the fifth force-induced orbital precession rate $\delta\omega/\omega$ and the parameterized post-Newtonian parameters $\gamma$ and $\beta$, all of which depend on the fifth force parameters and the interaction distance. It turns out that, due to the attractive fifth force, $\delta\omega/\omega$ is always greater than zero, $\gamma$ is always less than one, $\beta$ is greater than one at large distances, and additionally this class of theories is ruled out as an alternative theory to dark matter. We place stringent constraints on the fifth force parameters by combining the lunar laser ranging (LLR), Cassini, and Mercury precession experiments, and derive the upper bounds on the strength ratio of the fifth force to gravitational force at different scales from the LLR observation. We find that the Mercury constraint is not competitive with the LLR and Cassini constraints and the LLR observation imposes much more stringent bounds on the strength ratio on large scales than on small scales. Our results show that this theory is sufficiently close to GR for a small enough fifth force strength and can reduce to GR with a minimally coupled scalar field in the absence of fifth force.


Introduction
Although Einstein's general relativity (GR) has been very successful at interpreting gravity, it has been plagued by the problems of quantization [1,2] as well as dark matter and dark energy [3,4].It is commonly believed by scientists that GR is not the final theory of gravitational interaction.Therefore, studies of alternative theories of gravity play a crucial role in testing GR [5,6].It is now known that GR is the unique interacting theory of a Lorentz invariant massless spin-2 particle [7,8].Based on this, a natural way to extend GR is to add extra degrees of freedom.Scalar fields are widely used in the fundamental physics [9] and in cosmology [10,11,12,13], and additionally the discovery of Higgs boson [14] has shown that scalar particles really exist in nature.Based on these considerations, scalar degrees of freedom are introduced into most alternative theories of gravity, and the simplest one is the scalar-tensor gravity [15,16], which can generally generate a fifth force [17,18] to mediate gravitational interactions.
Scalar-tensor theories of gravity have been extensively studied on a huge range of scales, from laboratory tests a e-mail: zhxing@nwu.edu.cn[17,18,19], to solar system [20,21,22,23,24,25,26] and binary pulsar [27,28,29,30] tests, to galactic [31,32,33,34] and cosmological [35,36,37] tests.However, the literature cited above only studied the scalar-tensor gravity with one or two coupling functions, which is not the most general case.In this paper, we focus on the general free scalartensor gravity, which contains three free coupling functions, namely a free coupling function controlling the strength of the interaction, a free scalar field kinetic coupling function, and a free scalar potential function.
We first consider the simple case of a vanishing scalar potential.In the case, for a perfect fluid as a matter source, we derive the 2 post-Newtonian (PN), 3PN and 4PN equations for the massless tensor and massless scalar fields, respectively, and calculate in detail all ten parameterized post-Newtonian (PPN) parameters by solving these PN equations.We find that in the case, the fifth force of Yukawa-type is absent, the PPN parameters γ and β are model-dependent constants, and the remaining eight PPN parameters are all zero.
In particular, we focus on the theory of a non-vanishing free scalar potential.In the theory, by solving the 2PN and 4PN equations of the massless tensor and massive arXiv:2305.06752v2[gr-qc] 10 Apr 2024 scalar fields for a point-like matter source, we derive the effective gravitational constant G eff and the PPN parameters γ and β.We find that in the theory the fifth force of Yukawa-type appears and the three parameters mentioned above are no longer constants but depend on the interaction distance from the test mass to the gravitational source.The fifth force affects the orbital dynamics of the system, slightly breaking the inverse-square law of Newtonian gravity and thus inducing the perihelion precession of the orbit.We derive the perihelion precession induced by any deviation from the inverse-square law by starting from the Lagrangian of the orbital motion and using perturbation analyses.We obtain the fifth force-induced perihelion precession by reducing any deviation to the Yukawa-type fifth force.We find that the coupling strength α and mediator mass m 0 of the Yukawa-type fifth force can be represented as the certain combinations of three free functions of the theory and can be used to more concisely re-express the PPN parameters γ and β.By analyzing the functional form of the coupling strength α, we find that in the theory the fifth force of Yukawa-type is surprisingly always attractive.This leads to the following interesting conclusions for the general free scalar-tensor gravity.First, this theory (with an attractive Yukawa-type fifth force) is insufficient to explain the flat rotation curves of galaxies without the need for dark matter (building on previous works in [38,39,40,41]).Second, the fifth force-induced perihelion precession is always in the same direction as the orbital motion.Third, the PPN parameter γ is always less than one, while β is greater than one at large distances, and they both tend to 1 in the limit of large distances.In short, although not strictly true, this is a good physical picture: in the near-field region, the general free scalar-tensor gravity looks like GR plus an attractive fifth force of Yukawa-type.
Finally, we place stringent constraints on the fifth force parameter space (α, m 0 ) by combining the experiments of lunar laser ranging (LLR) [42], Cassini [43], and Mercury perihelion precession [44].We find that all values of the coupling strength are allowed when the mediator mass is greater than a certain value that is different for different experimental constraints.The constraint from Mercury precession is the weakest one of all.As the mediator mass decreases, the Cassini constraint becomes tighter and tighter and tends to a constant.For small mediator masses the LLR constraint is not competitive with the Cassini constraint, but as the mediator mass increases the LLR constraint becomes tighter and tighter until it reaches the tightest constraint.Furthermore, from the LLR observation, we derive the upper bounds on the fifth force to gravitational force strength ratio f 5 (r) at different scales.It is an interesting result that the upper bounds on the strength ratio f 5 (r) become tighter and tighter as the interaction distance increases.
The organization of this paper is as follows.In Sec. 2, we display the action and field equations for the general free scalar-tensor gravity.In Sec. 3, we calculate in detail all ten PPN parameters for the theory in the simple case of a vanishing scalar potential.In Sec. 4, we focus on the theory with a non-vanishing free scalar potential and calculate in detail the effective gravitational constant and the PPN parameters γ and β for the theory.In Sec. 5, we study the Yukawa-type fifth force in the theory and derive in detail the fifth force-induced perihelion precession.In Sec. 6, we discuss some interesting properties of the general free scalar-tensor gravity.In Sec. 7, we place stringent constraints on the parameter space of the attractive fifth force by solar system experiments and derive the upper bounds on the scalar to tensor force strength ratio at different distances from the LLR observation.We conclude in Sec. 8.

General free scalar-tensor gravity
A general free scalar-tensor gravity with three free coupling functions is given by the following action [16,15]: where G is the bare gravitational constant, g is the determinant of the metric g µν , R is the Ricci scalar, and S m is the action of matter fields ψ m .F (ϕ) is the coupling function that controls the local value of the gravitational coupling constant, W (ϕ) is the coupling function that describes the scalar field kinetic interaction, and V (ϕ) is the scalar field potential function that endows the scalar field with mass.Here, the three coupling functions F (ϕ), W (ϕ) and V (ϕ) are all free functions, and they or their combinations are not constrained by theoretical considerations.The variation of the action (1) with respect to the metric and scalar fields yields the dynamical equations of g µν and ϕ: where a prime ( ′ ) denotes differentiation with respect to ϕ, T µν ≡ (−2/ √ −g)δS m /δg µν is the energy-momentum tensor of matter fields, and T is its the trace.
3 All ten PPN parameters for V = 0 In order to distinguish different theories of gravity and to compare their predictions with a large number of experiments, the PPN formalism has been developed [45,46,47,48,49,50,44].In the PPN formalism, the PN metrics for gravity theories are parameterized by a set of PPN parameters, which have been measured by various high precision experiments [44].
In this section, we use the PPN formalism to derive in detail all ten PPN parameters of the general free scalartensor gravity with a vanishing scalar potential for a perfect fluid as a matter source.In the PPN formalism, the gravitational field of the matter source is weak (i.e., GM/r ≪ 1) and its typical velocity v is small (i.e., v 2 ∼ GM/r ≪ 1), and thus the PN metrics can be calculated by expanding the field equations to O(v n ) and solving them.
In the weak-field, the tensor and scalar fields can be expanded around a flat Minkowski metric η µν = diag(−1, 1, 1, 1) and a scalar background ϕ 0 as follows: where the perturbations h µν and φ can be expanded to O(v n ) and are denoted by superscript (n) below.For a perfect fluid as a matter source, its energy-momentum tensor is given by [50,44] where u µ is the fluid four-velocity, ρ is the rest mass density of the fluid, p is the fluid pressure, and Π is the internal energy per unit rest mass.
In order to solve the field equations in the PPN formalism, we use the PN gauge introduced in [50,44], By using Eqs.( 4), ( 5) and ( 6) and expanding the dynamical equations ( 2) and (3) (set V = 0) to O(v 2 ), we can obtain the 2PN equations of the fields φ, h 00 and h ij in the form (2) (2) where the subscript '0' indicates the value of a quantity at ϕ 0 .These are Poisson's equations and the solutions are (2) (2) where U is the Newtonian potential (see Eqs. ( 61)).Here, the quantity can be interpreted as the Newtonian gravitational constant and is set to one in the subsequent calculations.From the above equation, we see that the fifth force mediated by the scalar field is absent in the general free scalar-tensor gravity with a vanishing scalar potential.
In the PN gauge (6), expanding up to O(v 3 ), the timespace component of the metric field equation ( 2) is expressed in the form (2) using Eqs.( 11) and ( 68) the solution is (3) where V j and W j are the metric potentials, defined in Eqs.(61).
Using the PN gauge in Eq. ( 6) and T µν in Eq. ( 5), and expanding the time component of the metric field equation (2) to O(v 4 ), we obtain the 4PN equations of the field h 00 , ∇ φ) Substituting Eqs. ( 10), ( 11) and ( 12) into the above equation and using the relations (69) and (70), we obtain the solution in the form (4) where Φ 1 , Φ 2 , Φ 3 and Φ 4 are the metric potentials, defined in Eqs.(61).By comparing the metric perturbations in Eqs. ( 11), ( 12), ( 15) and ( 17) with the standard form of the PPN metric in Eqs.(60), we can identify all ten PPN parame-ters in the form where the significances of these parameters are given in Appendix A.
As we can see, in the case of a vanishing scalar potential, the PPN parameters γ and β are model-dependent constants, and the remaining PPN parameters are all zero.However, in the theory with a non-vanishing free scalar potential, due to the emergence of the Yukawa-type fifth force, these two parameters γ and β are no longer constants but distance-dependent, as we will discuss in the next section.
4 Effective gravitational constant and PPN parameters γ and β for V ̸ = 0 Now let us focus on the general free scalar-tensor gravity with three free functions in the case of a point mass M as a matter source, which can be described by the energymomentum tensor (5) with ρ = M δ(r) and vanishing v i , p and Π .

Effective gravitational constant G eff and PPN parameter γ
In the PPN formalism, imposing the PN gauge (6) and expanding the dynamical equations ( 2) and (3) to O(v 2 ), we obtain the 2PN equations of the fields φ, h 00 and h ij in the form (2) ∇ 2 (2) where the quantity m s can be interpreted as the effective mass of the scalar field, and given by Here, the first is a screened Poisson equation, and the last two can be simplified as the Poisson equations by using the first.By setting ρ = M δ(r) and the solutions of these equations are given by (2) (2) The effective gravitational constant G eff and the PPN parameter γ are defined in the form [44], (2) Using these, from Eqs. ( 26) and ( 27), we obtain where G eff and γ both are distance-dependent and can reduce to G/F 0 and 1 in the limit of m s r → ∞ and to Eqs. ( 13) and ( 18) in the limit of m s r → 0. It is clear that the Yukawa exponential factor has been restored since the scalar field is massive.This means that the Yukawa-type fifth force appears in the theory, as we will see in more detail in the next section.

PPN parameter β
By imposing the PN gauge (6) and the scalar field equation (3) is expanded up to O(v 4 ) in the form φ ,00 + 8πGρ φ , (31) where φ ,00 can be dropped for static solutions.The pressure p and specific internal energy Π can be dropped for a point-like matter source.The term ρ φ corresponds to gravitational self-energies and can be neglected.Using (2) φ and (2) h ij in Eqs.(25) and (27), the solution of the remaining equation is given by ×[e 3msr Ei(−3m s r)−e msr Ei(−m s r)] +2 [e 3msr Ei(−2m s r)−e msr ln(m s r)], (32) where Ei(x) is the exponential integral function, defined by Ei In the PN gauge (6), expanding the metric field equation (2) to O(v 4 ), we obtain the 4PN equation of the field h 00 in the form h 00 φ .(34) Here, the velocity v, pressure p and specific internal energy Π can be dropped, since the calculation is performed in the rest frame of the point-like matter source.The terms ρ h 00 and ρ φ can be neglected, since they correspond to gravitational self-energies and do not affect the calculation of the PPN parameter β.Using φ in Eqs. ( 25), ( 26), ( 27) and (32), the remaining equation is solved by Using the relationship h 00 = −2β(G eff M /r) 2 and G eff in Eq. ( 29), the PPN parameter β is given by where this parameter is also distance-dependent and can reduce to 1 and Eq. ( 19) in the cases of m s r → ∞ and m s r → 0, respectively.
We will see in the next two sections that these parameters G eff , γ and β in Eqs. ( 29), ( 30) and ( 36) can be expressed more concisely as functions of the fifth force parameters.

Yukawa-type fifth force and fifth force-induced perihelion precession
In this section we study the Yukawa-type fifth force and the fifth force-induced perihelion precession in the general free scalar-tensor gravity.

Yukawa-type fifth force
In order to more clearly describe the fifth force in the theory, by using the combinations of three functions F , W and V , we define the following two quantities: As we will see later, these two quantities are the coupling strength and mediator mass of the Yukawa-type fifth force.
Here, F 0 > 0 is a consequence of the fact that gravity is attractive (see Eq. ( 29) or ( 43)), and V ′′ 0 ≥ 0 is the non-negative mass squared of the scalar field.Additionally, the effective mass in Eq. ( 24) is non-negative, i.e., F 0 V ′′ 0 /(3F ′ 2 + 2F W ) 0 ≥ 0. These above relations imply Making use of the definitions of α and m 0 in Eqs. ( 37) and (38), the effective mass in Eq. ( 24) is expressed as and the effective gravitational constant in Eq. ( 29) is rewritten as Using this, the potential energy describing the gravitational interaction between two point-like objects with masses m and M is given by Differentiating this above, the force is where the first term is the gravitational force (tensor force) and the second term is the Yukawa-type fifth force (scalar force).

Fifth force-induced perihelion precession
Because of the presence of the Yukawa-type fifth force, the net force experienced by a planet does not vary exactly as inverse-square, which induces the perihelion precession of the planet orbit in which the major axis of the planet orbit slowly rotates in the orbital plane.In the weak-field limit, the Lagrangian describing the orbital dynamics of two point-like objects with masses m and M (m << M ) in the polar coordinates (r, θ) is given by where an overdot denotes differentiation with respect to time.By varying the Lagrangian and performing the replacement r = u −1 , the equation of motion is where l = r 2 θ is the (conserved) orbital angular momentum per unit mass.The force F (r) ≡ −dV (r)/dr can usually be rewritten in the following form where f (r) is the sum of various deviations from the inversesquare law of Newtonian gravity.Substituting the above equation into Eq.( 45) yields where p = l 2 /(G N M ) is the semi-latus rectum of the unperturbed orbit and is linked to the the semi-major axis a and eccentricity e of the unperturbed orbit by p = a(1 − e 2 ).In the above equation, f (u −1 ) can be regarded as a small perturbation, then the solution should lie near the unperturbed value.For this reason, expanding the right hand side of Eq. ( 47) in Taylor's series at p −1 , we obtain the solution in the form with where u e = eu p for a quasi-elliptical orbit.From Eq. (48) we see that the perihelion occurs at ω(θ − θ 0 ) ≡ ωθ n = 2πn , n ∈ Z.Using this, the precession angle per orbit, δθ ≡ (θ n+1 − θ n ) − 2π is given by Substituting Eq. ( 49) into the above equation and keeping only the leading order, we obtain the perihelion precession rate Note that the above relationship is also obeyed between each component of f (r) and its induced precession.By comparing Eq. ( 43) and Eq. ( 46), we obtain the strength ratio of the fifth force to gravitational force, and this is also the strength ratio of the scalar to tensor force in the general free scalar-tensor gravity.Substituting this into Eq.( 52), the fifth force-induced perihelion precession rate is given by Here, we note that in the presence of the Yukawa-type fifth force, as m 0 → 0 (see Eq. ( 40)), the induced precession goes to zero.

Properties of general free scalar-tensor gravity
In this section we summarize and discuss some interesting properties of the general free scalar-tensor gravity.Making use of the fifth force parameters α and m 0 in Eqs.(37) and (38), the PPN parameter γ in Eq. ( 30) can be rewritten as and when m s r > 1 the PPN parameter β in Eq. ( 36) can be approximated as We now get the four important Eqs.( 53), ( 54), ( 55) and (56) in this paper.Our main results are also summarized in these equations.According to Eq. ( 39), the coupling strength α is always non-negative.Now let us consider the two cases of α = 0 and α > 0.
In the case of α = 0 (i.e., F ′ = 0), the coupling function F is a constant, the non-minimal coupling disappears, and we have the following results: (i) f 5 (r) = 0 and δω 5 /ω = 0, i.e., the Yukawa-type fifth force and its induced precession are absent; (ii) the PPN parameters γ = β = 1.These indicate that the theory in the case reduces to GR with a minimally coupled scalar field.
In the case of α > 0 (i.e., F ′ ̸ = 0), the presence of the non-minimal coupling makes the theory non-trivial, and we find the following interesting results for the general free scalar-tensor gravity: (i) f 5 (r) > 0, i.e., the Yukawa-type fifth force in the theory is always attractive; (ii) δω 5 /ω > 0, i.e., the precession is in the same direction as the orbital motion; (iii) γ < 1, and β > 1 for m s r > 1.
Note that, according to the previous works 1 , the result (i) means that the general free scalar-tensor gravity is not sufficient to reproduce the flat rotation curves of galaxies without the need for dark matter, that is, the theory cannot be used as an alternative theory to dark matter.The results (ii) and (iii) show that the experimental constraints on the theory come only from the upper bound on δω/ω, the lower bound on γ and the upper bound on β.

Solar system constraints
In this section we place constraints on the parameter space of the attractive fifth force in the theory by solar system experiments involving the LLR measurement [42], the Cassini satellite [43], and the perihelion precession of Mercury [44].
The perihelion precession is one of the four classic solar system tests of GR.Bounds on anomalous precession deviating from GR in the solar system provide extremely stringent constraints on the fifth force, in which the most stringent bound comes from the LLR observation [42] searching for the anomalous precession of the Moon.The Moon moves on an elliptical orbit with a small eccentricity of 0.0549 and a semi-major axis of 0.00257 astronomical unit (AU ≃ 1.496 × 10 8 km).In the Earth-Moon system, the characteristic interaction distance is the semilatus rectum of the orbit, i.e., 0.00256AU.The observed quantities relevant to placing constraints are presented in Table 1.
The PPN parameters have been measured accurately by a large number of experiments [44].The current most stringent constraints on the PPN parameters γ and β come from the Cassini measurements of the Shapiro timedelay [43] and the perihelion precession of Mercury [44], respectively.These two experiments also provide a characteristic interaction distance, respectively.In the Cassini measurements [43], the radio signals sent between the Cassini spacecraft and the Earth were passing by the Sun at a distance of 0.00744AU, which is the characteristic interaction distance in the experiment.In the perihelion precession of Mercury [44], Mercury moves on an elliptical orbit with a moderate eccentricity of 0.2056 and a semi-major axis of 0.3871AU, and the characteristic interaction distance is the semi-latus rectum, i.e., 0.3707AU.The measured values of the relevant parameters are presented in Table 1.
From Table 1, we obtain the upper bounds δω ω ≤ 0.5 × 10 −11 and δω ω ≤ 1.21 × 10 −11 (57) 1 Refs.[38,39,40,41] have shown that an almost constant profile of galactic rotation curve is recovered by adding a repulsive Yukawa-type fifth force to gravitational force.55) and (56), we derive the observationally allowed range of the coupling strength α for each value of the mediator mass m 0 , shown in Fig. 1.Using Eqs. ( 53) and ( 54), the upper bounds on δω/ω in Eq. ( 57) can be translated into the upper bounds on f 5 (r) for each value of the effective mass m s , shown in Fig. 2. Note that, here, the upper bound of γ and the lower bounds of δω/ω and β cannot be used to constrain the theory, because of the conclusions (ii) and (iii).In Fig. 1, the yellow, cyan and green regions are excluded by the LLR measurement, the Cassini experiment,   1) are the interaction distances of the LLR, Cassini, and Mercury experiments, respectively.In the yellow regions, the dashed and solid lines are obtained by equating the 1σ and 2σ limits of Eqs.(57) with Eq. ( 54), respectively, and plotting the corresponding contour in the parameter space (α, m 0 ).The dashed and solid lines of the cyan regions are obtained in a similar way using Eqs.( 58) and (55).In the green regions, the dashed and solid lines are obtained in a similar way using Eqs.( 59) and ( 56), and the upper right regions of the dotted line indicate m s d i > 1 (d i = 0.3707AU), in which Eq. ( 56) holds.Clearly, the Mercury constraint on the parameter space (α, m 0 ) is the weakest one of all.Table 2 records several special points in Fig. 1.From Table 2 we obtain the upper bound of α ≤ 1.25 × 10 −5 at 2σ CL for m 0 = 0.This upper bound on α can be translated into the familiar bound ω ≥ 40000 in massless Brans-Dicke theory [15,44] by setting F (ϕ) = ϕ and W (ϕ) = ω(ϕ)/ϕ and substituting them into Eq.(37).We can see that at 2σ CL all values of α are allowed for m 0 ≥ 10 −9.65 eV in the LLR constraint, for m 0 ≥ 10 −13.4 eV in the Cassini constraint, and for m 0 ≥ 10 −15.4 eV in the Mercury constraint.At 2σ CL, the horizontal coordinate of the intersection between δω 5 /ω and γ in Eqs.(54) and ( 55) is m 0 = 10 −15.70 eV.We can see that in the current precision the range of values of α for m 0 ≤ 10 −15.70 eV is exactly the same as that of m 0 = 0.In other words, the upper bound on α is no longer sensitive to values of m 0 , when m 0 ≤ 10 −15.70 eV or m −1 0 ≥ 10 9 m, which is comparable to the interaction distance of 0.00744AU in Table 1.
Fig. 2 shows the constraints on the strength ratio f 5 (r) of the fifth force to gravitational force from the LLR observation searching for the anomalous precession of the Moon.The curves represent the 2σ CL upper bounds on the strength ratio f 5 (r) as a function of effective mass m s at different distances r/d M = 100(red), 10(orange), 2 (yellow), 1 (green), 0.99 (cyan), 0.9 (blue), 0.2 (purple), and ∈ [0, 0.1] (grey).Here, the effective mass m s is defined in Eq. (40), and d M = 0.00256AU is the semi-latus rectum of the Moon's orbit (see Table 1).We find that the 2σ CL upper bounds of f 5 (r ≤ 0.1d M ) are approximately the same as that of f 5 (0) (see the grey curve).As the distance r decreases, the excluded parameter space decreases, and then the constraints on the strength ratio f 5 (r) become  weak.Therefore, the excluded parameter space of r = 0 is located entirely in the excluded parameter space of r > 0. The LLR observation places much tighter fifth force constraint on large scales than on small scales.

Conclusions
General free scalar-tensor gravity with three free coupling functions is a natural extension to GR.One of the most important consequences of this extension to GR is the emergence of a Yukawa-type fifth force.
In this paper, we explored the impact of the Yukawatype fifth force in the theory on the PN metric and the orbital precession and derived in detail the PPN parameters γ and β and the fifth force-induced perihelion precession rate δω/ω, respectively.These parameters depend not only on the interaction distance from the test mass to the gravitational source, but also on the coupling strength and mediator mass of the Yukawa-type fifth force.We found that if the fifth force vanishes then all these parameters reduce to their values in GR.Our further analysis showed that the Yukawa-type fifth force in the theory is surprisingly always attractive, which leads to δω/ω > 0, γ < 1, and at large distances β > 1.In addition, this class of theories cannot completely substitute the role of dark matter in describing galactic dynamics.
Finally, we used solar system experiments (LLR, Cassini, and Mercury precession) to place stringent constraints on the parameter space of the attractive fifth force, and derived the upper bounds on the fifth force to gravitational force (scalar to tensor force) strength ratio at different scales with the LLR observation.It turned out that the Mercury constraint is the weakest one of all, for small mediator masses the LLR constraint is not competitive with the Cassini constraint, but as the mediator mass increases the LLR constraint becomes tighter and tighter until it reaches the tightest constraint.Furthermore, we found that the upper bounds on the scalar to tensor force strength ratio become tighter as its interaction distance increases.
In this paper, all coupling functions are free, but the real physical world is usually governed by symmetry constraints.So, in future work, we will study the general unfree scalar-tensor gravity constrained by theoretical considerations, which gives rise to extremely interesting results.

A PPN formalism
The PPN formalism was inspired by the earlier work of Eddington [45], Robertson [46] and Schiff [47].Nordtvedt [48] and Will [49] generalized the framework to a system of gravitating point masses and a perfect fluid, respectively.
In the PPN formalism, the metric for a perfect fluid is parametrized as follows [50,44]: ), where w i is coordinate velocity of the PPN coordinate system relative to the mean rest-frame of the universe.The coefficients γ, β, ξ, α 1 , α 2 , α The metric potentials U , U ij , Φ W , A, Φ 1 , Φ 2 , Φ 3 , Φ 4 , V i and W i are constructed as follows [50,44]: where ρ, p and Π are respectively the rest mass density, pressure and specific internal energy of the perfect fluid (see the energy-momentum tensor of the perfect fluid in Eq. ( 5)).The metric potentials satisfy the following simple and useful relationships [50,44]

Fig. 1 .
Fig. 1.Solar system constraints on the coupling strength α of the fifth force as a function of its mediator mass m0.The shaded regions are excluded by various experiments.The dashed and solid lines represent 1σ CL and 2σ CL, respectively.
α ≤ 10 −10.4 m 0 (eV) ∈ [0, +∞) LLR and the Mercury anomalous precession, respectively.The three parallel dotted lines from right to left indicate the curves m s d i = 1 where d i = 0.00256AU, 0.00744AU, and 0.3707AU (see Table

Fig. 2 .
Fig. 2. 2σ CL upper bounds from the LLR observation on the strength ratio f5(r) at different distances as a function of effective mass ms.

3 , ζ 1 , ζ 2 ,
ζ 3 and ζ 4 are the PPN parameters, whose significances are as follows: γ measures the amount of space-curvature produced by unit rest mass, β measures the non-linearity in the superposition law for gravity, ξ measures the violation of local position invariance, α 1 , α 2 and α 3 measure the violation of local Lorentz invariance, and α 3 , ζ 1 , ζ 2 , ζ 3 and ζ 4 measure the violation of total energy-momentum conservation.

Table 1 .
Current experimental constraints on the Moon perihelion precession rate δω/ω and the PPN parameters γ and β.

Table 2 .
2σ CL bounds on the coupling strength α of the fifth force for different values of its mediator mass m0.