# Radial velocity eclipse mapping of exoplanets

## Abstract

Planetary rotation rates and obliquities provide information regarding the history of planet formation, but have not yet been measured for evolved extrasolar planets. Here we investigate the theoretical and observational perspective of the Rossiter-McLauglin effect during secondary eclipse (RMse) ingress and egress for transiting exoplanets. Near secondary eclipse, when the planet passes behind the parent star, the star sequentially obscures light from the approaching and receding parts of the rotating planetary surface. The temporal block of light emerging from the approaching (blue-shifted) or receding (red- shifted) parts of the planet causes a temporal distortion in the planet’s spectral line profiles resulting in an anomaly in the planet’s radial velocity curve. We demonstrate that the shape and the ratio of the ingress-to-egress radial velocity amplitudes depends on the planetary rotational rate, axial tilt and impact factor (i.e. sky-projected planet spin-orbital alignment). In addition, line asymmetries originating from different layers in the atmosphere of the planet could provide information regarding zonal atmospheric winds and constraints on the hot spot shape for giant irradiated exoplanets. The effect is expected to be most-pronounced at near-infrared wavelengths, where the planet-to-star contrasts are large. We create synthetic near-infrared, high-dispersion spectroscopic data and demonstrate how the sky-projected spin axis orientation and equatorial velocity of the planet can be estimated. We conclude that the RMse effect could be a powerful method to measure exoplanet spins.

###### Subject headings:

planets and satellites: atmospheres, fundamental parameters — stars: planetary systems — techniques: spectroscopic — infrared: planetary systems^{1}

## 1. Introduction

Planetary rotation rates and axial tilts are critical parameters in determining seasonal climate variations (Williams & Kasting, 1997; Cowan et al., 2012) and are relevant to the planet formation and evolution history (Agnor, 1999). By definition the planet rotational rate and axial tilt (e.g. obliquity) are respectively the time it takes for a complete revolution and the angle between the planetary spin angular momentum and the planet’s orbital angular momentum. Planets with obliquity and are considered to have prograde and retrograde rotation, respectively. Planet spins are well-constrained for the Solar system planets, spanning a wide range of rotational rates and axial tilts (Cox & Pilachowski, 2000) and are considered to reflect the unique formation and evolutionary history of each planet (Laskar & Robutel, 1993). Planets accumulate rotational angular momentum from the relative motions of accreted material. The stochastic nature of planetary accretion from planetesimals allows for a random component to the net spin angular momentum. Prograde angular momentum (i.e. spin) could be accumulated by a planet on a circular orbit within a uniform surface density disk of small planetesimals. The planet clears a gap and thus accretes a larger fraction of material from the edges of its accretion zone. Retrograde spins are considered to originate from giant impacts during the early stages of planet formation. Therefore constraints on the planet spins are of high scientific interest, because deviations of the axes of rotation may have been caused by impacts of large bodies during their early history. Even the axes of rotation of the gaseous planets may have been affected by impacts on their rocky cores before these planets accumulated their large atmospheres of hydrogen and helium. Such atmospheres are considered to be accreted hydrodynamically, in flows quite different from those which govern the dynamics of planetesimals and lead to prograde rotation (de Pater & Lissauer, 2001; Faure, 2007).

Exoplanet tidal theory predicts obliquity erosion, on time scales Gyr, for planets on day orbits around normal low-mass stars, i.e. preventing seasonal variations (Heller et al., 2011a, b). Observational methods have been proposed to probe the rotation rates and obliquities for exoplanets from oblateness measurements and variability due to surface inhomogeneity e.g. Hui & Seager (2002); Seager & Hui (2002); Barnes & Fortney (2003); Pallé et al. (2008); Kawahara & Fujii (2010); Fujii & Kawahara (2012), but require precisions in excess of mag.

