BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

David Nesvorny

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

ResearcherUnclaimed source record
astro-ph.EPastro-ph.SRastro-ph.IMphysics.geo-phphysics.space-ph

Catalog footprint

What is connected

44works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

44 published item(s)

preprint2022arXiv

Asteroid Families: properties, recent advances and future opportunities

Collisions are one of the key processes shaping planetary systems. Asteroid families are outcomes of such collisions still identifiable across our solar system. The families provide a unique view of catastrophic disruption phenomena and have been in the focus of planetary scientists for more than a century. Most of them are located in the main belt, a ring of asteroids between Mars and Jupiter. Here we review the basic properties of the families, discuss some recent advances, and anticipate future challenges. This review pays more attention to dynamic aspects such as family identification, age determination, and long-term evolution. The text, however, goes beyond that. Especially, we cover the details of young families that see the major advances in the last years, and we anticipate it will develop even faster in the future. We also discuss the relevance of asteroid families for water-ice content in the asteroid belt and our current knowledge on links between families and main-belt comets.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Dynamical Implantation of Blue Binaries in the Cold Classical Kuiper Belt

Colors and binarity provide important constraints on the Kuiper belt formation. The cold classical objects at radial distance r=42-47 au from the Sun are predominantly very red (spectral slope s>17%) and often exist as equal-size binaries (~30% observed binary fraction). This has been taken as evidence for the in-situ formation of cold classicals. Interestingly, a small fraction (~10%) of cold classicals is less red with s<17%, and these "blue" bodies are often found in wide binaries. Here we study the dynamical implantation of blue binaries from r<42 au. We find that they can be implanted into the cold classical belt from a wide range of initial radial distances, but the survival of the widest blue binaries -- 2001 QW322 and 2003 UN284 -- implies formation at r>30 au. This would be consistent with the hypothesized less-red to very-red transition at 30<r<40 au. For any reasonable choice of parameters (Neptune's migration history, initial disk profile, etc.), however, our model predicts a predominance of blue singles, rather than blue binaries, which contradicts existing observations. We suggest that wide blue binaries formed in situ at r=42-47 au and their color reflects early formation in a protoplanetary gas disk. The predominantly VR colors of cold classicals may be related to the production of methanol and other hydrocarbons during the late stages of the disk, when the temperature at 45 au dropped to 20 K and carbon monoxide was hydrogenated.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Origin and dynamical evolution of the asteroid belt

The asteroid belt was dynamically shaped during and after planet formation. Despite representing a broad ring of stable orbits, the belt contains less than one one-thousandth of an Earth mass. The asteroid orbits are dynamically excited with a wide range in eccentricity and inclination and their compositions are diverse, with a general trend toward dry objects in the inner belt and more water-rich objects in the outer belt. Here we review models of the asteroid belt's origins and dynamical history. The classical view is that the belt was born with several Earth masses in planetesimals, then strongly depleted. However, it is possible that very few planetesimals ever formed in the asteroid region and that the belt's story is one of implantation rather than depletion. A number of processes may have implanted asteroids from different regions of the Solar System, dynamically removed them, and excited their orbits. During the gaseous disk phase these include the effects of giant planet growth and migration and sweeping secular resonances. After the gaseous disk phase these include scattering from resident planetary embryos, chaos in the giant planets' orbits, the giant planet instability, and long-term dynamical evolution. Different global models for Solar System formation imply contrasting dynamical histories of the asteroid belt. Vesta and Ceres may have been implanted from opposite regions of the Solar System -- Ceres from the Jupiter-Saturn region and Vesta from the terrestrial planet region -- and could therefore represent very different formation conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

TOI-216: Resonant Constraints on Planet Migration

TOI-216 is a pair of close-in planets with orbits deep in the 2:1 mean motion resonance. The inner, Neptune-class planet (TOI-216b) is near 0.12 au (orbital period $P_{\rm b} \simeq 17$ d) and has a substantial orbital eccentricity ($e_{\rm b} \simeq 0.16$), and large libration amplitude ($A_ψ\simeq 60^\circ$) in the resonance. The outer planet (TOI-216c) is a gas giant on a nearly circular orbit. We carry out $N$-body simulations of planet migration in a protoplanetary gas disk to explain the orbital configuration of TOI-216 planets. We find that TOI-216b's migration must have been halted near its current orbital radius to allow for a convergent migration of the two planets into the resonance. For the inferred damping-to-migration timescale ratio $τ_e/τ_a \simeq 0.02$, overstable librations in the resonance lead to a limit cycle with $A_ψ\simeq 80^\circ$ and $e_{\rm b}<0.1$. The system could have remained in this configuration for the greater part of the protoplanetary disk lifetime. If the gas disk was removed from inside out, this would have reduced the libration amplitude to $A_ψ\simeq 60^\circ$ and boosted $e_{\rm b}$ via the resonant interaction with TOI-216c. Our results suggest a relatively fast inner disk removal ($\sim 10^5$ yr). Another means of explaining the large libration amplitude is stochastic stirring from a (turbulent) gas disk. For that to work, overstable librations would need to be suppressed, $τ_e/τ_a \simeq 0.05$, and very strong turbulent stirring (or some other source of large stochastic forcing) would need to overcome the damping effects of gas. Hydrodynamical simulations can be performed to test these models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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. The best fits were obtained with the initial ejection velocities < ~20 m/s 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 56+/-6 Myr for the assumed asteroid density 1.5 g/cm3. 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Eccentric Early Migration of Neptune

The dynamical structure of the Kuiper belt can be used as a clue to the formation and evolution of the Solar System, planetary systems in general, and Neptune's early orbital history in particular. The problem is best addressed by forward modeling where different initial conditions and Neptune's orbital evolutions are tested, and the model predictions are compared to orbits of known Kuiper belt objects (KBOs). It has previously been established that Neptune radially migrated, by gravitationally interacting with an outer disk of planetesimals, from the original radial distance r < 25 au to its current orbit at 30 au. Here we show that the migration models with a very low orbital eccentricity of Neptune e_N < 0.05 do not explain KBOs with semimajor axes 50 < a < 60 au, perihelion distances q > 35 au and inclinations i < 10 deg. If e_N < 0.05 at all times, the Kozai cycles control the implantation process and the orbits with q > 35 au end up having, due to the angular momentum's z-component conservation, i > 10 deg. Better results are obtained when Neptune's eccentricity is excited to e_N ~ 0.1 and subsequently damped by dynamical friction. The low-e & low-i orbits at 50-60 au are produced in this model when KBOs are lifted from the scattered disk by secular cycles -- mainly the apsidal resonance nu_8 -- near various mean motion resonances. These results give support to a (mild) dynamical instability that presumably excited the orbits of giant planets during Neptune's early migration.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

OSSOS XX: The Meaning of Kuiper Belt Colors

Observations show that 100-km-class Kuiper belt objects (KBOs) can be divided in (at least) two color groups, hereafter red (R, g-i<1.2) and very red (VR, g-i>1.2), reflecting a difference in their surface composition. This is thought to imply that KBOs formed over a relatively wide range of radial distance, r. The cold classicals at 42<r<47 au are predominantly VR and known Neptune Trojans at r=30 au are mostly R. Intriguingly, however, the dynamically hot KBOs show a mix of R and VR colors and no correlation of color with r. Here we perform migration/instability simulations where the Kuiper belt is populated from an extended planetesimal disk. We find that the color observations can be best understood if R objects formed at r<r* and VR objects at r>r*, with 30<r*<40 au. The proposed transition at 30<r*<40 au would explain why the VR objects in the dynamically hot population have smaller orbital inclinations than the R objects, because the orbital excitation from Neptune weakens for orbits starting beyond 30 au. Possible causes of the R-VR color bimodality are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Stability of Jovian Trojans and their collisional families

The Jovian Trojans are two swarms of objects located around the L$_4$ and L$_5$ Lagrange points. The population is thought to have been captured by Jupiter during the Solar system's youth. Within the swarms, six collisional families have been identified in previous work, with four in the L$_4$ swarm, and two in the L$_5$. Our aim is to investigate the stability of the two Trojan swarms, with a particular focus on these collisional families. We find that the members of Trojan swarms escape the population at a linear rate, with the primordial L$_4$ (23.35% escape) and L$_5$ (24.89% escape) population sizes likely 1.31 and 1.35 times larger than today. Given that the escape rates were approximately equal between the two Trojan swarms, our results do not explain the observed asymmetry between the two groups, suggesting that the numerical differences are primordial in nature, supporting previous studies. Upon leaving the Trojan population, the escaped objects move onto orbits that resemble those of the Centaur and short-period comet populations. Within the Trojan collisional families, the 1996 RJ and 2001 UV$_{209}$ families are found to be dynamically stable over the lifetime of the Solar system, whilst the Hektor, Arkesilos and Ennomos families exhibit various degrees of instability. The larger Eurybates family shows 18.81\% of simulated members escaping the Trojan population. Unlike the L4 swarm, the escape rate from the Eurybates family is found to increase as a function of time, allowing an age estimation of approximately $1.045\pm 0.364 \times 10^9$ years.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Superparticle Method for Simulating Collisions

For problems in astrophysics, planetary science and beyond, numerical simulations are often limited to simulating fewer particles than in the real system. To model collisions, the simulated particles (aka superparticles) need to be inflated to represent a collectively large collisional cross section of real particles. Here we develop a superparticle-based method that replicates the kinetic energy loss during real-world collisions, implement it in an $N$-body code and test it. The tests provide interesting insights into dynamics of self gravitating collisional systems. They show how particle systems evolve over several free fall timescales to form central concentrations and equilibrated outer shells. The superparticle method can be extended to account for the accretional growth of objects during inelastic mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Role of Early Giant Planet Instability in the 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_0). The solar system giant planets rapidly formed during the protoplanetary disk stage and, after t_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 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 the 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_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 radially spread 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 Mars mass. This suggests efficient accretion of the terrestrial protoplanets during the first ~10 Myr of the solar system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Trans-Neptunian binaries (2018)

