Source author record

Kevin J. Walsh

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2023arXiv

Tidal disruption of near-Earth asteroids during close encounters with terrestrial planets

Numerical modeling has long suggested that gravitationally-bound (or so-called rubble-pile) near-Earth asteroids (NEAs) can be destroyed by tidal forces during close and slow encounters with terrestrial planets. However, tidal disruptions of NEAs have never been directly observed nor have they been directly attributed to any families of NEAs. Here we show population-level evidence for the tidal disruption of NEAs during close encounters with the Earth and Venus. Debiased model distributions of NEA orbits and absolute magnitudes based on observations by the Catalina Sky Survey during 2005--2012 underpredict the number of NEAs with perihelion distances coinciding with the semimajor axes of Venus and the Earth. A detailed analysis of the orbital distributions of the excess NEAs shows that their characteristics agree with the prediction for tidal disruptions, and they cannot be explained by observational selection effects or orbital dynamics. Accounting for tidal disruptions in evolutionary models of the NEA population partly bridges the gap between the predicted rate of impacts by asteroids with diameters of tens of meters and observed statistics of fireballs in the same size range.

preprint2022arXiv

Chemical Habitability: Supply and Retention of Life's Essential Elements During Planet Formation

Carbon, Hydrogen, Nitrogen, Oxygen, Phosphorus and Sulfur (CHNOPS) play key roles in the origin and proliferation of life on Earth. Given the universality of physics and chemistry, not least the ubiquity of water as a solvent and carbon as a backbone of complex molecules, CHNOPS are likely crucial to most habitable worlds. To help guide and inform the search for potentially habitable and ultimately inhabited environments, we begin by summarizing the CHNOPS budget of various reservoirs on Earth, their role in shaping our biosphere, and their origins in the Solar Nebula. We then synthesize our current understanding of how these elements behave and are distributed in diverse astrophysical settings, tracing their journeys from synthesis in dying stars to molecular clouds, protoplanetary settings, and ultimately temperate rocky planets around main sequence stars. We end by identifying key branching points during this journey, highlighting instances where a forming planets' distribution of CHNOPS can be altered dramatically, and speculating about the consequences for the chemical habitability of these worlds.

preprint2020arXiv

Charon: A brief history of tides

In 2015, the New Horizons spacecraft flew past Pluto and its moon Charon, providing the first clear look at the surface of Charon. New Horizons images revealed an ancient surface, a large, intricate canyon system, and many fractures, among other geologic features. Here, we assess whether tidal stresses played a significant role in the formation of tensile fractures on Charon. Although presently in a circular orbit, most scenarios for the orbital evolution of Charon include an eccentric orbit for some period of time and possibly an internal ocean. Past work has shown that these conditions could have generated stresses comparable in magnitude to other tidally fractured moons, such as Europa and Enceladus. However, we find no correlation between observed fracture orientations and those predicted to form due to eccentricity-driven tidal stress. It thus seems more likely that the orbit of Charon circularized before its ocean froze, and that either tidal stresses alone were insufficient to fracture the surface or subsequent resurfacing remove these ancient fractures.

preprint2020arXiv

Potential Themis Family Asteroid Contribution to the Jupiter-Family Comet Population

Recent dynamical analyses suggest that some Jupiter family comets (JFCs) may originate in the main asteroid belt instead of the outer solar system. This possibility is particularly interesting given evidence that icy main-belt objects are known to be present in the Themis asteroid family. We report results from dynamical analyses specifically investigating the possibility that icy Themis family members could contribute to the observed population of JFCs. Numerical integrations show that such dynamical evolution is indeed possible via a combination of eccentricity excitation apparently driven by the nearby 2:1 mean-motion resonance with Jupiter, gravitational interactions with planets other than Jupiter, and the Yarkovsky effect. We estimate that, at any given time, there may be tens of objects from the Themis family on JFC-like orbits with the potential to mimic active JFCs from the outer solar system, although not all, or even any, may necessarily be observably active. We find that dynamically evolved Themis family objects on JFC-like orbits have semimajor axes between 3.15 au and 3.40 au for the vast majority of their time on such orbits, consistent with the strong role that the 2:1 mean-motion resonance with Jupiter likely plays in their dynamical evolution. We conclude that a contribution from the Themis family to the active JFC population is plausible, although further work is needed to better characterize this contribution.

