Showing 1–50 of 58 results for author: Levison, H F

A Contact Binary Satellite of the Asteroid (152830) Dinkinesh

Authors: Harold F. Levison, Simone Marchi, Keith S. Noll, John R. Spencer, Thomas S. Statler, the Lucy mission team

Abstract: Asteroids with diameters less than about 5 km have complex histories because they are small enough for radiative torques, YORP, to be a notable factor in their evolution. (152830) Dinkinesh is a small asteroid orbiting the Sun near the inner edge of the Main Asteroid Belt with a heliocentric semimajor axis of 2.19 AU; its S type spectrum is typical of bodies in this part of the Main Belt. Here we… ▽ More Asteroids with diameters less than about 5 km have complex histories because they are small enough for radiative torques, YORP, to be a notable factor in their evolution. (152830) Dinkinesh is a small asteroid orbiting the Sun near the inner edge of the Main Asteroid Belt with a heliocentric semimajor axis of 2.19 AU; its S type spectrum is typical of bodies in this part of the Main Belt. Here we report observations by the Lucy spacecraft as it passed within 431 km of Dinkinesh. Lucy revealed Dinkinesh, which has an effective diameter of only $\sim$720 m, to be unexpectedly complex. Of particular note is the presence of a prominent longitudinal trough overlain by a substantial equatorial ridge, and the discovery of the first confirmed contact binary satellite, now named (152830) Dinkinesh I Selam. Selam consists of two near-equal sized lobes with diameters of $\sim$210 m and $\sim$230 m. It orbits Dinkinesh at a distance of 3.1 km with an orbital period of about 52.7 hr, and is tidally locked. The dynamical state, angular momentum, and geomorphologic observations of the system lead us to infer that the ridge and trough of Dinkinesh are probably the result of mass failure resulting from spin-up by YORP followed by the partial reaccretion of the shed material. Selam probably accreted from material shed by this event. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: 25 pages, 5 figures

Implications for the collisional strength of Jupiter Trojans from the Eurybates family

Authors: Raphael Marschall, David Nesvorný, Rogerio Deienno, Ian Wong, Harold F. Levison, William F. Bottke

