Showing 1–50 of 70 results for author: Boue, G

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

Abstract: PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observati… ▽ More PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases. △ Less

Submitted 8 June, 2024; originally announced June 2024.

The tidal deformation and atmosphere of WASP-12b from its phase curve

Authors: B. Akinsanmi, S. C. C. Barros, M. Lendl, L. Carone, P. E. Cubillos, A. Bekkelien, A. Fortier, H. -G. Florén, A. Collier Cameron, G. Boué, G. Bruno, B. -O. Demory, A. Brandeker, S. G. Sousa, T. G. Wilson, A. Deline, A. Bonfanti, G. Scandariato, M. J. Hooton, A. C. M. Correia, O. D. S. Demangeon, A. M. S. Smith, V. Singh, Y. Alibert, R. Alonso , et al. (63 additional authors not shown)

Abstract: Ultra-hot Jupiters present a unique opportunity to understand the physics and chemistry of planets at extreme conditions. WASP-12b stands out as an archetype of this class of exoplanets. We performed comprehensive analyses of the transits, occultations, and phase curves of WASP-12b by combining new CHEOPS observations with previous TESS and Spitzer data to measure the planet's tidal deformation, a… ▽ More Ultra-hot Jupiters present a unique opportunity to understand the physics and chemistry of planets at extreme conditions. WASP-12b stands out as an archetype of this class of exoplanets. We performed comprehensive analyses of the transits, occultations, and phase curves of WASP-12b by combining new CHEOPS observations with previous TESS and Spitzer data to measure the planet's tidal deformation, atmospheric properties, and orbital decay rate. The planet was modeled as a triaxial ellipsoid parameterized by the second-order fluid Love number, $h_2$, which quantifies its radial deformation and provides insight into the interior structure. We measured the tidal deformation of WASP-12b and estimated a Love number of $h_2=1.55_{-0.49}^{+0.45}$ (at 3.2$σ$) from its phase curve. We measured occultation depths of $333\pm24$ppm and $493\pm29$ppm in the CHEOPS and TESS bands, respectively, while the dayside emission spectrum indicates that CHEOPS and TESS probe similar pressure levels in the atmosphere at a temperature of 2900K. We also estimated low geometric albedos of $0.086\pm0.017$ and $0.01\pm0.023$ in the CHEOPS and TESS passbands, respectively, suggesting the absence of reflective clouds in the dayside of the WASP-12b. The CHEOPS occultations do not show strong evidence for variability in the dayside atmosphere of the planet. Finally, we refine the orbital decay rate by 12% to a value of -30.23$\pm$0.82 ms/yr. WASP-12b becomes the second exoplanet, after WASP-103b, for which the Love number has been measured (at 3$sigma$) from the effect of tidal deformation in the light curve. However, constraining the core mass fraction of the planet requires measuring $h_2$ with a higher precision. This can be achieved with high signal-to-noise observations with JWST since the phase curve amplitude, and consequently the induced tidal deformation effect, is higher in the infrared. △ Less

Submitted 20 February, 2024; v1 submitted 16 February, 2024; originally announced February 2024.

Comments: accepted for publication in A&A

Can one hear supercontinents in the tides of ocean planets?

Authors: Pierre Auclair-Desrotour, Mohammad Farhat, Gwenaël Boué, Mickaël Gastineau, Jacques Laskar

Abstract: Recent observations and theoretical progress made about the history of the Earth-Moon system suggest that tidal dissipation in oceans primarily drives the long term evolution of orbital systems hosting ocean planets. Particularly, they emphasise the key role played by the geometry of land-ocean distributions in this mechanism. However, the complex way continents affect oceanic tides still remains… ▽ More Recent observations and theoretical progress made about the history of the Earth-Moon system suggest that tidal dissipation in oceans primarily drives the long term evolution of orbital systems hosting ocean planets. Particularly, they emphasise the key role played by the geometry of land-ocean distributions in this mechanism. However, the complex way continents affect oceanic tides still remains to be elucidated. In the present study, we investigate the impact of a single supercontinent on the tidal response of an ocean planet and the induced tidally dissipated energy. The adopted approach is based on the linear tidal theory. By simplifying the continent to a spherical cap of given angular radius and position on the globe, we proceed to a harmonic analysis of the whole planet's tidal response including the coupling with the solid part due to ocean loading and self-attraction variations. In this framework, tidal flows are formulated analytically in terms of explicitly defined oceanic eigenmodes, as well as the resulting tidal Love numbers, dissipated power, and torque. The analysis highlights the symmetry breaking effect of the continent, which makes the dependence of tidal quantities on the tidal frequency become highly irregular. The metric introduced to quantify this continentality effect reveals abrupt transitions between polar and non-polar configurations, and between small-sized and medium-sized continents. Additionally, it predicts that a continent similar to South America or smaller (30°-angular radius) does not alter qualitatively the tidal response of a global ocean whatever its position on the planet. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: 35 pages, 13 figures, 5 tables. Accepted for publication in Astronomy & Astrophysics

MSC Class: 85-02

On the effect of tidal deformation on planetary phase curves

Authors: Babatunde Akinsanmi, Monika Lendl, Gwenael Boue, Susana C. C. Barros

Abstract: With the continuous improvement in the precision of exoplanet observations, it has become feasible to probe for subtle effects that can enable a more comprehensive characterization of exoplanets. A notable example is the tidal deformation of ultra-hot Jupiters by their host stars, whose detection can provide valuable insights into the planetary interior structure. In this work, we extend previous… ▽ More With the continuous improvement in the precision of exoplanet observations, it has become feasible to probe for subtle effects that can enable a more comprehensive characterization of exoplanets. A notable example is the tidal deformation of ultra-hot Jupiters by their host stars, whose detection can provide valuable insights into the planetary interior structure. In this work, we extend previous research on modeling deformation in transit light curves by proposing a simple approach to account for tidal deformation in phase curves. The planetary shape is modeled as a function of the second fluid Love number for radial deformation $h_{2f}$. We show that the effect of tidal deformation manifests across the full orbit of the planet as its projected area varies with phase, thereby allowing us to better probe the planet's shape in phase curves than in transits. Comparing the effects and detectability of deformation by different space-based instruments, we find that the effect of deformation is more prominent in infrared observations where the phase curve amplitude is the largest. A single JWST phase curve observation of a deformed planet, such as WASP-12b, can allow up to 17$σ$ measurement of $h_{2f}$ compared to 4$σ$ from transit-only observation. Such high precision $h_{2f}$ measurement can constrain the core mass of the planet to within 19\% of the total mass, thus providing unprecedented constraints on the interior structure. Due to the lower phase curve amplitudes in the optical, the other instruments provide $\leq4σ$ precision on $h_{2f}$ depending on the number of phase curves observed. We also find that detecting deformation from infrared phase curves is less affected by uncertainty in limb darkening, unlike detection in transits. Finally, the assumption of sphericity when analyzing the phase curve of deformed planets can lead to biases in several system parameters △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 15 pages, 14 figures, Accepted for publication in Astronomy and Astrophysics

Thermal tides in neutrally stratified atmospheres: Revisiting the Earth's Precambrian rotational equilibrium

Authors: Mohammad Farhat, Pierre Auclair-Desrotour, Gwenaël Boué, Russell Deitrick, Jacques Laskar

Abstract: Rotational dynamics of the Earth, over geological timescales, have profoundly affected local and global climatic evolution, probably contributing to the evolution of life. To better retrieve the Earth's rotational history, and motivated by the published hypothesis of a stabilized length of day during the Precambrian, we examine the effect of thermal tides on the evolution of planetary rotational m… ▽ More Rotational dynamics of the Earth, over geological timescales, have profoundly affected local and global climatic evolution, probably contributing to the evolution of life. To better retrieve the Earth's rotational history, and motivated by the published hypothesis of a stabilized length of day during the Precambrian, we examine the effect of thermal tides on the evolution of planetary rotational motion. The hypothesized scenario is contingent upon encountering a resonance in atmospheric Lamb waves, whereby an amplified thermotidal torque cancels the opposing torque of the oceans and solid interior, driving the Earth into a rotational equilibrium. With this scenario in mind, we construct an ab initio model of thermal tides on rocky planets describing a neutrally stratified atmosphere. The model takes into account dissipative processes with Newtonian cooling and diffusive processes in the planetary boundary layer. We retrieve from this model a closed-form solution for the frequency-dependent tidal torque which captures the main spectral features previously computed using 3D general circulation models. In particular, under longwave heating, diffusive processes near the surface and the delayed thermal response of the ground prove to be responsible for attenuating, and possibly annihilating, the accelerating effect of the thermotidal torque at the resonance. When applied to the Earth, our model prediction suggests the occurrence of the Lamb resonance in the Phanerozoic, but with an amplitude that is insufficient for the rotational equilibrium. Interestingly, though our study was motivated by the Earth's history, the generic tidal solution can be straightforwardly and efficiently applied in exoplanetary settings. △ Less

Submitted 22 September, 2023; v1 submitted 21 September, 2023; originally announced September 2023.

Comments: 20 pages (+14 for appendices), 6 figures

Did atmospheric thermal tides cause a daylength locking in the Precambrian? A review on recent results

Authors: Jacques Laskar, Mohammad Farhat, Margriet L. Lantink, Pierre Auclair-Desrotour, Gwenaël Boué, Matthias Sinnesael

Abstract: After the initial suggestion by Zahnle and Walker (1987) that the torque accelerating the spin rate of the Earth and produced by the heating of the atmosphere by the Sun could counteract the braking lunir-solar gravitational torque in the Precambrian, several authors have recently revisited this hypothesis. In these studies, it is argued that the geological evidences of the past spin state of the… ▽ More After the initial suggestion by Zahnle and Walker (1987) that the torque accelerating the spin rate of the Earth and produced by the heating of the atmosphere by the Sun could counteract the braking lunir-solar gravitational torque in the Precambrian, several authors have recently revisited this hypothesis. In these studies, it is argued that the geological evidences of the past spin state of the Earth play in favor of this atmospheric tidal locking of the length of the day (LOD). In the present review of the recent literature, we show that the drawn conclusions depend crucially on the consideration of the stromatolite geological LOD estimates obtained by Pannella at 1.88 and 2.0 Ga, which are subject to large uncertainties. When only the most robust cyclostatigraphic estimates of the LOD are retained, the LOD locking hypothesis is not supported. Moreover, the consideration of the published General Circulation Model numerical simulations and of new analytical models for the thermal atmospheric tides suggest that the atmospheric tidal resonance, which is the crucial ingredient for the LOD locking in the Precambrian, was never of sufficiently large amplitude to allow for this tidal LOD lock. △ Less

Submitted 21 September, 2023; v1 submitted 20 September, 2023; originally announced September 2023.

Comments: 16 pages, 9 figures

TOI-4860 b, a short-period giant planet transiting an M3.5 dwarf

Authors: J. M. Almenara, X. Bonfils, E. M. Bryant, A. Jordán, G. Hébrard, E. Martioli, A. C. M. Correia, N. Astudillo-Defru, C. Cadieux, L. Arnold, É. Artigau, G. Á. Bakos, S. C. C. Barros, D. Bayliss, F. Bouchy, G. Boué, R. Brahm, A. Carmona, D. Charbonneau, D. R. Ciardi, R. Cloutier, M. Cointepas, N. J. Cook, N. B. Cowan, X. Delfosse , et al. (25 additional authors not shown)

