The solid accretion rate, necessary to grow gas giant planetary cores within the disk lifetime, has been a major constraint for theories of planet formation. We tested the solid accretion rate efficiency on planetary cores of different masses embedded in their birth disk, by means of 3D radiation-hydrodynamics, where we followed the evolution of a swarm of embedded solids of different sizes. We found that using a realistic equation of state and radiative cooling, the disk at 5 au is able to cool efficiently and reduce its aspect ratio. As a result, the pebble isolation mass is reached before the core grows to 10 Earth masses, stopping efficiently the pebble flux and creating a transition disk. Moreover, the reduced isolation mass halts the solid accretion before the core reaches the critical mass, leading to a barrier to giant planet formation, and it explains the large abundance of super-Earth planets in the observed population.

A large fraction of stars is found to be part of binary or higher-order multiple systems. The ubiquity of planets found around single stars raises the question if and how planets in binary systems may form. Protoplanetary disks are the birthplaces of planets, and their characterization is crucial in order to understand the planet formation process. Our aim is to characterize the morphology of the GG Tau A disk, one of the largest and most massive circumbinary disks, and trace evidence for binary-disk interactions. We obtained observations in polarized scattered light of GG Tau A using the SPHERE/IRDIS instrument in the H-band filter. We analyze the observed disk morphology and substructures. We run 2D hydrodynamical models simulating the evolution of the circumbinary ring over the lifetime of the disk. The disk, as well as the cavity and the inner region are highly structured with several shadowed regions, spiral structures, and streamer-like filaments, some of them detected for the first time. The streamer-like filaments appear to connect the outer ring with the northern arc. Their azimuthal spacing suggests that they may be generated by periodic perturbations by the binary, tearing off material from the inner edge of the outer disk once during each orbit. By comparing observations to hydrodynamical simulations we find that the main features, in particular the gap size, as well as the spiral and streamer filaments, can be qualitatively explained by the gravitational interactions of a binary with semi-major axis of $\sim$35 au on an orbit coplanar with the circumbinary ring.

Hydrodynamical simulations of planet-disk interactions suggest that planets may be responsible for a number of the sub-structures frequently observed in disks in both scattered light and dust thermal emission. Despite the ubiquity of these features, direct evidence of planets embedded in disks and of the specific interaction features like spiral arms within planetary gaps still remain rare. In this study we discuss recent observational results in the context of hydrodynamical simulations in order to infer the properties of a putative embedded planet in the cavity of a transition disk. We imaged the transition disk SR 21 in H-band in scattered light with SPHERE/IRDIS and in thermal dust emission with ALMA band 3 (3mm) observations at a spatial resolution of 0.1". We combine these datasets with existing band 9 (430um) and band 7 (870um) ALMA continuum data. The Band 3 continuum data reveals a large cavity and a bright ring peaking at 53 au strongly suggestive of dust trapping.The ring shows a pronounced azimuthal asymmetry, with a bright region in the north-west that we interpret as a dust over-density. A similarly-asymmetric ring is revealed at the same location in polarized scattered light, in addition to a set of bright spirals inside the mm cavity and a fainter spiral bridging the gap to the outer ring. These features are consistent with a number of previous hydrodynamical models of planet-disk interactions, and suggest the presence of a ~1 MJup planet at 44 au and PA=11{\deg}. This makes SR21 the first disk showing spiral arms inside the mm cavity, as well as one for which the location of a putative planet can be precisely inferred. With the location of a possible planet being well-constrained by observations, it is an ideal candidate for follow-up observations to search for direct evidence of a planetary companion still embedded in its disk.

Recent ALMA observations revealed concentric annular structures in several young class-II objects. In an attempt to produce the rings and gaps in some of these systems, they have been modeled numerically with a single embedded planet assuming a locally isothermal equation of state. This is often justified by observations targeting the irradiation-dominated outer regions of disks (approximately 100 au). We test this assumption by conducting hydrodynamics simulations of embedded planets in thin locally isothermal and radiative disks that mimic the systems HD 163296 and AS 209 in order to examine the effect of including the energy equation in a seemingly locally isothermal environment as far as planet-disk interaction is concerned. We find that modeling such disks with an ideal equation of state makes a difference in terms of the number of produced rings and the spiral arm contrast in the disk. Locally isothermal disks produce sharper annular or azimuthal features and overestimate a single planet's gap-opening capabilities by producing multiple gaps. In contrast, planets in radiative disks carve a single gap for typical disk parameters. Consequently, for accurate modeling of planets with semimajor axes up to about 100 au, radiative effects should be taken into account even in seemingly locally isothermal disks. In addition, for the case of AS 209, we find that the primary gap is significantly different between locally isothermal and radiative models. Our results suggest that multiple planets are required to explain the ring-rich structures in such systems.

The most accurate method for modelling planetary migration and hence the formation of resonant systems is using hydrodynamical simulations. Usually, the force (torque) acting on a planet is calculated using the forces from the gas disc and the star, while the gas accelerations are computed using the pressure gradient, the star, and the planet's gravity, ignoring its own gravity. For the non-migrating the neglect of the disc gravity results in a consistent torque calculation while for the migrating case it is inconsistent. We aim to study how much this inconsistent torque calculation can affect the final configuration of a two-planet system. Our focus will be on low-mass planets because most of the multi-planetary systems, discovered by the Kepler survey, have masses around 10 Earth masses. Performing hydrodynamical simulations of planet-disc interaction, we measure the torques on non-migrating and migrating planets for various disc masses as well as density and temperature slopes with and without considering the disc self-gravity. Using this data, we find a relation that quantifies the inconsistency, use it in an N-body code, and perform an extended parameter study modelling the migration of a planetary system with different planet mass ratios and disc surface densities, in order to investigate the impact of the torque inconsistency on the architecture of the planetary system. Not considering disc self-gravity produces an artificially larger torque on the migrating planet that can result in tighter planetary systems. The deviation of this torque from the correct value is larger in discs with steeper surface density profiles. In hydrodynamical modelling of multi-planetary systems, it is crucial to account for the torque correction, otherwise the results favour more packed systems.

Planetary migration is a key link between planet formation models and observed exoplanet statistics. So far the theory of migration has focused on the interaction of planets with an inviscid or viscously evolving disk. Turbulent viscosity is thought to be the main driver of disk evolution and is known to affect the migration process. Recently, the topic of wind-driven accretion is experiencing a renaissance, as evidence is mounting that PPDs may be less turbulent than previously thought, and 3-D non-ideal MHD modeling of the wind-launching process is maturing. Aim: We wish to investigate how wind-driven accretion may affect migration. We aim for a qualitative exploration of the main effects, rather than a quantitative prediction. Methods: We perform 2-D hydrodynamic simulations with the FARGO3D code in the $(r,\phi)$-plane. The vertical coordinate and the launching of the wind are not treated explicitly. Instead, the torque of the wind onto the disk is treated using a simple 2-parameter formula treating the wind mass loss rate and the lever arm. Results: We find that the wind-driven accretion process has a different way of replenishing the co-orbital region than the viscous accretion. The former always injects mass from the outer edge of the co-orbital region and removes mass from the inner edge, while the latter injects or removes mass from the co-orbital region depending on the radial density gradients in the disk. The migration behavior can differ very much and under certain conditions it can drive rapid type-III-like outward migration. We derive an analytic expression for the parameters under which this outward migration occurs. Conclusion: If wind-driven accretion plays a role in the secular evolution of PPDs, migration studies have to include this process as well, because it can strongly affect the resulting migration rate and direction.

Context: Planets in accretion disks can excite spiral shocks, and---if massive enough---open gaps in their vicinity. Both of these effects can influence the overall disk thermal structure. Aims: We model planets of different masses and semimajor axes in disks of various viscosities and accretion rates to examine their impact on disk thermodynamics and highlight the mutable, non-axisymmetric nature of icelines in systems with massive planets. Methods: We conduct a parameter study using numerical hydrodynamics simulations where we treat viscous heating, thermal cooling and stellar irradiation as additional source terms in the energy equation, with some runs including radiative diffusion. Our parameter space consists of a grid containing different combinations of planet and disk parameters. Results: Both gap opening and shock heating can displace the iceline, with the effects being amplified for massive planets in optically thick disks. The gap region can split an initially hot (T>170 K) disk into a hot inner disk and a hot ring just outside of the planet's location, while shock heating can reshape the originally axisymmetric iceline into water-poor islands along spirals. We also find that radiative diffusion does not alter the picture significantly in this context. Conclusions: Shock heating and gap opening by a planet can effectively heat up optically thick disks and in general move and/or reshape the water iceline. This can affect the gap structure and migration torques. It can also produce azimuthal features that follow the trajectory of spiral arms, creating hot zones, "islands" of vapor and ice around spirals which could affect the accretion or growth of icy aggregates.

The recent detection of the third planet in Kepler-47 has shown that binary stars can host several planets in circumbinary orbits. To understand the evolution of such systems we have performed two-dimensional hydrodynamic simulations of the circumbinary disc with two embedded planets for several Kepler systems. In two cases, Kepler-47 and -413, the planets are captured in a 1:1 mean-motion resonance at the planet parking position near the inner edge of the disc. The orbits are fully aligned, have mean eccentricities of about 0.25 to 0.30, and the planets are entangled in a horseshoe type of motion. Subsequent n-body simulations without the disc show that the configurations are stable. Our results point to the existence of a new class of stable resonant orbits around binary stars. It remains to be seen if such orbits exist in reality.

As of today ten circumbinary planets orbiting solar type main sequence stars have been discovered. Nearly all of them orbit around the central binary very closely to the region of instability where it is difficult to form them in situ. It is assumed that they formed further out and migrated to their observed position. We extend previous studies to a more realistic thermal disc structure and determine what parameter influence the final parking location of a planet around a binary star. We perform two-dimensional numerical simulations of viscous accretion discs around a central binary that include viscous heating and radiative cooling from the disc surfaces. We vary the binary eccentricity as well as disc viscosity and mass. Concerning the disc evolution we find that it can take well over 100000 binary orbits until an equilibrium state is reached. As seen previously, we find that the central cavity opened by the binary becomes eccentric and precesses slowly in a prograde sense. Embedded planets migrate to the inner edge of the disc. In cases of lower disc viscosity they migrate further in maintaining a circular orbit, while for high viscosity they are parked further out on an eccentric orbit. The final location of an embedded planet is linked to its ability to open a gap in the disc. Gap opening planets separate inner from outer disc, preventing eccentricity excitation in the latter and making it more circular. This allows embedded planets to migrate closer to the binary, in agreement with the observations. The necessary condition for gap opening and the final planet position depend on the planet mass and disc viscosity.

The discovery of planets in close orbits around binary stars raises questions about their formation. It is believed that these planets formed in the outer regions of the disc and then migrated through planet-disc interaction to their current location. Considering five different systems (Kepler-16, -34, -35, -38 and -413) we model planet migration through the disc, with special focus on the final orbital elements of the planets. We investigate how the final orbital parameters are influenced by the disc and planet masses. Using 2D, locally isothermal, viscous hydrodynamical simulations we first model the disc dynamics for all five systems, followed by a study of the migration properties of embedded planets with different masses. For all systems we find that the discs become eccentric and precess slowly. We confirm the bifurcation feature in the precession period -- gap-size diagram for different binary mass ratios. The Kepler-16, -35, -38 and -413 systems lie on the lower branch and Kepler-34 on the upper one. For systems with small binary eccentricity we find a new non-monotonic, loop-like feature. In all systems the planets migrate to the inner edge of the disc cavity. Depending on the planet-disc mass ratio we observe two different regimes. Massive planets can significantly alter the disc structure by compressing and circularising the inner cavity and they remain on nearly circular orbits. Lower mass planets are strongly influenced by the disc, their eccentricity is exited to high values, and their orbits are aligned with the inner disc in a state of apsidal corotation. In our simulations the final locations of the planets are typically too large with respect to the observations because of the large inner gaps of the discs. The migrating planets in the most eccentric discs (around Kepler-34 and -413) show the largest final eccentricity in agreement with the observations.

Orbital mean motion resonances in planetary systems originate from dissipative processes in disk-planet interactions that lead to orbital migration. In multi-planet systems that host giant planets, the perturbation of the protoplanetary disk strongly affects the migration of companion planets. By studying the well-characterized resonant planetary system around GJ 876 we aim to explore which effects shape disk-driven migration in such a multi-planet system to form resonant chains. We modelled the orbital migration of three planets embedded in a protoplanetary disk using two-dimensional locally isothermal hydrodynamical simulations. We performed a parameter study by varying the disk thickness, $\alpha$ viscosity, mass as well as the initial position of the planets. Moreover, we have analysed and compared simulations with various boundary conditions at the disk's inner rim. We find that due to the high masses of the giant planets, substantial eccentricity can be excited in the disk. This results in large variations of the torque acting on the outer lower-mass planet, which we attribute to a shift of Lindblad and corotation resonances due to disk eccentricity. Depending on disk parameters, the migration of the outer planet can be stopped in a non-resonant state. In other models, the outer planet is able to open a partial gap and to circularize the disk again, later entering a 2:1 resonance with the most massive planet in the system to complete the observed Laplace resonance. Disk-mediated interactions between planets due to spiral waves and excitation of disk eccentricity cause deviations from smooth inward migration of exterior lower mass planets. Self-consistent modelling of the disk-driven migration of multi-planet systems is thus mandatory. Our results are compatible with a late migration of the outermost planet into the resonant chain, when the giant planet pair already is in resonance.

We study the evolution of circumbinary disks under the gravitational influence of the binary using two-dimensional hydrodynamical simulations to investigate the impact of disk and binary parameters on the dynamical aspects of the disk. To distinguish between physical and numerical effects we apply three hydrodynamical codes. First we analyse in detail numerical issues concerning the conditions at the boundaries and grid resolution. We then perform a series of simulations with different binary (eccentricity, mass ratio) and disk parameters (viscosity, aspect ratio) starting from a reference model with Kepler-16 parameters. Concerning the numerical aspects we find that the inner grid radius must be of the order of the binary semi-major axis, with free outflow conditions applied such that mass can flow onto the central binary. A closed inner boundary leads to unstable evolutions. We find that the inner disk turns eccentric and precesses for all investigated physical parameters. The precession rate is slow with periods ($T_\mathrm{prec}$) starting at around 500 binary orbits ($T_\mathrm{bin}$) for high viscosity and large $H/R$ where the inner hole is smaller and more circular. Reducing $\alpha$ and $H/R$ increases the gap size and $T_\mathrm{prec}$ reaches 2500 $T_\mathrm{bin}$. For varying binary mass ratios $q_\mathrm{bin}$ the gap size remains constant whereas $T_\mathrm{prec}$ decreases for increasing $q_\mathrm{bin}$. For varying binary eccentricities $e_\mathrm{bin}$ we find two separate branches in the gap size and eccentricity diagram. The bifurcation occurs at around $e_\mathrm{crit} \approx 0.18$ where the gap is smallest with the shortest $T_\mathrm{prec}$. For $e_\mathrm{bin}$ smaller and larger than $e_\mathrm{crit}$ the gap size and $T_\mathrm{prec}$ increase. Circular binaries create the most eccentric disks.

