The cores in molecular clouds are the densest and coldest regions of the interstellar medium (ISM). In these regions ISM-dust grains have the potential to coagulate. This study investigates the collisional evolution of the dust population by combining two models: a binary model that simulates the collision between two aggregates and a coagulation model that computes the dust size distribution with time. In the first, results from a parameter study quantify the outcome of the collision -sticking, fragmentation (shattering, breakage, and erosion) -and the effects on the internal structure of the particles in tabular format. These tables are then used as input for the dust evolution model, which is applied to an homogeneous and static cloud of temperature 10 K and gas densities between 10 3 and 10 7 cm −3 . The coagulation is followed locally on timescales of ∼10 7 yr. We find that the growth can be divided into two stages: a growth dominated phase and a fragmentation dominated phase. Initially, the mass distribution is relatively narrow and shifts to larger sizes with time. At a certain point, dependent on the material properties of the grains as well as on the gas density, collision velocities will become sufficiently energetic to fragment particles, halting the growth and replenishing particles of lower mass. Eventually, a steady state is reached, where the mass distribution is characterized by a mass spectrum of approximately equal amount of mass per logarithmic size bin. The amount of growth that is achieved depends on the cloud's lifetime. If clouds exist on free-fall timescales the effects of coagulation on the dust size distribution are very minor. On the other hand, if clouds have long-term support mechanisms, the impact of coagulation is important, resulting in a significant decrease of the opacity on timescales longer than the initial collision timescale between big grains.

The dust size distribution in molecular clouds can be strongly affected by ice-mantle formation and (subsequent) grain coagulation. Following previous work where the dust size distribution has been calculated from a state-of-the art collision model for dust aggregates that involves both coagulation and fragmentation (Paper I), the corresponding opacities are presented in this study. The opacities are calculated by applying the effective medium theory assuming that the dust aggregates are a mix of 0.1 μm silicate and graphite grains and vacuum. In particular, we explore how the coagulation affects the near-IR opacities and the opacity in the 9.7 μm silicate feature. We find that as dust aggregates grow to μm-sizes both the near-IR color excess and the opacity in the 9.7 μm feature increases. Despite their coagulation, porous aggregates help to prolong the presence of the 9.7 μm feature. We find that the ratio between the opacity in the silicate feature and the near-IR color excess becomes lower with respect to the ISM, in accordance with many observations of dark clouds. However, this trend is primarily a result of ice mantle formation and the mixed material composition of the aggregates, rather than being driven by coagulation. With stronger growth, when most of the dust mass resides in particles of size ∼10 μm or larger, both the near-IR color excess and the 9.7 μm silicate feature significantly diminish. Observations at additional wavelengths, in particular in the sub-mm range, are essential to provide quantitative constraints on the dust size distribution within dense cores. Our results indicate that the sub-mm index β will increase appreciably, if aggregates grow to ∼100 μm in size.

No abstract

The coagulation of dust aggregates occurs in various astrophysical environments. Each one is characterized by different conditions that influence the growth, e.g., relative velocities, composition, and size of the smallest constituents (monomers). Here we study the microphysics of collisions of dust aggregates in a four-dimensional parameter space. The parameters are the collision energy, the initial compactness of agglomerates, the mass ratio of collision partners, and the impact parameter. For this purpose we employ a state of the art molecular dynamics type of model that has been extensively and successfully tested against laboratory experiments. It simulates the motion of individual monomers interacting dynamically via van der Waals surface forces. The structure of aggregates is quantified by the filling factor that provides information about the internal structure, the packing density of monomers, and the projected surface area of aggregates. Our results show the importance of the impact parameter that causes formation of elongated particles, due to tensile forces acting in offset collisions. In head-on impacts, aggregates are compacted at lower energies. A sufficiently high energy causes restructuring to reach maximum compaction. If more energy is provided, pancake-like structures are formed. We find that the outcome of collisions can be represented in a simple way. A highly pronounced large fragment component coexists with a power-law distribution of small fragments. The structural parameter of these small fragments is described very well by a simple relation, largely independent of the initial compactness, impact energy or impact parameter. The simulations show that erosion by collisions with high mass-ratio can be significant. The ejected mass can be several orders of magnitude higher than the impactor mass. This contrasts with collisions of equal mass aggregates, where the same impact energy can lead to perfect sticking. These findings are summarized in the form of a simple collision recipe. The recipe specifies the outcome of a collision, averaged over the impact parameter. It is provided in tabular form for a range of physical parameters such as impact energy and pre-collision filling factor. The dependence on the mass ratio of impactor and target is taken into account by providing both a local and a global branch of the recipe.