Circumplanetary discs (CPDs) form around young gas giants and are thought to be the sites of moon formation as well as an intermediate reservoir of gas that feeds the growth of the gas giant. How the physical properties of such CPDs are affected by the planetary mass and the overall opacity is relatively poorly understood. In order to clarify this, we use the global radiation hydrodynamics code FARGOCA, with a grid structure that allows resolving the planetary gravitational potential sufficiently well for a CPD to form. We then study the gas flows and density/temperature structures that emerge as a function of planet mass, opacity and potential depth. Our results indicate interesting structure formation for Jupiter-mass planets at low opacities, which we subsequently analyse in detail. Using an opacity level that is 100 times lower than that of ISM dust, our Jupiter-mass protoplanet features an envelope that is sufficiently cold for a CPD to form, and a free-fall region separating the CPD and the circumstellar disc emerges. Interestingly, this free-fall region appears to be a result of supersonic erosion of outer envelope material, as opposed to the static structure formation that one would expect at low opacities. Our analysis reveals that the planetary spiral arms seem to pose a significant pressure barrier that needs to be overcome through radiative cooling in order for gas to free-fall onto the CPD. The circulation inside the CPD is near-keplerian and modified by the presence of CPD spiral arms. For high opacities we recover results from the literature, finding an essentially featureless hot envelope. With this work, we demonstrate the first simulation and analysis of a complete detachment process of a protoplanet from its parent disc in a 3D radiation hydrodynamics setting.

The full spatial structure and temporal evolution of the accretion flow into the envelopes of growing gas giants in their nascent discs is only accessible in simulations. Such simulations are constrained in their approach of computing the formation of gas giants by dimensionality, resolution, consideration of self-gravity, energy treatment and the adopted opacity law. Our study explores how a number of these parameters affect that measured accretion rate of a Saturn-mass planet. We present a global 3D radiative hydrodynamics framework using the FARGOCA-code. The planet is represented by a gravitational potential with a smoothing length at the location of the planet. No mass or energy sink is used, instead luminosity and gas accretion rates are self-consistently computed. We find that the gravitational smoothing length must be resolved by at least 10 grid cells to obtain converged measurements of the gas accretion rates. Secondly, we find gas accretion rates into planetary envelopes compatible with previous studies, and continue to explain those via the structure of our planetary envelopes and their luminosities. Our measured gas accretion rates are formally in the stage of Kelvin-Helmholtz contraction due to the modest entropy loss that can be obtained over the simulation time-scale, but our accretion rates are compatible with those expected during late run-away accretion. Our detailed simulations of the gas flow into the envelope of a Saturn-mass planet provides a framework for understanding the general problem of gas accretion during planet formation and highlights circulation features that develop inside the planetary envelopes. Those circulation features feedback into the envelope energetics and can have further implications for transporting dust into the inner regions of the envelope.

Gas-giant planets, like Jupiter and Saturn, acquire massive gaseous envelopes during the approximately 3 Myr-long lifetimes of protoplanetary discs. In the core accretion scenario, the formation of a solid core of around 10 Earth masses triggers a phase of rapid gas accretion. Previous 3D grid-based hydrodynamical simulations found runaway gas accretion rates corresponding to approximately 10 to 100 Jupiter masses per Myr. Such high accretion rates would result in all planets with larger-than-10-Earth-mass cores forming Jupiter-like planets, in clear contrast to the ice giants in the Solar System and the observed exoplanet population. In this work, we use 3D hydrodynamical simulations, that include radiative transfer, to model the growth of the envelope on planets with different masses. We find that gas flows rapidly through the outer part of the envelope, but this flow does not drive accretion. Instead, gas accretion is the result of quasi-static contraction of the inner envelope, which can be orders of magnitude smaller than the mass flow through the outer atmosphere. For planets smaller than Saturn, we measure moderate gas accretion rates that are below 1 Jupiter mass per Myr. Higher mass planets, however, accrete up to 10 times faster and do not reveal a self-driven mechanism that can halt gas accretion. Therefore, the reason for the final masses of Saturn and Jupiter remains difficult to understand, unless their completion coincided with the dissipation of the Solar Nebula.