The growth process of proto-planets can be sped-up by accreting a large number of solid, pebble-sized objects that are still present in the protoplanetary disc. It is still an open question on how efficient this process works in realistic turbulent discs. Here, we investigate the accretion of pebbles in turbulent discs that are driven by the purely hydrodynamical vertical shear instability (VSI). For this purpose, we perform global three-dimensional simulations of locally isothermal, VSI turbulent discs with embedded protoplanetary cores from 5 to 100 $M_\oplus$ that are placed at 5.2 au distance from the star. In addition, we follow the evolution of a swarm of embedded pebbles of different size under the action of drag forces between gas and particles in this turbulent flow. Simultaneously, we perform a set of comparison simulations for laminar viscous discs where the particles experience stochastic kicks. For both cases, we measure the accretion rate onto the cores as a function of core mass and Stokes number ($\tau_s$) of the particles and compare it to recent MRI turbulence simulations. Overall the dynamic is very similar for the particles in the VSI turbulent disc and the laminar case with stochastic kicks. For the small mass planets (i.e. 5 and 10 $M_\oplus$), well-coupled particles with $\tau_s = 1$, which have a size of about one meter at this location, we find an accretion efficiency (rate of particles accreted over drifting inward) of about 1.6-3%. For smaller and larger particles this efficiency is higher. However, the fast inward drift for $\tau_s = 1$ particles makes them the most effective for rapid growth, leading to mass doubling times of about 20,000 yr. For masses between 10 and 30 $M_\oplus$ the core reaches the pebble isolation mass and the particles are trapped at the pressure maximum just outside of the planet, shutting off further particle accretion.

The formation of satellite systems in circum-planetary discs is considered to be similar to the formation of rocky planets in a proto-planetary disc, especially Super-Earths. Thus, it is possible to use systems with large satellites to test formation theories that are also applicable to extrasolar planets. Furthermore, a better understanding of the origin of satellites might yield important information about the environment near the growing planet during the last stages of planet formation. In this work we investigate the formation and migration of the Jovian satellites through N-body simulations. We simulated a massive, static, low viscosity, circum-planetary disc in agreement with the minimum mass sub-nebula model prescriptions for its total mass. In hydrodynamic simulations we found no signs of gaps, therefore type II migration is not expected. Hence, we used analytic prescriptions for type I migration, eccentricity and inclination damping, and performed N-body simulations with damping forces added. Detailed parameter studies showed that the number of final satellites is strong influenced by the initial distribution of embryos, the disc temperature, and the initial gas density profile. For steeper initial density profiles it is possible to form systems with multiple satellites in resonance while a flatter profile favours the formation of satellites close to the region of the Galilean satellites. We show that the formation of massive satellites such as Ganymede and Callisto can be achieved for hotter discs with an aspect ratio of H/r ~ 0.15 for which the ice line was located around 30 R_J.

Using the HPC ressources of the state of Baden-W\"urttemberg, we modelled for the first time the luminous burst from a young massive star by accretion of material from its close environment. We found that the surroundings of young massive stars are shaped as a clumpy disk whose fragments provoke outbursts once they fall onto the protostar and concluded that similar strong luminous events observed in high-mass star forming regions may be a signature of the presence of such disks.

The surroundings of massive protostars constitute an accretion disc which has numerically been shown to be subject to fragmentation and responsible for luminous accretion-driven outbursts. Moreover, it is suspected to produce close binary companions which will later strongly influence the star's future evolution in the Hertzsprung-Russel diagram. We present three-dimensional gravitation-radiation-hydrodynamic numerical simulations of 100 Mo pre-stellar cores. We find that accretion discs of young massive stars violently fragment without preventing the (highly variable) accretion of gaseous clumps onto the protostars. While acquiring the characteristics of a nascent low-mass companion, some disc fragments migrate onto the central massive protostar with dynamical properties showing that its final Keplerian orbit is close enough to constitute a close massive proto-binary system, having a young high-mass and a low-mass component. We conclude on the viability of the disc fragmentation channel for the formation of such short-period binaries, and that both processes -close massive binary formation and accretion bursts- may happen at the same time. FU-Orionis-type bursts, such as observed in the young high-mass star S255IR-NIRS3, may not only indicate ongoing disc fragmentation, but also be considered as a tracer for the formation of close massive binaries - progenitors of the subsequent massive spectroscopic binaries - once the high-mass component of the system will enter the main-sequence phase of its evolution. Finally, we investigate the ALMA-observability of the disc fragments.

This review is based on lectures given at the 45th Saas-Fee Advanced Course 'From Protoplanetary Disks to Planet Formation' held in March 2015 in Les Diablerets, Switzerland. Starting with an overview of the main characterictics of the Solar System and extrasolar planets, we describe the planet formation process in terms of the sequential accretion scenario. First the growth processes of dust particles to planetesimals and subsequently to terrestrial planets or planetary cores are presented. This is followed by the formation process of the giant planets either by core accretion or gravitational instability. Finally, the dynamical evolution of the orbital elements as driven by disk-planet interaction and the overall evolution of multi-object systems is presented.

In weakly ionized discs turbulence can be generated through the vertical shear instability (VSI). Embedded planets feel a stochastic component in the torques acting on them which can impact their migration. In this work we study the interplay between a growing planet embedded in a protoplanetary disc and the VSI-turbulence. We performed a series of three-dimensional hydrodynamical simulations for locally isothermal discs with embedded planets in the mass range from 5 to 100 Earth masses. We study planets embedded in an inviscid disc that is VSI unstable, becomes turbulent and generates angular momentum transport with an effective $\alpha = 5 \cdot 10^{-4}$. This is compared to the corresponding viscous disc using exactly this $\alpha$-value. In general we find that the planets have only a weak impact on the disc turbulence. Only for the largest planet ($100 M_\oplus$) the turbulent activity becomes enhanced inside of the planet. The depth and width of a gap created by the more massive planets ($30, 100 M_\oplus$) in the turbulent disc equal exactly that of the corresponding viscous case, leading to very similar torque strengths acting on the planet, with small stochastic fluctuations for the VSI disc. At the gap edges vortices are generated that are stronger and longer lived in the VSI disc. Low mass planets (with $M_p \leq 10 M_\oplus$) do not open gaps in the disc in both cases but generate for the turbulent disc an over-density behind the planet that exerts a significant negative torque. This can boost the inward migration in VSI turbulent discs well above the Type I rate. Due to the finite turbulence level in realistic three-dimensional discs the gap depth will always be limited and migration will not stall in inviscid discs.

Context. Mass transfer due to Roche lobe overflow leads to the formation of an accretion disk around a weakly magnetized white dwarf (WD) in cataclysmic variables. At the inner edge of the disk, the gas comes upon the surface of the WD and has to get rid of its excess kinetic energy in order to settle down on the more slowly rotating outer stellar layers. This region is known as the boundary layer (BL). Aims. In this work we investigate the vertical structure of the BL, which is still poorly understood. We shall provide details of the basic structure of the two-dimensional (2D) BL and how it depends on parameters such as stellar mass and rotation rate, as well as the mass-accretion rate. We further investigate the destination of the disk material and compare our results with previous one-dimensional (1D) simulations. Methods. We solve the 2D equations of radiation hydrodynamics in a spherical geometry using a parallel grid-based code that employs a Riemann solver. The radiation energy is considered in the two-temperature approach with a radiative flux given by the flux-limited diffusion approximation. Results. The BL around a non-rotating WD is characterized by a steep drop in angular velocity over a width of only 1% of the stellar radius, a heavy depletion of mass, and a high temperature (500 000 K) as a consequence of the strong shear. Variations in Om_star, M_star, and M_dot influence the extent of the changes of the variables in the BL but not the general structure. Depending on Om_star, the disk material travels up to the poles or is halted at a certain latitude. The extent of mixing with the stellar material also depends on Om_star. We find that the 1D approximation matches the 2D data well, apart from an underestimated temperature.

We present recent improvements in the simulation of regolith sampling processes in microgravity using the numerical particle method smooth particle hydrodynamics (SPH). We use an elastic-plastic soil constitutive model for large deformation and failure flows for dynamical behaviour of regolith. In the context of projected small body (asteroid or small moons) sample return missions, we investigate the efficiency and feasibility of a particular material sampling method: Brushes sweep material from the asteroid's surface into a collecting tray. We analyze the influence of different material parameters of regolith such as cohesion and angle of internal friction on the sampling rate. Furthermore, we study the sampling process in two environments by varying the surface gravity (Earth's and Phobos') and we apply different rotation rates for the brushes. We find good agreement of our sampling simulations on Earth with experiments and provide estimations for the influence of the material properties on the collecting rate.

Recently, the vertical shear instability (VSI) has become an attractive purely hydrodynamic candidate for the anomalous angular momentum transport required for weakly ionized accretion disks. In direct three-dimensional numerical simulations of VSI turbulence in disks, a meridional circulation pattern was observed that is opposite to the usual viscous flow behavior. Here, we investigate whether this feature can possibly be explained by an anisotropy of the VSI turbulence. Using three-dimensional hydrodynamical simulations, we calculate the turbulent Reynolds stresses relevant for angular momentum transport for a representative section of a disk. We find that the vertical stress is significantly stronger than the radial stress. Using our results in viscous disk simulations with different viscosity coefficients for the radial and vertical direction, we find good agreement with the VSI turbulence for the stresses and meridional flow; this provides additional evidence for the anisotropy. The results are important with respect to the transport of small embedded particles in disks.

Most studies concerning the growth and evolution of massive planets focus either on their accretion or their migration only. In this work we study both processes concurrently to investigate how they might mutually affect each other. We modeled a 2-dimensional disk with a steady accretion flow onto the central star and embed a Jupiter mass planet at 5.2 au. The disk is locally isothermal and viscosity is modeled using a constant $\alpha$. The planet is held on a fixed orbit for a few hundred orbits to allow the disk to adapt and carve a gap. After this period, the planet is released and free to move according to the gravitational interaction with the gas disk. The mass accretion onto the planet is modeled by removing a fraction of gas from the inner Hill sphere, and the removed mass and momentum can be added to the planet. Our results show that a fast migrating planet is able to accrete more gas than a slower migrating planet. Utilizing a tracer fluid we analyzed the origin of the accreted gas which comes predominantly originating from the inner disk for a fast migrating planet. In case of slower migration the fraction of gas from the outer disk increases. We also found that even for very high accretion rates in some cases gas crosses the planetary gap from the inner to the outer disk. Our simulations show that the crossing of gas changes during the migration process as the migration rate slows down. Therefore classical type II migration where the planet migrates with the viscous drift rate and no gas crosses the gap is no general process but may only occur for special parameters and at a certain time during the orbital evolution of the planet.

A significant fraction of OB-type, main-sequence massive stars are classified as runaway and move supersonically through the interstellar medium (ISM). Their strong stellar winds interact with their surroundings where the typical strength of the local ISM magnetic field is about 3.5-7 micro-G, which can result in the formation of bow shock nebulae. We investigate the effects of such magnetic fields, aligned with the motion of the flow, on the formation and emission properties of these circumstellar structures. Our axisymmetric, magneto-hydrodynamical simulations with optically-thin radiative cooling, heating and anisotropic thermal conduction show that the presence of the background ISM magnetic field affects the projected optical emission our bow shocks at Ha and [OIII] lambda 5007 which become fainter by about 1-2 orders of magnitude, respectively. Radiative transfer calculations against dust opacity indicate that the magnetic field slightly diminishes their projected infrared emission and that our bow shocks emit brightly at 60 micron. This may explain why the bow shocks generated by ionizing runaway massive stars are often difficult to identify. Finally, we discuss our results in the context of the bow shock of Zeta Ophiuchi and we support the interpretation of its imperfect morphology as an evidence of the presence of an ISM magnetic field not aligned with the motion of its driving star.

Accretion-driven luminosity outbursts are a vivid manifestation of variable mass accretion onto protostars. They are known as the so-called FU Orionis phenomenon in the context of low-mass protostars. More recently, this process has been found in models of primordial star formation. Using numerical radiation hydrodynamics simulations, we stress that present-day forming massive stars also experience variable accretion and show that this process is accompanied by luminous outbursts induced by the episodic accretion of gaseous clumps falling from the circumstellar disk onto the protostar. Consequently, the process of accretion-induced luminous flares is also conceivable in the high-mass regime of star formation and we propose to regard this phenomenon as a general mechanism that can affect protostars regardless of their mass and/or the chemical properties of the parent environment in which they form. In addition to the commonness of accretion-driven outbursts in the star formation machinery, we conjecture that luminous flares from regions hosting forming high-mass star may be an observational implication of the fragmentation of their accretion disks.

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

Among the candidates for generating turbulence in accretion discs in situations with low intrinsic ionization the vertical shear instability (VSI) has become an interesting candidate, as it relies purely on a vertical gradient in the angular velocity. Existing simulations have shown that $\alpha$-values a few times $10^{-4}$ can be generated. The particle growth in the early planet formation phase is determined by the dynamics of dust particles. Here, we address in particular the efficiency of VSI-turbulence in concentrating particles in order to generate overdensities and low collision velocities. We perform 3D numerical hydrodynamical simulations of accretion discs around young stars that include radiative transport and irradiation from the central star. The motion of particles within a size range of a fraction of mm up to several m is followed using standard drag formula. We confirm that under realistic conditions the VSI is able to generate turbulence in full 3D protoplanetary discs. The irradiated disc shows turbulence within 10 to 60au. The mean radial motion of the gas is such that it is directed inward near the midplane and outward in the surface layers. We find that large particles drift inward with the expected speed, while small particles can experience phases of outward drift. Additionally, the particles show bunching behaviour with overdensities reaching 5 times the average value, which is strongest for dimensionless stopping times around unity. Particles in a VSI-turbulent discs are concentrated in large scale turbulent eddies and show low relative speeds that allow for growing collisions. The reached overdensities will also allow for the onset streaming instabilities further enhancing particle growth. The outward drift for small particles at higher disk elevations allows for the transport of processed high temperature material in the Solar System to larger distances.

Modern graphics processing units (GPUs) lead to a major increase in the performance of the computation of astrophysical simulations. Owing to the different nature of GPU architecture compared to traditional central processing units (CPUs) such as x86 architecture, existing numerical codes cannot be easily migrated to run on GPU. Here, we present a new implementation of the numerical method smooth particle hydrodynamics (SPH) using CUDA and the first astrophysical application of the new code: the collision between Ceres-sized objects. The new code allows for a tremendous increase in speed of astrophysical simulations with SPH and self-gravity at low costs for new hardware. We have implemented the SPH equations to model gas, liquids and elastic, and plastic solid bodies and added a fragmentation model for brittle materials. Self-gravity may be optionally included in the simulations and is treated by the use of a Barnes-Hut tree. We find an impressive performance gain using NVIDIA consumer devices compared to our existing OpenMP code. The new code is freely available to the community upon request.

Massive stars that have been ejected from their parent cluster and supersonically sailing away through the interstellar medium (ISM) are classified as exiled. They generate circumstellar bow shock nebulae that can be observed. We present two-dimensional, axisymmetric hydrodynamical simulations of a representative sample of stellar wind bow shocks from Galactic OB stars in an ambient medium of densities ranging from n_ISM=0.01 up to 10.0/cm3. Independently of their location in the Galaxy, we confirm that the infrared is the most appropriated waveband to search for bow shocks from massive stars. Their spectral energy distribution is the convenient tool to analyze them since their emission does not depend on the temporary effects which could affect unstable, thin-shelled bow shocks. Our numerical models of Galactic bow shocks generated by high-mass (~40 Mo) runaway stars yield H$\alpha$ fluxes which could be observed by facilities such as the SuperCOSMOS H-Alpha Survey. The brightest bow shock nebulae are produced in the denser regions of the ISM. We predict that bow shocks in the field observed at Ha by means of Rayleigh-sensitive facilities are formed around stars of initial mass larger than about 20 Mo. Our models of bow shocks from OB stars have the emission maximum in the wavelength range 3 <= lambda <= 50 micrometer which can be up to several orders of magnitude brighter than the runaway stars themselves, particularly for stars of initial mass larger than 20 Mo.

