Metallicity and Far-Infrared Luminosity of High Redshift Quasars
We present the results of an exploratory study of broad line region (BLR) metallicity in 34 2.2 z 4.6 quasars with far-infrared (FIR) luminosities (L) from 10 to 10 L. Quasar samples sorted by L might represent an evolutionary sequence if the star formation rates (SFRs) in quasar hosts generally diminish across quasar lifetimes. We use rest-frame ultraviolet spectra from the Sloan Digital Sky Survey to construct three composite spectra sorted by L, corresponding to average SFRs of 4980, 2130 and M yr after correcting for a nominal quasar FIR contribution. The measured N v 1240/C iv 1550 and Si iv 1397+O iv] 1402/C iv 1550 emission line ratios indicate super-solar BLR metallicities in all three composites, with no evidence for a trend with the star formation rate. The formal derived metallicities, Z 5–9 Z, are similar to those derived for the BLRs of other quasars at similar redshifts and luminosities. These results suggest that the ongoing star formation in the host is not responsible for the metal enrichment of the BLR gas. Instead, the BLR gas must have been enriched before the visible quasar phase. These results for high quasar metallicities, regardless of L, are consistent with evolution scenarios wherein visibly bright quasars appear after the main episode(s) of star formation and metal enrichment in the host galaxies. Finally, young quasars, those more closely associated with a recent merger or a blowout of gas and dust, may exhibit tracers of these events, such as redder continuum slopes and higher incidence of narrow absorption lines. With the caveat of small sample sizes, we find no relation between L and the reddening or the incidence of absorption lines.
keywords:quasars: general — quasars: emission lines — galaxies: evolution
Numerous studies have shown a strong correlation between supermassive black hole (SMBH) mass and host galaxy mass over broad mass ranges, implying that there is a close evolutionary relationship between them (Gebhardt et al., 2000; Merritt & Ferrarese, 2001; Tremaine et al., 2002; Marconi & Hunt, 2003; Häring & Rix, 2004; Shields et al., 2006). Luminous quasars at high redshifts, which represent the most massive SMBHs, are experiencing vigorous accretion of a significant portion of the final SMBH mass. This accretion activity is believed to be triggered by global processes (galaxy mergers, interactions or perhaps secular evolution) that also trigger major episodes of star formation in the massive host galaxies, which are also rapidly being assembled at high redshifts. Therefore, the processes of SMBH formation and growth resulting in the quasar phenomenon are directly linked to the birth of massive galaxies (Haehnelt et al., 1998; Richstone et al., 1998; Omont et al., 2001, 2003; Beelen et al., 2006; Cox et al., 2006). However, the nature of the relationship between quasars and galaxy formation is not well understood.
It is widely believed that strong interactions between gas rich galaxies can result in ultra-luminous infrared galaxies (ULIRGs), defined by L10, which have star formation rates (SFRs) of (Houck et al., 1985; Omont et al., 2001, 2003; Flores et al., 2004; Cox et al., 2005; Beelen et al., 2006; Daddi et al., 2007; Cao et al., 2008). A fraction of ULIRGs have been found to contain dust-enshrouded quasars (Sanders & Mirabel, 1996; Lonsdale et al., 2006). There is evidence at both low and high redshift that these embedded quasars are precursors to optically luminous quasars. For example, Cao et al. (2008) find that low redshift quasars with IR luminosities in the ULIRG range fall between optically selected PG quasars and ULIRGs (starbursts) for a variety of mid-IR spectroscopic indicators of starburst and active galactic nucleus (AGN) contributions including: polycyclic aromatic hydrocarbon (PAH) luminosities, fine structure emission line strengths, silicate absorption strengths, spectral slope and mid-IR colour indices. These findings suggest that IR luminous quasars could represent an intermediate stage between a dominant starburst and a dominant AGN phase.
Quasars with high SFRs therefore may be at an intermediate stage between an embedded quasar in a star forming galaxy and a visible quasar in a galaxy where most star formation has ceased (see also Sanders et al. (1988)). At low redshifts, numerous star formation indicators in quasar host galaxies are seen, e.g. high IR luminosities, PAH emission, strong (sub)-mm emission and strong CO emission (Hao et al., 2005; Schweitzer et al., 2006; Farrah et al., 2007a, b; Netzer et al., 2007). At high redshifts (z 2-6), observations at sub-mm and mm wavelengths (rest-frame far-infrared (FIR) to sub-mm, depending on wavelength and redshift range) of optically luminous quasars suggest that up to 30% of quasars also fall within the ULIRG range, similar to the IR quasars studied by Cao et al. (2008) (Carilli et al., 2001; Omont et al., 2001, 2003; Cox et al., 2005; Beelen et al., 2006; Hao et al., 2008). Coppin et al. (2008) compare dynamical, gas and SMBH masses of ten z sub-mm detected quasars and z sub-mm galaxies (SMGs), which are high redshift counterparts to ULIRGs, though less extreme with more evenly distributed star formation instead of a localised starburst (Menéndez-Delmestre et al., 2009). They find that the fainter half of their quasar sample could likely be ‘transition objects’ between SMGs and luminous quasars, based on their SMG-like surface densities and their proximity to the local M/M relation, since luminous quasars tend to lie above this relation and typical SMGs tend to lie below. The highest redshift sub-mm-selected source currently known has also been observed to posses similar stellar and gas masses to this z sample of transition objects, as well as a small AGN contribution (Coppin et al., 2009). These quasars eventually may produce enough energy from accretion to effectively blow out the surrounding gas, halting star formation and clearing the view to an optically bright quasar before halting their own growth (Wyithe & Loeb, 2003; Di Matteo et al., 2005; Hopkins et al., 2008). This regulation of star formation by the accreting SMBH could naturally produce the observed black hole-galaxy mass correlation, and different stages of the process should present different SFRs, with initially high SFRs declining as the AGN becomes more dominant, culminating in a phase of strong SMBH accretion and little or no ongoing star formation (Kauffmann & Haehnelt, 2000; Granato et al., 2004).
FIR luminosities can be used to determine the star formation rates in quasar host galaxies. The FIR emission is due to dust heated either by star formation or quasar emission (see Haas et al. (2003) for discussion). Lutz et al. (2007, 2008) find evidence that FIR emission in quasar host galaxies at all redshifts is caused by dust heated by star formation and not the quasar. The strength of PAH emission, which is found almost exclusively in star forming regions, tightly corresponds to the strength of FIR emission in the same hosts (Dale et al., 2001; Calzetti et al., 2007; Lutz et al., 2008). Beelen et al. (2006) measure FIR emission from six high redshift quasars and derive FIR to radio spectral indexes consistent with local star forming galaxies without AGN. The FIR emission is seemingly uncontaminated by hotter dust potentially heated by the AGN. These results along with others (see for example, Efstathiou & Rowan-Robinson (1995); Serjeant & Hatziminaoglou (2009)) further substantiate the claim that the FIR luminosity in quasars is dominated by star formation and not by the AGN.
FIR luminosity traces ongoing star formation, but past star formation also can be observed indirectly by measuring chemical abundances. Several studies have found that high-redshift quasars typically have metallicities greater than or equal to solar metallicity in the broad emission line region (BLR), which requires significant previous star formation in the host (Hamann & Ferland, 1999; Dietrich et al., 2003; Warner et al., 2004; Nagao et al., 2006). These BLR studies rely on emission line ratios such as N v/C iv and Si iv+O iv]/C iv that are sensitive to the selective enrichment of nitrogen as a secondary element and/or to the decreasing temperatures and increasing metal-line saturations that occur in metal-rich gas. The main result for typically solar or higher metallicities near quasars has been corroborated by independent studies of the narrow emission lines (Groves et al., 2006; Nagao et al., 2006) and narrow absorption lines (Hamann & Ferland, 1999; D’Odorico et al., 2004; Gabel et al., 2006; Simon & Hamann, 2010) in quasar spectra. The metal-rich BLR result is true even for the highest redshifts studied, e.g. Pentericci et al. (2002); Jiang et al. (2007); Juarez et al. (2009), with redshifts out to z = 6.4. There is no known change in metallicity with redshift (see Matsuoka et al. (2009) and references above). Simple chemical evolution models for quasars and elliptical galaxies find that galactic centres tend to be more metal-rich than their halos, and a centralised starburst can enrich the galactic centre to supersolar abundances in a short time () (Friaca & Terlevich, 1998; Hamann & Ferland, 1999; Granato et al., 2001; Hamann et al., 2002; Granato et al., 2004; Hamann et al., 2007; Juarez et al., 2009). These models, combined with the consistently super-solar gas abundances observed in quasar environments, imply that quasars tend to emerge after or near the end of (potentially) short centralised star formation epochs. If the quasar phase emerges when star formation is on the decline, less advanced environments might have higher SFRs and lower abundances (Georgakakis et al., 2009).
The absorption lines in quasar spectra provide additional information about outflows that might be related to the blowout of gas from the host galaxies, and perhaps, about the gaseous remnants of recent galaxy mergers. These absorption lines may be more common in quasars with higher SFRs if mergers and interactions trigger the star formation as in ULIRGs (Weymann et al., 1991; Becker et al., 2000; Richards, 2001; Rupke et al., 2005; Georgakakis et al., 2009). Quasar ‘associated’ C iv absorption lines (AALs), near the emission redshifts with velocity widths less than 500 km sand broad absorption lines (BALs), with velocity widths greater than a few thousand km s, are examples of these potential merger, interaction and outflow signatures. In systems with recent mergers, remnants from the interaction may manifest as a greater incidence of low-velocity ( km s) AALs. About 25% of bright quasars contain AALs with rest equivalent widths (REWs) of 0.3 Å, and around 10% of quasars have BALs (Ganguly et al., 2001; Vestergaard, 2003; Misawa et al., 2003; Trump et al., 2006; Nestor et al., 2008; Wild et al., 2008; Gibson et al., 2008; Rodríguez Hidalgo et al., 2010). A higher incidence of AALs or BALs may occur in the quasars with higher SFRs if these quasars have more recently experienced an interaction and/or there is a progression in quasar outflow characteristics with time.
In this paper, we present an exploratory observational study examining whether SFR in the host galaxies correlates with metallicity in the near-quasar environment. We measure metallicity in the quasar BLR from the rest frame ultraviolet (UV) spectrum and estimate galactic SFRs from the FIR luminosities.
We select quasars for this study from the sample observed with MAMBO at IRAM at 1.2 mm by Carilli et al. (2001) and Omont et al. (2001, 2003) and SCUBA at JCMT at 850 by Isaak et al. (2002), McMahon et al. (1999) and Priddey et al. (2003) and all compiled by Hao et al. (2008). The quasars all were selected to be optically bright with absolute B-band Magnitude M for the Carilli et al. objects, M for the Omont et al. objects and M111They use an Einstein deSitter cosmology, km s. for the McMahon et al. (1999); Isaak et al. (2002) and Priddey et al. (2003) objects. Carilli et al. (2001) observed a representative sample of 41 out of the more than 100 quasars with z 3.6 found as part of the Sloan Digital Sky Survey (SDSS) Galactic Cap and Southern Equatorial Stripe survey. Omont et al. (2001, 2003) observed a random selection of 97 radio quiet, optically luminous sources from the multicolour Palomar Digital Sky Survey available from G. Djorgovski’s web page222http://astro.caltech.edu/~george/z4.qsos with 3.9 z 4.5 and the Veron-Cetty & Veron (2000) catalogue with 1.8 z 2.8. McMahon et al. (1999) selected a small sample of 6 bright radio quiet quasars from the APM survey (Storrie-Lombardi et al., 1996). Isaak et al. (2002) selected a larger sample of the 76 most UV luminous radio quiet quasars known at the time of observation, and Priddey et al. (2003) selected a complimentary sample of 57 quasars from various large surveys.
We cross-reference this sample with the optical quasar spectra in the SDSS data-release 6, finding 116 objects with available spectra. The SDSS spectra from data release 6 have resolution R 2000 and wavelength coverage (Adelman-McCarthy & et al., 2008). We consider only those objects with redshifts between 2.17 and 4.75, compatible with SDSS spectral coverage of the full 1200 – 1600 Å rest-frame wavelength range. Spectra missing regions within the specified wavelength range are further excluded from the analysis. Five spectra with very low signal to noise ratios (S/N) also are rejected. The final sample consists of 34 optical SDSS spectra with a range of sub-mm brightnesses and a redshift range of 2.2 4.6.
The MAMBO 1.2 mm and SCUBA 850 m observations correspond to 160 – 400 m rest wavelengths, depending on quasar redshift. We follow the prescription in Hao et al. (2008) to convert from observed flux into FIR luminosity at 60 , L, in which we assume that a grey-body spectrum describes the rest frame FIR spectral energy distribution (SED) with a dust temperature of 41 K and dust emissivity index of 1.95, as detailed in Priddey & McMahon (2001). We find agreement with Hao et al. (2008) to within 10%.
The L, absolute B-band magnitude (M) and bolometric luminosity (L) for each object are shown in Figure 1, and listed in Table 1. The FIR luminosity, log(L, ranges from 12.1 to 13.4, is dominated by star formation and falls roughly within the ULIRG range. We estimate the star formation rate of each quasar host from L, corrected to exclude the small quasar contribution by assuming the hosts follow the same regression line as typical quasars for L vs. L as shown in Figure 1 of Hao et al. (2008). This corrected L is then used in Hao et al.’s equation 2, which is derived from the Kennicutt star formation rate law (Kennicutt, 1998), and the estimated star formation rates for each object are also listed in Table 1 . We convert the M from the observation papers to our adopted cosmology, and obtain M spanning the range of -26.7 to -29.4, which corresponds to log(L) of 47.2 to 48.3 log(erg s). L measures the quasar black hole accretion luminosity, and is estimated using a bolometric correction factor of 9.74 applied to the monochromatic continuum luminosity L(4400 Å) () following Vestergaard (2004).
We separate the sample into three bins based on L such that each bin contains a similar number of objects: FIR bright quasars with log(L/L) 13.17, FIR intermediate quasars with 12.8 log(L/L) 13.1 and FIR faint quasars with log(L/L) 12.75, in which the FIR faint luminosities are upper limits. Any objects with log(L/L) upper limits above 12.75 were not included in the study. The FIR-bright and FIR-faint bins each span the redshift range from 2.2 to 4.4. The FIR bright bin also includes one object with z 4.6, and the FIR intermediate bin spans a smaller redshift range from 3.7 to 4.4. The divisions between the L bins are represented in Figure 1 by horizontal dashed lines and average values in each bin for L and L are denoted by filled triangles. These average values are listed in Table 2 for each luminosity bin in columns 2 and 3, along with the corresponding star formation rates in column 4, and several abundance indicators, discussed in §’s 3.3 and 4.
|Composite spectra||log L||log L||Flux Ratio||Z/Z||Flux Ratio||Z/Z|
|log(erg s)||log(erg s)||Myr||N v/C iv||N v/C iv||Si iv+O iv]/C iv||Si iv+O iv]/C iv|
3.1 Composite Spectra
To make comparisons between different L, we create one composite spectrum from the SDSS spectra comprising each of the three FIR bins in the sample described in § 2. To create the composites, we begin by shifting each spectrum to its rest wavelength using the redshifts provided by SDSS (SubbaRao et al., 2002). We manually inspect each individual spectrum and remove absorption features by interpolating across the affected regions. The presence of narrow absorption lines and their relationship to L is discussed in § 3.4 below. Then we average together the spectra in each bin to create the final composites. To ensure that no single spectrum is dominating the composite we also compute the median spectra for each bin, compare them to the average spectra used throughout the rest of the analysis, and confirm that the two composite types are well-matched.
Our analysis is limited to the spectral region between Ly 1216 Å and C iv 1550 Å. These limits are imposed by the incomplete spectral coverage at longer wavelengths and by suppression in the Ly forest starting in the blue wing of the Ly emission line and extending to shorter wavelengths. Before normalising the composite spectra, we visually compare the continuum slopes of the three composites. We find no significant trend for reddening with FIR luminosity in the continuum slopes of the composite spectra. Differences between the composite slopes are negligible compared to the dispersion of slopes among the individual quasars making up each composite.
We fit a power-law continuum to the FIR-intermediate composite spectrum using wavelength regions devoid of absorption or emission, 25 Å wide and centred at 1460, and 1770, following Warner et al. (2003). Unfortunately, the FIR-faint and FIR-bright composites have incomplete spectral coverage past 1700 Å, so we fit the continuum using the available wavelength region around 1460 Å, plus the region from 1335-1355 Å, which is also used by Juarez et al. (2009), and which we determine by visual inspection to be devoid of absorption or any detectable emission. The normalised average composite spectra are shown in Figure 2.
3.2 Emission Line Flux Ratios
We measure several quasar broad emission line fluxes to calculate line flux ratios, which we then use to estimate the gas phase metallicity in the near-quasar environment. Flux ratios are measured for the emission lines that have been shown by e.g. Hamann et al. (2002); Juarez et al. (2009), to give good abundance estimates within the limited wavelength range from 1216 Å to 1600 Å: N v 1240/C iv 1549 and Si iv 1397+O iv] 1402/C iv 1549. These emission lines are labelled in Figure 2. Other ratios such as N III] 1750/O III] 1664, N v 1240/He II 1640 or N III 991/C III 977 are too weak given the S/N in the spectrum and/or they fall in a poorly characterised region of the spectrum. The measured line flux ratios are listed in columns 5 and 7 of Table 2.
Because the N v emission is strongly blended with the Ly emission, which is itself significantly degraded by absorption in the Ly forest, C iv is the only strong, non-blended emission line in the spectra. Thus, the C iv emission lines are fit with gaussians and used as templates for the other emission lines. The fits to C iv and the scaled fits to Ly and N v are shown for the three composites in Figure 3. The fits are performed using GATORPLOT, an IDL program written by C. Warner333http://www.astro.ufl.edu/~warner/GatorPlot/, and are carried out using the smallest number of gaussians possible to provide a good match to the data: one broad and one narrow component for each of the two doublet lines at 1548 and 1551 Å. The two gaussians comprising the 1548 emission line fit are held at fixed full width at half maximum (FWHM) and are scaled in total flux at the fixed central wavelength of the Ly 1216 Å and N v 1239, 1243 Å doublet emission to match these line strengths. Because C iv is known to be blue-shifted relative to Ly and the nominal quasar redshift by 310 km s (Tytler & Fan, 1992), the locations of the N v and Ly centroids are shifted from their laboratory values by this amount relative to C iv.
To accurately determine the N v flux, the red Ly wing beneath the N v emission must be characterised. To do this precisely, we scale the C iv fit to align with the non-absorbed region of Ly, which in some cases is quite small ( Å). The scaled fit may even rise above the data where heavy Ly forest absorption has eaten into the emission, as is the case for all of the three composites, or miss the peak of Ly as is the case in the FIR faint composite; both cases are clearly shown in Figure 3. In order to better match the Ly emission peak in the FIR faint composite, we allow the Ly gaussian to be narrower than the C iv gaussian. Narrowing the gaussian for the Ly emission line does not have a noticeable effect on the resulting N v strength. Regardless, a precise match to the Ly emission line is not required for this analysis, and does not effect the final N v emission results (see e.g. Baldwin et al. (2003) for further discussion).
The C iv and N v doublets are fit at their respective fixed separations and with a 3:2 intensity ratio (halfway between the allowed 2:1 and 1:1 intensity ratios) for the shorter and longer wavelength respectively, as in Baldwin et al. (2003). Possible broadening in N v is predicted if the N v emission forms at a larger distance from the quasar than the C iv emission. Peterson et al. (2004) and Peterson (2008) find that N v forms twice as close to the quasar as C iv, so if the line widths are controlled by virial motions, N v could be up to times broader than C iv. Ly, on the other hand, forms at about the same distance from the quasar as C iv and should be no broader (and probably narrower) than C iv (Peterson et al., 2004; Peterson, 2008). Therefore, the N v profiles we adopt, which have the same FWHM’s as C iv, may be narrower than the actual N v FWHM’s, while the Ly FWHM could be slightly broader. We estimate the uncertainties in the N v line strengths due to uncertainty in the continuum placement by repeating the line fits with different reasonable continuum heights. We also include estimates of N v line strength uncertainties due to blending with Ly and the poorly defined N v FWHM by repeating the line fits with a range of different N v and Ly FWHM’s. The overall uncertainties, combining these effects, are less than a factor of 1.5, so that the measurement errors in the N v/C iv ratios are also less than a factor of 1.5.
The Si iv+O iv] emission line is weaker than the C iv and N v emission lines in the composite spectra, however it also is isolated from other emission. This isolation means that measuring the equivalent width of the emission without fitting, using IRAF444IRAF is the Image Reduction and Analysis Facility, supported by NOAO and AURA Inc., produces accurate fluxes. We confirm the accuracy of the equivalent width measurements by comparing the previously measured fluxes from gaussian fits to C iv to the fluxes from C iv equivalent widths, and find compatible results. We measure Si iv+O iv] and C iv equivalent widths to calculate the Si iv+O iv]/C iv line flux ratios listed in Table 2. To estimate uncertainties in the Si iv+O iv] line strengths, we vary the equivalent width measurements by using a range of continuum heights and wavelength cutoffs for the edge of the emission lines. The overall uncertainties in the Si iv+O iv] line strengths are less than a factor of 1.2, so that the measurement errors in the Si iv+O iv]/C iv ratios are also less than a factor of 1.2.
We convert the N v/C iv flux ratios into metallicities using the theoretical relationship in which nitrogen abundance increases relative to carbon as metallicity increases because of secondary enrichment (CNO nucleosynthesis) processes, and determine the average metallicity for each L composite (Shields, 1976; Hamann & Ferland, 1993, 1999). The correlation is characterised by Hamann et al. (2002) with recent updates to the solar abundance ratios by Dhanda et al. (2007). We find supersolar metallicities for all three composite spectra, as listed in column 6 of Table 2.
Similar characterisations for the Si iv+O iv]/C iv flux ratio-metallicity correlation are performed by Nagao et al. (2006), using the older solar abundances also used by Hamann et al. (2002). We apply the latest corrections for the solar abundance ratios for this characterisation and determine that the metallicities are also supersolar according to this ratio, as shown in column 8 of Table 2.
There is no significant trend among the three composites in the line ratios, based both on our measurements of these ratios and on a visual inspection of the stacked composites in Figure 2, and correspondingly there is no significant trend in metallicity. The average metallicities across the three composites are Z9.5 Z for the N v/C iv ratio and 4.2 Z for the Si iv+O iv]/C iv ratio.
We note that the metallicities inferred from the two line ratios can differ by as much as a factor of 2 in the same spectrum. This is typical of the dispersion found between different line ratios in other studies, and it might be a good indication of the theoretical uncertainties (Dietrich et al., 2003; Nagao et al., 2006; Hamann & Simon, 2009). We compare our metallicities to the metallicities in the L range 10 to 10 erg s sampled by a large emission line study by Warner et al. (2004), who find metallicities of Z 3-6 Z, broadly consistent with the metallicities found in this sample. We also note that the results based on N v/C iv(and N v/He ii) tend to be higher than other line ratios, and thus our best guess at the metallicity in the current sample overall would be near the lower end of the range 4–9 Z (Hamann et al., 2002; Baldwin et al., 2003; Nagao et al., 2006). The relative metallicities between the three composites, which smooth over object-to-object scatter, are robust for differences greater than a factor of 2-3, based on the uncertainties in the line ratio measurements plus the theoretical uncertainties, and therefore useful for spotting strong metallicity trends among the different L bins.
3.4 Absorption Lines
We inspect each SDSS spectrum used to make the composites and count the number of C iv AALs with REW Å per spectrum. We find that 17/34 quasars have at least one AAL and 6/34 have absorption too broad to be classified as an AAL, yet too narrow to be classified as a BAL, with widths more than twice the thermal width and velocity dispersion expected for unrelated gas in the line of sight, which is indicative of gas forming in an outflow (Bahcall et al., 1967; Young et al., 1982). These statistics are broadly consistent with previous work, given the small numbers involved in this study (Vestergaard, 2003; Trump et al., 2006; Nestor et al., 2008; Wild et al., 2008; Gibson et al., 2008; Rodríguez Hidalgo et al., 2010). AALs and outflow lines appear in each of the three L bins, with no significant differences in their occurrence fractions among the bins that would indicate a trend with L, particularly given the small number of objects in each bin.
The average line flux ratios N v/C iv and Si iv+O iv]/C iv across the three composites, as shown in columns 5 and 7 of Table 2, correspond to average metallicities of Z 9.5 and 4.2 Z respectively, shown in columns 6 and 8 of Table 2. There is no significant trend in metallicity with L in this sample. We note that the FIR bright bin has the most luminous average L, while the FIR faint and intermediate bins have roughly equal, and less luminous average L. This difference in L could appear in the relative metallicities among the bins because more luminous quasars tend to be more metal-rich than less luminous quasars (Hamann & Ferland, 1999; Warner et al., 2004; Shemmer et al., 2004). However, the absence of a corresponding gradient in the metallicity is not surprising given the relatively high average of the total L, the range in average L among the composites of less than 0.4 dex (Table 2), which corresponds to a range in metallicity of no more than 0.2 dex in Warner et al. (2004), and the sensitivity of the data to trends in metallicity greater than 0.3 dex between the three L bins.
The fact that the data show no significant metallicity evolution with changing SFR (L) could be due to the small sample sizes being affected by object-to-object scatter or by the small dynamic range in L, with values that are all within the luminosity/SFR range of powerful ULIRGs.
Absorption features such as AALs and BALs and reddened continua may be more numerous among quasar environments that have undergone a recent merger or are participating in the blowout of gas that is revealing the visible quasar (Becker et al., 2000; Richards, 2001; D’Odorico et al., 2004; Georgakakis et al., 2009). Although these are optically selected quasars chosen to be bright in the rest-frame UV and strongly reddened sources would be excluded, there could still be a more subtle trend in sub-samples by L. We do not find a trend in the number of AALs per quasar or in the continuum slope per L bin. Each bin appears to have roughly the same percentage of quasars with AALs. The absence of a trend could be due simply to the small numbers in our sample. Or, if the consistency across bins is real, it could indicate that there is no progression in quasar outflow characteristics with changing SFRs. A larger sample size is needed to resolve this question.
The overall conclusion from this and other studies is that quasars appear to be metal-rich at all redshifts, with only small variations due to SMBH mass and luminosity dependencies (Warner et al., 2004; Shemmer et al., 2004; Nagao et al., 2006; Jiang et al., 2007; Juarez et al., 2009). The consistently high metallicities suggest that the star formation producing these FIR-bright quasars is not the star formation that determines the metallicity of the BLR. Instead, the BLR gas must have been enriched prior to these FIR-producing star formation episodes, possibly shortly after a merger or some other star formation trigger when the starburst first began. There are no significant differences among the L bins in the outflow fraction or in the amount of debris leftover from mergers in the form of absorption, which might correlate with L if this is indeed an evolutionary sequence. This scenario is consistent with models of galaxy and quasar evolution, which predict a lack of metallicity evolution. The models also find that quasar activity, though caused by a disruption of some sort, is significantly delayed after the disruption event and most of the star formation and enrichment in the host is complete even as the quasar phase first begins to emerge (Di Matteo et al., 2005; Hopkins et al., 2008; Li et al., 2008). The FIR-bright quasars may be in some sense less mature than the FIR-faint quasars, however, this is not manifest in significant metallicity evolution over the timescales probed by this sample, which are shorter than a quasar lifetime.
The reality of quasar-host galaxy formation could also be more complicated than any simple monotonic sequence in which high SFRs correlate cleanly with lower metallicities and an earlier stage of evolution. Veilleux et al. (2009) suggest a ‘softer’ ULIRG-quasar evolution paradigm with more scatter in the path of an individual galaxy’s evolution from ULIRG to quasar. They measure neon abundances of Z2.9 Z in the nuclear regions of ULIRGs. If these high metallicities are ubiquitous in ULIRGs, we will never observe metal-poor quasars, regardless of their stage of evolution. Other plausible scenarios suggest that AGN and star formation activity in galaxies may be episodic in nature (Di Matteo et al., 2008). Thus, as a quasar emerges at the tail end of a star-forming phase, there is no guarantee that its BLR metallicity is linked to this recent and/or ongoing star formation phase, because several previous episodes of star formation and AGN activity could already have occurred and enriched the BLR gas to supersolar values.
An alternative interpretation is provided by the works of Davé & Oppenheimer (2007) and Finlator & Davé (2008). They postulate a scenario where a star-forming galaxy could quickly reach an equilibrium metallicity early in its formation history. Higher mass galaxies would have higher abundances, but for a given galaxy, the abundance would not change significantly after an initial balance of in-flowing metal-poor and out-flowing metal-rich gas was reached. The lack of BLR metallicity evolution with star formation rate found in this study is consistent with these objects having previously reached an equilibrium metallicity, unaffected by the observed ongoing star formation.
We create composite rest-frame UV spectra for a sample of 34 quasars observed in the FIR and measure N v/C iv and Si iv+O iv]/C iv flux ratios. We convert these flux ratios into a metallicity for each of three L bins, which we interpret as SFR bins. We find that all three bins have supersolar metallicities of the order 4 to 9 times solar, and find no metallicity trend with L. We investigate the amount of outflowing gas, as manifest by both AALs and somewhat broader absorption lines, in each L bin and find no difference in number between the L bins, neither do we find a change in continuum slope with changing L, as would be consistent with more dusty star formation in hosts with higher SFRs. Together, these results suggest that the ongoing star formation in the host is not responsible for the metal enrichment of the BLR gas. Instead, the BLR gas must have been enriched before the visible quasar phase. We note that we require a substantially larger sample size with a broader range in L to confirm that these results are not affected by unusual individual objects or biased by the narrow L range.
- Adelman-McCarthy & et al. (2008) Adelman-McCarthy J. K., et al. 2008, ApJS, 175, 297
- Bahcall et al. (1967) Bahcall J. N., Sargent W. L. W., Schmidt M., 1967, ApJ, 149, L11+
- Baldwin et al. (2003) Baldwin J. A., Hamann F., Korista K. T., Ferland G. J., Dietrich M., Warner C., 2003, ApJ, 583, 649
- Becker et al. (2000) Becker R. H., White R. L., Gregg M. D., Brotherton M. S., Laurent-Muehleisen S. A., Arav N., 2000, ApJ, 538, 72
- Beelen et al. (2006) Beelen A., Cox P., Benford D. J., Dowell C. D., Kovács A., Bertoldi F., Omont A., Carilli C. L., 2006, ApJ, 642, 694
- Calzetti et al. (2007) Calzetti D., Kennicutt R. C., Engelbracht C. W., Leitherer C., Draine B. T., Kewley L., Moustakas J., Sosey M., Dale D. A., Gordon K. D., Helou G. X., Hollenbach D. J., Armus L., Bendo G., Bot C., Buckalew B. e. a., 2007, ApJ, 666, 870
- Cao et al. (2008) Cao C., Xia X. Y., Wu H., Mao S., Hao C. N., Deng Z. G., 2008, MNRAS, 390, 336
- Carilli et al. (2001) Carilli C. L., Bertoldi F., Rupen M. P., Fan X., Strauss M. A., Menten K. M., Kreysa E., Schneider D. P., Bertarini A., Yun M. S., Zylka R., 2001, ApJ, 555, 625
- Coppin et al. (2009) Coppin K. E. K., Smail I., Alexander D. M., Weiss A., Walter F., Swinbank A. M., Greve T. R., Kovacs A., De Breuck C., Dickinson M., Ibar E., Ivison R. J., Reddy N., Spinrad H., Stern D., Brandt W. N., Chapman S. C. e. a., 2009, MNRAS, 395, 1905
- Coppin et al. (2008) Coppin K. E. K., Swinbank A. M., Neri R., Cox P., Alexander D. M., Smail I., Page M. J., Stevens J. A., Knudsen K. K., Ivison R. J., Beelen A., Bertoldi F., Omont A., 2008, MNRAS, 389, 45
- Cox et al. (2005) Cox P., Beelen A., Bertoldi F., Omont A., Carilli C. L., Walter F., 2005, in Wilson A., ed., ESA Special Publication Vol. 577 of ESA Special Publication, Gas and dust in high redshift quasars. pp 115–120
- Cox et al. (2006) Cox T. J., Jonsson P., Primack J. R., Somerville R. S., 2006, MNRAS, 373, 1013
- Daddi et al. (2007) Daddi E., Dickinson M., Morrison G., Chary R., Cimatti A., Elbaz D., Frayer D., Renzini A., Pope A., Alexander D. M., Bauer F. E., Giavalisco M., Huynh M., Kurk J., Mignoli M., 2007, ApJ, 670, 156
- Dale et al. (2001) Dale D. A., Helou G., Contursi A., Silbermann N. A., Kolhatkar S., 2001, ApJ, 549, 215
- Davé & Oppenheimer (2007) Davé R., Oppenheimer B. D., 2007, MNRAS, 374, 427
- Dhanda et al. (2007) Dhanda N., Baldwin J. A., Bentz M. C., Osmer P. S., 2007, ApJ, 658, 804
- Di Matteo et al. (2008) Di Matteo T., Colberg J., Springel V., Hernquist L., Sijacki D., 2008, ApJ, 676, 33
- Di Matteo et al. (2005) Di Matteo T., Springel V., Hernquist L., 2005, Nature, 433, 604
- Dietrich et al. (2003) Dietrich M., Hamann F., Shields J. C., Constantin A., Heidt J., Jäger K., Vestergaard M., Wagner S. J., 2003, ApJ, 589, 722
- D’Odorico et al. (2004) D’Odorico V., Cristiani S., Romano D., Granato G. L., Danese L., 2004, MNRAS, 351, 976
- Efstathiou & Rowan-Robinson (1995) Efstathiou A., Rowan-Robinson M., 1995, MNRAS, 273, 649
- Farrah et al. (2007a) Farrah D., Bernard-Salas J., Spoon H. W. W., Soifer B. T., Armus L., Brandl B., Charmandaris V., Desai V., Higdon S., Devost D., Houck J., 2007, ApJ, 667, 149
- Farrah et al. (2007b) Farrah D., Lacy M., Priddey R., Borys C., Afonso J., 2007, ApJ, 662, L59
- Finlator & Davé (2008) Finlator K., Davé R., 2008, MNRAS, 385, 2181
- Flores et al. (2004) Flores H., Hammer F., Elbaz D., Cesarsky C. J., Liang Y. C., Fadda D., Gruel N., 2004, A&A, 415, 885
- Friaca & Terlevich (1998) Friaca A. C. S., Terlevich R. J., 1998, MNRAS, 298, 399
- Gabel et al. (2006) Gabel J. R., Arav N., Kim T.-S., 2006, ApJ, 646, 742
- Ganguly et al. (2001) Ganguly R., Bond N. A., Charlton J. C., Eracleous M., Brandt W. N., Churchill C. W., 2001, ApJ, 549, 133
- Gebhardt et al. (2000) Gebhardt K., Bender R., Bower G., Dressler A., Faber S. M., Filippenko A. V., Green R., Grillmair C., Ho L. C., Kormendy J., Lauer T. R., Magorrian J., Pinkney J., Richstone D., Tremaine S., 2000, ApJ, 539, L13
- Georgakakis et al. (2009) Georgakakis A., Clements D. L., Bendo G., Rowan-Robinson M., Nandra K., Brotherton M. S., 2009, MNRAS, 394, 533
- Gibson et al. (2008) Gibson R. R., Brandt W. N., Schneider D. P., Gallagher S. C., 2008, ApJ, 675, 985
- Granato et al. (2004) Granato G. L., De Zotti G., Silva L., Bressan A., Danese L., 2004, ApJ, 600, 580
- Granato et al. (2001) Granato G. L., Silva L., Monaco P., Panuzzo P., Salucci P., De Zotti G., Danese L., 2001, MNRAS, 324, 757
- Groves et al. (2006) Groves B. A., Heckman T. M., Kauffmann G., 2006, MNRAS, 371, 1559
- Haas et al. (2003) Haas M., Klaas U., Müller S. A. H., Bertoldi F., Camenzind M., Chini R., Krause O., Lemke D., Meisenheimer K., Richards P. J., Wilkes B. J., 2003, A&A, 402, 87
- Haehnelt et al. (1998) Haehnelt M. G., Natarajan P., Rees M. J., 1998, MNRAS, 300, 817
- Hamann & Ferland (1993) Hamann F., Ferland G., 1993, ApJ, 418, 11
- Hamann & Ferland (1999) Hamann F., Ferland G., 1999, ARA&A, 37, 487
- Hamann et al. (2002) Hamann F., Korista K. T., Ferland G. J., Warner C., Baldwin J., 2002, ApJ, 564, 592
- Hamann & Simon (2009) Hamann F., Simon L., 2009, in Cunha K., Spite M., Barbuy B., eds, Chemical Abundances in the Universe: Connecting First Stars to Planets Proceedings IAU Symposium No. 265, Quasar Metal Abundances & Host Galaxy Evolution
- Hamann et al. (2007) Hamann F., Warner C., Dietrich M., Ferland G., 2007, in Ho L. C., Wang J.-W., eds, The Central Engine of Active Galactic Nuclei Vol. 373 of Astronomical Society of the Pacific Conference Series, Quasar Metal Abundances and Host Galaxy Evolution. pp 653–+
- Hao et al. (2005) Hao C. N., Xia X. Y., Mao S., Wu H., Deng Z. G., 2005, ApJ, 625, 78
- Hao et al. (2008) Hao C.-N., Xia X.-Y., Shu-DeMao Deng Z.-G., Wu H., 2008, Chinese Journal of Astronomy and Astrophysics, 8, 12
- Häring & Rix (2004) Häring N., Rix H.-W., 2004, ApJ, 604, L89
- Hopkins et al. (2008) Hopkins P. F., Hernquist L., Cox T. J., Kereš D., 2008, ApJS, 175, 356
- Houck et al. (1985) Houck J. R., Schneider D. P., Danielson G. E., Neugebauer G., Soifer B. T., Beichman C. A., Lonsdale C. J., 1985, ApJ, 290, L5
- Isaak et al. (2002) Isaak K. G., Priddey R. S., McMahon R. G., Omont A., Peroux C., Sharp R. G., Withington S., 2002, MNRAS, 329, 149
- Jiang et al. (2007) Jiang L., Fan X., Vestergaard M., Kurk J. D., Walter F., Kelly B. C., Strauss M. A., 2007, AJ, 134, 1150
- Juarez et al. (2009) Juarez Y., Maiolino R., Mujica R., Pedani M., Marinoni S., Nagao T., Marconi A., Oliva E., 2009, A&A, 494, L25
- Kauffmann & Haehnelt (2000) Kauffmann G., Haehnelt M., 2000, MNRAS, 311, 576
- Kennicutt (1998) Kennicutt Jr. R. C., 1998, ARA&A, 36, 189
- Li et al. (2008) Li C., Kauffmann G., Heckman T. M., White S. D. M., Jing Y. P., 2008, MNRAS, 385, 1915
- Lonsdale et al. (2006) Lonsdale C. J., Farrah D., Smith H. E., 2006, Ultraluminous Infrared Galaxies. Springer Verlag, pp 285–+
- Lutz et al. (2007) Lutz D., Sturm E., Tacconi L. J., Valiante E., Schweitzer M., Netzer H., Maiolino R., Andreani P., Shemmer O., Veilleux S., 2007, ApJ, 661, L25
- Lutz et al. (2008) Lutz D., Sturm E., Tacconi L. J., Valiante E., Schweitzer M., Netzer H., Maiolino R., Andreani P., Shemmer O., Veilleux S., 2008, ApJ, 684, 853
- Marconi & Hunt (2003) Marconi A., Hunt L. K., 2003, ApJ, 589, L21
- Matsuoka et al. (2009) Matsuoka K., Nagao T., Maiolino R., Marconi A., Taniguchi Y., 2009, ArXiv e-prints
- McMahon et al. (1999) McMahon R. G., Priddey R. S., Omont A., Snellen I., Withington S., 1999, MNRAS, 309, L1
- Menéndez-Delmestre et al. (2009) Menéndez-Delmestre K., Blain A. W., Smail I., Alexander D. M., Chapman S. C., Armus L., Frayer D., Ivison R. J., Teplitz H., 2009, ApJ, 699, 667
- Merritt & Ferrarese (2001) Merritt D., Ferrarese L., 2001, ApJ, 547, 140
- Misawa et al. (2003) Misawa T., Eracleous M., Charlton J., Ganguly R., Tytler D., Kirkman D., O’Meara J., Suzuki N., Lubin D., 2003, in Bulletin of the American Astronomical Society Vol. 35 of Bulletin of the American Astronomical Society, A Census of Intrinsic Narrow Absorption Lines at z 3.0. pp 1329–+
- Nagao et al. (2006) Nagao T., Maiolino R., Marconi A., 2006, A&A, 447, 863
- Nagao et al. (2006) Nagao T., Marconi A., Maiolino R., 2006, A&A, 447, 157
- Nestor et al. (2008) Nestor D., Hamann F., Rodriguez Hidalgo P., 2008, MNRAS, 386, 2055
- Netzer et al. (2007) Netzer H., Lutz D., Schweitzer M., Contursi A., Sturm E., Tacconi L. J., Veilleux S., Kim D.-C., Rupke D., Baker A. J., Dasyra K., Mazzarella J., Lord S., 2007, ApJ, 666, 806
- Omont et al. (2003) Omont A., Beelen A., Bertoldi F., Cox P., Carilli C. L., Priddey R. S., McMahon R. G., Isaak K. G., 2003, A&A, 398, 857
- Omont et al. (2001) Omont A., Cox P., Bertoldi F., McMahon R. G., Carilli C., Isaak K. G., 2001, A&A, 374, 371
- Pentericci et al. (2002) Pentericci L., Fan X., Rix H.-W., Strauss M. A., Narayanan V. K., Richards G. T., Schneider D. P., Krolik J., Heckman T., Brinkmann J., Lamb D. Q., Szokoly G. P., 2002, AJ, 123, 2151
- Peterson (2008) Peterson B. M., 2008, New Astronomy Review, 52, 240
- Peterson et al. (2004) Peterson B. M., Ferrarese L., Gilbert K. M., Kaspi S., Malkan M. A., Maoz D., Merritt D., Netzer H., Onken C. A., Pogge R. W., Vestergaard M., Wandel A., 2004, ApJ, 613, 682
- Priddey et al. (2003) Priddey R. S., Isaak K. G., McMahon R. G., Omont A., 2003, MNRAS, 339, 1183
- Priddey & McMahon (2001) Priddey R. S., McMahon R. G., 2001, MNRAS, 324, L17
- Richards (2001) Richards G. T., 2001, ApJS, 133, 53
- Richstone et al. (1998) Richstone D., Ajhar E. A., Bender R., Bower G., Dressler A., Faber S. M., Filippenko A. V., Gebhardt K., Green R., Ho L. C., Kormendy J., Lauer T. R., Magorrian J., Tremaine S., 1998, Nature, 395, A14+
- Rodríguez Hidalgo et al. (2010) Rodríguez Hidalgo P., Hamann F., Nestor D. B., Shields J., 2010, in prep.
- Rupke et al. (2005) Rupke D. S., Veilleux S., Sanders D. B., 2005, ApJ, 632, 751
- Sanders & Mirabel (1996) Sanders D. B., Mirabel I. F., 1996, ARA&A, 34, 749
- Sanders et al. (1988) Sanders D. B., Soifer B. T., Elias J. H., Neugebauer G., Matthews K., 1988, ApJ, 328, L35
- Schweitzer et al. (2006) Schweitzer M., Lutz D., Sturm E., Contursi A., Tacconi L. J., Lehnert M. D., Dasyra K. M., Genzel R., Veilleux S., Rupke D., Kim D.-C., Baker A. J., Netzer H., Sternberg A., Mazzarella J., Lord S., 2006, ApJ, 649, 79
- Serjeant & Hatziminaoglou (2009) Serjeant S., Hatziminaoglou E., 2009, ArXiv e-prints
- Shemmer et al. (2004) Shemmer O., Netzer H., Maiolino R., Oliva E., Croom S., Corbett E., di Fabrizio L., 2004, ApJ, 614, 547
- Shields (1976) Shields G. A., 1976, ApJ, 204, 330
- Shields et al. (2006) Shields G. A., Salviander S., Bonning E. W., 2006, New Astronomy Review, 50, 809
- Simon & Hamann (2010) Simon L. E., Hamann F., 2010, submitted to MNRAS
- Storrie-Lombardi et al. (1996) Storrie-Lombardi L. J., McMahon R. G., Irwin M. J., Hazard C., 1996, ApJ, 468, 121
- SubbaRao et al. (2002) SubbaRao M., Frieman J., Bernardi M., Loveday J., Nichol B., Castander F., Meiksin A., 2002, in Starck J.-L., Murtagh F. D., eds, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series Vol. 4847 of Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, The Sloan Digital Sky Survey 1-Dimensional Spectroscopic Pipeline. pp 452–460
- Tremaine et al. (2002) Tremaine S., Gebhardt K., Bender R., Bower G., Dressler A., Faber S. M., Filippenko A. V., Green R., Grillmair C., Ho L. C., Kormendy J., Lauer T. R., Magorrian J., Pinkney J., Richstone D., 2002, ApJ, 574, 740
- Trump et al. (2006) Trump J. R., Hall P. B., Reichard T. A., Richards G. T., Schneider D. P., Vanden Berk D. E., Knapp G. R., Anderson S. F., Fan X., Brinkman J., Kleinman S. J., Nitta A., 2006, ApJS, 165, 1
- Tytler & Fan (1992) Tytler D., Fan X.-M., 1992, ApJS, 79, 1
- Veilleux et al. (2009) Veilleux S., Rupke D. S. N., Kim D., Genzel R., Sturm E., Lutz D., Contursi A., Schweitzer M., Tacconi L. J., Netzer H., Sternberg A., Mihos J. C., Baker A. J., Mazzarella J. M., Lord S., Sanders D. B., Stockton A., Joseph R. D., Barnes J. E., 2009, ApJS, 182, 628
- Veron-Cetty & Veron (2000) Veron-Cetty M. P., Veron P., 2000, VizieR Online Data Catalog, 7215, 0
- Vestergaard (2003) Vestergaard M., 2003, ApJ, 599, 116
- Vestergaard (2004) Vestergaard M., 2004, ApJ, 601, 676
- Warner et al. (2003) Warner C., Hamann F., Dietrich M., 2003, ApJ, 596, 72
- Warner et al. (2004) Warner C., Hamann F., Dietrich M., 2004, ApJ, 608, 136
- Weymann et al. (1991) Weymann R. J., Morris S. L., Foltz C. B., Hewett P. C., 1991, ApJ, 373, 23
- Wild et al. (2008) Wild V., Kauffmann G., White S., York D., Lehnert M., Heckman T., Hall P. B., Khare P., Lundgren B., Schneider D. P., vanden Berk D., 2008, MNRAS, 388, 227
- Wyithe & Loeb (2003) Wyithe J. S. B., Loeb A., 2003, ApJ, 595, 614
- Young et al. (1982) Young P., Sargent W. L. W., Boksenberg A., 1982, ApJS, 48, 455