Currently a rotational rate has been probed for only one extrasolar planet - the young fast rotator Pic b, (Snellen et al., 2014) while axial tilts have not yet been measured. Ground-based high-dispersion spectroscopy () in the near-infrared has recently become successful in characterising the atmospheres of hot Jupiters (Snellen et al., 2010; Brogi et al., 2012; Birkby et al., 2013; Brogi et al., 2014; Schwarz et al., 2015). At high spectral resolution molecular absorption bands are resolved into individual lines allowing their robust identification by line matching with model templates. As the planet orbits its star, the radial component of the planet orbital velocity changes by tens of km s, enabling a discrimination of the Doppler shifted planet spectrum from the steady telluric contamination. The planet signal is then extracted by cross-correlating the data with model spectra obtained by mixing the expected spectroscopically active trace gases in hot-Jupiter atmospheres and assuming a range of vertical temperature pressure profiles.

Kawahara (2012) considered the effect of planetary spin on the planetary radial velocity in dayside spectra of exoplanets, simulated the effect and concluded that planetary radial velocity could be a powerful means for constraining planet spins. In this paper we describe the potential of the Rossiter-McLauglin effect during secondary eclipse (RMse) combined with near-infrared high-dispersion spectroscopy to provide constraints on exoplanet rotation and obliquity.

The rest of the paper is organised as follows. Section 2 describes the rotational effect in transiting exoplanets. In Section 3 we derive the planet radial velocity anomaly due to the RMse effect. We choose the methodology of Ohta et al. (2005) who found an analytic solution for the RM effect at primary transits and show that their solution can equivalently be applied to the problem discussed here. Section 4 details the amplitude and shape of the RMse effect and discusses potential targets. Section 5 discusses the potential application of the method in light of the available and future instrumentation. Finally Section 6 is devoted to our conclusions.

## 2. RM effect at planet secondary eclipse

In its nature, the RM effect arises because of stellar rotation (Rossiter, 1924; McLaughlin, 1924). During a stellar occultation or a planetary transit, portions of the rotating star surface are temporarily obscured, causing the removal of particular radial velocity components from the stellar broadened absorption lines leading to a temporal radial velocity anomaly (Winn et al., 2005). The same effect is expected to take place with a rotating planet at secondary eclipse. When the planet passes behind the parent star the light from the approaching and receding parts of the rotating planetary surface sequentially enter/exit the geometric shadow of the star. The temporal block of light emerging from the approaching (blue-shifted) or receding (red-shifted) parts of the planet causes temporal distortion in the line profiles of the planet’s spectrum leading to an anomaly in the planet radial velocity curve. The shape and amplitude of this anomaly depend on the planet rotation rate, axial tilt and spatial orientation (i.e. planet spin-orbital alignment). Measurements of rotational rates and axial tilts are of high scientific interest as they provide constraints on exoplanet formation and evolution.

Orbital Parameters

& § 3 & Star mass

& § 3 & Planet mass

& § 3 & Orbital period

& Fig. 1 & Semimajor axis

& Fig. 1 & Planet orbital eccentricity

& Fig. 1 & Argument of periastron

& Fig. 1 & Eccentric anomaly

& Fig. 1 & Mean motion

& Fig. 1 & Mean anomaly

& Fig. 1 & time of pericentre passage

& Fig. 1 a & orbital inclination

& Fig. 1 a & True anomaly

& Fig. 1 a & planet to star distance

Internal Parameters of Star and Planet

& Fig. 1 a & Planet spin-to--axis angle

& Fig. 1 a & Planet angular velocity

& Fig. 1 b & Sky-projected spin-orbit angle

& § 3 & Stellar radius

& § 3 & Planet radius

& § 3 & Planet surface velocity,

Mathematical notation

& Fig. 1 a, b & Normal vector to the planet orbit

x & § 3 & Position of the star

x & § 3 & Position of the planet

& § 3 & Star to planet ratio

& Eq 7 & See Fig. 1 c

& Eq 12 & See Fig. 1 c

& Eq 13 & See Fig. 1 c