The supersonic motion of gravitating objects through a gaseous medium constitutes a classical problem in theoretical astrophysics. Its application covers a broad range of objects and scales from planets up to galaxies. Especially the dynamical friction, caused by the forming wake behind the object, plays an important role for the dynamics of the system. To calculate the dynamical friction, standard formulae, based on linear theory are often used. It is our goal to check the general validity of these formulae and provide suitable expressions for the dynamical friction acting on the moving object, based on the basic physical parameters of the problem. We perform sequences of high resolution numerical studies of rigid bodies moving supersonically through a homogeneous medium, and calculate the total drag acting on the object, which is the sum of gravitational and hydro drag. We study cases without gravity with purely hydrodynamical drag, as well as gravitating objects. From the final equilibrium state of the simulations, we compute for gravitating objects the dynamical friction by direct numerical integration of the gravitational pull acting on the embedded object. The numerical experiments confirm the known scaling laws for the dependence of the dynamical friction on the basic physical parameters as derived in earlier semi-analytical studies. As a new important result we find that the shock's stand-off distance is revealed as the minimum spatial interaction scale of dynamical friction. Below this radius, the gas settles into a hydrostatic state, which causes no net gravitational pull onto the moving body. Finally, we derive an analytic estimate for the stand-off distance that can be used to calculate the dynamical friction force.

We are observing, thanks to ALMA, the dust distribution in the region of active planet formation around young stars. This is a powerful tool to connect observations with theoretical models and improve our understandings of the processes at play. We want to test how a multi-planetary system shapes its birth disk and study the influence of the planetary masses and particle sizes on the final dust distribution. Moreover, we apply our model to the HL Tau system in order to obtain some insights on the physical parameters of the planets that are able to create the observed features. We follow the evolution of a population of dust particles, treated as Lagrangian particles, in two-dimensional, locally isothermal disks where two equal mass planets are present. The planets are kept in fixed orbits and they do not accrete mass. The outer planet plays a major role removing the dust particles in the co-orbital region of the inner planet and forming a particle ring which promotes the development of vortices respect to the single planetary case. The ring and gaps width depends strongly on the planetary mass and particle stopping times, and for the more massive cases the ring clumps in few stable points that are able to collect a high mass fraction. The features observed in the HL Tau system can be explained through the presence of several massive cores that shape the dust disk, where the inner planet(s) should have a mass on the order of 0.07 Jupiter masses and the outer one(s) on the order of 0.35 Jupiter masses. These values can be significantly lower if the disk mass turns out to be less than previously estimated. Decreasing the disk mass by a factor 10 we obtain similar gap widths for planets with a mass of 10 and 20 Earth masses respectively. Although the particle gaps are prominent, the expected gaseous gaps would be barely visible.

The existence of planets orbiting a central binary star system immediately raises questions regarding their formation and dynamical evolution. Recent discoveries of circumbinary planets by the Kepler space telescope has shown that some of these planets reside close to the dynamical stability limit where it is very difficult to form planets in situ. For binary systems with nearly circular orbits, such as Kepler-38, the observed proximity of planetary orbits to the stability limit can be understood by an evolutionary process in which planets form farther out in the disk and migrate inward to their observed position. The Kepler-34 system has a high orbital eccentricity of 0.52. Here, we analyse evolutionary scenarios for the planet observed around this system using two-dimensional hydrodynamical simulations. The highly eccentric binary opens a wide inner hole in the disk which is also eccentric, and displays a slow prograde precession. As a result of the large, eccentric inner gap, an embedded planet settles in a final equilibrium position that lies beyond the observed location of Kepler-34 b, but has the correct eccentricity. In this configuration the planetary orbit is aligned with the disk in a state of apsidal corotation.To account for the closer orbit of Kepler-34 b to the central binary, we considered a two-planet scenario and examined the evolution of the system through joint inward migration and capture into mean-motion resonances. When the inner planet orbits inside the gap of the disk, planet-planet scattering ensues. While often one object is thrown into a large, highly eccentric orbit, at times the system is left with a planet close to the observed orbit, suggesting that Kepler 34 might have had two circumbinary planets where one might have been scattered out of the system or into an orbit where it did not transit the central binary during the operation of Kepler.

Context. Disk accretion onto weakly magnetized stars leads to the formation of a boundary layer (BL) where the gas loses its excess kinetic energy and settles onto the star. There are still many open questions concerning the BL, for instance the transport of angular momentum (AM) or the vertical structure. Aims. It is the aim of this work to investigate the AM transport in the BL where the magneto-rotational instability (MRI) is not operating owing to the increasing angular velocity $\Omega(r)$ with radius. We will therefore search for an appropriate mechanism and examine its efficiency and implications. Methods. We perform 2D numerical hydrodynamical simulations in a cylindrical coordinate system $(r, \varphi)$ for a thin, vertically inte- grated accretion disk around a young star. We employ a realistic equation of state and include both cooling from the disk surfaces and radiation transport in radial and azimuthal direction. The viscosity in the disk is treated by the {\alpha}-model; in the BL there is no viscosity term included. Results. We find that our setup is unstable to the sonic instability which sets in shortly after the simulations have been started. Acoustic waves are generated and traverse the domain, developing weak shocks in the vicinity of the BL. Furthermore, the system undergoes recurrent outbursts where the activity in the disk increases strongly. The instability and the waves do not die out for over 2000 orbits. Conclusions. There is indeed a purely hydrodynamical mechanism that enables AM transport in the BL. It is efficient and wave mediated; however, this renders it a non-local transport method, which means that models of a effective local viscosity like the {\alpha}-viscosity are probably not applicable in the BL. A variety of further implications of the non-local AM transport are discussed.

Massive planets that open a gap in the accretion disk are believed to migrate with exactly the viscous speed of the disk, a regime termed type II migration. Population synthesis models indicate that standard type II migration is too rapid to be in agreement with the observations. We study the migration of massive planets between $2\times10^{-4}$ and $2\times10^{-3} M_\odot$ corresponding to 0.2 to 2 Jupiter masses $M_J$. in order to estimate the migration rate in comparison to type II migration. We follow the evolution of planets embedded in two-dimensional, locally isothermal disks with non-zero mass accretion which is explicitly modelled using suitable in- and outflow boundary conditions to ensure a specific accretion rate. After a certain relaxation time we release the planet and measure its migration through the disk and the dependence on parameters such as viscosity, accretion rate and planet mass. We study accreting and non-accretion planets. The inferred migration rate of the planet is determined entirely by the disk torques acting on it and is completely independent of the viscous inflow velocity, so there is no classical type II migration regime. Depending on the local disk mass the migration rate can be faster or slower than type II migration. From the torques and the accretion rate profile in the disk we see that the gap formed by the planet does not separate the inner from the outer disk as necessary for type II migration, rather gas crosses the gap or is accreted onto the planet.

Context: Anticyclonic vortices are considered as a favourable places for trapping dust and forming planetary embryos. On the other hand, they are massive blobs that can interact gravitationally with the planets in the disc. Aims: We aim to study how a vortex interacts gravitationally with a planet which migrates toward it or a planet which is created inside the vortex. Methods: We performed hydrodynamical simulations of a viscous locally isothermal disc using GFARGO and FARGO-ADSG. We set a stationary Gaussian pressure bump in the disc in a way that RWI is triggered. After a large vortex is established, we implanted a low mass planet in the outer disc or inside the vortex and allowed it to migrate. We also examined the effect of vortex strength on the planet migration and checked the validity of the final result in the presence of self-gravity. Results: We noticed regardless of the planet's initial position, the planet is finally locked to the vortex or its migration is stopped in a farther orbital distance in case of a stronger vortex. For the model with the weaker vortex, we studied the effect of different parameters such as background viscosity, background surface density, mass of the planet and different planet positions. In these models, while the trapping time and locking angle of the planet vary for different parameters, the main result, which is the planet-vortex locking, remains valid. We discovered that even a planet with a mass less than 5 * 10^{-7} M_{\star} comes out from the vortex and is locked to it at the same orbital distance. For a stronger vortex, both in non-self-gravitated and self-gravitating models, the planet migration is stopped far away from the radial position of the vortex. This effect can make the vortices a suitable place for continual planet formation under the condition that they save their shape during the planetary growth.

The origin of turbulence in accretion discs is still not fully understood. While the magneto-rotational instability is considered to operate in sufficiently ionized discs, its role in the poorly ionized protoplanetary disc is questionable. Recently, the vertical shear instability (VSI) has been suggested as a possible alternative. Our goal is to study the characteristics of this instability and the efficiency of angular momentum transport, in extended discs, under the influence of radiative transport and irradiation from the central star. We use multi-dimensional hydrodynamic simulations to model a larger section of an accretion disc. First we study inviscid and weakly viscous discs using a fixed radial temperature profile in two and three spatial dimensions. The simulations are then extended to include radiative transport and irradiation from the central star. In agreement with previous studies we find for the isothermal disc a sustained unstable state with a weak positive angular momentum transport of the order of $\alpha \approx 10^{-4}$. Under the inclusion of radiative transport the disc cools off and the turbulence terminates. For discs irradiated from the central star we find again a persistent instability with a similar $\alpha$ value as for the isothermal case. We find that the VSI can indeed generate sustained turbulence in discs albeit at a relatively low level with $\alpha$ about few times $10^{-4}$

Here we compute detailed model spectra of recently published optically thick one-dimensional radial baundary layer (BL) models in cataclysmic variables and compare them with observed soft X-ray/extreme ultraviolet (EUV) spectra of dwarf novae in outburst. Every considered BL model is divided into a number of rings, and for each ring, a structure model along the vertical direction is computed using the stellar-atmosphere method. The ring spectra are then combined into a BL spectrum taking Doppler broadening and limb darkening into account. Two sets of model BL spectra are computed, the first of them consists of BL models with fixed white dwarf (WD) mass (1 M_sun) and various relative WD angular velocities (0.2, 0.4, 0.6 and 0.8 break-up velocities), while the other deals with a fixed relative angular velocity (0.8 break-up velocity) and various WD masses (0.8, 1, and 1.2 M_sun). The model spectra show broad absorption features because of blending of numerous absorption lines, and emission-like features at spectral regions with only a few strong absorption lines. The model spectra are very similar to observed soft X-ray/EUV spectra of SS Cyg and U Gem in outburst. The observed SS Cyg spectrum could be fitted by BL model spectra with WD masses 0.8 - 1 M_sun and relative angular velocities 0.6 - 0.8 break up velocities. These BL models also reproduce the observed ratio of BL luminosity and disk luminosity. The difference between the observed and the BL model spectra is similar to a hot optically thin plasma spectrum and could be associated with the spectrum of outflowing plasma with a mass loss rate compatible with the BL mass accretion rate. The suggested method of computing BL spectra seems very promising and can be applied to other BL models for comparison with EUV spectra of dwarf novae in outburst.

PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2 with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars in wide fields to detect and characterize planets down to Earth-size by photometric transits, whose masses can then be determined by ground-based radial-velocity follow-up measurements. Asteroseismology will be performed for these bright stars to obtain highly accurate stellar parameters, including masses and ages. The combination of bright targets and asteroseismology results in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet radii, masses and ages, respectively. The planned baseline observing strategy includes two long pointings (2-3 years) to detect and bulk characterize planets reaching into the habitable zone (HZ) of solar-like stars and an additional step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will observe up to 1,000,000 stars and detect and characterize hundreds of small planets, and thousands of planets in the Neptune to gas giant regime out to the HZ. It will therefore provide the first large-scale catalogue of bulk characterized planets with accurate radii, masses, mean densities and ages. This catalogue will include terrestrial planets at intermediate orbital distances, where surface temperatures are moderate. Coverage of this parameter range with statistical numbers of bulk characterized planets is unique to PLATO 2.0.

Recently, a number of planets orbiting binary stars have been discovered by the Kepler space telescope. In a few systems the planets reside close to the dynamical stability limit. Due to the difficulty of forming planets in such close orbits, it is believed that they have formed further out in the disk and migrated to their present locations. Our goal is to construct more realistic models of planet migration in circumbinary disks, and to determine the final position of these planets. In our work, we focus on the system Kepler-38. The evolution of the circumbinary disk is studied using two-dimensional hydrodynamical simulations. We study locally isothermal disks as well as more realistic models with viscous heating, radiative cooling from the disk surfaces, and radiative diffusion in the disk mid plane. After the disk has been brought into equilibrium, a 115 Earth-mass planet is embedded and its evolution is followed. In all cases the planets stop inward migration near the inner edge of the disk. In isothermal disks with a typical disk scale height of H/r = 0.05, the final outcome agrees very well with the observed location of planet Kepler-38b. For the radiative models, the disk thickness and location of the inner edge is determined by the mass in the system. For surface densities in the order of 3000 g/cm^2 at 1 AU, the inner gap lies close to the binary and planets stop in the region between the 5:1 and 4:1 mean-motion resonances with the binary. A model with a disk with approximately a quarter of the mass yields a final position very close to the observed one. For planets migrating in circumbinary disks, the final position is dictated by the structure of the disk. Knowing the observed orbits of circumbinary planets, radiative disk simulations with embedded planets can provide important information on the physical state of the system during the final stages of its evolution.

The great diversity of extrasolar planetary systems has challenged our understanding of how planets form, and how their orbits evolve as they form. Among the various processes that may account for this diversity, the gravitational interaction between planets and their parent protoplanetary disc plays a prominent role in shaping young planetary systems. Planet-disc forces are large, and the characteristic times for the evolution of planets orbital elements are much shorter than the lifetime of protoplanetary discs. The determination of such forces is challenging, because it involves many physical mechanisms and it requires a detailed knowledge of the disc structure. Yet, the intense research of the past few years, with the exploration of many new avenues, represents a very significant improvement on the state of the discipline. This chapter reviews current understanding of planet-disc interactions, and highlights their role in setting the properties and architecture of observed planetary systems.

Disk accretion onto stars leads to the formation of a Boundary Layer (BL) near the stellar surface where the disk makes contact with the star. Albeit a large fraction of the total luminosity of the system originates from this tiny layer connecting the accretion disk and the accreting object, its structure has not been fully understood yet. It is the aim of this work, to obtain more insight into the Boundary Layer around the white dwarf in compact binary systems. There are still many uncertainties concerning the extent and temperature of the BL and the rotation rate of the white dwarf. We perform numerical hydrodynamical simulations, where the problem is treated in a one-dimensional, radial approximation (slim disk). The turbulence is described by the alpha parameter viscosity. We include both cooling from the disk surfaces and radial radiation transport. The radiation energy is treated in the one-temperature approximation. For a given M_dot our results show a strong dependence on the stellar mass and rotation rate. Both, the midplane and the effective temperature rise considerably with increasing stellar mass or decreasing stellar rotation rate. Our simulations further show, that the radiation energy and pressure are indeed important in the BL. However, some models show a low optical depth in the BL, making it necessary to find a better representation for optically thin regions. The combination of a high mass and a small radius, such as in white dwarfs, can lead to an enormous energy release in the BL, provided the WD rotates slowly. Since the radial extent of BLs is typically very small (about 0.02 to 0.05 R_star), this leads to surface temperatures of a few hundred thousand Kelvin. All of our models showed subsonic infall velocities with Mach numbers of < 0.4 at most.