The discovery and characteristics of transneptunian binaries are reviewed. In the 20 years since their first discovery, a wealth of detail has emerged including the frequency of binaries in different populations, their relative sizes and separations, and colors. Taken globally, these properties give strong clues to the origin and evolution of the populations where these binaries are found. In the last ten years and increasing number of binary orbits have been determined which yields a new trove of information on their masses and densities as well as details of their orbits including inclination, eccentricity and the timing of mutual events. In 2018, the study of transneptunian binaries remains one of the most active areas of progress in understanding the Solar System beyond Neptune.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Kuiper belt: formation and evolution

This chapter reviews accretion models for Kuiper belt objects (KBOs), discussing in particular the compatibility of the observed properties of the KBO population with the streaming instability paradigm. Then it discusses how the dynamical structure of the KBO population, including the formation of its 5 sub-components (cold, hot, resonant, scattered and fossilized), can be quantitatively understood in the framewok of the giant planet instability. We also establish the connections between the KBO population and the Trojans of Jupiter and Neptune, the irregular satellites of all giant planets, the Oort cloud and the D-type main belt asteroids. Finally, we discuss the collisional evolution of the KBO population, arguing that the current size-frequency distribution below 100~km in size has been achieved as a collisional equilibrium in a few tens of My inside the original massive trans-Neptunain disk, possibly with the exception of the cold population sub-component.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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 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. 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, or ~1000 bodies twice as massive as Pluto. We find that the probability for an outer disk object to end up on a stable orbit in the Kuiper belt is ~10^{-3}. Together, these results are consistent with having only two (known) Pluto-mass bodies in the Kuiper belt today (Pluto and Eris). The combined mass of Pluto-class objects in the original disk was ~2-8 Earth masses (M_Earth), which represents 10-40% of the estimated disk mass (M_disk ~ 20 M_Earth).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Asteroid Belt as a Relic From a Chaotic Early Solar System

