Issue 
A&A
Volume 628, August 2019



Article Number  A80  
Number of page(s)  10  
Section  Celestial mechanics and astrometry  
DOI  https://doi.org/10.1051/00046361/201935705  
Published online  09 August 2019 
Detecting the general relativistic orbital precession of the exoplanet HD 80606b
Institut d’Astrophysique de Paris, UMR 7095, CNRS, Sorbonne Université, 98 bis boulevard Arago, 75014 Paris, France
email: luc.blanchet@iap.fr, hebrard@iap.fr, francois.larrouturou@iap.fr
Received:
16
April
2019
Accepted:
28
May
2019
We investigate the relativistic effects in the orbital motion of the exoplanet HD 80606b with a high eccentricity of e ≃ 0.93. We propose a method to detect these effects (notably the orbital precession) based on measuring the successive eclipse and transit times of the exoplanet. In the case of HD 80606b, we find that in ten years (after approximately 33 periods) the instants of transits and eclipses are delayed with respect to the Newtonian prediction by about three minutes due to relativistic effects. These effects can be detected by comparing at different epochs the time difference between a transit and the preceding eclipse, and should be measurable by comparing events already observed on HD 80606 in 2010 with the Spitzer satellite together with those to be observed in the future with the James Webb Space Telescope.
Key words: relativistic processes / celestial mechanics / planetary systems / methods: observational
© L. Blanchet et al. 2019
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
Whereas more than 4000 exoplanets have been discovered so far, the planet orbiting the G5 star HD 80606 remains a remarkable and unique case. Its eccentricity is extremely high: e ≃ 0.93. HD 20782b is the only planet reported to have a higher eccentricity (O’Toole et al. 2009) but its high evalue stands on one single measurement and has not been confirmed up to now. The high eccentricity of HD 80606b was well established as soon as the planet was discovered through radial velocity measurements (Naef et al. 2001) and has been largely confirmed by subsequent observations. Today the most accurate system parameters are those determined by Hébrard et al. (2011), who provide e = 0.9330 ± 0.0005.
The long orbital period of the planet (P = 111.4367 ± 0.0004 days) implied a low probability of the orbital plan to be aligned with the line of sight to the Earth. Still, HD 80606b was successively discovered to pass behind its parent star (planetary eclipses) by Laughlin et al. (2009), then in front of it (planetary transits) by Moutou et al. (2009; see also GarciaMelendo & McCullough 2009; Fossey et al. 2009; Winn et al. 2009; Hidas et al. 2010). That configuration allows the planetary radius to be measured (R_{p} = 0.981 ± 0.023 R_{Jup}); the inclination of the orbit is also known (i = 89.269 ° ±0.018°), which provides the mass M_{p} = 4.08 ± 0.14 M_{Jup} from the skyprojected mass derived with radial velocities. The obliquity of the system could also be measured from the Rossiter–McLaughlin anomaly observed during transits, which revealed that the orbit of HD 80606b is prograde but inclined (Moutou et al. 2009; Pont et al. 2009; Winn et al. 2009); Hébrard et al. (2011) reported an obliquity of λ = 42 ° ±8°. Such a peculiar orbit may be explained by the influence of the companion star HD 80607, located 1200 AU further away, through a Kozai–Lidov mechanism (Pont et al. 2009; Correia et al. 2011).
The particularly high eccentricity of the planetary orbit raises the question of the feasibility of detecting the general relativistic (GR) precession of its periastron. Indeed it should be enhanced to
as compared to 43 arcsec/century for Mercury (see also Correia et al. 2011).
Besides Mercury, the relativistic precession has been measured in the solar system for some asteroids like Icarus (Gilvarry 1953). The tests done in the solar system agree with GR to within 10^{−4} (Will 1993). Outside the solar system the effect is well known in binary pulsars like the Hulse–Taylor pulsar (Hulse & Taylor 1975) where it has been measured with high precision (see Taylor & Weisberg 1982), and in other binary pulsars such as the double pulsar (Kramer et al. 2006). To a further distance, there are now attempts to measure the relativistic precession of the star S2 close to the galactic centre (Parsa et al. 2017). The measure of the relativistic periastron advance in stellar systems was proposed by Gimenez (1985), but it is quite complicated to discriminate it from the tidal precession (see e.g. Wolf et al. 2010).
Figueira et al. (2016) reported that the relativistic precession was not detectable in the HD 80606 system in spite of the large dataset secured over several years. They measured the variation of the longitude of the periastron using the whole available datasets: , thus in a nonsignificant way. Indeed, the accuracy reached on ω mainly from radial velocities is ±540 arcsec (Hébrard et al. 2011), which is not accurate enough.
The transiting nature of HD 80606b allows a better accuracy to be reached thanks to the precise timing that can be measured on the midpoint of each transit or eclipse. Typical accuracies of ±85 s (Hébrard et al. 2011) and ±260 s (Laughlin et al. 2009) were obtained on the midpoint of transits and eclipses, respectively, using the Spitzer Space Telescope and its IRAC instruments. Indeed, only a satellite on an Earthtrailing heliocentric orbit such as Spitzer could allow the whole 12hourlong transit of HD 80606b to be continuously observed together with sufficiently long offtransit references immediately before and after the event.
In this paper we propose to detect the GR precession of the periastron (and the relativistic orbital motion) for the exoplanet HD 80606b by using a method based on the transit times, computing a small drift in the successive instants of transits due to relativity. We find that in ten years (after 33 periods), due to the relativistic effects, the instants of transits are delayed with respect to the Newtonian predicted ones by about three minutes. Whereas the transit midtime could be measured with such an accuracy, the absolute value of that shift could not be detected because the system’s orbital parameters are not known today with a high enough accuracy.
However, we argue that the effect could actually be measurable from the observation of the change of the elapsed time t_{tr − ec} = t_{tr} − t_{ec} between an eclipse and the next transit^{1}. The time between an eclipse and a transit was measured in January 2010 using Spitzer to be
with a good accuracy of ±275 s. Due to the GR precession of the periastron, and the modification of the period due to relativity, we find that this time difference will be reduced by 183 s in 2020 (after 33 orbits), and by 271 s in 2024 (49 orbits). Such shift would be difficult to detect in 2020 using IRAC on Spitzer, but should be detectable a few years later using NIRCam or MIRI on James Webb Space Telescope (JWST).
Most of the paper is devoted to theoretical aspects, and it ends with a conclusion and discussion in Sect. 6 where we give further comments on previous observations and on the Newtonian perturbing effects. In Sect. 2 we introduce the geometrical conventions used to compute the times of transits and eclipses. Those are computed in two different ways: a postKeplerian parametrisation of the orbit in Sect. 3, and a Hamiltonian method using DelaunayPoincaré canonical variables in Sect. 4. Those methods give the exact solution of the first postNewtonian equations of motion. Our final results for the transit times are given in Tables 1–3 below. Furthermore, as a check of the latter methods, we also computed the transit times with a method relying on celestialmechanics perturbation theory in Sect. 5. As we used it, this procedure is valid on average over one orbit, and takes only into account the secular effects. We find that it is in good agreement with the postKeplerian and the Hamiltonian methods. A quick but rough estimate of the GR shift of transit times is also given in Appendix A.
Differences (in seconds) between the instants of transit computed with the relativistic (1PN) model and with the Newtonian model, taking as reference point (δt ≡ 0) the passage at periastron just before the transit of Jan. 2010.
Differences (in seconds) between the instants of eclipse computed with the relativistic (1PN) model and with the Newtonian model, taking as reference point () the passage at periastron just after the eclipse of Jan. 2010.
Differences (in seconds) between the time interval t_{tr − ec}(N) between the passage at the minimum approach point during the Nth eclipse and the passage at the minimum approach point during the Nth transit, and the time interval t_{tr − ec}(0) at the reference period, when the eclipse and transit where measured in January 2010 with the Spitzer satellite (Hébrard et al. 2011).
2. Geometry of the planetary transits
We introduce a frame {x, y, z} with origin at the centre of the star. The observer is in the direction x (Earth’s direction), while the projected angular momentum of the star is by convention along the direction z. Hence the plane {y, z} is the plane of the sky. The orbital plane of the planet is defined by its inclination I with respect to the plane {x, y}, and by its longitude Ω with respect to the plane {x, z}, following our conventions described in Fig. 1a. As depicted in Fig. 1c, the motion of the planet in the orbital plane is parametrised by polar coordinates (r, φ), where φ denotes the sum of the true anomaly and of ω_{0}, the latter denoting a constant angle, which can be viewed as defining the initial argument of the periastron at a reference time. Then the coordinates (x, y, z) of the planet in this frame read (see Figs. 1a and c)
Fig. 1. Geometry of the planetary transit: (a) Motion of the planet with respect to the observer, located at infinity in the direction of x of the frame {x, y, z}. The centre of the star is denoted O, P is the periastron at reference time t_{P}, and J_{⋆} is the angular momentum of the star projected on the plane of the sky {y, z}. (b) The transit as seen by the observer and the definition of the points {T_{i}}; here are only shown T_{1} and T_{3}. (c) The trajectory of the planet and the points {T_{i}} (transit) and (eclipse) as seen in the orbital plane {X, Y}. Up to an angle Ω, the observer is in the direction of X. (d) View of the transit from above. The direction (ON) is the line of “nodes”, corresponding to (OX) in (a). 

