A&A 452, 869-874 (2006)
DOI: 10.1051/0004-6361:20054578
L. Christensen1,2 - K. Jahnke1,3 - L. Wisotzki1 - S. F. Sánchez1,4 - K. Exter5 - M. M. Roth1
1 - Astrophysikalisches Institut Potsdam, An der Sternwarte 16,
14482 Potsdam, Germany
2 - European Southern Observatory, Casilla 19001, Santiago 19, Chile
3 - Max Planck Institut für Astronomie, Königstuhl 17, 69117
Heidelberg, Germany
4 - Centro Astronomico Hispano Aleman de Calar Alto, Spain
5 - Instituto de Astrofísica de Canarias, La Laguna, Tenerife, Spain
Received 23 November 2005 / Accepted 7 March 2006
Abstract
Powerful radio galaxies and radio-loud quasars at high redshifts are
frequently associated with extended emission-line regions
(EELRs). Here we investigate the [O II] EELR around the quasar
3C 196 at z=0.871 using integral field spectroscopy. We also detect
extended [Ne II] emission at a distance of about 30 kpc from
the core. The emission is aligned with the radio hot spots and shows a
redshifted and a blueshifted component with a velocity difference of
800 km s-1. The alignment effect and large velocities
support the hypothesis that the EELR is caused by a jet-cloud
interaction, which is furthermore indicated by the presence of a
pronounced bend in the radio emission at the location of the radio hot
spots. We also report observations of two other systems which do not
show as clear indications of interactions. We find a weaker alignment
of an [O II] EELR from the z=0.927 quasar 3C 336, while no
EELR is found around the core-dominated quasar OI 363 at z=0.63.
Key words: galaxies: kinematics and dynamics - quasars: emission lines - quasars: individual: 3C 196
Interactions of powerful radio emission with surrounding gas can explain the presence of extended optical emission line regions around high redshift, radio-loud AGN. However, most studies have focused on radio galaxies (e.g. Villar-Martin et al. 1997), rather than radio-loud quasars (RLQs), because there is no bright optical glare from the central AGN itself. In the case of broad absorption-line quasars, measurements of the blue shifted absorption lines have shown large outflow velocities (Turnshek 1984), but only a few observations have revealed a direct connection between the more extended gas and the central engine. Outflows driven by AGN activity are potentially an effective method for quenching star formation in the host galaxy (Di Matteo et al. 2005).
Luminous extended narrow emission-line regions (EELRs), reaching more than 100 kpc from the QSO nucleus, have been detected around quasars at redshifts between 1 and 4 (e.g. Bremer et al. 1992b,a; Heckman et al. 1991). At low redshift, detailed analyses of the host galaxies are enhanced by observations of the stellar continuum emission, and in several cases the line emission extends much farther than the continuum emission. Aside from interaction with the radio emission, possible explanations include photoionisation from the central source, or ionising radiation from massive stars in star-forming regions in the host galaxy. The interaction scenario is supported by observations of alignment of the radio jets with the optical continuum and extended emission line regions in a sample of radio galaxies (McCarthy et al. 1987). Such interactions result in shocks which compress and heat the gas, and subsequent recombination line emission will cool the material. As alternative explanations for EELRs, remnants from galaxy mergers (Stockton & MacKenty 1987), material falling into dark matter potentials (Haiman & Rees 2001), or cooling flows in massive galaxy clusters (e.g. Fabian & Crawford 1990) have been suggested.
Table 1: Log of the observations. The observations for 3C 196 were obtained in non-photometric conditions and all derived fluxes are relative.
These scenarios make different predictions for the surrounding material, and so the morphologies and kinematics of the EELRs can be used to discriminate between them. To investigate the structure and kinematics of these EELRs, narrow- and broad-band images, and slit spectroscopy have traditionally been used. Integral field spectroscopy (IFS) presents an alternative technique that allows imaging and spectroscopy simultaneously. IFS of six RLQs at redshift 0.26<z<0.60 showed that EELRs are common, although not always aligned with the radio axis (Crawford & Vanderriest 2000). An alignment between the radio axis and the extended emission around one lobe-dominated quasar was reported (Bremer 1997), and later IFS gave evidence for a jet-cloud interaction (Crawford & Vanderriest 1997). At z>2, detections of extended [O II] and [O III] emission lines have indicated no strong evolution with redshift, and a tendency for stronger line emission to be spatially coincident with the stronger radio emission (Wilman et al. 2000).
The AGN unification scheme states that radio galaxies and RLQs are the
same objects viewed at different angles relative to their jets
(Barthel 1989). IFS of radio galaxies have also indicated
jet-cloud interactions through the alignment effect and kinematics
(Sánchez et al. 2004; Márquez et al. 2000; Solórzano-Iñarrea & Tadhunter 2003). Like the EELRs from RLQs,
the alignment is roughly consistent with the radio morphology
(McCarthy et al. 1987). This alignment effect for RLQs and the EELRs is
also found at higher redshifts, but it is not as well-determined at
z>2 (Hutchings 1992; Heckman et al. 1991), possibly due to resonance
scattering of Ly
photons.
This paper presents observations of three RLQs at 0.6<z<0.9 and the first analysis of the alignment effect of the EELRs in this redshift range using IFS.
The data were obtained in connection with a project aimed at detecting
emission from galaxies associated with intervening damped
Lyman-
absorption lines (Christensen et al. 2005,2004).
From the original sample of seven QSOs at z<2, only three RLQs have
redshifted strong optical emission lines within the wavelength range
of the observations.
The observations were carried out with two integral field spectrographs; INTEGRAL (Arribas et al. 1998) mounted on the 4 m William Herschel Telescope, La Palma, and the Potsdam Multi Aperture Spectrophotometer (PMAS) mounted on the 3.5 m telescope at Calar Alto (Roth et al. 2000,2005). Table 1 presents a log of the observations.
The observations were obtained with the SB2 fibre bundle, which
consists of 189 object fibres plus 30 sky fibres arranged in a
ring with a diameter of 90
.
Each fibre has a diameter of 0
9 on the sky, giving a main field of view of
with a non-contiguous sampling and a
filling factor of about 67%. We used a 600 lines mm-1 grating
with a dispersion of 3 Å pixel-1 and a spectral resolution of
6 Å measured from the width of sky emission lines.
Data reduction was performed using IRAF tasks modified specifically
for the reduction of INTEGRAL data (see Garcia-Lorenzo et al. 2005). Bias
frames were subtracted and the frames were cleaned for cosmic ray hits
using the algorithm described in Pych (2004). All spectra were
extracted using the trace of the 219 spectra on the CCD found from an
exposure of a continuum lamp obtained at the beginning of the night.
Wavelength calibration was done using spectra of emission line lamps
also extracted for each fibre. The rms deviation of the wavelength
for sky emission lines for each spectrum was determined to be <0.2 Å. Differences between the individual fibre transmissions as a
function of wavelength were corrected for by modeling the transmission
of sky flat frames obtained at twilight. Two regions on the CCD were
affected by scattered light but there only in the blue end of the
spectral range and only affecting some of the spectra. Otherwise,
scattered light was found to be negligible compared to the overall
read noise of the CCD. Some sky fibres are contaminated by the QSO
flux because they are located next to the QSO spectra on the CCD, and
therefore affected by cross talk. The sky fibres were examined, and
those uncontaminated by the QSO flux were averaged and subtracted from
each object spectrum. We do not expect that any extended line
emission region is present at a distance of 45
.
To facilitate
inspection and visualisation, the data were interpolated onto a cube
of square spaxels (spatial elements) with sizes of
.
This interpolation and further
visualization of the data cubes was done with the Euro3D visualization
tool (Sánchez 2004). When one-dimensional spectra are extracted,
we retain the original spaxels sizes, instead using the interpolated
data cube.
The method for reducing PMAS data was essentially the same as for INTEGRAL data and was done with our own IDL-based software package P3d (Becker 2002). The main difference was the sky subtraction because PMAS does not have allocated sky fibres. Instead, an average sky spectrum was created from spectra at the edge of the field of view uncontaminated by the QSO or the nebulae. Because the spatial position of the nebulae was previously unknown, several different selections of sky spectra were examined before selecting an appropriate sky background spectrum. For the final sky background spectrum we selected spaxels where no emission lines were detected visually around the wavelength range of interest, but we did not detect line emission for the fibres at the edge of the field. Hence, an over-subtraction of the extended emission should be a small effect.
Both data sets were flux calibrated by comparison with observations of spectrophotometric standard stars observed at the same nights and with the same setup as used for the objects.
![]() |
Figure 1: One-dimensional spectra of the quasars extracted from the data cubes, by co-adding spectra from 20-30 spaxels. Emission lines have been indicated. |
Open with DEXTER |
To detect the EELRs it is first necessary to subtract the nuclear QSO emission because this is by far the dominant contribution. To get a clean spectrum of the EELR emission we use the approach described in Sánchez et al. (2004) and Sánchez et al. (2006). First, a two-dimensional model of the QSO PSF is made for each slice in the data cube with a one-pixel width using the program GALFIT (Peng et al. 2002). We allow both a point source representing the continuum (Gaussian component), and an extended host galaxy component to be present, but find that the host component is not detected in any of the objects. Assuming that the spatial location of the QSO emission varies smoothly with wavelength we make a model PSF in the form of a data cube. This model data cube is subtracted from the original date to create a residual data cube.
The INTEGRAL data cube of this object clearly shows extended
[O II] emission, even in the frame where the QSO emission is
not subtracted. After subtracting the QSO emission, residuals from
the QSO are still present due to uncertainties in the wavelength
dependent PSF determination, but this effect is strong for only the
central region as shown in Fig. 2. The two radio
hot-spots found at 5 GHz with a separation of 5
8 at
PA
120
(Pooley & Henbest 1974) are indicated by the
"+'' signs. The line emission is predominantly located in two regions
as noted before, and the morphology is similar to the [O II]
narrow band image presented in Ridgway & Stockton (1997). The north-western
part of the nebula appears to have a continuum counterpart which
extends to about 2
from the QSO (Boisse & Boulade 1990; Ridgway & Stockton 1997),
but the IFS data is not sensitive to the continuum. North of the
quasar we find that the line emission extends out to 3
(22 kpc), while a brighter region extends 5
(36 kpc) to the
south. Here and throughout we have assumed a flat cosmology with
H0=70 km s-1 Mpc-1,
and
= 0.7.
![]() |
Figure 2:
A narrow-band image of the [O II] emission at
z=0.871 (6954-6990 Å) from 3C 196 interpolated to a grid
scale of 0
![]() |
Open with DEXTER |
Co-adding all the spectra associated with the extended emission after
the QSO emission is subtracted results in an emission line with a
total flux of
erg cm-2 s-1 with a FWHM of
km s-1. The width [O II] does not vary
significantly over the face of the nebula within the uncertainties of
the emission line FWHM. However, it is expected that a
jet-cloud interaction would result in an increased velocity dispersion
at the location of the hot spots. With the current data we can not
detect this. Higher spatial resolution data are needed to determine
if this is the case.
Besides [O II] emission, we also detect extended
[Ne III] 3869 emission associated with the brightest
[O II] region to the south of the QSO, as shown in the
extracted one-dimensional spectrum in Fig. 3.
This is similar to the long-slit spectrum presented in
Steidel et al. (1997). The extension of the [Ne III] nebula appears
smaller than for the [O II] region, but the signal-to-noise
level does not allow for a detailed analysis. The emission line ratio
log([Ne III]/[O II])
is larger
than observed for H II regions in the Large Magellanic Cloud
(Oey et al. 2000), where values between -1.6 and -0.7 are found,
indicating that the ionising flux is harder around the RLQ. Compared
to a few other radio-loud objects (Stockton et al. 2002; Solórzano-Iñarrea & Tadhunter 2003), it
shows a similar line ratio. We note that the two emission line
FWHM are very different;
Å for [O II],
while [Ne III] is barely resolved (
Å when
corrected for the instrument resolution). This could signify that the
emission originates in different volumes, in which case the line ratio
has no physical meaning. Alternatively, a jet-cloud interaction could
significantly increase the FWHM of the [O II] line at
the location of the hot spots, while the [Ne III] only arises
from the brighter part of the nebula.
To check for the presence of fainter extended emission farther away
than 5
from the quasar centre, 30 spectra are co-added at
several spatial locations. From the non-detections of emission lines,
we derive an upper limit for emission in the surrounding field of
erg cm-2 s-1 arcsec-2 in the
INTEGRAL data. Not surprisingly, we neither find continuum emission
farther than 5
from the QSO.
![]() |
Figure 3: One-dimensional spectrum of the brightest region to the south of 3C 196 where both [O II] and [Ne III] emission appear extended. These lines are indicated by the vertical lines. The lower panel shows a one-dimensional spectrum of the northern region where no [Ne III] emission is detected. |
Open with DEXTER |
The spectra from various fibres show that the extended emission
around the QSO is progressively shifted in wavelength.
Fitting Gaussian profiles to the individual spectra gives
evidence for systematic velocity structure over the nebula, as shown in
the upper panel of Fig. 4. Only spectra with
[O II] detections larger than 2
above the background
are shown in colour. Compared to the systemic redshift of the QSO,
the bright southern region has a relative blueshift of -360 km s-1, whereas the fainter north western region is redshifted by
up to +500 km s-1. Another representation in the lower panel in
Fig. 4 shows a grey scale map of the intensity in
the [O II] line emission region derived by fitting emission
lines in individual spectra, and the velocity structures of the
emission are indicated by contours.
![]() |
Figure 4: Upper panel: velocity structures of the extended [O II] emission around 3C 196 determined after the QSO emission was subtracted. The zero-point of the velocity is the QSO redshift. Outlined circles indicate the spectra where no [O II] emission lines could be fit due to their low signal-to-noise levels. One spaxel close to the QSO centre has not been fitted, because the emission line is affected by QSO subtraction residuals. (See the online edition for a colour version of this figure.) Lower panel: image of the [O II] line emission intensity around 3C 196 with velocity contours overlaid. Contour levels are separated by 100 km s-1 and negative velocities are shown by the dashed lines. The "+'' signs indicate the positions of the two radio hot-spots which are roughly spatially coincident with the brightest emission line regions. |
Open with DEXTER |
Because our observing campaign was targeted towards finding emission
from a z=0.4367 galaxy responsible for the damped Ly
(DLA), the
wavelength coverage is not optimal compared to similar studies of
other QSOs, which use the [O II]/[O III] ratio to
derive an internal gas pressure (e.g. Crawford & Vanderriest 2000). Therefore
the present data cube does not allow modeling of the ionising
conditions in the nebulae.
The position of the EELR is not exactly aligned with the hot spots to the one arcsec accuracy. However, when we consider that the spectra do not sample the emission contiguously, both the position of the northern and the southern radio hot-spots are spatially coincident with the brightest emission line regions. The southern lobe is associated with the blueshifted emission, and the northern one with the redshifted part. This presents strong indications for an interaction with the radio jets. However, without polarization measurements from each hot spot we do not have the information on which is the nearest, and hence not completely rule out in-flow of the surrounding gas.
A galaxy responsible for the DLA line at z=0.4367 has been detected
1
5 the the south-east of the QSO
(Boisse & Boulade 1990; Chen et al. 2005; Le Brun et al. 1997). At this redshift the H
emission line falls close to the [O II] line at the QSO
redshift, and so potentially the south-western emission could
originate from this DLA galaxy. However, it is unlikely that the
brightest region of emission to the south is contaminated; the most
likely source for the emission is the QSO environment.
The [O II] emission line is within the wavelength range of the
IFS observations, but a subtraction of the QSO emission from the IFS
data cubes reveals no extended emission in either data sets. A
narrow-band image created from the INTEGRAL data cube at the
wavelength region around [O II] at z=0.63 with a width of
20 Å is shown in Fig. 5. No bright regions
with associated emission lines are found. Other narrow-band images
offset by up to 1000 km s-1 relative to the quasar redshift show
no line emission regions either. We also note that the nuclear
[O II] emission line shown in Fig. 1
appears weak. We estimate an upper detection limit for extended
[O II] emission of about
erg cm-2 s-1 arcsec-1in the INTEGRAL data set and
erg cm-2 s-1 arcsec-1 in
the PMAS data set, based on experiments with artificial emission
lines.
Hutchings (1992) obtained broad-band and narrow-band images
centered on the [O II] emission line at the quasar redshift and
found evidence for a detached region of [O II] emission
approximately 2
in size and offset by
3
to the
north-east. By creating a narrow-band image from the data cube
with the same wavelength range as in Hutchings (1992) we search for
this detached object, but it is not detected in our data.
![]() |
Figure 5:
Narrow-band image of OI 363 created from the INTEGRAL data at
![]() ![]() |
Open with DEXTER |
Figure 6 shows a narrow-band image created from the PMAS data cube. The morphology of the emission region is similar to the 28 Å wide [O II] narrow-band image presented in Ridgway & Stockton (1997), but like these authors we cannot confirm the extended emission reaching towards the north. The extended emission region appears to the south and north-west of the QSO, with the brightest region towards the south. The latter is spatially coincident with the radio jet, whereas the north-west emission line region is not aligned.
The PMAS data only allows detection of the emission to 2
from the QSO, and a detailed analysis of the kinematics of the
emission line region is not possible with the signal-to-noise level in
each individual spectrum from the data cube. In the co-added
one-dimensional spectrum for the whole emission line region the total
line flux measured is
erg cm-2 s-1 after a
correction for Galactic extinction is applied
(Schlegel et al. 1998). After correcting for the instrumental resolution,
the line width is
km s-1.
A narrow-band image from the data cube, offset to slightly longer
wavelengths (7190-7206 Å) shows that the emission feature extends
further to the north-west and is spatially coincident with the
emission feature denoted "2'' in a continuum image of
Steidel et al. (1997), which is marked by "+'' in
Fig. 6. A one-dimensional spectrum of that region in
the data cube indicates
,
in agreement with
z=0.931 measured by Steidel et al. (1997). Furthermore, the data cube
allows a confirmation of their object "3'' at z=0.892, 3
to
the east of the QSO. Possibly, the emission to the north-west is
unrelated to the QSO and only belongs to the galaxy "2''; thus we
determine the properties of the southern region only. Here the flux is
erg cm-2 s-1 and the FWHM is
km s-1. The redshift measured is
,
corresponding to a velocity difference of
km s-1relative to the QSO redshift.
![]() |
Figure 6:
Narrow-band image of the [O II] line emission region
from 3C 336 in the region 7167-7193 Å. The image is
interpolated to a pixel scale of 0
![]() ![]() ![]() |
Open with DEXTER |
Compared to the other RLQs where an alignment has been detected
(e.g. Crawford & Vanderriest 2000), the EELR around 3C 196 suggests the presence
of a jet-cloud interaction with strong velocity differences. The
emission line nebula has both a red- and a blueshifted component
aligned with the radio hot spots, and the velocity is larger than
expected for a gravitational origin. Rotation velocities for the most
rapidly rotating massive spiral disks are around 300 km s-1 at a
distance of about 10 kpc. If the EELR is a rotating disk, its mass
would have to be at least twice that of such a massive spiral. The
bend of the radio emission at the location of the hot spots noted in
Brown et al. (1986) supports the interpretation of an interaction.
The [O II] emission from the galaxy is spatially offset from
the intervening galaxy found by Boisse & Boulade (1990), which at a redshift
of 0.87 has H
emission lines very close to the [O II] at
the QSO redshift. With emission line fitting we are able to separate
the H
emission line from the [O II] emission, because
there is a 10 Å difference between the two lines, and the
positions are not exactly spatially coincident.
Extended [Ne III] emission is detected for 3C 196 at a distance of 30 kpc from the quasar nucleus, but only in the southern region where the [O II] emission is strongest. Because the [Ne III]/[O II] line flux ratio is larger than observed in H II regions, harder ionising radiation must be present if the emission originates in the same volume. For a detailed modeling more emission lines need to be observed. Ionisation conditions can in principle be analysed in the EELRs provided that the [O III] emission line is detected too. Combining the observed emission line ratios with knowledge of the quasar ionising flux allows determination of the nebular pressure through photoionisation modeling. However, this approach is only justified when emission lines arise in the same medium. This is sometimes neglected, as pointed out by Stockton et al. (2002), who find that the [O II] emission arises in a much denser medium than [O III] in one RLQ. Measurements of other line ratios could help to determine the nature of the nebulae and the presence of interaction signatures (e.g. Garcia-Lorenzo et al. 2005).
The alignment of the nebular emission line regions with the radio jets indicates that some interaction could be present. To determine whether the line emission is created by the ionising quasar flux which is collimated along the radio jet, it is necessary to determine the underlying continuum flux. If a significant stellar emission component is present then one can rule out this scenario. No extended continuum emission is detected in either object analysed here, but if it were present, the bright nuclear flux complicates the detection. Deeper (imaging) data could resolve this issue. Furthermore, if the extended emission is caused by an interaction with the radio jet, larger field of view observations of 3C 336 are necessary to analyse the effects at larger distances corresponding to the widely separated radio-lobes.
In the unification scheme, the core-dominated RLQs have their radio jet orientation closely aligned with the sight line. If radio jets cause interaction with the surrounding material and give rise to line emission, a core-dominated RLQ would have less extended emission than lobe-dominated RLQs. No EELR is detected for the one core-dominated RLQ analysed here.
An alternative scenario where the extended emission is dominated by the QSO ionising radiation directed in a cone (Haiman & Rees 2001; Weidinger et al. 2005) can not be rejected completely on the basis of this investigation. However, the velocity structure of the nebulae around 3C 196 supports the hypothesis that it is the interaction of the material with the radio jet which is the dominant effect and causes an outflow of the material.
Acknowledgements
L. Christensen acknowledges support by the German Verbundforschung associated with the ULTROS project, grant No. 05AE2BAA/4. S. F. Sánchez and K. Exter acknowledge the support from the Euro3D Research Training Network, grant No. HPRN-CT2002-00305. K. Jahnke acknowledges support from DLR. project No. 50 OR 0404.