BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Harold F. Levison

Harold F. Levison contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.EPastro-phastro-ph.IM

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

13 published item(s)

preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Formation and evolution of Pluto's small satellites

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Lagrangian Integrator for Planetary Accretion and Dynamics (LIPAD)

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Evidence from the asteroid belt for a violent past evolution of Jupiter&#39;s orbit

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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&#39;s mass - provided the planet&#39;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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Asteroids Were Born Big

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Cometary Origin of the Zodiacal Cloud and Carbonaceous Micrometeorites

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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

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 &#39;jumping Jupiter&#39; 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.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

The Calibration of the HST Kuiper Belt Object Search: Setting the Record Straight

The limiting magnitude of the HST data set used by Cochran et al. (1995) to detect small objects in the Kuiper belt is reevaluated, and the methods used are described in detail. It is shown, by implanting artificial objects in the original HST images, and re-reducing the images using our original algorithm, that the limiting magnitude of our images (as defined by the 50% detectability limit) is $V=28.4$. This value is statistically the same as the value found in the original analysis. We find that $\sim50%$ of the moving Kuiper belt objects with $V=27.9$ are detected when trailing losses are included. In the same data in which these faint objects are detected, we find that the number of false detections brighter than $V=28.8$ is less than one per WFPC2 image. We show that, primarily due to a zero-point calibration error, but partly due to inadequacies in modeling the HST&#39;S data noise characteristics and Cochran et al.&#39;s reduction techniques, Brown et al. 1997 underestimate the SNR of objects in the HST dataset by over a factor of 2, and their conclusions are therefore invalid.

0 citations0 reviewsNo code linkedOpen access