preprint2020arXiv

Spin-driven evolution of asteroids' top-shapes at fast and slow spins seen from (101955) Bennu and (162173) Ryugu

Proximity observations by OSIRIS-REx and Hayabusa2 provided clues on the shape evolution processes of the target asteroids, (101955) Bennu and (162173) Ryugu. Their oblate shapes with equatorial ridges, or the so-called top shapes, may have evolved due to their rotational conditions at present and in the past. Different shape evolution scenarios were previously proposed; Bennu's top shape may have been driven by surface processing, while Ryugu's may have been developed due to large deformation. These two scenarios seem to be inconsistent. Here, we revisit the structural analyses in earlier works and fill a gap to connect these explanations. We also apply a semi-analytical technique for computing the cohesive strength distribution in a uniformly rotating triaxial ellipsoid to characterize the global failure of top-shaped bodies. Assuming that the structure is uniform, our semi-analytical approach describes the spatial variations in failed regions at different spin periods; surface regions are the most sensitive at longer spin periods, while interiors fail structurally at shorter spin periods. This finding suggests that the shape evolution of a top shape may vary due to rotation and internal structure, which can explain the different evolution scenarios of Bennu's and Ryugu's top shapes. We interpret our results as the indications of top shapes' various evolution processes.

preprint2016arXiv

Orbital Simulations on Deflecting Near-Earth Objects by Directed Energy

Laser ablation of a Near-Earth Object (NEO) on a collision course with Earth produces a cloud of ejecta which exerts a thrust on the NEO, deflecting it from its original trajectory. Ablation may be performed from afar by illuminating an Earth-targeting asteroid or comet with a stand-off "DE- STAR" system consisting of a large phased-array laser in Earth orbit. Alternatively, a much smaller stand-on "DE-STARLITE" system may travel alongside the target, slowly deflecting it from nearby over a long period. This paper presents orbital simulations comparing the effectiveness of both systems across a range of laser and NEO parameters. Simulated parameters include magnitude, duration and, for the stand-on system, direction of the thrust, as well as the type, size and orbital characteristics of the target NEO. These simulations indicate that deflection distance is approximately proportional to the magnitude of thrust and to the square of the duration of ablation, and is inversely proportional to the mass. Furthermore, deflection distance shows strong dependence on thrust direction with the optimal direction of thrust varying with the duration of laser activity. As one example, consider a typical 325 m asteroid: beginning 15 yr in advance, just 2 N of thrust from a ~20 kW stand-on DE-STARLITE system is sufficient to deflect the asteroid by 2 R_e. Numerous scenarios are discussed as is a practical implementation of such a system consistent with current launch vehicle capabilities.

preprint2016arXiv

Terrestrial Planet Formation from an Annulus

It has been shown that some aspects of the terrestrial planets can be explained, particularly the Earth/Mars mass ratio, when they form from a truncated disk with an outer edge near 1.0 au (Hansen 2009). This has been previously modeled starting from an intermediate stage of growth utilizing pre-formed planetary embryos. We present simulations that were designed to test this idea by following the growth process from km-sized objects located between 0.7 to 1.0 au up to terrestrial planets. The simulations explore initial conditions where the solids in the disk are planetesimals with radii initially between 3 and 300 km, alternately including effects from a dissipating gaseous solar nebula and collisional fragmentation. We use a new Lagrangian code known as LIPAD (Levison et al. 2012), which is a particle-based code that models the fragmentation, accretion and dynamical evolution of a large number of planetesimals, and can model the entire growth process from km-sizes up to planets. A suite of large (Mars mass) planetary embryos is complete in only 1 Myr, containing most of the system mass. A quiescent period then persists for 10-20 Myr characterized by slow diffusion of the orbits and continued accretion of the remaining planetesimals. This is interupted by an instability that leads to embryos crossing orbits and embyro-embryo impacts that eventually produce the final set of planets. While this evolution is different than that found in other works exploring an annulus, the final planetary systems are similar, with roughly the correct number of planets and good Mars-analogs.

