Recent observational findings have suggested a positive correlation between the occurrence rates of inner super-Earths and outer giant planets. These results raise the question if the trend can be reproduced and explained by planet formation theory. Here, we investigate the properties of inner super-Earths and outer giant planets that form according to the core accretion picture. The aim of this research was to study mutual relations between these planet species in synthetic planet populations and to compare them to the observed exoplanet population. We invoke the Generation 3 Bern model of planet formation and evolution to simulate 1000 multi-planet systems. We then confront these synthetic systems with the observed sample, taking into account the detection bias that distorts the observed demographics. The formation of warm super-Earths and cold Jupiters in the same system is enhanced compared to the individual appearances, albeit weaker than proposed by observations. We attribute the discrepancy to dynamically active giant planets, which frequently disrupt inner super-Earth systems in the synthetic population. This result could be interpreted as an overestimation of the migration rates predicted by giant planet migration theory. We further find differences in the volatile content of planets in different system architectures and predict that high-density super-Earths are more likely to host an outer giant. This correlation can be tested observationally. Finally, super-Earth occurrence shows a negative correlation with metallicity in giant-hosting systems. The reason for this is that giant planets around high-metallicity stars are frequently dynamically active, which poses a threat to inner planetary systems.

If we want to understand planetesimal formation, the only data set we have is our own Solar System. It is particularly interesting as it is so far the only planetary system we know of that developed life. Understanding the conditions under which the Solar Nebula evolved is crucial in order to understand the different processes in the disk and the subsequent dynamical interaction between (proto-)planets, once the gas disk is gone. Protoplanetary disks provide a plethora of different parameters to explore. The question is whether this parameter space can be constrained, allowing simulations to reproduce the Solar System. Models and observations of planet formation provide constraints on the initial planetesimal mass in certain regions of the Solar Nebula. By making use of pebble flux-regulated planetesimal formation, we perform a parameter study with nine different disk parameters like the initial disk mass, initial disk size, initial dust-to-gas ratio, turbulence level, and more. We find that the distribution of mass in planetesimals in the disk depends on the planetesimal formation timescale and the pebbles' drift timescale. Multiple disk parameters can influence pebble properties and thus planetesimal formation. However, it is still possible to draw some conclusions on potential parameter ranges. Pebble flux-regulated planetesimal formation seems to be very robust, allowing simulations with a wide range of parameters to meet the initial planetesimal constraints for the Solar Nebula. I.e., it does not require a lot of fine tuning.

Dusty debris disks around pre- and main-sequence stars are potential signposts for the existence of planetesimals and exoplanets. Giant planet formation is therefore expected to play a key role in the evolution of the disk. This is indirectly confirmed by extant sub-millimeter near-infrared images of young protoplanetary and cool dusty debris disks around main sequence stars usually showing substantial spatial structures. A majority of recent discoveries of imaged giant planets have been obtained around young, early-type stars hosting a circumstellar disk. In this context, we have carried out a direct imaging program designed to maximize our chances of giant planet discovery and targeting twenty-two young, early-type stars. About half of them show indication of multi-belt architectures. Using the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE to acquire high-constrast coronagraphic differential near-infrared images, we have conducted a systematic search in the close environment of these young, dusty and early-type stars. We confirmed that companions detected around HIP 34276, HIP 101800 and HIP 117452 are stationary background sources and binary companions. The companion candidates around HIP 8832, HIP 16095 and HIP 95619 are determined as background contamination. For stars for which we infer the presence of debris belts, a theoretical minimum mass for planets required to clear the debris gaps can be calculated . The dynamical mass limit is at least $0.1 M_J$ and can exceed $1 M_J$. Direct imaging data is typically sensitive to planets down to $\sim 3.6 M_J$ at 1 $''$, and $1.7 M_J$ in the best case. These two limits tightly constrain the possible planetary systems present around each target. These systems will be probably detectable with the next generation of planet imagers.

Theoretical models of the ionization state in protoplanetary disks suggest the existence of large areas with low ionization and weak coupling between the gas and magnetic fields. In this regime hydrodynamical instabilities may become important. In this work we investigate the gas and dust structure and dynamics for a typical T Tauri system under the influence of the vertical shear instability (VSI). We use global 3D radiation hydrodynamics simulations covering all $360^\circ$ of azimuth with embedded particles of 0.1 and 1mm size, evolved for 400 orbits. Stellar irradiation heating is included with opacities for 0.1- to 10-$\mu$m-sized dust. Saturated VSI turbulence produces a stress-to-pressure ratio of $\alpha \simeq 10^{-4}$. The value of $\alpha$ is lowest within 30~au of the star, where thermal relaxation is slower relative to the orbital period and approaches the rate below which VSI is cut off. The rise in $\alpha$ from 20 to 30~au causes a dip in the surface density near 35~au, leading to Rossby wave instability and the generation of a stationary, long-lived vortex spanning about 4~au in radius and 40~au in azimuth. Our results confirm previous findings that mm size grains are strongly vertically mixed by the VSI. The scale height aspect ratio for 1mm grains is determined to be 0.037, much higher than the value $H/r=0.007$ obtained from millimeter-wave observations of the HL~Tau system. The measured aspect ratio is better fit by non-ideal MHD models. In our VSI turbulence model, the mm grains drift radially inwards and many are trapped and concentrated inside the vortex. The turbulence induces a velocity dispersion of $\sim 12$~m/s for the mm grains, indicating that grain-grain collisions could lead to fragmentation.

We present a numerical code for radiation hydrodynamics designed as a module for the freely available PLUTO code. We adopt a gray approximation and include radiative transfer following a two-moment approach by imposing the M1 closure to the radiation fields. This closure allows for a description of radiative transport in both the diffusion and free-streaming limits, and is able to describe highly anisotropic radiation transport as can be expected in the vicinity of an accreting planet in a protoplanetary disk. To reduce the computational cost caused by the timescale disparity between radiation and matter fields, we integrate their evolution equations separately in an operator-split way, using substepping to evolve the radiation equations. We further increase the code's efficiency by adopting the reduced speed of light approximation (RSLA). Our integration scheme for the evolution equations of radiation fields relies on implicit-explicit schemes, in which radiation-matter interaction terms are integrated implicitly while fluxes are integrated via Godunov-type solvers. The module is suitable for general astrophysical computations in 1, 2, and 3 dimensions in Cartesian, spherical and cylindrical coordinates, and can be implemented on rotating frames. We demonstrate the algorithm performance on different numerical benchmarks, paying particular attention to the applicability of the RSLA for computations of physical processes in protoplanetary disks. We show 2D simulations of vertical convection in disks and 3D simulations of gas accretion by planetary cores, which are the first of their kind to be solved with a two-moment approach.

Forming gas giant planets by the accretion of 100 km diameter planetesimals, a typical size that results from self-gravity assisted planetesimal formation, is often thought to be inefficient. Many models therefore use small km-sized planetesimals, or invoke the accretion of pebbles. Furthermore, models based on planetesimal accretion often use the ad hoc assumption of planetesimals distributed radially in a minimum mass solar nebula fashion. We wish to investigate the impact of various initial radial density distributions in planetesimals with a dynamical model for the formation of planetesimals on the resulting population of planets. In doing so, we highlight the directive role of the early stages of dust evolution into pebbles and planetesimals in the circumstellar disk on the following planetary formation. We have implemented a two population model for solid evolution and a pebble flux regulated model for planetesimal formation into our global model for planet population synthesis. This framework is used to study the global effect of planetesimal formation on planet formation. As reference, we compare our dynamically formed planetesimal surface densities with ad-hoc set distributions of different radial density slopes of planetesimals. Even though required, it is not solely the total planetesimal disk mass, but the planetesimal surface density slope and subsequently the formation mechanism of planetesimals, that enables planetary growth via planetesimal accretion. Highly condensed regions of only 100 km sized planetesimals in the inner regions of circumstellar disks can lead to gas giant growth. Pebble flux regulated planetesimal formation strongly boosts planet formation, because it is a highly effective mechanism to create a steep planetesimal density profile. We find this to lead to the formation of giant planets inside 1 au by 100 km already by pure planetesimal accretion.

The formation of planetesimals is an unsolved problem in planet formation theory. A prominent scenario for overcoming dust growth barriers in dead zones is the gravitational collapse of locally over-dense regions, shown to robustly produce $\sim$100 km sized objects. Still, the conditions under which planetesimal formation occurs remain unclear. For collapse to proceed, the self-gravity of an overdensity must overcome stellar tidal disruption on large scales and turbulent diffusion on small scales. Here, we relate the scales of streaming and Kelvin-Helmholtz instability, which both regulate particle densities on the scales of gravitational collapse, directly to planetesimal formation. We support our analytic findings by performing 3D hydrodynamical simulations of streaming and Kelvin-Helmholtz instability and planetesimal formation. We find that the vertical extent of the particle mid-plane layer and the radial width of streaming instability filaments are set by the same characteristic length scale, thus governing the strength of turbulent diffusion on the scales of planetesimal formation. We present and successfully test a collapse criterion: $0.1 Q \beta \epsilon^{-1}Z^{-1} \lesssim 1$ and show that even for Solar metallicities, planetesimals can form in dead zones of sufficiently massive disks. For a given gas Toomre-parameter $Q$, pressure gradient $\beta$, metallicity $Z$ and local particle enhancement $\epsilon$, the collapse criterion also provides a range of unstable scales, instituting a promising path for studying initial planetesimal mass distributions. Streaming instability is not required for planetesimal collapse, but by increasing $\epsilon$, can evolve a system to instability.

Photoevaporation and magnetically driven winds are two independent mechanisms to remove mass from protoplanetary disks. In addition to accretion, the effect of these two principles acting concurrently could be significant and the transition between those two has not been extensively studied and quantified in the literature yet. In order to contribute to the understanding of disk winds, we present the phenomena emerging in the framework of two-dimensional axisymmetric, non-ideal magnetohydrodynamic simulations including EUV-/ X-ray driven photoevaporation. Of particular interest are the examination of the transition region between photoevaporation and magnetically driven wind, the possibility of emerging magneto-centrifugal wind effects, as well as the morphology of the wind itself depending on the strength of the magnetic field. We use the PLUTO code in a 2.5D axisymmetric configuration with additional treatment of EUV-/ X-ray heating and dynamic ohmic diffusion based on a semi-analytical chemical model. We identify the transition between both outflow types to occur for values of the initial plasma beta $\beta \geq 10^7$, while magnetically driven winds generally outperform photoevaporation for stronger fields. In our simulations we observe irregular and asymmetric outflows for stronger magnetic fields. In the weak field regime the photoevaporation rates are slightly lowered by perturbations of the gas density in the inner regions of the disk. Overall, our results predict a wind with a lever arm smaller than 1.5, consistent with a hot magneto-thermal wind. Stronger accretion flows are present for values of $\beta < 10^7$.

Statistical analyses from exoplanet surveys around low-mass stars indicate that super-Earth and Neptune-mass planets are more frequent than gas giants around such stars, in agreement with core accretion theory of planet formation. Using precise radial velocities derived from visual and near-infrared spectra, we report the discovery of a giant planet with a minimum mass of 0.46 Jupiter masses in an eccentric 204-day orbit around the very low-mass star GJ 3512. Dynamical models show that the high eccentricity of the orbit is most likely explained from planet-planet interactions. The reported planetary system challenges current formation theories and puts stringent constraints on the accretion and migration rates of planet formation and evolution models, indicating that disc instability may be more efficient in forming planets than previously thought.

Context. Teegarden's Star is the brightest and one of the nearest ultra-cool dwarfs in the solar neighbourhood. For its late spectral type (M7.0V), the star shows relatively little activity and is a prime target for near-infrared radial velocity surveys such as CARMENES. Aims. As part of the CARMENES search for exoplanets around M dwarfs, we obtained more than 200 radial-velocity measurements of Teegarden's Star and analysed them for planetary signals. Methods. We find periodic variability in the radial velocities of Teegarden's Star. We also studied photometric measurements to rule out stellar brightness variations mimicking planetary signals. Results. We find evidence for two planet candidates, each with $1.1M_\oplus$ minimum mass, orbiting at periods of 4.91 and 11.4 d, respectively. No evidence for planetary transits could be found in archival and follow-up photometry. Small photometric variability is suggestive of slow rotation and old age. Conclusions. The two planets are among the lowest-mass planets discovered so far, and they are the first Earth-mass planets around an ultra-cool dwarf for which the masses have been determined using radial velocities.

Recent high-resolution interferometric observations of protoplanetary disks at (sub-)millimeter wavelengths reveal omnipresent substructures, such as rings, spirals, and asymmetries. A detailed investigation of eight rings detected in five disks by the DSHARP survey came to the conclusion that all rings are just marginally optically thick with optical depths between 0.2 and 0.5 at a wavelength of 1.25 mm. This surprising result could either be coincidental or indicate that the optical depth in all of the rings is regulated by the same process. We investigated if ongoing planetesimal formation could explain the "fine-tuned" optical depths in the DSHARP rings by removing dust and transforming it into "invisible" planetesimals. We performed a one-dimensional simulation of dust evolution in the second dust ring of the protoplanetary disk around HD 163296, including radial transport of gas and dust, dust growth and fragmentation, and planetesimal formation via gravitational collapse of sufficiently dense pebble concentrations. We show that planetesimal formation can naturally explain the observed optical depths if streaming instability regulates the midplane dust-to-gas ratio to unity. Furthermore, our simple monodisperse analytical model supports the hypothesis that planetesimal formation in dust rings should universally limit their optical depth to the observed range.

If planetesimal formation is an efficient process, as suggested by several models involving gravitational collapse of pebble clouds, then, before long, a significant part of the primordial dust mass should be absorbed in many km sized objects. A good understanding of the total amount of solids in the disk around a young star is crucial for planet formation theory. But as the mass of particles above the mm size cannot be assessed observationally, one must ask how much mass is hidden in bigger objects. We perform 0-d local simulations to study how the planetesimal to dust and pebble ratio is evolving in time and to develop an understanding of the potentially existing mass in planetesimals for a certain amount of dust and pebbles at a given disk age. We perform a parameter study based on a model considering dust growth, planetesimal formation and collisional fragmentation of planetesimals, while neglecting radial transport processes. While at early times, dust is the dominant solid particle species, there is a phase during which planetesimals make up a significant portion of the total mass starting at approximately $10^4 - 10^6$ yr. The time of this phase and the maximal total planetesimal mass strongly depend on the distance to the star $R$, the initial disk mass, and the efficiency of planetesimal formation $\epsilon$. After approximately $10^6$ yr, our model predicts planetesimal collisions to dominate, which resupplies small particles. In our model, planetesimals form fast and everywhere in the disk. For a given $\epsilon$, we were able to relate the dust content and mass of a given disk to its planetesimal content, providing us with some helpful basic intuition about mass distribution of solids and its dependence on underlying physical processes.

