Showing 1–50 of 141 results for author: Hirashita, H

JWST MIRI and NIRCam observations of NGC 891 and its circumgalactic medium

Authors: Jérémy Chastenet, Ilse De Looze, Monica Relaño, Daniel A. Dale, Thomas G. Williams, Simone Bianchi, Emmanuel M. Xilouris, Maarten Baes, Alberto D. Bolatto, Martha L. Boyer, Viviana Casasola, Christopher J. R. Clark, Filippo Fraternali, Jacopo Fritz, Frédéric Galliano, Simon C. O. Glover, Karl D. Gordon, Hiroyuki Hirashita, Robert Kennicutt, Kentaro Nagamine, Florian Kirchschlager, Ralf S. Klessen, Eric W. Koch, Rebecca C. Levy, Lewis McCallum , et al. (15 additional authors not shown)

Abstract: We present new JWST observations of the nearby, prototypical edge-on, spiral galaxy NGC 891. The northern half of the disk was observed with NIRCam in its F150W and F277W filters. Absorption is clearly visible in the mid-plane of the F150W image, along with vertical dusty plumes that closely resemble the ones seen in the optical. A $\sim 10 \times 3~{\rm kpc}^2$ area of the lower circumgalactic me… ▽ More We present new JWST observations of the nearby, prototypical edge-on, spiral galaxy NGC 891. The northern half of the disk was observed with NIRCam in its F150W and F277W filters. Absorption is clearly visible in the mid-plane of the F150W image, along with vertical dusty plumes that closely resemble the ones seen in the optical. A $\sim 10 \times 3~{\rm kpc}^2$ area of the lower circumgalactic medium (CGM) was mapped with MIRI F770W at 12 pc scales. Thanks to the sensitivity and resolution of JWST, we detect dust emission out to $\sim 4$ kpc from the disk, in the form of filaments, arcs, and super-bubbles. Some of these filaments can be traced back to regions with recent star formation activity, suggesting that feedback-driven galactic winds play an important role in regulating baryonic cycling. The presence of dust at these altitudes raises questions about the transport mechanisms at play and suggests that small dust grains are able to survive for several tens of million years after having been ejected by galactic winds in the disk-halo interface. We lay out several scenarios that could explain this emission: dust grains may be shielded in the outer layers of cool dense clouds expelled from the galaxy disk, and/or the emission comes from the mixing layers around these cool clumps where material from the hot gas is able to cool down and mix with these cool cloudlets. This first set of data and upcoming spectroscopy will be very helpful to understand the survival of dust grains in energetic environments, and their contribution to recycling baryonic material in the mid-plane of galaxies. △ Less

Submitted 15 August, 2024; originally announced August 2024.

Comments: Accepted for publication in Astronomy & Astrophysics; 16 pages, 8 figures

Accurate Simultaneous Constraints on the Dust Mass, Temperature and Emissivity Index of a Galaxy at Redshift 7.31

Authors: Hiddo Algera, Hanae Inami, Ilse De Looze, Andrea Ferrara, Hiroyuki Hirashita, Manuel Aravena, Tom Bakx, Rychard Bouwens, Rebecca Bowler, Elisabete Da Cunha, Pratika Dayal, Yoshinobu Fudamoto, Jacqueline Hodge, Alexander Hygate, Ivana van Leeuwen, Themiya Nanayakkara, Marco Palla, Andrea Pallottini, Lucie Rowland, Renske Smit, Laura Sommovigo, Mauro Stefanon, Aswin Vijayan, Paul van der Werf

Abstract: We present new multi-frequency ALMA continuum observations of the massive [$\log_{10}(M_\star/M_\odot) = 10.3_{-0.2}^{+0.1}$], UV-luminous [$M_\mathrm{UV} = -21.7 \pm 0.2$] $z=7.31$ galaxy REBELS-25 in Bands 3, 4, 5, and 9. Combining the new observations with previously-taken data in Bands 6 and 8, we cover the dust continuum emission of the galaxy in six distinct bands -- spanning rest-frame… ▽ More We present new multi-frequency ALMA continuum observations of the massive [$\log_{10}(M_\star/M_\odot) = 10.3_{-0.2}^{+0.1}$], UV-luminous [$M_\mathrm{UV} = -21.7 \pm 0.2$] $z=7.31$ galaxy REBELS-25 in Bands 3, 4, 5, and 9. Combining the new observations with previously-taken data in Bands 6 and 8, we cover the dust continuum emission of the galaxy in six distinct bands -- spanning rest-frame $50-350\,μ$m -- enabling simultaneous constraints on its dust mass ($M_\mathrm{dust}$), temperature ($T_\mathrm{dust}$) and emissivity index ($β_\mathrm{IR}$) via modified blackbody fitting. Given a fiducial model of optically thin emission, we infer a cold dust temperature of $T_\mathrm{dust} = 32_{-6}^{+9}\,$K and a high dust mass of $\log_{10}(M_\mathrm{dust}/M_\odot) = 8.2_{-0.4}^{+0.6}$, and moderately optically thick dust does not significantly alter these estimates. If we assume dust production is solely through supernovae (SNe), the inferred dust yield would be high, $y = 0.7_{-0.4}^{+2.3}\,M_\odot$ per SN. Consequently, we argue grain growth in the interstellar medium of REBELS-25 also contributes to its dust build-up. This is supported by the steep dust emissivity index $β_\mathrm{IR} = 2.5 \pm 0.4$ we measure for REBELS-25, as well as by its high stellar mass, dense interstellar medium, and metal-rich nature. Our results suggest that constraining the dust emissivity indices of high-redshift galaxies is important not only to mitigate systematic uncertainties in their dust masses and obscured star formation rates, but also to assess if dust properties evolve across cosmic time. We present an efficient observing setup to do so with ALMA, combining observations of the peak and Rayleigh-Jeans tail of the dust emission. △ Less

Submitted 5 August, 2024; originally announced August 2024.

Comments: 14 pages main text, 6 figures + appendices; re-submitted to MNRAS after a positive referee report

Evolution of grain size distribution in the circum-galactic medium

Authors: Hiroyuki Hirashita

Abstract: In order to theoretically understand dust properties in the circum-galactic medium (CGM), we construct a dust evolution model that incorporates the evolution of grain size distribution. We treat each of the galaxy and the CGM as a one-zone object, and consider the mass exchange between them. We take into account dust production and interstellar dust processing for the galaxy based on our previous… ▽ More In order to theoretically understand dust properties in the circum-galactic medium (CGM), we construct a dust evolution model that incorporates the evolution of grain size distribution. We treat each of the galaxy and the CGM as a one-zone object, and consider the mass exchange between them. We take into account dust production and interstellar dust processing for the galaxy based on our previous models, and newly incorporate sputtering in the hot phase and shattering in the cool phase for the CGM. We find that shattering increases the dust destruction (sputtering) efficiency in the CGM. The functional shape of the grain size distribution in the CGM evolves following that in the galaxy, but it is sensitive to the balance between sputtering and shattering in the CGM. For an observational test, we discuss the wavelength dependence of the reddening in the CGM traced by background quasar colors, arguing that, in order to explain the observed reddening level, a rapid inflow from the CGM to the galaxy is favored because of quick dust/metal enrichment. Small grain production by shattering in the CGM also helps to explain the rise of dust extinction toward short wavelengths. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 15 pages, 7 figures, accepted for publication in PASJ

Observational signatures of the dust size evolution in isolated galaxy simulations

Authors: Kosei Matsumoto, Hiroyuki Hirashita, Kentaro Nagamine, Stefan van der Giessen, Leonard E. C. Romano, Monica Relaño, Ilse De Looze, Maarten Baes, Angelos Nersesian, Peter Camps, Kuan-chou Hou, Yuri Oku