preprint2016arXiv

Yarkovsky V-shape identification of asteroid families

There are only a few known main belt (MB) asteroid families with ages greater than 2 Gyr (Brož et al. 2013, Spoto et al. 2015). Estimates based on the family producing collision rate suggest that the lack of >2 Gyr-old families may be due to a selection bias in current techniques used to identify families. Family fragments disperse in their orbital elements, semi-major axis, $a$, eccentricity, $e$, and inclination, $i$, due to secular resonances, close encounters with massive asteroids and the non-gravitational Yarkovsky force. This causes the family fragments to be indistinguishable from the background of the main belt making them more difficult to identify with the hierarchical clustering method (HCM) with increasing family age. The discovery of the Eulalia and new Polana families in the inner belt relied on new techniques because Yarkovsky spreading made them too disperse to be identified using the classical HCM. The technique used to discover the new Polana and Eulalia families is modified to identify asteroid families by searching for correlations between $a$ and asteroid diameter, $D$, or absolute magnitude, $H$. A group of asteroids is identified as a collisional family if its boundary in the $a$ vs. 1/$D$ or $a$ vs. $H$ planes has a characteristic V-shape which is due to the size dependent Yarkovsky spreading. The V-shape boundary is identified with two separate techniques. The first identifies a border by measuring the drop between the number of objects inside and outside of the border. The second identifies the V-shape border by measuring the number density of objects in $a$ vs. 1/$D$,$H$ space. The V-shape techniques are demonstrated on known families Erigone, Vesta, Koronis, and families difficult to identify by HCM such as Flora, Baptistina, new Polana, Eulalia and Karin. Application of the technique in a search for >2 Gyr-old families throughout the MB is discussed.

preprint2015arXiv

Earth and Terrestrial Planet Formation

The growth and composition of Earth is a direct consequence of planet formation throughout the Solar System. We discuss the known history of the Solar System, the proposed stages of growth and how the early stages of planet formation may be dominated by pebble growth processes. Pebbles are small bodies whose strong interactions with the nebula gas lead to remarkable new accretion mechanisms for the formation of planetesimals and the growth of planetary embryos. Many of the popular models for the later stages of planet formation are presented. The classical models with the giant planets on fixed orbits are not consistent with the known history of the Solar System, fail to create a high Earth/Mars mass ratio, and, in many cases, are also internally inconsistent. The successful Grand Tack model creates a small Mars, a wet Earth, a realistic asteroid belt and the mass-orbit structure of the terrestrial planets. In the Grand Tack scenario, growth curves for Earth most closely match a Weibull model. The feeding zones, which determine the compositions of Earth and Venus follow a particular pattern determined by Jupiter, while the feeding zones of Mars and Theia, the last giant impactor on Earth, appear to randomly sample the terrestrial disk. The late accreted mass samples the disk nearly evenly.

preprint2015arXiv

Formation and Evolution of Binary Asteroids

