This white paper proposes that AMBITION, a Comet Nucleus Sample Return mission, be a cornerstone of ESA's Voyage 2050 programme. We summarise some of the most important questions still open in cometary science after the successes of the Rosetta mission, many of which require sample analysis using techniques that are only possible in laboratories on Earth. We then summarise measurements, instrumentation and mission scenarios that can address these questions, with a recommendation that ESA select an ambitious cryogenic sample return mission. Rendezvous missions to Main Belt comets and Centaurs are compelling cases for M-class missions, expanding our knowledge by exploring new classes of comets. AMBITION would engage a wide community, drawing expertise from a vast range of disciplines within planetary science and astrophysics. With AMBITION, Europe will continue its leadership in the exploration of the most primitive Solar System bodies.

On 27 Apr 2015, when 67P/C-G was at 1.76 au from the Sun and moving towards perihelion, the OSIRIS and VIRTIS-M instruments on Rosetta observed the evolving dust and gas coma during a complete rotation of the comet. We aim to characterize the dust, H2O and CO2 gas spatial distribution in the inner coma. To do this we performed a quantitative analysis of the release of dust and gas and compared the observed H2O production rate with the one calculated using a thermo-physical model. For this study we selected OSIRIS WAC images at 612 nm (dust) and VIRTIS-M image cubes at 612 nm, 2700 nm (H2O) and 4200 nm (CO2). We measured the average signal in a circular annulus, to study spatial variation around the comet, and in a sector of the annulus, to study temporal variation in the sunward direction with comet rotation, both at a fixed distance of 3.1 km from the comet centre. The spatial correlation between dust and water, both coming from the sun-lit side of the comet, shows that water is the main driver of dust activity in this time period. The spatial distribution of CO2 is not correlated with water and dust. There is no strong temporal correlation between the dust brightness and water production rate as the comet rotates. The dust brightness shows a peak at 0deg sub-solar longitude, which is not pronounced in the water production. At the same epoch, there is also a maximum in CO2 production. An excess of measured water production, with respect to the value calculated using a simple thermo-physical model, is observed when the head lobe and regions of the Southern hemisphere with strong seasonal variations are illuminated. A drastic decrease in dust production, when the water production (both measured and from the model) displays a maximum, happens when typical Northern consolidated regions are illuminated and the Southern hemisphere regions with strong seasonal variations are instead in shadow.

The southern hemisphere of comet 67P/Churyumov-Gerasimenko (67P) became observable by the Rosetta mission in March 2015, a few months before cometary southern vernal equinox. The Anhur region in the southern part of the comet's larger lobe was found to be highly eroded, enriched in volatiles, and highly active. We analyze high-resolution images of the Anhur region pre- and post-perihelion acquired by the OSIRIS imaging system on board the Rosetta mission. The Narrow Angle Camera is particularly useful for studying the evolution in Anhur in terms of morphological changes and color variations.}{Radiance factor images processed by the OSIRIS pipeline were coregistered, reprojected onto the 3D shape model of the comet, and corrected for the illumination conditions. We find a number of morphological changes in the Anhur region that are related to formation of new scarps; removal of dust coatings; localized resurfacing in some areas, including boulders displacements; and vanishing structures, which implies localized mass loss that we estimate to be higher than 50 million kg. The strongest changes took place in and nearby the Anhur canyon-like structure, where significant dust cover was removed, an entire structure vanished, and many boulders were rearranged. All such changes are potentially associated with one of the most intense outbursts registered by Rosetta during its observations, which occurred one day before perihelion passage. Moreover, in the niche at the foot of a new observed scarp, we also see evidence of water ice exposure that persisted for at least six months. The abundance of water ice, evaluated from a linear mixing model, is relatively high (> 20%). Our results confirm that the Anhur region is volatile-rich and probably is the area on 67P with the most pristine exposures near perihelion.

Before Rosetta, the space missions Giotto and Stardust shaped our view on cometary dust, supported by plentiful data from Earth based observations and interplanetary dust particles collected in the Earth's atmosphere. The Rosetta mission at comet 67P/Churyumov-Gerasimenko was equipped with a multitude of instruments designed to study cometary dust. While an abundant amount of data was presented in several individual papers, many focused on a dedicated measurement or topic. Different instruments, methods, and data sources provide different measurement parameters and potentially introduce different biases. This can be an advantage if the complementary aspect of such a complex data set can be exploited. However, it also poses a challenge in the comparison of results in the first place. The aim of this work therefore is to summarise dust results from Rosetta and before. We establish a simple classification as a common framework for inter-comparison. This classification is based on a dust particle's structure, porosity, and strength as well as its size. Depending on the instrumentation, these are not direct measurement parameters but we chose them as they were the most reliable to derive our model. The proposed classification already proved helpful in the Rosetta dust community and we propose to take it into consideration also beyond. In this manner we hope to better identify synergies between different instruments and methods in the future.

