Dust plays a key role during star, disk and planet formation. Yet, its dynamics during the protostellar collapse remains a poorly investigated field. Recent studies seem to indicate that dust may decouple efficiently from the gas during these early stages. We aim to understand how much and in which regions dust grains concentrate during the early phases of the protostellar collapse, and see how it depends on the properties of the initial cloud and of the solid particles. We use the multiple species dust dynamics solver of the grid-based code RAMSES to perform various simulations of dusty collapses. We perform hydrodynamical and MHD simulations where we vary the maximum grain size, the thermal-to-gravitational energy ratio and the magnetic properties of the cloud. We simulate the simultaneous evolution of ten neutral dust grains species with grain sizes varying from a few nm to a few hundredth of microns. We obtain a significant decoupling between the gas and the dust for grains of typical sizes a few 10 microns. This decoupling strongly depends on the thermal-to-gravitational energy ratio, the grain sizes or the inclusion of a magnetic field. With a semi-analytic model calibrated on our results, we show that the dust ratio mostly varies exponentially with the initial Stokes number at a rate that depends on the local cloud properties. We find that larger grains tend to settle and drift efficiently in the first-core and in the newly formed disk. This can produce dust-to-gas ratios of several times the initial value. Dust concentrates in high density regions and is depleted in low density regions. The size at which grains decouple from the gas depends on the initial properties of the clouds. Since dust can not necessarily be used as a proxy for gas during the collapse, we emphasize on the necessity of including the treatment of its dynamics in collapse simulations.

Dusty debris disks around pre- and main-sequence stars are potential signposts for the existence of planetesimals and exoplanets. Giant planet formation is therefore expected to play a key role in the evolution of the disk. This is indirectly confirmed by extant sub-millimeter near-infrared images of young protoplanetary and cool dusty debris disks around main sequence stars usually showing substantial spatial structures. A majority of recent discoveries of imaged giant planets have been obtained around young, early-type stars hosting a circumstellar disk. In this context, we have carried out a direct imaging program designed to maximize our chances of giant planet discovery and targeting twenty-two young, early-type stars. About half of them show indication of multi-belt architectures. Using the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE to acquire high-constrast coronagraphic differential near-infrared images, we have conducted a systematic search in the close environment of these young, dusty and early-type stars. We confirmed that companions detected around HIP 34276, HIP 101800 and HIP 117452 are stationary background sources and binary companions. The companion candidates around HIP 8832, HIP 16095 and HIP 95619 are determined as background contamination. For stars for which we infer the presence of debris belts, a theoretical minimum mass for planets required to clear the debris gaps can be calculated . The dynamical mass limit is at least $0.1 M_J$ and can exceed $1 M_J$. Direct imaging data is typically sensitive to planets down to $\sim 3.6 M_J$ at 1 $''$, and $1.7 M_J$ in the best case. These two limits tightly constrain the possible planetary systems present around each target. These systems will be probably detectable with the next generation of planet imagers.

We present a fix to the overdamping problem found by Laibe & Price (2012) when simulating strongly coupled dust-gas mixtures using two different sets of particles using smoothed particle hydrodynamics. Our solution is to compute the drag at the barycentre between gas and dust particle pairs when computing the drag force by reconstructing the velocity field, similar to the procedure in Godunov-type solvers. This fixes the overdamping problem at negligible computational cost, but with additional memory required to store velocity derivatives. We employ slope limiters to avoid spurious oscillations at shocks, finding the van Leer Monotonized Central limiter most effective.

Instruments achieve sharper and finer observations of micron-in-size dust grains in the top layers of young stellar discs. To provide accurate models, we revisit the theory of dust settling for small grains, when gas stratification, dust inertia and finite correlation times for the turbulence should be handled simultaneously. We start from a balance of forces and derive distributions at steady-state. Asymptotic expansions require caution since limits do not commute. In particular, non-physical bumpy distributions appear when turbulence is purely diffusive. This excludes very short correlation times for real discs, as predicted by numerical simulations.

SPHERE (Beuzit et al,. 2019) has now been in operation at the VLT for more than 5 years, demonstrating a high level of performance. SPHERE has produced outstanding results using a variety of operating modes, primarily in the field of direct imaging of exoplanetary systems, focusing on exoplanets as point sources and circumstellar disks as extended objects. The achievements obtained thus far with SPHERE (~200 refereed publications) in different areas (exoplanets, disks, solar system, stellar physics...) have motivated a large consortium to propose an even more ambitious set of science cases, and its corresponding technical implementation in the form of an upgrade. The SPHERE+ project capitalizes on the expertise and lessons learned from SPHERE to push high contrast imaging performance to its limits on the VLT 8m-telescope. The scientific program of SPHERE+ described in this document will open a new and compelling scientific window for the upcoming decade in strong synergy with ground-based facilities (VLT/I, ELT, ALMA, and SKA) and space missions (Gaia, JWST, PLATO and WFIRST). While SPHERE has sampled the outer parts of planetary systems beyond a few tens of AU, SPHERE+ will dig into the inner regions around stars to reveal and characterize by mean of spectroscopy the giant planet population down to the snow line. Building on SPHERE's scientific heritage and resounding success, SPHERE+ will be a dedicated survey instrument which will strengthen the leadership of ESO and the European community in the very competitive field of direct imaging of exoplanetary systems. With enhanced capabilities, it will enable an even broader diversity of science cases including the study of the solar system, the birth and death of stars and the exploration of the inner regions of active galactic nuclei.