& Eq 14 & See Fig. 1 c

Note. – see Fig. 1 in Ohta et al. (2005) for all quantities marked with () and are defined as expected.

## 3. Formalism of the effect

To describe quantitatively the radial velocity anomaly caused by the RMse effect we assume a two-body problem with a central star and a planet of masses and , respectively. We refer the reader to Fig. 1 in Ohta et al. (2005) for a schematic illustration of the top view of the planetary orbit and their equations (1) to (7). The orbital velocity of the planet as a function of time with respect to the star, up to as detailed in Murray & Dermott (1999) is

(1) |

with all quantities defined in Table 1.

An eclipse or occultation of a part of the rotating planetary surface causes a time-dependant asymmetry in the absorption/emission line profiles. These asymmetries result in an apparent shift of the central spectral line positions when the lines are unresolved.

To describe the radial velocity anomaly caused by the planet rotation, similar to Ohta et al. (2005), we initially set the coordinate system at the star centre and its -axis to coincide with the observer’s line of sight (Fig. 1 a). The planet position is described with the coordinates , corresponding to the orbit plane position and the planet impact parameter.

For simplicity of the mathematical description of the problem, we choose a reference system () centred on the planet and rotated such that the -axis is parallel to the rotation axis of the planet (i.e. parallel to , see Fig. 1, b) and the rotation axis lies in the - plane. We define an angle between the sky-projected rotational angular velocity and the normal unit vector of the planet orbit, , see Fig. 1, b. This differs from the definition of Ohta et al. (2005), who assume to be the angle between the sky-projected stellar rotation axis and the normal vector of the planetary orbit .

In all calculations we ignore differential rotation of the planet surface as well as motions associated with atmospheric dynamics. A point on the surface of the rotating planet with coordinates will move with a velocity given by

(2) |

where is the angular velocity of the planet. The associated radiation will exhibit a Doppler shift defined as

(3) |

with respect to the observer along the -axis (i.e. the line of sight). We refer the reader to Sec. 3, Equations (14-19) in Ohta et al. (2005) for a derivation of the radial velocity profile for a star and adopt their expression (20) rewritten for the planet:

(4) |

Equation 4 relates the radial velocity change and the line intensity . Fig. 2 illustrates the different cases of the RM effect during secondary eclipse.

We evaluate the integrals assuming uniform model of the planet surface intensity . We ignore the role of planet limb-darkening as our goal is to estimate the first order rotational effect, and leave inclusion of the limb-darkening for future investigations. We also consider the star to be completely optically thick.

At ingress and egress the position of the stellar disc satisfies the relation . In order to simplify the computational task we rotate the coordinates in a time-dependent manner so that the stellar centre is always located along the new -axis, as in Ohta et al. (2005), see Fig. 1 b:

(5) |

The position of the star is given by

(6) |

where is defined as

(7) |

The intensity on the uniform planet surface at () is given by

(8) |

where , . The moments of the intensity then are

(9) |

and

(10) |

The range of the integrals is denoted with and is defined as the star-planet overlapping region (shaded area) in Fig. 1 c:

(11) |

The star and planet circles intersect at (), where

(12) |

(13) |

These calculations are equivalent to the calculations presented in Ohta et al. (2005) and to simplify the final result we introduce

(14) |

Then equations 9 and 10 are analytically integrated as

(15) |

and

(16) |

Finally combining these two results we find the radial velocity anomaly of the planet (during ingress or egress) as a function of the star position ():

(17) |

The final result is equivalent to the result of Ohta et al. (2005) which is expected given the identical overlapping area between the planet and star circles.

## 4. RMse effect amplitude and shape