The orbital structure of the asteroid belt holds a record of the Solar System's dynamical history. The current belt only contains ${\rm \sim 10^{-3}}$ Earth masses yet the asteroids' orbits are dynamically excited, with a large spread in eccentricity and inclination. In the context of models of terrestrial planet formation, the belt may have been excited by Jupiter's orbital migration. The terrestrial planets can also be reproduced without invoking a migrating Jupiter; however, as it requires a severe mass deficit beyond Earth's orbit, this model systematically under-excites the asteroid belt. Here we show that the orbits of the asteroids may have been excited to their current state if Jupiter and Saturn's early orbits were chaotic. Stochastic variations in the gas giants' orbits cause resonances to continually jump across the main belt and excite the asteroids' orbits on a timescale of tens of millions of years. While hydrodynamical simulations show that the gas giants were likely in mean motion resonance at the end of the gaseous disk phase, small perturbations could have driven them into a chaotic but stable state. The gas giants' current orbits were achieved later, during an instability in the outer Solar System. Although it is well known that the present-day Solar System exhibits chaotic behavior, our results suggest that the early Solar System may also have been chaotic.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The orbital distribution of trans-Neptunian objects beyond 50 au

The dynamical structure of the Kuiper belt beyond 50 au is not well understood. Here we report results of a numerical model with long-range, slow and grainy migration of Neptune. The model implies that bodies scattered outward by Neptune to semimajor axes a>50 au often evolve into resonances which subsequently act to raise the perihelion distances of orbits to q>40 au. The implication of the model is that the orbits with 50<a<100 au and q>40 au should cluster near (but not in) the resonances with Neptune (3:1 at a=62.6 au, 4:1 at a=75.9 au, 5:1 at a=88.0 au, etc.). The recent detection of several distant Kuiper Belt Objects (KBOs) near resonances is consistent with this prediction, but it is not yet clear whether the orbits are really non-resonant as our model predicts. We estimate from the model that there should presently be ~1600-2400 bodies at the 3:1 resonance and ~1000-1400 bodies at the 4:1 resonance (for q>40 au and diameters D>100 km). These results favorably compare with the population census of distant KBOs inferred from existing observations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Jumping Neptune Can Explain the Kuiper Belt Kernel

