We study the formation of dust in the expanding gas ejected as a result of a common envelope binary interaction. In our novel approach, we apply the dust formation model of Nozawa et al. to the outputs of the 3D hydrodynamic SPH simulation performed by Iaconi et al., that involves a giant of 0.88~\ms \ and 83~\rs, with a companion of 0.6~\ms \ placed on the surface of the giant in circular orbit. After simulating the dynamic in-spiral phase we follow the expansion of the ejecta for $\simeq 18\,000$~days. During this period the gas is able to cool down enough to reach dust formation temperatures. Our results show that dust forms efficiently in the window between $\simeq 300$~days (the end of the dynamic in-spiral) and $\simeq 5000$~days. The dust forms in two separate populations; an outer one in the material ejected during the first few orbits of the companion inside the primary's envelope and an inner one in the rest of the ejected material. The inner dust population dominates the grain size distribution at the end of the simulation. We are able to fit the grain size distribution at the end of the simulation with a double power law. The slope of the power law for smaller grains is flatter than that for larger grains, creating a knee-shaped distribution. The power law indexes are however different from the classical values determined for the interstellar medium. We also estimate that the contribution to cosmic dust by common envelope events is not negligible and comparable to that of novae and supernovae.

We investigate the common envelope binary interaction, that leads to the formation of compact binaries, such as the progenitor of Type Ia supernovae or of mergers that emit detectable gravitational waves. In this work we diverge from the classic numerical approach that models the dynamic in-spiral. We focus instead on the asymptotic behaviour of the common envelope expansion after the dynamic in-spiral terminates. We use the SPH code {\sc phantom} to simulate one of the setups from Passy et al., with a 0.88~\ms, 83~\rs \ RGB primary and a 0.6~\ms \ companion, then we follow the ejecta expansion for $\simeq 50$~yr. Additionally, we utilise a tabulated equation of state including the envelope recombination energy in the simulation (Reichardt et al.), achieving a full unbinding. We show that, as time passes, the envelope's radial velocities dominate over the tangential ones, hence allowing us to apply an homologous expansion kinematic model to the ejecta. The external layers of the envelope become homologous as soon as they are ejected, but it takes $\simeq 5000$~days ($\simeq 14$~yr) for the bulk of the unbound gas to achieve an homologous expanding regime. We observe that the complex distribution generated by the dynamic in-spiral evolves into a more ordered, ring-like shaped one in the asymptotic regime. We show that the thermodynamics of the expanding envelope are in very good agreement with those expected for an adiabatically expanding sphere under the homologous condition and give a prediction for the location and temperature of the photosphere assuming dust to be the main source of opacity.

Chondrites are one of the most primitive objects in the solar system, and keep the record of the degree of thermal metamorphism experienced in their parent bodies. This thermal history can be classified by the petrologic type. We investigate the thermal evolution of planetesimals to account for the current abundances (known as the fall statistics) of petrologic types 3 - 6 ordinary chondrites. We carry out a number of numerical calculations in which formation times and sizes of planetesimals are taken as parameters. We find that planetesimals that form within 2.0 Myr after the formation of Ca-Al-rich inclusions (CAIs) can contain all petrologic types of ordinary chondrites. Our results also indicate that plausible scenarios of planetesimal formation, which are consistent with the fall statistics, are that planetesimals with radii larger than 60 km start to form around 2.0 Myr after CAIs and/or that ones with radii less than 50 km should be formed within 1.5 Myr after CAIs. Thus, thermal modelling of planetesimals is important for revealing the occurrence and amount of metamorphosed chondrites, and for providing invaluable insights into planetesimal formation.

Extremely metal-poor (EMP) stars are the living fossils with records of chemical enrichment history at the early epoch of galaxy formation. By the recent large observation campaigns, statistical samples of EMP stars have been obtained. This motivates us to reconsider their classification and formation conditions. From the observed lower-limits of carbon and iron abundances of Acr(C) ~ 6 and [Fe/H]cr ~ -5 for C-enhanced EMP (CE-EMP) and C-normal EMP (CN-EMP) stars, we confirm that gas cooling by dust thermal emission is indispensable for the fragmentation of their parent clouds to form such low-mass, i.e., long-lived stars, and that the dominant grain species are carbon and silicate, respectively. We constrain the grain radius r_i^{cool} of a species i and condensation efficiency f_{ij} of a key element j as r_C^{cool}/f_{C, C} = 10 um and r_Sil^{cool}/f_{Sil, Mg} = 0.1 um to reproduce Acr(C) and [Fe/H]cr, which give a universal condition 10^{[C/H]-2.30} + 10^[Fe/H] > 10^{-5.07} for the formation of every EMP star. Instead of the conventional boundary [C/Fe] = 0.7 between CE- and CN-EMP stars, this condition suggests a physically meaningful boundary [C/Fe]_{b} = 2.30 above and below which carbon and silicate grains are dominant coolants, respectively.

The abundant forms in which the major elements in the universe exist have been determined from numerous astronomical observations and meteoritic analyses. Iron (Fe) is an exception, in that only depletion of gaseous Fe has been detected in the interstellar medium, suggesting that Fe is condensed into a solid, possibly the astronomically invisible metal. To determine the primary form of Fe, we replicated the formation of Fe grains in gaseous ejecta of evolved stars by means of microgravity experiments. We found that the sticking probability for formation of Fe grains is extremely small; only several atoms will stick per hundred thousand collisions, so that homogeneous nucleation of metallic Fe grains is highly ineffective, even in the Fe-rich ejecta of Type Ia supernovae. This implies that most Fe is locked up as grains of Fe compounds or as impurities accreted onto other grains in the interstellar medium.

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.

Presolar grains are small particles found in meteorites through their isotopic compositions which are considerably different from those of materials in the Solar System. If some isotopes in presolar grains diffused out beyond their grain sizes when they were embedded in parent bodies of meteorites, their isotopic compositions could be washed out, and hence the grains cannot be identified as presolar grains any more. We explore this possibility for the first time by self-consistently simulating the thermal evolution of planetesimals and the diffusion length of $^{18}$O in presolar silicate grains. Our results show that presolar silicate grains smaller than $\sim$ 0.03 $\mu m$ cannot keep their original isotopic compositions even if the host planetesimals experienced maximum temperature as low as 600 $^{\circ}$C. Since this temperature corresponds to the one experienced by petrologic type 3 chondrites, the isotopic diffusion can constrain the size of presolar silicate grains discovered in such chondrites to be larger than $\sim$ 0.03 $\mu m$. We also find that the diffusion lengths of $^{18}$O reach $\sim$ 0.3-2 $\mu m$ in planetesimals that were heated up to 700-800 $^{\circ}$C. This indicates that, if the original size of presolar grains spans a range from $\sim$ 0.001 $\mu m$ to $\sim$ 0.3 $\mu m$ like that in the interstellar medium, the isotopic records of the presolar grains may be almost completely lost in such highly thermalized parent bodies. We propose that isotopic diffusion could be a key process to control the size distribution and abundance of presolar grains in some types of chondrites.