In April 2016, the Rosetta spacecraft performed a low-altitude low-phase-angle flyby over the Imhotep-Khepry transition of 67P/Churyumov-Gerasimenko's nucleus. The OSIRIS/Narrow-Angle-Camera (NAC) acquired 112 images with mainly 3 broadband filters in the visible at a resolution of up to 0.53 m/px and for phase angles between 0.095{\deg} and 62{\deg}. Using those images, we have investigated the morphological and spectrophotometrical properties of this area. We assembled the images into coregistered color cubes. Using a 3D shape model, we produced the illumination conditions and georeference for each image. We projected the observations on a map to investigate its geomorphology. Observations were photometrically corrected using the Lommel-Seeliger disk law. Spectrophotometric analyses were performed on the coregistered color cubes. These data were used to estimate the local phase reddening. This region of the nucleus hosts numerous and varied types of terrains and features. We observe an association between a feature's nature, its reflectance, and its spectral slope. Fine material deposits exhibit an average reflectance and spectral slope, while terrains with diamictons, consolidated material, degraded outcrops, or features such as somber boulders, present a lower-than-average reflectance and higher-than-average spectral slope. Bright surfaces present here a spectral behavior consistent with terrains enriched in water-ice. We find a phase-reddening slope of 0.064{\pm}0.001{\%}/100nm/{\deg} at 2.7 au outbound, similarly to the one obtained at 2.3 au inbound during the February 2015 flyby. Identified as the source region of multiple jets and a host of water-ice material, the Imhotep-Khepry transition appeared in April 2016, close to the frost line, to further harbor several potential locations with exposed water-ice material among its numerous different morphological terrain units.

We directly measure twenty overhanging cliffs on the surface of comet 67P/Churyumov-Gerasimenko extracted from the latest shape model and estimate the minimum tensile strengths needed to support them against collapse under the comet's gravity. We find extremely low strengths of around one Pa or less (one to five Pa, when scaled to a metre length). The presence of eroded material at the base of most overhangs, as well as the observed collapse of two features and implied previous collapse of another, suggests that they are prone to failure and that true material strengths are close to these lower limits (although we only consider static stresses and not dynamic stress from, for example, cometary activity). Thus, a tensile strength of a few pascals is a good approximation for the tensile strength of 67P's nucleus material, which is in agreement with previous work. We find no particular trends in overhang properties with size, over the $\sim10-100$ m range studied here, or location on the nucleus. There are no obvious differences, in terms of strength, height or evidence of collapse, between the populations of overhangs on the two cometary lobes, suggesting that 67P is relatively homogenous in terms of tensile strength. Low material strengths are supportive of cometary formation as a primordial rubble pile or by collisional fragmentation of a small (tens of km) body.

Porosity is a key characteristic of the regolith on the surface of small bodies. The porosity of the regolith on the surface of asteroids is changed by applied pressure, and the relationship between pressure and porosity depends on the particle properties of the regolith. We performed compression measurements on samples of different materials, particle size distributions, and shapes to examine the relationship between particle size and pressure required for compression. We used a centrifuge and a compression testing machine for the experiments. The applied pressure for the centrifuge and the compression testing machine experiments ranged from 10^2 to 5x10^3 Pa and from 10^4 to 5x10^6 Pa, respectively. The initial porosity before compression was generally higher for samples with smaller particles and narrower particle size distributions. A sample compressed more easily when it consisted of smaller particles, probably due to smaller frictional forces between particles. We estimated the porosity of granular asteroids based on the results of our experiments. The estimated porosity at the center of a homogeneous asteroid with a radius of 50 km is 0.43. This porosity is typical of the porosity of asteroids of similar size.