The shape and amplitude of the RMse effect are illustrated in Fig. 2 for nine representative configurations of the planet spin-orbital alignments (i.e. impact parameter and ), corresponding to prograde rotation. In case of retrograde rotation (i. e. ) the radial velocity curves are inverted. The RMse computation in Fig. 2 assumes a Jupiter-like planet, with at the equator of km s, on a day orbit around a Sun-like star. We choose a period rather longer than the currently typical hot-Jupiter orbital periods because close-in exoplanets are expected to rapidly synchronise their rotation with the orbital period due to tides raised by the star on the planet. The time to spin-down the planet rotation is given by

(18) |

where , , , , , and are the planet’s tidal dissipation factor, radius, mass, rotational angular velocity, and synchronous (or orbital) angular velocity (Goldreich & Soter, 1966; Guillot et al., 1996). Assuming (Correia & Laskar, 2010) and orbital periods of 4 and 20 d, we find and yr, respectively implying that a d orbit required to have spin-down times comparable to the ages of the currently observed transiting hot Jupiters.

Our choice for the value of is driven by the correlation of the equatorial rotational velocities of solar system planets and their masses, suggesting that more massive planets rotate faster (Hughes, 2003). Currently a rotational rate has been constrained only for the massive and young planet Pic b ( km s), which is in line with the spin velocity-mass relation of the solar system (Snellen et al., 2014).

The maximum amplitude of the RMse effect is expected for central eclipses (i. e. ) and . In the case considered in Fig. 2, the amplitude is . Assuming the of Pic b we find rotation effect with an amplitude of . No effect is expected in cases where equal portions red- and blue-shifted portions of the planet are eclipsed (i.e. and ). This is also the case where the planet rotation axis is normal to the plane of the sky.

It should be pointed out that in all cases that include , is degenerate with , i.e. various combinations of these two parameters can produce the same radial velocity amplitude. Importantly, it is the amplitude difference of the ingress and egress RMse, their signs and the shape of the radial velocity curve that constrain the planet obliquity and break this degeneracy.

## 5. An observational perspective for RMse

We produce mock data to illustrate the planetary RMse effect assuming an 8 and 40 m class telescopes, equipped with a near-infrared, high spectral resolution (i.e. ) spectrometer. We assume that the planet orbit is constrained with enough precision to be subtracted prior to the search for the RMse effect. We calibrate our simulation adopting a precision of km s as achieved by Birkby et al. (2013) at m for water detection in the day side atmosphere of HD 189733b. We scale this precision by the square root of the number of collected spectra (i.e. 48), the ratio of the employed spectral coverage (assuming a hypothetical high-resolution spectrometer with a wavelength coverage of 400 nm) and a square root of the ratio of the target planet-to-star flux ratio and the planet-to-star flux ratio of HD 189733b, i.e. from Birkby et al. (2013). This factor plays a critical role as it accounts for the strength of the planet signal, which is determined by the planet to star flux contrast. The contrast decreased from to when increasing the planet orbital periods from day to day. We also factor in the sampling rate with the telescope diameter and the brightness of the target host star. An account for the change of the planet flux during a secondary eclipse is also incorporated by factoring the radial velocity uncertainties by square root of the flux of a secondary eclipse Mandel & Agol (2002) model with unity out-of eclipse baseline and zero in-eclipse flux.

Fig. 3 shows mock data generated assuming a 40 m-class telescope and a hypothetical planet host star of brightness mag (i.e. similar to the currently brightest transiting hot-Jupiter host star HD 189733). Fig. 3 shows the results for a planet with the physical properties of WASP-19 b, which is considered to be tidally locked. We assumed a prograde () and a retrograde () rotation with km s. We perform a RMse fit to the mock data utilising a Levenberg-Marquardt least-squares minimization and estimate and , (Markwardt, 2009). We find the planet signal to be well-detected (-) when combining nine secondary eclipses.

We also explored the potential of an 40 m-class telescope to detect the RMse effect in the currently brightest transiting exoplanet host-star HD 189733 b. Assuming a synchronised planet rotation (i.e. RMse amplitude of km s) and we find that secondary eclipses need to be co-added to detect the spin of HD 189733 b.