It is well known that the extinction properties along lines of sight to Type Ia supernovae are described by steep extinction curves with unusually low total-to-selective extinction ratios of Rv = 1.0-2.0. In order to reveal the properties of interstellar dust that causes such peculiar extinction laws, we perform the fitting calculations to the measured extinction curves by applying a two-component dust model composed of graphite and silicate. As for the size distribution of grains, we consider two function forms of the power-law and lognormal distributions. We find that the steep extinction curves derived from the one-parameter formula by Cardelli et al. (1989) with Rv = 2.0, 1.5, and 1.0 can be reasonably explained even by the simple power-law dust model that has a fixed power index of -3.5 with the maximum cut-off radii of a_{max} = 0.13 um, 0.094 um, and 0.057 um, respectively. These maximum cut-off radii are smaller than a_{max} ~ 0.24 um considered to be valid in the Milky Way, clearly demonstrating that the interstellar dust responsible for steep extinction curves is highly biased to smaller sizes. We show that the lognomal size distribution can also lead to good fits to the extinction curves with Rv = 1.0-3.1 by taking the appropriate combinations of the relevant parameters. We discuss that the extinction data at ultraviolet wavelengths are essential for constraining the composition and size distribution of interstellar dust.

Many astronomical objects are surrounded by dusty environments. In such dusty objects, multiple scattering processes of photons by circumstellar (CS) dust grains can effectively alter extinction properties. In this paper, we systematically investigate effects of multiple scattering on extinction laws for steady-emission sources surrounded by the dusty CS medium, using a radiation transfer simulation based on the Monte Carlo technique. In particular, we focus on whether and how the extinction properties are affected by properties of CS dust grains, adopting various dust grain models. We {\bf confirm} that behaviors of the (effective) extinction laws are highly dependent on the properties of CS grains. Especially, the total-to-selective extinction ratio $R_{V}$, which characterizes the extinction law, can be either increased or decreased, compared to the case without multiple scattering. We find that the criterion for this behavior is given by a ratio of albedos in the $B$ and $V$ bands. We also find that either small silicate grains or polycyclic aromatic hydrocarbons (PAHs) are necessary for realizing a low value of $R_{V}$ as often measured toward Type Ia supernovae, if the multiple scattering by CS dust is responsible for their non-standard extinction laws. Using the derived relations between the properties of dust grains and the resulting effective extinction laws, we propose that the extinction laws toward dusty objects could be used to constrain the properties of dust grains in CS environments.

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.

We present infrared multi-epoch observations of the dust forming nova V1280 Sco over $\sim$2000 days from the outburst. The temporal evolution of the infrared spectral energy distributions at 1272, 1616 and 1947 days can be explained by the emissions produced by amorphous carbon dust of mass (6.6--8.7)$\times$10$^{-8}$M$_{\odot}$ with a representative grain size of 0.01$~\mu$m and astronomical silicate dust of mass (3.4--4.3)$\times$10$^{-7}$M$_{\odot}$ with a representative grain size of 0.3--0.5$~\mu$m. Both of these dust species travel farther away from the white dwarf without an apparent mass evolution throughout those later epochs. The dust formation scenario around V1280 Sco suggested from our analyses is that the amorphous carbon dust is formed in the nova ejecta followed by the formation of silicate dust in the expanding nova ejecta or as a result of the interaction between the nova wind and the circumstellar medium.

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.

A few particles of presolar Al2O3 grains with sizes above 0.5 mum are believed to have been produced in the ejecta of core-collapse supernovae (SNe). In order to clarify the formation condition of such large Al2O3 grains, we investigate the condensation of Al2O3 grains for wide ranges of the gas density and cooling rate. We first show that the average radius and condensation efficiency of newly formed Al2O3 grains are successfully described by a non-dimensional quantity "Lambda_on" defined as the ratio of the timescale with which the supersaturation ratio increases to the collision timescale of reactant gas species at dust formation. Then, we find that the formation of submicron-sized Al2O3 grains requires at least ten times higher gas densities than those presented by one-dimensional SN models. This indicates that presolar Al2O3 grains identified as a SN origin might be formed in dense gas clumps, allowing us to propose that the measured sizes of presolar grains can be a powerful tool to constrain the physical conditions in which they formed. We also briefly discuss the survival of newly formed Al2O3 grains against the destruction in the shocked gas within the SN remnants.

The circumstellar (CS) environment is key to understanding progenitors of type Ia supernovae (SNe Ia), as well as the origin of a peculiar extinction property toward SNe Ia for cosmological application. It has been suggested that multiple scatterings of SN photons by CS dust may explain the non-standard reddening law. In this paper, we examine the effect of re-emission of SN photons by CS dust in the infrared (IR) wavelength regime. This effect allows the observed IR light curves to be used as a constraint on the position/size and the amount of CS dust. The method was applied to observed near-infrared (NIR) SN Ia samples; meaningful upper limits on the CS dust mass were derived even under conservative assumptions. We thereby clarify a difficulty associated with the CS dust scattering model as a general explanation for the peculiar reddening law, while it may still apply to a sub-sample of highly reddened SNe Ia. For SNe Ia in general, the environment at the interstellar scale appears to be responsible for the non-standard extinction law. Furthermore, deeper limits can be obtained using the standard nature of SN Ia NIR light curves. In this application, an upper limit of Mdot ~10^{-8}-10^{-7} Msun/yr (for the wind velocity of ~10 km/s) is obtained for a mass loss rate from a progenitor up to ~0.01 pc, and Mdot ~10^{-7}-10^{-6} Msun/yr up to ~0.1 pc.