Streaming instability is a privileged channel to bridge the gap between collisional growth of dust grains and planetesimal formation triggered by gravity. This instability is thought to develop through its secular mode, which is long-time growing and may not develop easily in real discs. We address this point by revisiting its perturbation analysis. A third-order expansion with respect to the Stokes number reveals important features over-looked so far. The secular mode can be stable. Epicycles can be unstable, more resistant to viscosity and are identified by Green's function analysis as promising channels for planetesimals formation.

Context. Atmospheric superrotating flows at the equator are an almost ubiquitous result of simulations of hot Jupiters, and a theory explaining how this zonally coherent flow reaches an equilibrium has been developed in the literature. However, this understanding relies on the existence of either an initial superrotating or a sheared flow, coupled with a slow evolution such that a linear steady state can be reached. Aims. A consistent physical understanding of superrotation is needed for arbitrary drag and radiative timescales, and the relevance of considering linear steady states needs to be assessed. Methods. We obtain an analytical expression for the structure, frequency and decay rate of propagating waves in hot Jupiter atmospheres around a state at rest in the 2D shallow water beta plane limit. We solve this expression numerically and confirm the robustness of our results with a 3D linear wave algorithm. We then compare with 3D simulations of hot Jupiter atmospheres and study the non linear momentum fluxes. Results. We show that under strong day night heating the dynamics does not transit through a linear steady state when starting from an initial atmosphere in solid body rotation. We further show that non linear effects favour the initial spin up of superrotation and that the acceleration due to the vertical component of the eddy momentum flux is critical to the initial development of superrotation. Conclusions. Overall, we describe the initial phases of the acceleration of superrotation, including consideration of differing radiative and drag timescales, and conclude that eddy-momentum driven superrotating equatorial jets are robust, physical phenomena in simulations of hot Jupiter atmospheres.

Tidal encounters in star clusters perturb discs around young protostars. In Cuello et al. (2019a, Paper I) we detailed the dynamical signatures of a stellar flyby in both gas and dust. Flybys produce warped discs, spirals with evolving pitch angles, increasing accretion rates, and disc truncation. Here we present the corresponding observational signatures of these features in optical/near-infrared scattered light and (sub-) millimeter continuum and CO line emission. Using representative prograde and retrograde encounters for direct comparison, we post-process hydrodynamical simulations with radiative transfer methods to generate a catalogue of multi-wavelength observations. This provides a reference to identify flybys in recent near-infrared and sub-millimetre observations (e.g., RW Aur, AS 205, HV Tau & DO Tau, FU Ori, V2775 Ori, and Z CMa).

Small dust grains are essential ingredients of star, disk and planet formation. We present an Eulerian numerical approach to study small dust grains dynamics in the context of star and protoplanetary disk formation. It is designed for finite volume codes. We use it to investigate dust dynamics during the protostellar collapse. We present a method to solve the monofluid equations of gas and dust mixtures with several dust species in the diffusion approximation implemented in the adaptive-mesh-refinement code RAMSES. It uses a finite volume second-order Godunov method with a predictor-corrector MUSCL scheme to estimate the fluxes between the grid cells. We benchmark our method against six distinct tests, dustyadvect, dustydiffuse, dustyshock, dustywave, settling and dustycollapse. We show that the scheme is second-order accurate in space on uniform grids and intermediate between second and first-order on non-uniform grids. We apply our method on various dustycollapse simulations of 1 solar mass cores composed of gas and dust. We developed an efficient approach to treat gas and dust dynamics in the diffusion regime on grid based codes. The canonical tests have been successfully passed. In the context of protostellar collapse, we show that dust is less coupled to the gas in the outer regions of the collapse where grains larger than approximately 100 micrometers fall significantly faster than the gas.

We present 3D smoothed particle hydrodynamics simulations of protoplanetary discs undergoing a flyby by a stellar perturber on a parabolic orbit lying in a plane inclined relative to the disc mid-plane. We model the disc as a mixture of gas and dust, with grains ranging from 1 {\mu}m to 10 cm in size. Exploring different orbital inclinations, periastron distances and mass ratios, we investigate the disc dynamical response during and after the flyby. We find that flybys induce evolving spiral structure in both gas and dust which can persist for thousands of years after periastron. Gas and dust structures induced by the flyby differ because of drag-induced effects on the dust grains. Variations in the accretion rate by up to an order of magnitude occur over a time-scale of order 10 years or less, inducing FU Orionis-like outbursts. The remnant discs are truncated and warped. The dust disc is left more compact than the gas disc, both because of disc truncation and accelerated radial drift of grains induced by the flyby.

We study the dynamics of a viscous protoplanetary disc hosting a population of dust grains with a range of sizes. We compute steady-state solutions, and show that the radial motion of both the gas and the dust can deviate substantially from those for a single-size dust population. Although the aerodynamic drag from the dust on the gas is weaker than in the case where all grains are optimally coupled to the gas, the cumulative "back-reaction" of the dust particles can still alter the gas dynamics significantly. In typical protoplanetary discs, the net effect of the dust back-reaction decreases the gas accretion flow compared to the dust-free (viscous) case, even for dust-to-gas ratios of order $1\%$. In the outer disc, where dust grains are typically less strongly coupled to the gas and settle towards the midplane, the dust back-reaction can even drive outward gas flow. Moreover, the radial inward drift of large grains is reduced below the gas motion in the inner disc regions, while small dust grains follow the gas dynamics over all the disc extent. The resulting dust and gas dynamics can give rise to observable structures, such as gas and dust cavities. Our results show that the dust back-reaction can play a major role in both the dynamics and observational appearance of protoplanetary discs, and cannot be ignored in models of protoplanetary disc evolution.

