The Effect of Tidal Interaction with a Gas Disk on Formation of Terrestrial Planets

We have performed N-body simulation on final accretion stage of terrestrial planets, including the effect of damping of eccentricity and inclination caused by tidal interaction with a remnant gas disk. As a result of runway and oligarchic accretion, about 20 Mars-sized protoplanets would be formed in nearly circular orbits with orbital separation of several to ten Hill radius. The orbits of the protoplanets would be eventually destabilized by long-term mutual gravity and/or secular resonance of giant gaseous planets. The protoplanets would coalesce with each other to form terrestrial planets through the orbital crossing. Previous N-body simulations, however, showed that the final eccentricities of planets are around 0.1, which are about 10 times higher than the present eccentricities of Earth and Venus. The obtained high eccentricities are the remnant of orbital crossing. We included the effect of eccentricity damping caused by gravitational interaction with disk gas as a drag force (“gravitational drag”) and carried out N-body simulation of accretion of protoplanets. We start with 15 protoplanets with 0.2M⊕ and integrate the orbits for 10⁷ years, which is consistent with the observationally inferred disk lifetime (in some runs, we start with 30 protoplanets with 0.1M⊕). In most runs, the damping time scale, which is equivalent to the strength of the drag force, is kept constant throughout each run in order to clarify the effects of the damping. We found that the planets' final mass, spatial distribution, and eccentricities depend on the damping time scale. If the damping time scale for a 0.2M⊕ mass planet at 1 AU is longer than 10⁸ years, planets grow to Earth's size, but the final eccentricities are too high as in gas-free cases. If it is shorter than 10⁶ years, the eccentricities of the protoplanets cannot be pumped up, resulting in not enough orbital crossing to make Earth-sized planets. Small planets with low eccentricities are formed with small orbital separation. On the other hand, if it is between 10⁶ and 10⁸ years, which may correspond to a mostly depleted disk (0.01-0.1% of surface density of the minimum mass model), some protoplanets can grow to about the size of Earth and Venus, and the eccentricities of such surviving planets can be diminished within the disk lifetime. Furthermore, in innermost and outermost regions in the same system, we often find planets with smaller size and larger eccentricities too, which could be analogous to Mars and Mercury. This is partly because the gravitational drag is less effective for smaller mass planets, and partly due to the “edge effect,” which means the innermost and outermost planets tend to remain without collision. We also carried out several runs with time-dependent drag force according to depletion of a gas disk. In these runs, we used exponential decay model with e-folding time of 3×10⁶ years. The orbits of protoplanets are stablized by the eccentricity damping in the early time. When disk surface density decays to ?1% of the minimum mass disk model, the damping force is no longer strong enough to inhibit the increase of the eccentricity by distant perturbations among protoplanets so that the orbital crossing starts. In this disk decay model, a gas disk with 10⁻⁴-10⁻³ times the minimum mass model still remains after the orbital crossing and accretional events, which is enough to damp the eccentricities of the Earth-sized planets to the order of 0.01. Using these results, we discuss a possible scenario for the last stage of terrestrial planet formation.     

Tidal Evolution of Low-mass Kepler Planets

The planets with a radius < 4 R_⊕ observed by the Kepler mission exhibit a unique feature, and propose a challenge for current planetary formation models. The tidal effect between a planet and its host star plays an essential role in reconfiguring the final orbits of the short-period planets. In this work, based on various initial Rayleigh distributions of the orbital elements, the final semi-major axis distributions of the planets with a radius < 4 R_⊕ after suffering tidal evolutions are investigated. Our simulations have qualitatively revealed some statistical properties: the semi-major axis and its peak value all increase with the increase of the initial semi-major axis and eccentricity. For the case that the initial mean semi-major axis is less than 0.1 au and the mean eccentricity is larger than 0.25, the results of numerical simulation are approximately consistent with the observation. In addition, the effects of other parameters, such as the tidal dissipation coefficient, stellar mass and planetary mass, etc., on the final semi-major axis distribution after tidal evolution are all relatively small. Based on the simulation results, we have tried to find some clues for the formation mechanism of low-mass planets. We speculate that these low-mass planets probably form in the far place of protoplanetary disk with a moderate eccentricity via the type I migration, and it is also possible to form in situ.     

