On the scaling of fragmentation and energy dissipation in collisions of dust aggregates

Fragmentation of granular clusters may be studied by experiments and by granular mechanics simulation. When comparing results, it is often assumed that results can be compared when scaled to the same value of E/Esep\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E/E_{\mathrm{sep}}$$\end{document}, where E denotes the collision energy and Esep\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_{\mathrm{sep}}$$\end{document} is the energy needed to break every contact in the granular clusters. The ratio E/Esep∝v2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E/E_{\mathrm{sep}}\propto v^2$$\end{document} depends on the collision velocity v but not on the number of grains per cluster, N. We test this hypothesis using granular-mechanics simulations on silica clusters containing a few thousand grains in the velocity range where fragmentation starts. We find that a good parameter to compare different systems is given by E/(NαEsep)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E/(N^{\alpha }E_{\mathrm{sep}})$$\end{document}, where α∼2/3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha \sim 2/3$$\end{document}. The occurrence of the extra factor Nα\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N^{\alpha }$$\end{document} is caused by energy dissipation during the collision such that large clusters request a higher impact energy for reaching the same level of fragmentation than small clusters. Energy is dissipated during the collision mainly by normal and tangential (sliding) forces between grains. For large values of the viscoelastic friction parameter, we find smaller cluster fragmentation, since fragment velocities are smaller and allow for fragment recombination.


Introduction
The understanding of the dynamics of dust particles and of agglomerates that consist of dust grains is of fundamental importance in the field of planet formation [1]. The earliest stage of the formation process is believed to be the agglomeration of dust grains in accretion disks around protostars as already discussed by Weidenschilling [2] who argues that relative velocities in such disks are sufficiently high to enable collisions, as well as by Weidenschilling and Cuzzi [3] who discuss turbulence as the driving force for particle collisions.
In such studies, the question whether agglomerates will fragment during collisions or stick and enable grain growth is of prime importance. This issue has been addressed by a variety of experimental [4][5][6][7][8][9] and theoretical/computational approaches [10][11][12][13][14], as reviewed in Refs. [15,16]. The results of such studies can be used to provide basic information for use in models that describe the temporal evolution of the agglomerate size distribution in dust disks [17,18].
Such models need a description of the fragmentation outcome of collisions spanning a wide range of cluster sizes and velocities, and for various relevant materials such as silica and ice. Also the comparison of experimental and computational data on cluster collisions-which have often been obtained in different cluster size regimes-needs a common base for discussion. Already early, Dominik and Tielens [10] argued that the prime quantity for cluster fragmentation in collisions is given by the ratio of the collision energy E and the separation energy E sep , which is the energy needed to break every contact in the aggregate; as soon as E exceeds E sep by a factor of around 10, catastrophic disruption of the cluster was predicted. This argument was later used repeatedly in order to discuss experimental and computational results of cluster collisions [4,13,19].
In the present paper, we analyze the validity of the argument that E∕E sep describes the fragmentation behavior in cluster collisions using granular mechanics simulations. To this end we vary the collision energy systematically by changing either the cluster size or the cluster velocity; we also change the energy dissipation during the collision to 1 3 33 Page 2 of 13 obtain a greater variety of collision results under otherwise identical conditions. The results will show that E∕E sep is not a good scaling parameter to describe fragmentation. However, we will demonstrate that E∕(N E sep ) , where N is the number of grains per cluster and ∼ 2∕3 , describes fragmentation well, and will discuss this finding.

Granular mechanics
Granular mechanics is the theory describing interactions of grains of a finite size. A theoretical description of grain interactions needs to know whether or not two grains are in contact and can then derive forces and torques based on the relative motion of the grains.
We use spherical grains, where each grain has the same radius, r, described by the materials parameters M and . Here M = Y∕(1 − ) 2 is the indentation modulus with Young's modulus Y and Poisson number , and is the specific surface energy. The model used in our work has already been used previously [19][20][21][22], and we give a short summary here for convenience.
Let us denote grain positions by ⃗ x , velocities by ⃗ v and angular velocities by ⃗ . Indices i and j refer to the respective grains. Two grains i and j interact with each other if their overlap is greater than 0. The contact area of two spherical grains is a circle of radius The relative motion can be split into normal and tangential motion as follows. With the normal unit vector the contact velocity is defined as The normal component of the velocity, ⃗ v n , is obtained by projecting ⃗ v c on n , and the tangential contribution of the velocity is then given by The unit tangent vector is defined as The normal force is given as a sum of the adhesive force and a viscoelastic term, and vanishes for negative overlap . The parameter A in the viscoelastic term controls the amount of friction in normal motion and will be denoted as friction parameter.
In a similar way, the relative angular velocity ⃗ = ⃗ i − ⃗ j can be split up into two contributions, twisting (twist), and rolling (roll), Tangential forces and torques for sliding (slid), rolling and twisting motion can be expressed as respectively, where the hat ( ̂ ) indicates a unit vector. yield denotes the critical distance over which two spheres can roll without breaking bonds.
Further details of the algorithm are provided in Ref. [20].

