Showing 1–50 of 67 results for author: Masset, F

A thermodynamic criterion for the formation of Circumplanetary Disks

Authors: Leonardo Krapp, Kaitlin M. Kratter, Andrew N. Youdin, Pablo Benítez-Llambay, Frédéric Masset, Philip J. Armitage

Abstract: The formation of circumplanetary disks is central to our understanding of giant planet formation, influencing their growth rate during the post-runaway phase and observability while embedded in protoplanetary disks. We use 3D global multifluid radiation hydrodynamics simulations with the FARGO3D code to define the thermodynamic conditions that enable circumplanetary disk formation around Jovian pl… ▽ More The formation of circumplanetary disks is central to our understanding of giant planet formation, influencing their growth rate during the post-runaway phase and observability while embedded in protoplanetary disks. We use 3D global multifluid radiation hydrodynamics simulations with the FARGO3D code to define the thermodynamic conditions that enable circumplanetary disk formation around Jovian planets on wide orbits. Our simulations include stellar irradiation, viscous heating, static mesh refinement, and active calculation of opacity based on evolving dust fluids. We find a necessary condition for the formation of circumplanetary disks in terms of a mean cooling time: when the cooling time is at least one order of magnitude shorter than the orbital time scale, the specific angular momentum of the gas is nearly Keplerian at scales of $R_{\rm{Hill}}/3$. We show that the inclusion of multifluid dust dynamics favors rotational support because dust settling produces an anisotropic opacity distribution that favors rapid cooling. In all our models with radiation hydrodynamics, specific angular momentum decreases as time evolves in agreement with the formation of an inner isentropic envelope due to compressional heating. The isentropic envelope can extend up to $R_{\rm{Hill}}/3$ and shows negligible rotational support. Thus, our results imply that young gas giant planets may host spherical isentropic envelopes, rather than circumplanetary disks. △ Less

Submitted 26 February, 2024; v1 submitted 22 February, 2024; originally announced February 2024.

Comments: 17 pages, 10 figures, submitted to ApJ. Comments are welcome

On the interaction of pebble accreting embryos with the gaseous disc: importance of thermal forces

Authors: S. Cornejo, F. S. Masset, F. J. Sánchez-Salcedo

Abstract: A planetary embryo embedded in a gaseous disc can grow by pebble accretion while subjected to a gravitational force from the disc that changes its orbital elements. Usually, that force is considered to arise from the Lindblad and corotation resonances with the embryo. However, more important contributions exist for low-mass planets. Radiative thermal diffusion in the vicinity of embryos yields an… ▽ More A planetary embryo embedded in a gaseous disc can grow by pebble accretion while subjected to a gravitational force from the disc that changes its orbital elements. Usually, that force is considered to arise from the Lindblad and corotation resonances with the embryo. However, more important contributions exist for low-mass planets. Radiative thermal diffusion in the vicinity of embryos yields an additional contribution to the disc's force that damps the eccentricity and inclination much more vigorously than the resonant interaction with the disc, and that in general induces fast inward migration. In addition, the irradiation of the disc by a hot embryo gives rise to an additional contribution that excites eccentricity and inclination, and induces outward migration. Which of the two contributions dominates depends on the embryo's luminosity. We assess the importance of these contributions (termed thermal forces) on the dynamics and growth of a set of pebble-accreting embryos initially of Martian mass, by means of N-body simulations that include analytic expressions for the disc's force. We find very different outcomes for the embryos subjected to thermal forces and those subjected only to resonant forces. Importantly, we find that the median final mass of the embryos subjected to thermal forces is nearly independent of the metallicity, whereas this mass roughly scales with the metallicity when they are subjected only to resonant forces. These results can be explained by the strong damping of eccentricity and inclination at low metallicity, which enhances the embryos' accretion efficiency. △ Less

Submitted 15 May, 2023; originally announced May 2023.

Comments: Accepted for publication in MNRAS

Evolution of the eccentricity and inclination of low-mass planets subjected to thermal forces: a numerical study

Authors: S. Cornejo, F. S. Masset, R. O. Chametla, S. Fromenteau