Giant planet migration through the action of disk torques and planet-planet scattering

This paper presents a parametric study of giant planet migration through the combined action of disk torques and planet-planet scattering. The torques exerted on planets during Type II migration in circumstellar disks readily decrease the semi-major axes a, whereas scattering between planets increases the orbital eccentricities ?. This paper presents a parametric exploration of the possible parameter space for this migration scenario using two (initial) planetary mass distributions and a range of values for the time scale of eccentricity damping (due to the disk). For each class of systems, many realizations of the simulations are performed in order to determine the distributions of the resulting orbital elements of the surviving planets; this paper presents the results of ∼8500 numerical experiments. Our goal is to study the physics of this particular migration mechanism and to test it against observations of extrasolar planets. The action of disk torques and planet-planet scattering results in a distribution of final orbital elements that fills the a-? plane, in rough agreement with the orbital elements of observed extrasolar planets. In addition to specifying the orbital elements, we characterize this migration mechanism by finding the percentages of ejected and accreted planets, the number of collisions, the dependence of outcomes on planetary masses, the time spent in 2:1 and 3:1 resonances, and the effects of the planetary IMF. We also determine the distribution of inclination angles of surviving planets and the distribution of ejection speeds for exiled planets.     

Study of dynamical processes at the final stages of planetary system formation: Resonance motion of giant planets

According to current observational data, planets of many exoplanetary systems have resonant motion. The formation of resonance configurations is studied within a unified model of planetary migration. Planets in the observed systems 24 Sex, HD 37124, HD 73526, HD 82943, HD 128311, HD 160691, Kepler 9, NN Ser, which are moving in the 2: 1 resonance, could have been captured into this resonance due to both the Type I and II migration with a wide range of parameters. The migration conditions are defined for the formation of HD 45364 and HD 200964 that are in the 3: 2 and 4: 3 first-order resonances, correspondingly. The results obtained for HD 200964 show that planets can be captured in the first-order resonances, when the outer-to-inner orbital period ratios for the planets are less than 3: 2, only if Type I migration rates are large, and the mass of at least one planet is substantially less than the modern masses of the observed giant planets. The formation of the HD 102272, HD 108874, HD 181433 and HD 202206 systems with planets in high-order resonances is considered. The capture into these resonances can be realized with very slow Type II migration. Possible bounds for migration parameters are considered. In particular, it has been found that the capture of HD 108874 into the 4: 1 resonance is possible only if the angle between the plane of planetary orbits and the plane of sky is appreciably less than 90°, i.e., the planetary masses are a few times larger than the minimum values. The capture of HD 202206 into the 5: 1 resonance is possible at low migration rates; however, another mechanism is required to explain the high observed eccentricity of the inner planet (for example, strong gravitational interaction between the planets). Resonant configurations can be disrupted due to the interaction between planets and remaining fragments of the planetesimal disk as, for example, may occur in the three-planet system 47 UMa. The specific orbital features observed for this system are explained.     

Eccentricity evolution of giant planet orbits due to circumstellar disk torques

The extrasolar planets discovered to date possess unexpected orbital elements. Most orbit their host stars with larger eccentricities and smaller semi-major axes than similarly sized planets in our own Solar System do. It is generally agreed that the interaction between giant planets and circumstellar disks (Type II migration) drives these planets inward to small radii, but the effect of these same disks on orbital eccentricity, ?, is controversial. Several recent analytic calculations suggest that disk-planet interactions can excite eccentricity, while numerical studies generally produce eccentricity damping. This paper addresses this controversy using a quasi-analytic approach, drawing on several preceding analytic studies. This work refines the current treatment of eccentricity evolution by removing several approximations from the calculation of disk torques. We encounter neither uniform damping nor uniform excitation of orbital eccentricity, but rather a function d?/dt that varies in both sign and magnitude depending on eccentricity and other Solar System properties. Most significantly, we find that for every combination of disk and planet properties investigated herein, corotation torques produce negative values of d?/dt for some range in ? within the interval [0.1, 0.5]. If corotation torques are saturated, this region of eccentricity damping disappears, and excitation occurs on a short timescale of less than 0.08 Myr. Thus, our study does not produce eccentricity excitation on a timescale of a few Myr—we obtain either eccentricity excitation on a short time scale, or eccentricity damping on a longer time scale. Finally, we discuss the implications of this result for producing the observed range in extrasolar planet eccentricity.     

