Source author record

David Nesvorný

David Nesvorný appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

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

20 published item(s)

preprint2026arXiv

A young progenitor for the most common planetary systems in the Galaxy

The Galaxy's most common known planetary systems have several Earth-to-Neptune-size planets in compact orbits. At small orbital separations, larger planets are less common than their smaller counterparts by an order of magnitude. The young star V1298 Tau hosts one such compact planetary system, albeit with four planets that are uncommonly large (5 to 10 Earth radii). The planets form a chain of near-resonances that result in transit-timing variations of several hours. Here we present a multi-year campaign to characterize this system with transit-timing variations, a method insensitive to the intense magnetic activity of the star. Through targeted observations, we first resolved the previously unknown orbital period of the outermost planet. The full 9-year baseline from these and archival data then enabled robust determination of the masses and orbital parameters for all four planets. We find the planets have low, sub-Neptune masses and nearly circular orbits, implying a dynamically tranquil history. Their low masses and large radii indicate that the inner planets underwent a period of rapid cooling immediately after dispersal of the protoplanetary disk. Still, they are much less dense than mature planets of comparable size. We predict the planets will contract to 1.5-4.0 Earth radii and join the population of super-Earths and sub-Neptunes that nature produces in abundance.

preprint2026arXiv

Two-source terrestrial planet formation with a sweeping secular resonance

The models that most successfully reproduce the orbital architecture of the Solar System terrestrial planets start from a narrow annulus of material that grows into embryos and then planets. However, it is not clear how this ring model can be made consistent with the chemical structure of the inner Solar System, which shows a reduced-to-oxidized gradient from Mercury to Mars and a parallel gradient in the asteroid belt. We propose that there were two primary reservoirs in the early inner Solar System: a narrow, refractory enriched ring inside of 1 au, and a less massive, extended planetesimal disk outside of 1 au with oxidation states ranging from enstatite chondrites to ordinary chondrites. We show through a suite of N-body simulations that an inwardly sweeping secular resonance, caused by aerodynamic drag and perturbations from a mean-motion resonant Jupiter and Saturn, gathers the outer planetesimal disk into a narrow ring that migrates radially, forms Mars, and contributes oxidized material to proto-Earth. Remaining unaccreted planetesimals can be implanted into the asteroid belt as the parent bodies of aubrites and non-carbonaceous iron meteorites, while the most reduced material is not implanted and thus unsampled in the meteorite collection. This model explains the oxidation and isotopic gradients within the inner Solar System within the context of a low-viscosity, magnetic wind-driven disk.

preprint2022arXiv

Exogenous delivery of water to Mercury

Radar and spacecraft observations show the permanently shadowed regions around Mercury's North Pole to contain water ice and complex organic material. One possible source of this material are impacts by interplanetary dust particles (IDPs), asteroids, and comets. We have performed numerical simulations of the dynamical evolution of asteroids and comets over the few Myr and checked for their impacts with Mercury. We use the N-body integrator RMVS/Swifter to propagate the Sun and the eight planets from their current positions. We add comets and asteroids to the simulations as massless test particles, based on their current orbital distributions. Asteroid impactors are assigned a probability of being water-rich (C-class) based on the measured distribution of taxonomic types. For comets, we assume a constant water fraction. For IDPs, we use a dynamical meteoroid model to compute the dust flux on Mercury. Relative to previous work on asteroid and comet impacts (Moses et al. 1999), we leverage 20 years of progress in minor body surveys. Immediate post-impact ejection of impactor material into outer space is taken into account as is the migration efficiency of water across Mercury's surface to the polar cold traps. We find that asteroids deliver $\sim 1 \times 10^{3}$ kg/yr of water to Mercury, comets deliver $\sim 1 \times 10^{3}$ kg/yr and IDPs deliver $\sim 16 \times 10^{3}$ kg/yr within a factor of several. Over a timescale of $\sim 1$ Gyr, this is enough to deliver the minimum amount of water required by the radar and MESSENGER observations. While other sources of water on Mercury are not ruled out by our analysis, we show that they are not required to explain the currently available observational lower limits.

preprint2022arXiv

HD 28109 hosts a trio of transiting Neptunian planets including a near-resonant pair, confirmed by ASTEP from Antarctica

