Issue |
A&A
Volume 593, September 2016
|
|
---|---|---|
Article Number | A3 | |
Number of page(s) | 21 | |
Section | Planets and planetary systems | |
DOI | https://doi.org/10.1051/0004-6361/201527288 | |
Published online | 29 August 2016 |
Submillimetre-sized dust aggregate collision and growth properties
Experimental study of a multi-particle system on a suborbital rocket
1 Institut für Geophysik und
extraterrestrische Physik, Technische Universität
Braunschweig, Mendelssohnstr.
3, 38106 Braunschweig, Germany
e-mail: j.brisset@tu-bs.de
2 Max Planck Institute for Solar System
Research, Justus-von-Liebig-Weg
3, 37077
Göttingen,
Germany
Received:
1
September
2015
Accepted:
25
May
2016
Context. In the very first steps of the formation of a new planetary system, dust agglomerates grow inside the protoplanetary disk that rotates around the newly formed star. In this disk, collisions between the dust particles, induced by interactions with the surrounding gas, lead to sticking. Aggregates start growing until their sizes and relative velocities are high enough for collisions to result in bouncing or fragmentation. With the aim of investigating the transitions between sticking and bouncing regimes for colliding dust aggregates and the formation of clusters from multiple aggregates, the Suborbital Particle and Aggregation Experiment (SPACE) was flown on the REXUS 12 suborbital rocket.
Aims. The collisional and sticking properties of sub-mm-sized aggregates composed of protoplanetary dust analogue material are measured, including the statistical threshold velocity between sticking and bouncing, their surface energy and tensile strength within aggregate clusters.
Methods. We performed an experiment on the REXUS 12 suborbital rocket. The protoplanetary dust analogue materials were micrometre-sized monodisperse and polydisperse SiO2 particles prepared into aggregates with sizes around 120 μm and 330 μm, respectively and volume filling factors around 0.37. During the experimental run of 150 s under reduced gravity conditions, the sticking of aggregates and the formation and fragmentation of clusters of up to a few millimetres in size was observed.
Results. The sticking probability of the sub-mm-sized dust aggregates could be derived for velocities decreasing from ~22 to 3 cm s-1. The transition from bouncing to sticking collisions happened at 12.7+2.1-1.4 cm s-1 for the smaller aggregates composed of monodisperse particles and at 11.5+1.9-1.3 and 11.7+1.9-1.3 cm s-1 for the larger aggregates composed of mono- and polydisperse dust particles, respectively. Using the pull-off force of sub-mm-sized dust aggregates from the clusters, the surface energy of the aggregates composed of monodisperse dust was derived to be 1.6 × 10-5 J m-2, which can be scaled down to 1.7 × 10-2 J m-2 for the micrometre-sized monomer particles and is in good agreement with previous measurements for silica particles. The tensile strengths of these aggregates within the clusters were derived to be 1.9+2.2-1.2 Pa and 1.6+0.7-0.6 Pa for the small and large dust aggregates, respectively. These values are in good agreement with recent tensile strength measurements for ~mm-sized silica aggregates.
Conclusions. Using our data on the sticking-bouncing threshold, estimates of the maximum aggregate size can be given. For a minimum mass solar nebula model, aggregates can reach sizes of ~1 cm.
Key words: protoplanetary disks / accretion, accretion disks / planets and satellites: formation
© ESO, 2016
1. Introduction
The formation of planets around a young star takes place in the remains of accreting clouds of gas and dust, which form disks around their pre-main-sequence stars (Weidenschilling & Cuzzi 1993; Weidenschilling 2000; Dominik et al. 2007). In these so-called protoplanetary disks (PPDs), the dust has condensed into solid grains that range in size from sub-μm- to μm (Bouwman et al. 2008). Several stochastic (Brownian motion) and aerodynamic processes (e.g. gas drag and turbulence) influence the dynamics of these grains in the nebula and induce velocity differences between them, which lead to collisions and initiate the growth of aggregates (e.g. see the reviews by Johansen et al. 2014; Testi et al. 2014; Blum & Wurm 2008).
For grain sizes of about one micrometre in these PPDs, collisions always lead to grain growth through van der Waals forces (Blum & Wurm 2008). This is the so-called “hit-and-stick” behaviour. The growing dust aggregates, however, soon become big enough to leave this regime. As their sizes increase, their relative velocities increase accordingly (Weidenschilling 1977a) and collisions then lead to restructuring, bouncing, or even fragmentation, rendering processes that allow further grain growth more complex (Blum & Wurm 2008; Güttler et al. 2010; Zsom et al. 2010). Reproducing the complete growth from μm-sized dust grains to km-sized planetesimals is an ongoing endeavour. Amongst the many processes to account for in PPDs, e.g. concentration of dust aggregates in streaming instabilities (Youdin & Goodman 2005) or gravitational instabilities of dense dust layers (Johansen et al. 2009, 2012), the collision behaviour of growing dust aggregates also serves as a crucial input to PPD models (Windmark et al. 2012a,b; Garaud et al. 2013). Knowing under which conditions and by which processes dust growth is possible is of utmost importance for understanding and simulating the first stages of the formation of planetary bodies.
While the exact conditions present in PPDs are still under discussion, investigating the collision behaviour of dust aggregates, without any influence of gas, temperature nor other disturbances is already very useful for developing dust evolution models. In order to further our understanding of the processes underlying the dust collision behaviour, many experiments as well as numerical simulations have been (and still are) performed (see reviews by Dominik et al. 2007; Blum & Wurm 2008). Güttler et al. (2010) compiled the results of several laboratory and microgravity dust collision experiments into a model predicting the outcome of a collision between two arbitrary dust aggregates, depending on their mass, porosity, and relative velocity. In particular, the transition between sticking and bouncing collisions is an area of interest because the number of experimental investigations for macroscopic dust aggregates probing into this transitions is very small. Indeed, depending on the collision outcomes, simulations of dust aggregate growth in PPDs can have very different results, even leading to stalling growth at centimetre sizes (Wada et al. 2011; Zsom et al. 2010). As far as we know, only two dust-aggregate collision experiments were conducted with particle sizes and at relative velocities adequate to observe this sticking-bouncing transition. Weidling et al. (2012) observed mm-sized silicate dust aggregates colliding at velocities of ~1 cm s-1 and Kothe et al. (2013) observed 100 μm-sized aggregates at ~10 cm s-1 (see Kothe et al. 2013, for an update of the collision model). These experiments were performed with many-particle systems, which also display “group” effects like clustering (Kothe et al. 2013).
In this context, the experimental investigation presented in this paper was designed to gather additional data on collisions between sub-mm-sized dust aggregates. One of the challenges to observing dust collisions in many-particle systems in a range of velocities down to 10 cm s-1 and below is the necessity to conduct these experiments under reduced gravity conditions. The flight profile of suborbital rockets allows for an experiment to remain in such reduced gravity conditions (down to 10-3g) for several minutes. This paper presents the Suborbital Particle Aggregation and Collision Experiment (SPACE) that flew on the REXUS 12 rocket (Rocket EXperiments for University Students, a joint project between DLR – Deutsches Zentrum für Luft-und Raumfahrt – and SNSB – Swedish National Space Board) in March 2012 and the analysis of the gathered data. In Sect. 2 we present the experimental setup and flight performance as well as the dust aggregates we chose to observe. Section 3 presents the results of the analysis of the gathered data. In Sect. 4 we discuss these results and their applicability to dust collisions in protoplanetary disks.
2. Experimental setup
This section presents the experiment hardware setup and the scientific data gathered by the SPACE experiment, which flew onboard the REXUS 12 suborbital rocket in March 2012. More details about the experimental hardware and rocket can be found in Brisset et al. (2013).
2.1. Experiment hardware
The purpose of the experiment being to observe dust-aggregate collisions in a multi-particle system, we placed sub-mm-sized dust aggregates into an evacuated glass container and recorded their behaviour during the extended microgravity phase with a high-speed camera. The particle container was divided in three cells, two smaller ones of dimensions 11 × 10 × 15 mm3 and a bigger cell with dimensions of 24 × 10 × 15 mm3 (see Fig. 1). These cells allowed for the observation of different types of dust aggregates during a single rocket flight. The inner glass walls of these cells were coated with an anti-adhesive nano-layer to reduce their sticking efficiency with the dust aggregates (see Fig. 5 in Brisset et al. 2013, for more details).
![]() |
Fig. 1 Experiment cells inside the rocket payload. a) Schematic of the REXUS 12 rocket with the position of the SPACE module marked by a red asterisk. b) CAD drawing of the SPACE payload including the vacuum chamber (1) and the camera (2) observing the interior of the chamber through a mirror (3). c) Schematic view of the glass experiment cells (5) installed on the shaking mechanism (6). The gears of the shaking mechanism activated by a motor induced the circular motion of the cells indicated on the figure by blue arrows. The circle indicates the limits of the vacuum chamber viewport. d) Photograph of the SPACE experiment cells inside the vacuum chamber and seen through the view port. The types of dust inserted into each cell were ~120 μm-sized monodisperse in the larger cell (I) and ~330 μm-sized poly- and monodisperse in each of the smaller cells (II and III, resp.). |
In these types of microgravity dust-aggregate collision experiments, experience has shown that in order to obtain and maintain a uniform spatial distribution of the aggregates in the cell volume (e.g. against microgravity disturbances by air drag acting on the rocket), a shaking of some kind has to take place. Indeed, if a many-particle system is left on its own in microgravity, it kinetically “cools” down quite rapidly, losing energy with each collision (Haff 1983) and leading to a lower collision frequency among the aggregates. Shaking also allows for the deagglomeration of clusters formed during the experiment preparation as well as for a control over the mean collision velocity between the aggregates. As suborbital rockets are subject to residual accelerations (see Appendix C.1), shaking of the experiment was an essential part of the hardware setup. The two main microgravity disturbances during flight were the atmospheric drag on the rocket body, which acts in the direction of flight, and the spinning centrifugal force, which acts in the outward radial direction and perpendicular to the direction of flight. Therefore, we decided to shake the dust cells in a circular motion along these two directions (see Fig. 1). The radius of the cell rotation being 1 mm, shaking frequencies between 5 and 30 Hz led to linear wall velocities between about 3 and 22 cm s-1.
The shaking-velocity profile included three cycles, each composed of a fast shaking period of 5 to 10 s and a longer slow shaking period of 15 to 25 s with transition ramps between them (see Fig. 2). This shaking profile was designed to observe three agglomeration phases (during ramping down the rotation speed of the cells) and three de-agglomeration phases (during ramping up the rotation speed of the cells). The slowest and longest shaking phase was scheduled around the apogee of the rocket trajectory, when the quality of the microgravity was at its best (i.e. the residual accelerations were lowest). As shown in Fig. 2, the shaking of the experiment kept the particle systems in a relatively high acceleration environment. However, this profile choice allowed for the observation of the dust aggregates over a continuous range of speeds, both decreasing and increasing. In addition, the slow shaking speeds of cycle 2 kept the particle systems under 1g for a significant amount of time (~30 s) during which the transition from a bouncing to a sticking regime were observed.
![]() |
Fig. 2 Measured rotation frequency of the particle cells during the experimental run of SPACE onboard the REXUS 12 rocket. Time is measured from the first recorded frame (~100 s after lift-off of the rocket). The corresponding wall centrifugal acceleration and maximal linear speed are indicated. The apogee of the rocket trajectory is indicated by a “+” sign. |
To be able to observe the aggregate collisions without the influence of gas, the glass containers were placed in a vacuum chamber, which was evacuated during the experiment. The pressure was logged during the rocket flight and remained below 10-4 mbar (10-2 Pa) for the duration of the experiment.
The evolution of the dust-aggregate ensembles in the three glass cells were simultaneously and continuously recorded with an autonomous high-speed camera with a chip size of 640 × 480 pixels, resulting in a spatial resolution of 56.6 μm/px, at a recording speed of 170 frames per second.
2.2. Dust aggregate properties
In this paper, we will refer to the single dust aggregates (a few 100 μm in size) introduced into the experiment glass cells before the beginning of the experiment as aggregates, dust aggregates, or monomer aggregates. These aggregates themselves consist of smaller dust grains of ~1 μm in size, which we will call monomer particles. In their storage container, these monomer particles form aggregates, which we then sieved to the desired size distributions. When several of these dust aggregates stick together during the experiment to form a bigger agglomerate, we will refer to it as a cluster.
The dust analogue material we used in the SPACE experiment was SiO2, identical to that used by Weidling et al. (2012) and Kothe et al. (2013), and frequently used in dust collision experiments (see Blum & Wurm 2008, and references therein). The material density of SiO2 is ρ = 2000 kg m-3 for the (slightly nano-porous) spherical monomer particles (monodisperse dust) and ρ = 2600 kg m-3 for the irregular monomer particles (polydisperse dust). Further properties can be retrieved from Blum et al. (2006b). Kothe et al. (2013) determined the average volume filling factor of the sieved clusters to be , in agreement with the measurements of Weidling et al. (2012). Building on Güttler et al. (2009) and Zsom et al. (2010), we consider this value of the aggregate filling factor to be representative of protoplanetary dust aggregates at the sizes of aggregates studied here. It should be mentioned that a filling factor of φ = 0.4 as found by Zsom et al. (2010) is representative of aggregates consisting of ~1 μm silicate particles. However, Kataoka et al. (2013) found that 0.1 μm ice particles form aggregates that are fluffier. Kothe et al. (2013) also investigated the inner structure of the aggregates with the result that the standard preparation process of the aggregates does not induce any kind of compacted shell, meaning that a homogeneous structure of the aggregates can be assumed (see Fig. 4 in Kothe et al. 2013).
Before the REXUS flight, the two SiO2 samples (monodisperse, see Fig. 3b., and polydisperse, see Fig. 3a.) were sieved into two different size distributions, one between 100 and 250 μm and the other one between 250 and 500 μm. The mean aggregate size was then measured to be 120 μm and 330 μm, respectively (see Fig. 4). The dust samples were introduced into the experiment cells as follows (see Fig. 1b.): 19.2 mg of the spherical monodisperse SiO2 dust and aggregate sizes between 100 and 250 μm were placed into the bigger glass cell, 15.3 mg of the spherical monodisperse SiO2 dust and aggregate sizes between 250 and 500 μm were placed into the right (seen from the camera location) smaller cell, and 19.8 mg of the irregular polydisperse SiO2 dust and aggregate sizes between 250 and 500 μm were placed into the left (seen from the camera location) smaller cell. These quantities are the total masses of dust inserted into each experiment cell.
![]() |
Fig. 3 SiO2 dust used in the SPACE experiment: a) SEM image of the polydisperse SiO2 monomer particles (image credit: E. Beitz); b) SEM image of monodisperse SiO2 monomer particles (image credit: E. Beitz). c) Microscope image of aggregates composed of irregular polydisperse SiO2, sieved between 250 and 500 μm. |
![]() |
Fig. 4 Size and mass distributions of the aggregate samples used in this study. a) Normalized cumulative number of aggregates with a size smaller than the size indicated on the x-axis. b) Normalized cumulative mass of aggregates smaller than the size indicated on the x-axis for the three different SiO2 aggregates used in the SPACE experiment. Black diamonds denote aggregates sieved between 100 and 250 μm and composed of monodisperse dust particles. Black triangles denote aggregates sieved between 250 and 500 μm and composed of monodisperse dust particles. Green asterisks denote aggregates sieved between 250 and 500 μm and composed of polydisperse dust. Dashed lines denote aggregate sizes measured during the experiment run on the rocket. |
2.3. Scientific data
![]() |
Fig. 5 Two image frames recorded during the SPACE experiment. a) At the very beginning of the experimental run, when all the dust aggregates were free-flying. b) During a slow shaking phase of the experiment (see Fig. 2). In slow shaking phases, the preferential accumulation of aggregates in the cells’ left upper corner due to microgravity disturbances (see C.1 for details) is visible. The flight direction of the rocket (i.e. direction of gas drag) and the outward radial direction (i.e. direction of centrifugal force) are indicated. |
The SPACE experiment performed nominally during the REXUS 12 flight. Two example frames of the gathered scientific data can be seen in Fig. 5. During the experiment run, the observed collisions between dust aggregates formed clusters on the inner walls of the glass containers. We show in Appendix B that this was due to the non-negligible residual acceleration increasing the number of collisions between the dust aggregates and the container walls. In the two directions in which the experiment cells were shaken, this had no major influence on the dust-glass interactions, as the shaking accelerations were almost 3 orders of magnitude stronger than the flight residual accelerations (see Appendix C.1). However, no shaking was applied in the line of sight of the camera. The increased number of dust-glass collisions due to the residual flight accelerations led to the sticking of a small number of aggregates to the cell walls with the very low sticking probability of about 0.5% (the cell walls were coated anti-adhesively, see Appendix B). These aggregates sticking to the walls acted as seeds for further cluster growth (see Appendix A).
During the first and third shaking phases (see Fig. 2), the aggregate speeds were too high (>10 cm s-1) to detect a significant change in collision behaviour. In the second phase, however, a clear transition from the bouncing to the sticking regime could be observed while the shaking speed was being ramped down (<5 cm s-1). During the speed ramp-up of this same phase, the dismantling of the formed clusters could be observed.
3. Results
3.1. Growing clusters on the cell walls
This section describes how cluster formation on the glass walls of the experiment cell can deliver information on dust aggregate properties and growth. To support the fact that clusters grew through aggregate-cluster collisions instead of aggregate-glass collisions we performed a series of Monte Carlo simulations in which we randomly deposited particles on a surface, representing the cell glass walls (see Appendix B). As detailed in Appendix B, we performed numerous runs with varying the number of aggregates as well as their sticking probability to the glass walls. In each case, the sticking probability in an aggregate-cluster collision was set to unity. We then compared the number and morphology of the forming clusters to those found in the experiments.
The result of this investigation was that a good match between simulation and experiment was found when the probability of a dust aggregate to stick to the glass upon collision was as low as 0.5%, which shows that the applied anti-adhesive coating was intact over 99.5% of the surface but showed defects at places where the impinging aggregates could stick. This means that it was not an enhanced glass sticking efficiency that led to clusters on the cell walls but indeed the enhanced number of collisions of the aggregates on the glass, induced by residual accelerations during the rocket flight (see Appendix C.1). It can be concluded that an individual dust cluster was not growing on the glass directly but on a single dust aggregates that stuck to the glass and served as a seed for further impinging aggregates. Accordingly, clusters in the SPACE experiment grew by collisions between dust aggregates and the clusters on the wall so that the further data analysis will deliver information on the collision and sticking properties of these aggregates.
The difference between collisions among free-flying aggregates and those between aggregates and clusters sticking to a glass wall is that the cluster mass in the latter case can be considered as infinite compared to a finite mass of free-flying clusters. Aggregate collisions simulated by Wada et al. (2011) show for example that a low-velocity collision (<2 m s-1 for their icy aggregates) only locally affects the involved aggregates (see their Fig. 5) so that their surface properties are unchanged elsewhere. A collision of a free-flying aggregate with a cluster or an aggregate on the wall is therefore equivalent to a collision with a free-flying cluster or aggregate, except for a slight difference (up to a factor 2 in the extreme case of two equal-size aggregates colliding) in the reduced mass, which is a negligible effect for clusters consisting of dozens of aggregates or more.
3.2. Determination of the aggregate collision velocities
During most phases of the experiment run (fast cell shaking), the optical depth in the SPACE suborbital flight experiment was too high and the image-recording frame rate too low to track aggregates individually. In these phases, it was not possible to determine collision velocities by following the aggregate positions over time. When decreasing the shaking speed under the sticking-bouncing transition, the experiment cells almost instantly got depleted in free-flying aggregates, which were then not available for individual aggregate tracking. Hence, only a very narrow span of time (~5 s during the shaking speed ramp-down of phase 2) allowed for aggregate tracking.
We therefore adopted a statistical approach to derive the aggregate velocities according to the shaking frequency of the experiment cells. The mean velocities derived from this statistical analysis for free-flying aggregates were then compared to a discrete set of aggregate velocities directly measured from the experiment data, when aggregate tracking was possible (shaking ramp-down of phase 2).
For this statistical approach, the residual accelerations during flight (see Appendix C.1) were neglected and the velocity of the dust aggregates was assumed to be solely induced by collisions with the container walls. The mean free path of the aggregates inside their cell volume was calculated by λ = 1/(nσ) with n being the number density of particles in the experiment cell and σ their mean mutual collision cross section, respectively. For the smaller size distribution counting ~4240 aggregates of a mean diametre of 120 μm and a cell volume of 15 × 10 × 24 mm3, the mean free path was calculated to be λ = 16 mm. For the larger size distributions counting ~375 aggregates of a mean diametre of 320 μm and a cell volume of 15 × 10 × 11 mm3, the mean free path was calculated to be 13 mm. Thus, in both cases the mean free path of the aggregates was comparable with the size of the cell. Considering collisions between aggregates and the cell walls to be perfectly inelastic, and knowing that the residual gas pressure in the cells was below 10-4 mbar during the experiment run (i.e. the influence of gas drag on the aggregate trajectories can be neglected), the free-flying aggregates were assumed to obtain and retain the maximum linear wall velocity until they collided with another aggregate or cluster on the cell walls.
To estimate the mean aggregate speed in the experiment cells during fast shaking phases, we considered that the configuration of the SPACE experiment, with clusters rotating at the speed of the cell walls (see Fig. 2 for the shaking profile), was similar to the one studied by Weidling et al. (2009) where an aggregate of known velocity collides with a wall moving vertically and sinusoidally with time. Therefore, we applied their method to statistically determine the probability of a certain relative velocity between the incoming free-flying aggregate and the rotating wall cluster. This method relies on the increased probability of aggregates to collide with clusters that are coming towards them in their current circular motion (clusters are following the motion of the experiment cells, see Brisset (2014) for details). The resulting probability distribution function can be seen in Fig. 6 (cumulative, red dashed line). The most frequent collision velocity was at twice the maximum linear wall velocity, 2vmax.
To confirm this modelled probability distribution, we measured the velocities of 51 free-flying aggregates by direct tracking during the short time interval in which this was possible (shaking speed ramp-down of phase 2). The obtained relative velocity distribution for collisions with clusters on the cell walls is shown in Fig. 6 (“+” signs). About 15% of the relative speeds were higher than 2vmax. This can be attributed either to semi-elastic collisions with the cell walls, or to extra kinetic energy stored in aggregate rotation that was released during the aggregate-wall collision. For the rest of the measured aggregate speeds, the velocity distribution is in very good agreement with the simple model. We can therefore assume a shaking speed dependant velocity distribution in the experiment cells at any time during the experiment run, as shown in Fig. 6. The Gaussian mean (location of the Gaussian distribution peak) and standard deviations of the velocity distribution measured on free-flying aggregates lie at 1.84. This value of the mean aggregate velocity will be used in the further data analysis.
![]() |
Fig. 6 Resulting relative collision velocity probability distribution function between an in-flying aggregate at vmax and a rotating cluster on the wall with linear velocity vmax (red dashed line). The cumulative distribution of 51 relative aggregate velocities measured during the shaking speed ramp-down of phase 2 is shown (“+” signs). The Gaussian mean (location of the Gaussian distribution peak) and standard deviations of the measured relative velocities are at 1.84 |
3.3. Number density of free-flying aggregates
The data processing and analysis steps are described in Appendix A. This data analysis showed that the background greyscale of the averaged recorded frame provides information on the number of free-flying aggregates in the inner cell volume at each moment. In the following, this is quantified to determine the number density of free-flying aggregates during the different phases of the experiment.
The background greyscale value of each averaged frame was determined as being the maximum of the normalized histogram of this frame (see Fig. A.2c., solid line). When no aggregates were flying in the cell volume, the background had the brightest possible value which is unity in the normalized greyscales shown in Fig. A.2c. When all aggregates were free-flying, like at the very beginning of the experiment (see Fig. A.2c., dashed red line), the maximum normalized greyscale value was only around 0.7, due to the shadowing effect of the free-floating aggregates.
This normalized background greyscale value, G, is related to the optical depth of the aggregate system τ, by G = e− τ. The optical density is in turn related to the aggregate number density in the experiment cell volume n, by τ = nσmonoL, where σmono is the projected cross section of one monomer aggregate and L is the depth of the experiment cell. Accordingly, the number density of free-flying aggregates is (1)G is measured as described above and L is known. σmono is calculated from the mean aggregate radius inside the experiment cell. The aggregate size distribution was measured in the laboratory (see Fig. 4, solid lines) and from the recorded experiment data during high velocity shaking phases (see Fig. 4, dashed lines). In these phases, individual aggregates cannot be tracked. However, the size distribution in still images can be measured. The data shows that the size distribution did not noticeably change during the flight: both cycle 1 and 3 final fast shaking phases show similar size distributions. As both the laboratory and the recorded size distributions are similar, the laboratory measurements can be used to asses σmono =
, Rmono being the measured aggregate radius (see Table 1 for measured radius values and their +/–1σ standard deviation).
The resulting number density of free-flying aggregate is shown in Fig. 7 as a function of time after start of the data recording. The slow shaking phase of the second cycle (see Fig. 2) can clearly be recognized at around 50 s, when no aggregates were free-flying in the volume. The two other slow shaking phases (cycles 1 and 3) were less pronounced at around 15 and 110 s, as both still had a high number of free-flying aggregates.
![]() |
Fig. 7 Retrieved number density of free-flying aggregates in the three cells of the SPACE experiment. The size distributions of large aggregates composed of polydisperse and monodisperse dust are shown as black solid and green dashed curves, respectively. The size distribution of small aggregates composed of monodisperse dust is shown as the red dash-dotted curve. The origin of time is at the start of data recording. The shaking cycles are indicated by dashed lines. |
3.4. Growth of clusters
When averaging the corrected frames (see Appendix A), the clusters growing on the glass walls become visible (see Fig. A.2b.). Fig. 8 shows the evolution of the total area covered by these clusters during the experiment run. During fast shaking phases, the clusters covered no or only small regions of the glass walls, while the maximum value of the cluster area was reached during the very slow shaking phase of cycle 2 at around 50 s after start of data recording.
![]() |
Fig. 8 Cell wall surface area covered by clusters during the SPACE experiment for the three dust types. The large size distribution of aggregates composed of poly- and monodisperse dust are shown by the black solid and green dashed curve, respectively. The small size distribution of aggregates composed of monodisperse dust is shown by the red dash-dotted curve. The origin of time is at the start of data recording. The shaking cycles are indicated by dashed lines. |
The absolute value of the area covered by clusters growing on the cell walls appears to be higher for the large size distribution of aggregates composed of monodisperse dust than for the two other types of dust. However, this has no influence on the further results of our analysis as only derivative values were used (see below). As this difference is not reflected in the background greyscale analysis, it is likely that more growth seeds were created on the glass surfaces perpendicular to the line of sight, compared to the two other cells. Surfaces parallel to the line of sight, which do not get accounted for in the absolute value of wall cluster surfaces, probably gathered more of the flying aggregates.
3.5. Aggregate pull-off forces
By analysing the evolution of the clusters on the experiment cell walls during the fragmentation phases (17 to 32 s, 62 to 76 s and 112 to 132 s, see Fig. 2), it was possible to determine the normalized cluster fragmentation rate (2)with n and
being the number density of free-flying aggregates in the experiment cell and its time derivate, respectively. Aggregates were indeed detaching directly from the clusters, as no aggregate rolling on the cell walls was observed. Sticking tests of dust on the coated cell glass before the flight, as well as the very low sticking probability of the aggregates on the cell glass during the flight (~0.5%, see Appendix B), showed that the pull-off force of aggregates from the cell walls is lower than the pull-off force between two or more dust aggregates inside a cluster. In the high acceleration environment experienced by the aggregates during the suborbital flight, aggregates detaching from the walls instead of from a dust cluster would first be rolling on the wall before detaching. This behaviour was not observed.
![]() |
Fig. 9 Normalized fragmentation rate of the clusters on the cell walls during the fragmentation phase of the second shaking cycle for the smaller and larger size distributions of aggregates composed of monodisperse SiO2 (red dash-dotted for small aggregates and green dashed line for large aggregates, respectively) and the larger size distribution of aggregates composed of polydisperse SiO2 (black solid line), plotted as a function of the centrifugal acceleration induced on the aggregates by the rotation of the cell wall. The curves are terminated at the acceleration where they reach the saturation of free-flying aggregates in the cell volume. The fragmentation rate maxima are shown by plus signs. |
Figure 9 shows the measured normalized fragmentation rate as a function of the centrifugal acceleration induced by the cell wall rotation during the second shaking cycle. The curves are terminated at the acceleration where a saturation of free-flying aggregates in the cell volume was reached. The normalized fragmentation rate increases with increasing centrifugal acceleration, indicating that more and more aggregates were pulled off clusters on the wall. The maxima (marked in Fig. 9) of these reduced fragmentation rates indicate the accelerations at which most of the aggregates composing the clusters on the cell walls were pulled off their parent clusters and, thus, allow the determination of a pull-off force for each investigated aggregate size and type. For the smaller aggregates composed of monodisperse dust, the linear wall acceleration at maximum fragmentation was ac = 16.6 m s-2. The corresponding pull-off force is N (see Table 1), with m being the mass of the aggregates pulled off the cluster. For the larger aggregates composed of mono- and polydisperse dust, the accelerations at maximum fragmentation were ac = 9.8 m s-2 and 9.1 m s-2, which corresponds to pull-off forces of
N and 1.6
N, respectively.
![]() |
Fig. 10 Monomer aggregate loss per cluster per frame (corresponding to 5.88 ms) for clusters on the cell walls during the fragmentation phase of the second shaking cycle for a) the small size distribution of aggregates composed of monodisperse SiO2; b) the large size distribution of aggregates composed of monodisperse SiO2; and c) the large size distribution of aggregates composed of polydisperse SiO2, respectively. |
To verify whether the calculated pull-off forces are realistic, we determined the number of monomer aggregates that each cluster lost between two consecutive image frames during the fragmentation phase. The result of this investigation is shown in Fig. 10. During the fragmentation phase, the shaking frequency of the experiment cells was increased over time, implying an increasing centrifugal acceleration level for the clusters on the wall. At centrifugal accelerations around the values for which the fragmentation rate of clusters reached its maximum (~63 s after start of data recording, see Fig. 10b.), the larger size distribution of aggregates composed of monodisperse dust displays a clear peak of monomer losses. This peak is also visible, although less marked, for the other two dust types. The aggregates composed of polydisperse dust lost up to 0.4 monomers per cluster and time interval between two frames (5.88 ms) at ~65 s (Fig. 10c.), meaning that not all clusters lost monomers at the same time. The small aggregates composed of monodisperse SiO2, however, lost up to about 4 monomers per cluster per frame at ~62 s (Fig. 10a.). This indicates that the small aggregates might have been pulled off the clusters in groups of four monomers instead of individually. The corrected aggregate radius and pull-off forces are listed in the last row of Table 1. Since the aggregates pulled off the clusters were four times larger in surface at identical acceleration, the force required to pull them off was about one order of magnitude higher than for single aggregates.
Accelerations at maximum fragmentation rate and pull-off forces Fpo for the three investigated types of dust.
3.6. Sticking probability
The analysis of the aggregate behaviour for increasing centrifugal accelerations shows that for accelerations under ~4 m s-2, fragmentation plays no role in the evolution of clusters of aggregates (see Fig. 9). During the second shaking cycle, aggregate sticking and cluster formation was observed at cell wall accelerations below 4 m s-2 (see Fig. 2). This means that the relatively high acceleration environment due to the shaking of the experiment cells did not influence the sticking behaviour of the aggregates. Its effects were only noticeable in the morphology of the formed clusters (see more details in Sect. 4.2 and Fig. 11). We therefore proceeded with the analysis of the sticking behaviour of the studied dust aggregates and measured their sticking probability as described in this section.
The number density profile determined in Sect. 3.3 can be used to derive the growth rate dn/ dt of clusters during shaking velocity ramp-down phases (6 to 15 s, 35 to 46 s and 85 to 112 s, see Fig. 2), which in turn can be used to calculate the sticking probability β of the coagulating system of aggregates.
Smoluchowski’s equation for particle coagulation can be adapted to the assumption that the growing clusters in the SPACE experiment were composed of a number i of single monomer aggregates (see Smoluchowski 1916; Blum 2006): (3)Here, n(i,t) is the number density of clusters composed of i monomer aggregates at a time t and Ka is the collision kernel for ballistic collisions between aggregates/clusters composed of i and j monomer aggregates, respectively. In this equation, the first term on the right-hand side accounts for the creation of clusters composed of i monomers by differently-sized aggregates, while the second term accounts for the depletion of clusters composed of i monomers by sticking collisions with other aggregates. The collisions taking place in the SPACE experiment cells during growth phases were mainly between monomer aggregates (j = 1) and clusters growing on the cell walls, thus depleting the inner cell volume of free-flying aggregates. Assuming that all clusters contain the same number of monomer aggregates at any given time, this implies that the first term on the right-hand side of Eq. (3) vanishes.
In the second term on the right-hand side, the collision kernel is defined as (4)Here, v(j,i) is the relative velocity between the clusters and the monomers, σ(j,i) the collision cross section between the two, and β(j,i;v) the sticking probability of their collision, respectively. For the collisions considered here, the sticking probability and relative velocity are considered to be the same for all aggregates at each moment. Hence, Eq. (3) becomes
(5)where n(t) is the number density of free-flying aggregates in the cell volume at time t determined in Sect. 3.3, dn(t) / dt is the derivative of the number density of aggregates incorporated in clusters on the cell walls determined in Sect. 3.4, β(t) and v(t) are the sticking probability and relative velocity of the collisions between the monomer aggregates and the clusters on the cell walls, σcross(t) is the mean collision cross section between a free-flying aggregate and a cluster on the wall, and n′(t) is the number density of clusters on the cell walls.
![]() |
Fig. 11 Typical wall cluster during the experiment from the large size distribution of aggregates composed of monodisperse dust. a) Original greyscale image. b) Contrast enhanced image. The red line in b) indicates the delimitation of the total cluster area. The packing density φp = 0.91 is determined by dividing the dark area by the area enclosed by the red line (see Sect. 4.2 for more details on the cluster morphology). |
The relative velocity v between the free-flying monomer aggregates and the rotating clusters on the cell walls has been determined in Sect. 3.2. The mean collision cross section σcross was calculated at each moment t of the experiment run by averaging over the collision cross sections of all clusters on the cell walls, (6)where N(t) = n′(t)V is the total number of clusters on the cell walls, V being the cell volume, and Rk and r are the radii of the wall cluster k and a monomer aggregate, respectively. In Eq. (6), the collision cross section σcross is taken as 2Rk(t)·2r. This assumes that the clusters on the cell walls have a thickness of one monomer aggregate and that the monomer aggregates approach them perpendicular to the observation direction, due to their motion induced by wall collisions (this collision configuration is shown by an artificially drawn incoming aggregate in the bottom left corner of Fig. 11a). To support this assumption, a typical wall cluster was chosen in the cell containing the larger size distribution of aggregates composed of monodisperse dust (see Fig. 11). In this cluster, the constituting monomer aggregates can very well be distinguished. This indicates that the cluster is very thin (most likely only one layer thick) and its packing density φp can be determined as the ratio between the area covered by aggregates to the total area of the cluster. For this cluster, the packing density is φp = 0.91. As the densest packing of circles on a surface is reached for hexagonal packing, where φp = 0.9069, it seems an appropriate approximation to assume that the clusters growing on the wall during the SPACE experiment possess a thickness of about one monomer (see Sect. 4.2 for more details on the cluster morphology). Thus, the quantity
(7)can be measured from the analysis of clusters growing on the wall during the SPACE experiment (see Sect. 3.4).
Finally, the sticking probability can be calculated as (8)All the quantities on the right hand side of this equation are measured at each frame during the experiment run. The resulting sticking probability between the monomer aggregates and the clusters on the cell walls calculated with Eq. (8) is shown in Fig. 12 as a function of the mean collision velocity and for all three aggregate types. During the growth phases of cycles 1 and 3, the sticking probability was very close to β = 0 because the wall speed was ~13 cm s-1 (see Fig. 2) so that the mean collision speed of 24 cm s-1 was much larger than the velocity range for sticking (see Fig. 12). Accordingly, only the sticking probability during the growth phase of cycle 2 is shown. The curves that can be distinguished growing in sticking probability with decreasing impact velocity are for the three different dust types during the growth phase of cycle 2. Below a certain velocity, the sticking probability of the dust aggregates rises very steeply and approaches β = 1. At this point, no free-flying aggregates are left in the cell volume. This velocity was determined to be the maximum velocity at which all aggregates are incorporated in clusters, leaving none free-flying in the cell volume. Table 2 lists the determined values for all three investigated types of aggregates. Both, the larger distributions of aggregates composed of mono- and polydisperse SiO2, showed perfect sticking at similar velocities of 11.5
and 11.7
cm s-1, respectively, while the smaller distribution of aggregates composed of monodisperse SiO2 stuck perfectly at 12.7
cm s-1 (error values are implied by the +/−1σ standard deviation on the collision velocities, see Sect. 3.2). The minimum velocity for which β< 0.05 was also determined. The smaller aggregates composed of monodisperse dust reached β = 0.05 at a mean collision velocity of 13.4
cm s-1 and the larger aggregates composed of mono- and polydisperse dust reached β = 0.05 at 12.2
and 12.6
cm s-1, respectively.
4. Discussion
During the suborbital flight, the dust aggregates inside the SPACE experiment cells were subjected to residual accelerations and potentially also to magnetic or electrostatic effects. We investigated the influence of these forces on the dust aggregates and present those in Appendix C. The result of this analysis is that aggregate motion induced by rocket residual accelerations, magnetic, or electrostatic effects can all be neglected compared to the accelerations induced by the cell shaking.
![]() |
Fig. 12 Sticking probability of the three types of dust aggregates investigated in the SPACE experiment plotted as a function of the statistical relative collision velocity computed in Sect. 3.2 (1.84vmax). The large size distribution of aggregates composed of poly- and monodisperse dust are shown by the black solid curve with squares and the green dashed curve with triangles, and the small size distribution of aggregates composed of monodisperse dust is shown by the red dash-dotted curve with diamonds. |
Perfect sticking velocities for all three aggregate types during the SPACE suborbital flight experiment.
In this section, we will start by comparing the present experiment results to previous work in Sect. 4.1, discuss the morphology of the formed aggregate clusters in Sect. 4.2 and derive the tensile strength of clusters and the surface energy of their aggregates in Sect. 4.3. Finally, conclusions will be derived for the behaviour of dust aggregates and clusters in protoplanetary disks in Sect. 4.4.
4.1. Comparing the SPACE results with previous work
The aggregate collisions investigated in this paper are different from the ones observed in previous work performed by Weidling et al. (2012) and Kothe et al. (2013). In the present work, aggregates did not collide while free-flying. Instead, one of the collision partners was attached to the glass cell wall and moving in a circular motion. Unlike in Weidling et al. (2012) and Kothe et al. (2013), this collision configuration prevented the formation of fractal clusters and collision energy did not dissipate through cluster restructuration or vibration. Instead, the relatively high acceleration environment at which the dust aggregates evolved during the SPACE experiment run imposed restructuring of the growing aggregates into compact mono-layered clusters (see Sect. 4.2 below). The sticking and fragmenting behaviour of the aggregates, however, was not influenced as we showed in Sects. 3.5 and 3.6. We could therefore deduce the aggregate and cluster properties described below in Sect. 4.3.
![]() |
Fig. 13 Flight collision results (dotted boxes) according to the aggregate size distributions (box heights) and their mean collision velocities during the experiment (box widths). As the larger aggregates composed of mono- and polydisperse SiO2 grains had very similar sticking properties, only the results for the aggregates composed of monodisperse SiO2 are plotted (upper box). The background colours correspond to the dust collision model developed by Güttler et al. (2010) with updates by Kothe et al. (2013). Green colour represents sticking collisions, yellow bouncing, and red fragmentation, respectively. The lines represent the transition zones between sticking and bouncing (solid line for a 50% sticking probability, dashed lines for 0, 25, 75 and 100% sticking probabilities) as well as the transition between bouncing and fragmentation (dash-dotted line for the onset of fragmentation), computed by Kothe et al. (in prep.). Mind that the collision model by Güttler et al. (2010) and Kothe et al. (2013) was developed for collisions between equal-mass aggregates, whereas the results presented here are dominated by single aggregates colliding with more massive clusters. |
Figure 13 compares the SPACE results with the dust collision model published by Güttler et al. (2010) with updates by Kothe et al. (2013). The collisions observed during the experiment run are represented by boxes with the box height given by their initial size distribution and the box width given by the range in mean relative collision velocity. The velocities for the onset of sticking (sticking probability β = 0.05) and for perfect sticking (β ≃ 1) listed in Table 2 are added to the boxes (vertical dotted lines), with the green and yellow regions denoting sticking and bouncing, respectively. As the larger aggregates composed of mono- and polydisperse dust had very similar sticking velocities, only the results for the monodisperse aggregates are plotted (upper box). The dust collision model is displayed in the background, delimiting the regions of the parametre range where sticking (green), bouncing (yellow), and fragmentation (red), respectively, dominate for collisions between similar-sized dust aggregates.
It can be noted that many sticking collisions occur in parametre ranges where bouncing is predicted by the current model. One explanation for this unexpected collision behaviour is the fact that the observed events are not aggregate-aggregate collisions, for which the model was developed, but aggregate-cluster collisions between free-flying aggregates and clusters growing on the cell walls. The mass ratio between colliding aggregates during the SPACE experiment covered up to 4 orders of magnitude. At the highest mass ratio between target cluster and projectile aggregate of ~104, only sticking collisions were observed. Kothe et al. (2013; see their Fig. 8) observed the enhanced sticking probability between clusters composed of a large number of aggregates, compared with the results of Weidling et al. (2012) who analysed collisions between aggregates of similar mass that resulted in bouncing for more than 90% of the collisions at velocities between 0.2 and 50 cm s-1. However, unlike both Kothe et al. (2013) and Weidling et al. (2012) who analysed collisions between same-sized aggregates, the data presented here reveals the role of the mass ratio between the colliding aggregates.
In Weidling et al. (2012) and Kothe et al. (2013), the coexistence of sticking and bouncing collisions in an extended region of the parametre space was observed. However, the rocket flight data display a very sharp transition from one regime to the other. This can be attributed to the statistical nature of the measurement: as the collision outcome is averaged over a high number of collisions, the transition appears sharper than for the individual collisions observed in the drop tower experiments.
4.2. Cluster morphology
The flat cluster structures found in our experiments (see Fig. 11) is due to an immediate restructuring of the clusters following a sticking collision. Due to the relatively high acceleration environment (≳1 m s-2) during the agglomeration phases (see Fig. 2), the centrifugal force experienced by the newly-acquired dust aggregate can potentially force individual aggregates to roll radially outward until no further rolling motion is possible. To evaluate whether this was the case in the SPACE experiment, we compare the centrifugal force with the rolling-friction force that counteracts any rolling motion. The latter was observed by Beitz et al. (2012; their Fig. 4), who witnessed the rolling of a dust aggregate above another under microgravity conditions and found that the initial velocity v ~ 2 cm s-1 was reduced to zero within a length of h ~ 4 mm. From these data we derive a rolling-friction acceleration of a ~ v2/ 2h = 0.1 m s-2, which is below the minimum centrifugal acceleration of 1 m s-2. Thus, we expect all aggregates to roll (but not to detach!) after sticking to the cluster.
The flattened, quasi-two-dimensional nature of the clusters, for which all constituent aggregates are in contact with the wall, considerably enhances the rigidity of the cluster. This, in turn, has consequences for the impinging dust aggregates. In contrast to the cluster-cluster collisions observed by Kothe et al. (2013) in which the weakly bound clusters could dissipate a much larger share of the impact energy, an impacting aggregate in our experimental configuration meets an aggregate that is attached to a stiff cluster. Hence, the physical situation is similar to an aggregate-aggregate collision, in which the impact energy is dissipated within the two aggregates in contact, with the exception that one of the aggregates has infinite mass (see also Brisset 2014, for more details). Therefore, the experiments performed here under low residual acceleration conditions are capable of deriving the sticking threshold of dust aggregates in mutual collisions but not the morphology of the clusters that would form under a perfect zero-gravity condition.
From numerical and experimental work by e.g. Blum et al. (2000), Krause & Blum (2004), Blum et al. (2006a), Blum (2006), Kothe et al. (2013), we know that growth through aggregate-aggregate and aggregate-cluster collisions in the sticking regime leads to very fluffy and fractal cluster structures. However, the collision environment of the protoplanetary nebula is not restricted to the hit-and-stick regime and can also lead to the formation of clusters that were rendered compact and stiff by bouncing collisions (Weidling et al. 2009; Zsom et al. 2010; Güttler et al. 2010). The structure of the clusters formed in the SPACE experiment is therefore, as well as the fractal aggregates observed in Kothe et al. (2013), representative of a particle population of the protoplanetary nebula.
4.3. Properties of clusters composed of sub-mm-sized aggregates
Here, the clusters built during the growth phases of the SPACE experiment will be investigated and some of their properties derived. Their inner cohesion is studied by investigating the tensile strength and the surface energy of their sub-mm-sized constituents.
Tensile strengths of aggregate clusters.
The pull-off forces required to detach a monomer aggregate from a wall cluster built during the slow shaking phase of the suborbital flight experiment are listed in Table 1. From these pull-off forces, the corresponding tensile strengths of clusters composed of monomer aggregates could be determined through T = Fpo/σagg, with σagg being the cross section of the detaching aggregate. We obtained 2.0 Pa and 2.1
Pa for the larger size distributions of aggregates composed of monodisperse and polydisperse dust, respectively. The tensile strength of agglomerates of four monomer aggregates composed of monodisperse dust is 2.7
Pa (see Sect. 3.5 for details).
The three values of the tensile strength are plotted in Fig. 14 as a function of the size of the detaching aggregates (diamonds and triangle for aggregates composed of mono- and polydisperse SiO2, respectively). The lines show the predicted values of the tensile strength of aggregates following the model by Skorov & Blum (2012), i.e. (9)where r and φ are the radius and the packing density of the aggregate. Skorov & Blum (2012) derived for the parametre T1 = 1.6 Pa. In Fig. 14, the dashed line is the model tensile strength computed for an aggregate packing density of φ = 0.3 as used by Skorov & Blum (2012) and the solid line for a packing density of φ = 0.37, the value assumed for the aggregates used in the SPACE experiment (see Sect. 2.2). Keeping in mind that the model predictions have no free parametre, the data from the SPACE experiment are in very good agreement with the predictions by Skorov & Blum (2012).
![]() |
Fig. 14 Tensile strength of aggregates in clusters investigated in the SPACE experiment for aggregates composed of mono- and polydisperse dust (diamonds and triangle, respectively) as well as of mm-sized dust aggregates measured by Blum et al. (2014; squares) in comparison to the model prediction by Skorov & Blum (2012). The model is applied for aggregate filling factors of 0.3 (dashed line) and 0.37 (solid line), respectively. |
Furthermore, Blum et al. (2014) also measured the tensile strength of SiO2 aggregate clusters, for aggregate diametres of 0.6 and 1.2 mm, respectively. Their results are plotted as squares in Fig. 14 and are in excellent agreement with the model by Skorov & Blum (2012) and the SPACE data. The concurrence of the theoretical model by Skorov & Blum (2012) with the results of two independent experimental methods is a strong support of the validity of the model for the derivation of the tensile strength of aggregate clusters. Please mind that the internal tensile strength of the dust-aggregate material with an individual aggregate is typically on the order of a few kPa (Blum et al. 2006b). The much lower tensile strength of adhering dust aggregates results from the relatively small contact area between the aggregates.
![]() |
Fig. 15 Aggregate pull-off forces as a function of their size for the two types of dust investigated in the SPACE suborbital flight experiment, i.e. monodisperse (diamonds) and polydisperse SiO2 (triangle). Linear fits through the data and the origin are also shown, with slopes of 7.6 × 10-4 (dashed line) and 9.8 × 10-4 (dotted line) for the two dust types, respectively. The smaller aggregates composed of monodisperse dust are assumed to detach in clumps of four aggregates. |
Estimation of the surface energy of macroscopic aggregates.
The pull-off forces listed in Table 1 are plotted in Fig. 15 as a function of the monomer aggregate size. For the smaller size distribution of aggregates composed of monodisperse dust, the detaching of agglomerates of four aggregates was assumed, as described in Sect. 3.5. The Johnson-Kendall-Roberts theory (Johnson et al. 1971) for the contact between soft spheres predicts a pull-off force (10)of one sphere from another as a function of the radius r. Here, γ is the surface energy of the material that forms the contact between the two spheres. Although this applies to perfect macroscopic soft spheres, the above relationship has been observed to be valid for μm-sized SiO2 particles as well (Heim et al. 1999). If this relationship is assumed to be valid for the monomer aggregates investigated by the SPACE experiment (mean diametres of ~120 and ~330 μm), an effective surface energy γeff must be introduced, as done in Weidling et al. (2012). Following their method, the effective surface energy can be scaled to the surface energy γ of the constituent micrometre-sized monomer particles via the aggregate filling factor φ and the Hertz factor
, where a0 is the radius of a monomer particle and a is the radius of the contact surface between two such monomer particles. This scaled effective surface energy can be written as
(11)where N is the number of connections from the separating aggregate to the cluster aggregates. The factor 2 accounts for the two sticking aggregates each having a surface energy of γ. As derived by Johnson et al. (1971; see Weidling et al. 2012, for more details),
(12)where E0 and ν are the monomer-particles’ Young’s modulus and Poisson number, respectively. Therefore,
(13)For the SiO2 dust used in the SPACE experiment, we assume ν = 0.17, E0 = 5.4 × 1010 Pa, a0 = 7.6 × 10-7 m and φ = 0.37. The number of connections from the separating aggregate to the cluster aggregates was determined from a typical cluster (Fig. 11) to be 2.5. In this cluster, each aggregate on the rim had between 2 and 3 neighbours.
Effective surface energies for the dust investigated in the SPACE rocket flight experiment.
The effective surface energy γeff was measured, to be 8.0 × 10-5 and 10.3 × 10-5 J m-2 for aggregates composed of mono- and polydisperse dust, respectively (see Fig. 15 and Table 3). For these values, the surface energy of a monomer particle was calculated to be γ = 1.7 × 10-2 J m-2 for the monodisperse and 2.0 × 10-2 J m-2 for the polydisperse SiO2 particles. These values are in very good agreement with 2.5 × 10-2 J m-2 measured for silica powder by Kendall et al. (1987) and 1.86 × 10-2 J m-2 measured for ~1 μm SiO2 particles by Heim et al. (1999). This confirms the validity of the scaling model developed in Weidling et al. (2012). This possible scaling from particles to aggregates could, for instance, be used in molecular dynamics simulations to investigate clusters composed of 100 μm-sized aggregates. Furthermore, it seems reasonable to try to adapt a collision recipe developed for clusters composed of μm-sized particles (Dominik & Tielens 1997) to clusters composed of ~100 μm-sized aggregates.
4.4. Application to protoplanetary disks
Three commonly used protoplanetary disk models were derived by Weidenschilling (1977b), Andrews & Williams (2007), and Desch (2007), respectively. These models are based on surface density (Σ) and temperature (T) profiles of the form (14)and
(15)where
is the distance to the central star, Σ0 and T0 are the values of the surface density and temperature at
1 AU, and δ and ϵ are the respective exponents of the power laws. All three models assume T0 = 280 K and ϵ = 0.5. Table 4 lists the surface density parametres. The Minimum Mass Solar Nebula (MMSN) model developed by Weidenschilling (1977b) is based on the minimum mass required in an original protoplanetary disk to build the Solar System as it exists today (Σ0 = 1700 g/cm2). Andrews & Williams (2007) derived a low density nebula model from the astronomical observations of circumstellar dust disks in Taurus-Auriga (Σ0 = 20 g/cm2) and Desch (2007) a high density model by considering planet migration following their formation in a compact configuration around the newly formed star (Σ0 = 50 500 g/cm2).
Model parametres for surface density of three protoplanetary disk models.
Turbulence in the disk is commonly assumed to follow a Kolmogorov cascade and the so-called α-prescription (Shakura & Sunyaev 1973), where the dimensionless parametre describes the strength of the turbulence, cs being the sound speed in the disk, and
and H the disk viscosity and height, respectively. Typical values for α are 10-5 in quiet regions of the disk, e.g. in a “dead-zone” that could develop around 1 AU (Turner et al. 2007; Brauer et al. 2008b), and up to 10-2 in more turbulent regions.
![]() |
Fig. 16 Relative velocities between dust aggregates in a protoplanetary disk computed according to a) Weidenschilling (1977b); b)Andrews & Williams (2007); and c)Desch (2007) at 1 AU and a turbulence parametre of α = 10-5. The velocity profiles are labelled in units of m s-1. The hashed boxes represent the collisions observed during the SPACE experiment (see text for more details). The dashed contour displays the minimum sticking velocity measured during the SPACE experiment for both size distributions. The dotted contour displays the same minimum sticking velocity in the case of aggregate-aggregate collisions instead of aggregate-cluster collisions (see text for details). |
To compute the dust aggregate relative velocities for different aggregate sizes in the PPD midplane, an adaptation of a code developed by Brauer et al. (2008a) was used. The effects considered to induce relative velocities were Brownian motion (Einstein 1905), radial drift (Whipple 1972; Weidenschilling 1977a) and turbulence (Ormel & Cuzzi 2007).
The computed relative velocity profiles for a low turbulence region (α = 10-5) can be seen in Fig. 16 for all three models described above at 1 AU. The aggregate and cluster size ranges investigated in the SPACE experiment are represented by hashed boxes, to indicate which “natural” relative velocities the involved aggregates would possess in a protoplanetary disk. The heights of the horizontal boxes represent the initial size distributions of aggregates inserted into the experiment cells and the widths of the horizontal boxes span from the monomer aggregate size to the mean size of clusters during slow shaking phases (see Sect. 4.1 for more details). In order to compare the relative velocities of aggregates in PPDs to the ones induced in the SPACE experiment, the minimum velocity for perfect aggregate sticking determined for both size distributions (see values in Table 2, dashed contour) was added to the figure.
As the collisions observed in the SPACE experiment occurred between individual aggregates and stiff aggregate clusters of infinite mass (see Sect. 4.2), the reduced mass of these collisions is the same than the mass of an individual aggregate. The change in momentum during the collision is therefore Δp = mv. This value of Δp could also be obtained by two aggregates of the same mass m (implying ) flying towards each other with the same velocity v (the relative velocity then would be 2v, see Brisset 2014, for more details). In Fig. 16, we also show the minimum sticking velocity between two same-size aggregates computed from our measurement of the sticking velocity between aggregates and stiff clusters (dotted line).
We can now see that for the MMSN (Weidenschilling 1977b, Fig. 16a.) and the compact (Desch 2007, Fig. 16c.) nebula models, the expected relative velocities for aggregate-aggregate and aggregate-cluster collisions are lower than the measured velocity for perfect aggregate sticking for the aggregates (both dashed and dotted lines). This indicates that at 1 AU and in a low-turbulence environment the investigated collisions would always lead to cluster growth. In the MMSN model, stiff clusters would stop growing through collisions with sub- to mm-sized particles when they reach about 4-mm in size (dotted line). The growth through same-sized aggregate-aggregate collisions, however, continues up to sizes of around 5 cm (dotted line in Fig. 16a). For the compact nebula model (Desch 2007, Fig. 16c.), growth of stiff clusters through aggregate sweep-up would only lead to cluster sizes of up to 10 cm, while aggregate-aggregate collisions could grow clusters up to about 50 cm in size. In the low density model (Andrews & Williams 2007, Fig. 16b.), growth of stiff clusters through aggregate sweep-up would only lead to cluster sizes between 0.1 and 0.2 mm, while aggregate-aggregate collisions could grow clusters up to about 1 mm in size. Above these maximum size listed, bouncing is expected.
5. Conclusion
In this paper we presented the results of a multi-particle collision experiment conducted under reduced gravity conditions. The particles investigated were sub-mm-sized aggregates composed of micrometre-sized spherical and irregular SiO2 dust particles. The experiment run on the REXUS 12 suborbital rocket provided 150 consecutive seconds of data observing the collision behaviour of these aggregates. The reduced gravity environment of the suborbital rocket was necessary to reach relative velocities between aggregates as low as <1 cm s-1 and to investigate the aggregate collision behaviour over a continuous range of relative velocities. During the SPACE experiment run, clusters of up to a few millimetres in size were formed. Owing to the rather high acceleration environment on the suborbital flight (≳1 m s-2, owing to the shaking of the experiment cell) the formed clusters restructured into monolayers through rolling of the aggregates upon sticking. Data analysis showed that this behaviour did not affect the measurements of the aggregate sticking probability and their pull-off force.
A direct measurement of the sticking probability for impact velocities from about 3 to 22 cm s-1 was possible (Fig. 12). The sticking probability of the dust aggregates was shown to be velocity dependant and the transition from a bouncing regime (sticking probability close to 0) to a sticking regime (sticking probability at 1) was very sharp. This indicates the existence of a sticking threshold for relative velocities between aggregates.
The suborbital flight also allowed for the study of the properties of sub-mm-sized aggregates built into clusters. The sequenced growth and fragmentation phases made it possible to measure the tensile strength and the effective surface energy of the aggregate clusters. The tensile strength values measured during the SPACE experiment are on the order of a few Pa, in agreement with independent measurements (Blum et al. 2014) and model predictions (Skorov & Blum 2012). The value for the effective surface energy of the aggregates is in very good agreement with the surface energy measured for ~1 μm particles (Heim et al. 1999), combined with the scaling model of Weidling et al. (2012). This indicates that dust aggregates can be used as monomer particles in molecular dynamics simulations, with a proper scaling of the particle parametres. This allows for simulations of clusters up to 10 cm in size with the computational capacities currently used.
As shown in Fig. 16, the measured sticking threshold for the aggregates studied in this work leads to different dust growth behaviour predictions depending on the protoplanetary disk model considered. For a disk dust surface density at 1 AU of 1700 g cm-2 (Weidenschilling 1977b), the results of the experiment presented in this work indicate a dust grain growth to sizes of about 1 cm at 1 AU. In the case of a lower density model (Andrews & Williams 2007), dust grains would grow only to a few 100 μm, while a higher density model (Desch 2007) would allow for uninterrupted dust grain growth.
Acknowledgments
We thank the REXUS/BEXUS project of the Deutsches Zentrum für Luft- und Raumfahrt (DLR) for the flight on the REXUS 12 rocket. This work was supported by the Interaction in Cosmic and Atmospheric Particle Systems (ICAPS) project of DLR (grants 50WM0936 and 50WM1236) and a fellowship from the International Max Planck Research School on Physical Processes in the Solar System and Beyond (IMPRS) at the Universities of Braunschweig and Göttingen. We also thank Dipl. Phys. Oliver Lenck from the Frauenhofer Institute for Surface Engineering and Thin Films of Braunschweig for the anti-adhesive glass coating of the particle containers.
References
- Anderson, M., Orsan, H., Koivusalo, R., Laukkanen, M., & Rigaldo, A. 2009, Reverse Engineering the REXUS Rocket, http://perso.centrale-marseille.fr/~arigaldo/visible/KTH_projects/SPACE_REXUS_Rocket.pdf [Google Scholar]
- Andrews, S. M., & Williams, J. P. 2007, ApJ, 659, 705 [NASA ADS] [CrossRef] [Google Scholar]
- Beitz, E., Güttler, C., Weidling, R., & Blum, J. 2012, Icarus, 218, 701 [NASA ADS] [CrossRef] [Google Scholar]
- Blum, J. 2006, Adv. Phys., 55, 881 [NASA ADS] [CrossRef] [Google Scholar]
- Blum, J., & Wurm, G. 2008, ARA&A, 46, 21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Blum, J., Wurm, G., Kempf, S., et al. 2000, Phys. Rev. Lett., 85, 2426 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Blum, J., Bruns, S., Rademacher, D., et al. 2006a, Phys. Rev. Lett., 97, 230601 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Blum, J., Schräpler, R., Davidsson, B. J. R., & Trigo-Rodriguez, J. M. 2006b, ApJ, 652, 1768 [NASA ADS] [CrossRef] [Google Scholar]
- Blum, J., Gundlach, B., Mühle, S., & Trigo-Rodriguez, J. M. 2014, Icarus, 235, 156 [NASA ADS] [CrossRef] [Google Scholar]
- Bouwman, J., Henning, T., Hillenbrand, L. A., et al. 2008, ApJ, 683, 479 [NASA ADS] [CrossRef] [Google Scholar]
- Brauer, F., Dullemond, C. P., & Henning, T. 2008a, A&A, 480, 859 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Brauer, F., Henning, T., & Dullemond, C. P. 2008b, A&A, 487, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Brisset, J. 2014, Ph.D. Thesis, Technische Universität Braunschweig, http://www.mps.mpg.de/3718395/Dissertation_2014_Brisset__Julie1.pdf [Google Scholar]
- Brisset, J., Heißelmann, D., Kothe, S., Weidling, R., & Blum, J. 2013, Rev. Sci. Instr., 84, 094501 [Google Scholar]
- Desch, S. J. 2007, ApJ, 671, 878 [NASA ADS] [CrossRef] [Google Scholar]
- Dominik, C., & Tielens, A. G. G. M. 1997, ApJ, 480, 647 [NASA ADS] [CrossRef] [Google Scholar]
- Dominik, C., Blum, J., Cuzzi, J., & Wurm, G. 2007, Protostars and Planets V, 783 [Google Scholar]
- Einstein, A. 1905, Annalen der Physik, 322, 549 [Google Scholar]
- Garaud, P., Meru, F., Galvagni, M., & Olczak, C. 2013, ApJ, 764, 146 [NASA ADS] [CrossRef] [Google Scholar]
- Güttler, C., Krause, M., Geretshauser, R., Speith, R., & Blum, J. 2009, ApJ, 701, 130 [NASA ADS] [CrossRef] [Google Scholar]
- Güttler, C., Blum, J., Zsom, A., Ormel, C. W., & Dullemond, C. P. 2010, A&A, 513, A56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Haff, P. K. 1983, J. Fluid Mech., 134, 401 [NASA ADS] [CrossRef] [Google Scholar]
- Heim, L.-O., Blum, J., Preuss, M., & Butt, H.-J. 1999, Phys. Rev. Lett., 83, 3328 [NASA ADS] [CrossRef] [Google Scholar]
- Johansen, A., Youdin, A., & Klahr, H. 2009, ApJ, 697, 1269 [NASA ADS] [CrossRef] [Google Scholar]
- Johansen, A., Youdin, A. N., & Lithwick, Y. 2012, A&A, 537, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Johansen, A., Blum, J., Tanaka, H., et al. 2014, Protostars and Planets VI, 547 [Google Scholar]
- Johnson, K. L., Kendall, K., & Roberts, A. D. 1971, Proc. R. Soc. Lond., 324, 301 [Google Scholar]
- Kataoka, A., Tanaka, H., Okuzumi, S., & Wada, K. 2013, A&A, 557, L4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kendall, K., McN. Alford, N., & Birchall, J. D. 1987, Nature, 325, 794 [NASA ADS] [CrossRef] [Google Scholar]
- Kothe, S., Blum, J., Weidling, R., & Güttler, C. 2013, Icarus, 225, 75 [NASA ADS] [CrossRef] [Google Scholar]
- Krause, M., & Blum, J. 2004, Phys. Rev. Lett., 93, 021103 [NASA ADS] [CrossRef] [MathSciNet] [PubMed] [Google Scholar]
- Nübold, H., Poppe, T., Rost, M., Dominik, C., & Glassmeier, K. 2003, Icarus, 165, 195 [NASA ADS] [CrossRef] [Google Scholar]
- Ormel, C. W., & Cuzzi, J. N. 2007, A&A, 466, 413 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Poppe, T., Blum, J., & Henning, T. 2000, ApJ, 533, 472 [NASA ADS] [CrossRef] [Google Scholar]
- Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 [NASA ADS] [Google Scholar]
- Skorov, Y., & Blum, J. 2012, Icarus, 221, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Smoluchowski, M. V. 1916, Physik. Z., 17, 557 [Google Scholar]
- Testi, L., Birnstiel, T., Ricci, L., et al. 2014, Protostars and Planets VI, 339 [Google Scholar]
- Turner, N. J., Sano, T., & Dziourkevitch, N. 2007, ApJ, 659, 729 [NASA ADS] [CrossRef] [Google Scholar]
- Wada, K., Tanaka, H., Suyama, T., Kimura, H., & Yamamoto, T. 2011, ApJ, 737 [Google Scholar]
- Weidenschilling, S. J. 1977a, MNRAS, 180, 57 [NASA ADS] [CrossRef] [Google Scholar]
- Weidenschilling, S. J. 1977b, Ap&SS, 51, 153 [Google Scholar]
- Weidenschilling, S. J. 2000, A&A, 92, 295 [Google Scholar]
- Weidenschilling, S. J., & Cuzzi, J. N. 1993, in Protostars and Planets III, eds. E. H. Levy, & J. I. Lunine, 1031 [Google Scholar]
- Weidling, R., Güttler, C., Blum, J., & Brauer, F. 2009, ApJ, 696, 2036 [NASA ADS] [CrossRef] [Google Scholar]
- Weidling, R., Güttler, C., & Blum, J. 2012, Icarus, 218, 688 [NASA ADS] [CrossRef] [Google Scholar]
- Whipple, F. L. 1972, in From Plasma to Planet (New-York: Wiley), 211 [Google Scholar]
- Windmark, F., Birnstiel, T., Güttler, C., et al. 2012a, A&A, 540 [Google Scholar]
- Windmark, F., Birnstiel, T., Ormel, C. W., & Dullemond, C. P. 2012b, A&A, 544, L16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Youdin, A. N., & Goodman, J. 2005, ApJ, 620, 459 [NASA ADS] [CrossRef] [Google Scholar]
- Zsom, A., Ormel, C. W., Güttler, C., Blum, J., & Dullemond, C. P. 2010, A&A, 513, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Appendix A: Analysis methods
Frame correction.
At the beginning of the data analysis sequence, the recorded frames were corrected as follows:
-
Elimination of the rotational movement of the experiment cells. This was done by detecting the cell limits in each frame and shifting the frame so that these limits were always in the same position. In these corrected frames, the free-floating particles possess cycloid motions instead of linear trajectories. The cells had an additional movement along the line of sight of the camera because they were slightly wobbling due to friction between the cog wheels of the shaking mechanism. This produced a non-zero noise level in the greyscale histogram of averaged frames (see Fig. A.2c.) that could be disregarded as not being related to the number of free-flying aggregates.
-
Correction of spatial variations in background illumination. The background illumination was not perfectly even over the field of view of the camera. To rectify this, the average background was determined over all the recorded frames and the mean value then subtracted from each frame (see Fig. A.1b.).
-
Checking for temporal variations in the transmission of the glass windows. As the number of collisions between the free-flying dust aggregates and the glass walls was very high and such collisions could pollute the glass walls with small dust particles, it was necessary to check for the glass cleanliness over the experiment run. This was done by picking a few random positions on a frame and checking their greyscale values at several moments in time during the experimental run when they were not covered by an aggregate. The result of this investigation was that the glass stayed clean in each of the three experiment cells confirming the low glass sticking efficiency determined in Appendix B. This also means that the anti-adhesive coating worked well on more than 99% of the glass surfaces and that the deposition of seed aggregates on the glass walls was probably due to small failures of the anti-adhesive coating.
The thus corrected frames were then used for the analysis of aggregate collisions and growth.
Frame averaging.
Because of the high optical depth of the free-flying aggregates in the experiment cells, it was not possible to measure their number and individual sizes directly. A way of still accessing this information was to perform a running average of the image frames. Because the free-floating aggregates were constantly moving, the total cross-sectional area of each aggregate was thus taken into account by contributing to the background greyscale. Furthermore, non-moving clusters growing on the glass walls became fully visible. Thus, each frame of the SPACE experiment was averaged over 201 frames, i.e. over a time of 1.2 s. An example of the frames obtained by this averaging can be seen in Fig. A.1c. The clusters on the glass walls are now apparent as dark spots and can be monitored. Furthermore, the background greyscale value of cluster-free parts of the frames delivers information on the number of dust aggregates that are free-flying, i.e. the darker the background, the more aggregates are free-flying, and the brighter the background, the more aggregates are incorporated in clusters on the cell walls.
![]() |
Fig. A.1 Image processing steps. a. Recorded SPACE data frame at 11.8 s after beginning of frame recording, showing one of the smaller experiment cells. b. The same frame after correction of background illumination. c. The same frame after background correction and averaging over 201 frames. The clusters on the cell walls become apparent. |
Statistical analysis methods.
The analysis of the frame background greyscale and of the clusters on the cell walls are in fact two distinct ways of retrieving information on the number and mass of the aggregates and clusters in the experiment. Figure A.2 shows histograms of two averaged data frames. In the left frame a., all of the dust aggregates present in the cell are free-flying. Its greyscale histogram (c, dashed red line) shows a rather diffuse distribution of medium-dark greyscale values. In the right frame b, however, all aggregates are incorporated in clusters and none are free-flying. The histogram now shows two distinct peaks (c, solid black line), a strong and narrow peak for high greyscale values accounting for the now very bright (corrected) background, and a smaller but also quite narrow peak at very low greyscale values accounting for the dark visible clusters on the cell walls. As the glass remained perfectly clean during the experiment run, these two frames can be used to calibrate the background greyscale to the fraction of aggregates free-flying in the cell volume at this time (see Sect. 3.3).
Another analysis method is to binarize the averaged frames and monitor the area covered by clusters at a certain point in time. Figure A.3 illustrates this by presenting two averaged and binarized frames at the beginning and the end of a fragmentation phase. Clusters (black areas) have partly dissolved over time. The total surface covered by clusters can be tracked during the experiment run (see Sect. 3.4).
![]() |
Fig. A.2 Analysis method used on the SPACE experiment data to determine the number of free-flying aggregates in the experiment cells at each moment in time. a) Region of interest inside an averaged data frame of the larger experiment cell (see Fig. 1b) at the beginning of the experiment run (1.2 s after start of recording), when all aggregates were free-flying. Image processing as shown in Fig. 7c was applied. b) Region of interest inside an averaged data frame of the same experiment cell during the slowest shaking phase (47.1 s after start of recording), when all aggregates are incorporated in clusters sticking to the glass walls. Image processing as shown in Fig. 7c was applied. c) greyscale histograms of frames a) and b) When all aggregates are free-flying a), the background is uniformly grey (dashed red line), whereas the dark clusters (small peak at lower greyscale values) and the bright background (high and narrow peak at higher greyscale values) can clearly be distinguished when the aggregates are all incorporated in clusters (solid black line). |
![]() |
Fig. A.3 Analysis of the clusters on the cell walls inside the SPACE experiment during the suborbital flight. a) Region of interest inside an averaged and binarized data frame of the bigger experiment cell at the beginning of a fragmentation phase (61.8 s after beginning of the recording). b) Region of interest inside an averaged and binarized data frame of the same cell at the end of this fragmentation phase (76.5 s after beginning of the recording). The area covered by clusters (black areas) has visibly decreased. |
Appendix B: Investigation of the cluster growth by aggregate-aggregate collisions
This appendix describes the Monte Carlo simulations performed depositing particles on a surface representing the glass walls of the SPACE experiment cell. The purpose of these simulations was to reproduce the aggregate deposition pattern seen on the cell walls and to investigate whether all aggregates stuck directly to the glass or tended to form clusters through dust-dust collisions. The parametres that were varied between the simulations were the initial number of particles to be deposited on the glass surface, Nsim, and a sticking probability for aggregate-glass collisions, κ.
Before the simulations are started, no particles are deposited on the surface. A number Nsim of particles are then created with a size distribution similar to the size distribution of the dust aggregates and inserted into the experiment cells. At each step of the simulation, one of the Nsim particles collides with a randomly chosen location on the surface. A collision with the bare glass surface either leads to a deposition of the aggregate according to the sticking probability κ, or to bouncing back into the cell volume. In the latter case, the particle collides with the surface again, on a new random position. A collision with another particle already sticking to the surface always leads to deposition (a sticking probability of unity is considered for aggregates impinging those sticking to the glass wall). The simulation stops when all Nsim particles are deposited.
As the frames recorded during the experiment had a resolution of 56.5 μm/pixel, it was possible to simulate the particle growth on a one-to-one scale. The smaller sieved aggregates observed in the SPACE experiment had a mean size of 120 μm ranging from ~20 to 300 μm in diametre and the bigger aggregates had a mean size of 330 μm ranging from ~50 to 600 μm (see Fig. 4). Hence, 1 to 81 pixel particles inserted into the simulation correspond to the size range of dust aggregates observed in the SPACE experiments. Simulating the wall surfaces by 221 × 362 pixel and 192 × 141 pixel corresponds to regions of interest inside the 15 × 24 mm2 and 15 × 11 mm2 visible surfaces of the larger and smaller experiment cells.
Monte Carlo simulation results: the best match between simulation and data for both dust distributions are presented.
For the small monodisperse SiO2 dust distribution, the real inserted number of aggregates was Nexp ≃ 4000. Thus, the initial number of particles created in the simulations was varied from Nsim = 3000 to Nsim = 5000, with a step size of 100 particles. For the larger size distributions of aggregates composed of mono-and polydisperse dust, the real number of inserted aggregates was Nexp ≃ 400. Here, the initial number of particles was varied from Nsim = 300 to Nsim = 500 with a step size of 10 particles. For both distributions, the glass sticking parametre κ was varied from κ = 0.0001 to κ = 0.1 with a linear step size of 0.0001.
A few examples of the simulation results can be seen in Figs. B.1 and B.2. For low sticking probabilities κ, only very few seeds deposited on the surface and clusters grow from there, leading to a few large clusters. For higher values of κ, many more particles deposit on the glass surface and the simulation results in a higher number of smaller clusters.
![]() |
Fig. B.1 Summary of the Monte Carlo simulation results for the agglomeration of the small monodisperse dust size distribution used in the SPACE experiment rocket flight (mean size of 120 μm). The simulation surface of 221 × 362 pixel represents a region of interest inside the bigger experiment cell. Each column has the same initial number of particles inserted into the simulation Nsim. Each row has the same glass sticking probability κ. |
![]() |
Fig. B.2 Monte Carlo simulation results for the agglomeration of the large size distribution of aggregates composed of monodisperse dust used in the SPACE experiment rocket flight (mean size of 330 μm). The simulation surface of 192 × 141 pixel represents a region of interest inside one of the smaller experiment cells. Each column has the same initial number of particles inserted into the simulation Nsim. Each row has the same glass sticking probability κ. |
To quantitatively compare the results of the simulation with the experimental data, a data frame taken at 47 s after start of recording was used. This frame was recorded during a slow shaking phase (see Fig. 2), while all particles were incorporated in clusters and none were free-floating. This frame was binarized and the resulting number of clusters and their size distribution were determined. The root mean square of the bin-by-bin difference between the measured cluster size distribution and the size distribution of the simulation result, σdist, was used as the comparison criterion. To smooth out binarization effects on the size distribution of the data frame, the bin sizes were chosen at 4 pixel for the small and 9 pixel for the larger dust distribution. The simulations were run 5000 times for each set of parametres and the mean σdist determined.
The best matches (minimum σdist) were found for Nsim = 4600 and Nsim = 360 and κ = 0.0055 and κ = 0.0045 for the smaller and larger aggregate distributions, respectively (Table B.1). Figure B.3 illustrates the best match between data and simulations for the smaller aggregate size distribution. The size distributions of 5000 simulation runs and of the measured data frame are plotted in grey and in light red, respectively. It can be seen that the size distribution of the data frame lies within the size distribution range obtained for the simulations. As the number Nsim of inserted particles is in good agreement with the number of aggregates inserted into the SPACE experiment cells, glass sticking values around 0.5% can be assumed.
It can be concluded that it was not an enhanced glass sticking efficiency that led to clusters on the cell walls but indeed the enhanced number of collisions of the aggregates on the glass, induced by residual accelerations during the rocket flight (see Table C.1). This is also supported by the fact that the glass walls did not get dirty during the experiment run despite the very high number of collisions. It can also be concluded that dust clusters were not growing on the glass directly but on a few dust aggregates that had stuck to failure points of the anti-cohesive coating of the glass walls and then served as seeds.
![]() |
Fig. B.3 Results of a Monte Carlo simulation for the smaller aggregate size distribution used in SPACE during the suborbital flight: a) binarized SPACE data frame at 47 s after start of data recording. At this point of the experiment run, all dust aggregates are clustered on the cell wall. b) Monte Carlo simulation depositing 4600 particles on a surface with a glass sticking parametre κ = 0.0055, c) size distribution of the data frame seen in a. (light red line, triangles) and 5000 runs of Monte Carlo simulations with the parametres of image b) (grey lines). The bin size for these distributions is 4 pixel (1.23 × 10-8 m2). |
Appendix C: Calculating perturbation accelerations during the suborbital flight
Appendix C.1: Residual accelerations during the rocket flight
During the suborbital flight of SPACE, the dust aggregates in the containers were subjected to two types of residual acceleration, caused by atmospheric drag and rocket spin.
The residual atmospheric drag accelerations were due to the fact that the REXUS 12 rocket had an apogee of 82 km, an altitude at which the aerodynamic effects of the remaining atmosphere on the rocket are still influential. The drag force induced along the direction of flight is , with ρ being the local air pressure, v, A and Cd the rocket’s flight velocity, cross section and drag coefficient, respectively. Comparing the rocket’s flight altitude profile with the recorded SPACE data, it could be determined that the streaming of the image frames started at an altitude of about 70 km, lasted through the passing of apogee (82 km) and stopped during descent at about 50 km altitude. Hence, it can be assumed that the relevant data to used for this analysis were recorded above an altitude of 70 km (the last seconds of recording were not relevant for the data analysis). The standard atmosphere model (1976 US Standard Atmosphere) at this altitude yields ρair = 7.42 × 10-5 kg/m3. The cylindrical rocket had a diametre of d = 0.356 m and, hence, a cross section of
m2. The drag coefficient Cd = 0.341 was calculated from measurements made on a REXUS flight in 2009 made by Anderson et al. (2009). The maximum rocket speed was v = 562 m s-1 at 70 km altitude, which leads to a maximum drag force on the rocket of F = 0.398 N. For a rocket mass of mr = 515 kg, the maximum residual drag acceleration was adrag = Fdrag/mr = 7.72×10-4 m s-2 = 7.86 × 10-5g, with g being the Earth’s gravitational acceleration. As the dust aggregates were free-flying in vacuum inside of the experiment cells, which are directly coupled to the rocket itself, this is the maximum residual drag acceleration they experienced relative to the cell walls.
The residual angular velocity of the rocket after de-spin was measured by the Service Module to be ωr = 0.19 rad/s. The resulting centrifugal acceleration on one of the SPACE dust aggregates was, hence, with R being the distance of the aggregate to the roll axis of the rocket. The way the experiment was designed, the distance to the rocket’s roll axis could be approximated to R = 0.04 m. Thus, the residual spin acceleration on the aggregates was about acentrif = 1.47 × 10-4g.
The effects of these residual accelerations on the behaviour of the dust aggregates can be seen in Fig. 5b. The tendency of the aggregates to gather in the upper left corner of the frame reveals the combination of a residual acceleration acting on them in the direction of flight (vertical on the image frames) and in the radially outward direction (horizontal on the image frames), increasing the number of their collisions in this corner of the cell.
To investigate whether these two residual accelerations had an influence on the collision behaviour of the observed aggregates, the acceleration values were compared to the those induced by the rotating cell walls. The minimum rotation frequency of the cells was fmin = 4.78 Hz, inducing a minimum wall acceleration of m s-2 = 9.2 × 10-2g, with
mm being the rotation radius of the experiment cells. This is almost 3 orders of magnitude stronger than the residual accelerations calculated above. Therefore, it is reasonable to neglect the effects of residual accelerations in the data results.
Appendix C.2: Magnetic and electrostatic effects
To validate our results as being meaningful for protoplanetary dust aggregation, it is also important to rule out magnetic or electrostatic effects on the collision properties. Both magnetic and electrostatic effects are known to play an important role on dust particle collisions (e.g. Poppe et al. 2000; Nübold et al. 2003). However, in this series of experiments with SiO2 dust aggregates, the intention was to investigate aggregate properties free of these influences.
Ruling out magnetic effects is trivial as none of the materials used for this experiment were magnetic.
Furthermore, to avoid electrostatic charging of any kind, the glass cells holding the particles were built inside of a frame acting as a Faraday cage and preventing external electric fields penetrating inside the cells. Before inserting the dust aggregates into the containers, they were grounded for several days to avoid initial particle charging.
![]() |
Fig. C.1 Reduced χ2 error statistics for linear and parabolic fits to the trajectories of 10 trackable aggregates in the cell containing the larger aggregates composed of monodisperse dust during the SPACE rocket flight experiment run, before (triangles, open for the linear fit and filled for the parabolic fit) and after impact with a glass wall (squares, open for the linear fit and filled for the parabolic fit). |
To make sure that these precautions guaranteed negligible electrostatic effects, the charge state of the aggregates observed during the SPACE flight was derived from the analysed data. During the SPACE experiment run, the relevant collisions were between free-flying aggregates and clusters forming on the cell walls. Electrostatic effects were ruled out by picking a few
random aggregates, which could be directly observed colliding with a container wall, and determining the form of their trajectories. In the case of electrostatic charging of the dust and the cell walls, the aggregates would behave like charged particles in an electric field: their trajectories approaching the cell wall would be non-linear with respect to time due to induced acceleration or deceleration. For the 10 analysed collisions, the reduced χ2 error statistics of both the linear and parabolic fit of their trajectories approaching and leaving the glass wall was measured (see Fig. C.1). The mean χ2 error statistics for a linear trajectory fit were 0.083 and 0.040 before and after the aggregate collision with the cell wall, respectively. For the parabolic fit, these values were 0.90 and 0.12, respectively, indicating that there was no non-linearity in the particle trajectories. The second degree coefficient of the formal parabolic fit, p2, can be related to an electrical charge of the corresponding aggregate. Each aggregate is considered as being a charged particle of charge q moving in a uniform electric field E, then , with m being the mass of the aggregate. The electric field is approximated to
with ϵ0 being the dielectric constant in vacuum and r the distance from the aggregate to the interacting charge (taken here as ~2 × aggregate radius). The mean number of elementary charges carried by each monomer particle (~1 μm in size) inside the aggregate is thus calculated to be 0.015. The induced accelerations are ael = 7.27 × 10-11 m s-2. These accelerations can be neglected compared to the ones induced by the cell shaking.
All Tables
Accelerations at maximum fragmentation rate and pull-off forces Fpo for the three investigated types of dust.
Perfect sticking velocities for all three aggregate types during the SPACE suborbital flight experiment.
Effective surface energies for the dust investigated in the SPACE rocket flight experiment.
Monte Carlo simulation results: the best match between simulation and data for both dust distributions are presented.
All Figures
![]() |
Fig. 1 Experiment cells inside the rocket payload. a) Schematic of the REXUS 12 rocket with the position of the SPACE module marked by a red asterisk. b) CAD drawing of the SPACE payload including the vacuum chamber (1) and the camera (2) observing the interior of the chamber through a mirror (3). c) Schematic view of the glass experiment cells (5) installed on the shaking mechanism (6). The gears of the shaking mechanism activated by a motor induced the circular motion of the cells indicated on the figure by blue arrows. The circle indicates the limits of the vacuum chamber viewport. d) Photograph of the SPACE experiment cells inside the vacuum chamber and seen through the view port. The types of dust inserted into each cell were ~120 μm-sized monodisperse in the larger cell (I) and ~330 μm-sized poly- and monodisperse in each of the smaller cells (II and III, resp.). |
In the text |
![]() |
Fig. 2 Measured rotation frequency of the particle cells during the experimental run of SPACE onboard the REXUS 12 rocket. Time is measured from the first recorded frame (~100 s after lift-off of the rocket). The corresponding wall centrifugal acceleration and maximal linear speed are indicated. The apogee of the rocket trajectory is indicated by a “+” sign. |
In the text |
![]() |
Fig. 3 SiO2 dust used in the SPACE experiment: a) SEM image of the polydisperse SiO2 monomer particles (image credit: E. Beitz); b) SEM image of monodisperse SiO2 monomer particles (image credit: E. Beitz). c) Microscope image of aggregates composed of irregular polydisperse SiO2, sieved between 250 and 500 μm. |
In the text |
![]() |
Fig. 4 Size and mass distributions of the aggregate samples used in this study. a) Normalized cumulative number of aggregates with a size smaller than the size indicated on the x-axis. b) Normalized cumulative mass of aggregates smaller than the size indicated on the x-axis for the three different SiO2 aggregates used in the SPACE experiment. Black diamonds denote aggregates sieved between 100 and 250 μm and composed of monodisperse dust particles. Black triangles denote aggregates sieved between 250 and 500 μm and composed of monodisperse dust particles. Green asterisks denote aggregates sieved between 250 and 500 μm and composed of polydisperse dust. Dashed lines denote aggregate sizes measured during the experiment run on the rocket. |
In the text |
![]() |
Fig. 5 Two image frames recorded during the SPACE experiment. a) At the very beginning of the experimental run, when all the dust aggregates were free-flying. b) During a slow shaking phase of the experiment (see Fig. 2). In slow shaking phases, the preferential accumulation of aggregates in the cells’ left upper corner due to microgravity disturbances (see C.1 for details) is visible. The flight direction of the rocket (i.e. direction of gas drag) and the outward radial direction (i.e. direction of centrifugal force) are indicated. |
In the text |
![]() |
Fig. 6 Resulting relative collision velocity probability distribution function between an in-flying aggregate at vmax and a rotating cluster on the wall with linear velocity vmax (red dashed line). The cumulative distribution of 51 relative aggregate velocities measured during the shaking speed ramp-down of phase 2 is shown (“+” signs). The Gaussian mean (location of the Gaussian distribution peak) and standard deviations of the measured relative velocities are at 1.84 |
In the text |
![]() |
Fig. 7 Retrieved number density of free-flying aggregates in the three cells of the SPACE experiment. The size distributions of large aggregates composed of polydisperse and monodisperse dust are shown as black solid and green dashed curves, respectively. The size distribution of small aggregates composed of monodisperse dust is shown as the red dash-dotted curve. The origin of time is at the start of data recording. The shaking cycles are indicated by dashed lines. |
In the text |
![]() |
Fig. 8 Cell wall surface area covered by clusters during the SPACE experiment for the three dust types. The large size distribution of aggregates composed of poly- and monodisperse dust are shown by the black solid and green dashed curve, respectively. The small size distribution of aggregates composed of monodisperse dust is shown by the red dash-dotted curve. The origin of time is at the start of data recording. The shaking cycles are indicated by dashed lines. |
In the text |
![]() |
Fig. 9 Normalized fragmentation rate of the clusters on the cell walls during the fragmentation phase of the second shaking cycle for the smaller and larger size distributions of aggregates composed of monodisperse SiO2 (red dash-dotted for small aggregates and green dashed line for large aggregates, respectively) and the larger size distribution of aggregates composed of polydisperse SiO2 (black solid line), plotted as a function of the centrifugal acceleration induced on the aggregates by the rotation of the cell wall. The curves are terminated at the acceleration where they reach the saturation of free-flying aggregates in the cell volume. The fragmentation rate maxima are shown by plus signs. |
In the text |
![]() |
Fig. 10 Monomer aggregate loss per cluster per frame (corresponding to 5.88 ms) for clusters on the cell walls during the fragmentation phase of the second shaking cycle for a) the small size distribution of aggregates composed of monodisperse SiO2; b) the large size distribution of aggregates composed of monodisperse SiO2; and c) the large size distribution of aggregates composed of polydisperse SiO2, respectively. |
In the text |
![]() |
Fig. 11 Typical wall cluster during the experiment from the large size distribution of aggregates composed of monodisperse dust. a) Original greyscale image. b) Contrast enhanced image. The red line in b) indicates the delimitation of the total cluster area. The packing density φp = 0.91 is determined by dividing the dark area by the area enclosed by the red line (see Sect. 4.2 for more details on the cluster morphology). |
In the text |
![]() |
Fig. 12 Sticking probability of the three types of dust aggregates investigated in the SPACE experiment plotted as a function of the statistical relative collision velocity computed in Sect. 3.2 (1.84vmax). The large size distribution of aggregates composed of poly- and monodisperse dust are shown by the black solid curve with squares and the green dashed curve with triangles, and the small size distribution of aggregates composed of monodisperse dust is shown by the red dash-dotted curve with diamonds. |
In the text |
![]() |
Fig. 13 Flight collision results (dotted boxes) according to the aggregate size distributions (box heights) and their mean collision velocities during the experiment (box widths). As the larger aggregates composed of mono- and polydisperse SiO2 grains had very similar sticking properties, only the results for the aggregates composed of monodisperse SiO2 are plotted (upper box). The background colours correspond to the dust collision model developed by Güttler et al. (2010) with updates by Kothe et al. (2013). Green colour represents sticking collisions, yellow bouncing, and red fragmentation, respectively. The lines represent the transition zones between sticking and bouncing (solid line for a 50% sticking probability, dashed lines for 0, 25, 75 and 100% sticking probabilities) as well as the transition between bouncing and fragmentation (dash-dotted line for the onset of fragmentation), computed by Kothe et al. (in prep.). Mind that the collision model by Güttler et al. (2010) and Kothe et al. (2013) was developed for collisions between equal-mass aggregates, whereas the results presented here are dominated by single aggregates colliding with more massive clusters. |
In the text |
![]() |
Fig. 14 Tensile strength of aggregates in clusters investigated in the SPACE experiment for aggregates composed of mono- and polydisperse dust (diamonds and triangle, respectively) as well as of mm-sized dust aggregates measured by Blum et al. (2014; squares) in comparison to the model prediction by Skorov & Blum (2012). The model is applied for aggregate filling factors of 0.3 (dashed line) and 0.37 (solid line), respectively. |
In the text |
![]() |
Fig. 15 Aggregate pull-off forces as a function of their size for the two types of dust investigated in the SPACE suborbital flight experiment, i.e. monodisperse (diamonds) and polydisperse SiO2 (triangle). Linear fits through the data and the origin are also shown, with slopes of 7.6 × 10-4 (dashed line) and 9.8 × 10-4 (dotted line) for the two dust types, respectively. The smaller aggregates composed of monodisperse dust are assumed to detach in clumps of four aggregates. |
In the text |
![]() |
Fig. 16 Relative velocities between dust aggregates in a protoplanetary disk computed according to a) Weidenschilling (1977b); b)Andrews & Williams (2007); and c)Desch (2007) at 1 AU and a turbulence parametre of α = 10-5. The velocity profiles are labelled in units of m s-1. The hashed boxes represent the collisions observed during the SPACE experiment (see text for more details). The dashed contour displays the minimum sticking velocity measured during the SPACE experiment for both size distributions. The dotted contour displays the same minimum sticking velocity in the case of aggregate-aggregate collisions instead of aggregate-cluster collisions (see text for details). |
In the text |
![]() |
Fig. A.1 Image processing steps. a. Recorded SPACE data frame at 11.8 s after beginning of frame recording, showing one of the smaller experiment cells. b. The same frame after correction of background illumination. c. The same frame after background correction and averaging over 201 frames. The clusters on the cell walls become apparent. |
In the text |
![]() |
Fig. A.2 Analysis method used on the SPACE experiment data to determine the number of free-flying aggregates in the experiment cells at each moment in time. a) Region of interest inside an averaged data frame of the larger experiment cell (see Fig. 1b) at the beginning of the experiment run (1.2 s after start of recording), when all aggregates were free-flying. Image processing as shown in Fig. 7c was applied. b) Region of interest inside an averaged data frame of the same experiment cell during the slowest shaking phase (47.1 s after start of recording), when all aggregates are incorporated in clusters sticking to the glass walls. Image processing as shown in Fig. 7c was applied. c) greyscale histograms of frames a) and b) When all aggregates are free-flying a), the background is uniformly grey (dashed red line), whereas the dark clusters (small peak at lower greyscale values) and the bright background (high and narrow peak at higher greyscale values) can clearly be distinguished when the aggregates are all incorporated in clusters (solid black line). |
In the text |
![]() |
Fig. A.3 Analysis of the clusters on the cell walls inside the SPACE experiment during the suborbital flight. a) Region of interest inside an averaged and binarized data frame of the bigger experiment cell at the beginning of a fragmentation phase (61.8 s after beginning of the recording). b) Region of interest inside an averaged and binarized data frame of the same cell at the end of this fragmentation phase (76.5 s after beginning of the recording). The area covered by clusters (black areas) has visibly decreased. |
In the text |
![]() |
Fig. B.1 Summary of the Monte Carlo simulation results for the agglomeration of the small monodisperse dust size distribution used in the SPACE experiment rocket flight (mean size of 120 μm). The simulation surface of 221 × 362 pixel represents a region of interest inside the bigger experiment cell. Each column has the same initial number of particles inserted into the simulation Nsim. Each row has the same glass sticking probability κ. |
In the text |
![]() |
Fig. B.2 Monte Carlo simulation results for the agglomeration of the large size distribution of aggregates composed of monodisperse dust used in the SPACE experiment rocket flight (mean size of 330 μm). The simulation surface of 192 × 141 pixel represents a region of interest inside one of the smaller experiment cells. Each column has the same initial number of particles inserted into the simulation Nsim. Each row has the same glass sticking probability κ. |
In the text |
![]() |
Fig. B.3 Results of a Monte Carlo simulation for the smaller aggregate size distribution used in SPACE during the suborbital flight: a) binarized SPACE data frame at 47 s after start of data recording. At this point of the experiment run, all dust aggregates are clustered on the cell wall. b) Monte Carlo simulation depositing 4600 particles on a surface with a glass sticking parametre κ = 0.0055, c) size distribution of the data frame seen in a. (light red line, triangles) and 5000 runs of Monte Carlo simulations with the parametres of image b) (grey lines). The bin size for these distributions is 4 pixel (1.23 × 10-8 m2). |
In the text |
![]() |
Fig. C.1 Reduced χ2 error statistics for linear and parabolic fits to the trajectories of 10 trackable aggregates in the cell containing the larger aggregates composed of monodisperse dust during the SPACE rocket flight experiment run, before (triangles, open for the linear fit and filled for the parabolic fit) and after impact with a glass wall (squares, open for the linear fit and filled for the parabolic fit). |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.