Transitional disks are protoplanetary disk around young stars that display inner holes in the dust distribution within a few AU, which is accompanied nevertheless by some gas accretion onto the central star. These cavities could possibly be created by the presence of one or more massive planets. If the gap is created by planets and gas is still present in it, then there should be a flow of gas past the planet into the inner region. It is our goal to study the mass accretion rate into the gap and in particular the dependency on the planet's mass and the thermodynamic properties of the disk. We performed 2D hydro simulations for disks with embedded planets. We added radiative cooling from the disk surfaces, radiative diffusion in the disk midplane, and stellar irradiation to the energy equation to have more realistic models. The mass flow rate into the gap region depends, for given disk thermodynamics, non-monotonically on the mass of the planet. Generally, more massive planets open wider and deeper gaps which would tend to reduce the mass accretion into the inner cavity. However, for larger mass planets the outer disk becomes eccentric and the mass flow rate is enhanced over the low mass cases. As a result, for the isothermal disks the mass flow is always comparable to the expected mass flow of unperturbed disks M_d, while for more realistic radiative disks the mass flow is very small for low mass planets (<= 4 M_jup) and about 50% for larger planet masses. For the radiative disks that critical planet mass for the disk to become eccentric is much larger that in the isothermal case. Massive embedded planets can reduce the mass flow across the gap considerably, to values of about an order of magnitude smaller than the standard disk accretion rate, and can be responsible for opening large cavities. The remaining mass flow into the central cavity is in good agreement with the observations.

Aims: The aim of this work is to gain a deeper insight into how much different aggregate types are affected by erosion. Especially, it is important to study the influence of the velocity of the impacting projectiles. We also want to provide models for dust growth in protoplanetary disks with simple recipes to account for erosion effects. Methods: To study the erosion of dust aggregates we employed a molecular dynamics approach that features a detailed micro-physical model of the interaction of spherical grains. For the first time, the model has been extended by introducing a new visco-elastic damping force which requires a proper calibration. Afterwards, different sample generation methods were used to cover a wide range of aggregate types. Results: The visco-elastic damping force introduced in this work turns out to be crucial to reproduce results obtained from laboratory experiments. After proper calibration, we find that erosion occurs for impact velocities of 5 m/s and above. Though fractal aggregates as formed during the first growth phase are most susceptible to erosion, we observe erosion of aggregates with rather compact surfaces as well. Conclusions: We find that bombarding a larger target aggregate with small projectiles results in erosion for impact velocities as low as a few m/s. More compact aggregates suffer less from erosion. With increasing projectile size the transition from accretion to erosion is shifted to higher velocities. This allows larger bodies to grow through high velocity collisions with smaller aggregates.

The transport of energy through radiation is very important in many astrophysical phenomena. In dynamical problems the time-dependent equations of radiation hydrodynamics have to be solved. We present a newly developed radiation-hydrodynamics module specifically designed for the versatile MHD code PLUTO. The solver is based on the flux-limited diffusion approximation in the two-temperature approach. All equations are solved in the co-moving frame in the frequency independent (grey) approximation. The hydrodynamics is solved by the different Godunov schemes implemented in PLUTO, and for the radiation transport we use a fully implicit scheme. The resulting system of linear equations is solved either using the successive over-relaxation (SOR) method (for testing purposes), or matrix solvers that are available in the PETSc library. We state in detail the methodology and describe several test cases in order to verify the correctness of our implementation. The solver works in standard coordinate systems, such as Cartesian, cylindrical and spherical, and also for non-equidistant grids. We have presented a new radiation-hydrodynamics solver coupled to the MHD-code \PLUTO that is a modern, versatile and efficient new module for treating complex radiation hydrodynamical problems in astrophysics. As test cases, either purely radiative situations, or full radiation-hydrodynamical setups (including radiative shocks and convection in accretion discs) have been studied successfully. The new module scales very well on parallel computers using MPI. For problems in star or planet formation, we have added the possibility of irradiation by a central source.

Context.Within the sequential accretion scenario of planet formation, planets are build up through a sequence sticking collisions. The outcome of collisions between porous dust aggregates is very important for the growth from very small dust particles to planetesimals. In this work we determine the necessary material properties of dust aggregates as a function the porosity. Aims: Continuum models such as SPH that are capable of simulating collisions of macroscopic dust aggregates require a set of material parameters. Some of them such as the tensile and shear strength are difficult to obtain from laboratory experiments. The aim of this work is to determine these parameters from ab-initio molecular dynamics simulations. Methods: We simulate the behavior of porous dust aggregates using a detailed micro-physical model of the interaction of spherical grains that includes adhesion forces, rolling, twisting, and sliding. Using different methods of preparing the samples we study the strength behavior of our samples with varying porosity and coordination number of the material. Results: For the tensile strength, we can reproduce data from laboratory experiments very well. For the shear strength, there are no experimental data available. The results from our simulations differ significantly from previous theoretical models, which indicates that the latter might not be sufficient to describe porous dust aggregates. Conclusions: We have provided functional behavior of tensile and shear strength of porous dust aggregates as a function of the porosity that can be directly applied in continuum simulations of these objects in planet formation scenarios.

We carry out three-dimensional Smoothed Particle Hydrodynamics simulations of spherical homogeneous SiO2 dust aggregates to investigate how the mass and the porosity of the aggregates affects their ability to survive an impact at various different collision velocities (between 1 - 27.5m/s). We explore how the threshold velocities for fragmentation vary with these parameters. Crucially, we find that the porosity plays a part of utmost importance in determining the outcome of collisions. In particular, we find that aggregates with filling factors >37% are significantly weakened and that the velocity regime in which the aggregates grow is reduced or even non-existent (instead, the aggregates either rebound off each other or break apart). At filling factors less than ~37% we find that more porous objects are weaker but not as weak as highly compact objects with filling factors >37%. In addition we find that (for a given aggregate density) collisions between very different mass objects have higher threshold velocities than those between very similar mass objects. We find that fragmentation velocities are higher than the typical values of 1m/s and that growth can even occur for velocities as high as 27.5m/s. Therefore, while the growth of aggregates is more likely if collisions between different sized objects occurs or if the aggregates are porous with filling factor <37%, it may also be hindered if the aggregates become too compact.

The Kepler-36 system consists of two planets that are spaced unusually close together, near the 7:6 mean motion resonance. While it is known that mean motion resonances can easily form by convergent migration, Kepler-36 is an extreme case due to the close spacing and the relatively high planet masses of 4 and 8 times that of the Earth. In this paper, we investigate whether such a system can be obtained by interactions with the protoplanetary disc. These discs are thought to be turbulent and exhibit density fluctuations which might originate from the magneto-rotational instability. We adopt a realistic description for stochastic forces due to these density fluctuations and perform both long term hydrodynamical and N-body simulations. Our results show that planets in the Kepler-36 mass range can be naturally assembled into a closely spaced planetary system for a wide range of migration parameters in a turbulent disc similar to the minimum mass solar nebula. The final orbits of our formation scenarios tend to be Lagrange stable, even though large parts of the parameter space are chaotic and unstable.

This study aims to determine the impact of stellar binary companions on the lifetime and evolution of circumstellar disks in the Chamaeleon I (Cha I) star-forming region by measuring the frequency and strength of accretion and circumstellar dust signatures around the individual components of T Tauri binary stars. We used high-angular resolution adaptive optics JHKL'-band photometry and 1.5-2.5mu spectroscopy of 19 visual binary and 7 triple stars in Cha I - including one newly discovered tertiary component - with separations between ~25 and ~1000au. The data allowed us to infer stellar component masses and ages and, from the detection of near-infrared excess emission and the strength of Brackett-gamma emission, the presence of ongoing accretion and hot circumstellar dust of the individual stellar component of each binary. Of all the stellar components in close binaries with separations of 25-100au, 10(+15-5)% show signs of accretion. This is less than half of the accretor fraction found in wider binaries, which itself appears significantly reduced (~44%) compared with previous measurements of single stars in Cha I. Hot dust was found around 50(+30-15)% of the target components, a value that is indistinguishable from that of Cha I single stars. Only the closest binaries (<25au) were inferred to have a significantly reduced fraction (<~25%) of components that harbor hot dust. Accretors were exclusively found in binary systems with unequal component masses M_secondary/M_primary < 0.8, implying that the detected accelerated disk dispersal is a function of mass-ratio. This agrees with the finding that only one accreting secondary star was found, which is also the weakest accretor in the sample. The results imply that disk dispersal is more accelerated the stronger the dynamical disk truncation, i.e., the smaller the inferred radius of the disk. (abridged)

The edges of magnetically-dead zones in protostellar disks have been proposed as locations where density bumps may arise, trapping planetesimals and helping form planets. Magneto-rotational turbulence in magnetically-active zones provides both accretion of gas on the star and transport of mass to the dead zone. We investigate the location of the magnetically-active regions in a protostellar disk around a solar-type star, varying the disk temperature, surface density profile, and dust-to-gas ratio. We also consider stellar masses between 0.4 and 2 $M_\odot$, with corresponding adjustments in the disk mass and temperature. The dead zone's size and shape are found using the Elsasser number criterion with conductivities including the contributions from ions, electrons, and charged fractal dust aggregates. The charged species' abundances are found using the approach proposed by S. Okuzumi. The dead zone is in most cases defined by the ambipolar diffusion. In our maps, the dead zone takes a variety of shapes, including a fish-tail pointing away from the star and islands located on and off the midplane. The corresponding accretion rates vary with radius, indicating locations where the surface density will increase over time, and others where it will decrease. We show that density bumps do not readily grow near the dead zone's outer edge, independently of the disk parameters and the dust properties. Instead, the accretion rate peaks at the radius where the gas-phase metals freeze out. This could lead to clearing a valley in the surface density, and to a trap for pebbles located just outside the metal freeze-out line.

Context: Bouncing collisions of dust aggregates within the protoplanetary may have a significant impact on the growth process of planetesimals. Yet, the conditions that result in bouncing are not very well understood. Existing simulations studying the bouncing behavior used aggregates with an artificial, very regular internal structure. Aims: Here, we study the bouncing behavior of sub-mm dust aggregates that are constructed applying different sample preparation methods. We analyze how the internal structure of the aggregate alters the collisional outcome and determine the influence of aggregate size, porosity, collision velocity, and impact parameter. Methods: We use molecular dynamics simulations where the individual aggregates are treated as spheres that are made up of several hundred thousand individual monomers. The simulations are run on GPUs. Results: Statistical bulk properties and thus bouncing behavior of sub-mm dust aggregates depend heavily on the preparation method. In particular, there is no unique relation between the average volume filling factor and the coordination number of the aggregate. Realistic aggregates bounce only if their volume filling factor exceeds 0.5 and collision velocities are below 0.1 m/s. Conclusions: For dust particles in the protoplanetary nebula we suggest that the bouncing barrier may not be such a strong handicap in the growth phase of dust agglomerates, at least in the size range of 0.1 mm.

The strength and direction of migration of low mass embedded planets depends on the disk's thermodynamic state, where the internal dissipation is balanced by radiative transport, and the migration can be directed outwards, a process which extends the lifetime of growing embryos. Very important parameters determining the structure of disks, and hence the direction of migration, are the viscosity and the adiabatic index. In this paper we investigate the influence of different viscosity prescriptions (alpha-type and constant) and adiabatic indices on disk structures and how this affects the migration rate of planets embedded in such disks. We perform 3D numerical simulations of accretion disks with embedded planets. We use the explicit/implicit hydrodynamical code NIRVANA that includes full tensor viscosity and radiation transport in the flux-limited diffusion approximation, as well as a proper equation of state for molecular hydrogen. The migration of embedded 20Earthmass planets is studied. Low-viscosity disks have cooler temperatures and the migration rates of embedded planets tend toward the isothermal limit. In these disks, planets migrate inwards even in the fully radiative case. The effect of outward migration can only be sustained if the viscosity in the disk is large. Overall, the differences between the treatments for the equation of state seem to play a more important role in disks with higher viscosity. A change in the adiabatic index and in the viscosity changes the zero-torque radius that separates inward from outward migration. For larger viscosities, temperatures in the disk become higher and the zero-torque radius moves to larger radii, allowing outward migration of a 20 Earth-mass planet to persist over an extended radial range. In combination with large disk masses, this may allow for an extended period of the outward migration of growing protoplanetary cores.

The strength and direction of migration of low mass planets depends on the disc's thermodynamics. In discs where the viscous heating is balanced by radiative transport, the migration can be directed outwards, a process which extends the lifetime of growing planetary embryos. We investigate the influence of opacity and stellar irradiation on the disc thermodynamics. Utilizing the resulting disc structure, we determine the regions of outward migration. We perform two-dimensional numerical simulations of equilibrium discs with viscous heating, radiative cooling and stellar irradiation. We use the hydrodynamical code NIRVANA that includes a full tensor viscosity and stellar irradiation, as well as a two temperature solver that includes radiation transport in the flux-limited diffusion approximation. The migration is studied by using torque formulae. In the constant opacity case, we reproduce the analytical results of a black-body disc: the stellar irradiation dominates in the outer regions -- leading to flaring -- while the viscous heating dominates close to the star. We find that the inner edge of the disc should not be significantly puffed-up by the stellar irradiation. If the opacity depends on the local density and temperature, the structure of the disc is different, and several bumps in the aspect ratio H/r appear, due to transitions between different opacity regimes. The bumps in the disc can shield the outer disc from stellar irradiation. Stellar irradiation is an important factor for determining the disc structure and has dramatic consequences for the migration of embedded planets. Compared to discs with only viscous heating, a stellar irradiated disc features a much smaller region of outward migration for a range of planetary masses. This suggests that the region where the formation of giant planet cores takes place is smaller, which in turn might lead to a shorter growth phase.

We performed micro-gravity collision experiments in our laboratory drop-tower using 5-cm-sized dust agglomerates with volume filling factors of 0.3 and 0.4, respectively. This work is an extension of our previous experiments reported in Beitz et al. (2011) to aggregates of more than one order of magnitude higher masses. The dust aggregates consisted of micrometer-sized silica particles and were macroscopically homogeneous. We measured the coefficient of restitution for collision velocities ranging from 1 cm/s to 0.5 m/s, and determined the fragmentation velocity. For low velocities, the coefficient of restitution decreases with increasing impact velocity, in contrast to findings by Beitz et al. (2011). At higher velocities, the value of the coefficient of restitution becomes constant, before the aggregates break at the onset of fragmentation. We interpret the qualitative change in the coefficient of restitution as the transition from a solid-body-dominated to a granular-medium-dominated behavior. We complement our experiments by molecular dynamics simulations of porous aggregates and obtain a reasonable match to the experimental data. We discuss the importance of our experiments for protoplanetary disks, debris disks, and planetary rings. The work is an extensional study to previous work of our group and gives a new insight in the velocity dependency of the coefficient of restitution due to improved measurements, better statistics and a theoretical approach.

