Anja C. Andersen

Stellar sources of dust in the high-redshift Universe

With the aim of investigating whether stellar sources can account for the ≥10 8 Mdust masses inferred from mm/sub-mm observations of samples of 5 <z< 6.4 quasars, we develop a chemical evolution model which follows the evolution of metals and dust on the stellar characteristic lifetimes, taking into account dust destruction mechanisms. Using a grid of stellar dust yields as a function of the initial mass and metallicity over the range 1-40 M � and 0-1 Z� , we show that the role of asymptotic giant branch (AGB) stars in cosmic dust evolution at high redshift might have been overlooked. In particular, we find that (i) for a stellar population forming according to a present-day Larson initial mass function (IMF) with mch = 0.35 M� , the characteristic time-scale at which AGB stars dominate dust production ranges between 150 and 500 Myr, depending both on the assumed star formation history and on the initial stellar metallicity; (ii) this result is only moderately dependent on the adopted stellar lifetimes, but it is significantly affected by variations of the IMF: for a mch = 5M � , dust from AGB starts to dominate only on time-scales larger than 1 Gyr and SNe are found to dominate dust evolution when mch ≥10 M� . We apply the chemical evolution model with dust to the host galaxy of the most distant quasar at z = 6.4, SDSS J1148+5251. Given the current uncertainties on the star formation history of the host galaxy, we have considered two models: (i) the star formation history obtained in a numerical simulation by Li et al. which predicts that a large stellar bulge is already formed at z = 6.4, and (ii) a constant star formation rate of 1000 Myr −1 , as suggested by the observations if most of the far-infrared luminosity is due to young stars. The total mass of dust predicted at z = 6.4 by the first model is 2 × 10 8 M� , within the range of values inferred by observations, with a substantial contribution (∼80 per cent) of AGB dust. When a constant star formation rate is adopted, the contribution of AGB dust decreases to ∼50 per cent but the total mass of dust formed is a factor of 2 smaller. Both models predict a rapid enrichment of the interstellar medium with metals and a relatively mild evolution of the carbon abundance, in agreement with observational constraints. This supports the idea that stellar sources can account for the dust observed but show that the contribution of AGB stars to dust production cannot be neglected, even at the most extreme redshifts currently accessible to observations.

Dust Extinction in High-z Galaxies with Gamma-Ray Burst Afterglow Spectroscopy : The 2175 Å Feature at z = 2.45

We report the clear detection of the 2175A dust absorption feature in the optical afterglow spectrum of the gamma-ray burst (GRB) GRB070802 at a redshift of z=2.45. This is the highest redshift for a detected 2175A dust bump to date, and it is the first clear detection of the 2175A bump in a GRB host galaxy, while several tens of optical afterglow spectra without the bump have been recorded in the past decade. The derived extinction curve gives A_V=0.8-1.5 depending on the assumed intrinsic slope. Of the three local extinction laws, an LMC type extinction gives the best fit to the extinction curve of the host of GRB070802. Besides the 2175A bump we find that the spectrum of GRB070802 is characterized by unusually strong low-ionization metal lines and possibly a high metallicity for a GRB sightline ([Si/H]=-0.46+/-0.38, [Zn/H]=-0.50+/-0.68). In particular, the spectrum of GRB070802 is unique for a GRB spectrum in that it shows clear CI absorption features, leading us to propose a correlation between the presence of the bump and CI. The gas to dust ratio for the host galaxy is found to be significantly lower than that of other GRB hosts with N(HI)/A_V=(2.4+/-1.0)x10^21 cm^-2 mag^-1, which lies between typical MW and LMC values. Our results are in agreement with the tentative conclusion reached by Gordon et al. 2003 that the shape of the extinction curve, in particular the presence of the bump, is affected by the UV flux density in the environment of the dust.

A log NH I = 22.6 Damped Lyα Absorber in a Dark Gamma-Ray Burst: The Environment of GRB 050401*

The optical afterglow spectrum of GRB 050401 (at z = 2.8992 ± 0.0004) shows the presence of a damped Lyα absorber (DLA), with log N = 22.6 ± 0.3. This is the highest column density ever observed in a DLA and is about 5 times larger than the strongest DLA detected so far in any QSO spectrum. From the optical spectrum, we also find a very large Zn column density, implying an abundance of [Zn/H] = -1.0 ± 0.4. These large columns are supported by the early X-ray spectrum from Swift XRT, which shows a column density (in excess of Galactic) of log NH = 22.21 assuming solar abundances (at z = 2.9). The comparison of this X-ray column density, which is dominated by absorption due to α-chain elements, and the H I column density derived from the Lyα absorption line allows us to derive a metallicity for the absorbing matter of [α/H] = -0.4 ± 0.3. The optical spectrum is reddened and can be well reproduced with a power law with SMC extinction, where AV = 0.62 ± 0.06. But the total optical extinction can also be constrained independent of the shape of the extinction curve: from the optical to X-ray spectral energy distribution, we find 0.5 AV 4.5. However, even this upper limit, independent of the shape of the extinction curve, is still well below the dust column that is inferred from the X-ray column density, i.e., AV = 9.1. This discrepancy might be explained by a small dust content with high metallicity (low dust-to-metals ratio). Gray extinction cannot explain the discrepancy, since we are comparing the metallicity to a measurement of the total extinction (without reference to the reddening). Little dust with high metallicity may be produced by sublimation of dust grains or may naturally exist in systems younger than a few hundred megayears.

