Planetesimal Formation Research Articles

The growth of super-large pre-planetary pebbles to an impact erosion limit

Abstract The early evolution of dust in protoplanetary disks is dominated by sticking collisions. However, this initial phase of particle growth faces constraints, notably from destructive encounters. To find the maximum particle size achievable, we studied collisional processes during a prolonged microgravity experiment aboard a suborbital flight. Specifically, we describe an impact erosion limit. We observed individual basalt beads, each measuring 0.5 mm in diameter, colliding with and either eroding or adhering to a cluster several centimetres in size. This cluster, formed from tribocharged particles, simulates an electrostatic growth phase that surpasses the classical bouncing barrier. We found a threshold velocity of about 0.5 m s−1, which separates additive and erosive impacts of individual beads. Numerical simulations of grains impacting clusters, for both low and high charge constituents, corroborate the experimental findings of surface erosion within the observed velocity range. This specific velocity threshold supports the potential formation of pebbles several centimetres in size within protoplanetary disks. Such dimensions place these pebbles well into a regime in which hydrodynamic interactions might facilitate the formation of planetesimals.

Clusters of tribocharged dust aggregates as pebbles in protoplanetary disks

In recent years, the tribocharging of colliding and bouncing submillimeter (submm) particles has been studied as a possible mechanism promoting the formation of large pebbles on centimeter (cm) to decimeter (dm) scales in protoplanetary disks. Here, we observe, for the first time, that it is not only monolithic, spherical particles, but also real dust aggregates, that become tribocharged and end up forming large clusters. For aggregates of ∼ 0.4,mm consisting of rm ∼ 1,micrometer (µm) sized dust, we determined net charge densities up to 10^-7,C/m² during our drop tower experiments. These charged aggregates form compact clusters up to 2,cm in size via collisions with other clusters and aggregates at collision velocities on the order of 1,cm/s. Size and speed are the only lower limits for growth, currently set by the limits of the experiment. However, these clusters already form under conditions that are well beyond the expected transition to bouncing for uncharged aggregates and clusters. Our findings further support the idea that collisional charging can leapfrog the traditional bouncing barrier and form larger clusters that then serve as large pebbles. These cm-sized clusters are more susceptible to further evolutionary steps via particle trapping, concentration, and planetesimal formation.

REsolved ALMA and SMA Observations of Nearby Stars (REASONS)

Context. Planetesimal belts are ubiquitous around nearby stars, and their spatial properties hold crucial information for planetesimal and planet formation models. Aims. We present resolved dust observations of 74 planetary systems as part of the REsolved ALMA and SMA Observations of Nearby Stars (REASONS) survey and archival reanalysis. Methods. We uniformly modelled interferometric visibilities for the entire sample to obtain the basic spatial properties of each belt, and combined these with constraints from multi-wavelength photometry. Results. We report key findings from a first exploration of this legacy dataset: (1) Belt dust masses are depleted over time in a radially dependent way, with dust being depleted faster in smaller belts, as predicted by collisional evolution. (2) Most belts are broad discs rather than narrow rings, with much broader fractional widths than rings in protoplanetary discs. We link broad belts to either unresolved substructure or broad planetesimal discs produced if protoplanetary rings migrate. (3) The vertical aspect ratios (h = H/R) of 24 belts indicate orbital inclinations of ~1–20º, implying relative particle velocities of ~0.1–4 km/s, and no clear evolution of heights with system age. This could be explained by early stirring within the belt by large bodies (with sizes of at least ~140 km to the size of the Moon), by inheritance of inclinations from the protoplanetary disc stage, or by a diversity in evolutionary pathways and gravitational stirring mechanisms. We release the REASONS legacy multidimensional sample of millimetre-resolved belts to the community as a valuable tool for follow-up multi-wavelength observations and population modelling studies.

Kinetic field theory applied to planetesimal formation I: Freely streaming dust particles

