Observational tests of asymptotically flat ${\cal R}^{2}$ spacetimes

A novel class of Buchdahl-inspired metrics with closed-form expressions was recently obtained based on Buchdahl's seminal work on searching for static, spherically symmetric metrics in ${\cal R}^{2}$ gravity in vacuo. Buchdahl-inspired spacetimes provide an interesting framework for testing predictions of ${\cal R}^{2}$ gravity models against observations. To test these Buchdahl-inspired spacetimes, we consider observational constraints imposed on the deviation parameter, which characterizes the deviation of the asymptotically flat Buchdahl-inspired metric from the Schwarzschild spacetime. We utilize several recent solar system experiments and observations of the S2 star in the Galactic center and the black hole shadow. By calculating the effects of Buchdahl-inspired spacetimes on astronomical observations both within and outside of the solar system, including the deflection angle of light by the Sun, gravitational time delay, perihelion advance, shadow, and geodetic precession, we determine observational constraints on the corresponding deviation parameters by comparing theoretical predictions with the most recent observations. Among these constraints, we find that the tightest one comes from the Cassini mission's measurement of gravitational time delay.


I. INTRODUCTION
Einstein's theory of general relativity (GR) has been the most successful theory describing the dynamics of massive objects under gravitational effects such as motions of binary stars and planetary motions near their host stars, predicting novel gravitational objects such as black holes, compact stars, and gravitational waves.It remains the only theory of gravity that passes all solar system and astronomical tests.In the last few years, two earliest predictions of GR have been tested and verified via black holes namely: the gravitational deflection of light passing by a black hole resulting in the formation of black hole shadow and the existence of gravitational waves.Recently the shadows of black holes at the center of M87 and Milky Way galaxies have been observed and analyzed in detail [1][2][3][4].Astronomers have also detected numerous gravitational wave signals generated by the merger of binary black holes of different masses [5,6].Despite these successes, some fundamental problems remain in the foundations of GR including its renormalization and establishing its unison with quantum mechanics thereby formulating a set of physical laws valid for both length scales, the very small and the very large.
The accelerated cosmic expansion observed in 1998 has spurred efforts to modify GR to account for the enigmatic "dark energy" component.Among the various modified theories of gravitation, the family of f (R) introduced by Buchdahl in early 1970s has become an active arena of investigation in the past 25 years [7][8][9][10].Within this ghost-free class of theories, pure R 2 gravity stands out for its scale invariance nature.Attempts to incorporate the matter sector, viz., the Glashow-Weinberg-Salam model of particle physics into pure quadratic gravity to form a renormalizable quantum gravity framework have been made in the form of adimensional gravity, or "agravity" [11,12].
can be recast using an auxiliary scalar field Φ as [8] resulting in the field equations in which Φ is equal to the Ricci scalar R. In vacuum, viz.T µν = 0, the terms 1 4 g µν Φ and 1 Φ (∇ µ ∇ ν Φ − g µν Φ) act as additional "matter" sources to the Einstein-Hilbert field equation, producing corrections beyond the vacuum solutions of GR.To date, there is yet empirical evidence of a scalar degree of freedom beside the established tensor ingredients.However, the dynamical nature of Φ, governed by the "harmonic" Equation ( 4), suggests potential manifestations for Φ near mass sources in the strong field regime.As demonstrated by one of the authors [13], the pure R 2 field equation admits a rich host of vacuum solutions with non-constant Ricci scalar.The Buchdahlinspired solutions found therein are asymptotically de Sitter (or anti-de Sitter) and entails a new (Buchdahl) parameter k, accounting for the term 1  Φ (∇ µ ∇ ν Φ − g µν Φ) and enabling spatial variations in Φ. Notably, considerations of a scalar degree of freedom have appeared in other contexts, such as the scalarization in neutron stars and black holes [14].
One important feature of pure R 2 gravity is its limit to GR.In particular, it concerns with the existence of a Newtonian behavior which is an essential requirement for any viable theory of gravitation.The limit is delicate [15].It is because whereas the vacuo of GR is Ricci flat, viz.R µν = 0, for pure R 2 gravity the vacuum background far way from mass sources is de Sitter (or anti-de Sitter) with a Ricci scalar R = 4Λ.In [15] it was found that, owing to the de Sitter background, the spin-2 tensor graviton excitations are massless instead of massive (which would have been the case if the background were locally flat).The massless modes should effectively carry a long-range interaction rather than a Yukawa short-range interaction.Based on this insight, one of the present authors [16] has recently established the emergence of a gravitational potential with the correct Newtonian tail on a de Sitter background for the pure R 2 theory.That is to say, pure R 2 gravity possesses a proper Newtonian limit, despite the absence of the Einstein-Hilbert term in its action (1).This finding strengthens the viability of pure R 2 as a candidate theory of gravitation.
An immediate implication is to find novel static and spherically symmetric spacetime or black hole solutions in R 2 theory which was pioneered by Buchdahl [17].His work culminated in the formulation of a second-order ordinary differential equation for finding metric coefficients which remained unsolved for a long time.Until recently, one of the coauthors (Nguyen) succeeded in obtaining vacuum solutions which are asymptotically de Sitter [13] or asymptotically flat as a special case [18] (with the axisymmetric extensions of the latter solution having been proposed in [19]).Our current study shall be concerned with this special Buchdahl-inspired solution.The metric involves a free parameter, the Buchdahl parameter k, which can be interpreted as a scalar hair, and setting k = 0 yields a Schwarzschild spacetime/black hole as a limiting case.In [18] it was shown that the respective static, spherically symmetric black hole solu-tion has different areas of event horizon depending on the chosen value of k.Specifically, the horizon area can be 0, 4πr 2 s , 16πr 2 s , and divergent for the following values k ∈ (−∞, r s ) ∪ (0, ∞), k = 0, k = −r s and k ∈ (−r s , 0), respectively.Here r s is the radius of the black hole horizon.A further investigation of this new solution is required to understand its phenomenology.
We need to clarify that the special Buchdahl-inspired solution is also Ricci scalar flat, namely, Φ → 0 everywhere in the vacuum exterior to a mass source.Whereas the term − 1 4 g µν Φ in Eq. ( 3) is negligible in this limit, contributions from the term 1 Φ (∇ µ ∇ ν Φ − g µν Φ) persist and are encoded by the (dimensionless) Buchdahl parameter k, defined to be k normalized by r s [20].As a higher-derivative characteristic, the value of k is systemdependent.It does not have a universal value but can vary from one system to another, depending on the composition of the matter source.In normal conditions, such as in the Solar System, k could be insignificant.In extreme conditions, such as around compact stars, intuitively, k may acquire large values.
In this article, we are motivated to analyze the experimental and observational implications of the special Buchdahl-inspired metric.We would like to see how much the new effects beyond the Schwarzschild case affect the dynamics of particles in geodesic motion in the special Buchdahl-inspired metric.To be more general, we consider a more general form of Buchdahl-inspired metric by treating several new parameters for describing the solution independently.By relying on the classical relativistic methods, we calculate perihelion shift, gravitational time delay, and geodesic precession of orbits, and test these results with the solar system experiments.We also attempt to consider the observational implications of the Buchdahl-inspired spacetime using the observational data of S2 star orbit about the Milky Way central black hole and to investigate the shadow of rotating solutions.
Our article is structured as follows.In Sec.II we provide a brief review of the general and special Buchdahlinspired metrics in different sets of coordinate systems.In Sec.III, we investigate the geodesics of both massless and massive objects in the general Buchdahl-inspired spacetime and derive in detail the effects of the Buchdahlinspired spacetime on observations in the Solar System experiments, black hole shadow of M87, and the orbit of S2 star at the Galactic center.The observational bounds on the deviation parameter, which characterizes the deviation of the asymptotically flat Buchdahl-inspired metric from the Schwarzschild spacetime, are obtained by comparing the theoretical predictions with observational data.Then, in Sec.IV, we study a spinning object in the general Buchdahl-inspired spacetime and derive the geodetic procession of its spin vector, from which we obtain the constraints on the corresponding deviation parameter by using the Gravity Probe B and lunar laser ranging data.A brief summary of our main results and some discussions are presented in Sec.V.

