Source author record

Bertram Bitsch

Bertram Bitsch appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2022arXiv

Emerging population of gap-opening planets around type-A stars -- Long-term evolution of the forming planets around HD 163296

The presence of forming planets embedded in their protoplanetary disks has been inferred from the detection of multiring structures in such disks. Most of these suspected planets are undetectable by direct imaging observations at current measurement sensitivities. Inward migration and accretion might make these putative planets accessible to the Doppler method, but the actual extent of growth and orbital evolution remains unconstrained. Under the premise that the gaps in the disk around HD 163296 originate from new-born planets, we investigate if and under which circumstances the gap-opening planets could represent progenitors of the exoplanet population detected around A-type stars. In particular, we study the dependence of final planetary masses and orbital parameters on the viscosity of the disk. The evolution of the embedded planets was simulated throughout the disk lifetime and up to 100 Myr after the dispersal of the disk, taking the evolving disk structure and a likely range of disk lifetimes into account. We find that the final configuration of the planets is largely determined by the $α$ viscosity parameter of the disk and less dependent on the choice for the disk lifetime and the initial planetary parameters. If we assume that planets such as those in HD 163296 evolve to form the observed exoplanet population of A-type stars, a $α$ parameter on the order of $3.16 \times 10^{-4} \lesssim α\lesssim 10^{-3}$ is required for the disks to induce sufficiently high migration rates. Depending on whether or not future direct imaging surveys will uncover a larger number of planets with $m_\mathrm{pl} \lesssim 3 M_\mathrm{Jup}$ and $a_\mathrm{pl} \gtrsim 10 \mathrm{AU}$ we expect the $α$ parameter to be at the lower or upper end of this range, always under the assumption that such disks indeed harbor wide orbit planets.

preprint2022arXiv

Interpreting the atmospheric composition of exoplanets: sensitivity to planet formation assumptions

Constraining planet formation based on the atmospheric composition of exoplanets is a fundamental goal of the exoplanet community. Existing studies commonly try to constrain atmospheric abundances, or to analyze what abundance patterns a given description of planet formation predicts. However, there is also a pressing need to develop methodologies that investigate how to transform atmospheric compositions into planetary formation inferences. In this study we summarize the complexities and uncertainties of state-of-the-art planet formation models and how they influence planetary atmospheric compositions. We introduce a methodology that explores the effect of different formation model assumptions when interpreting atmospheric compositions. We apply this framework to the directly imaged planet HR 8799e. Based on its atmospheric composition, this planet may have migrated significantly during its formation. We show that including the chemical evolution of the protoplanetary disk leads to a reduced need for migration. Moreover, we find that pebble accretion can reproduce the planet's composition, but some of our tested setups lead to too low atmospheric metallicities, even when considering that evaporating pebbles may enrich the disk gas. We conclude that the definitive inversion from atmospheric abundances to planet formation for a given planet may be challenging, but a qualitative understanding of the effects of different formation models is possible, opening up pathways for new investigations.

preprint2022arXiv

Shallowness of circulation in hot Jupiters -- Advancing the Ohmic dissipation model

The inflated radii of giant short-period extrasolar planets collectively indicate that the interiors of hot Jupiters are heated by some anomalous energy dissipation mechanism. Although a variety of physical processes have been proposed to explain this heating, recent statistical evidence points to the confirmation of explicit predictions of the Ohmic dissipation theory, elevating this mechanism as the most promising candidate for resolving the radius inflation problem. In this work, we present an analytic model for the dissipation rate and derive a simple scaling law that links the magnitude of energy dissipation to the thickness of the atmospheric weather layer. From this relation, we find that the penetration depth influences the Ohmic dissipation rate by an order of magnitude. We further investigate the weather layer depth of hot Jupiters from the extent of their inflation and show that, depending on the magnetic field strength, hot Jupiter radii can be maintained even if the circulation layer is relatively shallow. Additionally, we explore the evolution of zonal wind velocities with equilibrium temperature by matching our analytic model to statistically expected dissipation rates. From this analysis, we deduce that the wind speed scales approximately as $1/\sqrt{T_\mathrm{eq}-T_0}$, where $T_0$ is a constant that equals $T_0 \sim 1000~\mathrm{K}-1800~\mathrm{K}$ depending on planet-specific parameters (radius, mass, etc.). This work outlines inter-related constraints on the atmospheric flow and the magnetic field of hot Jupiters and provides a foundation for future work on the Ohmic heating mechanism.

preprint2021arXiv

How the planetary eccentricity influences the pebble isolation mass

