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.

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.

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-draginstabilities may also decrease considerably when multiple species are considered.

The migration of planetary cores embedded in a protoplanetary disk is an important mechanism within planetformation 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.

Lyman α observations of the transiting exoplanet HD 209458b enable the study of exoplanets exospheres exposed to stellar EUV fluxes, as well as the interacting stellar wind properties. In this study we present 3D hydrodynamical models for the stellarplanetary wind interaction including radiation pressure and charge exchange, together with photoionization, recombination and collisional ionization processes. Our models explore the contribution of the radiation pressure and charge exchange on the Lyα absorption profile in a hydrodynamical framework, and for a single set of stellar wind parameters appropriate for HD 209458. We find that most of the absorption is produced by the material from the planet, with a secondary contribution of neutralized stellar ions by charge exchange. At the same time, the hydrodynamic shock heats up the planetary material, resulting in a broad thermal profile. Meanwhile, the radiation pressure yielded a small velocity shift of the absorbing material. While neither charge exchange nor radiation pressure provide enough neutrals at the velocity needed to explain the observations at −100 km s −1 individually, we find that the two effects combined with the broad thermal profile are able to explain the observations.