Engineering topics

Nesvorný, David

Publications and source records attributed to Nesvorný, David.

17 records

DOE OSTI2022

Contribution of Ryugu-like material to Earth’s volatile inventory by Cu and Zn isotopic analysis

We report initial analyses showed that asteroid Ryugu’s composition is close to CI (Ivuna-like) carbonaceous chondrites (CCs) – the chemically most primitive meteorites, characterized by near-solar abundances for most elements. However, some isotopic signatures (for example, Ti, Cr) overlap with other CC groups, so the details of the link between Ryugu and the CI chondrites are not yet fully clear. Here we show that Ryugu and CI chondrites have the same zinc and copper isotopic composition. As the various chondrite groups have very distinct Zn and Cu isotopic signatures, our results point at a common genetic heritage between Ryugu and CI chondrites, ruling out any affinity with other CC groups. Since Ryugu’s pristine samples match the solar elemental composition for many elements, their Zn and Cu isotopic compositions likely represent the best estimates of the solar composition. Earth’s mass-independent Zn isotopic composition is intermediate between Ryugu/CC and non-carbonaceous chondrites (NCs), suggesting a contribution of Ryugu-like material to Earth’s budgets of Zn and other moderately volatile elements.

58 GEOSCIENCES↗

DOE OSTI2021

The Role of Early Giant-planet Instability in Terrestrial Planet Formation

The terrestrial planets are believed to have formed by violent collisions of tens of lunar- to Mars-size protoplanets at time t < 200 Myr after the protoplanetary gas disk dispersal (t {sub 0}). The solar system giant planets rapidly formed during the protoplanetary disk stage and, after t {sub 0}, radially migrated by interacting with outer disk planetesimals. An early (t < 100 Myr) dynamical instability is thought to have occurred with Jupiter having gravitational encounters with a planetary-size body, jumping inward by ∼0.2–0.5 au, and landing on its current, mildly eccentric orbit. Here we investigate how the giant-planet instability affected the formation of the terrestrial planets. We study several instability cases that were previously shown to match many solar system constraints. We find that resonances with giant planets help to remove solids available for accretion near ∼1.5 au, thus stalling the growth of Mars. It does not matter, however, whether the giant planets are placed on their current orbits at t {sub 0} or whether they realistically evolve in one of our instability models; the results are practically the same. The tight orbital spacing of Venus and Earth is difficult to reproduce in our simulations, including cases where bodies grow from a narrow annulus at 0.7–1 au, because protoplanets tend to spread radially during accretion. The best results are obtained in the narrow-annulus model when protoplanets emerging from the dispersing gas nebula are assumed to have (at least) the mass of Mars. This suggests efficient accretion of the terrestrial protoplanets during the first ∼10 Myr of the solar system.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Clarissa Family Age from the Yarkovsky Effect Chronology

The Clarissa family is a small collisional family composed of primitive C-type asteroids. It is located in a dynamically stable zone of the inner asteroid belt. In this work we determine the formation age of the Clarissa family by modeling planetary perturbations as well as thermal drift of family members due to the Yarkovsky effect. Simulations were carried out using the SWIFT-RMVS4 integrator modified to account for the Yarkovsky and Yarkovsky–O’Keefe–Radzievskii–Paddack (YORP) effects. We ran multiple simulations starting with different ejection velocity fields of fragments, varying proportion of initially retrograde spins, and also tested different Yarkovsky/YORP models. Our goal was to match the observed orbital structure of the Clarissa family which is notably asymmetrical in the proper semimajor axis, a {sub p}. The best fits were obtained with the initial ejection velocities ≲20 m s{sup −1} of diameter D ≃ 2 km fragments, ∼4:1 preference for spin-up by YORP, and assuming that ≃80% of small family members initially had retrograde rotation. The age of the Clarissa family was found to be t {sub age} = 56 ± 6 Myr for the assumed asteroid density ρ = 1.5 g cm{sup −3}. Small variation of density to smaller or larger value would lead to slightly younger or older age estimates. This is the first case where the Yarkovsky effect chronology has been successfully applied to an asteroid family younger than 100 Myr.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Bi-lobed Shape of Comet 67P from a Collapsed Binary