We present a methodology to interpret observations of protoplanetary discs where a flyby, also called a tidal encounter, is suspected. In case of a flyby, protoplanetary discs can be significantly disturbed. The resulting dynamical and kinematical signatures can last for several thousands of years after the flyby and hence deeply affect the evolution of the disc. These effects are stronger for closer encounters and more massive perturbers. For the very same flyby parameters, varying the inclination of the perturber's orbit produces a broad range of disc structures: spirals, bridges, warps and cavities. We study this kind of features both in the gas and in the dust for grains ranging from 1 {\mu}m to 10 cm in size. Interestingly, the dust exhibits a different dynamical behaviour compared to the gas because of gas-drag effects. Finally, flybys can also trigger high accretion events in the disc-hosting star, readily similar to FU Orionis-type outbursts. All this information can be used to infer the flyby parameters from an incomplete set of observations at different wavelengths. Therefore, the main scope of our flyby scene investigation (FSI) methodology is to help to interpret recent "puzzling" disc observations.

We present Phantom, a fast, parallel, modular and low-memory smoothed particle hydrodynamics and magnetohydrodynamics code developed over the last decade for astrophysical applications in three dimensions. The code has been developed with a focus on stellar, galactic, planetary and high energy astrophysics and has already been used widely for studies of accretion discs and turbulence, from the birth of planets to how black holes accrete. Here we describe and test the core algorithms as well as modules for magnetohydrodynamics, self-gravity, sink particles, H_2 chemistry, dust-gas mixtures, physical viscosity, external forces including numerous galactic potentials as well as implementations of Lense-Thirring precession, Poynting-Robertson drag and stochastic turbulent driving. Phantom is hereby made publicly available.

We present MULTIGRAIN, an algorithm for simulating multiple phases of small dust grains embedded in a gas, building on our earlier work in simulating two-phase mixtures of gas and dust in SPH (Laibe & Price 2012a,b; Price & Laibe 2015). The MULTIGRAIN method (Hutchison, Price & Laibe 2018) is more accurate than single-phase simulations because the gas experiences a backreaction from each dust phase and communicates this change to the other phases, thereby indirectly coupling the dust phases together. The MULTIGRAIN method is fast, explicit and low storage, requiring only an array of dust fractions and their derivatives defined for each resolution element. We demonstrate the MULTIGRAIN algorithm on test problems related to the settling of dust in the discs of gas around young stars, where solar systems are born. Finally I will discuss possible extensions of the method to incorporate both large and small grains, together with recent improvements in our numerical techniques for gas and dust mixtures. In particular, I will show how the 'overdamping' problem identified by Laibe & Price (2012a) can be solved.

We describe a new implementation of the one-fluid method in the SPH code Phantom to simulate the dynamics of dust grains in gas protoplanetary discs. We revise and extend previously developed algorithms by computing the evolution of a new fluid quantity that produces a more accurate and numerically controlled evolution of the dust dynamics. Moreover, by limiting the stopping time of uncoupled grains that violate the assumptions of the terminal velocity approximation, we avoid fatal numerical errors in mass conservation. We test and validate our new algorithm by running 3D SPH simulations of a large range of disc models with tightly- and marginally-coupled grains.

We present 3D hydrodynamical models of the HD142527 protoplanetary disc, a bright and well studied disc that shows spirals and shadows in scattered light around a 100 au gas cavity, a large horseshoe dust structure in mm continuum emission, together with mysterious fast radial flows and streamers seen in gas kinematics. By considering several possible orbits consistent with the observed arc, we show that all of the main observational features can be explained by one mechanism - the interaction between the disc and the observed binary companion. We find that the spirals, shadows and horseshoe are only produced in the correct position angles by a companion on an inclined and eccentric orbit approaching periastron - the 'red' family from Lacour et al. (2016). Dust-gas simulations show radial and azimuthal concentration of dust around the cavity, consistent with the observed horseshoe. The success of this model in the HD142527 disc suggests other mm-bright transition discs showing cavities, spirals and dust asymmetries may also be explained by the interaction with central companions.

We present a new algorithm, MULTIGRAIN, for modelling the dynamics of an entire population of small dust grains immersed in gas, typical of conditions that are found in molecular clouds and protoplanetary discs. The MULTIGRAIN method is more accurate than single-phase simulations because the gas experiences a backreaction from each dust phase and communicates this change to the other phases, thereby indirectly coupling the dust phases together. The MULTIGRAIN method is fast, explicit and low storage, requiring only an array of dust fractions and their derivatives defined for each resolution element.

We present Atacama Large Millimeter/sub-millimeter Array (ALMA) Cycle 2 observations of the 1.3 mm dust continuum emission of the protoplanetary disc surrounding the T Tauri star Elias 24 with an angular resolution of $\sim 0.2"$ ($\sim 28$ au). The dust continuum emission map reveals a dark ring at a radial distance of $0.47"$ ($\sim 65$ au) from the central star, surrounded by a bright ring at $0.58"$ ($\sim 81$ au). In the outer disc, the radial intensity profile shows two inflection points at $0.71"$ and $0.87"$ ($\sim 99$ and $121$ au respectively). We perform global three-dimensional smoothed particle hydrodynamic gas/dust simulations of discs hosting a migrating and accreting planet. Combining the dust density maps of small and large grains with three dimensional radiative transfer calculations, we produce synthetic ALMA observations of a variety of disc models in order to reproduce the gap- and ring-like features observed in Elias 24. We find that the dust emission across the disc is consistent with the presence of an embedded planet with a mass of $\sim 0.7\, \mathrm{M_{\mathrm{J}}}$ at an orbital radius of $\sim$ 60 au. Our model suggests that the two inflection points in the radial intensity profile are due to the inward radial motion of large dust grains from the outer disc. The surface brightness map of our disc model provides a reasonable match to the gap- and ring-like structures observed in Elias 24, with an average discrepancy of $\sim$ 5% of the observed fluxes around the gap region.