Abstract: By means of three dimensional, high resolution hydrodynamical simulations we study the orbital evolution of weakly eccentric or inclined low-mass protoplanets embedded in gaseous discs subject to thermal diffusion. We consider both non-luminous planets, and planets that also experience the radiative feedback from their own luminosity. We compare our results to previous analytical work, and find th… ▽ More By means of three dimensional, high resolution hydrodynamical simulations we study the orbital evolution of weakly eccentric or inclined low-mass protoplanets embedded in gaseous discs subject to thermal diffusion. We consider both non-luminous planets, and planets that also experience the radiative feedback from their own luminosity. We compare our results to previous analytical work, and find that thermal forces (the contribution to the disc's force arising from thermal effects) match those predicted by linear theory within $\sim 20$%. When the planet's luminosity exceeds a threshold found to be within $10$% of that predicted by linear theory, its eccentricity and inclination grow exponentially, whereas these quantities undergo a strong damping below this threshold. In this regime of low luminosity indeed, thermal diffusion cools the surroundings of the planet and allows gas to accumulate in its vicinity. It is the dynamics of this gas excess that contributes to damp eccentricity and inclination. The damping rates obtained can be up to $h^{-1}$ times larger than those due to the resonant interaction with the disc, where $h$ is the disc's aspect ratio. This suggests that models that incorporate planet-disc interactions using well-known formulae based on resonant wave-launching to describe the evolution of eccentricity and inclination underestimate the damping action of the disc on the eccentricity and inclination of low-mass planets by an order of magnitude. △ Less

Submitted 1 March, 2023; originally announced March 2023.

Comments: Accepted for publication in MNRAS

Planet-Disk Interactions

Authors: Sijme-Jan Paardekooper, Ruobing Dong, Paul Duffell, Jeffrey Fung, Frederic S. Masset, Gordon Ogilvie, Hidekazu Tanaka

Abstract: Planet-disk interactions, where an embedded massive body interacts gravitationally with the protoplanetary disk it was formed in, can play an important role in reshaping both the disk and the orbit of the planet. Spiral density waves are launched into the disk by the planet, which, if they are strong enough, can lead to the formation of a gap. Both effects are observable with current instruments.… ▽ More Planet-disk interactions, where an embedded massive body interacts gravitationally with the protoplanetary disk it was formed in, can play an important role in reshaping both the disk and the orbit of the planet. Spiral density waves are launched into the disk by the planet, which, if they are strong enough, can lead to the formation of a gap. Both effects are observable with current instruments. The back-reaction of perturbations induced in the disk, both wave-like and non-wavelike, is a change in orbital elements of the planet. The efficiency of orbital migration is a long-standing problem in planet formation theory. We discuss recent progress in planet-disk interactions for different planet masses and disk parameters, in particular the level of turbulence, and progress in modeling observational signatures of embedded planets. △ Less

Submitted 8 November, 2022; v1 submitted 17 March, 2022; originally announced March 2022.

Comments: 28 pages, 12 figures, chapter to appear in Protostars and Planets VII, eds. Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura. Updated version including community input

How the planetary eccentricity influences the pebble isolation mass

Authors: Raúl O. Chametla, Frédéric S. Masset, Clément Baruteau, Bertram Bitsch

Abstract: 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 minim… ▽ More 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. △ Less

Submitted 28 December, 2021; originally announced December 2021.

Comments: 9 pages, 7 figures. Accepted for publication in MNRAS

Eccentricity driving of pebble accreting low-mass planets

Authors: David A. Velasco-Romero, Frédéric S. Masset, Romain Teyssier

Abstract: By means of high resolution hydrodynamical, three-dimensional calculations with nested-meshes, we evaluate the eccentricity reached by a low-mass, luminous planet embedded in an inviscid disc with constant thermal diffusivity and subjected to thermal forces. We find that a cell size of at most one tenth of the size of the region heated by the planet is required to get converged results. When the p… ▽ More By means of high resolution hydrodynamical, three-dimensional calculations with nested-meshes, we evaluate the eccentricity reached by a low-mass, luminous planet embedded in an inviscid disc with constant thermal diffusivity and subjected to thermal forces. We find that a cell size of at most one tenth of the size of the region heated by the planet is required to get converged results. When the planet's luminosity is supercritical, we find that it reaches an eccentricity of order $10^{-2}$--$10^{-1}$, which increases with the luminosity and broadly scales with the disc's aspect ratio. Restricting our study to the case of pebble accretion, we incorporate to our model the dependence of the accretion rate of pebbles on the eccentricity. There is therefore a feedback between eccentricity, which determines the accretion rate and hence the planet's luminosity, and the luminosity, which yields the eccentricity attained through thermal forces. We solve for the steady state eccentricity and study how this quantity depends on the disc's turbulence strength parameter $α_z$, on the dimensionless stopping time of the pebbles $τ_s$, on the inward mass flux of pebbles and on the headwind (the difference between the gas velocity and the Keplerian velocity). We find that in general low-mass planets (up to a few Earth masses) reach eccentricities comparable to the disc's aspect ratio, or a sizeable fraction of the latter. Eccentric, low-mass protoplanets should therefore be the norm rather than the exception, even if they orbit far from other planets or from large scale disturbances in the disc. △ Less

Submitted 22 November, 2021; originally announced November 2021.

The importance of thermal torques on the migration of planets growing by pebble accretion

Authors: O. M. Guilera, M. M. Miller Bertolami, F. Masset, J. Cuadra, J. Venturini, M. P. Ronco

Abstract: A key process in planet formation is the exchange of angular momentum between a growing planet and the protoplanetary disc, which makes the planet migrate through the disc. Several works show that in general low-mass and intermediate-mass planets migrate towards the central star, unless corotation torques become dominant. Recently, a new kind of torque, called the thermal torque, was proposed as a… ▽ More A key process in planet formation is the exchange of angular momentum between a growing planet and the protoplanetary disc, which makes the planet migrate through the disc. Several works show that in general low-mass and intermediate-mass planets migrate towards the central star, unless corotation torques become dominant. Recently, a new kind of torque, called the thermal torque, was proposed as a new source that can generate outward migration of low-mass planets. While the Lindblad and corotation torques depend mostly on the properties of the protoplanetary disc and on the planet mass, the thermal torque depends also on the luminosity of the planet, arising mainly from the accretion of solids. Thus, the accretion of solids plays an important role not only in the formation of the planet but also in its migration process. In a previous work, we evaluated the thermal torque effects on planetary growth and migration mainly in the planetesimal accretion paradigm. In this new work, we study the role of the thermal torque within the pebble accretion paradigm. Computations are carried out consistently in the framework of a global model of planet formation that includes disc evolution, dust growth and evolution, and pebble formation. We also incorporate updated prescriptions of the thermal torque derived from high resolution hydrodynamical simulations. Our simulations show that the thermal torque generates extended regions of outward migration in low viscosity discs. This has a significant impact in the formation of the planets. △ Less

Submitted 10 August, 2021; originally announced August 2021.

Comments: Accepted for publication in MNRAS

Numerical study of coorbital thermal torques on cold or hot satellites

Authors: Raúl O. Chametla, Frédéric S. Masset

Abstract: We evaluate the thermal torques exerted on low-mass planets embedded in gaseous protoplanetary discs with thermal diffusion, by means of high-resolution three-dimensional hydrodynamics simulations. We confirm that thermal torques essentially depend on the offset between the planet and its corotation, and find a good agreement with analytic estimates when this offset is small compared to the size o… ▽ More We evaluate the thermal torques exerted on low-mass planets embedded in gaseous protoplanetary discs with thermal diffusion, by means of high-resolution three-dimensional hydrodynamics simulations. We confirm that thermal torques essentially depend on the offset between the planet and its corotation, and find a good agreement with analytic estimates when this offset is small compared to the size of the thermal disturbance. For larger offsets that may be attained in discs with a large pressure gradient or a small thermal diffusivity, thermal torques tend toward an asymptotic value broadly compatible with results from a dynamical friction calculation in an unsheared medium. We perform a convergence study and find that the thermal disturbance must be resolved over typically 10 zones for a decent agreement with analytic predictions. We find that the luminosity at which the net thermal torque changes sign matches that predicted by linear theory within a few percents. Our study confirms that thermal torques usually supersede Lindblad and corotation torques by almost an order of magnitude for low mass planets. As we increase the planetary mass, we find that the ratio of thermal torques to Lindblad and corotation torques is progressively reduced, and that the thermal disturbance is increasingly distorted by the horseshoe flow. Overall, we find that thermal torques are dominant for masses up to an order of magnitude larger than implemented in recent models of planetary population synthesis. We finally briefly discuss the case of stellar or intermediate-mass objects embedded in discs around AGNs. △ Less

Submitted 24 November, 2020; originally announced November 2020.

Comments: 12 pages, 9 figures. Accepted for publication in MNRAS

Dynamical friction with radiative feedback -- II. High resolution study of the subsonic regime

Authors: David A. Velasco-Romero, Frédéric S. Masset

Abstract: Recent work has suggested that the net gravitational force acting on a massive and luminous perturber travelling through a gaseous and opaque medium can have same direction as the perturber's motion (an effect sometimes called negative dynamical friction). Analytic results were obtained using a linear analysis and were later confirmed by means of non-linear numerical simulations which did not reso… ▽ More Recent work has suggested that the net gravitational force acting on a massive and luminous perturber travelling through a gaseous and opaque medium can have same direction as the perturber's motion (an effect sometimes called negative dynamical friction). Analytic results were obtained using a linear analysis and were later confirmed by means of non-linear numerical simulations which did not resolve the flow within the Bondi sphere of the perturber, hence effectively restricted to weakly perturbed regions of the flow (paper~I). Here we present high resolution simulations, using either 3D Cartesian or 2D cylindrical meshes that resolve the flow within the Bondi sphere. We perform a systematic study of the force as a function of the perturber's mass and luminosity, in the subsonic regime. We find that perturbers with mass $M$ smaller than a few $M_c\sim χc_s/G$ are subjected to a thermal force with a magnitude in good agreement with linear theory ($χ$ being the thermal diffusivity of the medium, $c_s$ the adiabatic sound speed and $G$ the gravitational constant), while for larger masses, the thermal forces are only a fraction of the linear estimate that decays as $M^{-1}$. Our analysis confirms the possibility of negative friction (hence a propulsion) on sufficiently luminous, low-mass embryos embedded in protoplanetary discs. Finally, we give an approximate expression of the total force at low Mach number, valid both for sub-critical ($M<M_c$) and super-critical ($M>M_c$) perturbers. △ Less

Submitted 28 April, 2020; originally announced April 2020.

Thermal torque effects on the migration of growing low-mass planets

Authors: O. M. Guilera, N. Cuello, M. Montesinos, M. M. Miller Bertolami, M. P. Ronco, J. Cuadra, F. S. Masset

Abstract: As planets grow the exchange of angular momentum with the gaseous component of the protoplanetary disc produces a net torque resulting in a variation of the semi-major axis of the planet. For low-mass planets not able to open a gap in the gaseous disc this regime is known as type I migration. Pioneer works studied this mechanism in isothermal discs finding fast inward type I migration rates that w… ▽ More As planets grow the exchange of angular momentum with the gaseous component of the protoplanetary disc produces a net torque resulting in a variation of the semi-major axis of the planet. For low-mass planets not able to open a gap in the gaseous disc this regime is known as type I migration. Pioneer works studied this mechanism in isothermal discs finding fast inward type I migration rates that were unable to reproduce the observed properties of extrasolar planets. In the last years, several improvements have been made in order to extend the study of type I migration rates to non-isothermal discs. Moreover, it was recently shown that if the planet's luminosity due to solid accretion is taken into account, inward migration could be slowed down and even reversed. In this work, we study the planet formation process incorporating, and comparing, updated type I migration rates for non-isothermal discs and the role of planet's luminosity over such rates. We find that the latter can have important effects on planetary evolution, producing a significant outward migration for the growing planets. △ Less

Submitted 24 April, 2019; originally announced April 2019.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society

Impact of thermal effects on the evolution of eccentricity and inclination of low-mass planets

Authors: Sebastien Fromenteau, Frederic Masset

Abstract: Using linear perturbation theory, we evaluate the time-dependent force exerted on an eccentric and inclined low-mass planet embedded in a gaseous protoplanetary disc with finite thermal diffusivity $χ$. We assume the eccentricity and inclination to be small compared to the size of the thermal lobes $λ\sim(χ/Ω)^{1/2}$, itself generally much smaller than the scalelength of pressure $H$. When the pla… ▽ More Using linear perturbation theory, we evaluate the time-dependent force exerted on an eccentric and inclined low-mass planet embedded in a gaseous protoplanetary disc with finite thermal diffusivity $χ$. We assume the eccentricity and inclination to be small compared to the size of the thermal lobes $λ\sim(χ/Ω)^{1/2}$, itself generally much smaller than the scalelength of pressure $H$. When the planet is non-luminous, we find that its eccentricity and inclination are vigorously damped by the disc, over a timescale shorter by a factor $H/λ$ than the damping timescale in adiabatic discs. On the contrary, when the luminosity-to-mass ratio of the planet exceeds a threshold that depends on the disc's properties, its eccentricity and inclination undergo an exponential growth. In the limit of a large luminosity, the growth rate of the eccentricity is 2.5~times larger than that of the inclination, in agreement with previous numerical work. Depending on their luminosity, planetary embryos therefore exhibit much more diverse behaviours than the mild damping of eccentricity and inclination considered hitherto. △ Less

Submitted 26 March, 2019; v1 submitted 11 March, 2019; originally announced March 2019.

Comments: 16 pages, 4 figures. Minor amendments to match the version published by MNRAS

Numerical study of dynamical friction with thermal effects -- I. Comparison to linear theory

Authors: David A. Velasco Romero, Frédéric S. Masset

Abstract: We investigate by means of numerical simulations with nested Cartesian meshes the force exerted on a massive and luminous perturber moving at constant speed through a homogeneous and opaque gas, taking into account thermal diffusion in the gas and the radiative feedback from the perturber. The force arising from the release of energy into the ambient medium by the perturber, or heating force, is d… ▽ More We investigate by means of numerical simulations with nested Cartesian meshes the force exerted on a massive and luminous perturber moving at constant speed through a homogeneous and opaque gas, taking into account thermal diffusion in the gas and the radiative feedback from the perturber. The force arising from the release of energy into the ambient medium by the perturber, or heating force, is directed along the direction of motion and induces an acceleration of the perturber. Its value is compared to analytic estimates in the low and high Mach number limits, and found to match those accurately. In addition, the drag exerted on a non-luminous perturber significantly departs from the adiabatic expression when thermal diffusion is taken into account. In the limit of a vanishing velocity, this drag tends to a finite value which we determine using linear perturbation theory and corroborate with numerical simulations. The drag on a non-luminous perturber in a non-adiabatic gas therefore behaves like dry or solid friction. We work out the luminosity threshold to get a net acceleration of the perturber and find it to be generally much smaller than the luminosity of accreting low-mass planetary embryos embedded in a gaseous protoplanetary disc at a few astronomical units. We also present in some detail our implementation of nested meshes, which runs in parallel over several \emph{Graphics Processing Units} (GPUs). △ Less

Submitted 7 December, 2018; originally announced December 2018.

Comments: 15 pages and 16 figures

Planet-disc interactions with Discontinuous Galerkin Methods using GPUs

Authors: David A. Velasco-Romero, Maria Han Veiga, Romain Teyssier, Frédéric S. Masset

Abstract: We present a two-dimensional Cartesian code based on high order discontinuous Galerkin methods, implemented to run in parallel over multiple GPUs. A simple planet-disc setup is used to compare the behaviour of our code against the behaviour found using the FARGO3D code with a polar mesh. We make use of the time dependence of the torque exerted by the disc on the planet as a mean to quantify the nu… ▽ More We present a two-dimensional Cartesian code based on high order discontinuous Galerkin methods, implemented to run in parallel over multiple GPUs. A simple planet-disc setup is used to compare the behaviour of our code against the behaviour found using the FARGO3D code with a polar mesh. We make use of the time dependence of the torque exerted by the disc on the planet as a mean to quantify the numerical viscosity of the code. We find that the numerical viscosity of the Keplerian flow can be as low as a few $10^{-8}r^2Ω$, $r$ and $Ω$ being respectively the local orbital radius and frequency, for fifth order schemes and resolution of $\sim 10^{-2}r$. Although for a single disc problem a solution of low numerical viscosity can be obtained at lower computational cost with FARGO3D (which is nearly an order of magnitude faster than a fifth order method), discontinuous Galerkin methods appear promising to obtain solutions of low numerical viscosity in more complex situations where the flow cannot be captured on a polar or spherical mesh concentric with the disc. △ Less

Submitted 3 May, 2018; originally announced May 2018.

Comments: Accepted for publication in MNRAS

Coorbital thermal torques on low-mass protoplanets

Authors: Frederic S. Masset

Abstract: Using linear perturbation theory, we investigate the torque exerted on a low-mass planet embedded in a gaseous protoplanetary disc with finite thermal diffusivity. When the planet does not release energy into the ambient disc, the main effect of thermal diffusion is the softening of the enthalpy peak near the planet, which results in the appearance of two cold and dense lobes on either side of the… ▽ More Using linear perturbation theory, we investigate the torque exerted on a low-mass planet embedded in a gaseous protoplanetary disc with finite thermal diffusivity. When the planet does not release energy into the ambient disc, the main effect of thermal diffusion is the softening of the enthalpy peak near the planet, which results in the appearance of two cold and dense lobes on either side of the orbit, of size smaller than the thickness of the disc. The lobes exert torques of opposite sign on the planet, each comparable in magnitude to the one-sided Lindblad torque. When the planet is offset from corotation, the lobes are asymmetric and the planet experiences a net torque, the `cold' thermal torque, which has a magnitude that depends on the relative value of the distance to corotation to the size of the lobes $\sim\sqrt{χ/Ω_p}$, $χ$ being the thermal diffusivity and $Ω_p$ the orbital frequency. We believe that this effect corresponds to the phenomenon named `cold finger' recently reported in numerical simulations, and we argue that it constitutes the dominant mode of migration of sub-Earth-mass objects. When the planet is luminous, the heat released into the ambient disc results in an additional disturbance that takes the form of hot, low-density lobes. They give a torque, named heating torque in previous work, that has an expression similar, but of opposite sign, to the cold thermal torque. △ Less

Submitted 14 September, 2017; v1 submitted 31 August, 2017; originally announced August 2017.

Comments: In press, updated to match MNRAS corrected proofs

Orbital advection with magnetohydrodynamics and vector potential

Authors: Wladimir Lyra, Colin P. McNally, Tobias Heinemann, Frederic Masset

Abstract: Orbital advection is a significant bottleneck in disk simulations, and a particularly tricky one when used in connection with magnetohydrodynamics. We have developed an orbital advection algorithm suitable for the induction equation with magnetic potential. The electromotive force is split into advection and shear term, and we find that we do not need an advective gauge as solving the orbital adve… ▽ More Orbital advection is a significant bottleneck in disk simulations, and a particularly tricky one when used in connection with magnetohydrodynamics. We have developed an orbital advection algorithm suitable for the induction equation with magnetic potential. The electromotive force is split into advection and shear term, and we find that we do not need an advective gauge as solving the orbital advection implicitly precludes the shear term from canceling the advection term. We prove and demonstrate the third order in time accuracy of the scheme. The algorithm is also suited for non-magnetic problems. Benchmarked results of (hydrodynamical) planet-disk interaction and the of the magnetorotational instability are reproduced. We include detailed descriptions of the construction and selection of stabilizing dissipations (or high frequency filters) needed to generate practical results. The scheme is self-consistent, accurate, and elegant in its simplicity, making it particularly efficiently in straightforward finite-difference methods. As a result of the work, the algorithm is incorporated in the public version of the Pencil Code, where it can be used by the community. △ Less

Submitted 3 September, 2017; v1 submitted 16 August, 2017; originally announced August 2017.

Comments: 11 pages + 4 appendix; 8+1 figures. Accepted by AJ; in press

Improved torque formula for low and intermediate mass planetary migration

Authors: María Alejandra Jiménez, Frédéric S. Masset

Abstract: The migration of planets on nearly circular, non-inclined orbits in protoplanetary discs is entirely described by the disc's torque. This torque is a complex function of the disc parameters, and essentially amounts to the sum of two components: the Lindblad torque and the corotation torque. Known torque formulae do not reproduce accurately the torque actually experienced in numerical simulations b… ▽ More The migration of planets on nearly circular, non-inclined orbits in protoplanetary discs is entirely described by the disc's torque. This torque is a complex function of the disc parameters, and essentially amounts to the sum of two components: the Lindblad torque and the corotation torque. Known torque formulae do not reproduce accurately the torque actually experienced in numerical simulations by low- and intermediate- mass planets in radiative discs. One of the main reasons for this inaccuracy is that these formulae have been worked out in two-dimensional analyses. Here we revisit the torque formula and update many of its dimensionless coefficients by means of tailored, three- dimensional numerical simulations. In particular, we derive the dependence of the Lindblad torque on the temperature gradient, the dependence of the corotation torque on the radial entropy gradient (and work out a suitable expression of this gradient in a three-dimensional disc). We also work out the dependence of the corotation torque on the radial temperature gradient, overlooked so far. Corotation torques are known to scale very steeply with the width of the horseshoe region. We extend the expression of this width to the domain of intermediate mass planets, so that our updated torque formula remains valid for planets up to typically several tens of Earth masses, provided these relatively massive planets do not significantly deplete their coorbital region. Our torque expression can be applied to low- and intermediate-mass planets in optically thick protoplanetary discs, as well as protomoons embedded in circumplanetary discs. △ Less

Submitted 27 July, 2017; originally announced July 2017.

Comments: 14 pages,9 figures, accepted for publication in MNRAS

Evolution of eccentricity and inclination of hot protoplanets embedded in radiative discs

Authors: Henrik Eklund, Frédéric S. Masset

Abstract: We study the evolution of the eccentricity and inclination of protoplanetary embryos and low-mass protoplanets (from a fraction of an Earth mass to a few Earth masses) embedded in a protoplanetary disc, by means of three dimensional hydrodynamics calculations with radiative transfer in the diffusion limit. When the protoplanets radiate in the surrounding disc the energy released by the accretion o… ▽ More We study the evolution of the eccentricity and inclination of protoplanetary embryos and low-mass protoplanets (from a fraction of an Earth mass to a few Earth masses) embedded in a protoplanetary disc, by means of three dimensional hydrodynamics calculations with radiative transfer in the diffusion limit. When the protoplanets radiate in the surrounding disc the energy released by the accretion of solids, their eccentricity and inclination experience a growth toward values which depend on the luminosity to mass ratio of the planet, which are comparable to the disc's aspect ratio and which are reached over timescales of a few thousand years. This growth is triggered by the appearance of a hot, under-dense region in the vicinity of the planet. The growth rate of the eccentricity is typically three times larger than that of the inclination. In long term calculations, we find that the excitation of eccentricity and the excitation of inclination are not independent. In the particular case in which a planet has initially a very small eccentricity and inclination, the eccentricity largely overruns the inclination. When the eccentricity reaches its asymptotic value, the growth of inclination is quenched, yielding an eccentric orbit with a very low inclination. As a side result, we find that the eccentricity and inclination of non-luminous planets are damped more vigorously in radiative discs than in isothermal discs. △ Less

Submitted 29 June, 2017; v1 submitted 6 April, 2017; originally announced April 2017.

Comments: 13 pages, 16 figures. The videos showing the temperature excess around the planet along the epicycle and inclined motion can be found at https://www.fis.unam.mx/~masset/EklundMasset2017/

Journal ref: Mon Not R Astron Soc (2017) 469 (1): 206-217

Gap formation by inclined massive planets in locally isothermal three-dimensional discs

Authors: R. O. Chametla, F. J. Sanchez-Salcedo, F. S. Masset, A. M. Hidalgo-Gamez

Abstract: We study gap formation in gaseous protoplanetary discs by a Jupiter mass planet. The planet's orbit is circular and inclined relative to the midplane of the disc. We use the impulse approximation to estimate the gravitational tidal torque between the planet and the disc, and infer the gap profile. For low-mass discs, we provide a criterion for gap opening when the orbital inclination is… ▽ More We study gap formation in gaseous protoplanetary discs by a Jupiter mass planet. The planet's orbit is circular and inclined relative to the midplane of the disc. We use the impulse approximation to estimate the gravitational tidal torque between the planet and the disc, and infer the gap profile. For low-mass discs, we provide a criterion for gap opening when the orbital inclination is $\leq 30^{\circ}$. Using the FARGO3D code, we simulate the disc response to an inclined massive planet. The dependence of the depth and width of the gap obtained in the simulations on the inclination of the planet is broadly consistent with the scaling laws derived in the impulse approximation. Although we mainly focus on planets kept on fixed orbits, the formalism permits to infer the temporal evolution of the gap profile in cases where the inclination of the planet changes with time. This study may be useful to understand the migration of massive planets on inclined orbit, because the strength of the interaction with the disc depends on whether a gap is opened or not. △ Less

Submitted 5 April, 2017; originally announced April 2017.

Comments: 16 pages, 13 figures, accepted for publication in MNRAS

Dynamical friction on hot bodies in opaque, gaseous media

Authors: Frédéric S. Masset, David A. Velasco Romero

Abstract: We consider the gravitational force exerted on a point-like perturber of mass $M$ travelling within a uniform gaseous, opaque medium at constant velocity $V$. The perturber irradiates the surrounding gas with luminosity $L$. The diffusion of the heat released is modelled with a uniform thermal diffusivity $χ$. Using linear perturbation theory, we show that the force exerted by the perturbed gas on… ▽ More We consider the gravitational force exerted on a point-like perturber of mass $M$ travelling within a uniform gaseous, opaque medium at constant velocity $V$. The perturber irradiates the surrounding gas with luminosity $L$. The diffusion of the heat released is modelled with a uniform thermal diffusivity $χ$. Using linear perturbation theory, we show that the force exerted by the perturbed gas on the perturber differs from the force without radiation (or standard dynamical friction). Hot, underdense gas trails the mass, which gives rise to a new force component, the heating force, with direction $+V$, thus opposed to the standard dynamical friction. In the limit of low Mach numbers, the heating force has expression $F_\mathrm{heat}=γ(γ-1)GML/(2χc_s^2)$, $c_s$ being the sound speed and $γ$ the ratio of specific heats. In the limit of large Mach numbers, $F_\mathrm{heat}=(γ-1)GML/(χV^2)f(r_\mathrm{min}V/4χ)$, where $f$ is a function that diverges logarithmically as $r_\mathrm{min}$ tends to zero. Remarkably, the force in the low Mach number limit does not depend on the velocity. The equilibrium speed, when it exists, is set by the cancellation of the standard dynamical friction and heating force. In the low Mach number limit, it scales with the luminosity to mass ratio of the perturber. Using the above results suggests that Mars- to Earth-sized planetary embryos heated by accretion in a gaseous protoplanetary disc should have eccentricities and inclinations that amount to a sizeable fraction of the disc's aspect ratio, for conditions thought to prevail at a few astronomical units. △ Less

Submitted 17 November, 2016; originally announced November 2016.

Comments: 11 pages, 3 figures, accepted for publication in MNRAS

Journal ref: Mon Not R Astron Soc (2017) 465 (3): 3175-3184

Planetary Torque in 3D Isentropic Disks

Authors: Jeffrey Fung, Frederic Masset, Elena Lega, David Velasco

Abstract: 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… ▽ More 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 ($γ$). 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 $γ$, up to supersonic speeds for the smallest $r_{\rm s}$ and $γ$ in our study. △ Less

Submitted 11 January, 2017; v1 submitted 28 October, 2016; originally announced October 2016.

Comments: Accepted to ApJ

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

Authors: Stefan Kraus, John D. Monnier, Michael J. Ireland, Gaspard Duchene, Catherine Espaillat, Sebastian Hoenig, Attila Juhasz, Chris Mordasini, Johan Olofsson, Claudia Paladini, Keivan Stassun, Neal Turner, Gautam Vasisht, Tim J. Harries, Matthew R. Bate, Jean-Francois Gonzalez, Alexis Matter, Zhaohuan Zhu, Olja Panic, Zsolt Regaly, Alessandro Morbidelli, Farzana Meru, Sebastian Wolf, John Ilee, Jean-Philippe Berger , et al. (53 additional authors not shown)

Abstract: 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 prot… ▽ More 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. △ Less

Submitted 16 August, 2016; v1 submitted 1 August, 2016; originally announced August 2016.

Comments: 13 pages, 5 figures, Proceedings of SPIE 2016

Circumplanetary disk or circumplanetary envelope?

Authors: J. Szulágyi, F. Masset, E. Lega, A. Crida, A. Morbidelli, T. Guillot

Abstract: 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 diamete… ▽ More 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. △ Less

Submitted 15 May, 2016; originally announced May 2016.

Comments: Accepted for publication at MNRAS

Long-term and large-scale hydrodynamical simulations of migrating planets

Authors: Pablo Benítez-Llambay, Ximena S. Ramos, Cristian Beaugé, Frédéric Masset

Abstract: We present a new method that allows long-term and large-scale hydrodynamical simulations of migrating planets over a grid-based Eulerian code. This technique, which consists in a remapping of the disk by tracking the planetary migration, enables runs of migrating planets over a time comparable to the age of protoplanetary disks. This method also has the potential to address efficiently problems re… ▽ More We present a new method that allows long-term and large-scale hydrodynamical simulations of migrating planets over a grid-based Eulerian code. This technique, which consists in a remapping of the disk by tracking the planetary migration, enables runs of migrating planets over a time comparable to the age of protoplanetary disks. This method also has the potential to address efficiently problems related with migration of multi-planet systems in gaseous disks, and to improve current results of migration of massive planets by including global viscous evolution as well as detailed studies of the co-orbital region during migration. We perform different tests using the public code FARGO3D to validate this method and compare its results with those obtained using a classical fixed grid. △ Less

Submitted 5 May, 2016; originally announced May 2016.

Comments: Accepted for publication in ApJ. For a movie describing the method, see https://youtu.be/66o0Z2lX8NU

FARGO3D: A new GPU-oriented MHD code

Authors: Pablo Benítez-Llambay, Frédéric Masset

Abstract: We present the FARGO3D code, recently publicly released. It is a magnetohydrodynamics code developed with special emphasis on protoplanetary disks physics and planet-disk interactions, and parallelized with MPI. The hydrodynamics algorithms are based on finite difference upwind, dimensionally split methods. The magnetohydrodynamics algorithms consist of the constrained transport method to preserve… ▽ More We present the FARGO3D code, recently publicly released. It is a magnetohydrodynamics code developed with special emphasis on protoplanetary disks physics and planet-disk interactions, and parallelized with MPI. The hydrodynamics algorithms are based on finite difference upwind, dimensionally split methods. The magnetohydrodynamics algorithms consist of the constrained transport method to preserve the divergence-free property of the magnetic field to machine accuracy, coupled to a method of characteristics for the evaluation of electromotive forces and Lorentz forces. Orbital advection is implemented, and an N-body solver is included to simulate planets or stars interacting with the gas. We present our implementation in detail and present a number of widely known tests for comparison purposes. One strength of FARGO3D is that it can run on both "Graphical Processing Units" (GPUs) or "Central Processing unit" (CPUs), achieving large speed up with respect to CPU cores. We describe our implementation choices, which allow a user with no prior knowledge of GPU programming to develop new routines for the CPU, and have them translated automatically for the GPU. △ Less

Submitted 7 February, 2016; originally announced February 2016.

Comments: Accepted for publication in ApJS

Horseshoe drag in three-dimensional globally isothermal disks

Authors: Frédéric Masset, Pablo Benítez-Llambay

Abstract: We study the horseshoe dynamics of a low-mass planet in a three-dimensional, globally isothermal, inviscid disk. We find, as reported in previous work, that the boundaries of the horseshoe region (separatrix sheets) have cylindrical symmetry about the disk's rotation axis. We interpret this feature as arising from the fact that the whole separatrix sheets have a unique value of Bernoulli's constan… ▽ More We study the horseshoe dynamics of a low-mass planet in a three-dimensional, globally isothermal, inviscid disk. We find, as reported in previous work, that the boundaries of the horseshoe region (separatrix sheets) have cylindrical symmetry about the disk's rotation axis. We interpret this feature as arising from the fact that the whole separatrix sheets have a unique value of Bernoulli's constant, and that this constant does not depend on altitude, but only on the cylindrical radius, in barotropic disks. We next derive an expression for the torque exerted by the horseshoe region onto the planet, or horseshoe drag. Potential vorticity is not materially conserved as in two-dimensional flows, but it obeys a slightly more general conservation law (Ertel's theorem) which allows to obtain an expression for the horseshoe drag identical to the expression in a two-dimensional disk. Our results are illustrated and validated by three-dimensional numerical simulations. The horseshoe region is found to be slightly more narrow than previously extrapolated from two-dimensional analyses with a suitable softening length of the potential. We discuss the implications of our results for the saturation of the corotation torque, and the possible connection to the flow at the Bondi scale, which the present analysis does not resolve. △ Less

Submitted 24 November, 2015; originally announced November 2015.

Comments: Accepted for publication in ApJ

Planet heating prevents inward migration of planetary cores

Authors: Pablo Benítez-Llambay, Frédéric Masset, Gloria Koenigsberger, Judit Szulágyi

Abstract: Planetary systems are born in the disks of gas, dust and rocky fragments that surround newly formed stars. Solid content assembles into ever-larger rocky fragments that eventually become planetary embryos. These then continue their growth by accreting leftover material in the disc. Concurrently, tidal effects in the disc cause a radial drift in the embryo orbits, a process known as migration. Fast… ▽ More Planetary systems are born in the disks of gas, dust and rocky fragments that surround newly formed stars. Solid content assembles into ever-larger rocky fragments that eventually become planetary embryos. These then continue their growth by accreting leftover material in the disc. Concurrently, tidal effects in the disc cause a radial drift in the embryo orbits, a process known as migration. Fast inward migration is predicted by theory for embryos smaller than three to five Earth masses. With only inward migration, these embryos can only rarely become giant planets located at Earth's distance from the Sun and beyond, in contrast with observations. Here we report that asymmetries in the temperature rise associated with accreting infalling material produce a force (which gives rise to an effect that we call "heating torque") that counteracts inward migration. This provides a channel for the formation of giant planets and also explains the strong planet-metallicity correlation found between the incidence of giant planets and the heavy-element abundance of the host stars. △ Less

Submitted 6 October, 2015; originally announced October 2015.

Comments: 19 pages, 4 figures

Journal ref: Nature 520, 63-65 (02 April 2015)

The explosion mechanism of core-collapse supernovae: progress in supernova theory and experiments

Authors: Thierry Foglizzo, Rémi Kazeroni, Jérôme Guilet, Frédéric Masset, Matthias González, Brendan K. Krueger, Jérôme Novak, Micaela Oertel, Jérôme Margueron, Julien Faure, Noël Martin, Patrick Blottiau, Bruno Peres, Gilles Durand

Abstract: The explosion of core-collapse supernova depends on a sequence of events taking place in less than a second in a region of a few hundred kilometers at the center of a supergiant star, after the stellar core approaches the Chandrasekhar mass and collapses into a proto-neutron star, and before a shock wave is launched across the stellar envelope. Theoretical efforts to understand stellar death focus… ▽ More The explosion of core-collapse supernova depends on a sequence of events taking place in less than a second in a region of a few hundred kilometers at the center of a supergiant star, after the stellar core approaches the Chandrasekhar mass and collapses into a proto-neutron star, and before a shock wave is launched across the stellar envelope. Theoretical efforts to understand stellar death focus on the mechanism which transforms the collapse into an explosion. Progress in understanding this mechanism is reviewed with particular attention to its asymmetric character. We highlight a series of successful studies connecting observations of supernova remnants and pulsars properties to the theory of core-collapse using numerical simulations. The encouraging results from first principles models in axisymmetric simulations is tempered by new puzzles in 3D. The diversity of explosion paths and the dependence on the pre-collapse stellar structure is stressed, as well as the need to gain a better understanding of hydrodynamical and MHD instabilities such as SASI and neutrino-driven convection. The shallow water analogy of shock dynamics is presented as a comparative system where buoyancy effects are absent. This dynamical system can be studied numerically and also experimentally with a water fountain. The potential of this complementary research tool for supernova theory is analyzed. We also review its potential for public outreach in science museums. △ Less

Submitted 6 February, 2015; v1 submitted 6 January, 2015; originally announced January 2015.

Comments: 19 pages, 6 figures, invited review accepted for publication in PASA

Accretion of Jupiter-mass Planets in the Limit of Vanishing Viscosity

Authors: J. Szulágyi, A. Morbidelli, A. Crida, F. Masset

Abstract: In the core-accretion model the nominal runaway gas-accretion phase brings most planets to multiple Jupiter masses. However, known giant planets are predominantly Jupiter-mass bodies. Obtaining longer timescales for gas accretion may require using realistic equations of states, or accounting for the dynamics of the circumplanetary disk (CPD) in low-viscosity regime, or both. Here we explore the se… ▽ More In the core-accretion model the nominal runaway gas-accretion phase brings most planets to multiple Jupiter masses. However, known giant planets are predominantly Jupiter-mass bodies. Obtaining longer timescales for gas accretion may require using realistic equations of states, or accounting for the dynamics of the circumplanetary disk (CPD) in low-viscosity regime, or both. Here we explore the second way using global, three-dimensional isothermal hydrodynamical simulations with 8 levels of nested grids around the planet. In our simulations the vertical inflow from the circumstellar disk (CSD) to the CPD determines the shape of the CPD and its accretion rate. Even without prescribed viscosity Jupiter's mass-doubling time is $\sim 10^4$ years, assuming the planet at 5.2 AU and a Minimum Mass Solar Nebula. However, we show that this high accretion rate is due to resolution-dependent numerical viscosity. Furthermore, we consider the scenario of a layered CSD, viscous only in its surface layer, and an inviscid CPD. We identify two planet-accretion mechanisms that are independent of the viscosity in the CPD: (i) the polar inflow -- defined as a part of the vertical inflow with a centrifugal radius smaller than 2 Jupiter-radii and (ii) the torque exerted by the star on the CPD. In the limit of zero effective viscosity, these two mechanisms would produce an accretion rate 40 times smaller than in the simulation. △ Less

Submitted 1 January, 2014; v1 submitted 21 December, 2013; originally announced December 2013.

Comments: 11 pages, 10 figures. Accepted for publication in Astrophysical Journal

Planet-Disc Interactions and Early Evolution of Planetary Systems

Authors: Clément Baruteau, Aurélien Crida, Sijme-Jan Paardekooper, Frédéric Masset, Jérôme Guilet, Bertram Bitsch, Richard P. Nelson, Wilhelm Kley, John C. B. Papaloizou

Abstract: 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 charact… ▽ More 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. △ Less

Submitted 16 December, 2013; originally announced December 2013.

Comments: 24 pages, 9 figures, refereed review chapter accepted for publication in Protostars and Planets VI, University of Arizona Press (2014), eds. H. Beuther, R. Klessen, C. Dullemond, Th. Henning. The talk presented at the Protostars and Planets VI conference is available at http://www.youtube.com/watch?v=_HMw4Lh7IOo

Low-mass planets in nearly inviscid disks: Numerical treatment

Authors: Wilhelm Kley, Tobias W. A. Mueller, Stefan M. Kolb, Pablo Benitez-Llambay, Frederic Masset

Abstract: 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 resolu… ▽ More 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. △ Less

Submitted 15 August, 2012; originally announced August 2012.

Comments: accepted by Astronomy & Astrophysics

On type-I migration near opacity transitions. A generalized Lindblad torque formula for planetary population synthesis

Authors: F. S. Masset

Abstract: We give an expression for the Lindblad torque acting on a low-mass planet embedded in a protoplanetary disk that is valid even at locations where the surface density or temperature profile cannot be approximated by a power law, such as an opacity transition. At such locations, the Lindblad torque is known to suffer strong deviation from its standard value, with potentially important implications f… ▽ More We give an expression for the Lindblad torque acting on a low-mass planet embedded in a protoplanetary disk that is valid even at locations where the surface density or temperature profile cannot be approximated by a power law, such as an opacity transition. At such locations, the Lindblad torque is known to suffer strong deviation from its standard value, with potentially important implications for type I migration, but the full treatment of the tidal interaction is cumbersome and not well suited to models of planetary population synthesis. The expression that we propose retains the simplicity of the standard Lindblad torque formula and gives results that accurately reproduce those of numerical simulations, even at locations where the disk temperature undergoes abrupt changes. Our study is conducted by means of customized numerical simulations in the low-mass regime, in locally isothermal disks, and compared to linear torque estimates obtained by summing fully analytic torque estimates at each Lindblad resonance. The functional dependence of our modified Lindblad torque expression is suggested by an estimate of the shift of the Lindblad resonances that mostly contribute to the torque, in a disk with sharp gradients of temperature or surface density, while the numerical coefficients of the new terms are adjusted to seek agreement with numerics. As side results, we find that the vortensity related corotation torque undergoes a boost at an opacity transition that can counteract migration, and we find evidence from numerical simulations that the linear corotation torque has a non-negligible dependency upon the temperature gradient, in a locally isothermal disk. △ Less

Submitted 13 June, 2012; originally announced June 2012.

Comments: Appeared in special issue of "Celestial Mechanics and Dynamical Astronomy" on Extrasolar Planetary Systems

Journal ref: 2011CeMDA, 111, 131

No snow-plough mechanism during the rapid hardening of supermassive black hole binaries

Authors: Clément Baruteau, Enrico Ramirez-Ruiz, Frédéric Masset

Abstract: We present two-dimensional hydrodynamical simulations of the tidal interaction between a supermassive black hole binary with moderate mass ratio, and the fossil gas disc where it is embedded. Our study extends previous one-dimensional height-integrated disc models, which predicted that the density of the gas disc between the primary and the secondary black holes should rise significantly during th… ▽ More We present two-dimensional hydrodynamical simulations of the tidal interaction between a supermassive black hole binary with moderate mass ratio, and the fossil gas disc where it is embedded. Our study extends previous one-dimensional height-integrated disc models, which predicted that the density of the gas disc between the primary and the secondary black holes should rise significantly during the ultimate stages of the binary's hardening driven by the gravitational radiation torque. This snow-plough mechanism, as we call it, would lead to an increase in the bolometric luminosity of the system prior to the binary merger, which could be detected in conjunction with the gravitational wave signal. We argue here that the snow-plough mechanism is unlikely to occur. In two-dimensions, when the binary's hardening timescale driven by gravitational radiation becomes shorter than the disc's viscous drift timescale, fluid elements in the inner disc get funneled to the outer disc through horseshoe trajectories with respect to the secondary. Mass leakage across the secondary's gap is thus found to be effective and, as a result, the predicted accretion disc luminosity will remain at roughly the same level prior to merger. △ Less

Submitted 19 March, 2012; originally announced March 2012.

Comments: 5 pages, 5 figures, accepted for publication in MNRAS

Recent developments in planet migration theory

Authors: Clément Baruteau, Frédéric Masset

Abstract: Planetary migration is the process by which a forming planet undergoes a drift of its semi-major axis caused by the tidal interaction with its parent protoplanetary disc. One of the key quantities to assess the migration of embedded planets is the tidal torque between the disc and planet, which has two components: the Lindblad torque and the corotation torque. We review the latest results on both… ▽ More Planetary migration is the process by which a forming planet undergoes a drift of its semi-major axis caused by the tidal interaction with its parent protoplanetary disc. One of the key quantities to assess the migration of embedded planets is the tidal torque between the disc and planet, which has two components: the Lindblad torque and the corotation torque. We review the latest results on both torque components for planets on circular orbits, with a special emphasis on the various processes that give rise to additional, large components of the corotation torque, and those contributing to the saturation of this torque. These additional components of the corotation torque could help address the shortcomings that have recently been exposed by models of planet population syntheses. We also review recent results concerning the migration of giant planets that carve gaps in the disc (type II migration) and the migration of sub-giant planets that open partial gaps in massive discs (type III migration). △ Less

Submitted 15 March, 2012; originally announced March 2012.

Comments: 52 pages, 18 figures. Review article to be published in "Tidal effects in Astronomy and Astrophysics", Lecture Notes in Physics

A Shallow Water Analogue of the Standing Accretion Shock Instability: Experimental Demonstration and Two-Dimensional Model

Authors: Thierry Foglizzo, Frédéric Masset, Jérôme Guilet, Gilles Durand

Abstract: Despite the sphericity of the collapsing stellar core, the birth conditions of neutron stars can be highly non spherical due to a hydrodynamical instability of the shocked accretion flow. Here we report the first laboratory experiment of a shallow water analogue, based on the physics of hydraulic jumps. Both the experiment and its shallow water modeling demonstrate a robust linear instability and… ▽ More Despite the sphericity of the collapsing stellar core, the birth conditions of neutron stars can be highly non spherical due to a hydrodynamical instability of the shocked accretion flow. Here we report the first laboratory experiment of a shallow water analogue, based on the physics of hydraulic jumps. Both the experiment and its shallow water modeling demonstrate a robust linear instability and nonlinear properties of symmetry breaking, in a system which is one million times smaller and about hundred times slower than its astrophysical analogue. △ Less

Submitted 15 December, 2011; originally announced December 2011.

Comments: 4 pages, 4 figures, accepted for publication in Phys. Rev. Letters. Supplementary Material (6 movies) available at http://irfu.cea.fr/Projets/SN2NS/outreach.html

Journal ref: Phys. Rev. Lett. 108, 051103 (2012)

Corotation torques experienced by planets embedded in weakly magnetized turbulent discs

Authors: C. Baruteau, S. Fromang, R. P. Nelson, F. Masset

Abstract: The migration of low-mass planets is driven by the differential Lindblad torque and the corotation torque in non-magnetic viscous models of protoplanetary discs. The corotation torque has recently received detailed attention as it may slow down, stall, or reverse migration. In laminar viscous disc models, the long-term evolution of the corotation torque is intimately related to viscous and thermal… ▽ More The migration of low-mass planets is driven by the differential Lindblad torque and the corotation torque in non-magnetic viscous models of protoplanetary discs. The corotation torque has recently received detailed attention as it may slow down, stall, or reverse migration. In laminar viscous disc models, the long-term evolution of the corotation torque is intimately related to viscous and thermal diffusion processes in the planet's horseshoe region. This paper examines the properties of the corotation torque in discs where MHD turbulence develops as a result of the magnetorotational instability, considering a weak initial toroidal magnetic field. We present results of 3D MHD simulations carried out with two different codes. Non-ideal MHD effects and the disc's vertical stratification are neglected, and locally isothermal disc models are considered. The running time-averaged torque exerted by the disc on a fixed planet is evaluated in three disc models. We first present results with an inner disc cavity (planet trap). As in viscous disc models, the planet is found to experience a positive running time-averaged torque over several hundred orbits, which highlights the existence of an unsaturated corotation torque maintained in the long term in MHD turbulent discs. Two disc models with initial power-law density and temperature profiles are also adopted, in which the time-averaged torque is found to be in decent agreement with its counterpart in laminar viscous disc models with similar viscosity at the planet location. Detailed analysis of the averaged torque density distributions indicates that the differential Lindblad torque takes very similar values in MHD turbulent and laminar viscous discs, and there exists an unsaturated corotation torque in MHD turbulent discs. This analysis also reveals the existence of an additional corotation torque in weakly magnetized discs. △ Less

Submitted 14 July, 2011; originally announced July 2011.

Comments: 16 pages, 13 figures, accepted for publication in Astronomy and Astrophysics

Meridional circulation in turbulent protoplanetary disks

Authors: Sebastien Fromang, Wladimir Lyra, Frederic Masset

Abstract: Based on the viscous disk theory, a number of recent studies have suggested there is large scale meridional circulation in protoplanetary disks. Such a flow could account for the presence of crystalline silicates, including calcium- and aluminum-rich inclusions (CAIs), at large distances from the sun. This paper aims at examining whether such large-scale flows exist in turbulent protoplanetary dis… ▽ More Based on the viscous disk theory, a number of recent studies have suggested there is large scale meridional circulation in protoplanetary disks. Such a flow could account for the presence of crystalline silicates, including calcium- and aluminum-rich inclusions (CAIs), at large distances from the sun. This paper aims at examining whether such large-scale flows exist in turbulent protoplanetary disks. High-resolution global hydrodynamical and magnetohydrodynamical (MHD) numerical simulations of turbulent protoplanetary disks were used to infer the properties of the flow in such disks. By performing hydrodynamic simulations using explicit viscosity, we demonstrate that our numerical setup does not suffer from any numerical artifact. The aforementioned meridional circulation is easily recovered in viscous and laminar disks and is quickly established. In MHD simulations, the magnetorotational instability drives turbulence in the disks. Averaging out the turbulent fluctuations on a long timescale, the results fail to show any large-scale meridional circulation. A detailed analysis of the simulations show that this lack of meridional circulation is due to the turbulent stress tensor having a vertical profile different from the viscous stress tensor. A simple model is provided that successfully accounts for the structure of the flow in the bulk of the disk. In addition to those results, possible deviations from standard vertically averaged alpha disk models are suggested by the simulations and should be the focus of future work. Global MHD numerical simulations of fully ionized and turbulent protoplanetary disks are not consistent with the existence of a large-scale meridional flow. As a consequence, the presence of crystalline silicates at large distance for the central star cannot be accounted for by that process as suggested by recent models based on viscous disk theory. △ Less

Submitted 24 August, 2011; v1 submitted 3 February, 2011; originally announced February 2011.

Comments: 16 pages, 13 figures, changes according to referee report. Accepted to Astronomy and Astrophysics

The mass-period distribution of close-in exoplanets

Authors: P. Benitez-Llambay, F. Masset, C. Beauge

Abstract: The lower limit to the distribution of orbital periods P for the current population of close-in exoplanets shows a distinctive discontinuity located at approximately one Jovian mass. Most smaller planets have orbital periods longer than P~2.5 days, while higher masses are found down to P~1 day. We analyze whether this observed mass-period distribution could be explained in terms of the combined… ▽ More The lower limit to the distribution of orbital periods P for the current population of close-in exoplanets shows a distinctive discontinuity located at approximately one Jovian mass. Most smaller planets have orbital periods longer than P~2.5 days, while higher masses are found down to P~1 day. We analyze whether this observed mass-period distribution could be explained in terms of the combined effects of stellar tides and the interactions of planets with an inner cavity in the gaseous disk. We performed a series of hydrodynamical simulations of the evolution of single-planet systems in a gaseous disk with an inner cavity mimicking the inner boundary of the disk. The subsequent tidal evolution is analyzed assuming that orbital eccentricities are small and stellar tides are dominant. We find that most of the close-in exoplanet population is consistent with an inner edge of the protoplanetary disk being located at approximately P>2 days for solar-type stars, in addition to orbital decay having been caused by stellar tides with a specific tidal parameter on the order of Q'*=10^7. The data is broadly consistent with planets more massive than one Jupiter mass undergoing type II migration, crossing the gap, and finally halting at the interior 2/1 mean-motion resonance with the disk edge. Smaller planets do not open a gap in the disk and remain trapped in the cavity edge. CoRoT-7b appears detached from the remaining exoplanet population, apparently requiring additional evolutionary effects to explain its current mass and semimajor axis. △ Less

Submitted 18 January, 2011; originally announced January 2011.

Comments: 8 Pages, 8 figures, accepted for publication in A&A

Saturated torque formula for planetary migration in viscous disks with thermal diffusion: recipe for protoplanet population synthesis

Authors: Frederic S. Masset, Jules Casoli

Abstract: We provide torque formulae for low mass planets undergoing type I migration in gaseous disks. These torque formulae put special emphasis on the horseshoe drag, which is prone to saturation: the asymptotic value reached by the horseshoe drag depends on a balance between coorbital dynamics (which tends to cancel out or saturate the torque) and diffusive processes (which tend to restore the unperturb… ▽ More We provide torque formulae for low mass planets undergoing type I migration in gaseous disks. These torque formulae put special emphasis on the horseshoe drag, which is prone to saturation: the asymptotic value reached by the horseshoe drag depends on a balance between coorbital dynamics (which tends to cancel out or saturate the torque) and diffusive processes (which tend to restore the unperturbed disk profiles, thereby desaturating the torque). We entertain here the question of this asymptotic value, and we derive torque formulae which give the total torque as a function of the disk's viscosity and thermal diffusivity. The horseshoe drag features two components: one which scales with the vortensity gradient, and one which scales with the entropy gradient, and which constitutes the most promising candidate for halting inward type I migration. Our analysis, which is complemented by numerical simulations, recovers characteristics already noted by numericists, namely that the viscous timescale across the horseshoe region must be shorter than the libration time in order to avoid saturation, and that, provided this condition is satisfied, the entropy related part of the horseshoe drag remains large if the thermal timescale is shorter than the libration time. Side results include a study of the Lindblad torque as a function of thermal diffusivity, and a contribution to the corotation torque arising from vortensity viscously created at the contact discontinuities that appear at the horseshoe separatrices. For the convenience of the reader mostly interested in the torque formulae, section 8 is self-contained. △ Less

Submitted 26 October, 2010; v1 submitted 9 September, 2010; originally announced September 2010.

Comments: Affiliation details changed. Fixed equation numbering issue. Biblio info added

Journal ref: Astrophysical Journal 723 (2010), 1393-1417

Long range outward migration of giant planets, with application to Fomalhaut b

Authors: A. Crida, F. Masset, A. Morbidelli

Abstract: Recent observations of exoplanets by direct imaging, reveal that giant planets orbit at a few dozens to more than a hundred of AU from their central star. The question of the origin of these planets challenges the standard theories of planet formation. We propose a new way of obtaining such far planets, by outward migration of a pair of planets formed in the 10 AU region. Two giant planets in me… ▽ More Recent observations of exoplanets by direct imaging, reveal that giant planets orbit at a few dozens to more than a hundred of AU from their central star. The question of the origin of these planets challenges the standard theories of planet formation. We propose a new way of obtaining such far planets, by outward migration of a pair of planets formed in the 10 AU region. Two giant planets in mean motion resonance in a common gap in the protoplanetary disk migrate outwards, if the inner one is significantly more massive than the outer one. Using hydrodynamical simulations, we show that their semi major axes can increase by almost one order of magnitude. In a flared disk, the pair of planets should reach an asymptotic radius. This mechanism could account for the presence of Fomalhaut b ; then, a second, more massive planet, should be orbiting Fomalhaut at about 75 AU. △ Less

Submitted 6 October, 2009; originally announced October 2009.

Comments: 6 pages, 4 figures, accepted for publication by ApJ Letters

Journal ref: Astrophys.J.705:L148-L152,2009

On the horseshoe drag of a low-mass planet. I - Migration in isothermal disks

Authors: J. Casoli, F. S. Masset

Abstract: We investigate the unsaturated horseshoe drag exerted on a low-mass planet by an isothermal gaseous disk. In the globally isothermal case, we use a formal- ism, based on the use of a Bernoulli invariant, that takes into account pressure effects, and that extends the torque estimate to a region wider than the horse- shoe region. We find a result that is strictly identical to the standard horsesho… ▽ More We investigate the unsaturated horseshoe drag exerted on a low-mass planet by an isothermal gaseous disk. In the globally isothermal case, we use a formal- ism, based on the use of a Bernoulli invariant, that takes into account pressure effects, and that extends the torque estimate to a region wider than the horse- shoe region. We find a result that is strictly identical to the standard horseshoe drag. This shows that the horseshoe drag accounts for the torque of the whole corotation region, and not only of the horseshoe region, thereby deserving to be called corotation torque. We find that evanescent waves launched downstream of the horseshoe U-turns by the perturbations of vortensity exert a feed-back on the upstream region, that render the horseshoe region asymmetric. This asymmetry scales with the vortensity gradient and with the disk's aspect ratio. It does not depend on the planetary mass, and it does not have any impact on the horseshoe drag. Since the horseshoe drag has a steep dependence on the width of the horseshoe region, we provide an adequate definition of the width that needs to be used in horseshoe drag estimates. We then consider the case of locally isothermal disks, in which the tempera- ture is constant in time but depends on the distance to the star. The horseshoe drag appears to be different from the case of a globally isothermal disk. The difference, which is due to the driving of vortensity in the vicinity of the planet, is intimately linked to the topology of the flow. We provide a descriptive inter- pretation of these effects, as well as a crude estimate of the dependency of the excess on the temperature gradient. △ Less

Submitted 27 July, 2009; originally announced July 2009.

Comments: Accepted for publication in ApJ

On the horseshoe drag of a low-mass planet. II Migration in adiabatic disks

Authors: F. S. Masset, J. Casoli

Abstract: We evaluate the horseshoe drag exerted on a low-mass planet embedded in a gaseous disk, assuming the disk's flow in the coorbital region to be adiabatic. We restrict this analysis to the case of a planet on a circular orbit, and we assume a steady flow in the corotating frame. We also assume that the corotational flow upstream of the U-turns is unperturbed, so that we discard saturation effects.… ▽ More We evaluate the horseshoe drag exerted on a low-mass planet embedded in a gaseous disk, assuming the disk's flow in the coorbital region to be adiabatic. We restrict this analysis to the case of a planet on a circular orbit, and we assume a steady flow in the corotating frame. We also assume that the corotational flow upstream of the U-turns is unperturbed, so that we discard saturation effects. In addition to the classical expression for the horseshoe drag in barotropic disks, which features the vortensity gradient across corotation, we find an additional term which scales with the entropy gradient, and whose amplitude depends on the perturbed pressure at the stagnation point of the horseshoe separatrices. This additional torque is exerted by evanescent waves launched at the horseshoe separatrices, as a consequence of an asymmetry of the horseshoe region. It has a steep dependence on the potential's softening length, suggesting that the effect can be extremely strong in the three dimensional case. We describe the main properties of the coorbital region (the production of vortensity during the U-turns, the appearance of vorticity sheets at the downstream separatrices, and the pressure response), and we give torque expressions suitable to this regime of migration. Side results include a weak, negative feed back on migration, due to the dependence of the location of the stagnation point on the migration rate, and a mild enhancement of the vortensity related torque at large entropy gradient. △ Less

Submitted 27 July, 2009; originally announced July 2009.

Comments: Accepted for publication in ApJ

The dynamical role of the circumplanetary disc in planetary migration

Authors: A. Crida, C. Baruteau, W. Kley, F. Masset

Abstract: 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… ▽ More 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. △ Less

Submitted 4 June, 2009; originally announced June 2009.

Comments: 15 pages, 19 figures (most in color), in press in Astronomy and Astrophysics

Journal ref: A&A 502, 679-693 (2009)

Type I planetary migration in a self-gravitating disk

Authors: C. Baruteau, F. Masset

Abstract: We investigate the tidal interaction between a low-mass planet and a self-gravitating protoplanetary disk, by means of two-dimensional hydrodynamic simulations. We first show that considering a planet freely migrating in a disk without self-gravity leads to a significant overestimate of the migration rate. The overestimate can reach a factor of two for a disk having three times the surface densi… ▽ More We investigate the tidal interaction between a low-mass planet and a self-gravitating protoplanetary disk, by means of two-dimensional hydrodynamic simulations. We first show that considering a planet freely migrating in a disk without self-gravity leads to a significant overestimate of the migration rate. The overestimate can reach a factor of two for a disk having three times the surface density of the minimum mass solar nebula. Unbiased drift rates may be obtained only by considering a planet and a disk orbiting within the same gravitational potential. In a second part, the disk self-gravity is taken into account. We confirm that the disk gravity enhances the differential Lindblad torque with respect to the situation where neither the planet nor the disk feels the disk gravity. This enhancement only depends on the Toomre parameter at the planet location. It is typically one order of magnitude smaller than the spurious one induced by assuming a planet migrating in a disk without self-gravity. We confirm that the torque enhancement due to the disk gravity can be entirely accounted for by a shift of Lindblad resonances, and can be reproduced by the use of an anisotropic pressure tensor. We do not find any significant impact of the disk gravity on the corotation torque. △ Less

Submitted 28 January, 2008; originally announced January 2008.

Comments: 42 pages, 17 figures, accepted for publication in ApJ

Building Giant-Planet Cores at a Planet Trap

Authors: Alessandro Morbidelli, Aurelien Crida, Frederic Masset, Richard P. Nelson

Abstract: A well-known bottleneck for the core-accretion model of giant-planet formation is the loss of the cores into the star by Type-I migration, due to the tidal interactions with the gas disk. It has been shown that a steep surface-density gradient in the disk, such as the one expected at the boundary between an active and a dead zone, acts as a planet trap and prevents isolated cores from migrating… ▽ More A well-known bottleneck for the core-accretion model of giant-planet formation is the loss of the cores into the star by Type-I migration, due to the tidal interactions with the gas disk. It has been shown that a steep surface-density gradient in the disk, such as the one expected at the boundary between an active and a dead zone, acts as a planet trap and prevents isolated cores from migrating down to the central star. We study the relevance of the planet trap concept for the accretion and evolution of systems of multiple planetary embryos/cores. We performed hydrodynamical simulations of the evolution of systems of multiple massive objects in the vicinity of a planet trap. The planetary embryos evolve in 3 dimensions, whereas the disk is modeled with a 2D grid. Synthetic forces are applied onto the embryos to mimic the damping effect that the disk has on their inclinations. Systems with two embryos tend to acquire stable, separated and non-migrating orbits, with the more massive embryo placed at the planet trap and the lighter one farther out in the disk. Systems of multiple embryos are intrinsically unstable. Consequently, a long phase of mutual scattering can lead to accreting collisions among embryos; some embryos are injected into the inner part of the disk, where they can be evacuated into the star by Type I migration. The system can resume a stable, non-migrating configuration only when the number of surviving embryos decreases to a small value (~2-4). This can explain the limited number of giant planets in our solar system. These results should apply in general to any case in which the Type-I migration of the inner embryo is prevented by some mechanism, and not solely to the planet trap scenario. △ Less

Submitted 15 November, 2007; originally announced November 2007.

Comments: in press in Astronomy and Astrophysics

On the corotation torque in a radiatively inefficient disk

Authors: C. Baruteau, F. Masset

Abstract: We consider the angular momentum exchange at the corotation resonance between a two-dimensional gaseous disk and a uniformly rotating external potential, assuming that the disk flow is adiabatic. We first consider the linear case for an isolated resonance, for which we give an expression of the corotation torque that involves the pressure perturbation, and which reduces to the usual dependence o… ▽ More We consider the angular momentum exchange at the corotation resonance between a two-dimensional gaseous disk and a uniformly rotating external potential, assuming that the disk flow is adiabatic. We first consider the linear case for an isolated resonance, for which we give an expression of the corotation torque that involves the pressure perturbation, and which reduces to the usual dependence on the vortensity gradient in the limit of a cold disk. Although this expression requires the solution of the hydrodynamic equations, it provides some insight into the dynamics of the corotation region. In the general case, we find an additional dependence on the entropy gradient at corotation. This dependence is associated to the advection of entropy perturbations. These are not associated to pressure perturbations. They remain confined to the corotation region, where they yield a singular contribution to the corotation torque. In a second part, we check our torque expression by means of customized two-dimensional hydrodynamical simulations. In a third part, we contemplate the case of a planet embedded in a Keplerian disk, assumed to be adiabatic. We find an excess of corotation torque that scales with the entropy gradient, and we check that the contribution of the entropy perturbation to the torque is in agreement with the expression obtained from the linear analysis. We finally discuss some implications of the corotation torque expression for the migration of low mass planets in the regions of protoplanetary disks where the flow is radiatively inefficient on the timescale of the horseshoe U-turns. △ Less

Submitted 17 September, 2007; originally announced September 2007.

Comments: 36 pages, 19 figures, accepted for publication in ApJ

Simulating planet migration in globally evolving disks

Authors: A. Crida, A. Morbidelli, F. Masset

Abstract: Numerical simulations of planet-disk interactions are usually performed with hydro-codes that -- because they consider only an annulus of the disk, over a 2D grid -- can not take into account the global evolution of the disk. However, the latter governs planetary migration of type II, so that the accuracy of the planetary evolution can be questioned. To develop an algorithm that models the loc… ▽ More Numerical simulations of planet-disk interactions are usually performed with hydro-codes that -- because they consider only an annulus of the disk, over a 2D grid -- can not take into account the global evolution of the disk. However, the latter governs planetary migration of type II, so that the accuracy of the planetary evolution can be questioned. To develop an algorithm that models the local planet-disk interactions together with the global viscous evolution of the disk, we surround the usual 2D grid with a 1D grid ranging over the real extension of the disk. The 1D and 2D grids are coupled at their common boundaries via ghost rings, paying particular attention to the fluxes at the interface, especially the flux of angular momentum carried by waves. The computation is done in the frame centered on the center of mass to ensure angular momentum conservation. The global evolution of the disk and the local planet-disk interactions are both well described and the feedback of one on the other can be studied with this algorithm, for a negligible additional computing cost with respect to usual algorithms. △ Less

Submitted 30 August, 2006; originally announced August 2006.

Comments: 12 pages, 11 figures, accepted for publication in A&A

Journal ref: Astron.Astrophys.461:1173-1183,2007

On the migration of protogiant solid cores

Authors: F. S. Masset, G. D'Angelo, W. Kley

Abstract: 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 calcul… ▽ More 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. △ Less

Submitted 7 July, 2006; originally announced July 2006.

Comments: Accepted for publication in ApJ

Journal ref: Astrophys.J.652:730-745,2006

A comparative study of disc-planet interaction

Authors: M. de Val-Borro, R. G. Edgar, P. Artymowicz, P. Ciecielag, P. Cresswell, G. D'Angelo, E. J. Delgado-Donate, G. Dirksen, S. Fromang, A. Gawryszczak, H. Klahr, W. Kley, W. Lyra, F. Masset, G. Mellema, R. P. Nelson, S. -J. Paardekooper, A. Peplinski, A. Pierens, T. Plewa, K. Rice, C. Schaefer, R. Speith

Abstract: 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 seve… ▽ More 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. △ Less

Submitted 1 June, 2006; v1 submitted 9 May, 2006; originally announced May 2006.

Comments: 32 pages, accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.370:529-558,2006

Disk-Planet Interactions During Planet Formation

Authors: J. C. B. Papaloizou, R. P. Nelson, W. Kley, F. S. Masset, P. Artymowicz

Abstract: 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 clos… ▽ More 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. △ Less

Submitted 8 March, 2006; originally announced March 2006.

Comments: 14 pages, 14 figures, to appear in "Protostars and Planets V". A version with high resolution figures is available from http://www.maths.qmul.ac.uk/~rpn/preprints/index.html

Surface density jumps as planet traps

Authors: F. Masset, A. Morbidelli, A. Crida, J. Ferreira

Abstract: [Abridged] The tidal torque exerted by a protoplanetary disk with power law surface density and temperature profiles onto an embedded protoplanetary embryo is generally a negative quantity that leads to the embryo inwards migration. Here we investigate how the tidal torque balance is affected at a disk surface density radial jump. The jump has two consequences : - it affects the differential L… ▽ More [Abridged] The tidal torque exerted by a protoplanetary disk with power law surface density and temperature profiles onto an embedded protoplanetary embryo is generally a negative quantity that leads to the embryo inwards migration. Here we investigate how the tidal torque balance is affected at a disk surface density radial jump. The jump has two consequences : - it affects the differential Lindblad torque. In particular if the disk is merely empty on the inner side, the differential Lindblad torque almost amounts to the large negative outer Lindblad torque. - It affects the corotation torque, which is a quantity very sensitive to the local gradient of the disk surface density. In particular if the disk is depleted on the inside and if the jump occurs radially over a few pressure scale-heights, the corotation torque is a positive quantity that is much larger than in a power-law disk. We show by means of customized numerical simulations of low mass planets embedded in protoplanetary nebulae with a surface density jump that the second effect is dominant, that is that the corotation torque largely dominates the differential Lindblad torque on the edge of a central depletion. As a consequence the type I migration of low mass objects reaching the jump should be halted, and all these objects should be trapped there. △ Less

Submitted 10 January, 2006; originally announced January 2006.

Comments: accepted for publication in ApJ