Orbital stability of systems of closely-spaced planets,II: configurations with coorbital planets

We numerically investigate the stability of systems of 1 \({{\rm M}_{\oplus}}\) planets orbiting a solar-mass star. The systems studied have either 2 or 42 planets per occupied semimajor axis, for a total of 6, 10, 126, or 210 planets, and the planets were started on coplanar, circular orbits with the semimajor axes of the innermost planets at 1 AU. For systems with two planets per occupied orbit, the longitudinal initial locations of planets on a given orbit were separated by either 60° (Trojan planets) or 180°. With 42 planets per semimajor axis, initial longitudes were uniformly spaced. The ratio of the semimajor axes of consecutive coorbital groups in each system was approximately uniform. The instability time for a system was taken to be the first time at which the orbits of two planets with different initial orbital distances crossed. Simulations spanned virtual times of up to 1 × 10⁸, 5 × 10⁵, and 2 × 10⁵ years for the 6- and 10-planet, 126-planet, and 210-planet systems, respectively. Our results show that, for a given class of system (e.g., five pairs of Trojan planets orbiting in the same direction), the relationship between orbit crossing times and planetary spacing is well fit by the functional form log(t _c /t ₀) = b β + c, where t _c is the crossing time, t ₀ = 1 year, β is the separation in initial orbital semimajor axis (in terms of the mutual Hill radii of the planets), and b and c are fitting constants. The same functional form was observed in the previous studies of single planets on nested orbits (Smith and Lissauer 2009). Pairs of Trojan planets are more stable than pairs initially separated by 180°. Systems with retrograde planets (i.e., some planets orbiting in the opposite sense from others) can be packed substantially more closely than can systems with all planets orbiting in the same sense. To have the same characteristic lifetime, systems with 2 or 42 planets per orbit typically need to have about 1.5 or 2 times the orbital separation as orbits occupied by single planets, respectively.     

On the stability of Earth-like planets in multi-planet systems

We present a continuation of our numerical study on planetary systems with similar characteristics to the Solar System. This time we examine the influence of three giant planets on the motion of terrestrial-like planets in the habitable zone (HZ). Using the Jupiter–Saturn–Uranus configuration we create similar fictitious systems by varying Saturn’s semi-major axis from 8 to 11 AU and increasing its mass by factors of 2–30. The analysis of the different systems shows the following interesting results: (i) Using the masses of the Solar System for the three giant planets, our study indicates a maximum eccentricity (max-e) of nearly 0.3 for a test-planet placed at the position of Venus. Such a high eccentricity was already found in our previous study of Jupiter–Saturn systems. Perturbations associated with the secular frequency g ₅ are again responsible for this high eccentricity. (ii) An increase of the Saturn-mass causes stronger perturbations around the position of the Earth and in the outer HZ. The latter is certainly due to gravitational interaction between Saturn and Uranus. (iii) The Saturn-mass increased by a factor 5 or higher indicates high eccentricities for a test-planet placed at the position of Mars. So that a crossing of the Earth’ orbit might occur in some cases. Furthermore, we present the maximum eccentricity of a test-planet placed in the Earth’ orbit for all positions (from 8 to 11 AU) and masses (increased up to a factor of 30) of Saturn. It can be seen that already a double-mass Saturn moving in its actual orbit causes an increase of the eccentricity up to 0.2 of a test-planet placed at Earth’s position. A more massive Saturn orbiting the Sun outside the 5:2 mean motion resonance (a _S  ≥9.7 AU) increases the eccentricity of a test-planet up to 0.4.     

Global dynamics of the Gliese 876 planetary system

