HST IMAGING OF FADING AGN CANDIDATES. I. HOST-GALAXY PROPERTIES AND ORIGIN OF THE EXTENDED GAS*

, , , , , , , , , , , , , and

Published 2015 April 14 © 2015. The American Astronomical Society. All rights reserved.
, , Citation William C. Keel et al 2015 AJ 149 155 DOI 10.1088/0004-6256/149/5/155

1538-3881/149/5/155

ABSTRACT

We present narrow- and medium-band Hubble Space Telescope imaging, with additional supporting ground-based imaging, spectrophotometry, and Fabry–Perot interferometric data, for eight galaxies identified as hosting a fading active galactic nucleus (AGN). These are selected to have AGN-ionized gas projected $\gt 10$ kpc from the nucleus and energy budgets with a significant shortfall of ionizing radiation between the requirement to ionize the distant gas and the AGN as observed directly, indicating fading of the AGN on ≈50,000 yr timescales. This paper focuses on the host-galaxy properties and origin of the gas. In every galaxy, we identify evidence of ongoing or past interactions, including tidal tails, shells, and warped or chaotic dust structures; a similarly selected sample of obscured AGNs with extended ionized clouds shares this high incidence of disturbed morphologies. Several systems show multiple dust lanes in different orientations, broadly fit by differentially precessing disks of accreted material viewed ∼1.5 Gyr after its initial arrival. The host systems are of early Hubble type; most show nearly pure de Vaucouleurs surface brightness profiles and Sérsic indices appropriate for classical bulges, with one S0 and one SB0 galaxy. The gas has a systematically lower metallicity than the nuclei; three systems have abundances uniformly well below solar, consistent with an origin in tidally disrupted low-luminosity galaxies, while some systems have more nearly solar abundances (accompanied by such signatures as multiple Doppler components), which may suggest redistribution of gas by outflows within the host galaxies themselves. These aspects are consistent with a tidal origin for the extended gas in most systems, although the ionized gas and stellar tidal features do not always match closely. Unlike extended emission regions around many radio-loud AGNs, these clouds are kinematically dominated by rotation, in some cases in warped disks. Outflows can play important kinematic roles only in localized regions near some of the AGNs. We find only a few sets of young star clusters potentially triggered by AGN outflows. In UGC 7342 and UGC 11185, multiple luminous star clusters are seen just within the projected ionization cones, potentially marking star formation triggered by outflows. As in the discovery example, Hanny's Voorwerp/IC 2497, there are regions in these clouds where the lack of a strong correlation between Hα surface brightness and ionization parameter indicates that there is unresolved fine structure in the clouds. Together with thin coherent filaments spanning several kpc, persistence of these structures over their orbital lifetimes may require a role for magnetic confinement. Overall, we find that the sample of fading AGNs occur in interacting and merging systems, that the very extended ionized gas is composed of tidal debris rather than galactic winds, and that these host systems are bulge-dominated and show no strong evidence of triggered star formation in luminous clusters.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

An important subset of active galactic nuclei (AGNs) is the small fraction with extended emission-line regions (EELRs), ionized gas extending as far as tens of kpc from the core. These regions provide information on the geometry of escaping ionizing radiation (as in ionization cones) on the history of the AGN via light-travel time to the clouds, and on the environment of the AGN. In radio-loud objects, these regions often have complex and large-amplitude velocity structures, indicating an important role for outflows. In contrast, radio-quiet objects with EELRs may show lower-amplitude and globally ordered velocity fields, consistent with tidal debris photoionized by the AGN but not otherwise disturbed.

Recent results have highlighted an important role for EELRs in probing the history of an AGN's radiative output on timescales up to 105 yr. An early discovery of the Galaxy Zoo citizen-science project (Lintott et al. 2008) was Hanny's Voorwerp, a galaxy-sized cloud of highly ionized gas associated with the spiral galaxy IC 2497 (Lintott et al. 2009). The object's spectrum clearly indicates photoionization by a luminous AGN, while direct observations of IC 2497 at all wavelengths show a deficit of at least two orders of magnitude in ionizing luminosity (Lintott et al. 2009; Schawinski et al. 2010). Furthermore, Hubble Space Telescope (HST) imaging and spectroscopy fail to show any high-ionization gas in IC 2497 which would result from an extant AGN (Keel et al. 2012b). These results are all explained if the AGN in IC 2497 was a QSO until $\approx {{10}^{5}}$ yr before our current view of it, and we see Hanny's Voorwerp as an ionization echo--the analog of a light echo, but expressed through recombination, which has a long timescale at low densities. Such dramatic changes in radiative output suggest correspondingly dramatic changes in the accretion onto the central black hole, providing a view on hitherto inaccessible timescales. Using the properties of Hanny's Voorwerp as a template and enabled by nearly 200 Galaxy Zoo volunteers, we searched for similar objects, using the distinct colors of nebulae with strong [O iii] $\lambda 5007$ emission when seen at low redshift in the SDSS filter set (Keel et al. 2012a); by analogy with Hanny's Voorwerp, these were known collectively as "Voorwerpjes" using the Dutch diminutive form.

A spectroscopic study of the most promising candidates from color and morphology yielded a sample of 19 additional AGNs with high-ionization gas extending more than 10 kpc in projection from the AGN, suitable for studying the history of these AGNs. Among these, eight show a strong deficit in ionizing luminosity based on both the nuclear emission-line fluxes and far-infrared luminosity, thereby excluding objects with strongly obscured AGNs from which radiation escapes in some directions to ionize surrounding gas. Fading AGNs are easier to find with this approach than brightening ones, since these have the same signature as a lack of surrounding cool gas at an appropriate position to be ionized. Broadly, the relative numbers of fading versus obscured AGNs in the Keel et al. (2012a) sample suggests that these AGNs are luminous in episodes lasting $0.2-2\times {{10}^{5}}$ yr, probably with many such episodes over much longer lifetimes. At least 3/4 of these are in interacting or merging systems, again consistent with ionizing gas-rich tidal debris at large radii in directions not blocked by a disk in the AGN host galaxy itself. Subsequent work by other groups has extended the case for long-timescale variations in AGNs. A set of objects with very similar properties at higher redshifts and higher luminosities has been reported by Schirmer et al. (2013), and deep [O iii] images of the merger remnant NGC 7252 by Schweizer et al. (2013) indicate a recent active episode of its nucleus (where there is no directly detected AGN at all). In a statistical context, Hickox et al. (2014) find that similarly dramatic AGN variations can account for the scatter in relationships between star formation indicators and AGN luminosity.

We have carried out an observational campaign to characterize the Voorwerpje systems in more detail, concentrating on the fading-AGN candidates, aiming both to probe the history of the AGN itself and to examine the properties of the host galaxies and their interaction and star-forming histories. We concentrate in this paper on the structure of the host galaxies, origin and fine structure of the ionized gas, and potential evidence for star formation triggered by outflow from the AGN. Companion papers deal with the implications of these data for these AGNs and their histories (C. J. Lintott et al. 2015, in preparation), and with the two sample AGNs whose double radio sources and giant ionized clouds violate the usual alignment pattern for radio jets and ionization cones (Mkn 1498 and NGC 5972; W. P. Maksym et al., in preparation). The high-resolution HST images support the interpretation that the AGNs in these objects are indeed fading, give hints as to the timescales involved, and in some cases suggest nuclear outflows possibly connected to the fading or mode change of the AGNs. We describe in turn the various kinds of observations we have carried out, the analysis techniques used in our measurements, and the results organized both by kind of data and (in an appendix) by galaxy.

In evaluating sizes and luminosities we compute distances using a Hubble constant of 73 km s−1 Mpc−1.

2. SAMPLE, OBSERVATIONS, AND DATA REDUCTION

2.1. Galaxy Sample

We observed the subset of EELR host galaxies from Keel et al. (2012a) which have ionizing-energy deficits—ratios of bolometric AGN luminosity to the required level for ionizing the most distant clouds—greater than about a factor of three. (The companion paper by Lintott et al. includes updated energy-budget calculations, incorporating WISE data to fill mid-IR gaps in the SED). Table 1 lists the target galaxies, AGN types, and redshifts.

Table 1.  Galaxy Properties and Filters

Galaxy AGN Type z [O iii] Filter Hα+[N ii] Filter Scale, kpc/''
Mkn 1498 Sy 1.9 0.0547 ACS FR551N ACS FR716N 1.11
NGC 5252 Sy 1.5 0.0228 WFPC2 FR533N WFPC2 F673N 0.50
NGC 5972 Sy 2 0.0297 ACS FR505N WFC3 F673N 0.63
SDSS 1510+07 Sy 2 0.0458 ACS FR505N ACS FR716N 0.89
SDSS 2201+11 Sy 2 0.0296 WFC3 FQ508N WFC3 F673N 0.59
Teacup AGN Sy 2 0.0852 ACS FR551N ACS FR716N 1.58
UGC 7342 Sy 2 0.0477 ACS FR505N ACS FR716N 0.98
UGC 11185 Sy 2 0.0412 ACS FR505N ACS FR656N 0.80

Download table as:  ASCIITypeset image

2.2.  HST Imaging

Our primary data set consists of HST images isolating [O iii] $\lambda 5007$ and Hα+[N ii] $\lambda 6583$ and neighboring continuum bands. For the starlight continuum, we used the F621M (center 6190 Å, FWHM = 700 Å) and F763M (center 7600 Å, FWHM = 780 Å) intermediate-band filters of the Wide-Field Camera 3 (WFC3) to combine reasonable surface brightness sensitivity with good rejection of the strong emission lines. For the emission lines, we used the ramp filters on the Advanced Camera for Surveys (ACS) except for a few instances where a galaxy's redshift put one of these lines in a narrow WFC3 filter. The ACS narrow ramp filters have a nominal FWHM of 2% of the center wavelength. Total exposures were 1250–1400 s in F621M, 1150–1260 s in F763M, and 2800–2900 s in the narrow bands, with the variations depending on viewing period per orbit. The observation sequences were set to give slightly longer exposures in Hα, known to be weaker than [O iii] emission. Each observation included two exposures with small 0farcs 145 dither offsets between them, to improve sampling of the final product and reduce the effects of hot or low pixels. For NGC 5972, already known to have very large emission regions, the dither step of 2farcs 41 was made large enough to fill the inter-chip gap with one exposure. The different cameras often required different guide stars, leading to small offsets ($\simeq 0\buildrel{\prime\prime}\over{.} 2$) in the world coordinate systems between them; we used the nuclei and bright stars to register the sets of images.

The CCDs in both ACS and WFC3 are affected by charge-transfer efficiency (CTE) losses, especially severe in the low-background regime of the narrowband data. Uncorrected, this leads to much of the charge originating in a pixel being spread along its readout column; not only does this lead to underestimating the flux of faint targets, but cosmic-ray events leave trails that are harder to reject cleanly. These are substantially improved through application of pixel-based reconstructions, essentially a deconvolution using physical properties of the CTE losses. We used the ACS routine of Anderson & Bedin (2010), which is now incorporated in the MAST archival processing, and the recently released WFC3 version described by Baggett et al. (2014). The output of these correction routines was combined for each filter using multidrizzle (Gonzaga et al. 2012). The LAcosmic routine (van Dokkum 2001) worked well to reject and patch cosmic-ray events for the WFC3 data; we used interactive removal with the IRAF imedit task for ACS images.