Abstract: We report the discovery and characterisation of a giant transiting planet orbiting a nearby M3.5V dwarf (d = 80.4 pc, $G$ = 15.1 mag, $K$=11.2 mag, R$_\star$ = 0.358 $\pm$ 0.015 R$_\odot$, M$_\star$ = 0.340 $\pm$ 0.009 M$_\odot$). Using the photometric time series from TESS sectors 10, 36, 46, and 63 and near-infrared spectrophotometry from ExTrA, we measured a planetary radius of 0.77 $\pm$ 0.03… ▽ More We report the discovery and characterisation of a giant transiting planet orbiting a nearby M3.5V dwarf (d = 80.4 pc, $G$ = 15.1 mag, $K$=11.2 mag, R$_\star$ = 0.358 $\pm$ 0.015 R$_\odot$, M$_\star$ = 0.340 $\pm$ 0.009 M$_\odot$). Using the photometric time series from TESS sectors 10, 36, 46, and 63 and near-infrared spectrophotometry from ExTrA, we measured a planetary radius of 0.77 $\pm$ 0.03 R$_J$ and an orbital period of 1.52 days. With high-resolution spectroscopy taken by the CFHT/SPIRou and ESO/ESPRESSO spectrographs, we refined the host star parameters ([Fe/H] = 0.27 $\pm$ 0.12) and measured the mass of the planet (0.273 $\pm$ 0.006 M$_J$). Based on these measurements, TOI-4860 b joins the small set of massive planets ($>$80 M$_E$) found around mid to late M dwarfs ($<$0.4 R$_\odot$), providing both an interesting challenge to planet formation theory and a favourable target for further atmospheric studies with transmission spectroscopy. We identified an additional signal in the radial velocity data that we attribute to an eccentric planet candidate ($e=0.66\pm0.09$) with an orbital period of $427\pm7$~days and a minimum mass of $1.66\pm 0.26$ M$_J$, but additional data would be needed to confirm this. △ Less

Submitted 12 January, 2024; v1 submitted 2 August, 2023; originally announced August 2023.

Comments: 16 pages, 14 figures, accepted for publication in A&A

The magnetic field and multiple planets of the young dwarf AU~Mic

Authors: J. -F. Donati, P. I. Cristofari, B. Finociety, B. Klein, C. Moutou, E. Gaidos, C. Cadieux, E. Artigau, A. C. M. Correia, G. Boué, N. J. Cook, A. Carmona, L. T. Lehmann, J. Bouvier, E. Martioli, J. Morin, P. Fouqué, X. Delfosse, R. Royon, G. Hébrard, S. H. P. Alencar, J. Laskar, L. Arnold, P. Petit, A. Kospal , et al. (3 additional authors not shown)

Abstract: In this paper we present an analysis of near-infrared spectropolarimetric and velocimetric data of the young M dwarf AU Mic, collected with SPIRou at the Canada-France-Hawaii telescope from 2019 to 2022, mostly within the SPIRou Legacy Survey. With these data, we study the large- and small-scale magnetic field of AU Mic, detected through the unpolarized and circularly-polarized Zeeman signatures o… ▽ More In this paper we present an analysis of near-infrared spectropolarimetric and velocimetric data of the young M dwarf AU Mic, collected with SPIRou at the Canada-France-Hawaii telescope from 2019 to 2022, mostly within the SPIRou Legacy Survey. With these data, we study the large- and small-scale magnetic field of AU Mic, detected through the unpolarized and circularly-polarized Zeeman signatures of spectral lines. We find that both are modulated with the stellar rotation period (4.86 d), and evolve on a timescale of months under differential rotation and intrinsic variability. The small-scale field, estimated from the broadening of spectral lines, reaches $2.61\pm0.05$ kG. The large-scale field, inferred with Zeeman-Doppler imaging from Least-Squares Deconvolved profiles of circularly-polarized and unpolarized spectral lines, is mostly poloidal and axisymmetric, with an average intensity of $550\pm30$ G. We also find that surface differential rotation, as derived from the large-scale field, is $\simeq$30% weaker than that of the Sun. We detect the radial velocity (RV) signatures of transiting planets b and c, although dwarfed by activity, and put an upper limit on that of candidate planet d, putatively causing the transit-timing variations of b and c. We also report the detection of the RV signature of a new candidate planet (e) orbiting further out with a period of $33.39\pm0.10$ d, i.e., near the 4:1 resonance with b. The RV signature of e is detected at 6.5$σ$ while those of b and c show up at $\simeq$4$σ$, yielding masses of $10.2^{+3.9}_{-2.7}$ and $14.2^{+4.8}_{-3.5}$ Earth masses for b and c, and a minimum mass of $35.2^{+6.7}_{-5.4}$ Earth masses for e. △ Less

Submitted 24 April, 2023; v1 submitted 19 April, 2023; originally announced April 2023.

Comments: MNRAS, in press (20 pages and 12 figures + 9 pages of supplementary material)

Examining the orbital decay targets KELT-9 b, KELT-16 b and WASP-4 b, and the transit-timing variations of HD 97658 b

Authors: J. -V. Harre, A. M. S. Smith, S. C. C. Barros, G. Boué, Sz. Csizmadia, D. Ehrenreich, H. -G. Florén, A. Fortier, P. F. L. Maxted, M. J. Hooton, B. Akinsanmi, L. M. Serrano, N. M. Rosário, B. -O. Demory, K. Jones, J. Laskar, V. Adibekyan, Y. Alibert, R. Alonso, D. R. Anderson, G. Anglada, J. Asquier, T. Bárczy, D. Barrado y Navascues, W. Baumjohann , et al. (65 additional authors not shown)

Abstract: Tidal orbital decay is suspected to occur especially for hot Jupiters, with the only observationally confirmed case of this being WASP-12 b. By examining this effect, information on the properties of the host star can be obtained using the so-called stellar modified tidal quality factor $Q_*'$, which describes the efficiency with which kinetic energy of the planet is dissipated within the star. Th… ▽ More Tidal orbital decay is suspected to occur especially for hot Jupiters, with the only observationally confirmed case of this being WASP-12 b. By examining this effect, information on the properties of the host star can be obtained using the so-called stellar modified tidal quality factor $Q_*'$, which describes the efficiency with which kinetic energy of the planet is dissipated within the star. This can help to get information about the interior of the star. In this study, we aim to improve constraints on the tidal decay of the KELT-9, KELT-16 and WASP-4 systems, to find evidence for or against the presence of this particular effect. With this, we want to constrain each star's respective $Q_*'$ value. In addition to that, we also aim to test the existence of the transit timing variations (TTVs) in the HD 97658 system, which previously favoured a quadratic trend with increasing orbital period. Making use of newly acquired photometric observations from CHEOPS and TESS, combined with archival transit and occultation data, we use Markov chain Monte Carlo (MCMC) algorithms to fit three models, a constant period model, an orbital decay model, and an apsidal precession model, to the data. We find that the KELT-9 system is best described by an apsidal precession model for now, with an orbital decay trend at over 2 $σ$ being a possible solution as well. A Keplerian orbit model with a constant orbital period fits the transit timings of KELT-16 b the best due to the scatter and scale of their error bars. The WASP-4 system is represented the best by an orbital decay model at a 5 $σ$ significance, although apsidal precession cannot be ruled out with the present data. For HD 97658 b, using recently acquired transit observations, we find no conclusive evidence for a previously suspected strong quadratic trend in the data. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Comments: 19 pages, 11 figures, 9 tables, accepted for publication in A&A

Journal ref: A&A 669, A124 (2023)

Tilting Uranus via the migration of an ancient satellite

Authors: Melaine Saillenfest, Zeeve Rogoszinski, Giacomo Lari, Kevin Baillié, Gwenaël Boué, Aurélien Crida, Valéry Lainey

Abstract: Context. The 98°-obliquity of Uranus is commonly attributed to giant impacts that occurred at the end of the planetary formation. This picture, however, is not devoid of weaknesses. Aims. On a billion-year timescale, the tidal migration of the satellites of Jupiter and Saturn has been shown to strongly affect their spin-axis dynamics. We aim to revisit the scenario of tilting Uranus in light of… ▽ More Context. The 98°-obliquity of Uranus is commonly attributed to giant impacts that occurred at the end of the planetary formation. This picture, however, is not devoid of weaknesses. Aims. On a billion-year timescale, the tidal migration of the satellites of Jupiter and Saturn has been shown to strongly affect their spin-axis dynamics. We aim to revisit the scenario of tilting Uranus in light of this mechanism. Methods. We analyse the precession spectrum of Uranus and identify the candidate secular spin-orbit resonances that could be responsible for the tilting. We determine the properties of the hypothetical ancient satellite required for a capture and explore the dynamics numerically. Results. If it migrates over 10 Uranus' radii, a single satellite with minimum mass 4e-4 Uranus' mass is able to tilt Uranus from a small obliquity and make it converge towards 90°. In order to achieve the tilting in less than the age of the Solar System, the mean drift rate of the satellite must be comparable to the Moon's current orbital expansion. Under these conditions, simulations show that Uranus is readily tilted over 80°. Beyond this point, the satellite is strongly destabilised and triggers a phase of chaotic motion for the planet's spin axis. The chaotic phase ends when the satellite collides into the planet, ultimately freezing the planet's obliquity in either a prograde, or plainly retrograde state (as Uranus today). Spin states resembling that of Uranus can be obtained with probabilities as large as 80%, but a bigger satellite is favoured, with mass 1.7e-3 Uranus' mass or more. Yet, a smaller ancient satellite is not categorically ruled out, and there is room for improving this basic scenario in future studies. Interactions among several pre-existing satellites is a promising possibility. △ Less

Submitted 21 September, 2022; originally announced September 2022.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 668, A108 (2022)

The Resonant Tidal Evolution of the Earth-Moon Distance

Authors: Mohammad Farhat, Pierre Auclair-Desrotour, Gwenaël Boué, Jacques Laskar

Abstract: Due to tidal interactions in the Earth-Moon system, the spin of the Earth slows down and the Moon drifts away. This recession of the Moon is now measured with great precision, but it has been realized, more than fifty years ago, that simple tidal models extrapolated back in time lead to an age of the Moon that is by far incompatible with the geochronological and geochemical evidence. In order to e… ▽ More Due to tidal interactions in the Earth-Moon system, the spin of the Earth slows down and the Moon drifts away. This recession of the Moon is now measured with great precision, but it has been realized, more than fifty years ago, that simple tidal models extrapolated back in time lead to an age of the Moon that is by far incompatible with the geochronological and geochemical evidence. In order to evade this problem, more elaborate models have been proposed, taking into account the oceanic tidal dissipation. However, these models did not fit both the estimated lunar age and the present rate of lunar recession simultaneously. Here we present a physical model that reconciles these two constraints and yields a unique solution of the tidal history. This solution fits well the available geological proxies for the history of the Earth-Moon system and consolidates the cyclostratigraphic method. The resulting evolution involves multiple crossings of resonances in the oceanic dissipation that are associated with significant and rapid variations in the lunar orbital distance, the Earth's length of the day, and the Earth's obliquity. △ Less

Submitted 1 July, 2022; originally announced July 2022.

Journal ref: A&A 665, L1 (2022)

Detection of the tidal deformation of WASP-103b at $3\,σ$ with CHEOPS

Authors: S. C. C. Barros, B. Akinsanm, G. Boué, A. M. S. Smith, J. Laskar, S. Ulmer-Moll, J. Lillo-Box, the CHEOPS team