We present a statistical analysis of the distribution of large scale topographic features on comet 67P/Churyumov-Gerasimenko. We observe that the cumulative cliff height distribution across the surface follows a power law with a slope equal to -1.69 +- 0.02. When this distribution is studied independently for each region, we find a good correlation between the slope of the power law and the orbital erosion rate of the surface. For instance, the northern hemisphere topography is dominated by structures on the 100~m scale while the southern hemisphere topography, illuminated at perihelion, is dominated by 10~m scale terrain features. Our study suggest that the current size of a cliff is controlled not only by material cohesion but by the dominant erosional process in each region. This observation can be generalized to other comets, where we argue that primitive nuclei are characterized by the presence of large cliffs with a cumulative height power index equal to or above -1.5, while older, eroded cometary surfaces have a power index equal to or below -2.3. In effect, our model shows that a measure of the topography provides a quantitative assessment of a comet's erosional history, i.e. its evolutionary age.

We collect observational evidence that supports the scheme of mass transfer on the nucleus of comet 67P/Churyumov-Gerasimenko. The obliquity of the rotation axis of 67P causes strong seasonal variations. During perihelion the southern hemisphere is four times more active than the north. Northern territories are widely covered by granular material that indicates back fall originating from the active south. Decimetre sized chunks contain water ice and their trajectories are influenced by an anti-solar force instigated by sublimation. OSIRIS observations suggest that up to 20 % of the particles directly return to the nucleus surface taking several hours of travel time. The back fall covered northern areas are active if illuminated but produce mainly water vapour. The decimetre chunks from the nucleus surface are too small to contain more volatile compounds such as CO 2 or CO. This causes a north-south dichotomy of the composition measurements in the coma. Active particles are trapped in the gravitational minimum of Hapi during northern winter. They are "shock frozen" and only reactivated when the comet approaches the sun after its aphelion passage. The insolation of the big cavity is enhanced by self-heating, i. e. reflection and IR radiation from the walls. This, together with the pristinity of the active back fall, explains the early observed activity of the Hapi region. Sobek may be a role model for the consolidated bottom of Hapi. Mass transfer in the case of 67P strongly influences the evolution of the nucleus and the interpretation of coma measurements.

The ESA Rosetta spacecraft has been tracking its target, the Jupiter-family comet 67P/Churyumov-Gerasimenko, in close vicinity for over two years. It hosts the OSIRIS instruments: the Optical, Spectroscopic, and Infrared Remote Imaging System composed of two cameras, see e.g. Keller et al. (2007). In some imaging sequences dedicated to observe dust particles in the comet's coma, the two cameras took images at the same time. The aim of this work is to use these simultaneous double camera observations to calculate the dust particles' distance to the spacecraft. As the two cameras are mounted on the spacecraft with an offset of 70 cm, the distance of particles observed by both cameras can be determined by a shift of the particles' apparent trails on the images. This paper presents first results of the ongoing work, introducing the distance determination method for the OSIRIS instrument and the analysis of an example particle. We note that this method works for particles in the range of about 500 m - 6000 m from the spacecraft.

We present the results of the photometric and spectrophotometric properties of the 67P/Churyumov-Gerasimenko nucleus derived with the OSIRIS instrument during the closest fly-by over the comet, which took place on 14 th February 2015 at a distance of {\~} 6 km from the surface. Several images covering the 0{\deg}-33{\deg} phase angle range were acquired, and the spatial resolution achieved was 11 cm/pxl. The flown-by region is located on the big lobe of the comet, near the borders of the Ash, Apis and Imhotep regions. Our analysis shows that this region features local heterogeneities at the decimetre scale. We observed difference of reflectance up to 40{\%} between bright spots and sombre regions, and spectral slope variations up to 50{\%}. The spectral reddening effect observed globally on the comet surface by Fornasier et al. (2015) is also observed locally on this region, but with a less steep behaviour. We note that numerous metre-sized boulders, which exhibit a smaller opposition effect, also appear spectrally redder than their surroundings. In this region, we found no evidence linking observed bright spots to exposed water-ice-rich material. We fitted our dataset using the Hapke 2008 photometric model. The region overflown is globally as dark as the whole nucleus (geometric albedo of 6.8{\%}) and it has a high porosity value in the uppermost-layers (86{\%}). These results of the photometric analysis at a decimetre scale indicate that the photometric properties of the flown-by region are similar to those previously found for the whole nucleus.