Materials parameters
We use silica grains with a set of standard parameters [20,23,24] as silica is often used in laboratory experiments as an analogue of interplanetary grain material [16,[25][26][27]. The grains have a radius of r = 0.76 μm and a mass density of = 2000 kg∕m 3 ; the mass of a grain is thus m = 3.68 ⋅ 10 −15 kg . The elastic properties are given by M = Y∕(1 − 2 ) = 55.6 GPa with Young's modulus Y = 54 GPa and Poisson number = 0.17 . The specific surface energy is assumed as = 0.05 J∕m 2 . We recently noted that our implementation of the granular-mechanics scheme uses only half of Young's modulus. The parameter A which controls the friction in normal motion has a standard value of A orig = 0.5 ns , as determined from experimental data [20]. We allow it here to vary between 0.5A orig to 1.00A orig in steps of 0.25A orig . This variation was done in order to allow for some variation in the fragmentation results for identical clusters (structure, v, N), as well as to study the influence of energy dissipation on the fragmentation behavior. The rolling-friction parameter yield is set to 3.2 nm as recommended by Heim et al. [25]. The collision energy E, i.e., the kinetic energy at the start of the simulation, where 2 clusters are shot towards each other with a relative velocity v, can be expressed as and is thus seen to be proportional to N and v 2 . For a velocity For studying fragmentation, the so-called break-up energy E break is introduced, which represents the energy necessary to break a contact. It is given by the work required to move two grains from their equilibrium overlap (i.e., the overlap, where adhesion and elastic repulsion cancel) to an overlap of 0; it is [28] Here K is a numerical constant equal to K = (12∕5) (3∕2) 2∕3 5∕3 = 21.19 [28]. In our case, it is E break = 68.46 aJ. Note that normal friction is neglected in the calculation of E break , Eq. (14), and will lead to an increase in this quantity. While the dependence of E break on the materials parameters is common knowledge [10], the prefactor K has been intensely discussed. Dominik and Tielens [10] advocate an increased value in order to take into account the energy dissipation by vibrations induced by the collision. From a consideration of experimental data [4], Ormel et al. [13] even argue for a value of K = 1760 . However, the magnitude of K does not influence the scaling behavior of fragmentation, which is the topic of the present study.
In addition, we define the total separation energy of a cluster (break-up of all contacts) as where N c denotes the number of contacts in each of the colliding clusters.

Simulations and their analysis
We build spherical granular clusters of radius R containing a specified number of grains, N, with a specified filling factor . In this work, we use = 0.28 . For the values of N = 1250 , 2500 and 3750 grains, the clusters thus have radii of R = 12.51 , 15.77 and 18.05 μm , respectively.
The clusters are constructed by an algorithm based on Ref. [21], which allows to produce clusters with pre-defined filling factor inside a spherical region. In short, new grains are added to the growing cluster at the positions of the smallest local filling factor until the prescribed number of grains have been put into the sphere. The local filling factor is calculated in a sphere of radius 3r around each grain. If a new grain has been rejected 9 times-grains are rejected if they overlap with other grains-it is attached to a randomly chosen grain; such situations may occur for large filling factors. Finally, the cluster is allowed to relax.
Since the grains are placed step by step by attaching one grain, the expected number of contacts is N c = N − 1 . In fact, we find 1251, 2506 and 3754 contacts, respectively. Thus, it holds to a good approximation N c = N . The separation energies thus amount to 85.6, 171.6 and 257.0 fJ for clusters with N = 1250 , 2500 and 3750 grains.
We study the collision of two clusters with relative velocity v in the range of 10-30 m/s; the velocity is varied in steps of 5 m/s. To this purpose, we duplicate the cluster and rotate it randomly. We perform only central collisions, i.e., with a vanishing impact parameter. In order to allow for statistical reliability of the results, each simulation is repeated 10 times; in each simulation, the relative orientations of the two colliding clusters are changed randomly. The results presented below are all averages over these 10 simulations. The setup is shown schematically in Fig. 1, while Fig. 2

displays snapshots of an actual collision.
This work uses the granular mechanics code LIGGGHTS [29], which implements the forces described in Sect. 2.1 above. Every simulation is performed with 1.2 ⋅ 10 6 time steps of 50 ps, leading to a total duration of 60 μs.
To analyze the collision results, we use the algorithm of cluster analysis proposed by Stoddard [30], which identifies the fragments generated by the collision. Let us denote by N 1 and N 2 the numbers of grains in the two largest fragments. We calculate the agglomeration parameter [31], and the fragmentation parameter [31],

Fig. 1 Schematic setup of a cluster-cluster collision
It is possible to calculate the grain density by using the methods of smoothed-particle dynamics [32,33]. For any point ⃗ r in space, the grain density can be calculated as Here the Gaussian serves as a a smoothing kernel. The matrix Σ is a 3 × 3 diagonal matrix, 4r 2 . Thus each grain is smeared over a width of 2r. Similarly, the density of grain contacts can be calculated by substituting ⃗ x i in Eq. (18) by the contact points, , if i and j are in contact, and summing over all contacts. The kinetic energy density is also given in analogy to Eq. (18), if the smoothing kernel is multiplied by the kinetic energy, E i , of grain i. We determine E i from both the translational and rotational energy of grain i; however, the rotational energy is always at least a factor of 10 −3 smaller than the translational energy. The kinetic energy of grain i is determined relative to the center of mass of its home cluster.

Representative case: time evolution
We show the time evolution of a prototypical collision-N = 3750 , A = A orig , v = 10 m∕s-in the snapshots provided in Figs. 3, 4 and 5. In these snapshots, the densities-as calculated by Eq. (18) and its analogs-are plotted in a plane containing the collision axis.
In the following, we denote as the contact plane the midplane between the two clusters, perpendicular to the collision axis. Here, the two clusters meet and the most energetic processes are going to occur in its neighborhood, which we denote as the collision zone, see Fig. 5a.
The grain density, Fig. 3, shows how during the collision the grain density in the collision zone strongly increases-approximately by a factor of 2. At around 3 μs , the shape of the merged cluster is almost spherical, but due to grain compression, its volume has considerably shrunk along the collision axis. In lateral direction, however, the merged cluster expanded slightly. The final . shape at the end of the collision is oblate, as the merged cluster has shrunk considerably along the collision axis. The density of contacts, Fig. 4, shows similar features. After 1 μs , many contacts in the central collision region are broken. Already at 2 μs , the destroyed contacts in the collision zone have not only been reformed but new contacts have been generated. Note that at this time, the highly coordinated core region is surrounded by a shell of low contact density, which indicates how the collision energy spreads out of the primary collision zone and destroys other grain contacts. At 3 μs , the contact density in the collision zone is even higher than the grain density, while initially, it had the same magnitude. This shows that during grain compaction, the number of contacts per grain, i.e., the coordination number, has increased above the value of 2 characteristic for clusters assembled by ballistic aggregation [14]. This effect has spread over the entire merged cluster at the end of the simulation.
The energy density, Fig. 5, shows most clearly how the energy slowly diffuses out of the contact region, see Fig. 5a at 1 μs , to the entire merged cluster. At 3 μs , the kinetic energy is only present in a bulge surrounding the merged cluster at the periphery of the collision zone, while all other energy has been dissipated away. We do not show the energy density at the end of the simulation, since it is close to zero everywhere.
Already in this simulation, around 52% of the grains were ejected from the merged cluster; this fraction keeps increasing with collision velocity. At 30 m/s, only 3.4% of the grains are retained in the largest fragment.
We conclude that these snapshots show that the collision is terminated within approximately 3 μs . The clusters do not interpenetrate each other at this filling factor, 0.28, but are stopped at the surface such that a thin collision zone develops, in which contacts are destroyed. The energy spreads out from this zone into the remaining parts of the clusters but has dissipated away at the end of the simulation. The collision zone has been compacted and grain coordination has been increased.
The results for other simulations look similar. The merged cluster becomes flattened along the collision axis for higher velocities such that it assumes a pancake shape, and more grains are ejected from it, see also Fig. 2. Further grain-based snapshots of such events can be seen in Figs. 7 and 9 of Ref. [19].

Scaling
The agglomeration parameter X, Eq. (16), is large if the largest fragment-i.e., the merged cluster-contains almost all the mass of the system. It is often assumed [19,31] that small small values X < 0.15 identify the fragmentation regime. On the other hand, the fragmentation parameter,  17), assumes large values when the mass contained in the two largest fragments is only a small fraction of the total mass. Hence these two parameters are often used to identify whether cluster collisions are in the sticking or in the fragmentation regime. We present these parameters in Fig. 6a and b for the set of simulations performed: velocities in the range of v = 10 -30 m/s, N = 1250-3750 grains per cluster, and friction coefficients varying between A∕A orig = 0.5-1.00. The figure shows the well known facts that with increasing velocity and decreasing grain numbers, more fragmentation occurs [19,22,34]. In addition, we observe that with increasing dissipation (here: normal friction) fragmentation is prevented. This appears to be a novel result; it was reached here because we use exactly the same clusters at the same velocities but with changed values of the energy dissipation. It appears plausible since dissipation destroys the kinetic energy that is available for grain bouncing and hence for cluster fragmentation. The effect is astonishingly strong; by reducing the friction coefficient by a factor of 0.5, nearly all collisions in the velocity and size window studied here become fragmenting. A reduction of A by a factor of 2 has a similar-but somewhat more pronounced-effect as a reduction of the collision energy by a factor of 2. This is also the reason why we did not change the friction coefficient in a wider range.
In Fig. 6, we used the scaling parameter E∕E sep , which compares the available collision energy to the total energy needed to break all intergranular contacts, E sep = NE break , Eq. (15). For fixed N and A, the energy is varied by the cluster velocity, v. Such a scaling has been proposed by Dominik and Tielens [10] and has been used repeatedly to compare fragmentation results among different clusters [4,13,19]. Figure 6 shows that this parameter does not allow to collapse all collision data on a single curve. However, this is possible, when changing the scaling parameter to explicitly contain the number of grains, N. Figure 7 shows that our data collapse on a single curve for all N, if we use the scaling parameter E∕(N E sep ) ; the best results were achieved with in the range 0.6-0.7. We will use = 2∕3 from now on. We note that for our clusters, E∕E sep ∝ v 2 , since the grain number drops out from the ratio, and the materials parameters appearing in E break are unchanged. A scaling parameter E∕(N E sep ) with ≠ 0 thus means that cluster size explicitly enters the description. Data collapse works better with the fragmentation parameter N s than with the agglomeration parameter X. This finding coincides with our experience that X has some ambiguity in defining fragmentation: small values of X might either mean strong fragmentation or intact bouncing of the clusters. We therefore propose that the use of N s is more adequate in quantifying cluster fragmentation.
Of course, for different dissipation constants A, different master curves are obtained, since the details of the collision and of the fragmentation change, see Sect. 3.4. It is possible to collapse all data to just one master curve by using a parameter E∕[N E sep (A∕A orig ] ) ; Fig. 8 demonstrates that such a data collapse works well with = 0.8 . The dependence on A means that for lower friction, the cluster will exhibit stronger fragmentation, which is expected.
Since the number of grains scales like the cluster volume, N ∝ R 3 , the value of = 2∕3 can be associated with a cluster-size dependence corresponding to that of the cluster cross section. A value of = 2∕3 can, however, also be discussed in terms of the cluster passage time since then the scaling parameter E∕(N 2∕3 E sep ) ∝ (v∕R) 2 = t −2 pass . Our result that the relevant scaling parameter for cluster fragmentation in collisions, E∕(N E sep ) with = 2∕3 , thus points at the fact that during cluster stopping in the collision, the collision energy cannot be used directly to break all clusters. Fragmentation first occurs in the shallow collision zone with an area ∝ N 2∕3 and then expands into the entire merged clusters (volume ∝ N ). These features explain the exponent of = 2∕3 . Evidently, energy dissipation in the collision dictates to what extent the entire collision energy is available to break all bonds, or is dissipated away, before it can reach each grain. For this reason, we now discuss the role of energy dissipation in the collisions studied here. Figure 9a and b show the fragmentation parameter plotted as a function of time for two collisions of clusters of N = 1250 grains each at collision velocity v = 20 m/s, with friction parameters A = 0.5A orig and A = A orig . We observe that for both collisions, N s goes up to almost 1 at times of 1.5-2.5 μs , indicating complete fragmentation. However, at later times, N s decreases again; this means that grains assemble again to form clusters. This process is particularly pronounced for the higher value of the friction parameter, Fig. 9b. The ensuing decrease in N s shows that fragmentation is not a simple process of breaking up several contacts just once. In fact, contacts can break, but contacts can also form (temporarily as well as permanently). In this specific case, contacts are more likely to form again for higher friction, since higher friction dissipates more kinetic energy and grains collide with lower relative velocities allowing for more new contacts to form. This is not only influenced by friction, but by the initial collision velocity as well, see Fig. 9c, since at lower collision velocities, contacts can form again.