Finally, we investigated the potential of an 8 m-class telescope to detect the RMse effect in K=5.5 host star with a WASP-19b-like planet and synchronised rotation. We find a number of secondary eclipses will have to be added in order to detect the RMse effect at - confidence.

## 6. Discussion and conclusions

We have discussed the Doppler signature of a spinning exoplanet during secondary eclipse and derived a relation for the radial velocity as a function of the planet’s rotational rate and axial tilt (Eq. 17). The RMse effect complements the tool-box offered by transiting exoplanets providing a proxy to exoplanet spins. We showed that the formalism of the RMse is equivalent to the formalism developed by Ohta et al. (2005) to describe the RM effect caused by a planet transiting its star just from a different perspective. We note that Hirano et al. (2011) derived an improved accuracy radial velocity anomaly curves for the RM effect during primary transits. However, the solution of Ohta et al. (2005) provides a conservative estimate of the radial velocity amplitude and is precise enough to illustrate the RMse effect.

We acknowledge that our assumption of a rigid rotating non-limb darkened planetary surface (i.e. ignoring differential rotation and atmospheric dynamics) is a crude approximation, but it illustrates the expected RMse radial velocity anomaly in the simplest case. The radial velocity originating from the RMse effect is expected to be degenerate with atmospheric dynamics. Exoplanets are expected to have atmospheric circulation with increasing dynamics at shorter orbital distances (Kataria et al., 2015). Therefore it is expected that the radial velocity curves of most close-in exoplanets may significantly differ from that assumed here, with super-rotating winds expected to enhance the radial velocity signals by factor of a few. The RMse effect presents the opportunity to probe atmospheric winds with altitude, when wide coverage high dispersion spectroscopy becomes capable of providing per-point precisions in the m s domain. In such cases it could be possible to measure longitudinally-integrated radial velocities and map the surface of a planet during the secondary eclipse ingress/egress phases, providing constraints on the hot-spot position. We postpone the derivation of the RMse effect in these cases for future theoretical global circulation model investigations and identify the top three planets expected to exhibit the strongest RMse effect assuming synchronised rotation. Applying Eq. 17 (assuming ) we find values of km s for targets WASP-103 b, WASP-12 b and WASP-19 b, respectively.

Kawahara (2012) demonstrated the rotational signature of an exoplanet in the planet radial velocity. The RMse effect discussed in this study is of similar nature. However the geometric shadow of the star rather than the changing planet phase is the factor removing velocity components from the visible planet surface. Importantly, because the RMse effect occurs only during ingress and egress, its detection requires significantly less observing time compared to the case in Kawahara (2012). RMse therefore could be much-well suited for planets with longer orbital periods.

We also demonstrated that to be put into practice the RMse effect requires exoplanet host stars brighter than the currently known exoplanet hosts, i.e. brighter than mag and large aperture telescopes (i e. 40 m). The strongest RMse effect for evolved stars is expected for non-synchronized rotation, i.e. planets on longer than 20 day orbits, with a typical km s or more. However, such planets exhibit lower temperatures (i.e. K), compared to the typical hot Jupiters giving a small planet-to-star flux contrasts (i.e. or lower at m). The detection of the RMse effect in such cases would be an extremely challenging if not impossible task, even for the upcoming 40 m class telescopes.

An interesting case could be young planets which still have not been synchronised. Observations using high-dispersion, near-infrared spectroscopy have currently constrained the rotational rate of the only one such case ( Pic b) from rotationally broadened absorption lines (Snellen et al., 2014). An important opportunity could exist if the orbit orientation of this planet allows transits and secondary eclipses, because the large brightness of the host star (i. e. mag) and the fast spin of the planet (i. e. km s) would produce a strong and detectable RMse signal. Another interesting opportunity could be offered by transiting brown dwarfs with orbital periods larger than 10-20 days. Such systems assume high flux contrasts between the star and the brown dwarf in the near-infrared which in case of a bright host star (i. e. mag) could provide an opportunity to probe the spins and latitudinal radial velocity maps of these objects.