The Rosetta spacecraft observations revealed that the nucleus of comet 67P/Churyumov–Gerasimenko consists of two similarly sized lobes connected by a narrow neck. Here, we evaluate the possibility that 67P is a collapsed binary. We assume that the progenitor of 67P was a binary and consider various physical mechanisms that could have brought the binary components together, including small-scale impacts and gravitational encounters with planets. We find that 67P could be a primordial body (i.e., not a collisional fragment) if the outer planetesimal disk lasted ≲10 Myr before it was dispersed by migrating Neptune. The probability of binary collapse by impact is ≃30% for tightly bound binaries. Most km-class binaries become collisionally dissolved. Roughly 10% of the surviving binaries later evolve to become contact binaries during the disk dispersal, when bodies suffer gravitational encounters with Neptune. Overall, the processes described in this work do not seem to be efficient enough to explain the large fraction (∼67%) of bi-lobed cometary nuclei inferred from spacecraft imaging.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Origin and Evolution of Short-period Comets

Comets are icy objects that orbitally evolve from the trans-Neptunian region into the inner solar system, where they are heated by solar radiation and become active due to the sublimation of water ice. Here we perform simulations in which cometary reservoirs are formed in the early solar system and evolved over 4.5 Gyr. The gravitational effects of Planet 9 (P9) are included in some simulations. Different models are considered for comets to be active, including a simple assumption that comets remain active for N{sub p}(q) perihelion passages with perihelion distance q<2.5 au. The orbital distribution and number of active comets produced in our model is compared to observations. The orbital distribution of ecliptic comets (ECs) is well reproduced in models with N{sub p}(2.5)≃500 and without P9. With P9, the inclination distribution of model ECs is wider than the observed one. We find that the known Halley-type comets (HTCs) have a nearly isotropic inclination distribution. The HTCs appear to be an extension of the population of returning Oort-cloud comets (OCCs) to shorter orbital periods. The inclination distribution of model HTCs becomes broader with increasing N{sub p}, but the existing data are not good enough to constrain N{sub p} from orbital fits. N{sub p}(2.5)>1000 is required to obtain a steady-state population of large active HTCs that is consistent with observations. To fit the ratio of the returning-to-new OCCs, by contrast, our model implies that N{sub p}(2.5)≲10, possibly because the detected long-period comets are smaller and much easier to disrupt than observed HTCs.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Neptune's orbital migration was grainy, not smooth

The Kuiper Belt is a population of icy bodies beyond the orbit of Neptune. The complex orbital structure of the Kuiper Belt, including several categories of objects inside and outside of resonances with Neptune, emerged as a result of Neptune’s migration into an outer planetesimal disk. An outstanding problem with the existing migration models is that they invariably predict excessively large resonant populations, while observations show that the non-resonant orbits are in fact common (e.g., the main belt population is ≃2–4 times larger than Plutinos in the 3:2 resonance). Here we show that this problem can be resolved if it is assumed that Neptune’s migration was grainy, as expected from scattering encounters of Neptune with massive planetesimals. The grainy migration acts to destabilize resonant bodies with large libration amplitudes, a fraction of which ends up on stable non-resonant orbits. Thus, the non-resonant-to-resonant ratio obtained with the grainy migration is higher, up to ∼10 times higher for the range of parameters investigated here, than in a model with smooth migration. In addition, the grainy migration leads to a narrower distribution of the libration amplitudes in the 3:2 resonance. The best fit to observations is obtained when it is assumed that the outer planetesimal disk below 30 au contained 1000–4000 Plutos. We estimate that the combined mass of Pluto-class objects in the original disk represented 10%–40% of the estimated disk mass (M{sub disk}≃20 M{sub Earth}). This constraint can be used to better understand the accretion processes in the outer solar system.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

