In this paper we extend our numerical method for simulating terrestrial planet formation from Leinhardt and Richardson (2005) to include dynamical friction from the unresolved debris component. In the previous work we implemented a rubble pile planetesimal collision model into direct N-body simulations of terrestrial planet formation. The new collision model treated both accretion and erosion of planetesimals but did not include dynamical friction from debris particles smaller than the resolution limit for the simulation. By extending our numerical model to include dynamical friction from the unresolved debris, we can simulate the dynamical effect of debris produced during collisions and can also investigate the effect of initial debris mass on terrestrial planet formation. We find that significant initial debris mass, 10% or more of the total disk mass, changes the mode of planetesimal growth. Specifically, planetesimals in this situation do not go through a runaway growth phase. Instead they grow concurrently, similar to oligarchic growth. In addition to including the dynamical friction from the unresolved debris, we have implemented particle tracking as a proxy for monitoring compositional mixing. Although there is much less mixing due to collisions and gravitational scattering when dynamical friction of the background debris is included, there is significant inward migration of the largest protoplanets in the most extreme initial conditions.

We report on the results of the first 3D SPH simulation of massive, gravitationally unstable protoplanetary disks with radiative transfer. We adopt a flux-limited diffusion scheme justified by the high opacity of most of the disk. The optically thin surface of the disk cools as a blackbody. The disks grow slowly in mass starting from a Toomre-stable initial condition to the point at which they become marginally unstable. We find that gravitationally bound clumps with masses close to a Jupiter mass can arise. Fragmentation appears to be driven by vertical convective-like motions capable of transporting the heat from the disk midplane to its surface on a timescale of only about 40 years at 10 AU. A larger or smaller cooling efficiency of the disk at the optically thin surface can promote or stifle fragmentation by affecting the vertical temperature profile, which determines whether convection can happen or not, and by regulating the accretion flow from optically thin regions towards overdense regions. We also find that the chances of fragmentation increase for a higher mean molecular weight $\mu$ since compressional heating is reduced. Around a star with mass $1 M_{\odot}$ only disks with $\mu \ge 2.4$, as expected for gas with a metallicity comparable to solar or higher, fragment. This suggests that disk instability, like core-accretion, should be more effective in forming gas giants at higher gas metallicities, consistent with the observed correlation between metallicity of the planet-hosting stars and frequency of planets

We present three-dimensional self-gravitating smoothed-particle hydrodynamics (SPH) simulations of an isothermal gaseous disc interacting with an embedded planet. Discs of varying stability are simulated with planets ranging from 10 Earth-masses to 2 Jupiter-masses. The SPH technique provides the large dynamic range needed to accurately capture the large scale behavior of the disc as well as the small scale interaction of the planet with surrounding material. Most runs used 10^5 gas particles, giving us the spatial resolution required to observe the formation of planets. We find four regions in parameter space: low-mass planets undergo Type I migration; higher-mass planets can form a gap; the gravitational instability mode of planet formation in marginally stable discs can be triggered by embedded planets; discs that are completely unstable can fragment to form many planets. The disc stability is the most important factor in determing which interaction a system will exhibit. For the stable disc cases, our migration and accretion time-scales are shorter and scale differently than previously suggested.

We present three-dimensional SPH calculations of giant planets embedded in gaseous disks. Our findings are collected into a map of parameter space, exhibiting four distinct regions: Type I migration, gap formation, triggered formation of more planets, and wholly unstable disks. For Type I migration of the planet due to secular interactions with the disk material, the migration rate depends linearly on the disk mass, and is independent of the initial planet mass. For more massive disks, the planet can disturb the disk strongly enough to trigger the collapse of gas into additional giant planets. When additional planets form, their interaction as point masses dominates the subsequent behavior of the system. This mechanism allows for the rapid formation of Jupiter-mass and higher planets. Migration due to interaction with the disk can significantly change the orbits of giant planets in gas disks.