Table 1 lists the filters used for the narrow emission-line bands for each galaxy. For brevity, we refer to SDSS J151004.01+074037.1 as SDSS 1510+07, to SDSS J220141.64+115124.3 as SDSS 2201+11, and note that the Teacup AGN is formally known as SDSS J143029.88+133912.0.

For most objects, line/continuum separation is very clean with these filters. For Mkn 1498 and the Teacup AGN, we used image arithmetic to correct for line emission in the nominal continuum bands, at levels of 10% in the narrowband count rate for [O iii] and 30% for Hα+[N ii] . For Mkn 1498, the scale factors were 0.10 for [O iii] and 0.15 for Hα+[N ii] . To improve continuum subtraction in [O iii] around the dust lanes of SDSS 2201+11, we used a synthetic continuum image extrapolated from the observed images for a closer match to the wavelength of the line emission, iteratively generated to minimize residuals in regions with weak or no line emission.

In the WFC3 images of UGC 7342, reflections from bright field stars produce bright elliptical rings superimposed on the outer parts of the host galaxy in smooth regions so that little information is lost.

NGC 5252, with its well-known ionization cones, satisfies our sample criteria relating nuclear luminosity and ionization level in the outer regions. We use the archival HST WFPC2 images described by Tsvetanov et al. (1996) and Morse et al. (1998), with analysis parallel to our handling of the new ACS and WFC3 images. The monochromatic field of the WFPC2 ramp filters is much smaller than the 40 × 80'' field of the ACS ramps, so the outermost filaments are progressively lost in the [O iii] image.

For the two largest systems, NGC 5972 and UGC 7342, the telescope orientation for the ACS observations was constrained to align their major axes nearly parallel to the orientation of the 40 × 80'' monochromatic field of view. The only known emission-line structure lost to the field size is an outlying, radial filament extending from 22–45'' to the south of NGC 5972; ground-based data give independent information on its [O iii] flux.

Figure 1 gives an overview of the target galaxies in our HST imagery, comparing continuum structures and [O iii] emission.

Figure 1.
Standard image High-resolution image
Figure 1.