In the near future, the Next Generation of Transit Surveys (NGTS), the Transiting Survey Satellite (TESS) and the Planetary Transits and Oscillations of stars (PLATO) projects are expected to significantly expand the sample of known planets hosted by bright stars (i.e. brighter than mag) and hence to provide more targets suitable for detection of the RMse effect.

The hypothetical radial velocity data in Section 5 assumes retrieval of the planet signal via cross-correlation with a model planet spectrum similar to Snellen et al. (2010); Brogi et al. (2012); Birkby et al. (2013); de Kok et al. (2013); Schwarz et al. (2015). Although currently the cross-correlation technique has not been demonstrated to provide per spectrum radial velocity measurements, in their Fig. 1, Snellen et al. (2010) presented results for the planet geocentric radial velocity as a function of the planet orbital phase at a high significance (-). It is expected that large-aperture near-future telescopes (i.e. 30-40 m) could provide even higher significance signals and planet radial velocity measurements from multiple spectra. In addition, cloud-free exoplanet atmospheres are expected to result in stronger cross-correlation signals compared to hazy or cloudy atmospheres for the same stellar and planetary physical properties. Cloud-free atmospheres are expected to have well-pronounced lines, i.e. non-muted absorption features which are expected for hazy atmospheres.

Stellar activity and pulsations are known to produce radial velocity variations that could in some cases mimic those induced by the orbital motion of exoplanets. This may lead to misinterpretations of radial velocity variations, especially when those variations have periods less or equal to the star rotational periods (Lagrange et al., 2010). Therefore, it is important to consider whether stellar activity, for instance the distortion to the stellar spectral lines induced by starspots on a rotating stellar surface, could induce residual signal misinterpreted as exoplanet signal. This is especially pertinent for co-added observations obtained over several nights.

A typical RMse observation would last for h on each of the ingress and egress phases of a WASP-19b-like planet (see Fig. 3). This time interval is significantly shorter than the rotation periods of the typical planet host stars ranging from 10 to 40 days (Paz-Chinchón et al., 2015). Furthermore, starspots are expected to produce sine-wave radial velocity variations (i.e. different from the RMse shape) with amplitudes of a few hundred m s in the optical and are expected to reduce by factor of at least a few at near infrared wavelengths, where the contrast between the photosphere and cool star spots (e.g. K) is significantly reduced.

Thus, if present, residual stellar signal caused by starspots would be a factor of a few tens smaller than the the expected 1-2 km s amplitude of the RMse effect for the closest synchronized planets and even more for non-synchronized planets on longer than 20 day period orbits. Only in the case of most evolved planets (i.e. typical hot Jupiters), if synchronized, with orbital/rotational periods from 4-5 days to day could a residual introduce larger radial velocity scatter when multiple observations are combined. However, such planets are expected to have RMse amplitudes much smaller than 1 km s and would be difficult targets by definition.

In the near infrared the planet spectrum is dominated by thermal radiation and is not expected to contain information from distortions caused by starspots. In case that the planet spectrum also contains a non-negligible component from a reflected star spectrum then the planet spectral lines will have distortions caused by the star spots. This would make the planet spectrum different and more difficult to cross correlate with model spectra. However, such an effect is not expected to mimic the RMse.

In conclusion, we have demonstrated the RMse effect and its potential to constrain planet rotational rates and sky projected spin-orbital alignments. We derived the radial velocity curve caused by the RM anomaly and estimated the amplitude of the effect for a hot Jupiter. Finally we discussed the prospects for detecting the effect and constraining planet spins and axial tilts from the current and upcoming instrumentation.

### Footnotes

- slugcomment: Received 2015 May 20; accepted 2015 June 12;

### References