Energy dissipation
For studying energy dissipation, we determine the time dependence of the power lost in the collision by all dissipation channels (sliding, twisting, rolling, normal motion). The power lost by dissipative forces is given by and the power lost by dissipative torques is given by (19) t pass = R∕v, where we sum over forces and torques exerted by grain j on grain i. In Fig. 10 we plot the time evolution of these powers (smoothed with a Gaussian filter of 15 ns width) for collisions of clusters with N = 1250 grains per cluster at a collision velocity of v = 20 m/s for two friction parameters, A = 0.5A orig and A = A orig . The figure shows that energy dissipation by forces (denoted by P f ) is dominant over energy dissipation by torques ( P t ). Energy dissipation due to normal-force dissipation ( P fn ) and sliding-force dissipation ( P fs ) is of similar order of magnitude. Since initially, the grains do not rotate, the collision excites grain rotation as can be seen by the positive sign in the sliding ( P ts ) rotational mode. The twisting mode ( P tt ) as well as (to a lesser extent) the rolling mode ( P tr ) are de-excited.  Table 1 assembles the energy balance at a time of 2 μs quantitatively. As Fig. 10 showed, dissipation data at 2 μs after the start of the collision is close to final. Table 1a quantifies that the energies dissipated by normal and sliding forces are dominant, while energies in rotational motion are minor; sliding induces rotation, and so this channel carries another sign than the other channels.
We compare the collision energy E = 459.7 pJ in Table 1b with the final energies. Around 20% of the initial energy are found at 2 s in translational motion of grains while their rotation carries only a minor fraction of energy. If we sum these kinetic energies with the dissipated energy, we obtain a value which agrees within better than 1% with the initial collision energy E. The small error is caused by the numerical integration of the power over time. The energy needed for the fragmentation of all grains, E sep , amounts to only 85.6 fJ in this case, and can be ignored in this calculation.
Fragment kinetic energies are higher for A = 0.5A orig , since less energy has been dissipated, while a large number of fragments has been created which then leave the collision zone without further interaction.
One can see that for the original value of A = A orig , the powers dissipated through normal and sliding motion are very close to each other, while the importance of normal motion for the dissipation decreases for lower A. This is clear from the fact that A enters into the normal force as a parameter.
A reduction of A by 50% does not lead to a similar reduction in the power dissipated in the normal-force channel ( P fn ), which is only reduced by 24%. In addition, this reduction is partly compensated by an increase in the power dissipated by sliding motion; this increase may be induced by the longer-lived contacts (reduced fragmentation, see Fig. 9) for large A. As a summary, we see that a reduction of A slightly (by 6%) reduces the energy dissipation during the collision, enhancing the kinetic energy going into the fragments. Figure 10 shows that the time scale for energy dissipation is around 2 μs . We may compare this time scale with the cluster passage time defined in Eq. (19). For a velocity of v = 20 m/s, t pass = 0.63 (0.79, 0.9) s for cluster sizes of N = 1250 (2500, 3750) grains. A comparison with Fig. 10 shows that the collisions are indeed over after around (2-3) t pass . Such short time scales for collisions are in agreement with previous results of cluster stopping in a granular environment [35], and can be justified by referring to energyproportional stopping as in Poncelet's law [36,37]. Figure 11 shows plots of the time evolution of the kinetic energy during the collision, averaged over 10 collisions. In Fig. 11a, we observe that for higher collision velocities, energy will initially be dissipated faster-roughly in agreement with the passage-time scale, Eq. (19)-but collisions at lower velocities dissipate energy for a longer time and, most strikingly, at the end of the simulation have dissipated a higher fraction of their starting energy. Figure 11b shows a similar effect for varied cluster size, N, such that largecluster collisions lead to longer-lasting energy dissipation. We argue that this feature is consistent with our observation that E∕E sep is not a good parameter to describe the collisions. As the number of grains and contacts effectively cancel out of E∕E sep , this parameter does not show any dependence on cluster size. It would describe a fragmentation process well, where impact energy is distributed through the cluster instantaneously. Realistic collisions, however, will initially store their energy in the collision zone and from there energy will gradually spread throughout the cluster. If contacts break too rapidly, there is no effective means for transporting the impact energy throughout the entire clusters, and contacts far away from the impact site are not affected by the collision. Collisions at lower velocities store less energy in the collision zone, break up less contacts and allow for energy to be distributed better throughout the cluster. This allows for more energy to be dissipated eventually.
The time evolution of the energy dissipation of colliding clusters is only negligibly affected by the dissipation constant A, Fig. 11c. This demonstrates that our choice of varying A to allow for a larger variety of collision situations affects the collision time scale only negligibly. However, energy dissipation proceeds for longer time scales for large values of A. This is caused by less fragmentation Table 1 Energy balance of the collision analyzed in Fig. 11 for two values of the friction constant A Top: Energies dissipated within 2 μs in the various channels, and their sum. Bottom: Comparison of the collision energy E with the kinetic and rotational energies of grains, E kin and E rot , respectively, at 2 μs , the total dissipated energy, E diss , which is given by the sum value of (a), and the final energy, occurring in this case, as discussed in connection with Fig. 9. Note that even though the energy dissipation processes do not appear to have completely stopped at the end of the simulation, the fragmentation patterns are final, see Fig. 9.
Minor relative motion in the grains still exists, but contacts do no longer change significantly.