Production of dust by massive stars at high redshift

The large amounts of dust detected in sub-millimeter galaxies and quasars at high redshift pose a challenge to galaxy formation models and theories of cosmic dust formation. At z>6 only stars of relatively high mass (>3 M⊙) are sufficiently short-lived to be potential stellar sources of dust. This review is devoted to identifying and quantifying the most important stellar channels of rapid dust formation. We ascertain the dust production efficiency of stars in the mass range 3–40 M⊙ using both observed and theoretical dust yields of evolved massive stars and supernovae (SNe) and provide analytical expressions for the dust production efficiencies in various scenarios. We also address the strong sensitivity of the total dust productivity to the initial mass function. From simple considerations, we find that, in the early Universe, high-mass (>3 M⊙) asymptotic giant branch stars can only be dominant dust producers if SNe generate ≲3×10−3 M⊙ of dust whereas SNe prevail if they are more efficient. We address the challenges in inferring dust masses and star-formation rates from observations of high-redshift galaxies. We conclude that significant SN dust production at high redshift is likely required to reproduce current dust mass estimates, possibly coupled with rapid dust grain growth in the interstellar medium.

Simulations of dust-trapping vortices in protoplanetary discs

Local three-dimensional shearing box simulations of the compressible coupled dust-gas equations are used in the fluid approximation to study the evolution of different initial vortex configurations in a protoplanetary disc and their dust-trapping capabilities. The initial conditions for the gas are derived from an analytic solution to the compressible Euler equation and the continuity equation. The solution is valid if there is a vacuum outside the vortex. In the simulations the vortex is either embedded in a hot corona, or it is extended in a cylindrical fashion in the vertical direction. Both configurations are found to survive for at least one orbit and lead to accumulation of dust inside the vortex. This confirms earlier findings that dust accumulates in anticyclonic vortices, indicating that this is a viable mechanism for planetesimal formation.

Lyα RADIATIVE TRANSFER WITH DUST: ESCAPE FRACTIONS FROM SIMULATED HIGH-REDSHIFT GALAXIES

The Lyα emission line is an essential diagnostic tool for probing galaxy formation and evolution. Not only is it commonly the strongest observable line from high-redshift galaxies, but from its shape detailed information about its host galaxy can be revealed. However, due to the scattering nature of Lyα photons increasing their path length in a nontrivial way, if dust is present in the galaxy, the line may be severely suppressed and its shape altered. In order to interpret observations correctly, it is thus of crucial significance to know how much of the emitted light actually escapes the galaxy. In the present work, using a combination of high-resolution cosmological hydrosimulations and an adaptively refinable Monte Carlo Lyα radiative transfer code including an environment dependent model of dust, the escape fractions f esc of Lyα radiation from high-redshift (z = 3.6) galaxies are calculated. In addition to the average escape fraction, the variation of f esc in different directions and from different parts of the galaxies is investigated, as well as the effect on the emergent spectrum. Escape fractions from a sample of simulated galaxies of representative physical properties are found to decrease for increasing galaxy virial mass M vir, from f esc approaching unity for M vir ~ 109 M ☉ to f esc less than 10% for M vir ~ 1012 M ☉. In spite of dust being almost gray, it is found that the emergent spectrum is affected nonuniformly, with the escape fraction of photons close to the line center being much higher than of those in the wings, thus effectively narrowing the Lyα line.

Genesis and evolution of dust in galaxies in the early Universe II. Rapid dust evolution in quasars at z 6

Aims. We intend to assess the most plausible scenarios for generating large amounts of dust in high-z quasars (QSOs) on the basis of observationally derived physical properties of QSOs at z 6. Methods. We use a chemical evolution model to compute the temporal progression of quantities such as the amount of dust and gas, stellar masses, star formation rates (SFRs) and the metallicity for various combinations of the initial mass function (IMF), the mass of the galaxy, dust production efficiencies, and the degree of dust destruction in the ISM. We investigate the influence of the SFR on the evolution of these quantities, and determine the earliest epochs at which agreement with observations can be achieved. We apply the obtained results to individual QSOs at z 6. Results. We find that large quantities of dust can be generated rapidly as early as 30 Myr after the onset of the starburst when the SFR .