- Agnor, C. B., Canup, R. M., & Levison, H. F. 1999, Icarus, 142, 219
- Barnes, J. W., & Fortney, J. J. 2003, ApJ, 588, 545
- Birkby, J. L., de Kok, R. J., Brogi, M., et al. 2013, MNRAS, 436, L35
- Brogi, M., de Kok, R. J., Birkby, J. L., Schwarz, H., & Snellen, I. A. G. 2014, A&A, 565, A124
- Brogi, M., Snellen, I. A. G., de Kok, R. J., et al. 2012, Nature, 486, 502
- Correia, A. C. M., & Laskar, J. 2010, Exoplanets, 239
- Cowan, N. B., Voigt, A., & Abbot, D. S. 2012, ApJ, 757, 80
- Cox, A. N., & Pilachowski, C. A. 2000, Physics Today, 53, 77
- de Kok, R. J., Brogi, M., Snellen, I. A. G., et al. 2013, A&A, 554, A82
- de Pater, I., & Lissauer, J. J. 2001, Planetary Sciences, by Imke de Pater and Jack J. Lissauer, pp. 544. ISBN 0521482194. Cambridge, UK: Cambridge University Press, December 2001.,
- Faure, G., Mensing, Teresa M. 2007, Introduction to Planetary Science: The Geological Perspective, by G. Faure and T.M. Mensing. Berlin: Springer, 2007.,
- Fujii, Y., & Kawahara, H. 2012, ApJ, 755, 101
- Goldreich, P., & Soter, S. 1966, Icarus, 5, 375
- Guillot, T., Burrows, A., Hubbard, W. B., Lunine, J. I., & Saumon, D. 1996, ApJ, 459, L35
- Heller, R., Barnes, R., & Leconte, J. 2011, Origins of Life and Evolution of the Biosphere, 41, 539
- Heller, R., Leconte, J., & Barnes, R. 2011, A&A, 528, A27
- Hirano, T., Suto, Y., Winn, J. N., et al. 2011, ApJ, 742, 69
- Hughes, D. W. 2003, Planet. Space Sci., 51, 517
- Hui, L., & Seager, S. 2002, ApJ, 572, 540
- Kataria, T., Showman, A. P., Fortney, J. J., et al. 2015, ApJ, 801, 86
- Kawahara, H. 2012, ApJ, 760, L13
- Kawahara, H., & Fujii, Y. 2010, ApJ, 720, 1333
- Lagrange, A.-M., Desort, M., & Meunier, N. 2010, A&A, 512, A38
- Laskar, J., & Robutel, P. 1993, Nature, 361, 608
- Mandel, K., & Agol, E. 2002, ApJ, 580, L171
- Markwardt, C. B. 2009, Astronomical Data Analysis Software and Systems XVIII, 411, 251
- McLaughlin, D. B. 1924, ApJ, 60, 22
- Murray, C. D., & Dermott, S. F. 1999, Solar system dynamics by Murray, C. D., 1999,
- Ohta, Y., Taruya, A., & Suto, Y. 2005, ApJ, 622, 1118
- Pallé, E., Ford, E. B., Seager, S., Montañés-Rodríguez, P., & Vazquez, M. 2008, ApJ, 676, 1319
- Paz-Chinchón, F., Leão, I. C., Bravo, J. P., et al. 2015, ApJ, 803, 69
- Rossiter, R. A. 1924, ApJ, 60, 15
- Schwarz, H., Brogi, M., de Kok, R., Birkby, J., & Snellen, I. 2015, A&A, 576, A111
- Seager, S., & Hui, L. 2002, ApJ, 574, 1004
- Snellen, I. A. G., de Kok, R. J., de Mooij, E. J. W., & Albrecht, S. 2010, Nature, 465, 1049
- Snellen, I. A. G., Brandl, B. R., de Kok, R. J., et al. 2014, Nature, 509, 63
- Williams, D. M., & Kasting, J. F. 1997, Icarus, 129, 254
- Winn, J. N., Noyes, R. W., Holman, M. J., et al. 2005, ApJ, 631, 1215