We investigate the condition for the formation of low-mass second-generation stars in the early universe. It has been proposed that gas cooling by dust thermal emission can trigger fragmentation of a low-metallicity star-forming gas cloud. In order to determine the critical condition in which dust cooling induces the formation of low-mass stars, we follow the thermal evolution of a collapsing cloud by a one-zone semi-analytic collapse model. Earlier studies assume the dust amount in the local universe, where all refractory elements are depleted onto grains, and/or assume the constant dust amount during gas collapse. In this paper, we employ the models of dust formation and destruction in early supernovae to derive the realistic dust compositions and size distributions for multiple species as the initial conditions of our collapse calculations. We also follow accretion of heavy elements in the gas phase onto dust grains, i.e., grain growth, during gas contraction. We find that grain growth well alters the fragmentation property of the clouds, and that this still does not approach to the value in the local universe. The critical conditions can be written by the gas metallicity Zcr and the initial depletion efficiency fdep,0 of gas-phase metal onto grains, or dust-to-metal mass ratio, as (Zcr/10^{-5.5} Zsun) = (fdep,0/0.18)^{-0.44} with small scatters in the range of Zcr = [0.06--3.2]x10^{-5} Zsun. We also show that the initial dust composition and size distribution are important to determine Zcr.

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 investigate the origin of carbon-enhanced metal-poor (CEMP) stars starting from the recently discovered $\rm [Fe/H]<-7.1$ star SMSS J031300 (Keller et al. 2014). We show that the elemental abundances observed on the surface of SMSS J031300 can be well fit by the yields of faint, metal free, supernovae. Using properly calibrated faint supernova explosion models, we study, for the first time, the formation of dust grains in such carbon-rich, iron-poor supernova ejecta. Calculations are performed assuming both unmixed and uniformly mixed ejecta and taking into account the partial destruction by the supernova reverse shock. We find that, due to the paucity of refractory elements beside carbon, amorphous carbon is the only grain species to form, with carbon condensation efficiencies that range between (0.15-0.84), resulting in dust yields in the range (0.025-2.25)M$_{\odot}$. We follow the collapse and fragmentation of a star forming cloud enriched by the products of these faint supernova explosions and we explore the role played by fine structure line cooling and dust cooling. We show that even if grain growth during the collapse has a minor effect of the dust-to-gas ratio, due to C depletion into CO molecules at an early stage of the collapse, the formation of CEMP low-mass stars, such as SMSS J031300, could be triggered by dust cooling and fragmentation. A comparison between model predictions and observations of a sample of C-normal and C-rich metal-poor stars supports the idea that a single common pathway may be responsible for the formation of the first low-mass stars.

We present Hubble Space Telescope (HST) observations of the field of the Type IIb supernova (SN) 2011dh in M51 performed at ~1161 rest-frame days after explosion using the Wide Field Camera 3 and near-UV filters F225W and F336W. A star-like object is detected in both bands and the photometry indicates it has negative (F225W - F336W) color. The observed object is compatible with the companion of the now-vanished yellow supergiant progenitor predicted in interacting binary models. We consider it unlikely that the SN is undergoing strong interaction and thus estimate that it makes a small contribution to the observed flux. The possibilities of having detected an unresolved light echo or an unrelated object are briefly discussed and judged unlikely. Adopting a possible range of extinction by dust, we constrain parameters of the proposed binary system. In particular, the efficiency of mass accretion onto the binary companion must be below 50%, if no significant extinction is produced by newly formed dust. Further multiband observations are required in order to confirm the identification of the object as the companion star. If confirmed, the companion star would already be dominant in the UV/optical regime, so it would readily provide a unique opportunity to perform a detailed study of its properties.

We examine dust formation in macronovae (as known as kilonovae), which are the bright ejecta of neutron star binary mergers and one of the leading sites of r-process nucleosynthesis. We find that dust grains of r-process elements are difficult to form because of the low number density of the r-process atoms, while carbon or elements lighter than irons can condense into dust if they are abundant, in light of the first macronova candidate associated with GRB 130603B. Dust grains absorb emission from ejecta with opacity even greater than that of the r-process elements, and re-emit photons at infrared wavelengths. Such dust emission can potentially account for the macronova without r-process nucleosynthesis as an alternative model. This dust scenario predicts a more featureless spectrum than the r-process model and day-scale optical-to-ultraviolet emission.

We investigate the formation of dust in a stellar wind during the red-supergiant (RSG) phase of a very massive Population III star with the zero-age main sequence mass of 500 M_sun. We show that, in a carbon-rich wind with a constant velocity, carbon grains can form with a lognormal-like size distribution, and that all of the carbon available for dust formation finally condense into dust for wide ranges of the mass-loss rate ((0.1-3)x10^{-3} M_sun yr^{-1}) and wind velocity (1-100 km s^{-1}). We also find that the acceleration of the wind driven by newly formed dust suppresses the grain growth but still allows more than half of gas-phase carbon to be finally locked up in dust grains. These results indicate that at most 1.7 M_sun of carbon grains can form in total during the RSG phase of 500 M_sun Population III stars. Such a high dust yield could place very massive primordial stars as important sources of dust at the very early epoch of the universe if the initial mass function of Population III stars was top-heavy. We also briefly discuss a new formation scenario of carbon-rich ultra-metal-poor stars considering the feedback from very massive Population III stars.

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.

Dust grains in low-metallicity star-forming regions may be responsible for the formation of the first low-mass stars. The minimal conditions to activate dust-induced fragmentation require the gas to be pre-enriched above a critical dust-to-gas mass ratio Dcr=[2.6--6.3]x10^-9 with the spread reflecting the dependence on the grain properties. The recently discovered Galactic halo star SDSS J102915+172927 has a stellar mass of 0.8 Msun and a metallicity of Z=4.5x10^-5 Zsun and represents an optimal candidate for the dust-induced low-mass star formation. Indeed, for the two most plausible Population III supernova progenitors, with 20 Msun and 35 Msun, the critical dust-to-gas mass ratio can be overcome provided that at least 0.4 Msun of dust condenses in the ejecta, allowing for moderate destruction by the reverse shock. Here we show that even if dust formation in the first supernovae is less efficient or strong dust destruction does occur, grain growth during the collapse of the parent gas cloud is sufficiently rapid to activate dust cooling and likely fragmentation into low-mass and long-lived stars. Silicates and magnetite grains can experience significant grain growth in the density range 10^9 /cc < nH<10^12 /cc by accreting gas-phase species (SiO, SiO2, and Fe) until their gas-phase abundance drops to zero, reaching condensation efficiencies =1. The corresponding increase in the dust-to-gas mass ratio allows dust-induced cooling and fragmentation to be activated at 10^12 /cc < nH < 10^14 /cc, before the collapsing cloud becomes optically thick to continuum radiation. We show that for all the initial conditions that apply to the parent cloud of SDSS J102915+172927, dust-driven fragmentation is able to account for the formation of the star.

