Search

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)--$\beta$ (ultraviolet spectral slope) relation, different dust sources predict different positions in the IRX--$\beta$ diagram. A larger $\beta$ is preferred under a similar IRX if dust growth is the dominant dust source., Comment: 11 pages, 4 figures, accepted for publication in MNRAS

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 $z2$, 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., 16 pages, 14 figures (13 in the main text, 1 in the appendix). Accepted for publication on MNRAS

Based on a one-zone evolution model of grain size distribution in a galaxy, we calculate the evolution of infrared spectral energy distribution (SED), considering silicate, carbonaceous dust, and polycyclic aromatic hydrocarbons (PAHs). The dense gas fraction ($\eta_\mathrm{dense}$) of the interstellar medium (ISM), the star formation time-scale ($\tau_\mathrm{SF}$), and the interstellar radiation field intensity normalized to the Milky Way value ($U$) are the main parameters. We find that the SED shape generally has weak mid-infrared (MIR) emission in the early phase of galaxy evolution because the dust abundance is dominated by large grains. At an intermediate stage ($t\sim 1$ Gyr for $\tau_\mathrm{SF}=5$ Gyr), the MIR emission grows rapidly because the abundance of small grains increases drastically by the accretion of gas-phase metals. We also compare our results with observational data of nearby and high-redshift ($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 $\eta_\mathrm{dense}$ (i.e.\ the ISM dominated by the diffuse phase) is favoured to reproduce the 8 $\mu$m intensity dominated by PAHs both for the nearby and the $z\sim 2$ samples. A long $\tau_\mathrm{SF}$ raises the 8 $\mu$m emission to a level consistent with the nearby low-metallicity galaxies. The broad match between the theoretical calculations and the observations supports our understanding of the grain size distribution, but the importance of the diffuse ISM for the PAH emission implies the necessity of spatially resolved treatment for the ISM., Comment: 15 pages, 10 figures, accepted for publication in MNRAS

Interstellar dust grains can be spun up by radiative torques, and the resulting centrifugal force may be strong enough to disrupt large dust grains. We examine the effect of this rotational disruption on the evolution of grain size distribution in galaxies. To this goal, we modify our previous model by assuming that rotational disruption is the major small-grain production mechanism. We find that rotational disruption can have a large influence on the evolution of grain size distribution in the following two aspects especially for composites and grain mantles (with tensile strength $\sim 10^7$ erg cm$^{-3}$). First, because of the short time-scale of rotational disruption, the small-grain production occurs even in the early phase of galaxy evolution. Therefore, even though stars produce large grains, the abundance of small grains can be large enough to steepen the extinction curve. Secondly, rotational disruption is important in determining the maximum grain radius, which regulates the steepness of the extinction curve. For compact grains with tensile strength $\gtrsim 10^9$ erg cm$^{-3}$, the size evolution is significantly affected by rotational disruption only if the radiation field is as strong as (or the dust temperature is as high as) expected for starburst galaxies. For compact grains, rotational disruption predicts that the maximum grain radius becomes less than 0.2 $\mu$m for galaxies with a dust temperature $\gtrsim 50$ K., Comment: 13 pages, 11 figures, Accepted for publication in MNRAS

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~\mu$m, where $a$ is the grain radius), while in (ii), the small-grain ($a\lesssim 0.01~\mu$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., Comment: 10 pages, 8 figures, accepted for publication in MNRAS

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., Comment: Accepted by Monthly Notices of the Royal Astronomical Society (2019/9/5) 7 pages 4 figures

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~\mu$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~\mu$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., Comment: 14 pages, 8 figures, accepted for publication in MNRAS

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--$\beta$ 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., Comment: 13 pages, 10 figures, published in MNRAS

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., 18 pages, 10 figures, accepted for publication in MNRAS

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 $\mu$m band is useful in further clarifying the dust temperatures. To overcome the current shallowness of 450 $\mu$m observations, we examine a future Antarctic 30-m class telescope with a suitable atmospheric condition for wavelengths $\lesssim 450~\mu$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, Comment: 11 pages, 6 figures, accepted for publication in MNRAS

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 ($\eta_\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 $\tau_\mathrm{SF}=5$ Gyr) and are coagulated. The filling factor drops down to 0.3 at grain radii $\sim 0.03~\mu$m for $\eta_\mathrm{dense}=0.5$. The grains are more porous for smaller $\eta_\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., Comment: 19 pages, 10 figures, accepted for publication in MNRAS

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., 28 pages, 17 figures, accepted for publication in A&A

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 \AA\ 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., Comment: 15 pages, 9 figures, accepted for publication in MNRAS

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., Comment: Paper accepted for publication in section 4. Extragalactic astronomy of Astronomy and Astrophysics (Feb. 2020). Correction typo in caption of Fig. 1

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.

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., 18 pages, 9 figures, accepted for publication in MNRAS

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., Comment: Astronomy & Astrophysics, accepted; 13 pages, 9 figures, 1 table

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~\mu$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., Comment: Accepted for publication into MNRAS, 16 pages, 8 figures

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~\mu$m for silicate, while graphite is not supported because of its strong 2175 \AA\ bump. Using amorphous carbon improves the fit to the reddening curves compared with graphite if the grain radius is $a\lesssim 0.03~\mu$m. For (ii), the cosmic extinction requires $\eta\gtrsim 10^{-2}$ ($\eta$ is the ratio of the halo dust mass to the stellar mass; the observationally suggested value is $\eta\sim 10^{-3}$) for silicate if $a\sim 0.03~\mu$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~$\mu$m) is accepted by the cosmic extinction; thus, we find that a grain radius range of $a\sim 0.01$--0.03 $\mu$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., Comment: 13 pages, 9 figures, accepted for publication in Planetary and Space Science (special issue for Cosmic Dust X)