Embedded planets disturb the density structure of the ambient disk and gravitational back-reaction will induce possibly a change in the planet's orbital elements. The accurate determination of the forces acting on the planet requires careful numerical analysis. Recently, the validity of the often used fast orbital advection algorithm (FARGO) has been put into question, and special numerical resolution and stability requirements have been suggested. In this paper we study the process of planet-disk interaction for small mass planets of a few Earth masses, and reanalyze the numerical requirements to obtain converged and stable results. One focus lies on the applicability of the FARGO-algorithm. Additionally, we study the difference of two and three-dimensional simulations, compare global with local setups, as well as isothermal and adiabatic conditions. We study the influence of the planet on the disk through two- and three-dimensional hydrodynamical simulations. To strengthen our conclusions we perform a detailed numerical comparison where several upwind and Riemann-solver based codes are used with and without the FARGO-algorithm. With respect to the wake structure and the torque density acting on the planet we demonstrate that the FARGO-algorithm yields correct results, and that at a fraction of the regular cpu-time. We find that the resolution requirements for achieving convergent results in unshocked regions are rather modest and depend on the pressure scale height of the disk. By comparing the torque densities of 2D and 3D simulations we show that a suitable vertical averaging procedure for the force gives an excellent agreement between the two. We show that isothermal and adiabatic runs can differ considerably, even for adiabatic indices very close to unity.

As planets form and grow within gaseous protoplanetary disks, the mutual gravitational interaction between the disk and planet leads to the exchange of angular momentum, and migration of the planet. We review current understanding of disk-planet interactions, focussing in particular on physical processes that determine the speed and direction of migration. We describe the evolution of low mass planets embedded in protoplanetary disks, and examine the influence of Lindblad and corotation torques as a function of the disk properties. The role of the disk in causing the evolution of eccentricities and inclinations is also discussed. We describe the rapid migration of intermediate mass planets that may occur as a runaway process, and examine the transition to gap formation and slower migration driven by the viscous evolution of the disk for massive planets. The roles and influence of disk self-gravity and magnetohydrodynamic turbulence are discussed in detail, as a function of the planet mass, as is the evolution of multiple planet systems. Finally, we address the question of how well global models of planetary formation that include migration are able to match observations of extrasolar planets.

In 2D simulations of thin gaseous disks with embedded planets or self-gravity the gravitational potential needs to be smoothed to avoid singularities in the numerical evaluation of the gravitational potential or force. In order to correctly resemble the realistic case of vertically extended 3D disks the softening prescription used in 2D needs to be adjusted properly. In this paper we analyze the embedded planet and the self-gravity case and provide a method to evaluate the required smoothing in 2D simulations of thin disks. Starting from the averaged hydrodynamic equations and using a vertically isothermal disk model, we calculate the force. We compare our results to the often used Plummer form of the potential which runs as \propto 1/(r^2+\epsilon^2)^{1/2}. For that purpose we compute the required smoothing length \epsilon as a function of distance r to the planet or to a disk element within a self-gravitating disk. We find that for larger distances \epsilon is determined solely by the vertical disk thickness H. For the planet case we find that outside r\simH a value of \epsilon = 0.7H describes the averaged force very well, while in the self-gravitating disk the value needs to be larger, \epsilon = 1.2H. For smaller distances the smoothing needs to be reduced significantly. Comparing torque densities of 3D and 2D simulations we show that the modification to the vertical density stratification as induced by an embedded planet needs to be taken into account to obtain agreeing results. In disk fragmentation simulations the choice of \epsilon can determine whether a disk will fragment or not. Because a wrong smoothing length can change even the direction of migration, it is very important to include the effect of the planet on the local scale height in 2D planet-disk simulations. We provide an approximate and fast method for this purpose which gives very good agreement to full 3D simulations.

The early planetesimal growth proceeds through a sequence of sticking collisions of dust agglomerates. Very uncertain is still the relative velocity regime in which growth rather than destruction can take place. The outcome of a collision depends on the bulk properties of the porous dust agglomerates. Continuum models of dust agglomerates require a set of material parameters that are often difficult to obtain from laboratory experiments. Here, we aim at determining those parameters from ab-initio molecular dynamics simulations. Our goal is to improveon the existing model that describe the interaction of individual monomers. We use a molecular dynamics approach featuring a detailed micro-physical model of the interaction of spherical grains. The model includes normal forces, rolling, twisting and sliding between the dust grains. We present a new treatment of wall-particle interaction that allows us to perform customized simulations that directly correspond to laboratory experiments. We find that the existing interaction model by Dominik & Tielens leads to a too soft compressive strength behavior for uni and omni-directional compression. Upon making the rolling and sliding coefficients stiffer we find excellent agreement in both cases. Additionally, we find that the compressive strength curve depends on the velocity with which the sample is compressed. The modified interaction strengths between two individual dust grains will lead to a different behaviour of the whole dust agglomerate. This will influences the sticking probabilities and hence the growth of planetesimals. The new parameter set might possibly lead to an enhanced sticking as more energy can be stored in the system before breakup.

In this paper we study the evolution of viscous and radiative circumstellar disks under the influence of a companion star. We focus on the eccentric {\gamma} Cephei and {\alpha} Centauri system as examples and compare the disk quantities such as disk eccentricity and precession rate to previous isothermal simulations. We perform two-dimensional hydrodynamical simulations of the binary star systems under the assumption of coplanarity of the disk, host star and binary companion. We use the grid-based, staggered mesh code FARGO with an additional energy equation to which we added radiative cooling based on opacity tables. The eccentric binary companion perturbs the disk around the primary star periodically. Upon passing periastron spirals arms are induced that wind from the outer disk towards the star. In isothermal simulations this results in disk eccentricities up to {\epsilon}_disk ~ 0.2, but in more realistic radiative models we obtain much smaller eccentricities of about {\epsilon}_disk ~ 0.04 - 0.06 with no real precession. Models with varying viscosity and disk mass indicate show that disks with less mass have lower temperatures and higher disk eccentricity. The rather large high disk eccentricities, as indicated in previous isothermal disk simulations, implied a more difficult planet formation in the {\gamma} Cephei system due to the enhanced collision velocities of planetesimals. We have shown that under more realistic conditions with radiative cooling the disk become less eccentric and thus planet formation may be made easier. However, we estimate that the viscosity in the disk has to very small, with {\alpha} \lesssim 0.001, because otherwise the disk's lifetime will be too short to allow planet formation to occur along the core instability scenario. We estimate that the periodic heating of the disk in eccentric binaries will be observable in the mid-IR regime.

The aim of the present paper is to investigate the spatial structure of a protoplanetary disc whose dynamics is governed by magnetorotational turbulence. We perform a series of local 3D chemo-radiative MHD simulations located at different radii of a disc which is twice as massive as the standard minimum mass solar nebula of Hayashi (1981). The ionisation state of the disc is calculated by including collisional ionisation, stellar X-rays, cosmic rays and the decay of radionuclides as ionisation sources, and by solving a simplified chemical network which includes the effect of the absorption of free charges by {\mu}m-sized dust grains. In the region where the ionisation is too low to assure good coupling between matter and magnetic fields, a non-turbulent central "dead zone" forms, which ranges approximately from a distance of 2 AU to 4 AU from the central star. The approach taken in the present work allows for the first time to derive the global spatial structure of a protoplanetary disc from a set of physically realistic numerical simulations.

In the framework of the coagulation scenario, kilometre-sized planetesimals form by subsequent collisions of pre-planetesimals of sizes from centimetre to hundreds of metres. Pre-planetesimals are fluffy, porous dust aggregates, which are inhomogeneous owing to their collisional history. Planetesimal growth can be prevented by catastrophic disruption in pre-planetesimal collisions above the destruction velocity threshold. We develop an inhomogeneity model based on the density distribution of dust aggregates, which is assumed to be a Gaussian distribution with a well-defined standard deviation. As a second input parameter, we consider the typical size of an inhomogeneous clump. These input parameters are easily accessible by laboratory experiments. For the simulation of the dust aggregates, we utilise a smoothed particle hydrodynamics (SPH) code with extensions for modelling porous solid bodies. The porosity model was previously calibrated for the simulation of silica dust, which commonly serves as an analogue for pre-planetesimal material. The inhomogeneity is imposed as an initial condition on the SPH particle distribution. We carry out collisions of centimetre-sized dust aggregates of intermediate porosity. We vary the standard deviation of the inhomogeneous distribution at fixed typical clump size. The collision outcome is categorised according to the four-population model. We show that inhomogeneous pre-planetesimals are more prone to destruction than homogeneous aggregates. Even slight inhomogeneities can lower the threshold for catastrophic disruption. For a fixed collision velocity, the sizes of the fragments decrease with increasing inhomogeneity. Pre-planetesimals with an active collisional history tend to be weaker. This is a possible obstacle to collisional growth and needs to be taken into account in future studies of the coagulation scenario.

The migration of planets plays an important role in the early planet-formation process. An important problem has been that standard migration theories predict very rapid inward migration, which poses problems for population synthesis models. However, it has been shown recently that low-mass planets (20-30 Earth Masses) that are still embedded in the protoplanetary disc can migrate outwards under certain conditions. Simulations have been performed mostly for planets at given radii for a particular disc model. Here, we plan to extend previous work and consider different masses of the disc to quantify the influence of the physical disc conditions on planetary migration. We perform three-dimensional (3D) radiation hydrodynamical simulations of embedded planets in protoplanteary discs. For planets on circular orbits at various locations we measure the radial dependece of the torques. For all considered planet masses (20-30 Earth masses) in this study we find outward migration within a limited radial range of the disc, typically from about 0.5 up to 1.5-2.5 a_Jup. Inside and outside this intervall, migration is inward and given by the Lindblad value for large radii. Because outward migration stops at a certain location in the disc, there exists a zero-torque distance for planetary embryos, where they can continue to grow without moving too fast. For higher disc masses (M_disc > 0.02 M_Sol) convection ensues, which changes the structure of the disc and therefore the torque on the planet as well. Outward migration stops at different points in the disc for different planetary masses, resulting in a quite extended region where the formation of larger cores might be easier. In higher mass discs, convection changes the disc's structure resulting in fluctuations in the surface density, which influence the torque acting on the planet, and therefore its migration rate.

Within the collision growth scenario for planetesimal formation, the growth step from centimetre sized pre-planetesimals to kilometre sized planetesimals is still unclear. The formation of larger objects from the highly porous pre-planetesimals may be halted by a combination of fragmentation in disruptive collisions and mutual rebound with compaction. However, the right amount of fragmentation is necessary to explain the observed dust features in late T Tauri discs. Therefore, detailed data on the outcome of pre-planetesimal collisions is required and has to be presented in a suitable and precise format. We propose and apply a new classification scheme for pre-planetesimal collisions based on the quantitative aspects of four fragment populations: the largest and second largest fragment, a power-law population, and a sub-resolution population. For the simulations of pre-planetesimal collisions, we adopt the SPH numerical scheme with extensions for the simulation of porous solid bodies. By means of laboratory benchmark experiments, this model was previously calibrated and tested for the correct simulation of the compaction, bouncing, and fragmentation behaviour of macroscopic highly porous silica dust aggregates. It is shown that previous attempts to map collision data were much too oriented on qualitatively categorising into sticking, bouncing, and fragmentation events. We show that the four-populations model encompasses all previous categorisations and in addition allows for transitions. This is because it is based on quantitative characteristic attributes of each population such as the mass, kinetic energy, and filling factor. As a demonstration of the applicability and the power of the four-populations model, we utilise it to present the results of a study on the influence of collision velocity in head-on collisions of intermediate porosity aggregates.

While planets in the solar system only have a low inclination with respect to the ecliptic there is mounting evidence that in extrasolar systems the inclination can be very high, at least for close-in planets. One process to alter the inclination of a planet is through planet-disc interactions. Recent simulations considering radiative transport have shown that the evolution of migration and eccentricity can strongly depend on the thermodynamic state of the disc. We extend previous studies to investigate the planet-disc interactions of fixed and moving planets on inclined and eccentric orbits. We also analyse the effect of the disc's thermodynamic properties on the orbital evolution of embedded planets in detail. The protoplanetary disc is modelled as a viscous gas where the internally produced dissipation is transported by radiation. For locally isothermal discs, we confirm previous results and find inclination damping and inward migration for planetary cores. For low inclinations i < 2 H/r, the damping is exponential, while di/dt is proportional to i^-2 for larger i. For radiative discs, the planetary migration is very limited, as long as their inclination exceeds a certain threshold. If the inclination is damped below this threshold, planetary cores with a mass up to approximately 33 Earth masses start to migrate outwards, while larger cores migrate inwards right from the start. The inclination is damped for all analysed planet masses. In a viscous disc an initial inclination of embedded planets will be damped for all planet masses. This damping occurs on timescales that are shorter than the migration time. If the inclination lies beneath a certain threshold, the outward migration in radiative discs is not handicapped. Outward migration is strongest for circular and non-inclined orbits.

Context. Abridged. Many stars are members of binary systems. During early phases when the stars are surrounded by discs, the binary orbit and disc midplane may be mutually inclined. The discs around T Tauri stars will become mildly warped and undergo solid body precession around the angular momentum vector of the binary system. It is unclear how planetesimals in such a disc will evolve and affect planet formation. Aims. We investigate the dynamics of planetesimals embedded in discs that are perturbed by a binary companion on a circular, inclined orbit. We examine collisional velocities of the planetesimals to determine when they can grow through accretion. We vary the binary inclination, binary separation, D, disc mass, and planetesimal radius. Our standard model has D=60 AU, inclination=45 deg, and a disc mass equivalent to the MMSN. Methods. We use a 3D hydrodynamics code to model the disc. Planetesimals are test particles which experience gas drag, the gravitational force of the disc, the companion star gravity. Planetesimal orbit crossing events are detected and used to estimate collisional velocities. Results. For binary systems with modest inclination (25 deg), disc gravity prevents planetesimal orbits from undergoing strong differential nodal precession (which occurs in absence of the disc), and forces planetesimals to precess with the disc on average. For bodies of different size the orbit planes become modestly mutually inclined, leading to collisional velocities that inhibit growth. For larger inclinations (45 degrees), the Kozai effect operates, leading to destructively large relative velocities. Conclusions. Planet formation via planetesimal accretion is difficult in an inclined binary system with parameters similar to those considered in this paper. For systems in which the Kozai mechanism operates, the prospects for forming planets are very remote.

Planet-disk interaction predicts a change in the orbital elements of an embedded planet. Through linear and fully hydrodynamical studies it has been found that migration is typically directed inwards. Hence, this migration process gives natural explanation for the presence of the 'hot' planets orbiting close to the parent star, and it plays a mayor role in explaining the formation of resonant planetary systems. However, standard migration models for locally isothermal disks indicate a too rapid inward migration for small mass planets, and a large number of massive planets are found very far away from the star. Recent studies, including more complete disk physics, have opened up new paths to slow down or even reverse migration. The new findings on migration are discussed and connected to the observational properties of planetary systems.

We have detected 2 circumbinary planets around the close binary system NN Serpentis using the orbital time delay effect measured via the sharp eclipses of the white dwarf primary. The present pre-cataclysmic binary was formed when the original - 2 M primary expanded into a red giant, causing the secondary star to drop from its original orbit at a separation of about 1.4 A.U. down to its current separation at 0.0043 A.U. A quasi-adiabatic evolution of the circumbinary planets' orbits during the common-envelope phase would have placed them in unstable configurations, suggesting that they may have suffered significant orbital drag effects and were originally in much larger orbits. Alternatively, they may have been created as 2nd-generation planets during the last million years from the substantial amount of material lost during the creation of the binary, making them the youngest planets known. Either solution shows how little we actually understand about planetary formation.

