Source author record

Hiroyuki Hirashita

Hiroyuki Hirashita appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

52works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

52 published item(s)

preprint2022arXiv

Analytic models of dust temperature in high-redshift galaxies

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}$.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

Dust diffusion in SPH simulations of an isolated galaxy

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2021arXiv

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

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.

preprint2021arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2019arXiv

Dust evolution in pre-stellar cores

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.

preprint2019arXiv

Galaxy simulation with the evolution of grain size distribution

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.

preprint2016arXiv

Dust evolution processes constrained by extinction curves in nearby galaxies

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.

preprint2016arXiv

Revisiting the lifetime estimate of large presolar grains in the interstellar medium

Some very large (>0.1 um) presolar grains are sampled in meteorites. We reconsider the lifetime of very large grains (VLGs) in the interstellar medium focusing on interstellar shattering caused by turbulence-induced large velocity dispersions. This path has never been noted as a dominant mechanism of destruction. We show that, if interstellar shattering is the main mechanism of destruction of VLGs, their lifetime is estimated to be $\gtrsim 10^8$ yr; in particular, very large SiC grains can survive cosmic-ray exposure time. However, most presolar SiC grains show residence times significantly shorter than 1 Gyr, which may indicate that there is a more efficient mechanism than shattering in destroying VLGs, or that VLGs have larger velocity dispersions than 10 km s$^{-1}$. We also argue that the enhanced lifetime of SiC relative to graphite can be the reason why we find SiC among $μ$m-sized presolar grains, while the abundance of SiC in the normal interstellar grains is much lower than graphite.

preprint2015arXiv

A new era of sub-millimeter GRB afterglow follow-ups with the Greenland Telescope

A planned rapid submillimeter (submm) Gamma Ray Burst (GRBs) follow-up observations conducted using the Greenland Telescope (GLT) is presented. The GLT is a 12-m submm telescope to be located at the top of the Greenland ice sheet, where the high-altitude and dry weather porvides excellent conditions for observations at submm wavelengths. With its combination of wavelength window and rapid responding system, the GLT will explore new insights on GRBs. Summarizing the current achievements of submm GRB follow-ups, we identify the following three scientific goals regarding GRBs: (1) systematic detection of bright submm emissions originating from reverse shock (RS) in the early afterglow phase, (2) characterization of forward shock and RS emissions by capturing their peak flux and frequencies and performing continuous monitoring, and (3) detections of GRBs as a result of the explosion of first-generation stars result of GRBs at a high redshift through systematic rapid follow ups. The light curves and spectra calculated by available theoretical models clearly show that the GLT could play a crucial role in these studies.

preprint2015arXiv

Dust Cooling in Supernova Remnants in the Large Magellanic Cloud

The infrared-to-X-ray (IRX) flux ratio traces the relative importance of dust cooling to gas cooling in astrophysical plasma such as supernova remnants (SNRs). We derive IRX ratios of SNRs in the LMC using Spitzer and Chandra SNR survey data and compare them with those of Galactic SNRs. IRX ratios of all the SNRs in the sample are found to be moderately greater than unity, indicating that dust grains are a more efficient coolant than gas although gas cooling may not be negligible. The IRX ratios of the LMC SNRs are systematically lower than those of the Galactic SNRs. As both dust cooling and gas cooling pertain to the properties of the interstellar medium, the lower IRX ratios of the LMC SNRs may reflect the characteristics of the LMC, and the lower dust-to- gas ratio (a quarter of the Galactic value) is likely to be the most significant factor. The observed IRX ratios are compared with theoretical predictions that yield IRX ratios an order of magnitude larger. This discrepancy may originate from the dearth of dust in the remnants due to either the local variation of the dust abundance in the preshock medium with respect to the canonical abundance or the dust destruction in the postshock medium. The non-equilibrium ionization cooling of hot gas, in particular for young SNRs, may also cause the discrepancy. Finally, we discuss implications for the dominant cooling mechanism of SNRs in low-metallicity galaxies.