We report on the discovery and characterisation of three planets orbiting the F8 star HD~28109, which sits comfortably in \tess's continuous viewing zone. The two outer planets have periods of $\rm 56.0067 \pm 0.0003~days$ and $\rm 84.2597_{-0.0008}^{+0.0010}~days$, which implies a period ratio very close to that of the first-order 3:2 mean motion resonance, exciting transit timing variations (TTVs) of up to $\rm 60\,mins$. These two planets were first identified by \tess, and we identified a third planet in the \textcolor{black}{\tess photometry} with a period of $\rm 22.8911 \pm 0.0004~days$. We confirm the planetary nature of all three planetary candidates using ground-based photometry from Hazelwood, ASTEP and LCO, including a full detection of the $\rm \sim9\,h$ transit of HD~28109 c from Antarctica. The radii of the three planets are \textcolor{black}{$\rm R_b=2.199_{-0.10}^{+0.098} ~R_{\oplus}$, $\rm R_c=4.23\pm0.11~ R_{\oplus}$ and $\rm R_d=3.25\pm0.11 ~R_{\oplus}$}; we characterise their masses using TTVs and precise radial velocities from ESPRESSO and HARPS, and find them to be $\rm M_b=18.5_{-7.6}^{+9.1}~M_{\oplus}$, $\rm M_c=7.9_{-3.0}^{+4.2}~M_{\oplus}$ and $\rm M_d=5.7_{-2.1}^{+2.7}~M_{\oplus}$, making planet b a dense, massive planet while c and d are both under-dense. We also demonstrate that the two outer planets are ripe for atmospheric characterisation using transmission spectroscopy, especially given their position in the CVZ of JWST. The data obtained to date are consistent with resonant (librating) and non-resonant (circulating) solutions; additional observations will show whether the pair is actually locked in resonance or just near-resonant.

preprint2022arXiv

Thermal processing of Jupiter Family Comets during their chaotic orbital evolution

Evidence for cometary activity beyond Jupiter and Saturn's orbits -- such as that observed for Centaurs and long period comets -- suggests that the thermal processing of comet nuclei starts long before they enter the inner Solar System, where they are typically observed and monitored. Such observations raise questions as to the depth of unprocessed material, and whether the activity of JFCs can be representative of any primitive material. Here we model the coupled thermal and dynamical evolution of Jupiter Family Comets (JFCs), from the moment they leave their outer Solar System reservoirs until their ejection into interstellar space. We apply a thermal evolution model to a sample of simulated JFCs obtained from dynamical simulations (arXiv:1706.07447) that successfully reproduce the orbital distribution of observed JFCs. We show that due to the stochastic nature of comet trajectories toward the inner solar system, all simulated JFCs undergo multiple heating episodes resulting in significant modifications of their initial volatile contents. A statistical analysis constrains the extent of such processing. We suggest that primordial condensed hypervolatile ices should be entirely lost from the layers that contribute to cometary activity observed today. Our results demonstrate that understanding the orbital (and thus, heating) history of JFCs is essential when putting observations in a broader context.

preprint2020arXiv

Modeling the chronologies and size distributions of Ceres and Vesta craters

We infer the crater chronologies of Ceres and Vesta from a self-consistent dynamical model of asteroid impactors. The model accounts for planetary migration/instability early in the solar system history and tracks asteroid orbits over 4.56 Gy. It is calibrated on the current population of the asteroid belt. The model provides the number of asteroid impacts on different worlds at any time throughout the solar system history. We combine the results with an impactor-crater scaling relationship to determine the crater distribution of Ceres and Vesta and compare these theoretical predictions with observations. We find that: (i) The Ceres and Vesta chronologies are similar, whereas they significantly differ from the lunar chronology. Therefore, using the lunar chronology for main belt asteroids, as often done in previous publications, is incorrect. (ii) The model results match the number and size distribution of large (diameter $>90$ km) craters observed on Vesta, but overestimate the number of large craters on Ceres. This implies that large crater erasure is required for Ceres. (iii) In a model where planetary migration/instability happens early, the probability to form the Rheasilvia basin on Vesta during the last 1 Gy is 10\%, a factor of $\sim1.5$ higher than for the late instability case and $\sim2.5$ times higher than found in previous studies. Thus, while the formation of the Rheasilvia at $\sim1$ Gy ago (Ga) would be somewhat unusual, it cannot be ruled out at more than $\simeq1.5σ$. In broader context, our work provides a self-consistent framework for modeling asteroid crater records.