Abstract: In this work, we model the collisional evolution of the Jupiter Trojans and determined under which conditions the Eurybates-Queta system survives. We show that the collisional strength of the Jupiter Trojans and the age of the Eurybates family and by extension Queta are correlated. The collisional grinding of the Jupiter Trojan population over 4.5 Gy results in a size-frequency distribution (SFD… ▽ More In this work, we model the collisional evolution of the Jupiter Trojans and determined under which conditions the Eurybates-Queta system survives. We show that the collisional strength of the Jupiter Trojans and the age of the Eurybates family and by extension Queta are correlated. The collisional grinding of the Jupiter Trojan population over 4.5 Gy results in a size-frequency distribution (SFD) that remains largely unaltered at large sizes (>10 km) but is depleted at small sizes (10 m to 1 km). This results in a turnover in the SFD, the location of which depends on the collisional strength of the material. It is to be expected that the Trojan SFD bends between 1 and 10 km. Based on the SFD of the Eurybates family, we find that the family was likely the result of a catastrophic impact onto a 100 km rubble pile target. This corresponds to objects with a rather low collisional strength (10 times weaker than that of basaltic material studied in Benz & Asphaug1999). Assuming this weak strength, and an initial cumulative slope of the size frequency distribution of 2.1 between diameters of 2 m and 100 km when the Trojans were captured, the existence of Queta, the satellite of Eurybates, implies an upper limit for the family age of 3.7 Gy. Alternatively, we demonstrate that an unconventional collisional strength with a minimum at 20 m is a plausible candidate to ensure the survival of Queta over the age of the Solar System. Finally, we show how different collisional histories change the expected number of craters on the targets of the Lucy mission and that Lucy will be able to differentiate between them. △ Less

Submitted 22 August, 2022; originally announced August 2022.

Comments: 19 pages, 13 figures

The orbit and density of the Jupiter Trojan satellite system Eurybates-Queta

Authors: M. E. Brown, H. F. Levison, K. S. Noll, R. Binzel, M. W. Buie, W. Grundy, S. Marchi, C. B. Olkin, J. Spencer, T. S. Statler, H. Weaver

Abstract: We report observations of the Jupiter Trojan asteroid (3548) Eurybates and its satellite Queta with the Hubble Space Telescope and use these observations to perform an orbital fit to the system. Queta orbits Eurybates with a semimajor axis of $2350\pm11$ km at a period of $82.46\pm0.06$ days and an eccentricity of $0.125\pm0.009$. From this orbit we derive a mass of Eurybates of… ▽ More We report observations of the Jupiter Trojan asteroid (3548) Eurybates and its satellite Queta with the Hubble Space Telescope and use these observations to perform an orbital fit to the system. Queta orbits Eurybates with a semimajor axis of $2350\pm11$ km at a period of $82.46\pm0.06$ days and an eccentricity of $0.125\pm0.009$. From this orbit we derive a mass of Eurybates of $1.51\pm0.03 \times 10^{17}$ kg, corresponding to an estimated density of $1.1\pm0.3$ g cm$^{-3}$, broadly consistent with densities measured for other Trojans, C-type asteroids in the outer main asteroid belt, and small icy objects from the Kuiper belt. Eurybates is the parent body of the only major collisional family among the Jupiter Trojans; its low density suggests that it is a typical member of the Trojan population. Detailed study of this system in 2027 with the Lucy spacecraft flyby should allow significant insight into collisional processes among what appear to be the icy bodies of the Trojan belt. △ Less

Submitted 3 June, 2021; originally announced June 2021.

Comments: Planetary Science Journal, in press

Lucy Mission to the Trojan Asteroids: Instrumentation and Encounter Concept of Operations

Authors: Catherine B. Olkin, Harold F. Levison, Michael Vincent, Keith S. Noll, John Andrews, Sheila Gray, Phil Good, Simone Marchi, Phil Christensen, Dennis Reuter, Harold Weaver, Martin Patzold, James F. Bell III, Victoria E. Hamilton, Neil Dello Russo, Amy Simon, Matt Beasley, Will Grundy, Carly Howett, John Spencer, Michael Ravine, Michael Caplinger

Abstract: The Lucy Mission accomplishes its science during a series of five flyby encounters with seven Trojan asteroid targets. This mission architecture drives a concept of operations design that maximizes science return, provides redundancy in observations where possible, features autonomous fault protection and utilizes onboard target tracking near closest approach. These design considerations reduce ri… ▽ More The Lucy Mission accomplishes its science during a series of five flyby encounters with seven Trojan asteroid targets. This mission architecture drives a concept of operations design that maximizes science return, provides redundancy in observations where possible, features autonomous fault protection and utilizes onboard target tracking near closest approach. These design considerations reduce risk during the relatively short time-critical periods when science data is collected. The payload suite consists of a color camera and infrared imaging spectrometer, a high-resolution panchromatic imager, and a thermal infrared spectrometer. The mission design allows for concurrent observations of all instruments. Additionally, two spacecraft subsystems will also contribute to the science investigations: the Terminal Tracking Cameras will obtain wide field-of-view imaging near closest approach to determine the shape of each of the Trojan targets and the telecommunication subsystem will carry out Doppler tracking of the spacecraft to determine the mass of each of the Trojan targets. △ Less

Submitted 9 April, 2021; originally announced April 2021.

Convex Shape and Rotation Model of Lucy Target (11351) Leucus from Lightcurves and Occultations

Authors: Stefano Mottola, Stephan Hellmich, Marc W. Buie, Amanda M. Zangari, Simone Marchi, Michael E. Brown, Harold F. Levison

Abstract: We report new photometric lightcurve observations of the Lucy Mission target (11351) Leucus acquired during the 2017, 2018 and 2019 apparitions. We use these data in combination with stellar occultations captured during five epochs (Buie et al. 2020) to determine the sidereal rotation period, the spin axis orientation, a convex shape model, the absolute scale of the object, its geometric albedo, a… ▽ More We report new photometric lightcurve observations of the Lucy Mission target (11351) Leucus acquired during the 2017, 2018 and 2019 apparitions. We use these data in combination with stellar occultations captured during five epochs (Buie et al. 2020) to determine the sidereal rotation period, the spin axis orientation, a convex shape model, the absolute scale of the object, its geometric albedo, and a model of the photometric properties of the target. We find that Leucus is a prograde rotator with a spin axis located within a sky-projected radius of 3° (1$σ$) from J2000 Ecliptic coordinates ($λ=208°$, $β=+77°$) or J2000 Equatorial Coordinates (RA=248$°$, Dec=+58$°$). The sidereal period is refined to $P_{sid}=445.683\pm0.007$ h. The convex shape model is irregular, with maximum dimensions of (60.8, 39.1, 27.8) km. The convex model accounts for global features of the occultation silhouettes, although minor deviations suggest that local and global concavities are present. We determine a geometric albedo $p_V=0.043\pm0.002$. The derived phase curve supports a D-type classification for Leucus. △ Less

Submitted 18 September, 2020; originally announced September 2020.

Detection of a Satellite of the Trojan Asteroid (3548) Eurybates -- A Lucy Mission Target

Authors: K. S. Noll, M. E. Brown, H. A. Weaver, W. M. Grundy, S. B. Porter, M. W. Buie, H. F. Levison, C. Olkin, J. R. Spencer, S. Marchi, T. S. Statler

Abstract: We describe the discovery of a satellite of the Trojan asteroid (3548) Eurybates in images obtained with the Hubble Space Telescope. The satellite was detected on three separate epochs, two in September 2018 and one in January 2020. The satellite has a brightness in all three epochs consistent with an effective diameter of d2 =1.2+/-0.4 km. The projected separation from Eurybates was s~1700-2300 k… ▽ More We describe the discovery of a satellite of the Trojan asteroid (3548) Eurybates in images obtained with the Hubble Space Telescope. The satellite was detected on three separate epochs, two in September 2018 and one in January 2020. The satellite has a brightness in all three epochs consistent with an effective diameter of d2 =1.2+/-0.4 km. The projected separation from Eurybates was s~1700-2300 km and varied in position, consistent with a large range of possible orbits. Eurybates is a target of the Lucy Discovery mission and the early detection of a satellite provides an opportunity for a significant expansion of the scientific return from this encounter. △ Less

Submitted 26 August, 2020; v1 submitted 4 August, 2020; originally announced August 2020.

Comments: 18 pages, 4 figures

Very Slow Rotators from Tidally Synchronized Binaries

Authors: D. Nesvorny, D. Vokrouhlicky, W. F. Bottke, H. F. Levison, W. M. Grundy

Abstract: A recent examination of K2 lightcurves indicates that ~15% of Jupiter Trojans have very slow rotation (spin periods Ps>100 h). Here we consider the possibility that these bodies formed as equal-size binaries in the massive outer disk at ~20-30 au. Prior to their implantation as Jupiter Trojans, tight binaries tidally evolved toward a synchronous state with Ps~Pb, where Pb is the binary orbit perio… ▽ More A recent examination of K2 lightcurves indicates that ~15% of Jupiter Trojans have very slow rotation (spin periods Ps>100 h). Here we consider the possibility that these bodies formed as equal-size binaries in the massive outer disk at ~20-30 au. Prior to their implantation as Jupiter Trojans, tight binaries tidally evolved toward a synchronous state with Ps~Pb, where Pb is the binary orbit period. They may have been subsequently dissociated by impacts and planetary encounters with at least one binary component retaining its slow rotation. Surviving binaries on Trojan orbits would continue to evolve by tides and spin-changing impacts over 4.5 Gyr. To explain the observed fraction of slow rotators, we find that at least ~15-20% of outer disk bodies with diameters 15<D<50 km would have to form as equal-size binaries with 12<ab/R<30, where ab is the binary semimajor axis and R=D/2. The mechanism proposed here could also explain very slow rotators found in other small body populations. △ Less

Submitted 31 March, 2020; v1 submitted 27 March, 2020; originally announced March 2020.

Comments: accepted in ApJL

Planetesimals to Terrestrial Planets: collisional evolution amidst a dissipating gas disk

Authors: Kevin J. Walsh, Harold F. Levison

Abstract: We present numerical simulations of terrestrial planet formation that examine the growth continuously from planetesimals to planets in the inner Solar System. Previous studies show that the growth will be inside-out, but it is still common practice to assume that the entire inner disk will eventually reach a bi-modal distribution of embryos and planetesimals. For the combinations of disk mass, ini… ▽ More We present numerical simulations of terrestrial planet formation that examine the growth continuously from planetesimals to planets in the inner Solar System. Previous studies show that the growth will be inside-out, but it is still common practice to assume that the entire inner disk will eventually reach a bi-modal distribution of embryos and planetesimals. For the combinations of disk mass, initial planetesimal radius and gas disk lifetime explored in this work the entire disk never reaches a simple bi-modal mass distribution. We find that the inside-out growth is amplified by the combined effects of collisional evolution of solid bodies and interactions with a dissipating gas disk. This leads to oligarchic growth never being achieved in different places of the disk at the same time, where in some cases the disk can simultaneoulsy support chaotic growth and giant impacts inside 1 au and runaway growth beyond 2 au. The planetesimal population is efficiently depleted in the inner disk where embryo growth primarily advances in the presence of a significant gas disk. Further out in the disk growth is slower relative to the gas disk dissipation, resulting in more excited planetesimals at the same stage of growth and less efficient accretion. This same effect drives mass loss due to collisional grinding strongly altering the surface density of the accreted planets relative to the initial mass distribution. This effect decreases the Mars-to-Earth mass ratios compared to previous works with no collisional grinding. Similar to some previous findings utilizing vastly different growth scenarios these simulations produce a first generation of planetary embryos that are stable for 10-20 Myr, or 5-10 e-folding times of the gas dissipation timescale, before having an instability and entering the chaotic growth stage. △ Less

Submitted 2 August, 2019; originally announced August 2019.

Comments: 17 figures

Journal ref: Icarus 329,1 (2019)

The Upcoming Mutual Event Season for the Patroclus - Menoetius Trojan Binary

Authors: W. M. Grundy, K. S. Noll, M. W. Buie, H. F. Levison

Abstract: We present new Hubble Space Telescope and ground-based Keck observations and new Keplerian orbit solutions for the mutual orbit of binary Jupiter Trojan asteroid (617) Patroclus and Menoetius, targets of NASA's Lucy mission. We predict event times for the upcoming mutual event season, which is anticipated to run from late 2017 through mid 2019. We present new Hubble Space Telescope and ground-based Keck observations and new Keplerian orbit solutions for the mutual orbit of binary Jupiter Trojan asteroid (617) Patroclus and Menoetius, targets of NASA's Lucy mission. We predict event times for the upcoming mutual event season, which is anticipated to run from late 2017 through mid 2019. △ Less

Submitted 8 March, 2019; originally announced March 2019.

Journal ref: Published in 2018 in Icarus 305, 198-202

Evidence for Very Early Migration of the Solar System Planets from the Patroclus-Menoetius binary Jupiter Trojan

Authors: David Nesvorny, David Vokrouhlicky, William F. Bottke, Harold F. Levison

Abstract: The orbital distribution of trans-Neptunian objects provides strong evidence for the radial migration of Neptune. The outer planets' orbits are thought to have become unstable during the early stages with Jupiter having scattering encounters with a Neptune-class planet. As a consequence, Jupiter jumped inward by a fraction of an au, as required from inner solar system constraints, and obtained its… ▽ More The orbital distribution of trans-Neptunian objects provides strong evidence for the radial migration of Neptune. The outer planets' orbits are thought to have become unstable during the early stages with Jupiter having scattering encounters with a Neptune-class planet. As a consequence, Jupiter jumped inward by a fraction of an au, as required from inner solar system constraints, and obtained its current orbital eccentricity. The timing of these events is often linked to the lunar Late Heavy Bombardment that ended ~700 Myr after the dispersal of the protosolar nebula (t_0). Here we show instead that planetary migration started shortly after t_0. Such early migration is inferred from the survival of the Patroclus-Menoetius binary Jupiter Trojan. The binary formed at t <~ t_0 within a massive planetesimal disk once located beyond Neptune. The longer the binary stayed in the disk, the greater the likelihood that collisions would strip its components from one another. The simulations of its survival indicate that the disk had to have been dispersed by migrating planets within <100 Myr of t_0. This constraint implies that the planetary migration is unrelated to the formation of the youngest lunar basins. △ Less

Submitted 11 September, 2018; originally announced September 2018.

Comments: published in Nature Astronomy

Opportunities for the Large Synoptic Survey Telescope to Find New L$_5$ Trojan and Hilda Lucy Encounter Targets

Authors: Megan E. Schwamb, Harold F. Levison, Marc W. Buie

Abstract: The Trojan clouds contain a surprisingly diverse population of objects. It is only by understanding this diversity can we unravel their history. The more objects Lucy visits the better we can constrain the formation and evolution of the giant planets. While Lucy will visit four Trojans in the L$_4$ cloud, its trajectory only passes one known Trojan system, the PM binary, in L$_5$. There is ample t… ▽ More The Trojan clouds contain a surprisingly diverse population of objects. It is only by understanding this diversity can we unravel their history. The more objects Lucy visits the better we can constrain the formation and evolution of the giant planets. While Lucy will visit four Trojans in the L$_4$ cloud, its trajectory only passes one known Trojan system, the PM binary, in L$_5$. There is ample time in the Lucy timeline to possibly add additional L$_5$ targets, thereby maximizing the mission's science return. In the next decade, the 8.4-m Large Synoptic Survey Telescope (LSST), is expected to detect nearly 300,000 Jupiter Trojans between approximately 16 and 24.5 r-band magnitudes in the $\sim$18,000 deg$^2$ main survey and the expected $\sim$4,000 deg$^2$ northern ecliptic extension. We explore potential opportunities for LSST to find new Lucy L$_5$ fly-by targets. △ Less

Submitted 5 September, 2018; v1 submitted 29 August, 2018; originally announced August 2018.

Comments: accepted to Research Notes of the AAS (RNAAS) - animated figure available in the online RNAAS published version (http://iopscience.iop.org/article/10.3847/2515-5172/aade00)

Origin and Evolution of Short-Period Comets

Authors: D. Nesvorny, D. Vokrouhlicky, L. Dones, H. F. Levison, N. Kaib, A. Morbidelli

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

Submitted 22 June, 2017; originally announced June 2017.

Terrestrial Planet Formation from an Annulus

Authors: Kevin J. Walsh, Hal F. Levison

Abstract: 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… ▽ More 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. △ Less

Submitted 21 September, 2016; originally announced September 2016.

Comments: 23 pages, 15 figs

Journal ref: Astronomical Journal, 152, 3 (2016)

Growing the terrestrial planets from the gradual accumulation of sub-meter sized objects

Authors: Harold F. Levison, Katherine A. Kretke, Kevin Walsh, William Bottke

Abstract: Building the terrestrial planets has been a challenge for planet formation models. In particular, classical theories have been unable to reproduce the small mass of Mars and instead predict that a planet near 1.5 AU should roughly be the same mass as the Earth. Recently, a new model called Viscous Stirred Pebble Accretion (VSPA) has been developed that can explain the formation of the gas giants.… ▽ More Building the terrestrial planets has been a challenge for planet formation models. In particular, classical theories have been unable to reproduce the small mass of Mars and instead predict that a planet near 1.5 AU should roughly be the same mass as the Earth. Recently, a new model called Viscous Stirred Pebble Accretion (VSPA) has been developed that can explain the formation of the gas giants. This model envisions that the cores of the giant planets formed from 100 to 1000 km bodies that directly accreted a population of pebbles --- sub-meter sized objects that slowly grew in the protoplanetary disk. Here we apply this model to the terrestrial planet region and find that it can reproduce the basic structure of the inner Solar System, including a small Mars and a low-mass asteroid belt. Our models show that for an initial population of planetesimals with sizes similar to those of the main belt asteroids, VSPA becomes inefficient beyond $\sim\!$1.5 AU. As a result, Mars's growth is stunted and nothing large in the asteroid belt can accumulate. △ Less

Submitted 7 October, 2015; originally announced October 2015.

Comments: Accepted for PNAS

Growing the gas-giant planets by the gradual accumulation of pebbles

Authors: Harold F. Levison, Katherine A. Kretke, Martin J. Duncan

Abstract: It is widely held that the first step in forming the gas giant planets, such as Jupiter and Saturn, is to form solid `cores' of roughly 10 M$_\oplus$. Getting the cores to form before the solar nebula dissipates ($\sim\!1-10\,$Myr) has been a major challenge for planet formation models. Recently models have emerged in which `pebbles' (centimeter- to meter-size objects) are first concentrated by ae… ▽ More It is widely held that the first step in forming the gas giant planets, such as Jupiter and Saturn, is to form solid `cores' of roughly 10 M$_\oplus$. Getting the cores to form before the solar nebula dissipates ($\sim\!1-10\,$Myr) has been a major challenge for planet formation models. Recently models have emerged in which `pebbles' (centimeter- to meter-size objects) are first concentrated by aerodynamic drag and then gravitationally collapse to form 100 --- 1000 km objects. These `planetesimals' can then efficiently accrete leftover pebbles and directly form the cores of giant planets. This model known as `pebble accretion', theoretically, can produce 10 M$_\oplus$ cores in only a few thousand years. Unfortunately, full simulations of this process show that, rather than creating a few 10 M$_\oplus$ cores, it produces a population of hundreds of Earth-mass objects that are inconsistent with the structure of the Solar System. Here we report that this difficulty can be overcome if pebbles form slowly enough to allow the planetesimals to gravitationally interact with one another. In this situation the largest planetesimals have time to scatter their smaller siblings out of the disk of pebbles, thereby stifling their growth. Our models show that, for a large, and physically reasonable region of parameter space, this typically leads to the formation of one to four gas giants between 5 and 15 AU in agreement with the observed structure of the Solar System. △ Less

Submitted 7 October, 2015; originally announced October 2015.

Journal ref: Nature, Volume 524, Issue 7565, pp. 322-324 (2015)

Evidence for Pebbles in Comets

Authors: K. A. Kretke, H. F. Levison

Abstract: When the EPOXI spacecraft flew by Comet 103P/Hartley 2, it observed large particles floating around the comet nucleus. These particles are likely low-density, centimeter- to decimeter-sized clumps of ice and dust. While the origin of these objects remains somewhat mysterious, it is possible that they are giving us important information about the earliest stages of our Solar System's formation. Rec… ▽ More When the EPOXI spacecraft flew by Comet 103P/Hartley 2, it observed large particles floating around the comet nucleus. These particles are likely low-density, centimeter- to decimeter-sized clumps of ice and dust. While the origin of these objects remains somewhat mysterious, it is possible that they are giving us important information about the earliest stages of our Solar System's formation. Recent advancements in planet formation theory suggest that planetesimals (or cometestimals) may grow directly from the gravitational collapse of aerodynamically concentrated small particles, often referred to as "pebbles." Here we show that the particles observed in the coma of 103P are consistent with the sizes of pebbles expected to efficiently form planetesimals in the region that this comet likely formed, while smaller pebbles are may be expected in the majority of comets, whose chemistry is often indicative of formation in the colder, outer regions of the protoplanetary disk. △ Less

Submitted 2 September, 2015; originally announced September 2015.

Journal ref: Icarus, Volume 262, December 2015, pp 9-13

Formation and evolution of Pluto's small satellites

Authors: Kevin J. Walsh, Harold F. Levison

Abstract: 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 th… ▽ More 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. △ Less

Submitted 5 May, 2015; originally announced May 2015.

Comments: 18 pages, 8 figures. Accepted in AJ

The Mutual Orbit, Mass, and Density of the Large Transneptunian Binary System Varda and Ilmarë

Authors: W. M. Grundy, S. B. Porter, S. D. Benecchi, H. G. Roe, K. S. Noll, C. A. Trujillo, A. Thirouin, J. A. Stansberry, E. Barker, H. F. Levison

Abstract: From observations by the Hubble Space Telescope, Keck II Telescope, and Gemini North Telescope, we have determined the mutual orbit of the large transneptunian object (174567) Varda and its satellite Ilmarë. These two objects orbit one another in a highly inclined, circular or near-circular orbit with a period of 5.75 days and a semimajor axis of 4810 km. This orbit reveals the system mass to be (… ▽ More From observations by the Hubble Space Telescope, Keck II Telescope, and Gemini North Telescope, we have determined the mutual orbit of the large transneptunian object (174567) Varda and its satellite Ilmarë. These two objects orbit one another in a highly inclined, circular or near-circular orbit with a period of 5.75 days and a semimajor axis of 4810 km. This orbit reveals the system mass to be (2.664 +/- 0.064) x 10^20 kg, slightly greater than the mass of the second most massive main-belt asteroid (4) Vesta. The dynamical mass can in turn be combined with estimates of the surface area of the system from Herschel Space Telescope thermal observations to estimate a bulk density of 1.24 +0.50 -0.35 g cm^-3. Varda and Ilmarë both have colors similar to the combined colors of the system, B-V = 0.886 +/- 0.025 and V-I = 1.156 +/- 0.029. △ Less

Submitted 3 May, 2015; originally announced May 2015.

Comments: 16 pages, 4 tables, 6 figures, in press in Icarus

Challenges in Forming the Solar System's Giant Planet Cores via Pebble Accretion

Authors: K. A. Kretke, H. F. Levison

Abstract: Though ~10 Earth mass rocky/icy cores are commonly held as a prerequisite for the formation of gas giants, theoretical models still struggle to explain how these embryos can form within the lifetimes of gaseous circumstellar disks. In recent years, aerodynamic-aided accretion of "pebbles," objects ranging from centimeters to meters in size, has been suggested as a potential solution to this long-s… ▽ More Though ~10 Earth mass rocky/icy cores are commonly held as a prerequisite for the formation of gas giants, theoretical models still struggle to explain how these embryos can form within the lifetimes of gaseous circumstellar disks. In recent years, aerodynamic-aided accretion of "pebbles," objects ranging from centimeters to meters in size, has been suggested as a potential solution to this long-standing problem. While pebble accretion has been demonstrated to be extremely effective in local simulations that look at the detailed behavior of these pebbles in the vicinity of a single planetary embryo, to date there have been no global simulations demonstrating the effectiveness of pebble accretion in a more complicated, multi-planet environment. Therefore, we have incorporated the aerodynamic-aided accretion physics into LIPAD, a Lagrangian code that can follow the collisional / accretional / dynamical evolution of a protoplanetary system, to investigate how pebble accretion manifests itself in the larger planet formation picture. We find that under generic circumstances, pebble accretion naturally leads to an "oligarchic" type of growth in which a large number of planetesimals grow to similar-sized planets. In particular, our simulations tend to form hundreds of Mars- and Earth-mass objects between 4 and 10 AU. While the merging of some oligarchs may allow them to grow massive enough to form giant planet cores, leftover oligarchs lead to planetary systems that cannot be consistent with our own solar system. We investigate various ideas presented in the literature (including evaporation fronts and planet traps) and find that none easily overcome this tendency toward oligarchic growth. △ Less

Submitted 13 October, 2014; v1 submitted 15 September, 2014; originally announced September 2014.

Comments: 20 pages, 33 figures, accepted in AJ

A Lagrangian Integrator for Planetary Accretion and Dynamics (LIPAD)

Authors: Harold F. Levison, Martin J. Duncan, Edward Thommes

Abstract: We presented the first particle based, Lagrangian code that can follow the collisional/accretional/dynamical evolution of a large number of km-sized planetesimals through the entire growth process to become planets. We refer to it as the 'Lagrangian Integrator for Planetary Accretion and Dynamics' or LIPAD. LIPAD is built on top of SyMBA, which is a symplectic $N$-body integrator. In order to hand… ▽ More We presented the first particle based, Lagrangian code that can follow the collisional/accretional/dynamical evolution of a large number of km-sized planetesimals through the entire growth process to become planets. We refer to it as the 'Lagrangian Integrator for Planetary Accretion and Dynamics' or LIPAD. LIPAD is built on top of SyMBA, which is a symplectic $N$-body integrator. In order to handle the very large number of planetesimals required by planet formation simulations, we introduce the concept of a `tracer' particle. Each tracer is intended to represent a large number of disk particles on roughly the same orbit and size as one another, and is characterized by three numbers: the physical radius, the bulk density, and the total mass of the disk particles represented by the tracer. We developed statistical algorithms that follow the dynamical and collisional evolution of the tracers due to the presence of one another. The tracers mainly dynamically interact with the larger objects (`planetary embryos') in the normal N-body way. LIPAD's greatest strength is that it can accurately model the wholesale redistribution of planetesimals due to gravitational interaction with the embryos, which has recently been shown to significantly affect the growth rate of planetary embryos . We verify the code via a comprehensive set of tests which compare our results with those of Eulerian and/or direct N-body codes. △ Less

Submitted 3 July, 2012; originally announced July 2012.

Comments: Accepted to the Astronomical Journal. See http://www.boulder.swri.edu/~hal/LIPAD.html for more detail including animations

Mutual Events in the Cold Classical Transneptunian Binary System Sila and Nunam

Authors: W. M. Grundy, S. D. Benecchi, D. L. Rabinowitz, S. B. Porter, L. H. Wasserman, B. A. Skiff, K. S. Noll, A. J. Verbiscer, M. W. Buie, S. W. Tourtellotte, D. C. Stephens, H. F. Levison

Abstract: Hubble Space Telescope observations between 2001 and 2010 resolved the binary components of the Cold Classical transneptunian object (79360) Sila-Nunam (provisionally designated 1997 CS29). From these observations we have determined the circular, retrograde mutual orbit of Nunam relative to Sila with a period of 12.50995 \pm 0.00036 days and a semimajor axis of 2777 \pm 19 km. A multi-year season… ▽ More Hubble Space Telescope observations between 2001 and 2010 resolved the binary components of the Cold Classical transneptunian object (79360) Sila-Nunam (provisionally designated 1997 CS29). From these observations we have determined the circular, retrograde mutual orbit of Nunam relative to Sila with a period of 12.50995 \pm 0.00036 days and a semimajor axis of 2777 \pm 19 km. A multi-year season of mutual events, in which the two near-equal brightness bodies alternate in passing in front of one another as seen from Earth, is in progress right now, and on 2011 Feb. 1 UT, one such event was observed from two different telescopes. The mutual event season offers a rich opportunity to learn much more about this barely-resolvable binary system, potentially including component sizes, colors, shapes, and albedo patterns. The low eccentricity of the orbit and a photometric lightcurve that appears to coincide with the orbital period are consistent with a system that is tidally locked and synchronized, like the Pluto-Charon system. The orbital period and semimajor axis imply a system mass of (10.84 \pm 0.22) \times 10^18 kg, which can be combined with a size estimate based on Spitzer and Herschel thermal infrared observations to infer an average bulk density of 0.72 +0.37 -0.23 g cm^-3, comparable to the very low bulk densities estimated for small transneptunian binaries of other dynamical classes. △ Less

Submitted 17 April, 2012; originally announced April 2012.

Comments: In press in Icarus

A Method to Constrain the Size of the Protosolar Nebula

Authors: K. A. Kretke, H. F. Levison, M. W. Buie, A. Morbidelli

Abstract: Observations indicate that the gaseous circumstellar disks around young stars vary significantly in size, ranging from tens to thousands of AU. Models of planet formation depend critically upon the properties of these primordial disks, yet in general it is impossible to connect an existing planetary system with an observed disk. We present a method by which we can constrain the size of our own pro… ▽ More Observations indicate that the gaseous circumstellar disks around young stars vary significantly in size, ranging from tens to thousands of AU. Models of planet formation depend critically upon the properties of these primordial disks, yet in general it is impossible to connect an existing planetary system with an observed disk. We present a method by which we can constrain the size of our own protosolar nebula using the properties of the small body reservoirs in the solar system. In standard planet formation theory, after Jupiter and Saturn formed they scattered a significant number of remnant planetesimals into highly eccentric orbits. In this paper, we show that if there had been a massive, extended protoplanetary disk at that time, then the disk would have excited Kozai oscillations in some of the scattered objects, driving them into high-inclination (i > 50 deg), low-eccentricity orbits (q > 30 AU). The dissipation of the gaseous disk would strand a subset of objects in these high-inclination orbits; orbits that are stable on Gyr time scales. To date, surveys have not detected any Kuiper Belt Objects with orbits consistent with this dynamical mechanism. Using these non-detections by the Deep Ecliptic Survey (DES) and the Palomar Distant Solar System Survey we are able to rule out an extended gaseous protoplanetary disk (R_D > 80 AU) in our solar system at the time of Jupiter's formation. Future deep all sky surveys such as the Large Synoptic Survey Telescope (LSST) will all us to further constrain the size of the protoplanetary disk. △ Less

Submitted 10 February, 2012; originally announced February 2012.

Comments: 10 pages, Accepted to AJ

Reassessing the formation of the inner Oort cloud in an embedded star cluster

Authors: R. Brasser, M. J. Duncan, H. F. Levison, M. E. Schwamb, M. E. Brown

Abstract: We re-examine the formation of the inner Oort comet cloud while the Sun was in its birth cluster with the aid of numerical simulations. This work is a continuation of an earlier study (Brasser et al., 2006) with several substantial modifications. First, the system consisting of stars, planets and comets is treated self-consistently in our N-body simulations, rather than approximating the stellar e… ▽ More We re-examine the formation of the inner Oort comet cloud while the Sun was in its birth cluster with the aid of numerical simulations. This work is a continuation of an earlier study (Brasser et al., 2006) with several substantial modifications. First, the system consisting of stars, planets and comets is treated self-consistently in our N-body simulations, rather than approximating the stellar encounters with the outer Solar System as hyperbolic fly-bys. Second, we have included the expulsion of the cluster gas, a feature that was absent previously. Third, we have used several models for the initial conditions and density profile of the cluster -- either a Hernquist or Plummer potential -- and chose other parameters based on the latest observations of embedded clusters from the literature. {These other parameters result in the stars being on radial orbits and the cluster collapses.} Similar to previous studies, in our simulations the inner Oort cloud is formed from comets being scattered by Jupiter and Saturn and having their pericentres decoupled from the planets by perturbations from the cluster gas and other stars. We find that all inner Oort clouds formed in these clusters have an inner edge ranging from 100 AU to a few hundred AU, and an outer edge at over 100\,000 AU, with little variation in these values for all clusters. All inner Oort clouds formed are consistent with the existence of (90377) Sedna, an inner Oort cloud dwarf planetoid, at the inner edge of the cloud: Sedna tends to be at the innermost 2% for Plummer models, while it is 5% for Hernquist models. We emphasise that the existence of Sedna is a generic outcome. We define a `concentration radius' for the inner Oort cloud and find that its value increases with increasing number of stars in the cluster, ranging from 600 AU to 1500 AU for Hernquist clusters and from 1500 AU to 4000 AU for Plummer clusters... △ Less

Submitted 23 October, 2011; originally announced October 2011.

Comments: Accepted Icarus 2011

Ridge Formation and De-Spinning of Iapetus via an Impact-Generated Satellite

Authors: H. F. Levison, K. J. Walsh, A. C. Barr, L. Dones

Abstract: We present a scenario for building the equatorial ridge and de-spinning Iapetus through an impact-generated disk and satellite. This impact puts debris into orbit, forming a ring inside the Roche limit and a satellite outside. This satellite rapidly pushes the ring material down to the surface of Iapetus, and then itself tidally evolves outward, thereby helping to de-spin Iapetus. This scenario ca… ▽ More We present a scenario for building the equatorial ridge and de-spinning Iapetus through an impact-generated disk and satellite. This impact puts debris into orbit, forming a ring inside the Roche limit and a satellite outside. This satellite rapidly pushes the ring material down to the surface of Iapetus, and then itself tidally evolves outward, thereby helping to de-spin Iapetus. This scenario can de-spin Iapetus an order of magnitude faster than when tides due to Saturn act alone, almost independently of its interior geophysical evolution. Eventually, the satellite is stripped from its orbit by Saturn. The range of satellite and impactor masses required is compatible with the estimated impact history of Iapetus. △ Less

Submitted 9 May, 2011; originally announced May 2011.

Comments: 19 pages, 3 figures; Icarus, in press

Journal ref: Icarus v. 214, p. 773-778, 2011

Five New and Three Improved Mutual Orbits of Transneptunian Binaries

Authors: W. M. Grundy, K. S. Noll, F. Nimmo, H. G. Roe, M. W. Buie, S. B. Porter, S. D. Benecchi, D. C. Stephens, H. F. Levison, J. A. Stansberry

Abstract: We present three improved and five new mutual orbits of transneptunian binary systems (58534) Logos-Zoe, (66652) Borasisi-Pabu, (88611) Teharonhiawako-Sawiskera, (123509) 2000 WK183, (149780) Altjira, 2001 QY297, 2003 QW111, and 2003 QY90 based on Hubble Space Telescope and Keck 2 laser guide star adaptive optics observations. Combining the five new orbit solutions with 17 previously known orbits… ▽ More We present three improved and five new mutual orbits of transneptunian binary systems (58534) Logos-Zoe, (66652) Borasisi-Pabu, (88611) Teharonhiawako-Sawiskera, (123509) 2000 WK183, (149780) Altjira, 2001 QY297, 2003 QW111, and 2003 QY90 based on Hubble Space Telescope and Keck 2 laser guide star adaptive optics observations. Combining the five new orbit solutions with 17 previously known orbits yields a sample of 22 mutual orbits for which the period P, semimajor axis a, and eccentricity e have been determined. These orbits have mutual periods ranging from 5 to over 800 days, semimajor axes ranging from 1,600 to 37,000 km, eccentricities ranging from 0 to 0.8, and system masses ranging from 2 x 10^17 to 2 x 10^22 kg. Based on the relative brightnesses of primaries and secondaries, most of these systems consist of near equal-sized pairs, although a few of the most massive systems are more lopsided. The observed distribution of orbital properties suggests that the most loosely-bound transneptunian binary systems are only found on dynamically cold heliocentric orbits. Of the 22 known binary mutual or-bits, orientation ambiguities are now resolved for 9, of which 7 are prograde and 2 are retro-grade, consistent with a random distribution of orbital orientations, but not with models predicting a strong preference for retrograde orbits. To the extent that other perturbations are not dominant, the binary systems undergo Kozai oscillations of their eccentricities and inclinations with periods of the order of tens of thousands to millions of years, some with strikingly high amplitudes. △ Less

Submitted 14 March, 2011; originally announced March 2011.

Comments: In press in Icarus, 29 pages, 5 figures, 13 tables

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

Authors: David Nesvorny, David Vokrouhlicky, William F. Bottke, Keith Noll, Harold F. Levison

Abstract: 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… ▽ More 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. △ Less

Submitted 28 February, 2011; originally announced February 2011.

Comments: The Astronomical Journal, in press

Planetesimal-driven planet migration in the presence of a gas disk

Authors: Christopher C. Capobianco, Martin Duncan, Harold F. Levison

Abstract: We report here on an extension of a previous study by Kirsh et al. (2009) of planetesimal-driven migration using our N-body code SyMBA (Duncan et al., 1998). The previous work focused on the case of a single planet of mass Mem, immersed in a planetesimal disk with a power-law surface density distribution and Rayleigh distributed eccentricities and inclinations. Typically 10^4-10^5 equal-mass plane… ▽ More We report here on an extension of a previous study by Kirsh et al. (2009) of planetesimal-driven migration using our N-body code SyMBA (Duncan et al., 1998). The previous work focused on the case of a single planet of mass Mem, immersed in a planetesimal disk with a power-law surface density distribution and Rayleigh distributed eccentricities and inclinations. Typically 10^4-10^5 equal-mass planetesimals were used, where the gravitational force (and the back-reaction) on each planetesimal by the Sun and planetwere included, while planetesimal-planetesimal interactions were neglected. The runs reported on here incorporate the dynamical effects of a gas disk, where the Adachi et al. (1976) prescription of aerodynamic gas drag is implemented for all bodies. In some cases the Papaloizou and Larwood (2000) prescription of Type-I migration for the planet are implemented, as well as a mass distribution. In the gas-free cases, rapid planet migration was observed - at a rate independent of the planet's mass - provided the planet's mass was not large compared to the mass in planetesimals capable of entering its Hill sphere. In such cases, both inward and outward migrations can be self-sustaining, but there is a strong propensity for inward migration. When a gas disk is present, aerodynamic drag can substantially modify the dynamics of scattered planetesimals. For sufficiently large or small mono-dispersed planetesimals, the planet typically migrates inward. However, for a range of plausible planetesimal sizes (i.e. 0.5-5.0 km at 5.0 AU in a minimum mass Hayashi disk) outward migration is usually triggered, often accompanied by substantial planetary mass accretion. The origins of this behaviour are explained in terms of a toy model. The effects of including a size distribution and torques associated with Type-I migration are also discussed. △ Less

Submitted 23 September, 2010; originally announced September 2010.

Comments: 37 pages, 17 figures, Accepted for publication in Icarus

Evidence from the asteroid belt for a violent past evolution of Jupiter's orbit

Authors: Alessandro Morbidelli, Ramon Brasser, Rodney Gomes, Harold F. Levison, Kleomenis Tsiganis

Abstract: We use the current orbital structure of large (>50km) asteroids in the main asteroid belt to constrain the evolution of the giant planets when they migrated from their primordial orbits to their current ones. Minton & Malhotra (2009) showed that the orbital distribution of large asteroids in the main belt can be reproduced by an exponentially-decaying migration of the giant planets on a time scale… ▽ More We use the current orbital structure of large (>50km) asteroids in the main asteroid belt to constrain the evolution of the giant planets when they migrated from their primordial orbits to their current ones. Minton & Malhotra (2009) showed that the orbital distribution of large asteroids in the main belt can be reproduced by an exponentially-decaying migration of the giant planets on a time scale of tau ~ 0.5My. However, self-consistent numerical simulations show that the planetesimal-driven migration of the giant planets is inconsistent with an exponential change in their semi major axes on such a short time scale (Hahn & Malhotra, 1999). In fact, the typical time scale is tau > 5My. When giant planet migration on this time scale is applied to the asteroid belt, the resulting orbital distribution is incompatible with the observed one. However, the planet migration can be significantly sped up by planet-planet encounters. Consider an evolution where both Jupiter and Saturn have close encounters with a Neptune-mass planet (presumably Uranus or Neptune themselves) and where this third planet, after being scattered inwards by Saturn, is scattered outwards by Jupiter. This scenario leads to a very rapid increase in the orbital separation between Jupiter and Saturn that we show here to have only mild effects on the structure of asteroid belt. This type of evolution is called a jumping-Jupiter case. Our results suggest that the total mass and dynamical excitation of the asteroid belt before migration were comparable to those currently observed. Moreover, they imply that, before migration, the orbits of Jupiter and Saturn were much less eccentric than the current ones. △ Less

Submitted 8 September, 2010; originally announced September 2010.

Comments: AJ, in press

Modeling the Formation of Giant Planet Cores I: Evaluating Key Processes

Authors: H. F. Levison, E. Thommes, M. J. Duncan

Abstract: One of the most challenging problems we face in our understanding of planet formation is how Jupiter and Saturn could have formed before the the solar nebula dispersed. The most popular model of giant planet formation is the so-called 'core accretion' model. In this model a large planetary embryo formed first, mainly by two-body accretion. This is then followed by a period of inflow of nebular g… ▽ More One of the most challenging problems we face in our understanding of planet formation is how Jupiter and Saturn could have formed before the the solar nebula dispersed. The most popular model of giant planet formation is the so-called 'core accretion' model. In this model a large planetary embryo formed first, mainly by two-body accretion. This is then followed by a period of inflow of nebular gas directly onto the growing planet. The core accretion model has an Achilles heel, namely the very first step. We have undertaken the most comprehensive study of this process to date. In this study we numerically integrate the orbits of a number of planetary embryos embedded in a swarm of planetesimals. In these experiments we have included: 1) aerodynamic gas drag, 2) collisional damping between planetesimals, 3) enhanced embryo cross-sections due to their atmospheres, 4) planetesimal fragmentation, and 5) planetesimal driven migration. We find that the gravitational interaction between the embryos and the planetesimals lead to the wholesale redistribution of material - regions are cleared of material and gaps open near the embryos. Indeed, in 90% of our simulations without fragmentation, the region near that embryos is cleared of planetesimals before much growth can occur. The remaining 10%, however, the embryos undergo a burst of outward migration that significantly increases growth. On timescales of ~100,000 years, the outer embryo can migrate ~6 AU and grow to roughly 30 Earth-masses. We also find that the inclusion of planetesimal fragmentation tends to inhibit growth. △ Less

Submitted 16 December, 2009; originally announced December 2009.

Comments: Accepted to AJ, 62 pages 11 figures

(47171) 1999 TC36, A Transneptunian Triple

Authors: S. D. Benecchi, K. S. Noll, W. M. Grundy, H. F. Levison

Abstract: We present new analysis of HST images of (47171) 1999 TC36 that confirm it as a triple system. Fits to the point-spread function consistently show that the apparent primary is itself composed of two similar-sized components. The two central components, A1 and A2, can be consistently identified in each of nine epochs spread over seven years of time. In each instance the component separation, rang… ▽ More We present new analysis of HST images of (47171) 1999 TC36 that confirm it as a triple system. Fits to the point-spread function consistently show that the apparent primary is itself composed of two similar-sized components. The two central components, A1 and A2, can be consistently identified in each of nine epochs spread over seven years of time. In each instance the component separation, ranging from 0.023+/-0.002 to 0.031+/-0.003 arcsec, is roughly one half of the Hubble Space Telescope's diffraction limit at 606 nm. The orbit of the central pair has a semi-major axis of a~867 km with a period of P~1.9 days. These orbital parameters yield a system mass that is consistent with Msys = 12.75+/-0.06 10^18 kg derived from the orbit of the more distant secondary, component B. The diameters of the three components are dA1= 286(+45,-38) km, dA2= 265(+41,-35 km and dB= 139(+22,-18) km. The relative sizes of these components are more similar than in any other known multiple in the solar system. Taken together, the diameters and system mass yield a bulk density of p=542(+317,-211) kg m^-3. HST Photometry shows that component B is variable with an amplitude of >=0.17+/-0.05 magnitudes. Components A1 and A2 do not show variability larger than 0.08+/-0.03 magnitudes approximately consistent with the orientation of the mutual orbit plane and tidally-distorted equilibrium shapes. The system has high specific angular momentum of J/J'=0.93, comparable to most of the known Transneptunian binaries. △ Less

Submitted 10 December, 2009; originally announced December 2009.

Comments: 16 pages, 8 figures, 6 tables. Accepted to Icarus

How Common are Extrasolar, Late Heavy Bombardments?

Authors: Mark Booth, Mark C. Wyatt, Alessandro Morbidelli, Amaya Moro-Martín, Harold F. Levison

Abstract: The habitability of planets is strongly affected by impacts from comets and asteroids. Indications from the ages of Moon rocks suggest that the inner Solar System experienced an increased rate of impacts roughly 3.8 Gya known as the Late Heavy Bombardment (LHB). Here we develop a model of how the Solar System would have appeared to a distant observer during its history based on the Nice model of… ▽ More The habitability of planets is strongly affected by impacts from comets and asteroids. Indications from the ages of Moon rocks suggest that the inner Solar System experienced an increased rate of impacts roughly 3.8 Gya known as the Late Heavy Bombardment (LHB). Here we develop a model of how the Solar System would have appeared to a distant observer during its history based on the Nice model of Gomes et al. (2005). We compare our results with observed debris discs. We show that the Solar System would have been amongst the brightest of these systems before the LHB. Comparison with the statistics of debris disc evolution shows that such heavy bombardment events must be rare occurring around less than 12% of Sun-like stars. △ Less

Submitted 17 November, 2009; originally announced November 2009.

Comments: 2 pages, 1 figure. Contribution to the conference proceedings for 'Pathways towards habitable planets'

Cometary Origin of the Zodiacal Cloud and Carbonaceous Micrometeorites

Authors: David Nesvorny, Peter Jenniskens, Harold F. Levison, William F. Bottke, David Vokrouhlicky

Abstract: 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… ▽ More 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. △ Less

Submitted 26 October, 2009; v1 submitted 23 September, 2009; originally announced September 2009.

Journal ref: Astrophys.J.713:816-836,2010

Constructing the secular architecture of the solar system II: The terrestrial planets

Authors: Ramon Brasser, Alessandro Morbidelli, Rodney Gomes, Kleomenis Tsiganis, Harold F. Levison

Abstract: We investigate the dynamical evolution of the terrestrial planets during the planetesimal-driven migration of the giant planets. A basic assumption of this work is that giant planet migration occurred after the completion of terrestrial planet formation, such as in the models that link the former to the origin of the Late Heavy Bombardment. The divergent migration of Jupiter and Saturn causes th… ▽ More We investigate the dynamical evolution of the terrestrial planets during the planetesimal-driven migration of the giant planets. A basic assumption of this work is that giant planet migration occurred after the completion of terrestrial planet formation, such as in the models that link the former to the origin of the Late Heavy Bombardment. The divergent migration of Jupiter and Saturn causes the g5 eigenfrequency to cross resonances of the form g5=gk with k ranging from 1 to 4. Consequently these secular resonances cause large-amplitude responses in the eccentricities of the terrestrial planets. We show that the resonances g5=g_4 and g5=g3 do not pose a problem if Jupiter and Saturn have a fast approach and departure from their mutual 2:1 mean motion resonance. On the other hand, the resonance crossings g5=g2 and g5=g1 are more of a concern as they tend to yield a terrestrial system incompatible with the current one. We offer two solutions to this problem. The first uses the fact that a secular resonance crossing can also damp the amplitude of a Fourier mode if the latter is large originally. A second scenario involves a 'jumping Jupiter' in which encounters between an ice giant and Jupiter, without ejection of the former, cause the latter to migrate away from Saturn much faster than if migration is driven solely by encounters with planetesimals. In this case, the g5=g2 and g5=g1 resonances can be jumped over, or occur very briefly. △ Less

Submitted 10 September, 2009; originally announced September 2009.

Comments: Astronomy & Astrophysics (2009) in press

Constructing the secular architecture of the solar system I: The giant planets

Authors: Morbidelli Alessandro, Ramon Brasser, Kleomenis Tsiganis, Rodney Gomes, Harold F. Levison

Abstract: Using numerical simulations, we show that smooth migration of the giant planets through a planetesimal disk leads to an orbital architecture that is inconsistent with the current one: the resulting eccentricities and inclinations of their orbits are too small. The crossing of mutual mean motion resonances by the planets would excite their orbital eccentricities but not their orbital inclinations… ▽ More Using numerical simulations, we show that smooth migration of the giant planets through a planetesimal disk leads to an orbital architecture that is inconsistent with the current one: the resulting eccentricities and inclinations of their orbits are too small. The crossing of mutual mean motion resonances by the planets would excite their orbital eccentricities but not their orbital inclinations. Moreover, the amplitudes of the eigenmodes characterising the current secular evolution of the eccentricities of Jupiter and Saturn would not be reproduced correctly; only one eigenmode is excited by resonance-crossing. We show that, at the very least, encounters between Saturn and one of the ice giants (Uranus or Neptune) need to have occurred, in order to reproduce the current secular properties of the giant planets, in particular the amplitude of the two strongest eigenmodes in the eccentricities of Jupiter and Saturn. △ Less

Submitted 10 September, 2009; originally announced September 2009.

Comments: Astronomy & Astrophysics (2009) in press

Asteroids Were Born Big

Authors: Alessandro Morbidelli, William Bottke, David Nesvorny, Harold F. Levison

Abstract: 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 SF… ▽ More 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. △ Less

Submitted 15 July, 2009; originally announced July 2009.

Comments: Icarus (2009) in press

The History of the Solar System's Debris Disc: Observable Properties of the Kuiper Belt

Authors: Mark Booth, Mark C. Wyatt, Alessandro Morbidelli, Amaya Moro-Martín, Harold F. Levison

Abstract: The Nice model of Gomes et al. (2005) suggests that the migration of the giant planets caused a planetesimal clearing event which led to the Late Heavy Bombardment (LHB) at 880 Myr. Here we investigate the IR emission from the Kuiper belt during the history of the Solar System as described by the Nice model. We describe a method for easily converting the results of n-body planetesimal simulation… ▽ More The Nice model of Gomes et al. (2005) suggests that the migration of the giant planets caused a planetesimal clearing event which led to the Late Heavy Bombardment (LHB) at 880 Myr. Here we investigate the IR emission from the Kuiper belt during the history of the Solar System as described by the Nice model. We describe a method for easily converting the results of n-body planetesimal simulations into observational properties (assuming black-body grains and a single size distribution) and further modify this method to improve its realism (using realistic grain properties and a three-phase size distribution). We compare our results with observed debris discs and evaluate the plausibility of detecting an LHB-like process in extrasolar systems. Recent surveys have shown that 4% of stars exhibit 24 um excess and 16% exhibit 70 um excess. We show that the Solar System would have been amongst the brightest of these systems before the LHB at both 24 and 70 um. We find a significant increase in 24 um emission during the LHB, which rapidly drops off and becomes undetectable within 30 Myr, whereas the 70 um emission remains detectable until 360 Myr after the LHB. Comparison with the statistics of debris disc evolution shows that such depletion events must be rare occurring around less than 12% of Sun-like stars and with this level of incidence we would expect approximately 1 of the 413 Sun-like, field stars so far detected to have a 24 um excess to be currently going through an LHB. We also find that collisional processes are important in the Solar System before the LHB and that parameters for weak Kuiper belt objects are inconsistent with the Nice model interpretation of the LHB. △ Less

Submitted 19 June, 2009; originally announced June 2009.

Comments: 16 pages, 13 figures, accepted for publication in MNRAS

Journal ref: MNRAS 399, 385 (2009)

Detectability of Terrestrial Planets in Multi-Planet Systems: Preliminary Report

Authors: Wesley A. Traub, Charles Beichman, Andrew F. Boden, Alan P. Boss, Stefano Casertano, Joseph Catanzarite, Debra Fischer, Eric. B. Ford, Andrew Gould, Sam Halverson, Andrew Howard, Shigeru Ida, N. Jeremy Kasdin, Gregory P. Laughlin, Harold F. Levison, Douglas Lin, Valeri Makarov, James Marr, Matthew Muterspaugh, Sean N. Raymond, Dmitry Savransky, Michael Shao, Alessandro Sozzetti, Cengxing Zhai

Abstract: We ask if Earth-like planets (terrestrial mass and habitable-zone orbit) can be detected in multi-planet systems, using astrometric and radial velocity observations. We report here the preliminary results of double-blind calculations designed to answer this question. We ask if Earth-like planets (terrestrial mass and habitable-zone orbit) can be detected in multi-planet systems, using astrometric and radial velocity observations. We report here the preliminary results of double-blind calculations designed to answer this question. △ Less

Submitted 5 April, 2009; originally announced April 2009.

Comments: 10 pages, 0 figures

Considerations on the magnitude distributions of the Kuiper belt and of the Jupiter Trojans

Authors: Morbidelli Alessandro, Harold F. Levison, William Bottke, Luke Dones, David Nesvorny

Abstract: By examining the absolute magnitude (H) distributions (hereafter HD) of the cold and hot populations in the Kuiper belt and of the Trojans of Jupiter, we find evidence that the Trojans have been captured from the outer part of the primordial trans-Neptunian planetesimal disk. We develop a sketch model of the HDs in the inner and outer parts of the disk that is consistent with the observed distri… ▽ More By examining the absolute magnitude (H) distributions (hereafter HD) of the cold and hot populations in the Kuiper belt and of the Trojans of Jupiter, we find evidence that the Trojans have been captured from the outer part of the primordial trans-Neptunian planetesimal disk. We develop a sketch model of the HDs in the inner and outer parts of the disk that is consistent with the observed distributions and with the dynamical evolution scenario known as the `Nice model'. This leads us to predict that the HD of hot population should have the same slope of the HD of the cold population for 6.5 < H < 9, both as steep as the slope of the Trojans' HD. Current data partially support this prediction, but future observations are needed to clarify this issue. Because the HD of the Trojans rolls over at H~9 to a collisional equilibrium slope that should have been acquired when the Trojans were still embedded in the primordial trans-Neptunian disk, our model implies that the same roll-over should characterize the HDs of the Kuiper belt populations, in agreement with the results of Bernstein et al. (2004) and Fuentes and Holman (2008). Finally, we show that the constraint on the total mass of the primordial trans-Neptunian disk imposed by the Nice model implies that it is unlikely that the cold population formed beyond 35 AU. △ Less

Submitted 5 March, 2009; originally announced March 2009.

Comments: Icarus (2009) in press

A Search for Multi-Planet Systems Using the Hobby-Eberly Telescope

Authors: Robert A. Wittenmyer, Michael Endl, William D. Cochran, Harold F. Levison, Gregory W. Henry

Abstract: Extrasolar multiple-planet systems provide valuable opportunities for testing theories of planet formation and evolution. The architectures of the known multiple-planet systems demonstrate a fascinating level of diversity, which motivates the search for additional examples of such systems in order to better constrain their formation and dynamical histories. Here we describe a comprehensive inves… ▽ More Extrasolar multiple-planet systems provide valuable opportunities for testing theories of planet formation and evolution. The architectures of the known multiple-planet systems demonstrate a fascinating level of diversity, which motivates the search for additional examples of such systems in order to better constrain their formation and dynamical histories. Here we describe a comprehensive investigation of 22 planetary systems in an effort to answer three questions: 1) Are there additional planets? 2) Where could additional planets reside in stable orbits? and 3) What limits can these observations place on such objects? We find no evidence for additional bodies in any of these systems; indeed, these new data do not support three previously announced planets (HD 20367b: Udry et al. 2003, HD 74156d: Bean et al. 2008, and 47 UMa c: Fischer et al. 2002). The dynamical simulations show that nearly all of the 22 systems have large regions in which additional planets could exist in stable orbits. The detection-limit computations indicate that this study is sensitive to close-in Neptune-mass planets for most of the systems targeted. We conclude with a discussion on the implications of these non-detections. △ Less

