We numerically explore the possibility that the large orbital inclination of the martian satellite Deimos originated in an orbital resonance with an ancient inner satellite of Mars more massive than Phobos. We find that Deimos's inclination can be reliably generated by outward evolution of a martian satellite that is about 20 times more massive than Phobos through the 3:1 mean-motion resonance with Deimos at 3.3 Mars radii. This outward migration, in the opposite direction from tidal evolution within the synchronous radius, requires interaction with a past massive ring of Mars. Our results therefore strongly support the cyclic martian ring-satellite hypothesis of Hesselbrock and Minton (2017). Our findings, combined with the model of Hesselbrock and Minton (2017), suggest that the age of the surface of Deimos is about 3.5-4 Gyr, and require Phobos to be significantly younger.

We numerically simulate the past tidal evolution of the five large moons of Uranus (Miranda, Ariel, Umbriel, Titania, and Oberon). We find that the most recent major mean-motion resonance (MMR) between any two moons, the Ariel-Umbriel 5:3 MMR, had a large effect on the whole system. Our results suggest that this resonance is responsible for the current 4.3$^{\circ}$ inclination of Miranda (instead of previously proposed 3:1 Miranda-Umbriel MMR), and that all five moons had their inclinations excited during this resonance. Miranda experienced significant tidal heating during the Ariel-Umbriel 5:3 MMR due to its eccentricity being excited by Ariel's secular perturbations. This tidal heating draws energy from shrinking of Miranda's orbit, rather than Ariel's outward evolution, and can generate heat flows in excess of 100 mW m$^{-2}$, sufficient to produce young coronae on Miranda. We find that this MMR was followed by a sequence of secular resonances, which reshuffled the moons' eccentricities and inclinations. We also find that the precession of Oberon's spin axis is close to a resonance with the precession of Umbriel's orbital plane, and that this spin-orbit resonance was likely excited during the Ariel-Umbriel 5:3 MMR. After the exit from the MMR, subsequent Ariel-Umbriel secular resonance and Oberon-Umbriel spin-orbit resonance may be able to explain the current low inclinations of Ariel and Umbriel. The age of Miranda's surface features tentatively suggests Uranian tidal $Q=15,000-20,000$, which can be further refined in future work.

The Moon is known to have a small liquid core, and it is thought that in the distant past the core may have produced strong magnetic fields recorded in lunar samples. Here we implement a numerical model of lunar orbital and rotational dynamics that includes the effects of a liquid core. In agreement with previous work, we find that the lunar core is dynamically decoupled from the lunar mantle, and that this decoupling happened very early in lunar history. Our model predicts that the lunar core rotates sub-synchronously, and the difference between the core and the mantle rotational rates was significant when the Moon had a high forced obliquity during and after the Cassini State transition. We find that the presence of the lunar liquid core further destabilizes synchronous rotation of the mantle for a wide range of semimajor axes centered around the Cassini State transition. CMB torques make it even more likely that the Moon experienced large-scale inclination damping during the Cassini State transition. We present estimates for the mutual core-mantle obliquity as a function of Earth-Moon distance, and we discuss plausible absolute time-lines for this evolution. We conclude that our results are consistent with the hypothesis of a precession-driven early lunar dynamo and may explain the variability of the inferred orientation of the past lunar dynamo.

Quantifying the energy budget of Earth in the first few million years following the Moon-forming giant impact is vital to understanding Earth's initial thermal state and the dynamics of lunar tidal evolution. After the impact, the body was substantially vaporized and rotating rapidly, very different from the planet we know today. The subsequent evolution of Earth's energy budget, as the body cooled and angular momentum was transferred during lunar tidal recession, has not been accurately calculated with all relevant energy components included. Here, we use giant impact simulations and planetary structure models to calculate the energy budget at stages in Earth's evolution. We show that the figure and internal structure of Earth changed substantially during its post-impact evolution and that changes in kinetic, potential, and internal energy were all significant. These changes have important implications for the dynamics of tidal recession and the thermal structure of early Earth.

Using numerical integrations, we find that the orbital eccentricity of Saturn's moon Iapetus undergoes prominent multi-Myr oscillations. We identify the responsible resonant argument to be $\varpi-\varpi_{g5}+\Omega-\Omega_{eq}$, with the terms being the longitudes of pericenter of Iapetus and planetary secular mode $g_5$, Iapetus's longitude of the node and Saturn's equinox. We find that this argument currently (on a $10^7$~yr timescale) appears to librate with a very large amplitude. On longer timescales, the behavior of this resonant angle is strongly dependent on the resonant interaction between Saturn's spin axis and the planetary mode $f_8$, with long-term secular resonance being possible if Saturn's equinox is librating relative to the node of the $f_8$ eigenmode. We present analytical estimates of the dependence of the resonant argument on the orbital elements of Iapetus. We find that this Iapetus-$g_5$ secular resonance could have been established only after the passage of Iapetus through the 5:1 mean-motion resonance with Titan, possibly in the last Gyr. Using numerical simulations, we show that the capture into the secular resonance appears to be a low-probability event. While the Iapetus-$g_5$ secular resonance can potentially help us put new constraints on the past dynamics of the Saturnian system, uncertainties in both the spin axis dynamics of Saturn and the tidal evolution rate of Titan make it impossible to make any firm conclusions about the resonance's longevity and origin.