The Kuiper belt is a population of icy bodies beyond the orbit of Neptune. A particularly puzzling and up-to-now unexplained feature of the Kuiper belt is the so-called `kernel', a concentration of orbits with semimajor axes a~44 AU, eccentricities e~0.05, and inclinations i<5 deg. Here we show that the Kuiper belt kernel can be explained if Neptune's otherwise smooth migration was interrupted by a discontinuous change of Neptune's semimajor axis when Neptune reached ~28 AU. Before the discontinuity happened, planetesimals located at ~40 AU were swept into Neptune's 2:1 resonance, and were carried with the migrating resonance outwards. The 2:1 resonance was at ~44 AU when Neptune reached ~28 AU. If Neptune's semimajor axis changed by fraction of AU at this point, perhaps because Neptune was scattered off of another planet, the 2:1 population would have been released at ~44 AU, and would remain there to this day. We show that the orbital distribution of bodies produced in this model provides a good match to the orbital properties of the kernel. If Neptune migration was conveniently slow after the jump, the sweeping 2:1 resonance would deplete the population of bodies at ~45-47 AU, thus contributing to the paucity of the low-inclination orbits in this region. Special provisions, probably related to inefficiencies in the accretional growth of sizable objects, are still needed to explain why only a few low-inclination bodies have been so far detected beyond ~47 AU.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Evidence for Slow Migration of Neptune from the Inclination Distribution of Kuiper Belt Objects

Much of the dynamical structure of the Kuiper belt can be explained if Neptune migrated over several AU, and/or if Neptune was scattered to an eccentric orbit during planetary instability. An outstanding problem with the existing formation models is that the distribution of orbital inclinations they predicted is narrower than the one inferred from observations. Here we perform numerical simulations of Kuiper belt formation starting from an initial state with Neptune at 20 < a_{N,0} < 30 AU and a dynamically cold outer disk extending from beyond a_{N,0} to 30 AU. Neptune's orbit is migrated into the disk on an e-folding timescale 1 <= tau <= 100 Myr. A small fraction (~10^{-3}) of the disk planetesimals become implanted into the Kuiper belt in the simulations. By analyzing the orbital distribution of the implanted bodies in different cases we find that the inclination constraint implies that tau >= 10 Myr and a_{N,0} <= 25 AU. The models with tau < 10 Myr do not satisfy the inclination constraint, because there is not enough time for various dynamical processes to raise inclinations. The slow migration of Neptune is consistent with other Kuiper belt constraints, and with recently developed models of planetary instability/migration. Neptune's eccentricity and inclination are never large in these models (e_N<0.1, i_N<2 deg), as required to avoid excessive orbital excitation in the >40 AU region, where the Cold Classicals presumably formed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Hunt for Exomoons with Kepler (HEK): V. A Survey of 41 Planetary Candidates for Exomoons

We present a survey of 41 Kepler Objects of Interest (KOIs) for exomoons using Bayesian photodynamics, more than tripling the number of KOIs surveyed with this technique. We find no compelling evidence for exomoons although thirteen KOIs yield spurious detections driven by instrumental artifacts, stellar activity and/or perturbations from unseen bodies. Regarding the latter, we find seven KOIs exhibiting >5 sigma evidence of transit timing variations, including the 'mega-Earth' Kepler-10c, likely indicating an additional planet in that system. We exploit the moderately large sample of 57 unique KOIs surveyed to date to infer several useful statistics. For example, although there is a diverse range in sensitivities, we find that we are sensitive to Pluto-Charon mass-ratio systems for ~40% of KOIs studied and Earth-Moon mass-ratios for 1 in 8 cases. In terms of absolute mass, our limits probe down to 1.7 Ganymede masses, with a sensitivity to Earth-mass moons for 1 in 3 cases studied and to the smallest moons capable of sustaining an Earth-like atmosphere (0.3 Earth masses) for 1 in 4. Despite the lack of positive detections to date, we caution against drawing conclusions yet, since our most interesting objects remain under analysis. Finally, we point out that had we searched for the photometric transit signals of exomoons alone, rather than using photodynamics, we estimate that 1 in 4 KOIs would have erroneously been concluded to harbor exomoons due to residual time correlated noise in the Kepler data, posing a serious problem for alternative methods.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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 Nesvorny & Morbidelli (2012), in which the obliquities of Jupiter and Saturn can be excited when the spin-orbit resonances occur. The most notable resonances occur when the $s_7$ and $s_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_8$ resonance. This would set strict constrains 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_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_8$ resonance. This probably happened during the late stage of planetary migration when the evolution of the $s_8$ frequency was very slow, and the conditions for capture into the spin-orbit resonance with $s_8$ were satisfied. However, whether or not Saturn is in the spin-orbit resonance with $s_8$ at the present time is not clear, because the existing observations of Saturn's spin precession and internal structure models have significant uncertainties.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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 deg, 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 four out of ten instability cases studied here, the orbital perturbations were too large to be plausible. In one case, Iapetus's orbit was practically unneffected. 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 some 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 present value (~30% of trials lead to >5 deg) while leaving its orbital eccentricity low.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Outward migration of Jupiter and Saturn in 3:2 or 2:1 resonance in radiative disks: implications for the Grand Tack and Nice models

Embedded in the gaseous protoplanetary disk, Jupiter and Saturn naturally become trapped in 3:2 resonance and migrate outward. This serves as the basis of the Grand Tack model. However, previous hydrodynamical simulations were restricted to isothermal disks, with moderate aspect ratio and viscosity. Here we simulate the orbital evolution of the gas giants in disks with viscous heating and radiative cooling. We find that Jupiter and Saturn migrate outward in 3:2 resonance in modest-mass ($M_{disk} \approx M_{MMSN}$, where MMSN is the "minimum-mass solar nebula") disks with viscous stress parameter $α$ between $10^{-3}$ and $10^{-2} $. In disks with relatively low-mass ($M_{disk} \lesssim M_{MMSN}$) , Jupiter and Saturn get captured in 2:1 resonance and can even migrate outward in low-viscosity disks ($α\le 10^{-4}$). Such disks have a very small aspect ratio ($h\sim 0.02-0.03$) that favors outward migration after capture in 2:1 resonance, as confirmed by isothermal runs which resulted in a similar outcome for $h \sim 0.02$ and $α\le 10^{-4}$. We also performed N-body runs of the outer Solar System starting from the results of our hydrodynamical simulations and including 2-3 ice giants. After dispersal of the gaseous disk, a Nice model instability starting with Jupiter and Saturn in 2:1 resonance results in good Solar Systems analogs. We conclude that in a cold Solar Nebula, the 2:1 resonance between Jupiter and Saturn can lead to outward migration of the system, and this may represent an alternative scenario for the evolution of the Solar System.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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_N) planet wide of the 5:3 resonance (period 80.1 d) or a 1 Saturn-mass planet wide of the 7:3 resonance (period 109.2 d). 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 hour 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_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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Hunt for Exomoons with Kepler (HEK): IV. A Search for Moons around Eight M-Dwarfs

With their smaller radii and high cosmic abundance, transiting planets around cool stars hold a unique appeal. As part of our on-going project to measure the occurrence rate of extrasolar moons, we here present results from a survey focussing on eight Kepler planetary candidates associated with M-dwarfs. Using photodynamical modeling and Bayesian multimodal nested sampling, we find no compelling evidence for an exomoon in these eight systems. Upper limits on the presence of such bodies probe down to $\sim0.4M_{\oplus}$ in the best case. For KOI-314, we are able to confirm the planetary nature of two out of the three known transiting candidates using transit timing variations. Of particular interest is KOI-314c, which is found to have a mass of $1.0_{-0.3}^{+0.4}M_{\oplus}$, making it the lowest mass transiting planet discovered to date. With a radius of $1.61_{-0.15}^{+0.16}R_{\oplus}$, this Earth-mass world is likely enveloped by a significant gaseous envelope comprising $\geq17_{-13}^{+12}$% of the planet by radius. We find evidence to support the planetary nature of KOI-784 too via transit timing, but we advocate further observations to verify the signals. In both systems, we infer that the inner planet has a higher density than the outer world, which may be indicative of photo-evaporation. These results highlight both the ability of Kepler to search for sub-Earth mass moons and the exciting ancillary science which often results from such efforts.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Possible Moon of Kepler-90g is a False Positive

The discovery of an exomoon would provide deep insights into planet formation and the habitability of planetary systems, with transiting examples being particularly sought after. Of the hundreds of Kepler planets now discovered, the seven-planet system Kepler-90 is unusual for exhibiting an unidentified transit-like signal in close proximity to one of the transits of the long-period gas-giant Kepler-90g, as noted by Cabrera et. al. (2014). As part of the 'Hunt for Exomoons with Kepler' (HEK) project, we investigate this possible exomoon signal and find it passes all conventional photometric, dynamical and centroid diagnostic tests. However, pixel-level light curves indicate that the moon-like signal occurs on nearly all of the target's pixels, which we confirm using a novel way of examining pixel-level data which we dub the 'transit centroid'. This test reveals that the possible exomoon to Kepler-90g is likely a false positive, perhaps due to a cosmic ray induced Sudden Pixel Sensitivity Dropout (SPSD). This work highlights the extreme care required for seeking non-periodic low-amplitude transit signals, such as exomoons.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Transit timing variations for planets coorbiting in the horseshoe regime

While not detected yet, pairs of exoplanets in the 1:1 mean motion resonance probably exist. Low eccentricity, near-planar orbits, which in the comoving frame follow the horseshoe trajectories, are one of the possible stable configurations. Here we study transit timing variations 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 transit timing variations can directly constrain the (i) ratio of planetary masses, and (ii) 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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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 (Nesvorny et al. 2013). 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Capture of Trojans by Jumping Jupiter

Jupiter Trojans are thought to be survivors of a much larger population of planetesimals that existed in the planetary region when planets formed. They can provide important constraints on the mass and properties of the planetesimal disk, and its dispersal during planet migration. Here we tested a possibility that the Trojans were captured during the early dynamical instability among the outer planets (aka the Nice model), when the semimajor axis of Jupiter was changing as a result of scattering encounters with an ice giant. The capture occurs in this model when Jupiter's orbit and its Lagrange points become radially displaced in a scattering event and fall into a region populated by planetesimals (that previously evolved from their natal transplanetary disk to ~5 AU during the instability). Our numerical simulations of the new capture model, hereafter jump capture, satisfactorily reproduce the orbital distribution of the Trojans and their total mass. The jump capture is potentially capable of explaining the observed asymmetry in the number of leading and trailing Trojans. We find that the capture probability is (6-8)x10^-7 for each particle in the original transplanetary disk, implying that the disk contained (3-4)x10^7 planetesimals with absolute magnitude H<9 (corresponding to diameter D=80 km for a 7% albedo). The disk mass inferred from this work, M_disk=14-28 M_Earth, is consistent with the mass deduced from recent dynamical simulations of the planetary instability.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

KOI-142, the King of Transit Variations, is a Pair of Planets near the 2:1 Resonance

The Transit Timing Variations (TTVs) can be used as a diagnostic of gravitational interactions between planets in a multi-planet system. Many Kepler Objects of Interest (KOIs) exhibit significant TTVs, but KOI-142.01 stands out among them with an unrivaled, 12-hour TTV amplitude. Here we report a thorough analysis of KOI-142.01's transits. We discover periodic Transit Duration Variations (TDVs) of KOI-142.01 that are nearly in phase with the observed TTVs. We show that KOI-142.01's TTVs and TDVs uniquely detect a non-transiting companion with a mass 0.7 that of Jupiter (KOI-142c). KOI-142.01's mass inferred from the transit variations is consistent with the measured transit depth, suggesting a Neptune class planet (KOI-142b). The orbital period ratio P_c/P_b=2.03 indicates that the two planets are just wide of the 2:1 resonance. The present dynamics of this system, characterized here in detail, can be used to test various formation theories that have been proposed to explain the near-resonant pairs of exoplanets.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Hunt for Exomoons with Kepler (HEK): II. Analysis of Seven Viable Satellite-Hosting Planet Candidates

From the list of 2321 transiting planet candidates announced by the Kepler Mission, we select seven targets with favorable properties for the capacity to dynamically maintain an exomoon and present a detectable signal. These seven candidates were identified through our automatic target selection (TSA) algorithm and target selection prioritization (TSP) filtering, whereby we excluded systems exhibiting significant time-correlated noise and focussed on those with a single transiting planet candidate of radius less than 6 Earth radii. We find no compelling evidence for an exomoon around any of the seven KOIs but constrain the satellite-to-planet mass ratios for each. For four of the seven KOIs, we estimate a 95% upper quantile of M_S/M_P<0.04, which given the radii of the candidates, likely probes down to sub-Earth masses. We also derive precise transit times and durations for each candidate and find no evidence for dynamical variations in any of the KOIs. With just a few systems analyzed thus far in the in-going HEK project, projections on eta-moon would be premature, but a high frequency of large moons around Super-Earths/Mini-Neptunes would appear to be incommensurable with our results so far.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Hunt for Exomoons with Kepler (HEK): III. The First Search for an Exomoon around a Habitable-Zone Planet

Kepler-22b is the first transiting planet to have been detected in the habitable-zone of its host star. At 2.4 Earth radii, Kepler-22b is too large to be considered an Earth-analog, but should the planet host a moon large enough to maintain an atmosphere, then the Kepler-22 system may yet possess a telluric world. Aside from being within the habitable-zone, the target is attractive due to the availability of previously measured precise radial velocities and low intrinsic photometric noise, which has also enabled asteroseismology studies of the star. For these reasons, Kepler-22b was selected as a target-of-opportunity by the 'Hunt for Exomoons with Kepler' (HEK) project. In this work, we conduct a photodynamical search for an exomoon around Kepler-22b leveraging the transits, radial velocities and asteroseismology plus several new tools developed by the HEK project to improve exomoon searches. We find no evidence for an exomoon around the planet and exclude moons of mass >0.5 Earth masses to 95% confidence. By signal injection and blind retrieval, we demonstrate that an Earth-like moon is easily detected for this planet even when the time-correlated noise of the data set is taken into account. We provide updated parameters for the planet Kepler-22b including a revised mass of <53 Earth masses to 95% confidence and an eccentricity of 0.13(-0.13)(+0.36) by exploiting Single-body Asterodensity Profiling (SAP). Finally, we show that Kepler-22b has a >95% probability of being within the empirical habitable-zone but a <5% probability of being within the conservative habitable-zone.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Delivery of Dark Material to Vesta via Carbonaceous Chondritic Impacts

NASA's Dawn spacecraft observations of asteroid (4) Vesta reveal a surface with the highest albedo and color variation of any asteroid we have observed so far. Terrains rich in low albedo dark material (DM) have been identified using Dawn Framing Camera (FC) 0.75 μm filter images in several geologic settings: associated with impact craters (in the ejecta blanket material and/or on the crater walls and rims); as flow-like deposits or rays commonly associated with topographic highs; and as dark spots (likely secondary impacts) nearby impact craters. This DM could be a relic of ancient volcanic activity or exogenic in origin. We report that the majority of the spectra of DM are similar to carbonaceous chondrite meteorites mixed with materials indigenous to Vesta. Using high-resolution seven color images we compared DM color properties (albedo, band depth) with laboratory measurements of possible analog materials. Band depth and albedo of DM are identical to those of carbonaceous chondrite xenolith-rich howardite Mt. Pratt (PRA) 04401. Laboratory mixtures of Murchison CM2 carbonaceous chondrite and basaltic eucrite Millbillillie also show band depth and albedo affinity to DM. Modeling of carbonaceous chondrite abundance in DM (1-6 vol%) is consistent with howardite meteorites. We find no evidence for large-scale volcanism (exposed dikes/pyroclastic falls) as the source of DM. Our modeling efforts using impact crater scaling laws and numerical models of ejecta reaccretion suggest the delivery and emplacement of this DM on Vesta during the formation of the ~400 km Veneneia basin by a low-velocity (<2 km/sec) carbonaceous impactor. This discovery is important because it strengthens the long-held idea that primitive bodies are the source of carbon and probably volatiles in the early Solar System.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dynamics of pebbles in the vicinity of a growing planetary embryo: hydro-dynamical simulations

Understanding the growth of the cores of giant planets is a difficult problem. Recently, Lambrechts and Johansen (2012; LJ12) proposed a new model in which the cores grow by the accretion of pebble-size objects, as the latter drift towards the star due to gas drag. Here, we investigate the dynamics of pebble-size objects in the vicinity of planetary embryos of 1 and 5 Earth masses and the resulting accretion rates. We use hydrodynamical simulations, in which the embryo influences the dynamics of the gas and the pebbles suffer gas drag according to the local gas density and velocities. The pebble dynamics in the vicinity of the planetary embryo is non-trivial, and it changes significantly with the pebble size. Nevertheless, the accretion rate of the embryo that we measure is within an order of magnitude of the rate estimated in LJ12 and tends to their value with increasing pebble-size. We conclude that the model by LJ12 has the potential to explain the rapid growth of giant planet cores. The actual accretion rates however, depend on the surface density of pebble size objects in the disk, which is unknown to date.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Statistical Study of the Early Solar System's Instability with 4, 5 and 6 Giant Planets

Several properties of the Solar System, including the wide radial spacing and orbital eccentricities of giant planets, can be explained if the early Solar System evolved through a dynamical instability followed by migration of planets in the planetesimal disk. Here we report the results of a statistical study, in which we performed nearly 10^4 numerical simulations of planetary instability starting from hundreds of different initial conditions. We found that the dynamical evolution is typically too violent, if Jupiter and Saturn start in the 3:2 resonance, leading to ejection of at least one ice giant from the Solar System. Planet ejection can be avoided if the mass of the transplanetary disk of planetesimals was large (M_disk>50 M_Earth), but we found that a massive disk would lead to excessive dynamical damping (e.g., final e_55 < 0.01 compared to present e_55=0.044, where e_55 is the amplitude of the fifth eccentric mode in the Jupiter's orbit), and to smooth migration that violates constraints from the survival of the terrestrial planets. Better results were obtained when the Solar System was assumed to have five giant planets initially and one ice giant, with the mass comparable to that of Uranus and Neptune, was ejected into interstellar space by Jupiter. The best results were obtained when the ejected planet was placed into the external 3:2 or 4:3 resonance with Saturn and M_disk ~ 20 M_Earth. The range of possible outcomes is rather broad in this case, indicating that the present Solar System is neither a typical nor expected result for a given initial state, and occurs, in best cases, with only a ~5% probability (as defined by the success criteria described in the main text). The case with six giant planets shows interesting dynamics but does offer significant advantages relative to the five planet case.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Detection and Characterization of a Nontransiting Planet by Transit Timing Variations

The Kepler Mission is monitoring the brightness of ~150,000 stars searching for evidence of planetary transits. As part of the "Hunt for Exomoons with Kepler" (HEK) project, we report a planetary system with two confirmed planets and one candidate planet discovered using the publicly available data for KOI-872. Planet b transits the host star with a period P_b=33.6d and exhibits large transit timing variations indicative of a perturber. Dynamical modeling uniquely detects an outer nontransiting planet c near the 5:3 resonance (P_c=57.0d) of mass 0.37 times that of Jupiter. Transits of a third planetary candidate are also found: a 1.7-Earth radius super-Earth with a 6.8d period. Our analysis indicates a system with nearly coplanar and circular orbits, reminiscent of the orderly arrangement within the solar system.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Hunt for Exomoons with Kepler (HEK): I. Description of a New Observational Project

Two decades ago, empirical evidence concerning the existence and frequency of planets around stars, other than our own, was absent. Since this time, the detection of extrasolar planets from Jupiter-sized to most recently Earth-sized worlds has blossomed and we are finally able to shed light on the plurality of Earth-like, habitable planets in the cosmos. Extrasolar moons may also be frequent habitable worlds but their detection or even systematic pursuit remains lacking in the current literature. Here, we present a description of the first systematic search for extrasolar moons as part of a new observational project called "The Hunt for Exomoons with Kepler" (HEK). The HEK project distills the entire list of known transiting planet candidates found by Kepler (2326 at the time of writing) down to the most promising candidates for hosting a moon. Selected targets are fitted using a multimodal nested sampling algorithm coupled with a planet-with-moon light curve modelling routine. By comparing the Bayesian evidence of a planet-only model to that of a planet-with-moon, the detection process is handled in a Bayesian framework. In the case of null detections, upper limits derived from posteriors marginalised over the entire prior volume will be provided to inform the frequency of large moons around viable planetary hosts, eta-moon. After discussing our methodologies for target selection, modelling, fitting and vetting, we provide two example analyses.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamical Model for the Zodiacal Cloud and Sporadic Meteors

The solar system is dusty, and would become dustier over time as asteroids collide and comets disintegrate, except that small debris particles in interplanetary space do not last long. They can be ejected from the solar system by Jupiter, thermally destroyed near the Sun, or physically disrupted by collisions. Also, some are swept by the Earth (and other planets), producing meteors. Here we develop a dynamical model for the solar system meteoroids and use it to explain meteor radar observations. We find that the Jupiter Family Comets (JFCs) are the main source of the prominent concentrations of meteors arriving to the Earth from the helion and antihelion directions. To match the radiant and orbit distributions, as measured by the Canadian Meteor Orbit Radar (CMOR) and Advanced Meteor Orbit Radar (AMOR), our model implies that comets, and JFCs in particular, must frequently disintegrate when reaching orbits with low perihelion distance. Also, the collisional lifetimes of millimeter particles may be longer (>10^5 yr at 1 AU) than postulated in the standard collisional models (10^4 yr at 1 AU), perhaps because these chondrule-sized meteoroids are stronger than thought before. Using observations of the Infrared Astronomical Satellite (IRAS) to calibrate the model, we find that the total cross section and mass of small meteoroids in the inner solar system are (1.7-3.5)x10^11 km^2 and 4x10^19 g, respectively, in a good agreement with previous studies. The mass input required to keep the Zodiacal Cloud (ZC) in a steady state is estimated to be 10^4-10^5 kg/s. The input is up to 10 times larger than found previously, mainly because particles released closer to the Sun have shorter collisional lifetimes, and need to be supplied at a faster rate.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamics of Dust Particles Released from Oort Cloud Comets and Their Contribution to Radar Meteors

The Oort Cloud Comets (OCCs), exemplified by the Great Comet of 1997 (Hale-Bopp), are occasional visitors from the heatless periphery of the solar system. Previous works hypothesized that a great majority of OCCs must physically disrupt after one or two passages through the inner solar system, where strong thermal gradients can cause phase transitions or volatile pressure buildup. Here we study the fate of small debris particles produced by OCC disruptions to determine whether the imprints of a hypothetical population of OCC meteoroids can be found in the existing meteor radar data. We find that OCC particles with diameters D<10 um are blown out from the solar system by radiation pressure, while those with D>1 um have a very low Earth-impact probability. The intermediate particle sizes, D=100 um, represent a sweet spot. About 1% of these particles orbitally evolve by Poynting-Robertson drag to reach orbits with semimajor axis a=1 AU. They are expected to produce meteors with radiants near the apex of the Earth's orbital motion. We find that the model distributions of their impact speeds and orbits provide a good match to radar observations of apex meteors, except for the eccentricity distribution, which is more skewed toward e=1 in our model. Finally, we propose an explanation for the long-standing problem in meteor science related to the relative strength of apex and helion/antihelion sources. As we show in detail, the observed trend, with the apex meteors being more prominent in observations of highly sensitive radars, can be related to orbital dynamics of particles released on the long-period orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Observed Binary Fraction Sets Limits on the Extent of Collisional Grinding in the Kuiper Belt

The size distribution in the cold classical Kuiper belt can be approximated by two idealized power laws: one with steep slope for radii R>R* and one with shallow slope for R<R*, where R*~25-50 km. Previous works suggested that the SFD roll-over at R* can be the result of extensive collisional grinding in the Kuiper belt that led to the catastrophic disruption of most bodies with R<R*. Here we use a new code to test the effect of collisions in the Kuiper belt. We find that the observed roll-over could indeed be explained by collisional grinding provided that the initial mass in large bodies was much larger than the one in the present Kuiper belt, and was dynamically depleted. In addition to the size distribution changes, our code also tracks the effects of collisions on binary systems. We find that it is generally easier to dissolve wide binary systems, such as the ones existing in the cold Kuiper belt today, than to catastrophically disrupt objects with R~R*. Thus, the binary survival sets important limits on the extent of collisional grinding in the Kuiper belt. We find that the extensive collisional grinding required to produce the SFD roll-over at R* would imply a strong gradient of the binary fraction with R and separation, because it is generally easier to dissolve binaries with small components and/or those with wide orbits. The expected binary fraction for R<R* is <0.1. The present observational data do not show such a gradient. Instead, they suggest a large binary fraction of ~0.4 for R=30-40 km. This may indicate that the roll-over was not produced by disruptive collisions, but is instead a fossil remnant of the KBO formation process.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Young Solar System's Fifth Giant Planet?

Recent studies of solar system formation suggest that the solar system's giant planets formed and migrated in the protoplanetary disk to reach resonant orbits with all planets inside 15 AU from the Sun. After the gas disk's dispersal, Uranus and Neptune were likely scattered by gas giants, and approached their current orbits while dispersing the transplanetary disk of planetesimals, whose remains survived to this time in the region known as the Kuiper belt. Here we performed N-body integrations of the scattering phase between giant planets in an attempt to determine which initial states are plausible. We found that the dynamical simulations starting with a resonant system of four giant planets have a low success rate in matching the present orbits of giant planets, and various other constraints (e.g., survival of the terrestrial planets). The dynamical evolution is typically too violent, if Jupiter and Saturn start in the 3:2 resonance, and leads to final systems with fewer than four planets. Several initial states stand out in that they show a relatively large likelihood of success in matching the constraints. Some of the statistically best results were obtained when assuming that the solar system initially had five giant planets and one ice giant, with the mass comparable to that of Uranus and Neptune, was ejected to interstellar space by Jupiter. This possibility appears to be conceivable in view of the recent discovery of a large number free-floating planets in interstellar space, which indicates that planet ejection should be common.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Do Planetary Encounters Reset Surfaces of Near Earth Asteroids?

Processes such as the solar wind sputtering and micrometeorite impacts can modify optical properties of surfaces of airless bodies. This explains why spectra of the main belt asteroids, exposed to these `space weathering' processes over eons, do not match the laboratory spectra of ordinary chondrite (OC) meteorites. In contrast, an important fraction of Near Earth Asteroids (NEAs), defined as Q-types in the asteroid taxonomy, display spectral attributes that are a good match to OCs. Here we study the possibility that the Q-type NEAs underwent recent encounters with the terrestrial planets and that the tidal gravity (or other effects) during these encounters exposed fresh OC material on the surface (thus giving it the Q-type spectral properties). We used numerical integrations to determine the statistics of encounters of NEAs to planets. The results were used to calculate the fraction and orbital distribution of Q-type asteroids expected in the model as a function of the space weathering timescale, t_sw (see main text for definition), and maximum distance, r*, at which planetary encounters can reset the surface. We found that t_sw ~ 1e6 yr (at 1 AU) and r* ~ 5 R_pl, where R_pl is the planetary radius, best fit the data. Values t_sw < 1e5 yr would require that r* > 20 R_pl, which is probably implausible because these very distant encounters should be irrelevant. Also, the fraction of Q-type NEAs would be probably much larger than the one observed if t_sw > 1e7 yr. We found that t_sw ~ q^2, where q is the perihelion distance, expected if the solar wind sputtering controls t_sw, provides a better match to the orbital distribution of Q-type NEAs than models with fixed t_sw. We also discuss how the Earth magnetosphere and radiation effects such as YORP can influence the spectral properties of NEAs.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Asteroids Were Born Big

