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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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

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.

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

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

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

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.

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

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.

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.

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.

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.

We present a theoretical model aimed at explaining the IRX-$\beta$ 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.

Grain growth during star formation affects the physical and chemical processes in the evolution of star-forming clouds. We investigate the origin of the millimeter (mm)-sized grains recently observed in Class I protostellar envelopes. We use the coagulation model developed in our previous paper and find that a hydrogen number density of as high as $10^{10}~{\rm cm^{-3}}$, instead of the typical density $10^5~{\rm cm^{-3}}$, is necessary for the formation of mm-sized grains. Thus, we test a hypothesis that such large grains are transported to the envelope from the inner, denser parts, finding that gas drag by outflow efficiently "launches" the large grains as long as the central object has not grown to $\gtrsim 0.1$ M$_{\odot}$. By investigating the shattering effect on the mm-sized grains, we ensure that the large grains are not significantly fragmented after being injected in the envelope. We conclude that the mm-sized grains observed in the protostellar envelopes are not formed in the envelopes but formed in the inner parts of the star-forming regions and transported to the envelopes before a significant mass growth of the central object, and that they survive in the envelopes.

Formation and evolution of galaxies have been a central driving force in the studies of galaxies and cosmology. Recent studies provided a global picture of cosmic star formation history. However, what drives the evolution of star formation activities in galaxies has long been a matter of debate. The key factor of the star formation is the transition of hydrogen from atomic to molecular state, since the star formation is associated with the molecular phase. This transition is also strongly coupled with chemical evolution, because dust grains, i.e., tiny solid particles of heavy elements, play a critical role in molecular formation. Therefore, a comprehensive understanding of neutral-molecular gas transition, star formation and chemical enrichment is necessary to clarify the galaxy formation and evolution. Here we present the activity of SKA-JP galaxy evolution sub-science working group (subSWG) Our activity is focused on three epochs: z \sim 0, 1, and z > 3. At z \sim 0, we try to construct a unified picture of atomic and molecular hydrogen through nearby galaxies in terms of metallicity and other various ISM properties. Up to intermediate redshifts z \sim 1, we explore scaling relations including gas and star formation properties, like the main sequence and the Kennicutt-Schmidt law of star forming galaxies. To connect the global studies with spatially-resolved investigations, such relations will be plausibly a viable way. For high redshift objects, the absorption lines of HI 21-cm line will be a very promising observable to explore the properties of gas in galaxies. By these studies, we will surely witness a real revolution in the studies of galaxies by SKA.

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 $\mu$m-sized presolar grains, while the abundance of SiC in the normal interstellar grains is much lower than graphite.

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.

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.

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.

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.

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.

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.

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 \mu 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$\sigma$ detection threshold is achieved for lines at 12.2 and 8.6${\,\rm \mu 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.

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.

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 $\lambda_{\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, $\lambda_{\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 $\lambda_{\max}$. It is significant that the evolutionary effect appears in the perpendicular direction to the effect of $r_{V,\rm cut}$ on the $K$ - $\lambda_{\max}$ diagram. Very narrow polarization curves can be reproduced if we change the type of particles (prolate/oblate) and/or vary $a/b$.

A 12-m diameter radio telescope will be deployed to the Summit Station in Greenland to provide direct confirmation of a Super Massive Black Hole (SMBH) by observing its shadow image in the active galaxy M87. The telescope (Greenland Telescope: GLT) is to become one of the Very Long Baseline Interferometry (VLBI) stations at sub-millimeter (submm) regime, providing the longest baseline > 9,000 km to achieve an exceptional angular resolution of 20 micro arc sec at 350 GHz, which will enable us to resolve the shadow size of ~40 micro arc sec. The triangle with the longest baselines formed by the GLT, the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile, and the Submillimeter Array (SMA) in Hawaii will play a key role for the M87 observations. We have been working on the image simulations based on realistic conditions for a better understanding of the possible observed images. In parallel, retrofitting of the telescope and the site developments are in progress. Based on three years of opacity monitoring at 225 GHz, our measurements indicate that the site is excellent for submm observations, comparable to the ALMA site. The GLT is also expected to make single-dish observations up to 1.5 THz.

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.

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.

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 <\epsilon$, where $\epsilon$ is the star formation efficiency of the dense clouds. However, if $\epsilon$ 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.

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.

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.

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.

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.

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.

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.

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.

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 \mu 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 \mu 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.

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.

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.

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.

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.

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.

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.

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.

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 ($\nu\la$ several GHz; $\lambda\ga$ several cm) and spatial extent of dust affect the radio-FIR relation.

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.

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.

We apply the supernova(SN) extinction curves to reproduce the observed properties of SST J1604+4304 which is a young infrared (IR) galaxy at z = 1. The SN extinction curves used in this work were obtained from models of unmixed ejecta of type II supernovae(SNe II) for the Salpeter initial mass function (IMF) with a mass range from 8 to 30 M_sun or 8 to 40 M_sun. The effect of dust distributions on the attenuation of starlight is investigated by performing the chi-square fitting method against various dust distributions. These are the commonly used uniform dust screen, the clumpy dust screen, and the internal dust geometry. We add to these geometries three scattering properties, namely, no-scattering, isotropic scattering, and forward-only scattering. Judging from the chi-square values, we find that the uniform screen models with any scattering property provide good approximations to the real dust geometry. Internal dust is inefficient to attenuate starlight and thus cannot be the dominant source of the extinction. We show that the SN extinction curves reproduce the data of SST J1604+4304 comparable to or better than the Calzetti extinction curve. The Milky Way extinction curve is not in satisfactory agreement with the data unless several dusty clumps are in the line of sight. This trend may be explained by the abundance of SN-origin dust in these galaxies; SN dust is the most abundant in the young IR galaxy at z = 1, abundant in local starbursts, and less abundant in the Galaxy. If dust in SST J1604+4304 is dominated by SN dust, the dust production rate is about 0.1 M_sun per SN.

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 {\mu}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 {\mu}m and 140 {\mu}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 {\mu}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).

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 \AA. 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.

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.

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 in extragalactic HII regions, with the aim of understanding the role of dust and different physical conditions in the ionized medium. First, we compiled several observational data sets for Galactic and extragalactic HII regions and confirm that extragalactic HII regions follow the same size (D)--density (n) relation as Galactic ones. Motivated by the inability of static models to explain this, we then modelled the evolution of the size--density relation of HII regions by considering their star formation history, the effects of dust, and pressure-driven expansion. The results are compared with our sample data whose size and density span roughly six orders of magnitude. We find that the extragalactic size--density relation does not result from an evolutionary sequence of HII regions but rather reflects a sequence with different initial gas densities (``density hierarchy''). Moreover, the size of many HII regions is limited by dust absorption of ionizing photons, rather than consumption by ionizing neutral hydrogen. Dust extinction of ionizing photons is particularly severe over the entire lifetime of compact HII regions with typical gas densities of greater than 10^3 cm^{-3}. Hence, as long as the number of ionizing photons is used to trace massive star formation, much star-formation activity could be missed. Such compact dense environments, the ones most profoundly obscured by dust, have properties similar to ``maximum--intensity starbursts''. This implies that submillimeter and infrared wavelengths may be necessary to accurately assess star formation in these extreme conditions both locally and at high redshift.

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.

We investigate the FIR properties of a sample of BCDs observed by AKARI. By utilizing the data at wavelengths of $\lambda =65 \mu$m, 90 $\mu$m, and 140 $\mu$m, we find that the FIR colours of the BCDs are located at the natural high-temperature extension of those of the Milky Way and the Magellanic Clouds. This implies that the optical properties of dust in BCDs are similar to those in the Milky Way. Indeed, we explain the FIR colours by assuming the same grain optical properties, which may be appropriate for amorphous dust grains, and the same size distribution as those adopted for the Milky Way dust. Since both interstellar radiation field and dust optical depth affect the dust temperature, it is difficult to distinguish which of these two physical properties is responsible for the change of FIR colours. Then, in order to examine if the dust optical depth plays an important role in determining the dust temperature, we investigate the correlation between FIR colour (dust temperature) and dust-to-gas ratio. We find that the dust temperature tends to be high as the dust-to-gas ratio decreases but that this trend cannot be explained by the effect of dust optical depth. Rather, it indicates a correlation between dust-to-gas ratio and interstellar radiation field. Although the metallicity may also play a role in this correlation, we suggest that the dust optical depth could regulate the star formation activities, which govern the interstellar radiation field. We also mention the importance of submillimetre data in tracing the emission from highly shielded low-temperature dust.

We present recent advances in theoretical studies of the formation and evolution of dust in primordial supernovae (SNe) that are considered to be the main sources of dust in the early universe. Being combined with the results of calculations of dust formation in the ejecta of Population III SNe, the investigations of the evolution of newly formed dust within supernova remnants (SNRs) show that smaller grains are predominantly destroyed by sputtering in the shocked gas, while larger grains are injected into the ambient medium. The mass of dust grains surviving the destruction in SNRs reaches up to 0.1--15 $M_\odot$, which is high enough to account for the content of dust observed for the host galaxies of quasars at $z > 5$. In addition, the transport of dust formed in the ejecta causes the formation of low-mass stars in the dense shell of primordial SNRs and affects the elemental composition of those stars. We also show that the flat extinction curve is expected in the high-redshift universe where SNe are the possible sources of dust.

We investigate shattering and coagulation of dust grains in turbulent interstellar medium (ISM). The typical velocity of dust grain as a function of grain size has been calculated for various ISM phases based on a theory of grain dynamics in compressible magnetohydrodynamic turbulence. In this paper, we develop a scheme of grain shattering and coagulation and apply it to turbulent ISM by using the grain velocities predicted by the above turbulence theory. Since large grains tend to acquire large velocity dispersions as shown by earlier studies, large grains tend to be shattered. Large shattering effects are indeed seen in warm ionized medium (WIM) within a few Myr for grains with radius $a\ga 10^{-6}$ cm. We also show that shattering in warm neutral medium (WNM) can limit the largest grain size in ISM ($a\sim 2\times 10^{-5} \mathrm{cm}$). On the other hand, coagulation tends to modify small grains since it only occurs when the grain velocity is small enough. Coagulation significantly modifies the grain size distribution in dense clouds (DC), where a large fraction of the grains with $a<10^{-6}$ cm coagulate in 10 Myr. In fact, the correlation among $R_V$, the carbon bump strength, and the ultraviolet slope in the observed Milky Way extinction curves can be explained by the coagulation in DC. It is possible that the grain size distribution in the Milky Way is determined by a combination of all the above effects of shattering and coagulation. Considering that shattering and coagulation in turbulence are effective if dust-to-gas ratio is typically more than $\sim 1/10$ of the Galactic value, the regulation mechanism of grain size distribution should be different between metal-poor and metal-rich environments.

We report basic far-infrared (FIR) properties of eight blue compact dwarf galaxies (BCDs) observed by AKARI. We measure the fluxes at the four FIS bands (wavelengths of 65 um, 90 um, 140 um, and 160 um). Based on these fluxes, we estimate basic quantities about dust: dust temperature, dust mass, and total FIR luminosity. We find that the typical dust temperature of the BCD sample is systematically higher than that of normal spiral galaxies, although there is a large variety. The interstellar radiation field estimated from the dust temperature ranges up to 100 times of the Galactic value. This confirms the concentrated star-forming activity in BCDs. The star formation rate can be evaluated from the FIR luminosity as 0.01--0.5 $M_\odot$ yr$^{-1}$. Combining this quantity with gas mass taken from the literature, we estimate the gas consumption timescales (gas mass divided by the star formation rate), which prove to span a wide range from 1 Gyr to 100 Gyr. A natural interpretation of this large variety can be provided by intermittent star formation activity. We finally show the relation between dust-to-gas ratio and metallicity (we utilize our estimate of dust mass, and take other necessary quantities from the literature). There is a positive correlation between dust-to-gas ratio and metallicity as expected from chemical evolution models.

We investigate the extinction curves of young galaxies in which dust is supplied from Type II supernovae (SNe II) and/or pair instability supernovae (PISNe). Since at high redshift (z>5), low-mass stars cannot be dominant sources for dust grains, SNe II and PISNe, whose progenitors are massive stars with short lifetimes, should govern the dust production. Here, we theoretically investigate the extinction curves of dust produced by SNe II and PISNe, taking into account reverse shock destruction induced by collision with ambient interstellar medium. We find that the extinction curve is sensitive to the ambient gas density around a SN, since the efficiency of reverse shock destruction strongly depends on it. The destruction is particularly efficient for small-sized grains, leading to a flat extinction curve in the optical and ultraviolet wavelengths. Such a large ambient density as n_H > 1 cm^{-3} produces too flat an extinction curve to be consistent with the observed extinction curve for SDSS J104845.05+463718.3 at z=6.2. Although the extinction curve is highly sensitive to the ambient density, the hypothesis that the dust is predominantly formed by SNe at z~6 is still allowed by the current observational constraints. For further quantification, the ambient density should be obtained by some other methods. Finally we also discuss the importance of our results for observations of high-z galaxies, stressing a possibility of flat extinction curves.

Aims:In some of the lensed quasars, color differences between multiple images are observed at optical/near-infrared wavelengths. There are three possible origins of the color differences: intrinsic variabilities of quasars, differential dust extinction, and quasar microlensing. We examine how these three possible scenarios can reproduce the observed chromaticity. Methods:We evaluate how much color difference between multiple images can be reproduced by the above three possible scenarios with realistic models; (i) an empirical relation for intrinsic variabilities of quasars, (ii) empirical relations for dust extinction and theoretically predicted inhomogeneity in galaxies, or (iii) a theoretical model for quasar accretion disks and magnification patterns in the vicinity of caustics. Results:We find that intrinsic variabilities of quasars cannot be a dominant source responsible for observed chromatic features in multiple quasars. In contrast, either dust extinction or quasar microlensing can nicely reproduce the observed color differences between multiple images in most of the lensed quasars. Taking into account the time interval between observations at different wavebands in our estimations, quasar microlensing is a more realistic scenario to reproduce the observed color differences than dust extinction. All the observed color differences presented in this paper can be explained by a combination of these two effects, but monitoring observations at multiple wavebands are necessary to disentangle these.

A recent data analysis of the far-infrared (FIR) map of the Galaxy and the Magellanic Clouds has shown that there is a tight correlation between two FIR colours: the 60 um-100 um and 100 um-140 um colours. This FIR colour relation called ``main correlation'' can be interpreted as indicative of a sequence of various interstellar radiation fields with a common FIR optical property of grains. In this paper, we constrain the FIR optical properties of grains by comparing the calculated FIR colours with the observational main correlation. We show that neither of the ``standard'' grain species (i.e. astronomical silicate and graphite grains) reproduces the main correlation. However, if the emissivity index at ~ 100--200 um is changed to ~ 1--1.5 (not ~ 2 as the above two species), the main correlation can be successfully explained. Thus, we propose that the FIR emissivity index is ~ 1--1.5 for the dust in the Galaxy and the Magellanic Clouds at ~ 100--200 um. We also consider the origin of the minor correlation called ``sub-correlation'', which can be used to estimate the Galactic star formation rate.

We investigate the radio spectral energy distributions (SEDs) of young star-forming galaxies and how they evolve with time. The duration and luminosity of the nonthermal radio emission from supernova remnants (SNRs) are constrained by using the observational radio SEDs of SBS 0335-052 and I Zw 18, which are the two lowest-metallicity blue compact dwarf galaxies in the nearby universe. The typical radio ``fluence'' for SNRs in SBS 0335-052, that is the radio energy emitted per SNR over its radiative lifetime, is estimated to be $\sim 6$--$22\times 10^{22} {\rm W Hz^{-1} yr}$ at 5 GHz. On the other hand, the radio fluence in I Zw 18 is $\sim 1$--$3\times 10^{22} {\rm W Hz^{-1} yr}$ at 5 GHz. We discuss the origin of this variation and propose scaling relations between synchrotron luminosity and gas density. We have also predicted the time dependence of the radio spectral index and of the spectrum itself, for both the ``active'' (SBS 0335-052) and ``passive'' (I Zw 18) cases. These models enable us to roughly age date and classify radio spectra of star-forming galaxies into active/passive classes. Implications for high-z galaxy evolution are also discussed.

Recent observations have revealed that damped Ly$\alpha$ clouds (DLAs) host star formation activity. In order to examine if such star formation activity can be triggered by ionization fronts, we perform high-resolution hydrodynamics and radiative transfer simulations of the effect of radiative feedback from propagating ionization fronts on high-density clumps. We examine two sources of ultraviolet (UV) radiation field to which high-redshift (z ~ 3) galaxies could be exposed: one corresponding to the UV radiation originating from stars within the DLA, itself, and the other corresponding to the UV background radiation. We find that, for larger clouds, the propagating I-fronts created by local stellar sources can trigger cooling instability and collapse of significant part, up to 85%, of the cloud, creating conditions for star formation in a timescale of a few Myr. The passage of the I-front also triggers collapse of smaller clumps (with radii below ~4 pc), but in these cases the resulting cold and dense gas does not reach conditions conducive to star formation. Assuming that 85% of the gas initially in the clump is converted into stars, we obtain a star formation rate of $\sim 0.25 M_\odot {yr}^{-1} {kpc}^{-2}$. This is somewhat higher than the value derived from recent observations. On the other hand, the background UV radiation which has harder spectrum fails to trigger cooling and collapse. Instead, the hard photons which have long mean-free-path heat the dense clumps, which as a result expand and essentially dissolve in the ambient medium. Therefore, the star formation activity in DLAs is strongly regulated by the radiative feedback, both from the external UV background and internal stellar sources and we predict quiescent evolution of DLAs (not starburst-like evolution).

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.

We investigate the possibility to search for metal-poor high column density (> 10^{23} cm^{-2}) clouds at high redshift (z) by using gamma ray burst (GRB) afterglows. Such clouds could be related to primeval galaxies which may cause a burst of star formation. We show that a large part of hydrogen is in molecular form in such a high column density environment. Therefore, hydrogen molecules (H2) rather than hydrogen atoms should be searched for. Then we show that infrared H2 lines are detectable for metal-poor (< 0.01 solar metallicity) high column density (log NH [cm^{-2}] > 23.5) systems at high-z without suffering dust extinction. Optical properties of dust in infrared could also be constrained by observations of high column density systems. Some possible scenarios of producing high column density systems are finally discussed in the context of galaxy evolution.

In the early epoch of galaxy evolution, dust is only supplied by supernovae (SNe). With the aid of a new physical model of dust production by SNe developed by Nozawa et al. (2003) (N03), we constructed a model of dust emission from forming galaxies on the basis of the theoretical framework of Takeuchi et al. (2003) (T03). N03 showed that the produced dust species depends strongly on the mixing within SNe. We treated both unmixed and mixed cases and calculated the infrared (IR) spectral energy distribution (SED) of forming galaxies for both cases. Our model SED is less luminous than the SED of T03 model by a factor of 2-3. The difference is due to our improved treatment of UV photon absorption cross section, as well as different grain size and species newly adopted in this work. The SED for the unmixed case is found to have an enhanced near to mid-IR (N-MIR) continuum radiation in its early phase of the evolution (age < 10^{7.25} yr) compared with that for the mixed case. The strong N--MIR continuum is due to the emission from Si grains, which only exist in the species of the unmixed dust production. We also calculated the IR extinction curves for forming galaxies. Then we calculated the SED of a local starbursting dwarf galaxy SBS 0335-052. Our present model SED naturally reproduced the strong N--MIR continuum and the lack of cold FIR emission of SBS 0335-052. We found that only the SED of unmixed case can reproduce the NIR continuum of this galaxy. We then made a prediction for the SED of another typical star-forming dwarf, I Zw 18. We also presented the evolution of the SED of LBGs. Finally, we discussed the possibility of observing forming galaxies at z > 5.

We investigate the extinction curves of young galaxies in which dust is supplied from Type II supernovae (SNe II) and/or pair instability supernovae (PISNe). We adopt Nozawa et al. (2003) for compositions and size distribution of grains formed in SNe II and PISNe. We find that the extinction curve is quite sensitive to internal metal mixing in supernovae (SNe). The extinction curves predicted from the mixed SNe are dominated by SiO2 and is characterised by steep rise from infrared to ultraviolet (UV). The dust from unmixed SNe shows shallower extinction curve, because of the contribution from large-sized (~ 0.1 um) Si grains. However, the progenitor mass is important in unmixed SNe II: If the progenitor mass is smaller than ~ 20 Msun, the extinction curve is flat in UV; otherwise, the extinction curve rises toward the short wavelength. The extinction curve observed in a high-redshift quasar (z=6.2) favours the dust production by unmixed SNe II. We also provide some useful observational quantities, so that our model might be compared with future high-z extinction curves.

In order to interpret H2 (molecular hydrogen) quasar absorption line observations of damped Ly-alpha systems (DLAs) and sub-DLAs, we model their H2 abundance as a function of dust-to-gas ratio, including H2 self-shielding and dust extinction against dissociating photons. Then, we constrain the physical state of gas by using H2 data. Using H2 excitation data for DLA with H2 detections, we derive a gas density 1.5 < log n [cm^-3] < 2.5, temperature 1.5 < log T [K] < 3, and internal UV radiation field (in units of the Galactic value) 0.5 < log \chi < 1.5. We then find that the observed relation between molecular fraction and dust-to-gas ratio of the sample is naturally explained by the above conditions. However, it is still possible that H2 deficient DLAs and sub-DLAs with H2 fractions less than ~ 10^-6 are in a more diffuse and warmer state. The efficient photodissociation by the internal UV radiation field explains the extremely small H2 fraction (< 10^-6) observed for \kappa < 1/30 (\kappa is the dust-to-gas ratio in units of the Galactic value); H2 self-shielding causes a rapid increase and the large variations of H2 abundance for \kappa > 1/30. We finally propose an independent method to estimate the star formation rates of DLAs from H2 abundances; such rates are then critically compared with those derived from other proposed methods. The implications for the contribution of DLAs to the cosmic star formation history are briefly discussed.

We investigate the role of dust in star formation activity of extremely metal-poor blue compact dwarf galaxies (BCDs). Observations suggest that star formation in BCDs occurs in two different regimes: "active" and "passive". The "active" BCDs host super star clusters (SSCs), and are characterised by compact size, rich H2 content, large dust optical depth, and high dust temperature; the "passive" BCDs are more diffuse with cooler dust, and lack SSCs and large amounts of H2. By treating physical processes concerning formation of stars and dust, we are able to simultaneously reproduce all the above properties of both modes of star formation (active and passive). We find that the difference between the two regimes can be understood through the variation of the "compactness" of the star-forming region: an "active" mode emerges if the region is compact (with radius $\la 50$ pc) and dense (with gas number density $\ga 500$ cm$^{-3}$). The dust, supplied from Type II supernovae in a compact star-forming region, effectively reprocesses the heating photons into the infrared and induces a rapid H2 formation over a period of several Myr. This explains the high infrared luminosity, high dust temperature, and large H2 content of active BCDs. Moreover, the gas in "active" galaxies cools ($\la 300$ K) on a few dynamical timescales, producing a "run-away" star formation episode because of the favourable (cool) conditions. The mild extinction and relatively low molecular content of passive BCDs can also be explained by the same model if we assume a diffuse region (with radius $\ga 100$ pc and gas number density $\la 100$ cm$^{-3}$). We finally discuss primordial star formation in high-redshift galaxies in the context of the "active" and "passive" star formation scenario.

Understanding how galaxies formed their first stars is a vital cosmological question, but the study of high-redshift objects, caught in the act of forming their first stars, is difficult. Here we argue that two extremely low-metallicity Blue Compact Dwarf galaxies (BCDs), IZw18 and SBS0335-052, could be local templates for primordial star formation, since both lack evolved ($> $1 Gyr) stellar populations; but they form stars differently.

Infrared (IR) luminosity of galaxies originating from dust emission can be used as an indicator of the star formation rate (SFR). Inoue et al. (2000, IHK) have derived a formula for the conversion from IR luminosity to SFR by using the following three quantities: the fraction of Lyman continuum luminosity absorbed by gas (f), the fraction of UV luminosity absorbed by dust (epsilon), and the fraction of dust heating from old (>10^8 yr) stellar populations (eta). We develop a method to estimate those three quantities based on the idea that the various way of SFR estimates should return the same SFR. After applying our method to samples of galaxies, the following results are obtained. First, our method is applied to star-forming galaxies, finding that f~0.6, epsilon~0.5, and eta~0.4 as representative values. Next, we apply the method to a starburst sample, which shows larger extinction than the star-forming galaxy sample. With the aid of f, epsilon, and eta, we estimate reliable SFRs. Moreover, the H-alpha luminosity, if the H-alpha extinction is corrected by using the Balmer decrement, is suitable for a statistical analysis of SFR, because the same correction factor for the Lyman continuum extinction is applicable to both normal and starburst galaxies over all the range of SFR. The metallicity dependence of f and epsilon is also tested: Only the latter proves to have a correlation with metallicity. As an extension of our result, we show that all UV, H-alpha, and IR comoving luminosity densities at z=0 give a consistent SFR (~ 3x10^{-2}h M_sun/Mpc^3). Useful formulae for SFR estimate are listed.

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.

To interpret H_2 quasar absorption line observations in Damped Ly-alpha clouds (DLAs), we model the H_2 spatial distribution within a DLA. Based on numerical simulations of disk structures with parameters similar to those derived for such absorbers, we calculate the H_2 distribution as a function of ultraviolet background (UVB) intensity and dust-to-gas ratio. For typical values of these two quantities we find that the area in which the H_2 fraction exceeds 10^{-6} (typical observational detection limit) only covers $\la 10$% of the disk surface, i.e. H_2 has a very inhomogeneous, clumpy distribution even at these low abundance levels. This explains the relative paucity of H_2 detections in DLAs. We also show the dependence of the covering fraction of H_2 on dust-to-gas ratio and UVB intensity and we comment on the physics governing the H_2 chemical network at high redshift.

Before the end of 2002 will be launched the GALEX satellite (a NASA/SMEX project) which will observe all the sky in Ultraviolet (UV) through filters at 1500 and 2300 A down to m(AB) 21. In 2004 will be launched the ASTRO-F satellite which will perform an all sky survey at Far-Infrared (FIR) wavelengths. The cross-correlation of both surveys will lead to very large samples of galaxies for which FIR and UV fluxes will be available. Using the FIR to UV flux ratio as a quantitative tracer of the dust extinction we will be able to measure the extinction in the nearby universe (z<0.2) and to perform a statistically significant analysis of the extinction as a function of galactic properties. Of particular interest is the construction of pure FIR and UV selected samples for which the extinction will be measured as templates for the observation of high redshift galaxies.

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 5<z<20. In particular, we focus on the interplay between dust formation in Type II supernova ejecta and H2 production on dust grains. Such effect causes an enhancement of star formation rate by an order of magnitude on a timescale (~3--5 galactic dynamical times) shorter than the Hubble timescale. We also find that about half of the radiative energy from stars is reprocessed by dust grains and is finally radiated in the far infrared (FIR). Typical star formation rates and luminosities (FIR, UV and metal-line luminosities) are calculated for a large set of (M_vir, z_vir). Using these results and the Press-Schechter formalism, we calculate galaxy number counts and integrated light from high-redshift (z>5) 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.

This paper investigates the origin of the observed large variety in dust-to-gas ratio among blue compact dwarf galaxies (BCDs). By applying our chemical evolution model, we find that the dust destruction can largely suppress the dust-to-gas ratio when the metallicity of a BCD reaches $12+\log{\rm (O/H)}\sim 8$, i.e., a typical metallicity level of BCDs. We also show that dust-to-gas ratio is largely varied owing to the change of dust destruction efficiency that has two effects: (i) a significant contribution of Type Ia supernovae to total supernova rate; (ii) variation of gas mass contained in a star-forming region. While mass loss from BCDs was previously thought to be the major cause for the variance of dust-to-gas ratio, we suggest that the other two effects are also important. We finally discuss the intermittent star formation history, which naturally explains the large dispersion of dust-to-gas ratio among BCDs.

During the early stages of galaxy evolution, the metallicity is generally low and nearby metal-poor star-forming galaxies may provide templates for primordial star formation. In particular, the dust content of such objects is of great importance since early molecular formation can take place on grains. To gain insight into primeval galaxies at high redshift, we examine the dust content of the nearby extremely low-metallicity galaxy SBS 0335-052 which hosts a very young starburst (< 10 Myr). In young galaxies, the dust formation rate in Type II supernovae governs the amount of dust, and by incorporating recent results on dust production in Type II supernovae we model the evolution of dust content. If the star-forming region is compact (< 100 pc), as suggested by observations of SBS 0335-052, our models consistently explain the quantity of dust, far-infrared luminosity, and dust temperature in this low-metallicity object. We also discuss the H_2 abundance. The compactness of the region is important to H_2 formation, because the optical depth of dust for UV photons becomes large and H_2 dissociation is suppressed. We finally focus on implications for damped Ly alpha systems.

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.

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.

The virial mass ($M_{\rm vir}$)-metallicity relation among the Local Group dwarf spheroidal galaxies (dSphs) is examined. Hirashita, Takeuchi, & Tamura showed that the dSphs can be divided into two distinct classes with respect to the relation between their virial masses and luminosities: low-mass ($M_{\rm vir} \la 10^8 M_\odot$) and high-mass ($M_{\rm vir} \ga 10^8 M_\odot$) groups. We see that both the mass-metallicity and the mass-luminosity relations of the high-mass dSphs are understood as a low-mass extension of giant ellipticals. On the contrary, we find that the classical galactic-wind model is problematic to apply to the low-mass dSphs, whose low binding energy is comparable to that released by several supernova explosions. A strongly regulated star formation in their formation phase is required to reproduce their observed metallicity. Such regulation is naturally expected in a gas cloud with the primordial elemental abundance according to Nishi & Tashiro. A significant scatter in the mass-metallicity relation for the low-mass dSphs is also successfully explained along with the scenario of Hirashita and coworkers. We not only propose a new picture for a chemical enrichment of the dSphs, but also suggest that the mass-metallicity and the mass-luminosity relations be understood in a consistent context.

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.

We investigate what we can learn about galaxy formation and evolution from the data which will be obtained by the forthcoming large submillimeter/millimeter facilities, mainly by the Atacama Submillimeter Telescope Experiment (ASTE) and the Atacama Large Millimeter Array/Large Millimeter and Submillimeter Array (ALMA/LMSA). We first calculate the source counts from 350 um to 3 mm using the empirical infrared galaxy number count model of Takeuchi et al. (2001). Based on the number counts, we evaluate the source confusion and determine the confusion limit at various wavebands as a function of the characteristic beam size. At submillimeter wavelengths, source confusion with the 10-15-m class facilities becomes severe at 0.1 to 1 mJy level, and astrometry and flux measurements will be difficult. However, we show that very a large-area survey of submillimeter sources brighter than 10-50 mJy can provide a unique constraint on infrared galaxy evolution at z = 1-2. We also find that the 5\sigma-confusion limit of the LMSA is fainter than 1 \mu Jy. The source counts at such faint flux levels give important information of the epoch of galaxy formation. We then show that multiband photometry from the infrared (by ASTRO-F) to the millimeter can be utilized as a redshift estimator. In addition, we compare the observed 1.4, 5, and 8-GHz source counts with our model counts to examine the contribution of star forming galaxies to the faint radio galaxies.

We model the star formation history (SFH) and the chemical evolution of the Galactic disk by combining an infall model and a limit-cycle model of the interstellar medium (ISM). Recent observations have shown that the SFH of the Galactic disk violently variates or oscillates. We model the oscillatory SFH based on the limit-cycle behavior of the fractional masses of three components of the ISM. The observed period of the oscillation ($\sim 1$ Gyr) is reproduced within the natural parameter range. This means that we can interpret the oscillatory SFH as the limit-cycle behavior of the ISM. We then test the chemical evolution of stars and gas in the framework of the limit-cycle model, since the oscillatory behavior of the SFH may cause an oscillatory evolution of the metallicity. We find however that the oscillatory behavior of metallicity is not prominent because the metallicity reflects the past integrated SFH. This indicates that the metallicity cannot be used to distinguish an oscillatory SFH from one without oscillations.

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.

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.

It has been reported recently that there are some early-type spiral (Sa--Sab) galaxies having evident star-forming regions which concentrate in their own central 1-kpc. In such central region, is the mechanism of the star formation distinct from that in disks of spiral galaxies? To reveal this, we estimate the star formation efficiency (SFE) in this central 1-kpc star-forming region of some early-type spiral galaxies, taking account of the condition for this 1-kpc region to be self-gravitating. Using two indicators of present star formation rate (H$\alpha$ and infrared luminosity), we estimate the SFE to be a few percents. This is equivalent to the observational SFE in the disks of late-type spiral (Sb--) galaxies. This coincidence may support the universality of the mean SFE of spiral galaxies reported in the recent studies. That is, we find no evidence of distinct mechanism of the star formation in the central 1-kpc region of early-type galaxies. Also, we examine the structure of the central star-forming region, and discuss a method for estimating the mass of star-forming regions.

We propose a limit-cycle scenario of star formation history for any morphological type of spiral galaxies. It is known observationally that the early-type spiral sample has a wider range of the present star formation rate (SFR) than the late-type sample. This tendency is understood in the framework of the limit-cycle model of the interstellar medium (ISM), in which the SFR cyclically changes in accordance with the temporal variation of the mass fraction of the three ISM components. When the limit-cycle model of the ISM is applied, the amplitude of variation of the SFR is expected to change with the supernova (SN) rate. Observational evidence indicates that the early-type spiral galaxies show smaller rates of present SN than late-type ones. Combining this evidence with the limit-cycle model of the ISM, we predict that the early-type spiral galaxies show larger amplitudes in their SFR variation than the late-types. Indeed, this prediction is consistent with the observed wider range of the SFR in the early-type sample than in the late-type sample. Thus, in the framework of the limit-cycle model of the ISM, we are able to interpret the difference in the amplitude of SFR variation among the morphological classes of spiral galaxies.

We analyze the star formation history (SFH) of the Galactic disk by using an infall model. Based on the observed SFH of the Galactic disk, we first determine the timescales of the gas infall into the Galactic disk ($t_{\rm in}$) and that of the gas consumption to form stars ($t_{\rm sf}$). Since each of the two timescales does not prove to be determined independently from the SFH, we first fix $t_{\rm sf}$. Then, $t_{\rm in}$ is determined so that we minimize $\chi^2$. Consequently, we choose three parameter sets: ($t_{\rm sf}$ [Gyr], $t_{\rm in}$ [Gyr]) = (6.0, 23), (11, 12), and (15, 9.0), where we set the Galactic age as 15 Gyr. All of the three cases predict almost identical star formation history. Next, we test the intermittency (or variability) of the star formation rate (SFR) along with the smooth SFH suggested from the infall model. The large value of the $\chi^2$ statistic supports the violent time-variation of the SFH. If we interpret the observed SFH with smooth and variable components, the amplitude of the variable component is comparable to the smooth component. Thus, intermittent SFH of the Galactic disk is strongly suggested. We also examined the metallicity distribution of G-dwarfs. We found that the true parameter set lies between ($t_{\rm sf}$ [Gyr], $t_{\rm in}$ [Gyr]) = (6, 23) and (11, 12), though we should need a more sophisticated model including the process of metal enrichment within the Galactic halo.

We construct a new algorithm for estimating the star formation rate (SFR) of galaxies from their infrared (IR) luminosity by developing the theory of the IR emission from a dusty H II region. The derived formula is ${\rm SFR}/(\MO {\rm yr^{-1}})=\{3.3\times 10^{-10}(1-\eta)/ (0.4-0.2f+0.6 \epsilon)\}(L{\rm ^{obs}_{IR}}/\LO)$, where $f$ is the fraction of ionizing photons absorbed by hydrogen, $\epsilon$ is the efficiency of dust absorption for nonionizing photons from OB stars, and $\eta$ is the cirrus fraction of observed IR luminosity. The previous conversion formulae of SFR from the IR luminosity is applicable to only the case where the observed IR luminosity is nearly equal to the bolometric luminosity (starburst galaxies etc.), except for some empirical formulae. On the other hand, our theoretical SFR is applicable to galaxies even with a moderate star formation activity. That is, our simple and convenient formula is significantly useful for estimating the SFR of various morphologies and types of galaxies --- from early elliptical to late spiral and irregular galaxies, or from active starburst to quiescent galaxies --- as far as they have neither an abnormal dust-to-gas ratio nor an evident active galactic nucleus.

We constructed a model of infrared and sub-mm (hereafter IR) galaxy number count and estimated history of the IR luminosity density. We treat the evolutionary change of galaxy luminosities as a stepwise nonparametric form, in order to explore the most suitable evolutionary history which reproduces the present observational results. We found the evolutionary patterns which satisfy both constraints required from Cosmic Infrared Background (CIRB) and IR galaxy number counts. One order of magnitude increase of luminosity at redshift $z=0.75 - 1.0$ was found in IR $60 \mu$m luminosity density evolution. We also found that a large number of galaxies ($\sim 10^7$ in the whole sky) will be detected in all-sky survey at far-infrared by Infrared Imaging Surveyor (IRIS); Japanese infrared satellite project Astro-F.

On the basis of the propagating star-formation model, we investigated the star-formation activities of dwarf irregular galaxies (dIrrs) by considering two processes: the heating from stellar feedback and the cooling of the heated gas. After examining the timescales of the two processes, we propose that continuous global star-formation activity is difficult in dIrrs, since their small sizes make the feedback efficient and their small metallicities prevent the cooling from becoming effective. Thus, the intermittent nature of the star-formation activities of dIrrs, which is due to the small metallicity as well as the small size, is predicted. We emphasize that the size of a galaxy is an important factor concerning star-formation activity. The intermittence of the star-formation activity is also supported by the observed scatter of the UBV colors and H$\alpha$ equivalent widths of the dIrr sample. However, we note that efficient interstellar mixing may make the cooling time much shorter.

In this paper, we examine whether tidal forces exerted by the Galaxy or M31 have an influence on the Local Group dwarf spheroidal galaxies (dSphs) which are their companions. We focus on the surface brightness profiles of the dSphs, especially their core radii because it is suggested based on the numerical simulations that tidal disturbance can make core radii extended. We examine the correlation for the dSphs between the distances from their parent galaxy (the Galaxy or M31) and the compactnesses of their surface brightness profiles by using a parameter ``C'' defined newly in this paper. Consequently, we find no significant correlation. We make some remarks on the origin of this result by considering three possible scenarios; tidal picture, dark matter picture, and heterogeneity of the group of dSphs, each of which has been often discussed to understand fundamental properties and formation processes of dSphs.

We discuss the time variation of dust-to-gas mass ratio in spiral galaxies using the multi-phase model of interstellar medium. The typical timescale of the phase change of an interstellar gas is $\sim 10^7$--$10^8$ yr in spiral galaxies. Since the phase transition changes the filling factor of the cold gas where the dust growth occurs, the dust growth rate varies on that timescale. In order to examine the response of the dust-to-gas ratio to the phase transition, we construct a model of the time evolution of the dust-to-gas ratio. We adopt the three phase model for the interstellar gas and the Ikeuchi-Tomita model for the mass exchange among the phases. According to the model, three types of solutions are possible: [1] all the gas is transformed to a hot gas; [2] a stable stationary state of three phases is realized; [3] the filling factor of each phase cyclically changes. For each of the three types of solutions, the dust-to-gas ratio behaves as follows: [1] almost all the dust is destroyed (the dust-to-gas ratio becomes $\sim 0$); [2] the dust-to-gas ratio converges to a stationary state; [3] the dust-to-gas ratio varies cyclically in response to the phase transition. In the case of [3], the amplitude of the variation of the dust-to-gas ratio is large (nearly an order of magnitude) if the dust growth timescale is shorter than the phase transition timescale. This condition is easily satisfied in spiral galaxies. However, it is difficult for dwarf galaxies to realize the condition because their small metallicity makes the dust growth timescale long.

We examine the correlation between physical quantities to explore the existence of dark matter (DM) in the Local Group dwarf spheroidal galaxies (dSphs). In order to clarify whether DM exists in the dSphs we compare two extreme models of their internal kinematics: [1] a tidal model (the mass of a dSph is estimated with the luminous mass, since the observed large velocity dispersions of the dSphs are attributed to the tidal force by the Galaxy), and [2] a DM model (the mass of a dSph is estimated with the virial mass, since the velocity dispersion is considered to reflect the DM potential). In both models, we find that the relation between the surface brightness of the dSphs and the tidal force by the Galaxy is consistently interpreted. This makes a critical remark about the previous studies that concluded that the tidal force is effective for the dSphs based only on model [1]. We also check the correlation between the Galactocentric distance and the tidal force. Consequently, both models are also supported in this test. Thus, we are unable to judge which of the two is more promising for the dSphs in correlation investigations. The physical process of tidal disruption is also considered. Since the timescale for the tidally perturbed stars to escape is much shorter than the orbital periods of the dSphs, it is difficult for the tidally perturbed dSphs to exist as an assembly. Thus, we suggest that the gravitational field of DM is necessary to bind the dSphs with a large velocity dispersion.

We examine the dust-to-gas ratio as a function of metallicity for dwarf galaxies [dwarf irregular galaxies (dIrrs) and blue compact dwarf galaxies (BCDGs)]. Using a one-zone model and adopting the instantaneous recycling approximation, we prepare a set of basic equations which describes processes of dust formation and destruction in a galaxy. Four terms are included for the processes: dust formation from heavy elements ejected by stellar mass loss, dust destruction in supernova remnants, dust destruction in star-forming regions, and accretion of heavy elements onto preexisting dust grains. Solving the equations, we compare the result with observational data of nearby dIrrs and BCDGs. The solution is consistent with the data within the reasonable ranges of model parameters constrained by the previous examinations. This means that the model is successful in understanding the dust amount of nearby galaxies. We also show that the accretion rate of heavy element onto preexisting dust grains is less effective than the condensation of heavy elements in dwarf galaxies.

We discuss the time evolution of dust-to-gas mass ratio in the context of multi-phase model of interstellar medium. Phase transition of interstellar gas is considered to occur on a timescale of $10^{7--8}$ yr, according to a nonlinear open system model of interstellar medium. Since the phase transition changes the dust formation and destruction rates, the dust-to-gas ratio also fluctuates on the same timescale. This explains the scatter of the dust-to-gas ratios of spiral galaxies quantitatively, though we should note the large observational uncertainty.

We consider evolution of interstellar clouds in Local Group dwarf spheroidal galaxies (dSphs) in the context of their observed star formation histories. The Local Group dSphs generally experienced initial bursts of star formations in their formation epochs ($\sim 15$ Gyr ago), when hot gas originating from the supernovae can make the cold interstellar clouds evaporate. We find that the maximum size of evaporating cloud is 10 pc. Thus, clouds larger than 10 pc can survive during the initial star formation. These surviving clouds can contribute to the second star formation to produce ``intermediate-age ($\sim$ 3--10 Gyr ago) stellar populations.'' Assuming that collisions between clouds induce star formation and that the timescale of the second star formation is a few Gyr, we estimate the total mass of the clouds. The total mass is about $10^{4}M_\odot$, which is 1--3 orders of magnitude smaller than the typical stellar mass of a present dSph. This implies that the initial star formation is dominant over the second star formation, which is broadly consistent with the observed star formation histories. However, the variety of the dSphs in their star formation histories suggests that the effects of environments on the dSphs may be important.

Infrared Imaging Surveyor (IRIS) is a satellite which will be launched at the beginning of 2003. One of the main purposes of the IRIS mission is an all-sky survey in the far-infrared region with a flux limit much deeper than that of IRAS. Detection of a large number of galaxies ($\sim several \times 10^6$ in the whole sky) is expected in this survey. We investigated the expected optical and near-infrared (NIR) number counts of galaxies detected by the far-infrared scanner (FIS) of IRIS (hereafter, IRIS galaxies) and possibility of optical and NIR follow-up of them. The spectral energy distribution and the luminosity function of the IRIS galaxies are modeled based on the properties of galaxies observed by IRAS. The IRIS galaxies are divided into two populations according to their infrared luminosities ($L_{IR}$); normal spirals ($L_{IR}<10^{10}\LO$) and starbursts ($L_{IR}>10^{10}\LO$). The expected number counts of IRIS galaxies for both of the populations are calculated in $B$ and $H$ bands. We show that about 60 normal galaxies and about 80 starburst galaxies are detected per square degree in both of the two bands, when galaxy evolution is not taken into account. All of the normal population of IRIS galaxies are located at the redshift $z\ltsim 0.1$. As for the starburst population, we also calculated number of galaxies with a simple model of evolution. The total number of starburst population predicted by the evolution model is larger by 20% than that expected from the non-evolution model. In the evolution model, the numbers of low-z ($z<1$), intermediate-z ($1<z<3$), and high-z ($z>3$) galaxies are 100, 20, and 0.2 per square degree, respectively.

Infrared Imaging Surveyor (IRIS) is a satellite which will be launched in the beginning of 2003. One of the main purposes of the IRIS mission is an all-sky survey at far-infrared (FIR) with a flux limit much deeper than that of IRAS. In order to examine the performance of the survey, we estimated the FIR galaxy counts in four (50, 70, 120, and 150 $\mu$m) bands based on some models. We adopted a multicomponent model which consists of cirrus and starburst components for galaxy spectra, and the nearby FIR luminosity function derived from that of IRAS galaxies. We derived the number counts, redshift distributions, and infrared diffuse background radiation spectra for i) no evolution, ii) pure luminosity evolution, iii) pure density evolution with $q_0 = 0.1$ and 0.5. We found that a large number of galaxies ($\sim {\rm a few} \times 10^6$ in the whole sky) will be detected in this survey. With the aid of a vast number of detection, we will detect the effect of galaxy evolution, and evaluate the amplitude of evolution at least in the nearby universe in the IRIS survey, though it will be still difficult to constrain which type of evolution takes place from the number count alone. We also studied the estimation of redshifts of detected galaxies by their infrared colors alone. Although significant contamination takes place among nearby faint galaxies and high-$z$ ones, we found that rough estimation of galaxy redshift can be practicable by jointly using present and future optical surveys.

