Issue |
A&A
Volume 524, December 2010
|
|
---|---|---|
Article Number | A33 | |
Number of page(s) | 20 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361/201015504 | |
Published online | 22 November 2010 |
Mid-infrared spectroscopy of Spitzer-selected ultra-luminous starbursts at z ~ 2⋆
1
UPMC Univ. Paris 06, UMR7095, Institut d’Astrophysique de
Paris, 75014
Paris, France
e-mail: fiolet@iap.fr
2
CNRS, UMR7095, Institut d’Astrophysique de Paris,
75014
Paris,
France
3
Univ. Paris-Sud, Institut d’Astrophysique Spatiale,
UMR8617, 91405
Orsay,
France
4
CNRS, 91405
Orsay,
France
5
Spitzer Science Center, California Institute of
Technology, MS
220-6, Pasadena,
CA
91125,
USA
6
Department of Physics and Astronomy, Rutgers, the State University
of New Jersey, 136 Frelinghuysen
Road, Piscataway
NJ
08854,
USA
7
Max-Planck Institut für extraterrestrische Physik,
Postfach 1312, 85741
Garching,
Germany
8
Department of Physics & Astronomy, University of
Sussex, Falmer,
Brighton
BN1 9RH,
UK
9
University of Minnesota, 116 Church St SE, Minneapolis, MN
55414,
USA
10
North American ALMA Science Center, NRAO,
Charlottesville,
USA
11
National Radio Astronomy Observatory, PO Box O, Socorro, NM
87801,
USA
12
INAF - IASF Milano, via E. Bassini 15, 20133, Italy
13
Herschel Science Center, California Institute of
Technology, 100-22,
Pasadena, CA
91125,
USA
Received:
30
July
2010
Accepted:
26
August
2010
Context.Spitzer’s wide-field surveys and followup capabilities have allowed a new breakthrough in mid-IR spectroscopy up to redshifts ≥ 2, especially for 24 μm detected sources.
Aims. We want to study the mid-infrared properties and the starburst and AGN contributions, of 24 μm sources at z ~ 2, through analysis of mid-infrared spectra combined with millimeter, radio, and infrared photometry. Mid-infrared spectroscopy allows us to recover accurate redshifts.
Methods. A complete sample of 16 Spitzer-selected
sources (ULIRGs) believed to be starbursts at z ~ 2
(“5.8 μm-peakers”) was selected in the (0.5 deg2) J1064+56
SWIRE Lockman Hole field (“Lockman-North”). These sources have
S24 μm > 0.5 mJy,
a stellar emission peak redshifted to 5.8 μm, and
. The entire sample was observed with the low
resolution units of the Spitzer/IRS infrared spectrograph. These sources
have 1.2 mm observations with IRAM 30 m/MAMBO and very deep 20 cm observations from the
VLA. Nine of our sources also benefit from 350 μm observation and
detection from CSO/SHARC-II. All these data were jointly analyzed.
Results. The entire sample shows good quality IRS spectra dominated by strong PAH features. The main PAH features at 6.2, 7.7, 8.6, and 11.3 μm have high S/N average luminosities of 2.90 ± 0.31, 10.38 ± 1.09, 3.62 ± 0.27, and 2.29 ± 0.26 × 1010 L⊙, respectively. Thanks to their PAH spectra, we derived accurate redshifts spanning from 1.750 to 2.284. The average of these redshifts is 2.017 ± 0.038. This result confirms that the selection criteria of “5.8 μm-peakers” associated with a strong detection at 24 μm are reliable to select sources at z ~ 2. We have analyzed the different correlations between PAH emission and infrared, millimeter, and radio emissions. Practically all our sources are strongly dominated by starburst emission, with only one source showing an important AGN contribution. We have also defined two subsamples based on the equivalent width at 7.7 μm to investigate AGN contributions.
Conclusions. Our sample contains strong starbursts and represents a particularly 24 μm-bright class of SMGs. The very good correlation between PAH and far-IR luminosities is now confirmed in high-z starburst ULIRGs. These sources show a small AGN contribution to the mid-IR, around ~20% or less in most cases.
Key words: galaxies: high-redshift / galaxies: starburst / galaxies: active / infrared: galaxies / submillimeter: galaxies / techniques: spectroscopic
Colour figures and Appendices are only available in the electronic form at http://www.aanda.org
© ESO, 2010
1. Introduction
The mid-infrared (MIR) regime is known to provide very rich diagnostics of the interstellar medium (ISM) of galaxies. Hot dust from starbursts and especially AGN contributes to continuum emission. Dust also produces various absorption features, most spectacularly from silicates. However, the MIR spectra of star-forming galaxies are dominated by strong emission features attributed to polycyclic aromatic hydrocarbons (PAHs). These large molecules, typically a few hundred carbon atoms in the ISM, are the carriers of the series of prominent emission bands observed between 3 and 19 μm, especially between 6 and 12 μm (Leger & Puget 1984; Allamandola et al. 1985). PAHs efficiently absorb UV and optical photons from young stars and reemit the energy mostly in the mid-infrared.
After their discovery in the local ISM and nearby star-forming regions, PAH features have proved to be common, not only in various regions of the Milky Way, but in other galaxies ranging from local systems to high-z starbursts (e.g., Tielens 2008). PAH features were the focus of a good fraction of the MIR spectroscopy programs of the infrared space observatories ISO and Spitzer. Important key projects of both missions were devoted to MIR spectroscopy of large samples of galaxies. The ISO Key Project on Nearby Galaxies (Helou et al. 2000) showed an amazing similarity of PAH spectra in various galaxy environments. More detailed studies with the SINGS spectral mapping program (Kennicutt et al. 2003) using the Spitzer Infrared Spectrograph (IRS: Houck et al. 2004) provided further evidence of interband PAH feature strength variations, especially in the presence of weak AGN (Smith et al. 2007).
The high sensitivity of Spitzer/IRS allowed the extension of PAH spectroscopy to z ≳ 2 for samples of hundreds of galaxies that could be observed in the ~5–35 μm spectral range. This led to the measurement of a large number of redshifts of dust-enshrouded Spitzer sources at 1.5 < z < 2.5 in the so-called “redshift desert” of optical spectroscopy. These MIR spectra allow one to distinguish power-law AGN-dominated spectra from PAH starburst-dominated spectra and composite systems. This long series of high-z IRS studies has thus established the criteria for the identification of various types of sources from Spitzer 3.5–24 μm photometry (Fadda et al. 2010; Houck et al. 2005; Yan et al. 2005, 2007; Sajina et al. 2007; Weedman et al. 2006; Dasyra et al. 2009; Hernán-Caballero et al. 2009; Bertincourt et al. 2009). Particular types of sources can be selected using various criteria: dust-enshrouded AGN from strong 24 μm fluxes and/or power-law IRAC 3.6–8.0 μm spectra (Dey et al. 2008; Dey & The Ndwfs/MIPS Collaboration 2009), starbursts from the presence of a visible redshifted 1.6 μm stellar bump (Farrah et al. 2006, 2008; Lonsdale et al. 2009; Huang et al. 2009; Desai et al. 2009), type 1 or type 2 QSOs (Lutz et al. 2008; Martínez-Sansigre et al. 2008), and submillimeter-selected galaxies (SMGs: Valiante et al. 2007; Menéndez-Delmestre et al. 2007, 2009; Pope et al. 2008; Coppin et al. 2010).
Indeed, strong PAH emission is so characteristic of starbursts that the strength of PAH features may be used as a star formation rate indicator. This is justified by the good correlation that has been found between the luminosity of PAH features (mostly at 6.2 and 7.7 μm) and the total IR or far-IR luminosity, which in the absence of an AGN is taken to be a good tracer of star-forming activity (Kennicutt 1998). Such a correlation has been checked for starburst galaxies at both low and high redshift (e.g., Brandl et al. 2006; Pope et al. 2008; Menéndez-Delmestre et al. 2009). It is the origin of the correlation found between the ISO 7.0 μm- or Spitzer 8.0 μm-band flux and the star formation rate (e.g., Elbaz et al. 2002; Wu et al. 2006; Alonso-Herrero et al. 2006; Farrah et al. 2007) in low-z luminous infrared galaxies. It is similarly related to the use of the single-band 24 μm luminosity for determining star formation rates of galaxies at high redshift, especially at z ~ 2 (e.g., Bavouzet et al. 2008; Rieke et al. 2009; Elbaz 2010,Fiolet et al. in prep.).
Such rich diagnostics of the properties of the ISM and star formation at high redshift provided by mid-IR PAH spectral features have justified the many studies performed with the unique capabilities of Spitzer/IRS. We have used this experience to focus one of the last observations of “cold” Spitzer to obtain high-quality mid-IR spectra of z ~ 2 starbursts. Our sample was Spitzer-selected in a field with very rich multiwavelength data. We aimed to select starburst-dominated ultra-luminous IR galaxies at z ~ 2 based on the presence of a rest-frame 1.6 μm stellar bump – caused by the photospheric opacity minimum in giant or supergiant stars (e.g. John 1988; Simpson & Eisenhardt 1999; Weedman et al. 2006) – in the 5.8 μm IRAC band, together with very strong PAH emission implied by a high 24 μm flux density. Using these “5.8 μm-peakers”, we sought to address in the best conditions the various goals allowed by such studies: determination of redshifts from PAH features, accurate enough for allowing CO line observations with current millimeter facilities; detailed examination of the properties of PAH features in such high-z sources, especially their luminosity ratios, and comparison to local systems; estimation of the ratio between PAH bands and the weak underlying continuum that may constrain the modest AGN contribution; and accurate assessment of the correlation of PAH luminosities with total and far-infrared luminosities (and with the radio luminosity). Silicate absorption at 9.7 μm may provide additional information about the distribution of interstellar dust.
In this paper, we present the observations and the results from IRS spectroscopy of a complete sample of 16 sources selected by their IRAC and MIPS fluxes. The sample selection and the observations are described in Sect. 2. The method adopted to reduce these observations is presented in Sect. 3. The results and their analysis are discussed in Sects. 4 and 5. In Sect. 6, we discuss the discrimination between starburst- and AGN-dominated sources, and we make comparisons with the other samples of SMGs. Throughout the paper, we adopt a standard flat cosmology with H0 = 71 kms-1Mpc-1, ΩM = 0.27, and ΩΛ = 0.73 (Spergel et al. 2003).
2. Sample selection and observations
2.1. Sample selection
Our sample is based on the Spitzer/SWIRE
“5.8 μm-peaker” sample of (Fiolet et al.
2009, hereafter F09). Using the 2007 SWIRE internal catalogue, we have selected
the sixteen 24 μm-brightest sources of this sample with
S24 μm > 0.5 mJy.
These sources obey the “5.8 μm-peaker” criteria with
S3.6 μm < S4.5 μm < S5.8 μm > S8.0 μm.
Some of our sources have no detection in the 8.0 μm band, but we assume
that these sources are also “5.8 μm-peakers” if their fluxes at
5.8 μm are greater than the detection limits at 8.0 μm
(~40 μJy). To remove low-redshift interlopers, we also require that
the sources be optically faint, i.e., (Lonsdale et al. 2006). Our sample is complete for these criteria over
0.5 deg2. This area includes the “Lockman-North” field observed with the VLA
(Owen & Morrison 2008; Owen et al. 2009) and Herschel (Oliver et al. 2010). For this sample from F09, we have
1.2 mm flux densities (S1.2 mm) obtained with IRAM/MAMBO
[Institut de Radioastronomie Millimétrique/Max-Planck Millimeter Bolometer] array (Kreysa et al. 1998). 62% (10/16) of our sources are
detected with
S/N > 3 at
1.2 mm (Table 1) and are thus submillimeter
galaxies (SMGs) with
S1.2 mm > 2 mJy
(⟨ S1.2 mm ⟩ = 2.72 ± 0.54 mJy). The six other sources
have an average ⟨ S1.2 mm ⟩ = 1.01 ± 0.68 mJy (see
Table 1 and F09). This sample is thus
representative of a Spitzer-selected subclass of powerful
z ~ 2 ULIRGs (SMGs and OFRGs, Chapman et al. in prep; Magdis et al. 2010), of which there are ~2000 in all
SWIRE fields that peak at 4.5 or 5.8 μm and have
S24 μm > 500 μJy.
Spitzer mid-IR, submillimeter, and radio data of the selected sample.
Nine of our sources detected at 1.2 mm were observed at 350 μm by Kovács et al. (2010). All these sources have been detected, and our 16 sources are also detected at 250–350 μm by Herschel (Fiolet et al. in prep; Magdis et al. 2010). The “5.8 μm-peaker” criterium also implies that our sample is a mass selected-sample and has z ~ 2.
All of our sources benefit from very deep radio data at 1400, 610, and 324 MHz from the VLA and GMRT (Owen & Morrison 2008; Owen et al. 2009,Owen et al. in prep.), with observed flux densities S20 cm > 50 μJy (Table 1 and F09).
Our sample differs from the other samples of ULIRGs essentially in its selection criteria and redshift range. Indeed, Farrah et al. (2008) and Sajina et al. (2007) have selected bright sources at 24 μm. However, the sample of Sajina et al. (2007) is built to select the reddest objects in the [24 μm − 8 μm] vs. [24 μm − R] color–color diagram and is composed of PAH-weak and PAH-strong sources. In the comparisons that follow, we consider only the PAH-strong sources of this sample. Pope et al. (2008), Menéndez-Delmestre et al. (2009), and Valiante et al. (2007) have built samples of SMGs detected at 850 μm. Nevertheless, the sources from Pope et al. (2008) selected with S24 μm > 200 μ Jy are fainter than our sources (S24 μm > 500 μJy). Menéndez-Delmestre et al. (2009) and Valiante et al. (2007) have also selected their sample based on 1.4 GHz flux density. The sample of Desai et al. (2009) is relatively similar to ours, with sources selected by a redshifted 1.6 μm stellar bump in the IRAC bands and S24 μm > 500 μJy. These sources are also optically faint and satisfy the definition of Dust-Obscured Galaxies (DOGs: Dey et al. 2008). The sample of Huang et al. (2009) is also selected with S24 μm > 500 μJy and with color–color criteria yielding star-forming galaxies. Coppin et al. (2010) have built a sample of SMGs selected to be AGN-dominated. While our sample has redshifts from ~1.7 to ~2.3, the other samples span larger redshift ranges. Shi et al. (2009) have selected sources at intermediate redshift around z ~ 1. Farrah et al. (2008) observed sources with z = 1.3−2.2 centered at z ~ 1.7 (“4.5 μm-peakers”). The samples of SMGs from Pope et al. (2008); Menéndez-Delmestre et al. (2009) and Valiante et al. (2007), and the DOG sample of Desai et al. (2009), span redshifts from ~1 to ~3 (0.6 to 3.6 for Menéndez-Delmestre et al. 2009). The sample of star-forming galaxies of Huang et al. (2009) has redshifts between 1.7 and 2.1 centered at z ~ 1.9. Finally, the sample of Coppin et al. (2010) ranges from z ~ 2.45 to ~3.4.
2.2. Observations
The photometric redshifts of our targets, based on the IRAC bands, lie in the range z ~ 1.5−2.7 with uncertainties of ± 0.5 (see F09). In this redshift range, the low-resolution (R ~ 60−127) Long-Low observing modules of IRS, LL1 (19.5−38.0 μm) and LL2 (14.0−21.3 μm), can cover most of the PAH emission features from 6.2 to 11.3 μm in the rest frame. Indeed, the actual redshifts have all proved to lie in the range 1.7 < z < 2.3 (see Sect. 4.1).
This project (GO50119, PI G. Lagache) was observed for a total time of 53 h (41 h on-source) during the last cold campaign of Spitzer, on January 20–22, 2009. We have divided our sample into four groups based on 24 μm flux density. For each group, the on-target exposure time is summarized in Table 2. We chose these exposure times following the guidelines from the IRS IST report, based on IRS low-resolution observations of two faint sources by Teplitz et al.1 and the Spectroscopy “Performance Estimation Tool” (SPEC-PET). We also used our own IRS observations of distant galaxies to check the expected S/N. The integration times have been fixed to achieve a S/N ~ 4–5 per wavelength element in LL1 and LL2. The observations were done in mapping mode to maximize the rejection of “rogue” pixels and (immediately after applying the “skydark” calibrations) to limit contamination by latent charge from bright objects in our “Long” spectra of faint sources.
Summary of Spitzer IRS observations.
3. Data reduction
The observed spectra were processed with version S18.5.0 of the Spitzer Science Center pipeline. The pipeline produces the basic calibrated data (BCD), which are corrected for different effects (e.g., cosmic ray removal, masking of saturated pixels, and linearization of the science data signal ramps).
Applying the IDL code IrsLow, developed by Fadda et al. (2010) and available on-line2, we have extracted the spectra from the BCD. IrsLow is a data reduction code built especially for low-resolution IRS observations. This code makes a background subtraction and extracts the spectra of the target. The background is estimated after the user has manually masked the regions of interest of the data. In practice, we mask the position of the sources on the spectra to estimate and subtract the background. The spectrum is then extracted using a PSF fitting. The code also allows the simultaneous extraction of two different spectra in the case of two nearby sources. The details of using IrsLow and some examples can be found in Fadda et al. (2010). The output of the program is an ASCII file containing the spectrum and its errors. To facilitate our analysis, we make some modifications to each of these files. First, we limit the wavelengths to 20.5 μm through 35 μm for the first order (i.e., LL1) – the errors on the flux outside this range becoming quite large – and to <20.5 μm for the second order (i.e., LL2). Second, we sort the data by wavelength. Finally, we smooth the spectrum in IDL. This filter makes a running mean over five points providing thus an effective spectral channel width of 0.10 μm. To avoid artificial improvement of S/N, the errors on the spectra are kept unchanged by the smooth filter.
4. Results: spectra
The spectra of our 16 sources are shown in Appendix A (see also Figs. 3 and B). These spectra show high-S/N (>5σ) detections of PAH features at 6.2, 7.7, 8.6, and (in many cases) 11.3 μm for z ~ 2 (Table 4).
Considering the uncertainties of our data, all our sources have MIPS S24 μm from the 2007 SWIRE internal catalogue that are compatible with the 24 μm flux density (S24 μm,spectra) that we have directly measured from the spectra after convolution with the MIPS 24 μm filter profile (Fig. A). For 8/16, the difference is less than 10%. These flux densities are reported in Table 1. A similar result was obtained by Pope et al. (2008).
![]() |
Fig. 1 Ratio of S24 μm,spectra deduced from our spectra using the MIPS 24 μm filter, to S24 μm from the 2007 SWIRE internal catalogue, as a function of S24 μm. The black square is the mean for our full sample. The error bars correspond to 1σ. |
Open with DEXTER |
4.1. PAH redshifts
Redshifts.
Applying the method described in Bertincourt et al. (2009), we have determined redshifts from PAH features (zPAH). This method is based on the cross-correlation of 21 templates with our IRS spectra. These templates are spectra of different kinds of sources. Five templates are dominated by PAHs (Smith et al. 2007). 13 are less PAH-dominated local ULIRGs from Armus et al. (2007), including among others Arp 220 and Mrk 231. The remaining three are a radio galaxy, a quasar, and a Wolf-Rayet galaxy. The best estimate of the redshift, zPAH, is the average of redshifts found for all templates after rejection in cases of very high χ2 (χ2/SNR > 8) or too low SNR(<1). These redshifts are reported in Table 3.
Fluxes and equivalent widths for PAH features from PAHFIT (Smith et al. 2007).
![]() |
Fig. 2 PAH spectroscopic redshift zPAH as a function of photometric redshift from F09. The dashed lines correspond to Δlog (1 + z) = ± 0.06 or ± 15%. |
Open with DEXTER |
4.2. Spectral decomposition
We perform spectral decompositions using the PAHFIT code developed by Smith et al. (2007). PAHFIT adjusts some PAH features and spectral lines as a function of rest wavelength (Fig. 3). The PAH features are represented by Drude profiles and the atomic and molecular spectral lines by Gaussians. The decomposition also tries to fit a stellar continuum component (T ⋆ = 5000 K) and dust components for different temperatures (35, 50, 135, 200, 400, 600, 800, 1000, and 1200 K) represented by modified blackbodies. Only the intermediate dust temperatures make significant contributions to the mid-IR continuum. The dust extinction, mostly due to silicate absorption, is also taken into account. In our decompositions, we have adopted a screen geometry, corresponding to a uniform mask of dust absorbing the emission, to model the silicate absorption (Smith et al. 2007). Figure 3 presents an example of the result of the spectral decomposition. The PAH complexes at 6.2, 7.7, 8.6, and 11.3 μm are clearly visible. As proposed by Smith et al. (2007), the 7.7, 8.6, and 11.3 μm features are better fitted by combinations of several features at 7.42, 7.60, and 7.85 μm for the 7.7 μm feature, at 8.33 and 8.61 μm for the 8.6 μm feature, and at 11.23 and 11.33 μm for the 11.3 μm feature. The results of our decompositions are somewhat uncertain, due to the difficulty of estimating the silicate absorption. However, this absorption is known to be weak to moderate in such starburst sources (e.g. Spoon et al. 2007; Farrah et al. 2008; Weedman & Houck 2008; Li & Draine 2001; Yan et al. 2005). Our targets appear to have especially weak absorption making them a good sample for analyzing with this decomposition approach. The spectral decompositions for the entire sample are reported in Appendix B for individual sources, and in Fig. 4 for the stacked spectrum of all 16 sources (see also Fig. 12). It is seen from individual spectra and the stack that the level of the continuum between 5 and 11 μm remains modest compared to the PAH features (except for L10, Sect. 6.1).
![]() |
Fig. 3 Spectral decomposition for source L17. The solid green line is the fitted model. Blue lines above the continuum are PAH features; violet lines are spectral lines. Thin grey lines represent thermal dust continuum components. The thick red line shows the total continuum (stars+dust). The dotted black line shows the extinction (e−τλ, = 1 for no extinction). In this example, the top axis is labelled in observed wavelength. The redshift is zPAH = 1.959. Spectral decompositions for the other sources are presented in Appendix B. |
Open with DEXTER |
In the stacked spectrum (Fig. 4), a weak line is visible at ~9.7 μm and identified as the 0–0 S(3) pure rotational transition of molecular hydrogen at 9.665 μm (J = 5−3) in the PAHFIT spectral decomposition. This line is also visible in the partial stacks of Fig. 12 and may be present in the spectral decompositions of 13/16 sources (Fig. B). The detection of this line is confirmed by the detailed data analysis that we perform in a separate paper (Fiolet et al. in prep.). The very large luminosity of this line, stronger than the CO luminosity, probably implies very strong shocks to excite its upper level ~2500 K above the ground state of H2, although H2 excitation could also be influenced by X-rays or UV from AGN. We also see hints of other lines, although not at significant levels: the S(5) and S(4) transitions of molecular hydrogen at 6.86 μm and 8.08 μm, respectively, as well as the ionized species [Ar III] at 8.98 μm and [S IV] at 10.52 μm.
![]() |
Fig. 4 Spectral decomposition for the stacked spectrum of all 16 sources. The solid green line is the fitted model. The blue lines above the continuum are the PAH features. The narrow violet lines are the spectral lines. The grey lines represent the thermal dust continuum components. The thick red line shows the total continuum (stars+dust). The dotted black line shows the extinction (e−τλ, = 1 if no extinction). |
Open with DEXTER |
For each spectrum, we evaluate the continuum flux density and luminosity density at 5.8 μm (νLν(5.8 μm)). These values are reported in Table 5. There is only a small dispersion of νLν(5.8 μm) for our sources (average = 1.87 × 1011L⊙, median = 1.77 × 1011L⊙, min = 0.71 × 1011L⊙, max = 2.99 × 1011L⊙). In Fig. 5, we plot PAH luminosity at 7.7 μm (L7.7 μm) as a function of νLν(5.8 μm). The PAH luminosity at 7.7 μm is calculated as the integral of the total flux between 7.3 μm and 7.9 μm after subtraction of the continuum contribution (see Eq. (1)). We have also plotted the sample of Sajina et al. (2007). To make a true comparison, we have applied a similar reduction and spectral decomposition to the PAH-rich sources from Sajina et al. (2007). The values of νLν(5.8 μm) found in this case show no significant difference from the published continuum luminosities. We see, in Fig. 5, that our sample shows somewhat weaker νLν(5.8 μm) than the PAH-rich sample of Sajina et al. (2007) or the AGN-dominated sample of Coppin et al. (2010), while the values of L7.7 μm stay in the same range. This is likely due to a stronger AGN contribution in the sample of Sajina et al. (2007). Note also that the ratio L7.7 μm/νLν(5.8 μm) is very consistent with that found for the local starbursts of Brandl et al. (2006).
Luminosities for different wavelengths.
4.3. PAH features
Thanks to the spectral decomposition, we have determined the fluxes (F)
and equivalent widths (EW) for the PAH complexes at 6.2, 7.7, 8.6, and
11.3 μm (Table 4). The flux is
calculated by integrating the spectra between 6.2 and 6.3 μm, 7.3 and
7.9 μm, 8.3 and 8.7 μm, and 11.2 and
11.4 μm, respectively, after subtracting the continuum flux. We
calculate the associated luminosities
L6.2 μm,
L7.7 μm,
L8.6 μm, and
L11.3 μm, reported in Table 5, using: (1)
where DL is the luminosity distance and
with
Sλ and Sc the
total and the continuum flux densities at the considered wavelengths.
The equivalent width is calculated following: (2)All our sources show
L7.7 μm ≳ 5 × 1010 L⊙
and
EW7.7 μm > 2 μm,
i.e., well above the
EW7.7 μm > 0.8 μm
limit for “strong” PAH sources defined by Sajina et al.
(2007). All our sources except L10 also show
EW6.2 μm ≫ 0.2 μm
and
L6.2 μm ≫ 2 × 1010 L⊙.
These values are usually used to separate starburst- from AGN-dominated sources (Armus et al. 2007).
![]() |
Fig. 5 L7.7 μm as a function of the continuum luminosity νLν(5.8 μm). The red diamonds are the sources of our sample. The open (full) black triangle is the stack of our sources with EW7.7 μm > (<)6μm. The full black circle is the composite spectrum of the SMG sample from Menéndez-Delmestre et al. (2009). The black crosses are the PAH-rich sources from Sajina et al. (2007). The black asterisks are the AGN-dominated sources from Coppin et al. (2010). The black line is the best fit of all these samples. The dashed black line is the best fit to the local starbursts from Brandl et al. (2006) (black squares). |
Open with DEXTER |
Our sample has PAH luminosities on average comparable to those of the PAH-rich sources of Sajina et al. (2007) and slightly larger than those of the SMGs observed by Pope et al. (2008) and Menéndez-Delmestre et al. (2009), and of the “4.5 μm-peakers” of Farrah et al. (2008). Such a behavior might be naturally explained by stronger S24 μm than SMG samples samples and higher redshift than Farrah et al. (2008)’s sample. The S/N of our average luminosities are higher than the other samples (Table 6). Indeed, the mean values of L6.2 μm, L7.7 μm, L8.6 μm, and L11.3 μm are 2.90 ± 0.31, 10.38 ± 1.09, 3.62 ± 0.27, and 2.29 ± 0.26 × 1010 L⊙, respectively, for our sample. The stated uncertainties are the standard deviations of the means. These values of luminosities are simple averages. The weighted averages and the measurements from the stacked spectra of our sources give luminosities consistent with the simple averages within the uncertainties. For example, the weighted average of L6.2 μm is (2.70 ± 0.09) × 1010 L⊙, and the value derived from the stacked spectrum is (2.81 ± 0.94) × 1010 L⊙.
As seen in Figs. 6–8, and in Table 7, our PAH luminosities L6.2 μm, L7.7 μm, L8.6 μm, and L11.3 μm have ratios similar to those of SMGs, as well as local starbursts and star-forming galaxies at z ~ 1. The AGN-dominated sample from Coppin et al. (2010) seems to have slightly different ratios, with a higher L7.7 μm/L8.6 μm.
Average luminosities for different samples.
The samples plotted in Figs. 6–8 follow the relations: Note
that these relations are just plain fits of the displayed samples of Figs. 6–8. As these
samples may be biased in various ways, and the properties of the galaxies may vary with
redshift, one should not expect a strong statistical consistency for these relations. The
relations for L6.2 μm and
L11.3 μm are consistent with the
relations found by Pope et al. (2008) for local
starbursts and SMGs.
These relations show that both L7.7 μm/L6.2 μm and L7.7 μm/L11.3 μm are ~4 for our sources at z ~ 2 (see also Figs. 13 and 14 and Table 7), i.e., about twice as large as for the local starbursts of Brandl et al. (2006) and the sources of Farrah et al. (2008) at z ~ 1.7. On the other hand, the ratio L7.7 μm/L8.6 μm is about twice as small in these luminous galaxies at z ~ 2 as in fainter local starbursts (Table 7). This behavior is similar to that of SMGs (Table 7 and Sect. 6.2). The difference with respect to local starbursts could be explained by a difference in the size or ionization of PAH grains (Sects. 6.1, 6.2 and Fig. 14), or by a different extinction due to silicates around 10 μm (affecting the 8.6 and 11.3 μm bands) or ice at ~6 μm (absorbing the 6.2 μm band) (Menéndez-Delmestre et al. 2009; Spoon et al. 2002). However, the effects of extinction seem somewhat contradictory for the different band ratios. The finding of a lower 7.7/8.6 ratio compared to normal, local starbursts, is consistent with a low silicate extinction in our sample, as does the generally low level of extinction in the spectral fits. However, the “larger” 7.7/11.3 ratio that we report is going in the other direction; similarly, our “larger” 7.7/6.2 ratio would imply a larger ice absorption than in local starbursts.
![]() |
Fig. 6 L7.7 μm as a function of L6.2 μm. The red diamonds are the sources of our sample. The open (full) black triangle is the stack of our subsample with EW7.7 μm > (<) 6 μm. The full black circle is the stacked spectrum of the SMG sample from Menéndez-Delmestre et al. (2009). The black circles are the intermediate-redshift star-forming galaxies from Shi et al. (2009). The black stars are the sample of SMGs from Pope et al. (2008). The black crosses are the PAH-rich sources from Sajina et al. (2007). The black asterisks are the AGN-dominated sources from Coppin et al. (2010). The full black squares are the sources from Farrah et al. (2008). The open black squares are the local star-forming galaxies from Brandl et al. (2006). The solid line is our best fit to all these samples. The dashed line is the best fit found by Pope et al. (2008) for their SMGs and the local starbursts. |
Open with DEXTER |
![]() |
Fig. 7 L7.7 μm as a function of L8.6 μm. The symbols and the line are the same as in Fig. 6. |
Open with DEXTER |
![]() |
Fig. 8 L7.7 μm as a function of L11.3 μm. The symbols and the lines are the same as in Fig. 6. |
Open with DEXTER |
L10 is the only source that seems not to follow the relation between L7.7 μm and L6.2 μm, showing too weak a feature at 6.2 μm. However, Table 5 and Fig. 6 show that the deviation of L10 from the correlation remains within the statistical errors of this noisy spectrum. The offset is probably due to higher extinction or a more energetically important AGN (Rigopoulou et al. 1999); see also Sect. 6.1.
4.4. Silicate absorption
The determination of the optical depth of the silicate absorption at 9.7 μm (τ9.7), reported in Table 4, is a product of the general fit of the spectra with PAHFIT. It is based on interstellar extinction similar to that of the Milky Way (Smith et al. 2007). We note that the silicate absorption is weak on average, ⟨ τ9.7 ⟩ = 0.28 (0.50 if we consider only the sources with τ9.7 ≠ 0) vs. 1.09 for the PAH-rich sample of Sajina et al. (2007). Indeed, for 44% (7/16) of our sources, PAHFIT provides τ9.7 = 0. It is clear that all values of extinction (τ9.7) are very uncertain because of the weakness of the continuum and its uncertainty, and because of the lack of measurements at long wavelengths. We thus state only the qualitative conclusion that the average extinction is likely weak in most of our sources. Considering the values of EW6.2 and τ9.7, most of our sources belong to class 1C of the diagram of Spoon et al. (2007, their Figs. 1 and 3. This class is defined by spectra dominated by PAH emission with weak silicate absorption.
5. PAH emission and star formation
As noted above, it is well known that infrared PAH emission is strongly correlated with star formation through excitation and IR fluorescence of PAHs induced by the UV radiation of young stars. The exceptional quality of our PAH spectra at z ~ 2 allows us to further assess the correlation between PAH luminosity and total IR luminosity, and to compare with other starburst samples at z ~ 2 and at lower redshift. The deep radio data allow a precise comparison with the radio luminosity which is also known to trace the starburst IR luminosity well.
5.1. The correlation between IR and PAH luminosities
The correlation between the strength of the PAH emission and the far-IR luminosity in starbursts is well established although not understood in detail. Remarkably, it has been verified among local starburst galaxies over two orders of magnitude in LIR, from 1010 to 1012 L⊙, by Brandl et al. (2006). The approximate proportionality between L7.7 μm and LIR has been extended through the ULIRG range up to 1013 L⊙ and z ≥ 2 by several studies (Shi et al. 2009; Pope et al. 2008; Menéndez-Delmestre et al. 2009). Our sample gives us an opportunity to check how tight the relation is among the strongest starbursts at z ~ 2.
Average luminosities ratios for different samples.
For the nine sources of our sample observed at 350 μm by Kovács et al. (2010) – all SMGs (F09) – we have an
accurate estimate of the infrared luminosity LIR computed by
Kovács et al. (2010), who fit a distribution of
grey-body models with different dust temperatures (multi-temperature model). All these
sources are powerful ULIRGs with LIR close to
1013 L⊙
(average = 8.87 × 1012L⊙,
⟨ ΔLIR ⟩ = 2.81 × 1012L⊙;
see Table 5). The average value of the ratio of IR
and PAH luminosities is (6)It
is seen in Table 7 and Fig. 9 that the
L7.7 μm/LIR
ratio is within a factor ~2 of those for the other high-z samples.
However, it is a factor ~3.5 larger than for less luminous local starbursts as yielded
by Eq. (6). Schweitzer et al. (2006) have also demonstrated the existence of a correlation
between the strength of the PAH emission and the far-IR luminosity for local
starburst-dominated ULIRGs and QSO hosts. Nevertheless, the
⟨ L7.7 μm/LIR ⟩
ratio that they found is probably lower because of the more compact nature of those
objects.
In Fig. 9, we plot
L7.7 μm as a function of
LIR for different samples. We can see that our sample has
LIR comparable to those of the other SMG samples (Pope et al. 2008; Menéndez-Delmestre et al. 2009), and the PAH-rich sources from Sajina et al. (2007). The star-forming galaxies of
Shi et al. (2009) have
LIR values slightly lower than our sample. We also note that
the best fit of all these samples is consistent with the relations found by Pope et al. (2008) and Menéndez-Delmestre et al. (2009) for SMGs and local starbursts (Brandl et al. 2006). The AGN-dominated sample of Coppin et al. (2010) shows
LIR and
L7.7 μm greater than our sources, but
their ratio is not very different from our sample. The sources of Shi et al. (2009) seem to follow a slightly different relation but with
the same slope as the fit found by Menéndez-Delmestre
et al. (2009). The best fit of all these samples, from local starbursts to SMGs,
is (7)This relation is
consistent with L7.7 μm being a good tracer
of high-z starbursts: a strong
L7.7 μm should be the sign of a high
SFR. If we use our relation between L7.7 μm
and LIR and assume the relation from Kennicutt (1998) between SFR and LIR based
on a 0.1–100 M⊙ Salpeter-like IMF, we obtain the following
relation
(8)where SFR is in
units of M⊙ yr-1.
This relation is slightly different from that found by Menéndez-Delmestre et al. (2009). The latter is based on the luminosity at 7.7 μm calculated from the total flux between 7.1 μm and 8.3 μm after continuum substraction, instead of between 7.3 μm and 7.9 μm as in our case. However, the values of L7.7 μm computed with the two definitions are very comparable. Applying Eq. (8), we find on average ⟨ SFR ⟩ = 1260 ± 470 M⊙ yr-1. This value is consistent with the average SFR of 1060 M⊙ yr-1 computed from the radio luminosity in F09, and close to the median SFR of ~1200 M⊙ yr-1 found for SMGs by Menéndez-Delmestre et al. (2009).
5.2. The correlation between 1.4 GHz and PAH luminosities
![]() |
Fig. 9 L7.7 μm as a function of LIR. The symbols are the same as in Fig. 6. The solid line is our best fit to all samples. The dashed line is the best fit found by Pope et al. (2008) from fitting both their SMGs and the local starbursts. The dotted line is the best fit found by Shi et al. (2009). The dot-dashed line is the best fit found by Menéndez-Delmestre et al. (2009). |
Open with DEXTER |
The correlation between radio luminosity (L1.4 GHz) and far-infrared luminosity in star-forming regions and in local starbursts is well known (e.g., Helou et al. 1985; Condon 1992; Crawford et al. 1996; Sanders & Mirabel 1996). This correlation has been confirmed at z ~ 2 by Chapman et al. (2005) (see also F09; Condon 1992; Smail et al. 2002; Yun & Carilli 2002; Ivison et al. 2002; Kovács et al. 2006; Sajina et al. 2008; Younger et al. 2009; Ivison et al. 2010; Magnelli et al. 2010).
![]() |
Fig. 10 L7.7 μm as a function of rest-frame νLν(1.4 GHz). The symbols are the same as in Fig. 6. The black line is our best fit to all these samples. |
Open with DEXTER |
From the correlation between LIR and L7.7 μm (Fig. 9, Eq. (7)), we might expect a very good correlation between PAH and radio luminosities, as previously observed by Pope et al. (2008). Since our sources benefit from exceptionally good radio data at 20 cm, 50 cm, and 90 cm (Sect. 2.1), it is worthwhile to directly check the tightness of this relation.
It is straightforward to compute the rest-frame luminosity at 1.4 GHz
(L1.4GHz = νLν(1.4 GHz))
using, e.g., the formula given by Kovács et al.
(2006): (9)where
DL is the luminosity distance and
α is the radio spectral index defined by the best power law fit,
Sν ∝ να,
between S20 cm, S50 cm, and
S90 cm. We have computed
νLν(1.4 GHz) for our
sample and those of Sajina et al. (2007), Pope et al. (2008), Shi
et al. (2009), and Brandl et al. (2006).
For our sample, we have used α as computed in F09 and reported in
Table 1. For the PAH-rich sample of Sajina et al. (2007), we have used the values published
in Sajina et al. (2008). For the other samples, we
use the average value of α, ⟨ α ⟩ = −0.64 found in
F09, which is consistent with the typical radio spectral index for star-forming galaxies
(e.g., Condon 1992). The values of
νLν(1.4 GHz) are reported
in Table 5.
As expected, we observe a correlation between L1.4GHz and
L7.7 μm. Considering our sample, but also
the local star-forming galaxy sample (Brandl et al.
2006), the sample of SMGs (Pope et al.
2008), the PAH-rich sources from Sajina et al.
(2007), and the intermediate-redshift starbursts (Shi et al. 2009), we obtain the following relation (Fig. 10): (10)The effect of the
positive k-correction and the value of α are
particularly important for the determination of L1.4GHz and
may significantly affect the reliability of the relation stated in Eq. (10). The correlation between
L1.4GHz and
L7.7 μm seems to be more scattered than
the correlation between LIR and
L7.7 μm. It is not surprising
considering the uncertainties on L1.4GHz and the scatter of
the correlation between LIR and
L1.4GHz as shown by e.g. Kovács et al. (2010).
5.3. PAH emission and stellar mass
![]() |
Fig. 11 L7.7 μm as a function of the stellar bump luminosity L1.6 μm. The symbols are the same as in Fig. 6. Only the samples at z ~ 2 are plotted. The line is the best fit found for the “normal” star-forming galaxies of Lu et al. (2003). |
Open with DEXTER |
The stellar bump luminosity at 1.6 μm (νLν(1.6 μm), reported in Table 5, gives an estimate of the stellar mass (M ⋆ ) of our sources (see, e.g., F09; Seymour et al. 2007; Lonsdale et al. 2009). Applying the same mass-to-light ratio as F09 (M ⋆ /νLν(1.6 μm) = 0.7M⊙/L⊙), M ⋆ derived from νLν(1.6 μm) ranges from 1.1 to 3.0 × 1011 M⊙ (average = 1.6 × 1011 M⊙, median = 1.4 × 1011 M⊙).
In Fig. 11, we plot L7.7 μm as a function of νLν(1.6 μm) for our sources and the PAH-rich sample from Sajina et al. (2007). Our sample shows no clear difference from the sources of Sajina et al. (2007) in terms of νLν(1.6 μm). Both samples follow the relation found for “normal” star-forming galaxies (Sajina et al. 2007; Lu et al. 2003). This result implies that L7.7 μm and νLν(1.6 μm) do not trace AGN (Sajina et al. 2007). Since the PAH luminosity is tightly correlated with the star formation rate, the relation between νLν(1.6 μm) and L7.7 μm seems to confirm the existence of a correlation between star formation rate and stellar mass at z ~ 2 as previously observed by Pannella et al. (2009) and Mobasher et al. (2009). However, note that the mass-to-light ratio remains uncertain because of uncertainty in the age of the dominant stellar population.
6. Discussion
6.1. Distinction between starburst and AGN
![]() |
Fig. 12 PAHFIT spectral decomposition of stacked spectra. Upper panel: 7 sources with EW7.7 μm > 6 μm. Lower panel: 9 sources with EW7.7 μm < 6 μm. The solid green line is the fitted model. The blue lines above the continuum are the PAH features. The narrow violet lines are the spectral lines. The grey lines represent the thermal dust continuum components. The thick red line shows the total continuum (stars+dust). The dotted black line shows the extinction (e−τλ = 1 if no extinction). |
Open with DEXTER |
As previously discussed, the relative strength of the PAH features and the mid-IR continuum is a very efficient criterion for distinguishing starburst- from AGN-dominated ULIRGs (e.g., Yan et al. 2007; Sajina et al. 2007; Weedman et al. 2006; Houck et al. 2005). Genzel et al. (1998) and Rigopoulou et al. (1999) use the parameter 7.7 μm-(L/C), defined as the ratio between the peak flux density at 7.7 μm and the continuum flux density at the same wavelength, to distinguish between AGN and starburst. AGN-dominated composite sources show values of 7.7 μm-(L/C) weaker than starburst-dominated systems. An equivalent straightforward indicator is the equivalent width of the PAH features, especially EW7.7 μm, since it is proportional to the ratio of L7.7 μm to the 7.7 μm continuum. The equivalent widths of our sources span a broad range, more than a factor of 5. Despite the large uncertainty in the continuum, e.g., as estimated by PAHFIT, it is clear that EW7.7 μm should be a relevant parameter for discussion of AGN contributions. For our entire sample, we find ⟨ EW7.7 μm ⟩ = 5.64 μm (median = 4.88) vs. 3.04 μm for the stacked spectrum of Menéndez-Delmestre et al. (2009).
Based on the values of EW7.7 μm (see Table 4), we have constructed two subsamples: the first has EW7.7 μm > 6 μm (7 sources, ⟨ EW7.7 μm ⟩ = 8.33 μm) and the second has EW7.7 μm < 6 μm (9 sources, ⟨ EW7.7 μm ⟩ = 3.55 μm). We have produced stacked spectra of these two subsamples, built from the weighted means of the individual spectra scaled to ⟨ F7.7 μm ⟩ . They are plotted in Fig. 12. It is obvious, as expected, that the second subsample has a significantly stronger continuum.
Applying the same method as Menéndez-Delmestre et al. (2009), we estimate an X-ray luminosity (LX) from the 10 μm continuum flux (Table 4) via the Krabbe et al. (2001) relation. Then assuming the LX/LIR found for AGN by Alexander et al. (2005), we estimate the contribution of AGN to LIR. We find that the AGN contribution is ~20% for our whole sample and ~10% and ~25% for the subsamples with EW7.7 μm > 6 μm and with EW7.7 μm < 6 μm, respectively. Menéndez-Delmestre et al. (2009) have found ~32% to be the AGN contribution in SMGs. The difference can be interpreted as a difference in the continuum flux density, with continuum larger by a factor ~2 for the sample of Menéndez-Delmestre et al. (2009).
In Fig. 9 (see also Table 6), it is seen that the two subsamples are not different in term of LIR and follow the correlation between L7.7 μm and LIR found in Sect. 5.1. This result tends to argue that the majority of the IR luminosity is associated with the starburst and not with the AGN. This is not the case in the sample of Menéndez-Delmestre et al. (2009), which shows a smaller L7.7 μm/LIR ratio, implying that the AGN contributes more to LIR (Table 7).
Another possible sign of a strong AGN contribution is a large L7.7 μm/L6.2 μm ratio, as discussed by, e.g., Rigopoulou et al. (1999) and Pope et al. (2008). They have observed such a rise in L7.7 μm versus L6.2 μm for some SMGs, which they attribute to greater silicate absorption in the presence of a more luminous AGN than is seen in most SMGs or starburst galaxies. Stronger silicate absorption significantly affects L8.6 μm and L11.3 μm, while it may lead to an overestimate of L7.7 μm. However, as seen in Figs. 6 and 8, there is no significant difference in the value of L7.7 μm/L6.2 μm between sources with large and small EWs, perhaps because the silicate absorption is small.
We see in Fig. 10 and Table 6 that sources with EW7.7 μm < 6 μm seem to have an average value of νLν(1.4 GHz) somewhat greater than that of sources with EW7.7 μm > 6 μm. Radio excess could be a good indication of AGN strength (Seymour et al. 2008; Archibald et al. 2001; Reuland et al. 2004). However, we note that practically all of this difference is due to two sources in the small-EW subsample, L11 and L17, which are close to the radio loudness limit at z ~ 2 as defined by Jiang et al. (2007) and Sajina et al. (2008): νLν(1.4 GHz) = 3.66 × 107 L⊙. They are thus good candidates for having significant AGN emission. Applying the same method as previously to determinate the AGN contribution, we find 26% and 19% for L11 and L17, respectively.
To summarize: although small, it seems that the AGN contribution is slightly greater in the sources with EW7.7 μm < 6 μm than the systems with EW7.7 μm > 6 μm. Nevertheless, all these sources show strong PAH emission and individually follow more or less the relations found for local starbursts (Brandl et al. 2006). They are probably all starburst-dominated or AGN-starburst composite sources, with at most a weak AGN contribution to the mid-infrared.
Our sample is particularly homogeneous in term of PAH luminosity (see Table 5). Only one source, L10, is different. Indeed, this source shows a spectrum slightly different from the others. Its continuum is stronger and has a larger slope (Fig. B). Moreover, L10 deviates from the correlation between L6.2 μm and L7.7 μm with a low L6.2 μm. L10 also presents the weakest EW6.2 μm, EW7.7 μm, and EW8.6 μm of our sample, with values of 0.10, 2.32, and 0.77 μm, respectively. Such a difference can be explained by a larger AGN contribution to the mid-infrared emission, and also by a low signal-to-noise ratio (see Fig. A). Nevertheless, the fact that L10 follows the other correlations – especially the correlation between νLν(1.4 GHz) and L7.7 μm – seems to prove that this source is starburst-dominated. The same kind of source was observed by Pope et al. (2008).
6.2. Comparison with other studies
The present work is based on a complete sample of 16 sources selected to be starburst “5.8 μm-peakers” that are bright at 24 μm, with S24 μm > 0.50 mJy. The average flux density of our sample is ⟨ S24 μm ⟩ = 0.63 mJy. This is a brighter sample of SMGs compared to that of Pope et al. (2008) (⟨ S24 μm ⟩ = 0.45 mJy) or Menéndez-Delmestre et al. (2009) (⟨ S24 μm ⟩ = 0.33 mJy). On the other hand, the PAH-rich sample of Sajina et al. (2007) is biased toward the brightest 24μm sources, with a mean ⟨ S24 μm ⟩ = 1.19 mJy. Our sample is also brighter than the ULIRGs at z ~ 2 of Fadda et al. (2010) (⟨ S24 μm ⟩ = 0.27 mJy). All these samples have an average redshift z ~ 2, but ours spans a narrower redshift range than the others. Indeed, our sample has zPAH between 1.75 and 2.28, while the samples of Pope et al. (2008) and Menéndez-Delmestre et al. (2009), the ULIRGs from Fadda et al. (2010) and the PAH-rich sample of Sajina et al. (2007) span z = 0.93−2.55, 0.69–3.62, 1.62–2.44, and 0.82–2.47, respectively.
![]() |
Fig. 13 Upper panel: L7.7 μm/L11.3 μm as a function of L6.2 μm/L11.3 μm. Lower panel: L7.7 μm/L11.3 μm as a function of L8.6 μm/L11.3 μm. In both panels, the red diamonds are the sources in our sample. The small black circles are the sample of O’Dowd et al. (2009). The open (full) black triangle is the stack of our sources with EW7.7 μm > (< )6 μm. The black cross shows the average for the PAH-rich sources from Sajina et al. (2007). The black star shows the average for the sample of SMGs from Pope et al. (2008). The open black square shows the average for the local star-forming galaxies from Brandl et al. (2006). The full black square is the average for the sample of “4.5 μm-peakers” from Farrah et al. (2008). The large black circle is based on the composite spectrum of the SMG sample from Menéndez-Delmestre et al. (2009). The black line and the dotted lines are the best linear regression and the 1σ dispersion found by Galliano et al. (2008) for Galactic HII regions, dwarf spirals, and starburst galaxies. |
Open with DEXTER |
The differences between PAH luminosities for all the samples are small, but our PAH luminosities show an higher S/N than the other high-redshift samples (Table 6). Our sample has PAH luminosities in all bands quite comparable to those of Sajina et al. (2007) and about twice as large as those of Pope et al. (2008) and Menéndez-Delmestre et al. (2009).
As regards the PAH-band luminosity ratios, we see in Table 7 that our sample is not very different in term of ratios from the other SMG samples (Pope et al. 2008; Menéndez-Delmestre et al. 2009; Sajina et al. 2007). Nevertheless, Menéndez-Delmestre et al. (2009) find lower ratios for L6.2 μm/LIR and L7.7 μm/LIR. This is probably due to a greater AGN contribution to LIR (Sect. 6.1).
Our sample also shows a higher L7.7 μm/LIR ratio than the lower-redshift sample of Brandl et al. (2006). This difference could be explained by a more extended PAH distribution in our sources. Indeed, if the PAH distribution is extended, a larger fraction of PAH can survive in strong UV radiation field from AGN or young stars. Then L7.7 μm/LIR becomes higher than in compact PAH distributions where PAH are easier destroyed by UV radiation field (e.g., Huang et al. 2009).
In Fig. 13, we plot the different PAH ratios of these samples following O’Dowd et al. (2009). We find that our sample and the other SMG samples show L7.7 μm/L11.3 μm vs. L6.2 μm/L11.3 μm ratios consistent with the sample of O’Dowd et al. (2009). The sample of Brandl et al. (2006) falls in the border of the dispersion for the PAH ratios of Galactic HII regions, dwarf spirals, and starburst galaxies found by Galliano et al. (2008). The plot of L7.7 μm/L11.3 μm vs. L8.6 μm/L11.3 μm shows that all these samples have ratios different from those of O’Dowd et al. (2009) because of stronger emission at 8.6 μm. Our sample does not follow the relation found by Galliano et al. (2008).
![]() |
Fig. 14 Adapted from O’Dowd et al. (2009): L11.3 μm/L7.7 μm as a function of L6.2 μm/L7.7 μm. The symbols are the same as in Fig. 13. The black lines represent the expected ratios for fully neutral or fully ionized PAH molecules of a given number of carbon atoms from the models of Draine & Li (2001). |
Open with DEXTER |
A likely explanation of these differences in PAH ratios is modification of the size or ionization of PAH grains. Figure 14, adapted from O’Dowd et al. (2009), seems to prove that our sample has PAH grains with sizes comparable to those in the other samples of SMGs, but with a slightly lower fraction of ionized grains than the samples of Pope et al. (2008) and (especially) Sajina et al. (2007). Brandl et al. (2006) find smaller and more neutral grains in local starbursts than our sample. The difference could arise due to the strongest ionizing UV radiation from young stars in the most powerful starbursts but less ionising than in the sample of Pope et al. (2008) and Menéndez-Delmestre et al. (2009) probably because of lower AGN contamination.
7. Summary and conclusions
In this paper, we have presented the results of Spitzer/IRS observations of a complete sample of 16 sources selected to be 24 μm-bright and “5.8 μm-peakers.” The spectra obtained show very strong PAH features at 6.2, 7.7, 8.6, and 11.3 μm, along with a weak continuum. Thanks to these exceptionally strong features, we have estimated accurate PAH redshifts that span the range z = 1.75−2.28. Our average PAH redshift is ⟨ zPAH ⟩ = 2.02, and the average error of the individual zPAH is ⟨ ΔzPAH ⟩ = 0.010 (median = 0.007; min = 0.004; max = 0.064). The stacked spectrum of our 16 sources displays evidence for the pure rotational 0–0 S(3) molecular hydrogen line at λ = 9.67 μm, which is confirmed and analyzed in a separate paper.
Thanks to the very good quality of our IRS spectra, we have calculated PAH luminosities L6.2 μm, L7.7 μm, L8.6 μm, and L11.3 μm to have average values 2.90 ± 0.31, 10.38 ± 1.09, 3.62 ± 0.27, and 2.29 ± 0.26 × 1010L⊙, respectively. These luminosities have a S/N higher than the other SMGs samples. We have studied the correlation between L7.7 μm and the other PAH luminosities. All our sources, except perhaps L10, follow a correlation similar to that found for local starbursts, and they are not very different from other samples of SMGs at z ~ 2. We have also confirmed the very good correlation, previously observed, between L7.7 μm and LIR at high redshift, and we have verified its extension to the radio luminosity νLν(1.4 GHz). It also extends more loosely to stellar mass through νLν(1.6 μm).
All these relations, the luminosity ratios, and the equivalent widths allow us to estimate the AGN contribution to the mid-IR luminosity. We conclude that our sources are starburst-dominated and that the AGN contribution is ~20%. This sample is the most pure selection of massive starbursts at z ~ 2 compared to the other z ~ 2 Spitzer-selected samples. The known fact that the equivalent width of PAH features is a good discriminant between starburst-dominated and AGN-dominated sources is reinforced by the study of two subsamples of sources with EW7.7 μm > 6 μm and EW7.7 μm < 6 μm. The subsample with EW7.7 μm < 6 μm shows a larger AGN contribution (~25%) than the EW7.7 μm > 6 μm subsample (~10%), although the AGN contribution is small in both subsamples.
Acknowledgments
We thank Karin Menéndez-Delmestre and Hélène Roussel for their helpful contributions. This work is based primarily on IRS observations and observations made within the context of the SWIRE survey with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under NASA contract. This work includes observations made with IRAM, which is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain). The VLA is operated by NRAO, the National Radio Astronomy Observatory, a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. The CSO is funded by NSF Cooperative Agreement AST-083826. G.L., B.B., A.B. and F.B. had support for this work provided by ANR “D-SIGALE” ANR-06– BLAN-0170 grant.
References
- Alexander, D. M., Bauer, F. E., Chapman, S. C., et al. 2005, ApJ, 632, 736 [NASA ADS] [CrossRef] [Google Scholar]
- Allamandola, L. J., Tielens, A. G. G. M., & Barker, J. R. 1985, ApJ, 290, L25 [NASA ADS] [CrossRef] [Google Scholar]
- Alonso-Herrero, A., Rieke, G. H., Rieke, M. J., et al. 2006, ApJ, 650, 835 [NASA ADS] [CrossRef] [Google Scholar]
- Archibald, E. N., Dunlop, J. S., Hughes, D. H., et al. 2001, MNRAS, 323, 417 [NASA ADS] [CrossRef] [Google Scholar]
- Armus, L., Charmandaris, V., Bernard-Salas, J., et al. 2007, ApJ, 656, 148 [NASA ADS] [CrossRef] [Google Scholar]
- Bavouzet, N., Dole, H., Le Floc’h, E., et al. 2008, A&A, 479, 83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bertincourt, B., Helou, G., Appleton, P., et al. 2009, ApJ, 705, 68 [NASA ADS] [CrossRef] [Google Scholar]
- Bolzonella, M., Miralles, J.-M., & Pelló, R. 2000, A&A, 363, 476 [NASA ADS] [Google Scholar]
- Brandl, B. R., Bernard-Salas, J., Spoon, H. W. W., et al. 2006, ApJ, 653, 1129 [NASA ADS] [CrossRef] [Google Scholar]
- Chapman, S. C., Blain, A. W., Smail, I., & Ivison, R. J. 2005, ApJ, 622, 772 [NASA ADS] [CrossRef] [Google Scholar]
- Condon, J. J. 1992, ARA&A, 30, 575 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- Coppin, K., Pope, A., Menéndez-Delmestre, K., et al. 2010, ApJ, 713, 503 [NASA ADS] [CrossRef] [Google Scholar]
- Crawford, T., Marr, J., Partridge, B., & Strauss, M. A. 1996, ApJ, 460, 225 [NASA ADS] [CrossRef] [Google Scholar]
- Dasyra, K. M., Yan, L., Helou, G., et al. 2009, ApJ, 701, 1123 [NASA ADS] [CrossRef] [Google Scholar]
- Desai, V., Soifer, B. T., Dey, A., et al. 2009, ApJ, 700, 1190 [NASA ADS] [CrossRef] [Google Scholar]
- Dey, A., Soifer, B. T., Desai, V., et al. 2008, ApJ, 677, 943 [NASA ADS] [CrossRef] [Google Scholar]
- Dey, A., & The Ndwfs/MIPS Collaboration. 2009, in ASP Conf. Ser. 408, ed. W. Wang, Z. Yang, Z. Luo, & Z. Chen, 411 [Google Scholar]
- Draine, B. T., & Li, A. 2001, ApJ, 551, 807 [NASA ADS] [CrossRef] [Google Scholar]
- Elbaz, D. 2010, in IAU Symp., 267, 17 [Google Scholar]
- Elbaz, D., Cesarsky, C. J., Chanial, P., et al. 2002, A&A, 384, 848 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Fadda, D., Yan, L., Lagache, G., et al. 2010, ApJ, 719, 425 [NASA ADS] [CrossRef] [Google Scholar]
- Farrah, D., Lonsdale, C. J., Borys, C., et al. 2006, ApJ, 641, L17 [NASA ADS] [CrossRef] [Google Scholar]
- Farrah, D., Bernard-Salas, J., Spoon, H. W. W., et al. 2007, ApJ, 667, 149 [NASA ADS] [CrossRef] [Google Scholar]
- Farrah, D., Lonsdale, C. J., Weedman, D. W., et al. 2008, ApJ, 677, 957 [NASA ADS] [CrossRef] [Google Scholar]
- Fiolet, N., Omont, A., Polletta, M., et al. 2009, A&A, 508, 117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Galliano, F., Madden, S. C., Tielens, A. G. G. M., Peeters, E., & Jones, A. P. 2008, ApJ, 679, 310 [NASA ADS] [CrossRef] [Google Scholar]
- Genzel, R., Lutz, D., Sturm, E., et al. 1998, ApJ, 498, 579 [NASA ADS] [CrossRef] [Google Scholar]
- Helou, G., Soifer, B. T., & Rowan-Robinson, M. 1985, ApJ, 298, L7 [NASA ADS] [CrossRef] [Google Scholar]
- Helou, G., Lu, N. Y., Werner, M. W., Malhotra, S., & Silbermann, N. 2000, ApJ, 532, L21 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
- Hernán-Caballero, A., Pérez-Fournon, I., Hatziminaoglou, E., et al. 2009, MNRAS, 395, 1695 [NASA ADS] [CrossRef] [Google Scholar]
- Houck, J. R., Roellig, T. L., van Cleve, J., et al. 2004, ApJS, 154, 18 [NASA ADS] [CrossRef] [Google Scholar]
- Houck, J. R., Soifer, B. T., Weedman, D., et al. 2005, ApJ, 622, L105 [NASA ADS] [CrossRef] [Google Scholar]
- Huang, J., Faber, S. M., Daddi, E., et al. 2009, ApJ, 700, 183 [NASA ADS] [CrossRef] [Google Scholar]
- Ivison, R. J., Greve, T. R., Smail, I., et al. 2002, MNRAS, 337, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Ivison, R. J., Magnelli, B., Ibar, E., et al. 2010, A&A, 518, L31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jiang, L., Fan, X., Ivezić, et al. 2007, ApJ, 656, 680 [NASA ADS] [CrossRef] [Google Scholar]
- John, T. L. 1988, A&A, 193, 189 [NASA ADS] [Google Scholar]
- Kennicutt, Jr., R. C. 1998, ARA&A, 36, 189 [NASA ADS] [CrossRef] [Google Scholar]
- Kennicutt, Jr., R. C.,Armus, L., Bendo, G., et al. 2003, PASP, 115, 928 [NASA ADS] [CrossRef] [Google Scholar]
- Kovács, A., Chapman, S. C., Dowell, C. D., et al. 2006, ApJ, 650, 592 [NASA ADS] [CrossRef] [Google Scholar]
- Kovács, A., Omont, A., Beelen, A., et al. 2010, ApJ, 717, 29 [NASA ADS] [CrossRef] [Google Scholar]
- Krabbe, A., Böker, T., & Maiolino, R. 2001, ApJ, 557, 626 [NASA ADS] [CrossRef] [Google Scholar]
- Kreysa, E., Gemuend, H.-P., Gromke, J., et al. 1998, in Advanced Technology MMW, Radio, and Terahertz Telescopes, ed. T. G. Phillips, SPIE Conf., 3357, 319 [Google Scholar]
- Leger, A., & Puget, J. L. 1984, A&A, 137, L5 [NASA ADS] [Google Scholar]
- Li, A., & Draine, B. T. 2001, ApJ, 554, 778 [NASA ADS] [CrossRef] [Google Scholar]
- Lonsdale, C. J., Omont, A., Polletta, M. d. C., et al. 2006, in BAAS, 38, 1171 [Google Scholar]
- Lonsdale, C. J., Polletta, M. d. C.,Omont, A., et al. 2009, ApJ, 692, 422 [NASA ADS] [CrossRef] [Google Scholar]
- Lu, N., Helou, G., Werner, M. W., et al. 2003, ApJ, 588, 199 [NASA ADS] [CrossRef] [Google Scholar]
- Lutz, D., Sturm, E., Tacconi, L. J., et al. 2008, ApJ, 684, 853 [NASA ADS] [CrossRef] [Google Scholar]
- Magdis, G. E., Hwang, H., Elbaz, D., et al. 2010, MNRAS, 409, 22 [NASA ADS] [CrossRef] [Google Scholar]
- Magnelli, B., Lutz, D., Berta, S., et al. 2010, A&A, 518, L28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Martínez-Sansigre, A., Lacy, M., Sajina, A., & Rawlings, S. 2008, ApJ, 674, 676 [NASA ADS] [CrossRef] [Google Scholar]
- Menéndez-Delmestre, K., Blain, A. W., Alexander, D. M., et al. 2007, ApJ, 655, L65 [NASA ADS] [CrossRef] [Google Scholar]
- Menéndez-Delmestre, K., Blain, A. W., Smail, I., et al. 2009, ApJ, 699, 667 [NASA ADS] [CrossRef] [Google Scholar]
- Mobasher, B., Dahlen, T., Hopkins, A., et al. 2009, ApJ, 690, 1074 [NASA ADS] [CrossRef] [Google Scholar]
- O’Dowd, M. J., Schiminovich, D., Johnson, B. D., et al. 2009, ApJ, 705, 885 [NASA ADS] [CrossRef] [Google Scholar]
- Oliver, S. J., Wang, L., Smith, A. J., et al. 2010, A&A, 518, L21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Owen, F. N., & Morrison, G. E. 2008, AJ, 136, 1889 [NASA ADS] [CrossRef] [Google Scholar]
- Owen, F. N., Morrison, G. E., Klimek, M. D., & Greisen, E. W. 2009, AJ, 137 [Google Scholar]
- Pannella, M., Carilli, C. L., Daddi, E., et al. 2009, ApJ, 698, L116 [NASA ADS] [CrossRef] [Google Scholar]
- Pope, A., Chary, R.-R., Alexander, D. M., et al. 2008, ApJ, 675, 1171 [NASA ADS] [CrossRef] [Google Scholar]
- Reuland, M., Röttgering, H., van Breugel, W., & De Breuck, C. 2004, MNRAS, 353, 377 [NASA ADS] [CrossRef] [Google Scholar]
- Rieke, G. H., Alonso-Herrero, A., Weiner, B. J., et al. 2009, ApJ, 692, 556 [NASA ADS] [CrossRef] [Google Scholar]
- Rigopoulou, D., Spoon, H. W. W., Genzel, R., et al. 1999, AJ, 118, 2625 [NASA ADS] [CrossRef] [Google Scholar]
- Rowan-Robinson, M., Babbedge, T., Oliver, S., et al. 2008, MNRAS, 386, 697 [NASA ADS] [CrossRef] [Google Scholar]
- Sajina, A., Yan, L., Armus, L., et al. 2007, ApJ, 664, 713 [NASA ADS] [CrossRef] [Google Scholar]
- Sajina, A., Yan, L., Lutz, D., et al. 2008, ApJ, 683, 659 [NASA ADS] [CrossRef] [Google Scholar]
- Sanders, D. B., & Mirabel, I. F. 1996, ARA&A, 34, 749 [NASA ADS] [CrossRef] [Google Scholar]
- Schweitzer, M., Lutz, D., Sturm, E., et al. 2006, ApJ, 649, 79 [NASA ADS] [CrossRef] [Google Scholar]
- Seymour, N., Stern, D., De Breuck, C., et al. 2007, ApJS, 171, 353 [NASA ADS] [CrossRef] [Google Scholar]
- Seymour, N., Dwelly, T., Moss, D., et al. 2008, MNRAS, 386, 1695 [NASA ADS] [CrossRef] [Google Scholar]
- Shi, Y., Rieke, G. H., Ogle, P., Jiang, L., & Diamond-Stanic, A. M. 2009, ApJ, 703, 1107 [NASA ADS] [CrossRef] [Google Scholar]
- Simpson, C., & Eisenhardt, P. 1999, PASP, 111, 691 [NASA ADS] [CrossRef] [Google Scholar]
- Smail, I., Ivison, R. J., Blain, A. W., & Kneib, J.-P. 2002, MNRAS, 331, 495 [NASA ADS] [CrossRef] [Google Scholar]
- Smith, J. D. T., Draine, B. T., Dale, D. A., et al. 2007, ApJ, 656, 770 [NASA ADS] [CrossRef] [Google Scholar]
- Spergel, D. N., Verde, L., Peiris, H. V., et al. 2003, ApJS, 148, 175 [NASA ADS] [CrossRef] [Google Scholar]
- Spoon, H. W. W., Keane, J. V., Tielens, A. G. G. M., et al. 2002, A&A, 385, 1022 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Spoon, H. W. W., Marshall, J. A., Houck, J. R., et al. 2007, ApJ, 654, L49 [NASA ADS] [CrossRef] [Google Scholar]
- Tielens, A. G. G. M. 2008, ARA&A, 46, 289 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Valiante, E., Lutz, D., Sturm, E., et al. 2007, ApJ, 660, 1060 [NASA ADS] [CrossRef] [Google Scholar]
- Weedman, D. W., & Houck, J. R. 2008, ApJ, 686, 127 [NASA ADS] [CrossRef] [Google Scholar]
- Weedman, D., Polletta, M., Lonsdale, C. J., et al. 2006, ApJ, 653, 101 [NASA ADS] [CrossRef] [Google Scholar]
- Wu, Y., Charmandaris, V., Hao, L., et al. 2006, ApJ, 639, 157 [NASA ADS] [CrossRef] [Google Scholar]
- Yan, L., Chary, R., Armus, L., et al. 2005, ApJ, 628, 604 [NASA ADS] [CrossRef] [Google Scholar]
- Yan, L., Sajina, A., Fadda, D., et al. 2007, ApJ, 658, 778 [NASA ADS] [CrossRef] [Google Scholar]
- Younger, J. D., Omont, A., Fiolet, N., et al. 2009, MNRAS, 394, 1685 [NASA ADS] [CrossRef] [Google Scholar]
- Yun, M. S., & Carilli, C. L. 2002, ApJ, 568, 88 [NASA ADS] [CrossRef] [Google Scholar]
Online material
Appendix A: Spectra and best fit-templates
In Fig. A, we present our spectra and the best-fit templates among the 21 templates described in Sect. 4.1 using the method described in Bertincourt et al. (2009). We also plot the 24 μm flux densities from the SWIRE catalogue and the values determined from the IRS spectra convolved with the MIPS 24 μm filter.
![]() |
Fig. A.1 IRS spectra of our sources. The black line is our data. In shaded grey, we plot the 1σ deviation. The green short-dashed line is the best-fit PAH template from Smith et al. (2007). The red diamond is the 24 μm flux density from the SWIRE catalogue. The red square is the 24 μm flux density extracted from the spectrum with the 24 μm filter profile. The pink shape shows the region where the transmission of the 24 μm is maximum. All these spectra are plotted in observed wavelength. The templates are scaled to the “(1 + z) × 7.7 μm” flux. The source ID, the names of the best-fit templates, and zPAH are reported in each panel. |
Open with DEXTER |
![]() |
Fig. A.1 continued. |
Open with DEXTER |
Appendix B: Spectral decomposition of the spectra
In Fig. B, we present the decompositions of our spectra made with the PAHFIT code developed by Smith et al. (2007).
![]() |
Fig. B.1 Decompositions of the spectra of our sources. The solid green line is the fitted model. The blue lines above the continuum are the PAH features. The violet lines are the spectral lines. The thin grey lines represent the thermal dust continuum components. The thick red line shows the total continuum (stars+dust). The dotted black line shows the extinction (e−τλ, = 1 if no extinction). The source ID is reported as the title of each panel. |
Open with DEXTER |
![]() |
Fig. B.1 continued. |
Open with DEXTER |
All Tables
All Figures
![]() |
Fig. 1 Ratio of S24 μm,spectra deduced from our spectra using the MIPS 24 μm filter, to S24 μm from the 2007 SWIRE internal catalogue, as a function of S24 μm. The black square is the mean for our full sample. The error bars correspond to 1σ. |
Open with DEXTER | |
In the text |
![]() |
Fig. 2 PAH spectroscopic redshift zPAH as a function of photometric redshift from F09. The dashed lines correspond to Δlog (1 + z) = ± 0.06 or ± 15%. |
Open with DEXTER | |
In the text |
![]() |
Fig. 3 Spectral decomposition for source L17. The solid green line is the fitted model. Blue lines above the continuum are PAH features; violet lines are spectral lines. Thin grey lines represent thermal dust continuum components. The thick red line shows the total continuum (stars+dust). The dotted black line shows the extinction (e−τλ, = 1 for no extinction). In this example, the top axis is labelled in observed wavelength. The redshift is zPAH = 1.959. Spectral decompositions for the other sources are presented in Appendix B. |
Open with DEXTER | |
In the text |
![]() |
Fig. 4 Spectral decomposition for the stacked spectrum of all 16 sources. The solid green line is the fitted model. The blue lines above the continuum are the PAH features. The narrow violet lines are the spectral lines. The grey lines represent the thermal dust continuum components. The thick red line shows the total continuum (stars+dust). The dotted black line shows the extinction (e−τλ, = 1 if no extinction). |
Open with DEXTER | |
In the text |
![]() |
Fig. 5 L7.7 μm as a function of the continuum luminosity νLν(5.8 μm). The red diamonds are the sources of our sample. The open (full) black triangle is the stack of our sources with EW7.7 μm > (<)6μm. The full black circle is the composite spectrum of the SMG sample from Menéndez-Delmestre et al. (2009). The black crosses are the PAH-rich sources from Sajina et al. (2007). The black asterisks are the AGN-dominated sources from Coppin et al. (2010). The black line is the best fit of all these samples. The dashed black line is the best fit to the local starbursts from Brandl et al. (2006) (black squares). |
Open with DEXTER | |
In the text |
![]() |
Fig. 6 L7.7 μm as a function of L6.2 μm. The red diamonds are the sources of our sample. The open (full) black triangle is the stack of our subsample with EW7.7 μm > (<) 6 μm. The full black circle is the stacked spectrum of the SMG sample from Menéndez-Delmestre et al. (2009). The black circles are the intermediate-redshift star-forming galaxies from Shi et al. (2009). The black stars are the sample of SMGs from Pope et al. (2008). The black crosses are the PAH-rich sources from Sajina et al. (2007). The black asterisks are the AGN-dominated sources from Coppin et al. (2010). The full black squares are the sources from Farrah et al. (2008). The open black squares are the local star-forming galaxies from Brandl et al. (2006). The solid line is our best fit to all these samples. The dashed line is the best fit found by Pope et al. (2008) for their SMGs and the local starbursts. |
Open with DEXTER | |
In the text |
![]() |
Fig. 7 L7.7 μm as a function of L8.6 μm. The symbols and the line are the same as in Fig. 6. |
Open with DEXTER | |
In the text |
![]() |
Fig. 8 L7.7 μm as a function of L11.3 μm. The symbols and the lines are the same as in Fig. 6. |
Open with DEXTER | |
In the text |
![]() |
Fig. 9 L7.7 μm as a function of LIR. The symbols are the same as in Fig. 6. The solid line is our best fit to all samples. The dashed line is the best fit found by Pope et al. (2008) from fitting both their SMGs and the local starbursts. The dotted line is the best fit found by Shi et al. (2009). The dot-dashed line is the best fit found by Menéndez-Delmestre et al. (2009). |
Open with DEXTER | |
In the text |
![]() |
Fig. 10 L7.7 μm as a function of rest-frame νLν(1.4 GHz). The symbols are the same as in Fig. 6. The black line is our best fit to all these samples. |
Open with DEXTER | |
In the text |
![]() |
Fig. 11 L7.7 μm as a function of the stellar bump luminosity L1.6 μm. The symbols are the same as in Fig. 6. Only the samples at z ~ 2 are plotted. The line is the best fit found for the “normal” star-forming galaxies of Lu et al. (2003). |
Open with DEXTER | |
In the text |
![]() |
Fig. 12 PAHFIT spectral decomposition of stacked spectra. Upper panel: 7 sources with EW7.7 μm > 6 μm. Lower panel: 9 sources with EW7.7 μm < 6 μm. The solid green line is the fitted model. The blue lines above the continuum are the PAH features. The narrow violet lines are the spectral lines. The grey lines represent the thermal dust continuum components. The thick red line shows the total continuum (stars+dust). The dotted black line shows the extinction (e−τλ = 1 if no extinction). |
Open with DEXTER | |
In the text |
![]() |
Fig. 13 Upper panel: L7.7 μm/L11.3 μm as a function of L6.2 μm/L11.3 μm. Lower panel: L7.7 μm/L11.3 μm as a function of L8.6 μm/L11.3 μm. In both panels, the red diamonds are the sources in our sample. The small black circles are the sample of O’Dowd et al. (2009). The open (full) black triangle is the stack of our sources with EW7.7 μm > (< )6 μm. The black cross shows the average for the PAH-rich sources from Sajina et al. (2007). The black star shows the average for the sample of SMGs from Pope et al. (2008). The open black square shows the average for the local star-forming galaxies from Brandl et al. (2006). The full black square is the average for the sample of “4.5 μm-peakers” from Farrah et al. (2008). The large black circle is based on the composite spectrum of the SMG sample from Menéndez-Delmestre et al. (2009). The black line and the dotted lines are the best linear regression and the 1σ dispersion found by Galliano et al. (2008) for Galactic HII regions, dwarf spirals, and starburst galaxies. |
Open with DEXTER | |
In the text |
![]() |
Fig. 14 Adapted from O’Dowd et al. (2009): L11.3 μm/L7.7 μm as a function of L6.2 μm/L7.7 μm. The symbols are the same as in Fig. 13. The black lines represent the expected ratios for fully neutral or fully ionized PAH molecules of a given number of carbon atoms from the models of Draine & Li (2001). |
Open with DEXTER | |
In the text |
![]() |
Fig. A.1 IRS spectra of our sources. The black line is our data. In shaded grey, we plot the 1σ deviation. The green short-dashed line is the best-fit PAH template from Smith et al. (2007). The red diamond is the 24 μm flux density from the SWIRE catalogue. The red square is the 24 μm flux density extracted from the spectrum with the 24 μm filter profile. The pink shape shows the region where the transmission of the 24 μm is maximum. All these spectra are plotted in observed wavelength. The templates are scaled to the “(1 + z) × 7.7 μm” flux. The source ID, the names of the best-fit templates, and zPAH are reported in each panel. |
Open with DEXTER | |
In the text |
![]() |
Fig. A.1 continued. |
Open with DEXTER | |
In the text |
![]() |
Fig. B.1 Decompositions of the spectra of our sources. The solid green line is the fitted model. The blue lines above the continuum are the PAH features. The violet lines are the spectral lines. The thin grey lines represent the thermal dust continuum components. The thick red line shows the total continuum (stars+dust). The dotted black line shows the extinction (e−τλ, = 1 if no extinction). The source ID is reported as the title of each panel. |
Open with DEXTER | |
In the text |
![]() |
Fig. B.1 continued. |
Open with DEXTER | |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.