X-ray photoevaporation’s limited success in the formation of planetesimals by the streaming instability

X-ray photoevaporation’s limited success in the formation of planetesimals by the streaming instability

Authors:

Ercolano et al

Abstract:

The streaming instability is often invoked as solution to the fragmentation and drift barriers in planetesimal formation, catalysing the aggregation of dust on kyr time-scales to grow km-sized cores. However, there remains a lack of consensus on the physical mechanism(s) responsible for initiating it. One potential avenue is disc photoevaporation, wherein the preferential removal of relatively dust-free gas increases the disc metallicity. Late in the disc lifetime, photoevaporation dominates viscous accretion, creating a gradient in the depleted gas surface density near the location of the gap. This induces a local pressure maximum that collects drifting dust particles, which may then become susceptible to the streaming instability. Using a one-dimensional viscous evolution model of a disc subject to internal X-ray photoevaporation, we explore the efficacy of this process to build planetesimals. Over a range of parameters, we find that the amount of dust mass converted into planetesimals is often <1 a="" across="" an="" and="" are="" at="" au.="" be="" best="" common="" conclude="" considered="" contrast="" cores.="" debris="" discs="" driven="" far-ultra-violet="" few="" for="" formation="" in="" investigation="" large="" m="" may="" mechanism="" model="" most="" of="" our="" photoevaporation="" planetary="" rather="" recent="" relevant="" reported="" results="" similar="" spread="" tens="" than="" that="" the="" to="" we="">100 au) disc radii. The discrepancies are primarily a consequence of the different photoevaporation profiles assumed. Until observations more tightly constrain photoevaporation models, the relevance of this process to the formation of planets remains uncertain.

Redox States of Initial Atmospheres Outgassed on Rocky Planets and Planetesimals

Redox States of Initial Atmospheres Outgassed on Rocky Planets and Planetesimals

Authors:

Schaefer et al

Abstract:

The Earth and other rocky planets and planetesimals in the solar system formed through the mixing of materials from various radial locations in the solar nebula. This primordial material likely had a range of oxidation states as well as bulk compositions and volatile abundances. We investigate the oxygen fugacity produced by the outgassing of mixtures of solid meteoritic material, which approximate the primitive nebular materials. We find that the gas composition and oxygen fugacity of binary and ternary mixtures of meteoritic materials vary depending on the proportion of reduced versus oxidized material, and also find that mixtures using differentiated materials do not show the same oxygen fugacity trends as those using similarly reduced but undifferentiated materials. We also find that simply mixing the gases produced by individual meteoritic materials together does not correctly reproduce the gas composition or oxygen fugacity of the binary and ternary mixtures. We provide tabulated fits for the oxygen fugacities of all of the individual materials and binary mixtures that we investigate. These values may be useful in planetary formation models, models of volatile transport on planetesimals or meteorite parent bodies, or models of trace element partitioning during metal-silicate fractionation.

Formation of Close-in SuperEarths by Giant Impacts

Formation of Close-in Super-Earths by Giant Impacts: Effects of Initial Eccentricities and Inclinations of Protoplanets

Authors:
Matsumoto et al

Abstract:
Recent observations have revealed the eccentricity and inclination distributions of close-in super-Earths. These distributions have the potential to constrain their formation processes. In the in-situ formation scenario, the eccentricities and inclinations of planets are determined by gravitational scattering and collisions between protoplanets on the giant impact stage. We investigate the effect of the initial eccentricities and inclinations of protoplanets on the formation of close-in super-Earths. We perform N-body simulations of protoplanets in gas-free disks, changing the initial eccentricities and inclinations systematically. We find that while the eccentricities of protoplanets are well relaxed through their evolution, the inclinations are not. When the initial inclinations are small, they are not generally pumped up since scattering is less effective and collisions occur immediately after orbital crossing. On the other hand, when the initial inclinations are large, they tend to be kept large since collisional damping is less effective. Not only the resultant inclinations of planets, but also their number, eccentricities, angular momentum deficit, and orbital separations are affected by the initial inclinations of protoplanets.

Initial mass function of planetesimals formed by the streaming instability

Initial mass function of planetesimals formed by the streaming instability

Authors:

Schäfer et al

Abstract:

The streaming instability is a mechanism to concentrate solid particles into overdense filaments that undergo gravitational collapse and form planetesimals. However, it remains unclear how the initial mass function of these planetesimals depends on the box dimensions of numerical simulations. To resolve this, we perform simulations of planetesimal formation with the largest box dimensions to date, allowing planetesimals to form simultaneously in multiple filaments that can only emerge within such large simulation boxes. In our simulations, planetesimals with sizes between 80 km and several hundred kilometers form. We find that a power law with a rather shallow exponential cutoff at the high-mass end represents the cumulative birth mass function better than an integrated power law. The steepness of the exponential cutoff is largely independent of box dimensions and resolution, while the exponent of the power law is not constrained at the resolutions we employ. Moreover, we find that the characteristic mass scale of the exponential cutoff correlates with the mass budget in each filament. Together with previous studies of high-resolution simulations with small box domains, our results therefore imply that the cumulative birth mass function of planetesimals is consistent with an exponentially tapered power law with a power-law exponent of approximately -1.6 and a steepness of the exponential cutoff in the range of 0.3-0.4.

Rocky Planetesimal Formation via Fluffy Aggregates of Nanograins

Rocky Planetesimal Formation via Fluffy Aggregates of Nanograins

Authors:

Arakawa et al

Abstract:

Several pieces of evidence suggest that silicate grains in primitive meteorites are not interstellar grains but condensates formed in the early solar system. Moreover, the size distribution of matrix grains in chondrites implies that these condensates might be formed as nanometer-sized grains. Therefore, we propose a novel scenario for rocky planetesimal formation in which nanometer-sized silicate grains are produced by evaporation and recondensation events in early solar nebula, and rocky planetesimals are formed via aggregation of these nanograins. We reveal that silicate nanograins can grow into rocky planetesimals via direct aggregation without catastrophic fragmentation and serious radial drift, and our results provide a suitable condition for protoplanet formation in our solar system.

Gas Giants Need Lots of Icy Planetesimals to Form

Planet formation with envelope enrichment: new insights on planetary diversity

Authors:

Venturini et al

Abstract:

We compute, for the first time, self-consistent models of planet growth including the effect of envelope enrichment. The change of envelope metallicity is assumed to be the result of planetesimal disruption or icy pebble sublimation. We solve internal structure equations taking into account global energy conservation for the envelope to compute in-situ planetary growth. We consider different opacities and equations of state suited for a wide range of metallicities. We find that envelope enrichment speeds up the formation of gas giants. It also explains naturally the formation of low and intermediate mass objects with large fractions of H-He (~ 20 - 30 % in mass). High opacity models explain well the metallicity of the giant planets of the solar system, whereas low opacity models are suited for forming small mass objects with thick H-He envelopes and gas giants with sub-solar envelope metallicities. We find good agreement between our models and the estimated water abundance for WASP-43b. For HD 189733b, HD 209458b and WASP-12b we predict fractions of water larger than what is estimated from observations, by at least a factor ~ 2. Envelope enrichment by icy planetesimals is the natural scenario to explain the formation of a large variety of objects, ranging from mini-Neptunes, to gas giants. We predict that the total and envelope metallicity decrease with planetary mass.

How Fast Meteorites Burn up in Gas Giant Envelopes

Efficiency of Planetesimal Ablation in Giant Planetary Envelopes

Authors:

Pinhas et al

Abstract:

Observations of exoplanetary spectra are leading to unprecedented constraints on their atmospheric elemental abundances, particularly O/H, C/H, and C/O ratios. Recent studies suggest that elemental ratios could provide important constraints on formation and migration mechanisms of giant exoplanets. A fundamental assumption in such studies is that the chemical composition of the planetary envelope represents the sum-total of compositions of the accreted gas and solids during the formation history of the planet. We investigate the efficiency with which accreted planetesimals ablate in a giant planetary envelope thereby contributing to its composition rather than sinking to the core. From considerations of aerodynamic drag causing ‘frictional ablation’ and the envelope temperature structure causing ‘thermal ablation’, we compute mass ablations for impacting planetesimals of radii 30 m to 1 km for different compositions (ice to iron) and a wide range of velocities and impact angles, assuming spherical symmetry. Icy impactors are fully ablated in the outer envelope for a wide range of parameters. Even for Fe impactors substantial ablation occurs in the envelope for a wide range of sizes and velocities. For example, iron impactors of sizes below ∼0.5 km and velocities above ∼30 km/s are found to ablate by ∼ 60-80% within the outer envelope at pressures below 103 bar due to frictional ablation alone. For deeper pressures (∼107 bar), substantial ablation happens over a wider range of parameters. Therefore, our exploratory study suggests that atmospheric abundances of volatile elements in giant planets reflect their accretion history during formation.

