The outcome of protoplanetary dust growth: pebbles, boulders, or planetesimals? I. Mapping the zoo of laboratory collision experiments

The growth processes from protoplanetary dust to planetesimals are not fully understood. Laboratory experiments and theoretical models have shown that collisions among the dust aggregates can lead to sticking, bouncing, and fragmentation. However, no systematic study on the collisional outcome of protoplanetary dust has been performed so far so that a physical model of the dust evolution in protoplanetary disks is still missing. We intend to map the parameter space for the collisional interaction of arbitrarily porous dust aggregates. This parameter space encompasses the dust-aggregate masses, their porosities and the collision velocity. With such a complete mapping of the collisional outcomes of protoplanetary dust aggregates, it will be possible to follow the collisional evolution of dust in a protoplanetary disk environment. We use literature data, perform own laboratory experiments, and apply simple physical models to get a complete picture of the collisional interaction of protoplanetary dust aggregates. In our study, we found four different types of sticking, two types of bouncing, and three types of fragmentation as possible outcomes in collisions among protoplanetary dust aggregates. We distinguish between eight combinations of porosity and mass ratio. For each of these cases, we present a complete collision model for dust-aggregate masses between 10^-12 and 10^2 g and collision velocities in the range 10^-4 to 10^4 cm/s for arbitrary porosities. This model comprises the collisional outcome, the mass(es) of the resulting aggregate(s) and their porosities. We present the first complete collision model for protoplanetary dust. This collision model can be used for the determination of the dust-growth rate in protoplanetary disks.Comment: accepted by Astronomy and Astrophysic

Collisions of inhomogeneous pre-planetesimals

In the framework of the coagulation scenario, kilometre-sized planetesimals form by subsequent collisions of pre-planetesimals of sizes from centimetre to hundreds of metres. Pre-planetesimals are fluffy, porous dust aggregates, which are inhomogeneous owing to their collisional history. Planetesimal growth can be prevented by catastrophic disruption in pre-planetesimal collisions above the destruction velocity threshold. We develop an inhomogeneity model based on the density distribution of dust aggregates, which is assumed to be a Gaussian distribution with a well-defined standard deviation. As a second input parameter, we consider the typical size of an inhomogeneous clump. These input parameters are easily accessible by laboratory experiments. For the simulation of the dust aggregates, we utilise a smoothed particle hydrodynamics (SPH) code with extensions for modelling porous solid bodies. The porosity model was previously calibrated for the simulation of silica dust, which commonly serves as an analogue for pre-planetesimal material. The inhomogeneity is imposed as an initial condition on the SPH particle distribution. We carry out collisions of centimetre-sized dust aggregates of intermediate porosity. We vary the standard deviation of the inhomogeneous distribution at fixed typical clump size. The collision outcome is categorised according to the four-population model. We show that inhomogeneous pre-planetesimals are more prone to destruction than homogeneous aggregates. Even slight inhomogeneities can lower the threshold for catastrophic disruption. For a fixed collision velocity, the sizes of the fragments decrease with increasing inhomogeneity. Pre-planetesimals with an active collisional history tend to be weaker. This is a possible obstacle to collisional growth and needs to be taken into account in future studies of the coagulation scenario.Comment: 12 pages, 9 figures, 4 table

The four-populations model: a new classification scheme for pre-planetesimal collisions

Within the collision growth scenario for planetesimal formation, the growth step from centimetre sized pre-planetesimals to kilometre sized planetesimals is still unclear. The formation of larger objects from the highly porous pre-planetesimals may be halted by a combination of fragmentation in disruptive collisions and mutual rebound with compaction. However, the right amount of fragmentation is necessary to explain the observed dust features in late T Tauri discs. Therefore, detailed data on the outcome of pre-planetesimal collisions is required and has to be presented in a suitable and precise format. We propose and apply a new classification scheme for pre-planetesimal collisions based on the quantitative aspects of four fragment populations: the largest and second largest fragment, a power-law population, and a sub-resolution population. For the simulations of pre-planetesimal collisions, we adopt the SPH numerical scheme with extensions for the simulation of porous solid bodies. By means of laboratory benchmark experiments, this model was previously calibrated and tested for the correct simulation of the compaction, bouncing, and fragmentation behaviour of macroscopic highly porous silica dust aggregates. It is shown that previous attempts to map collision data were much too oriented on qualitatively categorising into sticking, bouncing, and fragmentation events. We show that the four-populations model encompasses all previous categorisations and in addition allows for transitions. This is because it is based on quantitative characteristic attributes of each population such as the mass, kinetic energy, and filling factor. As a demonstration of the applicability and the power of the four-populations model, we utilise it to present the results of a study on the influence of collision velocity in head-on collisions of intermediate porosity aggregates.Comment: 14 pages, 11 figures, 5 tables, to be published in Astronomy and Astrophysic

