-
Dust rings trap protoplanets on eccentric orbits and get consumed by them
Authors:
David A. Velasco-Romero,
Frédéric S. Masset,
Alessandro Morbidelli,
Pablo Benítez-Llambay,
Leonardo Krapp,
Elena Lega
Abstract:
We study the orbital evolution and mass growth of protoplanets with masses $M \in [0.1-8]$~M$_\oplus$ in the vicinity of a dusty ring, using three-dimensional numerical simulations with a two-fluid model and nested-meshes. We find two stable, eccentric orbits that lock the planet in the ring vicinity, thereby inhibiting its migration and allowing it to accrete dust from the ring. One of these orbi…
▽ More
We study the orbital evolution and mass growth of protoplanets with masses $M \in [0.1-8]$~M$_\oplus$ in the vicinity of a dusty ring, using three-dimensional numerical simulations with a two-fluid model and nested-meshes. We find two stable, eccentric orbits that lock the planet in the ring vicinity, thereby inhibiting its migration and allowing it to accrete dust from the ring. One of these orbits has an eccentricity comparable to the aspect ratio of the gaseous disc and has its periastron within the ring, enabling intermittent accretion during each pass. The other orbit has a smaller eccentricity and an apoastron slightly inside the ring. A planet locked at the outer orbit efficiently accretes from the ring and can reach the critical mass for runaway gas accretion on timescales $\gtrsim 10^5$ yr (for a 10~M$_\oplus$ dust ring at 10~au) while a planet locked at the inner orbit has a slower growth and might not supersede the super-Earth stage over the disc lifetime. While in our runs a low-mass embryo forming within the ring eventually joins the outer orbit, it is likely that the path taken depends on the specific details of the ring. The trapping on the outer orbit arises from an intermittent, strong thermal force at each passage through the ring, where the accretion rate spikes. It is insensitive to uncertainties that plague models considering planets trapped on circular orbits in rings. It is highly robust and could allow a growing planet to follow an expanding ring over large distances.
△ Less
Submitted 23 July, 2024;
originally announced July 2024.
-
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 20 August, 2024; v1 submitted 22 February, 2024;
originally announced February 2024.
-
A Fast second-order solver for stiff multifluid dust and gas hydrodynamics
Authors:
Leonardo Krapp,
Juan Garrido-Deutelmoser,
Pablo Benítez-Llambay,
Kaitlin M. Kratter
Abstract:
We present MDIRK: a Multifluid second-order Diagonally-Implicit Runge-Kutta method to study momentum transfer between gas and an arbitrary number ($N$) of dust species. The method integrates the equations of hydrodynamics with an Implicit Explicit (IMEX) scheme and solves the stiff source term in the momentum equation with a diagonally-implicit asymptotically stable Runge-Kutta method (DIRK). In p…
▽ More
We present MDIRK: a Multifluid second-order Diagonally-Implicit Runge-Kutta method to study momentum transfer between gas and an arbitrary number ($N$) of dust species. The method integrates the equations of hydrodynamics with an Implicit Explicit (IMEX) scheme and solves the stiff source term in the momentum equation with a diagonally-implicit asymptotically stable Runge-Kutta method (DIRK). In particular, DIRK admits a simple analytical solution that can be evaluated with $\mathcal{O}(N)$ operations, instead of standard matrix inversion, which is $\mathcal{O}(N)^3$. Therefore the analytical solution significantly reduces the computational cost of the multifluid method, making it suitable for studying the dynamics of systems with particle-size distributions. We demonstrate that the method conserves momentum to machine precision and converges to the correct equilibrium solution with constant external acceleration. To validate our numerical method we present a series of simple hydrodynamic tests, including damping of sound waves, dusty shocks, a multi-fluid dusty Jeans instability, and a steady-state gas-dust drift calculation. The simplicity of MDIRK lays the groundwork to build fast high-order asymptotically stable multifluid methods.
△ Less
Submitted 26 February, 2024; v1 submitted 29 September, 2023;
originally announced October 2023.
-
RAM: Rapid Advection Algorithm on Arbitrary Meshes
Authors:
Pablo Benítez-Llambay,
Leonardo Krapp,
Ximena S. Ramos,
Kaitlin M. Kratter
Abstract:
The study of many astrophysical flows requires computational algorithms that can capture high Mach number flows, while resolving a large dynamic range in spatial and density scales. In this paper we present a novel method, RAM: Rapid Advection Algorithm on Arbitrary Meshes. RAM is a time-explicit method to solve the advection equation in problems with large bulk velocity on arbitrary computational…
▽ More
The study of many astrophysical flows requires computational algorithms that can capture high Mach number flows, while resolving a large dynamic range in spatial and density scales. In this paper we present a novel method, RAM: Rapid Advection Algorithm on Arbitrary Meshes. RAM is a time-explicit method to solve the advection equation in problems with large bulk velocity on arbitrary computational grids. In comparison with standard up-wind algorithms, RAM enables advection with larger time steps and lower truncation errors. Our method is based on the operator splitting technique and conservative interpolation. Depending on the bulk velocity and resolution, RAM can decrease the numerical cost of hydrodynamics by more than one order of magnitude. To quantify the truncation errors and speed-up with RAM, we perform one and two-dimensional hydrodynamics tests. We find that the order of our method is given by the order of the conservative interpolation and that the effective speed up is in agreement with the relative increment in time step. RAM will be especially useful for numerical studies of disk-satellite interaction, characterized by high bulk orbital velocities, and non-trivial geometries. Our method dramatically lowers the computational cost of simulations that simultaneously resolve the global disk and well inside the Hill radius of the secondary companion.
△ Less
Submitted 9 May, 2023;
originally announced May 2023.
-
Quantifying the Impact of the Dust Torque on the Migration of Low-mass Planets
Authors:
Octavio M. Guilera,
Pablo Benitez-Llambay,
Marcelo M. Miller Bertolami,
Martin E. Pessah
Abstract:
Disk solids are critical in many planet formation processes, however, their effect on planet migration remains largely unexplored. Here we assess for the first time this important issue by building on the systematic measurements of dust torques on an embedded planet by Benitez-Llambay & Pessah (2018). Adopting standard models for the gaseous disk and its solid content, we quantify the impact of th…
▽ More
Disk solids are critical in many planet formation processes, however, their effect on planet migration remains largely unexplored. Here we assess for the first time this important issue by building on the systematic measurements of dust torques on an embedded planet by Benitez-Llambay & Pessah (2018). Adopting standard models for the gaseous disk and its solid content, we quantify the impact of the dust torque for a wide range of conditions describing the disk/planet system. We show that the total torque can be positive and revert inward planet migration for planetary cores with $M_{\rm p} \lesssim 10 M_\oplus$. We compute formation tracks for low-mass embryos for conditions usually invoked when modeling planet formation processes. Our most important conclusion is that dust torques can have a significant impact on the migration and formation history of planetary embryos. The most important implications of our findings are: $\it{i})$ For nominal dust-to-gas mass ratios $ε\simeq 0.01$, low-mass planets migrate outwards beyond the water ice-line if most of the mass in solids is in particles with Stokes numbers St $\simeq 0.1$. $\it{ii})$. For $ε\gtrsim 0.02-0.05$, solids with small Stokes numbers, St $\simeq 0.01$, can play a dominant role if most of the mass is in those particles. $\it{iii})$ Dust torques have the potential to enable low-mass planetary cores formed in the inner disk to migrate outwards and act as the seed for massive planets at distances of tens of au.
△ Less
Submitted 3 May, 2023;
originally announced May 2023.
-
Hydro-, Magnetohydro-, and Dust-Gas Dynamics of Protoplanetary Disks
Authors:
G. Lesur,
B. Ercolano,
M. Flock,
M. -K. Lin,
C. -C. Yang,
J. A. Barranco,
P. Benitez-Llambay,
J. Goodman,
A. Johansen,
H. Klahr,
G. Laibe,
W. Lyra,
P. Marcus,
R. P. Nelson,
J. Squire,
J. B. Simon,
N. Turner,
O. M. Umurhan,
A. N. Youdin
Abstract:
The building of planetary systems is controlled by the gas and dust dynamics of protoplanetary disks. While the gas is simultaneously accreted onto the central star and dissipated away by winds, dust grains aggregate and collapse to form planetesimals and eventually planets. This dust and gas dynamics involves instabilities, turbulence and complex non-linear interactions which ultimately control t…
▽ More
The building of planetary systems is controlled by the gas and dust dynamics of protoplanetary disks. While the gas is simultaneously accreted onto the central star and dissipated away by winds, dust grains aggregate and collapse to form planetesimals and eventually planets. This dust and gas dynamics involves instabilities, turbulence and complex non-linear interactions which ultimately control the observational appearance and the secular evolution of these disks.
This chapter is dedicated to the most recent developments in our understanding of the dynamics of gaseous and dusty disks, covering hydrodynamic and magnetohydrodynamic turbulence, gas-dust instabilities, dust clumping and disk winds. We show how these physical processes have been tested from observations and highlight standing questions that should be addressed in the future.
△ Less
Submitted 18 March, 2022;
originally announced March 2022.
-
Dust Settling Instability in Protoplanetary Discs
Authors:
Leonardo Krapp,
Andrew N. Youdin,
Kaitlin M. Kratter,
Pablo Benítez-Llambay
Abstract:
The streaming instability (SI) has been extensively studied in the linear and non-linear regimes as a mechanism to concentrate solids and trigger planetesimal formation in the midplane of protoplanetary discs. A related dust settling instability (DSI) applies to particles while settling towards the midplane. The DSI has previously been studied in the linear regime, with predictions that it could t…
▽ More
The streaming instability (SI) has been extensively studied in the linear and non-linear regimes as a mechanism to concentrate solids and trigger planetesimal formation in the midplane of protoplanetary discs. A related dust settling instability (DSI) applies to particles while settling towards the midplane. The DSI has previously been studied in the linear regime, with predictions that it could trigger particle clumping away from the midplane. This work presents a range of linear calculations and non-linear simulations, performed with FARGO3D, to assess conditions for DSI growth. We expand on previous linear analyses by including particle size distributions and performing a detailed study of the amount of background turbulence needed to stabilize the DSI. When including binned size distributions, the DSI often produces converged growth rates with fewer bins than the standard SI. With background turbulence, we find that the most favorable conditions for DSI growth are weak turbulence, characterized by $α\lesssim 10^{-6}$ with intermediate-sized grains that settle from one gas scale-height. These conditions could arise during a sudden decrease in disc turbulence following an accretion outburst. Ignoring background turbulence, we performed a parameter survey of local 2D DSI simulations. Particle clumping was either weak or occurred slower than particles settle. Clumping was reduced by a factor of two in a comparison 3D simulation. Overall, our results strongly disfavor the hypothesis that the DSI significantly promotes planetesimal formation. Non-linear simulations of the DSI with different numerical methods could support or challenge these findings.
△ Less
Submitted 15 September, 2020; v1 submitted 9 April, 2020;
originally announced April 2020.
-
Low-mass planet migration in three dimensional wind-driven inviscid discs: a negative corotation torque
Authors:
Colin P. McNally,
Richard P. Nelson,
Sijme-Jan Paardekooper,
Pablo Benitez-Llambay,
Oliver Gressel
Abstract:
We present simulations of low-mass planet--disc interactions in inviscid three-dimensional discs. We show that a wind-driven laminar accretion flow through the surface layers of the disc does not significantly modify the migration torque experienced by embedded planets. More importantly, we find that 3D effects lead to a dramatic change in the behaviour of the dynamical corotation torque compared…
▽ More
We present simulations of low-mass planet--disc interactions in inviscid three-dimensional discs. We show that a wind-driven laminar accretion flow through the surface layers of the disc does not significantly modify the migration torque experienced by embedded planets. More importantly, we find that 3D effects lead to a dramatic change in the behaviour of the dynamical corotation torque compared to earlier 2D theory and simulations. Although it was previously shown that the dynamical corotation torque could act to slow and essentially stall the inward migration of a low-mass planet, our results in 3D show that the dynamical corotation torque has the complete opposite effect and speeds up inward migration. Our numerical experiments implicate buoyancy resonances as the cause. These have two effects: (i) they exert a direct torque on the planet, whose magnitude relative to the Lindblad torque is measured in our simulations to be small; (ii) they torque the gas librating on horseshoe orbits in the corotation region and drive evolution of its vortensity, leading to the negative dynamical corotation torque. This indicates that at low turbulent viscosity, the detailed vertical thermal structure of the protoplanetary disc plays an important role in determining the migration behaviour of embedded planets. If this result holds up under a more refined treatment of disc thermal evolution, then it has important implications for understanding the formation and early evolution of planetary systems.
△ Less
Submitted 13 March, 2020; v1 submitted 25 February, 2020;
originally announced February 2020.
-
Predicting the observational signature of migrating Neptune-sized planets in low-viscosity disks
Authors:
Philipp Weber,
Sebastián Pérez,
Pablo Benítez-Llambay,
Oliver Gressel,
Simon Casassus,
Leonardo Krapp
Abstract:
The migration of planetary cores embedded in a protoplanetary disk is an important mechanism within planet-formation theory, relevant for the architecture of planetary systems. Consequently, planet migration is actively discussed, yet often results of independent theoretical or numerical studies are unconstrained due to the lack of observational diagnostics designed in light of planet migration. I…
▽ More
The migration of planetary cores embedded in a protoplanetary disk is an important mechanism within planet-formation theory, relevant for the architecture of planetary systems. Consequently, planet migration is actively discussed, yet often results of independent theoretical or numerical studies are unconstrained due to the lack of observational diagnostics designed in light of planet migration. In this work we follow the idea of inferring the migration behavior of embedded planets by means of the characteristic radial structures that they imprint in the disk's dust density distribution. We run hydrodynamical multifluid simulations of gas and several dust species in a locally isothermal $α$-disk in the low-viscosity regime ($α=10^{-5}$) and investigate the obtained dust structures. In this framework, a planet of roughly Neptune mass can create three (or more) rings in which dust accumulates. We find that the relative spacing of these rings depends on the planet's migration speed and direction. By performing subsequent radiative transfer calculations and image synthesis we show that - always under the condition of a near-inviscid disk - different migration scenarios are, in principle, distinguishable by long-baseline, state-of-the-art ALMA observations.
△ Less
Submitted 28 October, 2019; v1 submitted 4 September, 2019;
originally announced September 2019.
-
Long Baseline Observations of the HD100546 Protoplanetary Disk with ALMA
Authors:
Sebastián Pérez,
Simon Casassus,
Antonio Hales,
Sebastián Marino,
Anthony Cheetham,
Alice Zurlo,
Lucas Cieza,
Ruobing Dong,
Felipe Alarcón,
Pablo Benítez-Llambay,
Ed Fomalont,
Henning Avenhaus
Abstract:
Using the Atacama Large Millimeter/submillimeter Array (ALMA), we observed the young Herbig star HD 100546, host to a prominent disk with a deep, wide gap in the dust. The high-resolution 1.3 mm continuum observation reveals fine radial and azimuthal substructures in the form of a complex maze of ridges and trenches sculpting a dust ring. The $^{12}$CO(2-1) channel maps are modulated by wiggles or…
▽ More
Using the Atacama Large Millimeter/submillimeter Array (ALMA), we observed the young Herbig star HD 100546, host to a prominent disk with a deep, wide gap in the dust. The high-resolution 1.3 mm continuum observation reveals fine radial and azimuthal substructures in the form of a complex maze of ridges and trenches sculpting a dust ring. The $^{12}$CO(2-1) channel maps are modulated by wiggles or kinks that deviate from Keplerian kinematics particularly over the continuum ring, where deviations span 90$^\circ$ in azimuth, covering 5 km s$^{-1}$. The most pronounced wiggle resembles the imprint of an embedded massive planet of at least 5 M$_{\rm Jup}$ predicted from previous hydrodynamical simulations (Perez, Casassus, & Benitez-Llambay 2018). Such planet is expected to open a deep gap in both gas and dust density fields within a few orbital timescales, yet the kinematic wiggles lie near ridges in the continuum. The lesser strength of the wiggles in the $^{13}$CO and C$^{18}$O isotopologues show that the kinematic signature weakens at lower disk heights, and suggests qualitatively that it is due to vertical flows in the disk surface. Within the gap, the velocity field transitions from Keplerian to strongly non-Keplerian via a twist in position angle, suggesting the presence of another perturber and/or an inner warp. We also present VLT/SPHERE sparse aperture masking data which recovers scattered light emission from the gap's edges but shows no evidence for signal within the gap, discarding a stellar binary origin for its opening.
△ Less
Submitted 16 January, 2020; v1 submitted 14 June, 2019;
originally announced June 2019.
-
Streaming Instability for Particle-Size Distributions
Authors:
Leonardo Krapp,
Pablo Benítez-Llambay,
Oliver Gressel,
Martin E. Pessah
Abstract:
The streaming instability is thought to play a central role in the early stages of planet formation by enabling the efficient bypass of a number of barriers hindering the formation of planetesimals. We present the first study exploring the efficiency of the linear streaming instability when a particle-size distribution is considered. We find that, for a given dust-to-gas mass ratio, the multi-spec…
▽ More
The streaming instability is thought to play a central role in the early stages of planet formation by enabling the efficient bypass of a number of barriers hindering the formation of planetesimals. We present the first study exploring the efficiency of the linear streaming instability when a particle-size distribution is considered. We find that, for a given dust-to-gas mass ratio, the multi-species streaming instability grows on timescales much longer than those expected when only one dust species is involved. In particular, distributions that contain close-to-order-unity dust-to-gas mass ratios lead to unstable modes that can grow on timescales comparable, or larger, with those of secular instabilities. We anticipate that processes leading to particle segregation and/or concentration can create favourable conditions for the instability to grow fast. Our findings may have important implications for a large number of processes in protoplanetary disks that rely on the streaming instability as usually envisioned for a unique dust species. Our results suggest that the growth rates of other resonant-drag-instabilities may also decrease considerably when multiple species are considered.
△ Less
Submitted 10 June, 2019; v1 submitted 30 May, 2019;
originally announced May 2019.
-
Probing the Protosolar Disk Using Dust Filtering at Gaps in the Early Solar System
Authors:
Troels Haugbølle,
Philipp Weber,
Daniel P. Wielandt,
Pablo Benítez-Llambay,
Martin Bizzarro,
Oliver Gressel,
Martin E. Pessah
Abstract:
Jupiter and Saturn formed early, before the gas disk dispersed. The presence of gap-opening planets affects the dynamics of the gas and embedded solids and halts the inward drift of grains above a certain size. A drift barrier can explain the absence of calcium aluminium rich inclusions (CAIs) in chondrites originating from parent bodies that accreted in the inner solar system. Employing an interd…
▽ More
Jupiter and Saturn formed early, before the gas disk dispersed. The presence of gap-opening planets affects the dynamics of the gas and embedded solids and halts the inward drift of grains above a certain size. A drift barrier can explain the absence of calcium aluminium rich inclusions (CAIs) in chondrites originating from parent bodies that accreted in the inner solar system. Employing an interdisciplinary approach, we use a $μ$-X-Ray-fluorescence scanner to search for large CAIs and a scanning electron microscope to search for small CAIs in the ordinary chondrite NWA 5697. We carry out long-term, two-dimensional simulations including gas, dust, and planets to characterize the transport of grains within the viscous $α$-disk framework exploring the scenarios of a stand-alone Jupiter, Jupiter and Saturn \textit{in situ}, or Jupiter and Saturn in a 3:2 resonance. In each case, we find a critical grain size above which drift is halted as a function of the physical conditions in the disk. From the laboratory search we find four CAIs with a largest size of $\approx$200$\,μ$m. \Combining models and data, we provide an estimate for the upper limit of the $α$-viscosity and the surface density at the location of Jupiter, using reasonable assumptions about the stellar accretion rate during inward transport of CAIs, and assuming angular momentum transport to happen exclusively through viscous effects. Moreover, we find that the compound gap structure in the presence of Saturn in a 3:2 resonance favors inward transport of grains larger than CAIs currently detected in ordinary chondrites.
△ Less
Submitted 17 July, 2019; v1 submitted 28 March, 2019;
originally announced March 2019.
-
Migrating super-Earths in low-viscosity discs: unveiling the roles of feedback, vortices, and laminar accretion flows
Authors:
Colin P. McNally,
Richard P. Nelson,
Sijme-Jan Paardekooper,
Pablo Benítez-Llambay
Abstract:
We present the highest resolution study to date of super-Earths migrating in inviscid and low-viscosity discs, motivated by the connection to laminar, wind-driven models of protoplanetary discs. Our models unveil the critical role of vortices in determining the migration behaviour for partial gap-opening planets. Vortices form in pressure maxima at gap edges, and prevent the disc-feedback stopping…
▽ More
We present the highest resolution study to date of super-Earths migrating in inviscid and low-viscosity discs, motivated by the connection to laminar, wind-driven models of protoplanetary discs. Our models unveil the critical role of vortices in determining the migration behaviour for partial gap-opening planets. Vortices form in pressure maxima at gap edges, and prevent the disc-feedback stopping of migration for intermediate planets in low-viscosity and inviscid discs, contrary to the concept of the `inertial limit' or `disc feedback' halting predicted from analytical models. Vortices may also form in the corotation region, and can dramatically modify migration behaviour through direct gravitational interaction with the planet. These features become apparent at high resolution, and for all but the highest viscosities there exist significant difficulties in obtaining numerically converged results. The migration of partial gap-opening planets, however, clearly becomes chaotic for sufficiently low viscosities. At moderate viscosity, a smooth disc-feedback regime is found in which migration can slow substantially, and the migration time-scale observed corresponds to migration being driven by diffusive relaxation of the gap edges. At high viscosity classical Type I migration is recovered. For Jupiter-analogue planets in inviscid discs, a wide, deep gap is formed. Transient Type II migration occurs over radial length-scales corresponding to the gap width, beyond which migration can stall. Finally, we examine the particle trapping driven by structures left in inviscid discs by a migrating planet, and find that particle traps in the form of multiple rings and vortices can persist long after the planet has passed. In this case, the observation of particle traps by submillimetre interferometers such as ALMA cannot be used to infer the current presence of an adjacent planet.
△ Less
Submitted 13 January, 2019; v1 submitted 30 November, 2018;
originally announced November 2018.
-
Asymptotically stable numerical method for multispecies momentum transfer: gas and multifluid-dust test suite and implementation in FARGO3D
Authors:
Pablo Benítez-Llambay,
Leonardo Krapp,
Martín E. Pessah
Abstract:
We present an asymptotically and unconditionally stable numerical method to account for the momentum transfer between multiple species. Momentum is conserved to machine precision. This implies that the asymptotic equilibrium corresponds to the velocity of the center of mass. Aimed at studying dust dynamics, we implement this numerical method in the publicly available code FARGO3D. To validate our…
▽ More
We present an asymptotically and unconditionally stable numerical method to account for the momentum transfer between multiple species. Momentum is conserved to machine precision. This implies that the asymptotic equilibrium corresponds to the velocity of the center of mass. Aimed at studying dust dynamics, we implement this numerical method in the publicly available code FARGO3D. To validate our implementation, we develop a test suite for an arbitrary number of species, based on analytical or exact solutions of problems related to perfect damping, damped sound waves, shocks, local and global gas-dust radial drift in a disk and, linear streaming instability. In particular, we obtain first-order, steady-state solutions for the radial drift of multiple dust species in protoplanetary disks, in which the pressure gradient is not necessarily small. We additionally present non-linear shearing-box simulations of the streaming instability and compare them with previous results obtained with Lagrangian particles. We successfully validate our implementation by recovering the solutions from the test suite to second- and first-order accuracy in space and time, respectively. From this, we conclude that our scheme is suitable, and very robust, to study the self-consistent dynamics of several fluids. In particular, it can be used for solving the collisions between gas and dust in protoplanetary disks, with any degree of coupling.
△ Less
Submitted 10 June, 2019; v1 submitted 19 November, 2018;
originally announced November 2018.
-
Dust segregation in Hall-dominated turbulent protoplanetary disks
Authors:
Leonardo Krapp,
Oliver Gressel,
Pablo Benítez-Llambay,
Turlough P. Downes,
Gopakumar Mohandas,
Martin E. Pessah
Abstract:
Imaging of the dust continuum emitted from disks around nearby protostars reveals diverse substructure. In recent years, theoretical efforts have been intensified to investigate how far the intrinsic dynamics of protoplanetary disks (PPDs) can lead to such features. Turbulence in the realm of non-ideal magnetohydrodynamics (MHD) is one candidate for explaining the generation of zonal flows which c…
▽ More
Imaging of the dust continuum emitted from disks around nearby protostars reveals diverse substructure. In recent years, theoretical efforts have been intensified to investigate how far the intrinsic dynamics of protoplanetary disks (PPDs) can lead to such features. Turbulence in the realm of non-ideal magnetohydrodynamics (MHD) is one candidate for explaining the generation of zonal flows which can lead to local dust enhancements. Adopting a radially varying cylindrical disk model, and considering combinations of vertical and azimuthal initial net flux, we perform 3D non-ideal MHD simulations aimed at studying self-organization induced by the Hall effect in turbulent PPDs. To this end, new modules have been incorporated into the NIRVANA-III and FARGO3D MHD codes. We moreover include dust grains, treated in the fluid approximation, in order to study their evolution subject to the emerging zonal flows. In the regime of a dominant Hall effect, we robustly obtain large-scale organized concentrations in the vertical magnetic field that remain stable for hundreds of orbits. For disks with vertical initial net flux alone, we confirm the presence of zonal flows and vortices that introduce regions of super-Keplerian gas flow. Including a moderately strong net-azimuthal magnetic flux can significantly alter the dynamics, partially preventing the self-organization of zonal flows. For plasma beta-parameters smaller than 50, large-scale, near-axisymmetric structures develop in the vertical magnetic flux. In all cases, we demonstrate that the emerging features are capable of accumulating dust grains for a range of Stokes numbers.
△ Less
Submitted 30 May, 2019; v1 submitted 23 August, 2018;
originally announced August 2018.
-
Observability of planet-disc interactions in CO kinematics
Authors:
Sebastian Perez,
Simon Casassus,
Pablo Benítez-Llambay
Abstract:
Empirical evidence of planets in gas-rich circumstellar discs is required to constrain giant planet formation theories. Here we study the kinematic patterns which arise from planet-disc interactions and their observability in CO rotational emission lines. We perform three-dimensional hydrodynamical simulations of single giant planets, and predict the emergent intensity field with radiative transfe…
▽ More
Empirical evidence of planets in gas-rich circumstellar discs is required to constrain giant planet formation theories. Here we study the kinematic patterns which arise from planet-disc interactions and their observability in CO rotational emission lines. We perform three-dimensional hydrodynamical simulations of single giant planets, and predict the emergent intensity field with radiative transfer. Pressure gradients at planet-carved gaps, spiral wakes and vortices bear strong kinematic counterparts. The iso-velocity contours in the CO(2-1) line centroids $v_\circ$ reveal large-scale perturbations, corresponding to abrupt transitions from below sub-Keplerian to super-Keplerian rotation along with radial and vertical flows. The increase in line optical depth at the edge of the gap also modulates $v_\circ$, but this is a mild effect compared to the dynamical imprint of the planet-disc interaction. The large-scale deviations from the Keplerian rotation thus allow the planets to be indirectly detected via the first moment maps of molecular gas tracers, at ALMA angular resolutions. The strength of these deviations depends on the mass of the perturber. This initial study paves the way to eventually determine the mass of the planet by comparison with more detailed models.
△ Less
Submitted 13 June, 2018;
originally announced June 2018.
-
Orbital migration and Resonance Offset of the Kepler-25 and K2-24 systems
Authors:
C. Charalambous,
X. S. Ramos,
P. Benítez-Llambay,
C. Beaugé
Abstract:
Based on the model described in Ramos et al., 2017, we present an analytical+numerical study of the resonance capture under Type-I migration for the Kepler-25 (Marcy et al., 2014) and K2-24 (Petigura et al., 2016) Kepler systems, both close to a 2/1 mean-motion resonance. We find that, depending on the flare index and the proximity to the central star, the average value of the period-ratio between…
▽ More
Based on the model described in Ramos et al., 2017, we present an analytical+numerical study of the resonance capture under Type-I migration for the Kepler-25 (Marcy et al., 2014) and K2-24 (Petigura et al., 2016) Kepler systems, both close to a 2/1 mean-motion resonance. We find that, depending on the flare index and the proximity to the central star, the average value of the period-ratio between two consecutive planets show a significant deviation with respect to the resonant nominal value, up to values well in agreement with the observations.
△ Less
Submitted 15 March, 2018;
originally announced March 2018.
-
Characterizing the variable dust permeability of planet-induced gaps
Authors:
Philipp Weber,
Pablo Benítez-Llambay,
Oliver Gressel,
Leonardo Krapp,
Martin E. Pessah
Abstract:
Aerodynamic theory predicts that dust grains in protoplanetary disks will drift radially inward on comparatively short timescales. In this context, it has long been known that the presence of a gap opened by a planet can alter the dust dynamics significantly. In this paper, we carry out a systematic study employing long-term numerical simulations aimed at characterizing the critical particle-size…
▽ More
Aerodynamic theory predicts that dust grains in protoplanetary disks will drift radially inward on comparatively short timescales. In this context, it has long been known that the presence of a gap opened by a planet can alter the dust dynamics significantly. In this paper, we carry out a systematic study employing long-term numerical simulations aimed at characterizing the critical particle-size for retention outside a gap as a function of particle size and for various key parameters defining the protoplanetary disk model. To this end, we perform multifluid hydrodynamical simulations in two dimensions, including different dust species, which we treat as pressureless fluids. We initialize the dust outside of the planet's orbit and study under which conditions dust grains are able to cross the gap carved by the planet. In agreement with previous work, we find that the permeability of the gap depends both on dust dynamical properties and the gas disk structure: while small dust follows the viscously accreting gas through the gap, dust grains approaching a critical size are progressively filtered out. Moreover, we introduce and compute a depletion factor that enables us to quantify the way in which higher viscosity, smaller planet mass, or a more massive disk can shift this critical size to larger values. Our results indicate that gap-opening planets may act to deplete the inner reaches of protoplanetary disks of large dust grains -- potentially limiting the accretion of solids onto forming terrestrial planets.
△ Less
Submitted 25 January, 2018; v1 submitted 24 January, 2018;
originally announced January 2018.
-
Torques Induced by Scattered Pebble-flow in Protoplanetary Disks
Authors:
Pablo Benítez-Llambay,
Martin E. Pessah
Abstract:
Fast inward migration of planetary cores is a common problem in the current planet formation paradigm. Even though dust is ubiquitous in protoplanetary disks, its dynamical role in the migration history of planetary embryos has not been assessed. In this Letter, we show that the scattered pebble-flow induced by a low-mass planetary embryo leads to an asymmetric dust-density distribution that is ab…
▽ More
Fast inward migration of planetary cores is a common problem in the current planet formation paradigm. Even though dust is ubiquitous in protoplanetary disks, its dynamical role in the migration history of planetary embryos has not been assessed. In this Letter, we show that the scattered pebble-flow induced by a low-mass planetary embryo leads to an asymmetric dust-density distribution that is able to exert a net torque. By analyzing a large suite of multifluid hydrodynamical simulations addressing the interaction between the disk and a low-mass planet on a fixed circular orbit, and neglecting dust feedback onto the gas, we identify two different regimes, gas- and gravity-dominated, where the scattered pebble-flow results in almost all cases in positive torques. We collect our measurements in a first torque map for dusty disks, which will enable the incorporation of the effect of dust dynamics on migration into population synthesis models. Depending on the dust drift speed, the dust-to-gas mass ratio/distribution and the embryo mass, the dust-induced torque has the potential to halt inward migration or even induce fast outward migration of planetary cores. We thus anticipate that dust-driven migration could play a dominant role during the formation history of planets. Because dust torques scale with disk metallicity, we propose that dust-driven outward migration may enhance the occurrence of distant giant planets in higher-metallicity systems.
△ Less
Submitted 22 November, 2018; v1 submitted 24 January, 2018;
originally announced January 2018.
-
Planetary migration and the origin of the 2:1 and 3:2 (near)-resonant population of close-in exoplanets
Authors:
X. S. Ramos,
C. Charalambous,
P. Benítez-Llambay,
C. Beaugé
Abstract:
We present an analytical and numerical study of the orbital migration and resonance capture of fictitious two-planet systems with masses in the super-Earth range undergoing Type-I migration. We find that, depending on the flare index and proximity to the central star, the average value of the period ratio, $P_2/P_1$, between both planets may show a significant deviation with respect to the nominal…
▽ More
We present an analytical and numerical study of the orbital migration and resonance capture of fictitious two-planet systems with masses in the super-Earth range undergoing Type-I migration. We find that, depending on the flare index and proximity to the central star, the average value of the period ratio, $P_2/P_1$, between both planets may show a significant deviation with respect to the nominal value. For planets trapped in the 2:1 commensurability, offsets may reach values on the order of $0.1$ for orbital periods on the order of $1$ day, while systems in the 3:2 mean-motion resonance (MMR) show much smaller offsets for all values of the semimajor axis. These properties are in good agreement with the observed distribution of near-resonant exoplanets, independent of their detection method. We show that 2:1-resonant systems far from the star, such as HD82943 and HR8799, are characterized by very small resonant offsets, while higher values are typical of systems discovered by Kepler with orbital periods approximately a few days. Conversely, planetary systems in the vicinity of the 3:2 MMR show little offset with no significant dependence on the orbital distance. In conclusion, our results indicate that the distribution of Kepler planetary systems around the 2:1 and 3:2 MMR are consistent with resonant configurations obtained as a consequence of a smooth migration in a laminar flared disk, and no external forces are required to induce the observed offset or its dependence with the commensurability or orbital distance from the star.
△ Less
Submitted 21 April, 2017;
originally announced April 2017.
-
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.
-
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.
-
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.
-
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.
-
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.
-
Origin and Detectability of coorbital planets from radial velocity data
Authors:
C. A. Giuppone,
P. Benitez-Llambay,
C. Beauge
Abstract:
We analyze the possibilities of detection of hypothetical exoplanets in coorbital motion from synthetic radial velocity (RV) signals, taking into account different types of stable planar configurations, orbital eccentricities and mass ratios. For each nominal solution corresponding to small-amplitude oscillations around the periodic solution, we generate a series of synthetic RV curves mimicking t…
▽ More
We analyze the possibilities of detection of hypothetical exoplanets in coorbital motion from synthetic radial velocity (RV) signals, taking into account different types of stable planar configurations, orbital eccentricities and mass ratios. For each nominal solution corresponding to small-amplitude oscillations around the periodic solution, we generate a series of synthetic RV curves mimicking the stellar motion around the barycenter of the system. We then fit the data sets obtained assuming three possible different orbital architectures: (a) two planets in coorbital motion, (b) two planets in a 2/1 mean-motion resonance, and (c) a single planet. We compare the resulting residuals and the estimated orbital parameters.
For synthetic data sets covering only a few orbital periods, we find that the discrete radial velocity signal generated by a coorbital configuration could be easily confused with other configurations/systems, and in many cases the best orbital fit corresponds to either a single planet or two bodies in a 2/1 resonance. However, most of the incorrect identifications are associated to dynamically unstable solutions.
We also compare the orbital parameters obtained with two different fitting strategies: a simultaneous fit of two planets and a nested multi-Keplerian model. We find that the nested models can yield incorrect orbital configurations (sometimes close to fictitious mean-motion resonances) that are nevertheless dynamically stable and with orbital eccentricities lower than the correct nominal solutions.
Finally, we discuss plausible mechanisms for the formation of coorbital configurations, by the interaction between two giant planets and an inner cavity in the gas disk. For equal mass planets, both Lagrangian and anti-Lagrangian configurations can be obtained from same initial condition depending on final time of integration.
△ Less
Submitted 16 March, 2012; v1 submitted 1 December, 2011;
originally announced December 2011.
-
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.