No Arabic abstract
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{mu}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{mu}m silicate feature. We find that as dust aggregates grow to {mu}m-sizes both the near-IR color excess and the opacity in the 9.7 {mu}m feature increases. Despite their coagulation, porous aggregates help to prolong the presence of the 9.7{mu}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{mu}m or larger, both the near-IR color excess and the 9.7{mu}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 {beta} will increase appreciably, if aggregates grow to ~100{mu}m in size.
Context. Grains in circumstellar disks are believed to grow by mutual collisions and subsequent sticking due to surface forces. Results of many fields of research involving circumstellar disks, such as radiative transfer calculations, disk chemistry, magneto-hydrodynamic simulations largely depend on the unknown grain size distribution. Aims. As detailed calculations of grain growth and fragmentation are both numerically challenging and computationally expensive, we aim to find simple recipes and analytical solutions for the grain size distribution in circumstellar disks for a scenario in which grain growth is limited by fragmentation and radial drift can be neglected. Methods. We generalize previous analytical work on self-similar steady-state grain distributions. Numerical simulations are carried out to identify under which conditions the grain size distributions can be understood in terms of a combination of power-law distributions. A physically motivated fitting formula for grain size distributions is derived using our analytical predictions and numerical simulations. Results. We find good agreement between analytical results and numerical solutions of the Smoluchowski equation for simple shapes of the kernel function. The results for more complicated and realistic cases can be fitted with a physically motivated black box recipe presented in this paper. Our results show that the shape of the dust distribution is mostly dominated by the gas surface density (not the dust-to-gas ratio), the turbulence strength and the temperature and does not obey an MRN type distribution.
Dust and gas energetics are incorporated into a cluster-scale simulation of star formation in order to study the effect of heating and cooling on the star formation process. We build on our previous work by calculating separately the dust and gas temperatures. The dust temperature is set by radiative equilibrium between heating by embedded stars and radiation from dust. The gas temperature is determined using an energy-rate balance algorithm which includes molecular cooling, dust-gas collisional energy transfer, and cosmic-ray ionization. The fragmentation proceeds roughly similarly to simulations in which the gas temperature is set to the dust temperature, but there are differences. The structure of regions around sink particles have properties similar to those of Class 0 objects, but the infall speeds and mass accretion rates were, on average, higher than those seen for regions forming only low-mass stars. The gas and dust temperature have complex distributions not well modeled by approximations that ignore the detailed thermal physics. There is no simple relationship between density and kinetic temperature. In particular, high density regions have a large range of temperatures, determined by their location relative to heating sources. The total luminosity underestimates the star formation rate at these early stages, before ionizing sources are included, by an order of magnitude. As predicted in our previous work, a larger number of intermediate mass objects form when improved thermal physics is included, but the resulting IMF still has too few low mass stars. However, if we consider recent evidence on core-to-star efficiencies, the match to the IMF is improved.
Density profiles of isolated cores derived from thermal dust continuum emission rely on models of dust properties, such as mass opacity, which are poorly constrained. With complementary measures from near-infrared extinction maps, we can assess the reliability of commonly-used dust models. In this work, we compare Herschel-derived maps of the optical depth with equivalent maps derived from CFHT WIRCAM near-infrared observations for three isolated cores: CB68, L429, and L1552. We assess the dust opacities provided from four models: OH1a, OH5a, Orm1, and Orm4. Although the consistency of the models differs between the three sources, the results suggest that the optical properties of dust in the envelopes of the cores are best described by either silicate and bare graphite grains (e.g., Orm1) or carbonaceous grains with some coagulation and either thin or no ice mantles (e.g., OH5a). None of the models, however, individually produced the most consistent optical depth maps for every source. The results suggest that either the dust in the cores is not well described by any one dust property model, the application of the dust models cannot be extended beyond the very center of the cores, or more complex SED fitting functions are necessary.
It has recently been shown that turbulence in the interstellar medium (ISM) can significantly accelerate the growth of dust grains by accretion of molecules, but the turbulent gas-density distribution also plays a crucial role in shaping the grain-size distribution. The growth velocity, i.e., the rate of change of the mean grain radius, is proportional to the local gas density if the growth species (molecules) are well-mixed in the gas. As a consequence, grain growth happens at vastly different rates in different locations, since the gas-density distribution of the ISM shows a considerable variance. Here, it is shown that grain-size distribution (GSD) rapidly becomes a reflection of the gas-density distribution, irrespective of the shape of the initial GSD. This result is obtained by modelling ISM turbulence as a Markov process, which in the special case of an Ornstein-Uhlenbeck process leads to a lognormal gas-density distribution, consistent with numerical simulations of isothermal compressible turbulence. This yields an approximately lognormal GSD; the sizes of dust grains in cold ISM clouds may thus not follow the commonly adopted power-law GSD with index -3.5, but corroborates the use of a log-nomral GSD for large grains, suggested by several studies. It is also concluded that the very wide range of gas densities obtained in the high Mach-number turbulence of molecular clouds must allow formation of a tail of very large grains reaching radii of several microns.
AGB stars are, together with supernovae, the main contributors of stellar dust to the interstellar medium (ISM). Dust grains formed by AGB stars are thought to be large. However, as dust nucleation and growth within their outflows are still not understood, the dust-grain size distribution (GSD) is unknown. This is an important uncertainty regarding our knowledge of the chemical and physical history of interstellar dust, as AGB dust forms $sim$ 70% of the starting point of its evolution. We expand on our chemical kinetics model, which uniquely includes a comprehensive dust-gas chemistry. The GSD is now allowed to deviate from the commonly assumed canonical Mathis et al. (1977) distribution. We find that the specific GSD can significantly influence the dust-gas chemistry within the outflow. Our results show that the level of depletion of gas-phase species depends on the average grain surface area of the GSD. Gas-phase abundance profiles and their possible depletions can be retrieved from observations of molecular emission lines when using a range of transitions. Due to degeneracies within the prescription of GSD, specific parameters cannot be retrieved, only (a lower limit to) the average grain surface area. Nonetheless, this can discriminate between dust composed of predominantly large or small grains. We show that when combined with other observables such as the spectral energy distribution and polarised light, depletion levels from molecular gas-phase abundance profiles can constrain the elusive GSD of the dust delivered to the ISM by AGB outflows.