BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Gwenaël Boué

Gwenaël Boué appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.EPastro-ph.IMmath.DSphysics.geo-phastro-ph.COastro-ph.SRnlin.CDastro-ph.GAmath.NANumerical Analysisphysics.ao-ph

Catalog footprint

What is connected

22works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2022arXiv

The Resonant Tidal Evolution of the Earth-Moon Distance

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Constraining the Earth's Dynamical Ellipticity from Ice Age Dynamics

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Cassini states of a rigid body with a liquid core

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Tidal evolution of the Keplerian elements

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Anisotropic q-Gaussian velocity distributions in LambdaCDM halos

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Radial Velocity Data Analysis with Compressed Sensing Techniques

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Secular and tidal evolution of circumbinary systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A simple model of the chaotic eccentricity of Mercury

Mercury's eccentricity is chaotic and can increase so much that collisions with Venus or the Sun become possible (Laskar, 1989, 1990, 1994, 2008, Batygin & Laughlin, 2008, Laskar & Gastineau, 2009). This chaotic behavior results from an intricate network of secular resonances, but in this paper, we show that a simple integrable model with only one degree of freedom is actually able to reproduce the large variations in Mercury's eccentricity, with the correct amplitude and timescale. We show that this behavior occurs in the vicinity of the separatrices of the resonance g1-g5 between the precession frequencies of Mercury and Jupiter. However, the main contribution does not come from the direct interaction between these two planets. It is due to the excitation of Venus' orbit at Jupiter's precession frequency g5. We use a multipolar model that is not expanded with respect to Mercury's eccentricity, but because of the proximity of Mercury and Venus, the Hamiltonian is expanded up to order 20 and more in the ratio of semimajor axis. When the effects of Venus' inclination are added, the system becomes nonintegrable and a chaotic zone appears in the vicinity of the separatrices. In that case, Mercury's eccentricity can chaotically switch between two regimes characterized by either low-amplitude circulations or high-amplitude librations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Degeneracy in the characterization of non-transiting planets from transit timing variations

The transit timing variation (TTV) method allows the detection of non-transiting planets through their gravitational perturbations. Since TTVs are strongly enhanced in systems close to mean-motion resonances (MMR), even a low mass planet can produce an observable signal. This technique has thus been proposed to detect terrestrial planets. In this letter, we analyse TTV signals for systems in or close to MMR in order to illustrate the difficulties arising in the determination of planetary parameters. TTVs are computed numerically with an n-body integrator for a variety of systems close to MMR. The main features of these TTVs are also derived analytically. Systems deeply inside MMR do not produce particularly strong TTVs, while those close to MMR generate quasiperiodic TTVs characterised by a dominant long period term and a low amplitude remainder. If the remainder is too weak to be detected, then the signal is strongly degenerate and this prevents the determination of the planetary parameters. Even though an Earth mass planet can be detected by the TTV method if it is close to a MMR, it may not be possible to assert that this planet is actually an Earth mass planet. On the other hand, if the system is right in the center of a MMR, the high amplitude oscillation of the TTV signal vanishes and the detection of the perturber becomes as difficult as it is far from MMR.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

New analytical expressions of the Rossiter-McLaughlin effect adapted to different observation techniques

The Rossiter-McLaughlin (hereafter RM) effect is a key tool for measuring the projected spin-orbit angle between stellar spin axes and orbits of transiting planets. However, the measured radial velocity (RV) anomalies produced by this effect are not intrinsic and depend on both instrumental resolution and data reduction routines. Using inappropriate formulas to model the RM effect introduces biases, at least in the projected velocity Vsin(i) compared to the spectroscopic value. Currently, only the iodine cell technique has been modeled, which corresponds to observations done by, e.g., the HIRES spectrograph of the Keck telescope. In this paper, we provide a simple expression of the RM effect specially designed to model observations done by the Gaussian fit of a cross-correlation function (CCF) as in the routines performed by the HARPS team. We derived also a new analytical formulation of the RV anomaly associated to the iodine cell technique. For both formulas, we modeled the subplanet mean velocity v_p and dispersion beta_p accurately taking the rotational broadening on the subplanet profile into account. We compare our formulas adapted to the CCF technique with simulated data generated with the numerical software SOAP-T and find good agreement up to Vsin(i) < 20 km/s. In contrast, the analytical models simulating the two different observation techniques can disagree by about 10 sigma in Vsin(i) for large spin-orbit misalignments. It is thus important to apply the adapted model when fitting data.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Orbital migration induced by anisotropic evaporation. Can hot Jupiters form hot Neptunes ?