Satellites of asteroids have been discovered in nearly every known small body population, and a remarkable aspect of the known satellites is the diversity of their properties. They tell a story of vast differences in formation and evolution mechanisms that act as a function of size, distance from the Sun, and the properties of their nebular environment at the beginning of Solar System history and their dynamical environment over the next 4.5 Gyr. The mere existence of these systems provides a laboratory to study numerous types of physical processes acting on asteroids and their dynamics provide a valuable probe of their physical properties otherwise possible only with spacecraft. Advances in understanding the formation and evolution of binary systems have been assisted by: 1) the growing catalog of known systems, increasing from 33 to nearly 250 between the Merline et al. (2002) Asteroids III chapter and now, 2) the detailed study and long-term monitoring of individual systems such as 1999 KW4 and 1996 FG3, 3) the discovery of new binary system morphologies and triple systems, 4) and the discovery of unbound systems that appear to be end-states of binary dynamical evolutionary paths. Specifically for small bodies (diameter smaller than 10 km), these observations and discoveries have motivated theoretical work finding that thermal forces can efficiently drive the rotational disruption of small asteroids. Long-term monitoring has allowed studies to constrain the system's dynamical evolution by the combination of tides, thermal forces and rigid body physics. The outliers and split pairs have pushed the theoretical work to explore a wide range of evolutionary end-states.

preprint2015arXiv

Formation and evolution of Pluto's small satellites

Pluto's system of 5 known satellites are in a puzzling orbital configuration. Each of the four small satellites are on low-eccentricity and low-inclination orbits situated near a mean motion resonance with the largest satellite Charon. The Pluto-Charon binary likely formed as a result of a giant impact and so the simplest explanation for the small satellites is that they accreted from debris of that collision. The Pluto-Charon binary has evolved outward since its formation due to tidal forces, which drove them into their current doubly synchronous state. Meanwhile, leftover debris from the formation of Charon was not initially distant enough from Pluto-Charon to explain the orbits of the current small satellites. The outstanding problems of the system are the movement of debris outward and the small satellites location near mean motion resonances with Charon. This work explores the dynamical behavior of collisionally interacting debris orbiting the Pluto-Charon system. While this work specifically tests initial disk and ring configurations designed to mimic the aftermath of the disruption of satellites by heliocentric impactors, we generally find that collisional interactions can help move material outwards and keep otherwise unstable material dynamically bound to the Pluto-Charon system. These processes can produce rings of debris whose orbits evolve rapidly due to collisional processes, with increasing pericenters and decreasing semimajor axes. While these rings and disks of debris eventually build satellites significantly further out than the initial locations of a disrupted satellite, they do not show a strong preference for building satellites in or near mean motion resonances with Charon under a wide array of tested conditions.

preprint2015arXiv

Highly Siderophile Elements in the Earth's Mantle as a Clock for the Moon-forming Impact

According to the generally accepted scenario, the last giant impact on the Earth formed the Moon and initiated the final phase of core formation by melting the Earth's mantle. A key goal of geochemistry is to date this event, but different ages have been proposed. Some argue for an early Moon-forming event, approximately 30 million years (Myr) after the condensation of the first solids in the Solar System, whereas others claim a date later than 50 Myr (and possibly as late as around 100 My) after condensation. Here we show that a Moon-forming event at 40 Myr after condensation, or earlier, is ruled out at a 99.9 per cent confidence level. We use a large number of N-body simulations to demonstrate a relationship between the time of the last giant impact on an Earth-like planet and the amount of mass subsequently added during the era known as Late Accretion. As the last giant impact is delayed, the late-accreted mass decreases in a predictable fashion. This relationship exists within both the classical scenario and the Grand Tack scenario of terrestrial planet formation, and it holds across a wide range of disk conditions. The concentration of highly siderophile elements (HSEs) in Earth's mantle constrains the mass of chondritic material added to Earth during Late Accretion. Using HSE abundance measurements, we determine a Moon-formation age of 95 +/- 32 Myr since condensation. The possibility exists that some late projectiles were differentiated and left an incomplete HSE record in Earth's mantle. Even in this case, various isotopic constraints strongly suggest that the late-accreted mass did not exceed 1 per cent of Earth's mass, and so the HSE clock still robustly limits the timing of the Moon-forming event to significantly later than 40 My after condensation.

preprint2015arXiv

The Compositional Structure of the Asteroid Belt

