Dust formation in primordial Type II supernovae

We have investigated the formation of dust in the ejecta of Type II supernovae (SNe), mostly of primordial composition, to answer the question of where the first solid particles are formed in the Universe. However, we have also considered non-zero progenitor metallicity values up to Z = Z_⊙ . The calculations are based on standard nucleation theory, and the scheme has been tested for the first time on the well-studied case of SN1987A, yielding results that are in agreement with the available data. We find that: (i) the first dust grains are predominantly made of silicates, amorphous carbon (AC), magnetite and corundum; and (ii) the largest grains are the AC ones, with sizes around 300 Å, whereas the other grain types have smaller radii, around 10–20 Å . The grain size distribution depends somewhat on the thermodynamics of the ejecta expansion, and variations in the results by a factor ≈2 might occur within reasonable estimates of the relevant parameters. Also, and for the same reason, the grain size distribution is essentially unaffected by metallicity changes. The predictions on the amount of dust formed are very robust: for Z =0 , we find that SNe with masses in the range (12–35) M_⊙ produce about 0.08 M_⊙≲ M _d≲0.3 M_⊙ of dust per supernova. The above range increases by roughly three times as the metallicity is increased to solar values. We discuss the implications and the cosmological consequences of the results.     

Joint formation of QSOs and spheroids: QSOs as clocks of star formation in spheroids

Direct and indirect observational evidence leads to the conclusion that high-redshift QSOs did shine in the core of early-type protogalaxies during their main episode of star formation. Exploiting this fact, we derive the rate of formation of this kind of stellar system at high redshift by using the QSO luminosity function. The elemental proportions in elliptical galaxies, the descendants of the QSO hosts, suggest that the star formation was more rapid in more massive objects. We show that this is expected to occur in dark matter haloes, when the processes of cooling and heating are considered. This is also confirmed by comparing the observed submm counts with those derived by coupling the formation rate and the star formation rate of the spheroidal galaxies with a detailed model for their SED evolution. In this scenario SCUBA galaxies and Lyman-break galaxies are early-type protogalaxies forming the bulk of their stars before the onset of QSO activity.     

Extinction curves flattened by reverse shocks in supernovae

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⁻³  produces too flat an extinction curve to be consistent with the observed extinction curve for SDSS J1048+4637 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.     

A strategy for finding gravitationally lensed distant supernovae

Distant Type Ia and II supernovae (SNe) can serve as valuable probes of the history of the cosmic expansion and star formation, and provide important information on their progenitor models. At present, however, there are few observational constraints on the abundance of SNe at high redshifts. A major science driver for the Next Generation Space Telescope is the study of such very distant SNe. In this paper we discuss strategies for finding and counting distant SNe by using repeat imaging of supercritical intermediate redshift clusters whose mass distributions are well constrained via modelling of strongly lensed features. For a variety of different models for the star formation history and supernova progenitors, we estimate the likelihood of detecting lensed SNe as a function of their redshift. In the case of a survey conducted with Hubble Space Telescope ( HST ), we predict a high probability of seeing a supernova in a single return visit with either Wide Field Planetary Camera 2 or Advanced Camera for Surveys, and a much higher probability of detecting examples with     in the lensed case. Most events would represent magnified SNe II at     and a fraction will be more distant examples. We discuss various ways to classify such events using ground-based infrared photometry. We demonstrate an application of the method using the HST archival data and discuss the case of a possible event found in the rich cluster AC 114     .     

Extinction correction for Type Ia supernova rates – I. The model

In this paper, we present and discuss a new Monte Carlo approach aimed at correcting the observed supernova (SN) rates for the effects of host galaxy dust extinction. The problem is addressed in a general way and the model includes SN position distributions, SN light-curve and spectral library and dust properties and distribution as input ingredients. Even though the recipe we propose is in principle applicable to all SN types, in this paper, we illustrate the use of our model only for Type Ia. These represent, in fact, the simplest test case, basically due to their spectroscopic homogeneity, which to a first approximation allows one to treat them all in the same way. This test case shows that the final results do not depend critically on the spiral arm dust geometry, while the total amount of dust, its properties and the size of the Galactic bulge do have a strong effect. With the availability of more complete spectral libraries and a more accurate knowledge of SN spatial distribution, the method we propose here can be easily extended to core collapse events.     

Dust coagulation in star formation with different metallicities

Dust grains coagulate into larger aggregates in dense gas. This changes their size distribution and possibly affects the thermal evolution of star-forming clouds. We here investigate dust coagulation in collapsing pre-stellar cores with different metallicities by considering the thermal motions of grains. We show that coagulation does occur even at low metallicity  ∼10⁻⁶ Z_⊙  . However, we also find (i) that the H₂ formation rate on dust grains is reduced only after the majority of H₂ 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.     

Constraints on cosmological anisotropy out to z= 1 from Type Ia supernovae

A combined sample of 79 high- and low-redshift Type Ia supernovae (SNe) is used to set constraints on the degree of anisotropy in the Universe out to z ≃1. First, we derive the global most probable values of matter density Ω_M, the cosmological constant Ω_Λ and the Hubble constant H ₀, and find them to be consistent with the published results from the two data sets of Riess et al. and Perlmutter et al. We then examine the Hubble diagram (HD, i.e., the luminosity–redshift relation) in different directions on the sky by utilizing spherical harmonic expansion. In particular, via the analysis of the dipole anisotropy, we divide the sky into the two hemispheres that yield the most discrepant of the three cosmological parameters, and the scatter χ _HD² in each case. The most discrepant values roughly move along the locus −4Ω_M+3Ω_Λ=1 (cf. Perlmutter et al.), but by no more than Δ≈2.5 along this line. For a perfect Friedmann–Robertson–Walker universe, Monte Carlo realizations that mimic the current set of SNe yield values higher than the measured Δ in ∼1/5 of the cases (for Ω_M). We discuss implications for the validity of the Cosmological Principle, and possible calibration problems in the SNe data sets.     

