# The Very Extended Ionized Nebula around the Quasar MR2251-178

## 1 Introduction

In the ongoing quest to better understand the luminous quasar population, one of the most fertile areas of study has been the investigation of their galaxy environments. The host galaxy not only must provide fueling for the central engine, but should also display the effects of the strong nuclear ionizing emission and the impact of dynamic activity, such as jets and winds.

Broadband observations of quasar environments suggest a surprising variety of host galaxy morphologies, including spirals and ellipticals, nearby companions, and tidal interactions. While early ground-based work indicated a large percentage of spirals amongst low-redshift quasar hosts ($`40`$%; Hutchings et al. 1984a ) and no concrete evidence for ellipticals, subsequent HST observations find an elliptical fraction of greater than half (Bahcall et al. (1997)). More recent HST work suggests in fact that nearly all radio-loud and radio-quiet quasars reside in massive ellipticals (McLure et al. (1999); see also Disney et al. (1995)). A substantial fraction of quasar host galaxies exhibit tidal tails and streamers indicative of interactions (e.g., Hutchings & Campbell (1983); Hutchings et al. 1984b ; Stockton & MacKenty (1987); Bahcall et al. (1997); Hutchings et al. (1999)). Evidence for extended gas from interactions has also been found at H I 21 cm (e.g., Lim & Ho (1999)). The presence of companion galaxies is common; the QSO/galaxy correlation function is evidently several times that for galaxies alone (Fisher et al. (1996)).

Given the strong line emission in the nuclei of active galaxies and quasars, as well as the prevalence of ionized gas in interacting systems, narrowband observations of quasar environments have proved interesting as well. The imaging survey of Stockton & MacKenty (1987), the largest to date, found highly-structured \[O III\] emission in a quarter of 47 luminous QSOs, with typical extents of a few tens of kpc. Similar spectroscopic observations have detected extended regions of line emission in half of the objects observed (Boroson et al. (1985)). The line ratios usually suggest photoionization by the nuclear power-law spectrum (e.g., Boroson et al. (1985); Boisson et al. (1994)), although stellar absorption lines from (presumably in-situ) stars are sometimes seen (e.g., Miller et al. (1996)). Kinematic studies of the ionized gas component are scarce, but the gas motions generally appear to be complex and chaotic (e.g., Durret et al. (1994)). The size and luminosity of the ionized component appear to be correlated with both the narrow-line nuclear luminosity and the radio power of the quasar (Boroson et al. (1985); Stockton & MacKenty (1987); Véron-Cetty & Woltjer (1990)).

The quasar MR 2251–178 is one of a few radio-quiet quasars which exhibit an extended gaseous envelope (e.g., Bergeron et al. (1983)). This quasar was first discovered as a strong and variable X-ray source by the Ariel V satellite (Cooke et al. (1978)). Subsequent observations identified the object as a quasar at a redshift of $`0.0638\pm 0.0015`$ (Ricker et al. (1978); Canizares et al. (1978)), residing in the outskirts of a small cluster (Phillips (1980)). In this Letter, we present deep H$`\alpha `$ observations of MR 2251–178 obtained with the TAURUS Tunable Filter (TTF), a new, etalon-based instrument which has been optimized for the detection of faint, extended emission-line gas. These new data allow us to better constrain the extent, velocity field, and origin of the ionized nebula around MR 2251–178.

## 2 Observations and Reductions