Voronoi grids have been successfully used to represent density structures of gas in astronomical hydrodynamics simulations. While some codes are explicitly built around using a Voronoi grid, others, such as Smoothed Particle Hydrodynamics (SPH), use particle-based representations and can benefit from constructing a Voronoi grid for post-processing their output. So far, calculating the density of each Voronoi cell from SPH data has been done numerically, which is both slow and potentially inaccurate. This paper proposes an alternative analytic method, which is fast and accurate. We derive an expression for the integral of a cubic spline kernel over the volume of a Voronoi cell and link it to the density of the cell. Mass conservation is ensured rigorously by the procedure. The method can be applied more broadly to integrate a spherically symmetric polynomial function over the volume of a random polyhedron.

Streaming instability is a powerful mechanism which concentrates dust grains in pro- toplanetary discs, eventually up to the stage where they collapse gravitationally and form planetesimals. Previous studies inferred that it should be ineffective in viscous discs, too efficient in inviscid discs, and may not operate in local pressure maxima where solids accumulate. From a linear analysis of stability, we show that streaming instability behaves differently inside local pressure maxima. Under the action of the strong differential advection imposed by the bump, a novel unstable mode develops and grows even when gas viscosity is large. Hence, pressure bumps are found to be the only places where streaming instability occurs in viscous discs. This offers a promising way to conciliate models of planet formation with recent observations of young discs.

Planet formation is thought to occur in discs around young stars by the aggregation of small dust grains into much larger objects. The growth from grains to pebbles and from planetesimals to planets is now fairly well understood. The intermediate stage has however been found to be hindered by the radial-drift and fragmentation barriers. We identify a powerful mechanism in which dust overcomes both barriers. Its key ingredients are i) backreaction from the dust onto the gas, ii) grain growth and fragmentation, and iii) large-scale gradients. The pile-up of growing and fragmenting grains modifies the gas structure on large scales and triggers the formation of pressure maxima, in which particles are trapped. We show that these self-induced dust traps are robust: they develop for a wide range of disc structures, fragmentation thresholds and initial dust-to-gas ratios. They are favored locations for pebbles to grow into planetesimals, thus opening new paths towards the formation of planets.

We perform numerical simulations of dusty, supersonic turbulence in molecular clouds. We model 0.1, 1 and 10 {\mu}m sized dust grains at an initial dust-to-gas mass ratio of 1:100, solving the equations of combined gas and dust dynamics where the dust is coupled to the gas through a drag term. We show that, for 0.1 and 1 {\mu}m grains, the dust-to-gas ratio deviates by typically 10-20% from the mean, since the stopping time of the dust due to gas drag is short compared to the dynamical time. Contrary to previous findings, we find no evidence for orders of magnitude fluctuation in the dust-to-gas ratio for 0.1 {\mu}m grains. Larger, 10 {\mu}m dust grains may have dust-to-gas ratios increased by up to an order of magnitude locally. Both small (0.1 {\mu}m) and large ($\gtrsim$ 1 {\mu}m) grains trace the large-scale morphology of the gas, however we find evidence for 'size-sorting' of grains, where turbulence preferentially concentrates larger grains into dense regions. Size-sorting may help to explain observations of 'coreshine' from dark clouds, and why extinction laws differ along lines of sight through molecular clouds in the Milky Way compared to the diffuse interstellar medium.

We aim to understand under which conditions a low mass planet can open a gap in viscous dusty protoplanetary discs. For this purpose, we extend the theory of dust radial drift to include the contribution from the tides of an embedded planet and from the gas viscous forces. From this formalism, we derive i) a grain size-dependent criterion for dust gap opening in discs, ii) an estimate of the location of the outer edge of the dust gap and iii) an estimate of the minimum Stokes number above which low-mass planets are able to carve gaps which appear only in the dust disc. These analytical estimates are particularly helpful to appraise the minimum mass of an hypothetical planet carving gaps in discs observed at long wavelengths and high resolution. We validate the theory against 3D SPH simulations of planet-disc interaction in a broad range of dusty protoplanetary discs. We find a remarkable agreement between the theoretical model and the numerical experiments.

We investigate whether the rings, lopsided features and horseshoes observed at millimetre wavelengths in transitional discs can be explained by the dynamics of gas and dust at the edge of the cavity in circumbinary discs. We use 3D dusty smoothed particle hydrodynamics calculations to show that binaries with mass ratio $q \gtrsim 0.04$ drive eccentricity in the central cavity, naturally leading to a crescent-like feature in the gas density, which is accentuated in the mm dust grain population with intensity contrasts in mm-continuum emission of 10 or higher. We perform mock observations to demonstrate that these features closely match those observed by ALMA, suggesting that the origin of rings, dust horseshoes and other non-axisymmetric structures in transition discs can be explained by the presence of massive companions.

The Protoplanetary Discussions conference --- held in Edinburgh, UK, from 7th --11th March 2016 --- included several open sessions led by participants. This paper reports on the discussions collectively concerned with the multiphysics modelling of protoplanetary discs, including the self-consistent calculation of gas and dust dynamics, radiative transfer and chemistry. After a short introduction to each of these disciplines in isolation, we identify a series of burning questions and grand challenges associated with their continuing development and integration. We then discuss potential pathways towards solving these challenges, grouped by strategical, technical and collaborative developments. This paper is not intended to be a review, but rather to motivate and direct future research and collaboration across typically distinct fields based on \textit{community driven input}, to encourage further progress in our understanding of circumstellar and protoplanetary discs.