The past decade has brought major improvements in large-scale asteroid discovery and characterization with over half a million known asteroids and over 100,000 with some measurement of physical characterization. This explosion of data has allowed us to create a new global picture of the Main Asteroid Belt. Put in context with meteorite measurements and dynamical models, a new and more complete picture of Solar System evolution has emerged. The question has changed from "What was the original compositional gradient of the Asteroid Belt?" to "What was the original compositional gradient of small bodies across the entire Solar System?" No longer is the leading theory that two belts of planetesimals are primordial, but instead those belts were formed and sculpted through evolutionary processes after Solar System formation. This article reviews the advancements on the fronts of asteroid compositional characterization, meteorite measurements, and dynamical theories in the context of the heliocentric distribution of asteroid compositions seen in the Main Belt today. This chapter also reviews the major outstanding questions relating to asteroid compositions and distributions and summarizes the progress and current state of understanding of these questions to form the big picture of the formation and evolution of asteroids in the Main Belt. Finally, we briefly review the relevance of asteroids and their compositions in their greater context within our Solar System and beyond.

preprint2015arXiv

The Dynamical Evolution of the Asteroid Belt

The asteroid belt is the leftover of the original planetesimal population in the inner solar system. However, currently the asteroids have orbits with all possible values of eccentricities and inclinations compatible with long-term dynamical stability, whereas the initial planetesimal orbits should have been quasi-circular and almost co-planar. The total mass in the asteroid population is a small fraction of that existing primordially. Also, asteroids with different chemical/mineralogical properties are not ranked in an orderly manner with mean heliocentric distance as one could expect from the existence of a radial gradient of the temperature in the proto-planetary disk, but they are partially mixed. These properties show that the asteroid belt has been severely sculpted by one or a series of processes during its lifetime. This paper reviews the processes that have been proposed so far, discussing the properties that they explain and the problems that they are confronted with. Emphasis is paid to the interplay between the dynamical and the collisional evolution of the asteroid population, which allows the use of the size distribution to constrain the dynamical models. We divide the asteroid belt evolution into three phases. The first phase started during the lifetime of the gaseous proto-planetary disk, when the giant planets formed and presumably experienced large-scale migrations, and continued after the removal of the gas, during the build-up of the terrestrial planets. The second phase occurred after the removal of the gaseous proto-planetary disk and it became particularly lively for the asteroid belt when the giant planets suddenly changed their orbits, as a result of a mutual dynamical instability and the interaction with the trans-Neptunian planetesimal disk. The third phase covers the aftermath of the giant planet instability, until today.

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.

preprint2014arXiv

Properties and evolution of NEO families created by tidal disruption at Earth

