Researcher profile

Martin J. Duncan

Martin J. Duncan contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

astro-ph.EP astro-ph astro-ph.SR

Trust snapshot

Quick read

Trust 19 - Baseline

5works

0followers

3topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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.

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.

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.

preprint2011arXiv

55 Cancri: A Coplanar Planetary System that is Likely Misaligned with its Star

Although the 55 Cnc system contains multiple, closely packed planets that are presumably in a coplanar configuration, we use numerical simulations to demonstrate that they are likely to be highly inclined to their parent star's spin axis. Due to perturbations from its distant binary companion, this planetary system precesses like a rigid body about its parent star. Consequently, the parent star's spin axis and the planetary orbit normal likely diverged long ago. Because only the projected separation of the binary is known, we study this effect statistically, assuming an isotropic distribution for wide binary orbits. We find that the most likely projected spin-orbit angle is ~50 degrees, with a ~30% chance of a retrograde configuration. Transit observations of the innermost planet - 55 Cnc e - may be used to verify these findings via the Rossiter-McLaughlin effect. 55 Cancri may thus represent a new class of planetary systems with well-ordered, coplanar orbits that are inclined with respect to the stellar equator.

preprint2010arXiv

Eccentricity Trap: Trapping of Resonantly Interacting Planets near the Disk Inner Edge

Using orbital integration and analytical arguments, we have found a new mechanism (an "eccentricity trap") to halt type I migration of planets near the inner edge of a protoplanetary disk. Because asymmetric eccentricity damping due to disk-planet interaction on the innermost planet at the disk edge plays a crucial role in the trap, this mechanism requires continuous eccentricity excitation and hence works for a resonantly interacting convoy of planets. This trap is so strong that the edge torque exerted on the innermost planet can completely halt type I migrations of many outer planets through mutual resonant perturbations. Consequently, the convoy stays outside the disk edge, as a whole. We have derived semi-analytical formula for the condition for the eccentricity trap and predict how many planets are likely to be trapped. We found that several planets or more should be trapped by this mechanism in protoplanetary disks that have cavities. It can be responsible for the formation of non-resonant, multiple, close-in super-Earth systems extending beyond 0.1AU. Such systems are being revealed by radial velocity observations to be quite common around solar-type stars.

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'S data noise characteristics and Cochran et al.'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.