No Arabic abstract
Dust is formed out of stellar material and is constantly affected by different mechanisms occurring in the ISM. Dust grains behave differently under these mechanisms depending on their sizes, and therefore the dust grain size distribution also evolves as part of the dust evolution itself. Following how the grain size distribution evolves is a difficult computing task that is just recently being overtaking. Smoothed particle hydrodynamic (SPH) simulations of a single galaxy as well as cosmological simulations are producing the first predictions of the evolution of the dust grain size distribution. We compare for the first time the evolution of the dust grain size distribution predicted by the SPH simulations with the results provided by the observations. We analyse how the radial distribution of the small to large grain mass ratio (D(S)/D(L)) changes over the whole discs in three galaxies: M 101, NGC 628 and M 33. We find good agreement between the observed radial distribution of D(S)/D(L) and what is obtained from the SPH simulations of a single galaxy. The central parts of NGC 628, at high metallicity and with a high molecular gas fraction, are mainly affected not only by accretion but also by coagulation of dust grains. The centre of M 33, having lower metallicity and lower molecular gas fraction, presents an increase of D(S)/D(L), showing that shattering is very effective in creating a large fraction of small grains. Observational results provided by our galaxies confirm the general relations predicted by the cosmological simulations based on the two grain size approximation. However, we present evidence that the simulations could be overestimating the amount of large grains in high massive galaxies.
We revisit the evolution model of grain size distribution in a galaxy for the ultimate purpose of implementing it in hydrodynamical simulations. We simplify the previous model in such a way that some model-dependent assumptions are replaced with simpler functional forms. For the first test of the developed framework, we apply it to a one-zone chemical evolution model of a galaxy, confirming that our new model satisfactorily reproduces the previous results and that efficient coagulation of small grains produced by shattering and accretion is essential in reproducing the so-called MRN grain size distribution. For the next step, in order to test if our model can be treated together with the hydrodynamical evolution of the interstellar medium (ISM), we post-process a hydrodynamical simulation of an isolated disc galaxy using the new grain evolution model. We sample hydrodynamical particles representing each of the dense and diffuse ISM phases. By this post-processing, we find that the processes occurring in the dense gas (grain growth by accretion and coagulation) are important in reproducing the grain size distribution consistent with the Milky Way extinction curve. In our model, the grain size distributions are similar between the dense and diffuse ISM, although we observe a larger dispersion in the dense ISM. Moreover, we also show that even if we degrade the grain radius resolution (with 16 grid points), the overall shape of grain size distribution (and of resulting extinction curve) can be captured.
Interstellar dust grains can be spun up by radiative torques, and the resulting centrifugal force may be strong enough to disrupt large dust grains. We examine the effect of this rotational disruption on the evolution of grain size distribution in galaxies. To this goal, we modify our previous model by assuming that rotational disruption is the major small-grain production mechanism. We find that rotational disruption can have a large influence on the evolution of grain size distribution in the following two aspects especially for composites and grain mantles (with tensile strength $sim 10^7$ erg cm$^{-3}$). First, because of the short time-scale of rotational disruption, the small-grain production occurs even in the early phase of galaxy evolution. Therefore, even though stars produce large grains, the abundance of small grains can be large enough to steepen the extinction curve. Secondly, rotational disruption is important in determining the maximum grain radius, which regulates the steepness of the extinction curve. For compact grains with tensile strength $gtrsim 10^9$ erg cm$^{-3}$, the size evolution is significantly affected by rotational disruption only if the radiation field is as strong as (or the dust temperature is as high as) expected for starburst galaxies. For compact grains, rotational disruption predicts that the maximum grain radius becomes less than 0.2 $mu$m for galaxies with a dust temperature $gtrsim 50$ K.
In most debris discs, dust grain dynamics is strongly affected by stellar radiation pressure. As this mechanism is size-dependent, we expect dust grains to be spatially segregated according to their sizes. However, because of the complex interplay between radiation pressure, collisions and dynamical perturbations, this spatial segregation of the particle size distribution (PSD) has proven difficult to investigate with numerical models. We propose to explore this issue using a new-generation code that can handle some of the coupling between dynamical and collisional effects. We investigate how PSDs behave in both unperturbed discs at rest and in discs pertubed by planetary objects. We use the DyCoSS code of Thebault(2012) to investigate the coupled effect of collisions, radiation pressure and dynamical perturbations in systems having reached a steady state. We consider 2 setups: a narrow ring perturbed by an exterior planet, and an extended disc into which a planet is embedded. For both setups we consider an additional unperturbed case with no planet. We also investigate how possible spatial size segregation affect disc images at different wavelengths. We find that PSDs are always strongly spatially segregated. The only case for which they follow a standard dn/dr = C.r**(-3.5) law is for an unperturbed narrow ring, but only within the parent body ring itself. For all other configurations, the PSD can strongly depart from such power laws and have strong spatial gradients. As an example, the geometrical cross section of the disc is rarely dominated by the smallest grains on bound orbits, as it is expected to be in standard PSDs in s**q with q<-3. Although the exact profiles and spatial variations of PSDs are a complex function of the considered set-up, we are however able to derive some robust results that should be useful for image-or-SED-fitting models of observed discs.
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.
Based on a one-zone evolution model of grain size distribution in a galaxy, we calculate the evolution of infrared spectral energy distribution (SED), considering silicate, carbonaceous dust, and polycyclic aromatic hydrocarbons (PAHs). The dense gas fraction ($eta_mathrm{dense}$) of the interstellar medium (ISM), the star formation time-scale ($tau_mathrm{SF}$), and the interstellar radiation field intensity normalized to the Milky Way value ($U$) are the main parameters. We find that the SED shape generally has weak mid-infrared (MIR) emission in the early phase of galaxy evolution because the dust abundance is dominated by large grains. At an intermediate stage ($tsim 1$ Gyr for $tau_mathrm{SF}=5$ Gyr), the MIR emission grows rapidly because the abundance of small grains increases drastically by the accretion of gas-phase metals. We also compare our results with observational data of nearby and high-redshift ($zsim 2$) galaxies taken by textit{Spitzer}. We broadly reproduce the flux ratios in various bands as a function of metallicity. We find that small $eta_mathrm{dense}$ (i.e. the ISM dominated by the diffuse phase) is favoured to reproduce the 8 $mu$m intensity dominated by PAHs both for the nearby and the $zsim 2$ samples. A long $tau_mathrm{SF}$ raises the 8 $mu$m emission to a level consistent with the nearby low-metallicity galaxies. The broad match between the theoretical calculations and the observations supports our understanding of the grain size distribution, but the importance of the diffuse ISM for the PAH emission implies the necessity of spatially resolved treatment for the ISM.