Abstract: 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… ▽ More 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 ($νI_ν(8 μ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 $νI_ν(8 μ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. △ Less

Submitted 25 July, 2024; v1 submitted 4 February, 2024; originally announced February 2024.

Comments: 16 pages, 12 figures, accepted by A&A

Dust enrichment in the circum-galactic medium

Authors: Mau Otsuki, Hiroyuki Hirashita

Abstract: To understand the origin of dust in the circum-galactic medium (CGM), we develop a dust enrichment model. We describe each of the central galaxy and its CGM as a single zone, and consider the mass exchange between them through galactic inflows and outflows. We calculate the evolution of the gas, metal, and dust masses in the galaxy and the CGM. In the galaxy, we include stellar dust production and… ▽ More To understand the origin of dust in the circum-galactic medium (CGM), we develop a dust enrichment model. We describe each of the central galaxy and its CGM as a single zone, and consider the mass exchange between them through galactic inflows and outflows. We calculate the evolution of the gas, metal, and dust masses in the galaxy and the CGM. In the galaxy, we include stellar dust production and interstellar dust processing following our previous models. The dust in the galaxy is transported to the CGM via galactic outflows, and it is further processed by dust destruction (sputtering) in the CGM. We parameterize the time-scale or efficiency of each process and investigate the effect on the dust abundance in the CGM. We find that the resulting dust mass is sensitive to the dust destruction in the CGM, and the dust supply from galactic outflows, both of which directly regulate the dust abundance in the CGM. The inflow time-scale also affects the dust abundance in the CGM because it determines the gas mass evolution (thus, the star formation history) in the galaxy. The dust abundance in the CGM, however, is insensitive to stellar dust formation in the galaxy at later epochs because the dust production is dominated by dust growth in the interstellar medium. We also find that the resulting dust mass in the CGM is consistent with the value derived from a large sample of SDSS galaxies. △ Less

Submitted 2 February, 2024; originally announced February 2024.

Submillimetre galaxies as laboratories for dust grain coagulation

Authors: Hiroyuki Hirashita, Chian-Chou Chen

Abstract: Coagulation in the dense interstellar medium (ISM) is an important process that determines the size of the largest grains. We use submillimetre galaxies (SMGs) as laboratories of grain coagulation, since some of them host the densest ISM on a galactic scale among various populations of galaxies known. We examine how large the grains can be in such dense environments based on the mean ISM density e… ▽ More Coagulation in the dense interstellar medium (ISM) is an important process that determines the size of the largest grains. We use submillimetre galaxies (SMGs) as laboratories of grain coagulation, since some of them host the densest ISM on a galactic scale among various populations of galaxies known. We examine how large the grains can be in such dense environments based on the mean ISM density estimated from the observed typical dust mass density in SMGs. We also consider local density enhancement based on a model of supersonic turbulence, which is expected from strong stellar feedback. In the unlimited coagulation model, in which we do not impose any coagulation threshold velocity, grains as large as $\sim 30~μ$m can form under the observationally estimated mean gas density if the Mach number of turbulence is $\mathcal{M}\gtrsim 3$. We exclude this possibility since the observed emissivity index $β\simeq 2$ in the far infrared (FIR) indicates that such large grains cannot actively form in SMGs. This means that coagulation does not proceed in an unlimited way: 30-$μ$m grains should have velocities larger than the coagulation threshold. If we use a coagulation threshold (upper limit) grain velocity ($\sim 0.08$ km s$^{-1}$) taken from a theoretical study, grains likely grow only up to $μ$m size, which is small enough not to affect the FIR emissivity index. The above results indicate that SMGs can be used to constrain the physical processes relevant to coagulation. △ Less

Submitted 6 October, 2023; originally announced October 2023.

Comments: 8 pages, 4 figures, accepted for publication in MNRAS

Effects of dust grain size distribution on the abundances of CO and H$_2$ in galaxy evolution

Authors: Hiroyuki Hirashita

Abstract: We model the effect of grain size distribution in a galaxy on the evolution of CO and H$_2$ abundances. The formation and dissociation of CO and H$_2$ in typical dense clouds are modelled in a manner consistent with the grain size distribution. The evolution of grain size distribution is calculated based on our previous model, which treats the galaxy as a one-zone object but includes various dust… ▽ More We model the effect of grain size distribution in a galaxy on the evolution of CO and H$_2$ abundances. The formation and dissociation of CO and H$_2$ in typical dense clouds are modelled in a manner consistent with the grain size distribution. The evolution of grain size distribution is calculated based on our previous model, which treats the galaxy as a one-zone object but includes various dust processing mechanisms in the interstellar medium (ISM). We find that typical dense clouds become fully molecular (H$_2$) when the dust surface area increases by shattering while an increase of dust abundance by dust growth in the ISM is necessary for a significant rise of the CO abundance. Accordingly, the metallicity dependence of the CO-to-H$_2$ conversion factor, $X_\mathrm{CO}$, is predominantly driven by dust growth. We also examine the effect of grain size distribution in the galaxy by changing the dense gas fraction, which controls the balance between coagulation and shattering, clarifying that the difference in the grain size distribution significantly affects $X_\mathrm{CO}$ even if the dust-to-gas ratio is the same. The star formation time-scale, which controls the speed of metal enrichment also affects the metallicity at which the CO abundance rapidly increases (or $X_\mathrm{CO}$ drops). We also propose dust-based formulae for $X_\mathrm{CO}$, which need further tests for establishing their usefulness. △ Less

Submitted 26 April, 2023; originally announced April 2023.

Comments: 14 pages, 7 figures, accepted for publication in MNRAS

Kpc-scale properties of dust temperature in terms of dust mass and star formation activity

Authors: I-Da Chiang, Hiroyuki Hirashita, Jeremy Chastenet, Eric W. Koch, Adam K. Leroy, Erik Rosolowsky, Karin M. Sandstrom, Amy Sardone, Jiayi Sun, Thomas G. Williams

Abstract: We investigate how the dust temperature is affected by local environmental quantities, especially dust surface density ($Σ_\mathrm{dust}$), dust-to-gas ratio (D/G) and interstellar radiation field. We compile multi-wavelength observations in 46 nearby galaxies, uniformly processed with a common physical resolution of $2~$kpc. A physical dust model is used to fit the infrared dust emission spectral… ▽ More We investigate how the dust temperature is affected by local environmental quantities, especially dust surface density ($Σ_\mathrm{dust}$), dust-to-gas ratio (D/G) and interstellar radiation field. We compile multi-wavelength observations in 46 nearby galaxies, uniformly processed with a common physical resolution of $2~$kpc. A physical dust model is used to fit the infrared dust emission spectral energy distribution (SED) observed with WISE and Herschel. The star formation rate (SFR) is traced with GALEX ultraviolet data corrected by WISE infrared. We find that the dust temperature correlates well with the SFR surface density ($Σ_{\rm SFR}$), which traces the radiation from young stars. The dust temperature decreases with increasing D/G at fixed $Σ_{\rm SFR}$ as expected from stronger dust shielding at high D/G, when $Σ_\mathrm{SFR}$ is higher than $\sim 2\times 10^{-3}~\rm M_\odot~yr^{-1}~kpc^{-2}$. These measurements are in good agreement with the dust temperature predicted by our proposed analytical model. Below this range of $Σ_\mathrm{SFR}$, the observed dust temperature is higher than the model prediction and is only weakly dependent on D/G, which is possibly due to the dust heating from old stellar population or the variation of SFR within the past $10^{10}~$yr. Overall, the dust temperature as a function of $Σ_\mathrm{SFR}$ and $Σ_\mathrm{dust}$ predicted by our analytical model is consistent with observations. We also notice that at fixed gas surface density, $Σ_{\rm SFR}$ tends to increase with D/G, i.e. we can empirically modify the Kennicutt-Schmidt law with a dependence on D/G to better match observations. △ Less

Submitted 1 February, 2023; originally announced February 2023.

Comments: 16 pages, 8 figures, accepted for publication in MNRAS

Effects of dust sources on dust attenuation properties in IllustrisTNG galaxies at $z\sim 7$

Authors: Yuan-Ming Hsu, Hiroyuki Hirashita, Yen-Hsing Lin, Peter Camps, Maarten Baes

Abstract: Dust emission from high-redshift galaxies gives us a clue to the origin and evolution of dust in the early Universe. Previous studies have shown that different sources of dust (stellar dust production and dust growth in dense clouds) predict different ultraviolet (UV) extinction curves for galaxies at $z\sim 7$ but that the observed attenuation curves depend strongly on the geometry of dust and st… ▽ More Dust emission from high-redshift galaxies gives us a clue to the origin and evolution of dust in the early Universe. Previous studies have shown that different sources of dust (stellar dust production and dust growth in dense clouds) predict different ultraviolet (UV) extinction curves for galaxies at $z\sim 7$ but that the observed attenuation curves depend strongly on the geometry of dust and star distributions. Thus, we perform radiative transfer calculations under the dust-stars geometries computed by a cosmological hydrodynamic simulation (IllustrisTNG). This serves to investigate the dust attenuation curves predicted from `realistic' geometries. We choose objects with stellar mass and star formation rate appropriate for Lyman break galaxies at $z\sim 7$. We find that the attenuation curves are very different from the original extinction curves in most of the galaxies. This makes it difficult to constrain the dominant dust sources from the observed attenuation curves. We further include infrared dust emission in the analysis and plot the infrared excess (IRX)-UV spectral slope ($β$) diagram. We find that different sources of dust cause different IRX-$β$ relations for the simulated galaxies. In particular, if dust growth is the main source of dust, a variation of dust-to-metal ratio causes a more extended sequence with smaller IRX in the IRX-$β$ diagram. Thus, the comprehensive analysis of the abundances of dust and metals, the UV slope, and the dust emission could provide a clue to the dominant dust sources in the Universe. △ Less

Submitted 9 December, 2022; originally announced December 2022.

Comments: 11 pages, 6 figures, accepted for publication in MNRAS

Evolution of grain size distribution with enhanced abundance of small carbonaceous grains in galactic environments

Authors: Hiroyuki Hirashita

Abstract: We propose an updated dust evolution model that focuses on the grain size distribution in a galaxy. We treat the galaxy as a one-zone object and include five main processes (stellar dust production, dust destruction in supernova shocks, grain growth by accretion and coagulation, and grain disruption by shattering). In this paper, we improve the predictions related to small carbonaceous grains, whi… ▽ More We propose an updated dust evolution model that focuses on the grain size distribution in a galaxy. We treat the galaxy as a one-zone object and include five main processes (stellar dust production, dust destruction in supernova shocks, grain growth by accretion and coagulation, and grain disruption by shattering). In this paper, we improve the predictions related to small carbonaceous grains, which are responsible for the 2175 Å bump in the extinction curve and the polycyclic aromatic hydrocarbon (PAH) emission features in the dust emission spectral energy distribution (SED), both of which were underpredicted in our previous model. In the new model, we hypothesize that small carbonaceous grains are not involved in interstellar processing. This avoids small carbonaceous grains being lost by coagulation. We find that this hypothetical model shows a much better match to the Milky Way (MW) extinction curve and dust emission SED than the previous one. The following two additional modifications further make the fit to the MW dust emission SED better: (i) The chemical enrichment model is adjusted to give a nearly solar metallicity in the present epoch, and the fraction of metals available for dust growth is limited to half. (ii) Aromatization for small carbonaceous grains is efficient, so that the aromatic fraction is unity at grain radii $\lesssim 20$ Å. As a consequence of our modelling, we succeed in obtaining a dust evolution model that explains the MW extinction curve and dust emission SED at the same time. △ Less

Submitted 17 November, 2022; originally announced November 2022.

Comments: 11 pages, 7 figures, accepted for publication in MNRAS

Cosmic evolution of grain size distribution in galaxies using the $ν^2$GC semi-analytic model

Authors: Ryu Makiya, Hiroyuki Hirashita

Abstract: We investigate the cosmological evolution of interstellar dust with a semi-analytical galaxy formation model ($ν^2$GC), focusing on the evolution of grain size distribution. The model predicts the statistical properties of dust mass and grain size distribution in galaxies across cosmic history. We confirm that the model reproduces the relation between dust-to-gas ratio and metallicity in the local… ▽ More We investigate the cosmological evolution of interstellar dust with a semi-analytical galaxy formation model ($ν^2$GC), focusing on the evolution of grain size distribution. The model predicts the statistical properties of dust mass and grain size distribution in galaxies across cosmic history. We confirm that the model reproduces the relation between dust-to-gas ratio and metallicity in the local Universe, and that the grain size distributions of the Milky Way (MW)-like sample become similar to the so-called MRN distribution that reproduces the observed MW extinction curve. Our model, however, tends to overpredict the dust mass function at the massive end at redshift $z\lesssim 0.8$ while it reproduces the abundance of dusty galaxies at higher redshifts. We also examine the correlation between grain size distribution and galaxy properties (metallicity, specific star formation rate, gas fraction, and stellar mass), and observe a clear trend of large-grain-dominated, small-grain-dominated, and MRN-like grain size distributions from unevolved to evolved stages. As a consequence, the extinction curve shapes are flat, steep, and intermediate (MW-like) from the unevolved to evolved phases. At a fixed metallicity, the grain size distribution tends to have larger fractions of small grains at lower redshift; accordingly, the extinction curve tends to be steeper at lower redshift. We also predict that supersolar-metallicity objects at high redshift have flat extinction curves with weak 2175 Åbump strength. △ Less

Submitted 12 October, 2022; originally announced October 2022.

Comments: 12 pages, 9 figures, accepted to MNRAS

Compact dust emission in a gravitationally lensed massive quiescent galaxy at z = 2.15 revealed in ~130 pc-resolution observations by ALMA

Authors: Takahiro Morishita, Abdurro'uf, Hiroyuki Hirashita, Andrew B. Newman, Massimo Stiavelli, Marco Chiaberge

Abstract: We present new observations of MRG-M2129, a quiescent galaxy at z = 2.15 with the Atacama Large Millimeter/submillimeter Array (ALMA). With the combination of the gravitational lensing effect by the foreground cluster and the angular resolution provided by ALMA, our data reveal 1.2 mm continuum emission at $\sim130$ pc angular resolution. Compact dust continuum is detected at 7.9 $σ$ in the target… ▽ More We present new observations of MRG-M2129, a quiescent galaxy at z = 2.15 with the Atacama Large Millimeter/submillimeter Array (ALMA). With the combination of the gravitational lensing effect by the foreground cluster and the angular resolution provided by ALMA, our data reveal 1.2 mm continuum emission at $\sim130$ pc angular resolution. Compact dust continuum is detected at 7.9 $σ$ in the target but displaced from its stellar peak position by $62 \pm 38$ mas, or $\sim169 \pm 105$ pc in the source plane. We find considerably high dust-to-stellar mass ratio, $4 \times 10^{-4}$. From non-detection of the [C i] 3P2 -> 3P1 line, we then derive $3 σ$ upper limits on the molecular gas-to-dust mass ratio $δ_\mathrm{GDR} < 60$ and the molecular gas-to-stellar mass ratio fH2 < 2.3%. The derived $δ_\mathrm{GDR}$ is $>2\times$ smaller than the typical value assumed for quiescent galaxies in the literature. Our study supports that there exists a broad range of $δ_\mathrm{GDR}$ and urges submillimeter follow-up observations of quenching/recently quenched galaxies at similar redshifts. Based on the derived low $δ_\mathrm{GDR}$ and other observed dust properties, we argue that the central black hole is still active and regulates star formation in the system. Our study exhibits a rare case of a gravitationally lensed type 2 QSO harbored by a quiescent galaxy. △ Less

Submitted 7 September, 2022; v1 submitted 22 August, 2022; originally announced August 2022.

Comments: Accepted for publication in ApJ

Analytic models of dust temperature in high-redshift galaxies

Authors: Hiroyuki Hirashita, I-Da Chiang

Abstract: We investigate physical reasons for high dust temperatures ($T_\mathrm{dust}\gtrsim 40$ K) observed in some high-redshift ($z>5$) galaxies using analytic models. We consider two models that can be treated analytically: the radiative transfer (RT) model, {where a broad distribution of values for $T_\mathrm{dust}$ is considered}, and the one-tempearture (one-$T$) model, which assumes {uniform… ▽ More We investigate physical reasons for high dust temperatures ($T_\mathrm{dust}\gtrsim 40$ K) observed in some high-redshift ($z>5$) galaxies using analytic models. We consider two models that can be treated analytically: the radiative transfer (RT) model, {where a broad distribution of values for $T_\mathrm{dust}$ is considered}, and the one-tempearture (one-$T$) model, which assumes {uniform $T_\mathrm{dust}$}. These two extremes {serve to bracket the most realistic scenario}. We adopt the Kennicutt--Schmidt (KS) law to relate stellar radiation field to gas surface density, and vary the dust-to-gas ratio. As a consequence, our model is capable of predicting the relation between the surface density of star formation rate ($Σ_\mathrm{SFR}$) or dust mass ($Σ_\mathrm{dust}$) and $T_\mathrm{dust}$. We show that the high $T_\mathrm{dust}$ observed at $z\gtrsim 5$ favour low dust-to-gas ratios ($\lesssim 10^{-3}$). An enhanced star formation compared with the KS law gives an alternative explanation for the high $T_\mathrm{dust}$. The dust temperatures are similar between the two (RT and one-$T$) models as long as we use ALMA Bands 6--8. We also examine the relation among $Σ_\mathrm{SFR}$, $Σ_\mathrm{dust}$ and $T_\mathrm{dust}$ without assuming the KS law, and confirm the consistency with the actual observational data at $z>5$. In the one-$T$ model, we also examine a clumpy dust distribution, which predicts lower $T_\mathrm{dust}$ because of the leakage of stellar radiation. This enhances the requirement of low dust abundance or high star formation efficiency to explain the observed high $T_\mathrm{dust}$. △ Less

Submitted 9 August, 2022; originally announced August 2022.

Comments: 12 pages, 5 figures, Accepted for publication in MNRAS

Dust grain size evolution in local galaxies: a comparison between observations and simulations

Authors: M. Relano, I. De Looze, A. Saintonge, K. -C. Hou, L. Romano, K. Nagamine, H. Hirashita, S. Aoyama, I. Lamperti, U. Lisenfeld, M. Smith, J. Chastenet, T. Xiao, Y. Gao, M. Sargent, S. A. van der Giessen

Abstract: The evolution of the dust grain size distribution has been studied in recent years with great detail in cosmological hydrodynamical simulations taking into account all the channels under which dust evolves in the interstellar medium. We present a systematic analysis of the observed spectral energy distribution of a large sample of galaxies in the local universe in order to derive not only the tota… ▽ More The evolution of the dust grain size distribution has been studied in recent years with great detail in cosmological hydrodynamical simulations taking into account all the channels under which dust evolves in the interstellar medium. We present a systematic analysis of the observed spectral energy distribution of a large sample of galaxies in the local universe in order to derive not only the total dust masses but also the relative mass fraction between small and large dust grains (DS/DL). Simulations reproduce fairly well the observations except for the high stellar mass regime where dust masses tend to be overestimated. We find that ~45% of galaxies exhibit DS/DL consistent with the expectations of simulations, while there is a sub-sample of massive galaxies presenting high DS/DL (log(DS/DL)~-0.5), and deviating from the prediction in simulations. For these galaxies, which also have high molecular gas mass fractions and metallicities, coagulation is not an important mechanism affecting the dust evolution. Including diffusion, transporting large grains from dense regions to a more diffuse medium where they can be easily shattered, would explain the observed high DS/DL values in these galaxies. With this study we reinforce the use of the small-to-large grain mass ratio to study the relative importance of the different mechanisms in the dust life cycle. Multi-phase hydrodynamical simulations with detailed feedback prescriptions and more realistic subgrid models for the dense phase could help to reproduce the evolution of the dust grain size distribution traced by observations. △ Less

Submitted 26 July, 2022; originally announced July 2022.

Comments: 32 pages, 22 figures, 4 tables. Accepted in MNRAS

Dissecting Nearby Galaxies with piXedfit: II. Spatially Resolved Scaling Relations Among Stars, Dust, and Gas

Authors: Abdurro'uf, Yen-Ting Lin, Hiroyuki Hirashita, Takahiro Morishita, Sandro Tacchella, Po-Feng Wu, Masayuki Akiyama, Tsutomu T. Takeuchi

Abstract: We study spatially resolved scaling relations among stars, dust, and gas in ten nearby spiral galaxies. In a preceding paper Abdurro'uf et al. (2022), we have derived spatially resolved properties of the stellar population and dust by panchromatic spectral energy distribution (SED) fitting using piXedfit. Now, we investigate resolved star formation ($Σ_{\rm H_{2}}$--$Σ_{\rm SFR}$--$Σ_{*}$) and dus… ▽ More We study spatially resolved scaling relations among stars, dust, and gas in ten nearby spiral galaxies. In a preceding paper Abdurro'uf et al. (2022), we have derived spatially resolved properties of the stellar population and dust by panchromatic spectral energy distribution (SED) fitting using piXedfit. Now, we investigate resolved star formation ($Σ_{\rm H_{2}}$--$Σ_{\rm SFR}$--$Σ_{*}$) and dust scaling relations. While the relations with all sub-galactic regions of the galaxies are reasonably tight ($σ\lesssim 0.3$ dex), we find that most of the scaling relations exhibit galaxy-to-galaxy variations in normalization and shape. Only two relations of $Σ_{\rm dust}$--$Σ_{\rm gas}$ and $Σ_{\rm dust}$--$Σ_{\rm H_{2}}$ do not show noticeable galaxy-to-galaxy variations among our sample galaxies. We further investigate correlations among the scaling relations. We find significant correlations among the normalization of the $Σ_{\rm H_{2}}$--$Σ_{\rm SFR}$--$Σ_{*}$ relations, which suggest that galaxies with higher levels of resolved $\text{H}_{2}$ fraction ($f_{\rm H_{2}}$) tend to have higher levels of resolved star formation efficiency (SFE) and specific star formation rate (sSFR). We also observe that galaxies with higher levels of resolved dust-to-stellar mass ratios tend to have higher levels of resolved sSFR, SFE, and $f_{\rm H_{2}}$. Moreover, we find that galaxies with higher global sSFR and less compact morphology tend to have higher levels of the resolved sSFR, SFE, and $f_{\rm H_{2}}$, which can explain the variations in the normalization of the $Σ_{\rm H_{2}}$--$Σ_{\rm SFR}$--$Σ_{*}$ relationships. Overall, we observe indications of the contributions of both global and local factors in governing the star formation process in galaxies. △ Less

Submitted 18 July, 2022; originally announced July 2022.

Comments: 28 pages, 13 figures, accepted for publication in ApJ. piXedfit v1.0 is publicly available at https://github.com/aabdurrouf/piXedfit. Documentation is available at https://pixedfit.readthedocs.io/en/latest/. Some tutorials of practical usages in jupyter notebooks can be found at https://github.com/aabdurrouf/piXedfit/tree/main/examples

Dust SEDs in Milky Way-like galaxies in the IllustrisTNG simulations based on the evolution of grain size distribution

Authors: Chiung-Yin Chang, Yu-Hsiu Huang, Hiroyuki Hirashita, Andrew P. Cooper

Abstract: To understand how the evolution of grain size distribution in galaxies affects observed dust properties, we apply a post-processing dust evolution model to galaxy merger trees from the IllustrisTNG cosmological hydrodynamical simulation. Our dust model includes stellar dust production, sputtering in hot gas, dust growth by accretion and coagulation in the dense interstellar medium (ISM), and shatt… ▽ More To understand how the evolution of grain size distribution in galaxies affects observed dust properties, we apply a post-processing dust evolution model to galaxy merger trees from the IllustrisTNG cosmological hydrodynamical simulation. Our dust model includes stellar dust production, sputtering in hot gas, dust growth by accretion and coagulation in the dense interstellar medium (ISM), and shattering in the diffuse ISM. We decompose the grain size distribution into different dust species depending on the elemental abundances and the dense ISM fraction given by the simulation. In our previous work, we focused on Milky Way (MW) analogs and reproduced the observed MW extinction curve. In this study, we compute dust spectral energy distributions (SEDs) for the MW analogues. Our simulated SEDs broadly reproduce the observed MW SED within their dispersion and so does the observational data of nearby galaxies, although they tend to underpredict the MW SED at short wavelengths where emission is dominated by polycyclic aromatic hydrocarbons (PAHs). We find that metallicity and dense gas fraction are the most critical factors for the SED shape, through their influence on coagulation and shattering.The overall success of our models in reproducing the MW SED further justifies the dust evolution processes included in the model and predicts the dispersion in the SEDs caused by the variety in the assembly history. We also show that the most significant increase in the dust SED occurs between redshifts $z\sim 3$ and 2 in the progenitors of the simulated MW-like galaxies. △ Less

Submitted 7 April, 2022; originally announced April 2022.

Comments: 10 pages, 8 figure, accepted for publication in MNRAS

The co-evolution of molecular hydrogen and the grain size distribution in an isolated galaxy

Authors: Leonard E. C. Romano, Kentaro Nagamine, Hiroyuki Hirashita

Abstract: Understanding the evolution of dust and molecular hydrogen (H$_2$) is a critical aspect of galaxy evolution, as they affect star formation and the spectral energy distribution of galaxies. We use the $N$-body/smoothed-particle-hydrodynamics code {\sc Gadget-4} to compute the evolution of dust and H$_2$ in a suite of numerical simulations of an isolated Milky-Way-like galaxy. The evolution of the f… ▽ More Understanding the evolution of dust and molecular hydrogen (H$_2$) is a critical aspect of galaxy evolution, as they affect star formation and the spectral energy distribution of galaxies. We use the $N$-body/smoothed-particle-hydrodynamics code {\sc Gadget-4} to compute the evolution of dust and H$_2$ in a suite of numerical simulations of an isolated Milky-Way-like galaxy. The evolution of the full grain size distribution (GSD) is solved by sampling the grain size on a logarithmically spaced grid with 30 bins. The evolution of a primordial chemistry network with twelve species is solved consistently with the hydrodynamic evolution of the system, including star formation, metal and energy ejections from stars into the interstellar medium through supernova feedback and stellar winds. The formation model for H$_2$ considers the GSD and photo-dissociation through the UV radiation of young stars. We identify the processes needed for producing a sizeable amount of H$_2$, verify that the resulting star formation law in the later stages of galaxy evolution is consistent with observations of local spirals, and show that our model manages to produce a galactic molecular gas fraction in line with observations of Milky-Way-like galaxies. We stress the importance of the co-evolution of the GSD and H$_2$, as models assuming a fixed MRN shape for the GSD overestimate the production of H$_2$ in regimes where the dust abundance is dominated by large grains and underestimate it in the regime where the dust is dominated by small grains, both of which are realized in simulations of dust evolution. △ Less

Submitted 17 May, 2022; v1 submitted 11 February, 2022; originally announced February 2022.

Comments: Accepted for publication in MNRAS. 16 pages, 13 figures

Dust diffusion in SPH simulations of an isolated galaxy

Authors: Leonard E. C. Romano, Kentaro Nagamine, Hiroyuki Hirashita

Abstract: We compute the evolution of the grain size distribution (GSD) in a suite of numerical simulations of an isolated Milky-Way-like galaxy using the $N$-body/smoothed-particle-hydrodynamics code {\sc Gadget-4}. The full GSD is sampled on a logarithmically spaced grid with 30 bins, and its evolution is calculated self-consistently with the hydrodynamical and chemical evolution of the galaxy using a sta… ▽ More We compute the evolution of the grain size distribution (GSD) in a suite of numerical simulations of an isolated Milky-Way-like galaxy using the $N$-body/smoothed-particle-hydrodynamics code {\sc Gadget-4}. The full GSD is sampled on a logarithmically spaced grid with 30 bins, and its evolution is calculated self-consistently with the hydrodynamical and chemical evolution of the galaxy using a state-of-the-art star formation and feedback model. In previous versions of this model, the GSD tended to be slightly biased towards larger grains and the extinction curve had a tendency to be flatter than the observations. This work addresses these issues by considering the diffusion of dust and metals through turbulence on subgrid scales and introducing a multi-phase subgrid model that enables a smoother transition from diffuse to dense gas. We show that diffusion can significantly enhance the production of small grains and improve the agreement with the observed dust extinction curve in the Milky Way. △ Less

Submitted 17 May, 2022; v1 submitted 10 February, 2022; originally announced February 2022.

Comments: Accepted for publication in MNRAS. 20 pages, 23 figures

Evolution of dust grain size distribution and grain porosity in galaxies

Authors: Hiroyuki Hirashita, Vladimir B. Il'in

Abstract: The radiative properties of interstellar dust are affected not only by the grain size distribution but also by the grain porosity. We develop a model for the evolution of size-dependent grain porosity and grain size distribution over the entire history of galaxy evolution. We include stellar dust production, supernova dust destruction, shattering, coagulation, and accretion. Coagulation is {assume… ▽ More The radiative properties of interstellar dust are affected not only by the grain size distribution but also by the grain porosity. We develop a model for the evolution of size-dependent grain porosity and grain size distribution over the entire history of galaxy evolution. We include stellar dust production, supernova dust destruction, shattering, coagulation, and accretion. Coagulation is {assumed to be} the source of grain porosity. We use a one-zone model with a constant dense gas fraction ($η_\mathrm{dense}$), which regulates the balance between shattering and coagulation. We find that porosity develops after small grains are sufficiently created by the interplay between shattering and accretion (at age $t\sim 1$ Gyr for star formation time-scale $τ_\mathrm{SF}=5$ Gyr) and are coagulated. The filling factor drops down to 0.3 at grain radii $\sim 0.03~μ$m for $η_\mathrm{dense}=0.5$. The grains are more porous for smaller $η_\mathrm{dense}$ because small grains, from which porous coagulated grains form, are more abundant. We also calculate the extinction curves based on the above results. The porosity steepens the extinction curve significantly for silicate, but not much for amorphous carbon. The porosity also increases the collisional cross-sections and produces slightly more large grains through the enhanced coagulation; however, the extinction curve does not necessarily become flatter because of the steepening effect by porosity. We also discuss the implication of our results for the Milky Way extinction curve. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Comments: 19 pages, 10 figures, accepted for publication in MNRAS

A simple numerical experiment on the dust temperature bias for Lyman break galaxies at $z\gtrsim 5$

Authors: Yung Ying Chen, Hiroyuki Hirashita, Wei-Hao Wang, Naomasa Nakai

Abstract: Some studies suggest that the dust temperatures ($T_\mathrm{d}$) in high-redshift ($z\gtrsim 5$) Lyman break galaxies (LBGs) are high. However, possible observational bias in $T_\mathrm{d}$ is yet to be understood. Thus, we perform a simple test using random realizations of LBGs with various stellar masses, dust temperatures, and dust-to-stellar mass ratios, and examine how the sample detected by… ▽ More Some studies suggest that the dust temperatures ($T_\mathrm{d}$) in high-redshift ($z\gtrsim 5$) Lyman break galaxies (LBGs) are high. However, possible observational bias in $T_\mathrm{d}$ is yet to be understood. Thus, we perform a simple test using random realizations of LBGs with various stellar masses, dust temperatures, and dust-to-stellar mass ratios, and examine how the sample detected by ALMA is biased in terms of $T_\mathrm{d}$. We show that ALMA tends to miss high-$T_\mathrm{d}$ objects even at total dust luminosity. LBGs are, however, basically selected by the stellar UV luminosity. The dust-temperature bias in a UV-selected sample is complicated because of the competing effects between high $T_\mathrm{d}$ and low dust abundance. For ALMA Band 6, there is no tendency of high-$T_\mathrm{d}$ LBGs being more easily detected in our experiment. Thus, we suggest that the observed trend of high $T_\mathrm{d}$ in $z\gtrsim 5$ LBGs is real. We also propose that the 450 $μ$m band is useful in further clarifying the dust temperatures. To overcome the current shallowness of 450 $μ$m observations, we examine a future Antarctic 30-m class telescope with a suitable atmospheric condition for wavelengths $\lesssim 450~μ$m, where the detection is not confusion-limited. We find that, with this telescope, an $L_\mathrm{IR}$-selected sample with $\log(L_\mathrm{IR}/\mathrm{L}_{\odot})>11$ is constructed for $z\gtrsim 5$, and detection in the intermediate-$M_\star$ (stellar mass) range [$9<\log (M_\star /M_{\odot})<9.5$] is much improved, especially at high $T_\mathrm{d}$. △ Less

Submitted 26 October, 2021; originally announced October 2021.

Comments: 11 pages, 6 figures, accepted for publication in MNRAS

Dissecting Nearby Galaxies with piXedfit: I. Spatially Resolved Properties of Stars, Dust, and Gas as Revealed by Panchromatic SED Fitting

Authors: Abdurro'uf, Yen-Ting Lin, Hiroyuki Hirashita, Takahiro Morishita, Sandro Tacchella, Masayuki Akiyama, Tsutomu T. Takeuchi, Po-Feng Wu

Abstract: We study spatially resolved properties (on spatial scales of $\sim 1-2$ kpc out to at least $3$ effective radii) of the stars, dust, and gas in ten nearby spiral galaxies. The properties of the stellar population and dust are derived by fitting the spatially resolved spectral energy distribution (SED) with more than 20 photometric bands ranging from far-ultraviolet to far-infrared. Our newly devel… ▽ More We study spatially resolved properties (on spatial scales of $\sim 1-2$ kpc out to at least $3$ effective radii) of the stars, dust, and gas in ten nearby spiral galaxies. The properties of the stellar population and dust are derived by fitting the spatially resolved spectral energy distribution (SED) with more than 20 photometric bands ranging from far-ultraviolet to far-infrared. Our newly developed software piXedfit performs point spread function matching of images, pixel binning, and models the stellar light, dust attenuation, dust emission, and emission from a dusty torus heated by an active galactic nucleus simultaneously through the energy balance approach. With this self-consistent analysis, we present the spatially resolved version of the IRX--$β$ relation, finding that it is consistent with the relationship from the integrated photometry. We show that the old stellar populations contribute to the dust heating, which causes an overestimation of star formation rate (SFR) derived from the total ultraviolet and infrared luminosities on kpc scales. With archival high-resolution maps of atomic and molecular gas, we study the radial variation of the properties of the stellar populations (including stellar mass, age, metallicity, and SFR), dust (including dust mass, dust temperature, and abundance of polycyclic aromatic hydrocarbon), gas, as well as dust-to-stellar mass and dust-to-gas mass ratios. We observe a depletion of molecular gas mass fraction in the central region of the majority of the galaxies, suggesting that the lack of available fuel is an important factor in suppressing the specific SFR at the center. △ Less

Submitted 14 December, 2021; v1 submitted 6 October, 2021; originally announced October 2021.

Comments: 35 pages, 19 figures, 4 tables. Accepted by ApJ (results remain the same). Complementary figures to Fig. 4 (for 10 galaxies) and high resolution version of Fig. 8 are available at https://github.com/aabdurrouf/maps-nearby-galaxies

Geometry effects on dust attenuation curves with different grain sources at high redshift

Authors: Yen-Hsing Lin, Hiroyuki Hirashita, Peter Camps, Maarten Baes

Abstract: Dust has been detected in high-redshift ($z>5$) galaxies but its origin is still being debated. Dust production in high-redshift galaxies could be dominated by stellar production or by accretion (dust growth) in the interstellar medium. Previous studies have shown that these two dust sources predict different grain size distributions, which lead to significantly different extinction curves. In thi… ▽ More Dust has been detected in high-redshift ($z>5$) galaxies but its origin is still being debated. Dust production in high-redshift galaxies could be dominated by stellar production or by accretion (dust growth) in the interstellar medium. Previous studies have shown that these two dust sources predict different grain size distributions, which lead to significantly different extinction curves. In this paper, we investigate how the difference in the extinction curves affects the dust attenuation properties of galaxies by performing radiative transfer calculations. To examine the major effects of the dust--stars distribution geometry, we adopt two representative cases in spherical symmetry: the well-mixed geometry (stars and dust are homogeneously mixed) and the two-layer geometry (young stars are more concentrated in the centre). In both cases, we confirm that the attenuation curve can be drastically steepened by scattering and by different optical depths between young and old stellar populations, and can be flattened by the existence of unobscured stellar populations. We can reproduce similar attenuation curves even with very different extinction curves. Thus, we conclude that it is difficult to distinguish the dust sources only with attenuation curves. However, if we include information on dust emission and plot the IRX (infrared excess)--$β$ (ultraviolet spectral slope) relation, different dust sources predict different positions in the IRX--$β$ diagram. A larger $β$ is preferred under a similar IRX if dust growth is the dominant dust source. △ Less

Submitted 7 September, 2021; originally announced September 2021.

Comments: 11 pages, 4 figures, accepted for publication in MNRAS

Journal ref: MNRAS, 507, 2755-2765 (2021)

Shattering as a source of small grains in the circum-galactic medium

Authors: Hiroyuki Hirashita, Ting-Wen Lan

Abstract: Observed reddening in the circum-galactic medium (CGM) indicates a significant abundance of small grains, of which the origin is still to be clarified. We examine a possible path of small-grain production through shattering of pre-existing large grains in the CGM. Possible sites where shattering occurs on a reasonable time-scale are cool clumps with hydrogen number density $n_\mathrm{H}\sim 0.1$ c… ▽ More Observed reddening in the circum-galactic medium (CGM) indicates a significant abundance of small grains, of which the origin is still to be clarified. We examine a possible path of small-grain production through shattering of pre-existing large grains in the CGM. Possible sites where shattering occurs on a reasonable time-scale are cool clumps with hydrogen number density $n_\mathrm{H}\sim 0.1$ cm$^{-3}$ and gas temperature $T_\mathrm{gas}\sim 10^4$ K, which are shown to exist through observations of Mg II absorbers. We calculate the evolution of grain size distribution in physical conditions appropriate for cool clumps in the CGM, starting from a large-grain-dominated distribution suggested from theoretical studies. With an appropriate gas turbulence model expected from the physical condition of cold clumps (maximum eddy size and velocity of $\sim$100 pc and 10 km s$^{-1}$, respectively), together with the above gas density and temperature and the dust-to-gas mass ratio inferred from observations (0.006), we find that small-grain production occurs on a time-scale (a few $\times 10^8$ yr) comparable to the lifetime of cool clumps derived in the literature. Thus, the physical conditions of the cool clouds are favrourable for small-grain production. We also confirm that the reddening becomes significant on the above time-scale. Therefore, we conclude that small-grain production by shattering is a probable cause for the observed reddening in the CGM. We also mention the effect of grain materials (or their mixtures) on the reddening at different redshifts (1 and 2). △ Less

Submitted 13 May, 2021; originally announced May 2021.

Comments: 13 pages, 8 figures, Accepted for publication in MNRAS

Evolution of dust porosity through coagulation and shattering in the interstellar medium

Authors: Hiroyuki Hirashita, Vladimir B. Il'in, Laurent Pagani, Charlène Lefèvre

Abstract: The properties of interstellar grains, such as grain size distribution and grain porosity, are affected by interstellar processing, in particular, coagulation and shattering, which take place in the dense and diffuse interstellar medium (ISM), respectively. In this paper, we formulate and calculate the evolution of grain size distribution and grain porosity through shattering and coagulation. For… ▽ More The properties of interstellar grains, such as grain size distribution and grain porosity, are affected by interstellar processing, in particular, coagulation and shattering, which take place in the dense and diffuse interstellar medium (ISM), respectively. In this paper, we formulate and calculate the evolution of grain size distribution and grain porosity through shattering and coagulation. For coagulation, we treat the grain evolution depending on the collision energy. Shattering is treated as a mechanism of forming small compact fragments. The balance between these processes are determined by the dense-gas mass fraction $η_\mathrm{dense}$, which determines the time fraction of coagulation relative to shattering. We find that the interplay between shattering supplying small grains and coagulation forming porous grains from shattered grains is fundamentally important in creating and maintaining porosity. The porosity rises to 0.7--0.9 (or the filling factor 0.3--0.1) around grain radii $a\sim 0.1~μ$m. We also find that, in the case of $η_\mathrm{dense}=0.1$ (very efficient shattering with weak coagulation) porosity significantly enhances coagulation, creating fluffy submicron grains with filling factors lower than 0.1. The porosity enhances the extinction by 10--20 per cent at all wavelengths for amorphous carbon and at ultraviolet wavelengths for silicate. The extinction curve shape of silicate becomes steeper if we take porosity into account. We conclude that the interplay between shattering and coagulation is essential in creating porous grains in the interstellar medium and that the resulting porosity can impact the grain size distributions and extinction curves. △ Less

Submitted 6 January, 2021; originally announced January 2021.

Comments: 17 pages, 9 figures, accepted for publication in MNRAS

Evolution of the grain size distribution in Milky Way-like galaxies in post-processed IllustrisTNG simulations

Authors: Yu-Hsiu Huang, Hiroyuki Hirashita, Yun-Hsin Hsu, Yen-Ting Lin, Dylan Nelson, Andrew P. Cooper

Abstract: We model dust evolution in Milky Way-like galaxies by post-processing the IllustrisTNG cosmological hydrodynamical simulations in order to predict dust-to-gas ratios and grain size distributions. We treat grain-size-dependent dust growth and destruction processes using a 64-bin discrete grain size evolution model without spatially resolving each galaxy. Our model broadly reproduces the observed du… ▽ More We model dust evolution in Milky Way-like galaxies by post-processing the IllustrisTNG cosmological hydrodynamical simulations in order to predict dust-to-gas ratios and grain size distributions. We treat grain-size-dependent dust growth and destruction processes using a 64-bin discrete grain size evolution model without spatially resolving each galaxy. Our model broadly reproduces the observed dust--metallicity scaling relation in nearby galaxies. The grain size distribution is dominated by large grains at $z\gtrsim 3$ and the small-grain abundance rapidly increases by shattering and accretion (dust growth) at $z\lesssim 2$. The grain size distribution approaches the so-called MRN distribution at $z\sim 1$, but a suppression of large-grain abundances occurs at $z<1$. Based on the computed grain size distributions and grain compositions, we also calculate the evolution of the extinction curve for each Milky Way analogue. Extinction curves are initially flat at $z>2$, and become consistent with the Milky Way extinction curve at $z\lesssim 1$ at $1/λ< 6~\rm μm^{-1}$. However, typical extinction curves predicted by our model have a steeper slope at short wavelengths than is observed in the Milky Way. This is due to the low-redshift decline of gas-phase metallicity and the dense gas fraction in our TNG Milky Way analogues that suppresses the formation of large grains through coagulation. △ Less

Submitted 27 November, 2020; originally announced November 2020.

Comments: 16 pages, 14 figures (13 in the main text, 1 in the appendix). Accepted for publication on MNRAS

Spectral energy distributions of dust and PAHs based on the evolution of grain size distribution in galaxies

Authors: Hiroyuki Hirashita, Weining Deng, Maria S. Murga

Abstract: 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 ($η_\mathrm{dense}$) of the interstellar medium (ISM), the star formation time-scale ($τ_\mathrm{SF}$), and the interstellar radiation field… ▽ More 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 ($η_\mathrm{dense}$) of the interstellar medium (ISM), the star formation time-scale ($τ_\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 ($t\sim 1$ Gyr for $τ_\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 ($z\sim 2$) galaxies taken by \textit{Spitzer}. We broadly reproduce the flux ratios in various bands as a function of metallicity. We find that small $η_\mathrm{dense}$ (i.e.\ the ISM dominated by the diffuse phase) is favoured to reproduce the 8 $μ$m intensity dominated by PAHs both for the nearby and the $z\sim 2$ samples. A long $τ_\mathrm{SF}$ raises the 8 $μ$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. △ Less

Submitted 2 October, 2020; originally announced October 2020.

Comments: 15 pages, 10 figures, accepted for publication in MNRAS

The gas, metal and dust evolution in low-metallicity local and high-redshift galaxies

Authors: Ambra Nanni, Denis Burgarella, Patrice Theulé, Benoit Côté, Hiroyuki Hirashita

Abstract: The chemical enrichment in the interstellar medium (ISM) of galaxies is regulated by several physical processes: stellar evolution, grain formation and destruction, galactic inflows and outflows. Understanding such processes is essential to follow the chemical enrichment of galaxies through the cosmic epochs, and to interpret the observations. Despite the importance of such topics, the efficiency… ▽ More The chemical enrichment in the interstellar medium (ISM) of galaxies is regulated by several physical processes: stellar evolution, grain formation and destruction, galactic inflows and outflows. Understanding such processes is essential to follow the chemical enrichment of galaxies through the cosmic epochs, and to interpret the observations. Despite the importance of such topics, the efficiency of the different processes driving the evolution of baryons in galaxies, remain controversial. We revise the current description of metal and dust evolution in local low-metallicity dwarf galaxies and we develop a description for Lyman Break Galaxies. Our main goal is to reproduce i) the peak in the mass of dust over the mass of stars (sMdust) observed within few hundred Myrs; ii) the decrease of the sMdust at later time. The spectral energy distribution of the galaxies is fitted with the "Code Investigating GALaxies Emission" (CIGALE), through which the stellar and dust masses, and the star formation rate are estimated. For some of the dwarf galaxies, the metal and gas content are also available. We run different calculations of chemical evolution in galaxies, and we fit the observed properties through the model predictions. We show that i) a top-heavy initial mass function that favours massive stars and a dust condensation fraction for Type II Supernovae (SNe II) of 50% or more help to reproduce the peak of sMdust observed after 100 Myrs since the beginning of the cycle; ii) galactic outflows play a crucial role in reproducing the decline in sMdust with age, and they are more efficient than grain destruction from SNe II; iii) a star formation efficiency (mass of gas converted into stars) of few per cent is required to explain the metallicity of local dwarf galaxies; iv) dust growth in the ISM is not necessary to reproduce the sMdust and, if present, its effect is erased by galactic outflows. △ Less

Submitted 26 June, 2020; originally announced June 2020.

Comments: 28 pages, 17 figures, accepted for publication in A&A

Journal ref: A&A 641, A168 (2020)

Effects of rotational disruption on the evolution of grain size distribution in galaxies

Authors: Hiroyuki Hirashita, Thiem Hoang

Abstract: 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… ▽ More 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 $μ$m for galaxies with a dust temperature $\gtrsim 50$ K. △ Less

Submitted 18 March, 2020; originally announced March 2020.

Comments: 13 pages, 11 figures, Accepted for publication in MNRAS

Observational and theoretical constraints on the formation and early evolution of the first dust grains in galaxies at 5 < z < 10

Authors: D. Burgarella, A. Nanni, H. Hirashita, P. Theule, A. K. Inoue, T. T. Takeuchi

Abstract: The first generation of stars were born a few hundred million years after the big bang. These stars synthesized elements heavier than H and He, that are later expelled into the interstellar medium, initiating the rise of metals. Within this enriched medium, the first dust grains formed. This event is cosmological crucial for molecule formation as dust plays a major role by cooling low-metallicity… ▽ More The first generation of stars were born a few hundred million years after the big bang. These stars synthesized elements heavier than H and He, that are later expelled into the interstellar medium, initiating the rise of metals. Within this enriched medium, the first dust grains formed. This event is cosmological crucial for molecule formation as dust plays a major role by cooling low-metallicity star-forming clouds which can fragment to create lower mass stars. Collecting information on these first dust grains is difficult because of the negative alliance of large distances and low dust masses. We combine the observational information from galaxies at redshifts 5 < z < 10 to constrain their dust emission and theoretically understand the first evolutionary phases of the dust cycle. Spectral energy distributions (SEDs) are fitted with CIGALE and the physical parameters and their evolution are modelled. From this SED fitting, we build a dust emission template for this population of galaxies in the epoch of reionization. Our new models explain why some early galaxies are observed and others are not. We follow in time the formation of the first grains by supernovae later destroyed by other supernova blasts and expelled in the circumgalactic and intergalactic media. We have found evidence for the first dust grains formed in the universe. But, above all, this letter underlines the need to collect more data and to develop new facilities to further constrain the dust cycle in galaxies in the epoch of reionization. △ Less

Submitted 11 February, 2020; v1 submitted 5 February, 2020; originally announced February 2020.

Comments: Paper accepted for publication in section 4. Extragalactic astronomy of Astronomy and Astrophysics (Feb. 2020). Correction typo in caption of Fig. 1

Journal ref: A&A 637, A32 (2020)

Self-consistent modelling of aromatic dust species and extinction curves in galaxy evolution

Authors: Hiroyuki Hirashita, Maria S. Murga

Abstract: We formulate and calculate the evolution of dust in a galaxy focusing on the distinction among various dust components -- silicate, aromatic carbon, and non-aromatic carbon. We treat the galaxy as a one-zone object and adopt the evolution model of grain size distribution developed in our previous work. We further include aromatization and aliphatization (inverse reaction of aromatization). We rega… ▽ More We formulate and calculate the evolution of dust in a galaxy focusing on the distinction among various dust components -- silicate, aromatic carbon, and non-aromatic carbon. We treat the galaxy as a one-zone object and adopt the evolution model of grain size distribution developed in our previous work. We further include aromatization and aliphatization (inverse reaction of aromatization). We regard small aromatic grains in a radius range of 3--50 Å as polycyclic aromatic hydrocarbons (PAHs). We also calculate extinction curves in a consistent manner with the abundances of silicate and aromatic and non-aromatic carbonaceous dust. Our model nicely explains the PAH abundance as a function of metallicity in nearby galaxies. The extinction curve become similar to the Milky Way curve at age $\sim$ 10 Gyr, in terms of the carbon bump strength and the far-ultraviolet slope. We also apply our model to starburst galaxies by shortening the star formation time-scale (0.5 Gyr) and increasing the dense-gas fraction (0.9), finding that the extinction curve maintains bumpless shapes (because of low aromatic fractions), which are similar to the extinction curves observed in the Small Magellanic Cloud and high-redshift quasars. Thus, our model successfully explains the variety in extinction curve shapes at low and high redshifts. △ Less

Submitted 7 January, 2020; originally announced January 2020.

Comments: 15 pages, 9 figures, accepted for publication in MNRAS

Dust evolution in pre-stellar cores

Authors: Charlène Lefèvre, Laurent Pagani, Bilal Ladjelate, Michiel Min, Hiroyuki Hirashita, Robert Zylka

Abstract: Dust grains are the building {blocks} of future planets. They evolve in size, shape and composition during the life cycle of the interstellar medium. We seek to understand the process which leads from diffuse medium grains to dust grains in the vicinity of protostars inside disks. As a first step, we propose to characterize the dust evolution inside pre-stellar cores thanks to multi-wavelength obs… ▽ More Dust grains are the building {blocks} of future planets. They evolve in size, shape and composition during the life cycle of the interstellar medium. We seek to understand the process which leads from diffuse medium grains to dust grains in the vicinity of protostars inside disks. As a first step, we propose to characterize the dust evolution inside pre-stellar cores thanks to multi-wavelength observations. We will present how NIKA2 maps are crucial to better constrain dust properties and {we will} introduce SIGMA: a new flexible dust model in open access. △ Less

Submitted 14 November, 2019; v1 submitted 8 November, 2019; originally announced November 2019.

Comments: 6 pages, 5 figures, to appear in the proceedings of the international conference entitled mm Universe @ NIKA2, Grenoble (France), June 2019, EPJ Web of conferences

The fate of the interstellar medium in early-type galaxies. I. First direct measurement of the timescale of dust removal

Authors: Michał J. Michałowski, J. Hjorth, C. Gall, D. T. Frayer, A. -L. Tsai, H. Hirashita, K. Rowlands, T. T. Takeuchi, A. Leśniewska, D. Behrendt, N. Bourne, D. H. Hughes, E. Spring, J. Zavala, P. Bartczak

Abstract: An important aspect of quenching star formation is the removal of the cold interstellar medium (ISM; non-ionised gas and dust) from a galaxy. In addition, dust grains can be destroyed in a hot or turbulent medium. The adopted timescale of dust removal usually relies on uncertain theoretical estimates. It is tricky to track the dust removal, because usually dust is constantly replenished by consecu… ▽ More An important aspect of quenching star formation is the removal of the cold interstellar medium (ISM; non-ionised gas and dust) from a galaxy. In addition, dust grains can be destroyed in a hot or turbulent medium. The adopted timescale of dust removal usually relies on uncertain theoretical estimates. It is tricky to track the dust removal, because usually dust is constantly replenished by consecutive generations of stars. Our objective is to measure observationally the timescale of dust removal. We here explore an approach to select galaxies which do have detectable amounts of dust and cold ISM but exhibit a low current dust production rate. Any decrease of the dust and gas content as a function of the age of such galaxies therefore must be attributed to processes governing the ISM removal. We used a sample of galaxies detected by Herschel in the far-infrared with visually assigned early-type morphology or spirals with red colours. We also obtained JCMT/SCUBA-2 observations for five of them. We discovered an exponential decline of the dust-to-stellar mass ratio with age, which we interpret as an evolutionary trend of dust removal from these galaxies. For the first time we directly measure the dust removal timescale in such galaxies to be tau=(2.5+-0.4) Gyr (the corresponding half-life time is (1.75+-0.25) Gyr). This quantity may be used in models in which it must be assumed a priori and cannot be derived. Any process which removes dust in these galaxies, such as dust grain destruction, cannot happen on shorter timescales. The timescale is comparable to the quenching timescales found in simulations for galaxies with similar stellar masses. The dust is likely of internal, not external origin. It was either formed in the past directly by supernovae, or from seeds produced by SNe and with grain growth in the ISM contributing substantially to the dust mass accumulation. △ Less

Submitted 14 October, 2019; originally announced October 2019.

Comments: Astronomy & Astrophysics, accepted; 13 pages, 9 figures, 1 table

Journal ref: A&A 632, A43 (2019)

Expected dust grain size distributions in galaxies detected by ALMA at $z>7$

Authors: Hsin-Min Liu, Hiroyuki Hirashita

Abstract: The dust properties in high-redshift galaxies provide clues to the origin of dust in the Universe. Although dust has been detected in galaxies at redshift $z>7$, it is difficult to constrain the dominant dust sources only from the total dust amount. Thus, we calculate the evolution of grain size distribution, expecting that different dust sources predict different grain size distributions. Using t… ▽ More The dust properties in high-redshift galaxies provide clues to the origin of dust in the Universe. Although dust has been detected in galaxies at redshift $z>7$, it is difficult to constrain the dominant dust sources only from the total dust amount. Thus, we calculate the evolution of grain size distribution, expecting that different dust sources predict different grain size distributions. Using the star formation time-scale and the total baryonic mass constrained by the data in the literature, we calculate the evolution of grain size distribution. To explain the total dust masses in ALMA-detected $z>7$ galaxies, the following two solutions are possible: (i) high dust condensation efficiency in stellar ejecta, and (ii) efficient accretion (dust growth by accreting the gas-phase metals in the interstellar medium). We find that these two scenarios predict significantly different grain size distributions: in (i), the dust is dominated by large grains ($a\gtrsim 0.1~μ$m, where $a$ is the grain radius), while in (ii), the small-grain ($a\lesssim 0.01~μ$m) abundance is significantly enhanced by accretion. Accordingly, extinction curves are expected to be much steeper in (ii) than in (i). Thus, we conclude that extinction curves provide a viable way to distinguish the dominant dust sources in the early phase of galaxy evolution. △ Less

Submitted 19 September, 2019; v1 submitted 18 September, 2019; originally announced September 2019.

Comments: 10 pages, 8 figures, accepted for publication in MNRAS

Modelling the evolution of PAH abundance in galaxies

Authors: Shiau-Jie Rau, Hiroyuki Hirashita, Maria Sergeevna Murga

Abstract: We investigate the evolution of polycyclic aromatic hydrocarbon (PAH) abundance in a galaxy, which is a crucial step to understand the evolution of bright emission features in the mid-infrared range. We calculate the evolution of dust grain size distribution in a manner consistent with the physical conditions of the interstellar medium by post-processing our previous hydrodynamical simulation of a… ▽ More We investigate the evolution of polycyclic aromatic hydrocarbon (PAH) abundance in a galaxy, which is a crucial step to understand the evolution of bright emission features in the mid-infrared range. We calculate the evolution of dust grain size distribution in a manner consistent with the physical conditions of the interstellar medium by post-processing our previous hydrodynamical simulation of an isolated disc galaxy. We also differentiate between aromatic and non-aromatic grains for carbonaceous dust species and explicitly considered the aromatization process. As a consequence, our model explains the metallicity dependence of PAH abundances in nearby galaxies well. The PAH abundance increase is driven particularly by the interplay between shattering and accretion (dust growth). The fast aromatization guarantees that the small carbonaceous grains trace PAHs very well. Since shattering and accretion are sensitive to the dust abundance, we predict that the PAH-to-dust abundance ratio increases as the metallicity increases. This is consistent with the observation data of nearby galaxies. △ Less

Submitted 6 September, 2019; originally announced September 2019.

Comments: Accepted by Monthly Notices of the Royal Astronomical Society (2019/9/5) 7 pages 4 figures

Galaxy simulation with the evolution of grain size distribution

Authors: Shohei Aoyama, Hiroyuki Hirashita, Kentaro Nagamine

Abstract: We compute the evolution of interstellar dust in a hydrodynamic simulation of an isolated disc galaxy. We newly implement the evolution of full grain size distribution by sampling 32 grid points on the axis of the grain radius. We solve it consistently with the chemical enrichment and hydrodynamic evolution of the galaxy. This enables us to theoretically investigate spatially resolved evolution of… ▽ More We compute the evolution of interstellar dust in a hydrodynamic simulation of an isolated disc galaxy. We newly implement the evolution of full grain size distribution by sampling 32 grid points on the axis of the grain radius. We solve it consistently with the chemical enrichment and hydrodynamic evolution of the galaxy. This enables us to theoretically investigate spatially resolved evolution of grain size distribution in a galaxy. The grain size distribution evolves from a large-grain-dominated ($\gtrsim 0.1~μ$m) phase to a small-grain production phase, eventually converging to a power-law-like grain size distribution similar to the so-called MRN distribution. We find that the small-grain abundance is higher in the dense ISM in the early epoch ($t\lesssim 1$ Gyr) because of efficient dust growth by accretion, while coagulation makes the small-grain abundance less enhanced in the dense ISM later. This leads to steeper extinction curves in the dense ISM than in the diffuse ISM in the early phase, while they show the opposite trend later. The radial trend is also described by faster evolution in the inner part. We also confirm that the simulation reproduces the observed trend in the relation between dust-to-gas ratio and metallicity, and in the radial gradients of dust-to-gas ratio and dust-to-metal ratio. Since the above change in the grain size distribution occurs at $t\sim 1$ Gyr, the age and density dependence of grain size distribution has a significant impact on the extinction curves even at high redshift. △ Less

Submitted 24 November, 2019; v1 submitted 5 June, 2019; originally announced June 2019.

Comments: Accepted for publication into MNRAS, 16 pages, 8 figures

Radiation-pressure-driven dust transport to galaxy halos at $z\sim 10$

Authors: Hiroyuki Hirashita, Akio K. Inoue

Abstract: The origin of dust in galaxy halos or in the circum-galactic medium (CGM) is still a mystery. We investigate if the radiation pressure in high-redshift ($z\sim 10$) galaxies can efficiently transport dust to halos. To clarify the first dust enrichment of galaxy halos in the early Universe, we solve the motion of a dust grain considering radiation pressure, gas drag, and gravity in the vertical dir… ▽ More The origin of dust in galaxy halos or in the circum-galactic medium (CGM) is still a mystery. We investigate if the radiation pressure in high-redshift ($z\sim 10$) galaxies can efficiently transport dust to halos. To clarify the first dust enrichment of galaxy halos in the early Universe, we solve the motion of a dust grain considering radiation pressure, gas drag, and gravity in the vertical direction of the galactic disc. Radiation pressure is estimated in a consistent manner with the stellar spectra and dust extinction. As a consequence, we find that dust grains with radii $a\sim 0.1~μ$m successfully escape from the galactic disc if the ongoing star formation episode converts more than 15 per cent of the baryon content into stars and lasts $\gtrsim 30$ Myr, while larger and smaller grains are trapped in the disc because of gravity and gas drag, respectively. We also show that grain charge significantly enhances gas drag at a few--10 scale heights of the galactic disc, where the grain velocities are suppressed to $\sim 1$ km s$^{-1}$. There is an optimum dust-to-gas ratio ($\sim 10^{-3}$) in the galactic disc and an optimum virial mass $\sim 10^{10}$--$10^{11}$ M$_{\odot}$ for the transport of $a\sim 0.1~μ$m grains to the halo. We conclude that early dust enrichment of galaxy halos at $z\gtrsim 10$ is important for the origin of dust in the CGM. △ Less

Submitted 14 May, 2019; originally announced May 2019.

Comments: 14 pages, 8 figures, accepted for publication in MNRAS

Dust scaling relations in a cosmological simulation

Authors: Kuan-Chou Hou, Shohei Aoyama, Hiroyuki Hirashita, Kentaro Nagamine, Ikkoh Shimizu

Abstract: To study the dust evolution in the cosmological structure formation history, we perform a smoothed particle hydrodynamic simulation with a dust enrichment model in a cosmological volume. We adopt the dust evolution model that represents the grain size distribution by two sizes and takes into account stellar dust production and interstellar dust processing. We examine the dust mass function and the… ▽ More To study the dust evolution in the cosmological structure formation history, we perform a smoothed particle hydrodynamic simulation with a dust enrichment model in a cosmological volume. We adopt the dust evolution model that represents the grain size distribution by two sizes and takes into account stellar dust production and interstellar dust processing. We examine the dust mass function and the scaling properties of dust in terms of the characteristics of galaxies. The simulation broadly reproduces the observed dust mass functions at redshift $z = 0$, except that it overproduces the massive end at dust mass $M_\mathrm{d} \gtrsim 10^{8}$ ${\rm M}_\odot$. This overabundance is due to overproducing massive gas/metal-rich systems, but we also note that the relation between stellar mass and gas-phase metallicity is reproduced fairly well by our recipe. The relation between dust-to-gas ratio and metallicity shows a good agreement with the observed one at $z=0$, which indicates successful implementation of dust evolution in our cosmological simulation. Star formation consumes not only gas but also dust, causing a decreasing trend of the dust-to-stellar mass ratio at the high-mass end of galaxies. We also examine the redshift evolution up to $z \sim~ 5$, and find that the galaxies have on average the highest dust mass at $z = 1-2$. For the grain size distribution, we find that galaxies with metallicity $\sim 0.3~ Z_\odot$ tend to have the highest small-to-large grain abundance ratio; consequently, the extinction curves in those galaxies have the steepest ultraviolet slopes. △ Less

Submitted 9 January, 2019; originally announced January 2019.

Comments: 18 pages, 9 figures, accepted for publication in MNRAS

Remodeling the evolution of grain size distribution in galaxies

Authors: Hiroyuki Hirashita, Shohei Aoyama

Abstract: 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 o… ▽ More 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. △ Less

Submitted 18 October, 2018; originally announced October 2018.

Comments: 18 pages, 10 figures, accepted for publication in MNRAS

Comparison of cosmological simulations and deep submillimetre galaxy surveys

Authors: Shohei Aoyama, Hiroyuki Hirashita, Chen-Fatt Lim, Yu-Yen Chang, Wei-Hao Wang, Kentaro Nagamine, Kuan-Chou Hou, Ikkoh Shimizu, Hui-Hsuan Chung, Chien-Hsiu Lee, Xian-Zhong Zheng

Abstract: Recent progress in submillimetre surveys by single-dish telescopes allows us to further challenge the consistency between cosmological simulations and observations. In particular, we compare our simulations that include dust formation and destruction with the recent SCUBA-2 surveys (`STUDIES') by putting emphases on basic observational properties of dust emission such as dust temperature, size of… ▽ More Recent progress in submillimetre surveys by single-dish telescopes allows us to further challenge the consistency between cosmological simulations and observations. In particular, we compare our simulations that include dust formation and destruction with the recent SCUBA-2 surveys (`STUDIES') by putting emphases on basic observational properties of dust emission such as dust temperature, size of infrared (IR)-emitting region, IR luminosity function and IRX--$β$ relation. After confirming that our models reproduce the local galaxy properties, we examine the STUDIES sample at $z\approx 1-4$, finding that the simulation reproduces the aforementioned quantities except for the $z\gtrsim 2$ IR luminosity function at the massive end ($\sim 10^{13}$ L$_{\odot}$). This means that the current simulation correctly reproduces the overall scaling between the size and luminosity (or star formation rate) of dusty region, but lacks extreme starburst phenomena at $z\gtrsim 2$. We also discuss extinction curves and possible AGN contribution. △ Less

Submitted 21 January, 2019; v1 submitted 27 September, 2018; originally announced September 2018.

Comments: 13 pages, 10 figures, published in MNRAS

Dust abundance and grain size in galaxy halos

Authors: Hiroyuki Hirashita, Chih-Yu Lin

Abstract: We investigate the abundance and properties (especially, grain size) of dust in galaxy halos using available observational data in the literature. There are two major sets of data. One is (i) the reddening curves at redshifts $z\sim 1$ and 2 derived for Mg II absorbers, which are assumed to trace the medium in galaxy halos. The other is (ii) the cosmic extinction up to $z\sim 2$ mainly traced by d… ▽ More We investigate the abundance and properties (especially, grain size) of dust in galaxy halos using available observational data in the literature. There are two major sets of data. One is (i) the reddening curves at redshifts $z\sim 1$ and 2 derived for Mg II absorbers, which are assumed to trace the medium in galaxy halos. The other is (ii) the cosmic extinction up to $z\sim 2$ mainly traced by distant background quasars. For (i), the observed reddening curves favor a grain radius of $a\sim 0.03~μ$m for silicate, while graphite is not supported because of its strong 2175 Å bump. Using amorphous carbon improves the fit to the reddening curves compared with graphite if the grain radius is $a\lesssim 0.03~μ$m. For (ii), the cosmic extinction requires $η\gtrsim 10^{-2}$ ($η$ is the ratio of the halo dust mass to the stellar mass; the observationally suggested value is $η\sim 10^{-3}$) for silicate if $a\sim 0.03~μ$m as suggested by the reddening curve constraint. Thus, for silicate, we do not find any grain radius that satisfies both (i) and (ii) unless the halo dust abundance is much larger than suggested by the observations. For amorphous carbon, in contrast, a wide range of grain radius ($a\sim 0.01$--0.3~$μ$m) is accepted by the cosmic extinction; thus, we find that a grain radius range of $a\sim 0.01$--0.03 $μ$m is supported by combining (i) and (ii). We also discuss the origin of dust in galaxy halos, focusing on the importance of grain size in the physical mechanism of dust supply to galaxy halos. △ Less

Submitted 3 April, 2018; originally announced April 2018.

Comments: 13 pages, 9 figures, accepted for publication in Planetary and Space Science (special issue for Cosmic Dust X)

Cosmological simulation with dust formation and destruction

Authors: Shohei Aoyama, Kuan-Chou Hou, Hiroyuki Hirashita, Kentaro Nagamine, Ikkoh Shimizu

Abstract: To investigate the evolution of dust in a cosmological volume, we perform hydrodynamic simulations, in which the enrichment of metals and dust is treated self-consistently with star formation and stellar feedback. We consider dust evolution driven by dust production in stellar ejecta, dust destruction by sputtering, grain growth by accretion and coagulation, and grain disruption by shattering, and… ▽ More To investigate the evolution of dust in a cosmological volume, we perform hydrodynamic simulations, in which the enrichment of metals and dust is treated self-consistently with star formation and stellar feedback. We consider dust evolution driven by dust production in stellar ejecta, dust destruction by sputtering, grain growth by accretion and coagulation, and grain disruption by shattering, and treat small and large grains separately to trace the grain size distribution. After confirming that our model nicely reproduces the observed relation between dust-to-gas ratio and metallicity for nearby galaxies, we concentrate on the dust abundance over the cosmological volume in this paper. The comoving dust mass density has a peak at redshift $z\sim 1$--2, coincident with the observationally suggested dustiest epoch in the Universe. {In the local Universe}, roughly 10 per cent of the dust is contained in the intergalactic medium (IGM), where only 1/3--1/4 of the dust survives against dust destruction by sputtering. We also show that the dust mass function is roughly reproduced at $\lesssim 10^8$ M$_\odot$, while the massive end still has a discrepancy, which indicates {the necessity of stronger feedback in massive galaxies}. %%The relation showed that accretion is essential for dusty galaxies. In addition, our model broadly reproduces the observed radial profile of dust surface density in the circum-galactic medium (CGM). While our model satisfies the observational constraints for the dust extinction {on cosmological scales}, it predicts that the dust in the CGM and IGM is dominated by large ($> 0.03~μ$m) grains, which is in tension with the steep reddening curves {observed} in the CGM. △ Less

Submitted 22 July, 2018; v1 submitted 12 February, 2018; originally announced February 2018.

Comments: 17 pages, 9 figures. published on MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, 2018, Volume 478, Issue 4, Pages 4905-4921

Exploring the Dust Population in Cold Diffuse Clouds

Authors: Steven J. Gibson, Hiroyuki Hirashita, Aaron C. Bell, Mary E. Spraggs, Alberto Noriega-Crespo, Sean J. Carey, William T. Reach, Christopher M. Brunt

Abstract: The formation and evolution of cold diffuse clouds (CDCs), the parent objects of dense molecular clouds, affects both the star formation process and that of larger-scale galactic evolution. We have begun a pilot study of one CDC's dust content, with the goal of quantifying the abundances of different types of dust and relating these to the relative abundance of molecular gas, the cloud's physical… ▽ More The formation and evolution of cold diffuse clouds (CDCs), the parent objects of dense molecular clouds, affects both the star formation process and that of larger-scale galactic evolution. We have begun a pilot study of one CDC's dust content, with the goal of quantifying the abundances of different types of dust and relating these to the relative abundance of molecular gas, the cloud's physical properties, and its general stage of development. Using photometry from AKARI and other surveys, we have extracted a sample spectral energy distribution (SED) of the CDC dust thermal emission over the near-infrared to submillimeter range. The extracted SED closely resembles others in the literature, confirming our isolation of the cloud emission from other sources along the sight line. We plan to fit this SED with dust models at each position in the cloud, automating our procedure to map out the structure of this CDC and others. △ Less

Submitted 4 December, 2017; originally announced December 2017.

Comments: 4 pages, 3 figures, to appear in proceedings of AKARI 2017 Conference, "The Cosmic Wheel and the Legacy of the AKARI Archive: From Galaxies and Stars to Planets and Life"

Populating H$_2$ and CO in galaxy simulation with dust evolution

Authors: Li-Hsin Chen, Hiroyuki Hirashita, Kuan-Chou Hou, Shohei Aoyama, Ikkoh Shimizu, Kentaro Nagamine

Abstract: There are two major theoretical issues for the star formation law (the relation between the surface densities of molecular gas and star formation rate on a galaxy scale): (i) At low metallicity, it is not obvious that star-forming regions are rich in H$_2$ because the H$_2$ formation rate depends on the dust abundance; and (ii) whether or not CO really traces H$_2$ is uncertain, especially at low… ▽ More There are two major theoretical issues for the star formation law (the relation between the surface densities of molecular gas and star formation rate on a galaxy scale): (i) At low metallicity, it is not obvious that star-forming regions are rich in H$_2$ because the H$_2$ formation rate depends on the dust abundance; and (ii) whether or not CO really traces H$_2$ is uncertain, especially at low metallicity. To clarify these issues, we use a hydrodynamic simulation of an isolated disc galaxy with a spatial resolution of a few tens parsecs. The evolution of dust abundance and grain size distribution is treated consistently with the metal enrichment and the physical state of the interstellar medium. We compute the H$_2$ and CO abundances using a subgrid post-processing model based on the dust abundance and the dissociating radiation field calculated in the simulation. We find that when the metallicity is $\lesssim 0.4$ Z$_\odot$ ($t<1$ Gyr), H$_2$ is not a good tracer of star formation rate because H$_2$-rich regions are limited to dense compact regions. At $Z\gtrsim 0.8$ Z$_\odot$, a tight star formation law is established for both H$_2$ and CO. At old ($t \sim 10$ Gyr) ages, we also find that adopting the so-called MRN grain size distribution with an appropriate dust-to-metal ratio over the entire disc gives reasonable estimates for the H$_2$ and CO abundances. For CO, improving the spatial resolution of the simulation is important while the H$_2$ abundance is not sensitive to sub-resolution structures at $Z\gtrsim 0.4$ Z$_\odot$. △ Less

Submitted 1 November, 2017; originally announced November 2017.

Comments: 21 pages, 14 figures, accepted for publication in MNRAS

Dust masses of $z>5$ galaxies from SED fitting and ALMA upper limits

Authors: Hiroyuki Hirashita, Denis Burgarella, Rychard J. Bouwens

Abstract: We aim at constraining the dust mass in high-redshift ($z\gtrsim 5$) galaxies using the upper limits obtained by ALMA in combination with the rest-frame UV--optical spectral energy distributions (SEDs). For SED fitting, because of degeneracy between dust extinction and stellar age, we focus on two extremes: continuous star formation (Model A) and instantaneous star formation (Model B). We apply th… ▽ More We aim at constraining the dust mass in high-redshift ($z\gtrsim 5$) galaxies using the upper limits obtained by ALMA in combination with the rest-frame UV--optical spectral energy distributions (SEDs). For SED fitting, because of degeneracy between dust extinction and stellar age, we focus on two extremes: continuous star formation (Model A) and instantaneous star formation (Model B). We apply these models to Himiko (as a representative UV-bright object) and a composite SED of $z>5$ Lyman break galaxies (LBGs). For Himiko, Model A requires a significant dust extinction, which leads to a high dust temperature $>70$ K for consistency with the ALMA upper limit. This high dust temperature puts a strong upper limit on the total dust mass $M_\mathrm{d}\lesssim 2\times 10^6$ M$_{\odot}$, and the dust mass produced per supernova (SN) $m_\mathrm{d,SN}\lesssim 0.1$ M$_{\odot}$. Such a low $m_\mathrm{d,SN}$ suggests significant loss of dust by reverse shock destruction or outflow, and implies that SNe are not the dominant source of dust at high $z$. Model B allows $M_\mathrm{d}\sim 2\times 10^7$ M$_{\odot}$ and $m_\mathrm{d,SN}\sim 0.3$ M$_{\odot}$.} We could distinguish between Models A and B if we observe Himiko at {wavelength $<$ 1.2 mm by ALMA. For the LBG sample, we obtain $M_\mathrm{d}\lesssim 2\times 10^6$ M$_{\odot}$ for a typical LBG at $z>5$, but this only puts an upper limit for $m_\mathrm{d,SN}$ as $\sim 2$ M$_{\odot}$. This clarifies the importance of observing UV-bright objects (like Himiko) to constrain the dust production by SNe. △ Less

Submitted 7 September, 2017; originally announced September 2017.

Comments: 12 pages, 6 figures, accepted for publication in MNRAS

Evolution of dust extinction curves in galaxy simulation

Authors: Kuan-Chou Hou, Hiroyuki Hirashita, Kentaro Nagamine, Shohei Aoyama, Ikkoh Shimizu

Abstract: To understand the evolution of extinction curve, we calculate the dust evolution in a galaxy using smoothed particle hydrodynamics simulations incorporating stellar dust production, dust destruction in supernova shocks, grain growth by accretion and coagulation, and grain disruption by shattering. The dust species are separated into carbonaceous dust and silicate. The evolution of grain size distr… ▽ More To understand the evolution of extinction curve, we calculate the dust evolution in a galaxy using smoothed particle hydrodynamics simulations incorporating stellar dust production, dust destruction in supernova shocks, grain growth by accretion and coagulation, and grain disruption by shattering. The dust species are separated into carbonaceous dust and silicate. The evolution of grain size distribution is considered by dividing grain population into large and small gains, which allows us to estimate extinction curves. We examine the dependence of extinction curves on the position, gas density, and metallicity in the galaxy, and find that extinction curves are flat at $t \lesssim 0.3$ Gyr because stellar dust production dominates the total dust abundance. The 2175 Å bump and far-ultraviolet (FUV) rise become prominent after dust growth by accretion. At $t \gtrsim 3$ Gyr, shattering works efficiently in the outer disc and low density regions, so extinction curves show a very strong 2175 Å bump and steep FUV rise. The extinction curves at $t\gtrsim 3$ Gyr are consistent with the Milky Way extinction curve, which implies that we successfully included the necessary dust processes in the model. The outer disc component caused by stellar feedback has an extinction curves with a weaker 2175 Å bump and flatter FUV slope. The strong contribution of carbonaceous dust tends to underproduce the FUV rise in the Small Magellanic Cloud extinction curve, which supports selective loss of small carbonaceous dust in the galaxy. The snapshot at young ages also explain the extinction curves in high-redshift quasars. △ Less

Submitted 6 April, 2017; originally announced April 2017.

Comments: 16 pages, 10 figures, accepted for publication in MNRAS

Dust evolution with active galactic nucleus feedback in elliptical galaxies

Authors: Hiroyuki Hirashita, Takaya Nozawa

Abstract: We have recently suggested that dust growth in the cold gas phase dominates the dust abundance in elliptical galaxies while dust is efficiently destroyed in the hot X-ray emitting plasma (hot gas). In order to understand the dust evolution in elliptical galaxies, we construct a simple model that includes dust growth in the cold gas and dust destruction in the hot gas. We also take into account the… ▽ More We have recently suggested that dust growth in the cold gas phase dominates the dust abundance in elliptical galaxies while dust is efficiently destroyed in the hot X-ray emitting plasma (hot gas). In order to understand the dust evolution in elliptical galaxies, we construct a simple model that includes dust growth in the cold gas and dust destruction in the hot gas. We also take into account the effect of mass exchange between these two gas components induced by active galactic nucleus (AGN) feedback. We survey reasonable ranges of the relevant parameters in the model and find that AGN feedback cycles actually produce a variety in cold gas mass and dust-to-gas ratio. By comparing with an observational sample of nearby elliptical galaxies, we find that, although the dust-to-gas ratio varies by an order of magnitude in our model, the entire range of the observed dust-to-gas ratios is difficult to be reproduced under a single parameter set. Variation of the dust growth efficiency is the most probable solution to explain the large variety in dust-to-gas ratio of the observational sample. Therefore, dust growth can play a central role in creating the variation in dust-to-gas ratio through the AGN feedback cycle and through the variation in dust growth efficiency. △ Less

Submitted 25 January, 2017; originally announced January 2017.

Comments: 14 pages, 9 figures, accepted for publication in Planetary and Space Science (special issue for Cosmic Dust IX)

Effects of dust evolution on the abundances of CO and H$_2$

Authors: Hiroyuki Hirashita, Nanase Harada

Abstract: The CO-to-H$_2$ conversion factor ($X_\mathrm{CO}$) is known to correlate with the metallicity ($Z$). The dust abundance, which is related to the metallicity, is responsible for this correlation through dust shielding of dissociating photons and H$_2$ formation on dust surfaces. In this paper, we investigate how the relation between dust-to-gas ratio and metallicity ($\mathcal{D}$--$Z$ relation) a… ▽ More The CO-to-H$_2$ conversion factor ($X_\mathrm{CO}$) is known to correlate with the metallicity ($Z$). The dust abundance, which is related to the metallicity, is responsible for this correlation through dust shielding of dissociating photons and H$_2$ formation on dust surfaces. In this paper, we investigate how the relation between dust-to-gas ratio and metallicity ($\mathcal{D}$--$Z$ relation) affects the H$_2$ and CO abundances (and $X_\mathrm{CO}$) of a `molecular' cloud. For the $\mathcal{D}$--$Z$ relation, we adopt a dust evolution model developed in our previous work, which treats the evolution of not only dust abundance but also grain sizes in a galaxy. Shielding of dissociating photons and H$_2$ formation on dust are solved consistently with the dust abundance and grain sizes. As a consequence, our models {predict consistent metallicity dependence of $X_\mathrm{CO}$ with observational data}. Among various processes driving dust evolution, grain growth by accretion has the largest impact on the $X_\mathrm{CO}$--$Z$ relation. The other processes also have some impacts on the $X_\mathrm{CO}$--$Z$ relation, but their effects are minor compared with the scatter of the observational data at the metallicity range ($Z\gtrsim 0.1$ Z$_\odot$) where CO could be detected. We also find that dust condensation in stellar ejecta has a dramatic impact on the H$_2$ abundance at low metallicities ($\lesssim 0.1$ Z$_\odot$), relevant for damped Lyman $α$ systems and nearby dwarf galaxies, and that the grain size dependence of H$_2$ formation rate is also important. △ Less

Submitted 17 January, 2017; originally announced January 2017.

Comments: 18 pages, 7 figures, accepted for publication in MNRAS

Constraint on dust evolution processes in normal galaxies at $z>6$ detected by ALMA

Authors: W. -C. Wang, H. Hirashita, K. -C. Hou

Abstract: Recent ALMA observations of high-redshift normal galaxies have been providing a great opportunity to clarify the general origin of dust in the Universe, not biased to very bright special objects even at $z>6$. To clarify what constraint we can get for the dust enrichment in normal galaxies detected by ALMA, we use a theoretical model that includes major processes driving dust evolution in a galaxy… ▽ More Recent ALMA observations of high-redshift normal galaxies have been providing a great opportunity to clarify the general origin of dust in the Universe, not biased to very bright special objects even at $z>6$. To clarify what constraint we can get for the dust enrichment in normal galaxies detected by ALMA, we use a theoretical model that includes major processes driving dust evolution in a galaxy; that is, dust condensation in stellar ejecta, dust growth by the accretion of gas-phase metals, and supernova destruction. Using the dust emission fluxes detected in two normal galaxies at $z>6$ by ALMA as a constraint, we can get the range of the time-scales (or efficiencies) of the above mentioned processes. We find that if we assume extremely high condensation efficiency in stellar ejecta ($f_{\mathrm{in}} \ga 0.5$), rapid dust enrichment by stellar sources in the early phase may be enough to explain the observed ALMA flux, unless dust destruction by supernovae in those galaxies is stronger than that in nearby galaxies. If we assume a condensation efficiency expected from theoretical calculations ($f_{\mathrm{in}} \la 0.1$), strong dust growth (even stronger than assumed for nearby galaxies if they are metal-poor galaxies) is required. These results indicate that the normal galaxies detected by ALMA at $z>6$ are biased to objects (i) with high dust condensation efficiency in stellar ejecta, (ii) with strong dust growth in very dense molecular clouds, or (iii) with efficient dust growth because of fast metal enrichment up to solar metallicity. A measurement of metallicity is crucial to distinguish among these possibilities. △ Less

Submitted 21 November, 2016; originally announced November 2016.

Comments: 8 pages, 4 figures. Accepted for publication in MNRAS

Galaxy Simulation with Dust Formation and Destruction

Authors: Shohei Aoyama, Kuan-Chou Hou, Ikkoh Shimizu, Hiroyuki Hirashita, Keita Todoroki, Jun-Hwan Choi, Kentaro Nagamine

Abstract: We perform smoothed particle hydrodynamics (SPH) simulations of an isolated galaxy with a new treatment for dust formation and destruction. To this aim, we treat dust and metal production self-consistently with star formation and supernova feedback. For dust, we consider a simplified model of grain size distribution by representing the entire range of grain sizes with large and small grains. We in… ▽ More We perform smoothed particle hydrodynamics (SPH) simulations of an isolated galaxy with a new treatment for dust formation and destruction. To this aim, we treat dust and metal production self-consistently with star formation and supernova feedback. For dust, we consider a simplified model of grain size distribution by representing the entire range of grain sizes with large and small grains. We include dust production in stellar ejecta, dust destruction by supernova (SN) shocks, grain growth by accretion and coagulation, and grain disruption by shattering. We find that the assumption of fixed dust-to-metal mass ratio becomes no longer valid when the galaxy is older than 0.2 Gyr, at which point the grain growth by accretion starts to contribute to the nonlinear rise of dust-to-gas ratio. As expected in our previous one-zone model, shattering triggers grain growth by accretion since it increases the total surface area of grains. Coagulation becomes significant when the galaxy age is greater than $\sim$ 1 Gyr: at this epoch the abundance of small grains becomes high enough to raise the coagulation rate of small grains. We further compare the radial profiles of dust-to-gas ratio $(\mathcal{D})$ and dust-to-metal ratio $(\mathcal{D}/Z)$ (i.e., depletion) at various ages with observational data. We find that our simulations broadly reproduce the radial gradients of dust-to-gas ratio and depletion. In the early epoch ($\lesssim 0.3$ Gyr), the radial gradient of $\mathcal{D}$ follows the metallicity gradient with $\mathcal{D}/Z$ determined by the dust condensation efficiency in stellar ejecta, while the $\mathcal{D}$ gradient is steeper than the $Z$ gradient at the later epochs because of grain growth by accretion. The framework developed in this paper is applicable to any SPH-based galaxy evolution simulations including cosmological ones. △ Less

Submitted 18 January, 2017; v1 submitted 23 September, 2016; originally announced September 2016.

Comments: 18 pages, 16 figures, published on MNRAS

Journal ref: MNRAS (2017) 466 (1): 105-121

Dust evolution processes constrained by extinction curves in nearby galaxies

Authors: Kuan-Chou Hou, Hiroyuki Hirashita, Michał J. Michałowski

Abstract: Extinction curves, especially those in the Milky Way (MW), the Large Magellanic Cloud (LMC), and the Small Magellanic Cloud (SMC), have provided us with a clue to the dust properties in the nearby Universe. We examine whether or not these extinction curves can be explained by well known dust evolution processes. We treat the dust production in stellar ejecta, destruction in supernova shocks, dust… ▽ More Extinction curves, especially those in the Milky Way (MW), the Large Magellanic Cloud (LMC), and the Small Magellanic Cloud (SMC), have provided us with a clue to the dust properties in the nearby Universe. We examine whether or not these extinction curves can be explained by well known dust evolution processes. We treat the dust production in stellar ejecta, destruction in supernova shocks, dust growth by accretion and coagulation, and dust disruption by shattering. To make a survey of the large parameter space possible, we simplify the treatment of the grain size distribution evolution by adopting the `two-size approximation', in which we divide the grain population into small ($\lesssim 0.03~μ$m) and large ($\gtrsim 0.03~μ$m) grains. It is confirmed that the MW extinction curve can be reproduced in reasonable ranges for the time-scale of the above processes with a silicate-graphite mixture. This indicates that the MW extinction curve is a natural consequence of the dust evolution through the above processes. We also find that the same models fail to reproduce the SMC/LMC extinction curves. Nevertheless, this failure can be remedied by giving higher supernova destruction rates for small carbonaceous dust and considering amorphous carbon for carbonaceous dust; these modification fall in fact in line with previous studies. Therefore, we conclude that the current dust evolution scenario composed of the aforementioned processes is successful in explaining the extinction curves. All the extinction curves favor efficient interstellar processing of dust, especially, strong grain growth by accretion and coagulation. △ Less

Submitted 22 August, 2016; originally announced August 2016.

Comments: 19 pages, 14 figures