We report the detection of a transiting Earth-size planet around GJ 357, a nearby M2.5V star, using data from the Transiting Exoplanet Survey Satellite (TESS). GJ 357 b (TOI-562.01) is a transiting, hot, Earth-sized planet (Teq=525+-11 K) with a radius of Rb=1.217+-0.084 Re and an orbital period of Pb=3.93 d. Precise stellar radial velocities from CARMENES and PFS, as well as archival data from HIRES, UVES, and HARPS also display a 3.93-day periodicity, confirming the planetary nature and leading to a planetary mass of Mb=1.84+-0.31 Me. In addition to the radial velocity signal for GJ 357 b, more periodicities are present in the data indicating the presence of two further planets in the system: GJ 357 c, with a minimum mass of Mc=3.40+-0.46 Me in a 9.12 d orbit, and GJ 357 d, with a minimum mass of Md=6.1+-1.0 Me in a 55.7 d orbit inside the habitable zone. The host is relatively inactive and exhibits a photometric rotation period of Prot=78+-2 d. GJ 357 b is to date the second closest transiting planet to the Sun, making it a prime target for further investigations such as transmission spectroscopy. Therefore, GJ 357 b represents one of the best terrestrial planets suitable for atmospheric characterization with the upcoming JWST and ground-based ELTs.

Due to the gas rich environments of early circumstellar disks, the gravitational collapse of cool, dense regions of the disk form fragments largely composed of gas. During formation, disk fragments may attain increased metallicities as they interact with the surrounding disk material, whether through particle migration to pressure maxima or through mutual gravitational interaction. In this paper, we investigate the ability of fragments to collect and retain a significant solid component through gas-particle interactions in high-resolution 3D self-gravitating shearing box simulations. The formation of axissymmetric perturbations associated with gravitational instabilities allows particles of intermediate sizes to concentrate through aerodynamic drag forces. By the onset of fragmentation, the mass of local particle concentrations within the fragment are comparable to that of the gas component and the sebsequent gravitational collapse results in the formation of a solid core. We find that these cores can be up to several tens of Earth masses, depending on grain size, before the fragment center reaches temperatures which would sublimate solids. The solid fraction and total mass of the fragment also depend on the metallicity of the young parent protoplanetary disk, with higher initial metallicities resulting in larger fragments and larger solid cores. Additionally, the extended atmospheres of these soon-to-be gas giants or brown dwarfs are occasionally enriched above the initial metallicity, provided no solid core forms in the center and are otherwise lacking in heavier elements when a core does form.

PDS 70b is the most robustly detected young planet imaged in the gap of a transition disk so far, found at a separation of ~195 mas (~22 au) from its host star and at a position angle of about 155 deg. This system is therefore a unique laboratory to characterize the properties of young planetary systems at the stage of their formation. We aim to trace direct and indirect imprints of PDS 70b on the gas and dust emission of the circumstellar disk in order to study the properties of this ~5 Myr young planetary system. We obtained ALMA band 7 observations of PDS 70 in dust continuum and $^{12}CO$ (3-2) and combined them with archival data resulting in an unprecedented angular resolution of about 70 mas (~8 au). We derive an upper limit on circumplanetary material at the location of PDS 70b of ~0.01 $M_{\oplus}$ and find a highly structured circumstellar disk in both dust and gas. The outer dust ring peaks at ~0.65" (74 au) and reveals a possible second unresolved peak at ~0.53" (60 au). The integrated intensity of CO also shows evidence of a depletion of emission at ~0.2" (23 au) with a width of ~0.1" (11 au). The gas kinematics show evidence of a deviation from Keplerian rotation inside $\lesssim$0.8" (91 au). This implies a pressure gradient that can account for the location of the dust ring well beyond the location of PDS 70b. Further in, we detect an inner disk which appears to be connected to the outer disk by a possible bridge feature in the North-West region in both gas and dust. We compare the observations to hydrodynamical simulations including a planet with different masses covering the estimated mass range previously derived from NIR photometry (~5-9 M$_{Jup}$). We find that even a planet with a mass of 10 M$_{Jup}$ may not be sufficient to explain the extent of the wide gap and an additional low-mass companion may be needed to account for the observed disk morphology.

We propose an expression for a local planetesimal formation rate proportional to the instantaneous radial pebble flux. The result --- a radial planetesimal distribution --- can be used as initial condition to study the formation of planetary embryos. We follow the idea that one needs particle traps to locally enhance the dust-to-gas ratio sufficiently such that particle gas interactions can no longer prevent planetesimal formation on small scales. The location of these traps can emerge everywhere in the disk. Their occurrence and lifetime is subject of ongoing research, thus they are implemented via free parameters. This enables us to study the influence of the disk properties on the formation of planetesimals, predicting their time dependent formation rates and location of primary pebble accretion. We show that large $\alpha$-values of $0.01$ (strong turbulence) prevent the formation of planetesimals in the inner part of the disk, arguing for lower values of around $0.001$ (moderate turbulence), at which planetesimals form quickly at all places where they are needed for proto-planets. Planetesimals form as soon as dust has grown to pebbles ($\sim\mathrm{mm}$ to $\mathrm{dm}$) and the pebble flux reaches a critical value, which is after a few thousand years at $2-3\,$AU and after a few hundred thousand years at $20-30\,$AU. Planetesimal formation lasts until the pebble supply has decreased below a critical value. The final spatial planetesimal distribution is steeper compared to the initial dust and gas distribution which helps to explain the discrepancy between the minimum mass solar nebula and viscous accretion disks.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

The streaming instability, as an example of instabilities driven by particle feedback on a gas flow, has been proven to have a major role in controlling the formation of planetesimals. Here, we present experiments to approach this situation in the laboratory for particles in the Knudsen flow regime. In these experiments, we observe a particle cloud trapped for about 30 s in a rotating system under Earth's gravity. For average dust-to-gas ratios up to 0.08, particles behave like individual test particles. Their sedimentation speed is identical to that of a single free-falling particle, even in locally denser regions. However, for higher dust-to-gas ratios, the motion of particles becomes sensitive to clumping. Particles in locally denser regions now sediment faster. Their sedimentation speed then depends linearly on the overall dust-to-gas ratio. This clearly shows a transition from tracerlike behavior to collective behavior. Beyond these findings, these types of experiments can now be used as a gauge to test the particle feedback models in astrophysical hydrocodes which are currently used for numerical simulations of streaming instabilities.

Hydrodynamical instabilities in disks around young stars depend on the thermodynamic stratification of the disk and on the local rate of thermal relaxation. Here, we map the spatial extent of unstable regions for the Vertical Shear Instability (VSI), the Convective OverStability (COS), and the amplification of vortices via the Subcritical Baroclinic Instability (SBI). We use steady state accretion disk models, including stellar irradiation, accretion heating and radiative transfer. We determine the local radial and vertical stratification and thermal relaxation rate in the disk, in dependence of the stellar mass, disk mass and mass accretion rate. We find that passive regions of disks - i.e. the midplane temperature dominated by irradiation - are COS unstable about one pressure scale height above the midplane and VSI unstable at radii $> 10 \, \text{au}$. Vortex amplification via SBI should operate in most parts of active and passive disks. For active parts of disks (midplane temperature determined by accretion power) COS can become active down to the midplane. Same is true for the VSI because of the vertically adiabatic stratification of an internally heated disk. If hydro instabilities or other non-ideal MHD processes are able to create $\alpha$-stresses ($> 10^{-5}$) and released accretion energy leads to internal heating of the disk, hydrodynamical instabilities are likely to operate in significant parts of the planet forming zones in disks around young stars, driving gas accretion and flow structure formation. Thus hydro-instabilities are viable candidates to explain the rings and vortices observed with ALMA and VLT.

This chapter highlights the properties of turbulence and meso-scale flow structures in protoplanetary disks and their role in the planet formation process. Here we focus on the formation of planetesimals from a gravitational collapse of a pebble cloud. Large scale and long lived flow structures - vortices and zonal flows - are a consequence of weak magneto and hydrodynamic instabilities in the pressure and entropy stratified quasi-Keplerian shear flow interacting with the fast rotation of the disk. The vortices and zonal flows on the other hand are particle traps tapping into the radial pebble flux of the disk, leading to locally sufficient accumulations to trigger gravitational collapse, directly converting pebbles to many kilometer sized planetesimals. This collapse is moderated by the streaming instability, which is a back-reaction from the particle accumulations onto the gas flow. Without trapping pebbles and increasing thus the local solid to gas ratio, this back-reaction would ultimately prevent the formation of planetsimals via turbulent diffusion. The formation of long lived flow structures is therefore a necessary condition for an efficient and fast formation of planetesimals.

The formation of planets remains one of the most challenging problems of contemporary astrophysics. Starting with micron-sized dust grains, coagulation models predict growth up to centimeter (pebbles), but growth beyond this size is difficult because of fragmentation and drift. Ways to bypass this problem have focused on inhomogeneities in the flow, be that zonal flows, streaming instability, or vortices. Because vortices are in equilibrium between the Coriolis and the pressure force, the pressureless grains will orbit along a vortex streamline experiencing a drag force. This is a very effective mechanism to concentrate pebbles as also seen in numerical simulations and possibly in ALMA observations. Yet, a high pebble load is dangerous for the vortex, and we showed that in two-dimensional simulations the backreaction eventually leads to vortex disruption. We investigate whether the same happens in three dimensions. We perform 3D simulations with pebbles in a local box finding that, although the pebbles disturb the vortex around the midplane, the column does not get destroyed. This result is important because, based on the previous 2D result suggesting complete disruption, the vortex interpretation of ALMA observations has been called into question. We show instead that the vortex behaves like a Taylor column, and the pebbles as obstacles to the flow. Pebble accumulation in the center of the vortices proceeds to roughly the same concentration as in the control run without backreaction.

Context. A low-mass brown dwarf has been recently imaged around HR 2562 (HD 50571), a star hosting a debris disk resolved in the far infrared. Interestingly, the companion location is compatible with an orbit coplanar with the disk and interior to the debris belt. This feature makes the system a valuable laboratory to analyze the formation of substellar companions in a circumstellar disk and potential disk-companion dynamical interactions. Aims. We aim to further characterize the orbital motion of HR 2562 B and its interactions with the host star debris disk. Methods. We performed a monitoring of the system over ~10 months in 2016 and 2017 with the VLT/SPHERE exoplanet imager. Results. We confirm that the companion is comoving with the star and detect for the first time an orbital motion at high significance, with a current orbital motion projected in the plane of the sky of 25 mas (~0.85 au) per year. No orbital curvature is seen in the measurements. An orbital fit of the SPHERE and literature astrometry of the companion without priors on the orbital plane clearly indicates that its orbit is (quasi-)coplanar with the disk. To further constrain the other orbital parameters, we used empirical laws for a companion chaotic zone validated by N-body simulations to test the orbital solutions that are compatible with the estimated disk cavity size. Non-zero eccentricities (>0.15) are allowed for orbital periods shorter than 100 yr, while only moderate eccentricities up to ~0.3 for orbital periods longer than 200 yr are compatible with the disk observations. A comparison of synthetic Herschel images to the real data does not allow us to constrain the upper eccentricity of the companion.

Several protoplanetary disks observed by ALMA show dust concentrations consistent with particle trapping in giant vortices. The formation and survival of vortices is of major importance for planet formation, because vortices act as particle traps and are therefore preferred locations of planetesimal formation. Recent studies showed that the vertical shear instability (VSI) is capable of generating turbulence and small vortices in protoplanetary disks that have the proper radial and vertical stratification and thermally relax on sufficiently short time scales. But the effect of the azimuthal extend of the disk is often neglected as the disks azimuth is limited to $\Delta \phi \leq \pi/2$. We aim to investigate the influence of the azimuthal extent of the disk on the long-term evolution of a protoplanetary disk and the possibility of large vortices forming. To this end, we perform 3-dimensional simulations for up to 1000 local orbits using different values of $\Delta \phi = \pi/2 $ to $2\pi$ for VSI in disks with a prescribed radial density and temperature gradient cooling on short timescales. We find the VSI capable of forming large vortices which can exist at least several hundred orbits in simulations covering a disk with $\Delta \phi \geq \pi$. This suggests the VSI to be capable to form vortices or at least to trigger vortex formation via a secondary instability, e.g. Rossby Wave Instability or Kelvin Helmholtz Instability.

The G-type star GJ504A is known to host a 3 to 35 MJup companion whose temperature, mass, and projected separation all contribute to make it a test case for the planet formation theories and for atmospheric models of giant planets and light brown dwarfs. We collected data from the CHARA interferometer, SOPHIE spectrograph, and VLT/SPHERE high contrast imager to revisit the properties of the system. We measure a radius of 1.35+/- 0.04Rsun for GJ504A which yields isochronal ages of 21+/-2Myr or 4.0+/-1.8Gyr for the system and line-of-sight stellar rotation axis inclination of $162.4_{-4.3}^{+3.8}$ degrees or $18.6_{-3.8}^{+4.3}$ degrees. We re-detect the companion in the Y2, Y3, J3, H2, and K1 dual band SPHERE images. The complete 1-4 $\mu$m SED shape of GJ504b is best reproduced by T8-T9.5 objects with intermediate ages ($\leq1.5$Gyr), and/or unusual dusty atmospheres and/or super-solar metallicities. All six atmospheric models used yield $\mathrm{T_{eff}=550 \pm 50}$K for GJ504b and point toward a low surface gravity (3.5-4.0 dex). The accuracy on the metallicity value is limited by model-to-model systematics. It is not degenerate with the C/O ratio. We derive $\mathrm{log\:L/L_{\odot}=-6.15\pm0.15}$ dex for the companion compatible with masses of $\mathrm{M=1.3^{+0.6}_{-0.3}M_{Jup}}$ and $\mathrm{M=23^{+10}_{-9} M_{Jup}}$ for the young and old age ranges, respectively. The semi-major axis (sma) is above 27.8 au and the eccentricity lower than 0.55. The posterior on GJ~504b's orbital inclination suggests a misalignment with GJ~504A rotation axis. We combine the radial velocity and multi-epoch imaging data to exclude additional objects (90\% prob.) more massive than 2.5 and 30 $\mathrm{M_{Jup}}$ with sma in the range 0.01-80 au for the young and old system ages, respectively. The companion is in the envelope of the population of planets synthetized with our core-accretion model.