Winds of M- and S-type AGB stars: An unorthodox suggestion for the driving mechanism

Context: Current knowledge suggests that the dust-driven wind scenario provides a realistic framework for understanding mass loss from C-rich AGB stars. For M-type objects, however, recent detailed models demonstrate that radiation pressure on silicate grains is not sufficient to drive the observed winds, contrary to previous expectations. Aims: In this paper, we suggest an alternative mechanism for the mass loss of M-type AGB stars, involving the formation of both carbon and silicate grains due to non-equilibrium effects, and we study the viability of this scenario. Methods: We model the dynamical atmospheres and winds of AGB stars by solving the coupled system of frequency-dependent radiation hydrodynamics and time-dependent dust formation, using a parameterized description of non-equilibrium effects in the gas phase. This approach allows us to assess under which circumstances it is possible to drive winds with small amounts of carbon dust and to get silicate grains forming in these outflows at the same time. Results: The properties of the resulting wind models, such as mass-loss rates and outflow velocities, are well within the observed limits for M-type AGB stars. Furthermore, according to our results, it is quite unlikely that significant amounts of silicate grains will condense in a wind driven by a force totally unrelated to dust formation, as the conditions in the upper atmosphere and wind acceleration region put strong constraints on grain growth. Conclusions: .The proposed scenario provides a natural explanation for the observed similarities in wind properties of M-type and C-type AGB stars and implies a smooth transition for stars with increasing carbon abundance, from solar-composition to C-rich AGB stars, possibly solving the longstanding problem of the driving mechanism for stars with a C/O close to one.

Dust and the Type II-Plateau Supernova 2004et

We present mid-infrared (MIR) observations of the Type II-plateau supernova (SN) 2004et, obtained with the {\it Spitzer Space Telescope} between days 64 and 1406 past explosion. Late-time optical spectra are also presented. For the period 300-795 days past explosion, we argue that the spectral energy distribution of SN 2004et comprises (a) a hot component due to emission from optically thick gas, as well as free-bound radiation, (b) a warm component due to newly formed, radioactively heated dust in the ejecta, and (c) a cold component due to an IR echo from the interstellar-medium dust of the host galaxy, NGC 6946. There may also have been a small contribution to the IR SED due to free-free emission from ionised gas in the ejecta. We reveal the first-ever spectroscopic evidence for silicate dust formed in the ejecta of a supernova. This is supported by our detection of a large, but progressively declining, mass of SiO. However, we conclude that the mass of directly detected ejecta dust grew to no more than a few times 10^(-4)Msun. We also provide evidence that the ejecta dust formed in comoving clumps of fixed size. We argue that, after about two years past explosion, the appearance of wide, box-shaped optical line profiles was due to the impact of the ejecta on the progenitor circumstellar medium and that the subsequent formation of a cool, dense shell was responsible for a later rise in the MIR flux. This study demonstrates the rich, multi-faceted ways in which a typical core-collapse supernova and its progenitor can produce and/or interact with dust grains. The work presented here adds to the growing number of studies which do not support the contention that SNe are responsible for the large mass of observed dust in high-redshift galaxies.

Genesis and evolution of dust in galaxies in the early Universe - I. Modelling dust evolution in starburst galaxies

Aims. The aim is to elucidate the astrophysical conditions required for generating large amounts of dust in massive starburst galaxies at high redshift. Methods. We have developed a numerical galactic chemical evolution model. The model is constructed such that the effect of a wide range of parameters can be investigated. It takes into account results from stellar evolution models, a differentiation between diverse types of core collapse supemovae (CCSN), and the contribution of asymptotic giant branch (AGB) stars in the mass range 3-8 M ⊙ . We consider the lifetime-dependent yield injection into the interstellar medium (ISM) by all sources, and dust destruction due to supernova (SN) shocks in the ISM. We ascertain the temporal progression of the dust mass and the dust-to-gas and dust-to-metal mass ratios, as well as other physical properties of a galaxy, and study their dependence on the mass of the galaxy, the initial mass function (IMF), dust production efficiencies, and dust destruction in the ISM. Results. The amount of dust and the physical properties of a galaxy strongly depend on the initial gas mass available. Overall, while the total amount of dust produced increases with galaxy mass, the detailed outcome depends on the SN dust production efficiency, the IMF, and the strength of dust destruction in the ISM. Dust masses are higher for IMFs biased towards higher stellar masses, even though these IMFs are more strongly affected by dust destruction in the ISM. The sensitivity to the IMF increases as the mass of the galaxy decreases. SNe are primarily responsible for a significant enrichment with dust at early epochs ( 10 8 M ⊙ . Our preferred scenario is dominated by SN dust production in combination with top-heavy IMFs and moderate dust destruction in the ISM.