Planets orbiting post-common envelope binaries provide fundamental information on planet formation and evolution. We searched for such planets in NN Ser ab, an eclipsing short-period binary that shows long-term eclipse time variations. Using published, reanalysed, and new mid-eclipse times of NN Ser ab obtained between 1988 and 2010, we find excellent agreement with the light-travel-time effect by two additional bodies superposed on the linear ephemeris of the binary. Our multi-parameter fits accompanied by N-body simulations yield a best fit for the objects NN Ser (ab)c and d locked in a 2:1 mean motion resonance, with orbital periods P_c=15.5 yrs and P_d=7.7 yrs, masses M_c sin i_c = 6.9 M_Jup and M_d sin i_d = 2.2 M_Jup, and eccentricities e_c=0 and e_d=0.20. A secondary chi**2 minimum corresponds to an alternative solution with a period ratio of 5:2. We estimate that the progenitor binary consisted of an A star with ~2 M_Sun and the present M dwarf secondary at an orbital separation of ~1.5 AU. The survival of two planets through the common-envelope phase that created the present white dwarf requires fine tuning between the gravitational force and the drag force experienced by them in the expanding envelope. The alternative is a second-generation origin in a circumbinary disk created at the end of this phase. In that case, the planets would be extremely young with ages not exceeding the cooling age of the white dwarf of 10**6 yrs.

With an average eccentricity of about 0.29, the eccentricity distribution of extrasolar planets is markedly different from the solar system. Among other scenarios considered, it has been proposed that eccentricity may grow through planet-disc interaction. Recently, it has been noticed that the thermodynamical state of the disc can significantly influence the migration properties of growing protoplanets. However, the evolution of planetary eccentricity in radiative discs has not been considered yet. In this paper we study the evolution of planets on eccentric orbits that are embedded in a three-dimensional viscous disc and analyse the disc's effect on the orbital evolution of the planet. We use the three-dimensional hydrodynamical code NIRVANA that includes full tensor viscosity and implicit radiation transport in the flux-limited diffusion approximation. The code uses the FARGO-algorithm to speed up the simulations. First we measure the torque and power exerted on the planet by the disc for fixed orbits, and then we let the planet start with initial eccentricity and evolve it in the disc. For locally isothermal we confirm previous results and find eccentricity damping and inward migration for planetary cores. In the case of radiative discs, the planets experience an inward migration as long as its eccentricity lies above a certain threshold. After the damping of eccentricity cores with masses below 33 Earthmasses begin to migrate outward in radiative discs, while higher mass cores always migrate inward. For all planetary masses studied (up to 200 Earthmasses) we find eccentricity damping. In viscous discs the orbital eccentricity of embedded planets is damped during the evolution independent of the mass. Hence, planet-disc interaction does not seem to be a viable mechanism to explain the observed high eccentricity of exoplanets.

We study the effects of diffusion on the non-linear corotation torque, or horseshoe drag, in the two-dimensional limit, focusing on low-mass planets for which the width of the horseshoe region is much smaller than the scale height of the disc. In the absence of diffusion, the non-linear corotation torque saturates, leaving only the Lindblad torque. Diffusion of heat and momentum can act to sustain the corotation torque. In the limit of very strong diffusion, the linear corotation torque is recovered. For the case of thermal diffusion, this limit corresponds to having a locally isothermal equation of state. We present some simple models that are able to capture the dependence of the torque on diffusive processes to within 20% of the numerical simulations.

We investigate a massive ($\varSigma \sim 10000 g cm^{-2}$ at 1 AU) protoplanetary disc model by means of 3D radiation magnetohydrodynamics simulations. The vertical structure of the disc is determined self-consistently by a balance between turbulent heating caused by the MRI and radiative cooling. Concerning the vertical structure, two different regions can be distinguished: A gas-pressure dominated, optically thick midplane region where most of the dissipation takes place, and a magnetically dominated, optically thin corona which is dominated by strong shocks. At the location of the photosphere, the turbulence is supersonic ($M \sim 2$), which is consistent with previous results obtained from the fitting of spectra of YSOs. It is known that the turbulent saturation level in simulations of MRI-induced turbulence does depend on numerical factors such as the numerical resolution and the box size. However, by performing a suite of runs at different resolutions (using up to 64x128x512 grid cells) and with varying box sizes (with up to 16 pressure scaleheights in the vertical direction), we find that both the saturation levels and the heating rates show a clear trend to converge once a sufficient resolution in the vertical direction has been achieved.

Among multi-planet planetary systems there are a large fraction of resonant systems. Studying the dynamics and formation of these systems can provide valuable informations on processes taking place in protoplanetary disks where the planets are thought have been formed. The recently discovered resonant system HD 60532 is the only confirmed case, in which the central star hosts a pair of giant planets in 3:1 mean motion resonance. We intend to provide a physical scenario for the formation of HD 60532, which is consistent with the orbital solutions derived from the radial velocity measurements. Observations indicate that the system is in an antisymmetric configuration, while previous theoretical investigations indicate an asymmetric equilibrium state. The paper aims at answering this discrepancy as well. We performed two-dimensional hydrodynamical simulations of thin disks with an embedded pair of massive planets. Additionally, migration and resonant capture are studied by gravitational N-body simulations that apply properly parametrized non-conservative forces. Our simulations suggest that the capture into the 3:1 mean motion resonance takes place only for higher planetary masses, thus favouring orbital solutions having relatively smaller inclination i=20 degrees. The system formed by numerical simulations qualitatively show the same behaviour as HD 60532. We also find that the presence of an inner disk (between the inner planet and the star) plays a very important role in determining the final configurations of resonant planetary systems. Its damping effect on the inner planet's eccentricity is responsible for the observed antisymmetric state of HD 60532.

(Note: this is a shortened version of the original A&A-style structured abstract). The physical nature of the strong photometric variability of T Tau Sa, the more massive member of the Southern "infrared companion" to T Tau, has long been debated. Intrinsic luminosity variations due to variable accretion were originally proposed but later challenged in favor of apparent fluctuations due to time-variable foreground extinction. In this paper we use the timescale of the variability as a diagnostic for the underlying physical mechanism. Because the IR emission emerging from Sa is dominantly thermal emission from circumstellar dust at <=1500K, we can derive a minimum size of the region responsible for the time-variable emission. In the context of the variable foreground extinction scenario, this region must be (un-) covered within the variability timescale, which implies a minimum velocity for the obscuring foreground material. If this velocity supercedes the local Kepler velocity we can reject foreground extinction as a valid variability mechanism. The variable accretion scenario allows for shorter variability timescales since the variations in luminosity occur on much smaller scales, essentially at the surface of the star, and the disk surface can react almost instantly on the changing irradiation with a higher or lower dust temperature and according brightness. We have detected substantial variations at long wavelengths in T Tau S: +26% within four days at 12.8 micron. We show that this short-term variability cannot be due to variable extinction and instead must be due to variable accretion. Using a radiative transfer model of the Sa disk we show that variable accretion can in principle also account for the much larger (several magnitude) variations observed on timescales of several years. For the long-term variability, however, also variable foreground extinction is a viable mechanism.

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.

The recently discovered planetary system HD45364 which consists of a Jupiter and Saturn mass planet is very likely in a 3:2 mean motion resonance. The standard scenario to form planetary commensurabilities is convergent migration of two planets embedded in a protoplanetary disc. When the planets are initially separated by a period ratio larger than two, convergent migration will most likely lead to a very stable 2:1 resonance for moderate migration rates. To avoid this fate, formation of the planets close enough to prevent this resonance may be proposed. However, such a simultaneous formation of the planets within a small annulus, seems to be very unlikely. Rapid type III migration of the outer planet crossing the 2:1 resonance is one possible way around this problem. In this paper, we investigate this idea in detail. We present an estimate for the required convergent migration rate and confirm this with N-body and hydrodynamical simulations. If the dynamical history of the planetary system had a phase of rapid inward migration that forms a resonant configuration, we predict that the orbital parameters of the two planets are always very similar and hence should show evidence of that. We use the orbital parameters from our simulation to calculate a radial velocity curve and compare it to observations. Our model can explain the observational data as good as the previously reported fit. The eccentricities of both planets are considerably smaller and the libration pattern is different. Within a few years, it will be possible to observe the planet-planet interaction directly and thus distinguish between these different dynamical states.

To explain important properties of extrasolar planetary systems (eg. close-in hot Jupiters, resonant planets) an evolutionary scenario which allows for radial migration of planets in disks is required. During their formation protoplanets undergo a phase in which they are embedded in the disk and interact gravitationally with it. This planet-disk interaction results in torques (through gravitational forces) acting on the planet that will change its angular momentum and result in a radial migration of the planet through the disk. To determine the outcome of this very important process for planet formation, dedicated high resolution numerical modeling is required. This contribution focusses on some important aspects of the numerical approach that we found essential for obtaining successful results. We specifically mention the treatment of Coriolis forces, Cartesian grids, and the FARGO method.

As of today over 40 planetary systems have been discovered in binary star systems. In all cases the configuration appears to be circumstellar, where the planets orbit around one of the stars, the secondary acting as a perturber. The formation of planets in binary star systems is more difficult than around single stars due to the gravitational action of the companion on the dynamics of the protoplanetary disk. In this contribution we first briefly present the relevant observational evidence for planets in binary systems. Then the dynamical influence that a secondary companion has on a circumstellar disk will be analyzed through fully hydrodynamical simulations. We demonstrate that the disk becomes eccentric and shows a coherent precession around the primary star. Finally, fully hydrodynamical simulations of evolving protoplanets embedded in disks in binary star systems are presented. We investigate how the orbital evolution of protoplanetary embryos and their mass growth from cores to massive planets might be affected in this very dynamical environment. We consider, in particular, the planet orbiting the primary in the system Gamma Cephei.

We study the torque on low-mass planets embedded in protoplanetary discs in the two-dimensional approximation, incorporating non-isothermal effects. We couple linear estimates of the Lindblad (or wave) torque to a simple, but non-linear, model of adiabatic corotation torques (or horseshoe drag), resulting in a simple formula that governs Type I migration in non-isothermal discs. This formula should apply in optically thick regions of the disc, where viscous and thermal diffusion act to keep the horseshoe drag unsaturated. We check this formula against numerical hydrodynamical simulations, using three independent numerical methods, and find good agreement.

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.

Numerical simulations of planets embedded in protoplanetary gaseous discs are a precious tool for studying the planetary migration ; however, some approximations have to be made. Most often, the selfgravity of the gas is neglected. In that case, it is not clear in the literature how the material inside the Roche lobe of the planet should be taken into account. Here, we want to address this issue by studying the influence of various methods so far used by different authors on the migration rate. We performed high-resolution numerical simulations of giant planets embedded in discs. We compared the migration rates with and without gas selfgravity, testing various ways of taking the circum-planetary disc (CPD) into account. Different methods lead to significantly different migration rates. Adding the mass of the CPD to the perturbing mass of the planet accelerates the migration. Excluding a part of the Hill sphere is a very touchy parameter that may lead to an artificial suppression of the type III, runaway migration. In fact, the CPD is smaller than the Hill sphere. We recommend excluding no more than a 0.6 Hill radius and using a smooth filter. Alternatively, the CPD can be given the acceleration felt by the planet from the rest of the protoplanetary disc. The gas inside the Roche lobe of the planet should be very carefully taken into account in numerical simulations without any selfgravity of the gas. The entire Hill sphere should not be excluded. The method used should be explicitly given. However, no method is equivalent to computing the full selfgravity of the gas.

In this paper we investigate the influence of radiative transport on the growth of the magnetorotational instability (MRI) in accretion discs. The analysis is performed by use of analytical and numerical means. We provide a general dispersion relation together with the corresponding eigenfunctions describing the growth rates of small disturbances on a homogeneous background shear flow. The dispersion relation includes compressibility and radiative effects in the flux-limited diffusion approximation. By introducing an effective speed of sound, all the effects of radiation transport can be subsumed into one single parameter. It can be shown that the growth rates of the vertical modes -- which are the fastest growing ones -- are reduced by radiative transport. For the case of non-vertical modes, the growth rates may instead be enhanced. We quantify the effects of compressibility and radiative diffusion on the growth rates for the gas-pressure dominated case. The analytical discussion is supplemented by numerical simulations, which are also used for a first investigation of the non-linear stage of the MRI in gas-pressure dominated accretion discs with radiation transport included.

The formation of resonant planets pairs in exoplanetary systems involves planetary migration inside the protoplanetary disc : an inwards migrating outer planet captures in Mean Motion Resonance an inner planet. During the migration of the resonant pair of planets, the eccentricities are expected to rise excessively, if no damping mechanism is applied on the inner planet. We express the required damping action to match the observations, and we show that the inner disk can play this role. This result applies for instance to the system GJ 876 : we reproduce the observed orbital elements through a fully hydrodynamical simulation of the evolution of the resonant planets.

We study the development of finite eccentricity in accretion disks in close binary systems using a two-dimensional grid-based numerical scheme. We perform detailed parameter studies to explore the dependence on viscosity, disk aspect ratio, the inclusion of a mass-transfer stream and the role of the boundary conditions. We consider mass ratios 0.05<q<0.3 appropriate to superoutbursting cataclysmic binary systems. Instability to the formation of a precessing eccentric disk that attains a quasi-steady state with mean eccentricity in the range 0.3-0.5 occurs readily. The shortest growth times are ~15 binary orbits for the largest viscosities and the instability mechanism is for the most part consistent with the mode-coupling mechanism associated with the 3:1 resonance proposed by Lubow. However, the results are sensitive to the treatment of the inner boundary and to the incorporation of the mass-transfer stream. In the presence of a stream we found a critical viscosity below which the disk remains circular. Incorporation of a mass-transfer stream tends to impart stability for small enough viscosity (or, equivalently, mass-transfer rate through the disk) and does assist in obtaining a prograde precession rate that is in agreement with observations. For the larger q the location of the 3:1 resonance is pushed outwards towards the Roche lobe where higher-order mode couplings and nonlinearity occur. It is likely that three-dimensional simulations that properly resolve the disk's vertical structure are required to make significant progress in this case.

In isothermal disks the migration of protoplanets is directed inward. For small planetary masses the standard type-I migration rates are so fast that this may result in an unrealistic loss of planets into the stars. We investigate the planet-disk interaction in non-isothermal disks and analyze the magnitude and direction of migration for an extended range of planet masses. We have performed detailed two-dimensional numerical simulations of embedded planets including heating/cooling effects as well as radiative diffusion for realistic opacities. In radiative disks, small planets with M_planet < 50 M_Earth do migrate outward with a rate comparable to absolute magnitude of standard type-I migration. For larger masses the migration is inward and approaches the isothermal, type-II migration rate. Our findings are particularly important for the first growth phase of planets and ease the problem of too rapid inward type-I migration.

Over 30 planetary systems have been discovered to reside in binary stars. For small separations gravitational perturbation of the secondary star has a strong influence on the planet formation process. It truncates the protoplanetary disk, may shortens its lifetime, and stirs up the embedded planetesimals. Due to its small semi-major axis (18.5 AU) and large eccentricity (e=0.35) the binary $\gamma$ Cephei represents a particularly challenging example. In the present study we model the orbital evolution and growth of embedded protoplanetary cores of about 30 earth masses in the putative protoplanetary disk surrounding the primary star in the $\gamma$ Cep system. We assume coplanarity of the disk, binary and planet and perform two-dimensional hydrodynamic simulations of embedded cores in a protoplanetary disk. The presence of the eccentric secondary star perturbs the disk periodically and generates strong spiral arms at periapse which propagate toward the disk centre. The disk also becomes slightly eccentric (with e_d = 0.1-0.15), and displays a slow retrograde precession in the inertial frame. For all initial separations (2.5 to 3.5 AU) we find inward migration of the cores. For initial semi-major axes (a_p \gsim 2.7), we find a strong increase in the planetary eccentricity despite the presence of inward migration. Only cores which are initially far from the disk outer edge have a bounded orbital eccentricity which converges, roughly to the value of the planet observed in the $\gamma$ Cep system. We have shown that under the condition protoplanetary cores can form at around 2.5 AU, it is possible to evolve and grow such a core to form a planet with final outcome similar to that observed.