Submitted 3 March, 2009; originally announced March 2009.

Comments: Accepted to ApJS. Includes 39 pages of radial-velocity data tables

Mutual Orbits and Masses of Six Transneptunian Binaries

Authors: W. M. Grundy, K. S. Noll, M. W. Buie, S. D. Benecchi, D. C. Stephens, H. F. Levison

Abstract: We present Hubble Space Telescope observations of six binary transneptunian systems: 2000 QL251, 2003 TJ58, 2001 XR254, 1999 OJ4, (134860) 2000 OJ67, and 2004 PB108. The mutual orbits of these systems are found to have periods ranging from 22 to 137 days, semimajor axes ranging from 2360 to 10500 km, and eccentricities ranging from 0.09 to 0.55. These orbital parameters enable estimation of syst… ▽ More We present Hubble Space Telescope observations of six binary transneptunian systems: 2000 QL251, 2003 TJ58, 2001 XR254, 1999 OJ4, (134860) 2000 OJ67, and 2004 PB108. The mutual orbits of these systems are found to have periods ranging from 22 to 137 days, semimajor axes ranging from 2360 to 10500 km, and eccentricities ranging from 0.09 to 0.55. These orbital parameters enable estimation of system masses ranging from 0.2 to 9.7 x 10+18 kg. For reasonable assumptions of bulk density (0.5 to 2.0 g cm-3), the masses can be combined with visible photometry to constrain sizes and albedos. The resulting albedos are consistent with an emerging picture of the dynamically "Cold" Classical subpopulation having relatively high albedos, compared with comparablysized objects on more dynamically excited orbits. △ Less