Context. Giant planets open gaps in their protoplanetary and subsequently suffer so-called type II migration. Schematically, planets are thought to be tightly locked within their surrounding disks, and forced to follow the viscous advection of gas onto the central star. This fundamental principle however has recently been questioned, as migrating planets were shown to decouple from the gas' radial drift. Aims. In this framework, we question whether the traditionally used linear scaling of migration rate of a giant planet with the disk's viscosity still holds. Additionally, we assess the role of orbit-crossing material as part of the decoupling mechanism. Methods. We have performed 2D (r, {\theta}) numerical simulations of point-mass planets embedded in locally isothermal {\alpha}-disks in steady-state accretion, with various values of {\alpha}. Arbitrary planetary accretion rates were used as a means to diminish or nullify orbit-crossing flows. Results. We confirm that the migration rate of a gap-opening planet is indeed proportional to the disk's viscosity, but is not equal to the gas drift speed in the unperturbed disk. We show that the role of gap-crossing flows is in fact negligible. Conclusions. From these observations, we propose a new paradigm for type II migration : a giant planet feels a torque from the disk that promotes its migration, while the gap profile relative to the planet is restored on a viscous timescale, thus limiting the planet migration rate to be proportional to the disk's viscosity. Hence, in disks with low viscosity in the planet region, type II migration should still be very slow. Key words. protoplanetary disks; planet-disk interactions; planets and satellites: formation

The growth of a planetary core by pebble accretion stops at the so called pebble isolation mass, when the core generates a pressure bump that traps drifting pebbles outside its orbit. If the isolation mass is very small, then gas accretion is protracted and the planet remains at a few Earth masses with a mainly solid composition. For larger values of the pebble isolation mass, the planet might be able to accrete gas from the protoplanetary disc and grow into a gas giant. Previous works have determined a scaling of the pebble isolation mass with cube of the disc aspect ratio. Here we expand on previous measurements and explore the dependency of the pebble isolation mass on all relevant parameters of the protoplanetary disc. We use 3D hydrodynamical simulations to measure the pebble isolation mass and derive a simple scaling law that captures the dependence on the local disc structure and the turbulent viscosity parameter $\alpha$. We find that small pebbles, coupled to the gas, with Stokes number $\tau_{\rm f}<0.005$ can drift through the partial gap at pebble isolation mass. However, as the planetary mass increases, particles must be decreasingly smaller to penetrate through the pressure bump. Turbulent diffusion of particles, however, can lead to an increase of the pebble isolation mass by a factor of two, depending on the strength of the background viscosity and on the pebble size. We finally explore the implications of the new scaling law of the pebble isolation mass on the formation of planetary systems by numerically integrating the growth and migration pathways of planets in evolving protoplanetary discs. Compared to models neglecting the dependence of the pebble isolation mass on the $\alpha$-viscosity, our models including this effect result in larger core masses for giant planets. These larger core masses are more akin to the core masses of the giant planets in the Solar System.

A large fraction of giant planets have gaseous envelopes that are limited to about 10 % of their total mass budget. Such planets are present in the Solar System (Uranus, Neptune) and are frequently observed in short periods around other stars (the so-called Super-Earths). In contrast to these observations, theoretical calculations based on the evolution of hydrostatic envelopes argue that such low mass envelopes cannot be maintained around cores exceeding five Earth masses. Instead, under nominal disc conditions, these planets would acquire massive envelopes through runaway gas accretion within the lifetime of the protoplanetary disc. In this work, we show that planetary envelopes are not in hydrostatic balance, which slows down envelope growth. A series of 3-dimensional, global, radiative hydrodynamical simulations reveal a steady state gas flow, which enters through the poles and exits in the disc midplane. Gas is pushed through the outer envelope in about 10 orbital timescales. In regions of the disc that are not significantly dust-depleted, envelope accretion onto cores of about five Earth masses can get stalled as the gas flow enters the deep interior. Accreted solids sublimate deep in the convective interior, but small opacity-providing grains are trapped in the flow and do not settle, which further prevents rapid envelope accretion. The transition to runaway gas accretion can however be reached when cores grow larger than typical Super-Earths, beyond 15 Earth masses, and preferably when disc opacities are below kappa=1 cm^2/g. These findings offer an explanation for the typical low-mass envelopes around the cores of Super-Earths.