The discoveries of huge amounts of dust and unusual extinction curves in high-redshift quasars (z > 4) cast challenging issues on the origin and properties of dust in the early universe. In this Letter, we investigate the evolutions of dust content and extinction curve in a high-z quasar, based on the dust evolution model taking account of grain size distribution. First, we show that the Milky-Way extinction curve is reproduced by introducing a moderate fraction (~0.2) of dense molecular-cloud phases in the interstellar medium for a graphite-silicate dust model. Then we show that the peculier extinction curves in high-z quasars can be explained by taking a much higher molecular-cloud fraction (>0.5), which leads to more efficient grain growth and coagulation, and by assuming amorphous carbon instead of graphite. The large dust content in high-z quasar hosts is also found to be a natural consequence of the enhanced dust growth. These results indicate that grain growth and coagulation in molecular clouds are key processes that can increase the dust mass and change the size distribution of dust in galaxies, and that, along with a different dust composition, can contribute to shape the extinction curve., 5 pages, 4 figures, Accepted for publication in MNRAS Letters

Coreshine in dense molecular cloud cores (dense cores) is interpreted as evidence for micrometer-sized grains (referred to as very large grains, VLGs). VLGs may have a significant influence on the total dust amount and the extinction curve. We estimate the total abundance of VLGs in the Galaxy, assuming that dense cores are the site of VLG formation. We find that the VLG abundance relative to the total dust mass is roughly $\phi_\mathrm{VLG}\sim 0.01(1-\epsilon )/\epsilon (\tau_\mathrm{SF}/5\times 10^9~\mathrm{yr})^{-1} (f_\mathrm{VLG}/0.5)(t_\mathrm{shat}/10^8~\mathrm{yr})$, where $\epsilon$ is the star formation efficiency in dense cores, $\tau_\mathrm{SF}$ the timescale of gas consumption by star formation, $f_\mathrm{VLG}$ the fraction of dust mass eventually coagulated into VLGs in dense cores, and $t_\mathrm{shat}$ the lifetime of VLGs (determined by shattering). Adopting their typical values for the Galaxy, we obtain $\phi_\mathrm{VLG}\sim 0.02$--0.09. This abundance is well below the value detected in the heliosphere by Ulysses and Galileo, which means that local enhancement of VLG abundance in the solar neighborhood is required if the VLGs originate from dense cores. We also show that the effects of VLGs on the extinction curve are negligible even with the upper value of the above range, $\phi_\mathrm{VLG}\sim 0.09$. If we adopt an extreme value, $\phi_\mathrm{VLG}\sim 0.5$, close to that inferred from the above spacecraft data, the extinction curve is still in the range of the variation in Galactic extinction curves, but is not typical of the diffuse ISM., Comment: 8 pages, 2 figures, Accepted for publication in Planetary and Space Science

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 \AA\ 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 \AA\ 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 \AA\ 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., Comment: 16 pages, 10 figures, accepted for publication in MNRAS

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 $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., Comment: 12 pages, 6 figures, accepted for publication in MNRAS

We investigate the evolution of extinction curves in galaxies based on our evolution model of grain size distribution. In this model, we considered various processes: dust formation by SNe II and AGB stars, dust destruction by SN shocks in the ISM, metal accretion onto the surface of grains (referred to as grain growth), shattering and coagulation. We find that the extinction curve is flat in the earliest stage of galaxy evolution. As the galaxy is enriched with dust, shattering becomes effective to produce a large abundance of small grains ($a \la 0.01\;\mu$m). Then, grain growth becomes effective at small grain radii, forming a bump at $a \sim 10^{-3}${--}$10^{-2}\;\mu$m on the grain size distribution. Consequently, the extinction curve at ultraviolet (UV) wavelengths becomes steep, and a bump at $1/\lambda \sim 4.5\;\mu{\rm m}^{-1}\;(\lambda: \mbox{wavelength})$ on the extinction curve becomes prominent. Once coagulation becomes effective, the extinction curves become flatter, but the UV extinction remains overproduced when compared with the Milky Way extinction curve. This discrepancy can be resolved by introducing a stronger contribution of coagulation. Consequently, an interplay between shattering and coagulation could be important to reproduce the Milky Way extinction curve., Comment: 10 pages, 8 figures, accepted for publication in MNRAS