Submitted 7 January, 2009; v1 submitted 16 December, 2008; originally announced December 2008.

Comments: 18 pages, 3 figures, 10 tables; in press in Icarus

The Correlated Colors of Transneptunian Binaries

Authors: S. D. Benecchi, K. S. Noll, W. M. Grundy, M. W. Buie, D. C. Stephens, H. F. Levison

Abstract: We report resolved photometry of the primary and secondary components of 23 transneptunian binaries obtained with the Hubble Space Telescope. V-I colors of the components range from 0.7 to 1.5 with a median uncertainty of 0.06 magnitudes. The colors of the primaries and secondaries are correlated with a Spearman rank correlation probability of 99.99991%, 5 sigma for a normal distribution. Fits t… ▽ More We report resolved photometry of the primary and secondary components of 23 transneptunian binaries obtained with the Hubble Space Telescope. V-I colors of the components range from 0.7 to 1.5 with a median uncertainty of 0.06 magnitudes. The colors of the primaries and secondaries are correlated with a Spearman rank correlation probability of 99.99991%, 5 sigma for a normal distribution. Fits to the primary vs. secondary colors are identical to within measurement uncertainties. The color range of binaries as a group is indistinguishable from that of the larger population of apparently single transneptunian objects. Whatever mechanism produced the colors of apparently single TNOs acted equally on binary systems. The most likely explanation is that the colors of transneptunian objects and binaries alike are primordial and indicative of their origin in a locally homogeneous, globally heterogeneous protoplanetary disk. △ Less