Planet migration is inherently a three-dimensional (3D) problem, because Earth-size planetary cores are deeply embedded in protoplanetary disks. Simulations of these 3D disks remain challenging due to the steep requirement in resolution. Using two different hydrodynamics code, FARGO3D and PEnGUIn, we simulate disk-planet interaction for a 1 to 5 Earth-mass planet embedded in an isentropic disk. We measure the torque on the planet and ensure that the measurements are converged both in resolution and between the two codes. We find that the torque is independent of the smoothing length of the planet's potential ($r_{\rm s}$), and that it has a weak dependence on the adiabatic index of the gaseous disk ($\gamma$). The torque values correspond to an inward migration rate qualitatively similar to previous linear calculations. We perform additional simulations with explicit radiative transfer using FARGOCA, and again find agreement between 3D simulations and existing torque formulae. We also present the flow pattern around the planets, and show that active flow is present within the planet's Hill sphere, and meridional vortices are shed downstream. The vertical flow speed near the planet is faster for a smaller $r_{\rm s}$ or $\gamma$, up to supersonic speeds for the smallest $r_{\rm s}$ and $\gamma$ in our study.

We present three-dimensional simulations with nested meshes of the dynamics of the gas around a Jupiter mass planet with the JUPITER and FARGOCA codes. We implemented a radiative transfer module into the JUPITER code to account for realistic heating and cooling of the gas. We focus on the circumplanetary gas flow, determining its characteristics at very high resolution ($80\%$ of Jupiter's diameter). In our nominal simulation where the temperature evolves freely by the radiative module and reaches 13000 K at the planet, a circumplanetary envelope was formed filling the entire Roche-lobe. Because of our equation of state is simplified and probably overestimates the temperature, we also performed simulations with limited maximal temperatures in the planet region (1000 K, 1500 K, and 2000 K). In these fixed temperature cases circumplanetary disks (CPDs) were formed. This suggests that the capability to form a circumplanetary disk is not simply linked to the mass of the planet and its ability to open a gap. Instead, the gas temperature at the planet's location, which depends on its accretion history, plays also fundamental role. The CPDs in the simulations are hot and cooling very slowly, they have very steep temperature and density profiles, and are strongly sub-Keplerian. Moreover, the CPDs are fed by a strong vertical influx, which shocks on the CPD surfaces creating a hot and luminous shock-front. In contrast, the pressure supported circumplanetary envelope is characterized by internal convection and almost stalled rotation.

The terrestrial planets and the asteroids dominant in the inner asteroid belt are water poor. However, in the protoplanetary disk the temperature should have decreased below water condensation level well before the disk was photoevaporated. Thus, the global water depletion of the inner Solar System is puzling. We show that, even if the inner disk becomes cold, there cannot be direct condensation of water. This is because the snowline moves towards the Sun more slowly than the gas itself. The appearance of ice in a range of heliocentric distances swept by the snowline can only be due to the radial drift of icy particles from the outer disk. However, if a sufficiently massive planet is present, the radial drift of particles is interrupted, because the disk acquires a superKeplerian rotation just outside of the planetary orbit. From this result, we propose that the precursor of Jupiter achieved about 20 Earth masses when the snowline was still around 3 AU. This effectively fossilized the snowline at that location. Although cooling, the disk inside of the Jovian orbit remained ice-depleted because the flow of icy particles from the outer system was intercepted by the planet. This scenario predicts that planetary systems without giant planets should be much more rich in water in their inner regions than our system. We also show that the inner edge of the planetesimal disk at 0.7AU, required in terrestrial planet formation models to explain the small mass of Mercury and the absence of planets inside of its orbit, could be due to the silicate condensation line, fossilized at the end of the phase of streaming instability that generated the planetesimal seeds. Thus, when the disk cooled, silicate particles started to drift inwards of 0.7AU without being sublimated, but they could not be accreted by any pre-existing planetesimals.

For the very first time we present 3D simulations of planets embedded in stellar irradiated discs. It is well known that thermal effects could reverse the direction of planetary migration from inwards to outwards, potentially saving planets in the inner, optically thick parts of the protoplanetary disc. When considering stellar irradiation in addition to viscous friction as a source of heating, the outer disc changes from a shadowed to a flared structure. Using a suited analytical formula it has been shown that in the flared part of the disc the migration is inwards; planets can migrate outwards only in shadowed regions of the disc, { because the radial gradient of entropy is stronger there}. In order to confirm this result numerically, we have computed the total torque acting on planets held on fixed orbits embedded in stellar irradiated 3D discs using the hydrodynamical code FARGOCA. We find qualitatively good agreement between the total torque obtained with numerical simulations and the one predicted by the analytical formula. For large masses (>20 Earth masses) we find quantitative agreement, and we obtain outwards migration regions for planets up to 60 Earth masses in the early stages of accretional discs. We find nevertheless that the agreement with the analytic formula is quite fortuitous because the formula underestimates the size of the horseshoe region; this error is compensated by imperfect estimates of other terms, most likely the cooling rate and saturation.

The migration strength and direction of embedded low-mass planets depends on the disc structure. In discs with an efficient radiative transport, the migration can be directed outwards for planets with more than 3-5 Earth masses. This is due to the entropy driven corotation torque, a process that extends the lifetimes of growing planetary embryos. We investigate the influence on the disc structure caused by a jump in the alpha parameter of the viscosity to model a dead-zone structure in the disc. We focus on M-dot discs, which have a constant net mass flux. Using the resulting disc structure, we investigate the consequences for the formation of planetesimals and determine the regions of outward migration for proto-planets. We performed numerical hydrosimulations of M-dot discs in the r-z-plane. We used the explicit/implicit hydrodynamical code FARGOCA 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. Viscosity transitions inside the disc create transitions in density that stop inward migration for small planets through the so-called "planet trap" mechanism. This mechanism also works for planets down to M_P > 0.5M_E. Additionally, the viscosity transitions change the pressure gradient in the disc, which facilitates planetesimal formation via the streaming instability. However, a very steep transition in viscosity is needed to achieve in a pressure bump in the disc. The transition in viscosity facilitates planetesimal formation and can stop the migration of small-mass planets (M_P>0.5M_E), but still does not halt inward migration of smaller planets and planetesimals that are affected by gas drag. A very steep, probably unrealistic viscosity gradient is needed to trap planets of smaller masses and halt gas-drag-driven planetesimal migration at a pressure bump.

In this paper, we address the migration of small mass planets in 3D radiative disks. Indeed, migration of small planets is known to be too fast inwards in locally isothermal conditions. However, thermal effects could reverse its direction, potentially saving planets in the inner, optically thick parts of the protoplanetary disc. This effect has been seen for masses larger than 5 Earth masses, but the minimum mass for this to happen has never been probed numerically, although it is of crucial importance for planet formation scenarios. We have extended the hydro-dynamical code FARGO to 3D, with thermal diffusion. With this code, we perform simulations of embedded planets down to 2 Earth masses. For a set of discs parameters for which outward migration has been shown in the range of $[5, 35]$ Earth masses, we find that the transition to inward migration occurs for masses in the range $[3, 5]$ Earth masses. The transition appears to be due to an unexpected phenomenon: the formation of an asymmetric cold and dense finger of gas driven by circulation and libration streamlines. We recover this phenomenon in 2D simulations where we control the cooling effects of the gas through a simple modeling of the energy equation.

The strength and direction of migration of embedded low mass planets depends on the disc's structure. It has been shown that, in discs with viscous heating and radiative transport, the migration can be directed outwards. In this paper we investigate the influence of a constant dM/dt-flux through the disc, as well as the influence of the disc's metallicity on the disc's thermodynamics. We focus on dM/dt discs, which have a net mass flux through them. Utilizing the resulting disc structure, we determine the regions of outward migration in the disc. We perform numerical hydrosimulations of dM/dt discs with viscous heating, radiative cooling and stellar irradiation in 2D in the r-z-plane. We use the explicit/implicit hydrodynamical code FARGOCA 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 of embedded planets is studied by using torque formulae. For a disc of gas surface density and viscosity, we find that the discs thermal structure depends on the product Sigma_G \nu and the amount of heavy elements, while the migration of planets additionally to the mentioned quantities, depends on the amount of viscosity itself. As a result of this, the disc structure can not be approximated by simple power laws. During the lifetime of the disc, the structure of the disc changes significantly in a non-linear way in the inner parts. In the late stages of the discs evolution, outward migration is only possible if the metallicity of the disc is high. For low metallicity, planets would migrate inwards and could potentially be lost. The presented disc structures and migration maps have important consequences on the formation of planets, as they can give hints on the different formation mechanisms for different types of planets as a function of metallicity.

We examine the gas circulation near a gap opened by a giant planet in a protoplanetary disk. We show with high resolution 3D simulations that the gas flows into the gap at high altitude over the mid-plane, at a rate dependent on viscosity. We explain this observation with a simple conceptual model. From this model we derive an estimate of the amount of gas flowing into a gap opened by a planet with Hill radius comparable to the scale-height of a layered disk (i. e. a disk with viscous upper layer and inviscid midplane). Our estimate agrees with modern MRI simulations(Gressel et al., 2013). We conclude that gap opening in a layered disk can not slow down significantly the runaway gas accretion of Saturn to Jupiter-mass planets.

In the Solar System, planets have a small inclination with respect to the equatorial plane of the Sun, but there is evidence that in extrasolar systems the inclination can be very high. This spin-orbit misalignment is unexpected, as planets form in a protoplanetary disc supposedly aligned with the stellar spin. Planet-planet interactions are supposed to lead to a mutual inclination, but the effects of the protoplanetary disc are still unknown. We investigate therefore planet-disc interactions for planets above 1M_Jup. We check the influence of the inclination i, eccentricity e, and mass M_p of the planet. We perform 3D numerical simulations of protoplanetary discs with embedded high-mass planets. We provide damping formulae for i and e as a function of i, e, and M_p that fit the numerical data. For highly inclined massive planets, the gap opening is reduced, and the damping of i occurs on time-scales of the order of 10^-4 deg/yr M_disc/(0.01 M_star) with the damping of e on a smaller time-scale. While the inclination of low planetary masses (<5M_Jup) is always damped, large planetary masses with large i can undergo a Kozai-cycle with the disc. These Kozai-cycles are damped in time. Eccentricity is generally damped, except for very massive planets (M_p = 5M_Jup) where eccentricity can increase for low inclinations. The dynamics tends to a final state: planets end up in midplane and can then, over time, increase their eccentricity as a result of interactions with the disc. The interactions with the disc lead to damping of i and e after a scattering event of high-mass planets. If i is sufficiently reduced, the eccentricity can be pumped up because of interactions with the disc. If the planet is scattered to high inclination, it can undergo a Kozai-cycle with the disc that makes it hard to predict the exact movement of the planet and its orbital parameters at the dispersal of the disc.

The observed wide eccentricity distribution of extrasolar giant planets is thought to be the result of dynamical instabilities and gravitational scattering among planets. Previously, it has been assumed that the orbits in giant planet systems become gravitationally unstable after the gas nebula dispersal. It was not well understood, however, how these unstable conditions were established in the first place. In this work we numerically simulate the evolution of systems of three planets as the planets sequentially grow to Jupiter's mass, and dynamically interact among themselves and with the gas disk. We use the hydro-dynamical code FARGO that we modified by implementing the $N$-body integrator SyMBA. The new code can handle close encounters and collisions between planets. To test their stability, the planetary systems were followed with SyMBA for up to $10^8$ yr after the gas disk dispersal. We find that dynamics of the growing planets is complex, because migration and resonances raise their orbital eccentricities, and cause dynamical instabilities when gas is still around. If the dynamical instabilities occur early, planets can be removed by collisions and ejections, and the system rearranges into a new, more stable configuration. In this case, the planetary systems emerging from the gas disks are expected to be stable, and would need to be destabilized by other means (low-mass planets, planetesimal disks, etc.). Alternatively, for the giant planet system to be intrinsically unstable upon the gas disk dispersal, a special timing would be required with the growth of (at least some of) the giant planets having to occur near the end of the gas disk lifetime.

An accelerated logarithmic potential models the mean motion of stars in a flat rotation curve galaxy that sustains a wind system. For stars outside the galactic wind launching region, the asymmetric removal of linear momentum by the wind is seen as a perturbing acceleration superimposed onto the galactic potential. We study the integrability of stellar motion in an accelerated logarithmic potential. We use surfaces of section of the dynamical system to probe the integrability of motion. We provide numerical evidence that motion in an accelerated logarithmic potential is non-integrable. Large scale chaotic diffusion occurs for lower values of the projected angular momentum along the direction of acceleration and persists at all values of the angular momentum in the outer part of the galaxy inside the truncation radius where the galactic acceleration balances the wind-induced acceleration.

We propose a morphological multi-scale analysis of large-scale structures obtained by computer simulations and by observations. Structures are obtained at different scales by applying a wavelet transform on the observed and simulated data. Application of a segmentation algorithm allows a quantitative morphological description of the structures at each scale. The morphological parameter which we propose represents the deviation of a structure from sphericity. The dependence on scale of this simple parameter is shown to characterize in simulations the underlying cosmological model. We compare the HDM, CDM and mixed models with the CfA catalogue. This comparison favours a mixed model containing $65\%$ of CDM, $30\%$ of HDM, and $5\%$ baryonic matter.