Terrestrial Planet Formation Surrounding Close Binary Stars

Disk material has been observed around both components of some young close binary star systems. It has been shown that if planets form at the right places within such disks, they can remain dynamically stable for very long times. Herein, we numerically simulate the late stages of terrestrial planet growth in circumbinary disks around 'close' binary star systems with stellar separations between 0.05 AU and 0.4 AU and binary eccentricities up to 0.8. In each simulation, the sum of the masses of the two stars is 1 solar mass, and giant planets are included. Our results are statistically compared to a set of planet formation simulations in the Sun-Jupiter-Saturn system that begin with essentially the same initial disk of protoplanets. The planetary systems formed around binaries with apastron distances less than ~ 0.2 AU are very similar to those around single stars, whereas those with larger maximum separations tend to be sparcer, with fewer planets, especially interior to 1 AU. We also provide formulae that can be used to scale results of planetary accretion simulations to various systems with different total stellar mass, disk sizes, and planetesimal masses and densities.Comment: 60 pages, 4 tables, and 11 low resolution eps figures. Article with high resolution figures is available at http://www-personal.umich.edu/~equintan/publications.html . Accepted for publication in Icaru

A scaling law for accretion zone sizes

Current theories of runaway planetary accretion require small random velocities of the accreted particles. Two body gravitational accretion cross sections which ignore tidal perturbations of the Sun are not valid for the slow encounters which occur at low relative velocities. Wetherill and Cox have studied accretion cross sections for rocky protoplanets orbiting at 1 AU. Using analytic methods based on Hill's lunar theory, one can scale these results for protoplanets that occupy the same fraction of their Hill sphere as does a rocky body at 1 AU. Generalization to bodies of different sizes is achieved here by numerical integrations of the three-body problem. Starting at initial positions far from the accreting body, test particles are allowed to encounter the body once, and the cross section is computed. A power law is found relating the cross section to the radius of the accreting body (of fixed mass)

Composition of Massive Giant Planets

The two current models for giant planet formation are core accretion and disk instability. We discuss the core masses and overall planetary enrichment in heavy elements predicted by the two formation models, and show that both models could lead to a large range of final compositions. For example, both can form giant planets with nearly stellar compositions. However, low-mass giant planets, enriched in heavy elements compared to their host stars, are more easily explained by the core accretion model. The final structure of the planets, i.e., the distribution of heavy elements, is not firmly constrained in either formation model.Comment: 6 pages, Proceedings of IAU Symposium 276 (Invited talk), The Astrophysics of Planetary Systems: Formation, Structure, and Dynamical Evolution. Turin, Italy, Oct. 201

Kepler's Multiple Planet Systems

More than one-third of the 4700 planet candidates found by NASAs Kepler spacecraft are associated with target stars that have more than one planet candidate, and such multis account for the vast majority of candidates that have been verified as true planets.The large number of multis tells us that flat multiplanet systems like our Solar System are common. Virtually all of the candidate planetary systems are stable, as tested by numerical integrations that assume a physically motivated mass-radius relationship. Statistical studies performed on these candidate systems reveal a great deal about the architecture of planetary systems, including the typical spacing of orbits and flatness. The characteristics of several of the most interesting confirmed Kepler & K2 multi-planet systems will also be discussed

Long-term evolution of a planetesimal swarm in the vicinity of a protoplanet

Many models of planet formation involve scenarios in which one or a few large protoplanets interact with a swarm of much smaller planetesimals. In such scenarios, three-body perturbations by the protoplanet as well as mutual collisions and gravitational interactions between the swarm bodies are important in determining the velocity distribution of the swarm. We are developing a model to examine the effects of these processes on the evolution of a planetesimal swarm. The model consists of a combination of numerical integrations of the gravitational influence of one (or a few) massive protoplanets on swarm bodies together with a statistical treatment of the interactions between the planetesimals. Integrating the planetesimal orbits allows us to take into account effects that are difficult to model analytically or statistically, such as three-body collision cross-sections and resonant perturbations by the protoplanet, while using a statistical treatment for the particle-particle interactions allows us to use a large enough sample to obtain meaningful results

Formation of Jupiter using opacities based on detailed grain physics

Numerical simulations, based on the core-nucleated accretion model, are presented for the formation of Jupiter at 5.2 AU in 3 primordial disks with three different assumed values of the surface density of solid particles. The grain opacities in the envelope of the protoplanet are computed using a detailed model that includes settling and coagulation of grains and that incorporates a recalculation of the grain size distribution at each point in time and space. We generally find lower opacities than the 2% of interstellar values used in previous calculations [Hubickyj, O., Bodenheimer, P., Lissauer, J. J., 2005. Icarus 179, 415--431; Lissauer, J. J., Hubickyj, O., D'Angelo, G., Bodenheimer, P., 2009. Icarus 199, 338-350]. These lower opacities result in more rapid heat loss from and more rapid contraction of the protoplanetary envelope. For a given surface density of solids, the new calculations result in a substantial speedup in formation time as compared with those previous calculations. Formation times are calculated to be 1.0, 1.9, and 4.0 Myr, and solid core masses are found to be 16.8, 8.9, and 4.7 Earth masses, for solid surface densities, sigma, of 10, 6, and 4 grams per squared centimeter, respectively. For sigma=10 and sigma=6 g/cm^2, respectively, these formation times are reduced by more than 50% and more than 80% compared with those in a previously published calculation with the old approximation to the opacity