We report on X-ray spectral evolution of the nearby Type IIn supernova (SN) 2005ip, based on Chandra and Swift observations covering from ~1 to 6 years after the explosion. X-ray spectra in all epochs are well fitted by a thermal emission model with kT > 7 keV. The somewhat high temperature suggests that the X-ray emission mainly arises from the circumstellar medium heated by the forward shock. We find that the spectra taken 2-3 years since the explosion are heavily absorbed N_H ~ 5e22 cm^{-2}, but the absorption gradually decreases to the level of the Galactic absorption N_H ~ 4e20 cm^{-2} at the final epoch. This indicates that the SN went off in a dense circumstellar medium and that the forward shock has overtaken it. The intrinsic X-ray luminosity stays constant until the final epoch when it drops by a factor of ~2. The intrinsic 0.2-10 keV luminosity during the plateau phase is measured to be ~1.5e41 erg/s, ranking SN 2005ip as one of the brightest X-ray SNe. Based on the column density, we derive a lower-limit of a mass-loss rate to be M_dot ~ 0.015 (V_w/100 km/s) M_sun/yr, which roughly agrees with that inferred from the X-ray luminosity, M_dot ~ 0.02 (V_w/100 km/s) M_sun/yr, where V_w is the circumstellar wind speed. Such a high mass-loss rate suggests that the progenitor star had eruptive mass ejections like a luminous blue variable star. The total mass ejected in the eruptive period is estimated to be ~15 M_sun, indicating that the progenitor mass is greater than ~25 M_sun.

The non-steady-state formation of small clusters and the growth of grains accompanied by chemical reactions are formulated under the consideration that the collision of key gas species (key molecule) controls the kinetics of dust formation process. The formula allows us to evaluate the size distribution and condensation efficiency of dust formed in astrophysical environments. We apply the formulation to the formation of C and MgSiO3 grains in the ejecta of supernovae, as an example, to investigate how the non-steady effect influences the formation process, condensation efficiency f_{con}, and average radius a_{ave} of newly formed grains in comparison with the results calculated with the steady-state nucleation rate. We show that the steady-state nucleation rate is a good approximation if the collision timescale of key molecule tau_{coll} is much smaller than the timescale tau_{sat} with which the supersaturation ratio increases; otherwise the effect of the non-steady state becomes remarkable, leading to a lower f_{con} and a larger a_{ave}. Examining the results of calculations, we reveal that the steady-state nucleation rate is applicable if the cooling gas satisfies Lambda = tau_{sat}/tau_{coll} > 30 during the formation of dust, and find that f_{con} and a_{ave} are uniquely determined by Lambda_{on} at the onset time t_{on} of dust formation. The approximation formulae for f_{con} and a_{ave} as a function of Lambda_{on} could be useful in estimating the mass and typical size of newly formed grains from observed or model-predicted physical properties not only in supernova ejecta but also in mass-loss winds from evolved stars.

Supernovae (SNe) have been proposed to be the main production sites of dust grains in the Universe. Our knowledge on their importance to dust production is, however, limited by observationally poor constraints on the nature and amount of dust particles produced by individual SNe. In this paper, we present a spectrum covering optical through near-Infrared (NIR) light of the luminous Type IIn supernova (SN IIn) 2010jl around one and half years after the explosion. This unique data set reveals multiple signatures of newly formed dust particles. The NIR portion of the spectrum provides a rare example where thermal emission from newly formed hot dust grains is clearly detected. We determine the main population of the dust species to be carbon grains at a temperature of ~1,350 - 1,450K at this epoch. The mass of the dust grains is derived to be ~(7.5 - 8.5) x 10^{-4} Msun. Hydrogen emission lines show wavelength-dependent absorption, which provides a good estimate on the typical size of the newly formed dust grains (~0.1 micron, and most likely <~0.01 micron). We attribute the dust grains to have been formed in a dense cooling shell as a result of a strong SN-circumstellar media (CSM) interaction. The dust grains occupy ~10% of the emitting volume, suggesting an inhomogeneous, clumpy structure. The average CSM density is required to be >~3 x 10^{7} cm^{-3}, corresponding to a mass loss rate of >~0.02 Msun yr^{-1} (for a mass loss wind velocity of ~100 km s^{-1}). This strongly supports a scenario that SN 2010jl and probably other luminous SNe IIn are powered by strong interactions within very dense CSM, perhaps created by Luminous Blue Variable (LBV)-like eruptions within the last century before the explosion.

Modern data of the extinction curve from the ultraviolet to the near infrared are revisited to study the property of dust grains in the Milky Way (MW) and the Small Magellanic Cloud (SMC). We confirm that the graphite-silicate mixture of grains yields the observed extinction curve with the simple power-law distribution of the grain size but with a cutoff at some maximal size: the parameters are tightly constrained to be $q = 3.5 \pm 0.2$ for the size distribution $a^{-q}$ and the maximum radius $a_{max} = 0.24 \pm 0.05$ um, for both MW and SMC. The abundance of grains, and hence the elemental abundance, is constrained from the reddening versus hydrogen column density, E(B-V)/N_H. If we take the solar elemental abundance as the standard for the MW, >56 % of carbon should be in graphite dust, while it is <40 % in the SMC using its available abundance estimate. This disparity and the relative abundance of C to Si explain the difference of the two curves. We find that 50-60 % of carbon may not necessarily be in graphite but in the amorphous or the glassy phase. Iron may also be in the metallic phase or up to ~80 % in magnetite rather than in silicates, so that the Mg/Fe ratio in astronomical olivine is arbitrary. With these substitutions the parameters of the grain size remain unchanged. The mass density of dust grains relative to hydrogen is $\rho_{dust}/\rho_H = 1/(120 {+10 \atop -16})$ for the MW and $1/(760 {+70 \atop -90}) for the SMC under some abundance constraints. We underline the importance of the wavelength-dependence slope of the extinction curve in the near infrared in constructing the dust model: if $A_{\lambda} \propto \lambda^{-gamma}$ with gamma ~ 1.6, the power-law grain-size model fails, whereas it works if gamma ~ 1.8-2.0.

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.

We present a study of SN 2009js in NGC 918. Multi-band Kanata optical photometry covering the first ~120 days show the source to be a Type IIP SN. Reddening is dominated by that due to our Galaxy. One-year-post-explosion photometry with the NTT, and a Subaru optical spectrum 16 days post-discovery, both imply a good match with the well-studied subluminous SN 2005cs. The plateau phase luminosity of SN 2009js and its plateau duration are more similar to the intermediate luminosity IIP SN 2008in. Thus, SN 2009js shares characteristics with both subluminous and intermediate luminosity SNe. Its radioactive tail luminosity lies between SN 2005cs and SN 2008in, whereas its quasi-bolometric luminosity decline from peak to plateau (quantified by a newly-defined parameter Delta[logL] measuring adiabatic cooling following shock breakout) is much smaller than both the others. We estimate the ejected mass of 56Ni to be low (~0.007 Msun). The SN explosion energy appears to have been small, similar to SN 2005cs. SN 2009js is the first subluminous SN IIP to be studied in the mid-infrared. It was serendipitously caught by Spitzer at very early times. In addition, it was detected by WISE 105 days later with a significant 4.6 micron flux excess above the photosphere. The infrared excess luminosity relative to the photosphere is clearly smaller than that of SN 2004dj extensively studied in the mid-infrared. The excess may be tentatively assigned to heated dust with mass ~3e-5 Msun, or to CO fundamental emission as a precursor to dust formation.