The giant impact hypothesis remains the leading theory for lunar origin. However, current models struggle to explain the Moon's composition and isotopic similarity with Earth. Here we present a new lunar origin model. High-energy, high-angular momentum giant impacts can create a post-impact structure that exceeds the corotation limit (CoRoL), which defines the hottest thermal state and angular momentum possible for a corotating body. In a typical super-CoRoL body, traditional definitions of mantle, atmosphere and disk are not appropriate, and the body forms a new type of planetary structure, named a synestia. Using simulations of cooling synestias combined with dynamic, thermodynamic and geochemical calculations, we show that satellite formation from a synestia can produce the main features of our Moon. We find that cooling drives mixing of the structure, and condensation generates moonlets that orbit within the synestia, surrounded by tens of bars of bulk silicate Earth (BSE) vapor. The moonlets and growing moon are heated by the vapor until the first major element (Si) begins to vaporize and buffer the temperature. Moonlets equilibrate with BSE vapor at the temperature of silicate vaporization and the pressure of the structure, establishing the lunar isotopic composition and pattern of moderately volatile elements. Eventually, the cooling synestia recedes within the lunar orbit, terminating the main stage of lunar accretion. Our model shifts the paradigm for lunar origin from specifying a certain impact scenario to achieving a Moon-forming synestia. Giant impacts that produce potential Moon-forming synestias were common at the end of terrestrial planet formation.

In the giant impact hypothesis for lunar origin, the Moon accreted from an equatorial circum-terrestrial disk; however the current lunar orbital inclination of 5 degrees requires a subsequent dynamical process that is still debated. In addition, the giant impact theory has been challenged by the Moon's unexpectedly Earth-like isotopic composition. Here, we show that tidal dissipation due to lunar obliquity was an important effect during the Moon's tidal evolution, and the past lunar inclination must have been very large, defying theoretical explanations. We present a new tidal evolution model starting with the Moon in an equatorial orbit around an initially fast-spinning, high-obliquity Earth, which is a probable outcome of giant impacts. Using numerical modeling, we show that the solar perturbations on the Moon's orbit naturally induce a large lunar inclination and remove angular momentum from the Earth-Moon system. Our tidal evolution model supports recent high-angular momentum giant impact scenarios to explain the Moon's isotopic composition and provides a new pathway to reach Earth's climatically favorable low obliquity.

1I/`Oumuamua is the first known interstellar small body, probably being only about 100~m in size. Against expectations based on comets, `Oumuamua does not show any activity and has a very elongated figure, and also exhibits undamped rotational tumbling. In contrast, `Oumuamua's trajectory indicates that it was moving with the local stars, as expected from a low-velocity ejection from a relatively nearby system. Here I assume that `Oumuamua is typical of 100-m interstellar objects, and speculate on its origins. I find that giant planets are relatively inefficient at ejecting small bodies from inner solar systems of main-sequence stars, and that binary systems offer a much better opportunity for ejections of non-volatile bodies. I also conclude that `Oumuamua is not a member of a collisional population, which could explain its dramatic difference from small asteroids. I observe that 100-m small bodies are expected to carry little mass in realistic collisional populations, and that occasional events when whole planets are disrupted in catastrophic encounters may dominate interstellar population of 100-m fragments. Unlike the Sun or Jupiter, red dwarf stars are very dense and are capable of thoroughly tidally disrupting terrestrial planets. I conclude that the origin of `Oumuamua as a fragment from a planet that was tidally disrupted and then ejected by a dense member of a binary system could explain its peculiarities.