DYNAMICS AND TRANSIT VARIATIONS OF RESONANT EXOPLANETS

Transit timing variations (TTVs) are deviations of the measured midtransit times from the exact periodicity. One of the most interesting causes of TTVs is the gravitational interaction between planets. Here we consider a case of two planets in a mean motion resonance (orbital periods in a ratio of small integers). This case is important because the resonant interaction can amplify the TTV effect and allow planets to be detected more easily. We develop an analytic model of the resonant dynamics valid for small orbital eccentricities and use it to derive the principal TTV terms. We find that a resonant system should show TTV terms with two basic periods (and their harmonics). The resonant TTV period is proportional ( m / M {sub *}){sup −2/3}, where m and M {sub *} are the planetary and stellar masses. For m = 10{sup −4} M {sub *}, for example, the TTV period exceeds the orbital period by about two orders of magnitude. The amplitude of the resonant TTV terms scales linearly with the libration amplitude. The ratio of the TTV amplitudes of two resonant planets is inversely proportional to the ratio of their masses. These and other relationships discussed in the main text can be used to aid the interpretation of TTV observations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

THE SCHULHOF FAMILY: SOLVING THE AGE PUZZLE

The Schulhof family, a tight cluster of small asteroids around the central main belt body (2384) Schulhof, belongs to a so far rare class of very young families (estimated ages less than 1 Myr). Characterization of these asteroid clusters may provide important insights into the physics of the catastrophic disruption of their parent body. The case of the Schulhof family has been up to now complicated by the existence of two proposed epochs of its origin. In this paper, we first use our own photometric observations, as well as archival data, to determine the rotation rate and spin axis orientation of the largest fragment (2384) Schulhof. Our data also allow us to better constrain the absolute magnitude of this asteroid, and thus also improve the determination of its geometric albedo. Next, using the up-to-date catalog of asteroid orbits, we perform a new search of smaller members in the Schulhof family, increasing their number by 50%. Finally, the available data are used to access Schulhof's family age anew. We now find that the younger of the previously proposed two ages of this family is not correct, resulting from a large orbital uncertainty of single-opposition members. Our new runs reveal a single age solution of about 800 kyr with a realistic uncertainty of 200 kyr.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

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 three 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 semimajor 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 ∼10° in inclination and at ∼0.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.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2015

TILTING JUPITER (A BIT) AND SATURN (A LOT) DURING PLANETARY MIGRATION

We study the effects of planetary late migration on the gas giants’ obliquities. We consider the planetary instability models from Nesvorný and Morbidelli, in which the obliquities of Jupiter and Saturn can be excited when spin–orbit resonances occur. The most notable resonances occur when the s{sub 7} and s{sub 8} frequencies, changing as a result of planetary migration, become commensurate with the precession frequencies of Jupiter’s and Saturn’s spin vectors. We show that Jupiter may have obtained its present obliquity by crossing of the s{sub 8} resonance. This would set strict constraints on the character of migration during the early stage. Additional effects on Jupiter’s obliquity are expected during the last gasp of migration when the s{sub 7} resonance was approached. The magnitude of these effects depends on the precise value of the Jupiter’s precession constant. Saturn’s large obliquity was likely excited by capture into the s{sub 8} resonance. This probably happened during the late stage of planetary migration when the evolution of the s{sub 8} frequency was very slow, and the conditions for capture into the spin–orbit resonance with s{sub 8} were satisfied. However, whether or not Saturn is in the spin–orbit resonance with s{sub 8} at the present time is not clear because the existing observations of Saturn’s spin precession and internal structure models have significant uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Excitation of the orbital inclination of Iapetus during planetary encounters