The Gliese 876 planetary system consists of two Jupiter-like planets having a nearly commensurate 2:1 orbital periods ratio. Because the semimajor axes of the planets are very small (of the order 0.1 au and 0.2 au, respectively), and the eccentricity of the inner companion is ≃0.3, the mutual perturbations are extremely large. However, many authors claim the long-term orbital stability of the system, at least over 500 Myr for initial conditions found by Rivera & Lissauer. Results of investigations of a migration of initially separated planets into the close 2:1 mean motion resonance lock from Lee & Peale also support the conclusion that the system should be stable for the lifetime of the parent star. Initial conditions of the system, found from non-linear N -body fits by Laughlin & Chambers and Rivera & Lissauer, to the radial velocity curve, formally allow for a variety of orbital configurations of the GJ 876 system, e.g. coplanar, with planetary inclinations in the range [≃30°, 90°], and with relative inclinations of orbital planes as high as 80°. Our work is devoted to the stability investigation of the systems originating from the fitted initial conditions. We study neighbourhoods of these initial states in the orbital parameter space. We found estimations of the 2:1 mean motion resonance width and dynamical limitations on the planetary masses. We also obtain a global representation of the domains of the orbital parameters space in which initial conditions leading to stable evolutions can be found. Our results can be useful in localization of the best, stable fits to the observational data. In our investigations we use the MEGNO technique (the Mean Exponential Growth factor of Nearby Orbits) invented by Cincotta & Simó. It allows us to distinguish efficiently and precisely between chaotic and regular behaviour of a planetary system.     

On the Instability of Jupiter's Trojans

We have numerically explored the mechanisms that destabilize Jupiter's Trojan orbits outside the stability region defined by Levison et al. (1997, Nature385, 42-44). Different models have been exploited to test various possible sources of instability on timescales on the order of ∼10⁸ years.In the restricted three-body model, only a few Trojan orbits become unstable within 10⁸ years. This intrinsic instability contributes only marginally to the overall instability found by Levison et al.In a model where the orbital parameters of both Jupiter and Saturn are fixed, we have investigated the role of Saturn and its gravitational influence. We find that a large fraction of Trojan orbits become unstable because of the direct nonresonant perturbations by Saturn. By shifting its semimajor axis at constant intervals around its present value we find that the near 5:2 mean motion resonance between the two giant planets (the Great Inequality) is not responsible for the gross instability of Jupiter's Trojans since short-term perturbations by Saturn destabilize Trojans, even when the two planets are far out of the resonance.Secular resonances are an additional source of instability. In the full six-body model with the four major planets included in the numerical integration, we have analyzed the effects of secular resonances with the node of the planets. Trojan asteroids have relevant inclinations, and nodal secular resonances play an important role. When a Trojan orbit becomes unstable, in most cases the libration amplitude of the critical argument of the 1:1 mean motion resonance grows until the asteroid encounters the planet. Libration amplitude, eccentricity, and nodal rate are linked for Trojan orbits by an algebraic relation so that when one of the three parameters is perturbed, the other two are affected as well. There are numerous secular resonances with the nodal rate of Jupiter that fall inside the region of instability and contribute to destabilize Trojans, in particular the ν₁₆. Indeed, in the full model the escape rate over 50 Myr is higher compared to the fixed model.Some secular resonances even cross the stability region delimited by Levison et al. and cause instability. This is the case of the 3:2 and 1:2 nodal resonances with Jupiter. In particular the 1:2 is responsible for the instability of some clones of the L4 Trojan (3540) Protesilaos.     

On the Inclination Distribution of the Jovian Irregular Satellites