How big were the first planetesimals? We attempt to answer this question by conducting coagulation simulations in which the planetesimals grow by mutual collisions and form larger bodies and planetary embryos. The size frequency distribution (SFD) of the initial planetesimals is considered a free parameter in these simulations, and we search for the one that produces at the end objects with a SFD that is consistent with asteroid belt constraints. We find that, if the initial planetesimals were small (e.g. km-sized), the final SFD fails to fulfill these constraints. In particular, reproducing the bump observed at diameter D~100km in the current SFD of the asteroids requires that the minimal size of the initial planetesimals was also ~100km. This supports the idea that planetesimals formed big, namely that the size of solids in the proto-planetary disk ``jumped'' from sub-meter scale to multi-kilometer scale, without passing through intermediate values. Moreover, we find evidence that the initial planetesimals had to have sizes ranging from 100 to several 100km, probably even 1,000km, and that their SFD had to have a slope over this interval that was similar to the one characterizing the current asteroids in the same size-range. This result sets a new constraint on planetesimal formation models and opens new perspectives for the investigation of the collisional evolution in the asteroid and Kuiper belts as well as of the accretion of the cores of the giant planets.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Cometary Origin of the Zodiacal Cloud and Carbonaceous Micrometeorites

The zodiacal cloud is a thick circumsolar disk of small debris particles produced by asteroid collisions and comets. Here, we present a zodiacal cloud model based on the orbital properties and lifetimes of comets and asteroids, and on the dynamical evolution of dust after ejection. The model is quantitatively constrained by IRAS observations of thermal emission, but also qualitatively consistent with other zodiacal cloud observations. We find that 85-95% of the observed mid-infrared emission is produced by particles from the Jupiter-family comets (JFCs) and $<$10% by dust from long period comets. Asteroidal dust is found to be present at $<$10%. We suggest that spontaneous disruptions of JFCs, rather than the usual cometary activity driven by sublimating volatiles, is the main mechanism that librates cometary particles into the zodiacal cloud. Our results imply that JFC particles represent $\sim$85% of the total mass influx at Earth. Since their atmospheric entry speeds are typically low ($\approx$14.5 km s$^{-1}$ mean for D=100-200 $μ$m with $\approx$12 km s$^{-1}$ being the most common case), many JFC grains should survive frictional heating and land on the Earth's surface. This explains why most micrometeorites collected in antarctic ice have primitive carbonaceous composition. The present mass of the inner zodiacal cloud at $<$5 AU is estimated to be 1-$2\times10^{19}$ g, mainly in D=100-200 $μ$m particles. The inner zodiacal cloud should have been $>10^4$ times brighter during the Late Heavy Bombardment (LHB) epoch $\approx$3.8 Gyr ago, when the outer planets scattered numerous comets into the inner solar system. The bright debris disks with a large 24-$μ$m excess observed around mature stars may be an indication of massive cometary populations existing in those systems.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Fast Inversion Method for Determination of Planetary Parameters from Transit Timing Variations

The Transit Timing Variation (TTV) method relies on monitoring changes in timing of transits of known exoplanets. Non-transiting planets in the system can be inferred from TTVs by their gravitational interaction with the transiting planet. The TTV method is sensitive to low-mass planets that cannot be detected by other means. Here we describe a fast algorithm that can be used to determine the mass and orbit of the non-transiting planets from the TTV data. We apply our code, ttvim.f, to a wide variety of planetary systems to test the uniqueness of the TTV inversion problem and its dependence on the precision of TTV observations. We find that planetary parameters, including the mass and mutual orbital inclination of planets, can be determined from the TTV datasets that should become available in near future. Unlike the radial velocity technique, the TTV method can therefore be used to characterize the inclination distribution of multi-planet systems.

0 citations0 reviewsNo code linkedOpen access