Apparent disk-mass reduction and planetesimal formation in gravitationally unstable disks in Class 0/I YSOs

Apparent disk-mass reduction and planetesimal formation in gravitationally unstable disks in Class 0/I YSOs

Authors:

Tsukamoto et al

Abstract:

We investigate the dust structure of gravitationally unstable gas disks undergoing mass accretion from the envelope, envisioning application to Class 0/I YSOs. We compute evolution of the surface density and dust size by taking into account dust collisional growth and radial drift. We find that the dust disk quickly settles into the steady state and the dust mass in the steady-state disk decreases by a factor of 1/2 to 1/3, while the radiative flux of the dust thermal emission also decreases by a factor of 1/3 to 1/5, both compared to that for a disk with ISM dust-to-gas mass ratio and micron-sized dust. We suggest that the disk mass in the Class 0/I YSOs is underestimated by factor of 1/3 to 1/5 when it is calculated from the dust thermal emission assuming an ISM dust-to-gas mass ratio and micron-sized dust opacity, and that a larger fraction of the disks in Class 0/I YSOs is gravitationally unstable than previously considered. We derive an empirical formula for the disk-mass reduction rate, which can be used to test whether or not the disks observed in Class 0/I YSOs are gravitationally unstable. We also investigate the orbital radius of planetesimal formation rP and show that rP becomes ∼20 AU. Because rP increases as the gas surface density increases and a gravitationally unstable disk has a theoretical maximum gas surface density, rP∼20 AU is the theoretical maximum radius for planetesimal formation. We find that the dust particles migrate inwardly in the form of "pebble", and we suggest that planet formation via pebble accretion in Class 0/I phase is preferable to that in the Class II phase because the dust is supplied by envelope accretion and a significant amount of dust particles (Mdust∼10−2M⊙) pass through the disk during the Class 0/I phase.

FAILED GROWTH AT THE BOUNCING BARRIER IN PLANETESIMAL FORMATION

FAILED GROWTH AT THE BOUNCING BARRIER IN PLANETESIMAL FORMATION

Authors:

Kruss et al

Abstract:

In laboratory experiments, we studied collisions of ensembles of compact (filling factor = 0.33) millimeter dust aggregates composed of micrometer quartz grains. We used cylindrical aggregates, triangular aggregates, square aggregates, and rectangular aggregates. Ensembles of equal size aggregates as well as ensembles with embedded larger aggregates were studied. The typical collision velocities are 10–20 mm s−1. High spatial and temporal resolution imaging unambiguously shows that individual collisions lead to sticking with a high probability of 20%. This leads to connected clusters of aggregates. The contact areas between two aggregates increase with collision velocity. However, this cluster growth is only temporary, as subsequent collisions of aggregates and clusters eventually lead to the detachment of all aggregates from a cluster. The contacts are very fragile as aggregates cannot be compressed further or fragment under our experimental conditions to enhance the contact stability. Therefore, the evolution of the ensemble always leads back to a distribution of individual aggregates of initial size. This supports and extends earlier experiments showing that a bouncing barrier in planetesimal formation would be robust against shape and size variations.

Cometary ices in forming protoplanetary disc midplanes

Cometary ices in forming protoplanetary disc midplanes

Authors:

Drozdovskaya et al

Abstract:

Low-mass protostars are the extrasolar analogues of the natal Solar system. Sophisticated physicochemical models are used to simulate the formation of two protoplanetary discs from the initial prestellar phase, one dominated by viscous spreading and the other by pure infall. The results show that the volatile prestellar fingerprint is modified by the chemistry en route into the disc. This holds relatively independent of initial abundances and chemical parameters: physical conditions are more important. The amount of CO2 increases via the grain-surface reaction of OH with CO, which is enhanced by photodissociation of H2O ice. Complex organic molecules are produced during transport through the envelope at the expense of CH3OH ice. Their abundances can be comparable to that of methanol ice (few per cent of water ice) at large disc radii (R > 30 au). Current Class II disc models may be underestimating the complex organic content. Planet population synthesis models may underestimate the amount of CO2 and overestimate CH3OH ices in planetesimals by disregarding chemical processing between the cloud and disc phases. The overall C/O and C/N ratios differ between the gas and solid phases. The two ice ratios show little variation beyond the inner 10 au and both are nearly solar in the case of pure infall, but both are subsolar when viscous spreading dominates. Chemistry in the protostellar envelope en route to the protoplanetary disc sets the initial volatile and prebiotically significant content of icy planetesimals and cometary bodies. Comets are thus potentially reflecting the provenances of the midplane ices in the solar nebula.

Watching the Variation in Composition of the Debris Disk Around White Dwarf WD 1145+017

Spectroscopic Evolution of Disintegrating Planetesimals: Minutes to Months Variability in the Circumstellar Gas Associated with WD 1145+017

Authors:

Redfield et al

Abstract:

With the recent discovery of transiting planetary material around WD 1145+017, a critical target has been identified that links the evolution of planetary systems with debris disks and their accretion onto the star. We present a series of observations, five epochs over a year, taken with Keck and the VLT, which for the first time show variability of circumstellar absorption in the gas disk surrounding WD 1145+017 on timescales of minutes to months. Circumstellar absorption is measured in more than 250 lines of 14 ions among ten different elements associated with planetary composition, e.g., O, Mg, Ca, Ti, Cr, Mn, Fe, Ni. Broad circumstellar gas absorption with a velocity spread of 225 km/s is detected, but over the course of a year blue shifted absorption disappears while redshifted absorption systematically increases. A correlation of equivalent width and oscillator strength indicates that the gas is not highly optically thick (median tau approximately 2). We discuss simple models of an eccentric disk coupled with magnetospheric accretion to explain the basic observed characteristics of these high resolution and high signal-to-noise observations. Variability is detected on timescales of minutes in the two most recent observations, showing a loss of redshifted absorption for tens of minutes, coincident with major transit events and consistent with gas hidden behind opaque transiting material. This system currently presents a unique opportunity to learn how the gas causing the spectroscopic, circumstellar absorption is associated with the ongoing accretion evidenced by photospheric contamination, as well as the transiting planetary material detected in photometric observations.

Planetesimal Formation by Gravitational Instability of a Porous-Dust Disk

Planetesimal Formation by Gravitational Instability of a Porous-Dust Disk

Authors:

Michikoshi et al

Abstract:

Recently it is proposed that porous icy dust aggregates are formed by pairwise accretion of dust aggregates beyond the snowline. We calculate the equilibrium random velocity of porous dust aggregates taking into account mutual gravitational scattering, collisions, gas drag, and turbulent stirring and scattering. We find that the disk of porous dust aggregates becomes gravitationally unstable as they evolve through gravitational compression in the minimum-mass solar nebula model for a reasonable range of turbulence strength, which leads to rapid formation of planetesimals.

Three-lane signatures of planets in planetesimal disks

Three-lane signatures of planets in planetesimal disks

Authors:

Demidova et al

Abstract:

In massive numerical experiments we show that a planet embedded in a planetesimal disk induces a characteristic three-lane "planetosignature" representing a pattern of three stellar-centric rings, one bright coorbital with the planet, and two dark gaps in the radial distribution of the particles. The gaps correspond to orbital resonances 2/1 and 1/2 with the planet. This theoretical prediction is in perfect agreement with recent ALMA observations of the disk of HL Tau.

Dust Evolution and the Formation of Planetesimals

Dust Evolution and the Formation of Planetesimals

Authors:

Birnstiel et al

Abstract:

The solid content of circumstellar disks is inherited from the interstellar medium: dust particles of at most a micrometer in size. Protoplanetary disks are the environment where these dust grains need to grow at least 13 orders of magnitude in size. Our understanding of this growth process is far from complete, with different physics seemingly posing obstacles to this growth at various stages. Yet, the ubiquity of planets in our galaxy suggests that planet formation is a robust mechanism. This chapter focuses on the earliest stages of planet formation, the growth of small dust grains towards the gravitationally bound "planetesimals", the building blocks of planets. We will introduce some of the key physics involved in the growth processes and discuss how they are expected to shape the global behavior of the solid content of disks. We will consider possible pathways towards the formation of larger bodies and conclude by reviewing some of the recent observational advances in the field.