Open with DEXTER 
Since x is the direction of the observer, the condition for the (partial or complete) transit of the planet is that y^{2} + z^{2} ≤ (R_{⋆} + R_{p})^{2}, where R_{⋆} and R_{p} are the radius of the star and planet, hence
Since R_{⋆} + R_{p} ≪ r, it is easy to see that this condition can be satisfied when sin I ≪ 1 or sin φ ≪ 1, and cos I cos Ωsinφ + sinΩcosφ≪1. With the conventions of Fig. 1a we are interested in a transit for which φ ≃ π and Ω ≃ 0^{2}. On the other hand, for an eclipse during which the planet is passing beyond the star as seen from the observer, we have φ ≃ 0 and Ω ≃ 0. In an expansion to the first order in we find that the coordinates of the planet during and around the transit are
where the impact parameter b of the planet is the distance between the trajectory of the planet projected on the plane {y, z} and the centre of the star. Neglecting the bending of the trajectory during the transit and eclipses^{3}, the impact parameter is given by (see Fig. 1b)
where r(π) and r(0) are the values of the radial coordinate of the planet at φ = π and φ = 0, when the transit and eclipse occur. Since the eclipse will be seen from behind the star and so much farther away, we shall have to add to the instants of eclipse a Roemer time delay accounting for the propagation of light at finite velocity c, and which is a 0.5PN effect ∝1/c. Thus where
We define the positions of the planet {T_{1}, T_{2}, T_{3}, T_{4}} to be respectively the entry of the transit, the start of the full transit, and the exits of the full and partial transits. Similarly the points corresponding to the eclipse are shown in Fig. 1c. The transit and eclipse conditions for the polar coordinates (r, φ) at the passage of these points read
Given a model r(φ) for the planetary orbit, this determines the true anomaly φ for each of these points, where the vertical coordinate Y(b) is given in terms of the impact parameter b by (see Figs. 1b and c)
The cases for which (R_{⋆} ± R_{p})^{2} = b^{2} correspond to the planet just grazing tangentially the star, either from the exterior or from the interior of the stellar disc. In addition to Eqs. (8) and (9) we compute also the points corresponding to the minimal approach to the center of the star, T_{m} and , whose transit conditions are obviously (see Fig. 1b)
Having determined the polar coordinates of the planet we shall deduce the instants t_{i} and of passage at each of these points by using the law of motion of the planet on the trajectory.
3. PostNewtonian motion of the planet
In this section we compute the times of transits and eclipses, that is, the instants of the successive passages at the positions T_{i} and defined by the geometric conditions in Eqs. (8)–(10), in a relativistic model, including the first postNewtonian (1PN) corrections beyond a Keplerian orbit. More precisely, we compute the difference between the transit times predicted by the relativistic 1PN model and those predicted by the Newtonian model. This will permit us to assess whether the latter relativistic corrections could be detectable, and we propose an observable quantity, directly measurable in future observations of the planet HD 80606b.
3.1. Keplerian parametrisation of the orbit
In the Keplerian model the orbit of the planet and its motion on the orbit are given by the usual six orbital elements {a, e, I, ℓ, ω_{0}, Ω}, where I and Ω specify the orbital plane (see Fig. 1a), a is the semimajor axis, e is the eccentricity, ω_{0} denotes the constant angular position of the periastron, and the mean anomaly ℓ is defined by ℓ = n_{0}(t − t_{0, P}), where t_{0, P} is the instant of passage at the periastron and n_{0} = 2π/P_{0} is the mean motion, with P_{0} the period. The Keplerian relative orbit in polar coordinates (r, φ), is conveniently parametrised by the eccentric anomaly ψ as
The relative orbit in polar coordinates reads
from which we deduce the impact parameter in the Newtonian model, say b_{0}, to be
In the Newtonian model we have Kepler’s third law n_{0} = (GM/a^{3})^{1/2}, where M = M_{⋆} + M_{p} is the total mass of the star plus planet system. Furthermore the energy and the angular momentum of the orbit are given by the usual Keplerian formulas
where μ = M_{⋆}M_{p}/M is the reduced mass.
The conditions for the transits and eclipses are given by Eqs. (8) with (9), and determine the values of the eccentric anomaly at these points. In the Newtonian analysis the eccentric anomaly ψ_{0} is obtained by solving
where Y_{0} = Y(b_{0}) denotes the vertical coordinate (Eq. (9)) computed with the Newtonian impact parameter given by Eq. (13).
3.2. Relativistic corrections to the Keplerian motion
Next we consider the relativistic 1PN model see (e.g. Wagoner & Will 1976). For this we shall use an explicit solution of the equations of motion that is particularly elegant, called the quasiKeplerian representation of the 1PN motion (Damour & Deruelle 1985). In this representation the orbit and motion are parametrised by the eccentric anomaly ψ in a Keplerianlooking form as
where ℓ = n(t − t_{p}) and n = 2π/P, with now P being the period of the relativistic model. The motion is parametrised by six constants: the mean motion n, a particular definition of semimajor axis a_{r} differing from a in the Newtonian model by a small 1PN term, the precession of the orbit K including the relativistic precession, and three types of eccentricities e_{r}, e_{t}, and e_{φ} differing from e and from each other by small 1PN terms. We shall pose
The constant k denotes the relativistic advance of the periastron per orbital revolution, which is where Δ is the angle of return to the periastron, and whose GR value has been given in Eq. (1). The effect of precession can be seen clearly with the expression of the relative orbit in polar coordinates, which takes the more involved form (Damour & Deruelle 1985)
Besides the usual contribution of the precession , we note the presence of a term vanishing in the small mass ratio limit, when ν = μ/M → 0.
Crucial to the postKeplerian formalism are the explicit expressions of all these constants in terms of the conserved energy E and angular momentum J of the orbit. Since we want to compare the relativistic model of Eq. (16) to the Newtonian one of Eq. (11) we will compare, for convenience, orbits that have the same energy E and angular momentum J. Thus we need to express the relations in terms of the “Newtonian” semimajor axis a and eccentricity e defined in terms of E and J by the Newtonian formulas in Eq. (14). We have (Damour & Deruelle 1985; see also Blanchet 2014)
Among the effects described by these constants, the main ones are those associated with the mean motion n and the periastron precession K. The expressions for n and K in terms of E and J are invariant, meaning they do not depend on the coordinate system. Since
we expect the effect of the precession of the orbit to be dominant over the shift in the orbital period. The other relations, linking the semimajor axis and eccentricities to E and J, depend on the coordinate system; they are given here in harmonic coordinates. In this paper we compute an observable quantity (defined by Eq. (30) below), which thus does not depend on the choice of coordinate system.
We now investigate the difference between the analysis using the relativistic model and the Newtonian one concerning the planetary transits. The conditions for the transits and eclipses are still given by Eqs. (8) with (9), with the impact parameter given by Eq. (6). However, the impact parameter will differ from the Newtonian value by an amount δb. We readily find this amount thanks to the explicit formula Eq. (18), using a first order expansion in the small parameters given by Eq. (19),
where φ_{N} = (2N + 1)π − ω_{0} and . Here N is the number of orbits since the reference time, with N = 0 corresponding to the transit observed in January 2010. Then the coordinate Y(b) is modified by and we find (see Eq. (9))^{4},
Naturally for T_{m} and , Eq. (10), we just have δY = δbcotI.
We are now in a position to solve the transit conditions and find the modification of the eccentric anomaly δψ with respect to the Newtonian model, that is, the solution of Eq. (15). Among all the relativistic effects contributing to δψ, only the one associated with the relativistic precession k is “secular”, and therefore grows when ψ_{0} → ψ_{0} + 2πN proportionally to the number of orbits N. The other effects, associated with the relativistic parameters ξ, ε_{r}, and ε_{φ}, are periodic and therefore average to zero in the long term. Accordingly we split
Similarly, we distinguish δY = δY_{secular} + δY_{periodic}, where the secular part is given by the Ndependent terms in Eq. (21). Perturbing Eq. (8) around the solution ψ_{0} of the Newtonian condition of Eq. (15), we explicitly find
Finally, having computed δψ with the transit condition we readily obtain the time difference between the instants of transits (and eclipses) in the relativistic model and in the Newtonian model by using the equation for the mean anomaly ℓ = n(t − t_{P})=ψ − e_{t} sin ψ. We notice that there is another secular term associated with the correction in the period, n = n_{0}(1 + ζ), where ζ is the second invariant of the motion besides the relativistic precession k. We thus write similarly δt = δt_{secular} + δt_{periodic} and find
When computing , the 1.5PN shift in the Roemer effect also has to be taken in account, although we find that this effect is very small. With the help of Eq. (7), it comes after N orbits
Finally let us comment on the approximation we have made here, namely the expansion to first order in the relativistic parameters k, ζ, ξ, and so on. This approximation is valid as long as the number of orbits N over which we integrate is much smaller than the inverse of these relativistic effects, for example as long as N ≪ 1/k. In the case of HD 80606b we have k ≃ 5 × 10^{−7} and we are considering N = 33 orbits between the transit of January 2010 and that of 2020 (see the Introduction), so our approximation is justified.
We present our results^{5} in Tables 1–3. From the above analysis we compute for each orbit N the quantities (with i ∈ {1, 2, m, 3, 4})
where t_{0, i}(N) and denote the Newtonian values. Obviously we have t_{0, i}(N)=t_{0, i}(0)+NP_{0} and since the Newtonian orbit is fixed, and the motion periodic with period P_{0} = 2π/n_{0}. Together with the shift of instants of transit, we compute also the relativistic effect on the total duration and the entrance of the Nth transit, as defined by
together with similar quantities and for the Nth eclipse. The results are given in Tables 1 and 2.
As we see in Table 1, the relativistic model predicts that after ten years the planet will begin the transit about three minutes earlier than the Newtonian prediction, and about 4.5 min earlier after 15 years. Given that the precision on the measured instants of transits can be as good as 85 s (Hébrard et al. 2011), this is already an indication that the relativistic effects might be measurable if the orbital parameters were particularly well known. However, the tenyearslong relativistic effect on the transit duration δt_{14} and transit entrance δt_{12} are on the order of seconds, and seem a priori hardly detectable. Concerning the eclipse in Table 2, the effect is much smaller, of the order of ten seconds after ten years with respect to the Newtonian model. As the effect on the transit and the eclipse timings is different, this provides a way to directly detect it.
Next we define t_{tr − ec}(N) to be the time interval between the passage at the minimum approach point during the Nth eclipse and the passage at the minimum approach point during the Nth transit:
By a straightforward calculation, using the fact that the Newtonian orbit is fixed, we obtain
The important point about this result is that the righthand side of the equation represents the relativistic prediction, which is tabulated in Tables 1 and 2, while the lefthand side is directly measurable: namely t_{tr − ec}(0) has been measured with good precision between the eclipse and the transit of 2010, while future observation campaigns could allow t_{tr − ec}(N) to be measured ten to fifteen years later. From Table 3 we predict that the observations in 2020 will measure that the time difference between the transit and eclipse (with N = 33) is shorter by 182.8 s (approximately three minutes), because of the relativistic effect, as compared to the one which was measured ten years ago in 2010. The values of this observable quantity given by Eq. (30) are also presented for the transits of the next seven years in Table 3; in December 2024, the difference between the eclipse and transit (with N = 49) is shorter by 271.4 s (4.5 min) as compared to the one measured in 2010.
The 2010 measurement yielded (Laughlin et al. 2009; Hébrard et al. 2011)
with good published uncertainty in the measurement, of the order of 4.5 min, comparable to the relativistic effect we want to detect. Therefore we conclude that for the next observations of the transit/eclipse in the coming years, we should already be on the verge of detecting the relativistic effect.
4. Hamiltonian integration of the motion
At the 1PN order, the Hamiltonian of the relative motion of two point masses reads (Blanchet 2014)
where is the reduced linear momentum, which is conjugate to X, and we pose P_{R} = P ⋅ X/R with R = X, and P^{2} = P^{2}. We shall use the DelaunayPoincaré canonical variables. We start with the usual orbital elements {a, e, ℓ, ω}, where a and e are the semimajor axis and eccentricity, ω is the argument of the periastron, and ℓ = n(t − t_{P}) denotes the mean anomaly, where by definition^{6}. Since the motion takes place in a fixed orbital plane, we do not need to consider the inclination I and the longitude of the node Ω. Thus we have
together with ℓ = ψ − e sin ψ, where ψ is the eccentric anomaly. Then Eq. (32) reads
where we introduced 𝒳 = 1 − e cos ψ for convenience. The conjugate DelaunayPoincaré variables {λ, Λ, h, ℋ} are then defined by
and the Hamiltonian equations take the canonical form
With the explicit expression of the 1PN Hamiltonian in Eq. (34) we obtain
Let us now turn to the integration of the latter equations. As we see from Eq. (37) we need to obtain the indefinite integrals and . Looking at Gradshteyn & Ryzhik (1980), one finds all the needed integrals:
With these elementary integrals it is straightforward to obtain the complete solution of the motion. We first compute Λ(ψ) and ℋ(ψ). To the zeroth order their values Λ_{0} and ℋ_{0} are constant, so the solution is described by constant orbital elements a_{0}, e_{0}, and , and we can integrate using n_{0}dt = 𝒳_{0}dψ to the requested order, where 𝒳_{0} = 1 − e_{0} cos ψ. Using the definitions of Λ and ℋ in Eq. (35) we obtain the solution as a = a_{0} + δa and e = e_{0} + δe, where
In exactly the same way, dh/dt can be integrated to give ω = ω_{0} + δω, with
Finally, we integrate the mean anomaly ℓ = ψ − e sin ψ = λ + h. The point is that in the equation for λ, we first need to express the mean motion as n = n_{0} + δn, where and δa has been found in Eq. (39). We thus obtain ℓ = ℓ_{0} + δℓ, where ℓ_{0} = n_{0}(t − t_{0, P})=ψ − e_{0} sin ψ (with t_{0, P} the constant passage at periastron), and
We note that we could separate, as in Sect. 3, the secular contributions given by the first terms in both δω and δℓ from the manifestly periodic ones.
Combining all those equations, we can solve the transit condition to Newtonian order (Eq. (15)), which gives ψ_{0}, and then the correction δψ coming from Eq. (22), where the modification of the impact parameter is given by perturbing Eq. (13):
Finally this gives the modification of the instants of transits due to the 1PN corrections as
where now 𝒳_{0} = 1 − e_{0} cos ψ_{0}. This method gives the same results as the postKeplerian ones, reported in Tables 1 and 2, up to 0.5% after 33 cycles. This difference can be explained by comparing Eqs. (21) and (42): the computation of δb slightly differs depending on the method used. Moreover by replacing roughly φ_{N} → φ_{N}/(1 + k) in Eq. (21), the second order perturbations can be estimated to play a role at the level of ∼5 × 10^{−6}.
5. Lagrangian perturbation theory
It is instructive to work out the problem by means of a Lagrangian perturbation method. In this case we shall focus on secular effects (on a timescale much longer than the orbital period) after averaging the perturbation equations. The Lagrangian of the relative motion of two point masses at the 1PN order is given by
where the 1PN perturbation function ℛ depends on the relative position x and velocity v = dx/dt of the particles (we pose r = x and ṙ = dr/dt), and reads (see e.g. Blanchet & Iyer 2003),
In the Lagrangian perturbation formalism, when looking at secular effects, we can apply the usual perturbation equations of celestial mechanics directly with the perturbation function ℛ in the Lagrangian (see e.g. Damour & EspositoFarèse 1994). As in Sect. 4 we choose the independent orbital elements to be {a, e, ℓ, ω}, and we pose n = (GM/a^{3})^{1/2}. The perturbation (Eq. (45)) depends only on a, e, and ℓ, with the dependence on ℓ being implicit in the eccentric anomaly ψ(e, ℓ), which is the solution of the Kepler equation, Eq. (11b). The four left perturbation equations read (see e.g. Brouwer & Clemence 1961; Blanchet & Novak 2011)
It must be remembered that in perturbation theory n is not equal to 2π/P with P the period. Indeed we have the usual definition ℓ = n(t − t_{P}) but the instant of passage at periastron depends on time: t_{P} = t_{P}(t). Instead the period will be given by averaging Eq. (46c) for ℓ.
As in Sect. 4 we denote the constant orbital elements to zeroth order by a_{0}, e_{0}, and pose . To first order the equations for a and e can readily be integrated as
where the perturbation function reduces in this case to (with 𝒳_{0} = 1 − e_{0} cos ψ)
The results of Eqs. (47) and (48) coincide exactly with those of the Hamiltonian formalism, Eq. (39). Next, the equations for ℓ and ω to first order become
Finally we have to perform the orbital average of these equations, . Either we can compute the average directly from Eq. (49), or in a simpler way we can substitute ℛ by its average ⟨ℛ⟩ in Eq. (49). We have
hence we end up with the two main results
These are in agreement with the relative modification of the mean motion or period P given by the relativistic parameter ζ in Eq. (19a), and of course with the relativistic precession parameter k in Eq. (19b). Thus the two averaged Eqs. (51) reproduce the two gaugeinvariant equations of the quasiKeplerian formalism in Sect. 3. They can also be obtained by averaging the equations of the Hamiltonian formalism in Sect. 4.
The times of transit and eclipse are straightforwardly computed in perturbation theory using Eq. (51) and the transit conditions of Eq. (8). In this way we have confirmed the results reported in Tables 1–3, but there is a typical difference of nearly 2% after 33 cycles, due to the fact that in this perturbative approach we neglect the periodic terms. As we discussed in Sect. 3, the relativistic effect on the times of transit after a large number of orbits essentially depends on the precession k, and the modification of the orbital period ζ, with the latter effect being smaller (see Eq. (20)).
6. Conclusion and discussion
In this paper we investigated the relativistic effects in the orbital motion of the high eccentricity exoplanet HD 80606b orbiting around the G5 star HD 80606. We proposed a method to detect these effects, based on the accurate measurement of the elapsed time between a midtransit instant of the planet passing in front of the parent star, and the preceding mideclipse instant when it passes behind the star.
We presented different computations of the relativistic effects on the transit and eclipse times. One is based on the postKeplerian parametrisation of the orbit in Sect. 3, another one is a Hamiltonian method with DelaunayPoincaré canonical variables in Sect. 4, and finally a Lagrangian perturbation approach for computing the secular effects in Sect. 5. A crude understanding and estimation of the main effect is also provided in Appendix A. These methods gave consistent results, which are reported in Tables 1–3.
We found that in ten to fifteen years, corresponding to 33–49 orbital periods, the time difference between the eclipse and the next transit is reduced by approximately 3–4.5 min due to the relativistic effects for this planet. We conclude that these effects should be detectable for the next observations of the full transit and eclipse of HD 80606b in coming years by comparing to previous observations done in 2010.
By comparison to shortperiod exoplanets, only a few longduration eclipses and transits of HD 80606b have been observed today. Laughlin et al. (2009) reported the whole eclipse measured with Spitzer in November 2007 with a midtime accuracy of ±260 s. de Wit et al. (2016) secured a similar observation with Spitzer of the whole eclipse of January 2010. A few days after again with Spitzer, Hébrard et al. (2011) observed the whole transit of January 2010 with a midtime accuracy of ±85 s. Roberts et al. (2013) observed the same transit with the MOST satellite and reached a lower accuracy of ±294 s. All the other available observations were made from the ground in February 2009 (Moutou et al. 2009; GarciaMelendo & McCullough 2009; Fossey et al. 2009; Hidas et al. 2010), June 2009 (Winn et al. 2009), or January 2010 (Shporer et al. 2010) and could only observe partial portions of transits with a given telescope; so their accuracies on the midtransit times were ±310 s or even poorer.
Thus, the best available measurement of t_{tr − ec} was obtained in January 2010 with Spitzer with an accuracy of ±275 s. A similar accuracy could be reached in 2020 with Spitzer, but it would be probably insufficient to detect the effect, slightly below 200 s at that time. Thanks to its 6.5m aperture diameter, the James Webb Space Telescope (JWST) would provide an accuracy on t_{tr − ec} a few times better than Spitzer (0.85m diameter). So after JWST has started its operation, hopefully in 2021, it should be feasible to detect the effect with its NIRCam or MIRI instruments, while the effect will approach 300 s by comparison to 2010. Spitzer and JWST are among the rare telescopes able to continuously observe the 12hour duration transit of HD 80606b, which is mandatory to reach a good timing accuracy. Other current or future telescopes such as TESS, CHEOPS, or PLATO could also sample longduration transits; however, they are smaller and less sensitive than Spitzer or JWST, so are not able to reach their accuracy on timing. In addition they operate in the optical domain, whereas Spitzer and JWST observe in infrared where eclipses are deeper and thus more accurately measured. JWST clearly is the best facility to attempt that detection.
Finally we quickly discuss the other perturbing effects that may affect the measurement. HD 80606b is the only planetary companion detected today in that system. More than 15 years of radialvelocity monitoring^{7} did not provide a hint of any additional companion (Hébrard et al. 2011; Figueira et al. 2016). With the available data, one can put a 3σ upper limit of 0.4 m/s/yr of a possible linear drift in addition to the radial velocity signature of HD 80606b; this allows any additional planetary companion with a skyprojected mass larger than 0.5 Jupiter mass to be excluded with an orbital period shorter than 40 years around HD 80606. On the other hand, the companion star HD 80607 is too far away (1200 AU) to produce a noticeable effect on the orbital precession of HD 80606b.
Besides the possibility that there might be some disturbing other bodies, we have to consider the precession of the orbit caused by the oblateness of the star, the tidal effects between the star and the planet, and the Lense–Thirring effect. The orbital precession rate due to the quadrupolar deformation J_{2} of the star is
Assuming the solar value J_{2} ∼ 10^{−7} we obtain Δ_{J2} ∼ 0.4 arcsec/century for HD 80606b, which is negligible for our purpose. As for the orbital precession induced by the tidal interaction of the planet with its parent star, it produces a supplementary orbital precession at the rate (see e.g. Jordàn & Bakos 2008; Fabrycky 2010)
where k_{p} denotes the Love number of the planet, which is approximately k_{p} ∼ 0.25 for a hot Jupiter, and k_{⋆} that of the star, expected to be approximately k_{⋆} ∼ 0.01 (Jordàn & Bakos 2008). We find that the first term, due to the tidal deformation of the planet, gives a rather large contribution of about 32 arcsec/century, while the second term in Eq. (53), due to the deformation of the star, is much smaller, approximately 2 arcsec/century. Hence the tidal interaction represents a nonnegligible effect, but nevertheless one that is smaller than the GR relativistic precession:
Hence we conclude that the tidal interaction should be taken into account in the precise data analysis of the transit times aimed at measuring the GR effect. The Lense–Thirring effect, due to the angular momentum S_{⋆} of the star, induces a precession of the line of nodes with a rate of
when M_{p} ≪ M_{⋆}. Considering the extreme case of a spinorbit alignment and assuming a solar value S_{⋆} ∼ 10^{42} kg m^{2} s^{−1}, it is negligible, as expected: Δ_{LT} ∼ 0.07 arcsec/century.
To compare with the solar system, those three effects are negligible for Mercury, for which they are on the order of Δ_{J2} ≃ 0.2 arcsec/century, Δ_{T} ≃ 2 × 10^{−6} arcsec/century, and Δ_{LT} ∼ 2 × 10^{−3} arcsec/century. In the case of Mercury the periastron advance is largely dominated by the Newtonian effect of the other planets, mainly Venus and Jupiter: Δ_{planets} = 532.3 arcsec/century ≃12.4 Δ_{GR}. From this point of view, HD 80606/HD 80606b represents a cleaner system than the solar system.
While in the final stage of the preparation of this paper, we found that this method was already suggested by Jordàn & Bakos (2008), Pàl & Kocsis (2008) and later by Iorio (2011). Those works assume a central transit and use only the effect of the periastron shift. Furthermore they are generic while here we also investigate in detail the remarkable case of HD 80606b.
For HD 80606b, Hébrard et al. (2011) gave sinΩ ≃ 0.019.
When (R_{⋆} ± R_{p})^{2} = b^{2} the 1PN perturbation of the Newtonian transit condition breaks down since the denominator of Eq. (22) is zero. This corresponds to the case where the planet is just grazing the disc of the star either from inside or outside. For HD 8606b we have b_{0} ≃ 0.73 (R_{⋆} + R_{p}) and b_{0} ≃ 0.90 (R_{⋆} − R_{p}) so we are far from this case.
We took our values in Table 1 from Hébrard et al. (2011). We note that our angular conventions in Fig. 1a are different: I = λ (hence I = 42^{ °} adopted here) and .
Heyl & Gladman (2007) proposed a method to detect nontransiting planetary companions via the Newtonian secular precession they induce on an (observed) transiting planet, over and above the effect of general relativity.
Acknowledgments
The authors would like to thank Alain Lecavelier des Étangs and Gilles EspositoFarèse for motivating discussions at an early stage of the work.
References
 Blanchet, L. 2014, Liv. Rev. Relativ., 17, 2 [NASA ADS] [CrossRef] [Google Scholar]
 Blanchet, L., & Iyer, B. R. 2003, Class. Quant. Grav., 20, 755 [NASA ADS] [CrossRef] [Google Scholar]
 Blanchet, L., & Novak, J. 2011, MNRAS, 412, 2530 [NASA ADS] [CrossRef] [Google Scholar]
 Brouwer, D., & Clemence, G. 1961, Methods of Celestial Mechanics (New York: Academic Press) [Google Scholar]
 Correia, A. C. M., Laskar, J., Farago, F., & Boué, G. 2011, Celest. Mech. Dyn. Astron., 111, 105 [NASA ADS] [CrossRef] [Google Scholar]
 Damour, T., & Deruelle, N. 1985, Ann. Inst. H. Poincaré Phys. Théor., 43, 107 [Google Scholar]
 Damour, T., & EspositoFarèse, G. 1994, Phys. Rev. D, 49, 1693 [NASA ADS] [CrossRef] [Google Scholar]
 de Wit, J., Lewis, N. K., Langton, J., et al. 2016, ApJ, 820, L33 [NASA ADS] [CrossRef] [Google Scholar]
 Fabrycky, D. 2010, in NonKeplerian Dynamics of Exoplanets, ed. S. Seager (Tucson: The University of Arizona Press), 217 [Google Scholar]
 Figueira, P., Santerne, A., Suárez Mascareño, A., et al. 2016, A&A, 592, A143 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Fossey, S. J., Waldmann, I. P., & Kipping, D. M. 2009, MNRAS, 396, L16 [NASA ADS] [CrossRef] [Google Scholar]
 GarciaMelendo, E., & McCullough, P. R. 2009, ApJ, 698, 558 [NASA ADS] [CrossRef] [Google Scholar]
 Gilvarry, J. J. 1953, Phys. Rev., 89, 1046 [NASA ADS] [CrossRef] [Google Scholar]
 Gimenez, A. 1985, ApJ, 405, 167 [Google Scholar]
 Gradshteyn, I., & Ryzhik, I. 1980, Table of Integrals, Series and Products (New York: Academic Press) [Google Scholar]
 Hébrard, G., Desert, J.M., Diaz, R. F., et al. 2011, A&A, 516, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Heyl, J. S., & Gladman, B. J. 2007, MNRAS, 377, 1511 [NASA ADS] [CrossRef] [Google Scholar]
 Hidas, M. G., Tsapras, Y., Street, R. A., et al. 2010, MNRAS, 406, 1146 [NASA ADS] [Google Scholar]
 Hulse, R., & Taylor, J. 1975, ApJ, 195, L51 [NASA ADS] [CrossRef] [Google Scholar]
 Iorio, L. 2011, MNRAS, 411, 167 [NASA ADS] [CrossRef] [Google Scholar]
 Jordàn, A., & Bakos, G. 2008, ApJ, 685, 543 [NASA ADS] [CrossRef] [Google Scholar]
 Kramer, M., Stairs, I. H., Manchester, R. N., et al. 2006, Science, 314, 97 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Laughlin, G., Deming, D., Langton, J., et al. 2009, Nature, 457, 562 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Moutou, C., Hébrard, G., Bouchy, F., et al. 2009, A&A, 498, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Naef, D., Latham, D. W., Mayor, M., et al. 2001, A&A, 375, L27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 O’Toole, S. J., Tinney, C. G., Jones, H. R. A., et al. 2009, MNRAS, 392, 641 [NASA ADS] [CrossRef] [Google Scholar]
 Pàl, A., & Kocsis, B. 2008, MNRAS, 389, 191 [NASA ADS] [CrossRef] [Google Scholar]
 Parsa, M., Eckart, A., Shahzamanian, B., et al. 2017, ApJ, 845, 22 [NASA ADS] [CrossRef] [Google Scholar]
 Pont, F., Hebrard, G., Irwin, J. M., et al. 2009, A&A, 502, 695 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Roberts, J. E., Barnes, J. W., Rowe, J. F., & Fortney, J. F. 2013, ApJ, 762, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Shporer, A., Winn, J. N., Dreizler, S., et al. 2010, ApJ, 722, 880 [NASA ADS] [CrossRef] [Google Scholar]
 Taylor, J., & Weisberg, J. 1982, ApJ, 253, 908 [NASA ADS] [CrossRef] [Google Scholar]
 Wagoner, R., & Will, C. 1976, ApJ, 210, 764 [NASA ADS] [CrossRef] [Google Scholar]
 Will, C. M. 1993, Theory and Experiment in Gravitational Physics (Cambridge: Cambridge University Press) [CrossRef] [Google Scholar]
 Winn, J. N., Howard, A. W., Johnson, J. A., et al. 2009, ApJ, 703, 2091 [NASA ADS] [CrossRef] [Google Scholar]
 Wolf, M., Claret, A., Kotková, L., et al. 2010, A&A, 509, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A: Rough estimation of the effect
