A&A 371, 97-106 (2001)
DOI: 10.1051/0004-6361:20010359
I. Márquez1 - P. Petitjean2,3 - B. Théodore4,5 - M. Bremer6 - G. Monnet7 - J.-L. Beuzit8
1 - Instituto de Astrofísica de Andalucía (CSIC),
Apdo. 3004, 18080 Granada, Spain
2 -
Institut d'Astrophysique de Paris, CNRS, 98bis Bd Arago,
75014 Paris, France
3 -
UA CNRS 173 - DAEC, Observatoire de Paris-Meudon, 92195 Meudon
Cedex, France
4 -
Service d'Aéronomie du CNRS, BP 3, 91371 Verrière le Buisson,
France
5 -
ACRI, 260 route du Pin Montard, BP 234, 06904 Sophia-Antipolis,
France
6 -
Bristol Univ. (Dept. of Physics)
H H Wills Physics Laboratory, Tyndall Av, Bristol BS8 1TL, UK
7 -
European Southern Observatory, Karl Schwarschild Strae 2,
85748 Garching-bei-München, Germany
8 -
Canada-France-Hawaii Telescope Corporation, 65-1238 Mamaloha Highway, Kamuela, HI 96743, USA
Received 26 January 2001 / Accepted 8 March 2001
Abstract
We present the results of adaptive-optics imaging in the H and K bands
of 12 low and intermediate redshift (z < 0.6) quasars using the
PUEO system mounted on the Canada-France-Hawaii telescope. Five
quasars are radio-quiet and seven are radio-loud. The images,
obtained under poor seeing conditions, and with the QSOs (mV > 15.0) themselves as reference for the correction, have typical
spatial resolution of FWHM 0.3 arcsec before deconvolution. The
deconvolved H-band image of PG 1700+514 has a spatial resolution
of 0.16 arcsec and reveals a wealth of details on the companion and
the host-galaxy.
Four out of the twelve quasars have close companions and obvious signs
of interactions. The two-dimensional images of three of the host-galaxies unambiguously reveal bars and spiral arms. The morphology of
the other objects are difficult to determine from one dimensional surface brightness profile and deeper images are needed.
Analysis of mock data shows that elliptical galaxies are always
recognized as such, whereas disk hosts can be missed for small disk scale
lengths and large QSO contributions.
Key words: galaxies: active - galaxies: quasars - galaxies: fundamental parameters - galaxies: photometry - infrared: galaxies
Evidence that nearby bright galaxies contain massive dark objects in their center has become increasingly compelling over the last few years and early suggestions that a tight correlation exists between the mass of the dark object and the mass of the bulge (Kormendy & Richstone 1995) have been convincingly corroborated (Magorrian et al. 1998; Ferrarese & Merrit 2000). It is thus possible that AGN activity is a usual episode of the history of most, if not all, present-day bright galaxies. One way to investigate this is to determine the luminosity and morphology of galaxies hosting quasars. In addition, this gives clues on the range of conditions needed for strong nuclear activity to occur.
Recent observations of host-galaxies with HST have questioned the previous belief that radio-loud and radio-quiet quasars are found preferentially in the center of, respectively, elliptical and spiral galaxies (Bahcall et al. 1997; McLure et al. 1999). This supports previous conclusions from IR ground-based observations (Dunlop et al. 1993; McLeod & Rieke 1995; Taylor et al. 1996) that the distinction in host-galaxy characteristics between the two classes of QSOs is subtle (see also Hooper et al. 1997, for larger redshift). More surprising is the finding by McLure et al. (1999) that host-galaxies of luminous AGNs (both radio-quiet and radio-loud) are all massive ellipticals whereas Bahcall et al. (1997) found that, in the same redshift range, host-galaxies are of various morphologies but that all radio-loud quasars have bright elliptical hosts or occur in interacting systems (Kirkhakos et al. 1999).
The detection and analysis of host-galaxies is difficult even from space (Bahcall et al. 1994; Bahcall et al. 1995; McLeod & Rieke 1995; McLeod & McLeod 2001). Indeed, the determination of the PSF and the subtraction of the point source image are crucial in this work. Differentiation between the two classical profiles, either an exponential disk or a de Vaucouleurs power-law, is effective only in the regions close to the center, or in the far-wings of the PSF (see Fig. 1 of McLure et al. 2000). This is the main limiting factor in determining host-galaxy morphologies from the ground. With the advent of adaptive optics, it will be possible to alleviate the seeing limitation (Stockton et al. 1998; Aretxaga et al. 1998; Hutchings et al. 1999). Observing in the infrared minimizes the difference in luminosity between the host and nucleus again improving our ability to determine the host morphology.
PSF determination is still a major problem but the difficulties are balanced by the prospect of using very soon 10 m class telescopes which will provide higher sensitivity and better spatial resolution.
In this paper we present the result of a pilot programme aiming at testing the capabilities of adaptive optics in this field. We present the data in Sect. 2, discuss each object in Sect. 3, analyze the results in Sect. 4 and conclude in Sect. 5.
Name | Other names | mv | z | run1 | Exposure | FWHM2 |
(minutes) | (arcsec) | |||||
Q 0955+326 | 3C 232 | 15.8 | 0.530 | 1 | 20 | 0.36 |
PG 1001+291 | 16.0 | 0.329 | 1 | 60 | 0.24 | |
PG 1012+008 | 15.6 | 0.185 | 1 | 20 | 0.26 | |
PKS 1302-102 | 15.2 | 0.286 | 1, 2 | 20, 48 | 0.26, 0.42 | |
PG 1402+261 | 15.5 | 0.164 | 1, 2 | 40, 64 | 0.26, 0.34 | |
B2 1425+26 | 15.7 | 0.366 | 1, 2 | 35, 64 | 0.24, 0.34 | |
Q 1618+177 | 3C 334 | 16.0 | 0.555 | 1 | 30 | 0.36 |
PG 1700+518 | 15.1 | 0.290 | 1, 2 | 40, 40 | 0.26, 0.30 | |
Q 1704+608 | 3C 351 | 15.3 | 0.371 | 1, 2 | 40, 64 | 0.26, 0.43 |
B2 1721+34 | 15.4 | 0.206 | 1 | 75 | 0.24 | |
PG 2112+059 | 15.6 | 0.466 | 1, 2 | 30, 32 | 0.24, 0.34 | |
PKS 2128-12 | 16.2 | 0.501 | 2 | 24a | 0.61, 0.48 |
Name | z | Companions | Extension | ( ) | ( ) | mr1/4H( mr1/4K) | Host type | ||
5 arcsec | 10 arcsec | arcsec | kpca | ||||||
Q 0955+326 | 0.530 | 2.2 | 7.9 | 14.18 | 14.25 | 12.56 | E | ||
PG 1001+291 | 0.329 | 2.9 | 8.1 | 13.87 | 14.10 | 14.84 | SBa | ||
PG 1012+008 | 0.185 | 1 | 1 | 3.0 | 6.3 | 16.71 | 16.97 | 17.19 | E-Sa |
PKS 1302-102 | 0.286 | 2 | 2.6 | 7.1 | 13.45 (12.53) | 13.84 (12.89) | 13.84 (13.08) | E | |
PG 1402+261 | 0.164 | 3.0 | 5.3 | 13.23 (11.91) | 14.07 (12.86) | 14.07 (11.67) | SBa | ||
B2 1425+26 | 0.366 | 2 | 2.3 | 7.5 | 14.26 (13.30) | 14.69 (14.01) | 15.01 (14.23) | E-Sa | |
Q 1618+177 | 0.555 | 1 | 2 | 1.5 | 5.4 | 14.68 | 15.18 | 15.48 | |
PG 1700+518 | 0.290 | 1 | 2.5 | 6.7 | 12.90 (11.75) | 13.89 (12.50) | 13.89 (12.15) | Sa? | |
Q 1704+608 | 0.371 | 1 | 2.3 | 6.9 | 13.30 (12.21) | 14.37 (12.89) | 14.30 (12.84) | E | |
B2 1721+34 | 0.206 | 2 | 1.5 | 3.2 | 13.95 | 14.13 | 15.07 | ||
PG 2112+059 | 0.466 | 1 | 1 | 2.0 | 6.6 | 13.64 (12.64) | 14.10 (12.96) | 14.71 (13.51) | E? |
PKS 2128-12 | 0.501 | 1 | 1.7 | 6.1 | 13.42 (12.77) | 14.39 (13.33) | 14.50 (13.21) | ||
a H0 = 100 km s-1 Mpc-1; = 0.5. |
In order to use adaptive optics correction quasars were selected such that the nuclei were bright enough to be used as the wavefront reference point source. The sample of radio-quiet quasars were all PG quasars with and with redshift less than 0.6. The radio-loud objects were selected from 3C, 4C, B2 and PKS catalogues with the same magnitude and z criteria. The final objects observed (see Table 1) were selected based upon the suitability for the observing conditions on the observing runs.
We used the CFHT adaptive optics bonnette (PUEO) and the IR camera KIR on May 1998 (run 1) and May 1999 (run 2). The weather conditions were poor during both runs and the FWHM of the seeing PSF was never better than 0.8 arcsec. The adaptative-optics correction was performed on the QSOs themselves. The quasar was centered successively in the center of the four quadrants of the detector. The exposure time for individual images was two minutes. The background was determined by median-averaging the frames and the flat-field was taken to be the normalized dark-substracted background. The images were then aligned and added. The final images have a typical resolution of FWHM 0.3 arcsec. After each science observation an image of a star with similar magnitude as the QSO was taken in order to determine the PSF and use it to deconvolve the images. Due to rapid variations in the wheather conditions however, it was not always possible to follow this predefined procedure.
A synthetic PSF function, derived from the stellar images was used to deconvolve each of the images. As it was not always possible to apply a standard procedure due to fluctuating seeing conditions, a careful, although somewhat arbitrary choice of the PSF had to be done. In Fig. 2 we show the images of PKS 1700+514 obtained using, for the deconvolution, three different PSFs from stars observed during the same night. These have respectively, FWHM = 0.30, 0.42 and 0.48 arcsec. The initial image of the object has a resolution of FWHM = 0.26 arcsec and the star observed just after the science exposure has FWHM = 0.48 arcsec. It is apparent that the best result is obtained using the star with the FWHM closest to that of the science exposure. Here, we were guided in the exercice by the existence of the HST image by Hines et al. (1999). In general, this illustrates the crucial role played by a careful PSF determination in AO observations.
Results are summarized in Table 2. Columns and give, respectively, the number of objects (probably companions) found within 5 and 10 arcsec from the quasar down to mH = 20.5; columns and give the maximum radial distance (in arcsec and kilo-parsec) to which the host is detected at a significance level of 3 above the background; column gives the total magnitude of the object in the H-band and columns and those of the host-galaxy as derived from the PSF subtraction and profile fitting respectively (see Sect. 4); the assigned morphology, which comes from the 2D brightness distribution and the comparison of the two profile fittings, is given in column .
Figure 1: Images of PKS 1302-102 in the H-band. The top-left panel and bottom panel correspond to the image, respectively, before and after deconvolution. The spatial resolutions are respectively FWHM 0.32 and 0.24 arcsec. The top right panel shows the two companions after substraction of the PSF and a model for the host-galaxy. The inset in the bottom panel corresponds to a higher contrast version of the inner 4 arcsec 4 arcsec | |
Open with DEXTER |
Figure 2: Images of PG 1700+514 after deconvolution using the PSF given by three stars observed during the same night. The best deconvolution is obtained using the star with the FWHM closest to that of the quasar. The resulting image (bottom) has a final resolution of FWHM = 0.16 arcsec | |
Open with DEXTER |
The image obtained at CFHT is shown in Fig. 2. We confirm the findings by Stockton et al. (1998) that the companion has the appearance of an arc with several condensations. We used different PSF to deconvolve the image. The best deconvolution is obtained using the star with the FWHM closest to that of the AGN (0.30 arcsec, see Sect. 2). The image has a final resolution of 0.16 arcsec and is probably the best image obtained yet on this object. The companion is seen as a highly disturbed system with a bright nucleus and a ring-like structure; the nucleus being decentered with respect to the ring. The host-galaxy is clearly seen around the quasar with a bright extension to the south-west, first noted by Stickel et al. (1995) and clearly visible in the optical images by Stockton et al. (1998). In addition, we detect a bright knot to the south-east which is not seen in the NICMOS data probably because of the presence of residuals in the PSF subtraction. The comparison between the HST and CFHT images of PG 1700+514 shows how powerful AO can be, and bodes well for the use of the technique on 10 m-class telescopes. No obvious relation is found between the near-IR image and the radio map (Hutchings et al. 1992).
3C 232 - TON 0469 This object has attracted much interest because it is located 2 arcmin (8h-1 kpc at z = 0.0049) north of the nearby (z = 0.0049) galaxy NGC 3067 and has been considered as the prototype of the galaxy-quasar physical association as a probe for anomalous redshifts (Hoyle & Burbidge 1996). The HST spectrum shows a strong Mg II system at = 0.0049 (Tumlinson et al. 1999) which has been shown to be due to an H I tidal tail originating in the disk of the foreground galaxy (Carilli et al. 1989). The MERLIN map at 1.6 GHz shows a slight extension to the east (Akujor et al. 1994). The host-galaxy is detected in the CFHT image with an extension to the south-west up to 5 arcsec from the quasar). From both the image (Fig. 3) and the profile fitting (Fig. 4) we derive that the host is better described as an elliptical galaxy.
PG 1001+291 - Ton 0028 Only one image from the HST archive has
been published yet (Boyce et al. 1999). These authors favor the
interpretation that the object is an interacting object with two
nuclei, one 1.9 arcsec to the south-west and the other 2.3 arcsec to
the north-east. We find no evidence for the presence of two nuclei. In
particular there is no residual after subtraction of a
host-galaxy model. The profile is not well fitted by a single disk
component, due to the presence of a bump which is characteristic of
the bar component in the profiles of barred spirals. We therefore
conclude that the elongated NE-SW feature seen in our image probably
reveals the presence of a strong bar in an otherwise spiral galaxy.
PG 1012+008 This radio-quiet quasar interacts with two
companions located to the north and east (Heckman et al. 1984). It has
been imaged with HST by Bahcall et al. (1997) who concluded from
2D-modelling that the host is a spiral galaxy with
r1/2 = 10.8 kpc. They also fit the one-dimensional profile as
an elliptical with r1/2 = 24.5 kpc. The best 2D fit
by McLure et al. (1999)
is a large elliptical with r1/2 = 23 kpc. We
find that the best fit to the H-band profile is a r1/4 law with a
scalelength of 2.9 arcsec or 6h-1 kpc (see Fig. 3).
PG 1402+261 - TON 0182 Bahcall et al. (1997) determined from HST
imaging that the host of this quasar is a bright spiral galaxy with
prominent H II regions located along the spiral arms. The CFHT
image shows that there is a bright elongation extending
2 arcsec away from the nucleus on both sides of it in the
NW-SE direction, probably tracing a strong bar also seen
as a bump in the surface brightness profile, see Fig. 4.
B2 1425+26 - TON 0202 This quasar is hosted by a bright galaxy
probably in mild interaction with one of the three companions located
within 10 arcsec from the point source and listed by
Hutchings et al. (1984) and Block & Stockton (1991). Two of these
companions are detected in our H and K-band images. The host-galaxy
has been classified as elliptical by previous investigators
(Malkan 1984; Kirkhakos et al. 1999). However, the isophotes in our
image are distorted and an arc-like feature is clearly seen to the
west (see also the F555W image of 1999). Moreover, an
emission line nebulosity is detected by H,
H,
[O
III] emission up to 3 arcsec from the nucleus
(Boroson et al. 1984; Stockton & MacKenty 1987). Our fit shows that an elliptical galaxy
(or an early spiral) is a good description for the host-galaxy (see
Fig. 4).
3C 334 The host-galaxy has twisted isophotes. As noted by
Lehnert et al. (1999), the galaxy appears elongated approximatively in
the same direction as the radio structure (Hutchings et al. 1998). Note
that the arc-like structure seen to the south in the deconvolved image
was detected by its [OII] emission (Hes et al. 1996). Comparison
of the new AO corrected images with previous ground-based images
(e.g. Hes et al. 1996; Márquez et al. 1999) clearly shows the power
of AO technique. We note that three Lyman-
absorptions are
detected in the quasar spectrum with
= 0.37, 1.00 and
0.21 Å at redshifts
= 0.5387, 0.5449 and 0.5491,
slightly smaller than the emission redshift
= 0.555
(Jannuzi et al. 1998). Spectroscopy of the companions is required to
determine if they are somehow associated with these absorption
systems. The profile fitting does not allow us to distinguish between a
disk-like or an elliptical host.
3C 351 We observed this field in the H and K-bands (see
Fig. 3). The two objects closest to the line of sight in
our image are situated 7 arcsec north-east of the quasar. One is a
bright spiral galaxy, well resolved in our images, the other is very
faint (it is barely visible in Fig. 3 but it is more
clearly seen in the deconvolved image) and could be a companion of the
former. Lanzetta et al. (1995) and Le Brun et al. (1996) have searched
the field around 3C 351 for galaxies responsible for Lyman-and C IV absorption observed by Bahcall et al. (1993) in the HST
spectrum of the quasar. There are two Lyman-
absorption lines
at
= 0.2216 and 0.2229, the former showing C IV
absorption as well. Le Brun et al. (1996) identified the metal line
system with a galaxy 710
h-150 kpc away from the line of
sight. The redshift of the two objects closest to the quasar have not been
determined however and it is possible that these two objects are responsible
for the absorptions. If they are at z = 0.222, the impact parameter
is of the order of only 15h-1 kpc. It is very important to
confirm this for our understanding of the nature of H I halos
around low-z galaxies. Boyce et al. (1998) detect a companion at
about 3.3 arcsec east to the quasar. We do not detect
this companion in any of our images although it is within the
possibilities of our imaging. This may indicate that the flux in the
HST-F702W filter is dominated by line emission or that the object is very
blue. As there is no strong intervening absorption in the spectrum of
the quasar, it is probable that this companion is at the redshift of
the quasar. Note that 3C 351 exhibits a strong associated system with
H I Lyman-,
C IV, N V and O VI
absorptions. The presence of the associated system does not seem to
be related to any other imaging property of the QSO.
B2 1721+34 The host-galaxy is detected up to 1.5 arcsec from the
central point-source in the H-band CFHT image (see
Fig. 3). The one-dimensional profile shows that a r1/4law fits the galaxy profile better than an exponential fit, but images
with better S/N ratio are needed to confirm this at higher
significance.
Figure 3: Images of objects described in Sect. 3.3 and Tables 1 and 2. They are K band images for PKS 1302-102 and PKS 2128-12 and H band images for the rest. The scale is marked with the horizontal bar in the top-left figure and is the same for all the figures | |
Open with DEXTER |
PG 2112+059 Two companions are detected within 10 arcsec of
the nucleus. An elliptical profile
better describes the light distribution of the host-galaxy.
PKS 2128-12 Disney et al. (1995) fitted the one dimensional
profile with a de Vaucouleurs law of radius = 37.4 kpc.
The images of this object have the poorest resolution in our sample.
We are therefore unable to determine a reliable fit to the host-galaxy
profile. We detect a companion at 7 arcsec north-east of the
nucleus.
Figure 4: Top: surface brightness profiles (dots and error bars) obtained with the undeconvolved images in the H-band (except for PKS 2128-12, for which we used the K-band image because of its better PSF). Fits by a disk (dashed line) and a r1/4 law (dotted-dashed line) are overplotted. Bottom: residuals from the disk (stars) and r1/4 (crosses) fits in percentage. The vertical lines are drawn at respectively two and three times the PSF FWHM | |
Open with DEXTER |
Figure 5: Illustrative examples of mocked images and their fits. The top panel is the surface brightness profiles with models overplotted (exponential disk as a solid line and r1/4 law as a dashed line). a) The input is a disk galaxy and the mocked image is better fitted by a disk host (the quasar contributes 25% to the total light). b) The input is a disk galaxy but the mocked image is equally well fitted by the two laws. The disk parameters are the same as is a), but the quasar contributes 70% to the total light. c) When the input is an elliptical galaxy, the mocked image is always better fitted by a r1/4 profile | |
Open with DEXTER |
The magnitudes of the hosts have been derived by integrating the r1/4 profiles for all the objects. Indeed, it can be seen in Fig. 4 that there is an excess of light at small radii compared to the disk profile for all objects. This suggests that in our sample, the disk galaxies also have a strong bulge and/or a strong bar. This is confirmed by the 2D luminosity distribution. Results are given in Table 2. We have also subtracted a scaled version of the most suitable PSF for each nucleus (imposing a non-negative profile in the center, see Márquez et al. 1999) The resulting host magnitudes, computed by integrating the PSF-subtracted images, are in good agreement with those obtained from the profile fitting (see Table 2).
In order to test our fitting procedure, we have generated images of model elliptical and disk galaxies with scale-lengths and effective surface brightness within the range derived from the data. The same orientation and axis ratio is given to all of them. A point source is added in the center of the galaxy to mimic the quasar. An appropriate amount of noise is added, and then the images are convolved with a typical observed PSF. The mock data images are analyzed in the same way as real data.
Figure 6: Test of the fitting procedure on disk-galaxies. We plot here the difference between the input and output magnitudes versus the input magnitude in the case the host-galaxy has the same magnitude as the quasar. A dot surrounded by a circle means that the fit by an r1/4 law is at least as good as the fit by an exponential disk | |
Open with DEXTER |
Figure 7: Test of the fitting procedure on elliptical galaxies. We plot here the difference between the input and output magnitudes versus the input magnitude. Triangles correspond to the systems where the QSO contributes to the total luminosity as the host-galaxy does, squares represent those cases in which the host luminosity is half that of the QSO | |
Open with DEXTER |
We first note that an elliptical galaxy is always recognized as an elliptical galaxy by the fitting procedure, whereas a disk-galaxy is better fitted by a r1/4 law when the unresolved point-source contributes more than half the total light. This is illustrated in Fig. 5. This means that, at least with data of similar quality to those presented here, the fraction of elliptical galaxies in the sample may be overpredicted. Going deeper, at least 0.5 to 1 mag, should help solve this problem as it is apparent that the distinction between spiral and elliptical profiles is easier when the galaxy is detected at larger distances from the central point-source.
It is interesting to note that the output magnitudes are brighter than the input in both cases, elliptical or disk galaxies (see Figs. 6 and 7). The reason for this is probably the difficulty in determining the extension of the PSF wings which, if not subtracted properly, will artificially increase the flux of the host-galaxy. In the case of spirals, the difference is as large as 0.6 mag when the contribution of the point-source is the same as the contribution of the host-galaxy (see Fig. 6). For the ellipticals, the difference is less but still important when the QSO dominates the total flux.
Note that the ratio between the QSO and the host-galaxy luminosities is expected to increase with redshift. The above bias tends to imply that host-galaxy luminosities could be overestimated.
Adaptive-optics imaging in the H and K bands has been used to study the morphology of QSO host-galaxies at low and intermediate redshifts (z < 0.6). We detect the host-galaxies in 11 out of 12 quasars, of which 5 are radio-quiet and 7 are radio-loud. The images, obtained under poor seeing conditions, and with the QSOs themselves as reference for the correction, have typical spatial resolution of FWHM 0.3 arcsec before deconvolution. In the best case, the deconvolved H-band image of PG 1700+514 (with a spatial resolution of 0.16 arcsec) reveals a wealth of detail on the companion and the host-galaxy, and is probably the best-quality image of this object thus far.
Four of the quasars in our sample have close companions and show obvious signs of interactions. The two-dimensional images of three of the host-galaxies unambiguously reveal bars and spiral arms. For the other objects, it is difficult to determine the host-galaxy morphology on the basis of one dimensional surface brightness fits alone.
We have simulated mocked images of host-galaxies, both spirals and ellipticals, and applied the same analysis as to the data. Disk hosts can be missed for small disk scalelengths and large QSO contributions. In this case, the host-galaxy can be misidentified as an elliptical galaxy. Elliptical galaxies are always recognized as such, but with a luminosity which can be overestimated by up to 0.5 mag. The reason for this is that the method used here tends to attribute some of the QSO light to the host. This is also the case for disk galaxies with a strong contribution of the unresolved component.
Acknowledgements
We acknowledge the suggestions by the referee, Alan Stockton, that helped to improve the presentation. I. Márquez acknowledges financial support from the Spanish Ministerio de Educación y Ciencia (EX94-8826734). This work is financed by DGICyT grants PB93-0139 and PB96-0921. Financial support to develop the present investigation has been obtained through the Junta de Andalucía TIC-114.