Searching For Rapid Orbital Decay of WASP-18b
Ashlee N. Wilkins, Laetitia Delrez, Adrian J. Barker, Drake Deming,, Douglas Hamilton, Michael Gillon, and Emmanuel Jehin

TL;DR
This study investigates the orbital decay of exoplanet WASP-18b to measure the star's tidal dissipation efficiency, finding no rapid decay and suggesting the star's Q' is higher than previously inferred, indicating less tidal dissipation.
Contribution
The paper provides the first nine-year observational analysis of WASP-18b's orbital decay, constraining the star's tidal quality factor Q' to be greater than 10^6, challenging prior expectations.
Findings
No evidence of rapid orbital decay over nine years
Q' for the host star is at least 10^6 at 95% confidence
Supports the idea that F-stars are less tidally dissipative than solar-type stars
Abstract
The WASP-18 system, with its massive and extremely close-in planet, WASP-18b (M = 10.3M, a = 0.02 AU, P = 22.6 hours), is one of the best known exoplanet laboratories to directly measure Q', the modified tidal quality factor and proxy for efficiency of tidal dissipation, of the host star. Previous analysis predicted a rapid orbital decay of the planet toward its host star that should be measurable on the time scale of a few years, if the star is as dissipative as is inferred from the circularization of close-in solar-type binary stars. We have compiled published transit and secondary eclipse timing (as observed by WASP, TRAPPIST, and Spitzer) with more recent unpublished light curves (as observed by TRAPPIST and HST) with coverage spanning nine years. We find no signature of a rapid decay. We conclude that the absence of rapid orbital decay most likely derives from Q' being…
| Parameter | Value | Average | Reference |
|---|---|---|---|
| The Star: WASP-18 | |||
| 1.29 0.16 | Doyle et al. 2013 | ||
| Radius (R⊙) | 1.22 0.11 | ||
| 1.15 0.02 | Bonfanti et al. 2016 | ||
| 6400 100 | Hellier et al. 2009 | ||
| Teff (K) | 6400 75 | 6322 72 | Doyle et al. 2013 |
| 6167 7 | Bonfanti et al. 2016 | ||
| 4.4 0.15 | Hellier et al. 2009 | ||
| log | 4.32 0.09 | 4.32 0.10 | Doyle et al. 2013 |
| 4.39 0.01 | Bonfanti et al. 2016 | ||
| < 2.0 | Southworth et al. 2009 | ||
| Age (Gyr) | 0.5 - 1.5 | Hellier et al. 2009 | |
| 0.9 0.2 | Bonfanti et al. 2016 | ||
| 1.281 | Southworth et al. 2009 | ||
| M∗ (M⊙) | 1.24 0.04 | 1.25 0.04 | Triaud et al. 2010 |
| 1.22 0.03 | Enoch et al. 2010 | ||
| The Planet: WASP-18b | |||
| P (days) | 0.94145299 8.7 | Hellier et al. 2009 | |
| 0.02047 | Southworth et al. 2009 | ||
| a (AU) | 0.02034 | ||
| 0.02020 | Triaud et al. 2010 | ||
| 86 2.5 | Hellier et al. 2009 | ||
| (∘) | 83.3 | ||
| 80.6 | Triaud et al. 2010 | ||
| 10.43 | Southworth et al. 2009 | ||
| Mp (MJ) | 10.27 | ||
| 10.11 | Triaud et al. 2010 | ||
| Facility | Date | Original Reference(s) | Orbit | BJD (TDB) |
|---|---|---|---|---|
| WASP | May - December 2006 | Hellier et al. (2009) | 0 | 2454221.481630.00038 |
| Spitzer | December 20, 2008 | Nymeyer et al. (2011); Maxted et al. (2013) | 636.5 | 2454820.71680.0007 |
| Spitzer | December 24, 2008 | Nymeyer et al. (2011); Maxted et al. (2013) | 640.5 | 2454824.48150.0006 |
| Warm Spitzer | January 23, 2010 | Maxted et al. (2013) | 1061.5 | 2455220.83370.0006 |
| Warm Spitzer | January 24, 2010 | Maxted et al. (2013) | 1062 | 2455221.30420.0001 |
| Warm Spitzer | August 23, 2010 | Maxted et al. (2013) | 1285.5 | 2455431.71910.0003 |
| Warm Spitzer | August 24, 2010 | Maxted et al. (2013) | 1286 | 2455432.18970.0001 |
| TRAPPIST | September 30, 2010 | Maxted et al. (2013) | 1327 | 2455470.78850.00040 |
| TRAPPIST | October 2, 2010 | Maxted et al. (2013) | 1330 | 2455473.61440.00090 |
| TRAPPIST | December 23, 2010 | Maxted et al. (2013) | 1416 | 2455554.57860.00050 |
| TRAPPIST | January 8, 2011 | Maxted et al. (2013) | 1433 | 2455570.5840 |
| TRAPPIST | November 11, 2011 | Maxted et al. (2013) | 1758 | 2455876.55590.0013 |
| HST | April 22, 2014 | This work | 2840.5 | 2456895.67730.0006 |
| HST | April 22, 2014 | This work | 2841 | 2456896.14780.0008 |
| TRAPPIST | August 20, 2015 | This work | 3223 | 2457255.7832 |
| TRAPPIST | October 21, 2015 | This work | 3291 | 2457319.8010 |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Searching For Rapid Orbital Decay of WASP-18b
Ashlee N. Wilkins11affiliation: Department of Astronomy, University of Maryland, College Park, MD 20742-2421, [email protected] , Laetitia Delrez22affiliation: Cavendish Laboratory, University of Cambridge, J J Thomson Avenue, Cambridge CB3 0HE , Adrian J. Barker33affiliation: Department of Applied Mathematics, School of Mathematics, University of Leeds, Leeds, LS2 9JT, UK , Drake Deming11affiliation: Department of Astronomy, University of Maryland, College Park, MD 20742-2421, [email protected] ,
Douglas Hamilton11affiliation: Department of Astronomy, University of Maryland, College Park, MD 20742-2421, [email protected] , Michael Gillon44affiliation: Space Sciences, Technologies and Astrophysics Research (STAR) Institute, Université de Liège, allée du 6 Août 19C, 4000 Liège, Belgium , and Emmanuel Jehin44affiliation: Space Sciences, Technologies and Astrophysics Research (STAR) Institute, Université de Liège, allée du 6 Août 19C, 4000 Liège, Belgium
Abstract
The WASP-18 system, with its massive and extremely close-in planet, WASP-18b (Mp = 10.3MJ, a = 0.02 AU, P = 22.6 hours), is one of the best known exoplanet laboratories to directly measure Q’, the modified tidal quality factor and proxy for efficiency of tidal dissipation, of the host star. Previous analysis predicted a rapid orbital decay of the planet toward its host star that should be measurable on the time scale of a few years, if the star is as dissipative as is inferred from the circularization of close-in solar-type binary stars. We have compiled published transit and secondary eclipse timing (as observed by WASP, TRAPPIST, and Spitzer) with more recent unpublished light curves (as observed by TRAPPIST and HST) with coverage spanning nine years. We find no signature of a rapid decay. We conclude that the absence of rapid orbital decay most likely derives from Q’ being larger than was inferred from solar-type stars, and find that Q’ 1, at 95 % confidence; this supports previous work suggesting that F-stars, with their convective cores and thin convective envelopes, are significantly less tidally dissipative than solar-type stars, with radiative cores and large convective envelopes.
stars: individual (WASP-18) – planets and satellites: atmospheres – techniques: photometric – techniques: spectroscopic
††slugcomment: Draft for the Astrophysical Journal
1 Introduction
The discovery of WASP-18b (Hellier et al., 2009; Southworth et al., 2009), with its large mass (10.3 MJ) and small orbit (22.6 hours) (see Table LABEL:tab:system_details for other parameters), elicited one primary question: how could it exist? A planet of that mass and proximity should raise a substantial tidal distortion (tidal bulge) in the central star. Because the star is not a perfectly elastic body, and because the planet orbits more quickly than the star rotates, the tidal bulge would lag behind the planet, causing the planet’s orbital motion to accelerate and the orbit to shrink (Goldreich & Soter, 1966). Hellier et al. (2009) calculated a 0.65 Myr future lifetime for the planet, assuming that the star is as dissipative as is inferred from the circularization of close solar-type binary stars (Meibom & Mathieu, 2005; Ogilvie & Lin, 2007). The estimated age of the star is a few hundred million years (Bonfanti et al., 2016; Hellier et al., 2009) to 2 Gyr (Southworth et al., 2009); finding a planet with a lifetime that is such a small fraction of the system’s age is extremely improbable. Hamilton (2009) discusses several alternative explanations, ranging from an overestimation of the decay rate (due to unmodeled nuances of tidal physics leading to an underestimation of the tidal Q’ parameter) to a non-tidal mechanism holding the planet in place (e.g. influence of another body in the system). Barker & Ogilvie (2009) investigate the efficiency of tidal dissipation in the convective envelopes of F-stars, which have both convective cores and convective envelopes; G stars, on which most studies of exoplanetary tidal decay focus (e.g. Jackson et al. 2009; Birkby et al. 2014), have radiative cores and thicker convective envelopes. The Barker & Ogilvie (2009) calculations reveal that tidal dissipation within F stars is generally much less efficient than within G-stars, and therefore that planetary tidal decay around stars like WASP-18 would be imperceptibly low over a decadel timespan (Barker & Ogilvie, 2010; Barker, 2011). If, however, tidal dissipation within WASP-18 behaved as is usually inferred for solar-type stars, Birkby et al. (2014) predict that its transit should occur progressively earlier at each observation, accumulating to a measurable shift of nearly six minutes over ten years. This is the largest predicted shift of any planet, making the WASP-18 system possibly the best known laboratory for direct measurements of the stellar tidal Q’ parameter. Maciejewski et al. (2016) potentially measured the tidal decay of WASP-12b, but Hoyer et al. (2016) ruled out the orbital decay of WASP-43b proposed by Jiang et al. (2016).
In this Letter, we bring together published measurements of transit and secondary eclipse timing from discovery (Hellier et al., 2009), Spitzer (Nymeyer et al., 2011; Maxted et al., 2013), and ground-based TRAPPIST (Maxted et al., 2013) observations, and new analyses of unpublished archival (HST), and recent TRAPPIST data. We place strong limits on the maximum rate of the system’s orbital decay and discuss the implications.
2 New Observations
2.1 TRAPPIST
The TRAnsiting Planets and PlanetesImals Small Telescope - South (TRAPPIST-S, Jehin et al. 2011; Gillon et al. 2011) is a ground-based, 60-cm robotic telescope based at the La Silla Observatory used to study both exoplanets and small bodies in the Solar System. TRAPPIST observed two WASP-18b photometric transits in the fall of 2015 in the broad-band Sloan-z filter, centered at 0.9134 m.
2.2 Hubble Space Telescope
The Hubble Space Telescope (HST) observed WASP-18b in 2014 in spatial scan mode (Deming et al., 2013) over its full phase (PID 13467, PI Bean), including one full transit, one full secondary eclipse, and one extra eclipse ingress. All observations were made with the Wide Field Camera 3 (WFC3) G141 infrared grism, covering 1.1 - 1.7 m. While the primary deliverable from such observations is the spectrum, we sum over wavelength to extract a photometric light curve. To maximize observing efficiency, the scan reverses direction, rather than taking the time to reset to the starting point, at the end of each scan. This introduces a non-constant offset requiring separate analysis of the forward and reverse scans.
3 Analysis: Deriving the new White Light Curves
Table LABEL:tab:all_obs includes the transit and secondary eclipse times used in this analysis. We describe here how we generated white light curves and transit fits to the new TRAPPIST and HST data.
3.1 TRAPPIST Light Curves
We reduce our TRAPPIST data in the methods described by Gillon et al. (2012). We calculated the best-fit transit curve for each observation using the TRAPPIST MCMC procedure (Gillon et al. 2009 and references therein), executing the Mandel & Agol (2002) algorithm to find the new best-fit light curve parameters. We generated the curve plotted in Figure 1 with the BATMAN procedure (Kreidberg, 2015), given those orbital parameters.
3.2 HST White Light Curves
As has been studied extensively (e.g. Sing et al. 2016), the HST WFC3 camera, while improved over its predecessor NICMOS, has persistent systematic errors that seem to be a function of incident flux (Wilkins et al., 2014), with three distinctive effects: a visit-long ramp, an orbit-long ramp, and a “hook” within orbits (Berta et al., 2012; Wilkins et al., 2014). We reduce the WFC3 data and mitigate systematics in a modified divide-oot method – a method of averaging out all three systematic effects (Wilkins et al., 2014; Deming et al., 2013), including the correction to the STScI wavelength calibrations found in Wilkins et al. (2014). To fit the transit, we use the non-linear, fourth-order limb darkening coefficients from Claret (2000) in the Mandel & Agol (2002) light curve models, and derive “prayer-bead” error bars as in Gillon et al. (2009). To fit the the secondary eclipse, we use the same procedure in the limit of no limb darkening, such that the shape is that of a trapezoid. We analyze the forward and reverse scans independently, as mentioned in § 2, due to a non-linear offset between the two; the final timing results agree and are thus shown as an average in the table.
4 Results: Transit Timing Evolution over Nine Years
We have compiled all published transit and secondary eclipse observations of WASP-18b and added them to the new observations obtained by HST and TRAPPIST to produce a data set spanning more than nine years. The full data set is found in Table LABEL:tab:all_obs. Published light curve solutions came from four observing campaigns:
WASP: The Wide-Angle Search for Planets (WASP) Project (Pollacco et al., 2006) announced the discovery and initial orbital solution of WASP-18b as observed in transit by the WASP-South Survey and in radial velocity with the CORALIE spectrograph (Hellier et al., 2009), and confirmed with the Danish 1.5m telescope at ESO (Southworth et al., 2009). The Southworth et al. (2009) ephemeris was later found to be erroneous (Southworth et al., 2010); we use only the Hellier et al. (2009) ephemeris.
Spitzer: Nymeyer et al. (2011) observed two secondary eclipses of WASP-18b via the Spitzer Exoplanet Target of Opportunity Program with the Infrared Array Camera (IRAC, PID 50517). The first secondary eclipse was observed in the 3.6 m and 5.8 m channels on December 20, 2008, the second in the 4.5 m and 8.0 m channels on December 24, 2008. Maxted et al. (2013) reanalyzed the Nymeyer et al. (2011) points.
Warm Spitzer: Maxted et al. (2013) observed two full phases of WASP-18b’s orbit with warm Spitzer, one with the 3.6 m channel on January 23, 2010, and the other with the 4.5 m channel on August 23, 2010.
TRAPPIST: In addition to the unpublished, new transit curves presented as part of this work, TRAPPIST also observed WASP-18b five times in transit in late 2010 and early 2011, also in the Sloan-z’ filter Maxted et al. (2013).
To search for tidal decay, we study the correlation between the number of orbits since discovery ephemeris and transit (or eclipse) arrival time. In the case of no orbital evolution, this correlation would be linear, and the slope of the line would be the planetary orbital period. We allow for the possibility of decay by including a second-order term that is dependent on the rate of any orbital evolution. We first perform a multivariate linear regression and find a plausible fit (). To explore the trade-off between the linear and quadratic terms of the fits, we also perform a Markov Chain Monte Carlo (MCMC) quadratic fit using (Foreman-Mackey et al., 2013); the results of both fits, which are in excellent agreement, are in Figure 2. With , we find the period P = days, in agreement with the Hellier et al. 2009 P = 0.94145299 8.7 days. If WASP-18 were as tidally dissipative as is inferred from the circularization of solar-type close binary stars, there should be a definitive deviation from linear behavior, i.e., the quadratic term should be nonzero. We measure an upper limit for the magnitude of the quadratic term, and we therefore find no confirmation of rapid tidal decay for the WASP-18 system. Indeed, as discussed in the next section, we should not have expected to find evidence of rapid decay.
5 Discussion: Implications of the Absence of Rapid Tidal Decay
suggested by Hellier et al. (2009); Birkby et al. (2014), we turn instead to the predictions of Barker & Ogilvie (2009, 2010); Barker (2011); Lanza et al. (2011). We first briefly review the discussion of these predictions as they apply to WASP-18, and then calculate a constraint on the Q’ of WASP-18.
5.1 Tidal Dissipation in G vs. F Stars
Tides raised within a central star by a planetary companion are dissipated within the star and angular momentum is transferred between the stellar spin and planetary orbit in the process (e.g. Ogilvie 2014). For short-period planets (orbiting sub-synchronously rotating stars), like WASP-18b, that also have approximately circular orbits, tidal dissipation in the star causes the planet to lose angular momentum and spiral inward, because the tidal bulge raised in the star lags the planet when the planet’s orbital period is less than the star’s rotational period (i.e. Porb < Prot). This is the opposite of the Earth-Moon system, in which the Moon recedes from the Earth because the bulge leads the Moon (since Porb >Prot). The rate of change of the orbit depends on the efficiency of tidal dissipation within the host star; this is where stellar structure becomes important.
The tide in the star is often decomposed into two contributions: an equilibrium tide and a dynamical tide (e.g. Zahn 1977). Dissipation of both components is expected to become less efficient in stars slightly more massive than the Sun (i.e. F stars). While we often generalize Sun-like stars (typically defined as 0.5M*⊙* M*∗* 1.3 M*⊙) to have radiative cores and convective envelopes and more massive stars to have the opposite, development of convective cores and radiative envelopes is actually a continuum. WASP-18, for example, is a 1.2-M⊙* F6 star, and, according to MESA stellar structure models (Paxton et al., 2011), should have a convective core within the innermost 6 % of the stellar radius, and a convective envelope in the outer 15 %; it is therefore intermediate between an solar-mass and high-mass star. For tidal dissipation, therefore, an F star like WASP-18 is not ”sun-like.”
We quantify the efficiency of tidal dissipation using the tidal quality factor, Q, defined as (Goldreich, 1963):
[TABLE]
where is the maximum energy stored in the tidal bulge and , intrinsically negative, is the energy dissipated in one tidal period. We use the modified Q (i.e. Q’) convention throughout this Letter:
[TABLE]
where is the tidal Love number. Q’ is almost certainly not a single constant number for all stars (even of the same spectral type), but is instead a complicated function of the stellar mass, structure, rotation, and tidal periods, as well as the planetary properties (e.g. Ogilvie 2014). Q’ is the Q of an equivalent homogeneous body ( = 3/2). A large Q’ corresponds to weak or inefficient tidal dissipation, and a smaller Q’ corresponds to strong or efficient dissipation. We investigate here the Q’ of the star (WASP-18), not the planet (WASP-18b); the planet’s tidal Q’ is relevant for its own tidal evolution, and leads to synchronization of its rotation and circularization of its orbit.
The equilibrium tide is dissipated within the convective envelope of the star by the effective viscosity of convective turbulence (Zahn, 1966; Goldreich & Nicholson, 1977); however, the effective viscosity may be significantly reduced in the case of a short-period planet (Penev & Sasselov, 2011; Ogilvie & Lesur, 2012). In addition, in F stars, the outer convection zone is thin and of very low mass, so it is expected to be much less dissipative than in G-stars; the effective tidal Q’ could be as high as 1011 (Barker & Ogilvie, 2009) for a star like WASP-18 at the tidal frequencies of interest. Dissipation in the convective core of an F star is also likely weak (e.g. Zahn 1977).
The dynamical tide primarily consists of internal gravity (g-mode) waves that are tidally excited at the convective-envelope-radiative-core boundary and propagate inwards to the center of the star. These waves are thought to be damped by radiative diffusion or nonlinear effects. If they can reach the center, they become geometrically focused and, if the planet exciting them is sufficiently massive – like WASP-18b – they may reach sufficiently large amplitudes such that they break, leading to significantly enhanced tidal dissipation (Goodman & Dickson, 1998; Ogilvie & Lin, 2007; Barker & Ogilvie, 2010; Barker, 2011). This process deposits angular momentum into the star, thereby removing angular momentum from the planet’s orbit; the star’s rotation gets faster (“spin-up”), while the planet’s orbit shrinks. If WASP-18 were Sun-like, WASP-18b would be sufficiently massive to cause wave breaking, and we would expect the planet to rapidly spiral into its star. However, in the case of an F-star like WASP-18, the convective core prevents the tidally-excited gravity waves from reaching the center where they would be focused, so that they may never reach such large amplitudes to break, though they may be subject to weaker nonlinear effects (eg. Barker & Ogilvie 2011; Weinberg et al. 2012; Essick & Weinberg 2016). The dissipation would be significantly reduced, save for select resonant tidal frequencies, so that we would expect the planet to remain in the orbit in which it was discovered (Barker & Ogilvie, 2009). Furthermore, the lingering thin outer convective envelope in an F star of WASP-18’s mass would inhibit radiative damping of the waves near the top of the radiative zone (relative to more massive A-stars). The dissipation that Valsecchi & Rasio (2014) find for WASP-71 may be moderately higher than we would expect in WASP-18 precisely because it is a more massive (1.5 versus 1.2 ) star, and therefore has a thinner outer convective envelope than WASP-18, but what they obtain is still very weak.
Were a resonance present, Q’ could indeed be very low, and therefore the star could be quite dissipative. However, the above arguments and those of, e.g. Lanza et al. (2011); Barker & Ogilvie (2009), support a high-Q’, generally minimally-dissipative scenario for a star like WASP-18.
5.2 Estimating the Tidal Q’ for WASP-18
When a planet transits its host star, we have a convenient time point from which to measure any changes in the orbit, which we infer through a shift in the transit (or secondary eclipse) arrival time. Birkby et al. (2014) show that the expected shift can be reduced to:
[TABLE]
where is the planet-to-star mass ratio (for WASP-18, = 0.007843), is the stellar-radius-to-semi-major-axis ratio (for WASP-18, = 0.2789), T is the elapsed time, and P is the orbital period of the planet. Therefore, in a quadratic fit of the form
[TABLE]
where the linear coefficient corresponds to the period of the planet’s orbit, the quadratic term is defined by Equation 3. Rearranging, we find that Q’ depends on the quadratic coefficient as:
[TABLE]
We fit for the coefficients in Equation 4, and thus the period and Q’, as discussed in § 4 and shown in Figure 2.
Equation 3 makes clear that a planet must be close-in and massive (relative to the radius and mass of the host star), its orbital period must be short, and it must orbit a star with a favorable , in order to produce any discernible shift in time. Currently, in the NASA Exoplanet Archive, only eight confirmed planets have both masses larger than 1.0 and orbital periods of roughly one day or less. The addition of recently-announced KELT-16b (Oberst et al., 2016) makes nine. Of those, one is around a pulsar and four (WASP-18b, KELT-16b, WASP-12b, and WASP-103b) are around stars more massive than 1.2 and therefore likely possessing convective cores that preclude tidal wave breaking at the center. Of the remaining four, WASP-43b has already demonstrated no rapid tidal decay (Hoyer et al., 2016), but WASP-19, WTS-2, and K2-22 all orbit stars less massive than the sun, and may be reasonable testbeds for dissipation within a star with a larger convective envelope and a smaller radiative core; Birkby et al. (2014) has already suggested that WTS-2 should have a barely-discernible shift for Q’=106 (17 s over 16 years).
Equation 3 assumes a stellar obliquity of zero and neglects tidal dissipation in the planet, assuming its orbit to be circularized and its spin to be synchronized and aligned with the orbit. The canonical value of Q’ is 106, as derived for stars from measurements of the orbits of binary star systems (e.g. Meibom & Mathieu 2005) and for Solar System giant planets from the orbits of their satellites (Zhang & Hamilton, 2008)).
We return to Figure 2, as we can now interpret the findings for (and therefore Q’) physically. The 95th percentile posterior probability distribution for q is effectively zero; given that Q’ , this means we only can definitively extract a lower limit, Q’1, taken at the 5th percentile of the distribution. Continued monitoring of this system should further constrain WASP-18’s Q’, and it follows from the discussion above that we will continue to find an increasing lower limit, i.e. no evidence of rapid tidal decay.
6 Conclusion
We have combined previously published and new data to find no conclusive evidence of rapid tidal decay of the orbit of WASP-18b, supporting predictions of little to no tidal decay for a short-period planet around an F star (Barker & Ogilvie, 2009, 2010; Barker, 2011), given our current understanding of the physics of tidal dissipation in F stars. We find for WASP-18 that Q’ 1at 95 % confidence. Further observations of WASP-18b and similar monitoring of planets like WASP-19b, WTS-2b, and K2-22b would add tighter observational constraints on stellar Q’ for various stellar types, and allow us to further probe the mechanisms of stellar tidal dissipation.
Acknowledgements
The authors thank the referee for a careful eye and helpful insights. TRAPPIST is a project funded by the Belgian Fund for Scientific Research (F.R.S.-FNRS) under grant FRFC 2.5.594.09.F, with the participation of the Swiss National Science Foundation (SNF). L.D. acknowledges support of the F.R.I.A. fund of the F.R.S-FNRS. M. Gillon and E. Jehin are F.R.S.-FNRS Research Associates. AJB was supported by The Leverhulme Trust through the award of an Early Career Fellowship.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1Barker (2011) Barker, A. J. 2011, MNRAS, 414, 1365
- 2Barker & Ogilvie (2009) Barker, A. J., & Ogilvie, G. I. 2009, MNRAS, 395, 2268
- 3Barker & Ogilvie (2010) —. 2010, MNRAS, 404, 1849
- 4Barker & Ogilvie (2011) —. 2011, MNRAS, 417, 745
- 5Berta et al. (2012) Berta, Z. K., Charbonneau, D., Désert, J.-M., et al. 2012, Ap J, 747, 35
- 6Birkby et al. (2014) Birkby, J. L., Cappetta, M., Cruz, P., et al. 2014, MNRAS, 440, 1470
- 7Bonfanti et al. (2016) Bonfanti, A., Ortolani, S., & Nascimbeni, V. 2016, A&A, 585, A 5
- 8Claret (2000) Claret, A. 2000, A&A, 363, 1081