Irregular satellites—moons that occupy large orbits of significant eccentricity e and/or inclination I—circle each of the giant planets. The irregulars often extend close to the orbital stability limit, about 1/3-1/2 of the way to the edge of their planet's Hill sphere. The distant, elongated, and inclined orbits suggest capture, which presumably would give a random distribution of inclinations. Yet, no known irregulars have inclinations (relative to the ecliptic) between 47 and 141°.This paper shows that many high-I orbits are unstable due to secular solar perturbations. High-inclination orbits suffer appreciable periodic changes in eccentricity; large eccentricities can either drive particles with ∼70°<I<110° deep into the realm of the regular satellites (where collisions and scatterings are likely to remove them from planetocentric orbits on a timescale of 10⁷-10⁹ years) or expel them from the Hill sphere of the planet.By carrying out long-term (10⁹ years) orbital integrations for a variety of hypothetical satellites, we demonstrate that solar and planetary perturbations, by causing particles to strike (or to escape) their planet, considerably broaden this zone of avoidance. It grows to at least 55°<I<130° for orbits whose pericenters freely oscillate from 0 to 360°, while particles whose pericenters are locked at ±90° (Kozai mechanism) can remain for longer times.We estimate that the stable phase space (over 10 Myr) for satellites trapped in the Kozai resonance contains ∼10% of all stable orbits, suggesting the possible existence of a family of undiscovered objects at higher inclinations than those currently known.     

Eccentricities of Planets in Binary Systems

The most puzzling property of the extrasolar planets discovered by recent radial velocity surveys is their high orbital eccentricities, which are very difficult to explain within our current theoretical paradigm for planet formation. Current data reveal that at least 25% of these planets, including some with particularly high eccentricities, are orbiting a component of a binary star system. The presence of a distant companion can cause significant secular perturbations in the orbit of a planet. At high relative inclinations, large-amplitude, periodic eccentricity perturbations can occur. These are known as “Kozai cycles” and their amplitude is purely dependent on the relative orbital inclination. Assuming that every planet host star also has a (possibly unseen, e.g., substellar) distant companion, with reasonable distributions of orbital parameters and masses, we determine the resulting eccentricity distribution of planets and compare it to observations? We find that perturbations from a binary companion always appear to produce an excess of planets with both very high (?0.6) and very low (e ? 0.1) eccentricities. The paucity of near-circular orbits in the observed sample implies that at least one additional mechanism must be increasing eccentricities. On the other hand, the overproduction of very high eccentricities observed in our models could be combined with plausible circularization mechanisms (e.g., friction from residual gas) to create more planets with intermediate eccentricities (e? 0.1–0.6).     

Trajectories of satellites under the combined influences of earth oblateness and air drag

Trajectories of satellites under the influences of earth oblateness and air drag are derived by the asymptotic method in nonlinear mechanics. Based on the assumptions: (1) the dominant oblateness factor of the earth is the second harmonic (J ₂), (2) a non-rotating, spherically symmetric atmosphere and an exponential distribution of atmospheric density, (3) original elliptical orbits being of small eccentricity, closed-form solutions for the improved first order approximation are obtained. After finding the osculating orbital elements of the resulting trajectories, we expose the behavior of osculating orbits at various inclinations.     

The origin of TNO 2004 XR190 as a primordial scattered object

Numerical integrations of the equations of motion of the giant planets and scattering particles show that there is a possible orbital itinerary that a particle may follow from a scattering mode up to a stable position near the orbit of 2004 XR₁₉₀. This orbital evolution requires that the particle gets trapped in a mean motion resonance with Neptune coupled with the Kozai resonance. Imposing migration on Neptune while a particle is experiencing both resonances can entail an escape from resonance at a low particle’s eccentricity. This eccentricity and the associated inclination are always similar to those of 2004 XR₁₉₀. I conclude that 2004 XR₁₉₀ was most likely a scattered object that went through those resonance processes and was eventually deposited at its current position. By the same argument, it is expected that there must exist several other objects with similar semimajor axis, eccentricity and inclination as those of 2004 XR₁₉₀.     

Vertical instability and inclination excitation during planetary migration

We consider a two-planet system migrating under the influence of dissipative forces that mimic the effects of gas-driven (Type II) migration. It has been shown that, in the planar case, migration leads to resonant capture after an evolution that forces the system to follow families of periodic orbits. Starting with planets that differ slightly from a coplanar configuration, capture can, also, occur and, additionally, excitation of planetary inclinations has been observed in some cases. We show that excitation of inclinations occurs, when the planar families of periodic orbits, which are followed during the initial stages of planetary migration, become vertically unstable. At these points, vertical critical orbits may give rise to generating stable families of \(3D\) periodic orbits, which drive the evolution of the migrating planets to non-coplanar motion. We have computed and present here the vertical critical orbits of the \(2/1\) and \(3/1\) resonances, for various values of the planetary mass ratio. Moreover, we determine the limiting values of eccentricity for which the “inclination resonance” occurs.     