We model the behaviour of dust grains entrained by photoevaporation-driven winds from protoplanetary discs assuming a non-rotating, plane-parallel disc. We obtain an analytic expression for the maximum entrainable grain size in extreme-UV radiation-driven winds, which we demonstrate to be proportional to the mass loss rate of the disc. When compared with our hydrodynamic simulations, the model reproduces almost all of the wind properties for the gas and dust. In typical turbulent discs, the entrained grain sizes in the wind are smaller than the theoretical maximum everywhere but the inner disc due to dust settling.

We investigate dust entrainment by photoevaporative winds in protoplanetary discs using dusty smoothed particle hydrodrodynamics (SPH). We use unequal-mass particles to resolve more than five orders of magnitude in disc/outflow density and a one-fluid formulation to efficiently simulate an equivalent magnitude range in drag stopping time. We find that only micron sized dust grains and smaller can be entrained in EUV driven winds. The maximum grain size is set by dust settling in the disc rather than aerodynamic drag in the wind. More generally, there is a linear relationship between the base flow density and the maximum entrainable grain size in the wind. A pileup of micron sized dust grains can occur in the upper atmosphere at critical radii in the disc as grains decouple from the low-density wind. Entrainment is a strong function of location in the disc, resulting in a size sorting of grains in the outflow---the largest grain being carried out between $10$--$20\,$AU. The peak dust density for each grain size occurs at the inner edge of its own entrainment region.

Here we derive a Parker-wind like solution for a stratified, plane-parallel atmosphere undergoing photoionisation. The difference compared to the standard Parker solar wind is that the sonic point is crossed only at infinity. The simplicity of the analytic solution makes it a convenient test problem for numerical simulations of photoevaporation in protoplanetary discs.

We identify two distinct physical mechanisms for dust gap opening by embedded planets in protoplanetary discs based on the symmetry of the drag-induced motion around the planet: I) A mechanism where low mass planets, that do not disturb the gas, open gaps in dust by tidal torques assisted by drag in the inner disc, but resisted by drag in the outer disc; and II) The usual, drag assisted, mechanism where higher mass planets create pressure maxima in the gas disc which the drag torque then acts to evacuate further in the dust. The first mechanism produces gaps in dust but not gas, while the second produces partial or total gas gaps which are deeper in the dust phase. Dust gaps do not necessarily indicate gas gaps.

We describe a simple method for simulating the dynamics of small grains in a dusty gas, relevant to micron-sized grains in the interstellar medium and grains of centimetre size and smaller in protoplanetary discs. The method involves solving one extra diffusion equation for the dust fraction in addition to the usual equations of hydrodynamics. This "diffusion approximation for dust" is valid when the dust stopping time is smaller than the computational timestep. We present a numerical implementation using Smoothed Particle Hydrodynamics (SPH) that is conservative, accurate and fast. It does not require any implicit timestepping and can be straightforwardly ported into existing 3D codes.

Protoplanetary discs are now routinely observed and exoplanets, after the numerous indirect discoveries, are starting to be directly imaged. To better understand the planet formation process, the next step is the detection of forming planets or of signposts of young planets still in their disc, such as gaps. A spectacular example is the ALMA science verification image of HL Tau showing numerous gaps and rings in its disc. To study the observability of planet gaps, we ran 3D hydrodynamical simulations of a gas and dust disc containing a 5 M J gap-opening planet and characterised the spatial distribution of migrating, growing and fragmenting dust grains. We then computed the corresponding synthetic images for ALMA. For a value of the dust fragmentation threshold of 15 m s --1 for the collisional velocity, we identify for the first time a self-induced dust pile up in simulations taking fragmentation into account. This feature, in addition to the easily detected planet gap, causes a second apparent gap that could be mistaken for the signature of a second planet. It is therefore essential to be cautious in the interpretation of gap detections.

We explain the axisymmetric gaps seen in recent long-baseline observations of the HL Tau protoplanetary disc with the Atacama Large Millimetre/Submillimetre Array (ALMA) as being due to the different response of gas and dust to embedded planets in protoplanetary discs. We perform global, three dimensional dusty smoothed particle hydrodynamics calculations of multiple planets embedded in dust/gas discs which successfully reproduce most of the structures seen in the ALMA image. We find a best match to the observations using three embedded planets with masses of 0.2, 0.27 and 0.55 $M_{\rm J}$ in the three main gaps observed by ALMA, though there remain uncertainties in the exact planet masses from the disc model.

We model the dust evolution in protoplanetary disks with full 3D, Smoothed Particle Hydrodynamics (SPH), two-phase (gas+dust) hydrodynamical simulations. The gas+dust dynamics, where aerodynamic drag leads to the vertical settling and radial migration of grains, is consistently treated. In a previous work, we characterized the spatial distribution of non-growing dust grains of different sizes in a disk containing a gap-opening planet and investigated the gap's detectability with the Atacama Large Millimeter/submillimeter Array (ALMA). Here we take into account the effects of grain growth and fragmentation and study their impact on the distribution of solids in the disk. We show that rapid grain growth in the two accumulation zones around planet gaps is strongly affected by fragmentation. We discuss the consequences for ALMA observations.

