On the Origin of Dust Structures in Protoplanetary Disks: Constraints from the Rossby Wave Instability

High resolution sub-mm observations of protoplanetary disks with ALMA have revealed that dust rings are common in large, bright disks. The leading explanation for these structures is dust-trapping in a local gas pressure maximum, caused by an embedded planet or other dynamical process. Independent of origin, such dust traps should be stable for many orbits to collect significant dust. However, ring-like perturbations in gas disks are also known to trigger the Rossby Wave Instability (RWI). We investigate whether axisymmetric pressure bumps can simultaneously trap dust and remain stable to the RWI. The answer depends on the thermodynamic properties of pressure bumps. For isothermal bumps, dust traps are RWI-stable for widths from ${\sim}1$ to several gas scale-heights. Adiabatic dust traps are stable over a smaller range of widths. For temperature bumps with no surface density component, however, all dust traps tend to be unstable. Smaller values of disk aspect ratio allow stable dust trapping at lower bump amplitudes and over a larger range of widths. We also report a new approximate criterion for RWI. Instability occurs when the radial oscillation frequency is $\lesssim75$\% of the Keplerian frequency, which differs from the well-known Lovelace necessary (but not sufficient) criterion for instability. Our results can guide ALMA observations of molecular gas by constraining the resolution and sensitivity needed to identify the pressure bumps thought to be responsible for dust rings.Comment: 11 pages, 5 figures, accepted to ApJ

Efectos magnetohidrodinámicos no ideales en discos protoplanetarios

En este trabajo realizamos un estudio de los efectos no ideales de la MHD, conocidos como efecto Hall, difusión Ambipolar y difusión de Ohm. Éstos son capaces de suprimir la inestabilidad Magneto Rotacional presente en discos de acreción magnetizados y generar nuevos mecanismos de acreción. Para el estudio de estos efectos es necesario el uso de códigos MHD, a fin de resolver la compleja dinámica involucrada. Así, el principal objetivo de este trabajo es la utilización, explotación y extensión del código público magnetohidrodinámico FARGO3D mediante el desarrollo de un nuevo módulo físico que incluya los efectos no ideales de la MHD conocidos como efecto Hall y difusión Ambipolar. Presentamos aquí una serie de simulaciones que demuestran la efectiva implementación numérica de estos módulos. Esta implementación abre una puerta hacia el estudio de nuevos efectos en la interacción planeta-disco

Asymptotically Stable Numerical Method for Multispecies Momentum Transfer:Gas and Multifluid Dust Test Suite and Implementation in FARGO3D

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.Comment: 27 pages, 12 figures; accepted for publication in ApJ

A Fast second-order solver for stiff multifluid dust and gas hydrodynamics

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.Comment: 21 pages, 7 figures, accepted for publication ApJ

RAM: Rapid Advection Algorithm on Arbitrary Meshes

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.Comment: 15 pages, 7 figures. Submitted to ApJ. Comments are welcom

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 $\alpha \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

Characterizing the variable dust permeability of planet-induced gaps

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.Comment: 18 pages, 14 figures, accepted for publication in Ap

Predicting the observational signature of migrating Neptune-sized planets in low-viscosity disks

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 $\alpha$-disk in the low-viscosity regime ($\alpha=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.Comment: 19 pages, 13 figures, accepted for publication in Astrophysical Journa

Determination of antibodies to Flavivirus in wild birds from Buenos Aires City

Saint Louis encephalitis virus (SLEV), West Nile virus (WNV) and Ilheus virus (ILHV) are flaviviruses maintained in nature by enzootic transmission networks between mosquitoes and birds. These viruses have been detected in South America and identified as a cause of neurological diseases in humans. There is a record of activity in monkeys and birds for ILHV in Chaco, Corrientes and Misiones; of WNV in horses and birds in the center and east of the country, while SLEV is widely distributed in Argentina and there is virological and serological evidence in birds, humans and various mammals. In order to detect viral circulation in free-living birds in the Autonomous City of Buenos Aires (CABA), seasonal sampling was carried out in the Costanera Sur Ecological Reserve.Para acceder a la videoconferencia completa, hacer clic en "Enlace externo".Sociedad Latinoamericana de Ecología de Vectore

Hydrodynamical interaction of stellar and planetary winds: Effects of charge exchange and radiation pressure on the observed Ly α absorption

Lyman α observations of the transiting exoplanet HD 209458b enable the study of exoplanet exospheres exposed to stellar extreme ultraviolet (EUV) fluxes, as well as the interacting stellar wind properties. In this study we present 3D hydrodynamical models for the stellar-planetary 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 to 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 yields a small velocity shift of the absorbing material. While neither charge exchange nor radiation pressure provides 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.Fil: Esquivel, A.. Universidad Nacional Autónoma de México. Instituto de Ciencias Nucleares; México. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Astronomía Teórica y Experimental. Universidad Nacional de Córdoba. Observatorio Astronómico de Córdoba. Instituto de Astronomía Teórica y Experimental; ArgentinaFil: Schneiter, Ernesto Matías. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Astronomía Teórica y Experimental. Universidad Nacional de Córdoba. Observatorio Astronómico de Córdoba. Instituto de Astronomía Teórica y Experimental; Argentina. Stockholms Universitet; SueciaFil: Villarreal D'angelo, Carolina Susana. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Astronomía Teórica y Experimental. Universidad Nacional de Córdoba. Observatorio Astronómico de Córdoba. Instituto de Astronomía Teórica y Experimental; ArgentinaFil: Sgró, Mario Agustín. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Astronomía Teórica y Experimental. Universidad Nacional de Córdoba. Observatorio Astronómico de Córdoba. Instituto de Astronomía Teórica y Experimental; ArgentinaFil: Krapp, Leonardo Javier. Stockholms Universitet; Suecia. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Astronomía Teórica y Experimental. Universidad Nacional de Córdoba. Observatorio Astronómico de Córdoba. Instituto de Astronomía Teórica y Experimental; Argentin