We have calculated the coherence and detectable lifetimes of synthetic near-Earth object (NEO) families created by catastrophic disruption of a progenitor as it suffers a very close Earth approach. The closest or slowest approaches yield the most violent `s-class' disruption events. We found that the average slope of the absolute magnitude (H) distribution, $N(H)\propto10^{(0.55\pm0.04)\,H}$, for the fragments in the s-class families is steeper than the slope of the NEO population \citep{mainzer2011} in the same size range. The families remain coherent as statistically significant clusters of orbits within the NEO population for an average of $\barτ_c = (14.7\pm0.6)\times10^3$ years after disruption. The s-class families are detectable with the techniques developed by \citet{fu2005} and \citet{Schunova2012} for an average duration ($\barτ_{det}$) ranging from about 2,000 to about 12,000 years for progenitors in the absolute magnitude ($H_p$) range from 20 to 13 corresponding to diameters in the range from about 0.5 to 10$\km$ respectively. The short detectability lifetime explains why zero NEO families have been discovered to-date. Nonetheless, every tidal disruption event of a progenitor with D$>0.5\km$ is capable of producing several million fragments in the $1\meter$ to $10\meter$ diameter range that can contribute to temporary local density enhancements of small NEOs in Earth's vicinity. We expect that there are about 1,200 objects in the steady state NEO population in this size range due to tidal disruption assuming that one $1\km$ diameter NEO tidally disrupts at Earth every 2,500 years. These objects may be suitable targets for asteroid retrieval missions due to their Earth-like orbits with corresponding low $v_{\infty}$. The fragments from the tidal disruptions at Earth have $\sim5\times$ the collision probability with Earth compared to the background NEO population.

preprint2014arXiv

Terrestrial Planet Formation at Home and Abroad

We review the state of the field of terrestrial planet formation with the goal of understanding the formation of the inner Solar System and low-mass exoplanets. We review the dynamics and timescales of accretion from planetesimals to planetary embryos and from embryos to terrestrial planets. We discuss radial mixing and water delivery, planetary spins and the importance of parameters regarding the disk and embryo properties. Next, we connect accretion models to exoplanets. We first explain why the observed hot Super Earths probably formed by in situ accretion or inward migration. We show how terrestrial planet formation is altered in systems with gas giants by the mechanisms of giant planet migration and dynamical instabilities. Standard models of terrestrial accretion fail to reproduce the inner Solar System. The "Grand Tack" model solves this problem using ideas first developed to explain the giant exoplanets. Finally, we discuss whether most terrestrial planet systems form in the same way as ours, and highlight the key ingredients missing in the current generation of simulations.

preprint2014arXiv

Water Delivery and Giant Impacts in the 'Grand Tack' Scenario

A new model for terrestrial planet formation (Hansen 2009, Walsh et al. 2011) has explored accretion in a truncated protoplanetary disk, and found that such a configuration is able to reproduce the distribution of mass among the planets in the Solar System, especially the Earth/Mars mass ratio, which earlier simulations have generally not been able to match. Walsh et al. tested a possible mechanism to truncate the disk--a two-stage, inward-then-outward migration of Jupiter and Saturn, as found in numerous hydrodynamical simulations of giant planet formation. In addition to truncating the disk and producing a more realistic Earth/Mars mass ratio, the migration of the giant planets also populates the asteroid belt with two distinct populations of bodies--the inner belt is filled by bodies originating inside of 3 AU, and the outer belt is filled with bodies originating from between and beyond the giant planets (which are hereafter referred to as `primitive' bodies). We find here that the planets will accrete on order 1-2% of their total mass from primitive planetesimals scattered onto planet-crossing orbits during the formation of the planets. For an assumed value of 10% for the water mass fraction of the primitive planetesimals, this model delivers a total amount of water comparable to that estimated to be on the Earth today. While the radial distribution of the planetary masses and the dynamical excitation of their orbits are a good match to the observed system, we find that the last giant impact is typically earlier than 20 Myr, and a substantial amount of mass is accreted after that event. However, 5 of the 27 planets larger than half an Earth mass formed in all simulations do experience large late impacts and subsequent accretion consistent with the dating of the Moon-forming impact and the estimated amount of mass accreted by Earth following that event.

preprint2013arXiv

Introducing the Eulalia and new Polana asteroid families: re-assessing primitive asteroid families in the inner Main Belt

The so-called Nysa-Polana complex of asteroids is a diverse and widespread group. It appears to be two overlapping families of different asteroid taxonomies: (44) Nysa is an E-type asteroid with the lowest number in the midst of a predominantly S-type cluster and (142) Polana is a B-type asteroid near the low-albedo B- and C-type cluster. Using the data from the Wide-field Infrared Survey Explorer (WISE) mission we have re-analyzed the region around the Nysa-Polana complex in the inner Main Belt, focusing on the low-albedo population. (142) Polana does not appear to be a member of the family of low-albedo asteroids in the Nysa-Polana complex. Rather, the largest is asteroid (495) Eulalia. This asteroid has never before been linked to this complex for an important dynamical reason: it currently has a proper eccentricity slightly below the range of most of the family members. However, its orbit is very close to the 3:1 mean motion resonance with Jupiter and is in a weak secular resonance. We show that its osculating eccentricity varies widely on short timescales and the averaged value diffuses over long timescales. The diffusive orbit, low-albedo, taxonomic similarity and semimajor axis strongly suggests that despite its current proper eccentricity, (495) Eulalia could have recently been at an orbit very central to the family. Hierarchical Clustering Method tests confirm that at an eccentricity of e=0.15, (495) Eulalia could be the parent of the family. The ``Eulalia family'' was formed between 900--1500 Myr ago, and likely resulted from the breakup of a 100--160 km parent body. There is also compelling evidence for an older and more widespread primitive family in the same region of the asteroid belt parented by asteroid (142) Polana.

