The asteroids (more precisely: objects of the main asteroid belt) and Kuiper Belt objects (more precisely: objects of the cold classical Kuiper Belt) are leftovers of the building material for our earth and all other planets in our solar system from more than 4.5 billion years ago. At the time of their formation those were typically 100 km large objects. They were called planetesimals, built up from icy and dusty grains. In our current paradigm of planet formation it was turbulent flows and metastable flow patterns, like zonal flows and vortices, that concentrated mm to cm sized icy dust grains in sufficient numbers that a streaming instability followed by a gravitational collapse of these particle clump was triggered. The entire picture is sometimes referred to as gravoturbulent formation of planetesimals. What was missing until recently, was a physically motivated prediction on the typical sizes at which planetesimals should form via this process. Our numerical simulations in the past had only shown a correlation between numerical resolution and planetesimal size and thus no answer was possible (Johansen et al.2011). But with the lastest series of simulations on JUQUEEN (Stephan & Doctor 2015), covering all the length scales down to the physical size of actual planetesimals, we were able to obtain values for the turbulent particle diffusion as a function of the particle load in the gas. Thus, we have all necessary data at hand to feed a 'back of the envelope' calculation that predicts the size of planetesimals as result of a competition between gravitational concentration and turbulent diffusion. Using the diffusion values obtained in the numerical simulations it predicts planetesimal sizes on the order of 100 km, which suprisingly coincides with the measured data from both asteroids (Bottke et al.2005) as well from Kuiper Belt objects (Nesvorny et al.2011).

We run global two dimensional hydrodynamical simulations, using the PLUTO code and the planet-disk model of Uribe et al. 2011, to investigate the effect of the convective overstability (CO) on planet-disk interactions. First, we study the long-term evolution of planet-induced vortices. We found that the CO leads to smoother planetary gap edges, thus weaker planet-induced vortices. The main result was the observation of two generation of vortices, which can pose an explanation for the location of the vortex in the Oph IRS48 system. The lifetime of the primary vortices, as well as the birth time of the secondary vortices are shown to be highly dependent on the thermal relaxation timescale. Second, we study the long-term evolution of the migration of low mass planets and assess whether the CO can prevent the saturation of the horseshoe drag. We found that the disk parameters that favour slow inward or outward migration oppose the amplification of vortices, meaning that the CO does not seem to be a good mechanism to prevent the saturation of the horseshoe drag. On the other hand, we observed a planetary trap, caused by vortices formed in the horseshoe region. This trap may be an alternative mechanism to prevent the fast type I migration rates.

The young (12+8−4 Myr) and nearby (19.44±0.05 pc) star β Pictoris is considered one of the best laboratories for the study of early phases of planetary systems formation since the identification of an extended debris disk surrounding the star in 1984. In 2009, we imaged at 3.8 μm with NaCo at VLT a gas giant planet around β Pictoris, roughly along the disk mid-plane, with a semi-major axis between 8 and 14 AU. We present here the first images of the planet in the J (1.265 μm), H (1.66 μm), and M' (4.78 μm) bands obtained between 2011 and 2012. We used these data to build the 1-5 μm spectral energy distribution (SED) of the companion, and to consolidate previous semi-major axis (8-10 AU) estimates. We compared the SED to seven atmospheric models to derive Teff = 1700 ± 100 K. We used the temperature and the luminosity of β Pictoris b to estimate new masses for the companion. We compared these masses to independent constraints set by the orbital parameters and the radial velocities and use them to discuss the formation history of the object.

Recent simulations show long -lived sub- and super-Keplerian flows in protoplanetary disks. These so-called zonal flows are found in local as well as global simulations of magneto-rotationally unstable disks. We investigate the strength and life-time of the resulting long-lived gas over- and under-densities as well as particle concentrations function of the azimuthal and radial size of the local shearing box. Changes in the azimuthal extent do not affect the zonal flow features. However, strength and life-time of zonal flows increase with increasing radial box sizes. Our simulations show indications, and support earlier results, that zonal flows have a natural length scale of approximately 5 pressure scale heights. For the first time, the reaction of dust particles in boxes with zonal flows are studied. We show that particles of some centimeters in size reach a hundred-fold higher density than initially, without any self-gravitating forces acting on the point masses. We further investigate collision velocities of dust grains in a turbulent medium.