Abstract Planet formation in the solar system was started when the first planetesimals were formed from the gravitational collapse of pebble clouds. Numerical simulations of this process, especially in the framework of streaming instability, produce various power laws for the initial mass function for planetesimals. While recent advances have shed light on turbulence and its role in particle clustering, a comprehensive theoretical framework linking turbulence characteristics to particle cluster properties and planetesimal mass function remains incomplete. Recently, a kinetic field theory for ensembles of point-like classical particles in or out of equilibrium has been applied to cosmic structure formation. This theory encodes the dynamics of a classical particle ensemble by a generating functional specified by the initial probability distribution of particles in phase space and their equations of motion. Here, we apply kinetic field theory to planetesimal formation. A model for the initial probability distribution of dust particles in phase space is obtained from a quasi-initial state for a three-dimensional streaming-instability simulation that is a particle distribution with velocities for gas and particles from the Nakagawa relations. The equations of motion are chosen for the simplest case of freely streaming particles. We calculate the non-linearly evolved density power spectrum of dust particles and find that it develops a universal k−3 tail at small scales, suggesting scale-invariant structure formation below a characteristic and time-dependent length scale. Thus, the KFT analysis indicates that the initial state for streaming instability simulations does not impose a constraint on structure evolution during planetesimal formation.

Streaming Instabilities in Accreting Protoplanetary Disks: A Parameter Study

The streaming instability (SI) is currently the leading candidate for triggering planetesimal formation in protoplanetary disks. Recently, a novel variation, the “azimuthal-drift” streaming instability (AdSI), was discovered in disks exhibiting laminar gas accretion. Unlike the classical SI, the AdSI does not require pressure gradients and can concentrate dust even at low abundances. We extend previous simulations of the AdSI to explore the impact of dust abundance, accretion-flow strength, pressure gradients, and grain size. For a dimensionless accretion-flow strength α M = 0.1 and particle Stokes number St = 0.1, we find the AdSI produces dust filaments for initial dust-to-gas ratios as low as ϵ = 0.01. For ϵ ≳ 1, maximum dust-to-gas ratios of order 100 are attained, which can be expected to undergo gravitational collapse. Furthermore, even in systems dominated by the classical SI, an accretion flow drives filament formation, without which the disk remains in a state of small-scale turbulence. Our results suggest that an underlying accretion flow facilitates dust concentration and may thus promote planetesimal formation.

Particle fragmentation inside planet-induced spiral waves

ABSTRACT Growing planets interact with their surrounding protoplanetary disc, generating feedback effects that may promote or suppress nearby planet formation. We study how spiral waves launched by planets affect the motion and collisional evolution of particles in the disc. To this end, we perform local 2D hydrodynamical simulations that include a gap-opening planet and integrate particle trajectories within the gas field. Our results show that particle trajectories bend at the location of the spiral wave, and collisions occurring within the spiral exhibit significantly enhanced collisional velocities compared to elsewhere. To quantify this effect, we ran simulations with varying planetary masses and particle sizes. The resulting collisional velocities within the spiral far exceed the typical fragmentation threshold, even for collisions between particles of relatively similar sizes and for planetary masses below the pebble isolation mass. If collisions within the spiral are frequent, this effect could lead to progressively smaller particle sizes as the radial distance from the planet decreases, impacting processes such as gap filtering, pebble accretion, and planetesimal formation.

Planetesimal gravitational collapse in a gaseous environment: Thermal and dynamic evolution

Planetesimal formation models often invoke the gravitational collapse of pebble clouds to overcome various barriers to grain growth and propose processes to concentrate particles sufficiently to trigger this collapse. On the other hand, the geochemical approach for planet formation constrains the conditions for planetesimal formation and evolution by providing temperatures that should be reached to explain the final composition of planetesimals, the building blocks of planets. To elucidate the thermal evolution during gravitational collapse, we used numerical simulations of a self-gravitating cloud of particles and gas coupled with gas drag. Our goal is to determine how the gravitational energy relaxed during the contraction is distributed among the different energy components of the system, and how this constrains a thermal and dynamical planetesimal's history. We identify the conditions necessary to achieve a temperature increase of several hundred kelvins, and as much as 1600 K. Our results emphasise the key role of the gas during the collapse.

High stability of charged particle clusters in protoplanetary disks