preprint2016arXiv

Characterizing the original ejection velocity field of the Koronis family

An asteroid family forms as a result of a collision between an impactor and a parent body. The fragments with ejection speeds higher than the escape velocity from the parent body can escape its gravitational pull. The cloud of escaping debris can be identified by the proximity of orbits in proper element, or frequency, domains. Obtaining estimates of the original ejection speed can provide valuable constraints on the physical processes occurring during collision, and used to calibrate impact simulations. Unfortunately, proper elements of asteroids families are modified by gravitational and non-gravitational effects, such as resonant dynamics, encounters with massive bodies, and the Yarkovsky effect, such that information on the original ejection speeds is often lost, especially for older, more evolved families. It has been recently suggested that the distribution in proper inclination of the Koronis family may have not been significantly perturbed by local dynamics, and that information on the component of the ejection velocity that is perpendicular to the orbital plane ($v_W$), may still be available, at least in part. In this work we estimate the magnitude of the original ejection velocity speeds of Koronis members using the observed distribution in proper eccentricity and inclination, and accounting for the spread caused by dynamical effects. Our results show that i) the spread in the original ejection speeds is, to within a 15% error, inversely proportional to the fragment size, and ii) the minimum ejection velocity is of the order of 50 m/s, with larger values possible depending on the orbital configuration at the break-up.

preprint2016arXiv

Constraints on the original ejection velocity fields of asteroid families

Asteroid families form as a result of large-scale collisions among main belt asteroids. The orbital distribution of fragments after a family-forming impact could inform us about their ejection velocities. Unfortunately, however, orbits dynamically evolve by a number of effects, including the Yarkovsky drift, chaotic diffusion, and gravitational encounters with massive asteroids, such that it is difficult to infer the ejection velocities eons after each family's formation. Here we analyze the inclination distribution of asteroid families, because proper inclination can remain constant over long time intervals, and could help us to understand the distribution of the component of the ejection velocity that is perpendicular to the orbital plane ($v_{W}$). From modeling the initial breakup, we find that the distribution of $v_{W}$ of the fragments, which manage to escape the parent body's gravity, should be more peaked than a Gaussian distribution (i.e., be leptokurtic) even if the initial distribution was Gaussian. We surveyed known asteroid families for signs of a peaked distribution of $v_{W}$ using a statistical measure of the distribution peakedness or flatness known as kurtosis. We identified eight families whose $v_{W}$ distribution is significantly leptokurtic. These cases (e.g. the Koronis family) are located in dynamically quiet regions of the main belt, where, presumably, the initial distribution of $v_{W}$ was not modified by subsequent orbital evolution. We suggest that, in these cases, the inclination distribution can be used to obtain interesting information about the original ejection velocity field.

preprint2016arXiv

Detection of the YORP Effect for Small Asteroids in the Karin Cluster

The Karin cluster is a young asteroid family thought to have formed only $\simeq 5.75$~My ago. The young age can be demonstrated by numerically integrating the orbits of Karin cluster members backward in time and showing the convergence of the perihelion and nodal longitudes (as well as other orbital elements). Previous work has pointed out that the convergence is not ideal if the backward integration only accounts for the gravitational perturbations from the Solar System planets. It improves when the thermal radiation force known as the Yarkovsky effect it is accounted for. This argument can be used to estimate the spin obliquities of the Karin cluster members. Here we take advantage of the fast growing membership of the Karin cluster and show that the obliquity distribution of diameter $D\simeq 1-2$ km Karin asteroids is bimodal, as expected if the YORP effect acted to move obliquities toward the extreme values ($0^\circ$ or $180^\circ$). The measured magnitude of the effect is consistent with the standard YORP model. The surface thermal conductivity is inferred to be $0.07$-0.2 W m$^{-1}$ K$^{-1}$ (thermal inertia $\simeq 300-500$ J m$^{-2}$ K$^{-1}$s$^{-1/2}$). We find that the strength of the YORP effect is roughly $\simeq 0.7$ of the nominal strength obtained for a collection of random Gaussian spheroids. These results are consistent with a surface composed of rough, rocky regolith. The obliquity values predicted here for 480 members of the Karin cluster can be validated by the light-curve inversion method.