Shattering of dust grains in the interstellar medium is a viable mechanism of small grain production in galaxies. We examine the robustness or uncertainty in the theoretical predictions of shattering. We identify $P_1$ (the critical pressure above which the deformation destroys the original lattice structures) as the most important quantity in determining the timescale of small grain production, and confirm that the same $P_1/t$ ($t$ is the duration of shattering) gives the same grain size distribution [$n(a)$, where $a$ is the grain radius] after shattering within a factor of 3. The uncertainty in the fraction of shocked material that is eventually ejected as fragments causes uncertainties in $n(a)$ by a factor of 1.3 and 1.6 for silicate and carbonaceous dust, respectively. The size distribution of shattered fragments have minor effects as long as the power index of the fragment size distribution is less than ~ 3.5, since the slope of grain size distribution $n(a)$ continuously change by shattering and becomes consistent with $n(a)\propto a^{-3.5}$. The grain velocities as a function of grain radius can have an imprint in the grain size distribution especially for carbonaceous dust. We also show that the formulation of shattering can be simplified without losing sufficient precision., Comment: 12 pages, 7 figures, Accepted for publication in Earth, Planets, and Space (Special Issue: Cosmic Dust V)

The central star-forming regions in three blue compact dwarf galaxies (He 2-10, NGC 5253, and II Zw 40) were observed in the 340 GHz (880 micron) band at 5 arcsec resolution with the Submillimetre Array (SMA). Continuum emission associated with the central star-forming complex was detected in all these galaxies. The SMA 880 micron flux is decomposed into free-free emission and dust emission by using centimetre-wavelength data in the literature. We find that free-free emission contributes half or more of the SMA 880 micron flux in the central starbursts in those three galaxies. In spite of the dominance of free-free emission at 880 micron, the radio-to-far infrared (FIR) ratios in the central star-forming regions are not significantly higher than those of the entire systems, showing the robustness of radio-FIR relation. Based on the robustness of the radio-FIR relation, we argue that the free--free fraction in the 880 micron emission is regulated by the dust temperature. We also analyze the CO (J = 3--2) emission data. We find that CO is a good tracer of the total gas mass in solar-metallicity object He 2-10. Low-metallicity objects, NGC 5253 and II Zw 40, have apparently high star formation efficiencies; however, this may be an artifact of significant dissociation of CO in the low-metallicity environments. We also point out a potential underestimate of dust mass, since the dust traced by emission is biased to the most luminous high-temperature regions, particularly when a system hosts a compact star-forming region where the dust temperature is high., Comment: 13 pages, 6 figures, accepted for publication in MNRAS

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., 18 pages, 16 figures, published on MNRAS