The collapse of dust particle clouds directly to km-sized planetesimals is a promising way to explain the formation of planetesimals, asteroids and comets. In the past, this collapse has been studied in stratified shearing box simulations with super-solar dust-to-gas ratio \epsilon, allowing for streaming instability (SI) and gravitational collapse. This paper studies the non-stratified SI under dust-to-gas ratios from \epsilon=0.1 up to \epsilon=1000 without self-gravity. The study covers domain sizes of L=0.1 H, 0.01 H and 0.001 H, in terms of gas disk scale height H, using the PencilCode. They are performed in radial-azimuthal (2-d) and radial-vertical (2.5-d) extent. The used particles of St=0.01 and 0.1 mark the upper end of the expected dust growth. SI-activity is found up to very high dust-to-gas ratios, providing fluctuations in the local dust-to-gas ratios and turbulent particle diffusion \delta. We find an SI-like instability that operates in r-\varphi even when vertical modes are suppressed. This new azimuthal streaming instability (aSI) shows similar properties and appearance as the SI. Both, SI and aSI, show diffusivity at \epsilon=100 only to be two orders of magnitude lower than at \epsilon=1, suggesting a \delta ~ \epsilon^{-1} relation that is shallow around \epsilon = 1. The (a)SI ability to concentrate particles is found to be uncorrelated with its strength in particle turbulence. Finally, we performed a resolution study to test our findings of the aSI. This paper stresses out the importance of properly resolving the (a)SI at high dust-to-gas ratios and planetesimal collapse simulations, leading else wise to potentially incomplete results.