Abstract: WASP-103b is the exoplanet with the highest expected deformation signature in its transit light curve and one of the shortest expected spiral-in times. Measuring the tidal deformation of the planet would allow us to estimate the second degree fluid Love number and gain insight into the planet's internal structure. Moreover, measuring the tidal decay timescale would allow us to estimate the stellar… ▽ More WASP-103b is the exoplanet with the highest expected deformation signature in its transit light curve and one of the shortest expected spiral-in times. Measuring the tidal deformation of the planet would allow us to estimate the second degree fluid Love number and gain insight into the planet's internal structure. Moreover, measuring the tidal decay timescale would allow us to estimate the stellar tidal quality factor, which is key to constraining stellar physics. We obtained 12 transit light curves of WASP-103b with the CHEOPS to estimate the tidal deformation and tidal decay of this extreme system. We modelled the high-precision CHEOPS transit light curves together with systematic instrumental noise using multi-dimensional Gaussian process regression informed by a set of instrumental parameters. To model the tidal deformation, we used a parametrisation model which allowed us to determine the second degree fluid Love number of the planet. We combined our light curves with previously observed transits of WASP-103b with HST and Spitzer. We estimate the radial Love number of WASP-103b to be $h_f = 1.59^{+0.45}_{-0.53}$. This is the first time that the tidal deformation is directly detected (at $3\, σ$) from the transit light curve of an exoplanet. Combining the transit times derived from CHEOPS, with the other transit times available in the literature, we find no significant orbital period variation for WASP-103b. However, the data show a hint of an orbital period increase instead of a decrease, as is expected for tidal decay. This could be either due to a visual companion star if this star is bound, the Applegate effect, or a statistical artefact. The estimated Love number of WASP-103b is similar to Jupiter's. This will allow us to constrain the internal structure and composition of WASP-103b, which could provide clues on the inflation of hot Jupiters. (Abridged) △ Less

Submitted 10 January, 2022; originally announced January 2022.

Comments: 21 pages, 9 figures, 11 tables, CHEOPS light curves available online, Published A&A

Ellipsoidal equilibrium figure and Cassini states of rotating planets and satellites deformed by a tidal potential in the spatial case

Authors: Hugo A. Folonier, Gwenaël Boué, Sylvio Ferraz-Mello

Abstract: The equilibrium figure of an inviscid tidally deformed body is the starting point for the construction of many tidal theories such as Darwinian tidal theories or the hydrodynamical Creep tide theory. This paper presents the ellipsoidal equilibrium figure when the spin rate vector of the deformed body is not perpendicular to the plane of motion of the companion. We obtain the equatorial and the pol… ▽ More The equilibrium figure of an inviscid tidally deformed body is the starting point for the construction of many tidal theories such as Darwinian tidal theories or the hydrodynamical Creep tide theory. This paper presents the ellipsoidal equilibrium figure when the spin rate vector of the deformed body is not perpendicular to the plane of motion of the companion. We obtain the equatorial and the polar flattenings as functions of the Jeans and the Maclaurin flattenings, and of the angle $θ$ between the spin rate vector and the radius vector. The equatorial vertex of the equilibrium ellipsoid does not point toward the companion, which produces a torque perpendicular to the rotation vector, which introduces terms of precession and nutation. We find that the direction of spin may differ significantly from the direction of the principal axis of inertia $C$, so the classical approximation $\mathsf{I}\vecω\approx C\vecω$ only makes sense in the neighborhood of the planar problem. We also study the so-called Cassini states. Neglecting the short-period terms in the differential equation for the spin direction and assuming a uniform precession of the line of the orbital ascending node, we obtain the same differential equation as that found by Colombo (1966). That is, a tidally deformed inviscid body has exactly the same Cassini states as a rotating axisymmetric rigid body, the tidal bulge having no secular effect at first order. △ Less

Submitted 27 December, 2021; originally announced December 2021.

Comments: 17 pages, 7 figures

The large obliquity of Saturn explained by the fast migration of Titan

Authors: Melaine Saillenfest, Giacomo Lari, Gwenaël Boué

Abstract: The obliquity of a planet is the tilt between its equator and its orbital plane. Giant planets are expected to form with near-zero obliquities. After its formation, some dynamical mechanism must therefore have tilted Saturn up to its current obliquity of 26.7°. This event is traditionally thought to have happened more than 4 Gyrs ago during the late planetary migration because of the crossing of a… ▽ More The obliquity of a planet is the tilt between its equator and its orbital plane. Giant planets are expected to form with near-zero obliquities. After its formation, some dynamical mechanism must therefore have tilted Saturn up to its current obliquity of 26.7°. This event is traditionally thought to have happened more than 4 Gyrs ago during the late planetary migration because of the crossing of a resonance between the spin-axis precession of Saturn and the nodal orbital precession mode of Neptune. Here, we show that the fast tidal migration of Titan measured by Lainey et al. (2020) is incompatible with this scenario, and that it offers a new explanation for Saturn's current obliquity. A significant migration of Titan would prevent any early resonance, invalidating previous constraints on the late planetary migration set by the tilting of Saturn. We propose instead that the resonance was encountered recently, about 1 Gyr ago, forcing Saturn's obliquity to increase from a small value (possibly less than 3°), up to its current state. This scenario suggests that Saturn's normalised polar moment of inertia lies between 0.224 and 0.237. Our findings bring out a new paradigm for the spin-axis evolution of Saturn, Jupiter, and possibly giant exoplanets in multiple systems, whereby obliquities are not settled once for all, but continuously evolve as a result of the migration of their satellites. △ Less

Submitted 8 October, 2021; originally announced October 2021.

Comments: Author's manuscript of the letter published in Nature Astronomy on 18-01-2021. The published version can be found at https://www.nature.com/articles/s41550-020-01284-x

Journal ref: Nature Astronomy 5, 345-349 (2021)

Librations of a body composed of a deformable mantle and a fluid core

Authors: Clodoaldo Ragazzo, Gwenaël Boué, Yeva Gevorgyan, Lucas S. Ruiz