Context. The initial particle growth in protoplanetary disks is limited by a bouncing barrier at submillimeter wavelengths. Bouncing leads to tribocharging and the electrostatic attraction of tribocharged aggregates may eventually draw them into large clusters. A charge- mediated growth phase allows for the formation of larger entities, namely, clusters of aggregates that are more prone to further particle concentrations, such as the streaming instability. Aims. We aim to quantify the strength of the electrostatic forces. Methods. In laboratory experiments, we used an acoustic trap to levitate small aggregates of tribocharged submm grains. These aggregates spin up within the trap until they lose grains. Thus, we used the centrifugal force as a measure of the local force. Results. Grains are regularly bound strongly to their neighbors. In comparison, the force at ejection can be stronger than the attractive scattering forces of the trap and can therefore be several orders of magnitude larger than expected. We note that these forces are long- ranging, compared to van der Waals forces. Thus, charged aggregates are much more stable than uncharged ones. Conclusions. Particle aggregates in disks might grow to centimeter clusters or larger as tribocharging increases the effective binding forces. This allows for hydrodynamic concentration and planetesimal formation to eventually take place throughout a wide part of the disk.

Effects of Dust Coagulation on Streaming Instability

Streaming instability (SI) in dust has long been thought to be a promising process in triggering planetesimal formation in the protoplanetary disks (PPDs). In this study, we present the first numerical investigation that models the SI in the vertically stratified disk together with the dust coagulation process. Our simulations reveal that, even with the initially small dust sizes, because dust coagulation promotes dust size growth, SI can eventually still be triggered. Specifically, dust coagulation, limited only by dust fragmentation, broadens the parameter boundaries obtained from previous SI studies using single dust species. We describe the various stages of dust dynamics along with their size evolution and explore the impact of different dust fragmentation velocities. Implications of these results for realistic PPDs are also discussed.

High-resolution ALMA Observations of Richly Structured Protoplanetary Disks in σ Orionis

The Atacama Large Millimeter/submillimeter Array (ALMA) has detected substructures in numerous protoplanetary disks at radii from a few to over 100 au. These substructures are commonly thought to be associated with planet formation, either by serving as sites fostering planetesimal formation or by arising as a consequence of planet–disk interactions. Our current understanding of substructures, though, is primarily based on observations of nearby star-forming regions with mild UV environments, whereas stars are typically born in much harsher UV environments, which may inhibit planet formation in the outer disk through external photoevaporation. We present high-resolution (∼8 au) ALMA 1.3 mm continuum images of eight disks in σ Orionis, a cluster irradiated by an O9.5 star. Gaps and rings are resolved in the images of five disks. The most striking of these is SO 1274, which features five gaps that appear to be arranged nearly in a resonant chain. In addition, we infer the presence of gap or shoulder-like structures in the other three disks through visibility modeling. These observations indicate that substructures robustly form and survive at semimajor axes of several tens of au or less in disks exposed to intermediate levels of external UV radiation as well as in compact disks. However, our observations also suggest that disks in σ Orionis are mostly small, and thus millimeter continuum gaps beyond a disk radius of 50 au are rare in this region, possibly due to either external photoevaporation or age effects.

Semi-quantitative μLIBS mapping of germanium diffusion between metal and silicate during planetary core–mantle segregation

Understanding the distribution of elements between the cores and mantles of planetary bodies is essential for elucidating the geological and geochemical processes operating during their formation. Because the semi-metallic element Ge shows complex siderophile, chalcophile, and lithophile properties, it is particularly significant for investigating core–mantle segregation and ore deposit formation. Recent advancements in in situ microanalysis techniques, such as μLIBS, have enabled high-resolution elemental mapping of Ge and other trace elements.This study explores the application of μLIBS imaging to investigate the high-temperature/high-pressure distribution of Ge between metallic and silicate matrices analogous to planetary cores and mantles. By cross-calibrating the μLIBS intensities with electron microprobe profile analyses carried out on piston cylinder experimental samples, we produced semi-quantitative Ge maps from which we extracted Ge concentration profiles to assess Ge diffusion from the silicate to the metal. We determined, for the first time, the diffusion coefficient of Ge between metal and silicate to be D(Ge)metal = 4.03E−13 ± 3.6E−14 m2/s at 1 GPa and 1350 °C. Our results demonstrate the capability and efficiency of μLIBS for providing detailed, high-resolution (15 μm) trace element concentration data at few ppm levels, offering a time- and cost-effective alternative to conventional techniques.These findings bring insights on planetesimal differentiation and on chemical exchanges between metal and silicate phases in a magma ocean and between two different metal phases occurring at the bottom of the magma ocean during planetesimal formation.