Using a one-zone model for the evolution of dust in spiral galaxies and applying the instantaneous recycling approximation to the model equations, we investigate dust-to-gas ratio of spiral galaxies. Four processes are considered; dust formation from heavy elements ejected by stellar mass loss, dust destruction in supernova remnants, dust destruction in star-forming regions, and accretion of heavy elements onto preexisting dust grains. The equations describe a simplified relation between dust-to-gas ratio and metallicity. The relation is independent of star formation rate. By comparing the theoretical relation with the observational data of nearby spiral galaxies, we show (i) that the accretion process onto the preexisting dust particles must be taken into account for the spiral galaxies; (ii) that the efficiency of dust production from heavy elements ejected by stars can be constrained by the spiral galaxies with low metallicity; (iii) that the Salpeter and Scalo initial mass functions are both consistent with the data of the spiral galaxies.

We analyzed 10 dwarf spheroidal galaxies (dSphs) in the Local Group, and found two distinct sequences on the $M_{\rm vir} / L$ -- $M_{\rm vir}$ plane: $M_{\rm vir} / L \propto M_{\rm vir}^{1.6}$ for $M_{\rm vir} < 10^8 M_\odot$ whereas $M_{\rm vir} / L \simeq {\rm const.}$ for $M_{\rm vir} > 10^8 M_\odot$ ($M_{\rm vir}$ and $L$ are the virial mass and the total luminosity of a dSph, respectively). We interpret the disconuity as the threshold for the gas in dSphs to be blown away by successive supernovae. We succeeded to give a quantitative explanation of the discontinuity mass of $10^8M_\odot$ as blow-away condition. We further derived the above relation for the low-mass dSphs, assuming that the initial star formation rate of the dSphs is proportional to the inverse of the cooling time. The relation of high-mass dSphs is also explained along with the same consideration, with the condition that the gas cannot been blown away effectively.