We investigate the pebble isolation mass for a planet on a fixed eccentric orbit in its protoplanetary disc by conducting a set of 2D hydrodynamical simulations including dust turbulent diffusion. A range of planet eccentricities up to $e=0.2$ is adopted. Our simulations also cover a range of $α-$turbulent viscosities, and for each pair $\{α,e\}$ the pebble isolation mass is estimated as the minimum planet mass in our simulations such that solids with a Stokes number $\gtrsim 0.05$ do not flow across the planet orbit and remain trapped around a pressure bump outside the planet gap. For $α<10^{-3}$, we find that eccentric planets reach a well-defined pebble isolation mass, which can be smaller than for planets on circular orbits when the eccentricity remains smaller than the disc's aspect ratio. We provide a fitting formula for how the pebble isolation mass depends on planet eccentricity. However, for $α> 10^{-3}$, eccentric planets cannot fully stall the pebbles flow, and thus do not reach a well-defined pebble isolation mass. Our results suggest that the maximum mass reached by rocky cores should exhibit a dichotomy depending on the disc turbulent viscosity. While being limited to ${\cal O}(10\,M_\oplus)$ in low-viscosity discs, this maximum mass could reach much larger values in discs with a high turbulent viscosity in the planet vicinity. Our results further highlight that pebble filtering by growing planets might not be as effective as previously thought, especially in high-viscosity discs, with important implications to protoplanetary discs observations.

preprint2021arXiv

Planetesimal rings as the cause of the Solar System's planetary architecture

Astronomical observations reveal that protoplanetary disks around young stars commonly have ring- and gap-like structures in their dust distributions. These features are associated with pressure bumps trapping dust particles at specific locations, which simulations show are ideal sites for planetesimal formation. Here we show that our Solar System may have formed from rings of planetesimals -- created by pressure bumps -- rather than a continuous disk. We model the gaseous disk phase assuming the existence of pressure bumps near the silicate sublimation line (at $T \sim$1400~K), water snowline (at $T \sim$170~K), and CO-snowline (at $T \sim$30~K). Our simulations show that dust piles up at the bumps and forms up to three rings of planetesimals: a narrow ring near 1~au, a wide ring between $\sim$3-4~au and $\sim$10-20~au, and a distant ring between $\sim$20~au and $\sim$45~au. We use a series of simulations to follow the evolution of the innermost ring and show how it can explain the orbital structure of the inner Solar System and provides a framework to explain the origins of isotopic signatures of Earth, Mars and different classes of meteorites. The central ring contains enough mass to explain the rapid growth of the giant planets' cores. The outermost ring is consistent with dynamical models of Solar System evolution proposing that the early Solar System had a primordial planetesimal disk beyond the current orbit of Uranus.

