Source author record

Matija Ćuk

Matija Ćuk 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

15works
4topics
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

15 published item(s)

preprint2022arXiv

Cupid Is Not Doomed Yet: On the Stability of the Inner Moons of Uranus

Some of the small inner moons of Uranus have very closely-spaced orbits. Multiple numerical studies have found that the moons Cressida and Desdemona, within the Portia sub-group, are likely to collide in less than 100 Myr. The subsequent discovery of three new moons (Cupid, Perdita, and Mab) made the system even more crowded. In particular, it has been suggested that the Belinda group (Cupid, Belinda, and Perdita) will become unstable in as little as 10$^5$ years. Here we revisit the issue of the stability of the inner moons of Uranus using updated orbital elements and considering tidal dissipation. We find that the Belinda group can be stable on $10^8$-year timescales due to an orbital resonance between Belinda and Perdita. We find that tidal evolution cannot form the Belinda-Perdita resonance, but convergent migration could contribute to the long-term instability of the Portia group. We propose that Belinda captured Perdita into the resonance during the last episode of disruption and re-accretion among the inner moons, possibly hundreds of Myr ago.

preprint2022arXiv

Three-Body Resonances in the Saturnian System

Saturn has a dynamically rich satellite system, which includes at least three orbital resonances between three pairs of moons: Mimas-Tethys 4:2, Enceladus-Dione 2:1, and Titan-Hyperion 4:3 mean-motion resonances. Studies of the orbital history of Saturn's moons usually assume that their past dynamics was also dominated solely by two-body resonances. Using direct numerical integrations, we find that three-body resonances among Saturnian satellites were quite common in the past, and could result in a relatively long-term, but finite capture time (10 Myr or longer). We find that these three-body resonances are invariably of the eccentricity type, and do not appear to affect the moons' inclinations. While some three-body resonances are located close to two-body resonances (but involve the orbital precession of the third body), others are isolated, with no two-body arguments being near resonance. We conclude that future studies of the system's past must take full account of three-body resonances, which have been overlooked in the past work.

preprint2021arXiv

Origin of the Moon

The Earth-Moon system is unusual in several respects. The Moon is roughly 1/4 the radius of the Earth - a larger satellite-to-planet size ratio than all known satellites other than Pluto's Charon. The Moon has a tiny core, perhaps with only ~1% of its mass, in contrast to Earth whose core contains nearly 30% of its mass. The Earth-Moon system has a high total angular momentum, implying a rapidly spinning Earth when the Moon formed. In addition, the early Moon was hot and at least partially molten with a deep magma ocean. Identification of a model for lunar origin that can satisfactorily explain all of these features has been the focus of decades of research.

preprint2020arXiv

Dynamical History of the Uranian System

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.

preprint2020arXiv

Evidence for a Past Martian Ring from the Orbital Inclination of Deimos

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.

preprint2020arXiv

Pathways to Sustainable Planetary Science

Climate change is a major impending threat to the future of humanity. According to the International Panel on Climate Change (IPCC), our emissions are estimated to have caused 0.8 deg C-1.2 deg C of anthropogenic global warming (AGW) above pre-industrial levels. AGW is likely to reach 1.5 degrees C between 2030 and 2052 if it continues to increase at the current rate. As the climate change is driven by the release of carbon dioxide and other greenhouse gases (GHG) into the atmosphere, there is a broad consensus that the mitigation of climate change requires transition to low GHG emission energy sources, technologies and practices. Implementing such changes systematically from individual to community-wide scales together with the resulting cultural changes and leadership towards environmental consciousness and responsibility are crucial to mitigate the looming damage of AGW. Given planetary scientists' wide recognition of the realities of climate change, and the need for us to maintain credibility by leading by example, it is appropriate to make own professional behavior more environmentally responsible. While scientists are few in numbers, and planetary scientists far fewer, high volumes of academic travel to conferences, panels, colloquia, and research collaboration visits together with extensive use of large, energetically demanding infrastructures make the "carbon footprint" of scientists much higher than that of an average citizen. This White Paper focuses on how modifying our activities, particularly associated with academic travel, can affect the carbon footprint of the planetary science community, and it makes recommendations on how the community and the funding agencies could best participate in the cultural change required to mitigate the damage that AGW will cause.

preprint2019arXiv

Early Dynamics of the Lunar Core

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.

preprint2019arXiv

The energy budget and figure of Earth during recovery from the Moon-forming giant impact

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.

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.

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.

preprint2014arXiv

Yarkovsky-Driven Spreading of the Eureka Family of Mars Trojans

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.

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.

preprint2012arXiv

Long-Term Stability of Horseshoe Orbits

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.

preprint2011arXiv

Chronology and Sources of Lunar Impact Bombardment

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]