We report the discovery of a bright, brown dwarf companion to the star HIP 64892, imaged with VLT/SPHERE during the SHINE exoplanet survey. The host is a B9.5V member of the Lower-Centaurus-Crux subgroup of the Scorpius Centaurus OB association. The measured angular separation of the companion ($1.2705\pm0.0023$") corresponds to a projected distance of $159\pm12$ AU. We observed the target with the dual-band imaging and long-slit spectroscopy modes of the IRDIS imager to obtain its SED and astrometry. In addition, we reprocessed archival NACO L-band data, from which we also recover the companion. Its SED is consistent with a young (<30 Myr), low surface gravity object with a spectral type of M9$_{\gamma}\pm1$. From comparison with the BT-Settl atmospheric models we estimate an effective temperature of $T_{\textrm{eff}}=2600 \pm 100$ K, and comparison of the companion photometry to the COND evolutionary models yields a mass of $\sim29-37$ M$_{\text{J}}$ at the estimated age of $16^{+15}_{-7}$ Myr for the system. HIP 64892 is a rare example of an extreme-mass ratio system ($q\sim0.01$) and will be useful for testing models relating to the formation and evolution of such low-mass objects.

Transition disks are considered sites of ongoing planet formation, and their dust and gas distributions could be signposts of embedded planets. The transition disk around the T Tauri star RY Lup has an inner dust cavity and displays a strong silicate emission feature. Using high-resolution imaging we study the disk geometry, including non-axisymmetric features, and its surface dust grain, to gain a better understanding of the disk evolutionary process. Moreover, we search for companion candidates, possibly connected to the disk. We obtained high-contrast and high angular resolution data in the near-infrared with the VLT/SPHERE extreme adaptive optics instrument whose goal is to study the planet formation by detecting and characterizing these planets and their formation environments through direct imaging. We performed polarimetric imaging of the RY~Lup disk with IRDIS (at 1.6 microns), and obtained intensity images with the IRDIS dual-band imaging camera simultaneously with the IFS spectro-imager (0.9-1.3 microns). We resolved for the first time the scattered light from the nearly edge-on circumstellar disk around RY~Lup, at projected separations in the 100 \,au range. The shape of the disk and its sharp features are clearly detectable at wavelengths ranging from 0.9 to 1.6 microns. We show that the observed morphology can be interpreted as spiral arms in the disk. This interpretation is supported by in-depth numerical simulations. We also demonstrate that these features can be produced by one planet interacting with the disk. We also detect several point sources which are classified as probable background objects.

Massive circumstellar disks are prone to gravitational instabilities, which trigger the formation of spiral arms that can fragment into bound clumps under the right conditions. Two dimensional simulations of self-gravitating disks are useful starting points for studying fragmentation, allowing for high-resolution simulations of thin disks. However, convergence issues can arise in 2D from various sources. One of these sources is the 2D approximation of self-gravity, which exaggerates the effect of self-gravity on small scales when the potential is not smoothed to account for the assumed vertical extent of the disk. This effect is enhanced by increased resolution, resulting in fragmentation at longer cooling timescales $\beta$. If true, it suggests that the 3D simulations of disk fragmentation may not have the same convergence problem and could be used to examine the nature of fragmentation without smoothing self-gravity on scales similar to the disk scale height. To that end, we have carried out local 3D self-gravitating disk simulations with simple $\beta$ cooling with fixed background irradiation to determine if 3D is necessary to properly describe disk fragmentation. Above a resolution of $\sim 40$ grid cells per scale height, we find that our simulations converge with respect to the cooling timescale. This result converges in agreement with analytic expectations which place a fragmentation boundary at $\beta_\mathrm{crit} = 3$.

Hydrodynamic, non-magnetic instabilities can provide turbulent stress in the regions of protoplanetary discs, where the MRI can not develop. The induced motions influence the grain growth, from which formation of planetesimals begins. Thermal relaxation of the gas constrains origins of the identified hydrodynamic sources of turbulence in discs. We estimate the radiative relaxation timescale of temperature perturbations and study the dependence of this timescale on the perturbation wavelength, the location within the disc, the disc mass, and the dust-to-gas mass ratio. We then apply thermal relaxation criteria to localise modes of the convective overstability, the vertical shear instability, and the zombie vortex instability. Our calculations employed the latest tabulated dust and gas mean opacities and we account for the collisional coupling to the emitting species. The relaxation criterion defines the bulk of a typical T Tauri disc as unstable to the development of linear hydrodynamic instabilities. The midplane is unstable to the convective overstability from at most $2\mbox{ au}$ and up to $40\mbox{ au}$, as well as beyond $140\mbox{ au}$. The vertical shear instability can develop between $15\mbox{ au}$ and $180\mbox{ au}$. The successive generation of (zombie) vortices from a seeded noise can work within the inner $0{.}8\mbox{ au}$. Dynamic disc modelling with the evolution of dust and gas opacities is required to clearly localise the hydrodynamic turbulence, and especially its non-linear phase.

The key aspect determining the post-formation luminosity of gas giants has long been considered to be the energetics of the accretion shock at the planetary surface. We use 1D radiation-hydrodynamical simulations to study the radiative loss efficiency and to obtain post-shock temperatures and pressures and thus entropies. The efficiency is defined as the fraction of the total incoming energy flux which escapes the system (roughly the Hill sphere), taking into account the energy recycling which occurs ahead of the shock in a radiative precursor. We focus here on a constant equation of state to isolate the shock physics but use constant and tabulated opacities. While robust quantitative results will require a self-consistent treatment including hydrogen dissocation and ionization, the results show the correct qualitative behavior and can be understood semi-analytically. The shock is found to be isothermal and supercritical for a range of conditions relevant to core accretion (CA), with Mach numbers greater than ca. 3. Across the shock, the entropy decreases significantly, by a few entropy units (k_B/baryon). While nearly 100 percent of the incoming kinetic energy is converted to radiation locally, the efficiencies are found to be as low as roughly 40 percent, implying that a meaningful fraction of the total accretion energy is brought into the planet. For realistic parameter combinations in the CA scenario, a non-zero fraction of the luminosity always escapes the system. This luminosity could explain, at least in part, recent observations in the LkCa 15 and HD 100546 systems.

The first long-baseline ALMA campaign resolved the disk around the young star HL Tau into a number of axisymmetric bright and dark rings. Despite the very young age of HL Tau these structures have been interpreted as signatures for the presence of (proto)planets. The ALMA images triggered numerous theoretical studies based on disk-planet interactions, magnetically driven disk structures, and grain evolution. Of special interest are the inner parts of disks, where terrestrial planets are expected to form. However, the emission from these regions in HL Tau turned out to be optically thick at all ALMA wavelengths, preventing the derivation of surface density profiles and grain size distributions. Here, we present the most sensitive images of HL Tau obtained to date with the Karl G. Jansky Very Large Array at 7.0 mm wavelength with a spatial resolution comparable to the ALMA images. At this long wavelength the dust emission from HL Tau is optically thin, allowing a comprehensive study of the inner disk. We obtain a total disk dust mass of 0.001 - 0.003 Msun, depending on the assumed opacity and disk temperature. Our optically thin data also indicate fast grain growth, fragmentation, and formation of dense clumps in the inner densest parts of the disk. Our results suggest that the HL Tau disk may be actually in a very early stage of planetary formation, with planets not already formed in the gaps but in the process of future formation in the bright rings.

Dust grains with sizes around (sub)mm are expected to couple only weakly to the gas motion in regions beyond 10 au of circumstellar disks. In this work, we investigate the influence of the spatial distribution of such grains on the (sub)mm appearance of magnetized protoplanetary disks. We perform non-ideal global 3D magneto-hydrodynamic stratified disk simulations including particles of different sizes (50 micron to 1 cm), using a Lagrangian particle solver. We calculate the spatial dust temperature distribution, including the dynamically coupled submicron-sized dust grains, and derive ideal continuum re-emission maps of the disk through radiative transfer simulations. Finally, we investigate the feasibility to observe specific structures in the thermal re-emission maps with the Atacama Large Millimeter/submillimeter Array (ALMA). The pressure bump close to the outer edge of the dead-zone leads to particle trapping in ring structures. More specifically, vortices in the disk concentrate the dust and create an inhomogeneous distribution of the solid material in the azimuthal direction. The large-scale disk perturbations are preserved in the (sub)mm re-emission maps. The observable structures are very similar to those expected to result from planet-disk interaction. The larger dust particles increase the brightness contrast between gap and ring structures. We find that rings, gaps and the dust accumulation in the vortex could be traced with ALMA down to a scale of a few astronomical units in circumstellar disks located in nearby star-forming regions. Finally, we present a brief comparison of these structures with those recently found with ALMA in the young circumstellar disks of HL Tau and Oph IRS 48.

Protoplanetary disks fragment due to gravitational instability when there is enough mass for self-gravitation, described by the Toomre parameter, and when heat can be lost at a rate comparable to the local dynamical timescale, described by t_c=beta Omega^-1. Simulations of self-gravitating disks show that the cooling parameter has a rough critical value at beta_crit=3. When below beta_crit, gas overdensities will contract under their own gravity and fragment into bound objects while otherwise maintaining a steady state of gravitoturbulence. However, previous studies of the critical cooling parameter have found dependence on simulation resolution, indicating that the simulation of self-gravitating protoplanetary disks is not so straightforward. In particular, the simplicity of the cooling timescale t_c prevents fragments from being disrupted by pressure support as temperatures rise. We alter the cooling law so that the cooling timescale is dependent on local surface density fluctuations, a means of incorporating optical depth effects into the local cooling of an object. For lower resolution simulations, this results in a lower critical cooling parameter and a disk more stable to gravitational stresses suggesting the formation of large gas giants planets in large, cool disks is generally suppressed by more realistic cooling. At our highest resolution however, the model becomes unstable to fragmentation for cooling timescales up to beta = 10.

We study the evolution of planet-induced vortices in radially stratified disks, with initial conditions allowing for radial buoyancy. For this purpose we run global two dimensional hydrodynamical simulations, using the PLUTO code. Planet-induced vortices are a product of the Rossby wave instability (RWI) triggered in the edges of a planetary gap. In this work we assess the influence of radial buoyancy for the development of the vortices. We found that radial buoyancy leads to smoother planetary gaps, which generates weaker vortices. This effect is less pronounced for locally isothermal and quasi-isothermal (very small cooling rate) disks. We observed the formation of two generations of vortices. The first generation of vortices is formed in the outer wall of the planetary gap. The merged primary vortex induces accretion, depleting the mass on its orbit. This process creates a surface density enhancement beyond the primary vortex position. The second generation of vortices arise in this surface density enhancement, indicating that the bump in this region is sufficient to trigger the RWI. The merged secondary vortex is a promising explanation for the location of the vortex in the Oph IRS48 system. Finally, we observed a nonmonotonic behavior for the vortex lifetimes as a function of the thermal relaxation timescale, agreeing with previous studies. The birth times of the secondary vortices also display a nonmonotonic behavior, which is correlated with the growth time of the primary vortex and its induced accretion.

Circumstellar disks are considered to be the birthplace of planets. Specific structures like spiral arms, gaps, and cavities are characteristic indicators of planet-disk interaction. Investigating these structures can provide insights into the growth of protoplanets and the physical properties of the disk. We investigate the feasibility of using molecular lines to trace planet-induced structures in circumstellar disks. Based on 3D hydrodynamic simulations of planet-disk interactions, we perform self-consistent temperature calculations and produce N-LTE molecular line velocity-channel maps and spectra of these disks using our new N-LTE line radiative transfer code Mol3D. Subsequently, we simulate ALMA observations using the CASA simulator. We consider two nearly face-on inclinations, 5 disk masses, 7 disk radii, and 2 different typical pre-main-sequence host stars (T Tauri, Herbig Ae). We calculate up to 141 individual velocity-channel maps for five molecules/isotopoloques in a total of 32 rotational transitions to investigate the frequency dependence of the structures indicated above. We find that the majority of protoplanetary disks in our parameter space could be detected in the molecular lines considered. However, unlike the continuum case, gap detection is not straightforward in lines. For example, gaps are not seen in symmetric rings but are masked by the pattern caused by the global (Keplerian) velocity field. We identify specific regions in the velocity-channel maps that are characteristic of planet-induced structures. Simulations of high angular resolution molecular line observations demonstrate the potential of ALMA to provide complementary information about the planet-disk interaction as compared to continuum observations. In particular, the detection of planet-induced gaps is possible under certain conditions.(abridged)

We analyse the concentration of solid particles in vortices created and sustained by radial buoyancy in protoplanetary disks, i.e. baroclinic vortex growth. Besides the gas drag acting on particles we also allow for back-reaction from dust onto the gas. This becomes important when the local dust-to-gas ratio approaches unity. In our 2D, local, shearing sheet simulations we see high concentrations of grains inside the vortices for a broad range of Stokes numbers, ${\rm St}$. An initial dust-to-gas ratio of 1:100 can easily be reversed to 100:1 for ${\rm St}=1$. The increased dust-to-gas ratio triggers the streaming instability, thus counter-intuitively limiting the maximal achievable overdensities. We find that particle trapping inside vortices opens the possibility for gravity-assisted planetesimal formation even for small particles ($\rm{St}=0.01$) and low initial dust-to-gas ratios (1:$10^4$)

Recent observations by the Atacama Large Millimeter/submillimeter Array (ALMA) of disks around young stars revealed distinct asymmetries in the dust continuum emission. In this work we want to study axisymmetric and non-axisymmetric structures, evocated by the magneto-rotational instability in the outer regions of protoplanetary disks. We combine the results of state-of-the-art numerical simulations with post-processing radiative transfer (RT) to generate synthetic maps and predictions for ALMA. We performed non-ideal global 3D MHD stratified simulations of the dead-zone outer edge using the FARGO MHD code PLUTO. The stellar and disk parameters are taken from a parameterized disk model applied for fitting high-angular resolution multi-wavelength observations of circumstellar disks. The 2D temperature and density profiles are calculated consistently from a given surface density profile and Monte-Carlo radiative transfer. The 2D Ohmic resistivity profile is calculated using a dust chemistry model. The magnetic field is a vertical net flux field. The resulting dust reemission provides the basis for the simulation of observations with ALMA. The fiducial model develops a large gap followed by a jump in surface density located at the dead-zone outer edge. The jump in density and pressure is strong enough to stop the radial drift of particles. In addition, we observe the generation of vortices by the Rossby wave instability (RWI) at the jumps location close to 60 AU. The vortices are steadily generated and destroyed at a cycle of 40 local orbits. The RT results and simulated ALMA observations predict the feasibility to observe such large scale structures appearing in magnetized disks without having a planet.

Young giant planets, which are embedded in a circumstellar disk, will significantly perturb the disk density distribution. This effect can potentially be used as an indirect tracer for planets. We investigate the feasibility of observing planet-induced gaps in circumstellar disks in scattered light. We perform 3D hydrodynamical disk simulations combined with subsequent radiative transfer calculations in scattered light for different star, disk, and planet configurations. The results are compared to those of a corresponding study for the (sub)mm thermal re-emission. The feasibility of detecting planet-induced gaps in scattered light is mainly influenced by the optical depth of the disk and therefore by the disk size and mass. Planet-induced gaps are in general only detectable if the photosphere of the disks is sufficiently disturbed. Within the limitations given by the parameter space here considered, we find that gap detection is possible in the case of disks with masses below $\sim 10^{-4\dots-3} \, \rm M_\odot$. Compared to the disk mass that marks the lower Atacama Large (Sub)Millimeter Array (ALMA) detection limit for the thermal radiation re-emitted by the disk, it is possible to detect the same gap both in re-emission and scattered light only in a narrow range of disk masses around $\sim 10^{-4} \, \rm M_\odot$, corresponding to $16\%$ of cases considered in our study.

Direct imaging surveys for exoplanets commonly exclude binary stars from their target lists, leaving a large part of the overall planet demography unexplored. To address this gap in our understanding of planet formation and evolution, we have launched the first direct imaging survey dedicated to circumbinary planets: SPOTS, the Search for Planets Orbiting Two Stars. In this paper, we discuss the theoretical context, scientific merit, and technical feasibility of such observations, describe the target sample and observational strategy of our survey, and report on the first results from our pilot survey of 26 targets with the VLT NaCo facility. While we have not found any confirmed substellar companions to date, a number of promising candidate companions remain to be tested for common proper motion in upcoming follow-up observations. We also report on the astrometry of the three resolved binaries in our target sample. This pilot survey constitutes a successful proof of concept for our survey strategy and paves the way for a second stage of exploratory observations with VLT SPHERE.

In a molecular cloud dust opacity typically dominates over gas opacity, yet in the vicinities of forming stars dust is depleted, and gas is the sole provider of opacity. In the optically thin circumstellar environments the radiation temperature cannot be assumed to be equal to the gas temperature, hence the two-temperature Planck means are necessary to calculate the radiative equilibrium. By using the two-temperature mean opacity one does obtain the proper equilibrium gas temperature in a circumstellar environment, which is in a chemical equilibrium. A careful consideration of a radiative transfer problem reveals that the equilibrium temperature solution can be degenerate in an optically thin gaseous environment. We compute mean gas opacities based on the publicly available code DFSYNTHE by Kurucz and Castelli. We performed the calculations assuming local thermodynamic equilibrium and an ideal gas equation of state. The values were derived by direct integration of the high-resolution opacity spectrum. We produced two sets of gas opacity tables: Rosseland means and two-temperature Planck means (the tables available via http://cdsweb.u-strasbg.fr/ as well as via http://www.mpia-hd.mpg.de/homes/malygin). For three metallicities [Me/H] = 0.0,+/-0.3 we covered the parameter range 3.48 <= log T_rad[K] <= 4.48 in radiation temperature, 2.8 <= log T_gas[K]} <= 6.0 in gas temperature, and -10 <= log P[dyn cm^-2] <= 6 in gas pressure. We show that in the optically thin circumstellar environment for a given stellar radiation field and local gas density there are several equilibrium gas temperatures possible. We conclude that, in general, equilibrium gas temperature cannot be determined without treating the temperature evolution.

This letter expands the stability criterion for radially stratified, vertically {unstratified} accretion disks incorporating thermal relaxation. We find a linear amplification of epicyclic oscillations in these disks that depends on the effective cooling time, i.e. an overstability. The growth rates of the overstability vanish for both extreme cases, e.g. infinite cooling time and instantaneous cooling, i.e. the adiabatic and fully isothermal cases. However, for thermal relaxation times $\tau$ on the order of the orbital frequency, $\tau\Omega \sim 1$, modes grow at a rate proportional to the square of the Brunt-V\"ais\"al\"a frequency. The overstability is based on epicyclic motions, with the thermal relaxation causing gas to heat while radially displaced inwards, and cool while radially displaced outwards. This causes the gas to have a lower density when moving outwards compared to when it moves inwards, so it feels the outwards directed pressure force more strongly on that leg of the journey. We suggest the term ``Convective Overstability" for the phenomenon that has already been numerically studied in the non-linear regime in the context of amplifying vortices in disks, under the name ``Subcritical Baroclinic Instability". The point of the present paper is to make clear that vortex formation in three-dimensional disks is neither subcritical, i.e. does not need a finite perturbation, nor is it baroclinic in the sense of geophysical fluid dynamics, which requires on vertical shear. We find that Convective Overstability is a linear instability that will operate under a wide range of physical conditions for circumstellar disks.

The opacity due to grains in the envelope of a protoplanet regulates the accretion rate of gas during formation, thus the final bulk composition of planets with primordial H/He is a function of it. Observationally, for exoplanets with known mass and radius it is possible to estimate the bulk composition via internal structure models. We first determine the reduction factor of the ISM grain opacity f_opa that leads to gas accretion rates consistent with grain evolution models. We then compare the bulk composition of synthetic low-mass and giant planets at different f_opa with observations. For f_opa=1 (full ISM opacity) the synthetic low-mass planets have too small radii, i.e., too low envelope masses compared to observations. At f_opa=0.003, the value calibrated with the grain evolution models, synthetic and actual planets occupy similar mass-radius loci. The mean enrichment of giant planets relative to the host star as a function of planet mass M can be approximated as Z_p/Z_star = beta*(M/M_Jup)^alpha. We find alpha=-0.7 independent of f_opa in synthetic populations in agreement with the observational result (-0.71+-0.10). The absolute enrichment level decreases from beta=8.5 at f_opa=1 to 3.5 at f_opa=0. At f_opa=0.003 one finds beta=7.2 which is similar to the observational result (6.3+-1.0). We thus find observational hints that the opacity in protoplanetary atmospheres is much smaller than in the ISM even if the specific value of the grain opacity cannot be constrained here. The result for the enrichment of giant planets helps to distinguish core accretion and gravitational instability. In the simplest picture of core accretion where first a critical core forms and afterwards only gas is added, alpha=-1. If a core accretes all planetesimals inside the feeding zone, alpha=-2/3. The observational result lies between these values, pointing to core accretion as the formation mechanism.

Context: Several recent studies have found that planet migration in adiabatic discs differs significantly from migration in isothermal discs. Depending on the thermodynamic conditions, i.e., the effectiveness of radiative cooling, and the radial surface density profile, planets migrate inward or outward. Clearly, this will influence the semimajor axis - mass distribution of planets as predicted by population synthesis simulations. Aims: Our goal is to study the global effects of radiative cooling, viscous torque desaturation and gap opening as well as stellar irradiation on the tidal migration of a synthetic planet population. Methods: We combine results from several analytical studies and 3D hydrodynamic simulations in a new semi-analytical migration model for the application in our planet population synthesis calculations. Results: We find a good agreement of our model with torques obtained in a 3D radiative hydrodynamic simulations. We find three convergence zones in a typical disc, towards which planets migrate from the in- and outside, affecting strongly the migration behavior of low-mass planets. Interestingly, this leads to slow type II like migration behavior for low-mass planets captured in those zones even without an ad hoc migration rate reduction factor or a yet to be defined halting mechanism. This means that the new prescription of migration including non-isothermal effects makes the preciously widely used artificial migration rate reduction factor obsolete. Conclusions: Outward migration in parts of a disc makes some planets survive long enough to become massive. The convergence zones lead to a potentially observable accumulations of low-mass planets at certain semimajor axes. Our results indicate that further studies of the mass where the corotation torque saturates will be needed since its value has a major impact on the properties of planet populations.

Planets appear to form in environments shaped by the gas flowing through protostellar disks to the central young stars. The flows in turn are governed by orbital angular momentum transfer. In this chapter we summarize current understanding of the transfer processes best able to account for the flows, including magneto-rotational turbulence, magnetically-launched winds, self-gravitational instability and vortices driven by hydrodynamical instabilities. For each in turn we outline the major achievements of the past few years and the outstanding questions. We underscore the requirements for operation, especially ionization for the magnetic processes and heating and cooling for the others. We describe the distribution and strength of the resulting flows and compare with the long-used phenomenological $\alpha$-picture, highlighting issues where the fuller physical picture yields substantially different answers. We also discuss the links between magnetized turbulence and magnetically-launched outflows, and between magnetized turbulence and hydrodynamical vortices. We end with a summary of the status of efforts to detect specific signatures of the flows.

Direct imaging has just started the inventory of the population of gas giant planets on wide-orbits around young stars in the solar neighborhood. Following this approach, we carried out a deep imaging survey in the near-infrared using VLT/NaCo to search for substellar companions. We report here the discovery in L' (3.8 microns) images of a probable companion orbiting at 56 AU the young (10-17 Myr), dusty, and early-type (A8) star HD 95086. This discovery is based on observations with more than a year-time-lapse. Our first epoch clearly revealed the source at 10 sigma while our second epoch lacked good observing conditions hence yielding a 3 sigma detection. Various tests were thus made to rule out possible artifacts. This recovery is consistent with the signal at the first epoch but requires cleaner confirmation. Nevertheless, our astrometric precision suggests the companion to be comoving with the star, with a 3 sigma confidence level. The planetary nature of the source is reinforced by a non-detection in Ks-band (2.18 microns) images according to its possible extremely red Ks - L' color. Conversely, background contamination is rejected with good confidence level. The luminosity yields a predicted mass of about 4-5MJup (at 10-17 Myr) using "hot-start" evolutionary models, making HD 95086 b the exoplanet with the lowest mass ever imaged around a star.

We previously reported the direct detection of a low mass companion at a projected separation of 55+-2 AU around the B9 type star {\kappa} Andromedae. The properties of the system (mass ratio, separation) make it a benchmark for the understanding of the formation and evolution of gas giant planets and brown dwarfs on wide-orbits. We present new angular differential imaging (ADI) images of the Kappa Andromedae system at 2.146 (Ks), 3.776 (L'), 4.052 (NB 4.05) and 4.78 {\mu}m (M') obtained with Keck/NIRC2 and LBTI/LMIRCam, as well as more accurate near-infrared photometry of the star with the MIMIR instrument. We derive a more accurate J = 15.86 +- 0.21, H = 14.95 +- 0.13, Ks = 14.32 +- 0.09 mag for {\kappa} And b. We redetect the companion in all our high contrast observations. We confirm previous contrasts obtained at Ks and L' band. We derive NB 4.05 = 13.0 +- 0.2 and M' = 13.3 +- 0.3 mag and estimate Log10(L/Lsun) = -3.76 +- 0.06. We build the 1-5 microns spectral energy distribution of the companion and compare it to seven PHOENIX-based atmospheric models in order to derive Teff = 1900+100-200 K. Models do not set constrains on the surface gravity. ``Hot-start" evolutionary models predict masses of 14+25-2 MJup based on the luminosity and temperature estimates, and considering a conservative age range for the system (30+120-10 Myr). ``warm-start" evolutionary tracks constrain the mass to M >= 11 MJup. Therefore, the mass of {\kappa} Andromedae b mostly falls in the brown-dwarf regime, due to remaining uncertainties in age and mass-luminosity models. According to the formation models, disk instability in a primordial disk could account for the position and a wide range of plausible masses of {\kappa} And b.

We have performed three-dimensional magneto-hydrodynamical simulations of stellar accretion disks, using the PLUTO code, and studied the accretion of gas onto a Jupiter-mass planet and the structure of the circumplanetary gas flow after opening a gap in the disk. We compare our results with simulations of laminar, yet viscous disks with different levels of an $\alpha$-type viscosity. In all cases, we find that the accretion flow across the surface of the Hill sphere of the planet is not spherically or azimuthally symmetric, and is predominantly restricted to the mid-plane region of the disk. Even in the turbulent case, we find no significant vertical flow of mass into the Hill sphere. The outer parts of the circumplanetary disk are shown to rotate significantly below Keplerian speed, independent of viscosity, while the circumplanetary disk density (therefore the angular momentum) increases with viscosity. For a simulation of a magnetized turbulent disk, where the global averaged alpha stress is $\alpha_{MHD}=10^{-3}$, we find the accretion rate onto the planet to be $\dot{M}\approx 2\times10^{-6}M_{J}yr^{-1}$ for a gap surface density of $12 g cm^{-2}$. This is about a third of the accretion rate obtained in a laminar viscous simulation with equivalent $\alpha$ parameter.

Direct imaging has confirmed the existence of substellar companions on wide orbits. To understand the formation and evolution mechanisms of these companions, the full population properties must be characterized. We aim at detecting giant planet and/or brown dwarf companions around young, nearby, and dusty stars. Our goal is also to provide statistics on the population of giant planets at wide-orbits and discuss planet formation models. We report a deep survey of 59 stars, members of young stellar associations. The observations were conducted with VLT/NaCo at L'-band (3.8 micron). We used angular differential imaging to reach optimal detection performance. A statistical analysis of about 60 % of the young and southern A-F stars closer than 65 pc allows us to derive the fraction of giant planets on wide orbits. We use gravitational instability models and planet population synthesis models following the core-accretion scenario to discuss the occurrence of these companions. We resolve and characterize new visual binaries and do not detect any new substellar companion. The survey's median detection performance reaches contrasts of 10 mag at 0.5as and 11.5 mag at 1as. We find the occurrence of planets to be between 10.8-24.8 % at 68 % confidence level assuming a uniform distribution of planets in the interval 1-13 Mj and 1-1000 AU. Considering the predictions of formation models, we set important constraints on the occurrence of massive planets and brown dwarf companions that would have formed by GI. We show that this mechanism favors the formation of rather massive clump (Mclump > 30 Mj) at wide (a > 40 AU) orbits which might evolve dynamically and/or fragment. For the population of close-in giant planets that would have formed by CA, our survey marginally explore physical separations (<20 AU) and cannot constrain this population.

Recent studies have shown that baroclinic vortex amplification is strongly dependent on certain factors, namely, the global entropy gradient, the efficiency of thermal diffusion and/or relaxation as well as numerical resolution. We conduct a comprehensive study of a broad range and combination of various entropy gradients, thermal diffusion and thermal relaxation time-scales via local shearing sheet simulations covering the parameter space relevant for protoplanetary disks. We measure the Reynolds stresses as a function of our control parameters and see that there is angular momentum transport even for entropy gradients as low as $\beta=-{d\ln s}/{d\ln r}={1}/{2}$, which corresponds to values observed in protoplanetary accretion disks. The amplification-rate of the perturbations, $\Gamma$, appears to be proportional to $\beta^2$ and thus proportional to the square of the \BV ($\Gamma \propto \beta^2 \propto N^2$). The saturation level of Reynolds stresses on the other hand seems to be proportional to $\beta^{1/2}$. This highlights the importance of baroclinic effects even for the low entropy gradients expected in protoplanetary disks.

Recent numerical simulations have shown long-lived axisymmetric sub- and super-Keplerian flows in protoplanetary disks. These zonal flows are found in local as well as global simulations of disks unstable to the magnetorotational instability. This paper covers our study of the strength and lifetime of zonal flows and the resulting long-lived gas over- and underdensities as functions of the azimuthal and radial size of the local shearing box. We further investigate dust particle concentrations without feedback on the gas and without self-gravity. Strength and lifetime of zonal flows increase with the radial extent of the simulation box, but decrease with the azimuthal box size. Our simulations support earlier results that zonal flows have a natural radial length scale of 5-7 gas pressure scale heights. This is the first study that combines three-dimensional MHD simulations of zonal flows and dust particles feeling the gas pressure. The pressure bumps trap particles with $\textrm{St} = 1$ very efficiently. We show that $\textrm{St} = 0.1$ particles (of some centimeters in size if at $5\textrm{AU}$ in an MMSN) reach a hundred-fold higher density than initially. This opens the path for particles of $\textrm{St} = 0.1$ and dust-to-gas ratio of 0.01 or for particles of $\textrm{St} \geq 0.5$ and dust-to-gas ratio $10^{-4}$ to still reach densities that potentially trigger the streaming instability and thus gravoturbulent formation of planetesimals.

Planets are supposed to form in circumstellar disks. The gravitational potential of a planet perturbs the disk and leads to characteristic structures, i.e. spiral waves and gaps, in the disk's density profile. We perform a large-scale parameter study of the observability of these planet-induced structures in circumstellar disks with ALMA. On the basis of HD and MHD simulations, we calculated the disk temperature structure and (sub)mm images of these systems. These were used to derive simulated ALMA images. Because appropriate objects are frequent in Taurus, we focused on a distance of 140pc and a declination of 20{\deg}. The explored range of star-disk-planet configurations consists of 6 HD simulations (including magnetic fields and different planet masses), 9 disk sizes, 15 total disk masses, 6 different central stars, and two different grain size distributions. On almost all scales and in particular down to a scale of a few AU, ALMA is able to trace disk structures induced by planet-disk interaction or by the influence of magnetic fields on the wavelength range between 0.4 and 2.0mm. In most cases, the optimum angular resolution is limited by the sensitivity. However, within the range of typical masses of protoplanetary disks (0.1-0.001Msun) the disk mass has a minor impact on the observability. It is possible to resolve disks down to 2.67e-6Msun and trace gaps induced by a planet with M_p/M_s = 0.001 in disks with 2.67e-4Msun with a signal-to-noise ratio greater than three. The central star has a major impact on the observability of gaps, as well as the considered maximum grainsize of the dust in the disk. In general, it is more likely to trace planet-induced gaps in our MHD models, because gaps are wider in the presence of magnetic fields. We also find that zonal flows resulting from MRI create gap-like structures in the disk's re-emission radiation, which are observable with ALMA.

We review our recent studies demonstrating that the radiation pressure problem in the formation of massive stars can be circumvented via an anisotropy of the thermal radiation field. Such an anisotropy naturally establishes with the formation of a circumstellar disk. The required angular momentum transport within the disk can be provided by developing gravitational torques. Radiative Rayleigh-Taylor instabilities in the cavity regions - as previously suggested in the literature - are not required and are shown to be not occurring in the context of massive star formation.

We investigate the characteristic properties of self-sustained MRI turbulence in low-ionized proto-planetary disks. We study the transition regime between active and dead-zone, performing 3D global non-ideal MHD simulations of stratified disk covering range of magnetic Reynolds number between 2700 < Rm < 6600. We found converged and saturated MRI turbulence for Rm > 5000 with a strength of alpha ~ 0.01. Below Rm < 5000 the MRI starts to decay at the midplane, having Elsasser numbers below one. We find a transition regime between 3300 < Rm < 5000 where the MRI turbulence is still sustained but damped. At around Rm < 3000 the MRI turbulence decays but could reestablished due to the accumulation of toroidal magnetic field or the radial transport of magnetic field from the active region. Below Rm < 3000 the MRI cannot be sustained and is decaying. Here hydro-dynamical motions, like density waves dominate. We observe anti-cyclonic vortices in the transition between dead-zone and active zone.

The research of exoplanets has entered an era in which we characterize extrasolar planets. This has become possible with measurements of radii and luminosities. Meanwhile, radial velocity surveys discover also very low-mass planets. Uniting all this observational data into one coherent picture to better understand planet formation is an important, but difficult undertaking. Our approach is to develop a model which can make testable predictions for all these observational methods. We continue to describe how we have extended our formation model into a self-consistently coupled formation and evolution model. We show how we calculate the internal structure of the solid core and radiogenic heating. We also improve the protoplanetary disk model. Finally, we conduct population synthesis calculations. We present how the planetary mass-radius relationship of planets with primordial H/He envelopes forms and evolves in time. The basic shape of the M-R relation can be understood from the core accretion model. Low-mass planets cannot bind massive envelopes, while super-critical cores necessarily trigger runway gas accretion, leading to "forbidden" zones in the M-R plane. For a given mass, there is a considerable diversity of radii. We compare the synthetic M-R relation with the observed one, finding good agreement for a>0.1 AU. The synthetic radius distribution is characterized by a strong increase towards small R, and a second, lower local maximum at ~1 Jovian radius. The increase towards small radii reflects the increase of the mass function towards low M. The second local maximum is due to the fact that radii are nearly independent of mass for giant planets. A comparison of the synthetic radius distribution with Kepler data shows agreement for R>2 Earth radii, but divergence for smaller radii. We predict that in the next few years, Kepler should find the second, local maximum at ~1 Jovian radius.

A first characterization of many exoplanets has recently been achieved by the observational determination of their radius. For some planets, a measurement of the luminosity has also been possible, with many more directly imaged planets expected in the future. The statistical characterization of exoplanets through their mass-radius and mass-luminosity diagram is thus becoming possible. This is for planet formation and evolution theory of similar importance as the mass-distance diagram. Our aim in this and a companion paper is to extend our formation model into a coupled formation and evolution model. We want to calculate in a self-consistent way all basic characteristics (M,a,R,L) of a planet and use the model for population synthesis calculations. Here we show how we solve the structure equations describing the gaseous envelope not only during the early formation phase, but also during gas runaway accretion, and during the evolutionary phase at constant mass on Gyr timescales. We then study the in situ formation and evolution of Jupiter, the mass-radius relationship of giants, the influence of the core mass on the radius and the luminosity both in the "hot start" and the "cold start" scenario. We put special emphasis on the comparison with other models. We find that our results agree very well with those of more complex models, despite a number of simplifications. The upgraded model yields the most important characteristics of a planet from its beginning as a seed embryo to a Gyr old planet. This is the case for all planets in a synthetic planetary population. Therefore, we can now use self-consistently the statistical constraints coming from all major observational techniques. This is important in a time where different techniques yield constraints on very diverse sub-populations of planets, and where its is challenging to put all these constraints together in one coherent picture.

Context: The global size and spatial distribution of dust is an important ingredient in the structure and evolution of protoplanetary disks and in the formation of larger bodies, such as planetesimals. Aims: We aim to derive simple equations that explain the global evolution of the dust surface density profile and the upper limit of the grain size distribution and which can readily be used for further modeling or for interpreting of observational data. Methods: We have developed a simple model that follows the upper end of the dust size distribution and the evolution of the dust surface density profile. This model is calibrated with state-of-the-art simulations of dust evolution, which treat dust growth, fragmentation, and transport in viscously evolving gas disks. Results: We find very good agreement between the full dust-evolution code and the toy model presented in this paper. We derive analytical profiles that describe the dust-to-gas ratios and the dust surface density profiles well in protoplanetary disks, as well as the radial flux by solid material "rain out", which is crucial for triggering any gravity assisted formation of planetesimals. We show that fragmentation is the dominating effect in the inner regions of the disk leading to a dust surface density exponent of -1.5, while the outer regions at later times can become drift-dominated, yielding a dust surface density exponent of -0.75. Our results show that radial drift is not efficient in fragmenting dust grains. This supports the theory that small dust grains are resupplied by fragmentation due to the turbulent state of the disk.

Context. This is the fourth paper in a series showing the results of planet population synthesis calculations. Aims. Our goal in this paper is to systematically study the effects of important disk properties, namely disk metallicity, mass and lifetime on fundamental planetary properties. Methods. For a large number of protoplanetary disks we calculate a population of planets with our core accretion formation model including planet migration and disk evolution. Results. We find a large number of correlations: Regarding the planetary initial mass function, metallicity, disk mass and disk lifetime have different roles: For high [Fe/H], giant planets are more frequent. For high disk masses, giant planets are more massive. For long disk lifetimes, giant planets are both more frequent and massive. At low metallicities, very massive giant planets cannot form, but otherwise giant planet mass and metallicity are uncorrelated. In contrast, planet masses and disk gas masses are correlated. The sweet spot for giant planet formation is at 5 AU. In- and outside this distance, higher planetesimals surface densities are necessary. Low metallicities can be compensated by high disk masses, and vice versa, but not ad infinitum. At low metallicities, giant planets only form outside the ice line, while at high metallicities, giant planet formation occurs throughout the disk. The extent of migration increases with disk mass and lifetime and usually decreases with metallicity. No clear correlation of metallicity and the semimajor axis of giant planets exists because in low [Fe/H] disks, planets start further out, but migrate more, whereas for high [Fe/H] they start further in, but migrate less. Close-in low mass planets have a lower mean metallicity than Hot Jupiters. Conclusions. The properties of protoplanetary disks are decisive for the properties of planets, and leave many imprints.

Core accretion and disk instability have traditionally been regarded as the two competing possible paths of planet formation. In recent years, evidence have accumulated in favor of core accretion as the dominant mode, at least for close-in planets. However, it might be hypothesized that a significant population of wide planets formed by disk instabilities could exist at large separations, forming an invisible majority. In previous work, we addressed this issue through a direct imaging survey of B2--A0-type stars, and concluded that <30% of such stars form and retain planets and brown dwarfs through disk instability, leaving core accretion as the likely dominant mechanism. In this paper, we extend this analysis to FGKM-type stars by applying a similar analysis to the Gemini Deep Planet Survey (GDPS) sample. The results strengthen the conclusion that substellar companions formed and retained around their parent stars by disk instabilities are rare. Specifically, we find that the frequency of such companions is <8% for FGKM-type stars under our most conservative assumptions, for an outer disk radius of 300 AU, at 99% confidence. Furthermore, we find that the frequency is always <10% at 99% confidence independently of outer disk radius, for any radius from 5 to 500 AU. We also simulate migration at a wide range of rates, and find that the conclusions hold even if the companions move substantially after formation. Hence, core accretion remains the likely dominant formation mechanism for the total planet population, for every type of star from M-type through B-type.

Context: When massive stars exert a radiation pressure onto their environment that is higher than their gravitational attraction, they launch a radiation-pressure-driven outflow. It has been claimed that a radiative Rayleigh-Taylor instability should lead to the collapse of the outflow cavity and foster the growth of massive stars. Aims: We investigate the stability of idealized radiation-pressure-dominated cavities, focusing on its dependence on the radiation transport approach for the stellar radiation feedback. Methods: We compare two different methods for stellar radiation feedback: gray flux-limited diffusion (FLD) and ray-tracing (RT). We also derive simple analytical models to support our findings. Results: Only the FLD cases lead to prominent instability in the cavity shell. The RT cases do not show such instability. The gray FLD method underestimates the opacity at the location of the cavity shell and leads to extended epochs of marginal Eddington equilibrium in the cavity shell, making them prone to the radiative Rayleigh-Taylor instability. In the RT cases, the radiation pressure exceeds gravity by 1-2 orders of magnitude. The radiative Rayleigh-Taylor instability is then consequently suppressed. Conclusions: Treating the stellar irradiation in the gray FLD approximation underestimates the radiative forces acting on the cavity shell. This can lead artificially to situations that are affected by the radiative Rayleigh-Taylor instability. The proper treatment of direct stellar irradiation by massive stars is crucial for the stability of radiation-pressure-dominated cavities.

We investigate the significance of large scale azimuthal, magnetic and velocity modes for the MRI turbulence in accretion disks. We perform 3D global ideal MHD simulations of global stratified proto-planetary disk models. Our domains span azimuthal angles of \pi/4, \pi/2, \pi and 2\pi. We observe up to 100% stronger magnetic fields and stronger turbulence for the restricted azimuthal domain models \pi/2 and \pi/4 compared to the full 2\pi model. We show that for those models, the Maxwell Stress is larger due to strong axisymmetric magnetic fields, generated by the \alpha \Omega dynamo. Large radial extended axisymmetric toroidal fields trigger temporal magnification of accretion stress. All models display a positive dynamo-\alpha in the northern hemisphere (upper disk). The parity is distinct in each model and changes on timescales of 40 local orbits. In model 2\pi, the toroidal field is mostly antisymmetric in respect to the midplane. The eddies of the MRI turbulence are highly anisotropic. The major wavelengths of the turbulent velocity and magnetic fields are between one and two disk scale heights. At the midplane, we find magnetic tilt angles around 8-9 degree increasing up to 12-13 degree in the corona. We conclude that an azimuthal extent of \pi is sufficient to reproduce most turbulent properties in 3D global stratified simulations of magnetised accretion disks.

There has been a long-standing discussion in the literature as to whether core accretion or disk instability is the dominant mode of planet formation. Over the last decade, several lines of evidence have been presented showing that core accretion is most likely the dominant mechanism for the close-in population of planets probed by radial velocity and transits. However, this does not by itself prove that core accretion is the dominant mode for the total planet population, since disk instability might conceivably produce and retain large numbers of planets in the far-out regions of the disk. If this is a relevant scenario, then the outer massive disks of B-stars should be among the best places for massive planets and brown dwarfs to form and reside. In this study, we present high-contrast imaging of 18 nearby massive stars, of which 15 are in the B2--A0 spectral type range and provide excellent sensitivity to wide companions. By comparing our sensitivities to model predictions of disk instability based on physical criteria for fragmentation and cooling, and using Monte-Carlo simulations for orbital distributions, we find that ~85 % of such companions should have been detected in our images on average. Given this high degree of completeness, stringent statistical limits can be set from the null-detection result, even with the limited sample size. We find that <32 % of massive stars form and retain disk instability planets, brown dwarfs and very low-mass stars of <100 Mjup within 300 AU, at 99 % confidence. These results, combined with previous findings in the literature, lead to the conclusion that core accretion is likely the dominant mode of planet formation.

We performed 3D MHD simulations of planet migration in stratified disks using the Godunov code PLUTO, where the disk is turbulent due to the magnetorotational instability. We study the migration for planets with different planet-star mass ratios $q=M_{p}/M_{s}$. In agreement with previous studies, for the low-mass planet cases ($q=5\times10^{-6}$ and $10^{-5}$), migration is dominated by random fluctuations in the torque. For a Jupiter-mass planet $(q=M_{p}/M_{s}=10^{-3}$ for $M_{s}=1M_{\odot})$, we find a reduction of the magnetic stress inside the orbit of the planet and around the gap region. After an initial stage where the torque on the planet is positive, it reverses and we recover migration rates similar to those found in disks where the turbulent viscosity is modelled by an $\alpha$ viscosity. For the intermediate-mass planets ($q=5\times10^{-5}, 10^{-4}$ and $2\times10^{-4}$) we find a new and so far unexpected behavior. In some cases they experience sustained and systematic outwards migration for the entire duration of the simulation. For this case, the horseshoe region is resolved and torques coming from the corotation region can remain unsaturated due to the stresses in the disk. These stresses are generated directly by the magnetic field. The magnitude of the horseshoe drag can overcome the negative Lindblad contribution when the local surface density profile is flat or increasing outwards, which we see in certain locations in our simulations due to the presence of a zonal flow. The intermediate-mass planet is migrating radially outwards in locations where there is a positive gradient of a pressure bump (zonal flow).

We present full 2 Pi global 3-D stratified MHD simulations of accretion disks. We interpret our results in the context of proto-planetary disks. We investigate the turbulence driven by the magneto-rotational instability (MRI) using the PLUTO Godunov code in spherical coordinates with the accurate and robust HLLD Riemann solver. We follow the turbulence for more than 1500 orbits at the innermost radius of the domain to measure the overall strength of turbulent motions and the detailed accretion flow pattern. We find that regions within two scale heights of the midplane have a turbulent Mach number of about 0.1 and a magnetic pressure two to three orders of magnitude less than the gas pressure, while outside three scale heights the magnetic pressure equals or exceeds the gas pressure and the turbulence is transonic, leading to large density fluctuations. The strongest large-scale density disturbances are spiral density waves, and the strongest of these waves has m=5. No clear meridional circulation appears in the calculations because fluctuating radial pressure gradients lead to changes in the orbital frequency, comparable in importance to the stress gradients that drive the meridional flows in viscous models. The net mass flow rate is well-reproduced by a viscous model using the mean stress distribution taken from the MHD calculation. The strength of the mean turbulent magnetic field is inversely proportional to the radius, so the fields are approximately force-free on the largest scales. Consequently the accretion stress falls off as the inverse square of the radius.

We present high-resolution computer simulations of dust dynamics and planetesimal formation in turbulence generated by the magnetorotational instability. We show that the turbulent viscosity associated with magnetorotational turbulence in a non-stratified shearing box increases when going from 256^3 to 512^3 grid points in the presence of a weak imposed magnetic field, yielding a turbulent viscosity of $\alpha\approx0.003$ at high resolution. Particles representing approximately meter-sized boulders concentrate in large-scale high-pressure regions in the simulation box. The appearance of zonal flows and particle concentration in pressure bumps is relatively similar at moderate (256^3) and high (512^3) resolution. In the moderate-resolution simulation we activate particle self-gravity at a time when there is little particle concentration, in contrast with previous simulations where particle self-gravity was activated during a concentration event. We observe that bound clumps form over the next ten orbits, with initial birth masses of a few times the dwarf planet Ceres. At high resolution we activate self-gravity during a particle concentration event, leading to a burst of planetesimal formation, with clump masses ranging from a significant fraction of to several times the mass of Ceres. We present a new domain decomposition algorithm for particle-mesh schemes. Particles are spread evenly among the processors and the local gas velocity field and assigned drag forces are exchanged between a domain-decomposed mesh and discrete blocks of particles. We obtain good load balancing on up to 4096 cores even in simulations where particles sediment to the mid-plane and concentrate in pressure bumps.

The most massive stars can form via standard disk accretion - despite of the radiation pressure generated - due to the fact that the massive accretion disk yields a strong anisotropy in the radiation field, releasing most of the radiation pressure perpendicular to the disk accretion flow. Here, we analyze the self-gravity of the forming circumstellar disk as the potential major driver of the angular momentum transport in such massive disks responsible for the high accretion rates needed for the formation of massive stars. For this purpose, we perform self-gravity radiation hydrodynamics simulations of the collapse of massive pre-stellar cores. The formation and evolution of the resulting circumstellar disk is investigated in 1.) axially symmetric simulations using an alpha-shear-viscosity prescription and 2.) a three-dimensional simulation, in which the angular momentum transport is provided self-consistently by developing gravitational torques in the self-gravitating accretion disk. The simulation series of different strength of the alpha-viscosity shows that the accretion history of the forming star is mostly independent of the alpha-viscosity-parameter. The accretion history of the three-dimensional run driven by self-gravity is more time-dependent than the viscous disk evolution in axial symmetry. The mean accretion rate, i.e. the stellar mass growth, is nearly identical to the alpha-viscosity models. We conclude that the development of gravitational torques in self-gravitating disks around forming massive stars provides a self-consistent mechanism to efficiently transport the angular momentum to outer disk radii. Also the formation of the most massive stars can therefore be understood in the standard accretion disk scenario.