We present a new approach to simulating mixtures of gas and dust in smoothed particle hydrodynamics (SPH). We show how the two-fluid equations can be rewritten to describe a single-fluid 'mixture' moving with the barycentric velocity, with each particle carrying a dust fraction. We show how this formulation can be implemented in SPH while preserving the conservation properties (i.e. conservation of mass of each phase, momentum and energy). We also show that the method solves two key issues with the two fluid approach: it avoids over-damping of the mixture when the drag is strong and prevents a problem with dust particles becoming trapped below the resolution of the gas. We also show how the general one-fluid formulation can be simplified in the limit of strong drag (i.e. small grains) to the usual SPH equations plus a diffusion equation for the evolution of the dust fraction that can be evolved explicitly and does not require any implicit timestepping. We present tests of the simplified formulation showing that it is accurate in the small grain/strong drag limit. We discuss some of the issues we have had to solve while developing this method and finally present a preliminary application to dust settling in protoplanetary discs.

Dust-gas mixtures are the simplest example of a two fluid mixture. We show that when simulating such mixtures with particles or with particles coupled to grids a problem arises due to the need to resolve a very small length scale when the coupling is strong. Since this is occurs in the limit when the fluids are well coupled, we show how the dust-gas equations can be reformulated to describe a single fluid mixture. The equations are similar to the usual fluid equations supplemented by a diffusion equation for the dust-to-gas ratio or alternatively the dust fraction. This solves a number of numerical problems as well as making the physics clear.

We derive the single-fluid evolution equations describing a mixture made of a gas phase and an arbitrary number of dust phases, generalising the approach developed in Laibe & Price (2014a). A generalisation for continuous dust distributions as well as analytic approximations for strong drag regimes are also provided. This formalism lays the foundation for numeri- cal simulations of dust populations in a wide range of astrophysical systems while avoiding limitations associated with a multiple-fluid treatment. The usefulness of the formalism is illustrated on a series of analytical problems, namely the dustybox, dustyshock and dustywave problems as well as the radial drift of grains and the streaming instability in protoplanetary discs. We find physical effects specific to the presence of several dust phases and multiple drag timescales, including non-monotonic evolution of the differential velocity between phases and increased efficiency of the linear growth of the streaming instability. Interestingly, it is found that under certain conditions, large grains can migrate outwards in protoplanetary discs. This may explain the presence of small pebbles at several hundreds of astronomical units from their central star.

The growth of dust particles into planet embryos needs to circumvent the radial-drift barrier, i.e. the accretion of dust particles onto the central star by radial migration. The outcome of the dust radial migration is governed by simple criteria between the dust-to-gas ratio and the exponents p and q of the surface density and temperature power laws. The transfer of radiation provides an additional constraint between these quantities because the disc thermal structure is fixed by the dust spatial distribution. To assess which discs are primarily affected by the radial-drift barrier, we used the radiative transfer code MCFOST to compute the temperature structure of a wide range of disc models, stressing the particular effects of grain size distributions and vertical settling. We find that the outcome of the dust migration process is very sensitive to the physical conditions within the disc. For high dust-to-gas ratios (> 0.01) or flattened disc structures (H/R < 0.05), growing dust grains can efficiently decouple from the gas, leading to a high concentration of grains at a critical radius of a few AU. Decoupling of grains can occur at a large fraction (> 0.1) of the initial radius, for a dust-to-gas ratio greater than ~ 0.05. The exact value of the required dust-to-gas ratio for dust to stop its migration is strongly dependent on the disc temperature structure. Non growing dust grains are accreted for discs with flat surface density profiles (p<0.7) while they always remain in the disc if the surface density is steep enough (p>1.2). Both the presence of large grains and vertical settling tend to favour the accretion of non growing dust grains onto the central object, but it slows down the migration of growing dust grains. All the disc configurations are found to have favourable temperature profiles over most of the disc to retain their planetesimals.

In this paper, we show how the two-fluid equations describing the evolution of a dust and gas mixture can be reformulated to describe a single fluid moving with the barycentric velocity of the mixture. This leads to evolution equations for the total density, momentum, the differential velocity between the dust and the gas phases and either the dust-to-gas ratio or the dust fraction. The equations are similar to the usual equations of gas dynamics, providing a convenient way to extend existing codes to simulate two-fluid mixtures without modifying the code architecture. Our approach avoids the inherent difficulties related to the standard approach where the two phases are separate and coupled via a drag term. In particular, the requirements of infinite spatial and temporal resolution as the stopping time tends to zero are no longer necessary. This means that both small and large grains can be straightforwardly treated with the same method, with no need for complicated implicit schemes. Since there is only one resolution scale the method also avoids the problem of unphysical trapping of one fluid (e.g. dust) below the resolution of the other. We also derive a simplified set of equations applicable to the case of strong drag/small grains, consisting of the standard fluid equations with a modified sound speed, plus an advection-diffusion equation for the dust-to-gas ratio. This provides a simple and fast way to evolve the mixture when the stopping time is smaller than the Courant timestep. We present a Smoothed Particle Hydrodynamics implementation in a companion paper.

In a companion paper we have shown how the equations describing gas and dust as two fluids coupled by a drag term can be reformulated to describe the system as a single fluid mixture. Here we present a numerical implementation of the one-fluid dusty gas algorithm using Smoothed Particle Hydrodynamics (SPH). The algorithm preserves the conservation properties of the SPH formalism. In particular, the total gas and dust mass, momentum, angular momentum and energy are all exactly conserved. Shock viscosity and conductivity terms are generalised to handle the two-phase mixture accordingly. The algorithm is benchmarked against a comprehensive suit of problems: dustybox, dustywave, dustyshock and dustyoscill, each of them addressing different properties of the method. We compare the performance of the one-fluid algorithm to the standard two-fluid approach. The one-fluid algorithm is found to solve both of the fundamental limitations of the two- fluid algorithm: it is no longer possible to concentrate dust below the resolution of the gas (they have the same resolution by definition), and the spatial resolution criterion h < csts, required in two-fluid codes to avoid over-damping of kinetic energy, is unnecessary. Implicit time stepping is straightforward. As a result, the algorithm is up to ten billion times more efficient for 3D simulations of small grains. Additional benefits include the use of half as many particles, a single kernel and fewer SPH interpolations. The only limitation is that it does not capture multi-streaming of dust in the limit of zero coupling, suggesting that in this case a hybrid approach may be required.