preprint2015arXiv

Dust processing in elliptical galaxies

We reconsider the origin and processing of dust in elliptical galaxies. We theoretically formulate the evolution of grain size distribution, taking into account dust supply from asymptotic giant branch (AGB) stars and dust destruction by sputtering in the hot interstellar medium (ISM), whose temperature evolution is treated by including two cooling paths: gas emission and dust emission (i.e. gas cooling and dust cooling). With our new full treatment of grain size distribution, we confirm that dust destruction by sputtering is too efficient to explain the observed dust abundance even if AGB stars continue to supply dust grains, and that, except for the case where the initial dust-to-gas ratio in the hot gas is as high as $\sim 0.01$, dust cooling is negligible compared with gas cooling. However, we show that, contrary to previous expectations, cooling does not help to protect the dust; rather, the sputtering efficiency is raised by the gas compression as a result of cooling. We additionally consider grain growth after the gas cools down. Dust growth by the accretion of gas-phase metals in the cold medium increase the dust-to-gas ratio up to $\sim 10^{-3}$ if this process lasts >10/(n_H/10^3 cm^{-3}) Myr, where $n_\mathrm{H}$ is the number density of hydrogen nuclei. We show that the accretion of gas-phase metals is a viable mechanism of increasing the dust abundance in elliptical galaxies to a level consistent with observations, and that the steepness of observed extinction curves is better explained with grain growth by accretion.

preprint2015arXiv

The origin of dust extinction curves with or without the 2175 A bump in galaxies: The case of the Magellanic Clouds

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.

preprint2014arXiv

Constraining dust formation in high-redshift young galaxies

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 $α$ 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~μ$Jy ($3σ$) 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.

preprint2014arXiv

Dense molecular cloud cores as a source of micrometer-sized grains in galaxies

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 $ϕ_\mathrm{VLG}\sim 0.01(1-ε)/ε(τ_\mathrm{SF}/5\times 10^9~\mathrm{yr})^{-1} (f_\mathrm{VLG}/0.5)(t_\mathrm{shat}/10^8~\mathrm{yr})$, where $ε$ is the star formation efficiency in dense cores, $τ_\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 $ϕ_\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, $ϕ_\mathrm{VLG}\sim 0.09$. If we adopt an extreme value, $ϕ_\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.

preprint2014arXiv

Discovery of an Outstanding Disk in the cD Galaxy of the Hydra A Cluster

The central cD galaxy of the Hydra A cluster has one of the most powerful active galactic nuclei (AGNs) in the nearby Universe (z <0.2). We report on the discovery of a dust lane in the cD galaxy using Subaru telescope. The i'-band image shows the existence of a dark band of the size of 3.6"x 0.7" (4 kpc x 0.8 kpc), which appears to be quite similar to the dust lane observed in Centaurus A. The morphology indicates that the cold disk that seen as the dust lane is almost edge-on and rotates around the AGN. Since the minor axis of the dust lane is nearly parallel to the radio jets emerging from the AGN, the disk is probably feeding its gas into the central black hole. From the absorption, we estimate the hydrogen column density of the lane is N_H=2.0 x 10^21 cm^-2, and the mass of the disk is ~8 x 10^7 M_sun. The column density is consistent with constraints obtained from Chandra X-ray observations. The age of the disk is >~ 4 x 10^7 yr. The position angle of the disk and the galaxy's photometric axis are misaligned, which may imply that the cold gas in the disk is brought via galaxy mergers. Our observations may indicate that the supply of cold gas by galaxy mergers is required for the most intensive feedback from AGNs.

preprint2014arXiv

Effects of grain growth on the interstellar polarization curve