Submitted 13 November, 2008; originally announced November 2008.

Comments: 28 pages, 4 figure, 4 tables. accepted to Icarus

On a Scattered-Disk Origin for the 2003 El61 Collisional Family - an Example of the Importance of Collisions on the Dynamics of Small Bodies

Authors: Harold F. Levison, Morbidelli Alessandro, David Vokrouhlicky, William Bottke

Abstract: The recent discovery of the 2003EL61 collisional family in the Kuiper belt (Brown et al., 2007) is surprising because the formation of such a family is a highly improbable event in today's belt. Assuming Brown et al.'s estimate of the size of progenitors, we find that the probability that a Kuiper belt object was involved in such a collision since primordial times is less than roughly 0.001. In… ▽ More The recent discovery of the 2003EL61 collisional family in the Kuiper belt (Brown et al., 2007) is surprising because the formation of such a family is a highly improbable event in today's belt. Assuming Brown et al.'s estimate of the size of progenitors, we find that the probability that a Kuiper belt object was involved in such a collision since primordial times is less than roughly 0.001. In addition, it is not possible for the collision to have occurred in a massive primordial Kuiper belt because the dynamical coherence of the family would not have survived whatever event produced the currently observed orbital excitation. Here we suggest that the family is the result of a collision between two scattered disk objects. We show that the probability that a collision occurred between two such objects with sizes similar to those advocated in Brown et al. (2007) and that the center of mass of the resulting family is on an orbit typical of the Kuiper belt can be as large as 47%. Given the large uncertainties involved in this estimate, this result is consistent with the existence of one such family. If true, this result has implications far beyond the origin of a single collisional family, because it shows that collisions played an important role in shaping the dynamical structure of the small body populations that we see today. △ Less