We aim to derive a simple analytic model to understand the essential properties of vertically settling growing dust grains in laminar protoplanetary discs. Separating the vertical dynamics from the motion in the disc midplane, we integrate the equations of motion for both a linear and an exponential grain growth rate. Numerical integrations are performed for more complex growth models. We find that the settling efficiency depends on the value of the dimensionless parameter gamma, which characterises the relative efficiency of grain growth with respect to the gas drag. Since gamma is expected to be of order as the initial dust-to-gas ratio in the disc (of order 10^-2), grain growth enhances the energy dissipation of the dust particles and improve the settling efficiency in protoplanetary discs. This behaviour is mostly independent of the growth model considered as well as of the radial drift of the particles.

We aim to study the migration of growing dust grains in protoplanetary discs, where growth and migration are tightly coupled. This includes the crucial issue of the radial-drift barrier for growing dust grains. We therefore extend the study performed in Paper I, considering models for grain growth and grain dynamics where both the migration and growth rate depend on the grain size and the location in the disc. The parameter space of disc profiles and growth models is exhaustively explored. In doing so, interpretations for the grain motion found in numerical simulations are also provided. We find that a large number of cases is required to characterise entirely the grains radial motion, providing a large number of possible outcomes. Some of them lead dust particles to be accreted onto the central star and some of them don't. We find then that q<1 is required for discs to retain their growing particles, where q is the exponent of the radial temperature profile T(R) proportional to R^-q. Additionally, the initial dust-to gas ratio has to exceed a critical value for grains to pile up efficiently, thus avoiding being accreted onto the central star. Discs are also found to retain efficiently small dust grains regenerated by fragmentation. We show how those results are sensitive to the turbulent model considered. Even though some physical processes have been neglected, this study allows to sketch a scenario in which grains can survive the radial-drift barrier in protoplanetary discs as they grow.

In a series of papers, we present a comprehensive analytic study of the global motion of growing dust grains in protoplanetary discs, addressing both the radial drift and the vertical settling of the particles. Here we study how the radial drift of dust particles is affected by grain growth. In a first step, toy models in which grain growth can either be constant, accelerate or decelerate are introduced. The equations of motion are analytically integrable and therefore the grains dynamics is easy to understand. The radial motion of growing grains is governed by the relative efficiency of the growth and migration processes which is expressed by the dimensionless parameter Lambda, as well as the exponents for the gas surface density and temperature profiles, denoted p and q respectively. When Lambda is of order unity, growth and migration are strongly coupled, providing the most efficient radial drift. For the toy models considered, grains pile up when -p+q+1/2<0. Importantly, we show the existence of a second process which can help discs to retain their solid materials. For accelerating growth, grains end up their migration at a finite radius, thus avoiding being accreted onto the central star.

We have performed three-dimensional two-fluid (gas-dust) hydrodynamical models of circumstellar discs with embedded protoplanets (3 - 333 M\oplu) and small solid bodies (radii 10cm to 10m). We find that high mass planets (\gtrsim Saturn mass) open sufficiently deep gaps in the gas disc such that the density maximum at the outer edge of the gap can very efficiently trap metre-sized solid bodies. This allows the accumulation of solids at the outer edge of the gap as solids from large radii spiral inwards to the trapping region. This process of accumulation occurs fastest for those bodies that spiral inwards most rapidly, typically metre-sized boulders, whilst smaller and larger objects will not migrate sufficiently rapidly in the discs lifetime to benefit from the process. Around a Jupiter mass planet we find that bound clumps of solid material, as large as several Earth masses, may form, potentially collapsing under self-gravity to form planets or planetesimals. These results are in agreement with Lyra et al. (2009), supporting their finding that the formation of a second generation of planetesimals or of terrestrial mass planets may be triggered by the presence of a high mass planet.

We present a new algorithm for simulating two-fluid gas and dust mixtures in Smoothed Particle Hydrodynamics (SPH), systematically addressing a number of key issues including the generalised SPH density estimate in multi-fluid systems, the consistent treatment of variable smoothing length terms, finite particle size, time step stability, thermal coupling terms and the choice of kernel and smoothing length used in the drag operator. We find that using double-hump shaped kernels improves the accuracy of the drag interpolation by a factor of several hundred compared to the use of standard SPH bell-shaped kernels, at no additional computational expense. In order to benchmark our algorithm, we have developed a comprehensive suite of standardised, simple test problems for gas and dust mixtures: dustybox, dustywave, dustyshock, dustysedov and dustydisc, the first three of which have known analytic solutions. We present the validation of our algorithm against all of these tests. In doing so, we show that the spatial resolution criterion \Delta < cs ts is a necessary condition in all gas+dust codes that becomes critical at high drag (i.e. small stopping time ts) in order to correctly predict the dynamics. Implicit timestepping and the implementation of realistic astrophysical drag regimes are addressed in a companion paper.