Carbon to oxygen ratios in extrasolar planetesimals

Carbon to oxygen ratios in extrasolar planetesimals

Authors:

Wilson et al

Abstract:

Observations of small extrasolar planets with a wide range of densities imply a variety of planetary compositions and structures. Currently, the only technique to measure the bulk composition of extrasolar planetary systems is the analysis of planetary debris accreting onto white dwarfs, analogous to abundance studies of meteorites. We present measurements of the carbon and oxygen abundances in the debris of planetesimals at ten white dwarfs observed with the Hubble Space Telescope, along with C/O ratios of debris in six systems with previously reported abundances. We find no evidence for carbon-rich planetesimals, with C/O less than 0.8 by number in all 16 systems. Our results place an upper limit on the occurrence of carbon-rich systems at less than 17 percent with a 2 sigma confidence level. The range of C/O of the planetesimals is consistent with that found in the Solar System, and appears to follow a bimodal distribution: a group similar to the CI chondrites, with log(less than C/O greater than)=-0.92, and oxygen-rich objects with C/O less than or equal to that of the bulk Earth. The latter group may have a higher mass fraction of water than the Earth, increasing their relative oxygen abundance.

Extraordinary Claim: the Earth Formed Volatile Free

United Theory of Planet Formation (I): Tandem Regime

Authors:

Ebisuzaki et al

Abstract:

We have obtained a steady-state, 1-D model of the accretion disk of a protostar taking into account the magneto-rotational instability (MRI). We find that the disk is divided into an outer turbulent region (OTR), a MRI suppressed region (MSR), and an inner turbulent region (ITR). The outer turbulent region is fully turbulent because of MRI. However, in the range, r_{out} (= 8 - 60 AU) from the central star, MRI is suppressed around the midplane of the gas disk and a quiet area without turbulence appears, because the degree of ionization of gas becomes low enough. The disk becomes fully turbulent again in the range r in (= 0.2 - 1 AU), which is called the inner turbulent region, because the midplane temperature become high enough (\gt 1000 K) due to gravitational energy release.

Planetesimals are formed through gravitational instability at the two distinct sites, outer and inner MRI fronts (the boundaries between the MRI suppressed region (MSR) and the outer and inner turbulent regions), because of the radial concentration of the solid particles. At the outer MRI front, icy particles grow through low-velocity collisions into porous aggregates with low densities. They eventually undergo gravitational instability to form icy planetesimals. On the other hand, rocky particles accumulate at the inner MRI front, since their drift velocities turn outward due to the local maximum in gas pressure. They undergo gravitational instability in a sub-disk of pebbles to form rocky planetesimals at the inner MRI front.

The tandem regime is consistent with the ABEL model, in which the Earth was initially formed as a completely volatile-free planet. The water and other volatile elements came later through the accretion of icy particles by the occasional scatterings in the outer regions.

Insights into planet formation from debris disks

Insights into planet formation from debris disks: I. The solar system as an archetype for planetesimal evolution

Authors:

Matthews et al

Abstract:

Circumstellar disks have long been regarded as windows into planetary systems. The advent of high sensitivity, high resolution imaging in the submillimetre where both the solid and gas components of disks can be detected opens up new possibilities for understanding the dynamical histories of these systems and therefore, a better ability to place our own solar system, which hosts a highly evolved debris disk, in context. Comparisons of dust masses from protoplanetary and debris disks have revealed a stark downturn in mass in millimetre-sized grains around a stellar age of 10 Myr, ostensibly in the "transition disk" phase, suggesting a period of rapid accretion of such grains onto planetesimals. This rapid formation phase is in keeping with radionucleide studies of Kuiper Belt Objects in the solar system. Importantly, this suggests that any thermal gradients in the gas of disks of this era will be "frozen in" to the planetesimals as they rapidly accrete from the solids and ices in their vicinity. Measurements of radial gradients in thermal tracers such as DHO, DCN and other tracers can therefore provide insight into the nascent solar system's abudances. In studies of dynamical evolution of the solar system, it is tacitly assumed that such abundances can reveal the location of formation for bodies now found in the asteroid belt and Kuiper belt. Similarly, evidence of gas detected from collisional evolution in young debris disks could potentially reveal how rapidly objects have dynamically evolved in those systems, most of which will be significantly younger than the solar system.