Saturn's moon, Iapetus, has an orbit in a transition region where the Laplace surface is bending from the equator to the orbital plane of Saturn. The orbital inclination of Iapetus to the local Laplace plane is ≅ 8°, which is unexpected because the inclination should be ≅ 0 if Iapetus formed from a circumplanetary disk on the Laplace surface. It thus appears that some process has pumped up Iapetus's inclination while leaving its eccentricity near zero (e ≅ 0.03 at present). Here, we examined the possibility that Iapetus's inclination was excited during the early solar system instability when encounters between Saturn and ice giants occurred. We found that the dynamical effects of planetary encounters on Iapetus's orbit sensitively depend on the distance of the few closest encounters. In 4 out of 10 instability cases studied here, the orbital perturbations were too large to be plausible. In one case, Iapetus's orbit was practically unaffected. In the remaining five cases, the perturbations of Iapetus's inclination were adequate to explain its present value. In three of these cases, however, Iapetus's eccentricity was excited to >0.1-0.25, and it is not clear whether it could have been damped to its present value (≅ 0.03) by a subsequent process (e.g., tides and dynamical friction from captured irregular satellites do not seem to be strong enough). Our results therefore imply that only 2 out of 10 instability cases (∼20%) can excite Iapetus's inclination to its current value (∼30% of trials lead to >5°) while leaving its orbital eccentricity low.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Transit timing variations for planets co-orbiting in the horseshoe regime

Although not yet detected, pairs of exoplanets in 1:1 mean motion resonance probably exist. Low eccentricity, near-planar orbits, which in the comoving frame follow horseshoe trajectories, are one of the possible stable configurations. Here we study transit timing variations (TTVs) produced by mutual gravitational interaction of planets in this orbital architecture, with the goal to develop methods that can be used to recognize this case in observational data. In particular, we use a semi-analytic model to derive parametric constraints that should facilitate data analysis. We show that characteristic traits of the TTVs can directly constrain the (1) ratio of planetary masses and (2) their total mass (divided by that of the central star) as a function of the minimum angular separation as seen from the star. In an ideal case, when transits of both planets are observed and well characterized, the minimum angular separation can also be inferred from the data. As a result, parameters derived from the observed transit timing series alone can directly provide both planetary masses scaled to the central star mass.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Orbital perturbations of the Galilean satellites during planetary encounters

The Nice model of the dynamical instability and migration of the giant planets can explain many properties of the present solar system, and can be used to constrain its early architecture. In the jumping-Jupiter version of the Nice model, required from the terrestrial planet constraint and dynamical structure of the asteroid belt, Jupiter has encounters with an ice giant. Here, we study the survival of the Galilean satellites in the jumping-Jupiter model. This is an important concern because the ice-giant encounters, if deep enough, could dynamically perturb the orbits of the Galilean satellites and lead to implausible results. We performed numerical integrations where we tracked the effect of planetary encounters on the Galilean moons. We considered three instability cases from Nesvorný and Morbidelli that differed in the number and distribution of encounters. We found that in one case, where the number of close encounters was relatively small, the Galilean satellite orbits were not significantly affected. In the other two, the orbital eccentricities of all moons were excited by encounters, Callisto's semimajor axis changed, and, in a large fraction of trials, the Laplace resonance of the inner three moons was disrupted. The subsequent evolution by tides damps eccentricities and can recapture the moons in the Laplace resonance. A more important constraint is represented by the orbital inclinations of the moons, which can be excited during the encounters and not appreciably damped by tides. We find that one instability case taken from Nesvorný and Morbidelli clearly does not meet this constraint. This shows how the regular satellites of Jupiter can be used to set limits on the properties of encounters in the jumping-Jupiter model, and help us to better understand how the early solar system evolved.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Photo-dynamical analysis of three Kepler objects of interest with significant transit timing variations