In a low-metallicity gas, rapid cooling by dust thermal emission is considered to induce cloud fragmentation and play a vital role in the formation of low-mass stars (<~ 1 M_sun) in metal-poor environments. We investigate how the growth of dust grains through accretion of heavy elements in the gas phase onto grain surfaces alters the thermal evolution and fragmentation properties of a collapsing gas cloud. We calculate directly grain growth and dust emission cooling in a self-consistent manner. We show that MgSiO3 grains grow sufficiently at gas densities nH = 10^{10}, 10^{12}, and 10^{14} /cc for metallicities Z = 10^{-4}, 10^{-5}, and 10^{-6} Zsun, respectively, where the cooling of the collapsing gas cloud is enhanced. The condition for efficient dust cooling is insensitive to the initial condensation factor of pre-existing grains within the realistic range of 0.001--0.1, but sensitive to metallicity. The critical metallicity is Zcrit ~ 10^{-5.5} Zsun for the initial grain radius r_{MgSiO3,0} <~ 0.01 um and Zcrit ~ 10^{-4.5} Zsun for r_{MgSiO3,0} >~ 0.1 um. The formation of a recently discovered low-mass star with extremely low metallicity (<= 4.5x10^{-5} Zsun) could have been triggered by grain growth.

The discovery of a low-mass star with such low metallicity as < 4.5x10^{-5} Z_sun reveals the critical role of dust in the formation of extremely metal-poor stars. In this paper we explore the effect of the growth of dust grains through accretion of gaseous refractory elements in very low-metallicity pre-stellar cores on the cloud fragmentation induced by the dust emission cooling. Employing a simple model of grain growth in a gravitationally collapsing gas, we show that Fe and Si grains can grow efficiently at hydrogen densities of ~10^{10}-10^{14} cm^{-3} in the clouds with metal abundances of -5 <~ [Fe, Si/H] <~ -3. The critical metal number abundances, above which the grain growth could induce the fragmentation of the gas clouds, are estimated to be A_{crit} ~ 10^{-9}-10^{-8}, unless the initial grain radius is too large (>~ 1 um) or the sticking probability is too small (<~ 0.01). We find that even if the initial dust-to-gas mass ratio is well below the minimum value required for the dust-induced fragmentation, the grain growth increases the dust mass high enough to cause the gas fragmentation into sub-solar mass clumps. We suggest that as long as the critical metal abundance is satisfied, the grain growth could play an important role in the formation of low-mass stars with metallicity as low as 10^{-5} Z_sun.

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.

Most of the observational studies of supernova (SN) explosions are limited to early phases (< a few yr after the explosion) of extragalactic SNe and observations of SN remnants (> 100 yr) in our Galaxy or very nearby galaxies. SNe at the epoch between these two, which we call "transitional" phase, have not been explored in detail except for several extragalactic SNe including SN 1987A in the Large Magellanic Cloud. We present theoretical predictions for the infrared (IR) dust emissions by several mechanisms; emission from dust formed in the SN ejecta, light echo by circumstellar and interstellar dust, and emission from shocked circumstellar dust. We search for IR emission from 6 core-collapse SNe at the transitional phase in the nearby galaxies NGC 1313, NGC 6946, and M101 by using the data taken with the AKARI satellite and Spitzer. Among 6 targets, we detect the emission from SN 1978K in NGC 1313. SN 1978K is associated with 1.3 x 10^{-3} Msun of silicate dust. We show that, among several mechanisms, the shocked circumstellar dust is the most probable emission source to explain the IR emission observed for CSM-rich SN 1978K. IR emission from the other 5 objects is not detected. Our current observations are sensitive to IR luminosity of > 10^{38} erg s^{-1}, and the non-detection of SN 1962M excludes the existence of the shocked circumstellar dust for a high gas mass-loss rate of sim 10^{-4} Msun yr^{-1}. Observations of SNe at the transitional phase with future IR satellites will fill the gap of IR observations of SNe with the age of 10-100 years, and give a new opportunity to study the circumstellar and interstellar environments of the progenitor, and possibly dust formation in SNe.

The origin of interstellar dust in galaxies is poorly understood, particularly the relative contributions from supernovae and the cool stellar winds of low-intermediate mass stars. Here, we present Herschel PACS and SPIRE photometry at 70-500um of the historical young supernova remnants: Kepler and Tycho; both thought to be the remnants of Type Ia explosion events. We detect a warm dust component in Kepler's remnant with T = 82K and mass 0.0031Msun; this is spatially coincident with thermal X-ray emission optical knots and filaments, consistent with the warm dust originating in the circumstellar material swept up by the primary blast wave of the remnant. Similarly for Tycho's remnant, we detect warm dust at 90K with mass 0.0086Msun. Comparing the spatial distribution of the warm dust with X-rays from the ejecta and swept-up medium, and Ha emission arising from the post-shock edge, we show that the warm dust is swept up interstellar material. We find no evidence of a cool (25-50 K) component of dust with mass >0.07Msun as observed in core-collapse remnants of massive stars. Neither the warm or cold dust components detected here are spatially coincident with supernova ejecta material. We compare the lack of observed supernova dust with a theoretical model of dust formation in Type Ia remnants which predicts dust masses of 0.088(0.017)Msun for ejecta expanding into surrounding densities of 1(5)cm-3. The model predicts that silicon- and carbon-rich dust grains will encounter the interior edge of the observed dust emission at 400 years confirming that the majority of the warm dust originates from swept up circumstellar or interstellar grains (for Kepler and Tycho respectively). The lack of cold dust grains in the ejecta suggests that Type Ia remnants do not produce substantial quantities of iron-rich dust grains and has important consequences for the 'missing' iron mass observed in ejecta.

We analyze the multi-wavelength photometric and spectroscopic data of 12 ultraluminous infrared galaxies (ULIRGs) at z ~ 1 and compare them with models of stars and dust in order to study the extinction law and star formation in young infrared (IR) galaxies. Five extinction curves, namely, the Milky Way (MW), the pseudo MW which is MW-like without the 2175 Angstrom feature, the Calzetti, and two SN dust curves, are applied to the data, by combining with various dust distributions, namely, the uniform dust screen, the clumpy dust screen, the internal dust geometry, and the composite geometry with a combination of dust screen and internal dust. Employing a minimum chi square method, we find that the foreground dust screen geometry, especially combined with the 8 - 40 M_sun SN extinction curve, provides a good approximation to the real dust geometry, whereas internal dust is only significant in 2 galaxies. The SN extinction curves, which are flatter than the others, reproduce the data of 8(67%) galaxies better. Dust masses are estimated to be in excess of ~ 10^8 M_sun. Inferred ages of the galaxies are very young, 8 of which range from 10 to 650 Myr. The SN-origin dust is the most plausible to account for the vast amount of dust masses and the flat slope of the observed extinction law. The inferred dust mass per SN ranges from 0.01 to 0.4 M_sun/SN.

