Observational signatures of the dust size evolution in isolated galaxy simulations

We aim to provide observational signatures of the dust size evolution in the ISM. In particular, we explore indicators of the polycyclic aromatic hydrocarbon (PAH) mass fraction ($q_{PAH}$), defined as the mass fraction of PAHs relative to total dust grains. In addition, we validate our dust evolution model by comparing the observational signatures from our simulations to observations. We used the hydrodynamic simulation code, GADGET4-OSAKA to model the dust properties of Milky Way-like and NGC 628-like galaxies representing star-forming galaxies. This code incorporates the evolution of grain size distributions driven by dust production and interstellar processing. Furthermore, we performed post-processing dust radiative transfer with SKIRT based on the simulations to predict the observational properties. We find that the intensity ratio between 8 um and 24 um is correlated with $q_{PAH}$ and can be used as an indicator of PAH mass fraction. However, this ratio is influenced by the radiation field. We suggest the 8 um-to-total infrared intensity ratio ($\nu I_\nu(8 \mu m)/I$(TIR)) as another indicator, since it is tightly correlated with $q_{PAH}$. Furthermore, we explored the spatially resolved $q_{PAH}$ in the simulated Milky Way-like galaxy using $\nu I_\nu(8 \mu m)/I$(TIR). We find that the spatially resolved $q_{PAH}$ increases with metallicity at metallicity at Z<0.2 Zsun due to the interplay between accretion and shattering while it decreases at Z>0.2 Zsun because of coagulation. Finally, we compared the above indicators in the NGC 628-like simulation with those observed in NGC 628 by recent observations. Consequently, we find that our simulation underestimates the PAH mass fraction throughout the entire galaxy by a factor of $\sim 8$ on average. This could be due to the efficient loss of PAHs by coagulation in our model.