KOI-227, KOI-319 and KOI-884 are identified here as (at least) two planet systems. For KOI-319 and KOI-884, the observed Transit Timing Variations (TTVs) of the inner transiting planet are used to detect an outer non-transiting planet. The outer planet in KOI-884 is ≅2.6 Jupiter masses and has the orbital period just narrow of the 3:1 resonance with the inner planet (orbital period ratio 2.93). The distribution of parameters inferred from KOI-319.01's TTVs is bimodal with either a ≅1.6 Neptune-mass (M{sub N}) planet wide of the 5:3 resonance (period 80.1 days) or a ≅1 Saturn-mass planet wide of the 7:3 resonance (period 109.2 days). The radial velocity measurements can be used in this case to determine which of these parameter modes is correct. KOI-227.01's TTVs with large ≅10 hr amplitude can be obtained for planetary-mass companions in various major resonances. Based on the Bayesian evidence, the current TTV data favor the outer 2:1 resonance with a companion mass ≅1.5 M{sub N}, but this solution implies a very large density of KOI-227.01. The inner and outer 3:2 resonance solutions with sub-Neptune-mass companions are physically more plausible, but will need to be verified.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

The effect of conjunctions on the transit timing variations of exoplanets

We develop an analytic model for transit timing variations produced by orbital conjunctions between gravitationally interacting planets. If the planetary orbits have tight orbital spacing, which is a common case among the Kepler planets, the effect of a single conjunction can be best described as: (1) a step-like change of the transit timing ephemeris with subsequent transits of the inner planet being delayed and those of the outer planet being sped up, and (2) a discrete change in sampling of the underlying oscillations from eccentricity-related interaction terms. In the limit of small orbital eccentricities, our analytic model gives explicit equations for these effects as a function of the mass and orbital separation of planets. We point out that a detection of the conjunction effect in real data is of crucial importance for the physical characterization of planetary systems from transit timing variations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Dynamical model for the toroidal sporadic meteors

More than a decade of radar operations by the Canadian Meteor Orbit Radar have allowed both young and moderately old streams to be distinguished from the dispersed sporadic background component. The latter has been categorized according to broad radiant regions visible to Earth-based observers into three broad classes: the helion and anti-helion source, the north and south apex sources, and the north and south toroidal sources (and a related arc structure). The first two are populated mainly by dust released from Jupiter-family comets and new comets. Proper modeling of the toroidal sources has not to date been accomplished. Here, we develop a steady-state model for the toroidal source of the sporadic meteoroid complex, compare our model with the available radar measurements, and investigate a contribution of dust particles from our model to the whole population of sporadic meteoroids. We find that the long-term stable part of the toroidal particles is mainly fed by dust released by Halley type (long period) comets (HTCs). Our synthetic model reproduces most of the observed features of the toroidal particles, including the most troublesome low-eccentricity component, which is due to a combination of two effects: particles' ability to decouple from Jupiter and circularize by the Poynting-Robertson effect, and large collision probability for orbits similar to that of the Earth. Our calibrated model also allows us to estimate the total mass of the HTC-released dust in space and check the flux necessary to maintain the cloud in a steady state.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

TTVFast: An efficient and accurate code for transit timing inversion problems

Transit timing variations (TTVs) have proven to be a powerful technique for confirming Kepler planet candidates, for detecting non-transiting planets, and for constraining the masses and orbital elements of multi-planet systems. These TTV applications often require the numerical integration of orbits for computation of transit times (as well as impact parameters and durations); frequently tens of millions to billions of simulations are required when running statistical analyses of the planetary system properties. We have created a fast code for transit timing computation, TTVFast, which uses a symplectic integrator with a Keplerian interpolator for the calculation of transit times. The speed comes at the expense of accuracy in the calculated times, but the accuracy lost is largely unnecessary, as transit times do not need to be calculated to accuracies significantly smaller than the measurement uncertainties on the times. The time step can be tuned to give sufficient precision for any particular system. We find a speed-up of at least an order of magnitude relative to dynamical integrations with high precision using a Bulirsch-Stoer integrator.

79 ASTRONOMY AND ASTROPHYSICS↗