Polar alignment of a dusty circumbinary disc–II. Application to 99 Herculis

ABSTRACT We investigate the formation of dust traffic jams in polar-aligning circumbinary discs. In our first paper, we found as the circumbinary disc evolves towards a polar configuration perpendicular to the binary orbital plane, the differential precession between the gas and dust components leads to multiple dust traffic jams. These dust traffic jams evolve to form a coherent dust ring. In part two, we use 3D smoothed particle hydrodynamical simulations of gas and dust to model an initially highly misaligned circumbinary disc around the 99 Herculis (99 Her) binary system. Our results reveal that the formation of these dust rings is observed across various disc parameters, including the disc aspect ratio, viscosity, surface density power-law index, and temperature power-law index. The dust traffic jams are long-lived and persist even when the disc is fully aligned polar. The midplane dust-to-gas ratio within the rings can surpass unity, which may be a favourable environment for planetesimal formation. Using 2D inviscid shearing box calculations with parameters from our 3D simulations, we find streaming instability modes with significant growth rates. The streaming instability growth time-scale is less than the tilt oscillation time-scale during the alignment process. Therefore, the dust ring will survive once the gas disc aligns polar, suggesting that the streaming instability may aid in forming polar planets around 99 Her.

The Bulk Densities of Small Solar System Bodies as a Probe of Planetesimal Formation

Constraining the formation processes of small solar system bodies is crucial for gaining insights into planetesimal formation. Their bulk densities, determined by their compressive strengths, offer valuable information about their formation history. In this paper, we utilize a formulation of the compressive strength of dust aggregates obtained from dust N-body simulations to establish the relation between the bulk density and diameter. We find that this relation can be effectively approximated by a polytrope with an index of 0.5, coupled with a formulation of the compressive strength of dust aggregates. The lowest-density trans-Neptunian objects (TNOs) and main-belt asteroids (MBAs) are well reproduced by dust aggregates composed of 0.1 μm sized grains. However, most TNOs, MBAs, comets, and near-Earth asteroids (NEAs) exhibit higher densities, suggesting the influence of compaction mechanisms such as collision, dust grain disruption, sintering, or melting, leading to further growth. We speculate that there are two potential formation paths for small solar system bodies. One involves the direct coagulation of primordial dust grains, resulting in the formation of first-generation planetesimals, including the lowest-density TNOs, MBAs, and the parent bodies of comets and NEAs. In this case, comets and NEAs are fragments or rubble piles of first-generation planetesimals, and the objects themselves or the rubble are composed of 0.1 μm sized grains. The other path involves the further potential fragmentation of first-generation planetesimals into the compact dust aggregates observed in protoplanetary disks, resulting in the formation of second-generation planetesimals composed of compact dust aggregates, which may contribute to explaining another formation process of comets and NEAs.

The Ni isotopic composition of Ryugu reveals a common accretion region for carbonaceous chondrites.

The isotopic compositions of samples returned from Cb-type asteroid Ryugu and Ivuna-type (CI) chondrites are distinct from other carbonaceous chondrites, which has led to the suggestion that Ryugu/CI chondrites formed in a different region of the accretion disk, possibly around the orbits of Uranus and Neptune. We show that, like for Fe, Ryugu and CI chondrites also have indistinguishable Ni isotope anomalies, which differ from those of other carbonaceous chondrites. We propose that this unique Fe and Ni isotopic composition reflects different accretion efficiencies of small FeNi metal grains among the carbonaceous chondrite parent bodies. The CI chondrites incorporated these grains more efficiently, possibly because they formed at the end of the disk's lifetime, when planetesimal formation was also triggered by photoevaporation of the disk. Isotopic variations among carbonaceous chondrites may thus reflect fractionation of distinct dust components from a common reservoir, implying CI chondrites/Ryugu may have formed in the same region of the accretion disk as other carbonaceous chondrites.

Dust Growth and Evolution in Protoplanetary Disks