Compression Behaviour of Porous Dust Agglomerates

The early planetesimal growth proceeds through a sequence of sticking collisions of dust agglomerates. Very uncertain is still the relative velocity regime in which growth rather than destruction can take place. The outcome of a collision depends on the bulk properties of the porous dust agglomerates. Continuum models of dust agglomerates require a set of material parameters that are often difficult to obtain from laboratory experiments. Here, we aim at determining those parameters from ab-initio molecular dynamics simulations. Our goal is to improveon the existing model that describe the interaction of individual monomers. We use a molecular dynamics approach featuring a detailed micro-physical model of the interaction of spherical grains. The model includes normal forces, rolling, twisting and sliding between the dust grains. We present a new treatment of wall-particle interaction that allows us to perform customized simulations that directly correspond to laboratory experiments. We find that the existing interaction model by Dominik & Tielens leads to a too soft compressive strength behavior for uni and omni-directional compression. Upon making the rolling and sliding coefficients stiffer we find excellent agreement in both cases. Additionally, we find that the compressive strength curve depends on the velocity with which the sample is compressed. The modified interaction strengths between two individual dust grains will lead to a different behaviour of the whole dust agglomerate. This will influences the sticking probabilities and hence the growth of planetesimals. The new parameter set might possibly lead to an enhanced sticking as more energy can be stored in the system before breakup.Comment: 11 pages, 14 figures, accepted for publication in A&

Breaking through: The effects of a velocity distribution on barriers to dust growth

It is unknown how far dust growth can proceed by coagulation. Obstacles to collisional growth are the fragmentation and bouncing barriers. However, in all previous simulations of the dust-size evolution in protoplanetary disks, only the mean collision velocity has been considered, neglecting that a small but possibly important fraction of the collisions will occur at both much lower and higher velocities. We study the effect of the probability distribution of impact velocities on the collisional dust growth barriers. Assuming a Maxwellian velocity distribution for colliding particles to determine the fraction of sticking, bouncing, and fragmentation, we implement this in a dust-size evolution code. We also calculate the probability of growing through the barriers and the growth timescale in these regimes. We find that the collisional growth barriers are not as sharp as previously thought. With the existence of low-velocity collisions, a small fraction of the particles manage to grow to masses orders of magnitude above the main population. A particle velocity distribution softens the fragmentation barrier and removes the bouncing barrier. It broadens the size distribution in a natural way, allowing the largest particles to become the first seeds that initiate sweep-up growth towards planetesimal sizes.Comment: 4 pages, 3 figures. Accepted for publication as a Letter in Astronomy and Astrophysic

Dust size distributions in coagulation/fragmentation equilibrium: Numerical solutions and analytical fits

Context. Grains in circumstellar disks are believed to grow by mutual collisions and subsequent sticking due to surface forces. Results of many fields of research involving circumstellar disks, such as radiative transfer calculations, disk chemistry, magneto-hydrodynamic simulations largely depend on the unknown grain size distribution. Aims. As detailed calculations of grain growth and fragmentation are both numerically challenging and computationally expensive, we aim to find simple recipes and analytical solutions for the grain size distribution in circumstellar disks for a scenario in which grain growth is limited by fragmentation and radial drift can be neglected. Methods. We generalize previous analytical work on self-similar steady-state grain distributions. Numerical simulations are carried out to identify under which conditions the grain size distributions can be understood in terms of a combination of power-law distributions. A physically motivated fitting formula for grain size distributions is derived using our analytical predictions and numerical simulations. Results. We find good agreement between analytical results and numerical solutions of the Smoluchowski equation for simple shapes of the kernel function. The results for more complicated and realistic cases can be fitted with a physically motivated "black box" recipe presented in this paper. Our results show that the shape of the dust distribution is mostly dominated by the gas surface density (not the dust-to-gas ratio), the turbulence strength and the temperature and does not obey an MRN type distribution.Comment: 16 pages, 9 figures, accepted for publication in A&

Crossing barriers in planetesimal formation: The growth of mm-dust aggregates with large constituent grains