II. SUMMARY OF BUCHDAHL-INSPIRED VACUUM SOLUTIONS
In a pioneering Nuovo Cimento work in 1962 [17], Buchdahl developed-but prematurely abandoned-a program to find vacuum configurations for pure R 2 gravity.Subsequent advancements made by one of us, documented in Refs.[13,18], completed his program and derives an exhaustive class of metric in a compact form, to be summarized in this Section.
The field equation in vacuum is which contains fourth derivatives of the metric components g µν in R and ∇ µ ∇ ν R. The solution in general thus involves two additional parameters.As we shall see, they are the scalar curvature, 4Λ, at spatial infinity, and a new (Buchdahl) parameter k which is of the dimension of length.The case of Λ = 0 is an asymptotically (anti-)de Sitter vacuum solution, whereas the case of Λ = 0 is an asymptotically flat vacuum solution.
The body of our orbital motion study in this paper is on the asymptotically flat vacuum solution outside of a static and spherical symmetric mass source.Nevertheless, for completeness, we shall expose the available representations of the metrics in this Section.

A. Asymptotically de Sitter vacuum solution in standard coordinates
In [13], a general Buchdahl-inspired metric was determined to be in a compact form (with dΩ 2 := dθ 2 + sin 2 θ dφ 2 ) The variables p and q obey first-order "evolution" rules whereas the Ricci scalar is The metric involves Λ, representing the scalar curvature at spatial infinity, and k, the Buchdahl parameter.When k = 0, metric (6) duly recovers the de Sitter metric [13].