The formation of resonant pairs of planets in exoplanetary systems involves planetary migration in the protoplanetary disc. After a resonant capture, the subsequent migration in this configuration leads to a large increase of planetary eccentricities if no damping mechanism is applied. This has led to the conclusion that the migration of resonant planetary systems cannot occur over large radial distances and has to be terminated sufficiently rapidly through disc dissipation. In this study, we investigate whether the presence of an inner disc might supply an eccentricity damping of the inner planet, and if this effect could explain the observed eccentricities in some systems. To investigate the influence of an inner disc, we first compute hydrodynamic simulations of giant planets orbiting with a given eccentricity around an inner gas disc, and measure the effect of the latter on the planetary orbital parameters. We then perform detailed long term calculations of the GJ 876 system. We also run N-body simulations with artificial forces on the planets mimicking the effects of the inner and outer discs. We find that the influence of the inner disc can not be neglected, and that it might be responsible for the observed eccentricities. In particular, we reproduce quite well the orbital parameters of a few systems engaged in 2:1 mean motion resonances: GJ 876, HD 73526, HD 82943 and HD 128311. Finally, we derive analytically the effect that the inner disc should have on the inner planet to reach a specific orbital configuration with a given damping effect of the outer disc on the outer planet.

Young planets embedded in their protoplanetary disk interact gravitationally with it leading to energy and angular momentum exchange. This interaction determines the evolution of the planet through changes to the orbital parameters. We investigate changes in the orbital elements of a 20 Earth--mass planet due to the torques from the disk. We focus on the non-linear evolution of initially non-vanishing eccentricity $e$ and/or inclination $i$. We treat the disk as a two- or three-dimensional viscous fluid and perform hydrodynamical simulations with an embedded planet. We find rapid exponential decay of the planet orbital eccentricity and inclination for small initial values of $e$ and $i$, in agreement with linear theory. For larger values of $e > 0.1$ the decay time increases and the decay rate scales as $\dot{e} \propto e^{-2}$, consistent with existing theoretical models. For large inclinations ($i$ > 6 deg) the inclination decay rate shows an identical scaling $di/dt \propto i^{-2}$. We find an interesting dependence of the migration on the eccentricity. In a disk with aspect ratio $H/r=0.05$ the migration rate is enhanced for small non-zero eccentricities ($e < 0.1$), while for larger values we see a significant reduction by a factor of $\sim 4$. We find no indication for a reversal of the migration for large $e$, although the torque experienced by the planet becomes positive when $e \simeq 0.3$. This inward migration is caused by the persisting energy loss of the planet. For non gap forming planets, eccentricity and inclination damping occurs on a time scale that is very much shorter than the migration time scale. The results of non linear hydrodynamic simulations are in very good agreement with linear theory for small $e$ and $i$.

Based on radial velocity measurements it has been found recently that the two giant planets detected around the star HD 73526 are in 2:1 resonance. However, as our numerical integration shows, the derived orbital data for this system result in chaotic behavior of the giant planets, which is uncommon among the resonant extrasolar planetary systems. We intend to present regular (non-chaotic) orbital solutions for the giant planets in the system HD 73526 and offer formation scenarios based on combining planetary migration and sudden perturbative effects such as planet-planet scattering or rapid dispersal of the protoplanetary disk. A comparison with the already studied resonant system HD 128311, exhibiting similar behavior, is also done. The new sets of orbital solutions have been derived by the Systemic Console (www.oklo.org). The stability of these solutions has been investigated by the Relative Lyapunov indicator, while the migration and scattering effects are studied by gravitational N-body simulations applying non-conservative forces as well. Additionally, hydrodynamic simulations of embedded planets in protoplanetary disks are performed to follow the capture into resonance. For the system HD 73526 we demonstrate that the observational radial velocity data are consistent with a coplanar planetary system engaged in a stable 2:1 resonance exhibiting apsidal corotation. We have shown that, similarly to the system HD 128311, the present dynamical state of HD 73526 could be the result of a mixed evolutionary process melting together planetary migration and a perturbative event.

Among the extrasolar planetary systems about 30 are located in a stellar binary orbiting one of the stars, preferably the more massive primary. The dynamical influence of the second companion alters firstly the orbital elements of the forming protoplanet directly and secondly the structure of the disk from which the planet formed which in turn will modify the planet's evolution. We present detailed analysis of these effects and present new hydrodynamical simulations of the evolution of protoplanets embedded in circumstellar disks in the presence of a companion star, and compare our results to the system $\gamma$ Cep. To analyse the early formation of planetary embryos, we follow the evolution of a swarm of planetesimals embedded in a circumstellar disk. Finally, we study the evolution of planets embedded in circumbinary disks.

Following the recent insight in the material structure of comets, protoplanetesimals are assumed to have low densities and to be highly porous agglomerates. It is still unclear if planetesimals can be formed from these objects by collisional growth. Therefore, it is important to study numerically the collisional outcome from low velocity impacts of equal sized porous agglomerates which are too large to be examined in a laboratory experiment. We use the Lagrangian particle method Smooth Particle Hydrodynamics to solve the equations that describe the dynamics of elastic and plastic bodies. Additionally, to account for the influence of porosity, we follow a previous developed equation of state and certain relations between the material strength and the relative density. Collisional growth seems possible for rather low collision velocities and particular material strengths. The remnants of collisions with impact parameters that are larger than 50% of the radius of the colliding objects tend to rotate. For small impact parameters, the colliding objects are effectively slowed down without a prominent compaction of the porous structure, which probably increases the possibility for growth. The protoplanetesimals, however, do not stick together for the most part of the employed material strengths. An important issue in subsequent studies has to be the influence of rotation to collisional growth. Moreover, for realistic simulations of protoplanetesimals it is crucial to know the correct material parameters in more detail.

The increase of computational resources has recently allowed high resolution, three dimensional calculations of planets embedded in gaseous protoplanetary disks. They provide estimates of the planet migration timescale that can be compared to analytical predictions. While these predictions can result in extremely short migration timescales for cores of a few Earth masses, recent numerical calculations have given an unexpected outcome: the torque acting on planets with masses between 5 M_Earth and 20 M_Earth is considerably smaller than the analytic, linear estimate. These findings motivated the present work, which investigates existence and origin of this discrepancy or ``offset'', as we shall call it, by means of two and three dimensional numerical calculations. We show that the offset is indeed physical and arises from the coorbital corotation torque, since (i) it scales with the disk vortensity gradient, (ii) its asymptotic value depends on the disk viscosity, (iii) it is associated to an excess of the horseshoe zone width. We show that the offset corresponds to the onset of non-linearities of the flow around the planet, which alter the streamline topology as the planet mass increases: at low mass the flow non-linearities are confined to the planet's Bondi sphere whereas at larger mass the streamlines display a classical picture reminiscent of the restricted three body problem, with a prograde circumplanetary disk inside a ``Roche lobe''. This behavior is of particular importance for the sub-critical solid cores (M <~ 15 M_Earth) in thin (H/r <~0.06) protoplanetary disks. Their migration could be significantly slowed down, or reversed, in disks with shallow surface density profiles.

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.

A significant number of the known multiple exoplanetary systems are containing a pair of giant planets engaged in a low order mean motion resonance. Such a resonant condition protects the dynamics of these planets resulting in very stable orbits. According to recent studies the capture into a resonance is the result of a planetary migration process induced by the interaction of the planets with a protoplanetary disk. If the migration is slow enough (adiabatic) next to a mean motion resonance, the two planets will also be in apsidal corotation. The recently refined orbital parameters of the system HD 128311 suggest that the two giant planets are in a 2:1 mean motion resonance, however without exhibiting apsidal corotation. Thus the evolution of this system can not be described by an adiabatic migration process alone. We present possible evolution scenarios of this system combining migration processes and sudden perturbations. We model migration scenarios through numerical integration of the gravitational N-body problem with additional non-conservative forces. Planet-planet scattering has been investigated by N-body simulations. We show that the present dynamical state of the system HD128311 may be explained by such evolutionary processes.

The discovery of close orbiting extrasolar giant planets led to extensive studies of disk planet interactions and the forms of migration that can result as a means of accounting for their location. Early work established the type I and type II migration regimes for low mass embedded planets and high mass gap forming planets respectively. While providing an attractive means of accounting for close orbiting planets intially formed at several AU, inward migration times for objects in the earth mass range were found to be disturbingly short, making the survival of giant planet cores an issue. Recent progress in this area has come from the application of modern numerical techniques which make use of up to date supercomputer resources. These have enabled higher resolution studies of the regions close to the planet and the initiation of studies of planets interacting with disks undergoing MHD turbulence. This work has led to indications of how the inward migration of low to intermediate mass planets could be slowed down or reversed. In addition, the possibility of a new very fast type III migration regime, that can be directed inwards or outwards, that is relevant to partial gap forming planets in massive disks has been investigated.

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 investigate the response of an accretion disk to the presence of a perturbing protoplanet embedded in the disk through time dependent hydrodynamical simulations. The disk is treated as a two-dimensional viscous fluid and the planet is kept on a fixed orbit. We run a set of simulations varying the planet mass, and the viscosity and temperature of the disk. All runs are followed until they reach a quasi-equilibrium state. We find that for planetary masses above a certain minimum mass, already 3 M_Jup for a viscosity of nu = 10^{-5}, the disk makes a transition from a nearly circular state into an eccentric state. Increasing the planetary mass leads to a saturation of disk eccentricity with a maximum value of around 0.25. The transition to the eccentric state is driven by the excitation of an m=2 spiral wave at the outer 1:3 Lindblad resonance. The effect occurs only if the planetary masses are large enough to clear a sufficiently wide and deep gap to reduce the damping effect of the outer 1:2 Lindblad resonance. An increase in viscosity and temperature in the disk, which both tend to close the gap, have an adverse influence on the disk eccentricity. In the eccentric state the mass accretion rate onto the planet is greatly enhanced, an effect that may ease the formation of massive planets beyond about 5 M_Jup that are otherwise difficult to reach.

The two planets about the star GJ 876 appear to have undergone extensive migration from their point of origin in the protoplanetary disk -- both because of their close proximity to the star (30 and 60 day orbital periods) and because of their occupying three stable orbital resonances at the 2:1 mean-motion commensurability. The resonances were most likely established by converging differential migration of the planets leading to capture into the resonances. A problem with this scenario is that continued migration of the system while it is trapped in the resonances leads to orbital eccentricities that rapidly exceed the observational upper limits of e_1 = 0.31 and e_2 = 0.05. As seen in forced 3-body simulations, lower eccentricities would persist during migration only for an applied eccentricity damping. Here we explore the evolution of the GJ 876 system using two-dimensional hydrodynamical simulations that include viscous heating and radiative effects. We find that a hydrodynamic evolution within the resonance, where only the outer planet interacts with the disk, always rapidly leads to large values of eccentricities that exceed those observed. Only if mass is removed from the disk on a time scale of the order of the migration time scale (before there has been extensive migration after capture), as might occur for photoevaporation in the late phases of planet formation, can we end up with eccentricities that are consistent with the observations.

We study the evolution and emission of circumbinary disks around close classical T Tauri binary systems. High resolution numerical hydrodynamical simulations are employed to model a system consisting of a central eccentric binary star within an irradiated accretion disk. A detailed energy balance including viscous heating, radiative cooling and irradiation from the central star is applied to calculate accurately the emitted spectral energy distribution. Numerical simulations using two different methods, the previously developed Dual-Grid technique with a finite difference discretization, and the Smoothed Particle Hydrodynamics method are employed to compare the hydrodynamical features and strengthen our conclusions. Physical parameters of the setup are chosen to model the close systems of DQ Tau and AK Sco. Using the self-consistent models, we are able to fit the observed spectral energy distributions by constraining parameters such as disk mass, density profile and radial extension for those systems. We find that the incorporation of irradiation effects is necessary to obtain correct disk temperatures.

We have performed Smoothed Particle Hydrodynamics (SPH) simulations to study the time evolution of one and two protoplanets embedded in a protoplanetary accretion disc. We investigate accretion and migration rates of a single protoplanet depending on several parameters of the protoplanetary disc, mainly viscosity and scale height. Additionally, we consider the influence of a second protoplanet in a long time simulation and examine the migration of the two planets in the disc, especially the growth of eccentricity and chaotic behaviour. One aim of this work is to establish the feasibility of SPH for such calculations considering that usually only grid-based methods are adopted. To resolve shocks and to prevent particle penetration, we introduce a new approach for an artificial viscosity, which consists of an additional artificial bulk viscosity term in the SPH-representation of the Navier-Stokes equation. This allows for an accurate treatment of the physical kinematic viscosity to describe the shear, without the use of an artificial shear viscosity.

The main properties of the observed extrasolar planets are reviewed with respect to their relevance to the formation scenario of planetary systems. Results of numerical computations of embedded planets in viscously evolving disks are presented. Emphasis is given to the accretion and migration process. New calculations on inviscid disks are shown. The second part of the talk concentrates on resonant planetary systems. Among the observed extrasolar systems there are 3 confirmed cases, Gl 876, HD 82943 and 55 Cnc, where the planets orbit their central star in a low order mean motion resonance. Results of numerical simulations modeling the formation and evolution of such systems are presented.

We study the time evolution of two protoplanets still embedded in a protoplanetary disk. The results of two different numerical approaches are presented and compared. In the first approach, the motion of the disk material is computed with viscous hydrodynamical simulations, and the planetary motion is determined by N-body calculations including exactly the gravitational forces exerted by the disk material. In the second approach, only the N-body integration is performed but with additional dissipative forces included such as to mimic the effect of the disk torques acting on the disk. This type of modeling is much faster than the full hydrodynamical simulations, and gives comparative results provided that parameters are adjusted properly. Resonant capture of the planets is seen in both approaches, where the order of the resonance depends on the properties of the disk and the planets. Resonant capture leads to a rise in the eccentricity and to an alignment of the spatial orientation of orbits. The numerical results are compared with the observed planetary systems in mean motion resonance (Gl 867, HD 82943, and 55 Cnc). We find that the forcing together of two planets by their parent disk produces resonant configurations similar to those observed, but that eccentricity damping greater than that obtained in our hydrodynamic simulations is required to match the GJ 876 observations.

With a series of numerical simulations, we analyze the thermo-hydrodynamical evolution of circumstellar disks containing Jupiter-size protoplanets. In the framework of the two-dimensional approximation, we consider an energy equation that includes viscous heating and radiative effects in a simplified, yet consistent form. Multiple nested grids are used in order to study both global and local features around the planet. By means of different viscosity prescriptions, we investigate various temperature regimes. A planetary mass range from 0.1 to 1 Mj is examined. Computations show that gap formation is a general property which affects density, pressure, temperature, optical thickness, and radiated flux distributions. However, it remains a prominent feature only when the kinematic viscosity is on the order of 10^(15) cm^2/s or lower. Around accreting planets, a circumplanetary disk forms that has a surface density profile decaying exponentially with the distance and whose mass is 5-6 orders of magnitudes smaller than Jupiter's mass. Circumplanetary disk temperature profiles decline roughly as the inverse of the distance from the planet. Temperatures range from some 10 to ~1000 K. Planetary accretion and migration rates depend on the viscosity regime, with discrepancies within an order of magnitude. Estimates of growth and migration time scales inferred by these models are on the same orders of magnitude as those previously obtained with locally isothermal simulations.

