ABSTRACT
Analyses of high-redshift ultraluminous infrared (IR) galaxies traditionally use the observed optical to submillimeter spectral energy distribution (SED) and estimates of the dynamical mass as observational constraints to derive the star formation rate (SFR), the stellar mass, and age of these objects. An important observational constraint neglected in the analysis is the mass of dust giving rise to the IR emission. In this paper we add this constraint to the analysis of AzTEC–3. Adopting an upper limit to the mass of stars and a bolometric luminosity for this object, we construct different stellar and chemical evolutionary scenarios, constrained to produce the inferred dust mass and observed luminosity before the associated stellar mass exceeds the observational limit. We use the PÉGASE population synthesis code and a chemical evolution model to follow the evolution of the galaxy's SED and its stellar and dust masses as a function of galactic age for seven different stellar initial mass functions (IMFs). We find that the model with a Top Heavy IMF provided the most plausible scenario consistent with the observational constraints. In this scenario the dust formed over a period of ∼200 Myr, with an SFR of ∼500 M☉ yr−1. These values for the age and SFR in AzTEC–3 are significantly higher and lower, respectively, from those derived without the dust mass constraint. However, this scenario is not unique, and others cannot be completely ruled out because of the prevailing uncertainties in the age of the galaxy, its bolometric luminosity, and its stellar and dust masses. A robust result of our models is that all scenarios require most of the radiating dust mass to have been accreted in molecular clouds. Our new procedure highlights the importance of a multiwavelength approach, and of the use of dust evolution models in constraining the age and the star formation activity and history in galaxies.
Export citation and abstract BibTeX RIS
1. INTRODUCTION
The origin of the dust in ultraluminous high-redshift galaxies with redshifts z ≳ 4 is still unresolved, with supernovae (SNe), evolved asymptotic giant branch (AGB) stars, the winds around active galactic nuclei (AGNs), or growth in molecular clouds as possible sources (Elvis et al. 2002; Morgan & Edmunds 2003; Maiolino et al. 2004; Dwek et al. 2007; Valiante et al. 2009; Dwek & Cherchneff 2011; Michałowski et al. 2010b; Gall et al. 2011a, 2011b). Further complicating the origin of dust in some objects is the energy source powering the IR emission. An IR luminosity powered by luminous massive stars will suggest a high present star formation rate (SFR), favoring an SN origin for the dust. An IR luminosity that is powered by an AGN suggests a low present SFR and may favor the delayed injection of dust by AGB stars that may have formed during an earlier phase of much higher star formation activity (Valiante et al. 2009; Dwek & Cherchneff 2011).
No such ambiguity exists for AzTEC–3, an ultraluminous IR galaxy located at redshift z = 5.3, that was first detected at submillimeter wavelengths by Scott et al. (2008). It is a radio-dim source (Younger et al. 2007) with no X-ray counterpart (Capak et al. 2011), and part of a young cluster, representing an overdensity of early massive galaxies at the same redshift (Capak et al. 2011). Optical to submillimeter observations summarized in Capak et al. (2011) show that the energy output is dominated by its far-infrared (IR) emission with luminosities between ∼(0.6–2) × 1013 L☉. The lack of any radio or X-ray counterparts suggests that this luminosity is powered by a 1800 M☉ yr−1 starburst (Riechers et al. 2010; Capak et al. 2011). The total CO-inferred gas mass is ∼5 × 1010 M☉, which will be depleted in less than ∼30 Myr at this SFR. Such energy output suggests the presence of a large amount of dust, which had to be created (at least) in less than ∼1.1 Gyr, the age of the universe at that redshift. AzTEC–3 therefore falls in the category of extreme starburst galaxies with large dust masses (≳ 108 M☉) described by Michałowski et al. (2010a) and Robson et al. (2004).
In this paper, we present a new comprehensive approach for studying the star formation properties and the origin of dust in AzTEC–3. The approach is quite general and can be applied to any local or high-redshift galaxy. We first use the UV to far-IR spectral energy distribution (SED) of the galaxy to determine its bolometric luminosity and dust mass. Its SFR depends on the adopted initial mass function (IMF), and we generate a series of models relating the bolometric luminosity to the stellar mass and galactic age for different stellar IMFs (Section 2). In Section 3 we use the observed UV–optical (UVO) spectrum, the UVNIR SED, and dynamical limits on the stellar mass to constrain the age of the galaxy from observations. We then generate different scenarios for the evolution of the dust, characterized by different IMFs, and examine which scenario is capable of producing the inferred dust mass in AzTEC–3 within the allowed age limit, that is, before the concurrent stellar mass production exceeds the stellar mass limit (Section 4). Our results are discussed and summarized in Section 5.
In all our calculations we adopt a flat ΛCDM cosmology, with a baryonic density parameter Ωb = 0.044, a total matter (dark+baryonic) density parameter of Ωm = 0.27, a vacuum energy density ΩΛ = 0.73, and a Hubble constant of H0 = 70 km s−1 Mpc−1 (Spergel et al. 2007).
2. PROPERTIES OF AzTEC–3
2.1. Bolometric Luminosity and Dust Mass
Figure 1 depicts the visible to far-IR flux, ν Fν, from the galaxy calculated for a redshift of z = 5.3. The observed fluxes were taken from Capak et al. (2011).
The figure shows that the energy output from the galaxy is dominated by the stellar radiation that is absorbed and reradiated by dust at IR wavelengths. Observed IR flux densities and their 1σ uncertainties are listed in Table 1. Figure 2 shows the far-IR galaxy's SED in its rest frame (left panel) and the observer's frame (right panel). Also shown in the figure are fits of the spectra of individual dust species to the observations. In order to determine the uncertainties in the derived dust temperatures, dust masses, and IR luminosities, we constructed 1000 realizations of the observed spectrum. The flux at each wavelength was randomly chosen from a normal distribution centered around the nominal flux. The derived quantities for the different dust compositions are listed in Table 2. The Fe and graphite dust spectra are similar to that of the silicate dust and were left out of the figure for the sake of clarity. Dust masses vary greatly from 0.3 × 109 M☉, if all the dust in AzTEC–3 is in the form of amorphous carbon, to 3.2 × 109 M☉, if it consists of only silicate dust. Because of the skin effect, the absorptivity of metallic Fe depends on the grain radius, and the mass presented in the table was calculated for a radius of 0.33 μm for which the mass absorption coefficient, κ, was maximum at the wavelength of 174.6 μm. This choice of κ minimizes the mass of Fe dust required to fit the observed fluxes. Since AzTEC–3 is not likely to contain only pure Fe, silicate, or carbon-type dust, any dust consisting of a mixture of these elements will vary between these mass estimates.
Download figure:
Standard image High-resolution imageTable 1. Summary of Submillimeter Observations of AzTEC–3a
λobs | Fν(λ) | σ | λ0 | Lν(λ0) |
---|---|---|---|---|
(μm) | (mJy) | (mJy) | (μm)a | (L☉ Hz−1) |
890 | 8.7 | 1.5 | 141.3 | 1.158 |
1100 | 7.6 | 1.2 | 174.6 | 1.011 |
2000 | 3.7 | 1.4 | 317.5 | 0.492 |
2730 | 0.47b | ... | 433.3 | 0.0412 |
3280 | 0.20b | ... | 520.6 | 0.0173 |
8165 | 0.087b | ... | 1296.0 | 0.0116 |
Notes. aRest-frame wavelengths and specific luminosities were calculated for a redshift of z = 5.3. b3σ upper limits.
Download table as: ASCIITypeset image
Table 2. Derived Temperatures, Masses, and IR Luminosities for Different Dust Compositions and Galaxy Templates
Composition | Td | Md | LIR | κ |
---|---|---|---|---|
(K) | (109M☉) | (1013L☉) | (cm2 g−1)a | |
Feb | 32.4 ± 6.1 | 1.4+1.2− 0.7 | 0.6+0.4− 0.2 | 16.5 |
Graphite | 27.5 ± 4.5 | 2.2+4.6− 2.2 | 0.5+0.3− 0.2 | 17.6 |
Silicate | 28.5 ± 4.3 | 3.2+2.6− 1.4 | 0.6+0.4− 0.2 | 10.5 |
A-carbon | 46.9 ± 14.8 | 0.3+0.35− 0.17 | 1.1+1.7− 0.67 | 28.0 |
M 82 template | ... | 2.6 ± 0.9 | 1.7 ± 0.6 | ... |
Arp 220 template | ... | 1.4 ± 0.8 | 0.8 ± 0.5 | ... |
NGC 6946 template | ... | 18 ± 8 | 0.6 ± 0.3 | ... |
Notes. aThe mass absorption coefficient at wavelength λ = 174.6 μm. bFe mass was calculated for a grain radius of 0.33 μm. See the text for details.
Download table as: ASCIITypeset image
To get a more realistic presentation of the mixture of dust compositions and size distribution in AzTEC–3 we also fitted its spectrum with that of two local starburst galaxies: M82, and Arp220, and a normal star forming galaxy NGC6946. The spectra of these galaxies were obtained by using population synthesis models to derive their stellar emission, chemical evolution models to derive their metallicity and dust composition, and radiative transfer models including the stochastic heating of dust to derive their IR spectrum and dust mass (Galliano et al. 2008). Transferred to a redshift of z = 5.3, the IR luminosities and dust masses for all three galaxies were derived from their normalization factors, obtained from the least-squares fit of their spectra to the AzTEC–3 observations. The spectra depicted in the figure represent the total IR emission from dust residing in H ii and photodissociation regions. Both starburst galaxies, M82 and Arp220, provided a very good fit to the data. Dust masses varied from 1.4 to 2.6 × 109 M☉and IR luminosities (∼3–1000 μm) from 0.79 to 1.7 × 1013 L☉. In contrast, the dust in NGC6946 is too cold, and its far-IR emission exceeds the 3σ upper limits on the underlying continuum of the detected CO(J = 5 → 4) and CO(J = 6 → 5) lines (Riechers et al. 2010). This is evident in the inferred dust mass for this galaxy, which is in excess of ∼1010 M☉.
For the sake of being definitive we adopt throughout this paper a dust mass, Mdust, of (2.0 ± 1.0) × 109 M☉, representing the average value between M82 and Arp220. Our choice of average dust mass derived from these template galaxies is more realistic than any choice based on the assumption that all the IR emission is produced by a single dust specie. The average IR luminosity of these two galaxies is 1.0 × 1013 L☉ and since the bolometric luminosity of AzTEC–3 is dominated by the IR emission component we adopt a bolometric luminosity, Lbol, that is equal to 1.0 × 1013 L☉.
2.2. The Star Formation Rate
The far-IR luminosity provides a good measure of the star formation activities in galaxies (Kennicutt 1998). In a galaxy with ongoing star formation the bolometric luminosity and stellar mass increase with time because of the cumulative contribution from long-lived low-mass stars. Derived SFRs depend therefore on the assumed star formation history of the galaxy. Furthermore, the relation between the luminosity and SFR depends on the adopted stellar IMF. Massive stars radiate much more efficiently than low mass ones (L∝M3.3, in solar units) so that an IMF weighted more heavily toward massive stars will require a lower SFR for a given observed luminosity. To explore the dependence of the SFR on the stellar IMF, we used seven different functional forms to characterize the latter. The parameters characterizing the IMFs are given in Table 3, and Figure 3 depicts the functional form of four select IMFs. Table 3 also includes two characteristic values of the different IMFs: <m >, which is the IMF-averaged mass, and m⋆, which is the mass of stars that need to be formed to produce one Type II SN event. With these definitions, the stellar birthrate (yr−1) is given by SFR/<m >and the SN rate by SFR/m⋆.
Download figure:
Standard image High-resolution imageTable 3. Characteristics of the Stellar IMFsa
IMF | α1 | α2 | α3 | M1 | M2 | M3 | M4 | <m >b | m⋆c |
---|---|---|---|---|---|---|---|---|---|
Glazebrook | 1.5 | 2.15 | ... | 0.1 | 0.5 | 100 | ... | 0.72 | 70.0 |
Kroupa | 1.3 | 2.3 | ... | 0.1 | 0.5 | 100 | ... | 0.64 | 89.6 |
Mass Heavy | 2.35 | ... | ... | 1.0 | 100 | ... | ... | 3.1 | 52.7 |
Paunchy | 1.0 | 1.7 | 2.6 | 0.1 | 0.5 | 4 | 100 | 1.1 | 66.4 |
Salpeter | 2.35 | ... | ... | 0.1 | 100 | ... | ... | 0.35 | 134.7 |
Starburst | 2.35 | ... | ... | 5.0 | 100 | ... | ... | 12.8 | 24.4 |
Top Heavy | 1.95 | ... | ... | 0.1 | 100 | ... | ... | 0.78 | 55.4 |
Notes. aAll masses are in M☉. Stellar IMFs are given by for M1 < M < M2, for M2 < M < M3, and for M3 < M < M4. bThe IMF-averaged stellar mass. cThe mass of all stars born per SN event (Dwek & Cherchneff 2011).
Download table as: ASCIITypeset image
To determine the SFR and the mass of stellar objects as a function of galactic age we ran the PÉGASE population synthesis code developed by Fioc & Rocca-Volmerange (1997) for the seven different stellar IMFs, listed in Table 3, for a constant SFR of 1 M☉ yr−1. Table 4 lists the corresponding bolometric luminosities and stellar masses at different galactic ages. For example, for the adopted Lbol = 1 × 1013 L☉, the table shows that for a Glazebrook IMF the SFR is ∼735 M☉ yr−1 and the total mass of stars is ∼5.6 × 1010 M☉ when the galaxy reaches an age of 100 Myr. At the same age, a starburst IMF will require an SFR of only ∼270 M☉ yr−1 to produce the same luminosity and will have created ∼7.5 × 109 M☉ of stars. Figure 4 depicts these results in graphic form, showing the constant SFR required to produce a bolometric luminosity of 1 × 1013 L☉ at a given age for four different IMFs. The bold section in each line depicts the allowed range of galaxy ages before the model violates the constraint on the stellar mass, a topic discussed in detail in Section 3.3 below. We emphasize that this figure does not represent the evolution of the SFR with time. For each stellar IMF, it depicts the value of the constant SFR needed to produce a bolometric luminosity of 1 × 1013 L☉ at the assumed galactic age. The required SFR is lower for older systems because a larger fraction of Lbol arises from the accumulation of lower mass stars over the galaxy's lifetime.
Download figure:
Standard image High-resolution imageTable 4. Starburst Luminosities and Stellar Massesa
Starburst Age (yr) | 107 | 108 | 109 | 1010 |
---|---|---|---|---|
Glazebrook | ||||
Lbol (L☉) | 9.26e+09 | 1.36e+10 | 1.74e+10 | 2.20e+10 |
Mstar (M☉) | 1.03e+07 | 7.61e+07 | 5.89e+08 | 4.50e+09 |
Kroupa | ||||
Lbol (L☉) | 7.33e+09 | 1.13e+10 | 1.54e+10 | 2.08e+10 |
Mstar (M☉) | 1.05e+07 | 8.06e+07 | 6.40e+08 | 4.92e+09 |
Mass Heavy | ||||
Lbol (L☉) | 1.20e+10 | 1.90e+10 | 2.65e+10 | 3.25e+10 |
Mstar (M☉) | 1.02e+07 | 6.70e+07 | 3.80e+08 | 1.15e+09 |
Paunchy | ||||
Lbol (L☉) | 7.83e+09 | 1.43e+10 | 2.05e+10 | 2.53e+10 |
Mstar (M☉) | 1.05e+07 | 7.64e+07 | 4.92e+08 | 3.08e+09 |
Salpeter | ||||
Lbol (L☉) | 4.71e+09 | 7.49e+09 | 1.06e+10 | 1.49e+10 |
Mstar (M☉) | 1.07e+07 | 8.76e+07 | 7.56e+08 | 6.46e+09 |
Starburst | ||||
Lbol (L☉) | 2.55e+10 | 3.67e+10 | 3.62e+10 | 3.41e+10 |
Mstar (M☉) | 9.20e+06 | 2.76e+07 | 3.17e+07 | 7.18e+07 |
Top Heavy | ||||
Lbol (L☉) | 1.56e+10 | 2.04e+10 | 2.34e+10 | 2.56e+10 |
Mstar (M☉) | 9.72e+06 | 6.23e+07 | 4.45e+08 | 3.32e+09 |
Note. aEntries are calculated for a constant SFR of 1M☉yr−1.
Download table as: ASCIITypeset image
3. THE AGE OF AzTEC–3
3.1. Spectral Constraints on the Age
The preceding section illustrated the importance of knowing the galactic age for determining the star formation, and the stellar and elemental enrichment histories of a galaxy. The UVO spectrum can contain important information for determining this age. Figure 1 in Capak et al. (2011) presents the full resolution spectrum of AzTEC–3, obtained with the Keck II telescope. The spectrum shows Si ii and C iv stellar absorption features which are indicative of the presence of a population of young O and B stars. In a single burst of star formation the spectral contribution from young stars should diminish with time. The presence of these features therefore led them to the conclusion that AzTEC–3 should be less than ∼30 Myr old. However, in an ongoing star formation scenario, the stellar population will always contain young O and B stars, and their spectral features will only fade at much later times as they are gradually overwhelmed by the featureless continuum from the less massive stars.
To examine if the observed spectrum can be used to determine the age of the galaxy we constructed synthetic UV spectra for constant star formation scenarios and different stellar IMFs (TH = Top Heavy; Sp = Salpeter) using Starburst99 (Leitherer et al. 1999; Vázquez & Leitherer 2005; Leitherer et al. 2010). Figure 5 compares the synthetic spectra with the observed one. The figure shows that the TH model spectra have not evolved significantly in shape between the ages of 30 and 200 Myr. Their total intensity has, of course evolved, an effect not shown in the figure because of the offset. They are also identical to the Sp model at 30 Myr. The simulations show that when star formation is an ongoing process, the spectrum does not evolve significantly because even a very small number of O/B stars will dominate the UV emission. The AzTEC–3 spectrum shows some of the same stellar absorption features present in the synthetic spectra, but contains in addition features caused by interstellar absorption, a process that is not included in the Starburst99 models. The right panel of the figure shows the correlation between the observed spectrum and that of one of the models. Except for a constant offset caused by the evolution of the total intensity of the spectrum, the correlation plot is identical for all models.
Download figure:
Standard image High-resolution imageThe similarity of the AzTEC–3 spectrum to the model spectra, and the fact that the model spectra are essentially identical even though they represent different epochs and were generated with different IMFs, shows that the observed ∼1200–1600 Å spectrum cannot be used to constrain the age of the galaxy when star formation is an ongoing process.
3.2. SED Constraints on the Age: The Age Degeneracy
The UVNIR SED of a galaxy can provide independent constraints on its age. Using the results of the previous section we created a library of intrinsic stellar SEDs for the seven IMFs and for a grid of galaxy ages, ranging up to 800 Myr. The results are shown in Figures 6 and 7 for four select IMFs and for ages between 10 and 400 Myr. Each panel in the figure presents the intrinsic SED (black line) generated by a distinct evolutionary sequence that is terminated at the assumed galactic age presented in the panel. Each evolutionary sequence is characterized by a constant SFR which is distinct from other sequences and determined by the stellar IMF and the galactic age for the adopted bolometric luminosity of 1 × 1013 L☉ (see Table 4 and Figure 4). So the Salpeter 20 Myr panel represents the SED of a galaxy that underwent a constant SFR of ∼1900 M☉ yr−1 for that length of time. In each panel, the intrinsic SED is clearly in excess of the observed stellar flux (red diamonds), which has been significantly attenuated by dust. Lacking a definitive dust model for the IR emission, we used the Calzetti law (Calzetti et al. 2000) to characterize the extinction. For each evolutionary sequence we calculated the visual optical depth, τ(V), that provided the best least-squares fit to the observations. Values of τ(V) ranged from ∼2 to 4. The attenuated stellar spectra are presented by violet lines in the figures.
Download figure:
Standard image High-resolution imageDownload figure:
Standard image High-resolution imageFigure 8 depicts the minimum χ2 of the fit of the attenuated spectra to the data as a function of galaxy age and stellar IMF. The figure shows that for each IMF there exists an age degeneracy for which the value of χ2 is essentially unchanged. This degeneracy is caused by the fact that the galaxy's SED changes very little with time for the first ∼100 Myr. The small changes that do occur are compensated for with small changes in the magnitude of the attenuation. This effect is exhibited in both the figures of the SEDs and that depicting the value χ2 with galactic age. For example, Figure 8 shows that the Salpeter and Glazebrook IMFs provide equally good fits to the observed UVNIR SED for ages up to ∼100 Myr, and the TH IMF for ages up to ∼200 Myr. The fits become significantly worse thereafter, an effect that is also clearly seen in Figures 6 and 7. The starburst IMF provides an equally good fit to the observed spectrum at all epochs. Better fits to the observed spectra at later epochs can probably be obtained if we adopted a different extinction law, or a time-dependent extinction law that reflects the changing dust composition with time. Such changes would need to be consistent with the far-IR emission as well. Such detailed investigation into the extinction law (i.e., optical properties of the dust and the physical distribution of the dust with respect to the radiation sources) is currently beyond the scope of this paper.
Download figure:
Standard image High-resolution image3.3. Stellar Mass Constraints of the Age: A Partial Lifting of the Degeneracy
The ambiguity in stellar age can be partially lifted by considering the upper limits on the stellar mass in the galaxy. Capak et al. (2011) derived a stellar mass of ∼1010 M☉, using the Maraston library of stellar SEDs (Maraston 2005) and a galactic age of ∼30 Myr. This age was inferred from the Keck II spectrum, assuming that all the star formation occurred in a single burst. Since we have shown that the UV spectrum cannot constrain the age of AzTEC–3, we prefer to use a weaker, but observationally determined, constraint on the stellar mass. Using kinematic CO data and limits on the angular size of the galaxy, Riechers et al. (2010) derived the values of Mgas, the mass of the gas (dominated by the molecular component), and Mdyn, the dynamical mass of the galaxy. The dynamical mass is given by Mdyn ≃ Mgas+ Mstars + MDM, where Mstars and MDM are, respectively, the mass in stars and dark matter. The dynamical mass offers therefore a strict upper limit on the stellar mass. Adopting a gas fraction (≡ Mgas/Mdyn) of 0.5 (Riechers et al. 2010), we adopt an upper limit of Mstars ≃ Mgas ≃ 5 × 1010 M☉ on the stellar mass of the galaxy.
Since the mass of stars builds up at different rates for the different IMFs, this upper limit can provide useful constraints on the age of the galaxy. This effect is shown in Figure 9, which depicts the derived mass of stars as a function of the assumed age of AzTEC–3 for different stellar IMFs. Each symbol represents a distinct evolutionary sequence in which the SFR proceeded at a constant rate determined by the assumed galactic age and stellar IMF (see Figure 4). The Salpeter IMF requires the largest SFR and is therefore the first to exceed this mass limit after tage(max) ∼ 30 Myr. The Glazebrook IMF exceeds this mass limit after ∼80 Myr and the TH IMF after ∼200 Myr. The starburst IMF is the most efficient one in producing the observed bolometric luminosity with a low SFR and does not have a large population of low-mass stars. Consequently, it does not exceed the adopted stellar mass limit within ∼1 Gyr, the approximate age of the universe at that redshift.
Download figure:
Standard image High-resolution image4. DUST EVOLUTION: A NEW CONSTRAINT ON GALACTIC PROPERTIES
The presence of (2 ± 1) × 109 M☉ of dust provides new additional constraints on the SFR and history of AzTEC–3. To demonstrate this effect we used the chemical evolution model described in Dwek & Cherchneff (2011) to follow the evolution of the dust for the four different astrophysical scenarios determined by the stellar IMFs and the upper limit of the mass of stars in the galaxy. For each IMF we adopted a galactic age, tage(max), taken to be equal to the last epoch before the mass of stars produced by the model violates the stellar mass limit. The ages are ∼30, 80, and 200 Myr for the Salpeter, Glazebrook, and TH IMFs, respectively (see Figure 9). For the starburst IMF we chose, quite arbitrary, a value of tage(max) = 800 Myr. The models then assume that star formation proceeds at a constant rate of ∼1800, 700, 500, and 300 M☉ yr−1, respectively (see Figure 4).
We used a closed box model to follow the chemical evolution of AzTEC–3. The initial gas mass was chosen so that each model reproduced the inferred gas mass of ∼5 × 1010 M☉ at age t = tage(max). The chemical evolution model follows the changes in metallicities over timescales that are comparable to the main-sequence lifetimes of the massive SN progenitors. Consequently, in addition to the delayed recycling of matter by long-lived AGB stars, the model also includes the delayed recycling of matter by massive stars. AGB yields for stars with masses below 8 M☉ were taken from Karakas & Lattanzio (2007) for a metallicity of z = 0.008. SN dust yields of massive stars were taken as an average between the zero and solar metallicity yields given by Heger & Woosley (2010) and Woosley & Heger (2007), respectively (see Table 5 for select stellar masses). Dust condensation efficiencies were chosen to be unity for AGB stars, and 0.5 for Type II SNe.
Table 5. Yields of Condensable Elements and Maximum Dust Masses in Massive Starsa
Element/ | 20 M☉ | 50 M☉ | ||
---|---|---|---|---|
Metallicity | z = 0 | z = Z☉ | z = 0 | z = Z☉ |
12C | 0.21 | 0.23 | 1.86 | 1.61 |
16O | 1.34 | 1.38 | 11.2 | 9.7 |
24Mg | 0.065 | 0.065 | 0.31 | 0.33 |
28Si | 0.0025 | 0.19 | 0.27 | 0.21 |
56Fe | ∼0 | 0.12 | 0.41 | 0.28 |
Carbon dust | 0.21 | 0.23 | 1.86 | 1.61 |
Silicate dustb | 0.11 | 0.52 | 1.10 | 1.00 |
Iron dust | 0.0 | 0.12 | 0.41 | 0.28 |
Total dustc | 0.32 | 0.87 | 3.4 | 2.9 |
Notes. aZero metallicity yields were calculated for models with no mixing and with explosion energies of 1.2 × 1051 erg and 5 × 1051 erg, for the 20 and 50 M☉ stars, respectively (Heger & Woosley 2010); Z☉ metallicity yields were taken from Woosley & Heger (2007). All entries are in units of M☉. bThe silicate dust mass was calculated assuming that all the Mg and Si are locked up in MgO and SiO2 dust precursors. cTotal dust masses listed assume a 100% condensation efficiency.
Download table as: ASCIITypeset image
The results of our model calculations are shown in Figure 10. The green curve depicts the evolution of the gas mass which plummets when the constant SFR depletes all the gas (and dust) into stars. The blue curves depict the evolution of the total dust mass, which is the sum of the contribution from SNe and AGB stars. The AGB contribution is a small fraction of the total dust mass and is shown separately by red curves. Models were run for two cases: the first ignores the effect of grain destruction and is depicted by the solid blue curves; and the second takes the effect of grain destruction into account and is depicted by the dashed blue curves. The grain destruction rate was calculated by integrating the mass of dust destroyed as a function of shock velocity (Jones 2004), over the evolution of a supernova remnant. The mass of dust destroyed by a single SN remnant during its lifetime is about mdest ≈ 3 M☉ (Dwek et al. 2007) and the grain destruction rate is given by SFR × mdest/m⋆. The bold section of each curve depicts the evolution over the time span allowed before the galaxy violates the upper limit on the stellar mass (see Figure 9). The shaded region depicts the range of dust masses inferred from the far-IR observations.
Download figure:
Standard image High-resolution imageThe figure shows that when grain destruction is taken into account, none of the models are capable of producing the required dust mass at any time. With no grain destruction, all models depicted in the figure are capable of producing the required dust mass, albeit at different quantities. These idealized models can be considered as representatives of models in which rapid accretion in molecular clouds compensates for the mass of dust that is destroyed in the interstellar medium (ISM). So the models show that rapid grain growth in clouds is an essential process for producing the inferred dust mass in AzTEC–3.
The model that is clearly most capable of producing the inferred dust mass is the one with the TH IMF. At the age of 200 Myr it produces more than the inferred upper range of the dust mass, thus allowing for a modest amount of grain destruction. It therefore does not require that the accretion in molecular clouds reconstitute all the dust that was destroyed in the more diffuse ISM. The Starburst model is just capable of producing the lower range of inferred dust mass at t = 400 Myr and requires that all the refractory elements returned by grain processing into the gas are quickly reaccreted onto surviving grains in molecular clouds. The models with a Salpeter or Glazebrook IMF fall in between these two models.
Figure 11 depicts the spectrum produced by stars generated by the most plausible scenario that satisfies both the stellar mass and the dust mass constraints of AzTEC–3: a galaxy with a TH IMF undergoing a constant SFR at a rate of 300 M☉ yr−1 over a time span of 200 Myr, and allowing for over half of its dust to be destroyed. The stellar SED is attenuated by a Calzetti extinction law with τ(V) = 2.3. The intrinsic stellar spectrum is shown as the gray curve, and has a luminosity of 1 × 1013 L☉, which lies between the luminosity of the two galactic templates, M82 and Arp 220 (see Table 2). This scenario is not unique, but the most plausible one among the models. A TH IMF was also the preferred scenario for the rapid and efficient production of dust in the early universe by Gall et al. (2011a) and Dunne et al. (2010).
Download figure:
Standard image High-resolution image5. A POPULATION III ORIGIN FOR THE DUST?
Population III (Pop III) stars, the first population of stars being born out of a pristine, zero metallicity gas, may also be important contributors to the dust in these young galaxies (Nozawa et al. 2003, 2009; Todini & Ferrara 2001; Cherchneff & Dwek 2010). They are believed to be massive (Bromm & Larson 2004) and supported by radiation pressure at the Eddington luminosity:
Consider a population of massive Pop III stars all of identical mass m, and let Lbol be the total bolometric output of this stellar population. To maintain a steady abundance of such stars, they have to be created at a rate given by
where τMS, their main-sequence lifetime, is determined by the time it takes to radiate all the energy released in the production of helium and heavy elements at the Eddington luminosity. It is independent of stellar mass and given by
The SFR needed to maintain the observed luminosity is independent of the stellar mass and given by SFR = m (dN/dt) ≈ 10−12 × Lbol M☉ yr−1. The rate of dust production is given by
where Ydust is the amount of dust produced in these objects.
Dust formation in massive Pop III stars has only been calculated for 170 M☉ stars (Nozawa et al. 2003; Cherchneff & Dwek 2010). Such a Pop III star will have a luminosity of ∼5 × 106 L☉. The dust yield depends on the whether the ejecta maintains its stratified compositional structure or becomes completely mixed and ranges from ∼7 to ∼35 M☉, for unmixed and fully mixed ejecta, respectively. The total mass of dust produced at time t, assuming no grain destruction, is given by
where the range of ages reflects the uncertainty in the dust yield of such stars. The time needed to create an average Mdust = 2 × 109 M☉ of dust by 170 M☉ Pop III stars is therefore ∼100–400 Myr. This timescale is significantly longer than their main-sequence lifetime and too long to maintain a pristine reservoir of gas to feed the formation of metal-free stars. We therefore conclude that the massive amount of dust seen in AzTEC–3 cannot have been produced by a burst of Pop III stars that are currently giving rise to the observed far-IR luminosity of the galaxy.
6. DISCUSSION
Deriving the physical properties of high-redshift galaxies from their spectra is an important goal for understanding their origin and evolution. In the absence of perfect and complete information, these derived physical properties are not unique, allowing for different evolutionary scenarios to explain the observations.
In this paper we explore a new method that can, in principle, discriminate between different evolutionary scenarios. The method requires that the mass of dust inferred from the far-IR emission to be present in these galaxies be produced within an allotted time span imposed by other observations. We applied our method to the analysis of the ultraluminous submillimeter galaxy AzTEC–3. The basic observations available to constrain the models are the high-resolution UVO spectrum, the UV to millimeter SED, and the CO transitions that were detected from this object.
A fundamental observed quantity is the bolometric luminosity of the galaxy which in AzTEC–3 is dominated by the far-IR emission, and can be used to infer the SFR of the galaxy. However, converting Lbol to an SFR requires knowledge of the stellar IMF and the age of the galaxy, since for a constant SFR and given IMF, the bolometric luminosity increases with time because of the contribution of long-lived low-mass stars to the emission.
The IR region of the spectrum can be used to derive the mass of dust giving rise to the emission. The dust mass depends on the dust composition and the dust temperature distribution, both unknown quantities. Estimates of the dust mass can be derived for pure dust compositions (graphite, amorphous carbon, or silicates), based on dust models for the Milky Way. Alternatively, the dust mass can also be estimated by fitting the observed SED with that of known galaxies for which the dust mass has been derived from detailed models.
The CO emission spectrum can be used to derive the dynamical mass of the galaxy and the mass of molecular gas. Both quantities are uncertain. The dynamical mass estimate depends on the basic premise that the system has relaxed to a Keplerian disk, and the determination of the gas mass depends on the uncertain X-factor that converts the CO line intensity to an H2 mass. Without any further assumptions, only an upper limit can be obtained for the stellar mass in this galaxy.
The high-resolution UV spectrum can be used to determine the age of the system. Massive young stars exhibit characteristic P-Cygni profiles. The presence of these features in the spectrum can therefore be used to determine the age of the system. For an instantaneous burst of star formation, these features disappear as these massive stars die off within ∼30 Myr. So their presence can be used to estimate the age. When star formation is a continuous process, the massive stars are always present, and the strengths of the features relative to the continuum change then much more gradually. Therefore, the age inferred from such spectrum will be much longer than that inferred for an instantaneous burst. From our simulations, assuming a constant SFR, we find that the UV spectrum cannot discriminate between systems that are between ∼30 to 400 Myr old, the relevant age limits for AzTEC–3.
The bolometric luminosity of AzTEC–3 lies between 0.6 and 2 × 1013 L☉, and to be definitive we adopted a value of Lbol=1 × 1013 L☉. For the dust mass we adopted a value of (2 ± 1) × 109 M☉. We adopted a very conservative upper limit of 5 × 1010 M☉ on the stellar mass, based on the estimated dynamical and gas mass from this object.
Using these basic constraints, Lbol, Mdust, and Mstars, we constructed models to follow the evolution of the galaxy's SED, stellar, and dust content as a function of time for seven different stellar IMFs. Viable models were those that were able to produce the inferred mass of dust before the concurrently produced mass in stars exceeded the mass limit. The results of these models are depicted in Figure 10 and described in Section 4.
The results of the paper can be briefly summarized as follows.
- 1.Most of the dust giving rise to the far-IR emission from AzTEC–3 must have been grown in molecular clouds. In the Milky Way, the need for grain growth in molecular clouds to explain the abundance of dust was first pointed out by Dwek & Scalo (1980), and more recently by Zhukovska et al. (2008) and Jones & Nuth (2011). Michałowski et al. (2010b) argued that growth in molecular clouds is the main source of dust in high-redshift (5 < z < 6.5) quasars. In this paper, we have extended this point to the ultraluminous starburst AzTEC–3 using detailed chemical evolution models.
- 2.There is a correlation between grain destruction and the required growth in molecular clouds. At the peak of its dust production, which occurs at ∼200 Myr, the TH model can form the inferred dust mass in AzTEC–3 without resorting to grain growth in the ISM, but with a significantly reduced (compared to the Milky Way) efficiency of grain destruction. At earlier epochs grain growth in molecular clouds must account for a significantly larger fraction of the dust mass. In the Starburst model any grain destruction must be rapidly balanced by accretion in the ISM. The Salpeter and Galzebrook IMFs are less efficient in producing dust than the TH one. They can afford only little grain destruction, and therefore require a large fraction of their dust to be of interstellar origin.
- 3.The model with the TH IMF is the most efficient one for producing the dust in AzTEC–3, and therefore the more favorable model for the origin of the dust in this galaxy. Taking the galaxy's age to be 200 Myr, the longest time allowed by the stellar mass constraint, the constant SFR required to provide the observed bolometric luminosity and the required dust mass is ∼500 M☉ yr−1. The inferred constant SFR is significantly smaller and the age is significantly longer than the previous estimates of ∼1800 M☉ yr−1and ∼30 Myr, respectively, assuming a single burst of star formation.
- 4.The TH model is, however, not unique. A clear discrimination between the different models requires a significant reduction in the uncertainties of the measured quantities. For example, a better determination of the stellar mass may have resulted in a lower value than the upper limit adopted in this paper, ruling out the models using the Salpeter or the Glazebrook functions to characterize the stellar IMF.
- 5.Finally, Pop III stars could not have produced both the currently observed luminosity and dust mass in AzTEC–3.
The studies presented here illustrate the power of multiwavelength observations and the simultaneous use of spectral, stellar, and dust mass constraints for determining the properties of high-redshift galaxies. Our studies also highlight the need for more observational data and improved model input parameters.
E.D. acknowledges helpful discussions with Dominik Riechers. GISMO millimeter observations used in the analysis were supported through NSF grants AST-0705185 and AST-1020981.