We discuss a physical interpretation of the relation between $M_{\rm vir}/L$ and $M_{\rm vir}$ of dwarf spheroidal galaxies (dSphs), where $M_{\rm vir}$ and $L$ are the virial mass and the total luminosity of a dSph, respectively. We used 11 dSphs in the Local Group as the sample. We find two distinct sequences on the $M_{\rm vir}/L$--$M_{\rm vir}$ plane: $M_{\rm vir}/L\propto M_{\rm vir}^{2.0}$ for dSphs with $M_{\rm vir}<10^8M_\odot$, whereas $M_{\rm vir}/L\simeq$ constant for $M_{\rm vir}>10^8M_\odot$. A ``discontinuity'' is seen at $M_{\rm vir}\simeq 10^8M_\odot$. We interpret the ``discontinuity'' as the threshold for the gas in dSphs to be blown away by successive supernovae. If a dSph has virial mass (most of which is dark mass) less than $10^8M_\odot$, the gas is blown away, while in a dSph of larger mass, the deep potential well prevents the blow-away mechanism from working effectively. Thus, large mass ratio of dark matter (DM) to baryonic matter (i.e., large $M_{\rm vir}/L$) is realized in a low-mass ($M_{\rm vir}<10^8M_\odot$) dSph through the gas depletion, whereas $M_{\rm vir}/L$ becomes lower in high-mass ($M_{\rm vir}>10^8M_\odot$) dSphs. We further make an attempt to explain the above relation for the low-mass dSphs, $M_{\rm vir}/L\propto M_{\rm vir}^{2.0}$, based on estimate of cooling time, using the scaling laws of virial temperature, virial mass and radius of a dSph and assuming that the heating by OB-star radiation terminates the star formation activity. We succeed in deriving the above relation for the mass-to-light ratio and luminosity.

We examine the possibility that ram pressure exerted by the galactic wind from the Galaxy could have stripped gas from the Local Group dwarf galaxies, thereby affecting their star formation histories. Whether gas stripping occurs or not depends on the relative magnitudes of two counteracting forces acting on gas in a dwarf galaxy: ram pressure force by the wind and the gravitational binding force by the dwarf galaxy itself. We suggest that the galactic wind could have stripped gas in a dwarf galaxy located within the distance of $R_{c}\simeq 120(r_{s}/1 {kpc})^{3/2} ({\cal E}_{b}/10^{50} {erg})^{-1/2}$ kpc (where $r_{s}$ is the surface radius and ${\cal E}_{b}$ is the total binding energy of the dwarf galaxy, respectively) from the Galaxy within a timescale of Gyr, thereby preventing star formation there. Our result based on this Galactic wind model explains the recent observation that dwarfs located close to the Galaxy experienced star formation only in the early phase of their lifetimes, whereas distant dwarfs are still undergoing star formation. The present star formation in the Large Magellanic Cloud can also be explained through our Galactic wind model.