Parameterized post-Newtonian orbital effects in extrasolar planets

Perturbative post-Newtonian variations of the standard osculating orbital elements are obtained by using the two-body equations of motion in the parameterized post-Newtonian theoretical framework. The results obtained are applied to the Einstein and Brans–Dicke theories. As a results, the semi-major axis and eccentricity exhibit periodic variation, but no secular changes. The longitude of periastron and mean longitude at epoch experience both secular and periodic shifts. The post-Newtonian effects are calculated and discussed for six extrasolar planets.


Introduction
At present, the post-Newtonian effect has been exhibited gradually in the wake of unceasing development in the post Newtonian celestial mechanics and due to that the accurate degree of astronomical instruments is heightened unceasingly. Hence some authors devoted to the research on the subject and scopes, such as Estabrook (1969), Nordtvedt (1970), Rubincam (1977), Brumberg (1972Brumberg ( , 1985Brumberg ( , 2010, Damour and Deruelle (1985), Soffel et al. (1987), Soffel (1989), Klioner and Kopejkin (1992), Calura et al. (1997), Brumberg et al. (1995), Brumberg and Brumberg (2001), Iorio (2005aIorio ( , 2005bIorio ( , 2007aIorio ( , 2011a, Will (2008), Everitt et al. (2011), Kopeikin et al. (2011) et al. (2011). In the post-Newtonian celestial mechanics there are some best methods One of the best methods is the method of parameterized post-Newtonian Formalism (PPN method) because the theories include the various different gravitational theories with different parameters, such as Einstein, Brans-Dicke and other theories. Hence some authors devoted to the research on this scope, such as, Misner et al. (1973), Nordtvedt (1976), Sarmiento (1982), Will (1981. Moreover, the Asymptotic method Kopejkin 1989, 1990) and DSX method (Damour et al. 1991(Damour et al. , 1992 are also desirable. At present, some authors not only studied the post-Newtonian effect on the motion of celestial objects in the solar system, but also in the extrasolar planetary system. It is interesting and significant for studying the post-Newtonian effects on the extrasolar planets because in the extrasolar planetary system the separation between planets and primary star is nearer mutually and planet mass is nearly Jupiter mass. Hence the post-Newtonian effect on the orbital elements of extraplanets is larger. In the recent years some authors studied the post-Newtonian effect or the relativistic effect in the extrasolar planets (Calura and Montanari 1999;Miralda-Escudé 2002;Wittenmyer et al. 2005;Iorio 2006Iorio , 2011bIorio , 2011cAdams and Laughlin 2006a, 2006b, 2006cHeyl and Giadman 2007;Pál and Kocsis 2008;Jordán and Bakos 2008;Ragozzine and Wolf 2009). However, these authors used the method of the general relativity or the post-Newtonian approximation to study this problem. This paper used the parameterized post-Newtonian theories to study and calculate the parameterized post-Newtonian effect on the extrasolar planets with large eccentric orbit.

R, S and W components for the parameterized post-Newtonian perturbing acceleration in the two-body problem
The relative acceleration of two-body with the post-Newtonian parameters is given by Will (1981) Here Here f denotes the true anomaly.
↔ N is a unit vector in the radial direction and λ are unit vectors in the orbital plane. ↔ N is directed along the radial direction, while λ is perpendicular to N . In the equation m denotes Gm and the right side should multipled by c −2 . G is the gravitational constant and c is the speed of light. Equation (1) can be written as Here → denotes vector. We resolve the acceleration a into a radial component L (the unit vector normal to the orbital plane). On comparison with the expression (2), we get three scalar accelerative components R, S and W R = m r 2 (2γ + 2β) Substituting the following formulas of the problem of two body into the above formula (Smart 1953) We obtain μ m e sin f, where p = a(1 − e 2 ). We can write simply the above formulas as the following formulas: Substituting r = p/(1 + e cos f ) into the right hand sides of the expression (2), we obtain In it, it is Based on the post-Newtonian parameters (Will 1981), in the general relativity the post-Newtonian parameters α 1 = α 2 = α 3 = 0, β = 1, γ = 1, ζ 2 = 0. and in the Brans-Dicke gravitational theories α 1 = α 2 = α 3 = 0, ς 2 = 0, β = 1, γ = 1+ω 2+ω , ω is the dimensionless constant of the theory, ω = 5 (Estabrook 1969; Nordtvedt 1970).

The post-Newtonian perturbing equations and the perturbing variables
Substituting the perturbing accelerations R, S, W for the formulas (3), into the following Gaussian equations (Brouwer and Clemence 1961) Here u is the eccentric anomaly.ω is the longitude of periastron and ω is the argument of periastron. ε is the mean longitude at epoch.
We obtain the set of the post-Newtonian perturbing equations In the set of (6) we transform independent variable time t into independent variable anomaly f by using dt = r 2 df/na 2 √ 1 − e 2 and n 2 a 3 = m, and then, integrating the equations In it σ denotes arbitrary orbital elements from a, e, ω, i, Ω, and ε 0 . Substituting dσ dt for the set of (6) and dt df = r 2 /na × 2 √ 1 − e 2 into the above definite integral expressions and integrating, one obtain the perturbation variables δl = n(t − t 0 ) + δε. u is the eccentric anomaly Here l denotes the mean longitude of periastron The perturbative solutions (12) and (13) of Keplerian elements include over ten kinds of gravitational theories as shown in Table 5.1 in Will (1981). Hence the formals (12)-(13) are important and worth-while. In this paper we only select two kind theories of general relativity and Brans-Dicke.
In the above last integral expression, we have used already the next integral expressions Here u is the eccentric anomaly.

The secular variations of the orbital elements
It is seen from the results of the integration (8) that there exist the secular terms for δω and δε 0 Here u denotes the eccentric anomaly and u 0 is the value of u as t = 0. All other terms are the periodic terms for δa, δe, δω and δε. The coefficients of the periodic terms are the amplitudes of the periodic terms.
It is interesting for studying the secular terms. Hence we take the secular terms from the expressions (8) or integrating the definite integration (7) and taking the lower limit f 0 = 0 and upper limit f = 2π , the results of integration are that the periodic terms disappear and the secular terms appear per cycle by letting m = GM/c 2 , we get a = e = I = Ω = 0, The time variation of periastron passage, τ , can be derived from the following relation The secular rates per year are =ė =İ =Ω = 0, ω = ω/P (rad/yr), ε = ε/P (rad/yr), l = 2π/P + ε 0 /P (rad/yr), τ = τ/P (s/yr). (16) Here P denotes the orbital period, in yr. Substituting K 1 , K 2 , K 3 and K 4 for the expressions (9) into the expressions (14) and by replacing G and c 2 , then, we obtain the formulas for the secular variables and the variable rate in the general relativity × 9 1 − e 2 − 7 1 − e 2 (rad/cycle), In the Brans-Dicke theory

Numerical calculation for six exotrasolar planets
In this paper we choose six exoplanets: HD16871b, HD68988b, HD217107b, HD88133b, XO-3b and GJ-436b as an example for the former exoplanets, their P (d), M * (M ) are retrieved from Bodenheimer et al. (2003) and a(Au), e and m b (m J ) are cited from www.mpia.de/homes/ Lyra/planet_naming.html.; for latter three exoplanets, their P (d) a(Au), M * (M ) and e are retrieved from Jordán and Bakos (2008) and m b (m J ) is cited from http://www. exoplanet.eu/index.php. These data are listed in Table 5 of the Appendix. Substituting those data into formulas (17)-(20), we obtain the numerical results for the secular variation of the orbital element of six exoplanets in Table 1 and Table 2. The results of the numerical values of advance of periastron of HD88133b, XO-3b and GJ-436b in this paper as compared with that of three exoplanets in the other author's work (Jordán and Bakos 2008) are listed in Table 3. We can seen from the above Table 3 that both results are nearly approximate in the relativistic effect, but there are some different. The calculated values of this paper are some larger than that of Jordán and Bakos (2008) This difference results in that this paper calculatesω GR by using the mass of two-body (parent star and exoplanet) and Jordán and Bakos only consider the mass of the parent star and neglect the mass of exoplanet.

Comparison with the planets in solar system
Substituting the data of Mercury and Jupiter into the formulas (11) for ω E , we obtain the results for the comparison of the perihelion of Mercury and Jupiter per century with that of two exoplanets per century listed in Table 4.
We can see from Table 4 that the values of advance of the periastron of the exoplanets are largest than that of the planets in the solar system. Therefore, it is important and meaningful for studying the motion of the exoplanets.
6.3 On the possibility of observing these effects Let us discuss the possibility of observing these effects. In the solar system the advance of perihelion of Mercury is 42.91 per century. We may see from Table 3 that in the extrasolar planetary system the maximal value of advance of XO-3b is 14246 per century which correspond to 332 time (fold) value of advance of perihelion of the Mercury. At present, some authors applied the method of TTV (transit timing variation) or the method of TDV (transit duration variation). i.e., the secular precession can be detected through the long-term change in P obs or in T D (TDV) to the observation of the extrasolar planetary system (Agol et al. 2005;Rabus et al. 2009;Gibson et al. 2009;Iorio 2011b). Therefore, the possibility that the non-Newtonian advances of the periastra of the extrasolar planets considered can be observed is certainly interesting and deserves further studies.

Slouly orbiting planers
The author emphasezes that when we consider slouly orbiting planers, we could look like secular term over relatively short observational time interval, i.e., the relatively short arcs or the short term effect are available and important. Hence the author takes the time interval per year in the Table 1 and Table 2.  (2008) m b denotes exoplanet mass which is cited from http://www.mpia.de/homes/Lyra/planet_naming.html for the former three references, and the latter three references is cited from http://www.exoplanet.eu/index.php for m b 6.5 Prospect for further investigation (the new try) At present, the fifth force, Yukawa-like interaction has been investigated in our solar system (Iorio 2007b;Haranas et al. 2011;Tsang 2012). It may be predicted that extrasolar planets may well be used also for constraining putative fifth force, Yukawa-like interaction in the further investigation.

Conclusions
In this paper we worked out parameterized post-Newtonian effect on the orbits of celestial objects. The semi-major axis and eccentricity exhibit periodic variation, but no secular variation. The longitude of periastron and mean longitude at epoch exhibit secular and periodic variation. The inclination and the longitude of ascending node are unaffected. Such effects on the orbits of the extrasolar planets may be observed possibly because their effects of advance of periastron are large as in the calculation of this paper. The results of this paper based on the parameterized post-Newtonian gravitational metric by the work of C.M. Will, amplified and extended his work.
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.