We study the role of the unknown microphysical properties of carbonaceous dust particles in determining the amount and size distribution of carbonaceous dust condensed in type II supernova explosions. We parametrize the microphysical properties in terms of the shape factor of the grain and the sticking coefficient of gas-phase carbon atoms onto the grain surfaces. We find that the amount of dust formed is fairly independent of these properties, within the parameter range considered, though limited by the available amount of carbon atoms not locked in CO molecules. However, we find that the condensation times and size distributions of dust grains depend sensitively on the microphysical parameters, with the mass distributions being weighted toward larger effective radii for conditions considering grains with higher sticking coefficients and/or more aspherical shapes. We discuss that this leads to important consequences on the predicted extinction law of SN dust and on the survival rate of the formed grains as they pass through the reverse shock of the SN. We conclude that a more detailed understanding of the dust formation process and of the microphysical properties of each dust species needs to be achieved before robust prediction on the SN dust yields can be performed.

We report far-infrared and submillimeter observations of Supernova 1987A, the star that exploded on February 23, 1987 in the Large Magellanic Cloud, a galaxy located 160,000 lightyears away. The observations reveal the presence of a population of cold dust grains radiating with a temperature of ~17-23 K at a rate of about 220 solar luminosity. The intensity and spectral energy distribution of the emission suggests a dust mass of ~0.4-0.7 solar mass. The radiation must originate from the SN ejecta and requires the efficient precipitation of all refractory material into dust. Our observations imply that supernovae can produce the large dust masses detected in young galaxies at very high redshifts.

We investigate the formation of dust grains in the ejecta of Type Ia supernovae (SNe Ia), adopting the carbon-deflagration W7 model. In the calculations of dust formation, we apply the nucleation and grain growth theory and consider the two cases with and without formation of CO and SiO molecules. The results of the calculations show that for the sticking probability of alpha_j=1, C, silicate, Si, and FeS grains can condense at early times of ~100--300 days after the explosion, whereas Fe and SiC grains cannot form substantially. Due to the low gas density in SNe Ia with no H-envelope, the average radii of the newly formed grains are generally below 0.01 micron, being much smaller than those in Type II-P SNe. This supports our previous conclusion that the radius of dust formed in the ejecta is smaller in SNe with less massive envelopes. The total dust mass ranges from 3x10^{-4} M_sun to 0.2 M_sun for alpha_j=0.1--1, depending on whether or not CO and SiO molecules are formed. We also estimate the optical depths and thermal emission by the newly formed dust and compare to the relevant observations of SNe Ia. We find that the formation of C grains in SNe Ia is suppressed to be consistent with the observational constraints. This implies that energetic photons and electrons heavily depress the formation efficiency of C grains or that the outermost C-O layer of SNe Ia is almost fully burned. Finally, we perform the calculations of dust destruction in the SN remnants and find that dust grains formed in the ejecta of SNe Ia are almost completely destroyed in the shocked gas before being injected into the interstellar medium. This indicates that SNe Ia are unlikely to be the major sources of interstellar dust.

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.

We present the results of near-infrared (NIR) multi-epoch observations of the optical transient in the nearby galaxy NGC300 (NGC300-OT) at 398 and 582 days after the discovery with the Infrared Camera (IRC) onboard AKARI. NIR spectra (2--5 um) of NGC300-OT were obtained for the first time. They show no prominent emission nor absorption features, but are dominated by continuum thermal emission from the dust around NGC300-OT. NIR images were taken in the 2.4, 3.2, and 4.1 um bands. The spectral energy distributions (SED) of NGC300-OT indicate the dust temperature of 810 (+-14) K at 398 days and 670 (+-12) K at 582 days. We attribute the observed NIR emission to the thermal emission from dust grains formed in the ejecta of NGC300-OT. The multi-epoch observations enable us to estimate the dust optical depth as larger than about 12 at 398 days and larger than about 6 at 582 days at 2.4 um, by assuming an isothermal dust cloud. The observed NIR emission must be optically thick, unless the amount of dust grains increases with time. Little extinction at visible wavelengths reported in earlier observations suggests that the dust cloud around NGC300-OT should be distributed inhomogeneously so as to not screen the radiation from the ejecta gas and the central star. The present results suggest the dust grains are not formed in spherically symmetric geometry, but rather in a torus, a bipolar outflow, or clumpy cloudlets.

We investigate the formation of dust grains in the ejecta of a SN IIb and their evolution in the shocked gas in the SNR by considering the uniform and power-law density structures for the CSM. Based on these calculations, we also simulate the time evolution of thermal emission from the shock-heated dust in the SNR and compare the results with the observations of Cas A SNR. We find that the total mass of dust formed in the SN IIb is as large as 0.167 M_sun but the average radius of dust is smaller than 0.01 mum and is significantly different from those in SNe II-P with the massive H-envelope. In the explosion with the small-mass H-envelope, the expanding He core undergoes little deceleration, so that the gas density in the He core is too low for large-sized grains to form. In addition, the low-mass H-envelope of the SN IIb leads to the early arrival of the reverse shock at the dust-forming region. If the CSM is more or less spherical, therefore, the newly formed grains would be completely destroyed in the relatively dense shocked gas for the CSM density of n_H > 0.1 cm^{-3}. However, the actual CSM is likely to be non-spherical, so that a part of grains could be ejected into the ISM without being shocked. We demonstrate that the time evolution of the SED by thermal dust emission is sensitive to the ambient gas density and structure that affects the passage of the reverse shock into the ejecta. Thus, the SED evolution well reflects the evolution of dust through erosion by sputtering and stochastic heating. For Cas A, we consider the CSM produced by the steady mass loss of ~8x10^{-5} M_sun/yr during the supergiant phase. Then we find the infrared SED of Cas A is reasonably reproduced by thermal emission from the newly formed dust of 0.08 M_sun, which consists of 0.008 M_sun shocked warm dust and 0.072 M_sun unshocked cold dust.

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.