The planetary mass-radius diagram is an observational result of central importance to understand planet formation. We present an updated version of our planet formation model based on the core accretion paradigm which allows to calculate planetary radii and luminosities during the entire formation and evolution of the planets. We first study with it the formation of Jupiter, and compare with previous works. Then we conduct planetary population synthesis calculations to obtain a synthetic mass-radius diagram which we compare with the observed one. Except for bloated Hot Jupiters which can be explained only with additional mechanisms related to their proximity to the star, we find a good agreement of the general shape of the observed and the synthetic mass-radius diagram. This shape can be understood with basic concepts of the core accretion model.

We investigate the radiation pressure feedback in the formation of massive stars in 1, 2, and 3D radiation hydrodynamics simulations of the collapse of massive pre-stellar cores. In contrast to previous research, we consider frequency dependent stellar radiation feedback, resolve the dust sublimation front in the vicinity of the forming star down to 1.27 AU, compute the evolution for several 10^5 yrs covering the whole accretion phase of the forming star, and perform a comprehensive survey of the parameter space. The most fundamental result is that the formation of a massive accretion disk in slowly rotating cores preserves a high anisotropy in the radiation field. The thermal radiation escapes through the optically thin atmosphere, effectively diminishing the radiation pressure feedback onto the accretion flow. Gravitational torques in the self-gravitating disk drive a sufficiently high accretion rate to overcome the residual radiation pressure. Simultaneously, the radiation pressure launches an outflow in the bipolar direction, which grows in angle with time and releases a substantial fraction of the initial core mass from the star-disk system. Summarized, for an initial core mass of 60, 120, 240, and 480 Msol these mechanisms allow the star to grow up to 28.2, 56.5, 92.6, and at least 137.2 Msol respectively.

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.

