Influence of the water content in protoplanetary discs on planet migration and formation

The temperature and density profiles of protoplanetary discs depend crucially on the mass fraction of micrometre-sized dust grains and on their chemical composition. A larger abundance of micrometre-sized grains leads to an overall heating of the disc, so that the water ice line moves further away from the star. An increase in the water fraction inside the disc, maintaining a fixed dust abundance, increases the temperature in the icy regions of the disc and lowers the temperature in the inner regions. Discs with a larger silicate fraction have the opposite effect. Here we explore the consequence of the dust composition and abundance for the formation and migration of planets. We find that discs with low water content can only sustain outwards migration for planets up to 4 Earth masses, while outwards migration in discs with a larger water content persists up to 8 Earth masses in the late stages of the disc evolution. Icy planetary cores that do not reach run-away gas accretion can thus migrate to orbits close to the host star if the water abundance is low. Our results imply that hot and warm super-Earths found in exoplanet surveys could have formed beyond the ice line and thus contain a significant fraction in water. These water-rich super-Earths should orbit primarily around stars with a low oxygen abundance, where a low oxygen abundance is caused by either a low water-to-silicate ratio or by overall low metallicity.Comment: 15 pages, 12 figures, accepted by A&

Planet population synthesis driven by pebble accretion in cluster environments

The evolution of protoplanetary discs embedded in stellar clusters depends on the age and the stellar density in which they are embedded. Stellar clusters of young age and high stellar surface density destroy protoplanetary discs by external photoevaporation and stellar encounters. Here we consider the effect of background heating from newly formed stellar clusters on the structure of protoplanetary discs and how it affects the formation of planets in these discs. Our planet formation model is build on the core accretion scenario including pebble accretion. We synthesize planet populations that we compare to observations. The giant planets in our simulations migrate over large distances due to the fast type-II migration regime induced by a high disc viscosity ($\alpha=5.4 \times 10^{-3}$). Cold Jupiters (r>1 AU) originate preferably from the outer disc, while hot Jupiters (r<0.1 AU) preferably form in the inner disc. We find that the formation of gas giants via pebble accretion is in agreement with the metallicity correlation, meaning that more gas giants are formed at larger metallicity. However, our synthetic population of isolated stars host a significant amount of giant planets even at low metallicity, in contradiction to observations where giant planets are preferably found around high metallicity stars, indicating that pebble accretion is very efficient in the standard pebble accretion framework. On the other hand, discs around stars embedded in cluster environments hardly form any giant planets at low metallicity in agreement with observations, where these changes originate from the increased temperature in the outer parts of the disc, which prolongs the core accretion time-scale of the planet. We therefore conclude that the outer disc structure and the planet's formation location determines the giant planet occurrence rate and the formation efficiency of cold and hot Jupiters.Comment: 12 pages, accepted for publication in MNRA

The growth of planets by pebble accretion in evolving protoplanetary discs

The formation of planets depends on the underlying protoplanetary disc structure, which influences both the accretion and migration rates of embedded planets. The disc itself evolves on time-scales of several Myr during which both temperature and density profiles change as matter accretes onto the central star. Here we use a detailed model of an evolving disc to determine the growth of planets by pebble accretion and their migration through the disc. Cores that reach their pebble isolation mass accrete gas to finally form giant planets with extensive gas envelopes, while planets that do not reach pebble isolation mass are stranded as ice giants and ice planets containing only minor amounts of gas in their envelopes. Unlike earlier population synthesis models, our model works without any artificial reductions in migration speed and for protoplanetary discs with gas and dust column densities similar to those inferred from observations. We find that in our nominal disc model the emergence of planetary embryos preferably occurs after approximately 2 Myr in order to not exclusively form gas giants, but also ice giants and smaller planets. The high pebble accretion rates ensure that critical core masses for gas accretion can be reached at all orbital distances. Gas giant planets nevertheless experience significant reduction in semi-major axes by migration. Considering instead planetesimal accretion for planetary growth, we show that formation time-scales are too long to compete with the migration time-scales and the dissipation time of the protoplanetary disc. Altogether, we find that pebble accretion overcomes many of the challenges in the formation of ice and gas giants in evolving protoplanetary discs.Comment: Accepted by A&A, now with language editin

Planet migration in three-dimensional radiative discs

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.Comment: 16 pages, Astronomy&Astrophysics, in pres

Influence of grain growth on the thermal structure of protoplanetary discs

The thermal structure of a protoplanetary disc is regulated by the opacity that dust grains provide. However, previous works have often considered simplified prescriptions for the dust opacity in hydrodynamical disc simulations, e.g. by considering only a single particle size. In the present work we perform 2D hydrodynamical simulations of protoplanetary discs where the opacity is self-consistently calculated for the dust population, taking into account the particle size, composition and abundance. We first compare simulations using single grain sizes to two different multi-grain size distributions at different levels of turbulence strengths, parameterized through the $\alpha$-viscosity, and different gas surface densities. Assuming a single dust size leads to inaccurate calculations of the thermal structure of discs, because the grain size dominating the opacity increases with orbital radius. Overall the two grain size distributions, one limited by fragmentation only and the other determined from a more complete fragmentation-coagulation equilibrium, give similar results for the thermal structure. We find that both grain size distributions give less steep opacity gradients that result in less steep aspect ratio gradients, in comparison to discs with only micrometer sized dust. Moreover, in the discs with a grain size distribution, the innermost outward migration region is removed and planets embedded is such discs experience lower migration rates. We also investigate the dependency of the water iceline position on the alpha-viscosity, the initial gas surface density at 1 AU and the dust-to-gas ratio and find $r_{ice} \propto \alpha^{0.61} \Sigma_{g,0}^{0.8} f_{DG}^{0.37}$ independently of the distribution used. The inclusion of the feedback loop between grain growth, opacities and disc thermodynamics allows for more self-consistent simulations of accretion discs and planet formation.Comment: Accepted by A&A, 27 pages, 19 figure