preprint2016arXiv

Dynamical Evidence for a Late Formation of Saturn's Moons

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.

preprint2016arXiv

Hektor - an exceptional D-type family among Jovian Trojans

In this work, we analyze Jovian Trojans in the space of suitable resonant elements and we identify clusters of possible collisional origin by two independent methods: the hierarchical clustering and a so-called "randombox". Compared to our previous work (Brož and Rozehnal 2011), we study a twice larger sample. Apart from Eurybates, Ennomos and $1996\,\rm RJ$ families, we have found three more clusters --- namely families around asteroids (20961)~Arkesilaos, (624)~Hektor in the $L_4$ libration zone and (247341)~$2001\,\rm UV_{209}$ in $L_5$. The families fulfill our stringent criteria, i.e. a high statistical significance, an albedo homogeneity and a steeper size-frequency distribution than that of background. In order to understand their nature, we simulate their long term collisional evolution with the Boulder code (Morbidelli et al. 2009) and dynamical evolution using a modified SWIFT integrator (Levison and Duncan, 1994). Within the framework of our evolutionary model, we were able to constrain the the age of the Hektor family to be either 1 to 4 Gyr or, less likely, 0.1 to 2.5 Gyr, depending on initial impact geometry. Since (624) Hektor itself seems to be a bilobed--shape body with a satellite (Marchis et al. 2014), i.e. an exceptional object, we address its association with the D--type family and we demonstrate that the moon and family could be created during a single impact event. We simulated the cratering event using a Smoothed Particle Hydrodynamics (SPH, Benz and Asphaug, 1994). This is also the first case of a family associated with a D--type parent body.

preprint2016arXiv

Jumping Jupiter can explain Mercury's orbit

The orbit of Mercury has large values of eccentricity and inclination that cannot be easily explained if this planet formed on a circular and coplanar orbit. Here, we study the evolution of Mercury's orbit during the instability related to the migration of the giant planets in the framework of the jumping Jupiter model. We found that some instability models are able to produce the correct values of Mercury's eccentricity and inclination, provided that relativistic effects are included in the precession of Mercury's perihelion. The orbital excitation is driven by the fast change of the normal oscillation modes of the system corresponding to the perihelion precession of Jupiter (for the eccentricity), and the nodal regression of Uranus (for the inclination).

preprint2015arXiv

Dynamical evolution of the Cybele asteroids

The Cybele region, located between the 2J:-1A and 5J:-3A mean-motion resonances, is adjacent and exterior to the asteroid main belt. An increasing density of three-body resonances makes the region between the Cybele and Hilda populations dynamically unstable, so that the Cybele zone could be considered the last outpost of an extended main belt. The presence of binary asteroids with large primaries and small secondaries suggested that asteroid families should be found in this region, but only relatively recently the first dynamical groups were identified in this area. Among these, the Sylvia group has been proposed to be one of the oldest families in the extended main belt. In this work we identify families in the Cybele region in the context of the local dynamics and non-gravitational forces such as the Yarkovsky and stochastic YORP effects. We confirm the detection of the new Helga group at $\simeq$3.65~AU, that could extend the outer boundary of the Cybele region up to the 5J:-3A mean-motion resonance. We obtain age estimates for the four families, Sylvia, Huberta, Ulla and Helga, currently detectable in the Cybele region, using Monte Carlo methods that include the effects of stochastic YORP and variability of the Solar luminosity. The Sylvia family should be $T = 1220 \pm 40$ Myr old, with a possible older secondary solution. Any collisional Cybele group formed prior to the late heavy bombardment would have been most likely completely dispersed in the jumping Jupiter scenario of planetary migration.

preprint2015arXiv

The evolution of asteroids in the jumping-Jupiter migration model