Collisions of mm-size dust aggregates play a crucial role in the early phases of planet formation. We developed a laboratory setup to observe collisions of dust aggregates levitating at mbar pressures and elevated temperatures of 800 K. We report on collisions between basalt dust aggregates of from 0.3 to 5 mm in size at velocities between 0.1 and 15 cm/s. Individual grains are smaller than 25 \mum in size. We find that for all impact energies in the studied range sticking occurs at a probability of 32.1 \pm 2.5% on average. In general, the sticking probability decreases with increasing impact parameter. The sticking probability increases with energy density (impact energy per contact area). We also observe collisions of aggregates that were formed by a previous sticking of two larger aggregates. Partners of these aggregates can be detached by a second collision with a probability of on average 19.8 \pm 4.0%. The measured accretion efficiencies are remarkably high compared to other experimental results. We attribute this to the rel. large dust grains used in our experiments, which make aggregates more susceptible to restructuring and energy dissipation. Collisional hardening by compaction might not occur as the aggregates are already very compact with only 54 \pm 1% porosity. The disassembly of previously grown aggregates in collisions might stall further aggregate growth. However, owing to the levitation technique and the limited data statistics, no conclusive statement about this aspect can yet be given. We find that the detachment efficiency decreases with increasing velocities and accretion dominates in the higher velocity range. For high accretion efficiencies, our experiments suggest that continued growth in the mm-range with larger constituent grains would be a viable way to produce larger aggregates, which might in turn form the seeds to proceed to growing planetesimals.Comment: 9 pages, 20 figure

Thermal history modeling of the H chondrite parent body

The cooling histories of individual meteorites can be empirically reconstructed by using ages from different radioisotopic chronometers with distinct closure temperatures. For a group of meteorites derived from a single parent body such data permit the reconstruction of the cooling history and properties of that body. Particularly suited are H chondrites because precise radiometric ages over a wide range of closure temperatures are available. A thermal evolution model for the H chondrite parent body is constructed by using all H chondrites for which at least three different radiometric ages are available. Several key parameters determining the thermal evolution of the H chondrite parent body and the unknown burial depths of the H chondrites are varied until an optimal fit is obtained. The fit is performed by an 'evolution algorithm'. Empirical data for eight samples are used for which radiometric ages are available for at least three different closure temperatures. A set of parameters for the H chondrite parent body is found that yields excellent agreement (within error bounds) between the thermal evolution model and empirical data of six of the examined eight chondrites. The new thermal model constrains the radius and formation time of the H chondrite parent body (possibly (6) Hebe), the initial burial depths of the individual H chondrites, the average surface temperature of the body, the average initial porosity of the material the body accreted from, and the initial 60Fe content of the H chondrite parent body.Comment: 16 pages, 7 figure

The outcome of protoplanetary dust growth: pebbles, boulders, or planetesimals? II. Introducing the bouncing barrier

The sticking of micron sized dust particles due to surface forces in circumstellar disks is the first stage in the production of asteroids and planets. The key ingredients that drive this process are the relative velocity between the dust particles in this environment and the complex physics of dust aggregate collisions. Here we present the results of a collision model, which is based on laboratory experiments of these aggregates. We investigate the maximum aggregate size and mass that can be reached by coagulation in protoplanetary disks. We model the growth of dust aggregates at 1 AU at the midplane at three different gas densities. We find that the evolution of the dust does not follow the previously assumed growth-fragmentation cycles. Catastrophic fragmentation hardly occurs in the three disk models. Furthermore we see long lived, quasi-steady states in the distribution function of the aggregates due to bouncing. We explore how the mass and the porosity change upon varying the turbulence parameter and by varying the critical mass ratio of dust particles. Particles reach Stokes numbers of roughly 10^-4 during the simulations. The particle growth is stopped by bouncing rather than fragmentation in these models. The final Stokes number of the aggregates is rather insensitive to the variations of the gas density and the strength of turbulence. The maximum mass of the particles is limited to approximately 1 gram (chondrule-sized particles). Planetesimal formation can proceed via the turbulent concentration of these aerodynamically size-sorted chondrule-sized particles.Comment: accepted for publication in A&

Superior antigen cross-presentation and XCR1 expression define human CD11c+CD141+ cells as homologues of mouse CD8+ dendritic cells

In recent years, human dendritic cells (DCs) could be subdivided into CD304+ plasmacytoid DCs (pDCs) and conventional DCs (cDCs), the latter encompassing the CD1c+, CD16+, and CD141+ DC subsets. To date, the low frequency of these DCs in human blood has essentially prevented functional studies defining their specific contribution to antigen presentation. We have established a protocol for an effective isolation of pDC and cDC subsets to high purity. Using this approach, we show that CD141+ DCs are the only cells in human blood that express the chemokine receptor XCR1 and respond to the specific ligand XCL1 by Ca2+ mobilization and potent chemotaxis. More importantly, we demonstrate that CD141+ DCs excel in cross-presentation of soluble or cell-associated antigen to CD8+ T cells when directly compared with CD1c+ DCs, CD16+ DCs, and pDCs from the same donors. Both in their functional XCR1 expression and their effective processing and presentation of exogenous antigen in the context of major histocompatibility complex class I, human CD141+ DCs correspond to mouse CD8+ DCs, a subset known for superior antigen cross-presentation in vivo. These data define CD141+ DCs as professional antigen cross-presenting DCs in the human