Submitted 3 September, 2008; originally announced September 2008.

Comments: The Astronomical Journal (2008) in press

(42355) Typhon-Echidna: Scheduling Observations for Binary Orbit Determination

Authors: W. M. Grundy, K. S. Noll, J. Virtanen, K. Muinonen, S. D. Kern, D. C. Stephens, J. A. Stansberry, H. F. Levison, J. R. Spencer

Abstract: We describe a strategy for scheduling astrometric observations to minimize the number required to determine the mutual orbits of binary transneptunian systems. The method is illustrated by application to Hubble Space Telescope observations of (42355) Typhon-Echidna, revealing that Typhon and Echidna orbit one another with a period of 18.971 +/- 0.006 days and a semimajor axis of 1628 +/- 29 km,… ▽ More We describe a strategy for scheduling astrometric observations to minimize the number required to determine the mutual orbits of binary transneptunian systems. The method is illustrated by application to Hubble Space Telescope observations of (42355) Typhon-Echidna, revealing that Typhon and Echidna orbit one another with a period of 18.971 +/- 0.006 days and a semimajor axis of 1628 +/- 29 km, implying a system mass of (9.49 +/- 0.52) x 10^17 kg. The eccentricity of the orbit is 0.526 +/- 0.015. Combined with a radiometric size determined from Spitzer Space Telescope data and the assumption that Typhon and Echidna both have the same albedo, we estimate that their radii are 76 +14/-16 and 42 +8/-9 km, respectively. These numbers give an average bulk density of only 0.44 +0.44/-0.17 g cm^-3, consistent with very low bulk densities recently reported for two other small transneptunian binaries. △ Less