We present the results of AKARI observations of the O-rich supernova remnant G292.0+1.8 using six IRC and four FIS bands covering 2.7-26.5 um and 50-180 um, respectively. The AKARI images show two prominent structures; a bright equatorial ring structure and an outer elliptical shell structure. The equatorial ring structure is clumpy and incomplete with its western end opened. The outer shell is almost complete and slightly squeezed along the north-south direction. The central position of the outer shell is ~ 1' northwest from the embedded pulsar and coincides with the center of the equatorial ring structure. The equatorial ring and the elliptical shell structures were partly visible in optical and/or X-rays, but they are much more clearly revealed in our AKARI images. There is no evident difference in infrared colors of the two prominent structures, which is consistent with the previous proposition that both structures are of circumstellar origin. However, we have detected faint infrared emission of a considerably high 15 to 24 um ratio associated with the supernova ejecta in the southeastern and northwestern areas. Our IRC spectra show that the high ratio is at least partly due to the emission lines from Ne ions in the supernova ejecta material. In addition we detect a narrow, elongated feature outside the SNR shell. We derive the physical parameters of the infrared-emitting dust grains in the shocked circumstellar medium and compare the result with model calculations of dust destruction by a SN shock. The AKARI results suggest that the progenitor was at the center of the infrared circumstellar shell in red supergiant stage and that the observed asymmetry in the SN ejecta could be a result of either a dense circumstellar medium in the equatorial plane and/or an asymmetric explosion.

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.

Core--collapsed supernovae (CCSNe) have been considered to be one of sources of dust in the universe. What kind and how much mass of dust are formed in the ejecta and are injected into the interstellar medium (ISM) depend on the type of CCSNe, through the difference in the thickness (mass) of outer envelope. In this review, after summarizing the existing results of observations on dust formation in CCSNe, we investigate formation of dust in the ejecta and its evolution in the supernova remnants (SNRs) of Type II--P and Type IIb SNe. Then, the time evolution of thermal emission from dust in the SNR of Type IIb SN is demonstrated and compared with the observation of Cas A. We find that the total dust mass formed in the ejecta does not so much depend on the type; $\sim 0.3-0.7 M_{\odot}$ in Type II--P SNe and $\sim 0.13 M_{\odot}$ in Type IIb SN. However the size of dust sensitively depends on the type, being affected by the difference in the gas density in the ejecta: the dust mass is dominated by grains with radii larger than 0.03 $\mu$m in Type II-P, and less than 0.006 $\mu$m in Type IIb, which decides the fate of dust in the SNR. The surviving dust mass is $\sim 0.04-0.2 M_{\odot}$ in the SNRs of Type II--P SNe for the ambient hydrogen density of $n_{\rm H}=10.0-1.0$ cm$^{-3}$, while almost all dust grains are destroyed in the SNR of Type IIb. The spectral energy distribution (SED) of thermal emission from dust in SNR well reflects the evolution of dust grains in SNR through erosion by sputtering and stochastic heating. The observed SED of Cas A SNR is reasonably reproduced by the model of dust formation and evolution for Type IIb SN.

We investigate the evolution of dust formed in Population III supernovae (SNe) by considering its transport and processing by sputtering within the SN remnants (SNRs). We find that the fates of dust grains within SNRs heavily depend on their initial radii $a_{\rm ini}$. For Type II SNRs expanding into the ambient medium with density of $n_{\rm H,0} = 1$ cm$^{-3}$, grains of $a_{\rm ini} < 0.05$ $\mu$m are detained in the shocked hot gas and are completely destroyed, while grains of $a_{\rm ini} > 0.2$ $\mu$m are injected into the surrounding medium without being destroyed significantly. Grains with $a_{\rm ini}$ = 0.05-0.2 $\mu$m are finally trapped in the dense shell behind the forward shock. We show that the grains piled up in the dense shell enrich the gas up to 10$^{-6}-10^{-4}$ $Z_\odot$, high enough to form low-mass stars with 0.1-1 $M_\odot$. In addition, [Fe/H] in the dense shell ranges from -6 to -4.5, which is in good agreement with the ultra-metal-poor stars with [Fe/H] < -4. We suggest that newly formed dust in a Population III SN can have great impacts on the stellar mass and elemental composition of Population II.5 stars formed in the shell of the SNR.

We present a theoretical model for Type Ib supernova (SN) 2006jc. We calculate the evolution of the progenitor star, hydrodynamics and nucleosynthesis of the SN explosion, and the SN bolometric light curve (LC). The synthetic bolometric LC is compared with the observed bolometric LC constructed by integrating the UV, optical, near-infrared (NIR), and mid-infrared (MIR) fluxes. The progenitor is assumed to be as massive as $40M_\odot$ on the zero-age main-sequence. The star undergoes extensive mass loss to reduce its mass down to as small as $6.9M_\odot$, thus becoming a WCO Wolf-Rayet star. The WCO star model has a thick carbon-rich layer, in which amorphous carbon grains can be formed. This could explain the NIR brightening and the dust feature seen in the MIR spectrum. We suggest that the progenitor of SN 2006jc is a WCO Wolf-Rayet star having undergone strong mass loss and such massive stars are the important sites of dust formation. We derive the parameters of the explosion model in order to reproduce the bolometric LC of SN 2006jc by the radioactive decays: the ejecta mass $4.9M_\odot$, hypernova-like explosion energy $10^{52}$ ergs, and ejected $^{56}$Ni mass $0.22M_\odot$. We also calculate the circumstellar interaction and find that a CSM with a flat density structure is required to reproduce the X-ray LC of SN 2006jc. This suggests a drastic change of the mass-loss rate and/or the wind velocity that is consistent with the past luminous blue variable (LBV)-like event.

SN 2006jc is a peculiar supernova (SN), in which the formation of dust has been confirmed at an early epoch of ~50 days after the explosion. We investigate the possibility of such an earlier formation of dust grains in the expanding ejecta of SN 2006jc, applying the Type Ib SN model that is developed to reproduce the observed light curve. We find that the rapid decrease of the gas temperature in SN 2006jc enables the condensation of C grains in the C-rich layer at 40-60 days after the explosion, which is followed by the condensation of silicate and oxide grains until ~200 days. The average radius of each grain species is confined to be less than 0.01 micron due to the low gas density at the condensation time. The calculated total dust mass reaches ~1.5 Msun, of which C dust shares 0.7 Msun. On the other hand, based on the calculated dust temperature, we show that the dust species and mass evaluated to reproduce the spectral energy distribution observed by AKARI and MAGNUM at day 200 are different from those obtained by the dust formation calculations; the dust species contributing to the observed flux are hot C and FeS grains with masses of $5.6 \times 10^{-4}$ Msun and $2.0 \times 10^{-3}$ Msun, respectively, though we cannot defy the presence of a large amount of cold dust such as silicate and oxide grains up to 0.5 Msun. One of the physical processes responsible for the difference between calculated and evaluated masses of C and FeS grains could be considered to be the destruction of small-sized clusters by energetic photons and electrons prevailing within the ejecta at the earlier epoch.

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.