B. Asymptotically de Sitter vacuum solution in "canonical" coordinates
The metric expressed in (6) contains a conformal factor, which is inversely proportional to the Ricci scalar.In [19], we considered making r a function of R such that the proper part of the metric satisfies g tt g RR = −1.The resulting metric is given by with The "evolution" rules now involve three functions p(R), q(R) and r(R) with the Ricci scalar being given by

C. Asymptotically flat vacuum solution in standard coordinates
In Ref. [18] we further found an exact closed analytical solution corresponding to the case of Λ = 0, which was called the special Buchdahl-inspired metric.This metric is Ricci scalar flat, but not Ricci flat.It describes an asymptotically flat spacetime.Hence, it is also appropriate to call it the asymptotically flat Buchdahl-inspired metric.
For the region r > r s , the metric of the asymptotically flat Buchdahl-inspired metric read where the function ρ(r) is and the dimensionless parameters are

D. Asymptotically flat vacuum solution in "q" coordinate
In [18], one of us reported yet an alternative expression for the asymptotically flat Buchdahl-inspired metric, Eq. ( 16).For the region q q + (which are defined in Eq. ( 20)): with

E. Asymptotically flat vacuum solution in isotropic coordinate
The metric ( 16) can be transformed into an isotropic form [21], which is symmetric with respect to a reciprocal coordinate transformation, per 4r

F. Asymptotically flat vacuum solution in Morris-Thorne form
In Ref. [21], the metric ( 16)-( 18) was brought into the Morris-Thorne form [22,23]: with Φ(R) and b(R) being the redshift and shape functions, respectively.Note that the relation y(R) is implicit by inverting Eq. ( 26).