preprint2013arXiv

Numerical simulations of granular dynamics. I. Hard-sphere discrete element method and tests

We present a new particle-based (discrete element) numerical method for the simulation of granular dynamics, with application to motions of particles on small solar system body and planetary surfaces. The method employs the parallel N-body tree code pkdgrav to search for collisions and compute particle trajectories. Collisions are treated as instantaneous point-contact events between rigid spheres. Particle confinement is achieved by combining arbitrary combinations of four provided wall primitives, namely infinite plane, finite disk, infinite cylinder, and finite cylinder, and degenerate cases of these. Various wall movements, including translation, oscillation, and rotation, are supported. We provide full derivations of collision prediction and resolution equations for all geometries and motions. Several tests of the method are described, including a model granular "atmosphere" that achieves correct energy equipartition, and a series of tumbler simulations that show the expected transition from tumbling to centrifuging as a function of rotation rate.

preprint2013arXiv

Will Comet ISON (C/2012 S1) Survive Perihelion?

On 2013 November 28 Comet ISON (C/2012 S1) will pass by the Sun with a perihelion distance of 2.7 solar radii. Understanding the possible outcomes for the comet's response to such a close passage by the Sun is important for planning observational campaigns and for inferring ISON's physical properties. We present new numerical simulations and interpret them in context with the historical track record of comet disruptions and of sungrazing comet behavior. Historical data suggest that sizes below ~200 m are susceptible to destruction by sublimation driven mass loss, while we find that for ISON's perihelion distance, densities lower than 0.1 g cm^-3 are required to tidally disrupt a retrograde or non-spinning body. Such low densities are substantially below the range of the best-determined comet nucleus densities, though dynamically new comets such as ISON have few measurements of physical properties. Disruption may occur for prograde rotation at densities up to 0.7 g cm^-3, with the chances of disruption increasing for lower density, faster prograde rotation, and increasing elongation of the nucleus. Given current constraints on ISON's nucleus properties and the typically determined values for these properties among all comets, we find tidal disruption to be unlikely unless other factors (e.g., spin-up via torquing) affect ISON substantially. Whether or not disruption occurs, the largest remnant must be big enough to survive subsequent mass loss due to sublimation in order for ISON to remain a viable comet well after perihelion.

preprint2012arXiv

A low mass for Mars from Jupiter's early gas-driven migration

Jupiter and Saturn formed in a few million years (Haisch et al. 2001) from a gas-dominated protoplanetary disk, and were susceptible to gas-driven migration of their orbits on timescales of only ~100,000 years (Armitage 2007). Hydrodynamic simulations show that these giant planets can undergo a two-stage, inward-then-outward, migration (Masset & Snellgrove 2001, Morbidelli & Crida 2007, Pierens & Nelson 2008). The terrestrial planets finished accreting much later (Klein et al. 2009), and their characteristics, including Mars' small mass, are best reproduced by starting from a planetesimal disk with an outer edge at about one astronomical unit from the Sun (Wetherill 1978, Hansen 2009) (1 AU is the Earth-Sun distance). Here we report simulations of the early Solar System that show how the inward migration of Jupiter to 1.5 AU, and its subsequent outward migration, lead to a planetesimal disk truncated at 1 AU; the terrestrial planets then form from this disk over the next 30-50 million years, with an Earth/Mars mass ratio consistent with observations. Scattering by Jupiter initially empties but then repopulates the asteroid belt, with inner-belt bodies originating between 1 and 3 AU and outer-belt bodies originating between and beyond the giant planets. This explains the significant compositional differences across the asteroid belt. The key aspect missing from previous models of terrestrial planet formation is the substantial radial migration of the giant planets, which suggests that their behaviour is more similar to that inferred for extrasolar planets than previously thought.