On periodic orbits and resonance in extrasolar planetary systems

We study the evolution of an extrasolar planetary system with two planets, for planar motion, starting from an exact resonant periodic motion and increasing the deviation from the equilibrium solution. We keep the semimajor axes and the eccentricities of the two planets fixed and we change the initial conditions by rotating the orbit of the outer planet by Δω. In this way the resonance is preserved, but we deviate from the exact periodicity and there is a transition from order to chaos as the deviation increases. There are three different routes to chaos, as far as the evolution of (ω ₂ ? ω ₁) is concerned: (a) Libration → rotation → chaos, with intermittent transition from libration to rotation in between, (b) libration → chaos and (c) libration → intermittent interchange between libration and rotation → chaos. This indicates that resonant planetary systems where the angle (ω ₂ ? ω ₁) librates or rotates are not different, but are closely connected to the exact periodic motion.     

Accretion of terrestrial planets from oligarchs in a turbulent disk

We have investigated the final accretion stage of terrestrial planets from Mars-mass protoplanets that formed through oligarchic growth in a disk comparable to the minimum mass solar nebula (MMSN), through N-body simulation including random torques exerted by disk turbulence due to Magneto-Rotational Instability. For the torques, we used the semi-analytical formula developed by Laughlin et al. [Laughlin, G., Steinacker, A., Adams, F.C., 2004. Astrophys. J. 608, 489-496]. The damping of orbital eccentricities (in all runs) and type-I migration (in some runs) due to the tidal interactions with disk gas is also included. Without any effect of disk gas, Earth-mass planets are formed in terrestrial planet regions in a disk comparable to MMSN but with too large orbital eccentricities to be consistent with the present eccentricities of Earth and Venus in our Solar System. With the eccentricity damping caused by the tidal interaction with a remnant gas disk, Earth-mass planets with eccentricities consistent with those of Earth and Venus are formed in a limited range of disk gas surface density (∼10⁻⁴ times MMSN). However, in this case, on average, too many (?6) planets remain in terrestrial planet regions, because the damping leads to isolation between the planets. We have carried out a series of N-body simulations including the random torques with different disk surface density and strength of turbulence. We found that the orbital eccentricities pumped up by the turbulent torques and associated random walks in semimajor axes tend to delay isolation of planets, resulting in more coagulation of planets. The eccentricities are still damped after planets become isolated. As a result, the number of final planets decreases with increase in strength of the turbulence, while Earth-mass planets with small eccentricities are still formed. In the case of relatively strong turbulence, the number of final planets are 4-5 at 0.5-2 AU, which is more consistent with Solar System, for relatively wide range of disk gas surface density (∼¹⁰⁻⁴-¹⁰⁻² times MMSN).     

The large-scale structure of the asteroid belt

We examine the distributions of 2888 numbered minor planets over orbital inclination, eccentricity, and semimajor axis, and define 19 zones which we believe adequately to isolate the selection biases in survey programs of the physical properties of minor planets. Six numbered asteroids have exceptional orbits and fall into no zone. We also call attention to rather sharp upper limits, which become increasingly stringent at larger heliocentric distances, on orbital inclinations and eccentricity.     

Stability of inclined orbits of terrestrial planets in habitable zones

In long-term stability studies of terrestrial planets moving in the habitable zone (HZ) of a sun-like star, we distinguish four different configurations: (i) planets moving in binary star systems, (ii) the inner type (where the gas giant moves outside the HZ), (iii) the outer type (where the gas giant is closer to the star, than the HZ) and (iv) the Trojan type (where the gas giant moves in the HZ). Since earlier calculations indicated, that the stability of the motion in the HZ also depends on the inclination of the terrestrial planet orbits, we present a detailed numerical investigation to show correlations between the eccentricity, the mass and the distance of the giant planet for various inclinations of the terrestrial planets. The orbital stability of the HZ was examined for all four configurations stated above. While we could find hardly any stable orbits for the first three types for inclinations higher than 40^°, the Trojan planets can be stable up to an inclination of 60^°. Additionally, we could also find some stabilizing effects of the inclination for the first three types. As dynamical model we used the elliptic restricted three-body problem, which consists of two massive and one mass-less body. This allows an application to all detected and future extrasolar single planet systems.     