Submitted 15 April, 2008; originally announced April 2008.

Origin of the Structure of the Kuiper Belt during a Dynamical Instability in the Orbits of Uranus and Neptune

Authors: Harold F. Levison, Alessandro Morbidelli, Christa Van Laerhoven, Rodney Gomes, Kleomenis Tsiganis

Abstract: We explore the origin and orbital evolution of the Kuiper belt in the framework of a recent model of the dynamical evolution of the giant planets, sometimes known as the Nice model. This model is characterized by a short, but violent, instability phase, during which the planets were on large eccentricity orbits. One characteristic of this model is that the proto-planetary disk must have been tru… ▽ More We explore the origin and orbital evolution of the Kuiper belt in the framework of a recent model of the dynamical evolution of the giant planets, sometimes known as the Nice model. This model is characterized by a short, but violent, instability phase, during which the planets were on large eccentricity orbits. One characteristic of this model is that the proto-planetary disk must have been truncated at roughly 30 to 35 AU so that Neptune would stop migrating at its currently observed location. As a result, the Kuiper belt would have initially been empty. In this paper we present a new dynamical mechanism which can deliver objects from the region interior to ~35 AU to the Kuiper belt without excessive inclination excitation. Assuming that the last encounter with Uranus delivered Neptune onto a low-inclination orbit with a semi-major axis of ~27 AU and an eccentricity of ~0.3, and that subsequently Neptune's eccentricity damped in ~1 My, our simulations reproduce the main observed properties of the Kuiper belt at an unprecedented level. △ Less