Over the past decade, advancement of observational capabilities, specifically the Atacama Large Millimeter/submillimeter Array (ALMA) and Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instruments, alongside theoretical innovations like pebble accretion, have reshaped our understanding of planet formation and the physics of protoplanetary disks. Despite this progress, mysteries persist along the winded path of micrometer-sized dust, from the interstellar medium, through transport and growth in the protoplanetary disk, to becoming gravitationally bound bodies. This review outlines our current knowledge of dust evolution in circumstellar disks, yielding the following insights: ▪Theoretical and laboratory studies have accurately predicted the growth of dust particles to sizes that are susceptible to accumulation through transport processes like radial drift and settling.▪Critical uncertainties in that process remain the level of turbulence, the threshold collision velocities at which dust growth stalls, and the evolution of dust porosity.▪Symmetric and asymmetric substructures are widespread. Dust traps appear to be solving several long-standing issues in planet formation models, and they are observationally consistent with being sites of active planetesimal formation.▪In some instances, planets have been identified as the causes behind substructures. This underlines the need to study earlier stages of disks to understand how planets can form so rapidly. In the future, better probes of the physical conditions in optically thick regions, including densities, turbulence strength, kinematics, and particle properties, will be essential for unraveling the physical processes at play.

53Mn-53Cr chronometry of ureilites: Implications for the timing of parent body accretion, differentiation and secondary reduction

Establishing the temporal evolution of the ureilite parent body(ies) is crucial for understanding the quantitative timescale of planetesimal formation and evolution in the protoplanetary disk. In order to establish a timeline for these early processes, age constraints on the accretion, differentiation and secondary reduction were obtained with the short-lived 53Mn-53Cr chronometer to whole-rock and sequentially digested fractions of main group ureilites. A whole-rock isochron dates the reservoir-scale Mn-Cr fractionation in the ureilite parent body(ies), associated with magmatic differentiation, to 2.89-0.51+0.56 Ma after CAI formation. This age implies that the ureilite parent body(ies) accreted no later than ∼1.5 Ma after CAI formation, at a time when the NC-CC dichotomy was already established. The 53Mn-53Cr systematics of fractions from chromite-bearing ureilites yield an age of 4.29-0.45+0.49 Ma after CAI formation for a secondary reduction event on the parent body. This event is commonly associated with the catastrophic disruption of the ureilite parent body while still hot. The chromite model ages are consistent with the isochron ages obtained from chromite-bearing ureilites. Collectively these ages indicate that chemical differentiation processes were underway on the ureilite parent body(ies) during the time interval when undifferentiated meteorite parent bodies were forming, and may have paused at the peak of planetesimal formation when planetary collisions were common.

Sequential giant planet formation initiated by disc substructure

Context. Planet formation models are necessary to understand the origins of diverse planetary systems. Circumstellar disc substructures have been proposed as preferred locations of planet formation, but a complete formation scenario has not been covered by a single model so far. Aims. We aim to study the formation of giant planets facilitated by disc substructure and starting with sub-micron-sized dust. Methods. We connect dust coagulation and drift, planetesimal formation, N-body gravity, pebble accretion, planet migration, planetary gas accretion, and gap opening in one consistent modelling framework. Results. We find rapid formation of multiple gas giants from the initial disc substructure. The migration trap near the substructure allows for the formation of cold gas giants. A new pressure maximum is created at the outer edge of the planetary gap, which triggers the next generation of planet formation resulting in a compact chain of giant planets. A high planet formation efficiency is achieved, as the first gas giants are effective at preventing dust from drifting further inwards, which preserves material for planet formation. Conclusions. Sequential planet formation is a promising framework to explain the formation of chains of gas and ice giants.

Properties of outer solar system pebbles during planetesimal formation from meteor observations