We investigate the evolution of protoplanets with different masses embedded in an accretion disk, via global fully three-dimensional hydrodynamical simulations. We consider a range of planetary masses extending from one and a half Earth's masses up to one Jupiter's mass, and we take into account physically realistic gravitational potentials of forming planets. In order to calculate accurately the gravitational torques exerted by disk material and to investigate the accretion process onto the planet, the flow dynamics has to be thoroughly resolved on long as well as short length scales. We achieve this strict resolution requirement by applying a nested-grid refinement technique which allows to greatly enhance the local resolution. Our results from altogether 51 simulations show that for large planetary masses, approximately above a tenth of the Jupiter's mass, migration rates are relatively constant, as expected in type II migration regime and in good agreement with previous two-dimensional calculations. In a range between seven and fifteen Earth's masses, we find a dependency of the migration speed on the planetary mass that yields time scales considerably longer than those predicted by linear analytical theories. This property may be important in determining the overall orbital evolution of protoplanets. The growth time scale is minimum around twenty Earth-masses, but it rapidly increases for both smaller and larger mass values. Significant differences between two- and three-dimensional calculations are found in particular for objects with masses smaller than ten Earth-masses. We also derive an analytical approximation for the numerically computed mass growth rates.

Recently, 59 low-luminosity object transits were reported from the Optical Gravitational Lensing Experiment (OGLE). Our follow-up low-resolution spectroscopy of 16 candidates provided two objects, OGLE-TR-3 and OGLE-TR-10, which have companions with radii compatible with those of gas-giant planets. Further high-resolution spectroscopy revealed a very low velocity variation (<500m/s) of the host star OGLE-TR-3 which may be caused by its unseen companion. An analysis of the radial velocity and light curve results in M<2.5 M_jup, R<1.6 R_jup, and an orbital separation of about 5 R_sol, which makes it the planet with the shortest period known. This allows to identify the low-luminosity companion of OGLE-TR-3 as a possible new gas-giant planet. If confirmed, this makes OGLE-TR-3 together with OGLE-TR-56 the first extrasolar planets detected via their transit light curves.

We study the evolution of a system consisting of two protoplanets still embedded in a protoplanetary disk. Results of two different numerical approaches are presented. In the first kind of model the motion of the disk material is followed by fully viscous hydrodynamical simulations, and the planetary motion is determined by N-body calculations including exactly the gravitational potential from the disk material. In the second kind we only solve the N-body part and add additional analytically given forces which model the effect of the torques of the disk. This type of modeling is of course orders of magnitudes faster than the full hydro-model. Another advantage of this two-fold approach is the possibility of adjusting the otherwise unknown parameters of the simplified model. The results give very good agreement between the methods. Using two different initial setups for the planets and disk, we obtain in the first case a resonant trapping into the 3:1 resonance, and in the second case a trapping into the 2:1 resonance. Resonant capture leads to a rise in the eccentricity and to an alignment of the of the spatial orientation of orbits. The characteristics of the numerical results agree very favorably with those of 3 observed planetary systems (GJ 876, HD 82943, and 55 Cnc) known to be in mean motion resonances.

We study analytically and numerically the stability of the pressure-less, viscously spreading accretion ring. We show that the ring is unstable to small non-axisymmetric perturbations. To perform the perturbation analysis of the ring we use a stretching transformation of the time coordinate. We find that to 1st order, one-armed spiral structures, and to 2nd order additionally two-armed spiral features may appear. Furthermore, we identify a dispersion relation determining the instability of the ring. The theoretical results are confirmed in several simulations, using two different numerical methods. These computations prove independently the existence of a secular spiral instability driven by viscosity, which evolves into persisting leading and trailing spiral waves. Our results settle the question whether the spiral structures found in earlier simulations of the spreading ring are numerical artifacts or genuine instabilities.

Planet evolution is tightly connected to the dynamics of both distant and close disk material. Hence, an appropriate description of disk-planet interaction requires global and high resolution computations, which we accomplish by applying a Nested-Grid method. Through simulations in two and three dimensions, we investigate how migration and accretion are affected by long and short range interactions. For small mass objects, 3D models provide longer growth and migration time scales than 2D ones do, whereas time lengths are comparable for large mass planets.

We investigate the possibility to find evidence for planets in circumstellar disks by infrared and submillimeter interferometry. Hydrodynamical simulations of a circumstellar disk around a solar-type star with an embedded planet of 1 Jupiter mass are presented. On the basis of 3D radiative transfer simulations, images of this system are calculated. These intensity maps provide the basis for the simulation of the interferometers VLTI (equipped with the mid-infrared instrument MIDI) and ALMA. While ALMA will provide the necessary basis for a direct gap and therefore indirect planet detection, MIDI/VLTI will provide the possibility to distinguish between disks with or without accretion on the central star on the basis of visibility measurements.

Recently, 46 low-luminosity object transits were reported from the Optical Gravitational Lensing Experiment. Our follow-up spectroscopy of the 16 most promising candidates provides a spectral classification of the primary. Together with the radius ratio from the transit measurements, we derived the radii of the low-luminosity companions. This allows to examine the possible sub-stellar nature of these objects. Fourteen of them can be clearly identified as low-mass stars. Two objects, OGLE-TR-03 and OGLE-TR-10 have companions with radii of 0.15 R_sun which is very similar to the radius of the transiting planet HD209458B. The planetary nature of these two objects should therefore be confirmed by dynamical mass determinations.

We study the evolution of circumbinary disks surrounding classical T Tau stars. High resolution numerical simulations are employed to model a system consisting of a central eccentric binary star within an accretion disk. The disk is assumed to be infinitesimally thin, however a detailed energy balance including viscous heating and radiative cooling is applied. A novel numerical approach using a parallelized Dual-Grid technique on two different coordinate systems has been implemented. Physical parameters of the setup are chosen to model the close systems of DQ Tau and AK Sco, as well as the wider systems of GG Tau and UY Aur. Our main findings are for the tight binaries a substantial flow of material through the disk gap which is accreted onto the central stars in a phase dependent process. We are able to constrain the parameters of the systems by matching both accretion rates and derived spectral energy distributions to observational data where available.

We investigate the possibility to find evidence for planets in circumstellar disks by infrared and submillimeter interferometry. We present simulations of a circumstellar disk around a solar-type star with an embedded planet of 1 Jupiter mass. The three-dimensional (3D) density structure of the disk results from hydrodynamical simulations. On the basis of 3D radiative transfer simulations, images of this system were calculated. The intensity maps provide the basis for the simulation of the interferometers VLTI (equipped with the mid-infrared instrument MIDI) and ALMA. While MIDI/VLTI will not provide the possibility to distinguish between disks with or without a gap on the basis of visibility measurements, ALMA will provide the necessary basis for a direct gap detection.

We study the evolution of embedded protoplanets in a protostellar disk using very high resolution nested-grid computations. This method allows us to perform global simulations of planets orbiting in disks and, at the same time, to resolve in detail the dynamics of the flow inside the Roche lobe of the planet. The primary interest of this work lies in the analysis of the gravitational torque balance acting on the planet. For this purpose we study planets of different masses, ranging from one Earth-mass up to one Jupiter-mass, assuming typical parameters of the protostellar disk. The high resolution of the method allows a precise determination of the mass flow onto the planet and the resulting torques. The obtained migration time scales are in the range from few times 10^4 years, for intermediate mass planets, to 10^6 years, for very low and high mass planets. Typical growth time scales depend strongly on the planetary mass, ranging from a few hundred years, in the case of Earth-type planets, to several ten thousand years, in the case of Jupiter-type planets.

The orbital parameters of the observed extrasolar planets differ strongly from those of our own solar system. The differences include planets with high masses, small semi-major axis and large eccentricities. We performed numerical computations of embedded planets in disks and follow their mass growth and orbital evolution over several thousand periods. We find that planets do migrate inwards on timescales of about $10^5$ years on nearly circular orbits, during which they may grow up to about 5 Jupiter masses. The interaction of the disk with several planets may halt the migration process and lead to a system similar to the solar planetary system.

The star formation process in molecular clouds usually leads to the formation of multiple stellar systems, mostly binaries. Remaining disks around those stars may be located around individual stars (circumstellar disks) or around the entire binary system (circumbinary disk). We shall briefly review the present observational evidence for both types of disks in binary stars, in particular the properties of circumbinary disks. We then present recent results of the theoretical modeling of the collapse and fragmentation of gravitationally unstable molecular cloud cores and their implications for binary and disk formation, and discuss the dynamical influence of the binary companions on disk truncation and gap formation. The presence of binaries may have profound influence on the process of planet formation as well. We present results on the stability and evolution of orbits of planets in disks around binaries.

The evolution of a system consisting of a protoplanetary disc with two embedded Jupiter sized planets is studied numerically. The disc is assumed to be flat and non-self gravitating, which is modeled by the planar (two-dimensional) Navier-Stokes equations. The mutual gravitational interaction of the planets and the star, and the gravitational torques of the disc acting on the planets and the central star are included. The planets have an initial mass of one Jupiter mass $M_{Jup}$ each and the radial distances from the star are one and two semi-major axis of Jupiter, respectively. During the evolution both planets increase their mass due to accretion of gas from the disc; after about 2500 orbital periods of the inner planet it has reached a mass of 2.3 and the outer planet of 3.2~$M_{Jup}$. The net gravitational torques exerted by the disc on the planets result in an inward migration of the outer planet on time-scales comparable to the viscous evolution time of the disc, while the semi-major axis of the inner planet remains constant. When the distance of close approach eventually becomes smaller than their mutual Hill radius the eccentricities increase strongly and the system may turn unstable. The implications for the origin of the solar system and the extrasolar planets are discussed.

We investigate the gravitational interaction of a Jovian mass protoplanet with a gaseous disc with aspect ratio and kinematic viscosity expected for the protoplanetary disc from which it formed. Different disc surface density distributions have been investigated. We focus on the tidal interaction with the disc with the consequent gap formation and orbital migration of the protoplanet. Nonlinear hydrodynamic simulations are employed using three independent numerical codes. A principal result is that the direction of the orbital migration is always inwards and such that the protoplanet reaches the central star in a near circular orbit after a characteristic viscous time scale of approximately 10,000 initial orbital periods. This was found to be independent of whether the protoplanet was allowed to accrete mass or not. Inward migration is helped through the disappearance of the inner disc, and therefore the positive torque it would exert, because of accretion onto the central star.Our results indicate that a realistic upper limit for the masses of closely orbiting giant planets is approximately 5 Jupiter masses, because of the reduced accretion rates obtained for planets of increasing mass. Assuming some process such as termination of the inner disc through a magnetospheric cavity stops the migration, the range of masses estimated for a number of close orbiting giant planets (Marcy, Cochran, & Mayor 1999; Marcy & Butler 1998) as well as their inward orbital migration can be accounted for by consideration of disc--protoplanet interactions during the late stages of giant planet formation. Maximally accreting protoplanets reached about four Jovian masses on reaching the neighbourhood of the central star.

The discovery of now about 20 extrasolar planets orbiting solar-type stars with properties quite different from those in our Solar System has raised many questions about the formation and evolution of planets. The tidal interaction between the planet and the surrounding disk determines the orbital properties and the mass of the planet. We have performed numerical computations of a planet embedded in a protoplanetary disk and found that for typical values of the viscosity the planet may easily grow upto ten Jupiter masses. New results on the mass evolution and the migration of the planet are presented.

Treating problems in full general relativity is highly complex and frequently approximate methods are employed to simplify the solution. We present comparative solutions of a infinitesimally thin relativistic, stationary, rigidly rotating disk obtained using the full equations and the approximate approach suggested by Wilson & Mathews. We find that the Wilson-Mathews method has about the same accuracy as the first post-Newtonian approximation.

We study the structure and dynamics of the gap created by a protoplanet in an accretion disc. The hydrodynamic equations for a flat, two-dimensional, non-selfgravitating protostellar accretion disc with an embedded, Jupiter sized protoplanet on a circular orbit are solved. To simulate possible accretion of mass onto the protoplanet we continually remove mass from the interior of the planet's Roche lobe which is monitored. Firstly, it is shown that consistent results independent on numerical issues (such as boundary or initial conditions, artificial viscosity or resolution) can be obtained. Then, a detailed parameter study delineates the influence of the disc viscosity and pressure on the magnitude of the accretion rate. We find that, even after the formation of a gap in the disc, the planet is still able to accrete more mass from the disc. This accretion occurs from regions of the disc which are radially exterior and interior to the planet's orbital radius. The rate depends on the magnitude of the viscosity and vertical thickness of the disc. For a disc viscosity alpha=10^{-3} and vertical thickness H/r=0.05 we estimate the time scale for the accumulation of one Jupiter mass to be of order hundred thousand years. For a larger(smaller) viscosity and disc thickness this accretion rate is increasing(decreasing). For a very small viscosity (alpha < 5 10^{-4}) the mass accretion rate through the gap onto the planet is markedly reduced, and the corresponding accretion time scale becomes larger than the viscous evolution time of the disc.

In numerical computations of astrophysical flows one frequently has to deal with situations which require the use of a corotating coordinate system. Then, in addition to the centrifugal force, the Coriolis terms have to be included in the Navier-Stokes equations. In usual numerical treatments of solving the equations of motion these non-inertial forces are included explicitly as source terms. In this letter we show that this procedure, if applied to the angular momentum equation, may lead to erroneous results. Using the specific sample calculation of a protoplanet embedded in a protostellar disk, we then demonstrate that only a conservative treatment of the Coriolis term in the angular momentum equation yields results which are identical to those obtained by computations performed in the inertial frame.

The structure of the boundary layer region between the disc and a comparatively slowly rotating star is studied using a causal prescription for viscosity. The vertically integrated viscous stress relaxes towards its equilibrium value on a relaxation timescale $\tau$, which naturally yields a finite speed of propagation for viscous information. For a standard alpha prescription with alpha in the range 0.1-0.01, and ratio of viscous speed to sound speed in the range 0.02-0.5, details in the boundary layer are strongly affected by the causality constraint. We study both steady state polytropic models and time dependent models, taking into account energy dissipation and transport. Steady state solutions are always subviscous with a variety of $\Omega$ profiles which may exhibit near discontinuities. For alpha =0.01 and small viscous speeds, the boundary layer adjusted to a near steady state. A long wavelength oscillation generated by viscous overstability could be seen at times near the outer boundary. Being confined there, the boundary layer remained almost stationary. However, for alpha =0.1 and large viscous speeds, short wavelength disturbances were seen throughout which could significantly affect the power output in the boundary layer. This could be potentially important in producing time dependent behaviour in accreting systems such as CVs and protostars.

Sequences of infinitesimally thin, uniformly rotating, self-gravitating relativistic discs with internal two-dimensional pressure have been constructed. It is shown that in weaker relativistic configurations the sequences undergo a continuous bifurcation from a disc to a ring structure, while in stronger relativistic cases the sequences terminate at the mass-shed limit where gravitational forces are exactly balanced by centrifugal forces.