Context. To form metre-sized pre-planetesimals in protoplanetary discs, growing grains have to decouple from the gas before they are accreted onto the central star during their phase of fast radial migration and thus overcome the so-called "radial-drift barrier" (often inaccurately referred to as the "metre-size barrier"). Aims. To predict the outcome of the radial motion of dust grains in protoplanetary discs whose surface density and temperature follow power-law profiles, with exponent p and q respectively. We investigate both the Epstein and the Stokes drag regimes which govern the motion of the dust. Methods. We analytically integrate the equations of motion obtained from perturbation analysis. We compare these results with those from direct numerical integration of the equations of motion. Then, using data from observed discs, we predict the fate of dust grains in real discs. Results. When a dust grain reaches the inner regions of the disc, the acceleration due to the increase of the pressure gradient is counterbalanced by the increase of the gas drag. We find that most grains in the Epstein (resp. the Stokes) regime survive their radial migration if-p+q+1/2 \leq0 (resp. if q\leq 2/3). The majority of observed discs satisfies both-p+q+ 1/2 \leq0 and q\leq 2/3: a large fraction of both their small and large grains remain in the disc, for them the radial drift barrier does not exist.

In a companion paper (Laibe & Price 2011b), we have presented an algorithm for simulating two-fluid gas and dust mixtures in Smoothed Particle Hydrodynamics (SPH). In this paper, we develop an implicit timestepping method that preserves the exact conservation of the both linear and angular momentum in the underlying SPH algorithm, but unlike previous schemes, allows the iterations to converge to arbitrary accuracy and is suited to the treatment of non- linear drag regimes. The algorithm presented in Paper I is also extended to deal with realistic astrophysical drag regimes, including both linear and non-linear Epstein and Stokes drag. The scheme is benchmarked against the test suite presented in Paper I, including i) the analytic solutions of the dustybox problem and ii) solutions of the dustywave, dustyshock, dustysedov and dustydisc obtained with explicit timestepping. We find that the implicit method is 1- 10 times faster than the explicit temporal integration when the ratio r between the the timestep and the drag stopping time is 1 < r < 1000.

In this paper we present the analytic solutions for two test problems involving two-fluid mixtures of dust and gas in an astrophysical context. The solutions provide a means of benchmarking numerical codes designed to simulate the non-linear dynamics of dusty gas. The first problem, dustybox, consists of two interpenetrating homogeneous fluids moving with relative velocity difference. We provide exact solutions to the full non-linear problem for a range of drag formulations appropriate to astrophysical fluids (i.e., various prescriptions for Epstein and Stokes drag in different regimes). The second problem, dustywave consists of the propagation of linear acoustic waves in a two-fluid gas-dust mixture. We provide the analytic solution for the case when the two fluids are interacting via a linear drag term. Both test problems are simple to set up in any numerical code and can be run with periodic boundary conditions. The solutions we derive are completely general with respect to both the dust-to-gas ratio and the amplitude of the drag coefficient. A stability analysis of waves in a gas-dust system is also presented, showing that sound waves in an astrophysical dust-gas mixture are linearly stable.

Aims: In order to understand the first stages of planet formation, when tiny grains aggregate to form planetesimals, one needs to simultaneously model grain growth, vertical settling and radial migration of dust in protoplanetary disks. In this study, we implement an analytical prescription for grain growth into a 3D two-phase hydrodynamics code to understand its effects on the dust distribution in disks. Methods: Following the analytic derivation of Stepinski & Valageas (1997), which assumes that grains stick perfectly upon collision, we implement a convenient and fast method of following grain growth in our 3D, two-phase (gas+dust) SPH code. We then follow the evolution of the size and spatial distribution of a dust population in a classical T Tauri star disk. Results: We find that the grains go through various stages of growth due to the complex interplay between gas drag, dust dynamics, and growth. Grains initially grow rapidly as they settle to the mid-plane, then experience a fast radial migration with little growth through the bulk of the disk, and finally pile-up in the inner disk where they grow more efficiently. This results in a bimodal distribution of grain sizes. Using this simple prescription of grain growth, we find that grains reach decimetric sizes in 10^5 years in the inner disk and survive the fast migration phase.

We investigate the behaviour of dust in protoplanetary disks under the action of gas drag using our 3D, two-fluid (gas+dust) SPH code. We present the evolution of the dust spatial distribution in global simulations of planetless disks as well as of disks containing an already formed planet. The resulting dust structures vary strongly with particle size and planetary gaps are much sharper than in the gas phase, making them easier to detect with ALMA than anticipated. We also find that there is a range of masses where a planet can open a gap in the dust layer whereas it doesn't in the gas disk. Our dust distributions are fed to the radiative transfer code MCFOST to compute synthetic images, in order to derive constraints on the settling and growth of dust grains in observed disks.

We present the first results of the treatment of grain growth in our 3D, two-fluid (gas+dust) SPH code describing protoplanetary disks. We implement a scheme able to reproduce the variation of grain sizes caused by a variety of physical processes and test it with the analytical expression of grain growth given by Stepinski & Valageas (1997) in simulations of a typical T Tauri disk around a one solar mass star. The results are in agreement with a turbulent growing process and validate the method. We are now able to simulate the grain growth process in a protoplanetary disk given by a more realistic physical description, currently under development. We discuss the implications of the combined effect of grain growth and dust vertical settling and radial migration on subsequent planetesimal formation.

We present the first results of the treatment of grain growth in our 3D, two-fluid (gas+dust) SPH code describing protoplanetary disks. We implement a scheme able to reproduce the variation of grain sizes caused by a variety of physical processes and test it with the analytical expression of grain growth given by Stepinski & Valageas (1997) in simulations of a typical T Tauri disk around a one solar mass star. The results are in agreement with a turbulent growing process and validate the method. We are now able to simulate the grain growth process in a protoplanetary disk given by a more realistic physical description, currently under development. We discuss the implications of the combined effect of grain growth and dust vertical settling and radial migration on subsequent planetesimal formation.