Conclusions
We discuss the scaling of fragmentation in granular cluster collisions using as a basis a series of simulations for granular silica clusters, = 0.28 , containing several thousand grains at velocities where fragmentation starts. We obtain the following findings. Cluster fragmentation is accompanied by energy dissipation.
Here the potential energy spent in cluster break-up plays only a minor role in the examples studied here, since the parameter E∕E sep , which describes the ratio of the available collision energy to the energy needed to fragment the cluster completely, is of the order of 1000. We find that the collision energy is spent mostly by normal and tangential (sliding) forces between grains, while torques mediated by grain rotation are negligible in this respect. In order to understand the effect of energy dissipation on cluster fragmentation in more detail, we allowed the viscoelastic friction parameter A, which describes energy dissipation of the normal contact force, to vary. Besides the obvious effect that lower values of A lead to less energy dissipation and hence a higher remaining kinetic energy of the fragments, we also find that it induces higher fragmentation. This effect is counter-intuitive at first sight, but is caused by the fact that high energy dissipation leads to smaller velocities of the fragments that thus have a higher chance to recombine to clusters. Indeed we see that for high A, fragmentation has a non-monotonic time evolution, and thus features a prominent grain recombination regime. As a consequence, dissipation also lasts for a longer time than for small A where the cluster is more completely fragmented.
The time scale of energy dissipation amounts to around three times the cluster passage time, Eq. (19). However, changes in the fragmentation parameters can proceed for a longer time, in particular for large clusters and low velocities, as at low energies new contacts are able to form.
In our systems, the difference between the number of grains, N, and the number of contacts, N c , is negligible, since for clusters assembled by ballistic deposition, the local coordination of each grain is 2-except at the surface-such there is only a constant factor between N and N c , which is of no consequence for scaling arguments. Also for more compact clusters, where the coordination number is larger, the difference between N and N c will only introduce a constant factor in the scaling parameter.
One might speculate to what extent the cluster shape influences collisions and the resulting fragmentation patterns. Such an influence is indeed to be expected: If two oblate clusters meet with their flat sides, more fragmentation is to be expected than for two prolate clusters colliding with their tips. This argument follows the dependence of our scaling on the cluster cross section via N . For the dissipation of non-collisional energy input-such as by granular temperature or release of elastic energy-such an argument is not so easily available, as such an energy input might either be homogeneous in the entire cluster or follow a law not determined by the cluster cross section. We therefore do not expect our analysis to be easily extended to non-collisional energy transfers.
The fragmentation behavior during cluster collisions is relevant for instance in the field of dust agglomeration in protoplanetary disks, and our material parameters were chosen to apply for such situations. Here, scaling parameters are used for comparing experimental and simulational data, or for extrapolating available data to other size or velocity regimes. We expect that our argument that E∕(N E sep ) rather than E∕E sep is the appropriate scaling parameter for cluster fragmentation will prove useful in data comparison and extrapolation as well as in numerical schemes for calculating the temporal evolution of dust disks.