Dust jets, i.e. fuzzy collimated streams of cometary material arising from the nucleus, have been observed in-situ on all comets since the Giotto mission flew by comet 1P/Halley in 1986. Yet their formation mechanism remains unknown. Several solutions have been proposed, from localized physical mechanisms on the surface/sub-surface (see review in Belton (2010)) to purely dynamical processes involving the focusing of gas flows by the local topography (Crifo et al. 2002). While the latter seems to be responsible for the larger features, high resolution imagery has shown that broad streams are composed of many smaller features (a few meters wide) that connect directly to the nucleus surface. We monitored these jets at high resolution and over several months to understand what are the physical processes driving their formation, and how this affects the surface. Using many images of the same areas with different viewing angles, we performed a 3-dimensional reconstruction of collimated jets, and linked them precisely to their sources on the nucleus. Results.We show here observational evidence that the Northern hemisphere jets of comet 67P arise from areas with sharp topographic changes and describe the physical processes involved. We propose a model in which active cliffs are the main source of jet-like features, and therefore the regions eroding the fastest on comets. We suggest that this is a common mechanism taking place on all comets.

We study the link between gravitational slopes and the surface morphology on the nucleus of comet 67P/Churyumov-Gerasimenko and provide constraints on the mechanical properties of the cometary material. We computed the gravitational slopes for five regions on the nucleus that are representative of the different morphologies observed on the surface, using two shape models computed from OSIRIS images by the stereo-photoclinometry (SPC) and stereo-photogrammetry (SPG) techniques. We estimated the tensile, shear, and compressive strengths using different surface morphologies and mechanical considerations. The different regions show a similar general pattern in terms of the relation between gravitational slopes and terrain morphology: i) low-slope terrains (0-20 deg) are covered by a fine material and contain a few large ($>$10 m) and isolated boulders, ii) intermediate-slope terrains (20-45 deg) are mainly fallen consolidated materials and debris fields, with numerous intermediate-size boulders from $<$1 m to 10 m for the majority of them, and iii) high-slope terrains (45-90 deg) are cliffs that expose a consolidated material and do not show boulders or fine materials. The best range for the tensile strength of overhangs is 3-15 Pa (upper limit of 150 Pa), 4-30 Pa for the shear strength of fine surface materials and boulders, and 30-150 Pa for the compressive strength of overhangs (upper limit of 1500 Pa). The strength-to-gravity ratio is similar for 67P and weak rocks on Earth. As a result of the low compressive strength, the interior of the nucleus may have been compressed sufficiently to initiate diagenesis, which could have contributed to the formation of layers. Our value for the tensile strength is comparable to that of dust aggregates formed by gravitational instability and tends to favor a formation of comets by the accrection of pebbles at low velocities.

We report on the first major temporal morphological changes observed on the surface of the nucleus of comet 67P/Churyumov-Gerasimenko, in the smooth terrains of the Imhotep region. We use images of the OSIRIS cameras onboard Rosetta to follow the temporal changes from 24 May 2015 to 11 July 2015. The morphological changes observed on the surface are visible in the form of roundish features, which are growing in size from a given location in a preferential direction, at a rate of 5.6 - 8.1$\times$10$^{-5}$ m s$^{-1}$ during the observational period. The location where changes started and the contours of the expanding features are bluer than the surroundings, suggesting the presence of ices (H$_2$O and/or CO$_2$) exposed on the surface. However, sublimation of ices alone is not sufficient to explain the observed expanding features. No significant variations in the dust activity pattern are observed during the period of changes.

One of the main aims of the ESA Rosetta mission is to study the origin of the solar system by exploring comet 67P/Churyumov-Gerasimenko at close range. In this paper we discuss the origin and evolution of comet 67P/Churyumov-Gerasimenko in relation to that of comets in general and in the framework of current solar system formation models. We use data from the OSIRIS scientific cameras as basic constraints. In particular, we discuss the overall bi-lobate shape and the presence of key geological features, such as layers and fractures. We also treat the problem of collisional evolution of comet nuclei by a particle-in-a-box calculation for an estimate of the probability of survival for 67P/Churyumov-Gerasimenko during the early epochs of the solar system. We argue that the two lobes of the 67P/Churyumov-Gerasimenko nucleus are derived from two distinct objects that have formed a contact binary via a gentle merger. The lobes are separate bodies, though sufficiently similar to have formed in the same environment. An estimate of the collisional rate in the primordial, trans-planetary disk shows that most comets of similar size to 67P/Churyumov-Gerasimenko are likely collisional fragments, although survival of primordial planetesimals cannot be excluded. A collisional origin of the contact binary is suggested, and the low bulk density of the aggregate and abundance of volatile species show that a very gentle merger must have occurred. We thus consider two main scenarios: the primordial accretion of planetesimals, and the re-accretion of fragments after an energetic impact onto a larger parent body. We point to the primordial signatures exhibited by 67P/Churyumov-Gerasimenko and other comet nuclei as critical tests of the collisional evolution.