Submitted 4 December, 2007; originally announced December 2007.

Evidence for Two Populations of Classical Transneptunian Objects: The Strong Inclination Dependence of Classical Binaries

Authors: Keith S. Noll, William M. Grundy, Denise C. Stephens, Harold F. Levison, Susan D. Kern

Abstract: We have searched 101 Classical transneptunian objects for companions with the Hubble Space Telescope. Of these, at least 21 are binary. The heliocentric inclinations of the objects we observed range from 0.6-34 degrees. We find a very strong anticorrelation of binaries with inclination. Of the 58 targets that have inclinations of less than 5.5 degrees, 17 are binary, a binary fraction of 29 +7/-… ▽ More We have searched 101 Classical transneptunian objects for companions with the Hubble Space Telescope. Of these, at least 21 are binary. The heliocentric inclinations of the objects we observed range from 0.6-34 degrees. We find a very strong anticorrelation of binaries with inclination. Of the 58 targets that have inclinations of less than 5.5 degrees, 17 are binary, a binary fraction of 29 +7/-6%. All 17 are similar-brightness systems. On the contrary, only 4 of the 42 objects with inclinations greater than 5.5 degrees have satellites and only 1 of these is a similar-brightness binary. This striking dichotomy appears to agree with other indications that the low eccentricity, non-resonant Classical transneptunian objects include two overlapping populations with significantly different physical properties and dynamical histories. △ Less

Submitted 10 March, 2008; v1 submitted 9 November, 2007; originally announced November 2007.

Comments: 32 pages, 4 figures, 2 tables in press for Icarus volume 194 Note added in proof and reference updates in this corrected version

Dynamical and Observational Constraints on Additional Planets in Highly Eccentric Planetary Systems

Authors: Robert A. Wittenmyer, Michael Endl, William D. Cochran, Harold F. Levison

Abstract: Long time coverage and high radial velocity precision have allowed for the discovery of additional objects in known planetary systems. Many of the extrasolar planets detected have highly eccentric orbits, which raises the question of how likely those systems are to host additional planets. We investigate six systems which contain a very eccentric (e>0.6) planet: HD 3651, HD 37605, HD 45350, HD 8… ▽ More Long time coverage and high radial velocity precision have allowed for the discovery of additional objects in known planetary systems. Many of the extrasolar planets detected have highly eccentric orbits, which raises the question of how likely those systems are to host additional planets. We investigate six systems which contain a very eccentric (e>0.6) planet: HD 3651, HD 37605, HD 45350, HD 80606, HD 89744, and 16 Cyg B. We present updated radial-velocity observations and orbital solutions, search for additional planets, and perform test particle simulations to find regions of dynamical stability. The dynamical simulations show that short-period planets could exist in the HD 45350 and 16 Cyg B systems, and we use the observational data to set tight detection limits, which rule out additional planets down to a few Neptune masses in the HD 3651, HD 45350, and 16 Cyg B systems. △ Less

Submitted 13 June, 2007; originally announced June 2007.

Comments: Accepted for publication in AJ

Journal ref: Astron.J.134:1276-1284,2007

Dynamics of the giant planets of the solar system in the gaseous proto-planetary disk and relationship to the current orbital architecture

Authors: Alessandro Morbidelli, Kleomenis Tsiganis, Aurelien Crida, Harold F. Levison, Rodney Gomes

Abstract: We study the orbital evolution of the 4 giant planets of our solar system in a gas disk. Our investigation extends the previous works by Masset and Snellgrove (2001) and Morbidelli and Crida (2007, MC07), which focussed on the dynamics of the Jupiter-Saturn system. The only systems that we found to reach a steady state are those in which the planets are locked in a quadruple mean motion resonanc… ▽ More We study the orbital evolution of the 4 giant planets of our solar system in a gas disk. Our investigation extends the previous works by Masset and Snellgrove (2001) and Morbidelli and Crida (2007, MC07), which focussed on the dynamics of the Jupiter-Saturn system. The only systems that we found to reach a steady state are those in which the planets are locked in a quadruple mean motion resonance (i.e. each planet is in resonance with its neighbor). In total we found 6 such configurations. For the gas disk parameters found in MC07, these configurations are characterized by a negligible migration rate. After the disappearance of the gas, and in absence of planetesimals, only two of these six configurations (the least compact ones) are stable for a time of hundreds of millions of years or more. The others become unstable on a timescale of a few My. Our preliminary simulations show that, when a planetesimal disk is added beyond the orbit of the outermost planet, the planets can evolve from the most stable of these configurations to their current orbits in a fashion qualitatively similar to that described in Tsiganis et al. (2005). △ Less

Submitted 17 July, 2007; v1 submitted 12 June, 2007; originally announced June 2007.

Comments: The Astronomical Journal (17/07/2007) in press

Journal ref: Astron.J.134:1790-1798,2007

The orbit, mass, size, albedo, and density of (65489) Ceto/Phorcys: A tidally-evolved binary Centaur

Authors: W. M. Grundy, J. A. Stansberry, K. S. Noll, D. C. Stephens, D. E. Trilling, S. D. Kern, J. R. Spencer, D. P. Cruikshank, H. F. Levison

Abstract: Hubble Space Telescope observations of Uranus- and Neptune-crossing object (65489) Ceto/Phorcys (provisionally designated 2003 FX128) reveal it to be a close binary system. The mutual orbit has a period of 9.554 +/- 0.011 days and a semimajor axis of 1840 +/- 48 km. These values enable computation of a system mass of (5.41 +/- 0.42) 10^18 kg. Spitzer Space Telescope observations of thermal emiss… ▽ More Hubble Space Telescope observations of Uranus- and Neptune-crossing object (65489) Ceto/Phorcys (provisionally designated 2003 FX128) reveal it to be a close binary system. The mutual orbit has a period of 9.554 +/- 0.011 days and a semimajor axis of 1840 +/- 48 km. These values enable computation of a system mass of (5.41 +/- 0.42) 10^18 kg. Spitzer Space Telescope observations of thermal emission at 24 and 70 microns are combined with visible photometry to constrain the system's effective radius (109 +10/-11 km) and geometric albedo (0.084 +0.021/-0.014). We estimate the average bulk density to be 1.37 +0.66/-0.32 g cm^-3, consistent with ice plus rocky and/or carbonaceous materials. This density contrasts with lower densities recently measured with the same technique for three other comparably-sized outer Solar System binaries (617) Patroclus, (26308) 1998 SM165, and (47171) 1999 TC36, and is closer to the density of the saturnian irregular satellite Phoebe. The mutual orbit of Ceto and Phorcys is nearly circular, with an eccentricity <= 0.015. This observation is consistent with calculations suggesting that the system should tidally evolve on a timescale shorter than the age of the solar system. △ Less

Submitted 12 April, 2007; originally announced April 2007.

The Dynamical Structure of the Kuiper Belt and its Primordial Origin

Authors: A. Morbidelli, H. F. Levison, R. Gomes

Abstract: This chapter discusses the dynamical properties of the Kuiper belt population. Then, it focuses on the characteristics of the Kuiper belt that cannot be explained by its evolution in the framework of the current solar system. We review models of primordial solar system evolution that have been proposed to reproduce the Kuiper belt features, outlining advantages and problems of each of them. This chapter discusses the dynamical properties of the Kuiper belt population. Then, it focuses on the characteristics of the Kuiper belt that cannot be explained by its evolution in the framework of the current solar system. We review models of primordial solar system evolution that have been proposed to reproduce the Kuiper belt features, outlining advantages and problems of each of them. △ Less

Submitted 21 March, 2007; originally announced March 2007.

Comments: to be published in `Transneptunian Objects', Barucci et al. eds., University of Arizona Press

Models of the Collisional Damping Scenario for Ice Giant Planets and Kuiper Belt Formatio

Authors: Harold F. Levison, Alessandro Morbidelli

Abstract: Chiang et al. 2006, hereafter C06 have recently proposed that the observed structure of the Kuiper belt could be the result of a dynamical instability of a system of ~5 primordial ice giant planets in the outer Solar System. According to this scenario, before the instability occurred, these giants were growing in a highly collisionally damped environment according to the arguments in Goldreich e… ▽ More Chiang et al. 2006, hereafter C06 have recently proposed that the observed structure of the Kuiper belt could be the result of a dynamical instability of a system of ~5 primordial ice giant planets in the outer Solar System. According to this scenario, before the instability occurred, these giants were growing in a highly collisionally damped environment according to the arguments in Goldreich et al. (2004a,b, hereafter G04). Here we test this hypothesis with a series of numerical simulations using a new code designed to incorporate the dynamical effects of collisions. We find that we cannot reproduce the observed Solar System. In particular, G04 and C06 argue that during the instability, all but two of the ice giants would be ejected from the Solar System by Jupiter and Saturn, leaving Uranus and Neptune behind. We find that ejections are actually rare and that instead the systems spread outward. This always leads to a configuration with too many planets that are too far from the Sun. Thus, we conclude that both G04's scheme for the formation of Uranus and Neptune and C06's Kuiper belt formation scenario are not viable in their current forms. △ Less

Submitted 18 January, 2007; originally announced January 2007.