In this work, we investigate the evolution of a primordial belt of asteroids, represented by a large number of massless test particles, under the gravitational effect of migrating Jovian planets in the framework of the jumping-Jupiter model. We perform several simulations considering test particles distributed in the Main Belt, as well as in the Hilda and Trojan groups. The simulations start with Jupiter and Saturn locked in the mutual 3:2 mean motion resonance plus 3 Neptune-mass planets in a compact orbital configuration. Mutual planetary interactions during migration led one of the Neptunes to be ejected in less than 10 Myr of evolution, causing Jupiter to jump by about 0.3 au in semi-major axis. This introduces a large scale instability in the studied populations of small bodies. After the migration phase, the simulations are extended over 4 Gyr, and we compare the final orbital structure of the simulated test particles to the current Main Belt of asteroids with absolute magnitude $H<9.7$. The results indicate that, in order to reproduce the present Main Belt, the primordial belt should have had a distribution peaked at $\sim10^{\circ}$ in inclination and at $\sim0.1$ in eccentricity. We discuss the implications of this for the Grand Tack model. The results also indicate that neither primordial Hildas, nor Trojans, survive the instability, confirming the idea that such populations must have been implanted from other sources. In particular, we address the possibility of implantation of Hildas and Trojans from the Main Belt population, but find that this contribution should be minor.

preprint2014arXiv

Hungaria Asteroid Family as the Source of Aubrite Meteorites

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.

preprint2013arXiv

Early Dynamical Instabilities in the Giant Planet Systems

The observed wide eccentricity distribution of extrasolar giant planets is thought to be the result of dynamical instabilities and gravitational scattering among planets. Previously, it has been assumed that the orbits in giant planet systems become gravitationally unstable after the gas nebula dispersal. It was not well understood, however, how these unstable conditions were established in the first place. In this work we numerically simulate the evolution of systems of three planets as the planets sequentially grow to Jupiter's mass, and dynamically interact among themselves and with the gas disk. We use the hydro-dynamical code FARGO that we modified by implementing the $N$-body integrator SyMBA. The new code can handle close encounters and collisions between planets. To test their stability, the planetary systems were followed with SyMBA for up to $10^8$ yr after the gas disk dispersal. We find that dynamics of the growing planets is complex, because migration and resonances raise their orbital eccentricities, and cause dynamical instabilities when gas is still around. If the dynamical instabilities occur early, planets can be removed by collisions and ejections, and the system rearranges into a new, more stable configuration. In this case, the planetary systems emerging from the gas disks are expected to be stable, and would need to be destabilized by other means (low-mass planets, planetesimal disks, etc.). Alternatively, for the giant planet system to be intrinsically unstable upon the gas disk dispersal, a special timing would be required with the growth of (at least some of) the giant planets having to occur near the end of the gas disk lifetime.

preprint2013arXiv

Olivine-dominated asteroids: Mineralogy and origin

Olivine-dominated asteroids are a rare type of objects formed either in nebular processes or through magmatic differentiation. The analysis of meteorite samples suggest that at least 100 parent bodies in the main belt experienced partial or complete melting and differentiation before being disrupted. However, only a few olivine-dominated asteroids, representative of the mantle of disrupted differentiated bodies, are known to exist. Due to the paucity of these objects in the main belt their origin and evolution have been a matter of great debate over the years. In this work we present a detailed mineralogical analysis of twelve olivine-dominated asteroids. Within our sample we distinguish two classes, one that we call monomineralic-olivine asteroids and another referred to as olivine-rich asteroids. For the monomineralic-olivine asteroids the olivine chemistry was found to range from ~ Fo49 to Fo70, consistent with the values measured for brachinites and R chondrites. In the case of the olivine-rich asteroids we determined their olivine and low-Ca pyroxene abundance using a new set of spectral calibrations derived from the analysis of R chondrites spectra. We found that the olivine abundance for these asteroids varies from 0.68 to 0.93, while the fraction of low-Ca pyroxene to total pyroxene ranges from 0.6 to 0.9. A search for dynamical connections between the olivine-dominated asteroids and asteroid families found no genetic link (of the type core-mantel-crust) between these objects.

preprint2013arXiv

On the Dynamics and Origin of Haumea's Moons

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.

preprint2013arXiv

Titan-Hyperion Resonance and the Tidal Q of Saturn

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.