MR 2251–178 was observed on August 30, 1998 using the TTF at the f/8 Cassegrain focus of the 3.9-meter AAT. The TTF instrument consists of a pair of modified high-finesse ($`N40`$) Queensgate etalons (blue and red) which can be tuned to provide narrowband imaging anywhere within the wavelength range 400 to 960 nm, through an arbitrary bandpass, with resolving powers of 100 to 1000.<sup>1</sup><sup>1</sup>1Consult Bland-Hawthorn & Jones (1998) and the TTF web page (http://www.aao.gov.au/local/www/jbh/ttf/) for further details. The observations of MR 2251–178 were made with the red side of the TTF, using a mediumband ($`\mathrm{\Delta }\lambda =26.0`$ nm) blocking filter centered at 707 nm, tilted by 16°. Two 600-second exposures were obtained at each of two etalon spacings. The exposures were dithered amongst pointings on a 15″ grid. The average atmospheric seeing of 1.3″ was oversampled by the 0.37″ pixels. The night was photometric. MR 2251–178 was also observed on September 3, 1998 in direct imaging mode, using a standard I-band filter.

Fits to a number of emission lines from observations of a calibration lamp (CuAr) were used to determine the relationship between wavelength, spatial position, and etalon gap spacing. The free spectral range of the etalon, i.e., the distance between orders, was found to be 265.8 Å, well-matched to the bandwidth of the blocking filter. The system was used in the 26th order of interference, with a 9 micron etalon gap, at an effective finesse of 39.7. This translates to a spectral resolution of 6.8 Å, or $`R1040`$, with an effective bandpass of $`12`$ Å. The pair of etalon spacings produced imagery with central wavelengths on the optical axis of 6983 Å and 6986 Å, corresponding to redshifted velocities of 60 and 195 km s<sup>-1</sup> relative to the quasar. The field is essentially monochromatic: pixels $`1`$′ from the optical axis have a central wavelength only 1.3 Å ($`60`$ km s<sup>-1</sup>) bluer than pixels on-axis. The optical axis is located approximately 35″ northwest of the quasar.

The data frames were bias-subtracted and flatfielded in the typical manner. An azimuthally-symmetric sky frame was produced from a mean radial sky spectrum for each image and subsequently subtracted. Images at each of the two etalon spacings were then aligned and combined using 9 stellar objects in the field. An image mask was used to simultaneously remove three ghost reflections from each image. The data were flux calibrated using observations of the standard star EG 21 (Stone & Baldwin (1983)).

The I-band observations were reduced following standard CCD procedures. The summed continuum image was matched to the narrowband imagery using a number of stars and subsequently subtracted. A few remaining stellar residuals and cosmetic defects were identified and repaired by hand.

## 3 Discussion

Figure The Very Extended Ionized Nebula around the Quasar MR2251-178 presents our final H$`\alpha `$ images of MR 2251–178 at each of the two etalon spacings (panels $`a`$ and $`b`$), as well as the I-band continuum image (panel $`c`$) and a combined H$`\alpha `$ image (panel $`d`$). The 2-$`\sigma `$ detection limit of the H$`\alpha `$ images is $``$ 5 $`\times `$ 10<sup>-18</sup> erg s<sup>-1</sup> cm<sup>-2</sup> arcsec<sup>-2</sup>, more than an order of magnitude fainter than previously published data on this object. In the discussion that follows, we assume $`H_0=50`$ km s<sup>-1</sup> Mpc<sup>-1</sup> and a corresponding image scale of 1.9 kpc arcsec<sup>-1</sup>.

### 3.1 Morphology

Published optical imagery and spectroscopy have detected two ionized gas components around the quasar MR 2251–178: an elongated, highly-ionized circumnuclear component of diameter $`27`$ kpc, and an extended envelope of faint H$`\alpha `$\- and \[O III\]-emitting filaments out to a radius of $`110`$ kpc (Bergeron et al. (1983); Hansen et al. (1984); Alighieri et al. (1984); Macchetto et al. (1990)). Our observations confirm these findings and reveal a number of new features.

The circumnuclear component of ionized gas is evidenced by strong H$`\alpha `$ emission directly surrounding the quasar, extending out to a radius of $`20`$ kpc. This emission is slightly elongated in the east-west direction, i.e., along the axis of the quasar’s radio jet ($`PA102`$°; Macchetto et al. (1990)), and has been identified as an “extended emission-line region” (EELR) of the quasar host galaxy (Mulder & Valentijn (1992)). As others have reported on characteristics of the EELR in some detail, we will say nothing more of it herein.

The extended ionized gas component is observed as a much larger network of diffuse and filamentary emission exterior to the EELR. Above a flux level of $`1.8\times 10^{17}`$ erg s<sup>-1</sup> cm<sup>-2</sup> arcsec<sup>-2</sup>, diffuse H$`\alpha `$ emission is visible surrounding the galaxy to a radius of at least 50 kpc. Knots and filaments of ionized gas are visible to radii up to 120 kpc towards the northeast, southeast, and southwest. The filaments display a spiral morphology and exhibit a remarkable azimuthal symmetry that has not been previously seen. The total extent of this complex ($`200`$ kpc) makes it the largest known around a quasar. We derive a total flux from the extended ionized gas component of $`1.4\times 10^{13}`$ erg s<sup>-1</sup> cm<sup>-2</sup>, or a luminosity of $`2.5\times 10^{42}`$ erg s<sup>-1</sup> at the redshift distance of the quasar (384 Mpc). This value, as well as the corresponding upper limit on the gas mass of $`6\times 10^{10}`$ M, is comparable or slightly greater than previously published results (e.g., Nørgaard-Nielsen et al. (1986); Macchetto et al. (1990)).

A few of the individual knots and filaments have been previously noted; these are indicated in Figure The Very Extended Ionized Nebula around the Quasar MR2251-178$`d`$, using the denotations from Macchetto et al. (1990). Several additional knots also noted by these authors and others (e.g., Alighieri et al. (1984)) are conspicuous in their absence from our imagery. Most notable among these are the “E” knots interior to the eastern filament in Figure The Very Extended Ionized Nebula around the Quasar MR2251-178. As will be noted in the following section, this absence is probably due to the large range of velocity spanned by the filament emission, relative to our narrow bandpasses.

### 3.2 Kinematics

The spectroscopy which initially discovered the extended ionized component of MR 2251–178 also provided the first kinematical information on the nebula (Bergeron et al. (1983)). Those observations revealed a normal rotation curve along each of three position angles, with peak velocities of $`150`$ km s<sup>-1</sup>. Surprisingly, the extended nebula appears to be rotating in an opposite sense to the inner (EELR) regions of the galaxy (e.g., Nørgaard-Nielsen et al. (1986)). The only published two-dimensional kinematics of the system are narrow-field H$`\alpha `$ Fabry-Perot observations of the EELR (Mulder & Valentijn (1992)).

Comparison of our H$`\alpha `$ imagery at each of the two distinct etalon positions reveals a broad velocity gradient across the quasar nebula. The maximum gradient appears to be along a position angle of $`40`$°, with the eastern filaments generally less redshifted than the western filaments.

The eastern filaments are clearly detected in both the “blue” and “red” images. Since an upper limit of $`200`$ km s<sup>-1</sup> has been measured for their line widths (Bergeron et al. (1983)), the central velocity of these filaments must lie within the overlap region of our two velocity regimes. Accounting for its distance from the optical axis, the ionized gas east of the quasar must be near the systemic velocity of the quasar.

The western filaments are detected only in the “red” image and appear to be redshifted from the systemic velocity by a value of $`200300`$ km s<sup>-1</sup>. This correlates well with the imaging observations of Macchetto et al. (1990), whose broad ($`\mathrm{\Delta }\lambda =64`$ Å) H$`\alpha `$ filter barely detected the western filaments out to a velocity of $`200`$ km s<sup>-1</sup> relative to the quasar. Furthermore, these authors more readily detected the western filaments in their \[O III\] imagery, which encompassed a redder range of velocities.

We therefore confirm that the velocity structure of the extended ionized gas component does not appear to follow the rotation curve of the quasar (Nørgaard-Nielsen et al. (1986); Mulder & Valentijn (1992)). The broad velocity range of our images and the limitation of two distinct velocity samples restricts us from analyzing more precise velocity variations, e.g., amongst individual knots and filaments.

### 3.3 The Origin of the H II Envelope

Several scenarios have been put forth for the origin of the ionized gas envelope in MR 2251–178, including tidal debris from an interaction with galaxy G1 some $`2\times 10^8`$ years ago, similar debris from a merging event, clouds expelled from the quasar and/or host galaxy, a cooling flow, or the ionized portion of a large H I envelope that is gravitationally bound to the quasar (e.g., Nørgaard-Nielsen et al. (1986); Macchetto et al. (1990)). Our observations reveal a tremendous spatial extent for the emission in this system. We also observe a relatively high degree of symmetry in the envelope, including at least two “arms” of ionized gas and a pervasive diffuse component. This symmetry, together with the organized large-scale kinematics, casts doubt upon the interaction model for creation of the ionized gas. Furthermore, the galaxy G1, the usual culprit quoted as an interacting companion to MR 2251–178, has a velocity redshifted by 1246 km s<sup>-1</sup> relative to the quasar (Bergeron et al. (1983)), yet its spatial position corresponds to that of relatively blueshifted filaments. Published spectroscopy shows the emission lines from these filaments to be symmetric and significantly narrower than this (e.g., Bergeron et al. (1983)), although we cannot rule out the presence of a very faint, high-velocity streamer extending out toward G1.

Nevertheless, the coarse kinematic structure of the extended ionized gas indicates that it represents a distinct component from the inner EELR region, since the latter appears to rotate in an opposite sense. This implies that the envelope most probably did not originate within the host galaxy of the quasar. A retrograde merging event, in which a small, gas-rich galaxy has been subsumed by the quasar, remains a possibility, however again the relatively well-ordered large-scale kinematics and azimuthal symmetry of the ionized envelope are problematic.

The cooling flow hypothesis appears equally unlikely. As pointed out by Macchetto et al. (1990), the large size and ordered kinematics of the envelope, and the off-centered position and peculiar systemic velocity of the quasar with respect to the underlying cluster, all argue against this scenario. By elimination, our observations therefore favor a model in which the extended ionized envelope resides within a large complex of H I gas centered about the quasar. If this is correct, a deep search for 21 cm H I line emission should reveal the massive neutral envelope and constrain its origin (e.g., captured intergalactic H I clouds, remnant accreting gas from galaxy formation, etc.).

### 3.4 The Ionization of the H II Envelope

The high excitation level of the brighter knots around MR 2251–178 (e.g., Bergeron et al. (1983); Nørgaard-Nielsen et al. (1986); Macchetto et al. (1990)) requires an energetic source of ionization. In-situ ionization by a faint hot stellar component is unlikely to be dominant for a number of reasons. Deep broadband imagery of MR 2251–178 reveals a morphology very unlike that of our narrowband observations. The R-band image of Hutchings et al. (1999) exhibits a faint component extended along the N–S direction, markedly different from the symmetric spiral pattern that we observe in H$`\alpha `$. Moreover, Hutchings et al. (1999) do not detect an excess of continuum emission at the locations of any of the brighter H$`\alpha `$ knots in the envelope. Perhaps even more damaging to the stellar ionization scenario is our detection of a large amount of diffuse ionized gas outside of the brighter knots. We therefore favor an external source of ionization.

The presence of a significant amount of ionized gas at large angles relative to the jet axis (PA $`102`$°; Macchetto et al. (1990)) seems to rule out ionization mechanisms directly associated with the quasar jet (e.g., shocks). We are therefore left with the possibility that the ionization of the nebula is sustained by the quasar radiation field. As described in detail in previous studies (e.g., Bergeron et al. (1983); Nørgaard-Nielsen et al. (1986); Macchetto et al. (1990)), the power radiated by the quasar can easily account for the high ionization level of the nebula. However, the relative symmetry of the envelope requires that the ionizing radiation is escaping the quasar symmetrically with respect to our line of sight. The large-scale radiation field from MR 2251–178 therefore shows no sign of anisotropy or alignment with the radio axis, contrary to expectations from unified models of active galactic nuclei that rely on orientation effects (e.g., Barthel (1989); Antonucci (1993); Urry & Padovani (1995)). This symmetry of the radiation field, if typical of all quasars, may have important consequences on our understanding of the impact that quasars have on the intergalactic environment (e.g., proximity effect of Ly$`\alpha `$ clouds; Carswell et al. (1982); Murdoch et al. (1986)).

## 4 Summary and Future Directions

Our deep H$`\alpha `$ observations of MR 2251–178 reveal a spiral complex of ionized material that extends more or less symmetrically out to $``$ 120 kpc from the quasar. The coarse velocity field derived from our data shows a NE – SW velocity gradient that is opposite to that of the inner line-emitting region. The morphology and kinematics of the nebula cannot be easily explained by an interaction/merger event or a cooling flow. We favor the scenario in which the ionized material is part of a larger neutral envelope that is photoionized by the radiation field of the quasar. Deep H I 21-cm observations are needed to confirm this model.

We suspect that significant quantities of ionized gas may be present around luminous quasars, but have remained undetected with standard narrowband imaging techniques. The observations presented herein constitute only 40 minutes of integration, yet reach more than an order of magnitude deeper in flux than previous narrowband imagery. The TTF instrument allows the observer to tune the filter to a very precise bandpass in order to match the expected emission and avoid bright sky features. The TTF therefore promises to greatly impact a number of observational programs, such as those that aim at parameterizing ionized quasar envelopes and quantifying the impact of quasars on the intergalactic environment.

SV is grateful for support of this research by NSF/CAREER grant AST-9874973 and by a Cottrell Scholarship awarded by the Research Corporation.