We apply the time evolution of grain size distributions by accretion and coagulation found in our previous work to the modelling of the wavelength dependence of interstellar linear polarization. We especially focus on the parameters of the Serkowski curve $K$ and $λ_{\max}$ characterizing the width and the maximum wavelength of this curve, respectively. We use aligned silicate and non-aligned carbonaceous spheroidal particles with different aspect ratios $a/b$. The imperfect alignment of grains with sizes larger than a cut-off size $r_{V,\rm cut}$ is considered. We find that the evolutionary effects on the polarization curve are negligible in the original model with commonly used material parameters (hydrogen number density $n_\mathrm{H}=10^3$ cm$^{-3}$, gas temperature $T_\mathrm{gas}=10$~K, and the sticking probability for accretion $S_\mathrm{acc}=0.3$). Therefore, we apply the tuned model where the coagulation threshold of silicate is removed. In this model, $λ_{\max}$ displaces to the longer wavelengths and the polarization curve becomes wider ($K$ reduces) on time-scales $\sim (30 - 50) (n_\mathrm{H}/10^3 \mathrm{cm}^{-3})^{-1}$ Myr. The tuned models at $T < 30 (n_\mathrm{H}/10^3 \mathrm{cm}^{-3})^{-1} $ Myr and different values of the parameters $r_{V,\rm cut}$ can also explain the observed trend between $K$ and $λ_{\max}$. It is significant that the evolutionary effect appears in the perpendicular direction to the effect of $r_{V,\rm cut}$ on the $K$ - $λ_{\max}$ diagram. Very narrow polarization curves can be reproduced if we change the type of particles (prolate/oblate) and/or vary $a/b$.

preprint2014arXiv

Evolution of dust content in galaxies probed by gamma-ray burst afterglows

Because of their brightness, gamma-ray burst (GRB) afterglows are viable targets for investigating the dust content in their host galaxies. Simple intrinsic spectral shapes of GRB afterglows allow us to derive the dust extinction. Recently, the extinction data of GRB afterglows are compiled up to redshift $z=6.3$, in combination with hydrogen column densities and metallicities. This data set enables us to investigate the relation between dust-to-gas ratio and metallicity out to high redshift for a wide metallicity range. By applying our evolution models of dust content in galaxies, we find that the dust-to-gas ratio derived from GRB afterglow extinction data are excessively high such that they can be explained with a fraction of gas-phase metals condensed into dust ($f_\mathrm{in}$) $\sim 1$, while theoretical calculations on dust formation in the wind of asymptotic giant branch stars and in the ejecta of Type II supernovae suggest a much more moderate condensation efficiency ($f_\mathrm{in}\sim 0.1$). Efficient dust growth in dense clouds has difficulty in explaining the excessive dust-to-gas ratio at metallicities $Z/\mathrm{Z}_\odot <ε$, where $ε$ is the star formation efficiency of the dense clouds. However, if $ε$ is as small as 0.01, the dust-to-gas ratio at $Z\sim 10^{-2}$ Z$_\odot$ can be explained with $n_\mathrm{H}\gtrsim 10^6$ cm$^{-3}$. Therefore, a dense environment hosting dust growth is required to explain the large fraction of metals condensed into dust, but such clouds should have low star formation efficiencies to avoid rapid metal enrichment by stars.

preprint2014arXiv

Evolution of grain size distribution in high-redshift dusty quasars: Integrating large amounts of dust and unusual extinction curves

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.

preprint2014arXiv