Abstract: We present fully three-dimensional equations to describe the rotations of a body made of a deformable mantle and a fluid core. The model in its essence is similar to that used by INPOP (Integration Planétaire de l'Observatoire de Paris), e.g. Viswanathan et al. (2019), and by JPL (Jet Propulsion Laboratory), e.g. Folkner et al. (2014), to represent the Moon. The intended advantages of our model ar… ▽ More We present fully three-dimensional equations to describe the rotations of a body made of a deformable mantle and a fluid core. The model in its essence is similar to that used by INPOP (Integration Planétaire de l'Observatoire de Paris), e.g. Viswanathan et al. (2019), and by JPL (Jet Propulsion Laboratory), e.g. Folkner et al. (2014), to represent the Moon. The intended advantages of our model are: straightforward use of any linear-viscoelastic model for the rheology of the mantle; easy numerical implementation in time-domain (no time lags are necessary); all parameters, including those related to the "permanent deformation", have a physical interpretation. The paper also contains: 1) A physical model to explain the usual lack of hydrostaticity of the mantle (permanent deformation). 2) Formulas for free librations of bodies in and out-of spin-orbit resonance that are valid for any linear viscoelastic rheology of the mantle. 3) Formulas for the offset between the mantle and the idealized rigid-body motion (Peale's Cassini states). 4) Applications to the librations of Moon, Earth, and Mercury that are used for model validation. △ Less

Submitted 22 March, 2022; v1 submitted 17 August, 2021; originally announced August 2021.

Comments: 92 pages, 15 figures

Journal ref: Celest Mech Dyn Astr 134, 10 (2022)

Constraining the Earth's Dynamical Ellipticity from Ice Age Dynamics

Authors: Mohammad Farhat, Jacques Laskar, Gwenaël Boué

Abstract: The dynamical ellipticity of a planet expresses the departure of its mass distribution from spherical symmetry. It enters as a parameter in the description of a planet's precession and nutation, as well as other rotational normal modes. In the case of the Earth, uncertainties in this quantity's history produce an uncertainty in the solutions for the past evolution of the Earth-Moon system. Constra… ▽ More The dynamical ellipticity of a planet expresses the departure of its mass distribution from spherical symmetry. It enters as a parameter in the description of a planet's precession and nutation, as well as other rotational normal modes. In the case of the Earth, uncertainties in this quantity's history produce an uncertainty in the solutions for the past evolution of the Earth-Moon system. Constraining this history has been a target of interdisciplinary efforts as it represents an astro-geodetic parameter whose variation is shaped by geophysical processes, and whose imprints can be found in the geological signal. We revisit the classical problem of its variation during ice ages, where glacial cycles exerted a varying surface loading that had altered the shape of the geoid. In the framework of glacial isostatic adjustment, and with the help of a recent paleoclimatic proxy of ice volume, we present the evolution of the dynamical ellipticity over the Cenozoic ice ages. We map out the problem in full generality identifying major sensitivities to surface loading and internal variations in parameter space. This constrained evolution is aimed to be used in future astronomical computations of the orbital and insolation quantities of the Earth. △ Less

Submitted 26 March, 2021; originally announced March 2021.

The past and future obliquity of Saturn as Titan migrates

Authors: Melaine Saillenfest, Giacomo Lari, Gwenaël Boué, Ariane Courtot

Abstract: Aims: It has recently been shown that the fast migration of Titan could be responsible for the current 26.7°-tilt of Saturn's spin axis. We aim to quantify the level of generality of this result and to measure the most likely sets of parameters. We also aim to determine the obliquity that Saturn will reach in the future. Methods: We explore a broad range of parameters and propagate numerically t… ▽ More Aims: It has recently been shown that the fast migration of Titan could be responsible for the current 26.7°-tilt of Saturn's spin axis. We aim to quantify the level of generality of this result and to measure the most likely sets of parameters. We also aim to determine the obliquity that Saturn will reach in the future. Methods: We explore a broad range of parameters and propagate numerically the orientation of Saturn's spin axis both backward and forward in time. Results: In the adiabatic regime, the likelihood of reproducing Saturn's current spin-axis orientation is maximised for primordial obliquities between about 2° and 7°. For a slightly faster migration than expected from radio-science experiments, non-adiabatic effects even allow for exactly null primordial obliquities. Starting from such small tilts, Saturn's spin axis can evolve up to its current state provided that: i) the semi-major axis of Titan changed by more than 5% of its current value, and ii) its migration rate does not exceed 10 times the nominal measured rate. In comparison, observational data suggest that the increase of Titan's semi-major axis exceeded 50% over 4 Gyrs, and error bars imply that the current migration rate is unlikely to be larger than 1.5 times its nominal value. Conclusions: If Titan did migrate substantially before today, tilting Saturn from a small obliquity is not only possible, but it is the most likely scenario. Saturn's obliquity is expected to be still increasing today and could exceed 65° in the future. Maximising the likelihood would also put strict constraints on Saturn's polar moment of inertia. However, the possibility remains that Saturn's primordial obliquity was already large, for instance as a result of a massive collision. The unambiguous distinction between these two scenarios would be given by a precise measure of Saturn's polar moment of inertia. △ Less

Submitted 8 March, 2021; v1 submitted 25 January, 2021; originally announced January 2021.

Comments: Published in Astronomy and Astrophysics. Some typographical errors have been corrected

Journal ref: A&A 647, A92 (2021)

Six transiting planets and a chain of Laplace resonances in TOI-178

Authors: A. Leleu, Y. Alibert, N. C. Hara, M. J. Hooton, T. G. Wilson, P. Robutel, J. -B. Delisle, J. Laskar, S. Hoyer, C. Lovis, E. M. Bryant, E. Ducrot, J. Cabrera, L. Delrez, J. S. Acton, V. Adibekyan, R. Allart, C. Allende Prieto, R. Alonso, D. Alves, D. R. Anderson, D. Angerhausen, G. Anglada Escudé, J. Asquier, D. Barrado , et al. (130 additional authors not shown)

Abstract: Determining the architecture of multi-planetary systems is one of the cornerstones of understanding planet formation and evolution. Resonant systems are especially important as the fragility of their orbital configuration ensures that no significant scattering or collisional event has taken place since the earliest formation phase when the parent protoplanetary disc was still present. In this cont… ▽ More Determining the architecture of multi-planetary systems is one of the cornerstones of understanding planet formation and evolution. Resonant systems are especially important as the fragility of their orbital configuration ensures that no significant scattering or collisional event has taken place since the earliest formation phase when the parent protoplanetary disc was still present. In this context, TOI-178 has been the subject of particular attention since the first TESS observations hinted at a 2:3:3 resonant chain. Here we report the results of observations from CHEOPS, ESPRESSO, NGTS, and SPECULOOS with the aim of deciphering the peculiar orbital architecture of the system. We show that TOI-178 harbours at least six planets in the super-Earth to mini-Neptune regimes, with radii ranging from 1.152(-0.070/+0.073) to 2.87(-0.13/+0.14) Earth radii and periods of 1.91, 3.24, 6.56, 9.96, 15.23, and 20.71 days. All planets but the innermost one form a 2:4:6:9:12 chain of Laplace resonances, and the planetary densities show important variations from planet to planet, jumping from 1.02(+0.28/-0.23) to 0.177(+0.055/-0.061) times the Earth's density between planets c and d. Using Bayesian interior structure retrieval models, we show that the amount of gas in the planets does not vary in a monotonous way, contrary to what one would expect from simple formation and evolution models and unlike other known systems in a chain of Laplace resonances. The brightness of TOI-178 allows for a precise characterisation of its orbital architecture as well as of the physical nature of the six presently known transiting planets it harbours. The peculiar orbital configuration and the diversity in average density among the planets in the system will enable the study of interior planetary structures and atmospheric evolution, providing important clues on the formation of super-Earths and mini-Neptunes. △ Less

Submitted 22 January, 2021; originally announced January 2021.

Discovery and characterization of the exoplanets WASP-148b and c. A transiting system with two interacting giant planets

Authors: G. Hebrard, R. F. Diaz, A. C. M. Correia, A. Collier Cameron, J. Laskar, D. Pollacco, J. -M. Almenara, D. R. Anderson, S. C. C. Barros, I. Boisse, A. S. Bonomo, F. Bouchy, G. Boue, P. Boumis, D. J. A. Brown, S. Dalal, M. Deleuil, O. Demangeon, A. P. Doyle, C. A. Haswell, C. Hellier, H. Osborn, F. Kiefer, U. C. Kolb, K. Lam , et al. (17 additional authors not shown)

Abstract: We present the discovery and characterization of WASP-148, a new extrasolar system that includes at least two giant planets. The host star is a slowly rotating inactive late-G dwarf with a V=12 magnitude. The planet WASP-148b is a hot Jupiter of 0.72 R_Jup and 0.29 M_Jup that transits its host with an orbital period of 8.80 days. We found the planetary candidate with the SuperWASP photometric surv… ▽ More We present the discovery and characterization of WASP-148, a new extrasolar system that includes at least two giant planets. The host star is a slowly rotating inactive late-G dwarf with a V=12 magnitude. The planet WASP-148b is a hot Jupiter of 0.72 R_Jup and 0.29 M_Jup that transits its host with an orbital period of 8.80 days. We found the planetary candidate with the SuperWASP photometric survey, then characterized it with the SOPHIE spectrograph. Our radial velocity measurements subsequently revealed a second planet in the system, WASP-148c, with an orbital period of 34.5 days and a minimum mass of 0.40 M_Jup. No transits of this outer planet were detected. The orbits of both planets are eccentric and fall near the 4:1 mean-motion resonances. This configuration is stable on long timescales, but induces dynamical interactions so that the orbits differ slightly from purely Keplerian orbits. In particular, WASP-148b shows transit-timing variations of typically 15 minutes, making it the first interacting system with transit-timing variations that is detected on ground-based light curves. We establish that the mutual inclination of the orbital plane of the two planets cannot be higher than 35 degrees, and the true mass of WASP-148c is below 0.60 M_Jup. We present photometric and spectroscopic observations of this system that cover a time span of ten years. We also provide their Keplerian and Newtonian analyses; these analyses should be significantly improved through future TESS~observations. △ Less

Submitted 24 June, 2020; v1 submitted 30 April, 2020; originally announced April 2020.

Comments: 16 pages, 10 figures, 7 tables, A&A in press

Journal ref: A&A 640, A32 (2020)

Cassini states of a rigid body with a liquid core

Authors: Gwenaël Boué

Abstract: The purpose of this work is to determine the location and stability of the Cassini states of a celestial body with an inviscid fluid core surrounded by a perfectly rigid mantle. Both situations where the rotation speed is either non-resonant or trapped in a p:1 spin-orbit resonance where p is a half integer are addressed. The rotation dynamics is described by the Poincaré-Hough model which assumes… ▽ More The purpose of this work is to determine the location and stability of the Cassini states of a celestial body with an inviscid fluid core surrounded by a perfectly rigid mantle. Both situations where the rotation speed is either non-resonant or trapped in a p:1 spin-orbit resonance where p is a half integer are addressed. The rotation dynamics is described by the Poincaré-Hough model which assumes a simple motion of the core. The problem is written in a non-canonical Hamiltonian formalism. The secular evolution is obtained without any truncation in obliquity, eccentricity nor inclination. The condition for the body to be in a Cassini state is written as a set of two equations whose unknowns are the mantle obliquity and the tilt angle of the core spin-axis. Solving the system with Mercury's physical and orbital parameters leads to a maximum of 16 different equilibrium configurations, half of them being spectrally stable. In most of these solutions the core is highly tilted with respect to the mantle. The model is also applied to Io and the Moon. △ Less

Submitted 1 April, 2020; originally announced April 2020.

Comments: 25 pages, 3 figures. To be published in Celestial Mechanics and Dynamical Astronomy

Andrade rheology in time-domain. Application to Enceladus' dissipation of energy due to forced libration

Authors: Yeva Gevorgyan, Gwenaël Boué, Clodoaldo Ragazzo, Lucas S. Ruiz, Alexandre C. M. Correia

Abstract: The main purpose of this work is to present a time-domain implementation of the Andrade rheology, instead of the traditional expansion in terms of a Fourier series of the tidal potential. This approach can be used in any fully three dimensional numerical simulation of the dynamics of a system of many deformable bodies. In particular, it allows large eccentricities, large mutual inclinations, and i… ▽ More The main purpose of this work is to present a time-domain implementation of the Andrade rheology, instead of the traditional expansion in terms of a Fourier series of the tidal potential. This approach can be used in any fully three dimensional numerical simulation of the dynamics of a system of many deformable bodies. In particular, it allows large eccentricities, large mutual inclinations, and it is not limited to quasi-periodic perturbations. It can take into account an extended class of perturbations, such as chaotic motions, transient events, and resonant librations. The results are presented by means of a concrete application: the analysis of the libration of Enceladus. This is done by means of both analytic formulas in the frequency domain and direct numerical simulations. We do not a priori assume that Enceladus has a triaxial shape, the eventual triaxiality is a consequence of the satellite motion and its rheology. As a result we obtain an analytic formula for the amplitude of libration that incorporates a new correction due to the rheology. Our results provide an estimation of the amplitude of libration of the core of Enceladus as 0.6% of that of the shell. They also reproduce the observed 10 GW of tidal heat generated by Enceladus with a value of $0.17\times 10^{14}$Pa$\cdot$s for the global effective viscosity under both Maxwell and Andrade rheology. △ Less

Submitted 19 December, 2019; originally announced December 2019.

Comments: Preprint 'elsart' style, 22 pages, 9 multiple figures. Accepted for publication in Icarus

Bias and robustness of eccentricity estimates from radial velocity data

Authors: N. C. Hara, G. Boué, J. Laskar, J. B. Delisle, N. Unger

Abstract: Eccentricity is a parameter of particular interest as it is an informative indicator of the past of planetary systems. It is however not always clear whether the eccentricity fitted on radial velocity data is real or if it is an artefact of an inappropriate modelling. In this work, we address this question in two steps: we first assume that the model used for inference is correct and present inter… ▽ More Eccentricity is a parameter of particular interest as it is an informative indicator of the past of planetary systems. It is however not always clear whether the eccentricity fitted on radial velocity data is real or if it is an artefact of an inappropriate modelling. In this work, we address this question in two steps: we first assume that the model used for inference is correct and present interesting features of classical estimators. Secondly, we study whether the eccentricity estimates are to be trusted when the data contain incorrectly modelled signals, such as missed planetary companions, non Gaussian noises, correlated noises with unknown covariance, etc. Our main conclusion is that data analysis via posterior distributions, with a model including a free error term gives reliable results provided two conditions. First, convergence of the numerical methods needs to be ascertained. Secondly, the noise power spectrum should not have a particularly strong peak at the semi period of the planet of interest. As a consequence, it is difficult to determine if the signal of an apparently eccentric planet might be due to another inner companion in 2:1 mean motion resonance. We study the use of Bayes factors to disentangle these cases. Finally, we suggest methods to check if there are hints of an incorrect model in the residuals. We show on simulated data the performance of our methods and comment on the eccentricities of Proxima b and 55 Cnc f. △ Less

Submitted 3 July, 2019; originally announced July 2019.

Comments: Accepted for publication in MNRAS

High order regularised symplectic integrator for collisional planetary systems

Authors: Antoine C. Petit, Jacques Laskar, Gwenaël Boué, Mickaël Gastineau

Abstract: We present a new mixed variable symplectic (MVS) integrator for planetary systems, that fully resolve close encounters. The method is based on a time regularisation that allows keeping the stability properties of the symplectic integrators, while also reducing the effective step size whenever two planets encounter. We use a high order MVS scheme such that it is possible to integrate with large tim… ▽ More We present a new mixed variable symplectic (MVS) integrator for planetary systems, that fully resolve close encounters. The method is based on a time regularisation that allows keeping the stability properties of the symplectic integrators, while also reducing the effective step size whenever two planets encounter. We use a high order MVS scheme such that it is possible to integrate with large time steps far away from close encounters. We show that this algorithm is able to resolve almost exact collisions (i.e with a mutual separation of a fraction of the physical radius) while using the same time-step as in weakly perturbed problem such as the Solar System. We demonstrate the long term behaviour on systems of six super-Earths experiencing strong scattering for 50 kyr. We compare our algorithm to hybrid methods such as MERCURY and show that for an equivalent cost we obtain much better energy conservation. △ Less

Submitted 1 July, 2019; v1 submitted 30 May, 2019; originally announced May 2019.

Comments: Accepted in Astronomy and Astrophysics, updated to match the published version

Journal ref: A&A 628, A32 (2019)

Orbital relaxation and excitation of planets tidally interacting with white dwarfs

Authors: Dimitri Veras, Michael Efroimsky, Valeri V. Makarov, Gwenaël Boué, Vera Wolthoff, Sabine Reffert, Andreas Quirrenbach, Pier-Emmanuel Tremblay, Boris T. Gänsicke

Abstract: Observational evidence of white dwarf planetary systems is dominated by the remains of exo-asteroids through accreted metals, debris discs, and orbiting planetesimals. However, exo-planets in these systems play crucial roles as perturbing agents, and can themselves be perturbed close to the white dwarf Roche radius. Here, we illustrate a procedure for computing the tidal interaction between a whit… ▽ More Observational evidence of white dwarf planetary systems is dominated by the remains of exo-asteroids through accreted metals, debris discs, and orbiting planetesimals. However, exo-planets in these systems play crucial roles as perturbing agents, and can themselves be perturbed close to the white dwarf Roche radius. Here, we illustrate a procedure for computing the tidal interaction between a white dwarf and a near-spherical solid planet. This method determines the planet's inward and/or outward drift, and whether the planet will reach the Roche radius and be destroyed. We avoid constant tidal lag formulations and instead employ the self-consistent secular Darwin-Kaula expansions from Boué & Efroimsky (2019), which feature an arbitrary frequency dependence on the quality functions. We adopt wide ranges of dynamic viscosities and spin rates for the planet in order to straddle many possible outcomes, and provide a foundation for the future study of individual systems with known or assumed rheologies. We find that: (i) massive Super-Earths are destroyed more readily than minor planets (such as the ones orbiting WD 1145+017 and SDSS J1228+1040), (ii) low-viscosity planets are destroyed more easily than high-viscosity planets, and (iii) the boundary between survival and destruction is likely to be fractal and chaotic. △ Less

Submitted 5 April, 2019; originally announced April 2019.

Comments: Accepted for publication in MNRAS

Tidal evolution of the Keplerian elements

Authors: Gwenaël Boué, Michael Efroimsky

Abstract: We address the expressions for the rates of the Keplerian orbital elements within a two-body problem perturbed by the tides in both partners. The formulae for these rates have appeared in the literature in various forms, at times with errors. We reconsider, from scratch, the derivation of these rates and arrive at the Lagrange-type equations which, in some details, differ from the corresponding eq… ▽ More We address the expressions for the rates of the Keplerian orbital elements within a two-body problem perturbed by the tides in both partners. The formulae for these rates have appeared in the literature in various forms, at times with errors. We reconsider, from scratch, the derivation of these rates and arrive at the Lagrange-type equations which, in some details, differ from the corresponding equations obtained previously by Kaula (1964). We also write down detailed expressions for $da/dt$, $de/dt$ and $di/dt$, to order $e^4$. They differ from Kaula's expressions which contain a redundant factor of $M/(M+M^{\prime}),$ with $M$ and $M^{\prime}$ being the masses of the primary and the secondary. As Kaula was interested in the Earth-Moon system, this redundant factor was close to unity and was unimportant in his developments. This factor, however, must be reinstated when Kaula's theory is applied to a binary composed of partners of comparable masses. We have found that, while it is legitimate to simply sum the primary's and secondary's inputs in $da/dt$ or $de/dt$, this is not the case for $di/dt$. So our expression for $di/dt$ differs from that of Kaula in two regards. First, the contribution due to the dissipation in the secondary averages out when the apsidal precession is uniform. Second, we have obtained an additional term which emerges owing to the conservation of the angular momentum: a change in the inclination of the orbit causes a change of the primary's plane of equator. △ Less

Submitted 26 June, 2019; v1 submitted 3 April, 2019; originally announced April 2019.

Journal ref: Celestial Mechanics and Dynamical Astronomy 131:30 (2019)

Secular spin-axis dynamics of exoplanets

Authors: Melaine Saillenfest, Jacques Laskar, Gwenaël Boué

Abstract: Context: Seasonal variations and climate stability of a planet are very sensitive to the planet obliquity and its evolution. This is of particular interest for the emergence and sustainability of land-based life, but orbital and rotational parameters of exoplanets are still poorly constrained. Numerical explorations usually realised in this situation are thus in heavy contrast with the uncertain n… ▽ More Context: Seasonal variations and climate stability of a planet are very sensitive to the planet obliquity and its evolution. This is of particular interest for the emergence and sustainability of land-based life, but orbital and rotational parameters of exoplanets are still poorly constrained. Numerical explorations usually realised in this situation are thus in heavy contrast with the uncertain nature of the available data. Aims: We aim to provide an analytical formulation of the long-term spin-axis dynamics of exoplanets, linking it directly to physical and dynamical parameters, but still giving precise quantitative results if the parameters are well known. Together with bounds for the poorly constrained parameters of exoplanets, this analysis is designed to allow a quick and straightforward exploration of the spin-axis dynamics. Methods: The long-term orbital solution is decomposed in quasi-periodic series and the spin-axis Hamiltonian is expanded in powers of eccentricity and inclination. Chaotic zones are measured by the resonance overlap criterion. Bounds for the poorly known parameters of exoplanets are obtained from physical grounds (rotational breakup) and dynamical considerations (equipartition of AMD). Results: This method gives accurate results when the orbital evolution is well known. The chaotic zones for planets of the Solar System can be retrieved in details from simple analytical formulas. For less constrained planetary systems, the maximal extent of the chaotic regions can be computed, requiring only the mass, the semi-major axis and the eccentricity of the planets present in the system. Additionally, some estimated bounds of the precession constant allow to classify which observed exoplanets are necessarily out of major spin-orbit secular resonances (unless the precession rate is affected by the presence of massive satellites). △ Less

Submitted 9 January, 2019; originally announced January 2019.

Comments: accepted by Astronomy & Astrophysics

Journal ref: A&A 623, A4 (2019)

Detectability of shape deformation in short-period exoplanets

Authors: B. Akinsanmi, S. C. C. Barros, N. C. Santos, A. C. M. Correia, P. F. L. Maxted, G. Boué, J. Laskar

Abstract: Short-period planets suffer from extreme tidal forces from their parent stars causing them to deform and attain non-spherical shapes. The planet shapes, modeled here as triaxial ellipsoids, can have an impact on the observed transit light curves and the planetary parameters derived. We investigate the detectability of tidal deformation in short-period planets from their transit light curves and th… ▽ More Short-period planets suffer from extreme tidal forces from their parent stars causing them to deform and attain non-spherical shapes. The planet shapes, modeled here as triaxial ellipsoids, can have an impact on the observed transit light curves and the planetary parameters derived. We investigate the detectability of tidal deformation in short-period planets from their transit light curves and the instrumental precision needed. We show how the detection of deformation from the light curve allows us to obtain an observational estimate of the second fluid Love number which gives valuable insight about the planet's internal structure. We adopted a model to calculate the shape of a planet due to the external potentials acting on it and used this model to modify the ellc transit tool. Our model is parameterized by the Love number, hence for a given light curve we can derive the value of the Love number that best matches the observations. We simulated the known cases of WASP-103b and WASP-121b expected to be highly deformed. Our analyses showed that instrumental precision $\leq$50ppm/min is needed to reliably estimate the Love number and detect tidal deformation. This precision can be achieved for WASP-103b in ~40 transits using HST and in ~300 transits using the forthcoming CHEOPS. However, fewer transits will be required for short-period planets that may be found around bright stars in the TESS and PLATO survey missions. The unprecedented precisions expected from PLATO and JWST can permit the detection of deformation with a single transit. However, the effects of instrumental and astrophysical noise must be well-considered as they can increase the number of transits required to reach the 50 ppm/min detection limit. We also show that improper modeling of limb darkening can act to bury signals related to the planet's shape thereby leading us to infer sphericity for a deformed planet. △ Less

Submitted 11 December, 2018; originally announced December 2018.

Comments: 9 pages, 10 figures, 2 tables, accepted for publication in A&A

Journal ref: A&A 621, A117 (2019)

Activity induced variation in spin-orbit angles as derived from Rossiter-McLaughlin measurements

Authors: M. Oshagh, A. H. M. J. Triaud, A. Burdanov, P. Figueira, A. Reiners, N. C. Santos, J. Faria, G. Boue, R. F. Diaz, S. Dreizler, S. Boldt, L. Delrez, E. Ducrot, M. Gillon, A. Guzman Mesa, E. Jehin, S. Khalafinejad, S. Kohl, L. Serrano, S. Udry

Abstract: One of the most powerful methods used to estimate sky-projected spin-orbit angles of exoplanetary systems is through a spectroscopic transit observation known as the Rossiter-McLaughlin (RM) effect. So far mostly single RM observations have been used to estimate the spin-orbit angle, and thus there have been no studies regarding the variation of estimated spin-orbit angle from transit to transit.… ▽ More One of the most powerful methods used to estimate sky-projected spin-orbit angles of exoplanetary systems is through a spectroscopic transit observation known as the Rossiter-McLaughlin (RM) effect. So far mostly single RM observations have been used to estimate the spin-orbit angle, and thus there have been no studies regarding the variation of estimated spin-orbit angle from transit to transit. Stellar activity can alter the shape of photometric transit light curves and in a similar way they can deform the RM signal. In this paper we discuss several RM observations, obtained using the HARPS spectrograph, of known transiting planets that all transit extremely active stars, and by analyzing them individually we assess the variation in the estimated spin-orbit angle. Our results reveal that the estimated spin-orbit angle can vary significantly (up to 42 degrees) from transit to transit, due to variation in the configuration of stellar active regions over different nights. This finding is almost two times larger than the expected variation predicted from simulations. We could not identify any meaningful correlation between the variation of estimated spin-orbit angles and the stellar magnetic activity indicators. We also investigated two possible approaches to mitigate the stellar activity influence on RM observations. The first strategy was based on obtaining several RM observations and folding them to reduce the stellar activity noise. Our results demonstrated that this is a feasible and robust way to overcome this issue. The second approach is based on acquiring simultaneous high-precision short-cadence photometric transit light curves using TRAPPIST/SPECULOOS telescopes, which provide more information about the stellar active region's properties and allow a better RM modeling. △ Less

Submitted 4 September, 2018; originally announced September 2018.

Comments: 18 pages, 7 figures, 6 tables, accepted for publication in A&A

Journal ref: A&A 619, A150 (2018)

Hill stability in the AMD framework

Authors: Antoine C. Petit, Jacques Laskar, Gwenaël Boué

Abstract: In a two-planet system, due to Sundman (1912) inequality, a topological boundary can forbid close encounters between the two planets for infinite time. A system is said Hill stable if it verifies this topological condition. Hill stability is widely used in the study of extra solar planets dynamics. However people often use the coplanar and circular orbits approximation. In this paper, we explain… ▽ More In a two-planet system, due to Sundman (1912) inequality, a topological boundary can forbid close encounters between the two planets for infinite time. A system is said Hill stable if it verifies this topological condition. Hill stability is widely used in the study of extra solar planets dynamics. However people often use the coplanar and circular orbits approximation. In this paper, we explain how the Hill stability can be understood in the framework of Angular Momentum Deficit (AMD). In the secular approximation, the AMD allows to discriminate between a priori stable systems and systems for which a more in depth dynamical analysis is required. We show that the general Hill stability criterion can be expressed as a function of only the semi major axes, the masses and the total AMD of the system. The proposed criterion is only expanded in the planets-to-star mass ratio $ε$ and not in the semi-major axis ratio, in eccentricities nor in the mutual inclination. Moreover the expansion in $ε$ remains excellent up to values of about $10^{-3}$ even for two planets with very different mass values. We performed numerical simulations in order to highlight the sharp change of behaviour between Hill stable and Hill unstable systems. We show that Hill stable systems tend to be very regular whereas Hill unstable ones often lead to rapid planet collisions. We also remind that Hill stability does not protect from the ejection of the outer planet. △ Less

Submitted 22 June, 2018; originally announced June 2018.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 617, A93 (2018)

Rotation of a rigid satellite with a fluid component. A new light onto Titan's obliquity

Authors: Gwenaël Boué, Nicolas Rambaux, Andy Richard

Abstract: We revisit the rotation dynamics of a rigid satellite with either a liquid core or a global sub-surface ocean. In both problems, the flow of the fluid component is assumed inviscid. The study of a hollow satellite with a liquid core is based on the Poincare-Hough model which provides exact equations of motion. We introduce an approximation when the ellipticity of the cavity is low. This simplifica… ▽ More We revisit the rotation dynamics of a rigid satellite with either a liquid core or a global sub-surface ocean. In both problems, the flow of the fluid component is assumed inviscid. The study of a hollow satellite with a liquid core is based on the Poincare-Hough model which provides exact equations of motion. We introduce an approximation when the ellipticity of the cavity is low. This simplification allows to model both types of satellite in the same manner. The analysis of their rotation is done in a non-canonical Hamiltonian formalism closely related to Poincare's "forme nouvelle des equations de la mecanique". In the case of a satellite with a global ocean, we obtain a seven-degree of freedom system. Six of them account for the motion of the two rigid components, and the last one is associated with the fluid layer. We apply our model to Titan for which the origin of the obliquity is still a debated question. We show that the observed value is compatible with Titan slightly departing from the hydrostatic equilibrium and being in a Cassini equilibrium state. △ Less

Submitted 14 September, 2017; originally announced September 2017.

Comments: 39 pages, 6 figures, accepted in CeMDA

AMD-stability in presence of first order Mean Motion Resonances

Authors: Antoine C. Petit, Jacques Laskar, Gwenaël Boué

Abstract: The AMD-stability criterion allows to discriminate between a-priori stable planetary systems and systems for which the stability is not granted and needs further investigations. AMD-stability is based on the conservation of the Angular Momentum Deficit (AMD) in the averaged system at all orders of averaging. While the AMD criterion is rigorous, the conservation of the AMD is only granted in absenc… ▽ More The AMD-stability criterion allows to discriminate between a-priori stable planetary systems and systems for which the stability is not granted and needs further investigations. AMD-stability is based on the conservation of the Angular Momentum Deficit (AMD) in the averaged system at all orders of averaging. While the AMD criterion is rigorous, the conservation of the AMD is only granted in absence of mean-motion resonances (MMR). Here we extend the AMD-stability criterion to take into account mean-motion resonances, and more specifically the overlap of first order MMR. If the MMR islands overlap, the system will experience generalized chaos leading to instability. The Hamiltonian of two massive planets on coplanar quasi-circular orbits can be reduced to an integrable one degree of freedom problem for period ratios close to a first order MMR. We use the reduced Hamiltonian to derive a new overlap criterion for first order MMR. This stability criterion unifies the previous criteria proposed in the literature and admits the criteria obtained for initially circular and eccentric orbits as limit cases. We then improve the definition of AMD-stability to take into account the short term chaos generated by MMR overlap. We analyze the outcome of this improved definition of AMD-stability on selected multi-planet systems from the Extrasolar Planets Encyclopeadia. △ Less

Submitted 21 July, 2017; v1 submitted 18 May, 2017; originally announced May 2017.

Comments: Accepted by A and A 07/10/17

Journal ref: A&A 607, A35 (2017)

Complete spin and orbital evolution of close-in bodies using a Maxwell viscoelastic rheology

Authors: Gwenaël Boué, Alexandre C. M. Correia, Jacques Laskar

Abstract: In this paper, we present a formalism designed to model tidal interaction with a viscoelastic body made of Maxwell material. Our approach remains regular for any spin rate and orientation, and for any orbital configuration including high eccentricities and close encounters. The method is to integrate simultaneously the rotation and the position of the planet as well as its deformation. We provide… ▽ More In this paper, we present a formalism designed to model tidal interaction with a viscoelastic body made of Maxwell material. Our approach remains regular for any spin rate and orientation, and for any orbital configuration including high eccentricities and close encounters. The method is to integrate simultaneously the rotation and the position of the planet as well as its deformation. We provide the equations of motion both in the body frame and in the inertial frame. With this study, we generalize preexisting models to the spatial case and to arbitrary multipole orders using a formalism taken from quantum theory. We also provide the vectorial expression of the secular tidal torque expanded in Fourier series. Applying this model to close-in exoplanets, we observe that if the relaxation time is longer than the revolution period, the phase space of the system is characterized by the presence of several spin-orbit resonances, even in the circular case. As the system evolves, the planet spin can visit different spin-orbit configurations. The obliquity is decreasing along most of these resonances, but we observe a case where the planet tilt is instead growing. These conclusions derived from the secular torque are successfully tested with numerical integrations of the instantaneous equations of motion on HD 80606b. Our formalism is also well adapted to close-in super-Earths in multiplanet systems which are known to have non-zero mutual inclinations. △ Less

Submitted 8 December, 2016; originally announced December 2016.

Comments: 31 pages, 5 figures

Journal ref: Celestial Mechanics and Dynamical Astronomy, Volume 126, Issue 1-3, pp. 31-60 (2016)

The Two Rigid Body Interaction using Angular Momentum Theory Formulae

Authors: Gwenaël Boué

Abstract: This work presents an elegant formalism to model the evolution of the full two rigid body problem. The equations of motion, given in a Cartesian coordinate system, are expressed in terms of spherical harmonics and Wigner D-matrices. The algorithm benefits from the numerous recurrence relations satisfied by these functions allowing a fast evaluation of the mutual potential. Moreover, forces and tor… ▽ More This work presents an elegant formalism to model the evolution of the full two rigid body problem. The equations of motion, given in a Cartesian coordinate system, are expressed in terms of spherical harmonics and Wigner D-matrices. The algorithm benefits from the numerous recurrence relations satisfied by these functions allowing a fast evaluation of the mutual potential. Moreover, forces and torques are straightforwardly obtained by application of ladder operators taken from the angular momentum theory and commonly used in quantum mechanics. A numerical implementation of this algorithm is made. Tests show that the present code is significantly faster than those currently available in literature. △ Less

Submitted 8 December, 2016; originally announced December 2016.

Comments: 13 pages, 1 figure, submitted to CeMDA on September 28, 2016

Radial-Velocity Fitting Challenge. II. First results of the analysis of the data set

Authors: X. Dumusque, F. Borsa, M. Damasso, R. Diaz, P. C. Gregory, N. C. Hara, A. Hatzes, V. Rajpaul, M. Tuomi, S. Aigrain, G. Anglada-Escude, A. S. Bonomo, G. Boue, F. Dauvergne, G. Frustagli, P. Giacobbe, R. D. Haywood, H. R. A. Jones, M. Pinamonti, E. Poretti, M. Rainer, D. Segransan, A. Sozzetti, S. Udry

Abstract: Radial-velocity (RV) signals induce RV variations an order of magnitude larger than the signal created by the orbit of Earth-twins, thus preventing their detection. The goal of this paper is to compare the efficiency of the different methods used to deal with stellar signals to recover extremely low-mass planets despite. However, because observed RV variations at the m/s precision level or below i… ▽ More Radial-velocity (RV) signals induce RV variations an order of magnitude larger than the signal created by the orbit of Earth-twins, thus preventing their detection. The goal of this paper is to compare the efficiency of the different methods used to deal with stellar signals to recover extremely low-mass planets despite. However, because observed RV variations at the m/s precision level or below is a combination of signals induced by unresolved orbiting planets, by the star, and by the instrument, performing such a comparison using real data is extremely challenging. To circumvent this problem, we generated simulated RV measurements including realistic stellar and planetary signals. Different teams analyzed blindly those simulated RV measurements, using their own method to recover planetary signals despite stellar RV signals. By comparing the results obtained by the different teams with the planetary and stellar parameters used to generate the simulated RVs, it is therefore possible to compare the efficiency of these different methods. The most efficient methods to recover planetary signals {take into account the different activity indicators,} use red-noise models to account for stellar RV signals and a Bayesian framework to provide model comparison in a robust statistical approach. Using the most efficient methodology, planets can be found down to K/N= K_pl/RV_rms*sqrt{N_obs}=5 with a threshold of K/N=7.5 at the level of 80-90% recovery rate found for a number of methods. These recovery rates drop dramatically for K/N smaller than this threshold. In addition, for the best teams, no false positives with K/N > 7.5 were detected, while a non-negligible fraction of them appear for smaller K/N. A limit of K/N = 7.5 seems therefore a safe threshold to attest the veracity of planetary signals for RV measurements with similar properties to those of the different RV fitting challenge systems. △ Less

Submitted 13 September, 2016; originally announced September 2016.

Comments: 36 pages (including 10 pages of appendix), 23 figures, Accepted in A&A

Journal ref: A&A 598, A133 (2017)

Radial Velocity Data Analysis with Compressed Sensing Techniques

Authors: Nathan C. Hara, Gwenaël Boué, Jacques Laskar, Alexandre C. M. Correia

Abstract: We present a novel approach for analysing radial velocity data that combines two features: all the planets are searched at once and the algorithm is fast. This is achieved by utilizing compressed sensing techniques, which are modified to be compatible with the Gaussian processes framework. The resulting tool can be used like a Lomb-Scargle periodogram and has the same aspect but with much fewer pe… ▽ More We present a novel approach for analysing radial velocity data that combines two features: all the planets are searched at once and the algorithm is fast. This is achieved by utilizing compressed sensing techniques, which are modified to be compatible with the Gaussian processes framework. The resulting tool can be used like a Lomb-Scargle periodogram and has the same aspect but with much fewer peaks due to aliasing. The method is applied to five systems with published radial velocity data sets: HD 69830, HD 10180, 55 Cnc, GJ 876 and a simulated very active star. The results are fully compatible with previous analysis, though obtained more straightforwardly. We further show that 55 Cnc e and f could have been respectively detected and suspected in early measurements from the Lick observatory and Hobby-Eberly Telescope available in 2004, and that frequencies due to dynamical interactions in GJ 876 can be seen. △ Less

Submitted 6 September, 2016; originally announced September 2016.

Comments: Accepted for publication in MNRAS

MSC Class: 62G

Secular and tidal evolution of circumbinary systems

Authors: Alexandre C. M. Correia, Gwenaël Boué, Jacques Laskar

Abstract: We investigate the secular dynamics of three-body circumbinary systems under the effect of tides. We use the octupolar non-restricted approximation for the orbital interactions, general relativity corrections, the quadrupolar approximation for the spins, and the viscous linear model for tides. We derive the averaged equations of motion in a simplified vectorial formalism, which is suitable to mode… ▽ More We investigate the secular dynamics of three-body circumbinary systems under the effect of tides. We use the octupolar non-restricted approximation for the orbital interactions, general relativity corrections, the quadrupolar approximation for the spins, and the viscous linear model for tides. We derive the averaged equations of motion in a simplified vectorial formalism, which is suitable to model the long-term evolution of a wide variety of circumbinary systems in very eccentric and inclined orbits. In particular, this vectorial approach can be used to derive constraints for tidal migration, capture in Cassini states, and stellar spin-orbit misalignment. We show that circumbinary planets with initial arbitrary orbital inclination can become coplanar through a secular resonance between the precession of the orbit and the precession of the spin of one of the stars. We also show that circumbinary systems for which the pericenter of the inner orbit is initially in libration present chaotic motion for the spins and for the eccentricity of the outer orbit. Because our model is valid for the non-restricted problem, it can also be applied to any three-body hierarchical system such as star-planet-satellite systems and triple stellar systems. △ Less

Submitted 11 August, 2016; originally announced August 2016.

Comments: 40 pages, 18 figures

Journal ref: Celestial Mechanics and Dynamical Astronomy 2016

Detecting ring systems around exoplanets using high resolution spectroscopy: the case of 51Pegb

Authors: N. C. Santos, J. H. C. Martins, G. Boué, A. C. M. Correia, M. Oshagh, P. Figueira, A. Santerne, S. G. Sousa, C. Melo, M. Montalto, I. Boisse, D. Ehrenreich, C. Lovis, F. Pepe, S. Udry, A. Garcia Munoz

Abstract: In this paper we explore the possibility that the recently detected reflected light signal of 51\,Peg\,b could be caused by a ring system around the planet. We use a simple model to compare the observed signal with the expected signal from a short-period giant planet with rings. We also use simple dynamical arguments to understand the possible geometry of such a system. We provide evidence that, t… ▽ More In this paper we explore the possibility that the recently detected reflected light signal of 51\,Peg\,b could be caused by a ring system around the planet. We use a simple model to compare the observed signal with the expected signal from a short-period giant planet with rings. We also use simple dynamical arguments to understand the possible geometry of such a system. We provide evidence that, to a good approximation, the observations are compatible with the signal expected from a ringed planet, assuming that the rings are non-coplanar with the orbital plane. However, based on dynamical arguments, we also show that this configuration is unlikely. In the case of coplanar rings we then demonstrate that the incident flux on the ring surface is about 2\% the value received by the planet, a value that renders the ring explanation unlikely. The results suggest that the signal observed cannot in principle be explained by a planet+ring system. We discuss, however, the possibility of using reflected light spectra to detect and characterize the presence of rings around short-period planets. Finally, we show that ring systems could have already been detected by photometric transit campaigns, but their signal could have been easily misinterpreted by the expected light curve of an eclipsing binary. △ Less

Submitted 2 September, 2015; originally announced September 2015.

Comments: Astronomy & Astrophysics, in press

Deformation and tidal evolution of close-in planets and satellites using a Maxwell viscoelastic rheology

Authors: Alexandre C. M. Correia, Gwenaël Boué, Jacques Laskar, Adrián Rodríguez

Abstract: In this paper we present a new approach to tidal theory. Assuming a Maxwell viscoelastic rheology, we compute the instantaneous deformation of celestial bodies using a differential equation for the gravity field coefficients. This method allows large eccentricities and it is not limited to quasi-periodic perturbations. It can take into account an extended class of perturbations, including chaotic… ▽ More In this paper we present a new approach to tidal theory. Assuming a Maxwell viscoelastic rheology, we compute the instantaneous deformation of celestial bodies using a differential equation for the gravity field coefficients. This method allows large eccentricities and it is not limited to quasi-periodic perturbations. It can take into account an extended class of perturbations, including chaotic motions and transient events. We apply our model to some already detected eccentric hot Jupiters and super-Earths in planar configurations. We show that when the relaxation time of the deformation is larger than the orbital period, spin-orbit equilibria arise naturally at half-integers of the mean motion, even for gaseous planets. In the case of super-Earths, these equilibria can be maintained for very low values of eccentricity. Our method can also be used to study planets with complex internal structures and other rheologies. △ Less

Submitted 7 November, 2014; originally announced November 2014.

Comments: 16 pages, 13 figures, 2 tables

Journal ref: A&A 571, A50 (2014)

Spin-orbit angle in compact planetary systems perturbed by an inclined companion. Application to the 55 Cancri system

Authors: Gwenaël Boué, Daniel C. Fabrycky

Abstract: The stellar spin orientation relative to the orbital planes of multiplanet systems are becoming accessible to observations. For example, 55 Cancri is a system composed of 5 planets orbiting a member of a stellar binary for which a projected obliquity of 72+-12 deg relative to the orbit of the innermost planet has been reported (Bourrier & Hebrard 2014). This large obliquity has been attributed to… ▽ More The stellar spin orientation relative to the orbital planes of multiplanet systems are becoming accessible to observations. For example, 55 Cancri is a system composed of 5 planets orbiting a member of a stellar binary for which a projected obliquity of 72+-12 deg relative to the orbit of the innermost planet has been reported (Bourrier & Hebrard 2014). This large obliquity has been attributed to the perturbation induced by the binary. Here we describe the secular evolution of similar systems and we discuss the case of the 55 Cancri system more deeply. We provide two different orbital configurations compatible with the currently available observations. △ Less

Submitted 20 October, 2014; originally announced October 2014.

Comments: Accepted for publication in the proceedings of IAU Symposium 310: Complex Planetary Systems

Improvements on analytic modelling of stellar spots

Authors: M. Montalto, G. Boué, M. Oshagh, I. Boisse, G. Bruno, N. C. Santos

Abstract: In this work we present the solution of the stellar spot problem using the Kelvin-Stokes theorem. Our result is applicable for any given location and dimension of the spots on the stellar surface. We present explicitely the result up to the second degree in the limb darkening law. This technique can be used to calculate very efficiently mutual photometric effects produced by eclipsing bodies occul… ▽ More In this work we present the solution of the stellar spot problem using the Kelvin-Stokes theorem. Our result is applicable for any given location and dimension of the spots on the stellar surface. We present explicitely the result up to the second degree in the limb darkening law. This technique can be used to calculate very efficiently mutual photometric effects produced by eclipsing bodies occulting stellar spots and to construct complex spot shapes. △ Less

Submitted 8 July, 2014; originally announced July 2014.

Comments: Resubmitted to MNRAS after accounting for minor comments of second review, 9 pages, 5 figures, software available at http://eduscisoft.com/KSINT/

Compact planetary systems perturbed by an inclined companion: II. Stellar spin-orbit evolution

Authors: Gwenaël Boué, Daniel Fabrycky

Abstract: The stellar spin orientation relative to the orbital planes of multiplanet systems are becoming accessible to observations. Here, we analyze and classify different types of spin-orbit evolution in compact multiplanet systems perturbed by an inclined outer companion. Our study is based on classical secular theory, using a vectorial approach developed in a separate paper. When planet-planet perturba… ▽ More The stellar spin orientation relative to the orbital planes of multiplanet systems are becoming accessible to observations. Here, we analyze and classify different types of spin-orbit evolution in compact multiplanet systems perturbed by an inclined outer companion. Our study is based on classical secular theory, using a vectorial approach developed in a separate paper. When planet-planet perturbations are truncated at the second order in eccentricity and mutual inclination, and the planet-companion perturbations are developed at the quadrupole order, the problem becomes integrable. The motion is composed of a uniform precession of the whole system around the total angular momentum, and in the rotating frame, the evolution is periodic. Here, we focus on the relative motion associated to the oscillations of the inclination between the planet system and the outer orbit, and of the obliquities of the star with respect to the two orbital planes. The solution is obtained using a powerful geometric method. With this technique, we identify four different regimes characterized by the nutation amplitude of the stellar spin-axis relative to the orbital plane of the planets. In particular, the obliquity of the star reaches its maximum when the system is in the Cassini regime where planets have more angular momentum than the star, and where the precession rate of the star is similar to that of the planets induced by the companion. In that case, spin-orbit oscillations exceed twice the inclination between the planets and the companion. Even if mutual inclination is only ~ 20 deg, this resonant case can cause the spin-orbit angle to oscillate between perfectly aligned and retrograde values. △ Less

Submitted 29 May, 2014; originally announced May 2014.

Comments: 29 pages, 15 figures, accepted for publication in ApJ

Compact planetary systems perturbed by an inclined companion: I. Vectorial representation of the secular model

Authors: Gwenaël Boué, Daniel Fabrycky

Abstract: The non-resonant secular dynamics of compact planetary systems are modeled by a perturbing function which is usually expanded in eccentricity and absolute inclination with respect to the invariant plane. Here, the expressions are given in a vectorial form which naturally leads to an expansion in eccentricity and mutual inclination. The two approaches are equivalent in most cases, but the vectorial… ▽ More The non-resonant secular dynamics of compact planetary systems are modeled by a perturbing function which is usually expanded in eccentricity and absolute inclination with respect to the invariant plane. Here, the expressions are given in a vectorial form which naturally leads to an expansion in eccentricity and mutual inclination. The two approaches are equivalent in most cases, but the vectorial one is specially designed for those where a quasi-coplanar system tilts as a whole by a large amount. Moreover, the vectorial expressions of the Hamiltonian and of the equations of motion are slightly simpler than those given in terms of the usual elliptical elements. We also provide the secular perturbing function in vectorial form expanded in semimajor axis ratio allowing for arbitrary eccentricities and inclinations. The interaction between the equatorial bulge of a central star and its planets is also provided, as is the relativistic periapse precession of any planet induced by the central star. We illustrate the use of this representation for following the secular oscillations of the terrestrial planets of the solar system, and for Kozai cycles as may take place in exoplanetary systems. △ Less

Submitted 29 May, 2014; originally announced May 2014.

Comments: 21 pages, 5 figures, accepted for publication in ApJ

Strong Dependence of the Inner Edge of the Habitable Zone on Planetary Rotation Rate

Authors: Jun Yang, Gwenael Boue, Daniel C. Fabrycky, Dorian S. Abbot

Abstract: Planetary rotation rate is a key parameter in determining atmospheric circulation and hence the spatial pattern of clouds. Since clouds can exert a dominant control on planetary radiation balance, rotation rate could be critical for determining mean planetary climate. Here we investigate this idea using a three-dimensional general circulation model with a sophisticated cloud scheme. We find that s… ▽ More Planetary rotation rate is a key parameter in determining atmospheric circulation and hence the spatial pattern of clouds. Since clouds can exert a dominant control on planetary radiation balance, rotation rate could be critical for determining mean planetary climate. Here we investigate this idea using a three-dimensional general circulation model with a sophisticated cloud scheme. We find that slowly rotating planets (like Venus) can maintain an Earth-like climate at nearly twice the stellar flux as rapidly rotating planets (like Earth). This suggests that many exoplanets previously believed to be too hot may actually be habitable, depending on their rotation rate. The explanation for this behavior is that slowly rotating planets have a weak Coriolis force and long daytime illumination, which promotes strong convergence and convection in the substellar region. This produces a large area of optically thick clouds, which greatly increases the planetary albedo. In contrast, on rapidly rotating planets a much narrower belt of clouds form in the deep tropics, leading to a relatively low albedo. A particularly striking example of the importance of rotation rate suggested by our simulations is that a planet with modern Earth's atmosphere, in Venus' orbit, and with modern Venus' (slow) rotation rate would be habitable. This would imply that if Venus went through a runaway greenhouse, it had a higher rotation rate at that time. △ Less

Submitted 19 April, 2014; originally announced April 2014.

Comments: 7 pages, 4 figures, accepted at Astrophysical Journal Letters

Journal ref: Yang et al. 2014 ApJ 787:L2

Anisotropic q-Gaussian velocity distributions in LambdaCDM halos

Authors: Leandro J. Beraldo e Silva, Gary A. Mamon, Manuel Duarte, Radoslaw Wojtak, Sébastien Peirani, Gwenaël Boué

Abstract: The velocity distribution function (VDF) of dark matter (DM) halos in $Λ$CDM dissipationless cosmological simulations, which must be non-separable in its radial and tangential components, is still poorly known. We present the first single-parameter, non-separable, anisotropic model for the VDF in $Λ$CDM halos, built from an isotropic $q$-Gaussian (Tsallis) VDF of the isotropic set of dimensionless… ▽ More The velocity distribution function (VDF) of dark matter (DM) halos in $Λ$CDM dissipationless cosmological simulations, which must be non-separable in its radial and tangential components, is still poorly known. We present the first single-parameter, non-separable, anisotropic model for the VDF in $Λ$CDM halos, built from an isotropic $q$-Gaussian (Tsallis) VDF of the isotropic set of dimensionless spherical velocity components (after subtraction of streaming motions), normalized by the respective velocity dispersions. We test our VDF on 90 cluster-mass halos of a dissipationless cosmological simulation. Beyond the virial radius, $r_{\rm vir}$, our model VDF adequately reproduces that measured in the simulated halos, but no $q$-Gaussian model can adequately represent the VDF within $r_{\rm vir}$, as the speed distribution function is then flatter-topped than any $q$-Gaussian can allow. Nevertheless, our VDF fits significantly better the simulations than the commonly used Maxwellian (Gaussian) distribution, at virtually all radii within $5\,r_{\rm vir}$. Within 0.4 (1) $r_{\rm vir}$, the non-Gaussianity index $q$ is (roughly) linearly related to the slope of the density profile and also to the velocity anisotropy profile. We provide a parametrization of the modulation of $q$ with radius for both the median fits and the fit of the stacked halo. At radii of a few percent of $r_{\rm vir}$, corresponding to the Solar position in the Milky Way, our best-fit VDF, although fitting better the simulations than the Gaussian one, overproduces significantly the fraction of high velocity objects, indicating that one should not blindly use these $q$-Gaussian fits to make predictions on the direct detection rate of DM particles. △ Less

Submitted 20 December, 2016; v1 submitted 24 October, 2013; originally announced October 2013.

Comments: This version consolidates the published version and the Erratum (changes in red)

Journal ref: MNRAS 452, 944 (2015)

Probing the effect of gravitational microlensing on the measurements of the Rossiter-McLaughlin effect

Authors: M. Oshagh, G. Boué, P. Figueira, N. C. Santos, N. Haghighipour

Abstract: In general, in the studies of transit light-curves and the Rossiter-McLaughlin (RM), the contribution of the planet's gravitational microlensing is neglected. Theoretical studies, have, however shown that the planet's microlensing can affect the transit light-curve and in some extreme cases cause the transit depth to vanish. In this letter, we present the results of our quantitative analysis of mi… ▽ More In general, in the studies of transit light-curves and the Rossiter-McLaughlin (RM), the contribution of the planet's gravitational microlensing is neglected. Theoretical studies, have, however shown that the planet's microlensing can affect the transit light-curve and in some extreme cases cause the transit depth to vanish. In this letter, we present the results of our quantitative analysis of microlening on the RM effect. Results indicate that for massive planets in on long period orbits, the planet's microlensing will have considerable contribution to the star's RV measurements. We present the details of our study, and discuss our analysis and results. △ Less

Submitted 26 August, 2013; originally announced August 2013.

Comments: 6 pages, 3 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A, 558, 65-70 (2013)

Effect of stellar spots on high-precision transit light-curve

Authors: M. Oshagh, N. C. Santos, I. Boisse, G. Boué, M. Montalto, X. Dumusque, N. Haghighipour

Abstract: Stellar activity features such as spots can create complications in determining planetary parameters through spectroscopic and photometric observations. The overlap of a transiting planet and a stellar spot, for instance, can produce anomalies in the transit light curves that may lead to inaccurate estimation of the transit duration, depth and timing. For instance such inaccuracies can affect the… ▽ More Stellar activity features such as spots can create complications in determining planetary parameters through spectroscopic and photometric observations. The overlap of a transiting planet and a stellar spot, for instance, can produce anomalies in the transit light curves that may lead to inaccurate estimation of the transit duration, depth and timing. For instance such inaccuracies can affect the precise derivation of the planet's radius. In this paper, we present the results of a quantitative study on the effects of stellar spots on high precision transit light curves. We show that spot anomalies can lead to the estimate of a planet radius that is 4% smaller than the real value. The effects on the transit duration can also be of the order of 4%, longer or shorter. Depending on the size and distribution of spots, anomalies can also produce transit timing variations with significant amplitudes. For instance, TTVs with signal amplitudes of 200 seconds can be produced when the spot is completely dark and as large as the largest Sun spot. Our study also indicates that the minimum size of a stellar spot with detectable affects on the high precision transit light curve is around 0.03 time the stellar radius for typical Kepler telescope precision. We also show that the strategy of including more free parameters (such as transit depth and duration) in the fitting procedure to measure the transit time of each individual transit will not produce accurate results in the case of active stars. △ Less

Submitted 18 July, 2013; v1 submitted 4 June, 2013; originally announced June 2013.

Comments: 7 pages, 9 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A, 556, 19-24 (2013)

Tidal damping of the mutual inclination in hierachical systems

Authors: A. C. M. Correia, G. Boue, J. Laskar, M. H. M. Morais

Abstract: Hierarchical two-planet systems, in which the inner body's semi-major axis is between 0.1 and 0.5 AU, usually present high eccentricity values, at least for one of the orbits. As a result of the formation process, one may expect that planetary systems with high eccentricities also have high mutual inclinations. However, here we show that tidal effects combined with gravitational interactions damp… ▽ More Hierarchical two-planet systems, in which the inner body's semi-major axis is between 0.1 and 0.5 AU, usually present high eccentricity values, at least for one of the orbits. As a result of the formation process, one may expect that planetary systems with high eccentricities also have high mutual inclinations. However, here we show that tidal effects combined with gravitational interactions damp the initial mutual inclination to modest values in timescales that are shorter than the age of the system. This effect is not a direct consequence of tides on the orbits, but it results from a secular forcing of the inner planet's flattening. We then conclude that these hierarchical planetary systems are unlikely to present very high mutual inclinations, at least as long as the orbits remain outside the Lidov-Kozai libration areas. The present study can also be extended to systems of binary stars and to planet-satellite systems. △ Less

Submitted 5 May, 2013; v1 submitted 4 March, 2013; originally announced March 2013.

Comments: 16 pages, 13 figures

Journal ref: Astronomy and Astrophysics 553, A39 (2013)

Transit timing variations in WASP-10b induced by stellar activity?

Authors: S. C. C. Barros, G. Boue, N. P. Gibson, D. Pollacco, A. Santerne, F. P. Keenan, I. Skillen, R. A. Street

Abstract: The hot-Jupiter WASP-10b was reported by Maciejewski et al. (2011a,b) to show transit timing variations (TTV) with an amplitude of ~ 3.5 minutes. These authors proposed that the observed TTVs were caused by a 0.1 MJup perturbing companion with an orbital period of ~ 5.23 d, and hence, close to the outer 5:3 mean motion resonance with WASP-10b. To test this scenario, we present eight new transit li… ▽ More The hot-Jupiter WASP-10b was reported by Maciejewski et al. (2011a,b) to show transit timing variations (TTV) with an amplitude of ~ 3.5 minutes. These authors proposed that the observed TTVs were caused by a 0.1 MJup perturbing companion with an orbital period of ~ 5.23 d, and hence, close to the outer 5:3 mean motion resonance with WASP-10b. To test this scenario, we present eight new transit light curves of WASP-10b obtained with the Faulkes Telescope North and the Liverpool Telescope. The new light curves, together with 22 previously published ones, were modelled with a Markov-Chain Monte-Carlo transit fitting code. (...) Our homogeneously derived transit times do not support the previous claimed TTV signal, which was strongly dependent on 2 previously published transits that have been incorrectly normalised. Nevertheless, a linear ephemeris is not a statistically good fit to the transit times of WASP-10b. We show that the observed transit time variations are due to spot occultation features or systematics. We discuss and exemplify the effects of occultation spot features in the measured transit times and show that despite spot occultation during egress and ingress being difficult to distinguish in the transit light curves, they have a significant effect in the measured transit times. We conclude that if we account for spot features, the transit times of WASP-10 are consistent with a linear ephemeris with the exception of one transit (epoch 143) which is a partial transit. Therefore, there is currently no evidence for the existence of a companion to WASP-10b. Our results support the lack of TTVs of hot-Jupiters reported for the Kepler sample. △ Less

Submitted 16 January, 2013; originally announced January 2013.

Comments: 17 pages, 14 figures. Accepted for publication in MNRAS

Impact of stellar companions on precise radial velocities

Authors: D. Cunha, P. Figueira, N. C. Santos, C. Lovis, G. Boué

Abstract: Context: With the announced arrival of instruments such as ESPRESSO one can expect that several systematic noise sources on the measurement of precise radial velocity will become the limiting factor instead of photon noise. A stellar companion within the fiber is such a possible noise source. Aims: With this work we aim at characterizing the impact of a stellar companion within the fiber to radial… ▽ More Context: With the announced arrival of instruments such as ESPRESSO one can expect that several systematic noise sources on the measurement of precise radial velocity will become the limiting factor instead of photon noise. A stellar companion within the fiber is such a possible noise source. Aims: With this work we aim at characterizing the impact of a stellar companion within the fiber to radial velocity measurements made by fiber-fed spectrographs. We consider the contaminant star either to be part of a binary system whose primary star is the target star, or as a background/foreground star. Methods: To carry out our study, we used HARPS spectra, co-added the target with contaminant spectra, and then compared the resulting radial velocity with that obtained from the original target spectrum. We repeated this procedure and used different tunable knobs to reproduce the previously mentioned scenarios. Results: We find that the impact on the radial velocity calculation is a function of the difference between individual radial velocities, of the difference between target and contaminant magnitude, and also of their spectral types. For the worst-case scenario in which both target and contaminant star are well centered on the fiber, the maximum contamination for a G or K star may be higher than 10 cm/s, on average, if the difference between target and contaminant magnitude is $Δm$ < 10, and higher than 1 m/s if $Δm$ < 8. If the target star is of spectral type M, $Δm$ < 8 produces the same contamination of 10 cm/s, and a contamination may be higher than 1 m/s △ Less

Submitted 12 December, 2012; originally announced December 2012.

Comments: Accepted for publication in A&A on 29/12/2019 - 14 pages

MAMPOSSt: Modelling Anisotropy and Mass Profiles of Observed Spherical Systems. I. Gaussian 3D velocities

Authors: Gary A. Mamon, Andrea Biviano, Gwenaël Boué

Abstract: Mass modelling of spherical systems through internal motions is hampered by the mass/velocity anisotropy (VA) degeneracy inherent in the Jeans equation, as well as the lack of techniques that are both fast and adaptable to realistic systems. A new fast method, called MAMPOSSt, which performs a maximum likelihood fit of the distribution of observed tracers in projected phase space, is developed and… ▽ More Mass modelling of spherical systems through internal motions is hampered by the mass/velocity anisotropy (VA) degeneracy inherent in the Jeans equation, as well as the lack of techniques that are both fast and adaptable to realistic systems. A new fast method, called MAMPOSSt, which performs a maximum likelihood fit of the distribution of observed tracers in projected phase space, is developed and thoroughly tested. MAMPOSSt assumes a shape for the gravitational potential, but instead of postulating a shape for the distribution function in terms of energy and angular momentum, or supposing Gaussian line-of-sight velocity distributions, MAMPOSSt assumes a VA profile and a shape for the 3D velocity distribution, here Gaussian. MAMPOSSt requires no binning, differentiation, nor extrapolation of the observables. Tests on cluster-mass haloes from LambdaCDM cosmological simulations show that, with 500 tracers, MAMPOSSt is able to jointly recover the virial radius, tracer scale radius, dark matter scale radius and outer or constant VA with small bias (<10% on scale radii and <2% on the two other quantities) and inefficiencies of 10%, 27%, 48% and 20%, respectively. MAMPOSSt does not perform better when some parameters are frozen, and even worse when the virial radius is set to its true value, which appears to be the consequence of halo triaxiality. The accuracy of MAMPOSSt depends weakly on the adopted interloper removal scheme, including an efficient iterative Bayesian scheme that we introduce here, which can directly obtain the virial radius with as good precision as MAMPOSSt. Our tests show that MAMPOSSt with Gaussian 3D velocities is very competitive with, and up to 1000x faster than other methods. Hence, MAMPOSSt is a very powerful and rapid tool for the mass and anisotropy modeling of systems such as clusters and groups of galaxies, elliptical and dwarf spheroidal galaxies. △ Less

Submitted 4 January, 2013; v1 submitted 6 December, 2012; originally announced December 2012.

Comments: MNRAS in press (fixes to eqs [20] and [30] plus minor touch-ups). 21 pages, 10 figures, 6 Tables

Journal ref: MNRAS 429, 3079 (2013)