It is also possible to estimate the PN effects by considering a stepbystep shift of the orbit along an approximate sequence of Keplerian ellipses. As depicted in Fig. A.1, at each step we have to consider two time delays: the one induced by the new position of the transit on the orbit, δt_{k}, and the time to pass from the previous periastron to the current one, δt_{P}.
Fig. 2. Displacement of the trajectory and time shifts of the exoplanet inbetween two approximate successive Keplerian orbits N − 1 and N. The two main contributions δt_{k} and δt_{P} of the (rough) estimation of the time shift are indicated. The point is the point corresponding to T_{N − 1} reported in the next orbit N. 

Open with DEXTER 
The first delay is due to the different orientation of the Nth Keplerian ellipse with respect to the (N − 1)th one, shifted by the relativistic precession angle ω_{N} − ω_{N − 1} ≃ 2πk (see Fig. A.1). It simply reads n_{0} δt_{k} ≃ ℓ_{N} − ℓ_{N − 1}, where n_{0} is the Newtonian mean motion and ℓ_{N} the mean anomaly solving the transit condition of the Nth orbit (for a given transit T_{i} or eclipse ). Of course the mean anomaly is counted from the periastron P_{N} of the Nth orbit, and is given in terms of the eccentric anomaly by Kepler’s equation ℓ_{N} = ψ_{N} − e sin ψ_{N}. The condition relating the corresponding true anomalies of the successive ellipses is thus
As for the second time delay, it is due to the modification of the mean motion (and thus the period) as given by n = n_{0}(1 + ζ) (see Eq. (19a)). We have ℓ = n_{0}(t − t_{PN}), which we can roughly equate to n(t − t_{PN − 1}), where t_{PN − 1} and t_{PN} are the successive instants of passage at periastron. Thus t_{PN} − t_{PN − 1} ≃ −ζ(t − t_{PN − 1}), and we roughly approximate t − t_{PN − 1} by the period so that n_{0} δt_{P} ≃ −2πζ. The time delay between the Nth orbit and the reference one finally becomes
It is clear from Fig. A.1 that the two effects are negative, and add up negatively in the total delay. Even if this method seems quite crude, we find that it gives a reasonable estimate of the time shifts. Indeed, the computed times differ from the “exact” postKeplerian results reported in Tables 1–3 by only 9% after 33 cycles.
All Tables
Differences (in seconds) between the instants of transit computed with the relativistic (1PN) model and with the Newtonian model, taking as reference point (δt ≡ 0) the passage at periastron just before the transit of Jan. 2010.
Differences (in seconds) between the instants of eclipse computed with the relativistic (1PN) model and with the Newtonian model, taking as reference point () the passage at periastron just after the eclipse of Jan. 2010.
Differences (in seconds) between the time interval t_{tr − ec}(N) between the passage at the minimum approach point during the Nth eclipse and the passage at the minimum approach point during the Nth transit, and the time interval t_{tr − ec}(0) at the reference period, when the eclipse and transit where measured in January 2010 with the Spitzer satellite (Hébrard et al. 2011).
All Figures
Fig. 1. Geometry of the planetary transit: (a) Motion of the planet with respect to the observer, located at infinity in the direction of x of the frame {x, y, z}. The centre of the star is denoted O, P is the periastron at reference time t_{P}, and J_{⋆} is the angular momentum of the star projected on the plane of the sky {y, z}. (b) The transit as seen by the observer and the definition of the points {T_{i}}; here are only shown T_{1} and T_{3}. (c) The trajectory of the planet and the points {T_{i}} (transit) and (eclipse) as seen in the orbital plane {X, Y}. Up to an angle Ω, the observer is in the direction of X. (d) View of the transit from above. The direction (ON) is the line of “nodes”, corresponding to (OX) in (a). 

Open with DEXTER  
In the text 
Fig. 2. Displacement of the trajectory and time shifts of the exoplanet inbetween two approximate successive Keplerian orbits N − 1 and N. The two main contributions δt_{k} and δt_{P} of the (rough) estimation of the time shift are indicated. The point is the point corresponding to T_{N − 1} reported in the next orbit N. 

Open with DEXTER  
In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.