Source author record

Pablo Benítez-Llambay

Pablo Benítez-Llambay appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

6works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

6 published item(s)

preprint2020arXiv

Dust Settling Instability in Protoplanetary Discs

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.

preprint2020arXiv

Long Baseline Observations of the HD100546 Protoplanetary Disk with ALMA

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.

preprint2016arXiv

FARGO3D: A new GPU-oriented MHD code

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.

preprint2016arXiv

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

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.

preprint2015arXiv

Horseshoe drag in three-dimensional globally isothermal disks

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.

preprint2015arXiv

Planet heating prevents inward migration of planetary cores

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.