preprint2020arXiv

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 $α$-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 α^{0.61} Σ_{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.

preprint2020arXiv

The origins of nearly coplanar, non-resonant systems of close-in super-Earths

Some systems of close-in "super-Earths" contain five or more planets on non-resonant but compact and nearly coplanar orbits. The Kepler-11 system is an iconic representative of this class of system. It is challenging to explain their origins given that planet-disk interactions are thought to be essential to maintain such a high degree of coplanarity, yet these same interactions invariably cause planets to migrate into chains of mean motion resonances. Here we mine a large dataset of dynamical simulations of super-Earth formation by migration. These simulations match the observed period ratio distribution as long as the vast majority of planet pairs in resonance become dynamically unstable. When instabilities take place resonances are broken during a late phase of giant impacts, and typical surviving systems have planet pairs with significant mutual orbital inclinations. However, a subset of our unstable simulations matches the Kepler-11 system in terms of coplanarity, compactness, planet-multiplicity and non-resonant state. This subset have dynamical instability phases typically much shorter than ordinary systems. Unstable systems may keep a high degree of coplanarity post-instability if planets collide at very low orbital inclinations ($\lesssim1^\circ$) or if collisions promote efficient damping of orbital inclinations. If planetary scattering during the instability takes place at low orbital inclinations ($\text{i}\lesssim1^\circ$), orbital inclinations are barely increased by encounters before planets collide.When planetary scattering pumps orbital inclinations to higher values ($\gtrsim 1^\circ$) planets tend to collide at higher mutual orbital inclinations, but depending on the geometry of collisions mergers' orbital inclinations may be efficiently damped. Each of these formation pathways can produce analogues to the Kepler-11 system.

preprint2016arXiv

Did Jupiter's core form in the innermost parts of the Sun's protoplanetary disk?

Jupiter's core is generally assumed to have formed beyond the snow line. Here we consider an alternative scenario, that Jupiter's core may have accumulated in the innermost parts of the protoplanetary disk. A growing body of research suggests that small particles ("pebbles") continually drift inward through the disk. If a fraction of drifting pebbles is trapped at the inner edge of the disk a several Earth-mass core can quickly grow. Subsequently, the core may migrate outward beyond the snow line via planet-disk interactions. Of course, to reach the outer Solar System Jupiter's core must traverse the terrestrial planet-forming region. We use N-body simulations including synthetic forces from an underlying gaseous disk to study how the outward migration of Jupiter's core sculpts the terrestrial zone. If the outward migration is fast (Tmig~10^4 years), the core simply migrates past resident planetesimals and planetary embryos. However, if its migration is slower (Tmig~10^5 years) the core removes solids from the inner disk by shepherding objects in mean motion resonances. In many cases the disk interior to 0.5-1 AU is cleared of embryos and most planetesimals. By generating a mass deficit close to the Sun, the outward migration of Jupiter's core may thus explain the absence of terrestrial planets closer than Mercury. Jupiter's migrating core often stimulates the growth of another large (~Earth-mass) core -- that may provide a seed for Saturn's core -- trapped in exterior resonance. The migrating core also may transport a fraction of terrestrial planetesimals, such as the putative parent bodies of iron meteorites, to the asteroid belt.

preprint2016arXiv

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.

preprint2016arXiv

Planet Formation Imager (PFI): science vision and key requirements

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.

preprint2016arXiv

Runaway gas accretion and gap opening versus type~I migration

Growing planets interact with their natal protoplanetary disc, which exerts a torque onto them allowing them to migrate in the disc. Small mass planets do not affect the gas profile and migrate in the fast type~I migration. Although type~I migration can be directed outwards for planets smaller than $20-30 M_\oplus$, planets above this mass should be lost into the central star long before the disc disperses. Massive planets open a gap and subsequently migrate in the slower, type~II migration, which could save them from migrating all the way to the star. Hence, growing giant planets can be saved if and only if they can reach the gap opening mass, because this extends their migration time-scale, allowing them to eventually survive at large orbits until the disc itself disperses. However, most of the previous studies only measured the torques on planets with fixed masses and orbits to determine the migration rate. Additionally, the transition between type~I and type~II migration itself is not well studied. Here we use isothermal 2D disc simulations with FARGO-2D1D to study the migration behaviour of gas accreting protoplanets in discs. We find that migrating giant planets always open gaps in the disc. We further show analytically and numerically that planets growing in the runaway gas accretion regime will reach gap opening masses before migrating all the way to the central star if the disc is not extremely viscous and/or thick. An accretion rate limited to the radial gas flow in the disc, in contrast, is not fast enough. When gas accretion by the planet is taken into account, the gap opening process is accelerated because the planet accretes material originating from its horseshoe region. This allows an accreting planet to transition to type~II migration before being lost even if gas fails to be provided for a rapid enough growth and the gap opening mass is not reached.

preprint2014arXiv

Earth-like Habitats in Planetary Systems

Understanding the concept of habitability is related to an evolutionary knowledge of the particular planet-in-question. Additional indications so-called "systemic aspects" of the planetary system as a whole governs a particular planet's claim on habitability. Here we focus on such systemic aspects and discuss their relevance to the formation of an 'Earth-like' habitable planet. We summarize our results obtained by lunar sample work and numerical models within the framework of the Research Alliance "Planetary Evolution and Life". We consider various scenarios which simulate the dynamical evolution of the Solar System and discuss the likelihood of forming an Earth-like world orbiting another star. Our model approach is constrained by observations of the modern Solar System and the knowledge of its history. Results suggest that the long-term presence of terrestrial planets is jeopardized due to gravitational interactions if giant planets are present. But habitability of inner rocky planets may be supported in those planetary systems hosting giant planets. Gravitational interactions within a complex multiple-body structure including giant planets may supply terrestrial planets with materials which formed in the colder region of the proto-planetary disk. During these processes, water, the prime requisite for habitability, is delivered to the inner system. This may occur either during the main accretion phase of terrestrial planets or via impacts during a post-accretion bombardment. Results for both processes are summarized and discussed with reference to the lunar crater record. Starting from a scenario involving migration of the giant planets this contribution discusses the delivery of water to Earth, the modification of atmospheres by impacts in a planetary system context and the likelihood of the existence of extrasolar Earth-like habitable worlds.

preprint2014arXiv

Stellar irradiated discs and implications on migration of embedded planets II: accreting-discs

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 ν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.

preprint2014arXiv

Stellar irradiated discs and implications on migration of embedded planets III: viscosity transitions

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.

preprint2014arXiv

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.

preprint2013arXiv

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.

preprint2013arXiv

Planet-Disc Interactions and Early Evolution of Planetary Systems

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.

preprint2012arXiv

Stellar irradiated discs and implications on migration of embedded planets I: equilibrium discs

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.

preprint2011arXiv

Evolution of inclined planets in three-dimensional radiative discs

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.

preprint2011arXiv

Range of outward migration and influence of the disc's mass on the migration of giant planet cores

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.

preprint2010arXiv

Orbital evolution of eccentric planets in radiative discs

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.

preprint2009arXiv

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.