Circumstellar dust distribution in systems with two planets in resonance

We investigate via numerical modeling the effects of two planets locked in resonance, and migrating outward, on the dust distribution of the natal circumstellar disk. We aim to test whether the dust distribution exhibits peculiar features arising from the interplay among the gravitational perturbations of the planets in resonance, the evolution of the gas, and its influence on the dust grains' dynamics. We focus on the 3:2 and 2:1 resonance, where the trapping may be caused by the convergent migration of a Jupiter- and Saturn-mass planet, preceding the common gap formation and ensuing outward (or inward) migration. Models show that a common gap also forms in the dust component -- similarly to what a single, more massive planet would generate -- and that outward migration leads to a progressive widening of the dust gap and to a decoupling from the gas gap. As the system evolves, a significantly wider gap is observed in the dust distribution, which ceases to overlap with the gas gap in the inner disk regions. At the outer edge of the gas gap, outward migration of the planets produces an over-density of dust particles, which evolve differently in the 3:2 and 2:1 resonances. For the 3:2, the dust trap at the gap's outer edge is partly efficient and a significant fraction of the grains filters through the gap. For the 2:1 resonance, the trap is more efficient and very few grains cross the gap, while the vast majority accumulate at the outer edge of the gap.Comment: 16 pages, 12 figures. Accepted for publicatio

Dust distribution around low-mass planets on converging orbits

Super-Earths can form at large orbital radii and migrate inward due to tidal interactions with the circumstellar disk. In this scenario, convergent migration may occur and lead to the formation of resonant pairs of planets. We explore the conditions under which convergent migration and resonance capture take place, and what dynamical consequences can be expected on the dust distribution surrounding the resonant pair. We combine hydrodynamic planet--disk interaction models with dust evolution calculations to investigate the signatures produced in the dust distribution by a pair of planets in mean-motion resonances. We find that convergent migration takes place when the outer planet is the more massive. However, convergent migration also depends on the local properties of the disk, and divergent migration may result as well. For similar disk parameters, the capture in low degree resonances (e.g., 2:1 or 3:2) is preferred close to the star where the resonance strength can more easily overcome the tidal torques exerted by the gaseous disk. Farther away from the star, convergent migration may result in capture in high degree resonances. The dust distribution shows potentially observable features typically when the planets are trapped in a 2:1 resonance. In other cases, with higher degree resonances (e.g., 5:4 or 6:5) dust features may not be sufficiently pronounced to be easily observable. The degree of resonance established by a pair of super-Earths may be indicative of the location in the disk where capture occurred. There can be significant differences in the dust distribution around a single super-Earth and a pair of super-Earths in resonance.Comment: Accepted for publication on Astronomy and Astrophysic

New Formation Models for the Kepler-36 System

Formation of the planets in the Kepler-36 system is modeled by detailed numerical simulations according to the core-nucleated accretion scenario. The standard model is updated to include the dissolution of accreting rocky planetesimals in the gaseous envelope of the planet, leading to substantial enrichment of the envelope mass in heavy elements and a non-uniform composition with depth. For Kepler-36 c, models involving in situ formation and models involving orbital migration are considered. The results are compared with standard formation models. The calculations include the formation (accretion) phase, as well as the subsequent cooling phase, up to the age of Kepler-36 (7 Gyr). During the latter phase, mass loss induced by stellar XUV radiation is included. In all cases, the results fit the measured mass, 7.84 M$_\oplus$, and radius, 3.68 R$_\oplus$, of Kepler-36 c. Two parameters are varied to obtain these fits: the disk solid surface density at the formation location, and the "efficiency" factor in the XUV mass loss rate. The updated models are hotter and therefore less dense in the silicate portion of the planet and in the overlying layers of H/He, as compared with standard models. The lower densities mean that only about half as much H/He is needed to be accreted to fit the present-day mass and radius constraints. For Kepler-36 b, an updated in situ calculation shows that the entire H/He envelope is lost, early in the cooling phase, in agreement with observation.Comment: 21 pages, 18 figures, 1 table. Accepted for publication in The Astrophysical Journa

Three-dimensional Disk-Planet Torques in a Locally Isothermal Disk

We determine an expression for the Type I planet migration torque involving a locally isothermal disk, with moderate turbulent viscosity (~0.0005 < alpha < ~0.05), based on three-dimensional nonlinear hydrodynamical simulations. The radial gradients (in a dimensionless logarithmic form) of density and temperature are assumed to be constant near the planet. We find that the torque is roughly equally sensitive to the surface density and temperature radial gradients. Both gradients contribute to inward migration when they are negative. Our results indicate that two-dimensional calculations with a smoothed planet potential, used to account for the effects of the third dimension, do not accurately determine the effects of density and temperature gradients on the three-dimensional torque. The results suggest that substantially slowing or stopping planet migration by means of changes in disk opacity or shadowing is difficult and appears unlikely for a disk that is locally isothermal. The scalings of the torque and torque density with planet mass and gas sound speed follow the expectations of linear theory. We also determine an improved formula for the torque density distribution that can be used in one-dimensional long-term evolution studies of planets embedded in locally isothermal disks. This formula can be also applied in the presence of mildly varying radial gradients and of planets that open gaps. We illustrate its use in the case of migrating super-Earths and determine some conditions sufficient for survival.Comment: 20 pages, 16 figures. To appear in The Astrophysical Journal. Updated with changes and corrections added in proo

Recent Simulations of the Late Stages Growth of Jupiter

Presented by Lissauer et al. (2009, Icarus 199, 338) are used to test the model of capture of Jupiter's irregular satellites within proto-Jupiter's distended and thermally-supported envelope. We find such capture highly unlikely, since the envelope shrinks too slowly for a large number of moons to be retained, and many of those that would be retained would orbit closer to the planet than do the observed Jovian irregulars. Our calculations do not address (and therefore do not exclude) the possibility that the irregular satellites were captured as a result of gas drag within a circumjovian disk. Support for this research from NASA Outer Planets Research Program is gratefully acknowledged