Planetesimal Formation Processes Research Articles

A self-consistent model for dust settling and the vertical shear instability in protoplanetary disks

Abstract The spatial distribution of dust particles in protoplanetary disks affects dust evolution and planetesimal formation processes. The vertical shear instability (VSI) is one of the candidate hydrodynamic mechanisms that can generate turbulence in the outer disk region and affect dust diffusion. Turbulence driven by the VSI has a predominant vertical motion that can prevent dust settling. On the other hand, the dust distribution controls the spatial distribution of the gas cooling rate, thereby affecting the strength of VSI-driven turbulence. Here, we present a semi-analytic model that determines the vertical dust distribution and the strength of VSI-driven turbulence in a self-consistent manner. The model uses an empirical formula for the vertical diffusion coefficient in VSI-driven turbulence obtained from our recent hydrodynamical simulations. The formula returns the vertical diffusion coefficient as a function of the vertical profile of the cooling rate, which is determined by the vertical dust distribution. We use this model to search for an equilibrium vertical dust profile where settling balances with turbulent diffusion for a given maximum grain size. We find that if the grains are sufficiently small, there exists a stable equilibrium dust distribution where VSI-driven turbulence is sustained at a level of αz ∼ 10−3, where αz is the dimensionless vertical diffusion coefficient. However, as the maximum grain size increases, the equilibrium solution vanishes because the VSI can no longer stop the settling of the grains. This runaway settling may explain highly settled dust rings found in the outer part of some protoplanetary disks.

Streaming instability on different scales – I. Planetesimal mass distribution variability

ABSTRACT We present numerical simulations of dust clumping and planetesimal formation initiated by the streaming instability (SI) with self-gravity. We examine the variability in the planetesimal formation process by employing simulation domains with large radial and azimuthal extents and a novel approach of re-running otherwise identical simulations with different random initializations of the dust density field. We find that the planetesimal mass distribution and the total mass of dust that is converted into planetesimals can vary substantially between individual small simulations and within the domains of larger simulations. Our results show that the non-linear nature of the developed SI introduces substantial variability in the planetesimal formation process that has not been previously considered and suggests larger scale dynamics may affect the process.

Modelling the distributions of white dwarf atmospheric pollution: a low Mg abundance for accreted planetesimals?

Abstract The accretion of planetesimals onto white dwarf atmospheres allows determination of the composition of this polluting material. This composition is usually inferred from observed pollution levels by assuming it originated from a single body. This paper instead uses a stochastic model wherein polluting planetesimals are chosen randomly from a mass distribution, finding that the single body assumption is invalid in ${>20\%}$ of cases. Planetesimal compositions are modelled assuming parent bodies that differentiated into core, mantle and crust components. Atmospheric levels of Ca, Mg and Fe in the model are compared to a sample of 230 DZ white dwarfs for which such pollution is measured. A good fit is obtained when each planetesimal has its core, mantle and crust fractions chosen independently from logit-normal distributions which lead to average mass fractions of fCru = 0.15, fMan = 0.49 and fCor = 0.36. However, achieving this fit requires a factor 4 depletion of Mg relative to stellar material. This depletion is unlikely to originate in planetesimal formation processes, but might occur from heating while the star is on the giant branch. Alternatively the accreted material has stellar abundance, and either the inferred low Mg abundance was caused by an incorrect assumption that Mg sinks slower than Ca and Fe, or there are unmodelled biases in the observed sample. Finally, the model makes predictions for the timescale on which the observed pollutant composition varies, which should be the longer of the sinking and disc timescales, implying variability on decadal timescales for DA white dwarfs.

Crater Density Predictions for New Horizons Flyby Target 2014 MU69

Abstract In preparation for the 2019 January 1 encounter between the New Horizons spacecraft and the Kuiper Belt object 2014 MU69, we provide estimates of the expected impact crater surface density on the Kuiper Belt object. Using the observed crater fields on Charon and Pluto down to the resolution limit of the 2015 New Horizons flyby of those bodies and estimates of the orbital distribution of the crater-forming projectiles, we calculate the number of craters per unit area formed as a function of the time a surface on 2014 MU69 has been exposed to bombardment. We find that if the shallow crater size distribution from roughly 1–15 km exhibited on Pluto and Charon is indeed due to the sizes of Kuiper Belt projectiles, 2014 MU69 should exhibit a surface that is only lightly cratered below 1 km scale, despite being bombarded for ∼4 billion years. Its surface should therefore be more clearly indicative of its accretionary environment. In addition, this object may be the first observed for which the majority of the bombardment is from exogenic projectiles moving at less than or near the speed of sound in the target materials, implying morphologies more akin to secondary craters elsewhere in the solar system. Lastly, if the shallow Kuiper Belt size distribution implied from the Pluto and Charon imaging is confirmed at 2014 MU69, then we conclude that this size distribution is a preserved relic of its state ≃4.5 Gyr ago and provides a direct constraint on the planetesimal formation process itself.