Turbulence and angular momentum transport in accretion disks remains a topic of debate. With the realization that dead zones are robust features of protoplanetary disks, the search for hydrodynamical sources of turbulence continues. A possible source is the baroclinic instability (BI), which has been shown to exist in unmagnetized non-barotropic disks. We present shearing box simulations of baroclinicly unstable, magnetized, 3D disks, in order to assess the interplay between the BI and other instabilities, namely the magneto-rotational instability (MRI) and the magneto-elliptical instability. We find that the vortices generated and sustained by the baroclinic instability in the purely hydrodynamical regime do not survive when magnetic fields are included. The MRI by far supersedes the BI in growth rate and strength at saturation. The resulting turbulence is virtually identical to an MRI-only scenario. We measured the intrinsic vorticity profile of the vortex, finding little radial variation in the vortex core. Nevertheless, the core is disrupted by an MHD instability, which we identify with the magneto-elliptic instability. This instability has nearly the same range of unstable wavelengths as the MRI, but has higher growth rates. In fact, we identify the MRI as a limiting case of the magneto-elliptic instability, when the vortex aspect ratio tends to infinity (pure shear flow). We conclude that vortex excitation and self-sustenance by the baroclinic instability in protoplanetary disks is viable only in low ionization, i.e., the dead zone. Our results are thus in accordance with the layered accretion paradigm. A baroclinicly unstable dead zone should be characterized by the presence of large-scale vortices whose cores are elliptically unstable, yet sustained by the baroclinic feedback. As magnetic fields destroy the vortices and the MRI outweighs the BI, the active layers are unmodified.

We review the current theoretical understanding how growth from micro-meter sized dust to massive giant planets occurs in disks around young stars. After introducing a number of observational constraints from the solar system, from observed protoplanetary disks, and from the extrasolar planets, we simplify the problem by dividing it into a number of discrete stages which are assumed to occur in a sequential way. In the first stage - the growth from dust to kilometer sized planetesimals - the aerodynamics of the bodies are of central importance. We discuss both a purely coagulative growth mode, as well as a gravoturbulent mode involving a gravitational instability of the dust. In the next stage, planetesimals grow to protoplanets of roughly 1000 km in size. Gravity is now the dominant force. The mass accretion can be strongly non-linear, leading to the detachment of a few big bodies from the remaining planetesimals. In the outer planetary system (outside a few AU), some of these bodies can become so massive that they eventually accrete a large gaseous envelope. This is the stage of giant planet formation, as understood within the core accretion-gas capture paradigm. We also discuss the direct gravitational collapse model where giant planets are thought to form directly via a gravitational fragmentation of the gas disk. In the inner system, protoplanets collide in the last stage - probably after the dispersal of the gaseous disk - in giant impacts until the separations between the remaining terrestrial planets become large enough to allow long term stability. We finish the review with some selected questions.

We present radiation hydrodynamics simulations of the collapse of massive pre-stellar cores. We treat frequency dependent radiative feedback from stellar evolution and accretion luminosity at a numerical resolution down to 1.27 AU. In the 2D approximation of axially symmetric simulations, it is possible for the first time to simulate the whole accretion phase of several 10^5 yr for the forming massive star and to perform a comprehensive scan of the parameter space. Our simulation series show evidently the necessity to incorporate the dust sublimation front to preserve the high shielding property of massive accretion disks. Our disk accretion models show a persistent high anisotropy of the corresponding thermal radiation field, yielding to the growth of the highest-mass stars ever formed in multi-dimensional radiation hydrodynamics simulations. Non-axially symmetric effects are not necessary to sustain accretion. The radiation pressure launches a stable bipolar outflow, which grows in angle with time as presumed from observations. For an initial mass of the pre-stellar host core of 60, 120, 240, and 480 Msol the masses of the final stars formed in our simulations add up to 28.2, 56.5, 92.6, and at least 137.2 Msol respectively.

We present radiation hydrodynamics simulations of the collapse of massive pre-stellar cores. We treat frequency dependent radiative feedback from stellar evolution and accretion luminosity at a numerical resolution down to 1.27 AU. In the 2D approximation of axially symmetric simulations, it is possible for the first time to simulate the whole accretion phase (up to the end of the accretion disk epoch) for the forming massive star and to perform a broad scan of the parameter space. Our simulation series show evidently the necessity to incorporate the dust sublimation front to preserve the high shielding property of massive accretion disks. While confirming the upper mass limit of spherically symmetric accretion, our disk accretion models show a persistent high anisotropy of the corresponding thermal radiation field. This yields to the growth of the highest-mass stars ever formed in multi-dimensional radiation hydrodynamics simulations, far beyond the upper mass limit of spherical accretion. Non-axially symmetric effects are not necessary to sustain accretion. The radiation pressure launches a stable bipolar outflow, which grows in angle with time as presumed from observations. For an initial mass of the pre-stellar host core of 60, 120, 240, and 480 Msun the masses of the final stars formed in our simulations add up to 28.2, 56.5, 92.6, and at least 137.2 Msun respectively.