Titan’s rotational state

In Noyelles et al. (Astron. Astrophys. 478, 959–970 (2008)), a resonance involving the wobble of Titan is hinted at. This paper studies this scenario and its consequences. The first step is to build an accurate analytical model that would help to find the likely resonances in the rotation of every synchronous body. In this model, I take the orbital eccentricity of the body into account, as well as its variable inclination with respect to Saturn’s equator. Then an analytical study using the second fundamental model of the resonance is performed to study the resonance of interest. Finally, I study the dissipative consequences of this resonance. I find that this resonance may have increased the wobble of Titan by several degrees. For instance, if Titan’s polar momentum C is equal to 0.355MR _T ² (M and R _T being, respectively, Titan’s mass and radius), the wobble might be forced to 41 degrees. Thanks to an original formula, I find that the dissipation associated with the forced wobble might not be negligible compared to the contribution of the eccentricity. I also suspect that, due to the forced wobble, Titan’s period of rotation may be somewhat underestimated by observers. Finally, I use the analytical model presented in this paper to compute the periods of the free librations of the four Galilean satellites as well as the Saturnian satellite Rhea. For Io and Europa, the results are consistent with previous studies. For the other satellites, the periods of the free librations are, respectively, 186.37 d, 23.38 y and 30.08 y for Ganymede, 2.44 y, 209.32 y and 356.54 y for Callisto, and 51.84 d, 2.60 y and 3.59 y for Rhea.     

Survival of Trojan-type companions of Neptune during primordial planet migration

We investigate the survivability of Trojan-type companions of Neptune during primordial radial migration of the giant planets Jupiter, Saturn, Uranus, and Neptune. We adopt the usual planet migration model in which the migration speed decreases exponentially with a characteristic time scale τ (the e-folding time). We perform a series of numerical simulations, each involving the migrating giant planets plus ∼1000 test particle Neptune Trojans with initial distributions of orbital eccentricity, inclination, and libration amplitude similar to those of the known jovian Trojans asteroids. We analyze these simulations to measure the survivability of Neptune's Trojans as a function of migration rate. We find that orbital migration with the characteristic time scale τ=10⁶ years allows about 35% of preexisting Neptune Trojans to survive to 5τ, by which time the giant planets have essentially reached their final orbits. In contrast, slower migration with τ=10⁷ years yields only a ∼5% probability of Neptune Trojans surviving to a time of 5τ. Interestingly, we find that the loss of Neptune Trojans during planetary migration is not a random diffusion process. Rather, losses occur almost exclusively during discrete prolonged episodes when Trojan particles are swept by secondary resonances associated with mean-motion commensurabilities of Uranus with Neptune. These secondary resonances arise when the circulation frequencies, f, of critical arguments for Uranus-Neptune mean-motion near-resonances (e.g., f^UN_1:2, f^UN_4:7) are commensurate with harmonics of the libration frequency of the critical argument for the Neptune-Trojan 1:1 mean-motion resonance (f^NT_1:1). Trojans trapped in the secondary resonances typically have their libration amplitudes amplified until they escape the 1:1 resonance with Neptune. Trojans with large libration amplitudes are susceptible to loss during sweeping by numerous high-order secondary resonances (e.g., f^UN_1:2≈11f^NT_1:1). However, for the slower migration, with τ=10⁷ years, even tightly bound Neptune Trojans with libration amplitudes below 10° can be lost when they become trapped in 1:3 or 1:2 secondary resonances between f^UN_1:2 and f^NT_1:1. With τ=10⁷ years the 1:2 secondary resonance was responsible for the single greatest episode of loss, ejecting nearly 75% of existing Neptune Trojans. This episode occurred during the late stages of planetary migration when the remnant planetesimal disk would have been largely dissipated. We speculate that if the number of bodies liberated during this event was sufficiently high they could have caused a spike in the impact rate throughout the Solar System.