No Arabic abstract
Context. Circumstellar disks are known to contain a significant mass in dust ranging from micron to centimeter size. Meteorites are evidence that individual grains of those sizes were collected and assembled into planetesimals in the young solar system. Aims. We assess the efficiency of dust collection of a swarm of non-drifting planetesimals {rev with radii ranging from 1 to $10^3$,km and beyond. Methods. We calculate the collision probability of dust drifting in the disk due to gas drag by planetesimal accounting for several regimes depending on the size of the planetesimal, dust, and orbital distance: the geometric, Safronov, settling, and three-body regimes. We also include a hydrodynamical regime to account for the fact that small grains tend to be carried by the gas flow around planetesimals. Results. We provide expressions for the collision probability of dust by planetesimals and for the filtering efficiency by a swarm of planetesimals. For standard turbulence conditions (i.e., a turbulence parameter $alpha=10^{-2}$), filtering is found to be inefficient, meaning that when crossing a minimum-mass solar nebula (MMSN) belt of planetesimals extending between 0.1 AU and 35 AU most dust particles are eventually accreted by the central star rather than colliding with planetesimals. However, if the disk is weakly turbulent ($alpha=10^{-4}$) filtering becomes efficient in two regimes: (i) when planetesimals are all smaller than about 10 km in size, in which case collisions mostly take place in the geometric regime; and (ii) when planetary embryos larger than about 1000 km in size dominate the distribution, have a scale height smaller than one tenth of the gas scale height, and dust is of millimeter size or larger in which case most collisions take place in the settling regime. These two regimes have very different properties: we find that the local filtering efficiency $x_{filter,MMSN}$ scales with $r^{-7/4}$ (where $r$ is the orbital distance) in the geometric regime, but with $r^{-1/4}$ to $r^{1/4}$ in the settling regime. This implies that the filtering of dust by small planetesimals should occur close to the central star and with a short spread in orbital distances. On the other hand, the filtering by embryos in the settling regime is expected to be more gradual and determined by the extent of the disk of embryos. Dust particles much smaller than millimeter size tend only to be captured by the smallest planetesimals because they otherwise move on gas streamlines and their collisions take place in the hydrodynamical regime. Conclusions. Our results hint at an inside-out formation of planetesimals in the infant solar system because small planetesimals in the geometrical limit can filter dust much more efficiently close to the central star. However, even a fully-formed belt of planetesimals such as the MMSN only marginally captures inward-drifting dust and this seems to imply that dust in the protosolar disk has been filtered by planetesimals even smaller than 1 km (not included in this study) or that it has been assembled into planetesimals by other mechanisms (e.g., orderly growth, capture into vortexes). Further refinement of our work concerns, among other things: a quantitative description of the transition region between the hydro and settling regimes; an assessment of the role of disk turbulence for collisions, in particular in the hydro regime; and the coupling of our model to a planetesimal formation model.
In circumstellar discs, collisional grinding of planetesimals produces second-generation dust. While it remains unclear whether this ever becomes a major component of the total dust content, the presence of such dust, and potentially the substructure within, it can be used to explore a discs physical conditions. A perturbing planet produces nonaxisymmetric structures and gaps in the dust, regardless of its origin. The dynamics of planetesimals, however, will be very different than that of small dust grains due to weaker gas interactions. Therefore, planetesimal collisions could create dusty disc structures that would not exist otherwise. In this work, we use N-body simulations to investigate the collision rate profile of planetesimals near mean-motion resonances. We find that a distinct bump or dip feature is produced in the collision profile, the presence of which depends on the libration width of the resonance and the separation between the peri- and apocenter distances of the edges of the resonance. The presence of one of these two features depends on the mass and eccentricity of the planet. Assuming that the radial dust emission traces the planetesimal collision profile, the presence of a bump or dip feature in the dust emission at the 2:1 mean-motion resonance can constrain the orbital properties of the perturbing planet. This assumption is valid, so long as radial drift does not play a significant role during the collisional cascade process. Under this assumption, these features in the dust emission should be marginally observable in nearby protoplanetary disks with ALMA.
We present a study on the formation of planetary systems around low mass stars similar to Trappist-1, through the accretion of either planetesimals or pebbles. The aim is to determine if the currently observed systems around low mass stars favour one scenario over the other. We ran numerous N-body simulations, coupled to a thermally evolving viscous disc model, including prescriptions for planet migration and photoevaporation. We examine the differences between the pebble and planetesimal accretion scenarios, but also look at the influences of disc mass, planetesimal size, and the percentage of solids locked up within pebbles. When comparing the resulting planetary systems to Trappist-1, we find that a wide range of initial conditions for both accretion scenarios can form planetary systems similar to Trappist-1, in terms of planet mass, periods, and resonant configurations. Typically these planets formed exterior to the water iceline and migrated in resonant convoys to close to the central star. When comparing the planetary systems formed from pebbles to those formed from planetesimals, we find a large number of similarities, including average planet masses, eccentricities, inclinations and period ratios. One major difference was that of the water content of the planets. When including the effects of ablation and full recycling of the planets envelope with the disc, planets formed from pebbles were extremely dry, whilst those formed from planetesimals were extremely wet. If the water content is not fully recycled and instead falls to the planets core, or if ablation of the water is neglected, then the planets formed from pebbles are extremely wet, similar to those formed from planetesimals. Should the water content of the Trappist-1 planets be determined accurately, this could point to a preferred formation pathway for planetary systems, or to specific physics that may be at play.
Pebble accretion is a new mechanism to quickly grow the cores of planets. In pebble accretion, gravity and gas drag conspire to yield large collisional cross sections for small particles in protoplanetary disks. However, before pebble accretion commences, aerodynamical deflection may act to prevent planetesimals from becoming large, because particles tend to follow gas streamlines. We derive the planetesimal radius where pebble accretion is initiated and determine the growth timescales of planetesimals by sweepup of small particles. We obtain the collision efficiency factor as the ratio of the numerically-obtained collisional cross section to the planetesimal surface area, from which we obtain the growth timescales. Integrations are conducted in the potential flow limit (steady, inviscid) and in the Stokes flow regime (steady, viscid). Only particles of stopping time $t_s ll t_X$ where $t_Xapprox10^3$ s experience aerodynamic deflection. Even in that case, the planetesimals gravity always ensures positive collision factors. The maximum growth timescale occurs typically at around $Rapprox100 mathrm{km}$, but is less for colder disks, corresponding to interactions shifting to the Safronov focusing regime. For particles $t_s gg t_X$ pebble accretion commences only after this phase and is characterized by a steep drop in growth timescales. Consequently, at distances beyond ~10 AU sweepup growth timescales are always longer than $10$ Myr, while in the inner disk (~<3 AU) the viability of the sweepup scenario is determined by the outcome of pebble-planetesimal collisions in the geometric regime. We present analytical fits for the collision efficiency factors and the minimum planetesimal size needed for pebble accretion. (Abridged)
As a test bed for the growth of protoplanetary bodies in a turbulent circumstellar disk we examine the fate of a boulder using direct numerical simulations of particle seeded gas flowing around it. We provide an accurate description of the flow by imposing no-slip and non-penetrating boundary conditions on the boulder surface using the immersed boundary method pioneered by Peskin (2002). Advected by the turbulent disk flow, the dust grains collide with the boulder and we compute the probability density function (PDF) of the normal component of the collisional velocity. Through this examination of the statistics of collisional velocities we test the recently developed concept of collisional fusion which provides a physical basis for a range of collisional velocities exhibiting perfect sticking. A boulder can then grow sufficiently rapidly to settle into a Keplerian orbit on disk evolution time scales.
We investigate the influence of impacts of large planetesimals and small planetary embryos on the early Martian surface on the hydrodynamic escape of an early steam atmosphere that is exposed to the high soft X-ray and EUV flux of the young Sun. Impact statistics in terms of number, masses, velocities, and angles of asteroid impacts onto the early Mars are determined via n-body integrations. Based on these statistics, smoothed particle hydrodynamics (SPH) simulations result in estimates of energy transfer into the planetary surface material and according surface heating. For the estimation of the atmospheric escape rates we applied a soft X-ray and EUV absorption model and a 1-D upper atmosphere hydrodynamic model to a magma ocean-related catastrophically outgassed steam atmosphere with surface pressure values of 52 bar H2O and 11 bar CO2. The estimated impact rates and energy deposition onto an early Martian surface can account for substantial heating. The energy influx and conversion rate into internal energy is most likely sufficient to keep a shallow magma ocean liquid for an extended period of time. Higher surface temperatures keep the outgassed steam atmosphere longer in vapor form and therefore enhance its escape to space within about 0.6 Myr after its formation.