Planetary bodies form by accretion of smaller bodies. It has been suggested that a very efficient way to grow protoplanets is by accreting particles of size <<km (e.g., chondrules, boulders, or fragments of larger bodies) as they can be kept dynamically cold. We investigate the effects of gas drag on the impact radii and the accretion rates of these particles. As simplifying assumptions we restrict our analysis to 2D settings, a gas drag law linear in velocity, and a laminar disk characterized by a smooth (global) pressure gradient that causes particles to drift in radially. These approximations, however, enable us to cover an arbitrary large parameter space. The framework of the circularly restricted three body problem is used to numerically integrate particle trajectories and to derive their impact parameters. Three accretion modes can be distinguished: hyperbolic encounters, where the 2-body gravitational focusing enhances the impact parameter; three-body encounters, where gas drag enhances the capture probability; and settling encounters, where particles settle towards the protoplanet. An analysis of the observed behavior is presented; and we provide a recipe to analytically calculate the impact radius, which confirms the numerical findings. We apply our results to the sweepup of fragments by a protoplanet at a distance of 5 AU. Accretion of debris on small protoplanets (<50 km) is found to be slow, because the fragments are distributed over a rather thick layer. However, the newly found settling mechanism, which is characterized by much larger impact radii, becomes relevant for protoplanets of ~10^3 km in size and provides a much faster channel for growth.

Context: Radiative feedback plays a crucial role in the formation of massive stars. The implementation of a fast and accurate description of the proceeding thermodynamics in pre-stellar cores and evolving accretion disks is therefore a main effort in current hydrodynamics simulations. Aims: We introduce our newly implemented three-dimensional frequency dependent radiation transport algorithm for hydrodynamics simulations of spatial configurations with a dominant central source. Methods: The module combines the advantage of the speed of an approximate Flux Limited Diffusion (FLD) solver with the high accuracy of a frequency dependent first order ray-tracing routine. Results: We prove the viability of the scheme in a standard radiation benchmark test compared to a full frequency dependent Monte-Carlo based radiative transfer code. The setup includes a central star, a circumstellar flared disk, as well as an envelope. The test is performed for different optical depths. Considering the frequency dependence of the stellar irradiation, the temperature distributions can be described precisely in the optically thin, thick, and irradiated transition regions. Resulting radiative forces onto dust grains are reproduced with high accuracy. The achievable parallel speedup of the method imposes no restriction on further radiative (magneto-) hydrodynamics simulations. Conclusions: The proposed approximate radiation transport method enables frequency dependent radiation hydrodynamics studies of the evolution of pre-stellar cores and circumstellar accretion disks around an evolving massive star in a highly efficient and accurate manner.

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.