preprint2012arXiv

Building Terrestrial Planets

This paper reviews our current understanding of terrestrial planets formation. The focus is on computer simulations of the dynamical aspects of the accretion process. Throughout the chapter, we combine the results of these theoretical models with geochemical, cosmochemical and chronological constraints, in order to outline a comprehensive scenario of the early evolution of our Solar System. Given that the giant planets formed first in the protoplanetary disk, we stress the sensitive dependence of the terrestrial planet accretion process on the orbital architecture of the giant planets and on their evolution. This suggests a great diversity among the terrestrial planets populations in extrasolar systems. Issues such as the cause for the different masses and accretion timescales between Mars and the Earth and the origin of water (and other volatiles) on our planet are discussed at depth.

preprint2012arXiv

Physical characterization and origin of binary near-Earth asteroid (175706) 1996 FG3

The near-Earth asteroid (NEA) (175706) 1996 FG3 is a particularly interesting spacecraft target: a binary asteroid with a low-DeltaV heliocentric orbit. The orbit of its satellite has provided valuable information about its mass density while its albedo and colors suggest it is primitive or part of the C-complex taxonomic grouping. We extend the physical characterization of this object with new observations of its emission at mid-Infrared (IR) wavelengths and with near-IR reflection spectroscopy. We derive an area-equivalent system diameter of 1.90 \pm 0.28 km (corresponding to approximate component diameters of 1.83 km and 0.51 km, respectively) and a geometric albedo of 0.039 \pm 0.012. 1996 FG3 was previously classified as a C-type asteroid, though the combined 0.4--2.5 micron spectrum with thermal correction indicates classification as B-type; both are consistent with the low measured albedo. Dynamical studies show that 1996 FG3 has most probably originated in the inner main asteroid belt. Recent work has suggested the inner Main Belt (142) Polana family as the possible origin of another low-DeltaV B-type NEA, (101955) 1999 RQ36. A similar origin for 1996 FG3 would require delivery by the overlapping Jupiter 7:2 and Mars 5:9 mean motion resonances rather than the nu-6 resonance, and we find this to be a low probability, but possible, origin.

preprint2011arXiv

Stability analysis of the Martian obliquity during the Noachian era

We performed numerical simulations of the obliquity evolution of Mars during the Noachian era, at which time the giant planets were on drastically different orbits than today. For the preferred primordial configuration of the planets we find that there are two large zones where the Martian obliquity is stable and oscillates with an amplitude lower than 20$^\circ$. These zones occur at obliquities below 30$^\circ$ and above 60$^\circ$; intermediate values show either resonant or chaotic behaviour depending on the primordial orbits of the terrestrial planets.

preprint2011arXiv

The Effect of an Early Planetesimal-Driven Migration of the Giant Planets on Terrestrial Planet Formation

The migration of the giant planets due to the scattering of planetesimals causes powerful resonances to move through the asteroid belt and the terrestrial planet region. Exactly when and how the giant planets migrated is not well known. In this paper we present results of an investigation of the formation of the terrestrial planets during and after the migration of the giant planets. The latter is assumed to have occurred immediately after the dissipation of the nebular disk -- i.e. "early" with respect to the timing of the Late Heavy Bombardment (LHB). The presumed cause of our modeled early migration of the giant planets is angular mometum transfer between the planets and scattered planetesimals.