Figure 1. Comparison of continuum and emission-line structure of the sample galaxies. Blue and red are the continuum WFC3 images, and their average with [O iii] emission (from ACS or WFC3 narrowband data) added is mapped to green. The scaling of [O iii] to continuum is arbitrary but consistent for six objects; for SDSS 1510+07 and SDSS 2201+11, with lower [O iii] surface brightness, we used a scaling three times larger for the emission component. For each color band, we use an offset log intensity mapping (intensity $\propto \;{\rm log} ([I+{\rm constant}]$, similar to the sinh mapping used in SDSS composite images). Intensity scales are consistent between images except for Mkn 1498, where the maximum mapped intensity is twice as high as the rest to show detail near the AGN. Changes in the continuum color are shown, but very small in this visualization due to the modest wavelength baseline between the F621M and F763M filter bands. The angular scales are independently chosen to show the emission structure; scale bars are 5'' long with arrows pointing north. North is to the top and east to the left for all except UGC 7342.

Standard image High-resolution image

2.2.1. Surface Photometry

To assess their surface brightness profiles and trace dust features, we fit families of elliptical isophotes to the continuum images using the stsdas ellipse task (which implements the algorithm from Jedrzejewski 1987) and created smooth models following these profiles using the bmodel task. The isophote-tracing task was run in an interactive mode allowing the masking of obvious dust lanes or patches. In some regions where one of the ellipse parameters changes strongly with radius, the bmodel output had arcs of missing values that we replaced with the median of nonzero values over 7 × 7 pixel regions. For SDSS 2201+11, the dust is so widespread that our starting fit used the maximum intensity of the image after reflection about both axes relative to the large-scale major axis. Dust features were quantified from attenuation images, the ratio between data, and these models for each continuum image. Color-ratio images between the two continuum bands provide an additional measure of the location and extent of dust attenuation, although with lower signal-to-noise ratio (S/N) because of the small wavelength leverage between the F621M and F763M passbands. The color-ratio image, however, is not subject to the shortcomings of the elliptical-model fits in generating spurious structure in complex host galaxies, giving an independent view of the reality of dust features.

2.2.2. Identification and Photometry of Candidate Star Clusters

Although the continuum bands were selected for their proximity to wavelengths of strong emission lines, we also use them to characterize candidate star clusters. The color baseline between F621M and F763M is quite similar to $R-I$ in the Bessell (essentially Cousins) system. We have estimated the transformations via synthetic photometry on the Pickles (1998) spectral library using the response curves for these WFC3 filters and the RI curves given by Pickles. For $R-I\lt 0.8$, there are closely linear relations between the ${{m}_{621}}-{{m}_{763}}$ color and $R-I$; redder than this, for M stars, the F621M filter includes the TiO and CaH bands (e.g., Bochanski et al. 2007), and the relation becomes nonlinear. We ignore this regime since our interest is in young and luminous star clusters. Using Vega magnitudes where colors for A0 V stars are zero, and incorporating the median redshift z = 0.041 of our objects, we find ${{m}_{621}}-{{m}_{763}}=0.92{{(}_{R-I)0}}$ where the subscript denotes that $R-I$ is at redshift zero. The scatter about this linear fit is 0.026 magnitude for all stars earlier than type M. The redshift correction is quite small in ${{m}_{621}}-{{m}_{763}}$, introducing a scatter $\sigma =0.013$ magnitude, and significant systematic trends only for M stars. Therefore, to estimate the standard colors of star clusters, we invert the above expression and use ${{(}_{R-I)0}}=1.09({{m}_{621}}-{{m}_{763}})$.

In measuring magnitudes and colors of candidate star clusters, we normalize to a radius of r = 0farcs 4 (8 pixels), using the modest curve-of-growth corrections for a point source (0.045 magnitude for F631M, 0.033 for F763M). We correct for foreground Milky Way extinction using the values from Schlafly & Finkbeiner (2011). To reduce the effects of gradients in the local background, we performed aperture photometry after subtracting a smooth galaxy model, and quote errors based on statistics of local background regions after this subtraction.

2.3. Long-slit Spectra

We have obtained long-slit spectra of regions seen as particularly interesting in the HST images, using two facilities. At Lick Observatory, we used the Kast double spectrograph on the 3 m Shane telescope. Separate red and blue detectors, on optical trains fed by a dichroic beamsplitter, covered wavelength ranges 3400–4650 and 4900–7300 Å. The slit width was 2farcs 0.

At the 6 m telescope (BTA) of the Special Astrophysical Observatory of the Russian Academy of Sciences, the SCORPIO-2 multi-mode focal reducer (Afanasiev & Moiseev 2011) was used in the long-slit spectrograph mode, covering the wavelength range 3800–8100 Å at a resolution of typically 4.5 Å FWHM. The slit width was varied with the seeing, typically 1farcs 0. These spectra are detailed in Table 2.

Table 2.  Long-slit Spectra

Galaxy P.A. (°) Telescope Slit ('') Exposure (s) UT date
NGC 5972 165 BTA 1.0 900 2013 Apr 6
SDSS 2201+11 180 BTA 1.0 1800 2012 Aug 25
UGC 7342 180 BTA 0.7 2700 2013 Apr 5
UGC 7342 140 BTA 0.7 3600 2013 Apr 5
UGC 11185 90 BTA 1.0 1800 2012 Aug 25
Mkn 1498 90 Lick 2.0 1800 2013 Mar 13
NGC 5972 105 Lick 2.0 1800 2013 Jan 18
UGC 7342 154 Lick 2.0 1800 2013 Jan 15
UGC 7342 108 Lick 2.0 1800 2013 Jan 16

Download table as:  ASCIITypeset image

For both sets of spectra, absolute flux calibration used observations of standard stars, within the scaling errors set by changes in atmospheric seeing. These data, along with the long-slit spectra described by Keel et al. (2012a), provide information on the kinematics and chemical abundances of the gas and constrain ionization mechanisms in various regions (as described in Section 5, and C. J. Lintott et al. 2015, in preparation).

Some of the Lick spectra were targeted at possible outflow structures near the nuclei, and lie near the minor axes of the extended emission. The Lick PA 154° spectrum includes the projected neighboring galaxy SDSS J121820.18+291447.7, which proved to be a background object at z = 0.0674. For UGC 11185, the BTA spectrum resolved three velocity systems in the emission lines on the eastern side of the galaxy. The velocity results from these spectra are summarized in Figure 2, combining the redshift values with intensity slices for orientation with galaxy and emission-cloud features.

Figure 2.

Figure 2. Velocity slices from the long-slit spectra in Table 2, derived from Gaussian fits to [O iii] $\lambda 5007$ and shown as cz in the heliocentric frame. For spatial comparison, traces of continuum and [O iii] flux along the slit are shown at the bottom of each panel, as in Figure 12 of (Keel et al. 2012a). Their intensities are arbitrarily scaled for visibility, with [O iii] as the thicker lines. The offset spectrum of UGC 7342 does not have an associated continuum profile; at its distance from the core the continuum is essentially undetected. The spectrum of UGC 7342 along PA 154° crosses the background galaxy at z = 0.0674, which appears as a continuum enhancement at the indicated location. UGC 11185 shows three distinct velocity systems; the gap in the most extensive (low-velocity) system may be due to blending with the locally stronger intermediate component. The [O iii] profiles are shown separately for these components (at different intensity scales for visibility).

Standard image High-resolution image

2.3.1. Abundance Estimates

We estimate abundances in the extended clouds using the model grids from Storchi-Bergmann et al. (1998), which are computed for a typical AGN continuum shape and solar N/O ratio. They found that [N ii]/Hα is strongly abundance-dependent when either [O iii]/Hβ (we denote this combination Abund1) or [O ii]/[O iii] (for Abund2) is used as an excitation indicator. We employ both of these sets of line ratios (wherever our data cover the relevant lines) on both the spectra reported here and those in Keel et al. (2012a). (In comparisons, the Storchi-Bergmann et al. 1998 ratios include both members of the [O ii], [O iii], and [N ii] line pairs.) Their models express the results in O/H; our results are listed in Table 3. From the spectra in Table 2, we consider nuclei and distinct clouds separately, divided where the [O iii] emission shows local minima, giving up to five independent measurements of the EELRs in each galaxy. We include as well the abundances estimated from line ratios in Keel et al. (2012a) for nuclei and clouds, when red data covering [N ii] are available (including NGC 5252 from Acosta-Pulido et al. 1996, leaving only SDSS 1510+07 missing), as well as comparison values from data in Lintott et al. (2009) for IC 2497 and Hanny's Voorwerp. In deriving the line ratios, we correct for underlying Hα and Hβ absorption in the stellar popuations using the modeled equivalent widths from Keel (1983) as functions of continuum color.

Table 3.  Emission-line Ratiosa and Abundance Estimatesb

Galaxy Location [O iii] $\lambda 5007$ Fluxc [O ii]/[O iii] [O iii]/Hβ [N ii]/Hα [S II]/Hα [O I]/Hα Abund1 Abund2
IC 2497 Nucleus 21 2.16 3.6 1.53 0.54 8.69 8.72
Lintott et al. 2009 Voorwerp 496 0.110 13.97 0.22 0.165 0.028 8.51 8.68
 
Mkn 1498 Nucleus 490. 0.05 4.95 1.33 0.22 0.07 8.67 9.30
2012 data cloud 34. 0.15 16.35 0.19 0.20 0.06 8.57 8.56
NGC 5252 Nucleus 0.43 7.39 1.22 0.89 0.66 8.71 8.80
Acosta-Pulido et al. (1996) 2–5'' SSE 0.20 12.9 0.89 0.66 0.16 8.70 8.90
  6–11'' SSE 0.24 14.00 0.37 0.44 0.12 8.51 8.63
  12–17'' SSE 0.26 14.61 0.67 0.40 0.12 8.66 8.77
  5–13'' NNW 0.35 12.67 0.87 0.52 0.17 8.73 8.78
  17–26'' NNW 0.29 4.92 0.54 0.51 8.47 8.71
NGC 5972 Nucleus 45.6 10.2 0.94 0.60 0.075 8.69
BTA 8–22'' N 88.5 10.87 1.01 0.59 0.119 8.69
  46–55'' N 3.0 13.38 0.77 8.70
  2–11'' S 173 9.28 0.99 0.55 0.095 8.66
  16–22'' S 33.8 12.77 0.79 0.47 0.052 8.71
  38–54'' S 6.31 22.5 0.95 0.95 8.85:
NGC 5972 Nucleus 190. 0.22 3.40 0.92 0.49 0.07 8.57 8.90
2012 data cloud 95. 0.47 11.3 1.04 0.61 0.12 8.72 8.77
SDSS 2201 Nucleus 28.2 2.43 2.67 1.34 0.83 0.133 8.67 8.68:
BTA 14–20'' S 170 0.31 9.04 0.42 0.27 0.049 8.52 8.62
  10–19'' N 186 0.27 12.22 0.38 0.26 0.048 8.52 8.60
SDSS 2201 Nucleus 17. 1.24 7.28 1.61 1.02 0.17 8.77 8.62
2012 data cloud 78. 0.16 13.4 0.53 0.35 0.07 8.59 8.80
Teacup Nucleus 1300. 0.23 10.57 0.31 0.15 0.15 8.44 8.57
2012 data cloud 210. 0.17 10.73 0.45 0.36 0.11 8.52 8.73
UGC 7342 Nucleus 169. 0.23 10.85 0.73 0.62 0.083 8.63 8.82
BTA 8–16'' NW 50.0 0.174 14.45 0.26 0.27 0.055 8.54 8.58
  18–26'' NW 32.9 0.22 15.91 0.45 0.31 0.108 8.68 8.68
  11–21'' SE 47.6 0.157 12.30 0.25 0.29 0.039 8.48 8.58
  23–30'' SE 26.9 0.173 13.6 0.26 0.31 0.069 8.48 8.60
UGC 7342 E filament 66.0 0.20 14.39 0.36 0.34 0.084 8.56 8.64
UGC 7342 Nucleus 140. 0.73 10.95 0.73 0.69 0.085 8.64 8.63
Lick PA154 3–10'' NW 64.1 15.54 0.29 0.34 0.070 8.60
  10–23'' NW 35.1 0.37 18.09 0.60 0.44 8.78 8.61
  8–14'' SE 22.0 0.51 20.17 0.46 0.73 8.78 8.54
  15–22'' SE 12.3 16.04 0.24 8.48
UGC 7342 Nucleus 617. 0.15 10.23 0.72 0.70 0.084 8.61 8.82
Lick PA108 17–25'' W 57.8 0.15 9.51 0.43 0.23 8.50 8.57
  17–33'' E 216. 11.09 0.46 0.45 0.104 8.56
UGC 11185 Nucleus 1370 0.20 9.32 1.66 0.80 0.170 8.78 9.10
BTA 4–7'' E 212 0.40 9.56 1.45 0.75 8.77 9.05
  5–14'' W 27.0 5.95 1.41 0.83 0.115 8.72
UGC 11185 Nucleus 380. 0.28 11.82 1.71 0.88 0.19 8.89 9.10
2012 data cloud 51. 0.19 10.78 1.46 0.72 0.07 8.80 9.10

Notes.

aRatios involving [O ii] use sums of the $\lambda \lambda 3726,3729$ lines. Ratios using [O iii] give the sum of the $\lambda \lambda 4959,5007$ lines, and those with [N ii] sum the $\lambda \lambda 6548,6583$ lens for consistency with the calculations from Storchi-Bergmann et al. (1998). bAbundance estimates are 12+log (O/H), following the model grids presented by Storchi-Bergmann et al. (1998) using [O iii]/Hβ versus [N ii]/Hα (Abund1) and [O ii]/[O iii] versus [N ii]/Hα (Abund2). c[O iii] $\lambda 5007$ flux within the slit from the indicated region is given in units of 10−16 erg cm−2 s−1.

Download table as:  ASCIITypeset image

As might be expected from the broader density range present in the nuclei, the two abundance indicators agree more closely for the low-density clouds than in the nuclei. The scatter between these values for individual locations may be described by $\sigma =0.19$ dex for the nuclei and $\sigma =0.13$ dex for the extended clouds, with a systematic offset such that values using Abund2 are larger by 0.16 dex for nuclei and 0.09 for the clouds. We thus compare clouds and nuclei for each galaxy using single methods.

In interpreting the results, the internal consistency allows firmer statements about relative abundances than comparison with the solar value, in view of the uncertainties in the solar oxygen abundance. Using the solar value $12+{\rm log} [{\rm O}/{\rm H}]{\mkern 1mu} =8.76\pm 0.05$ following Caffau et al. (2008), most systems have clouds with subsolar metallicity, and the clouds are systematically lower in O/H than the nuclei. The mean values are $12+{\rm log} [{\rm O}/{\rm H}]=8.62$ from Abund1 ($\sigma =0.11$) and 8.70 ($\sigma =0.15$) from Abund2 for the extended clouds, versus 8.67, $\sigma =0.11$ and 8.84, $\sigma =0.22$ for the nuclei. (The samples are not identical for the two abundance indicators since we do not always have [O ii] data.) On this scale, several of the galaxies have extended clouds with abundances as low as 0.5 solar. In such cases as NGC 5252, where the data include multiple abundance estimates and we see systematic changes with radius, we may speculate about mixtures of nuclear outflow and more pristine, recently acquired gas in tidal features.

The uncertainty in these abundance estimates is almost certainly dominated by such details as ranges of density, which introduce ranges of ionization parameters and hence line ratios beyond the model assumptions. Within the low-density extended clouds, each indicator is internally consistent; for galaxies where our spectra sample distinct multiple emission stuctures (often on opposite sides of the galaxies), the scatter about the mean for each system is $\sigma =0.07$ for Abund1 and $\sigma =0.04$ for Abund2. At this level of confidence, these results show that some systems are significantly subsolar throughout (SDSS 2201+11, UGC 7342, IC 2497/Hanny's Voorwerp) while others are closer to solar abundances (Mkn 1498, NGC 5972, the Teacup, UGC 11185). This distinction may reflect different contributions from gas originally associated with the host galaxies versus external sources such as tidal disruption of low-mass, low-metallicity galaxies.

2.4. Velocity Mapping

For NGC 5972, SDSS 2201+11, UGC 7342, and UGC 11185, we obtained full velocity fields from [O iii] emission using the SCORPIO-2 system at the 6 m Large Altazimuth Telescope (BTA) in its scanning Fabry–Perot interferometer mode. The choice of targets was set both by weather and the availability of appropriate blocking filters to isolate [O iii] $\lambda 5007$ at the galaxy redshift. (UGC 7342 was observed in the $\lambda 4959$ line to use an available filter). In addition, we re-reduced the similar [O iii] data for NGC 5252 originally presented by Moiseev et al. (2002), Afanasiev et al. (2007), and Afanasiev et al. (2007). The observational parameters are listed in Table 4.

Table 4.  Fabry–Perot Velocity Maps

Galaxy Spectral FWHM (km s−1) Seeing FWHM ('') Exposure times (s) UT Date
NGC 5252 67 2.2 450 × 24 1998 Feb 24/25
NGC 5972 140 2.2 180 × 30 2014 May 21/22
SDSS 2201+11 42 2.0 300 × 30 2012 Aug 23/24
UGC 7342 137 2.8 240 × 30 2014 Apr 0/1
UGC 11185 140 2.7 210 × 30 2014 Oct 20/21

Download table as:  ASCIITypeset image

Data reduction followed the procedures in Moiseev & Egorov (2008). Voigt-profile fits of the data cube as sliced in the spectral direction were used to generate spatial distributions of continuum intensity, continuum-subtracted line intensity, mean radial velocity, and velocity dispersion. The Fabry–Perot results are summarized graphically in Figure 3, with the HST [O iii] images aligned for comparison. In UGC 11185, where the long-slit data show multiple [O iii] emission components, we have examined the data cube to investigate its two-dimensional extent on the sky; the component to the west of the nucleus at $cz\approx 12700$ km s−1 is associated with the roughly conical structure just east of the nucleus in the HST images.

Figure 3.
Standard image High-resolution image
Figure 3.

Figure 3. Summary of BTA Fabry–Perot observations, showing line images from the HST images and as derived from the Fabry–Perot data cubes, radial-velocity centroid, and velocity dispersion from Voigt-profile fitting, corrected for the instrumental resolution in each case. In addition to the centroid velocities shown here, NGC 5972 has regions with asymmetric or split line profiles to the north and south of the nucleus. The western component of UGC 11185 falls outside the monochromatic field of the ACS ramp fiters and is thus not well shown in the HST [O iii] image.

Standard image High-resolution image

Harrison et al. (2015) have recently presented integral-field spectroscopy of the Teacup system, giving velocity mapping of the entire emission region. They show that the large emission-line loop coincides with a radio-continuum structure, and identify a kpc-scale outflow region near the AGN itself. We incorporate these results in our later discussion.

2.5. Ground-based Imaging

The SCORPIO-2 system at the BTA was also used in its direct-imaging mode for V- and R-band images of SDSS 2201+11 and UGC 11185. At a pixel scale of 0farcs 357 pixel−1 covering a field of 6farcm1 square, these have excellent surface brightness sensitivity, tracing tidal features beyond the sensitivity of the HST intermediate-band data. Especially for SDSS 2201+11, the V images also show faint extended parts of the [O iii] clouds.

For NGC 5972, where both the emission-line regions and stellar tidal tails extend beyond the HST field, we supplement those data with images from the KPNO 2.1 m telescope, using a V filter for the continuum and a filter nominally centered at 5142 Å with FWHM 174 Å, for redshifted [O iii]. A TI CCD gave 0farcs 30 pixels. Total exposure times were 90 minutes in the [O iii] filter and 30 minutes in the V continuum.

To seek tidal features around SDSS 1510+07 beyond the SDSS surface brightness limit, we initially collected a total exposure of 80 minutes in R with the SARA remote 0.6 m telescope at Cerro Tololo, and followed up this detection using a total exposure of 5 hr in two sessions with the SARA 1 m telescope at Kitt Peak. Surface brightness calibration used standard stars from Landolt (1992).

As part of our examination of the interaction history of the EELR objects without energy-budget shortfalls, we use archival and new BTA images of Mkn 78 in the continuum (filters SED607 and FN641). The new FN641 exposure was obtained in 2014 October to test the possibility of faint reflections in the rich stellar field toward Mkn 78.

3. HOST-GALAXY STRUCTURE

Figure 4 shows the HST continuum images of the galaxies, displayed to span a wide range in surface brightness. Several show obvious tidal structures in these images, and all show evidence for disturbance via tidal tails or disrupted dust structure. We discuss these in turn and compare to the analogous set of AGNs with EELRs but energetics consistent with obscuration rather than fading.

Figure 4.

Figure 4. Host galaxy images in the HST F621M filter, shown in an offset log scaling, with brightness and contrast independently scaled for each. NGC 5252 is shown from a WFPC2 image in the F588M band. The Mkn 1498 and Teacup images have been corrected for emission-line contamination. Each image is shown with north to the left and east at the bottom; as in Figure 1, arrows are 5'' long and point north.

Standard image High-resolution image

3.1. Merger and Interaction Signatures

Mkn 1498 shows rich structure not resolved in earlier images. There are two broad spiral bands, which persist in residuals after both ellipse fitting (Section 2.2.1) and using a model with a fixed orientation (which is not influenced by these features). These may be the aftermath of a merging event.

NGC 5252 shows a symmetric, edge-on stellar disk in a luminous bulge. The evidence for past interaction here consists of the small-scale half-spiral of dust seen near the nucleus and the kinematically decoupled stars and gas (Morse et al. 1998).

SDSS 1510+07 has a small dust patch within the bar, showing the south to be the near side. This is detected in both continuum bands from independent fitting. This dust lane is closest to the nucleus on its western end. There is some evidence in the F621M image for a symmetric counterpart on the far side, with its deepest pixels showing a flux deficit of $\approx 10$% rather than $\approx 28$% on the near side. Our ground-based images show a tail connecting the AGN host galaxy to another galaxy just to the northwest, SDSS J151001.83+074058.4, at a typical surface brightness of R = 25.7 mag arcsec−2, as well as a more diffuse asymmetric structure outside the disk to the southeast (Figure 5). This tail is at roughly right angles to the axes of the ionized clouds.

Figure 5.

Figure 5. Deep R image of SDSS 1510+07 from the SARA 1 m telescope, smoothed by a Gaussian of FWHM = 2farcs 3. A tail connecting to the companion galaxy northwest is detected here, and a more diffuse asymmetric extension on the SE side of the AGN host galaxy. The field is 199'' square with north at the top, with an offset log intensity mapping.

Standard image High-resolution image

Beyond the twisted dust lanes apparent in the HST data, the stellar envelope of SDSS 2201+11 shows asymmetry and outer loops indicating an aging merger (Figure 6).

Figure 6.

Figure 6. BTA R-band image of SDSS 2201+11 from a 5 minute exposure, showing outer tidal features. The emission clouds show weakly via the Hα+[N ii] lines.The field shown spans 96'' × 109'', with north at the top, displayed using an offset logarithmic intensity mapping.

Standard image High-resolution image

The Teacup host galaxy clearly has an asymmetric structure, well described as a tidal tail with a shell-like feature at radius 10.0'' on the opposite side. These structures resemble stages in the merger of a small dynamically cold system with a bulge-dominated one, as in the modeling of shell galaxies by Hernquist & Quinn (1988). This configuration is seen in their models of radial encounters early in the postmerger phase when only a single shell appears. The tidal tail distinguishes this origin for the shell or arc from the weak-interaction case modeled by Thomson (1991). Generally, Feldmann et al. (2008) find that such narrow tidal features around elliptical galaxies must result from acquisition of cold disk systems with mass ratios $\approx 0.1$. and may remain distinct for 1–2 Gyr.

In interpreting the tails of UGC 7342, the merger is more advanced than it first appears; our spectra show that the apparent companion to the southeast (SDSS J121820.18+291447.7, seen at the right edge of the frame in Figure 1) is a background star-forming system at z = 0.0674, far beyond UGC 7342 at z = 0.0477; it shows positive residuals in our [O iii] continuum-subtracted image only because its Hβ and [O iii] $\lambda 4959$ lines are in opposite wings of the ramp-filter transmission at this setting. This identification reduces the computed far-IR luminosity of UGC 7342, since the WISE images (Wright et al. 2010) show that the background system strongly dominates the flux at 22 μm. The IRAS and Akari limits for the blended combination (Keel et al. 2012a) must be reduced to account for this unrelated contribution. Tidal tails are seen to 37'' (35 kpc), with several distinct shell-like surface brightness steps around a 10'' radius (9.5 kpc).

All our faded-AGN candidates are found in interacting, merging, or postmerger systems. Since the fading timescales we derive are of the order of the light-travel times across the systems, in contrast to the dynamical timescales which are orders of magnitude longer, we would expect the same to be true of the non-fading (that is, mostly obscured) systems in the parent sample of AGNs with large-scale ionized clouds from Keel et al. (2012a). Of the 12 obscured-AGN host galaxies in this parent sample, only one has no obvious source for extraplanar gas (merging, postmerger, interacting, or ram-pressure stripping in Virgo in the case of NGC 4388). This object is Mkn 78, in an elliptical host system.

This prompts us to re-examine Mkn 78 to seek potential subtle tidal features. Using pre-refurbishment HST images, Whittle & Wilson (2004) report a dust lane crossing the inner $\approx 1$'' at a skew angle to the major axis, which they suggest may result from an advanced accretion event, one old enough to leave little or no trace in the outer isophotes. For more insight into the structure of Mkn 78, we analyze the BTA images in BVR and red-continuum bands from Smirnova et al. (2010) and a new exposure in the FN641 filter centered near 6410 Å. The deep image in the SED607 filter, which lies between [O iii] and [O i] $\lambda 6300$ lines, is especially important, rejecting the bright emission-line clouds. As Smirnova et al. note, the images show no asymmetric structure or shells to deep levels, to a surface brightness $R\approx 26$ mag arcsec−2. A resolved companion or tidal plume appears in both these data sets and, in hindsight, the SDSS images, about 12'' W of the nucleus. We isolate it by fitting elliptical isophotes to the galaxy after masking this region and contaminating objects and examine the residuals after subtracting a model based on these isophotes. These also reveal an asymmetric feature oriented E–W, which could be interpreted as a small-scale stellar disk seen edge-on, and a light deficit crossing the nucleus near PA 35°, as well as the diffuse excess. This excess light, detected at 8–16'' from the core, is asymmetric with a tail to the north and involves roughly $6\times {{10}^{-4}}$ of the galaxy's red continuum light. A tidal origin is plausible, although not fully demonstrated; this would mean that all EELRs in our parent sample are in interacting, merging, or ram-pressure stripped systems.

3.2. Surface Brightness Profiles

Photometric profiles in most of these galaxies (derived as in Section 2.2.1) show them to be strongly bulge-dominated, in some cases with no significant disk component. The averaged brightness profile in the Teacup galaxy is close to a de Vaucouleurs ${{r}^{1/4}}$ form over nearly 10 magnitudes in surface brightness and a range of 3–370 pixels (0.15–18'') in radius.

To quantify the bulge properties and derive detailed information or limits on disk components, we used the GALFIT software (Peng et al. 2002, 2010), version 3.0.5. Dust masks were produced in a semi-automated way, starting with the dust maps in Figure 7, using attenuation thresholds depending on the local S/N and manually selecting the threshold appropriate to various parts of the image. Interloper star and galaxy images were simply masked from the fits. Central PSF components to accommodate an unresolved AGN were allowed, but were significant only in Mkn 1498, NGC 5252, and the Teacup AGN. For galaxies observed with WFC3, the F763M image was fit to reduce the effects of dust and young stellar populations. For NGC 5252, we fit the F606W WFPC2 image. The Sérsic indices of "bulge" components from these fits are listed in Table 5. Here, n is the single-component Sérsic index and nd is the index when an additional component fixed to be an exponential disk is included. Disk/bulge and AGN/bulge ratios are in intensity. All n values fall into the domain occupied by classical bulges, rather than pseudobulges or other structures (Fisher & Drory 2008).

Figure 7.

Figure 7. Attenuation images of target systems, using a common linear intensity display running from 0 to 1.5 in ratio of image to model. Ordered as in Table 1, the top row includes (left to right) Mkn 1498, NGC 5252, NGC 5972, and SDSS 1510+07; the bottom row has SDSS 2201+11, the Teacup, UGC 7342, and UGC 11185. The field is 18farcs 0 E–W, except for SDSS 1510+07 and UGC 7342 at 9farcs 05, and UGC 11185 at 15farcs 0. North is at the top and east to the left.

Standard image High-resolution image

Table 5.  Fitted Sérsic Indices

Galaxy n nd Disk/Bulge AGN/Bulge
Mkn 1498 6.04 6.14 0.41 0.06
NGC 5252 3.55 5.18 0.13 0.01
NGC 5972 4.26 3.71 0.51 ...
SDSS 1510+07 4.6 2.41 1.58
SDSS 2201+11 3.51 7.12 0.44
Teacup 5.94 5.65 0.25 0.02
UGC 7342 3.95 4.61 0.09
UGC 11185 7.31 8.87 0.02

Download table as:  ASCIITypeset image

In SDSS 1510+07, the stellar disk was also retrieved by GALFIT, although simultaneous fits to the bar and disk are poor with this set of components. In some of the highly disturbed systems, the physical significance of a formal exponential-disk fit is not clear.

Some potential explanations for the fading-AGN phenomenon involve interactions between members of binary black holes; several sets of simulations show the individual accretion disks having rapidly fluctuating luminosities (e.g., Bogdanović et al. 2008). One signature of such a process would be the "scouring" of stars in the core (Hopkins & Hernquist 2010; Rusli et al. 2013) In this sample, appropriately accurate measurements of the innermost light profile proved infeasible because of dust lanes or excess light from the AGN itself.

3.3. Dust Structures

Dust can show whether or not the cool ISM in a galaxy is relaxed and, since its attenuation effects strongly favor detection on the near side of a galaxy, can also provide some three-dimensional information on its location. We examine dust attenuation primarily through the ratio of continuum images to smooth elliptical models generated from the surface brightness profiles (Section 2.2.1). This procedure may produce artifacts, especially near the cores or dust lanes, so we also compare the results to ratio maps between the two continuum wavelengths. These color maps have lower S/Ns for detection, since the wavelength baseline is small, but do not suffer from the modeling uncertainties that can occur with galaxy intensity.

These attenuation maps in the bluer continuum image (WFC3 F621M or WFPC2 F588N) are shown in Figure 7, to a common scale in intensity ratio. Complex dust lanes and patches are ubiquitous, often with distributions suggesting the aftermath of strong interactions or mergers. Only in SDSS 1510+07 do we fail to detect such disorganized dust structures.

In Mkn 1498, we see an irregular and asymmetric dust lane at a sharp angle to the emission structures, detected to 5farcs 6 (6.2 kpc) toward the southeast.

In NGC 5252, we find a small half-spiral of dust seen to the northwest of the nucleus, strongly misaligned with the stellar disk, as noted by Tsvetanov et al. (1996). This disk has an orientation which could be circular in the plane defined by the kinematics of the extended ionized gas filaments (Morse et al. 1998) and, at modest resolution, H i (Prieto & Freudling 1996).

The dust in both NGC 5972 and SDSS 2201+11 crosses the central region multiple times on opposite sides of their nuclei. In both cases, this may represent a side view of a differentially precessing disk (Steiman-Cameron et al. 1992); such a disk could continue to the ionized clouds in SDSS 2201+11, where a dust lane connects to the northern cloud (and in order to show distinct absorption, this feature must lie on the near side of the galaxy). Our long-slit and Fabry–Perot data show the northern side of the system to be receding, helping to specify the geometric conditions of our view.

In the Teacup system, the ratio images show detail in the complex dust features seen even in the unprocessed data, especially to the southeast of the nucleus with a connection spiraling outward to projected distance 1farcs 9 (3.3 kpc). The deepest attenuation in this image, seen 0farcs 3 W of the nucleus, has 0.41 of the model intensity in the green band and 0.58 in the red. The inner part of the dust band runs roughly north - south, at about 45° to the greatest elongation of the ionized gas.

UGC 7342 shows a warped dust lane near the core, which may hide the AGN itself from direct optical view.

3.3.1. Precessing Disks

Several of these systems show dust lanes resembling the complex crossing geometry in NGC 4753 (Sandage 1961), which prompted us to consider them in light of the model for differentially precessing disks of material developed by Steiman-Cameron & Durisen (1988) and applied to NGC 4753 by Steiman-Cameron et al. (1992). Using their scale-free logarithmic potential model with a flat rotation curve, a ring of gas in an orbit of radius r at circular velocity vc and inclination i to the principal plane of a flattened potential well, will precess at and angular rate of

Equation (1)

where η depends on the shape of a (generally triaxial) potential: for major and minor axes $a,c$,

Equation (2)

In applying the model, we assume that the potential shape is oblate with the same axial ratio as the starlight distribution (that is, that we see the mass distribution edge-on), giving $\eta =0.113-0.128$. Since Steiman-Cameron et al. (1992) find that the damping term in i corresponding to an effective viscosity is small, we ignore this effect.

We implemented their model for NGC 5972 and SDSS 2201+11 to explore the ranges of geometries and times that could match the dust geometry we observe; these are the systems with well-resolved structure including features expected in the model. The starting point is an idealized introduction of a flat annulus of interstellar matter to the host galaxy at t = 0 at an inclination to the equatorial plane estimated from the crossing angles of dust lanes; 33° for NGC 5972 and 21° for SDSS 2201+11. The circular velocity is taken from averaging our long-slit data; 190 km s−1 for NGC 5972 and 255 km s−1 for SDSS 2201+11. This fixes the radial scale, leaving the age and viewing geometry as free parameters. Seeing discrete crossing dust lanes, rather than a broad projected distribution as in Centaurus A, would mean that the dust is confined to narrow radial zones, affording an additional freedom in matching predictions (and making our efforts more of an existence proof than a detailed fit) since there is little additional information on their line-of-sight locations. The model must be old enough for dust lanes to cross at sharp angles as we observe, a limit which in turn depends on the radial distance range between the dust features. For NGC 5972, dust is detected in our ratio images to about 9'' or 5.5 kpc. Since dust detection is strongly favored in front of the bright inner bulge, we consider dust in rings from 3 to 10 kpc in radius. From a typical direction along the galaxy's equatorial plane, this leads to the crossing of the inner rings as soon as 0.9 Gyr, and up to 3 Gyr for the outer rings (Figure 8). For dust absorption, we examine only the parts of these precessing rings on the near side, while for comparison with the ionized gas, we include the far side as well. In NGC 5972, the inner ionized gas regions might be interpreted as the intersection of ionization cones with such rings of material. Supporting this idea, there is no emission associated with the dust lanes in front of the core, with some emission regions plausibly connected with one such lane when projected $\approx 3$ kpc from the nucleus.

Figure 8.

Figure 8. Schematic precessing-disk model for dust in NGC 5972. These are edge-on views at a time 1.35 Gyr after introduction of an annulus of material inclined by $i=33{}^\circ $ to the principal plane of the host galaxy, a timestep which captures three important dust locations in the images. In comparison to images of NGC 5972, north would be to the left. The upper panels show only the material on the front side, which might be visible through attenuation of starlight, while the lower panels show the entire configuration. The left panels show rings of tracers spaced radially every 0.1 kpc, while the right panels have 0.5 kpc spacing to indicate the role of dust in distinct lanes. Scales show offsets in kpc. The upper right panel, with spaced rings shown where they cross in front of the system, is superimposed on the dust-attenuation image from Figure 7.

Standard image High-resolution image

In SDSS 2201+11, there is a direct connection between the gas and dust distributions; a dust lane leads from the main dust disk into the north-northwest emission-line cloud (which is thereby on the near side, and receding as shown by both the longslit and Fabry–Perot observations). In modeling this system, we take an axial ratio of 0.75 for the host starlight and inclination $i=21{}^\circ $. Dust is detected to a radius 16'' or 9.5 kpc; we consider annuli from 3 to 10 kpc. In this system we can evaluate the circular velocity from the full 2D BTA velocity map, adopting 255 km s−1. Crossing of dust lanes near the core starts around 1.4 Gyr, occurring later for features farther out. A sample model is shown in Figure 9.

Figure 9.

Figure 9. Schematic precessing-disk model for dust in SDSS 2201+11. These are edge-on views at a time 2.5 Gyr after introduction of an annulus of material inclined by $i=21{}^\circ $ to the principal plane of the host galaxy, a timestep which captures three important dust locations in the images. In comparison to images of SDSS 2201+11, north would be to the left. The upper panels show only the material on the front side, which might be visible through attenuation of starlight, while the lower panels show the entire configuration. The left panels show rings of tracers spaced radially every 0.1 kpc, while the right panels have 0.5 kpc spacing to indicate the role of dust in distinct lanes. Scales show offsets in kpc. The upper right panel, with spaced rings shown where they cross in front of the system, is superimposed on the dust-attenuation image from Figure 7.

Standard image High-resolution image

If the ionized gas forms part of the outermost ring matching this model, its actual cone angle for illumination from the nucleus is substantially larger than the 23° we see in projection, since we see the warped disk at a sharp angle.

4. STAR CLUSTERS

The evidence for a region of outflow-triggered star formation in Hanny's Voorwerp (Keel et al. 2012b), and more generally the signs of a role for AGN feedback in regulating star formation (e.g., Fabian 2012), give impetus to search for young stars and clusters associated with these AGNs and surrounding gas clouds. These may be detected either directly, when in luminous blue clusters, or secondhand, as localized regions with [O iii] /Hα much lower than in the surrounding AGN-ionized gas; in Hanny's Voorwerp, both signatures occur in the same places. We require cluster candidates to be detected in both of our continuum bands, to reject remaining cosmic-ray artifacts and reduce confusion with fine structure in dust lanes. Single clusters will be essentially pure PSF images at these redshifts, although blending of adjacent clusters (as well as incomplete correction of CTE effects in the images) leads us to consider slightly resolved objects as well (with measurements described in Section 2.2.2). None of the target galaxies have a rich cluster population. Our color information is somewhat limited in precision by the modest wavelength baseline between our continuum filters, which were chosen primarily for the best match to emission-line bands. Except in the inner regions of the galaxies where Poisson noise increases, the $5\sigma $ detection threshold is about ${{m}_{621}}=25.5$, ${{m}_{763}}=24.7$ for measurements within a 0farcs 4 radius. For faint clusters or those with close neighbors we use r = 0farcs 2 and a point-source correction to match the zero points given for 0farcs 4 (http://www.stsci.edu/hst/wfc3/phot_zp_lbn). To reduce contamination from foreground stars or background galaxies we accept only cluster candidates with $(R-I)\lt 0.9$ and a computed absolute magnitude M621 fainter than -15. No such objects are found in NGC 5252 and SDSS 1510+07. Photometric values for candidate clusters are given in Table 6.

Table 6.  Candidate Star Clusters

Galaxy ID $\alpha ,\delta $ (2000) m763 m632 ${{(}_{R-I)0}}$
Mkn 1498 1 16:28:03.675+51:46:26.43 25.26 ± 0.16 25.33 ± 0.10 0.08 ± 0.18
  2 16:28:03.484+51:46:29.17 25.10 ± 0.15 25.03 ± 0.08 −0.07 ± 0.17
  3 16:28:02.826+51:46:28.25 24.75 ± 0.10 25.33 ± 0.09 0.63 ± 0.15
NGC 5972 1 15:38:54.374+17:01:29.88 24.46 ± 0.11 24.61 ± 0.08 0.16 ± 0.14
  2 15:38:53.830+17:01:37.16 24.49 ± 0.11 25.05 ± 0.12 0.67 ± 0.16
  3 15:38:53.979+17:01:40.28 24.81 ± 0.15 25.10 ± 0.12 0.32 ± 0.19
  4 15:38:53.188,+17:01:40.49 24.46 ± 0.07 24.74 ± 0.07 0.30 ± 0.11
  5 15:38:53.315+17:01:34.56 24.87 ± 0.11 25.00 ± 0.07 0.14 ± 0.13
  6 15:38:52.760+17:01:39.35 24.22 ± 0.06 24.64 ± 0.05 0.45 ± 0.10
  7 15:38:52.876+17:01:43.77 24.94 ± 0.12 25.33 ± 0.10 0.42 ± 0.16
  8 15:38:54.893+17:01:37.06 23.43 ± 0.03 24.58 ± 0.05 0.17 ± 0.07
  9 15:38:54.819+17:01:20.50 23.85 ± 0.04 24.26 ± 0.04 0.45 ± 0.07
  10 15:38:54.518+17:01:09.83 25.69 ± 0.23 26.42 ± 0.26 0.79 ± 0.35
Teacup 1 14:30:29.910+13:39:14.43 23.47 ± 0.04 23.81 ± 0.03 0.37 ± 0.05
  2 14:30:29.857+13:39:09.02 23.57 ± 0.07 23.84 ± 0.05 0.30 ± 0.09
  3 14:30:29.907+13:39:09.95 24.06 ± 0.05 24.44 ± 0.08 0.42 ± 0.12
  4 14:30:30.137+13:39:11.51 22.93 ± 0.02 23.45 ± 0.02 0.57 ± 0.03
  5 14:30:30.131+13:39:11.31 23.47 ± 0.03 23.66 ± 0.03 0.21 ± 0.05
  6 14:30:30.155+13:39:11.44 24.27 ± 0.06 24.54 ± 0.06 0.30 ± 0.09
  7 14:30:30.151+13:39:11.90 23.48 ± 0.04 23.86 ± 0.03 0.42 ± 0.05
UGC 7342 1 12:18:16.507+29:15:25.87 25.12 ± 0.15 25.35 ± 0.09 0.25 ± 0.17
  2 12:18:17.298+29:15:20.77 24.05 ± 0.06 24.16 ± 0.04 0.12 ± 0.07
  3 12:18:17.517+29:15:19.31 24.87 ± 0.12 25.08 ± 0.08 0.21 ± 0.14
  4 12:18:17.592+29:15:17.80 24.55 ± 0.09 25.12 ± 0.08 0.63 ± 0.12
  5 12:18:17.859+29:15:15.93 24.01 ± 0.05 24.45 ± 0.05 0.48 ± 0.07
  6 12:18:19.215+29:15:27.95 24.60 ± 0.08 25.14 ± 0.07 0.59 ± 0.13
  7 12:18:19.224+29:15:28.41 24.04 ± 0.05 24.82 ± 0.06 0.78 ± 0.09
  8 12:18:19.335+29:15:28.77 25.12 ± 0.13 25.35 ± 0.09 0.25 ± 0.16
  9 12:18:20.935+29:15:31.25 24.50 ± 0.08 24.72 ± 0.05 0.24 ± 0.11
  10 12:18:20.688+29:15:06.16 25.06 ± 0.13 25.26 ± 0.08 0.22 ± 0.15
  11 12:18:20.298+29:15:02.67 24.89 ± 0.11 25.46 ± 0.10 0.63 ± 0.15
  12 12:18:21.752+29:14:53.63 24.65 ± 0.09 24.71 ± 0.05 0.07 ± 0.12
  13 12:18:20.671+29:14:52.95 24.70 ± 0.09 24.87 ± 0.06 0.19 ± 0.13
  14 12:18:20.687+29:14:52.43 24.91 ± 0.11 25.21 ± 0.08 0.33 ± 0.14
  15 12:18:21.168+29:14:48.42 24.65 ± 0.09 24.73 ± 0.05 0.09 ± 0.12
  16 12:18:21.203+29:14:47.07 24.46 ± 0.07 24.74 ± 0.05 0.31 ± 0.10
  17 12:18:19.065+29:14:40.70 24.43 ± 0.17 25.12 ± 0.07 0.20 ± 0.18
  18 12:18:18.851+29:14:59.02 24.44 ± 0.07 24.74 ± 0.05 0.33 ± 0.10
  19 12:18:17.793+29:14:58.70 24.08 ± 0.05 24.68 ± 0.05 0.66 ± 0.07
  20 12:18:19.048+29:15:29.21 23.56 ± 0.03 23.84 ± 0.01 0.30 ± 0.03
  21 12:18:19.650+29:15:12.25 24.39 ± 0.16 25.03 ± 0.11 0.70 ± 0.19
  22 12:18:19.101+29:15:13.86 24.56 ± 0.16 24.70 ± 0.10 0.16 ± 0.18
  23 12:18:18.926+29:15:14.85 25.15 ± 0.18 25.19 ± 0.10 0.04 ± 0.20
UGC 11185 1 18:16:13.763+42:39:44.60 22.06 ± 0.02 22.64 ± 0.02 0.53 ± 0.03
  2 18:16:13.811+42:39:43.45 22.58 ± 0.03 22.84 ± 0.03 0.29 ± 0.03
  3 18:16:09.773+42:39:41.80 22.37 ± 0.03 23.00 ± 0.02 0.69 ± 0.04
  4 18:16:10.575+42:39:42.21 22.87 ± 0.04 23.44 ± 0.04 0.68 ± 0.06

Download table as:  ASCIITypeset image

None of these cluster candidates show an [O iii] /Hα signature as we saw in Hanny's Voorwerp. Most of them are found too far from strong line emission (and those which are close to emission filaments do not show the ratio reversal that indicates a change to ionization by the UV continuum of hot stars). We do see a broad association between regions of high [O iii] /Hα ("ionization cones") and the cluster candidates, which we might speculate indicates a long-term connection between the AGNs and cluster formation, as shown in Figure 10. The cone angles in a projection range from 57–68°, so a random distribution would have about 1/3 of the points within the cones. Among all cluster candidates, we see a small excess over this (20/43), although the significance is marginal and could be changed by small variations in the boundaries in UGC 7342 and NGC 5972. Any outflow-induced star formation within the emission regions in these galaxies must be a minor effect or have been operating over long times compared to the light-travel times sampled by the observed gas.

Figure 10.

Figure 10. Spatial distributions of star-cluster candidates. North is at the top in each panel; to show both stellar and gaseous structures, the base images here are the continuum plus [O iii] mixes used in Figure 1. Dashed lines mark approximate boundaries of ionization cones as inferred from the [O iii] /Hα ratio. The large circle in the Teacup system marks a likely background galaxy, which might otherwise be interpreted as $\approx 6$ blue clusters surrounding a redder object.

Standard image High-resolution image

Mkn 1498 has three candidate clusters, all to the southwest of the nucleus and roughly perpendicular to the axis of the extended ionized gas. All are far from the ionized gas.

NGC 5972 shows nine cluster candidates within the projected radius of 30'' (17 kpc), with no obvious preferred alignment. Their absolute magnitudes range from −10.2 to −11.3, with a color range centered near $(R-I)=0.30$. Four of these nine lie projected within the ionization cones

Cluster candidates exist in SDSS 2201+11, but are so close to the dust lanes that precise photometry is greatly limited by fluctuations in background and attenuation, so they are not shown in the montage or listed in Table 6.

In the Teacup host galaxy, we find several candidate clusters within the inner 5'', mostly concentrated in a clump about 3farcs 5 E of the nucleus or roughly aligned with the dust lane. We measure their fluxes after subtracting a smooth model for the galaxy and use local statistics to evaluate photometric errors. The prominent clump may in fact be a background galaxy, since its structure of a redder core surrounded by blue knots is reminiscent of many disk galaxies seen in the emitted UV; we measure total values with r = 1'' of ${{m}_{763}}=21.05\pm 0.03$, ${{m}_{632}}=21.61\pm 0.03$, ${{(}_{R-I)0}}=0.61\pm 0.05$. All of these candidates occur in regions of lower-ionization gas outside the ionization-cone boundaries. The two bright objects outside the dust lane have ${{M}_{R}}\approx -14$, ${{(}_{R-I)0}}\approx 0.30$.

In UGC 7342, we find 23 candidate clusters in our color and luminosity range, from ${{M}_{R}}=-13$ and fainter. Five of these are in the prominent tidal tail to the northwest, and three are projected in the inner 6'' (5.5 kpc). Sixteen are quite blue, near $R-I=0.30$. As with UGC 11185, the five cluster candidates in the tidal tail lie near the projected boundary of the large-scale ionization cone. The clusters here are more luminous than in the other systems, with three brighter than ${{M}_{R}}=-12$ and all of the detections brighter than ${{M}_{R}}=-11$. There is marginal evidence for a bimodal cluster population with peaks near ${{(}_{R-I)0}}=0.3$ and 0.75. Half of the candidates (12/23) are within the ionization cones, albeit in many cases also within stellar tidal features.

UGC 11185 is seen at a low galactic latitude, so the potential contamination by bright foreground stars makes the bright magnitude cut for candidates more important. This leaves four distinct cluster candidates within 30'' (24 kpc) of the nucleus, one of which has at least three subcomponents at the 0farcs 15 level. This object and a close neighbor lie on the inner edge of the tidal loop to the east of the nucleus; the other two lie almost opposite this complex on the west side of the system. All four objects are projected just within the angular span of the ionization cone seen closer to the AGN--two in a stellar tidal stream and two possibly associated with the companion spiral galaxy to the west-southwest. These are potential cases of outflow-induced star formation.

We evaluate the surface density of foreground and background objects, which could potentially masquerade as star clusters, using the outer parts of the WFC3 fields (except for NGC 5972, this used the CCD, not including the galaxy center). Performing similar photometry on objects which are unresolved or marginally resolved, excluding obvious background galaxies, we derived object counts in bins of R621 and $(R-I)$ for ${{R}_{621}}\lt 25.0$ and requiring detections on both filters (as for the cluster candidates). Based on the area covered by our target galaxies, the expected number of unrelated objects within our color range is as large as 1, only for NGC 5972 and UGC 7342, in the range ${{R}_{621}}=24-25$. Contamination does not significantly affect our statistics of cluster candidates. Most potentially contaminating objects are quite red; about 75% are in the range $(R-I)=0.7-2.5$.

None of these galaxies are rich in young clusters despite the strong evidence of interactions and mergers, and also the substantial supplies of gas and dust. We might speculate that AGN episodes could have previously quenched star formation, or that there is some strong selection in favor of bulge-dominated systems to have the giant ionized clouds. In the former case, dust lanes might not only be plausible locations for molecular clouds, but could play a role in shielding such clouds from radiative interference by the AGNs.

Our color cut is generous for the ages and reddenings of candidate clusters; in comparison, the Starburst 99 models (Leitherer et al. 1999) show simple stellar populations briefly becoming as red as $(R-I)=0.9$ at ages near $8\times {{10}^{6}}$ yr for high upper mass limits, as red supergiants are most important in the integrated light (so this is the time of peak near-IR luminosity), but otherwise remaining bluer than $(R-I)=0.6$. With many candidates far from the galaxy cores, we would expect only local reddening except near the dust lanes.

5. GAS KINEMATICS

For the large-scale emission structures, the main feature of the velocity fields in these systems is that most appear rotationally dominated, following ordinary rotation curve patterns with only modest local displacements or velocity dispersions. This stands in strong contrast to the EELRs around radio-loud quasars (Stockton et al. 2006; Fu & Stockton 2009), which show emission features over ranges of 1000–2000 km s−1 (possibly resulting from interactions with the jets). For the objects with Fabry–Perot velocity fields, we can extend this conclusion to the full systems rather than sampling along single slits.

To characterize local departures from ordered motion, we consider residuals from median-smoothed velocity fields in the BTA data. For each system we construct a smoothed velocity field by taking the median of values within 9 × 10'' (13 × 25 pixels) boxes, with the short axis along the overall velocity gradient. In comparison of the velocity and velocity-dispersion maps (Figure 3), we find that almost all values ${{\sigma }_{v}}\gt 100$ km s−1 occur in regions of low surface brightness or strong gradients, where both beam smearing and seeing-match issues are strongest.

NGC 5252 remains an outlier in that its velocity gradient is almost perpendicular to the ionization cone (as seen in the H i data of Prieto & Freudling 1996 and the optical spectra of Acosta-Pulido et al. 1996 and Morse et al. 1998). The velocity field in this case is not purely planar rotation, as seen from the inflections in its contours. The only locally high velocity dispersion is at the AGN itself.

The velocity field appears continuous in NGC 5972 into the detached emission regions outside the ACS narrowband images. The line profiles are asymmetric, possibly reflecting projected material at very different line-of-sight locations, within 3'' north and south of the nucleus, in the sense that additional components are present closer to systemic velocity than the main rotation pattern. This is not seen along the minor axis, so this velocity structure is not due to seeing effects. A similar effect occurs in the longslit data along PA 165°: the velocity shows a displacement toward the systemic value 11'' N and S of the nucleus, where the dust attenuation is strongest. This spectrum shows stellar absorption in the Mg+MgH region in the inner 6'' sufficient to show that the stellar rotation is in the same sense as we see in the gas.

In SDSS 2201+11, the southern emission-line cloud may be two projected in overlap differing in radial velocity by ∼50 km s−1. There are distinct velocity components separated by about 7'' east-west, with a local increase in ${{\sigma }_{v}}$ at their projected interface. The higher-velocity component corresponds to distinctly more diffuse material on the western side of this structure. The velocity contours in this system are not purely tangential; the cloud edges away from the main dust disk plane have smaller rotational velocity, which may be explained by a warped disk model. The long-slit rotation slice shows a nearly undisturbed, typical rotation curve with circular velocity near 235 km s−1. This spectrum shows stellar absorption in the Mg+MgH region in the inner 20'' sufficient to show that the stellar rotation is in the same sense as we see in the gas.

The emission-line mapping of the Teacup system reported by Harrison et al. (2015) shows a localized region of large blueshifts within 1''SW of the AGN, with projected velocity offset 760 km s−1. Their radio mapping shows continuum emission coincident with part of the large loop of emission to the west, suggesting that this feature also marks an outflow (perhaps entrained material around the result of a previous radio outburst).

In UGC 7342, there is an outlying region 100 km s−1 redward of the smooth velocity field seen 8farcs 5 northwest of the nucleus. The gas distribution in this system is perpendicular to the circumnuclear dust lanes, and so is the velocity gradient; this might be an extreme example of a warped and precessing disk, with damping significant in bringing the material toward the plane of the main galaxy potential.

Important orientation angles are compared in Table 7. The kinematic major axis, perpendicular to the kinematic line of nodes (LON), is often (but not always) near the major axis of the emission clouds.

Table 7.  Structural and Kinematic axes

Galaxy Gas Major Axis (°) Kinematic LON (°) Kinematic Major Axis (°) Dust Lanes (°)
NGC 5252 161 6 104
NGC 5972 169 96 6 169
UGC 7342 131 38 128 29
SDSS 2201+11 21 118 28 173

Download table as:  ASCIITypeset image

In UGC 11185, there are three [O iii] velocity components on the east side of the system. The higher-velocity extensive system is associated in the Fabry–Perot data with the conical or triangular region bright in the HST [O iii] image. This component dominates the velocity fits where it is brightest, and blending in this region may account for some of the velocity structure seen by Keel et al. (2012a); on both sides the lower-velocity "systemic" curve shows a smoother, albeit still disturbed and asymmetric, form. The higher-velocity emission could trace an outflow, or ionization of a complex tidal tail projected against the inner part of the galaxy. The high abundances in the UGC 11185 extended clouds may argue for an outflow interpretation.

The ionized gas distributions in the tidal tails broadly match the stellar tails seen in continuum images, but there are small-scale offsets that might in principle either show genuine offsets in gas and star distributions (which naturally results from the different radial and velocity distributions of these components in the original disk before tidal disturbance) or show us three-dimensional structure in the tidal debris weighted by the distribution of ionizing radiation for the gas. These differences are most clearly seen in the outer parts of NGC 5972 (mostly beyond the field of the HST emission-line images), as shown in Figure 11.

Figure 11.

Figure 11. Stellar and gas tails in NGC 5972, from data with the KPNO 2.1 m telescope. Contours show logarithmic intensity of a V continuum image, corrected for the [O iii] line emission, and the gray scale shows continuum-subtracted [O iii]. Especially in the southern (lower) tail, the line emission broadly follows the stellar tail, but has an offset center line. Residuals from bright stars appear at opposite corners of the view. The scale bar shows 1'.

Standard image High-resolution image

6. GAS FINE STRUCTURE

As was done for IC 2497 and Hanny's Voorwerp by Keel et al. (2012b), we compare the pixel-by-pixel surface brightness in Hα to the ionization parameter U to test for the structure in the ionized clouds beyond the resolution limit of HST data. We can map the observed [O iii] /Hα ratio to changes in U because these clouds are photoionized by a distant source, so that the continuum shape is fixed and we can isolate regions with a small range in (projected) radius from the AGN. Under these conditions, if a resolution element samples a smooth distribution of matter, the surface brightness in a recombination line will scale as ${\rm SB}\propto n_{e}^{2}$, and U will scale as $n_{e}^{-1}\propto {\rm S}{{{\rm B}}^{-1/2}}$. At these ionization and density levels, [O iii] /Hα scales in a nearly linear way with U (e.g., Ferland & Netzer 1983). Our test thus uses the relation between Hα SB and [O iii] /Hα for individual features in the ionized clouds; significant departures from the ${\rm S}{{{\rm B}}^{-1/2}}$ scaling indicate that the conditions above are not satisfied. In Hanny's Voorwerp, the departures were in the sense that brighter pixels had U decrease more slowly than predicted; a straightforward interpretation is that there is unresolved structure, so that regions of higher surface brightness have a larger number of unresolved clumps and filaments rather than higher mean density.

In these systems, this test gives mixed results. As shown in the samples in Figure 12, some regions (each including pixels at nearly the same projected distance from the nuclei) are reasonable fits to the expected behavior if the HST images resolve most of the structure in the ionized gas. Others show a much weaker dependence of ionization parameter on Hα surface brightness than expected, suggesting that additional fine-scale structure exists beyond the HST resolution limit. The latter cases include the inner filaments in UGC 7342, radial samples in the loop of the Teacup system, the outer parts of the radial clouds in Mk 1498, and the inner bright filaments in NGC 5972. In a very broad way, clouds at the largest radii, often with lower mean surface brightness, are more likely to follow the simple ${{n}_{e}}-U$ behavior for a smooth distribution of material. Small features cannot be evaluated with this test if there are too few pixels to establish the relationship.

Figure 12.

Figure 12. Examples of poor and good fits of the Hα surface brightness—ionization-parameter relation to the expectation if Hα surface brightness represents constant ne within each resolution element, shown by the schematic curve for each panel. Left column, poor fits indicating significant unresolved density structure. Regions in the right column have relations that are well fit by the simple expectations. Points shown are masked at $2\sigma $ error is Hα, and a 3 × 3 pixel median filter was used to reject residual bad pixels. Each panel is labelled with the number of pixels shown. For each data set, a 51 pixel running median is shown to illustrate the behavior of the distributions.

Standard image High-resolution image

In several of these galaxies, the ionized gas shows coherent filaments extending several kpc, sometimes in parallel sets. Good examples are seen in SDSS 1510+07, extending as much as 7farcs 8 (7.5 kpc). These are barely resolved, with widths around 0farcs 15 (140 pc). This combination of length and narrowness may indicate confinement; even at the thermal velocity for H (17 km s−1 rms) associated with a typical electron temperature 12,000 K, such narrow filaments would disperse much faster than the dynamical timescales $\approx {{10}^{8}}$ yr given by the rotation velocity and size of the filaments. One could picture the filaments constantly being regenerated by instabilities or being preserved by confinement, both of which may suggest magnetic fields which are dynamically important on local scales for the diffuse gas and ordered over kpc scales. Such fields could shape the geometries of the clouds on large scales, both in filaments and the twisted "corkscrew" forms which are suggested in several of these galaxies (such as UGC 7342).

Some of these systems have extended radio emission (notably NGC 5972, Mkn 1498, and on smaller scales the Teacup). In these, our conjecture about magnetic confinement might become testable through sufficiently sensitive and high-resolution Faraday rotation measurements. The nature of the emission clouds—externally photoionized sections of tidal debris—introduces the possibility that such fine structure (and possibly embedded magnetic fields) could be a common feature of gas in tidal features, currently observable only when we see optical emission lines from the ionized portions.

7. ORIGIN OF THE IONIZED GAS

Since the ionized clouds we examine are (by selection) more than 10 kpc from the AGN, they are unlikely to be regions of ordinary disk interstellar gas. Such distant material might represent outflows driven by the AGN, or recently acquired gas in tidal tails or streams. The morphology, kinematics, and abundances in the gas all bear on this question.

Tidal distortions are present in the stellar morphologies of 7/8 galaxies in this sample; in the remaining one, NGC 5252, we see dust out of the galaxy plane, and H i data show a rough polar ring of gas with no detected stellar counterpart. Where the large-scale kinematic pattern is dominated by ordered rotation and the abundances are sigificantly lower than at the nucleus, the case is strong that we are seeing photoionized tidal debris (although it is possible that outflows may be important in a small region around the AGN itself). This applies to UGC 7342, with some of the lowest abundance values we observe all across the ionized clouds; to SDSS 2201+11; and to Hanny's Voorwerp itself (in which H i data from Józsa et al. (2009) show the ionized gas to be part of a much larger, kinematically coherent structure). NGC 5252 falls in this category as well; the outer emission regions have significantly lower abundances than seen in the inner few arcseconds. Since we trace kinematics using emission lines, we do not always know whether the sense of rotational order matches that of the stars; in NGC 5972, we detect absorption lines well across enough of the galaxy to know that the senses do match.

The more complex kinematics of the Teacup system and UGC 11185, with greater velocity ranges and localized multiple components, make a case for outflows playing a more important role. In these systems, the extended-gas abundances are close to those in the respective nuclei, likewise consistent with an outflow. Kinematically, it is also possible that the second emission system seen in the eastern part of UGC 11185 is tidal debris projected against the host galaxy, but then one might expect lower abundances.

An additional signature of a small-scale outflow would be a velocity mismatch between the AGN emission and the interpolated central value of the velocity field (or slice along the slit location). We can examine this point for some of our sample. NGC 5252 shows a rotation of velocity contours by $\approx 120{}^\circ $ in the inner 10'', with a separate region in the inner 2'' at about 90° to the outer velocity field. Each of these might be as well described as a counterrotating region or a radial perturbation of order 70 km s−1 to a global pattern. In NGC 5972 and SDSS 2201, any specifically nuclear velocity offset must be at a level $\lt 30$ km s−1.

8. SUMMARY

We have used HST images and a range of supporting data to study the host galaxies of EELRs whose AGNs are insufficient to photoionize them to the observed levels, implying significant fading of the AGNs radiation on timescales $\lt {{10}^{5}}$ yr. The kinematics of the ionized gas are dominated by rotation. Outflows are not globally important in these objects, unlike the EELRs in some radio-loud AGN hosts. We see several clear cases of warped and precessing disks of gas and dust, and the gas is broadly associated with stellar tidal tails and irregular dust distributions in others. Each of the eight objects in our sample shows evidence of a merger or strong interaction; these EELRs are largely photoionized tidal debris. With only the single possible exception of Mkn 78, the sample of comparable nearby EELRs where the AGNs are energetically capable of lighting them up are likewise in interacting systems; this makes sense, given that the light-travel timescales we can study are much shorter than the dynamical times for evolution of tidal distortions (broadly, 105 versus 109 yr). Both the ionizing radiation from the AGNs and the large reservoirs of distant neutral gas are required to see EELRs such as these, giving a strong selection effect in favor of AGNs with tidal debris; episodic behavior of the AGNs similar to that we infer for these objects could be equally common in noninteracting galaxies but we would seldom be in a position to observe it. One avenue might be use of neighboring galaxies' IGM to measure or limit the AGN radiation reaching them; we have conducted a pilot survey for such cases (W. C. Keel et al. 2015, in preparation).

The kinematics in the extended clouds show some instances (SDSS 2201+11, UGC 7342, Hanny's Voorwerp itself) of material strongly dominated by ordered rotation. These generally show low (nitrogen/oxygen) abundances in the gas, consistent with an external tidal origin in less massive and lower-abundance systems. In some systems (the Teacup AGN, UGC 11185) , we see multiple superimposed velocity systems and localized high velocities (although lower than previously found in many radio-loud AGNs), sometimes accompanied by near-solar abundances across the clouds; these are candidates for outflows at least affecting the clouds.

In detail, the gas distributions show a variety of forms—giant loops, long filaments, and partial shells. The long filaments, and in some structures a breakdown of the expected relation between Hα surface brightness and ionization parameter in the gas, may indicate a role for magnetic fields, possibly pre-existing in the tidal tails rather than linked to the AGNs. Since detection these structures requires the angular resolution of HST (and in some cases lies beyond this limit), other tidal tails may be similarly structured, although we cannot yet observe at these levels with other techniques (such as H i mapping).

We find only marginal evidence for association between the locations of emerging UV flux (as traced by the ionized gas) and young star clusters; any role for triggering star formation by outflows in these cases must be very limited. None of the galaxies has extensive star formation; one is an S0, one an SB0, and the rest strongly bulge-dominated even when tidal tails are present.

This work was supported by NASA through STScI grants HST-GO-12525.01-A and -B. Some of the data presented in this paper were obtained from the Mikulski Archive for Space Telescopes (MAST). This research has made use of NASA's Astrophysics Data System and the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. We thank Linda Dressel for advice on setting up the observations, especially reducing intrusive reflections. Giles Novak provided useful pointers on magnetic fields in the ISM. Identification of this galaxy sample was possible through the efforts of nearly 200 Galaxy Zoo volunteers; we are grateful for their contributions, and thank once more the list of participants in Keel et al. (2012a). This work was partly supported by the Active Processes in Galactic and Extragalactic Objects basic research programme of the Department of Physical Sciences of the RAS OFN-17. The observations obtained with the 6 m telescope of the SAO of the RAS were carried out with the financial support of the Ministry of Education and Science of the Russian Federation (contract no. N14.619.21.0004, the project RFMEFI61914X0004). W.P. Maksym is grateful for support by the University of Alabama Research Stimulation Program. V.N. Bennert acknowledges assistance from a National Science Foundation (NSF) Research at Undergraduate Institutions (RUI) grant AST-1312296. Note that findings and conclusions do not necessarily represent views of the NSF. C. J. Lintott acknowledges funding from The Leverhulme Trust and the STFC Science in Society Program. A. Moiseev is also grateful for the financial support of the Dynasty Foundation and a grant from the President of the Russian Federation (MD3623.2015.2). K. Schawinski was supported by a NASA Einstein Fellowship at Yale, and gratefully acknowledges support from Swiss National Science Foundation Grant PP00P2 138979/1. Galaxy Zoo was made possible by funding from a Jim Gray Research Fund from Microsoft and The Leverhulme Trust. This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration.

Facilities: HST (ACS, WFC3, WFPC2), BTA, Shane (Kast Double spectrograph), SARA, WISE, KPNO:2.1 m

APPENDIX: OBJECT PROPERTIES

For more convenient reference, we collect here descriptions of our results organized by galaxy, rather than by technique as detailed earlier.

Mkn 1498 has the most prominent broad-line region seen in our sample, perhaps a clue to the orientation of its AGN structures. There are multiple partial circumnuclear rings of ionized gas at radii 0.5–1farcs 6 (0.6–1.8 kpc). While the starlight distribution of Mkn 1498 is nearly symmetric, it does include low-contrast fans or spirals, and there are chaotic dust lanes extending outward to 7.5 kpc. The fine-structure comparison of Hα surface brightness and [O iii] /Hα indicator of the ionization parameter shows mismatches in the extended northern and southern emission filaments. There are three candidate luminous star clusters, all located far from the emission-line filaments and dust lanes.

NGC 5252 has long been perhaps the best-known example of an AGN with ionization cones. This edge-on S0 system has a half-spiral of dust (possibly only the front half of a more complete inner disk), whose orientation is consistent with the kinematics of the ionized gas but not the stellar disk. As shown previously in both optical emission and H i data, the gas kinematics are strongly decoupled from the stars. The gas, presumably externally acquired, has internally well-ordered kinematics; in our Fabry–Perot maps, local regions of $\simeq 2$'' extent differ from a smooth velocity field by no more than ${{\sigma }_{v}}=27$ km s−1. We find no mismatch between Hα surface brightness and ionization parameters in the outer filaments, suggesting that the WFPC2 images resolved most of the density contrast here; the inner filaments are too small for this statistical test. The derived abundances show a substatial gradient from the near-nuclear features to the outer filaments, which could indicate that both a central outflow and less-processed tidal mateiral are present.

NGC 5972 combines low redshift with intrinsically large EELRs detected to at least 50 kpc radius, revealing very rich structure. The starlight shows tidal tails and shells, crossed by dust lanes which span a 66° range in orientation. These dust lanes can be fit schematically by a differentially precessing, warped disk with a structured radial distribution; such a distribution is seen in models with ages ∼1.5 Gyr. The gas velocity field includes asymmetric or multiple line profiles north and south of the nucleus, consistent with seeing into a patchy distribution of dust so the data sample different depths on adjacent lines of sight. The velocity field shows local departures ${{\sigma }_{v}}\lt 39$ km s−1 from a smoothed distribution; these departures are broadly correlated with dust lanes, indicating a role for line-of-sight weighting. There is only a partial match between the dust lanes and ionized gas in front of the central bulge; projection effects and illumination of the gas only within ionization cones seen nearly side-on could explain this. An [O iii] /Hα comparison shows evidence for unresolved density structure in the inner emission regions but not the outer areas. There is a modest population of candidate star clusters, not spatially associated with the dust lanes or ionized gas filaments. The outermost filaments of ionized gas align roughly with stellar tails, but their centerlines differ, suggesting either star/gas sorting as a result of different initial velocity distributions during the interaction or three-dimensional differences highlighted by the AGN's radiation pattern. The velocity field, both sampled with a long slit and via Fabry–Perot mapping, shows ripples where dust attentuation is important along the major axis, broadly consistent with the warped-disk geometry found from the dust lanes themselves.

SDSS 1510+07 is an SB0 galaxy with a faint tail connecting it to a dim companion at projected distance 45'' (40 kpc), and extending beyond that companion. Its EELRs consist of two large filamentary clouds roughly along the galaxy's minor axis, with some coherent filaments stretching as much as 7.5 kpc. These structures follow the expected relation between ionization level and Hα surface brightness, suggesting that the HST data resolve most of the density contrast in them.

SDSS 2201+11 has a continuum morphology consistent with a pure bulge plus absorption in a warped (differentially precessing) disk. The EELRs are two distinct clouds, largely disjointed from the emission around the AGN, about 20° out of the midplane of the dust and starlight, with a dust lane connecting to the northern cloud showing it to be on the near side (and receding, from the velocity data). The inner parts of the starlight distribution share the rotation sense of the gas. Deep images show outer tidal loops. The southern cloud may have two kinematically separate components. The dust lanes are well understood as parts of a precessing disk of material which would have been externally acquired ∼1.7 Gyr ago; if they are part of this disk, the EELR clouds may span a much wider cone angle than we observe. The gas velocity field has scatter among observed regions ${{\sigma }_{v}}\lt 59$ km s−1 about a smooth model. There may be a few luminous star clusters, but the candidates are too close to the dust lanes for accurate photometry. Both clouds show evidence of unresolved fine structure, through mismatches between the observed and predicted relations between Hα surface brightness and [O iii] /Hα. The abundances are substantially below solar, consistent with an origin of the cloud material in tidal distortion of a low-mass, low-abundance former companion.

The Teacup AGN has an EELR dominated by a 13 kpc filamentary loop to the northwest. The continuum shows an early-type morphology with a probable shell plus an oppositely directed tail, as seen in simulation of minor mergers involving a small, dynamically cold intruder. There are a few candidate luminous star clusters; a clump of such candidates half-behind the gas loop may be a background disk galaxy with redder core. Chaotic dust lanes appear up to 2'' (3 kpc) from the center, adding to the evidence for a past merger. Hα and [O iii] /Hα pixel values suggest unresolved structure at increasing levels for projected distances $r\gt 3\buildrel{\prime\prime}\over{.} 3$ (5 kpc) from the AGN, becoming more pronounced in the loop at $r\gt 5$''(8 kpc). Recent results by Harrison et al. (2015) show strongly blueshifted line emission just southeast of the AGN, and a radio-continuum counterpart to the large emission loop to thenortheast, suggesting outflow on these two scales. Near-solar abundances in the gas would be broadly consistent with an important role for AGN-driven outflow.

UGC 7342 shows a two-sided, roughly conical set of emission filaments extending at least 37 kpc from the AGN. The host galaxy shows complex dust lanes and tidal tails. The apparent companion galaxy SDSS J121820.18+291447.7 lies in the background at z = 0.0674, so the merger is advanced enough for the second galaxy to no longer be distinct. This system has a rich set of candidate luminous star clusters, with the brightest close to the center line of the brightest tidal tail. The gas velocity field shows ordered rotation, with significant ripples crossing the tidal tails. Some ionized gas filaments roughly match stellar tidal tails, but others have no continuum counterparts, so these distributions are largely uncoupled. The inner filaments in this system show a mismatch with predictions of the Hα–[O iii]/Hα relation, indicating substantial levels of unresolved density structure. We find low abundances throughout the EELR in this system, suggesting, with the morphological data, that it consists of tidal debris from a low-abundance disrupted system.

UGC 11185 is part of a strongly interacting galaxy pair, with a prominent tidal loop to the east and chaotic dust structure. Two of the four candidate star clusters are found together in this loop. The emission-line data suggest that the HST images resolve most of the density structure in the gas throughout the EELR, which plausibly traces a large set of ionization cones. The longslit data show multiple velocity compones in [O iii] to the east of the nucleus, with the strongest associated morphogicaly with the triangular or conical bright regions seen in the HST [O iii] image. Along with the relatively high abundances, this suggests that a nuclear outflow has been important in shaping the EELR of this system.

Footnotes

  • Based on observations with the NASA/ESA Hubble Space Telescope obtained at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract No. NAS5-26555.

Please wait… references are loading.
10.1088/0004-6256/149/5/155