The correlation between black hole mass and bulge velocity dispersion in hierarchical galaxy formation models

Popular models for the origin of gamma-ray bursts (GRBs) include short-lived massive stars as the progenitors of the fireballs. Hence the redshift distribution of GRBs should track the cosmic star formation rate of massive stars accurately. A significant proportion of high-mass star formation activity appears to occur in regions that are obscured from view in the optical waveband by interstellar dust. The amount of dust-enshrouded star formation activity taking place has been estimated by observing the thermal radiation from the dust that has been heated by young stars in the far-infrared and submillimetre wavebands. Here we discuss an alternative probe – the redshift distribution of GRBs. GRBs are detectable at the highest redshifts, and because gamma-rays are not absorbed by dust, the redshift distribution of GRBs should therefore be unaffected by dust extinction. At present the redshifts of GRBs can only be determined from the associated optical transient emission; however, useful information about the prevalence of dust-obscured star formation can also be obtained from the ratio of GRBs with and without an associated optical transient. Eight GRBs currently have spectroscopic redshifts. Once about a hundred redshifts are known, the population of GRBs will provide an important test of different models of the star formation history of the Universe.     

Dust formation and survival in supernova ejecta

The presence of dust at high redshift requires efficient condensation of grains in supernova (SN) ejecta, in accordance with current theoretical models. Yet observations of the few well-studied supernovae (SNe) and supernova remnants (SNRs) imply condensation efficiencies which are about two orders of magnitude smaller. Motivated by this tension, we have (i) revisited the model of Todini & Ferrara for dust formation in the ejecta of core collapse SNe, and (ii) followed, for the first time, the evolution of newly condensed grains from the time of formation to their survival – through the passage of the reverse shock – in the SNR. We find that  0.1–0.6  M_⊙  of dust form in the ejecta of 12–40 M_⊙ stellar progenitors. Depending on the density of the surrounding interstellar medium, between 2 and 20 per cent of the initial dust mass survives the passage of the reverse shock, on time-scales of about  4–8 × 10⁴  yr  from the stellar explosion. Sputtering by the hot gas induces a shift of the dust size distribution towards smaller grains. The resulting dust extinction curve shows a good agreement with that derived by observations of a reddened QSO at   z = 6.2  . Stochastic heating of small grains leads to a wide distribution of dust temperatures. This supports the idea that large amounts (∼0.1 M_⊙) of cold dust  ( T ∼ 40   K)  can be present in SNRs, without being in conflict with the observed infrared emission.     

On core-collapse supernovae in normal and in Seyfert galaxies

This paper estimates the relative frequency of different types of core-collapse supernovae, in terms of the ratio between the number of Type Ib–Ic and of Type II supernovae. We estimate independently for all normal and Seyfert galaxies whose radial velocity is ≤14 000 km s⁻¹, and which had at least one supernova event recorded in the Asiago catalogue from 1986 January to 2000 August. We find that the ratio is  ≈0.23±0.05  in normal galaxies. This value is consistent with constant star formation rate and with a Salpeter initial mass function and an average binary rate ≈50 per cent. On the contrary, Seyfert galaxies exceed the ratio in normal galaxies by a factor ≈4 at a confidence level ≳2 σ . A caveat is that the numbers for Seyferts are still small (six of Type Ib–Ic and six of Type II supernovae discovered as yet). Assumed to be real, this excess of Type Ib/c supernovae may indicate a burst of low-age star formation  ( τ ≲20 Myr)  , a high incidence of binary systems in the inner regions  ( r ≲0.4 R ₂₅)  of Seyfert galaxies, or a top-loaded mass function.     

The escape of ionizing photons from supernova-dominated primordial galaxies

In order to assess the contribution of Lyman break galaxies (LBGs) and Lyman α emitters (LAEs) at redshifts  3 < z < 7  to the ionization of intergalactic medium (IGM), we investigate the escape fractions of ionizing photons from supernova-dominated primordial galaxies by solving the three-dimensional (3D) radiative transfer. The model galaxy is employed from an ultra-high-resolution chemodynamic simulation of a primordial galaxy by Mori & Umemura, which well reproduces the observed properties of LAEs and LBGs. The total mass of model galaxy is  10¹¹ M_⊙  . We solve not only photoionization but also collisional ionization by shocks. In addition, according to the chemical enrichment, we incorporate the effect of dust extinction, taking the size distributions of dust into account. As a result, we find that dust extinction reduces the escape fractions by a factor of 1.5–8.5 in the LAE phase and by a factor of 2.5–11 in the LBG phase, while the collisional ionization by shocks increases the escape fractions by a factor of  ≈2  . The resultant escape fractions are 0.07–0.47 in the LAE phase and 0.06–0.17 in the LBG phase. These results are well concordant with the recent estimations derived from the flux density ratio at 1500 to 900 Å of LAEs and LBGs. Combining the resultant escape fractions with the luminosity functions of LAEs and LBGs, we find that high- z LAEs and LBGs can ionize the IGM at   z = 3–5  . However, ionizing radiation from LAEs as well as LBGs falls short of ionizing the IGM at   z > 6  . That implies that additional ionization sources may be required at   z > 6  .