Rossby wave instability does not require sharp resistivity gradients

Rossby wave instability (RWI) at dead zone boundaries may play an important role in planet formation. Viscous hydrodynamics results suggest RWI is excited only when the viscosity changes over a radial distance less than two density scale heights. However in the disks around Solar-mass T Tauri stars, it is not viscosity but magnetic forces that provide the accretion stress beyond about 10 AU, where surface densities are low enough so stellar X-rays and interstellar cosmic rays can penetrate. Here we aim to explore the conditions for RWI in the smooth transition with increasing distance, from resistive and magnetically-dead to conducting and magnetically-active. We perform 3D unstratified MHD simulations with the Pencil Code, using static resistivity profiles. As a result, we find that in MHD, contrary to viscous models, the RWI is triggered even with a gradual change in resistivity extending from 10 to 40 AU (i.e., spanning 15 scale heights for aspect ratio 0.1). This is because magneto-rotational turbulence sets in abruptly when the resistivity reaches a threshold level. At higher resistivities the longest unstable wavelength is quenched, resulting in a sharp decline of the Maxwell stress towards the star. The sharp gradient in the magnetic forces leads to a localized density bump, that is in turn Rossby wave unstable. We conclude that even weak gradients in the resistivity can lead to sharp transitions in the Maxwell stress. The upshot is that the RWI is more easily activated in the outer disk than previously thought. Rossby vortices at the outer dead zone boundary thus could underlie the dust asymmetries seen in the outer reaches of transition disks.

Formation of sharp eccentric rings in debris disks with gas but without planets

Debris disks around young stars (analogues of the Kuiper Belt in our Solar System) show a variety of non-trivial structures attributed to planetary perturbations and used to constrain the properties of the planets. However, these analyses have largely ignored the fact that some debris disks are found to contain small quantities of gas, a component that all such disks should contain at some level. Several debris disks have been measured with a dust-to-gas ratio around unity at which the effect of hydrodynamics on the structure of the disk cannot be ignored. Here we report linear and nonlinear modelling that shows that dust-gas interactions can produce some of the key patterns attributed to planets. We find a robust clumping instability that organizes the dust into narrow, eccentric rings, similar to the Fomalhaut debris disk. The conclusion that such disks might contain planets is not necessarily required to explain these systems.

Orbital migration of low-mass planets in evolutionary radiative models: Avoiding catastrophic infall

Outward migration of low-mass planets has recently been shown to be a possibility in non-barotropic disks. We examine the consequences of this result in evolutionary models of protoplanetary disks. Planet migration occurs towards equilibrium radii with zero torque. These radii themselves migrate inwards because of viscous accretion and photoevaporation. We show that as the surface density and temperature fall, the planet orbital migration and disk depletion timescales eventually become comparable, with the precise timing depending on the mass of the planet. When this occurs, the planet decouples from the equilibrium radius. At this time, however, the gas surface density is already too low to drive substantial further migration. A higher mass planet, of 10 Earth masses, can open a gap during the late evolution of the disk, and stops migrating. Low mass planets, with 1 or 0.1 Earth masses, released beyond 1 AU in our models, avoid migrating into the star. Our results provide support for the reduced migration rates adopted in recent planet population synthesis models.

Standing on the shoulders of giants: Trojan Earths and vortex trapping in low mass self-gravitating protoplanetary disks of gas and solids

Centimeter and meter sized solid particles in protoplanetary disks are trapped within long lived high pressure regions, creating opportunities for collapse into planetesimals and planetary embryos. We study the accumulations in the stable Lagrangian points of a giant planet, as well as in the Rossby vortices launched at the edges of the gap it carves. We employ the Pencil Code, tracing the solids with a large number of interacting Lagrangian particles, usually 100,000. For particles of 1 cm to 10 cm radii, gravitational collapse occurs in the Lagrangian points in less than 200 orbits. For 5 cm particles, a 2 Earth mass planet is formed. For 10 cm, the final maximum collapsed mass is around 3 Earth masses. The collapse of the 1 cm particles is indirect, following the timescale of depletion of gas from the tadpole orbits. In the edges of the gap vortices are excited, trapping preferentially particles of 30 cm radii. The rocky planet that is formed is as massive as 17 Earth masses, constituting a Super-Earth. By using multiple particle species, we find that gas drag modifies the streamlines in the tadpole region around the classical L4 and L5 points. As a result, particles of different radii have their stable points shifted to different locations. Collapse therefore takes longer and produces planets of lower mass. Three super-Earths are formed in the vortices, the most massive having 4.4 Earth masses. We conclude that a Jupiter mass planet can induce the formation of other planetary embryos in the outer edge of its gas gap. Trojan Earth mass planets are readily formed, and although not existing in the solar system, might be common in the exoplanetary zoo.

Embryos grown in the dead zone: Assembling the first protoplanetary cores in low mass self-gravitating circumstellar disks of gas and solids

In the borders of the dead zones of protoplanetary disks, the inflow of gas produces a local density maximum that triggers the Rossby wave instability. The vortices that form are efficient in trapping solids. We aim to assess the possibility of gravitational collapse of the solids within the Rossby vortices. We perform global simulations of the dynamics of gas and solids in a low mass non-magnetized self-gravitating thin protoplanetary disk with the Pencil code. We use multiple particle species of radius 1, 10, 30, and 100 cm. The dead zone is modeled as a region of low viscosity. The Rossby vortices excited in the edges of the dead zone are very efficient particle traps. Within 5 orbits after their appearance, the solids achieve critical density and undergo gravitational collapse into Mars sized objects. The velocity dispersions are of the order of 10 m/s for newly formed embryos, later lowering to less than 1 m/s by drag force cooling. After 200 orbits, 38 gravitationally bound embryos were formed inside the vortices, half of them being more massive than Mars. The embryos are composed primarily of same-sized particles. We conclude that the presence of a dead zone naturally gives rise to a population of protoplanetary cores in the mass range of 0.1-0.6 Earth masses, on very short timescales.

Global magnetohydrodynamical models of turbulence in protoplanetary disks I. A cylindrical potential on a Cartesian grid and transport of solids

We present global 3D MHD simulations of disks of gas and solids, aiming at developing models that can be used to study various scenarios of planet formation and planet-disk interaction in turbulent accretion disks. A second goal is to show that Cartesian codes are comparable to cylindrical and spherical ones in handling the magnetohydrodynamics of the disk simulations, as the disk-in-a-box models presented here develop and sustain MHD turbulence. We investigate the dependence of the magnetorotational instability on disk scale height, finding evidence that the turbulence generated by the magnetorotational instability grows with thermal pressure. The turbulent stresses depend on the thermal pressure obeying a power law of 0.24+/-0.03, compatible with the value of 0.25 found in shearing box calculations. The ratio of stresses decreased with increasing temperature. We also study the dynamics of boulders in the hydromagnetic turbulence. The vertical turbulent diffusion of the embedded boulders is comparable to the turbulent viscosity of the flow. Significant overdensities arise in the solid component as boulders concentrate in high pressure regions.