Constraints on Planetesimal Collision Models in Debris Disks

Constraints on Planetesimal Collision Models in Debris Disks

Authors:

MacGregor et al

Abstract:

Observations of debris disks offer a window into the physical and dynamical properties of planetesimals in extrasolar systems through the size distribution of dust grains. In particular, the millimeter spectral index of thermal dust emission encodes information on the grain size distribution. We have made new VLA observations of a sample of seven nearby debris disks at 9 mm, with 3" resolution and ∼5 μJy/beam rms. We combine these with archival ATCA observations of eight additional debris disks observed at 7 mm, together with up-to-date observations of all disks at (sub)millimeter wavelengths from the literature to place tight constraints on the millimeter spectral indices and thus grain size distributions. The analysis gives a weighted mean for the slope of the power law grain size distribution, n(a)∝a−q, of ⟨q⟩=3.36±0.02, with a possible trend of decreasing q for later spectral type stars. We compare our results to a range of theoretical models of collisional cascades, from the standard self-similar, steady-state size distribution (q=3.5) to solutions that incorporate more realistic physics such as alternative velocity distributions and material strengths, the possibility of a cutoff at small dust sizes from radiation pressure, as well as results from detailed dynamical calculations of specific disks. Such effects can lead to size distributions consistent with the data, and plausibly the observed scatter in spectral indices. For the AU Mic system, the VLA observations show clear evidence of a highly variable stellar emission component; this stellar activity obviates the need to invoke the presence of an asteroid belt to explain the previously reported compact millimeter source in this system.

Resolving the Planetesimal Belt of HR 8799 with ALMA

Resolving the Planetesimal Belt of HR 8799 with ALMA

Authors:

Booth et al

Abstract:

The star HR 8799 hosts one of the largest known debris discs and at least four giant planets. Previous observations have found evidence for a warm belt within the orbits of the planets, a cold planetesimal belt beyond their orbits and a halo of small grains. With the infrared data, it is hard to distinguish the planetesimal belt emission from that of the grains in the halo. With this in mind, the system has been observed with ALMA in band 6 (1.34 mm) using a compact array format. These observations allow the inner edge of the planetesimal belt to be resolved for the first time. A radial distribution of dust grains is fitted to the data using an MCMC method. The disc is best fit by a broad ring between 145+12−12 AU and 429+37−32 AU at an inclination of 40+5−6{\deg} and a position angle of 51+8−8{\deg}. A disc edge at ~145 AU is too far out to be explained simply by interactions with planet b, requiring either a more complicated dynamical history or an extra planet beyond the orbit of planet b.

Collision velocity of dust grains in self-gravitating protoplanetary discs

Collision velocity of dust grains in self-gravitating protoplanetary discs

Authors:

Booth et al

Abstract:

We have conducted the first comprehensive numerical investigation of the relative velocity distribution of dust particles in self-gravitating protoplanetary discs with a view to assessing the viability of planetesimal formation via direct collapse in such environments. The viability depends crucially on the large sizes that are preferentially collected in pressure maxima produced by transient spiral features (Stokes numbers, St∼1); growth to these size scales requires that collision velocities remain low enough that grain growth is not reversed by fragmentation. We show that, for a single sized dust population, velocity driving by the disc's gravitational perturbations is only effective for St>3, while coupling to the gas velocity dominates otherwise. We develop a criterion for understanding this result in terms of the stopping distance being of order the disc scale height. Nevertheless, the relative velocities induced by differential radial drift in multi-sized dust populations are too high to allow the growth of silicate dust particles beyond St∼10−2 or 10−1 (10cm to m sizes at 30au), such Stokes numbers being insufficient to allow concentration of solids in spiral features. However, for icy solids (which may survive collisions up to several 10ms−1), growth to St∼1 (10m size) may be possible beyond 30au from the star. Such objects would be concentrated in spiral features and could potentially produce larger icy planetesimals/comets by gravitational collapse. These planetesimals would acquire moderate eccentricities and remain unmodified over the remaining lifetime of the disc.