The Rosetta mission of the European Space Agency has been orbiting the comet 67P/Churyumov-Gerasimenko (67P) since August 2014 and is now in its escort phase. A large complement of scientific experiments designed to complete the most detailed study of a comet ever attempted are onboard Rosetta. We present results for the photometric and spectrophotometric properties of the nucleus of 67P derived from the OSIRIS imaging system, which consists of a Wide Angle Camera (WAC) and a Narrow Angle Camera (NAC). The disk-averaged phase function of the nucleus of 67P shows a strong opposition surge with a G parameter value of -0.13$\pm$0.01 in the HG system formalism and an absolute magnitude $H_v(1,1,0)$ = 15.74$\pm$0.02 mag. The integrated spectrophotometry in 20 filters covering the 250-1000 nm wavelength range shows a red spectral behavior, without clear absorption bands except for a potential absorption centered at $\sim$ 290 nm that is possibly due to SO$_2$ ice. The nucleus shows strong phase reddening, with disk-averaged spectral slopes increasing from 11\%/(100 nm) to 16\%/(100 nm) in the 1.3$^{\circ}$--54$^{\circ}$ phase angle range. The geometric albedo of the comet is 6.5$\pm$0.2\% at 649 nm, with local variations of up to $\sim$ 16\% in the Hapi region. From the disk-resolved images we computed the spectral slope together with local spectrophotometry and identified three distinct groups of regions (blue, moderately red, and red). The Hapi region is the brightest, the bluest in term of spectral slope, and the most active surface on the comet. Local spectrophotometry shows an enhancement of the flux in the 700-750 nm that is associated with coma emissions.

All atmosphere-less planetary bodies are covered with a dust layer, the so-called regolith, which determines the optical, mechanical and thermal properties of their surface. These properties depend on the regolith material, the size distribution of the particles it consists of, and the porosity to which these particles are packed. We performed experiments in parabolic flights to determine the gravity dependency of the packing density of regolith for solid-particle sizes of 60 $\mu$m and 1 mm as well as for 100-250 $\mu$m-sized agglomerates of 1.5 $\mu$m-sized solid grains. We utilized g-levels between 0.7 m s$^{-2}$ and 18 m s$^{-2}$ and completed our measurements with experiments under normal gravity conditions. Based on previous experimental and theoretical literature and supported by our new experiments, we developed an analytical model to calculate the regolith stratification of celestial rocky and icy bodies and estimated the mechanical yields of the regolith under the weight of an astronaut and a spacecraft resting on these objects.

We present laboratory experiments of a vertically vibrated granular medium consisting of 1 mm diameter glass beads with embedded 8 mm diameter intruder glass beads. The experiments were performed in the laboratory as well as in a parabolic flight under reduced-gravity conditions (on Martian and Lunar gravity levels). We measured the mean rise velocity of the large glass beads and present its dependence on the fill height of the sample containers, the excitation acceleration, and the ambient gravity level. We find that the rise velocity scales in the same manner for all three gravity regimes and roughly linearly with gravity.