Short period planets are subject to intense energetic irradiations from their stars. It has been shown that this can lead to significant atmospheric mass-loss and create smaller mass planets. Here, we analyse whether the evaporation mechanism can affect the orbit of planets. The orbital evolution of a planet undergoing evaporation is derived analytically in a very general way. Analytical results are then compared with the period distribution of two classes of inner exoplanets: Jupiter-mass planets and Neptune-mass planets. These two populations have a very distinct period distribution, with a probability lower than 10^-4 that they were derived from the same parent distribution. We show that mass ejection can generate significant migration with an increase of orbital period that matches very well the difference of distribution of the two populations. This would happen if the evaporation emanates from above the hottest region of planet surface. Thus, migration induced by evaporation is an important mechanism that cannot be neglected.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Pumping the eccentricity of exoplanets by tidal effect

Planets close to their host stars are believed to undergo significant tidal interactions, leading to a progressive damping of the orbital eccentricity. Here we show that, when the orbit of the planet is excited by an outer companion, tidal effects combined with gravitational interactions may give rise to a secular increasing drift on the eccentricity. As long as this secular drift counterbalances the damping effect, the eccentricity can increase to high values. This mechanism may explain why some of the moderate close-in exoplanets are observed with substantial eccentricity values.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Tidal dissipation in multi-planet systems and constraints to orbit-fitting

We present here in full details the linear secular theory with tidal damping that was used to constraint the fit of the HD10180 planetary system in (Lovis et al. 2011). The theory is very general and can provide some intuitive understanding of the final state of a planetary system when one or more planets are close to their central star. We globally recover the results of (Mardling 2007), but we show that in the HD209458 planetary system, the consideration of the tides raised by the central star on the planet lead to believe that the eccentricity of HD209458b is most probably much smaller than 0.01.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Tidal evolution of hierarchical and inclined systems

We investigate the dynamical evolution of hierarchical three-body systems under the effect of tides, when the ratio of the orbital semi-major axes is small and the mutual inclination is relatively large (greater than 20 degrees). Using the quadrupolar non-restricted approximation for the gravitational interactions and the viscous linear model for tides, we derive the averaged equations of motion in a vectorial formalism which is suitable to model the long-term evolution of a large variety of exoplanetary systems in very eccentric and inclined orbits. In particular, it can be used to derive constraints for stellar spin-orbit misalignment, capture in Cassini states, tidal-Kozai migration, or damping of the mutual inclination. Because our model is valid for the non-restricted problem, it can be used to study systems of identical mass or for the outer restricted problem, such as the evolution of a planet around a binary of stars. Here, we apply our model to three distinct situations: 1) the HD80606 planetary system, for which we obtain the probability density function distribution for the misalignment angle, with two pronounced peaks of higher probability around 53 and 109 degrees; 2) the HD98800 binary system, for which we show that initial prograde orbits inside the observed disc may become retrograde and vice-versa, only because of tidal migration within the binary stars; 3) the HD11964 planetary system, for which we show that tidal dissipation combined with gravitational perturbations may lead to a decrease in the mutual inclination, and a fast circularization of the inner orbit.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A collisionless scenario for Uranus tilting

The origin of the high inclination of Uranus' spin-axis (Uranus' obliquity) is one of the great unanswered questions about the Solar system. Giant planets are believed to form with nearly zero obliquity, and it has been shown that the present behaviour of Uranus' spin is essentially stable. Several attempts were made in order to solve this problem. Here we report numerical simulations showing that Uranus' axis can be tilted during the planetary migration, without the need of a giant impact, provided that the planet had an additional satellite and a temporary large inclination. This might have happened during the giant planet instability phase described in the Nice model. In our scenario, the satellite is ejected after the tilt by a close encounter at the end of the migration. This model can both explain Uranus' large obliquity and bring new constraints on the planet orbital evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Speed limit on Neptune migration imposed by Saturn tilting

In this Letter, we give new constraints on planet migration. They were obtained under the assumption that Saturn's current obliquity is due to a capture in resonance with Neptune's ascending node. If planet migration is too fast, then Saturn crosses the resonance without being captured and it keeps a small obliquity. This scenario thus gives a lower limit on the migration time scale tau. We found that this boundary depends strongly on Neptune's initial inclination. For two different migration types, we found that tau should be at least greater than 7 Myr. This limit increases rapidly as Neptune's initial inclination decreases from 10 to 1 degree. We also give an algorithm to know if Saturn can be tilted for any migration law.

0 citations0 reviewsNo code linkedOpen access