We present a theoretical model aimed at explaining the IRX-$��$ relation for high redshift (z >5) galaxies. Recent observations (Capak+2015; Bouwens+2016) have shown that early Lyman Break Galaxies, although characterized by a large UV attenuation (e.g. flat UV beta slopes), show a striking FIR deficit, i.e. they are "infrared-dark". This marked deviation from the local IRX-beta relation can be explained by the larger molecular gas content of these systems. While dust in the diffuse ISM attains relatively high temperatures (Td = 45 K for typical size a=0.1 um; smaller grains can reach Td = 60 K), a sizable fraction of the dust mass is embedded in dense gas, and therefore remains cold. If confirmed, the FIR deficit might represent a novel, powerful indicator of the molecular content of high-z galaxies which can be used to pre-select candidates for follow-up deep CO observations. Thus, high-z CO line searches with ALMA might be much more promising than currently thought., 8 pages, 4 Figures, MNRAS Submitted

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~\mu$m) and large ($\gtrsim 0.03~\mu$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., Comment: 19 pages, 14 figures

In high-redshift ($z>5$) quasars, a large amount of dust ($\textstyle\sim 10^{8} \mathrm{M}_{\sun}$) has been observed. In order to explain the large dust content, we focus on a possibility that grain growth by the accretion of heavy elements is the dominant dust source. We adopt a chemical evolution model applicable to nearby galaxies but utilize parameters adequate to high-$z$ quasars. It is assumed that metals and dust are predominantly ejected by Type II supernovae (SNe). We have found that grain growth strongly depends on the grain size distribution. If we simply use the size distribution of grains ejected from SNe, grain growth is inefficient because of the lack of small grains (i.e.\ small surface-to-volume ratio of the dust grains). However, if we take small grain production by interstellar shattering into consideration, grain growth is efficient enough to account for the rich dust abundance in high-$z$ quasars. Our results not only confirm that grain growth is necessary to explain the large amount of dust in high-$z$ quasars, but also demonstrate that grain size distributions have a critical impact on grain growth.

Dust grains grow in interstellar clouds by accretion and coagulation. In this paper, we focus on these two grain growth processes and numerically investigate how they interplay to increase the grain radii. We show that accretion efficiently depletes grains with radii on a time-scale of 10 Myr in solar-metallicity molecular clouds. Coagulation also occurs on a similar time-scale, but accretion is more efficient in producing a large bump in the grain size distribution. Coagulation further pushes the grains to larger sizes after a major part of the gas-phase metals are used up. Similar grain sizes are achieved by coagulation regardless of whether accretion takes place or not; in this sense, accretion and coagulation modify the grain size distribution independently. The increase of the total dust mass in a cloud is also investigated. We show that coagulation slightly ‘suppresses’ dust mass growth by accretion but that this effect is slight enough to be neglected in considering the grain mass budget in galaxies. Finally, we examine how accretion and coagulation affect the extinction curve: the ultraviolet slope and the carbon bump are enhanced by accretion, while they are flattened by coagulation.

Ground-based observations at terahertz (THz) frequencies are a newly explorable area of astronomy for the next ten years. We discuss science cases for a first-generation 10-m class THz telescope, focusing on the Greenland Telescope as an example of such a facility. We propose science cases and provide quantitative estimates for each case. The largest advantage of ground-based THz telescopes is their higher angular resolution (~ 4 arcsec for a 10-m dish), as compared to space or airborne THz telescopes. Thus, high-resolution mapping is an important scientific argument. In particular, we can isolate zones of interest for Galactic and extragalactic star-forming regions. The THz windows are suitable for observations of high-excitation CO lines and [N II] 205 um lines, which are scientifically relevant tracers of star formation and stellar feedback. Those lines are the brightest lines in the THz windows, so that they are suitable for the initiation of ground-based THz observations. THz polarization of star-forming regions can also be explored since it traces the dust population contributing to the THz spectral peak. For survey-type observations, we focus on ``sub-THz'' extragalactic surveys, whose uniqueness is to detect galaxies at redshifts z ~ 1--2, where the dust emission per comoving volume is the largest in the history of the Universe. Finally we explore possibilities of flexible time scheduling, which enables us to monitor active galactic nuclei, and to target gamma-ray burst afterglows. For these objects, THz and submillimeter wavelength ranges have not yet been explored., 39 pages, 17 figures, accepted for publication in PASJ (as a review paper)

Grain growth by the accretion of metals in interstellar clou ds (called ‘grain growth’) could be one of the dominant processes that determine the dust content in galaxies. The importance of grain size distribution for the grain growth is demo nstrated in this paper. First, we derive an analytical formula that gives the grain size distr ibution after the grain growth in individual clouds for any initial grain size distribution. The time-scale of the grain growth is very sensitive to grain size distribution, since the grai n growth is mainly regulated by the surface to volume ratio of grains. Next, we implement the results of grain growth into dust enrichment models of entire galactic system along with the grain formation and destruction in the interstellar medium, finding that the grain growth in c louds governs the dust content in nearby galaxies unless the grain size is strongly biased to sizes larger than � 0.1 µm or the power index of the grain size distribution is shallower t han � 2.5. The grain growth in clouds contributes to the rapid increase of dust-to-gas r atio at a certain metallicity level (called critical metallicity in Asano et al. 2011 and Inoue 2 011), which we find to be sensitive to grain size distribution. Thus, the grain growth efficient ly increase the dust mass not only in nearby galaxies but also in high-redshift quasars, whose metallicities are larger than the critical value. Our recipe for the grain growth is applicabl e for any grain size distribution and easily implemented into any framework of dust enrichment in galaxies.

The importance of the far-infrared (FIR) mapping is demonstrated for a face-on spiral galaxy, M81, by analysing its imaging data at 65, 90 and 140 μm taken by AKARI. The basic products are the dust temperature map, the dust optical depth map and the colour-colour diagram. The main features are as follows. (i) The dust temperature derived from the total fluxes at 90 and 140 μm reflects the relatively low temperatures seen in the interarm and spiral arms excluding the warm spots rather than the high temperatures in warm spots and the centre. This indicates that the total FIR luminosity is dominated by the dust heated by the general interstellar radiation field. (ii) The galaxy is more extended at 140 μm than at the other shorter wavelengths, which reflects the radial dust temperature gradient. (iii) The dust optical depth derived from the FIR mapping is broadly consistent with that estimated from the FIR-to-ultraviolet luminosity ratio. (iv) The FIR colour-colour diagram is useful to identify a 'contamination' of warm dust. The existence of small-scale warm star-forming regions is supported in the bright spots along the spiral arms. This contamination also leads to an underestimate of dust optical depth (or dust column density).

We aim at understanding the statistical properties of luminous sub-millimeter (submm) galaxies (SMGs) in the context of cosmological structure formation. By utilizing a cosmological N-body simulation to calculate the distribution of dark halos in the Universe, we consider the dust enrichment in individual halos by Type II supernovae (SNe II). The SN II rate is estimated under a star formation activity which is assumed to occur on a dynamical timescale in the dark matter potential. Our simple framework successfully explains the luminosity function, the typical star formation rate, and the typical dust mass of an observational SMG sample at z∼3. We also examine the clustering properties of SMGs, since a positive cross correlation between SMGs and Lyα emitters (LAEs) is indeed observed by a recent observation. In the simulation, we select SMGs by FIR dust luminosity >1012L⊙, while LAEs are chosen such that the age and the virial mass are consistent with the observed LAE properties. The SMGs and LAEs selected in this way show a spatial cross correlation whose strength is consistent with the observation. This confirms that the SMGs really trace the most clustered regions at z∼3 and that their luminosities can be explained by the dust accumulation as a result of their star formation activities. We extend our prediction to higher redshifts, finding that a statistical sample of submm galaxies at z≥6 can be obtained by ALMA with a 100 arcmin2 survey. With the same survey, a few submm galaxies at z∼10 may be detected.

Using the results from a previously developed Lyman Alpha/continuum production/transmission and dust enrichment model for Lyman Alpha Emitters (LAEs), based on cosmological SPH simulations, we assess the detectability of their dust-reprocessed sub-millimeter (submm) radiation. As supernovae (rather than evolved stars) control dust formation and destruction processes, LAEs are relatively dust-poor with respect to local galaxies: they have low dust-to-gas ratios (0.05 times the dust-to-gas ratio of the Milky Way) in spite of their relatively high metallicity, which is approximately 10-50% of the solar value. Using the derived escape fraction of ultraviolet (UV) continuum photons we compute the UV luminosity absorbed by dust and re-emitted in the far infrared. The LAE submm fluxes correlate with their Lyman Alpha luminosity: about (3%, 1%) at z=(5.7, 6.6) of the LAEs in our simulated sample (those with with Log L_\alpha > 43.1) would have fluxes at 850 micrometer (the optimal band for detection) in excess of 0.12 mJy and will be therefore detectable at 5 sigma with ALMA with an integration time of only 1 hour. Such detections would open a new window on the physical conditions prevailing in these most distant galaxies.

We investigate the size-density relation of H ii regions in blue compact dwarf galaxies (BCDs) by compiling observational data of their size (Di) and electron density (ne). We find that the size-density relation follows a relation with constant column density (ne ∝ Di−1) rather than with constant luminosity (ne ∝ Di−1.5). Such behavior resembles that of Galactic H ii regions, and may imply an underlying “scale-free” connection. Because this size-density relation cannot be explained by static models, we model and examine the evolution of the size-density relation of H ii regions by considering the star formation history and pressure-driven expansion of H ii regions. We find that the size-density relation of the entire BCD sample does not result from an evolutionary sequence of H ii regions but rather reflects a sequence with different initial gas densities (or “hierarchy” of density). We also find that the dust extinction of ionizing photons is significant for the BCD sample, despite their blue optical colors. This means that as long as the emission from H ii regions is used to trace massive star formation, we would miss the star formation activity in dense environments even in low-metallicity galaxies such as BCDs.

The Large and Small Magellanic Clouds (LMC and SMC, respectively) are observed to have characteristic dust extinction curves that are quite different from those of the Galaxy (e.g., strength of the 2175 A bump). Although the dust composition and size distribution of the Magellanic Clouds (MCs) that can self-consistently explain their observed extinction curves have been already proposed, it remain unclear whether and how the required dust properties can be achieved in the formation histories of the MCs. We therefore investigate the time evolution of the dust properties of the MCs and thereby derive their extinction curves using one-zone chemical evolution models with formation and evolution of small and large silicate and carbonaceous dust grains and dusty winds associated with starburst events. We find that the observed SMC extinction curve without a conspicuous 2175 A bump can be reproduced well by our SMC model, if the small carbon grains can be selectively lost through the dust wind during the latest starburst about 0.2 Gyr ago. We also find that the LMC extinction curve with a weak 2175 A bump can be reproduced by our LMC model with less efficient removal of dust through dust wind. We discuss possible physical reasons for different dust wind efficiencies between silicate and graphite and among galaxies., Comment: 31 pages, 13 figures, accepted for the publication in ApJ

Full calculations of the evolution of grain size distribution in galaxies are in general computationally heavy. In this paper, we propose a simple model of dust enrichment in a galaxy with a simplified treatment of grain size distribution by imposing a `two-size approximation'; that is, all the grain population is represented by small (grain radius a < 0.03 micron) and large (a > 0.03 micron) grains. We include in the model dust supply from stellar ejecta, destruction in supernova shocks, dust growth by accretion, grain growth by coagulation and grain disruption by shattering, considering how these processes work on the small and large grains. We show that this simple framework reproduces the main features found in full calculations of grain size distributions as follows. The dust enrichment starts with the supply of large grains from stars. At a metallicity level referred to as the critical metallicity of accretion, the abundance of the small grains formed by shattering becomes large enough to rapidly increase the grain abundance by accretion. Associated with this epoch, the mass ratio of the small grains to the large grains reaches the maximum. After that, this ratio converges to the value determined by the balance between shattering and coagulation, and the dust-to-metal ratio is determined by the balance between accretion and shock destruction. With a Monte Carlo simulation, we demonstrate that the simplicity of our model has an advantage in predicting statistical properties. We also show some applications for predicting observational dust properties such as extinction curves., 14 pages, 12 figures, accepted for publication in MNRAS

The small grain sizes produced by Type II supernova (SN II) models in young, metal-poor galaxies make the appearance of their infrared (IR) spectral energy distribution (SED) quite different from that of nearby, older galaxies. To study this effect, we have developed a model for the evolution of dust content and the IR SED of low-metallicity, extremely young galaxies based on Hirashita et al. (2002). We find that, even in the intense ultraviolet (UV) radiation field of very young galaxies, small silicate grains are subject to stochastic heating resulting in a broad temperature distribution and substantial MIR continuum emission. Larger carbonaceous grains are in thermal equilibrium at T \simeq 50 - 100K, and they also contribute to the MIR. We present the evolution of SEDs and IR extinction of very young, low-metallicity galaxies. The IR extinction curve is also shown. In the first few Myrs, the emission peaks at \lambda \sim 30-50um at later times dust self-absorption decreases the apparent grain temperatures, shifting the bulk of the emission into the submillimetre band. We successfully apply the model to the IR SED of a low metallicity (1/41 Z_\odot) dwarf galaxy SBS0335-052. We find the SED, optical properties and extinction of the star forming region to be consistent with a very young and compact starburst. We also predict the SED of another extremely low-metallicity galaxy, I Zw 18, for future observational tests. Some prospects for future observations are discussed., Comment: MNRAS in press, pages, 6 figures, using mn2e.cls. Abstract abridged

Stars form out of molecular gas and supply dust grains during their last evolutionary stages; in turn hydrogen molecules (H2) are produced more efficiently on dust grains. Therefore, dust can drastically accelerate H2 formation, leading to an enhancement of star formation activity. In order to examine the first formation of stars and dust in galaxies, we model the evolution of galaxies in the redshift range of 55) galaxies in sub-millimetre and near-infrared bands. We find that: i) ALMA can detect dust emission from several thousands of galaxies per square degree, and ii) NGST can detect the stellar emission from 10^6 galaxies per square degree. Further observational checks of our predictions include the integrated flux of metal (oxygen and carbon) lines. We finally discuss possible color selection strategies for high-redshift galaxy searches., Comment: 19 pages, LaTeX, 10 figures, accepted to MNRAS

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 $\alpha$ systems and nearby dwarf galaxies, and that the grain size dependence of H$_2$ formation rate is also important., Comment: 18 pages, 7 figures, accepted for publication in MNRAS

The ubiquity of planets poses an interesting question: when first planets are formed in galaxies. We investigate this problem by adopting a theoretical model developed for understanding the statistical properties of exoplanets. Our model is constructed as the combination of planet traps with the standard core accretion scenario in which the efficiency of forming planetary cores directly relates to the dust density in disks or the metallicity ([Fe/H]). We statistically compute planet formation frequencies (PFFs) as well as the orbital radius ($$) within which gas accretion becomes efficient enough to form Jovian planets. The three characteristic exoplanetary populations are considered: hot Jupiters, exo-Jupiters densely populated around 1 AU, and low-mass planets such as super-Earths. We explore the behavior of the PFFs as well as $$ for the three different populations as a function of metallicity ($-2 \leq$[Fe/H]$\leq -0.6$). We show that the total PFFs increase steadily with metallicity, which is the direct outcome of the core accretion picture. For the entire range of the metallicity considered here, the population of the low-mass planets dominates over the Jovian planets. The Jovian planets contribute to the PFFs above [Fe/H]-1. We find that the hot Jupiters form at lower metallcities than the exo-Jupiters. This arises from the radially inward transport of planetary cores by their host traps, which is more effective for lower metallicity disks due to the slower growth of the cores. The PFFs for the exo-Jupiters exceed those for the hot Jupiters around [Fe/H]-0.7. Finally, we show that the critical metallicity for forming Jovian planets is [Fe/H]-1.2, which is evaluated by comparing the values of $$ between the hot Jupiters and the low-mass planets. The comparison intrinsically links to the different gas accretion efficiency between them., Comment: 10 pages, 3 figures, 5 tables, accepted for publication in ApJ

Polycyclic aromatic hydrocarbons (PAHs) are some of the major dust components in the interstellar medium (ISM). We present our evolution models for the abundance of PAHs in the ISM on a galaxy-evolution timescale. We consider shattering of carbonaceous dust grains in interstellar turbulence as the formation mechanism of PAHs while the PAH abundance can be reduced by coagulation onto dust grains, destruction by supernova shocks, and incorporation into stars. We implement these processes in a one-zone chemical evolution model to obtain the evolution of the PAH abundance in a galaxy. We find that PAH formation becomes accelerated above certain metallicity where shattering becomes efficient. For PAH destruction, while supernova shock is the primary mechanism in the metal-poor environment, coagulation is dominant in the metal-rich environment. We compare the evolution of the PAH abundances in our models with observed abundances in galaxies with a wide metallicity range. Our models reproduce both the paucity of PAH detection in low metallicity galaxies and the metallicity-dependence of the PAH abundance in high-metallicity galaxies. The strong metallicity dependence of PAH abundance appears as a result of the strong metallicity dependence of the dust mass increase by the accretion of metals onto dust grains, which are eventually shattered into PAHs. We conclude that the observational trend of the PAH abundance can be a natural consequence of shattering of carbonaceous grains being the source of PAHs. To establish our scenario of PAH formation, observational evidence of PAH formation by shattering would be crucial., 13 pages, 7 figures, 2 tables, accepted for publication in MNRAS

Core-collapse supernovae (SNe) are believed to be the first significant source of dust in the Universe. Such SNe are expected to be the main dust producers in young high-redshift Lyman $\alpha$ emitters (LAEs) given their young ages, providing an excellent testbed of SN dust formation models during the early stages of galaxy evolution. We focus on the dust enrichment of a specific, luminous LAE (Himiko, $z\simeq 6.6$) for which a stringent upper limit of $52.1~\mu$Jy ($3\sigma$) has recently been obtained from ALMA continuum observations at 1.2 mm. We predict its submillimetre dust emission using detailed models that follow SN dust enrichment and destruction and the equilibrium dust temperature, and obtain a plausible upper limit to the dust mass produced by a single SN: $m_\mathrm{d,SN} < 0.15$--0.45 M$_\odot$, depending on the adopted dust optical properties. These upper limits are smaller than the dust mass deduced for SN 1987A and that predicted by dust condensation theories, implying that dust produced in SNe are likely to be subject to reverse shock destruction before being injected into the interstellar medium. Finally, we provide a recipe for deriving $m_\mathrm{d,SN}$ from submillimetre observations of young, metal poor objects wherein condensation in SN ejecta is the dominant dust formation channel., Comment: 10 pages, 3 figures, accepted for publication in MNRAS

We re-examine the effect of Lyman continuum ($\lambda \leq 912$ \AA) extinction (LCE) by dust in H {\sc ii} regions in detail and discuss how it affects the estimation of the global star formation rate (SFR) of galaxies. To clarify the first issue, we establish two independent methods for estimating a parameter of LCE ($f$), which is defined as the fraction of Lyman continuum photons contributing to hydrogen ionization in an H {\sc ii} region. One of those methods determines $f$ from the set of Lyman continuum flux, electron density and metallicity. In the framework of this method, as the metallicity and/or the Lyman photon flux increase, $f$ is found to decrease. The other method determines $f$ from the ratio of infrared flux to Lyman continuum flux. Importantly, we show that $f \la 0.5$ via both methods in many H {\sc ii} regions of the Galaxy. Thus, it establishes that dust in such H {\sc ii} regions absorbs significant amount of Lyman continuum photons directly. To examine the second issue, we approximate $f$ to a function of only the dust-to-gas mass ratio (i.e., metallicity), assuming a parameter fit for the Galactic H {\sc ii} regions. We find that a characteristic $\hat{f}$, which is defined as $f$ averaged over a galaxy-wide scale, is 0.3 for the nearby spiral galaxies. This relatively small $\hat{f}$ indicates that a typical increment factor due to LCE for estimating the global SFR ($1/\hat{f}$) is large ($\sim 3$) for the nearby spiral galaxies. Therefore, we conclude that the effect of LCE is not negligible relative to other uncertainties of estimating the SFR of galaxies., Comment: 18 papges, 11 figures, accepted by ApJ

A possible formation mechanism of hydrogen molecules on a galactic scale is examined. We are interested especially in the role of hydrogen molecules for formation and evolution of primordial galaxies. Then, formation process of hydrogen molecules in a very low-metallicity galaxy (I Zw 18; the most typical metal-deficient galaxy) is studied. Adopting a recent observational result of absorption lines of hydrogen molecule in I Zw 18, we obtain the upper limit for the ionization degree in the case where hydrogen molecules can form via the H$^{-}$-process, although they are generally believed to form on the surface of dust grains. Furthermore, we present a critical ionization degree, above which H$^{-}$-process can be dominant over the formation process on the surface of grains. Interestingly, this critical ionization degree is comparable to the upper limit of the ionization degree for I Zw 18. For determination of the formation process of hydrogen molecules, future observational facilities can be useful. Thus, we examine the detectability in some wavelengths for the metal deficient galaxies. According to our estimate, near-infrared line emission of hydrogen molecule is observable in the level of 10 $\mu$Jy, the free-free radio emission is in the level of mJy, and the far-infrared emission from the dust on which hydrogen molecules form can also be detected with 10-mJy level with its temperature of 16 K. The near-infrared line and the far-infrared continuum are feasible for ASTRO-F observations., Comment: 12pages, no figure. PASJ (Accepted; tex-file is replaced)

Recently reported infrared (IR) galaxy number counts and cosmic infrared background (CIRB) all suggest that galaxies have experienced a strong evolution sometime in their lifetime. We statistically estimate the galaxy evolution history from these data. We find that an order of magnitude increase of the far-infrared (FIR) luminosity at redshift z = 0.5 - 1.0 is necessary to reproduce the very high CIRB intensity at 140 um reported by Hauser et al. (1998). z \sim 0.75 and decreases to, even at most, a factor of 10 toward z \sim 5, though many variants are allowed within these constraints. This evolution history also satisfies the constraints from the galaxy number counts obtained by IRAS, ISO and, roughly, SCUBA. The rapid evolution of the comoving IR luminosity density required from the CIRB well reproduces the very steep slope of galaxy number counts obtained by ISO. We also estimate the cosmic star formation history (SFH) from the obtained FIR luminosity density, considering the effect of the metal enrichment in galaxies. The derived SFH increases steeply with redshift in 0 < z < 0.75, and becomes flat or even declines at z > 0.75. This is consistent with the SFH estimated from the reported ultraviolet luminosity density. In addition, we present the performance of the Japanese ASTRO-F FIR galaxy survey. We show the expected number counts in the survey. We also evaluate how large a sky area is necessary to derive a secure information of galaxy evolution up to z \sim 1 from the survey, and find that at least 50 - 300 deg^2 is required., 18 pages LaTeX, PASJ in press. Abstract abridged

Infrared (IR) dust emission from galaxies is frequently used as an indicator of star formation rate (SFR). However, the effect of the dust-to-gas ratio (i.e., amount of the dust) on the conversion law from IR luminosity to SFR has not so far been considered. Then, in this paper, we present a convenient analytical formula including this effect. In order to obtain the dependence on the dust-to-gas ratio, we extend the formula derived in our previous paper, in which a theoretical formula converting IR luminosity to SFR was derived. That formula was expressed as ${\rm SFR}/(M_\odot~{\rm yr}^{-1})=\{3.3\times 10^{-10}(1- \eta)/(0.4-0.2f+0.6\epsilon)\} (L_{\rm IR}/L_\odot)$, where f is the fraction of ionizing photons absorbed by hydrogen, $\epsilon$ is the efficiency of dust absorption for nonionizing photons, $\eta$ is the cirrus fraction of observed dust luminosity, and $L_{\rm IR}$ is the observed luminosity of dust emission in the 8-1000-$\mu$m range. Our formula explains the IR excess of the Galaxy and the Large Magellanic Cloud. In the current paper, especially, we present the metallicity dependence of our conversion law between SFR and $L_{\rm IR}$. This is possible since both f and $\epsilon$ can be estimated via the dust-to-gas ratio, which is related to metallicity. We have confirmed that the relation between the metallicity and the dust-to-gas ratio is applied to both giant and dwarf galaxies. Finally, we apply the result to the cosmic star formation history. We find that the comoving SFR at z=3 calculated from previous empirical formulae is underestimated by a factor of 4-5., Comment: 8 pages LaTeX, to appear in A&A