We present our latest results on near- to mid- infrared observation of SN2006jc at 200 days after the discovery using the Infrared Camera (IRC) on board $AKARI$. The near-infrared (2--5$\mu$m) spectrum of SN2006jc is obtained for the first time and is found to be well interpreted in terms of the thermal emission from amorphous carbon of 800$\pm 10$K with the mass of $6.9\pm 0.5 \times 10^{-5}M_{\odot}$ that was formed in the supernova ejecta. This dust mass newly formed in the ejecta of SN 2006jc is in a range similar to those obtained for other several dust forming core collapse supernovae based on recent observations (i.e., $10^{-3}$--$10^{-5}$$M_{\odot}$). Mid-infrared photometric data with {\it{AKARI}}/IRC MIR-S/S7, S9W, and S11 bands have shown excess emission over the thermal emission by hot amorphous carbon of 800K. This mid-infrared excess emission is likely to be accounted for by the emission from warm amorphous carbon dust of 320$\pm 10$K with the mass of 2.7$^{+0.7}_{-0.5} \times 10^{-3}M_{\odot}$ rather than by the band emission of astronomical silicate and/or silica grains. This warm amorphous carbon dust is expected to have been formed in the mass loss wind associated with the Wolf-Rayet stellar activity before the SN explosion. Our result suggests that a significant amount of dust is condensed in the mass loss wind prior to the SN explosion. A possible contribution of emission bands by precursory SiO molecules in 7.5--9.5$\mu$m is also suggested.

We investigate the evolution of dust that formed at Population III supernova (SN) explosions and its processing through the collisions with the reverse shocks resulting from the interaction of the SN ejecta with the ambient medium. In particular, we investigate the transport of the shocked dust within the SN remnant (SNR), and its effect on the chemical composition, the size distribution, and the total mass of dust surviving in SNR. We find that the evolution of the reverse shock, and hence its effect on the processing of the dust depends on the thickness of the envelope retained by the progenitor star. Furthermore, the transport and survival of the dust grains depend on their initial radius, a_{ini}, and composition: For Type II SNRs expanding into the interstellar medium (ISM) with a density of n_{H,0}=1 cm^{-3}, small grains with a_{ini} < ~ 0.05 micron are completely destroyed by sputtering in the postshock flow, while grains with a_{ini}= 0.05--0.2 micron are trapped into the dense shell behind the forward shock. Very large grains of a_{ini} > ~ 0.2 micron are ejected into the ISM without decreasing their sizes significantly. We find that the total mass fraction of dust that is destroyed by the reverse shock ranges from 0.2 to 1.0, depending on the energy of the explosion and the density of the ambient ISM. The results of our calculations have significant impact on the abundance pattern of subsequent generation of stars that form in the dense shell of primordial SNRs.

We investigate the destruction of dust grains by sputtering in the high-velocity interstellar shocks driven by supernovae (SNe) in the early universe to reveal the dependence of the time-scale of dust destruction on the gas density $n_{{\rm H}, 0}$ in the interstellar medium (ISM) as well as on the progenitor mass $M_{\rm pr}$ and explosion energy $E_{\rm 51}$ of SN. The sputtering yields for the combinations of dust and ion species of interest to us are evaluated by applying the so-called universal relation with a slight modification. The dynamics of dust grains and their destruction by sputtering in shock are calculated by taking into account the size distribution of each dust species, together with the time evolution of temperature and density of gas in spherically symmetric shocks. The results of calculations show that the efficiency of dust destruction depends not only on the sputtering yield but also on the initial size distribution of each grain species. The efficiency of dust destruction increases with increasing $E_{\rm 51}$ and/or increasing $n_{{\rm H}, 0}$, but is almost independent of $M_{\rm pr}$ as long as $E_{\rm 51}$ is the same. The mass of gas swept up by shock is the increasing function of $E_{\rm 51}$ and the decreasing function of $n_{{\rm H}, 0}$. Combining these results, we present the approximation formula for the time-scale of destruction for each grain species in the early universe as a function of $E_{\rm 51}$ and $n_{{\rm H}, 0}$. This formula is applicable for investigating the evolution of dust grains at the early epoch of the universe with the metallicity of $Z \la 10^{-3}$ $Z_\odot$. The effects of the cooling processes of gas on the destruction of dust are briefly discussed.

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.

We investigate the formation of dust grains in the ejecta of population III supernovae including pair--instability supernovae, applying a theory of non-- steady state nucleation and grain growth. In the calculations, the time evolution of gas temperature in theejecta, which strongly affects the number density and size of newly formed grains, is calculated by solving the radiative transfer equation taking account of the energy deposition of radio active elements. Two extreme cases are considered for the mixing of elements in the ejecta; unmixed and uniformly mixed cases within the He--core. The results of calculations are summarized as the followings; in the unmixed ejecta, a variety of grain species condense, reflecting the difference of the elemental composition at the formation site in the ejecta, otherwise only oxide grains condense in the uniformly mixed ejecta. The average size of newly formed grains spans the range of three orders of magnitude, depending on the grain species and the formation condition, and the maximum radius is limited to less than 1 $\mu$m, which does not depend on the progenitor mass. The size distribution function summed up over all grain species is approximated by a power--law formula whose index is -3.5 for the larger radius and -2.5 for the smaller one; the radius at the crossover point ranges from 0.004 to 0.1 $\mu$m, depending on the model of supernovae. The fraction of mass locked into dust grains increases with increasing the progenitor mass; 2--5 % of the progenitor mass for core collapse supernovae and 15--30 % for pair--instability supernovae whose progenitor mass ranges from 140 to 260 $M_{\odot}$. Thus, if the very massive stars populate the first generation stars, a large amount of dust grains would be produced in the early universe.

We conduct a direct comparison of three different representative numerical codes for constructing models of rapidly rotating neutron stars in general relativity. Our aim is to evaluate the accuracy of the codes and to investigate how the accuracy is affected by the choice of interpolation, domain of integration and equation of state. In all three codes, the same physical parameters, equations of state and interpolation method are used. We construct 25 selected models for polytropic equations of state and 22 models with realistic neutron star matter equations of state. The three codes agree well with each other (typical agreement is better than 0.1% to 0.01%) for most models, except for the extreme assumption of uniform density stars. We conclude that the codes can be used for the construction of highly accurate initial data configurations for polytropes of index N > 0.5 (which typically correspond to realistic neutron stars), when the domain of integration includes all space and for realistic equations with no phase transitions. With the exception of the uniform density case, the obtained values of physical parameters for the models considered in this paper can be regarded as ``standard'' and we display them in detail for all models.