Observations of proto-planetary disks, as well as theoretical modeling, suggest that in the late stages of accretion leading up to the formation of planetesimals, particles grew to pebbles the size of 1-mm to tens of cm, depending on the location and ambient conditions in the disk. That is the same size range that dominates the present-day comet and primitive asteroid mass loss. Meteoroids that size cause visible meteors on Earth. Here, we hypothesize that the size distribution and the physical and chemical properties of young meteoroid streams still contain information about the conditions in the solar nebula during these late stages of accretion towards planetesimal formation. If so, they constrain where long-period (Oort Cloud) comets, Jupiter-family (Scattered Disk Kuiper Belt) comets, and primitive asteroids (Asteroid Belt) formed. From video and visual observations of 47 young meteor showers, we find that freshly ejected meteoroids from long-period comets tend to have low bulk density and are distributed with equal surface area per log-mass interval (magnitude distribution index χ ∼ 1.85), suggesting gentle accretion conditions. Jupiter-family comets, on the other hand, mostly produce meteoroids twice as dense and distributed with a steeper χ ∼ 2.15 or even χ ∼ 2.5, which implies that those pebbles grew from particles fragmenting in a collisional cascade or by catastrophic collisions, respectively. Some primitive asteroids show χ > 2.5, with most mass in small particles, indicating an even more aggressive fragmentation by processes other than mutual collisions. Both comet populations contain an admixture of compact materials that are sometimes sodium-poor, but Jupiter-family comets show a higher percentage (∼8% on average) than long-period comet showers (∼4%) and a wider range of percentages among comets. While there are exceptions in both groups, the implication is that most long-period comets formed under gentle particle growth conditions, possibly near the 30 AU edge of the Trans Neptunian Disk, while most Jupiter family comets formed closer to the Sun where pebbles reached or passed the fragmentation barrier, and primitive asteroids formed in the region where the cores of the giant planets formed. This is possible if the Scattered Disk represents all objects scattered by Neptune during its migration, while the present-day outer Oort cloud formed only during and after the time of the planet instability, well after the Sun had moved away from sibling stars.

Hydrogen sulfide and metal-enriched atmosphere for a Jupiter-mass exoplanet.

As the closest transiting hot Jupiter to Earth, HD 189733b has been the benchmark planet for atmospheric characterization1-3. It has also been the anchor point for much of our theoretical understanding of exoplanet atmospheres from composition4, chemistry5,6, aerosols7 to atmospheric dynamics8, escape9 and modelling techniques10,11. Previous studies of HD 189733b have detected carbon and oxygen-bearing molecules H2O and CO (refs. 12,13) in the atmosphere. The presence of CO2 and CH4 has been claimed14,15 but later disputed12,16,17. The inferred metallicity based on these measurements, a key parameter in tracing planet formation locations18, varies from depletion19,20 to enhancement21,22, hindered by limited wavelength coverage and precision of the observations. Here we report detections of H2O (13.4σ), CO2 (11.2σ), CO (5σ) and H2S (4.5σ) in the transmission spectrum (2.4-5.0 μm) of HD 189733b. With an equilibrium temperature of about 1,200 K, H2O, CO and H2S are the main reservoirs for oxygen, carbon and sulfur. Based on the measured abundances of these three main volatile elements, we infer an atmospheric metallicity of three to five times stellar. The upper limit on the methane abundance at 5σ is 0.1 ppm, which indicates a low carbon-to-oxygen ratio (<0.2), suggesting formation through the accretion of water-rich icy planetesimals. The low oxygen-to-sulfur and carbon-to-sulfur ratios also support the planetesimal accretion formation pathway23.

Streaming Instability and Turbulence: Conditions for Planetesimal Formation

The streaming instability (SI) is a leading candidate for planetesimal formation, which can concentrate solids through two-way aerodynamic interactions with the gas. The resulting concentrations can become sufficiently dense to collapse under particle self-gravity, forming planetesimals. Previous studies have carried out large parameter surveys to establish the critical particle to gas surface density ratio (Z), above which SI-induced concentration triggers planetesimal formation. The threshold Z depends on the dimensionless stopping time (τ s , a proxy for dust size). However, these studies neglected both particle self-gravity and external turbulence. Here, we perform 3D stratified shearing box simulations with both particle self-gravity and turbulent forcing, which we characterize via a turbulent diffusion parameter, α D. We find that forced turbulence, at amplitudes plausibly present in some protoplanetary disks, can increase the threshold Z by up to an order of magnitude. For example, for τ s = 0.01, planetesimal formation occurs when Z ≳ 0.06, ≳0.1, and ≳0.2 at α D = 10−4, 10−3.5, and 10−3, respectively. We provide a single fit to the critical Z required for the SI to work as a function of α D and τ s (although limited to the range τ s = 0.01–0.1). Our simulations also show that planetesimal formation requires a mid-plane particle-to-gas density ratio that exceeds unity, with the critical value being largely insensitive to α D. Finally, we provide an estimation of particle scale height that accounts for both particle feedback and external turbulence.