The structure of protoplanetary discs around evolving young stars

The formation of planets with gaseous envelopes takes place in protoplanetary accretion discs on time-scales of several millions of years. Small dust particles stick to each other to form pebbles, pebbles concentrate in the turbulent flow to form planetesimals and planetary embryos and grow to planets, which undergo substantial radial migration. All these processes are influenced by the underlying structure of the protoplanetary disc, specifically the profiles of temperature, gas scale height and density. The commonly used disc structure of the Minimum Mass Solar Nebular (MMSN) is a simple power law in all these quantities. However, protoplanetary disc models with both viscous and stellar heating show several bumps and dips in temperature, scale height and density caused by transitions in opacity, which are missing in the MMSN model. These play an important role in the formation of planets, as they can act as sweet spots for the formation of planetesimals via the streaming instability and affect the direction and magnitude of type-I-migration. We present 2D simulations of accretion discs that feature radiative cooling, viscous and stellar heating, and are linked to the observed evolutionary stages of protoplanetary discs and their host stars. These models allow us to identify preferred planetesimal and planet formation regions in the protoplanetary disc as a function of the disc's metallicity, accretion rate and lifetime. We derive simple fitting formulae that feature all structural characteristics of protoplanetary discs during the evolution of several Myr. These fits are straightforward to apply for modelling any growth stage of planets where detailed knowledge of the underlying disc structure is required.Comment: Accepted by A&A, v3 corrected small typo in the fitting formula

Stellar abundance of binary stars: their role in determining the formation location of super-Earths and ice giants

Binary stars form from the same parent molecular cloud and thus have the same chemical composition. Forming planets take building material (solids) away from the surrounding protoplanetary disc. Assuming that the disc's accretion onto the star is the main process that clears the disc, the atmosphere of the star will show abundance reductions caused by the material accreted by the forming planet(s). If planets are only forming around one star of a binary system, the planet formation process can result in abundance differences in wide binary stars, if their natal protoplanetary discs do not interact during planet formation. Abundance differences in the atmospheres of wide binaries hosting giant planets have already been observed and linked to the formation location of giant planets. Here, we model how much building material is taken away for super-Earth planets that form inside/outside of the water ice line as well as ice giants forming inside/outside of the CO ice line. Our model predicts a significant abundance difference $\Delta$[X/H] in the stellar atmospheres of the planet-hosting binary component. Our model predicts that super-Earths that form inside the water ice line ($r<r_{\rm H_2O}$) will result in an $\Delta$[Fe/H]/$\Delta$[O/H] abundance difference in the their host star that is a factor of 2 larger than for super-Earths formed outside the water ice line ($r>r_{\rm H_2O}$) in the water rich parts of the disc. Additionally, our model shows that the $\Delta$[Fe/H]/$\Delta$[C/H] abundance difference in the host star is at least a factor of 3 larger for ice giants formed at $r<r_{\rm CO}$ compared to ice giants formed far out in the protoplanetary disc ($r>r_{\rm CO}$). Future observations of wide binary star systems hosting super-Earths and ice giants could therefore help to constrain the migration pathway of these planets and thus constrain planet formation theories.Comment: accepted by MNRA

Enriching inner discs and giant planets with heavy elements

Giant exoplanets seem to have on average a much larger heavy element content than the solar system giants. Past attempts to explain these heavy element contents include collisions between planets, accretion of volatile rich gas and accretion of gas enriched in micro-metre sized solids. However, these different theories individually could not explain the heavy element content of giants and the volatile to refractory ratios in atmospheres of giant planets at the same time. Here we combine the approaches of gas accretion enhanced with vapor and small micro-meter sized dust grains. As pebbles drift inwards, the volatile component evaporates and enriches the disc, while the smaller silicate core of the pebble continues to move inwards. The smaller silicate pebbles drift slower, leading to a pile-up of material interior to the water ice line, increasing the dust-to-gas ratio interior to the ice line. Under the assumption that these small dust grains follow the motion of the gas, gas accreting giants accrete large fractions of small solids in addition to the volatile vapor. The effectiveness of the solid enrichment requires a large disc radius to maintain the pebble flux for a long time and a large viscosity that reduces the size and inward drift of the small dust grains. However, this process depends crucially on the debated size difference of the pebbles interior and exterior of the water ice line. On the other hand, the volatile component released by the inward drifting pebbles can lead to a large enrichment with heavy element vapor, independently of a size difference of pebbles interior and exterior to the water ice line. Our model stresses the importance of the disc's radius and viscosity on the enrichment of dust and vapor. Consequently we show how our model could explain the heavy element content of the majority of giant planets by using combined estimates of dust and vapor enrichment.Comment: Accepted by A&A, 10 pages, 9 figure

Highly inclined and eccentric massive planets I: Planet-disc interactions

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.Comment: accepted for publication in Astronomy and Astrophysic