GBT Detection of Polarization-dependent HI Absorption and HI Outflows in Local ULIRGs and Quasars
We present the results of a 21-cm HI survey of 27 local massive gas-rich late-stage mergers and merger remnants with the Robert C. Byrd Green Bank Telescope (GBT). These remnants were selected from the Quasar/ULIRG Evolution Study (QUEST) sample of ultraluminous infrared galaxies (ULIRGs; L L) and quasars; our targets are all bolometrically dominated by active galactic nuclei (AGN) and sample the later phases of the proposed ULIRG-to-quasar evolutionary sequence. We find the prevalence of HI absorption (emission) to be 100% (29%) in ULIRGs with HI detections, 100% (88%) in FIR-strong quasars, and 63% (100%) in FIR-weak quasars. The absorption features are associated with powerful neutral outflows that change from being mainly driven by star formation in ULIRGs to being driven by the AGN in the quasars. These outflows have velocities that exceed 1500 km s in some cases. Unexpectedly, we find polarization-dependent HI absorption in 57% of our spectra (88% and 63% of the FIR-strong and FIR-weak quasars, respectively). We attribute this result to absorption of polarized continuum emission from these sources by foreground HI clouds. About 60% of the quasars displaying polarized spectra are radio-loud, far higher than the 10% observed in the general AGN population. This discrepancy suggests that radio jets play an important role in shaping the environments in these galaxies. These systems may represent a transition phase in the evolution of gas-rich mergers into “mature” radio galaxies.
Subject headings:galaxies: active — galaxies: evolution — quasars: absorption lines — ISM: jets and outflows — radio lines: galaxies
Recent theoretical work suggests that feedback from AGN plays an important role in establishing the present-day appearance of massive galaxies. Major galaxy mergers at high redshifts may trigger powerful starbursts, lead to the formation of elliptical galaxies, and account for the growth of supermassive black holes (e.g., Sanders et al., 1988; di Matteo et al., 2005; Springel et al., 2005; Hopkins et al., 2005, 2008). In this merger-driven evolutionary scenario, two gas-rich progenitors collide to form a completely obscured ULIRG where substantial inflows of gas trigger starbursts and feed the central black hole(s). As the system evolves, the obscuring gas and dust begin to disperse, bringing about the formation of a dusty IR-excess quasar and eventually a “naked” optical quasar. Powerful winds, driven by either the starburst or the central black hole, seem to be required to quench star formation and to stop the growth of the black hole in order to explain the observed relation (e.g., Murray et al., 2005).
There is growing support for the existence of this type of “negative feedback” in luminous late-stage mergers. Most galaxies with star formation rate (SFR) densities 0.1 yr kpc show signs of outflow (e.g., Heckman, 2002; Veilleux et al., 2005). These starburst-driven winds are often seen in ULIRGs, with mass outflow rates of 0.1–5 times the SFR and velocities 100–400 km s (Rupke et al., 2002, 2005a, 2005b, 2005c; Martin, 2005, 2006). In late-stage mergers and pure Seyfert galaxies, however, these winds are sometimes driven by the AGN, reaching outflow velocities 1000 km s (Rupke et al., 2005c; Krug et al., 2010). Similar AGN-driven outflows are also found in recently (re)started radio galaxies (e.g., Oosterloo et al., 2000; Morganti et al., 2005b). Several studies have also shown that these outflows are often multi-phased, with powerful molecular, neutral, and ionized outflows having all been detected in the same object (e.g., Morganti et al., 2007; Fischer et al., 2010; Feruglio et al., 2010; Rupke & Veilleux, 2011; Sturm et al., 2011).
In the present paper, we focus on the HI properties of the QUEST ULIRGs and quasars, searching for the presence and the driving source (AGN or starburst) of neutral outflows in merging systems as a function of merger stage. This paper is organized as follows. In Section 2, we discuss how our HI sample is selected. In Section 3, we describe our observational setup, data reduction, and spectral analysis. We then examine the HI profiles of our targets in Section 4 and the detection of polarized absorption (including a more sophisticated spectral decomposition that boosts our number of absorption detections) in Section 5. Finally, in Section 6, we discuss the possible origins of the HI outflows. We conclude in Section 7 with a summary of our findings. Throughout this paper, we adopt = 71 km s Mpc, = 0.27, and = 0.73 (Hinshaw et al., 2009).
2. Sample Selection
In order to assess the physical conditions and the dynamical evolution of the neutral gas in late-stage merger systems, we selected targets from the QUEST sample of ULIRGs and quasars that are bolometrically dominated by AGN. The details of the QUEST sample selection are outlined in previous QUEST papers (e.g., Schweitzer et al., 2006; Netzer et al., 2007; Veilleux et al., 2009a). Briefly, the original QUEST sample of ULIRGs begins with the 1-Jy sample, a flux-limited sample of 118 ULIRGs selected from a redshift survey of the IRAS faint source catalog (Kim & Sanders, 1998). To these, five well-studied local ULIRGs from the Revised Bright Galaxy Survey (RBGS; Sanders et al., 2003) were added given the availability of archival data from multiple observatories (see Veilleux et al., 2009a; Veilleux, 2012, for details). The original quasar sample contains 25 quasars, including 24 Palomar-Green (PG) quasars from the Bright Quasar Sample (Schmidt & Green, 1983) and B2 2201+31A which satisfies the magnitude criterion of Schmidt & Green (1983). An additional nine PG quasars were added to the sample from archival data. The selected quasars fall at the low end of the quasar -band luminosity function.
The ULIRGs and quasars in the QUEST sample are well-matched in redshift () and in near-IR imaging show tidal features that indicate merger origins (Veilleux et al., 2006, 2009b). The entire QUEST sample spans the merger evolution sequence from binary starburst-dominated ULIRGs, to completely coalesced AGN-dominated ULIRGs, to dust-enshrouded quasars, to optical quasars whose dust has been cleared.
Our HI sample is a subset of the QUEST targets. To ensure quality spectra, we require ; our 27 targets evenly sample coalesced ULIRGs111This category includes Mrk 273, the only source in our sample to have binary nuclei (see Appendix). (9 targets), far-infrared(FIR)-strong quasars (; 10 targets), and FIR-weak quasars (8 targets). The two classes of quasars have similar spectral energy distributions representative of AGN, with the FIR-strong quasars having an additional contribution to their FIR emission from intense star formation (Netzer et al., 2007). The morphologies of the quasar hosts are in agreement with the picture that FIR-weak quasars are at a later evolutionary stage than FIR-strong quasars, following the (near) cessation of star formation (Veilleux et al., 2009b). Thus, our targets sample the later stages of merger evolution from coalescence to the aftermath of significant star formation. Table 1 summarizes the properties of our HI sample.
3. Observations and Data
3.1. Observational Setup
The targets were observed with the GBT Spectrometer in the L band between 2011 August and 2012 March during the 160 hours of time allocated for GBT programs 11B075 and 12A075 (PI: Teng). The data were collected in 10 minute on-off source pairs (5 minutes per position) with 12.5 MHz bandwidth, two spectral windows centered at the same frequency (1420.4058 MHz), two linear polarization modes (XX and YY), and nine sampling levels. To minimize the effects of radio frequency interference (RFI) on the data, we used the notch filter for sources with heliocentric recessional velocities below 17000 km s. The total on-source time for each galaxy was determined during the observing runs based on the strength of the spectral features relative to the noise or up to a maximum of three hours. Most sources were observed over multiple observing sessions. The circumpolar starburst galaxy M82 was observed at the beginning of each observing session as a check for instrumental variations.
We also measured the continuum flux densities of our sample at 1.42 GHz using the AutoPeak procedure in order to check for time variability. However, due to the plethora of continuum sources within the large GBT beam, the 5- confusion limit is 100 mJy at this frequency. Therefore, we were only able to obtain confident continuum measurements for five targets. These values are in agreement with older values obtained with the Very Large Array (VLA) for their nuclear core flux densities (Condon et al., 2002), despite the much larger beam of the GBT. For the remainder of sample, we use published values of the core continuum flux densities as measured with the VLA (e.g., Barvainis et al., 1996; Condon et al., 1998; Rafter et al., 2009).
3.2. Data Reduction
The spectra were reduced using GBTIDL (Marganian et al., 2006). Individual five-second records showing large broadband harmonic RFI were flagged, and persistent RFI spikes that appear in the data were replaced by interpolation across the affected channels for both spectral windows and polarizations. For each polarization, the data from the two spectral windows were averaged together, and a polynomial of order 5 was fitted over a range of 1000 km s on each side of a galaxy’s systemic redshift in the Hanning-smoothed and then decimated222The total number of channels in the spectrum is decreased, or decimated, by a factor equal to the order of smoothing performed. spectrum to subtract the baseline. The same baseline-fitting regions and order of polynomial fit were used for both polarizations.
Flux calibration was derived from simultaneous observations of an internal noise diode. The baseline-subtracted and flux-calibrated scans for each polarization were accumulated and averaged for each galaxy. A fourth order boxcar smoothing was then applied to the averaged data with decimation, resulting in a resolution of 6 km s per channel. The two polarizations were then averaged together to produce the full intensity spectra shown in Figure 1. We ensured that our flux calibration is correct by comparing the full intensity profile of Mrk 231 with that measured by the VLA (Carilli et al., 1998); both profiles are identical. Discussions of individual objects in the Appendix include comparisons with previous data, whose generally good agreement with our data in terms of velocity profile and flux validate our data reduction process. Table 1 also provides a list of some observational parameters for our galaxies.
Nearly all of the galaxies were detected. Two galaxies, PG 0804+761 and F15462–0450, had strong RFI spikes overlapping with the expected HI frequencies. At the limit of three hours on source, PG 1126–041 and F21219–1757 were not detected, and F04103–2838 and F15130–1958 were marginally detected. The detections of PG 1426+015 and PG 1617+175 are uncertain because their putative emission features have widths similar to those of broadness of the baseline ripples before baseline subtraction.
3.3. Spectral Analysis
The individual HI profiles in the spectra of our sample galaxy are complex combinations of absorption and emission. For simplicity, we have measured the emission and absorption features separately, but we note that this approach may slightly alter the final measured values of these quantities. More complex decomposition of the emission and absorption is only possible for a few sources because of the degeneracy between emission and absorption, as discussed in Section 5.2.
For galaxies with significant emission (peak-to-rms ), we measure the central velocity (), width (, ), and total velocity-integrated flux () using the GBTIDL task gmeasure. In these instances, since there are no obvious companions enclosed within the GBT beam, the total flux is integrated between velocities at which the flux density reaches zero. The central velocity is the heliocentric velocity measured at the midpoint of the 20% flux level. The width of emission is measured at the 20% and 50% levels ( and , respectively). The measurement errors in these quantities are derived using the expressions presented in Wei et al. (2010) following the methodologies of Schneider et al. (1986, 1990) and Fouque et al. (1990). The HI column density is calculated using expression 12.17 from Rohlfs & Wilson (2000) for optically thin radiation. Assuming the GBT antenna temperature to flux density gain conversion of 2 K Jy:
where is the velocity integrated HI flux in units of Jy km s. The neutral gas mass is calculated following Haynes & Giovanelli (1984):
where is the recessional velocity corrected for the influence of the Virgo Cluster, the Great Attractor, and the Shapley Supercluster as calculated by the NASA Extragalactic Database (NED). All of these measured and derived quantities are listed in Table 2.
Only a few absorption features are significantly detected (by the criterion that peak-to-rms ). However, many of the absorption features are very broad but shallow, so the peak-to-rms ratio is not necessarily the best determinant of a detection. We also consider a feature to be detected if its full width at half maximum (FWHM) spans at least 10 velocity channels and it has a peak-to-rms . These absorption features are quantified using the GBTIDL task fitgauss, with one Gaussian component at a time added to a spectrum until the residuals are statistically featureless. With the exception of PG 1119+120, where a narrow absorption feature is present within the emission profile, only the absorption features are modeled. For PG 1119+120, the emission was first modeled with Gaussians and subtracted from the data to define a baseline. The absorption was then modeled from the emission-subtracted spectrum. The Gaussian models are plotted (as components and as sums) in Figure 1. The GBTIDL task measures the central velocity (), FWHM, the “peak” flux density, and their measurement errors for each Gaussian component. The full velocity-integrated absorbed flux (), optical depth (), and column density () are then calculated from the best-fit models. We assume the HI absorbing screen covers the nuclear continuum source, in which case the optical depth profile is:
where is the HI line profile as a function of velocity and is the continuum flux density (Whiteoak & Gardner, 1972). The column density is then calculated as:
where is the spin temperature. In this paper, we will assume a spin temperature of 1000 K (see Section 6.1). For absorption features in PG 1119+120, PG 1244+026, and PG 1440+356, the absorption troughs are deeper than their nuclear continua, suggesting that sources of continuum emission extend beyond the nucleus. Table 3 summarizes the measured and derived values for all of the detected absorption features.
As a reality check, we inspected the data to confirm that the choice of polynomial order for the baseline fitting does not affect the detection of absorption or emission. Specifically, we examined whether polynomial orders of 3 or 5 can introduce ripples into the data that masquerade as absorption or emission features. After repeating the baseline subtraction for our data using a range of polynomial orders (1–5), we confirm that the same spectral features are still detected, even for marginal detections such as that in F04103–2838. The measured widths and velocities are in agreement with those measured with the best-fit baseline to within the errors. Therefore, we are confident that even the weak detections noted as ”marginal” in Tables 2 and 3 are real.
4. HI Profiles of Late-stage Mergers
The individual HI profiles of our targets are discussed in the Appendix. In general, we see HI trends along the merger sequence but no correlation with infrared luminosity.
All of the detected ULIRGs show HI absorption, typically with multiple velocity components, many of which are narrow (FWHM km s) and some of which are redshifted. Only 2/7 (29%) of the detected ULIRGs have HI emission. The different velocity absorption components have different origins and are not necessarily caused by outflows or infall; kinematics of circumnuclear gas disks may mimic the properties of outflows or infall in single dish data. In a couple of cases, such as UGC 05101 and Mrk 273, the absorption profiles are nearly symmetrical about zero velocity. For the Type 2 AGN UGC 05101, the profile is likely from absorption due to a slightly inclined rotating HI disk superposed on an obscured nucleus. This interpretation is consistent with CO maps of UGC 05101 showing smooth, regular rotation of the molecular gas (Genzel et al., 1998; Wilson et al., 2008). On the other hand, the CO velocity field of Mrk 273 is very distorted (Downes & Solomon, 1998; Wilson et al., 2008), so it is unlikely that the gas is in well-ordered disk-like rotation. The double absorption in HI seen in Mrk 273 probably corresponds to the binary obscured nuclei in this object. Thus, CO maps can help eliminate gas disks as the origin of the absorption features. However, only a few galaxies in our sample have readily available CO maps so direct comparisons are difficult. Comparison of single dish CO and HI profiles for the QUEST targets is on-going, but it is beyond the scope of the present paper. The more complex profiles, such as those seen in F15250+3608, F05189–2524, and F15130–1938, may be interpreted as multiple infalling or outflowing components. They may also result from the superposition of an emission profile and a broad absorption profile. The latter case is discussed in more detail in Section 5.2.
In contrast, the HI profiles of the FIR-strong quasars are generally detected not only in absorption (88% or 7 of the 8 detected sources) but also in emission (75%). The absorption features are on average 60 km s broader than those seen in ULIRGs or approximately three times greater than the typical FWHM measurement error in the ULIRGs. Both redshifted and blueshifted absorption features are seen, although the systemic velocities of these systems are uncertain. Optical emission lines may underestimate recessional velocity since some fraction of the line emission may arise from outflowing material. In the case of PG 0050+124 (I Zw 1), this effect causes a shift of 700 km s (see Appendix). A similar situation could also hold for PG 1440+356. In three cases, PG 1119+120, PG 1244+026, and PG 1440+356, the minima of the absorption troughs exceed the nuclear continuum flux densities. This discrepancy implies that off-nuclear emission, not accounted for in the nuclear fluxes of Table 1, also contributes to the continuum emission in these four sources. Given resolution limit of the VLA core measurements of PG 1119+120 and PG 1440+356, the off-nuclear emission must be at least 20–30 kpc from the core; for PG 1244+026, the separation from the core must be at least 3 kpc. While we cannot formally rule out the possibility that these absorption features are absorption of HI against continuum from background quasars, this possibility is statistically unlikely. Another possibility is that the flux densities of the cores in these galaxies are time variable as exemplified by the variability of the core in PG 1351+640 (see Appendix).
The HI profiles of the FIR-weak quasars are different from those of the ULIRGs and FIR-strong quasars. HI is seen mostly in emission (88% or 7 of 8 detected sources), and only 25% (2/8) show broad absorption. The emission profiles are very broad (500 km s) and bumpy, perhaps beginning to resemble those of regularly rotating disks that typically display Gaussian or double-horned profiles; the measured widths and velocities are consistent with previous measurements of quasars (e.g., Bothun et al., 1984; Condon et al., 1985; Hutchings et al., 1987; Ho et al., 2008; König et al., 2009). These properties will be revisited in Section 5.2.
Considering now our entire HI sample, the distribution of equivalent widths (EW; defined as the integrated flux density for each component divided by the continuum flux density of the core) is plotted in Figure 2. For this definition, the equivalent widths of the absorption components are negative. It is clear that the strengths of the emission and absorption features progress along the merger sequence, with ULIRGs having the weakest features, the FIR-strong quasars having the strongest absorption features, and the FIR-weak quasars having the strongest emission features. In Figure 3, we compare the FWHMs of emission and absorption with the star formation rates of their host galaxies. For emission, the FIR-weak sources have the broadest widths, the ULIRGs have the narrowest widths, and the FIR-strong quasars bridge the gap between those two classes. In the ULIRGs, there appears to be a weak correlation between the width of the absorption components and the SFR, reminiscent of a similar relationship seen in Na I (Rupke et al., 2005b). However, given that multiple absorption components are associated with a single system (and thus a single SFR value), the validity of the connection between the SFR and the FWHM of the absorption components in ULIRGs is uncertain. On the other hand, absorption FWHM is not at all related to SFR in the quasars, further suggesting that the AGN play a role in these systems.
5. Polarization-dependent HI Absorption
An unexpected outcome of our GBT program is that many of the absorption features show a dependence on polarization (Figure 4). Indeed, 13 of our 23 GBT-detected targets show some profile difference in their two polarizations (1 ULIRG, 7 IR-excess quasars, and 5 optical quasars). For GBT observations, it is expected that there is a 6% difference between polarizations in terms of normalization, but the observed spectral shapes should be the same in both polarizations (R. Maddalena, private comm.). This expectation is confirmed by the fact that our M82 data remain in agreement after baseline subtraction. Polarization variation in our QUEST sample is not RFI-related, as we have removed time intervals where the data are affected by broadband RFI. Since it is improbable that the HI screen itself is polarized, the underlying continuum emission must be polarized.
5.1. Time Variability
The most dramatic examples of polarized spectra are seen in PG 1211+143 and PG 1440+356. Not only are these spectra polarized, but they are also time variable on timescales of 30 minutes (Figures 5 and 6). In fact, of our radio-detected sample galaxies, time variability is only seen in sources that have significant absorption and spectra that are significantly polarized. A possible concern is that since we observed these sources over several sessions, variability could be an instrumental effect. However, comparing the different sessions’ data for M82, we find that time variability and polarization are not seen beyond what is expected of typical GBT observations. Moreover, session-to-session instrumental variability would not explain the 30-minute variability seen within the same observing session. Time varying polarizations are not seen only in the targets with the brightest continua, so the variations are not likely to be due to instrumental systematics from poor cancellation of baseline ripples in the position-switching scans (O’Neil, 2002).
A different sort of instrumental effect may be leakage of linear polarization from the linearly polarized background continuum into the total intensity of the absorption profiles. We cannot yet exclude this scenario since these observations were not designed for polarimetry and our current data lack the cross-correlated products (XY and YX) that would be needed to confirm these time variations.
If these time variations are real, they could be caused by one of several processes. First, they could be due to variations in the compact continuum sources behind the HI screen. As the X-ray light curves of PG 1211+143 and PG 1440+356 show (Figure 7), the X-ray emitting central region can be highly variable on a timescale of about 30 minutes. Although not contemporaneous with our GBT data, these X-ray data suggest that the AGN continuum can vary on a 30-minute time scale if the radio continuum is indeed related to the X-ray producing regions. However, if the polarization-dependent variations are simply a result of the changing compact source continuum, then we would expect the continuum to vary in block and would observe a change in the equivalent width of the HI absorption, not the frequency-dependent variations we see in these targets. Second, variability could be a result of changes in the structure of the HI screen. However, the timescale for the destruction of small HI clouds in the interstellar medium is on the order of a Myr, depending on the exact properties of the interstellar medium causing the absorption (e.g., McKee & Cowie, 1977; Nagashima et al., 2006); it is therefore unlikely that cloud destruction causes the observed short-term variations. Similarly, it is unlikely that the movement of HI clouds into and out of our line of sight to the continuum source can cause the very broad variations observed in these objects.
The time variability of the polarization-dependent HI absorption is easier to explain if the polarized continuum emission comes from both the nuclear source and an extended source (e.g., radio jets). In this picture (Figure 8), the HI screen is also extended, covering both the nuclear source and extended jets with kinematically distinct clouds. As the central source fluctuates, the ratio of central-to-jet emission also changes, resulting in the observed wavelength-dependent polarization variations. In this scenario, there is no need for the structure of the HI screen to change on a fast timescale.
5.2. Decomposing the Spectrum
Regardless of whether the time variability is real, it is certainly the case that the spectra shown in Figure 4 are polarized. Again, we checked for baseline-subtraction-related effects on the detection of polarized absorption and found that the choice of the polynomial order of the baseline fit matters little to the detection of these polarization-dependent absorption features. The most obvious example of polarization not due to observational systematics is seen in the spectra of PG 0050+124 (I Zw 1). This HI profile is not polarized at relative velocities of 500–1000 km s but becomes polarization-dependent between relative velocities –500 and +500 km s (Figure 4). This behavior argues against attribution of apparent polarization dependencies to poor baseline subtraction or session-to-session variations in our GBT observations. At least in this one case, the polarization dependence must be due to something physically happening at relative velocities between –500 and +500 km s. When we compare the HI profiles of each galaxy in the two separate polarizations, it appears that one polarization typically has higher flux than the other, with the exception of PG 1440+356, in which there is a velocity shift instead. The effect is not systematic in the sense that one polarization is always higher than the other; each is nearly equally likely to be the polarization with the higher intensity. In some cases, particularly in PG 1351+640 and PG 1411+442, the lower-intensity profile appears to be the same as the higher-intensity profile with the addition of a broad absorption component.
To follow up on this intriguing hint, we produced the difference spectra (the spectrum in the polarization with the more negative HI profile minus the spectrum in the polarization with the more positive HI profile) shown in Figure 9. The absorption features in these spectra are then modeled the same way as before, using the GBTIDL task fitgauss. When these absorption models are included, many of the full intensity spectra that need multiple components to explain the blue- and red-shifted components can be explained simply by HI emission that is combined with a single broad absorption component. This pattern is most obvious in PG 1351+640 and PG 1411+442. The results of this type of complex decomposition of the spectra are provided in Tables 4 and 5. Similar decompositions can explain the shapes of the full intensity spectra of the ULIRGs F05189–2524 and F15130–1958, and the FIR-strong quasar PG 1244+026. However, since those spectra are not polarized, we cannot create difference spectra to model the broad absorption component. The fitgauss task statistically prefers the multiple-component models to the single-absorption model for these spectra. Overall, with the exception of PG 1501+106, the decompositions result in stronger features in both absorption and emission.
Given the results of this exercise, we find the instances of HI absorption (emission) to be 100% (29%) in detected ULIRGs, 100% (88%) in FIR-strong quasars, and 63% (100%) in FIR-weak quasars. The more complex decomposition of the spectra does not drastically change the distribution of EWs (Figure 10) compared to Figure 2. Figure 11 provides a comparison of the FWHMs of the emission and absorption features from the complex modeling with the star formation rates of these galaxies. As in Figure 3, the FIR-weak quasars have the broadest emission features and the ULIRGs have the narrowest emission features. In terms of absorption, there is no obvious relationship between the absorption feature widths and the star formation rates. As before, the FIR-weak quasars have the broadest absorption features, with the FIR-strong quasars spanning the range of widths between FIR-weak quasars and ULIRGs.
6. The Properties of HI Outflows
As discussed in Section 4, the absorption components are not necessarily caused by infall or outflow. It is possible, as seen in UGC 05101 (edge on) and Mrk 231 (face on), that these features can be associated with the kinematics of gas disks. A circumnuclear gas disk seen nearly edge on with a non-isotropic gas mass (or background continuum) distribution can mimic characteristics of infall and outflow. Unfortunately, our single dish data cannot unambiguously exclude gas disks as the source of the features we observe, but it is unlikely that all of our targets with absorption features have gas disks that are seen precisely edge-on.
Keeping this caveat in mind, we assume in the following discussion an outflow origin only for blueshifted HI absorption features with HI-corrected heliocentric velocities (333These systemic velocities are relative to the HI emission rather than the optical emission. See Table 6 for more details. A remaining source of uncertainty in the interpretation of HI line profiles is the question of whether the midpoint of the 20% flux level is always a good approximation of the true systemic velocity. Emission from molecular lines is a more reliable tracer of systemic velocity than HI, so the ideal solution here would be acquisition of single-dish CO spectroscopy for the full QUEST sample. Unfortunately, these observations are still pending.) below –50 km s (Rupke et al., 2005a). A number of mechanisms have been proposed to explain such HI outflows. They include adiabatically expanding broad emission line clouds (e.g., Elvis et al., 2002), starburst- and AGN-driven galactic winds (e.g., Heckman et al., 1990; Veilleux et al., 2005), and AGN-driven jets interacting with the surrounding interstellar medium (ISM; e.g., Oosterloo et al., 2000). These mechanisms predict different outflow locations (from parsec to kiloparsec scales) and geometries (wide angled vs. collimated).
6.1. Mass Outflow Rates
assuming spherical geometry, mass conservation, and a wind that flows from a distance of to infinity. The outflow is also assumed to subtend a solid angle of within which the gas has a covering fraction of 1. To facilitate comparisons with the values derived from young radio galaxies by Morganti et al. (2005b), we use the same assumptions of and kpc. Applying these assumptions and Equation 5 to the blueshifted absorption components where , we find mass outflow rates (wide) up to more than 6000 yr (summarized in Table 6). These values are much larger than those of 60 yr measured by Morganti et al. (2005b) and the range measured in ULIRGs (13–133 yr; Rupke et al., 2002).
Many of the mass outflow rates we infer seem unrealistically large; such high rates have not been measured in other astrophysical objects. These values do, however, depend on the assumptions we have made above. The mass outflow rate depends linearly on the size scale of the outflow region and the derived HI column density. Recall from Equation 4 that the column density depends on the spin temperature and optical depth profile, which in turn depends on the brightness of the background continuum. There are thus three key factors that affect our estimates of a mass outflow rate: the spin temperature of the gas, the assumed distance, and the background continuum. The assumption of 1000 K for the HI spin temperature seems to be reasonably conservative. In the vicinity of powerful AGN, the extreme conditions of the gas and the presence of shocks can result in spin temperatures 3000 K (e.g., Holt et al., 2006), although we note that absorption-line opacity will be dominated by the coldest material in any multiphase medium. Spatially resolved interferometric maps of several quasars suggest that the outflow regions have size scales 0.2–1.6 kpc (Oosterloo et al., 2000; Morganti et al., 2003; Emonts et al., 2005; Morganti et al., 2005a). Adopting the smaller distance of 0.2 kpc instead of 1 kpc would reduce the derived mass outflow rates by a factor of 5, more consistent with those derived for ULIRGs.
Finally, the estimated HI column density also depends strongly on the geometry of the source of continuum emission. We have assumed that this emission comes entirely from the core. However, the absorption troughs in several quasars in our sample are deeper than the core fluxes. This excess implies that off-nuclear emission also contributes to the observed continuum. In this case, in Equation (3) is overestimated and so is , by a factor equal to the actual nuclear/total continuum flux ratio, although a larger source size may compensate to some extent.
6.2. Jet-induced Outflows
If the off-nuclear continuum source is indeed physically related to the target and not a background quasar that happens to lie within the GBT beam, our polarized spectra give more credence to the idea that radio jets play an important role in driving HI gas outflows (e.g., Morganti, 2011). Several sources (see Appendix) in our sample with polarized spectra have spatially resolved radio jets detected at higher frequencies. The resolutions of these VLA data imply that the jets are launched within 0.4–0.9 kpc of the core. Considering a jet-induced outflow scenario for the HI absorption features, we again apply Equation 5 to the blueshifted outflow components and assume kpc and a narrow opening angle of 5 for the jet. Under these new assumptions, the range of mass outflow rates is 0.2–31 yr (see col. (8) of Table 6). We consider these values lower bounds on the actual mass outflow rates, while those derived using and kpc (col. (7) of Table 6) are upper bounds.
Under the jet scenario, the kinetic energy release rate can be calculated using the equation from Heckman et al. (2000):
The kinetic energy in the outflow is on the order of to erg s for our sample.
The association of jets with the detected HI absorption features is also seen in the distribution of radio loudness of our sample. Only one quasar in our sample is considered radio-loud based on optical-to-radio flux ratios (see Appendix), but optical data may be subject to contamination from the host galaxy, resulting in an underestimate of radio loudness. Using the 1.4 GHz to 2–10 keV ratio as an indicator of radio loudness and a radio-loud cutoff value empirically derived from 1600 AGN (La Franca et al., 2010), we find that many objects in our sample are radio-loud (Table 1). These results are consistent with those found in a recent study by Chandola et al. (2012) that suggest higher detection rates of HI absorption toward extended compact steep-spectrum and gigahertz-peaked spectrum sources than toward the cores of radio galaxies. The environments in which the jets interact with the interstellar medium can be complex, and we may be seeing a strongly polarized continuum being absorbed selectively by multiple neutral hydrogen clouds, resulting in polarized HI absorption.
The fraction of radio-loud sources with associated HI absorption (58%) is far higher than the 10% observed in the general population of AGN. Given the high star formation rates of the ULIRGs, the circumnuclear starburst can contaminate the 1.4 GHz measurements. Using the 1.4 GHz to SFR correlation of Bell (2003), we corrected the 1.4 GHz flux for star formation. Of the three ULIRGs with blueshifted absorption and X-ray detections, all (100%) are radio-loud. Similarly, of the six FIR-strong quasars with detected HI outflows and X-ray emission, three (50%) are considered radio-loud; one of the two (50%) FIR-weak quasars with HI outflows and X-ray flux are radio-loud. This pattern recalls an idea put forward by Wilson & Colbert (1995) on the formation of radio galaxies to explain why radio-quiet AGN outnumber radio-loud AGN. The authors suggested that major mergers, such as the progenitors of our sample of galaxies, result in coalesced supermassive black holes with high spin values due to the conservation of angular momentum. Substantial energy is then extracted from these fast-spinning black holes to drive jets, yielding radio-loud AGN. Thus, younger radio galaxies may have more powerful jets and therefore stronger jet-induced outflows. This scenario is consistent with the observations of Morganti et al. (2005b), whose sample of young or recently restarted radio galaxies show similar HI outflows. The Wilson & Colbert (1995) picture is also consistent with theoretical models of Tchekhovskoy et al. (2010) that suggest the radio loud/quiet dichotomy in AGN is due to the spin of the central black hole.
6.3. HI Outflows Along the Merger Sequence
Even though our sample with mass outflow rates is relatively small, we see some trends in our outflow measurements (Table 6 and Figure 12). However, given the small number of detected HI outflows, particularly with only two measurements for the FIR-weak quasars, we cannot statically confirm that the distributions we see in Figure 12 is significant. There is also no obvious relationships between mass outflow rate and the bolometric luminosity, AGN luminosity, or radio loudness for our sample.
Systems with deep absorption and outflows become more prevalent as we progress along the merger sequence from the newly coalesced ULIRGs to the FIR-strong quasars. The trend between outflow kinematics and star formation rates among the ULIRGs (Figures 3 and 11) suggests that the HI outflows in these objects are driven mainly by the starbursts, with a median mass outflow rate in the jet scenario ( (jet) 1 yr) similar to their SFRs and previous Na I-based outflow measurements (Rupke et al., 2002, 2005a, 2005b, 2005c; Martin, 2005, 2006). ULIRGs also tend to have the narrowest absorption features. In the FIR-strong quasar systems, the median mass outflow rate is (jet) 5.8 yr, considerably larger than their median SFR of 1.0 yr. While the median mass outflow rate of the FIR-weak quasars is lower, (jet) 0.3 yr, so is their median SFR (0.6 yr). Thus, the mass-loading factor, i.e., the ratio of mass outflow rate to SFR, is much greater in the FIR-strong than the value of near unity derived in ULIRGs. The same value is much lower in the FIR-weak quasars, but this is difficult to compare given that this is based on a single measurement. However, the distinction between the ULIRGs and FIR-strong quasars further supports the idea that the AGN is playing a role in driving HI outflows in the quasars before the end of large-scale star formation.
HI in emission is much more prevalent among the quasars than the ULIRGs, and the HI emission profiles increasingly resemble those of regularly-rotating neutral gas disks as we move from FIR-strong to FIR-weak quasars. The amount of HI seen in emission increases modestly from 410 among the FIR-strong quasars to 710 among the FIR-weak quasars. These masses are consistent with those measured in low-redshift quasars (König et al., 2009) and the high end of the gas mass range seen in nearby active galaxies (Ho et al., 2008). We may be witnessing the creation of HI disks from residual neutral gas at the centers of merger remnants.
We have analyzed of GBT HI spectra of 27 ULIRGs and quasars selected from the QUEST sample. We find that:
As mergers evolve from fully coalesced ULIRGs to FIR-strong quasars to FIR-weak quasars, the characteristics of their HI profiles also change from mostly absorption in the ULIRGs to mostly emission in the FIR-weak quasars, with the FIR-bright quasars having intermediate properties. The kinematics of these systems change from narrow absorption features likely due to outflows driven by starbursts in ULIRGs, to broad absorption features from jet-driven outflows driven by AGN in FIR-strong quasars, and finally to the emission profiles more characteristic of gas in rotation in FIR-weak quasars.
A large fraction of the quasars exhibit polarization-dependent absorption features, with 78% of FIR-strong quasars and 63% of FIR-weak quasars showing these features compared with 13% in ULIRGs. Some of these features vary on short time scales (down to 30 minutes). Neither time variability nor polarization-dependent spectral features are seen in contemporaneous observations of M82 (as expected), reducing the likelihood that these effects are instrumental artifacts. These features are likely due to HI absorption of polarized nuclear and off-nuclear continuum emission. The large fraction (60%) of radio-loud sources with HI outflow-associated absorption suggests a close relationship between radio jets and the HI screen.
The polarization-dependent absorption features are most likely due to the emergence of radio jets in these systems, contributing to the continuum emission against which the absorption features are detected. The 60% radio-loud fraction of our sources with HI outflows is far higher than the 10% seen in the general population of AGN. This distinction suggests that radio jets play an important role in shaping the environments in these systems, which may tract a transition phase in the evolution of gas-rich mergers into “mature” radio galaxies.
- Barvainis et al. (1996) Barvainis, R., Lonsdale, C., & Antonucci, R. 1996, AJ, 111, 1431
- Bell (2003) Bell, E.F. 2003, ApJ, 586, 794
- Blundell & Beasley (1998) Blundell, K.M. & Beasley, A.J. 1998, MNRAS, 299, 165
- Bothun et al. (1984) Bothun, G.D., Heckman, T.M., Schommer, R.A., & Balick, B. 1984, AJ, 89, 1293
- Carilli et al. (1998) Carilli, C.L., Wrobel, J.M., & Ulvestad, J.S. 1998, AJ, 115, 928
- Chandola et al. (2012) Chandola, Y., Gupta, N., & Saika, D.J. 2012, MNRAS, in press (arXiv: 1211.6852)
- Condon et al. (1985) Condon, J.J., Hutchings, J.B., & Gower, A.C. 1985, AJ, 90, 1642
- Condon et al. (1998) Condon, J.J. et al. 1998, AJ, 115, 1693
- Condon et al. (2002) Condon, J.J., Cotton, W.D., & Broderick, J.J. AJ, 2002, 124, 675
- di Matteo et al. (2005) Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 604
- Downes & Solomon (1998) Downes, D. & Solomon, P.M. 1998, ApJ, 507, 615
- Emonts et al. (2005) Emonts, B.H.C. et al. 2005, MNRAS, 362, 931
- Elvis et al. (2002) Elvis, M., Marengo, M. & Karovska, M. 2002, ApJ, 567, L107
- Feruglio et al. (2010) Feruglio et al. A&A, 518, L155
- Fischer et al. (2010) Fischer et al. A&A, 518, L41
- Fouque et al. (1990) Fouque, P., Durand, N., Bottinelli, L., Gouguenheim, L., & Paturel, G. 1990, A&AS, 86, 473
- Genzel et al. (1998) Genzel, R. et al. 1998, ApJ, 498, 579
- Haynes & Giovanelli (1984) Haynes, M.P. & Giovanelli, R. 1984, AJ, 89, 758
- Heckman et al. (1990) Heckman, T.M., Armus, L. & Miley, G.K. 1990, ApJS, 74, 833
- Heckman et al. (2000) Heckman, T.M., Lehnert, M.D., Strickland, D.K., & Armus, L. 2000, ApJS, 129, 493
- Heckman (2002) Heckman, T.M. 2002, in ASP Conf. Proc. 254, Extragalactic Gas at Low Redshift, ed. J. S. Mulchaey & J. Stocke (San Francisco, CA: ASP), 292
- Hinshaw et al. (2009) Hinshaw, G. et al. 2009, ApJS, 180, 225
- Ho et al. (2008) Ho, L.C., Darling, J., & Greene, J.E. 2008, ApJS, 177, 103
- Ho & Kim (2009) Ho, L.C. & Kim, M. 2009, ApJS, 184, 398
- Holt et al. (2006) Holt, J. et al. 2006, MNRAS, 370, 1633
- Hopkins et al. (2005) Hopkins, P.F. et al. 2005, ApJ, 630, 705
- Hopkins et al. (2008) Hopkins, P.F. et al. 2008, ApJS, 175, 356
- Hutchings et al. (1987) Hutchings, J.B., Gower, A.C., & Price, R. 1987, AJ, 93, 6
- Iwasawa et al. (2011) Iwasawa, K. et al. 2011, A&A, 428, 137
- Kennicutt (1998) Kennicutt, R. C. 1998, ARA&A, 36, 189
- Kim & Sanders (1998) Kim, D.-C. & Sanders, D. B. 1998, ApJS, 119, 41
- König et al. (2009) König, S., Eckart, A., García-Marín, M., & Huchtmeier, W. K., 2009, A&A, 507, 757
- Krug et al. (2010) Krug, H.B. et al. 2010, ApJ, 708, 1145
- La Franca et al. (2010) La Franca, F., Melini, G., & Fiore, F. 2010, ApJ, 718, 368
- Leipski et al. (2006) Leipski, C., Falcke, H., Bennert, N., & Huttemeister, S. 2006, A&A, 161, 172
- Marganian et al. (2006) Marganian, P., Garwood, R.W., Braatz, J.A., Radziwill, N.M., & Maddalena, R.J. 2006, Astronomical Data Analysis Software and Systems XV (ASP Conf. Ser. 351), ed. C. Gabriel, C. Arviset, D. Ponz, & S. Enrique (San Francisco, CA: ASP), 512
- Martin (2005) Martin, C.L. 2005, ApJ, 621, 227
- Martin (2006) Martin, C.L. 2006, ApJ, 647, 222
- McKee & Cowie (1977) McKee, C.F. & Cowie, L.L. 1997, ApJ, 215, 213
- Morganti et al. (2003) Morganti, R., Oosterloo, T. A., Emonts, B. H. C., van der Hulst, J. M., & Tadhunter, C. N. 2003, ApJ, 593, L69
- Morganti et al. (2004) Morganti et al. 2004, A&A, 424, 119
- Morganti et al. (2005a) Morganti et al. 2005a, A&A, 439, 521
- Morganti et al. (2005b) Morganti et al. 2005b, A&A, 444, 9
- Morganti et al. (2007) Morganti et al. 2007, A&A, 476, 735
- Morganti (2011) Morganti, R. 2011, (arXiv: 1112.5093)
- Murray et al. (2005) Murray, N., Quataert, E., & Thompson, T.A. 2005, ApJ, 618, 569
- Nagashima et al. (2006) Nagashima, M., Inutsuka, S., & Koyama, H. 2006, ApJ, 652, L41
- Netzer et al. (2007) Netzer et al. 2007, ApJ, 666, 806
- O’Neil (2002) O’Neil, K. 2002, ASPC, 278, 293
- Oosterloo et al. (2000) Oosterloo et al. 2000, AJ, 119, 2085
- Piconcelli et al. (2005) Piconcelli, E. et al. 2005, A&A, 432, 15
- Rafter et al. (2009) Rafter, S.E., Crenshaw, D.M., & Wiita, P.J. 2009, AJ, 137, 42
- Rohlfs & Wilson (2000) Rohlfs, K., & Wilson, T. L. 2000, Tools of radio astronomy / K. Rohlfs, T.L. Wilson. New York : Springer, 2000. (Astronomy and astrophysics library,ISSN0941-7834)
- Rupke et al. (2002) Rupke, D. S. N. et al. 2002, ApJ, 570, 588
- Rupke et al. (2005a) Rupke, D. S. N., Veilleux, S., & Sanders, D. B. 2005a, ApJS, 160, 115
- Rupke et al. (2005b) Rupke, D. S. N., Veilleux, S., & Sanders, D. B. 2005b, ApJ, 632, 751
- Rupke et al. (2005c) Rupke, D.S.N. et al. 2005c, ApJ, 632, 751
- Rupke & Veilleux (2011) Rupke, D.S.N. & Veilleux, S. 2011, ApJ, 729, L27
- Sanders et al. (1988) Sanders, D.B. et al. 1998, ApJ, 328, L35
- Sanders et al. (2003) Sanders D.B. et al. 2003, AJ, 126, 1607
- Schmidt & Green (1983) Schmidt, M. & Green, R.F. 1983, ApJ, 269, 352
- Schneider et al. (1986) Schneider, S.E., Helou, G., Salpeter, E.E., & Terzian, Y. 1986, AJ, 92, 742
- Schneider et al. (1990) Schneider, S.E., Thuan, T.X., Magri, C., & Wadiak, J.E. 1990, ApJS, 72, 245
- Schweitzer et al. (2006) Schweitzer et al. 2006, ApJ, 649, 79
- Scoville et al. (2000) Scoville, N.Z. et al. 2000, AJ, 119, 991
- Springel et al. (2005) Springel, V. et al. 2005, Nature, 435, 629
- Sturm et al. (2006) Sturm et al. 2006, ApJ, 652, L13
- Sturm et al. (2011) Sturm et al. 2011, ApJ, 733, L16
- Tchekhovskoy et al. (2010) Tchekhovskoy, A., Narayan, R., & McKinny, J.C. 2010, ApJ, 711, 50
- Teng & Veilleux (2010) Teng, S.H. & Veilleux, S. 2010, ApJ, 725, 1848
- Ulvestad et al. (1999) Ulvestad, J.S., Wrobel, J.M., & Carilli, C.L. 1999, ApJ, 516, 127
- Ulvestad et al. (2005) Ulvestad, J.S., Antonucci, R.R.J., & Barvainis, R. 2005, ApJ, 621, 123
- Veilleux et al. (2002) Veilleux, S., Kim, D.C., & Sanders, D.B. 2002, ApJS, 143, 315
- Veilleux et al. (2005) Veilleux, S., Cecil, G., & Bland-Hawthorn, J. 2005, ARAA, 43, 769
- Veilleux et al. (2006) Veilleux, S. et al. 2006, ApJ, 643, 707
- Veilleux et al. (2009a) Veilleux, S. et al. 2009a, ApJS, 182, 628
- Veilleux et al. (2009b) Veilleux, S. et al. 2009b, ApJ, 701, 587
- Veilleux (2012) Veilleux, S. 2012, Journal of Physics: Conference Series, 372, 12001
- Whiteoak & Gardner (1972) Whiteoak, J.B. & Gardner, F.F. 1972, A&A, 21, 159
- Wei et al. (2010) Wei, L.H., Kannappan, S.J., Vogel, S.N., & Baker, A.J. 2010, ApJ, 708, 841
- Wilson & Colbert (1995) Wilson, A.S., & Colbert, E.J.M. 1995, ApJ, 438, 62
- Wilson et al. (2008) Wilson, C.D. et al. 2008, ApJS, 178, 189
Appendix A Notes on Individual Galaxies
In this appendix, we describe the HI spectral properties of our GBT sample. When appropriate, we compare the present data with previous observations of these galaxies.
- PG 0007+106
The only radio-loud object in our sample based on classic radio-to-optical criteria. This source is classified as a triple system in the third Zwicky catalog (III Zw 2), and the observed velocity range is consistent with the A ( km s) and C ( km s) components of the system. The GBT detection is assumed to be that of the brightest component, A. The HI emission profile of the full intensity spectrum agrees with the old Arecibo spectrum of Hutchings et al. (1987). The absorption feature, also noted by Hutchings et al. (1987), is more prominent in the XX spectrum of the GBT data.
- PG 0050+124 (I Zw 1)
The heliocentric radial velocity listed in Table 1 for this source is from optical spectroscopy by Ho & Kim (2009). The optically-derived velocity is 700 km s lower than the HI-derived value of 18285 km s (Bothun et al., 1984). The optical emission lines may arise in part from outflowing gas, so the optically-derived value may is underestimate the actual systemic velocity of this object. The shape of the full intensity profile (Figure 1) at 18285 km s agrees with those previously published by Bothun et al. (1984) and Hutchings et al. (1987), and matches expectations for a regularly rotating disk. However, at lower velocities, the profile becomes polarized and a broad absorption feature is centered at the optically-derived velocity (see Figure 9).
There is a marginal detection of an emission feature and a weak but broad blueshifted absorption feature present.
No emission is seen in the data. However, several absorption components are detected. An alternative explanation for the spectral shape is the superposition of an emission profile with a broad absorption component.
- PG 0804+761
Strong RFI is present at the expected frequency of this galaxy. It is therefore excluded from our analysis.
- PG 0844+249
The broad emission profile is asymmetric. Although the profiles in the two different polarizations differ slightly, they have shapes similar to that observed by Ho et al. (2008). Our integrated GBT flux density is larger than that measured by Ho et al. (2008). The emission feature redward of the peak may be due to a companion at km s.
- UGC 05101
Multiple absorption components are seen in the GBT data. We interpret the HI absorption profile as a highly absorbed core coincident with an inclined regularly rotating disk.
- PG 1119+120
This object was previously observed by Bothun et al. (1984) and Hutchings et al. (1987) with Arecibo. The profiles of the older spectra resemble our GBT observation, although the narrow absorption feature is not apparent in the lower-resolution Arecibo data. There is no obvious detection of a companion in the GBT data near 16000 km s as suggested by Hutchings et al. (1987). High resolution VLA maps at 4.8 GHz show a bent jet extending to the northeast (Leipski et al., 2006). At the resolution of these VLA maps (0.75), the jet is resolved at a distance of 0.72 kpc.
- PG 1126-041
This target was statistically undetected after three hours on-source, although there is a very weak broad absorption feature centered near its systemic velocity.
- PG 1211+143
- PG 1229+204
This object was observed by Ho et al. (2008). The profile of the GBT spectrum roughly agrees with the earlier data.
- PG 1244+026
A broad, blueshifted absorption feature is present. As for F05189–2524, an alternative explanation for the spectral shape may be emission absorbed by a single broad absorption feature.
- Mrk 231
The deep absorption trough is associated with the face-on disk of the galaxy (Carilli et al., 1998; Ulvestad et al., 1999). Our GBT full intensity HI profile of this sources mostly agrees with these previous observations, but we detect an additional absorption feature centered at relative velocity –300 km s that was first suspected by Morganti (2011). This feature is slightly polarized, with the XX spectrum having lower intensity. Although the polarization is not as high as those seen in Figure 4, if we give the spectra of Mrk 231 the same treatment detailed in Section 5.2, we find that this feature has FWHM 1000 km s corresponding to a mass outflow rate of 0.3 yr in the jet scenario. The polarized feature is thus likely associated with the radio jet (Carilli et al., 1998; Ulvestad et al., 1999). The velocity of this jet-enhanced outflow is similar to the value measured using Na I in the jet-disturbed region (–1400 km s; Rupke & Veilleux, 2011).
- Mrk 273
This source is the only binary in our sample; its nuclear separation is 750 kpc (Scoville et al., 2000). X-ray observations by Iwasawa et al. (2011) identified its southwestern nucleus as the Seyfert 2 nucleus and suggested that the northern nucleus may host a heavily obscured AGN. The double-peaked HI absorption profile is very symmetric and is likely due to the two obscured nuclei. There appears to be a weak absorption feature in the YY spectrum at –500 km s.
- PG 1351+640
The XX spectrum appears to show emission, while YY shows broad absorption. The HI absorption data may also be variable. Leipski et al. (2006) note a jet-like structure in their 4.8 GHz VLA maps and time variability in the core flux density. The resolution of the Leipski et al. (2006) data (0.56) corresponds to a distance of 0.91 kpc. Although on a different spatial scale, this structure is consistent (similar position angle) with a single-sided extension hinted at by the 8.4 GHz VLBA snapshot data presented by Blundell & Beasley (1998). This feature is not seen in deeper VLBA imaging at multiple frequencies, but it should be noted that the data from Ulvestad et al. (2005) are only in the left-circular polarization.
- PG 1411+442
The spectrum of this object is very similar to that of PG 1351+640. Emission is seen in the XX spectrum and a two-component absorption feature in the YY spectrum.
- PG 1426+015
Ho et al. (2008) detect a fairly narrow () and well-defined HI line centered at . We have examined our spectrum at the corresponding frequency and see hints of emission over a wider range of velocities; however, this feature has such low significance and is so sensitive to the details of baseline fitting we classify its detection status as ”uncertain.”
- PG 1440+356
The profile of the HI absorption feature depends slightly on the polarization. Exceptionally, the feature is redshifted by 300 km s with respect to systemic velocity rather than being blueshifted. A possible explanation is that the systemic velocity taken from NED and based on optical emission lines underestimates the actual systemic velocity because the line-emitting gas arises in part from outflowing material. Another possibility is that we may be detecting infalling HI. Emission is also present in the blue half of the spectrum, although this is only seen at one of several epochs.
- PG 1448+273
The weakly polarized spectra show blueshifted absorption.
- PG 1501+106
A weak emission feature is present in the data. A blueshifted absorption feature is also visible, particularly in the XX spectrum. A lot of the data for this source were affected by small ripples in the baseline. As a reality check, we compared our data with similar GBT data from program 10A-070 (PI: Kannappan). The ripples also appear to be present in a large fraction of their data. The persistent ripples could be due to harmonic waves from a strong RFI signal at a frequency that is slightly outside the detector frequency range. Due to the expected weakness of the features, we conservatively flagged all of the data that showed these ripples.
Multiple weak absorption components are present in this object.
This is the only ULIRG in our sample that shows distinct polarization-dependent HI absorption. The spectra are very similar to those of PG 1351+640 and PG 1411+442, where weak emission is seen in the XX spectrum and absorption in YY.
Strong RFI is present at the expected observed frequency of this galaxy. It is therefore excluded from our analysis.
- PG 1617+175
Broad emission is detected in this object. The profiles of the emission differ slightly in the two polarizations. Notably, the emission feature in the YY spectrum seems to be affected by an additional absorption feature (Figure 9).
This object was not detected.
- PG 2130+099
This target was observed by Ho et al. (2008). Our full intensity profile of this source agrees with the older data from Arecibo. A blueshifted absorption feature is also seen in the XX spectrum. Leipski et al. (2006) detected what looks like jet-associated hot spots emanating from the core in 4.8 GHz A-array VLA data. The resolution of these VLA data (0.35) implies that the jet can be resolved at a distance of 0.42 kpc.
- PG 2214+139
The HI spectrum of this object is dominated by a broad emission feature centered at systemic velocity.
|J2000||J2000||(km s)||Strength||( yr)||(Jy)||Observations||(s)||(km s)||(mJy)||Order|
|PG 0007+106||00:10:31.0||+10:58:30||26783||QSO||Weak||12.23||1.2||0.099||–2.4||2011 Sept – 2012 Jan||7894||6.11||0.77||3|
|PG 0050+124||00:53:34.9||+12:41:36||17658||QSO||Strong||12.07||10.4||0.009||–5.0||2011 Sept – 2012 Jan||4815||5.78||0.94||3|
|F04103–2838||04:12:19.4||–28:30:25||35215||ULIRG||Strong||12.30||154.2||0.013||–3.4||2011 Sept – 2012 Mar||7452||6.43||0.88||3|
|F05189–2524||05:21:01.4||–25:21:45||12760||ULIRG||Strong||12.22||71.6||0.029||2011 Sept – 2012 Jan||12330||5.60||0.58||3|
|PG 0804+761||08:10:58.6||+76:02:43||29979||QSO||Strong||12.08||0.3||0.004||–5.3||2011 Aug||860||6.23||2.97|
|PG 0844+349||08:47:42.4||+34:45:04||19187||QSO||Weak||11.44||0.8||0.0025||–3.3||2011 Nov||3153||5.83||0.96||3|
|UGC 05101||09:35:51.6||+61:21:11||11802||ULIRG||Strong||12.05||110.2||0.170||–3.0||2011 Aug – 2011 Sept||2151||5.57||1.46||3|
|PG 1119+120||11:21:47.1||+11:44:18||15050||QSO||Strong||11.33||0.4||0.004||2011 Oct – 2012 Mar||3837||5.68||0.83||3|
|PG 1126–041||11:29:16.6||–04:24:08||18575||QSO||Strong||11.52||2.5||0.0025||–4.6||2011 Aug – 2012 Jan||8871||5.81||0.54||3|
|PG 1211+143||12:14:17.7||+14:03:13||24253||QSO||Strong||11.96||0.0||0.600||–2.8||2011 Aug – 2012 Mar||5986||6.02||0.92||3|
|PG 1229+204||12:32:03.6||+20:09:29||18890||QSO||Weak||11.56||0.4||0.003||–4.9||2012 Mar||3411||5.82||0.92||3|
|PG 1244+026||12:46:35.2||+02:22:09||14443||QSO||Strong||11.02||1.5||0.002||–5.1||2011 Aug – 2012 Mar||3555||5.66||1.00||3|
|Mrk 231||12:56:14.2||+56:52:25||12642||ULIRG||Strong||12.60||174.2||0.309||–2.7||2011 Aug – 2011 Sept||1257||5.60||1.52||3|
|Mrk 273||13:44:42.1||+55:53:13||11326||ULIRG||Strong||12.24||141.6||0.136||–3.5||2011 Aug – 2011 Sept||1290||5.50||1.31||3|
|PG 1351+640||13:53:15.8||+63:45:46||26442||QSO||Strong||12.04||26.1||0.034||–3.3||2011 Sept – 2012 Feb||8153||6.10||0.79||3|
|PG 1411+442||14:13:48.3||+44:00:14||26861||QSO||Weak||11.78||0.0||0.0025||–4.6||2011 Aug – 2012 Feb||6346||6.12||0.92||3|
|PG 1426+015||14:29:06.6||+01:17:06||25953||QSO||Weak||11.92||1.9||0.004||–5.0||2012 Mar||3153||6.08||1.09||5|
|PG 1440+356||14:42:07.4||+35:26:23||23700||QSO||Strong||11.80||10.3||0.005||–4.8||2011 Aug – 2011 Dec||6135||6.00||0.63||3|
|PG 1448+273||14:51:08.7||+27:09:27||19487||QSO||Weak||11.43||0.1||0.003||–4.8||2011 Aug – 2011 Feb||8917||5.85||0.63||3|
|PG 1501+106||15:04:01.2||+10:26:16||10919||QSO||Strong||11.33||0.0||0.0025||–5.7||2012 Feb – 2012 Mar||6540||5.53||0.77||3|
|F15130–1958||15:15:55.2||–20:09:17||32585||ULIRG||Strong||12.23||108.8||0.010||–2.4||2011 Dec – 2012 Mar||2666||6.33||1.25||3|
|F15250+3608||15:26:59.4||+35:58:38||16535||ULIRG||Strong||12.12||109.6||0.018||2011 Aug – 2011 Nov||2007||5.74||1.16||3|
|PG 1617+175||16:20:11.3||+17:24:28||33708||QSO||Weak||11.74||4.0||0.0025||2011 Aug – 2012 Mar||2547||6.38||1.59||3|
|F21219–1757||21:24:41.6||–17:44:46||33526||ULIRG||Strong||12.17||48.3||0.015||2011 Oct – 2012 Mar||8529||6.37||0.70||3|
|PG 2130+099||21:32:27.8||+10:08:19||18880||QSO||Strong||11.77||0.2||0.007||–2.7||2012 Mar||2001||5.82||1.13||5|
|PG 2214+139||22:17:12.2||+14:14:21||19715||QSO||Weak||11.77||0.0||0.0025||–5.0||2011 Sept – 2011 Oct||2121||5.85||1.32||3|
References. – Condon et al. (NVSS; 1998); Rafter et al. (2009); Barvainis et al. (1996). Note. – Col.(1): Galaxy name. Coordinate-based names beginning with “F” are sources in the IRAS Faint Source Catalog; the “PG” sources are Palomar Green quasars. Col.(2): Right ascension in J2000 coordinates. Col. (3): Declination in J2000 coordinates. Col.(4): Heliocentric radial velocity given by NED. These velocities are measured via optical emission lines and may underestimate the true recessional velocities since some fraction of the line emission may arise from the outflowing material itself. All velocities are quoted in the optical convention where . Col.(5): ULIRG or QSO source. Col.(6): FIR strength as defined by Netzer et al. (2007) where FIR-strong sources have . Col.(7): Bolometric luminosity. For ULIRGs, we assume . For PG QSOs, we assume (Netzer et al., 2007). Col.(8): Star formation rates in solar masses per year derived using the average AGN contribution fraction to the total bolometric luminosity measured by Veilleux et al. (2009a) for each galaxy. For PG 0007+106, the AGN fraction is assumed to be the mean for all other FIR-weak quasars. It is assumed that all of the emission from the starburst is reprocessed by dust and re-emitted into the IR band. The SFR is derived using the relation given in Kennicutt (1998) for a 0.1–100 initial mass function. Col.(9): 1.4 GHz continuum flux densities. Values taken with the GBT unless noted. The upper limits are given as the detection limit of the NRAO VLA Sky Survey (NVSS). Col.(10): 1.4 GHz to 2–10 keV flux ratio, where , implies a radio-loud source (La Franca et al., 2010). We used the (average, if multiple observations were performed) absorption-corrected 2–10 keV flux from Teng & Veilleux (2010) except for PG 0007+106 and PG 0804+761 whose X-ray measurements were taken from Piconcelli et al. (2005). These ratios are corrected for the 1.4 GHz contribution from star formation based on the SFRs listed in Table 1 and the SFR-1.4 GHz relation from Bell (2003), which assumes the same 0.1–100 initial mass function as Kennicutt (1998). For ULIRGs, due to their high line-of-sight column densities, their measured 2–10 keV fluxes may be underestimated and thus the listed values may represent upper limits. Col.(11): Period over which the target was observed. Col.(12): Total on source time in seconds after flagging. Col.(13): Velocity resolution of the full intensity spectrum. Col.(14): The channel-to-channel rms of the full intensity spectrum. Col.(15): The order of polynomial over which the baseline was fitted and then subtracted.
|(km s)||(km s)||(km s)||(mJy)||(Jy km s)||(10 cm)||(10)|
|PG 1426+015 (U)||5.12||4.69|
|PG 1617+175 (U)||6.94||4.36|
Note. – Col.(1): Galaxy name where (M) and (U) indicate marginal and uncertain detections, respectively, as noted in Section 3.2. Col.(2): Heliocentric velocity measured from HI relative to the systematic velocity presented in Table 1, and its associated measurement error. Col.(3): Velocity width of the emission measured at the 20% level and its error. Col.(4): Velocity width of the emission measured at the 50% level and its error. Col.(5): Maximum intensity of the spectrum. Col.(6): Signal-to-noise ratio given by the peak intensity relative to the channel rms listed in Table 1. Col.(7): Velocity-integrated HI flux in emission and its error. Col.(8): Derived HI column density calculated using Equation 1. Col.(9): Derived HI gas mass from emission and its error calculated using Equation 2.
|(km s)||(km s)||(mJy)||(10 cm )||(Jy km s)|
|Strong Detections: peak-to-rms|
|Weak Detections: peak-to-rms|
Note. – Col.(1): Galaxy name. Col.(2): Heliocentric velocity measured from HI relative to the systematic velocity presented in Table 1, as modeled using a Gaussian model, and its associated measurement error. Col.(3): Full width of the absorption feature at the 50% flux level as measured using a Gaussian model and its error. Col.(4): Peak absorption flux and its error as measured using a Gaussian model. Col.(5): Signal-to-noise ratio measured as a ratio between the peak flux and channel rms listed in Table 1. Col.(6): Integrated HI column density based on the absorption model, assuming a spin temperature of 1000 K and using Equation 4 based on the optical depth profile derived using Equation 3. For PG 1411+442 and PG 1501+106, the 1.4 GHz core flux densities are only upper limits, so their derived values are lower limits. Col.(7): Peak optical depth based on Equation 3. Col.(8): Velocity-integrated HI flux in absorption measured from the best-fit Gaussian component models.
|(km s)||(km s)||(km s)||(mJy)||(Jy km s)||(10 cm)||(10)|
|PG 1617+175 (U)||9.23||5.81|
Note. – Col.(1): Galaxy name. Col.(2): Heliocentric velocity measured from HI relative to the systematic velocity presented in Table 1, and its associated measurement error. Col.(3): Velocity width of the emission measured at the 20% level and its error. Col.(4): Velocity width of the emission measured at the 50% level and its error. Col.(5): The maximum intensity of the spectrum. Col.(6): Signal-to-noise ratio given by the peak intensity relative to the channel rms listed in Table 1. Col.(7): Velocity-integrated HI flux in emission and its error. Col.(8): Derived HI column density calculated using Equation 1. Col.(9): Derived HI gas mass from emission and its error, calculated using Equation 2.
|(km s)||(km s)||(mJy)||(10 cm )||(Jy km s)|
|PG 1617+175 (U)||2.08||–0.67|
Note. – Col.(1): Galaxy name. Col.(2): Heliocentric velocity measured from HI relative to the systematic velocity presented in Table 1, as modeled using a Gaussian fit, and its associated measurement error. Col.(3): Full width of the absorption feature at the 50% flux level as measured using a Gaussian model and its error. Col.(4): Peak absorption flux and its error as measured using a Gaussian model. Col.(5): Signal-to-noise ratio measured as a ratio between the peak flux and channel rms listed in Table 1. Col.(6): Integrated HI column density based on the absorption model, assuming a spin temperature of 1000 K, and using Equation 4 based on the optical depth profile derived using Equation 3. For PG 1501+106, the 1.4 GHz core flux density is only an upper limit, so its derived value is a lower limit. Col.(7): Peak optical depth based on Equation 3. Col.(8): Velocity-integrated HI flux in absorption measured from the best-fit Gaussian component models.
|(km s)||(km s)||(10 cm )||(km s)||(km s)||(km s)||( yr)||( yr)|
|PG 1617+175 (U)||–95|
Note. – The sum of the mass outflow rates from the multiple Gaussian components in the galaxy. Col.(1): Galaxy name. Col.(2): Heliocentric velocity of the absorption component as tabulated in Tables 3 and 5. Col. (3): With the exception of PG 1211+143 and PG 1440+356, this is the outflow velocity measured relative to the HI systemic velocity measured from HI emission using values tabulated in column (2) of Table 2 or 4. if no emission is detected. For PG 1211+143 and PG 1440+356, this velocity is relative to the systemic velocity measured from optical emission lines, since the detection of HI emission is uncertain (Figures 5 and 6). Col. (4): Integrated HI column density based on the absorption model from Tables 3 and 5. Col. (5): Full-width at zero intensity based on the best-fit Gaussian absorption model and its associated error. Col. (6): The maximum outflow velocity where FWZI. Col. (7): The velocity at which 98% of the gas has lower velocity where . Col. (8): Derived mass outflow rate calculated using Equation 5 and assuming the outflow is in a wide-angled cone where the mass is flowing into a solid angle of from a radius of 1 kpc following the same assumptions as Morganti et al. (2005b). Col. (9): Same as column (8) but for a jet whose opening angle is assumed to be 5 launched at a distance of 0.5 kpc. Col. (10): Log of the kinetic energy outflow rate as calculated using Equation 6 for the jet scenario.