The poorly-ionized interior of the protoplanetary disk is the location where dust coagulation processes may be most efficient. However even here, planetesimal formation may be limited by the loss of solid material through radial drift, and by collisional fragmentation of the particles. Our aim is to investigate the possibility that solid particles are trapped at local pressure maxima in the dynamically evolving disk. We perform the first 3-D global non-ideal MHD calculations of the disk treating the turbulence driven by the magneto-rotational instability. The domain contains an inner MRI-active region near the young star and an outer midplane dead zone, with the transition between the two modeled by a sharp increase in the magnetic diffusivity. The azimuthal magnetic fields generated in the active zone oscillate over time, changing sign about every 150 years. We thus observe the radial structure of the `butterfly pattern' seen previously in local shearing-box simulations. The mean magnetic field diffuses from the active zone into the dead zone, where the Reynolds stress nevertheless dominates. The greater total accretion stress in the active zone leads to a net reduction in the surface density, so that after 800 years an approximate steady state is reached in which a local radial maximum in the midplane pressure lies near the transition radius. We also observe the formation of density ridges within the active zone. The dead zone in our models possesses a mean magnetic field, significant Reynolds stresses and a steady local pressure maximum at the inner edge, where the outward migration of planetary embryos and the efficient trapping of solid material are possible.

We present the discovery of a brown dwarf or possible planet at a projected separation of 1.9" = 29 AU around the star GJ 758, placing it between the separations at which substellar companions are expected to form by core accretion (~5 AU) or direct gravitational collapse (typically >100 AU). The object was detected by direct imaging of its thermal glow with Subaru/HiCIAO. At 10-40 times the mass of Jupiter and a temperature of 550-640 K, GJ 758 B constitutes one of the few known T-type companions, and the coldest ever to be imaged in thermal light around a Sun-like star. Its orbit is likely eccentric and of a size comparable to Pluto's orbit, possibly as a result of gravitational scattering or outward migration. A candidate second companion is detected at 1.2" at one epoch.

As accretion in protoplanetary disks is enabled by turbulent viscosity, the border between active and inactive (dead) zones constitutes a location where there is an abrupt change in the accretion flow. The gas accumulation that ensues triggers the Rossby wave instability, that in turn saturates into anticyclonic vortices. It was suggested that the trapping of solids within them leads to a burst of planet formation on very short timescales. We perform two-dimensional global simulations of the dynamics of gas and solids in a non-magnetized thin protoplanetary disk with the Pencil Code. We use multiple particle species of radius 1, 10, 30, and 100 cm, solving for the particles' gravitational interaction by a particle-mesh method. The dead zone is modeled as a region of low viscosity. Adiabatic and locally isothermal equations of state are used. We find that the Rossby wave instability is triggered under a variety of conditions, thus making vortex formation a robust process. Inside the vortices, fast accumulation of solids occurs and the particles collapse into objects of planetary mass in timescales as short as five orbits. Because the drag force is size-dependent, aerodynamical sorting ensues within the vortical motion, and the first bound structures formed are composed primarily of similarly-sized particles. In addition to erosion due to ram pressure, we identify gas tides from the massive vortices as a disrupting agent of formed protoplanetary embryos. We also estimate the collisional velocity history of the particles that compose the most massive embryo by the end of the simulation, finding that the vast majority of them never experienced a collision with another particle at speeds faster than 1 m/s.

The migration of growing protoplanets depends on the thermodynamics of the ambient disc. Standard modelling, using locally isothermal discs, indicate in the low planet mass regime an inward (type-I) migration. Taking into account non-isothermal effects, recent studies have shown that the direction of the type-I migration can change from inward to outward. In this paper we extend previous two-dimensional studies, and investigate the planet-disc interaction in viscous, radiative discs using fully three-dimensional radiation hydrodynamical simulations of protoplanetary accretion discs with embedded planets, for a range of planetary masses. We use an explicit three-dimensional (3D) hydrodynamical code NIRVANA that includes full tensor viscosity. We have added implicit radiation transport in the flux-limited diffusion approximation, and to speed up the simulations significantly we have newly adapted and implemented the FARGO-algorithm in a 3D context. First, we present results of test simulations that demonstrate the accuracy of the newly implemented FARGO-method in 3D. For a planet mass of 20 M_earth we then show that the inclusion of radiative effects yields a torque reversal also in full 3D. For the same opacity law used the effect is even stronger in 3D than in the corresponding 2D simulations, due to a slightly thinner disc. Finally, we demonstrate the extent of the torque reversal by calculating a sequence of planet masses. Through full 3D simulations of embedded planets in viscous, radiative discs we confirm that the migration can be directed outwards up to planet masses of about 33 M_earth. Hence, the effect may help to resolve the problem of too rapid inward migration of planets during their type-I phase.

We employ the PLUTO code for computational astrophysics to assess and compare the validity of different numerical algorithms on simulations of the magneto-rotational instability in 3D accretion disks. In particular we stress on the importance of using a consistent upwind reconstruction of the electro-motive force (EMF) when using the constrained transport (CT) method to avoid the onset of numerical instabilities. We show that the electro-motive force (EMF) reconstruction in the classical constrained transport (CT) method for Godunov schemes drives a numerical instability. The well-studied linear growth of magneto-rotational instability (MRI) is used as a benchmark for an inter-code comparison of PLUTO and ZeusMP. We reproduce the analytical results for linear MRI growth in 3D global MHD simulations and present a robust and accurate Godunov code which can be used for 3D accretion disk simulations in curvilinear coordinate systems.

We study the behavior of magnetorotational turbulence in shearing box simulations with a radial and azimuthal extent up to ten scale heights. Maxwell and Reynolds stresses are found to increase by more than a factor two when increasing the box size beyond two scale heights in the radial direction. Further increase of the box size has little or no effect on the statistical properties of the turbulence. An inverse cascade excites magnetic field structures at the largest scales of the box. The corresponding 10% variation in the Maxwell stress launches a zonal flow of alternating sub- and super-Keplerian velocity. This in turn generates a banded density structure in geostrophic balance between pressure and Coriolis forces. We present a simplified model for the appearance of zonal flows, in which stochastic forcing by the magnetic tension on short time-scales creates zonal flow structures with life-times of several tens of orbits. We experiment with various improved shearing box algorithms to reduce the numerical diffusivity introduced by the supersonic shear flow. While a standard finite difference advection scheme shows signs of a suppression of turbulent activity near the edges of the box, this problem is eliminated by a new method where the Keplerian shear advection is advanced in time by interpolation in Fourier space.

Recently, the existence of geometrically thick dust structures in Active Galactic Nuclei (AGN) has been directly proven with the help of mid-infrared interferometry. The observations are consistent with a two-component model made up of a geometrically thin and warm central disk, surrounded by a colder, fluffy torus component. In an exploratory study, we investigate one possible physical mechanism, which could produce such a structure, namely the effect of stellar feedback from a young nuclear star cluster on the interstellar medium in centres of AGN. The model is realised with the help of the hydrodynamics code TRAMP. We follow the evolution of the interstellar medium by taking discrete mass loss and energy ejection due to stellar processes, as well as optically thin radiative cooling into account. In a post-processing step, we calculate observable quantities (spectral energy distributions and images) with the help of the radiative transfer code MC3D. The interplay between injection of mass, supernova explosions and radiative cooling leads to a two-component structure made up of a cold geometrically thin, but optically thick and very turbulent disk residing in the vicinity of the angular momentum barrier, surrounded by a filamentary structure. The latter consists of cold long radial filaments flowing towards the disk and a hot tenuous medium in between, which shows both inwards and outwards directed motions. This modelling is able to reproduce the range of observed neutral hydrogen column densities of a sample of Seyfert galaxies as well as the relation between them and the strength of the silicate 10 micron spectral feature. Despite being quite crude, our mean Seyfert galaxy model is even able to describe the SEDs of two intermediate type Seyfert galaxies observed with the Spitzer Space Telescope.

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.

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.

Recently, the MID-infrared Interferometric instrument (MIDI) at the VLTI has shown that dust tori in the two nearby Seyfert galaxies NGC 1068 and the Circinus galaxy are geometrically thick and can be well described by a thin, warm central disk, surrounded by a colder and fluffy torus component. By carrying out hydrodynamical simulations with the help of the TRAMP code (Klahr et al. 1999), we follow the evolution of a young nuclear star cluster in terms of discrete mass-loss and energy injection from stellar processes. This naturally leads to a filamentary large scale torus component, where cold gas is able to flow radially inwards. The filaments open out into a dense and very turbulent disk structure. In a post-processing step, we calculate observable quantities like spectral energy distributions or images with the help of the 3D radiative transfer code MC3D (Wolf 2003). Good agreement is found in comparisons with data due to the existence of almost dust-free lines of sight through the large scale component and the large column densities caused by the dense disk.

To treat the problem of growing protoplanetary disc solids across the meter barrier, we consider a very simplified two-component coagulation-fragmentation model that consists of macroscopic boulders and smaller dust grains, the latter being the result of catastrophic collisions between the boulders. Boulders in turn increase their radii by sweeping up the dust fragments. An analytical solution to the dynamical equations predicts that growth by coagulation-fragmentation can be efficient and allow agglomeration of 10-meter-sized objects within the time-scale of the radial drift. These results are supported by computer simulations of the motion of boulders and fragments in 3-D time-dependent magnetorotational turbulence. Allowing however the fragments to diffuse freely out of the sedimentary layer of boulders reduces the density of both boulders and fragments in the mid-plane, and thus also the growth of the boulder radius, drastically. The reason is that the turbulent diffusion time-scale is so much shorter than the collisional time-scale that dust fragments leak out of the mid-plane layer before they can be swept up by the boulders there. Our conclusion that coagulation-fragmentation is not an efficient way to grow across the meter barrier in fully turbulent protoplanetary discs confirms recent results by Brauer, Dullemond, & Henning who solved the coagulation equation in a parameterised turbulence model with collisional fragmentation, cratering, radial drift, and a range of particle sizes. We find that a relatively small population of boulders in a sedimentary mid-plane layer can populate the entire vertical extent of the disc with small grains and that these grains are not first generation dust, but have been through several agglomeration-destruction cycles during the simulations.

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.

The initial stages of planet formation in circumstellar gas discs proceed via dust grains that collide and build up larger and larger bodies (Safronov 1969). How this process continues from metre-sized boulders to kilometre-scale planetesimals is a major unsolved problem (Dominik et al. 2007): boulders stick together poorly (Benz 2000), and spiral into the protostar in a few hundred orbits due to a head wind from the slower rotating gas (Weidenschilling 1977). Gravitational collapse of the solid component has been suggested to overcome this barrier (Safronov 1969, Goldreich & Ward 1973, Youdin & Shu 2002). Even low levels of turbulence, however, inhibit sedimentation of solids to a sufficiently dense midplane layer (Weidenschilling & Cuzzi 1993, Dominik et al. 2007), but turbulence must be present to explain observed gas accretion in protostellar discs (Hartmann 1998). Here we report the discovery of efficient gravitational collapse of boulders in locally overdense regions in the midplane. The boulders concentrate initially in transient high pressures in the turbulent gas (Johansen, Klahr, & Henning 2006), and these concentrations are augmented a further order of magnitude by a streaming instability (Youdin & Goodman 2005, Johansen, Henning, & Klahr 2006, Johansen & Youdin 2007) driven by the relative flow of gas and solids. We find that gravitationally bound clusters form with masses comparable to dwarf planets and containing a distribution of boulder sizes. Gravitational collapse happens much faster than radial drift, offering a possible path to planetesimal formation in accreting circumstellar discs.

This document contains refereed supplementary information for the paper ``Rapid planetesimal formation in turbulent circumstellar discs''. It contains 15 sections (\S1.1 -- \S1.15) that address a number of subjects related to the main paper. We describe in detail the Poisson solver used to find the self-potential of the solid particles, including a linear and a non-linear test problem (\S1.3). Dissipative collisions remove energy from the motion of the particles by collisional cooling (\S1.4), an effect that allows gravitational collapse to occur in somewhat less massive discs (\S1.7). A resolution study of the gravitational collapse of the boulders is presented in \S1.6. We find that gravitational collapse can occur in progressively less massive discs as the grid resolution is increased, likely due to the decreased smoothing of the particle-mesh self-gravity solver with increasing resolution. In \S1.10 we show that it is in good agreement with the Goldreich & Ward (1973) stability analysis to form several-hundred-km-sized bodies, when the analysis is applied to 5 AU and to regions of increased boulder column density. \S11 is devoted to the measurement of random speeds and collision speeds between boulders. We find good agreement between our measurements and analytical theory for the random speeds, but the measured collision speeds are 3 times lower than expected from analytical theory. Higher resolution studies, and an improved analytical theory of collision speeds that takes into account epicyclic motion, will be needed to determine whether collision speeds have converged. In \S1.12 we present models with no magnetic fields. The boulder layer still exhibits strong clumping, due to the streaming instability, if the global solids-to-gas ratio is increased by a factor 3. Gravitational collapse occurs as readily as in magnetised discs.

Millimeter interferometry provides evidence for the presence of mm to cm size "pebbles" in the outer parts of disks around pre-main-sequence stars. The observations suggest that large grains are produced relatively early in disk evolution (< 1 Myr) and remain at large radii for longer periods of time (5 to 10 Myr). Simple theoretical estimates of the radial drift time of solid particles, however, imply that they would drift inward over a time scale of less than 0.1 Myr. In this paper, we address this conflict between theory and observation, using more detailed theoretical models, including the effects of sedimentation, collective drag forces and turbulent viscosity. We find that, although these effects slow down the radial drift of the dust particles, this reduction is not sufficient to explain the observationally determined long survival time of mm/cm-sized grains in protoplanetary disks. However, if for some reason the gas to dust ratio in the disk is reduced by at least a factor of 20 from the canonical value of 100 (for instance through photoevaporation of the gas), then the radial drift time scales become sufficiently large to be in agreement with observations.

We study the effect of an imposed vertical magnetic field on the turbulent mass diffusion properties of magnetorotational turbulence in protoplanetary discs. It is well-known that the effective viscosity generated by the turbulence depends strongly on the magnitude of such an external field. In this letter we show that the turbulent diffusion of the flow also grows, but that the diffusion coefficient does not rise with increasing vertical field as fast as the viscosity does. The vertical Schmidt number, i.e. the ratio between viscosity and vertical diffusion, can be close to 20 for high field magnitudes, whereas the radial Schmidt number is increased from below unity to around 3.5. Our results may have consequences for the interpretation of observations of dust in protoplanetary discs and for chemical evolution modelling of these discs.

We perform numerical simulations of a disc-planet system using various grid-based and smoothed particle hydrodynamics (SPH) codes. The tests are run for a simple setup where Jupiter and Neptune mass planets on a circular orbit open a gap in a protoplanetary disc during a few hundred orbital periods. We compare the surface density contours, potential vorticity and smoothed radial profiles at several times. The disc mass and gravitational torque time evolution are analyzed with high temporal resolution. There is overall consistency between the codes. The density profiles agree within about 5% for the Eulerian simulations while the SPH results predict the correct shape of the gap although have less resolution in the low density regions and weaker planetary wakes. The disc masses after 200 orbital periods agree within 10%. The spread is larger in the tidal torques acting on the planet which agree within a factor 2 at the end of the simulation. In the Neptune case the dispersion in the torques is greater than for Jupiter, possibly owing to the contribution from the not completely cleared region close to the planet.

We perform numerical simulations of the Kelvin-Helmholtz instability in the mid-plane of a protoplanetary disk. A two-dimensional corotating slice in the azimuthal--vertical plane of the disk is considered where we include the Coriolis force and the radial advection of the Keplerian rotation flow. Dust grains, treated as individual particles, move under the influence of friction with the gas, while the gas is treated as a compressible fluid. The friction force from the dust grains on the gas leads to a vertical shear in the gas rotation velocity. As the particles settle around the mid-plane due to gravity, the shear increases, and eventually the flow becomes unstable to the Kelvin-Helmholtz instability. The Kelvin-Helmholtz turbulence saturates when the vertical settling of the dust is balanced by the turbulent diffusion away from the mid-plane. The azimuthally averaged state of the self-sustained Kelvin-Helmholtz turbulence is found to have a constant Richardson number in the region around the mid-plane where the dust-to-gas ratio is significant. Nevertheless the dust density has a strong non-axisymmetric component. We identify a powerful clumping mechanism, caused by the dependence of the rotation velocity of the dust grains on the dust-to-gas ratio, as the source of the non-axisymmetry. Our simulations confirm recent findings that the critical Richardson number for Kelvin-Helmholtz instability is around unity or larger, rather than the classical value of 1/4

We study axisymmetric models of layered protoplanetary discs taking radiative transfer effects into account, and allowing for a residual viscosity in the dead zone. We also explore the effect of different viscosity prescriptions. In addition to the ring instability reported in the first paper of the series we find an oscillatory instability of the dead zone, accompanied by variations of the accretion rate onto the central star. We provide a simplified analytical description explaining the mechanism of the oscillations. Finally, we find that the residual viscosity enables stationary accretion in large regions of layered discs. Based on results obtained with the help of a simple 1-D hydrocode we identify these regions, and discuss conditions in which layered discs can give rise to FU~Orionis phenomena.

We study the formation of a giant gas planet by the core--accretion gas--capture process, with numerical simulations, under the assumption that the planetary core forms in the center of an anti-cyclonic vortex. The presence of the vortex concentrates particles of centimeter to meter size from the surrounding disk, and speeds up the core formation process. Assuming that a planet of Jupiter mass is forming at 5 AU from the star, the vortex enhancement results in considerably shorter formation times than are found in standard core--accretion gas--capture simulations. Also, formation of a gas giant is possible in a disk with mass comparable to that of the minimum mass solar nebula.

Debris rings of dust are found around young luminous stars such as HR4796A and HD141569. Some of these entities have sharp edges and gaps which have been interpreted as evidence for the presence of shepherding and embedded planets. Here we show that gaps and sharp edges in the debris disks of dust can also be spontaneously self generated if they are embedded in optically thin regions of gaseous disks. This clumping instability arises in regions where an enhancement in the dust density leads to local gas temperature and pressure increases. Consequently, the relative motion between the gas and the dust is modified. The subsequent hydrodynamic drag on the dust particles leads to further enhancement of their concentration. We show that this process is linearly unstable and leads to the formation of ring-like structures within the estimated life time of such young objects. Once the gas is removed (e.g. by photo evaporation) the structures are ``frozen'' and will persist, even when the gas might not be observable anymore.

We study the evolution of an embedded protoplanet in a circumstellar disk using the 3D-Radiation Hydro code TRAMP, and treat the thermodynamics of the gas properly in three dimensions. The primary interest of this work lies in the demonstration and testing of the numerical method. We show how far numerical parameters can influence the simulations of gap opening. We study a standard reference model under various numerical approximations. Then we compare the commonly used locally isothermal approximation to the radiation hydro simulation using an equation for the internal energy. Models with different treatments of the mass accretion process are compared. Often mass accumulates in the Roche lobe of the planet creating a hydrostatic atmosphere around the planet. The gravitational torques induced by the spiral pattern of the disk onto the planet are not strongly affected in the average magnitude, but the short time scale fluctuations are stronger in the radiation hydro models. An interesting result of this work lies in the analysis of the temperature structure around the planet. The most striking effect of treating the thermodynamics properly is the formation of a hot pressure--supported bubble around the planet with a pressure scale height of H/R ~ 0.5 rather than a thin Keplerian circumplanetary accretion disk. We also observe an outflow of gas above and below the planet during the gap opening phase.

We explore the effect of magnetorotational turbulence on the dynamics and concentrations of boulders in local box simulations of a sub-Keplerian protoplanetary disc. The solids are treated as particles each with an independent space coordinate and velocity. We find that the turbulence has two effects on the solids. 1) Meter and decameter bodies are strongly concentrated, locally up to a factor 100 times the average dust density, whereas decimeter bodies only experience a moderate density increase. The concentrations are located in large scale radial gas density enhancements that arise from a combination of turbulence and shear. 2) For meter-sized boulders, the concentrations cause the average radial drift speed to be reduced by 40%. We find that the densest clumps of solids are gravitationally unstable under physically reasonable values for the gas column density and for the dust-to-gas ratio due to sedimentation. We speculate that planetesimals can form in a dust layer that is not in itself dense enough to undergo gravitational fragmentation, and that fragmentation happens in turbulent density fluctuations in this sublayer.

We measure the turbulent diffusion coefficient of dust grains embedded in magnetorotational turbulence in a protoplanetary disc directly from numerical simulations and compare it to the turbulent viscosity of the flow. The simulations are done in a local coordinate frame comoving with the gas in Keplerian rotation. Periodic boundary conditions are used in all directions, and vertical gravity is not applied to the gas. Using a two-fluid approach, small dust grains of various sizes (with friction times up to $\varOmega_0 \tau_{\rm f}=0.02$) are allowed to move under the influence of friction with the turbulent gas. We measure the turbulent diffusion coefficient of the dust grains by applying an external sinusoidal force field acting in the vertical direction on the dust component only. This concentrates the dust around the mid-plane of the disc, and an equilibrium distribution of the dust density is achieved when the vertical settling is counteracted by the turbulent diffusion away from the mid-plane. Comparing with analytical expressions for the equilibrium concentration we deduce the vertical turbulent diffusion coefficient. The vertical diffusion coefficient is found to be lower than the turbulent viscosity and to have an associated vertical diffusion Prandtl number of about 1.5. A similar radial force field also allows us to measure the radial turbulent diffusion coefficient. We find a radial diffusion Prandtl number of about 0.85 and also find that the radial turbulent diffusion coefficient is around 70% higher than the vertical. We also find evidence for trapping of dust grains of intermediate friction time in turbulent eddies.

In this work we use the radiation hydrodynamic code TRAMP to perform a two-dimensional axially symmetric model of the layered disc. Using this model we follow the accumulation of mass in the dead zone due to the radially varying accretion rate. We found a new type of instability which causes the dead zone to split into rings. This "ring instability" works due to the positive feedback between the thickness of the dead zone and the mass accumulation rate. We give an analytical description of this instability, taking into account non-zero thickness of the dead zone and deviations from the Keplerian rotational velocity. The analytical model agrees reasonably well with results of numerical simulations. Finally, we speculate about the possible role of the ring instability in protoplanetary discs and in the formation of planets.

This paper contains a local linear stability analysis for accretion disks under the influence of a global radial entropy gradient beta = - d log T / d log r for constant surface density. Numerical simulations suggested the existence of an instability in two- and three-dimensional models of the solar nebula. The present paper tries to clarify, quantify, and explain such a global baroclinic instability for two-dimensional flat accretion disk models. As a result linear theory predicts a transient linear instability that will amplify perturbations only for a limited time or up to a certain finite amplification. This can be understood as a result of the growth time of the instability being longer than the shear time which destroys the modes which are able to grow. So only non-linear effects can lead to a relevant amplification. Nevertheless, a lower limit on the entropy gradient ~beta = 0.22 for the transient linear instability is derived, which can be tested in future non-linear simulations. This would help to explain the observed instability in numerical simulations as an ultimate result of the transient linear instability, i.e. the Global Baroclinic Instability.

In this paper we present the global baroclinic instability as a source for vigorous turbulence leading to angular momentum transport in Keplerian accretion disks. We show by analytical considerations and three-dimensional radiation hydro simulations that, in particular, protoplanetary disks have a negative radial entropy gradient, which makes them baroclinic. Two-dimensional numerical simulations show that a baroclinic flow is unstable and produces turbulence. These findings are tested for numerical effects by performing a simulation with a barotropic initial condition which shows that imposed turbulence rapidly decays. The turbulence in baroclinic disks transports angular momentum outward and creates a radially inward bound accretion of matter. Potential energy is released and excess kinetic energy is dissipated. Finally the reheating of the gas supports the radial entropy gradient, forming a self consistent process. We measure accretion rates in our 2D and 3D simulations of dotM= - 1E-9 - 1E-7 Msun/yr and viscosity parameters of alpha = 1E-4 - 1E-2, which fit perfectly together and agree reasonably with observations. The turbulence creates pressure waves, Rossby waves, and vortices in the R-phi plane of the disk. We demonstrate in a global simulation that these vortices tend to form out of little background noise and to be long-lasting features, which have already been suggested to lead to the formation of planets.

We investigate the possibility of mapping large-scale anti-cyclonic vortices, resulting from a global baroclinic instability, as pre-cursors of planet formation in proto-planetary disks with the planned Atacama Large Millimeter Array (ALMA). On the basis of three-dimensional radiative transfer simulations, images of a hydrodynamically calculated disk are derived which provide the basis for the simulation of ALMA. We find that ALMA will be able to trace the theoretically predicted large-scale anti-cyclonic vortex and will therefore allow testing of existing models of this very early stage of planet formation in circumstellar disks.

There have been several model analyses of the near and mid IR flux from the circumstellar ring around HR4796A. In the vicinity of a young star, the possibility that the dust ring is embedded within a residual protostellar gas disk cannot be ruled out. In a gas-rich environment, larger sizes ($>100 \mu m$) are needed for the particles to survive the radiative blow out. The total dust mass required to account for the IR flux is $< 10^{-1} M_\oplus$. The combined influence of gas and stellar radiation may also account for the observed sharp inner boundary and rapidly fading outer boundary of the ring. The pressure gradient induced by a small (10%) amplitude variation in the surface density distribution of a low-mass gaseous disk would be sufficient to modify the rotation speed of the gas.