The Hungaria asteroid group is located interior to the main asteroid belt, with semimajor axes between 1.8 and 2 AU, low eccentricities and inclinations of 16-35 degrees. Recently, it has been proposed that Hungaria asteroids are a secularly declining population that may be related to the Late Heavy Bombardment (LHB) impactors (\'Cuk et al. 2012, Bottke et al. 2012). While \'Cuk et al. (2012) and Bottke et al. (2012) have reproduced a Hungaria-like population that declined exponentially, the real Hungarias were never confirmed to be unstable to the same degree. Here we find that the stability of Hungarias is strongly dependent on the evolution of the eccentricity of Mars, which is chaotic and unpredictable on Gyr timescales. We find that the high Martian eccentricity chiefly affects Hungarias through close approaches with Mars, rather than planetary secular modes. However, current minimum perihelia of Hungarias (over Myr timescales) are not diagnostic of their long-term stability due to a number of secular and mean motion resonances affecting the Hungaria region (Milani et al. 2010). We conclude that planetary chaos makes it impossible to determine the effective lifetimes of observed Hungarias. Furthermore, long-term changes of Martian eccentricity could lead to variable Hungaria loss over time. We speculate that some of the most stable Hungarias may have been placed in their present orbit when the eccentricity of Mars was significantly higher than today.

We explore the past evolution of Saturn's moons using direct numerical integrations. We find that the past Tethys-Dione 3:2 orbital resonance predicted in standard models likely did not occur, implying that the system is less evolved than previously thought. On the other hand, the orbital inclinations of Tethys, Dione and Rhea suggest that the system did cross the Dione-Rhea 5:3 resonance, which is closely followed by a Tethys-Dione secular resonance. A clear implication is that either the moons are significantly younger than the planet, or that their tidal evolution must be extremely slow (Q > 80,000). As an extremely slow-evolving system is incompatible with intense tidal heating of Enceladus, we conclude that the moons interior to Titan are not primordial, and we present a plausible scenario for the system's recent formation. We propose that the mid-sized moons re-accreted from a disk about 100 Myr ago, during which time Titan acquired its significant orbital eccentricity. We speculate that this disk has formed through orbital instability and massive collisions involving the previous generation of Saturn's mid-sized moons. We identify the solar evection resonance perturbing a pair of mid-sized moons as the most likely trigger of such an instability. This scenario implies that most craters on the moons interior to Titan must have been formed by planetocentric impactors.

Out of nine known stable Mars Trojans, seven appear to be members of an orbital grouping including the largest Trojan, Eureka. In order to test if this could be a genetic family, we simulated the long term evolution of a tight orbital cluster centered on Eureka. We explored two cases: cluster dispersal through planetary gravity alone over 1 Gyr, and a 1 Gyr evolution due to both gravity and the Yarkovsky effect. We find that the dispersal of the cluster in eccentricity is primarily due to dynamical chaos, while the inclinations and libration amplitudes are primarily changed by the Yarkovsky effect. Current distribution of the cluster members orbits are indicative of an initially tight orbital grouping that was affected by a negative acceleration (i.e. one against the orbital motion) consistent with the thermal Yarkovsky effect. We conclude that the cluster is a genetic family formed either in a collision or through multiple rotational fissions. The cluster's age is on the order of 1 Gyr, and its long-term orbital evolution is likely dominated by the seasonal, rather than diurnal, Yarkovsky effect. If confirmed, Gyr-scale dominance of the seasonal Yarkovsky effect may indicate suppression of the diurnal Yarkovsky drift by the related YORP effect. Further study of Mars Trojans is essential for understanding the long-term orbital and rotational dynamics of small bodies in the absence of frequent collisions.

The Hungaria asteroids are interior to the main asteroid belt, with semimajor axes between 1.8 and 2 AU, low eccentricities and inclinations of 16-35 degrees. Small asteroids in the Hungaria region are dominated by a collisional family associated with (434) Hungaria. The dominant spectral type of the Hungaria group is the E or X-type (Warner et al, 2009), mostly due to the E-type composition of Hungaria and its genetic family. It is widely believed the E-type asteroids are related to the aubrite meteorites, also known as enstatite achondrites (Gaffey et al, 1992). Here we explore the hypothesis that aubrites originate in the Hungaria family. In order to test this connection, we compare model Cosmic Ray Exposure ages from orbital integrations of model meteoroids with those of aubrites. We show that long CRE ages of aubrites (longest among stony meteorite groups) reflect the delivery route of meteoroids from Hungarias to Earth being different than those from main-belt asteroids. We find that the meteoroids from Hungarias predominantly reach Earth by Yarkovsky-drifting across the orbit of Mars, with no assistance from orbital resonances. We conclude that the CRE ages of aubrites are fully consistent with a dominant source at the inner boundary of the Hungaria family at 1.7 AU. From here, meteoroids reach Earth through the Mars-crossing region, with relatively quick delivery times favored due to collisions (with Hungarias and the inner main-belt objects). We find that, after Vesta, (434) Hungaria is the best candidate for an asteroidal source of an achondrite group.

Lainey et al. (2012), by re-analyzing long-baseline astrometry of Saturn's moons, have found that the moons' tidal evolution is much faster than previously thought, implying an order of magnitude stronger tidal dissipation within Saturn. This result is controversial and implies recent formation of at least some of the mid-sized icy moons of Saturn. Here we show that this more intensive tidal dissipation is in full agreement with the evolved state of the Titan-Hyperion resonance. This resonance was previously thought to be non-tidal in origin, as the amount of tidal evolution required for its assembly is beyond what is possible in models that assume that all the major moons are primordial. We find that the survival of the Titan-Hyperion resonance is in agreement with a past Titan-Iapetus 5:1 resonance, but not with unbroken tidal evolution of Rhea from the rings to its current distance.

The dwarf planet Haumea has two large satellites, Namaka and Hi'iaka, which orbit at relatively large separations. Both moons have significant eccentricities and inclinations, in a pattern that is consistent with a past orbital resonance (Ragozzine and Brown, 2009). Based on our analysis, we find that the present system is not consistent with satellite formation close to the primary and tidal evolution though mean-motion resonances. We propose that Namaka experienced only limited tidal evolution, leading to the mutual 8:3 mean-motion resonance which redistributed eccentricities and inclinations between the moons. This scenario requires that the original orbit of Hi'iaka was mildly eccentric; we propose that this eccentricity was either primordial or acquired though encounters with other TNOs. Both dynamical stability and our preferred tidal evolution model imply that the moons' masses are only about one half of previously estimated values, suggesting high albedos and low densities. As the present orbits of the moons strongly suggest formation from a flat disk close to their present locations, we conclude that Hi'iaka and Namaka may be second-generation moons, formed after the breakup of a past large moon, previously proposed as the parent body of the Haumea family (Schlichting and Sari, 2009). We derive plausible parameters of that moon, consistent with the current models of Haumea's formation (Leinhardt et al, 2010). An interesting implication of this hypothesis is that Hi'iaka and Namaka may orbit retrograde with respect to Haumea's spin. Retrograde orbits of Haumea's moons would be in full agreement with available observations and our dynamical analysis, and could provide a unique confirmation the "disrupted satellite" scenario for the origin of the family.

Unlike Trojans, horseshoe coorbitals are not generally considered to be long-term stable (Dermott and Murray, 1981; Murray and Dermott, 1999). As the lifetime of Earth's and Venus's horseshoe coorbitals is expected to be about a Gyr, we investigated the possible contribution of late-escaping inner planet coorbitals to the lunar Late Heavy Bombardment. Contrary to analytical estimates, we do not find many horseshoe objects escaping after first 100 Myr. In order to understand this behaviour, we ran a second set of simulations featuring idealized planets on circular orbits with a range of masses. We find that horseshoe coorbitals are generally long lived (and potentially stable) for systems with primary-to-secondary mass ratios larger than about 1200. This is consistent with results of Laughlin and Chambers (2002) for equal-mass pairs or coorbital planets and the instability of Jupiter's horseshoe companions (Stacey and Connors, 2008). Horseshoe orbits at smaller mass ratios are unstable because they must approach within 5 Hill radii of the secondary. In contrast, tadpole orbits are more robust and can remain stable even when approaching within 4 Hill radii of the secondary.

The Moon has suffered intense impact bombardment ending at 3.9 Gyr ago, and this bombardment probably affected all of the inner Solar System. Basin magnetization signatures and lunar crater size-distributions indicate that the last episode of bombardment at about 3.85 Gyr ago was less extensive than previously thought. We explore the contribution of the primordial Mars-crosser population to early lunar bombardment. We find that Mars-crosser population initially decays with a 80-Myr half-life, with the long tail of survivors clustering on temporarily non-Mars-crossing orbits between 1.8 and 2 AU. These survivors decay with half-life of about 600 Myr and are progenitors of the extant Hungaria asteroid group in the same region. We estimate the primordial Mars-crosser population contained about 0.01-0.02 Earth masses. Such initial population is consistent with no lunar basins forming after 3.8 Gya and the amount of mass in the Hungaria group. As they survive longer and in greater numbers than other primordial populations, Mars-crossers are the best candidate for forming the majority of lunar craters and basins, including most of the Nectarian system. However, this remnant population cannot produce Imbrium and Orientale basins, which formed too late and are too large to be part of a smooth bombardment. We propose that the Imbrian basins and craters formed in a discrete event, consistent with the basin magnetization signatures and crater size-distributions. This late "impactor shower" would be triggered by a collisional disruption of a Vesta-sized body from this primordial Mars-crossing population (Wetherill, 1975) that was still comparable to the present-day asteroid belt a 3.9 Gya. This tidal disruption lead to a short-lived spike in bombardment by non-chondritic impactors with a non-asteroidal size-frequency distribution, in agreement with available evidence. [abridged]