Spectroscopy and thermal modelling of the first interstellar object 1I/2017 U1 ‘Oumuamua

During the formation and evolution of the Solar System, significant numbers of cometary and asteroidal bodies were ejected into interstellar space$^{1,2}$. It can be reasonably expected that the same happened for planetary systems other than our own. Detection of such Inter- stellar Objects (ISOs) would allow us to probe the planetesimal formation processes around other stars, possibly together with the effects of long-term exposure to the interstellar medium. 1I/2017 U1 'Oumuamua is the first known ISO, discovered by the Pan-STARRS1 telescope in October 2017$^3$.The discovery epoch photometry implies a highly elongated body with radii of $\sim 200 \times 20$ m when a comet-like geometric albedo of 0.04 is assumed. Here we report spectroscopic characterisation of 'Oumuamua, finding it to be variable with time but similar to organically rich surfaces found in the outer Solar System. The observable ISO population is expected to be dominated by comet-like bodies in agreement with our spectra, yet the reported inactivity implies a lack of surface ice. We show this is consistent with predictions of an insulating mantle produced by long-term cosmic ray exposure. An internal icy composition cannot therefore be ruled out by the lack of activity, even though 'Oumuamua passed within 0.25 au of the Sun.

Interpreting the extended emission around three nearby debris disc host stars

