A&A 367, 51-58 (2001)
DOI: 10.1051/0004-6361:20000418
T. Yamada1 - S. Yamazaki1 - M. Hattori1 - G. Soucail2 - J.-P. Kneib2
1 - Astronomical Institute, Tohoku University, Aoba-ku,
Sendai 980-77, Japan
2 -
OMP, 14 Av. E. Belin, 31400 Toulouse, France
Received 5 August 1999 / Accepted 27 November 2000
Abstract
We present new high signal-to-noise ratio spectra of the components B and C of the gravitational lensing system MG 2016+112. We show that image C displays strong emission lines of Ly,
NV, CIV, HeII, and CIII] redshifted to z=3.27, similar to images A and B. We examine the emission-line flux ratios in order to put constraints on the lens models as well as to investigate the intrinsic nature of MG 2016+112. The observed line ratios of B and C are consistent with those expected in the simple photo-ionization models for narrow-line region of active galactic nuclei (AGN) except for the enhanced NV lines. The line ratios difference of components B and C can be interpreted as a difference in ionization parameters. This result is consistent with lens model prediction that C is a fold image of a slightly outer part of the nucleus. MG 2016+112 is known to be very unique among the high-redshift AGN; it is neither an ordinary broad-line quasar nor a powerful radio galaxy as indicated by the width and flux ratio of the emission lines.
Together with other observed properties discussed in literature, we argue that MG 2016+112 may be an obscured luminous radio-quiet quasar.
Key words: galaxies: active - galaxies: quasars: emission lines - galaxies: quasars: individual: MG 2016+112
MG 2016+112 is one of the first discovered gravitational-lens system. Lawrence et al. (1984) observed three distinct radio sources, A, B, and C. Further higher resolution radio map resolved C into multiple components (Garrett et al. 1994, 1996). At optical wavelength, the image A and B are point-sources but the image C is a fainter resolved object (Lawrence et al. 1984; Schneider et al. 1985, 1986). Spectroscopic observations revealed that A and B show very similar spectra dominated by strong ultra-violet (UV) emission-lines redshifted to z=3.27 (Lawrence et al. 1984; Schneider et al. 1987). Although the optical spectrum of C has not been fully published so far, it has always been considered at the same redshift as A and B after its detection in the redshifted Ly
narrow-band data (Schneider et al. 1986).
The MG 2016+112 lens system is enigmatic for 3 reasons. First, the lensed object itself is very unique among known high-redshift galaxies (Lawrence et al. 1984). It has been conventionally called a "quasar" since it is a fairly luminous point-like object. However, A and B show only narrow ultra-violet (UV) emission lines and thus MG 2016+112 is not an ordinary broad-line quasar. Although strong narrow emission lines are typically seen for high-redshift powerful radio galaxies (HzPRGs), MG 2016+112 appears somewhat different from the known HzPRGs. Indeed, while HzPRGs are typically extended in optical images and have lobe-dominated radio structures with a scale of at least a few tens of kpc, A and B show only point-like features even in high-resolution optical and radio images. It is important to investigate the true nature of the lensed object.
Second, image C is a very complex object: at optical and near-infrared (NIR) wavelength, C is resolved and has an arc-like morphology. Radio-to-optical flux ratio of C is several times larger than those of A and B and cannot be explained by variability and time delay. At radio wavelength, C is resolved into two components, C1 and C2. While A, B and C2 (sometimes referred as C
in literatures) are point-like objects even at 15-mas resolution, C1 has been further resolved into three chain-like components (Garrett et al. 1996). Radio spectral shape of C2 is similar to those of A and B (
)
but C1 has significantly flatter one (
). Whether C1 is a radio galaxy at different redshift or a lensed image of the outer structure of the radio source at z=3.27 is still in question.
Finally the nature of the lens producing this multiple image system has been enigmatic. Deep optical and NIR images have detected a red galaxy D amid images A, B and C (Schneider et al. 1986; Langston et al. 1991; Lawrence et al. 1993). Galaxy D seems to correspond to an evolved giant elliptical galaxy at
.
However, the image separation requires a mass-to-light ratio for D much larger than typical one. While lens models have assumed the existence of a high-redshift massive cluster as an additional source of lensing mass (Narasimha et al. 1984, 1987, 1989; Nair & Garrett 1997; Langston et al. 1991), at first no signature of such a cluster has been seen in the optical and NIR observations (Schneider et al. 1987; Langston et al. 1991).
Hattori et al. (1997) detected an extended X-ray emission in the direction of MG 2016+112, possibly emitted by the hot gas in the lensing cluster of galaxy. A strong emission-line-like feature consistent with iron lines redshifted to z
was detected. This discovery could in principle solve the "dark-lens" problem for MG 2016+112, but opened another question of a "dark cluster", namely, the lack of optical counterpart of the X-ray hot gas. Very recently, Soucail et al. (2000) and Kneib et al. (1997) spectroscopically detected several galaxies at
and Benítez et al. (1999) observed a possible color-magnitude sequence of faint red galaxies in this field arguing for the existence of a distant
massive structure at
.
We have obtained new spectra of components B, C, D of the MG 2016+112 system. In this paper, we analyze the emission lines of image B and C in order not only to investigate the intrinsic nature of MG 2016+112 but also to understand the component C in view of the lens-model prediction. Observations and data reduction are described in Sect. 2. In Sect. 3, we examine the observed emission line properties. The line flux ratios are compared with the prediction of the photo-ionization models as well as those of other high-redshift objects in various categories. In Sect. 4, the nature of image C is discussed in the context of lens models. We then argue that MG 2016+112 may be a radio-quiet obscured quasar based on the results obtained in Sect. 3 as well as those discussed in literatures. Throughout this paper we use H0=50 kms-1/Mpc,
and
.
Spectroscopic observations were made with the Canada-France-Hawaii Telescope (CFHT) using the Subarcsecond Imaging Spectrograph (SIS; Le Févre et al. 1994) in August 1997. The Stis2
CCD with 21
m pixel and the R150 grating were used and the resultant dispersion is 2.88 Å per pixel. The wavelength coverage extends from 4500 to 9000 Å. Two multi-slit masks were used to obtain spectra of candidate cluster galaxies in the field. Each slit has 7.8-arcsec-long and 0.78-arcsec-wide aperture, which gives
18 Å instrumental resolution. MG 2016+112 B, B1, and D were observed within one slit of the first mask (Mask1) and C with a slit in the second one (Mask2). In total, 6 and 5.8 hours exposures were obtained for Mask1 and Mask2, respectively.
The data was pre-reduced using standard IRAF tasks. We then follow the reduction procedure of multi-slit data described in Le Févre et al. (1995). Wavelength calibration was done with the arc-line spectra taken at the observation. Typical internal error of wavelength determination is
0.3 Å. The flux calibration was done using the spectroscopic standard star GD248 (Oke 1990). Note that the flux calibration becomes more uncertain above 7500-8000 Å because no order-separating filter was used and the flux is contaminated by the UV and blue contribution of the second order filter. We tried to carefully take into account this effect but flux calibration may remain uncertain above 8000 Å.
Figures 1 and 2 show the obtained spectra of image B and C smoothed with a 5
pixel Gaussian filter. Strong UV emission lines of Ly,
NV
Å, CIV
Å, HeII
Å, and CIII]
Å are seen in both of the spectra. SiIV and OV lines at 1400 Å are clearly seen in the B spectrum but only marginal in C. The redshift of C is the same as B within the uncertainties. There is no sign of contamination by a system at a different redshift.
![]() |
Figure 1: Observed spectra of MG 2016+112 B |
Open with DEXTER |
![]() |
Figure 2: Observed spectra of MG 2016+112 C |
Open with DEXTER |
We now describe the observed properties of the emission lines of image B and C. The velocity width (FWHM) corrected for instrumental resolution and the relative flux values measured by Gaussian fitting procedure are listed in Tables 1 and 2, respectively. The flux values are normalized to CIV lines. We corrected for the reddening by the Galaxy using the value of the extinction in this region, AV=0.67 and AI=0.36, derived by Benítez et al. (1999). No correction for the internal reddening was applied since the rest-frame wavelength of the emission lines concerned here are rather close. We treat the unresolved doublet lines as a single line.
MG 2016+112 A![]() |
MG 2016+112 B | MG 2016+112 C | F10214![]() |
SMM 02399![]() |
|
Ly![]() |
<1000 | 460 | 630 | 900 | 1850 |
N V | <1600 | 880 | 840 | 1700 | 1790 |
C IV | <1100 | 580 | 580 | 1200 | 1560 |
He II | <1200 | 460 | 440 | 1150 | 2800 |
C III] | -- | 750 | 1100 | 7700 | |
C III] narrow | -- | 360 | 380 | 1000 | |
C III] broad | -- | 1750 | 2550 | 3700 |
(a) FWHM in km s-1. (b) From Lawrence et al. (1984). (c) From Searjeant et al. (1998). (d) From Ivison et al. (1998). |
MG 2016+112 B | MG 2016+112 C | F10214![]() |
SMM 02399![]() |
HzPRG![]() |
Quasar BLR![]() |
Seyfert 2![]() |
|
Ly![]() |
2.73 (0.09) | 3.81 (0.43) | 0.53 | 3.39 | 2.17 | 8.52 | 4.52 |
N V | 0.66 (0.07) | 1.19 (0.48) | 0.74 | 2.20 | 0.57 | 0.42 | |
C IV | 1.00 | 1.00 | 1.00 | 1.00 | 1.00 | 1.00 | 1.00 |
He II | 0.28 (0.03) | 0.55 (0.10) | 0.53 | 0.36 | 0.12 | 0.87 | 0.17 |
C III] | 0.40 (0.04) | 1.32 (0.27) | 2.70 | 0.45 | 0.49 | 0.46 | |
C III] narrow | 0.17 (0.02) | 0.51 (0.10) | 0.34 | ||||
C III] broad | 0.29 (0.06) | 1.25 (0.51) | 0.28 |
(a) Normalized to CIV line. Numbers in parenthesis show the error values. (b) From Searjeant et al. (1998). (c) Evaluated by us from the values of equivalent width in Ivison et al. (1998). (d) From McCarthy et al. (1993). |
The velocity width of the lines are between 450-900 kms-1. These lines are much narrower than quasar broad lines which have typical width of 5000-10000 kms-1 but as narrow as those of HzPRG, 500-1000 kms-1 (e.g., McCarthy 1993). In Table 1, we compare the emission-line width to the two infrared-selected gravitationally-lensed type-2 AGN at high redshift, namely 10214 (e.g., Rowan-Robinson et al. 1991) and SMM 02399-0136 (Ivison et al. 1998).
CIII] line of both B and C seem to have broad wings. Figures 3a,b shows the results of the two-component Gaussian fitting of the CIII] lines for B and C and 3c and 3d shows those with single-component for comparison. They cannot be perfectly fitted with a single component but can be better fitted with a narrow plus broad components as shown in the figures. It seems strange that the CIII] line has a broad component while the CIV line, which comes from the gas at higher ionization stage and may be closer to the central engine, shows only a narrow one. It may be due to the extinction of the broad component of the CIV line. Another possibility is that the wings are the lines of other ions contaminating this spectral range.
The broad CIII] feature is not a unique characteristics of MG 2016+112 lensed object. 10214 (Searjeant et al. 1997) and SMM 02399-0136 (Ivison et al. 1998) also show similar property. Searjeant et al. (1997) considered the possibility of contamination by SiIII] lines to the blue-wing feature of the CIII] line in the spectrum of 10214 but the model cannot fully explain the observed feature (Searjeant et al. 1997). For the case of MG 2016+112, the broad component seems to extend both side of the narrow component (Figs. 3a,b) and thus it is also difficult to be explained by the contamination of SiIII] lines. It is interesting that those two gravitationally-lensed type-2 AGN as well as MG 2016+112 show evidence of a broad CIII] line.
Next, we investigate the line flux ratios. There are significant differences between line ratios of component B and C. While the HeII/CIV ratio of image B is 0.28, the ratio of C is 0.55. Also, the CIII]/CIV ratio of B and C is 0.17 and 0.51, respectively if we consider the narrow components. If we fit the lines with single Gaussian component, then the ratio is 0.40 and 1.32. The emission-line gas in image B seems to be at higher ionization stage than in C. In Fig 4, we compare these line ratios with those predicted by photoionization models (CLOUDY90, Ferland 1988).
We examined the cases with hydrogen density
and 1000 cm-3, power-low ionization continuum with
and 1.4, and a range of ionization parameter,
=10-2.5-10-1 and solar-abundance. The obtained line ratios are fairly consistent with the typical photo-ionization models considered for narrow-line region of AGN. The differences in line ratios between B and C may be interpreted as the difference in ionization parameters.
![]() |
Figure 3: Results of the two-components Gaussian fitting for the C III] lines of component B (panel a)) and C (panel b)) and of the one-component ones (panels c) and d)). The LDECONV task developed by Michitoshi Yoshida was used for the fitting |
Open with DEXTER |
There may be an effect of the reddening by dust in the object on the observed line ratio. According to Calzetti's reddening formula for starburst galaxies (Sawicki & Yee 1998), the reddening of
E(B-V) = 1 shifts the observed values with log(CIV/HeII) = -0.12
and
log(CIII]/HeII) = +0.30, respectively, which is shown by the arrow in Fig. 4. Even if there is a fairly large amount of reddening,
,
the discussion we give here is not so much affected.
There are weak but fairly significant NV lines in the spectra of B and C. In Fig. 5 we plotted the observed NV line ratios as well as those predicted by photoionization models. Clearly, photoionizatin models that can explain the line ratios at lower ionization state are not consistent with the observed NV flux. The observed NV lines are more than several times stronger than the predicted value. Reddening correction moves the points further away from the models. The origin of NV line may be different from those of other lines since the ionization potential of N++++ is 77.5 eV which is significantly higher than those of CIV (47.9 eV), CIII] (24.4 eV), and HeII (24.6 eV). Indeed, the observed line width of the NV line is larger than others (see Table 1). The enhancement of the NV line may also be due to the nitrogen over abundance. The excess of the NV line intensity is not special for MG 2016+112 B; Hamman & Ferland (1993) also compare their photoionization models with the observed lines of broad-line region (BLR) of quasars and found a large nitrogen over abundance. There is also similar nitrogen problem for the optical NII lines observed in the spectra of local AGN (e.g., Osterbrook 1986).
The fast shock models (Mouri & Taniguchi, private communication) can provide the NV/CIII] and NV/HeII ratios that match the observed values but then predict too strong CIV lines. The line ratios observed in image B and C are not perfectly understood with simple photoionization or shock models.
In Fig. 6, we also plot typical line ratios of various types of high-redshift AGN taken from McCarthy (1993), ultra-steep spectrum HzPRGs compiled by Röettgering et al. (1997), and the
infrared-selected type-2 quasars, 10214 (Serjeant et al. 1997). Many of the HzPRGs in Röettgering et al. (1997) are distributed at around log(CIV/HeII)
and log(CIII]/HeII)
,
which is consistent with the average flux ratio given by McCarthy (1993). The line ratios of image B is very different from typical HzPRGs. The line ratios of C and 10214 lies near the edge of the distribution of line ratios of HzPRGs.
We also plotted the line ratio of the broad lines observed in quasar spectra compiled by Baldwin et al. (1979). They are very different from both HzPRGs and image B and C and can be explained by the photoionization models with high density such as 109-10 cm-3.
Several lens models have been proposed for MG 2016+112 (Narasimha et al. 1984, 1987, 1989; Langston et al. 1991; Nair & Garrett 1997; Benítez et al. 1999). These models are constructed to match the optical and/or radio data such as image positions and flux ratios of the various components. The results of our spectroscopic observations can be used to put further constraints on the lens model.
As discussed in Sect. 3, the difference of CIII], CIV, and HeII line ratios in image B and C can be interpreted as a difference in ionization degree. Photoionization model predict smaller ionization parameter for C, which is a natural consequence if the component C is dominated by the light of a region at 1.5-2 times larger radius from the nucleus than the source region of the component B (assuming similar density).
![]() |
Figure 4:
Flux
ratios of the emission lines seen in the spectrum of image B and
C. Squares show those using only the narrow component of the
C III] line and circles those using the C III] flux obtained with
single-component fitting. The grids show the line flux ratio
predicted by photo-ionization models calculated using
CLOUDY90 (Ferland et al. 1988). The cases of hydrogen density of 100
cm-3 and 1000 cm-3 and input power-law spectra with
energy index -1.0 and -1.4 are plotted for ranges of
ionization parameter log![]() |
Open with DEXTER |
It is unlikely that the differences in line ratios of B and C is due to the contamination by a possible radio galaxy at the similar redshift assumed as the counterpart of the bright flat-spectrum radio component C1. If the high CIII]/CIV ratio observed in the spectrum of image C is due to such a contamination, for example, the assumed radio galaxy should have CIII]/CIV ratio larger than 1. Radio galaxies rarely show such high CIII]/CIV ratio.
Our spectroscopic results thus support lens models like the ones proposed by Langston et al. (1991) and Benítez et al. (1999). In these models, the AGN and the narrow-emission-line region is located somewhat outside the fold caustic in the source plane, which results in forming the images A and B. The outer regions of the source extends over the caustic and the region which is very close to the caustic is largely amplified and form an arc-like lensed images at the position of component C. Radio sources at position C are interpreted as a lensed image of a "jet" which extends inside the caustic. The counter images of the component C will exist at position A and B but is mostly swamped by the AGN light due to a smaller amplification factor (compared to position C).
![]() |
Figure 5: Same as Fig. 4, but for N V line |
Open with DEXTER |
The observed radio flux density of image B at 1.47 GHz is 61.7 mJy (Lawrence et al. 1984), which corresponds to a luminosity density
ergs-1 Hz-1
and
ergs-1 Hz-1
at 2.7 and 8.4 GHz, respectively, with a luminosity distance of 26.5 Gpc and a spectral index
(Lawrence et al. 1984).
is a gravitational lensing amplification factor. According to the recent lens model by Benítez et al. (1999), the amplification factor of image B is estimated to be
6.
The radio power of MG 2016+112 is fairly large even corrected from the lensing amplification (e.g., Lawrence et al. 1984). Danese et al. (1987) obtained radio
luminosity functions of at 2.4 GHz for various type of objects at .
Even the most luminous local objects have radio power
1032 ergs-1 Hz-1, which is more than a order of magnitude fainter than MG 2016+112. Dunlop & Peacock (1990) evaluated radio luminosity function of the radio-selected quasars and radio galaxies. Image B has radio power as large as those of very luminous radio sources at
which are definitely categorized as "quasars" or "powerful radio galaxies" (
1033-1034 erg s-1 Hz-1 at 2.7 GHz). Thus it is very certain that the object harbors AGN as powerful as luminous as quasars.
Bischef & Becker (1998) investigated radio emissions for a sample of 4079 known quasars based on the NVSS radio catalog. These quasars are selected from the Verón-Cetty & Verón (VCV) catalog and constitute the largest compilation so far to study the radio properties of radio- and optically-selected quasars based on homogeneous radio observations in one frequency. In their Fig. 4, they presented a distribution of radio power at 8.4 GHz with redshift. The bimodal distribution of the radio power is evident at least to .
At higher redshift, the number of known radio-luminous quasars is too small to draw firm conclusion, but the tendency seems to hold. The radio power of MG 2016+112 at 8.4 GHz is likely to lie on the extension of the lower sequence, if we adopt an amplification factor of
10. In Fig 7 we plot the observed and amplification-corrected radio power of MG 2016+112 B superposed on the figure of Bischef & Becker's (their Fig. 4).
![]() |
Figure 6: Line ratios for various kinds of high-redshift AGN. The solid squares are the average line ratios given in McCarthy (1993) and the open squares are those of USS HzPRGs in Röttgering et al. (1997). Triangles are the infrared-selected type-2 AGN, 10214, using the narrow-component (filled) and total (open) C III] flux. The ratios of broad lines quasar quoted from Baldwin (1979) are shown as crosses |
Open with DEXTER |
Radio loudness of AGN is conventionally defined by the radio to optical (or ultra violet) flux ratio. Bischef & Becker (1998) also obtained the distribution of the radio loudness for VCV quasars. According to their definition, radio-loud quasars have log
larger
than 1. The upper sequence in Fig. 4 of Bischef & Becker (1998) corresponds to radio-loud quasars and the lower to radio-quiet quasars. It is difficult, however, to evaluate the radio loudness of MG 2016+112 by using the radio to optical flux ratio since it is an obscured narrow-line object. We may see only a scattered light of the nucleus in the optical wavelength and the contamination by light from the host galaxy must be relatively large.
The nominal value, log
/
LB = 3.4,
thus does not mean that this object is a radio-loud AGN. Indeed, as mentioned
in the previous subsection, the observed radio power of MG 2016+112 after lensing amplification correction seems too faint to be classified as radio-loud, which suggests that this object is a radio-quiet AGN.
It is not surprising that the scattered component is more than hundred times smaller than the intrinsic luminosity and the intrinsic log
/LB value of the object can be
1.
![]() |
Figure 7: The observed (thick open circle) and the lensing amplification corrected (filled circle) radio luminosities of MG 2016+112 superposed in Fig. 4 of Bischof & Becker (1990) (scanned by us) which shows the distribution of radio luminosity at 8.4 GHz vs redshift for VCV quasars. Crosses are values for the broad lines of type-1 quasars |
Open with DEXTER |
How typical the radio power is if MG 2016+112 is a radio-quiet AGN? Figure 7 suggests that it may be one of the most radio active object among the radio-quiet quasars at this redshift. Kukula et al. (1998) recently investigated the correlation between radio and optical luminosity of a sample of nearby radio-quiet quasars and Seyfert galaxies
(z < 0.2). The most luminous nearby quasars with
have log
ergs-1 Hz-1. If the correlation holds for high-redshift quasars which are typically 50-100 times more luminous than nearby AGNs, MG 2016+112 (
erg s-1 Hz-1) may be one of the most luminous quasars in optical wavelength, too, even if there is some radio excess.
HST NICMOS observations (Falco et al.) revealed that the rest-frame optical light of MG 2016+112 A and B is dominated by the point sources. The resolution limit of NICMOS is 0.15 arcsec. If we adopt the lens model of Benítez et al. (1999), the size of the unresolved lensed object must be smaller than 0.03 arcsec, which corresponds to
200 pc at z=3.27. It may not be surprising that A and B are not resolved if we observe only the obscured nucleus or its scattered light which may come from inside the narrow-line region.
Many of HzPRGs show resolved faint extension in NIR images which may be star lights of their host galaxies. On the other hand, the known type-2 AGN or obscured quasars at 1-2.5 are not resolved or only marginally resolved (Ohta et al. 1995; Almani et al. 1995; Ivison et al. 1998). There could be some difference between observable properties of host galaxies of HzPRGs and radio-quiet type-2 AGN or obscured quasars.
The radio images of A and B are also point-like with
15 mas resolution (Garrett et al. 1996). 15 mas corresponds to about 40 pc in the physical scale at z=3.27.
Kukula et al. (1998) presented the 1.4 GHz maps of the nearby radio-quiet quasars with
0.5 arcsec resolution. They found a significant fraction of the radio emission in radio-quiet quasars originates in a compact nuclear source directly associated with the quasar. Therefore the point-like morphology of the radio emission of MG 2016+112 is not surprising if it is a radio-quiet quasar. At the same time, some outer structures are also seen in the radio maps of the nearby radio-quiet quasars. The maximum extent of the radio emission of resolved sources is typically a few kpc. It is possible that such a structure associated with the lensed object extends over the diamond caustic of the source plane to form the strongly amplified radio image at the position of component C.
Since only the narrow emission lines are observed in UV wavelength, MG 2016+112 is very likely to be an obscured luminous AGN. Although it is need to see the properties of lines in the rest-frame optical wavelength redshifted to near infrared before concluding whether it is really a luminous analogue of the Seyfert 2 galaxies (namely, type-2 quasars) or just a partially-obscured normal quasars, it is interesting to summarize the possible type-2 quasar nature of MG 2016+112 since the number of known radio-quiet luminous type-2 quasars and candidates are still very limited. Only several examples of type-2 AGN or candidates at high redshift which are serendipitously discovered in far-infrared, sub-mm, and X-ray source surveys. If MG 2016+112 at z=3.27 is really a type-2 quasar, it is the highest-redshift one identified so far and thus provides a unique and important example for further studies of type-2 quasars. The properties of MG 20116+112 are summarized as follows.
(1) The central engine of MG 2016+112 must be obscured since only the narrow emission lines are observed. (2) The radio and optical morphologies are both compact, which is not compatible with typical properties of HzPRGs but common for radio-quiet AGNs. (3) The radio power is much smaller than typical radio-loud quasars at
but consistent with luminous radio-quiet quasar. (4) Finally, the emission line flux ratios of image B are not compatible with typical HzPRGs.
We thus consider that MG 2016+112 is not a typical powerful radio galaxies but may be naturally classified as a rare example of a high-redshift radio-quiet type-2 quasar.
Since 10214 and SMM 02399-0136 are detected in CO in sub-millimeter (Rowan-Robinson et al. 1993; Ivison et al. 1998; Frayer et al. 1998), it will be interesting to investigate the cold gas and dust properties of MG 2016+112. As expected from the obscuration of the nucleus, large dust content may be a general property of type-2 quasars.
Acknowledgements
We thank Drs. H. Mouri and Y. Taniguchi for kindly providing us the results of the shock-model calculations. This research was partially supported by grants-in-aid for scientific research of the Japanese Ministry of Education, Science, Sports and Culture (09740168, 07055044). JPK thanks CNRS for support and the Yamada Science Foundation for fruitful visits in Japan. MH also thanks the financial supports of Yamada Science Foundation.