Previous models of the combined growth and migration of protoplanets needed large ad hoc reduction factors for the type I migration rate as found in the isothermal approximation. In order to eliminate these factors, a simple semi-analytical model is presented that incorporates recent results on the migration of low mass planets in non-isothermal disks. It allows for outward migration. The model is used to conduct planetary populations synthesis calculations. Two points with zero torque are found in the disks. Planets migrate both in- and outward towards these convergence zones. They could be important for accelerating planetary growth by concentrating matter in one point. We also find that the updated type I migration models allow the formation of both close-in low mass planets, but also of giant planets at large semimajor axes. The problem of too rapid migration is significantly mitigated.

We performed 3D MHD numerical simulations of planet migration in stratified disks using the Godunov code PLUTO (Mignone et al. 2007). The disk is invaded by turbulence generated by the magnetorotational instability (MRI). We study the migration for planets with different mass to primary mass ratio. The migration of the low-mass planet (q=Mp/Ms=10−5) is dominated by random fluctuations in the torque and there is no defined direction of migration on timescales of 100 orbits. The intermediate-mass planet (q=Mp/Ms=10−4) can experience systematic outwards migration that was sustained for the times we were able to simulate.

3D global MHD simulations of magneto-driven turbulence are performed for the disk of 100 AU with reduced amount of 10μm fluffy dust grains. We use X-ray and cosmic ray ionization, as well as simplified treatment of recombination on dust grains. The ionization of gas and charging of dust grains are dynamically evolving during the simulation, making the zone of high magnetic dissipation (’dead’ zone) variable. In our simulations, the jump in MRI-driven turbulent viscosity inside and outside of dead zone is insignificant. We find no hard edge, but rather a smooth transition between active and dead zone. Subsequently, there is no visible pressure bump at outer edge of the dead zone.

We present high resolution computer simulations of dust dynamics and planetesimal formation in turbulence triggered by the magnetorotational instability. Particles representing approximately meter-sized boulders clump in large scale overpressure regions in the simulation box. These overdensities readily contract due to the combined gravity of the particles to form gravitationally bound clusters with masses ranging from a few to several ten times the mass of the dwarf planet Ceres. Gravitationally bound clumps are observed to collide and merge at both moderate and high resolution. The collisional products form the top end of a distribution of planetesimal masses ranging from less than one Ceres mass to 35 Ceres masses. It remains uncertain whether collisions are driven by dynamical friction or underresolution of clumps.

We present 3D global non-ideal MHD simulations with a self consistent dynamic evolution of ionization fraction of the gas as result of reduced chemical network. We include X-ray ionization from the star as well as cosmic ray ionization. Based on local gas density and temperature in our chemical network, we determine the magnetic resistivity, which is fed back in MHD simulations. Parameters for dust size and abundance are chosen to have accreting layers and a laminar “dead” mid-plane.

We present 3D global MHD simulations of proto-planetary disks calculated with the ZeusMP code. We focus on gas dynamics; the magnetic diffusivity and temperature are fixed during the simulation. A zone with low gas ionization at the midplane is included within ±2 scale heights. We mimic the ‘snow’-line radius with one order of magnitude jump in magnetic diffusivity η. Resulting turbulent Maxwell and Reynolds stresses are present at the midplane despite of low ionization. We find no radial inhomogeneities in turbulent α stress for a mild ionization contrast at the ‘snow’-line. A smooth azimuthal magnetic field is produced in the dead zone which may be a driving force for a weak accretion flow.

We briefly overview our newly developed radiation transport module for MHD simulations and two actual applications. The method combines the advantage of the speed of the Flux-Limited Diffusion approximation and the high accuracy obtained in ray-tracing methods.

For a long time it was believed that thermal convection could serve as the driving mechanism for turbulence and angular momentum transport in accretion disks. Even it is meanwhile accepted that convection had to leave that role to the magneto rotational instability, it is still an important effect arising in a realistic treatment of accretion disks, i.e. with proper thermodynamics and radiation transport. We review the history of thermal convection in astrophysical disks and show the relevant analytic and numerical work, including energy transport by convection and the effect of “negative” Reynolds stresses. We will also place the convective instability into the context of the magnetorotational instability and planet–disk interaction.