Evolution of the extinction curves in galaxies

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\;μ$m). Then, grain growth becomes effective at small grain radii, forming a bump at $a \sim 10^{-3}${--}$10^{-2}\;μ$m on the grain size distribution. Consequently, the extinction curve at ultraviolet (UV) wavelengths becomes steep, and a bump at $1/λ\sim 4.5\;μ{\rm m}^{-1}\;(λ: \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.

preprint2014arXiv

Formation History of Polycyclic Aromatic Hydrocarbons in Galaxies

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.

preprint2014arXiv

Molecular hydrogen emission in the interstellar medium of the Large Magellanic Cloud

We present the detection and analysis of molecular hydrogen emission toward ten interstellar regions in the Large Magellanic Cloud. We examined low-resolution infrared spectral maps of twelve regions obtained with the Spitzer infrared spectrograph (IRS). The pure rotational 0--0 transitions of H$_2$ at 28.2 and 17.1${\,\rm μm}$ are detected in the IRS spectra for ten regions. The higher level transitions are mostly upper limit measurements except for three regions, where a 3$σ$ detection threshold is achieved for lines at 12.2 and 8.6${\,\rm μm}$. The excitation diagrams of the detected H$_2$ transitions are used to determine the warm H$_2$ gas column density and temperature. The single-temperature fits through the lower transition lines give temperatures in the range $86-137\,{\rm K}$. The bulk of the excited H$_2$ gas is found at these temperatures and contributes $\sim$5-17% to the total gas mass. We find a tight correlation of the H$_2$ surface brightness with polycyclic aromatic hydrocarbon and total infrared emission, which is a clear indication of photo-electric heating in photodissociation regions. We find the excitation of H$_2$ by this process is equally efficient in both atomic and molecular dominated regions. We also present the correlation of the warm H$_2$ physical conditions with dust properties. The warm H$_2$ mass fraction and excitation temperature show positive correlations with the average starlight intensity, again supporting H$_2$ excitation in photodissociation regions.

preprint2014arXiv

Planet Traps and First Planets: the Critical Metallicity for Gas Giant Formation

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 ($<R_{rapid}>$) 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 $<R_{rapid}>$ 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 $<R_{rapid}>$ between the hot Jupiters and the low-mass planets. The comparison intrinsically links to the different gas accretion efficiency between them.

preprint2014arXiv

The dust production rate of AGB stars in the Magellanic Clouds

We compare theoretical dust yields for stars with mass 1 Msun < mstar < 8 Msun, and metallicities 0.001 < Z < 0.008 with observed dust production rates (DPR) by carbon- rich and oxygen-rich Asymptotic Giant Branch (C-AGB and O-AGB) stars in the Large and Small Magellanic Clouds (LMC, SMC). The measured DPR of C-AGB in the LMC are reproduced only if the mass loss from AGB stars is very efficient during the carbon-star stage. The same yields over-predict the observed DPR in the SMC, suggesting a stronger metallicity dependence of the mass-loss rates during the carbon- star stage. DPR of O-AGB stars suggest that rapid silicate dust enrichment occurs due to efficient hot-bottom-burning if mstar > 3 Msun and Z > 0.001. When compared to the most recent observations, our models support a stellar origin for the existing dust mass, if no significant destruction in the ISM occurs, with a contribution from AGB stars of 70% in the LMC and 15% in the SMC.

preprint2014arXiv

Two-size approximation: a simple way of treating the evolution of grain size distribution in galaxies

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.

preprint2013arXiv

Condition for the formation of micron-sized dust grains in dense molecular cloud cores

We investigate the condition for the formation of micron-sized grains in dense cores of molecular clouds. This is motivated by the detection of the mid-infrared emission from deep inside a number of dense cores, the so-called `coreshine,' which is thought to come from scattering by micron-sized grains. Based on numerical calculations of coagulation starting from the typical grain size distribution in the diffuse interstellar medium, we obtain a conservative lower limit to the time $t$ to form micron-sized grains: $t/t_\mathrm{ff}>3 (5/S) (n_\mathrm{H}/10^5 \mathrm{cm}^{-3})^{-1/4}$ (where $t_\mathrm{ff}$ is the free-fall time at hydrogen number density $n_\mathrm{H}$ in the core, and $S$ the enhancement factor to the grain-grain collision cross-section to account for non-compact aggregates). At the typical core density $n_\mathrm{H}=10^5 \mathrm{cm}^{-3}$, it takes at least a few free-fall times to form the micron-sized grains responsible for coreshine. The implication is that those dense cores observed in coreshine are relatively long-lived entities in molecular clouds, rather than dynamically transient objects that last for one free-fall time or less.

preprint2013arXiv

Effects of grain growth mechanisms on the extinction curve and the metal depletion in the interstellar medium

Dust grains grow their sizes in the interstellar clouds (especially in molecular clouds) by accretion and coagulation. Here we model and test these processes by examining the consistency with the observed variation of the extinction curves in the Milky Way. We find that, if we simply use the parameters used in previous studies, the model fails to explain the flattening of far-UV extinction curve for large $R_V$ (flatness of optical extinction curve) and the existence of carbon bump even in flat extinction curves. This discrepancy is resolved by adopting a `tuned' model, in which coagulation of carbonaceous dust is less efficient (by a factor of 2) and that of silicate is more efficient with the coagulation threshold removed. The tuned model is also consistent with the relation between silicon depletion (indicator of accretion) and $R_V$ if the duration of accretion and coagulation is >100(n_H/10^3 cm^{-3})^{-1} Myr, where n_H is the number density of hydrogen nuclei in the cloud. We also examine the relations between each of the extinction curve features (UV slope, far-UV curvature, and carbon bump strength) and $R_V$. The correlation between UV slope and $R_V$, which is the strongest among the three correlations, is well reproduced by the tuned model. For far-UV curvature and carbon bump strength, the observational data are located between the tuned model and the original model without tuning, implying that the large scatters in the observational data can be explained by the sensitive response to the coagulation efficiency. The overall success of the tuned model indicates that accretion and coagulation are promising mechanisms of producing the variation of extinction curves in the Milky Way, although we do not exclude possibilities of other dust-processing mechanisms changing extinction curves.

preprint2013arXiv

Evolution of dust grain size distribution by shattering in the interstellar medium: robustness and uncertainty

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.

preprint2013arXiv

What determines the grain size distribution in galaxies?

We construct a dust evolution model taking into account the grain size distribution, and investigate what kind of dust processes determine the grain size distribution at each stage of galaxy evolution. In addition to the dust production by type II supernovae (SNeII) and asymptotic giant branch (AGB) stars, we consider three processes in the ISM: (i) dust destruction by SN shocks, (ii) metal accretion onto the surface of preexisting grains in the cold neutral medium (CNM) (called grain growth), and (iii) grain-grain collisions (shattering and coagulation) in the warm neutral medium (WNM) and CNM. We found that the grain size distribution in galaxies is controlled by stellar sources in the early stage of galaxy evolution, and that afterwards the main processes that govern the size distribution changes to those in the ISM. Since shattering produces a large abundance of small grains (consequently, the surface-to-volume ratio of grains increases), it enhances the efficiency of grain growth, contributing to the significant increase of the total dust mass. Grain growth creates a large bump in the grain size distribution around a ~ 0.01 μm. Coagulation occurs effectively after the number of small grains is enhanced by shattering, and the grain size distribution is deformed to have a bump at a ~ 0.03 - 0.05 μm at t ~ 10 Gyr. We conclude that the evolutions of the total dust mass and the grain size distribution in galaxies are closely related to each other, and the grain size distribution changes considerably through the galaxy evolution because the dominant dust processes which regulate the grain size distribution change.

preprint2012arXiv

Dust formation history of galaxies: a critical role of metallicity for the dust mass growth by accreting materials in the interstellar medium

This paper investigate what is the main driver of the dust mass growth in the interstellar medium (ISM) by using a chemical evolution model of galaxy with metals (elements heavier than helium) in dust phase in addition to the total amount of metals. We consider asymptotic giant branch (AGB) stars, type II supernovae (SNe II) and the dust mass growth in the ISM as the sources of dust, and SN shocks as the destruction mechanism of dust. Further, to describe the dust evolution precisely, our model takes into account the age and metallicity (the ratio of metal mass to ISM mass) dependence of the sources of dust. We particularly focused on the dust mass growth, and found that the dust mass growth in the ISM is regulated by the metallicity. To quantify this aspect, we introduce a "critical metallicity", which is a metallicity at which the contribution of stars (AGB stars and SNe II) equals that of the dust mass growth in the ISM. If the star formation timescale is shorter, the value of the critical metallicity is higher, but the galactic age at which the metallicity reaches the critical metallicity is shorter. From observations, it was expected that the dust mass growth was the dominant source of dust in the Milky Way and dusty QSOs at high redshifts. By introducing the critical metallicity, it is clearly shown that the dust mass growth is the main source of dust in such galaxies with various star formation timescales and ages. The dust mass growth in the ISM is regulated by metallicity, and we stress that the critical metallicity works as an indicator to judge whether the grain growth in the ISM is dominant source of dust in a galaxy, especially because of the strong and nonlinear dependence on the metallicity.

preprint2012arXiv

Dust growth in the interstellar medium: How do accretion and coagulation interplay?

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 $a\la 0.001 \micron$ on a time-scale of $\la 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 \textit{enhanced} by accretion, while they are flattened by coagulation.

preprint2012arXiv

Impact of grain size distributions on the dust enrichment in high-redshift quasars

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.

preprint2012arXiv

Properties of free-free, dust, and CO emissions in the starbursts of blue compact dwarf galaxies

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.

preprint2012arXiv

Synthesized grain size distribution in the interstellar medium

We examine a synthetic way of constructing the grain size distribution in the interstellar medium (ISM). First we formulate a synthetic grain size distribution composed of three grain size distributions processed with the following mechanisms that govern the grain size distribution in the Milky Way: (i) grain growth by accretion and coagulation in dense clouds, (ii) supernova shock destruction by sputtering in diffuse ISM, and (iii) shattering driven by turbulence in diffuse ISM. Then, we examine if the observational grain size distribution in the Milky Way (called MRN) is successfully synthesized or not. We find that the three components actually synthesize the MRN grain size distribution in the sense that the deficiency of small grains by (i) and (ii) is compensated by the production of small grains by (iii). The fraction of each {contribution} to the total grain processing of (i), (ii), and (iii) (i.e., the relative importance of the three {contributions} to all grain processing mechanisms) is 30-50%, 20-40%, and 10-40%, respectively. We also show that the Milky Way extinction curve is reproduced with the synthetic grain size distributions.

preprint2011arXiv

Central free-free-dominated 880 $\micron$ emission in II Zw 40

The central star-forming region in a blue compact dwarf galaxy, II Zw 40, was observed in the 340 GHz ($880 \micron$) band at $\sim 5$ arcsec (250 pc) resolution with the Submillimetre Array (SMA). A source associated with the central star-forming complex was detected with a flux of $13.6\pm 2.0$ mJy. The structure is more extended than the beam in the east-west direction. The SMA 880 $\micron$ flux is analyzed by using theoretical models of radio spectral energy distribution along with centimetre interferometric measurements in the literature. We find (i) that the SMA 880 $\micron$ flux is dominated ($\sim 75$ per cent) by free-free emission from the central compact star-forming region, and (ii) that the contribution from dust emission to the SMA 880 $\micron$ flux is at most $4\pm 2.5$ mJy. We also utilize our models to derive the radio--FIR relation of the II Zw 40 centre, suggesting that free-free absorption at low frequencies ($ν\la$ several GHz; $λ\ga$ several cm) and spatial extent of dust affect the radio-FIR relation.

preprint2011arXiv

Dust properties in afterglow of GRB 071025 at z~5

At high redshift, the universe was so young that core-collapse supernovae (SNe) are suspected to be the dominant source of dust production. However, some observations indicate that the dust production by SNe is an inefficient process, casting doubts on the existence of abundant SNe-dust in the early universe. Recently, Perley et al. (2010) reported that the afterglow of GRB 071025 - an unusually red GRB at z ~ 5 - shows evidence for the SNe-produced dust. Since this is perhaps the only high redshift GRB exhibiting compelling evidence for SNe-dust but the result could easily be affected by small systematics in photometry, we re-examined the extinction properties of GRB 071025 using our own optical/near-infrared data at a different epoch. In addition, we tested SNe-dust models with different progenitor masses and dust destruction efficiencies to constrain the dust formation mechanisms. By searching for the best-fit model of the afterglow spectral energy distribution, we confirm the previous claim that the dust in GRB 071025 is most likely to originate from SNe. We also find that the SNe-dust model of 13 or 25M_{\odot} without dust destruction fits the extinction property of GRB 071025 best, while pair-instability SNe (PISNe) models with a 170M_{\odot} progenitor poorly fit the data. Our results indicate that, at least in some systems at high redshift, SNe with intermediate masses within 10 - 30M_{\odot} were the main contributors for the dust enrichment, and the dust destruction effect due to reverse shock was negligible.

preprint2011arXiv

Effects of grain size distribution on the interstellar dust mass growth

Grain growth by the accretion of metals in interstellar clouds (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 demonstrated in this paper. First, we derive an analytical formula that gives the grain size distribution 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 grain 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 clouds governs the dust content in nearby galaxies unless the grain size is strongly biased to sizes larger than $\sim 0.1 \micron$ or the power index of the grain size distribution is shallower than $\sim -2.5$. The grain growth in clouds contributes to the rapid increase of dust-to-gas ratio at a certain metallicity level (called critical metallicity in Asano et al. 2011 and and Inoue 2011), which we find to be sensitive to grain size distribution. Thus, the grain growth efficiently 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 applicable for any grain size distribution and easily implemented into any framework of dust enrichment in galaxies.

preprint2011arXiv

The Role of Dust in the Early Universe I : Protogalaxy Evolution

We develop one-zone galaxy formation models in the early Universe, taking into account dust formation and evolution by supernova (SN) explosions. Especially we focus on the time evolution of dust size distribution, because ${\rm H}_{2}$ formation on the dust surface plays a critical role in the star formation process in the early Universe. In the model we assume that star formation rate (SFR) is proportional to the total amount of ${\rm H}_{2}$. We consistently treat (i) the formation and size evolution of dust, (ii) the chemical reaction networks including ${\rm H}_{2}$ formation both on the surface of dust and in gas phase, and (iii) the SFR in the model. First, we find that, because of dust destruction due to both reverse and forward shocks driven by SNe, H$_{2}$ formation is more suppressed than that without dust destruction. At the galaxy age of $sim0.8 {\rm Gyr}$, for galaxy models with virial mass $M_{\rm vir}=10^{9} M_{\odot}$ and formation redshift $z_{\rm vir}=10$, the molecular fraction is 2.5 orders of magnitude less in the model with dust destruction by both shocks than that in the model without dust destruction. Second, we show that the H$_{2}$ formation rate strongly depends on the ISM density around SN progenitors. The SFR in higher ISM density is lower, since dust destruction by reverse shocks is more effective in higher ISM density. We conclude that not only the amount but also the size distribution of dust being related with the star formation activity strongly affects the evolution of galaxies in the early Universe.

preprint2010arXiv

Effects of grain shattering by turbulence on extinction curves in starburst galaxies

Dust grains can be efficiently accelerated and shattered in warm ionized medium (WIM) because of the turbulent motion. This effect is enhanced in starburst galaxies, where gas is ionized and turbulence is sustained by massive stars. Moreover, dust production by Type II supernovae (SNe II) can be efficient in starburst galaxies. In this paper, we examine the effect of shattering in WIM on the dust grains produced by SNe II. We find that although the grains ejected from SNe II are expected to be biased to large sizes ($a\ga 0.1 \micron$, where $a$ is the grain radius) because of the shock destruction in supernova remnants, the shattering in WIM is efficient enough in $\sim 5$ Myr to produce small grains if the metallicity is nearly solar or more. The production of small grains by shattering steepens the extinction curve. Thus, steepening of extinction curves by shattering should always be taken into account for the system where the metallicity is solar and the starburst age is typically larger than 5 Myr. These conditions may be satisfied not only in nearby starbursts but also in high redshift ($z>5$) quasars.

preprint2010arXiv

Importance of far-infrared mapping in a spiral galaxy: AKARI observation of M81

The importance of the far-infrared (FIR) mapping is demonstrated for a face-on spiral galaxy, M81, by analyzing its imaging data at 65, 90, and 140 μm taken by AKARI. 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 μm 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 diagramis 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).

preprint2010arXiv

Shattering by turbulence as a production source of very small grains

The origin of grain size distribution in the interstellar medium is one of the most fundamental problems in the interstellar physics. In the Milky Way, smaller grains are more abundant in number, but their origins are not necessarily specified and quantified. One of the most efficient drivers of small grain production is interstellar turbulence, in which dust grains can acquire relative velocities large enough to be shattered. Applying the framework of shattering developed in previous papers, we show that small ($a\la 0.01~\micron$) grains reach the abundance level observed in the Milky Way in $\sim 10^8$ yr (i.e. within the grain lifetime) by shattering in warm neutral medium. We also show that if part of grains experience additional shattering in warm ionized medium, carbonaceous grains with $a\sim 0.01~\micron$ are redistributed into smaller sizes. This could explain the relative enhancement of very small carbonaceous grains with $a\sim 3$--100 Å. Our theory also explains the ubiquitous association between large grains and very small grains naturally. Some tests for our theory are proposed in terms of the metallicity dependence.

preprint2009arXiv

Detecting Lyman Alpha Emitters in the Sub-millimeter

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_α> 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.

preprint2009arXiv

Dust coagulation in star formation with different metallicities

Dust grains coagulate into larger aggregates in dense gas. This changes their size distribution and possibly affects the thermal evolution of star-forming clouds. We here investigate dust coagulation in collapsing pre-stellar cores with different metallicities by considering the thermal motions of grains. We show that coagulation does occur even at low metallicity $\sim 10^{-6} \Zsun$. However, we also find (i) that the H$_2$ formation rate on dust grains is reduced only after the majority of H$_2$ is formed; and (ii) that the dust opacity is modified only after the core becomes optically thick. Therefore, we conclude that the effects of dust coagulation can safely be neglected in discussing the temperature evolution of the pre-stellar cores for any metallicity as long as the grain motions are thermal.

preprint2006arXiv

Common Correlations between 60, 100 and 140 um Intensities in the Galactic Plane and Magellanic Clouds

We investigate the far-infrared SED of the Galaxy and the Magellanic Clouds by using the COBE (COsmic Background Explorer) / DIRBE (Diffuse InfraRed Background Experiment) ZSMA (Zodi - Subtracted Mission Average) maps at wavelengths of 60 um, 100 um and 140 um. We analyze three regions: the Galactic plane region with the Galactic latitude |b|<5 degree, the Large Magellanic Cloud (LMC) region, and the Small Magellanic Cloud (SMC) region. Because the dust optical depth is much smaller in the far-infrared than in the visible, we may observe cumulative far-infrared radiation from regions with various interstellar radiation field (IRSF) in a line of sight. As consequence of considering such an effect, we find a common far-infrared color correlation between the 140 - 100 um and 60 - 100 um intensity ratios in all the three galaxies. Although this color correlation cannot be explained by any existing model, it fits very well the far-infrared color of nearby star forming galaxies.

preprint2001arXiv

Absorption Measurement of Hydrogen Molecules in the Early Universe

We investigate the observability of hydrogen molecules in absorption. The absorption efficiency of the hydrogen molecules become comparable with or larger than that of the dust grains in the metal-poor condition expected in the early Universe. If we can use bright infrared continuum sources behind the molecular gas clouds, the absorption measurement of the hydrogen molecules will be an important technique to explore the primordial gas clouds that are contracting into first-generation objects.

preprint2001arXiv

H2 Formation in Low Metallicity Galaxies

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 $μ$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.

preprint2000arXiv

Exploring Galaxy Evolution from Infrared Number Counts and Cosmic Infrared Background

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.