We conducted micro-gravity experiments to study the outcome of collisions between sub-mm-sized dust agglomerates consisting of \mu m-sized SiO2 monomer grains at velocities of several cm/s. Prior to the experiments, we used X-ray computer tomography (nano-CT) imaging to study the internal structure of these dust agglomerates and found no rim compaction so that their collision behavior is not governed by preparation-caused artefacts. We found that collisions between these dust aggregates can lead either to sticking or to bouncing, depending mostly on the impact velocity. While previous collision models derived the transition between both regimes from contact physics, we used the available empirical data from these and earlier experiments to derive a power law relation between dust-aggregate mass and impact velocity for the threshold between the two collision outcomes. In agreement with earlier experiments, we show that the transition between both regimes is not sharp, but follows a shallower power law than predicted by previous models (G\"uttler et al. 2010). Furthermore, we find that sticking between dust aggregates can lead to the formation of larger structures. Collisions between aggregates-of-aggregates can lead to growth at higher velocities than homogeneous dust agglomerates.

Recently published crater statistics on the small asteroids 25143 Itokawa and 433 Eros show a significant depletion of craters below approx. 100 m in diameter. Possible mechanisms that were brought up to explain this lack of craters were seismic crater erasure and self armoring of a coarse, boulder covered asteroid surface. While seismic shaking has been studied in this context, the concept of armoring lacks a deeper inspection and an experimental ground truth. We therefore present cratering experiments of glass bead projectiles impacting into granular glass bead targets, where the grain sizes of projectile and target are in a similar range. The impact velocities are in the range of 200 to 300 m/s. We find that craters become fainter and irregular shaped as soon as the target grains are larger than the projectile sizes and that granular craters rarely form when the size ratio between projectile and target grain is around 1:10 or smaller. In that case, we observe a formation of a strength determined crater in the first struck target grain instead. We present a simple model based on the transfer of momentum from the projectile to this first target grain, which is capable to explain our results with only a single free parameter, which is moreover well determined by previous experiments. Based on estimates of typical projectile size and boulder size on Itokawa and Eros, given that our results are representative also for km/s impact velocities, armoring should play an important role for their evolution.

The formation of planetesimals is often accredited to collisional sticking of dust grains. The exact process is unknown, as collisions between larger aggregates tend to lead to fragmentation or bouncing rather than sticking. Recent laboratory experiments have however made great progress in the understanding and mapping of the complex physics involved in dust collisions. We want to study the possibility of planetesimal formation using the results from the latest laboratory experiments, particularly by including the fragmentation with mass transfer effect, which might lead to growth even at high impact velocities. We present a new experimentally and physically motivated dust collision model capable of predicting the outcome of a collision between two particles of arbitrary masses and velocities. It is used together with a continuum dust-size evolution code that is both fast in terms of execution time and able to resolve the dust well at all sizes, allowing for all types of interactions to be studied without biases. We find that for the general dust population, bouncing collisions prevent the growth above millimeter-sizes. However, if a small number of cm-sized particles are introduced, for example due to vertical mixing or radial drift, they can act as a catalyst and start to sweep up the smaller particles. At a distance of 3 AU, 100-meter-sized bodies are formed on a timescale of 1 Myr. We conclude that direct growth of planetesimals might be a possibility thanks to a combination of the existence of a bouncing barrier and the fragmentation with mass transfer effect. The bouncing barrier is here even beneficial, as it prevents the growth of too many large particles that would otherwise only fragment among each other, and creates a reservoir of small particles that can be swept up by larger bodies. However, for this process to work, a few seeds of cm in size or larger have to be introduced.

The formation of planetesimals in the early Solar System is hardly understood, and in particular the growth of dust aggregates above millimeter sizes has recently turned out to be a difficult task in our understanding [Zsom et al. 2010, A&A, 513, A57]. Laboratory experiments have shown that dust aggregates of these sizes stick to one another only at unreasonably low velocities. However, in the protoplanetary disk, millimeter-sized particles are known to have been ubiquitous. One can find relics of them in the form of solid chondrules as the main constituent of chondrites. Most of these chondrules were found to feature a fine-grained rim, which is hypothesized to have formed from accreting dust grains in the solar nebula. To study the influence of these dust-coated chondrules on the formation of chondrites and possibly planetesimals, we conducted collision experiments between millimeter-sized, dust-coated chondrule analogs at velocities of a few cm/s. For 2 and 3 mm diameter chondrule analogs covered by dusty rims of a volume filling factor of 0.18 and 0.35-0.58, we found sticking velocities of a few cm/s. This velocity is higher than the sticking velocity of dust aggregates of the same size. We therefore conclude that chondrules may be an important step towards a deeper understanding of the collisional growth of larger bodies. Moreover, we analyzed the collision behavior in an ensemble of dust aggregates and non-coated chondrule analogs. While neither the dust aggregates nor the solid chondrule analogs show sticking in collisions among their species, we found an enhanced sicking efficiency in collisions between the two constituents, which leads us to the conjecture that chondrules might act as "catalyzers" for the growth of larger bodies in the young Solar System.

Over the past years the processes involved in the growth of planetesimals have extensively been studied in the laboratory. Based on these experiments, a dust-aggregate collision model was developed upon which computer simulations were based to evaluate how big protoplanetary dust aggregates can grow and to analyze which kinds of collisions are relevant in the solar nebula and are worth further studies in the laboratory. The sticking threshold velocity of millimeter-sized dust aggregates is such a critical value that had so far only theoretically been derived, as the relevant velocities could not be reached in the laboratory. We developed a microgravity experiment that allows us for the first time to study free collisions of mm-sized dust aggregates down to velocities of ~0.1 cm/s to assess this part of the protoplanetary dust evolution model. Here, we present the results of 125 free collisions between dust aggregates of 0.5 to 2 mm diameter. Seven collisions with velocities between 0.2 and 3 cm/s led to sticking, suggesting a transition from perfect sticking to perfect bouncing with a certain sticking probability instead of a sharp velocity threshold. We developed a model to explain the physical processes involved in dust-aggregate sticking, derived dynamic material properties of the dust aggregates from the results of the collisions, and deduced the velocity below which aggregates always stick. For millimeter-sized porous dust aggregates this velocity is 8e-5 m/s.

Collisions between centimeter- to decimeter-sized dusty bodies are important to understand the mechanisms leading to the formation of planetesimals. We thus performed laboratory experiments to study the collisional behavior of dust aggregates in this size range at velocities below and around the fragmentation threshold. We developed two independent experimental setups with the same goal to study the effects of bouncing, fragmentation, and mass transfer in free particle-particle collisions. The first setup is an evacuated drop tower with a free-fall height of 1.5 m, providing us with 0.56 s of microgravity time so that we observed collisions with velocities between 8 mm/s and 2 m/s. The second setup is designed to study the effect of partial fragmentation (when only one of the two aggregates is destroyed) and mass transfer in more detail. It allows for the measurement of the accretion efficiency as the samples are safely recovered after the encounter. Our results are that for very low velocities we found bouncing as could be expected while the fragmentation velocity of 20 cm/s was significantly lower than expected. We present the critical energy for disruptive collisions Q*, which showed up to be at least two orders of magnitude lower than previous experiments in the literature. In the wide range between bouncing and disruptive collisions, only one of the samples fragmented in the encounter while the other gained mass. The accretion efficiency in the order of a few percent of the particle's mass is depending on the impact velocity and the sample porosity. Our results will have consequences for dust evolution models in protoplanetary disks as well as for the strength of large, porous planetesimal bodies.

In the last years, experiments have shown that collisions above the fragmentation threshold velocity are a potentially important growth process for protoplanatary dust aggregates. To obtain deeper understanding of this process, we performed laboratory and drop-tower experiments to study multiple impacts of small, porous dust-aggregate projectiles onto sintered dust targets. Projectile and target consisted of 1.5 micron monodisperse, spherical SiO2 monomers with volume filling factors of 0.15 (projectiles) and 0.45 (targets). The projectiles were accelerated by a solenoid magnet and combined with a magazine with which 25 impacts onto the same spot on the target could be performed in vacuum. We measured the mass-accretion efficiency and the volume filling factor for different impact velocities between 1.5 and 6.0 m/s. The experiments at the lowest impact speeds were performed in the Bremen drop-tower under microgravity conditions. Within this velocity range we found a linear increase of the accretion efficiency with increasing velocity. In the laboratory experiments, the accretion efficiency increases from 0.12 to 0.21 in units of the projectile mass. The recorded images of the impacts showed that the mass transfer from the projectile to the target leads to the growth of a conical structure on the target. From the images we also measured the volume filling factors of the grown structures, which ranged from 0.15 (uncompacted) to 0.40 (significantly compacted) with increasing impact speed. These results augment our knowledge of the aggregate growth in protoplanetary disks and should be taken into account for future models of protoplanetary dust growth.

A highly favoured mechanism of planetesimal formation is collisional growth. Single dust grains, which follow gas flows in the protoplanetary disc, hit each other, stick due to van der Waals forces and form fluffy aggregates up to centimetre size. The mechanism of further growth is unclear since the outcome of aggregate collisions in the relevant velocity and size regime cannot be investigated in the laboratory under protoplanetary disc conditions. Realistic statistics of the result of dust aggregate collisions beyond decimetre size is missing for a deeper understanding of planetary growth. Joining experimental and numerical efforts we want to calibrate and validate a computer program that is capable of a correct simulation of the macroscopic behaviour of highly porous dust aggregates. After testing its numerical limitations thoroughly we will check the program especially for a realistic reproduction of various benchmark experiments. We adopt the smooth particle hydrodynamics (SPH) numerical scheme with extensions for the simulation of solid bodies and a modified version of the Sirono porosity model. Experimentally measured macroscopic material properties of silica dust are implemented. We calibrate and test for the compressive strength relation and the bulk modulus. SPH has already proven to be a suitable tool to simulate collisions at rather high velocities. In this work we demonstrate that its area of application can not only be extended to low-velocity experiments and collisions. It can also be used to simulate the behaviour of highly porous objects in this velocity regime to a very high accuracy.The result of the calibration process in this work is an SPH code that can be utilised to investigate the collisional outcome of porous dust in the low-velocity regime.

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.

Recent years have shown many advances in our knowledge of the collisional evolution of protoplanetary dust. Based on a variety of dust-collision experiments in the laboratory, our view of the growth of dust aggregates in protoplanetary disks is now supported by a deeper understanding of the physics involved in the interaction between dust agglomerates. However, the parameter space, which determines the collisional outcome, is huge and sometimes inaccessible to laboratory experiments. Very large or fluffy dust aggregates and extremely low collision velocities are beyond the boundary of today's laboratories. It is therefore desirable to augment our empirical knowledge of dust-collision physics with a numerical method to treat arbitrary aggregate sizes, porosities and collision velocities. In this article, we implement experimentally-determined material parameters of highly porous dust aggregates into a Smooth Particle Hydrodynamics (SPH) code, in particular an omnidirectional compressive-strength and a tensile-strength relation. We also give a prescription of calibrating the SPH code with compression and low-velocity impact experiments. In the process of calibration, we developed a dynamic compressive-strength relation and estimated a relation for the shear strength. Finally, we defined and performed a series of benchmark tests and found the agreement between experimental results and numerical simulations to be very satisfactory. SPH codes have been used in the past to study collisions at rather high velocities. At the end of this work, we show examples of future applications in the low-velocity regime of collisional evolution.

In the recent years we have performed various experiments on the collision dynamics of highly porous dust aggregates and although we now have a comprehensive picture of the micromechanics of those aggregates, the macroscopic understanding is still lacking. We are therefore developing a mechanical model to describe dust aggregate collisions with macroscopic parameters like tensile strength, compressive strength and shear strength. For one well defined dust sample material, the tensile and compressive strength were measured in a static experiment and implemented in a Smoothed Particle Hydrodynamics (SPH) code. A laboratory experiment was designed to compare the laboratory results with the results of the SPH simulation. In this experiment, a mm-sized glass bead is dropped into a cm-sized dust aggregate with the previously measured strength parameters. We determine the deceleration of the glass bead by high-speed imaging and the compression of the dust aggregate by x-ray micro-tomography. The measured penetration depth, stopping time and compaction under the glass bead are utilized to calibrate and test the SPH code. We find that the statically measured compressive strength curve is only applicable if we adjust it to the dynamic situation with a 'softness' parameter. After determining this parameter, the SPH code is capable of reproducing experimental results, which have not been used for the calibration before.

To study the evolution of protoplanetary dust aggregates, we performed experiments with up to 2600 collisions between single, highly-porous dust aggregates and a solid plate. The dust aggregates consisted of spherical SiO$_2$ grains with 1.5$\mu$m diameter and had an initial volume filling factor (the volume fraction of material) of $\phi_0=0.15$. The aggregates were put onto a vibrating baseplate and, thus, performed multiple collisions with the plate at a mean velocity of 0.2 m s$^{-1}$. The dust aggregates were observed by a high-speed camera to measure their size which apparently decreased over time as a measure for their compaction. After 1000 collisions the volume filling factor was increased by a factor of two, while after $\sim2000$ collisions it converged to an equilibrium of $\phi\approx0.36$. In few experiments the aggregate fragmented, although the collision velocity was well below the canonical fragmentation threshold of $\sim1$ m s$^{-1}$. The compaction of the aggregate has an influence on the surface-to-mass ratio and thereby the dynamic behavior and relative velocities of dust aggregates in the protoplanetary nebula. Moreover, macroscopic material parameters, namely the tensile strength, shear strength, and compressive strength, are altered by the compaction of the aggregates, which has an influence on their further collisional behavior. The occurrence of fragmentation requires a reassessment of the fragmentation threshold velocity.

In order to investigate the hypothesis that dust aggregates were transformed to meteoritic chondrules by nebular lightning, we exposed silicatic and metallic dust samples to electric discharges with energies of 120 to 500 J in air at pressures between 10 and 10^5 Pa. The target charges consisted of powders of micrometer-sized particles and had dimensions of mm. The dust samples generally fragmented leaving the major fraction thermally unprocessed. A minor part formed sintered aggregates of 50 to 500 micrometer. In a few experiments melt spherules having sizes smaller than 180 micrometer in diameter (and, generally, interior voids) were formed; the highest spherule fraction was obtained with metallic Ni. Our experiments indicate that chondrule formation by electric current or by particle bombardment inside a discharge channel is unlikely.