Cool debris discs are a relic of the planetesimal formation process around their host star, analogous to the solar system's Edgeworth-Kuiper belt. As such, they can be used as a proxy to probe the origin and formation of planetary systems like our own. The Herschel Open Time Key Programmes "DUst around NEarby Stars" (DUNES) and "Disc Emission via a Bias-free Reconnaissance in the Infrared/Submillimetre" (DEBRIS) observed many nearby, sun-like stars at far-infrared wavelengths seeking to detect and characterize the emission from their circumstellar dust. Excess emission attributable to the presence of dust was identified from around $\sim$ 20% of stars. Herschel's high angular resolution ($\sim$ 7" FWHM at 100 $\mu$m) provided the capacity for resolving debris belts around nearby stars with radial extents comparable to the solar system (50 to 100 au). As part of the DUNES and DEBRIS surveys, we obtained observations of three debris disc stars, HIP 22263 (HD 30495), HIP 62207 (HD 110897), and HIP 72848 (HD 131511), at far-infrared wavelengths with the Herschel PACS instrument. Combining these new images and photometry with ancilliary data from the literature, we undertook simultaneous multi-wavelength modelling of the discs' radial profiles and spectral energy distributions using three different methodologies: single annulus, modified black body, and a radiative transfer code. We present the first far-infrared spatially resolved images of these discs and new single-component debris disc models. We characterize the capacity of the models to reproduce the disc parameters based on marginally resolved emission through analysis of two sets of simulated systems (based on the HIP 22263 and HIP 62207 data) with the noise levels typical of the Herschel images. We find that the input parameter values are recovered well at noise levels attained in the observations presented here.

N-BODY SIMULATION OF PLANETESIMAL FORMATION THROUGH GRAVITATIONAL INSTABILITY OF A DUST LAYER IN LAMINAR GAS DISK

We investigate the formation process of planetesimals from the dust layer by the gravitational instability in the gas disk using local $N$-body simulations. The gas is modeled as a background laminar flow. We study the formation process of planetesimals and its dependence on the strength of the gas drag. Our simulation results show that the formation process is divided into three stages qualitatively: the formation of wake-like density structures, the creation of planetesimal seeds, and their collisional growth. The linear analysis of the dissipative gravitational instability shows that the dust layer is secularly unstable although Toomre's $Q$ value is larger than unity. However, in the initial stage, the growth time of the gravitational instability is longer than that of the dust sedimentation and the decrease in the velocity dispersion. Thus, the velocity dispersion decreases and the disk shrinks vertically. As the velocity dispersion becomes sufficiently small, the gravitational instability finally becomes dominant. Then wake-like density structures are formed by the gravitational instability. These structures fragment into planetesimal seeds. The seeds grow rapidly owing to mutual collisions.

N-BODY SIMULATION OF PLANETESIMAL FORMATION THROUGH GRAVITATIONAL INSTABILITY AND COAGULATION. II. ACCRETION MODEL

The gravitational instability of a dust layer is one of the scenarios for planetesimal formation. If the density of a dust layer becomes sufficiently high as a result of the sedimentation of dust grains toward the midplane of a protoplanetary disk, the layer becomes gravitationally unstable and spontaneously fragments into planetesimals. Using a shearing box method, we performed local $N$-body simulations of gravitational instability of a dust layer and subsequent coagulation without gas and investigated the basic formation process of planetesimals. In this paper, we adopted the accretion model as a collision model. A gravitationally bound pair of particles is replaced by a single particle with the total mass of the pair. This accretion model enables us to perform long-term and large-scale calculations. We confirmed that the formation process of planetesimals is the same as that in the previous paper with the rubble pile models. The formation process is divided into three stages: the formation of non-axisymmetric structures, the creation of planetesimal seeds, and their collisional growth. We investigated the dependence of the planetesimal mass on the simulation domain size. We found that the mean mass of planetesimals formed in simulations is proportional to $L_y^{3/2}$, where $L_y$ is the size of the computational domain in the direction of rotation. However, the mean mass of planetesimals is independent of $L_x$, where $L_x$ is the size of the computational domain in the radial direction if $L_x$ is sufficiently large. We presented the estimation formula of the planetesimal mass taking into account the simulation domain size.

Detailed characterization of stellar high energy (FUV/EUV/X-ray) radiation fields during protoplanetary system formation

AbstractYoung stars undergoing the conversion of pre-main-sequence circumstellar disks into protoplanetary systems are strong sources of high energy (FUV/EUV/X-ray) radiation that controls the physical and chemical processes in their circumstellar environment out to hundreds of AU from the star. The high energy radiation resulting from magnetic activity and accretion onto the central star controls the thermal structure of disks, the formation process of planetesimals, and the photoexcitation and photoionization of protoplanets and young planetary atmospheres. Modeling of the dust and gas evolution requires an accurate understanding of the local radiation field throughout the ultraviolet (UV) and X-ray spectral regions, even those parts of the spectrum that are impossible to observe from Earth.Our current research efforts are directed towards developing a better understanding of UV (using HST and FUSE) and X-ray (using Chandra, XMM-Newton, and Swift) stellar activity and the resulting radiation fields during pre-main-sequence evolution from ages of a few to several hundred million years. These studies include extensive UV and X-ray spectral sampling of individual stars in nearby star formation regions and the various moving groups of the Local Association, including our HST Cycle 17 Large Project (GO-11616), which is using 111 HST orbits to observe 32 T Tauri stars with the COS UV spectrograph. Most young stars are well over 100 pc from the Sun and are consequently hard to observe in the UV and X-ray regions at even moderate spectral resolution. However, members of the Local Association, whose ages range from 7 Myr to a few hundred Myr, surround the Sun at distances of 50 pc or less and permit the detailed study of the later stages of the early evolution of stellar activity when gas giant and terrestrial protoplanets are forming. We illustrate our methodology using the 12 Myr old early-M dwarf AU Mic, which possesses a striking dust debris disk, as an example.

N‐Body Simulation of Planetesimal Formation through Gravitational Instability of a Dust Layer

We performed N-body simulations of a dust layer without a gas component and examined the formation process of planetesimals. We found that the formation process of planetesimals can be divided into three stages: the formation of non-axisymmetric wake-like structures, the creation of aggregates, and the collisional growth of the aggregates. Finally, a few large aggregates and many small aggregates are formed. The mass of the largest aggregate is larger than the mass predicted by the linear perturbation theory. We examined the dependence of system parameters on the planetesimal formation. We found that the mass of the largest aggregates increase as the size of the computational domain increases. However the ratio of the aggregate mass to the total mass $M_\mathrm{aggr}/M_\mathrm{total}$ is almost constant $0.8-0.9$. The mass of the largest aggregate increases with the optical depth and the Hill radius of particles.