G. A more generic Morris-Thorne form
In Ref. [24], we generalized the metric in ( 23)-( 27) by making two modifications: (i) replacing k with η in the redshift function (see below); (ii) treating k, η and ζ as independent parameters (in contrast to Eq. ( 27), where ζ = 1 + 3 k2 .The generalized metric is expressed as with Here we would like to present several remarks about the properties of the solution in the generic Morris-Thorne form.
Remark 2. Although the metric ( 28)-( 31) seems to have four parameters k, η, ζ, r s , it effectively depends on only three parameters: with the relation y(R) being implicit in Eq. ( 35).

III. GEODESICS AND CLASSICAL TESTS OF THE GENERAL BUCHDAHL-INSPIRED METRICS
In this section, we will present the geodesics evolution of a massive/massless particle orbiting the Buchdahlinspired metrics.From the geodesic evolution, we are able to calculate some observational quantities so that we can use data to constrain them.We claim that the data from the solar system, experiments can make stronger constraints than those obtained upon using the stellar stars orbiting the supermassive black hole in the galactic center, so we will consider the solar system experiments first.

A. Geodesics in the general Buchdahl-inspired metrics
In this subsection, we give the general Buchdahlinspired metrics in the standard coordinates and consider the geodesics of both massless and massive objects in this metric.Different from the metric given in ( 16), here we treat the parameters k, η, and ζ as independent parameters.Using the coordinates (t, R, θ, φ), the metric is given by In the weak field approximation, this metric can be expressed as The radius r s can be related to the ADM mass of the solution via r s = 2GM with M being the ADM mass and G being the gravitational constant.Comparing this weak field expansion with the Newtonian limit, one can relate G with the Newtonian gravitational constant Let us first consider the evolution of a particle in the general Buchdahl-inspired metric (36).For a massive/massless particle, if one ignores self-gravitational effects, its evolution is governed by the following geodesics where λ denotes the affine parameter of the geodesics and Γ µ νρ represent the Christoffel symbols of the general Buchdahl-inspired metric.Considering the general Buchdahl-inspired metric is static and spherically symmetric, it has two Killing vectors, ξ µ t = {∂ t , 0, 0, 0} and ξ µ φ = {0, 0, ∂ φ , 0}, which leads to two constants of motion E and L (conserved energy and angular momentum), i.e., For geodesics, we also have g µν dx µ dλ dx ν dλ = ε with ε = −1 for timelike geodesics which describes evolution of massive particle and ε = 0 for null geodesics which describes the evolution of massless particle.Then using (39) and (40), one obtains Without generality, we consider the evolution of the particle in the equatorial plane, i.e., we can set θ = π/2 and dθ/dλ = 0. Then one can simplify the above equation into the form where V eff (R) denotes the effective potential of the particle, For later convenience, it is useful to give the derivative of r with respect to φ, which can be obtained by using (40) in Eq. ( 42) and is given by This equation is the starting point for later calculations of the light deflection angle, gravitational time delay, and perihelion advance in the general Buchdahl-inspired metric in the following subsections.

B. Light deflection angle
The precise measurements of the deflection of the light passing by the Sun play an essential role in the establishment of GR.These data can also be used for constraining any possible derivation of the deflection angle in many modified gravities from that in GR.Here our purpose is to calculate the possible effects of the general Buchdahl-inspired metric on the deflection angle of the light and then constrain them by using the most recent measurements.
For the propagation of the light in the general Buchdahl-inspired metric (36), we have ε = 0. Introducing the impact parameter Eq. ( 44) can be transformed into where ± represents the cases with increasing and decreasing R, respectively.In general, for a bending light that does not fall into the object described by the general Buchdahl-inspired metric (36), the range of the allowed R is determined by the condition dR dλ ≥ 0. In the general Buchdahl-inspired metric, this implies the allowed range of R should be R 0 ≤ R < +∞ with R 0 denoting the closest approach of the light to the Sun.R 0 is a root of dR dλ = 0, thus one has Then the deflection of the angle of the light can be calculated by using Since we are considering the deflection of the light by the Sun, it is convenient to employ the weak field approximation by expanding the above integral in terms of r s /R, which gives where G N = (1+η)G is the Newtonian gravitational constant.For a special Buchdahl-inspired metric, as given in ( 16) with η = k and ζ = 1 + 3 k2 , one has where G N = (1 + k)G is the Newtonian gravitational constant for the Buchdahl-inspired metric.Now we consider the light deflected by the Sun.One can express the deflection angle ∆φ in terms of the ∆φ GR = 1.75 ′′ as ∆φ ∆φ The deflection angles of the lights from distance sources by the Sun have been measured in many experiments in the past 100 years.Currently, the most precise measurement was carried out by using the technique of the verylong baseline interferometry [29].By using the result of this measurement, one can constrain the parameter

C. Gravitational Time Delay
Gravitational time delay is an important phenomenon in that lights or radio waves can take more time to travel if they pass by a massive object, like the Sun or a planet.This phenomenon can be precisely measured by sending a radar signal from Earth or a spacecraft passing through the Sun and reflecting off another planet or spacecraft.The effects of the general Buchdahl-inspired metric (36) on the gravitational time delay can be derived from Eq. ( 46), from which one obtains Considering a radio wave travels from the Sun to the point R A , the time spent during this process can be calculated from the integral Here R 0 is the closest approach of the radio wave to the Sun which can be determined by Eq. (47).In the weak field approximation, one has The first term is the travel time of light in flat spacetime and the rest part contains both the contributions to the travel time in Schwarzschild metric and the new effects of the general Buchdahl-inspired metric.When k = 0 = η, the above expression exactly reduces to the Schwarzschild result.
There are two different cases in the experiments of measuring the gravitational time delay by sending a radar wave from Earth or spacecraft and then reflecting off another planet or spacecraft.One is the inferior conjunction case, in which the planet (or spacecraft, denoted by B), which reflects the radar signal, is located between the Earth (or spacecraft, denoted by A) and the Sun.The calculation of the time delay due to the general Buchdahlinspired metric of this case is simple and can be obtained by Another one is the superior conjunction case, in which the planet that reflects the radar signal and the Earth are on opposite sides of the Sun.For this case, the gravitational time delay is given by The most precise results related to the gravitational time delay were obtained from the Cassini experiments [30].This result ruled out many modified gravity theories that predicted larger deviations from GR.Here we would like to use its results to constrain the parameter k and η in the general Buchdahl-inspired metric.The Cassini experiment was conducted in June 2002, and the test of the gravitational time delay was achieved in the measurement of the frequency shift of radio waves to and from the Cassini spacecraft as they passed near the Sun.In the superior conjunction case, the relative change in the frequency is related to the time delay ∆t S via where ν 0 denotes the frequency of the radio waves emitted from the Earth and ν(t) is the frequency of the radio wave reflected back to the Earth at t.Using Eq. ( 59), one has Alternatively, Using the measurement performed in the Cassini experiment [30], one can constrain the parameters k and η in the general Buchdahl-inspired metric to be This constraint is stronger than that obtained by the observations of the deflection angle.For a special Buchdahl-inspired metric ( 16) with η = k and ζ = 1 + 3 k2 , the above constraint leads to a constraint on the parameter k as (64)

D. Perihelion Advance
Now we consider the orbit's perihelion advance for a massive particle moving in the general Buchdahl-inspired metric.For a massive particle, we have ε = −1.We still start with Eq. ( 44) by introducing a new variable x = 1/R, which leads to Taking the derivative on both sides of the above equation with respect to φ and expanding the equation about the small parameter r s (in the weak field approximation), one obtains The right-hand side of the above equation can be treated as perturbations to Newtonian gravity, which contains two parts.The first part 3G N M x 2 represents the correction from the Schwarzschild metric in GR and the last two terms (the second and the third term on the right-hand side of the above equation) are the new effects from the general Buchdahl-inspired metric.When the perturbations are absent, the above equation has an exact solution x 0 (φ) for a bounded orbit, which describes an elliptical orbit with the eccentricity e.
The perturbations from GR and the general Buchdahlinspired metric lead to the small derivation to the exact elliptical orbit.Thus we can write the orbit which contains the effects of the perturbations in the form of where x 0 is the elliptical orbit with the eccentricity e given by Eq. ( 67) and x 1 (φ) is the small corrections to the elliptical orbit, which satisfies Using x 0 (φ) in Eq. ( 67), one gets where The solution of x 1 is given by Only the last term contributes to the perihelion advance of a massive particle moving in the general Buchdahlinspired metric.For this reason, we can drop other terms and write the solution of x(φ) as where This expression represents the angular shift of the perihelia per orbit.For an ellipse described by Eq. ( 67), one can relate the angular momentum L of the massive particle to the semi-major axis a 0 of the ellipse as Thus one obtains It is evident that the above expression reduces to the Schwarzschild result by taking k = 0 = η.There are several observations of the perihelion advance that can be used to constrain the parameters η and k in the general Buchdahl-inspired metric.Now we consider three different observations related to the phenomenon of the perihelion advance in very different scales, i.e., the perihelion advances of the laser-ranged satellites orbiting the earth [31], of Mercury orbiting the Sun [32], and of the S2 star orbiting the supermassive black hole in the central region of our Milky Way galaxy [33].
Let us first consider the measured perihelion advance of LAGEOS satellites around the Earth.Using 13 years of tracking data of the LAGEOS satellites, the precession of the periapsis of the LAGEOS II satellite was measured to be [31] ∆φ ∆φ GR = 1 + (0.28 ± 2.14) × 10 −3 . (80) From this result, one obtains This bound corresponds to a constraint on the parameter k for the special Buchdahl-inspired metric with η = k and ζ = 1 + 3 k2 .We now turn to consider the observation of the anomalous perihelion advance for Mercury.The most accurate measurement of the perihelion advance was performed by the MESSENGER mission [32], which measures the perihelion advance for Mercury to be ∆φ = (42.9799± 0.009) ′′ /century. ( With this measurement, the bound on the η+ k 1+η arising from the general Buchdahl metric can be computed using the experimental error 0.009 ′′ /century, which yields This bound is better than that from the LAGEOS satellites by two orders of magnitude.Again, from the above constraint, one can get the constraint on the parameter k for the special Buchdahl-inspired metric with η = k and ζ = 1 + 3 k2 , i.e., −1.6 × 10 −5 < k < 1.6 × 10 −5 . (85) Here we would like to mention that the measurement of the anomalous perihelion advance for Mercury will be improved significantly in the near future from the joint European-Japanese BepiColombo project, which was launched in October 2018 [34,35].It is expected that this mission will improve the accuracy of the perihelion advance to be 10 −4 as/century.This accuracy is one order of magnitudes better than the current accuracy of MESSENGER mission [32].Thus, with the BepiColombo project, we expect to improve the constraints on the parameters η+ k 1+η arising from the general Buchdahl-inspired metric or the parameter k in the special Buchdahl-inspired spacetime to −10 −6 η+ k 1+η or k 10 −6 , which is much more restricted than that obtained from the Cassini experiment.
At last, let us consider the S2 star orbiting the central black hole of the Milky Way galaxy.Comparing the above two measurements, the observations of the S2 star provide a very different environment to test gravity in the strong gravity regime.The Schwarzschild precession of the S2 star was measured recently by the GRAVITY collaboration [33], which gives where per orbit period from the prediction of GR.For the effect of the general Buchdahl-inspired metric on the preces-sion, this observation leads to which corresponds to for the special Buchdahl-inspired metric with η = k and ζ = 1 + 3 k2 .

IV. GEODETIC PRECESSION OF SPINNING OBJECTS IN THE GENERAL BUCHDAHL-INSPIRED METRIC
In this section, we calculate the geodetic precession of spinning objects in the general Buchdahl-inspired metric.In curved spacetime, the evolution of a spinning particle follows two equations, the geodesics equation and the parallel transport equation, where u µ = dx µ /dλ is the four-velocity of the particle and s µ denotes the four-spin vector.u µ and s ν satisfy the following orthogonal condition The four-spin vector also satisfies the normalization condition Without generality, we consider the evolution of the particle in the equatorial plane, i.e., we can set θ = π/2 and dθ/dλ = 0. Let us study the test spinning particle moves in a circular orbit, i.e., Ṙ = 0, and its four-velocity u µ can be expressed as Then the angular velocity of the spinning particle is written as The stable circular orbit in the equatorial plane requires Solve these two equations, one has Along this circular orbit, the parallel transport equation ( 94) can be casted in the form of Differentiating (103) with respect to the affine parameter λ and converting λ → t using the relation dt = u t dλ one arrives at a second-order ordinary differential equation of s R , (107) This equation admits an exact solution where represents the oscillating frequency pertaining to the spin four-vector s µ .With the solution of s R , the other three components s t , s θ , and s φ can be immediately solved, giving In obtaining these solutions, we have used initial conditions, s t (0) = s θ (0) = s φ (0) = 0, which means spin vector s µ was initially directed along the radial direction.Comparing Eqs. ( 109) and (102), it is obvious to observe that the two frequencies, the oscillating frequency ω g of rotation of the spin vector and the orbital frequency Ω of a massive spinning particle along the circular orbit, are different.This difference leads to a precession of the spin vector.This is the phenomenon called geodetic precession.For one complete period of the circular orbit, the angle of the geodetic precession can be expressed as When k = 0 = η, the above result reduces to the geodetic precession of the Schwarzschild metric.
The geodetic precession can be tested and measured by using the gyroscopes in near-earth artificial satellites.One such experiment is performed by the Gravity probe B experiment, which was spaced at an altitude of 642 km and had an orbital time period of 97.65 min.According to GR, the geodetic effect induces a precession of the gyroscope spin axis by 6,606.1 milliarcseconds (mas) per year.Gravity Probe B measures this effect to be [36] ∆Θ = (6601.8± 18.3)mas/year, (114) which leads to For the special Buchdahl-inspired metric, the above result gives If one treats the Earth-Moon system as a gyroscope orbiting the Sun, its geodetic precession due to the gravitational field of the Sun has also been measured by using the Lunar laser ranging data.Recent measurement of the geodetic precession yields a relative deviation from GR as [37] ∆Θ − ∆Θ GR ∆Θ GR = −0.0019± 0.0064, (117) which gives Again, this bound corresponds to for the special Buchdahl-inspired metric with η = k and ζ = 1 + 3 k2 .

V. CONCLUSION
In this paper, we study the observational constraints that can be imposed on the asymptotically flat Buchdahlinspired solution.For this purpose, we calculate theoretically the effects of the parameter k on several solar system experiments and black hole observations.Specifically, we calculate in detail the deflection angle of light by the Sun, gravitational time delay, perihelion advance, and geodetic procession for massless and massive objects in the Buchdahl-inspired spacetime.With these theoretical predictions, we derive the constraints on the parameter k in the asymptotically flat Buchdahl-inspired spacetime by confronting our theoretical calculations with observations.Our results are summarized in Table I.In addition, we have provided different comparisons of parameters from modified gravity and the general relativity.For instance, Eq. (51) provides the comparison of deflection angle of light between two theories; Eq. (58) provides the comparison of gravitational time delay; Eqs. ( 79) and (86) provide similar comparison for perihelion and periastron precession in both theories while Eq.(117) gives a comparison of geodetic precession as predicted by two theories.
It is worth mentioning here that the measurement of the gravitational time delay by the Cassini experiment provides the most sensitive tool to constrain the parameter k in the Solar System.Another important constraint comes from observing the perihelion advance for Mercury by the MESSENGER mission.As we mentioned, the measurement of the anomalous perihelion advance for Mercury will be improved significantly in the near future from the joint European-Japanese BepiColombo project, which was launched in October 2018 [34,35].It is expected that this mission will improve the accuracy of the perihelion advance to be 10 −4 as/century, which can be used to improve the constraints on the parameters η+ k 1+η arising from the general Buchdahl-inspired metric or the parameter k in the special Buchdahl-inspired spacetime to −10 −6 η+ k 1+η or k 10 −6 , which is much more restricted than that obtained from the Cassini experiment.
In contrast to GR, pure R 2 gravity does not adhere to Birkhoff's theorem.As a higher-derivative characteristic, the Buchdahl parameter k of its vacuum solution exterior to a star is system-dependent.Therefore, our empirical tests as presented in this article were carried out under this premise.We have focused on the exterior vacuum solution, deferring the theoretical determination of k for future exploration [38].
The determination of k is, in principle, contingent on the composition-specifically, the equation of state and the distribution of matter within the host star.Typically, this inquiry involves matching the interior and exterior solutions across the star's surface.An alternative approach entails deriving a set of Tolman-Oppenheimer-Volkoff (TOV) equations governing the pressure and density of the star material.By numerically solving these equations in conjunction with the metric components, the exterior vacuum configuration of a star can be obtained based on a presumed equation of state and conditions at the star's center [39][40][41].Progresses on this front have recently been made by one of us in reducing the TOV equations for f (R) gravity to a single integrodifferential equation.This simplication has enabled our investigation into the interior-exterior matching for the Buchdahl-inspired solution, with detailed findings to be reported separately [38].
While most test results presented in Table I align with GR, the cases for S2 star and M87*, which may qualify as in a strong field regime, show large deviations for k from 0, albeit with too large error bars for definitive conclusions; see Eqs. (89) and (92).As investigated in [21], a k value in the range (−1, 0) has been associated with the potential formation of wormholes.Theoretically, such spacetime configurations could support the possibility of closed timelike curves, recently explored in [42].Consequently, future tests of the Buchdahl-inspired solution and pure R 2 gravity in strong field regimes may be warranted.

Table I .
(36)ary of estimates for bounds of the parameter η+ k 1+η arising in the general Buchdahl-inspired metric(36)and the parameter k in the special Buchdahl-inspired metric from several observations.