A&A 428, 817-821 (2004)
DOI: 10.1051/0004-6361:20041819
J. D. Kurk 1 - L. Pentericci 2 - R. A. Overzier 1 - H. J. A. Röttgering 1 - G. K. Miley 1
1 - Sterrewacht Leiden, PO Box 9513, 2300 RA, Leiden,
The Netherlands
2 -
Max-Planck-Institut für Astronomie, Königstuhl 17,
69117, Heidelberg, Germany
Received 10 August 2004 / Accepted 7 October 2004
Abstract
Radio galaxy PKS 1138-262 is a massive galaxy at z =
2.16, located in a dense environment. We have found an overdensity of Ly
emitting galaxies in this field, consistent with a proto-cluster
structure associated with the radio galaxy. Recently, we have
discovered forty candidate H
emitters by their excess near infrared
narrow band flux. Here, we present infrared spectroscopy of
nine of the brightest candidate H
emitters. All these candidates
show an emission line at the expected wavelength. The identification
of three of these lines with H
is confirmed by accompanying
[N II] emission. The spectra of the other candidates are
consistent with H
emission at
,
one being a QSO as
indicated by the broadness of its emission line. The velocity
dispersion of the emitters (360 km s-1) is significantly smaller
than that of the narrow band filter used for their selection (1600 km s-1). We therefore conclude that the emitters are associated
with the radio galaxy. The star formation rates (SFRs) deduced from
the H
flux are in the range 6-44
yr-1 and the SFR
density observed is 5-10 times higher than in the HDF-N at z =
2.23. The properties of the narrow emission lines indicate that the
emitters are powered by star formation and contain very young (<100 Myr) stellar populations with moderately high metallicities.
Key words: galaxies: active - galaxies: clusters: general - galaxies: evolution - cosmology: observations - cosmology: early Universe
High redshift clusters are prime subjects for the study of galaxy
formation and cosmology. The powerful radio galaxy PKS 1138-262 at z = 2.156appears to be the brightest galaxy in a high redshift cluster. We have
discovered an overdensity of Ly
emitters within 1.5 Mpc of 1138-262 (Kurk et al. 2000, Paper I). The redshifts of 14 emitters were
spectroscopically confirmed to be in the range
2.14 < z < 2.18(Pentericci et al. 2000, Paper II). In addition, we have carried out near
infrared imaging (Kurk et al. 2004, Paper III). The number of K band
galaxies and extremely red objects in this field is higher than in
blank fields. We found 40 objects with excess narrow band flux,
consistent with H
emission at
.
The surface density of
H
emitters increases towards the radio galaxy and their average Kmagnitude is lower and therefore their inferred stellar mass higher
than for the Ly
emitters. Here, we present infrared spectroscopy
of nine candidate H
emitters to confirm their redshift and
determine the velocity dispersion of the sample. We assume a flat
Universe with h0=0.65 and
.
![]() |
Figure 1:
The ten two dimensional spectra observed through four
slits. Skyline residuals are visible as vertical lines with a higher
noise level. Slit and spectrum number are indicated, as are the
object number from the H![]() ![]() |
Open with DEXTER |
With the aim of detecting H emitting galaxies in the proto-cluster
associated with PKS 1138-262, imaging in
and in a narrow band filter
(
,
)
was carried
out, employing two pointings covering a total field of 12.5 arcmin2. There are 40 objects with rest frame equivalent width
(
Å (see Paper III). From the list of 29 candidate
H
emitters within 1
3 from the radio galaxy, we
selected those with H
flux
erg cm-2 s-1. Furthermore, we selected those that were conveniently located
for placement in the slit for spectroscopic follow-up, which was
carried out with ISAAC at VLT Antu (UT1)
. The short wavelength camera
of ISAAC is equipped with a Rockwell Hawaii 10242 pixel Hg:Cd:Te
array which has a projected pixel scale of 0.147
.
We used the
medium resolution grating in second order resulting in a dispersion of 1.23 Å. The observations were carried out in the nights of March 23
and 25, 2002 under variable seeing, which was just below 1
for
most of the time. The
slit employed
resulted in a resolution of
.
During acquisition, the slit
was first positioned on a bright point source (a star or the radio
galaxy) and subsequently positioned at the midpoint between two
candidates, which was always within 32
.
We have employed four
slit positions, each targeted at two or three candidate emitters for 3.5 h (3.75 in one case). In total, nine candidates were observed,
one of which was included in two slits. Per slit, we obtained 14 (15)
frames of 15 min with offsets of 15, 18 or 20
in ABBA sequence with additional random jitter offsets up to 5, 2, or 1
,
respectively, to avoid recurrent registration of spectra on
bad pixels. Standard stars were observed with the same slit at a range
of airmasses during the nights to correct for telluric absorption and
to calibrate the data in flux. All observations were carried out at
airmass below 1.8.
Standard data reduction was carried out using pairwise frame subtraction, resulting in a final frame with positive spectra accompanied by negative spectra on both sides. Care was taken during setup to ensure that these negative spectra did not overlap with positive ones. One dimensional spectra of the candidate emitters were extracted from the positive two-dimensional spectra using the spatial profile of a standard star spectrum observed during that night and averaged per 3 pixels yielding bins of 3.7 Å. The spectral resolution is 7 Å, as measured from the FWHM of the skylines. The wavelength calibration is based on the OH skylines observed. The telluric standards observed at a range of airmasses show only small variations (less than a percent on average). An average absorption spectrum per night was used to correct for telluric absorption. There are no spectrophotometric flux standards in the infrared. We have therefore used one of the telluric standards (Hip043868) with spectral type B1. Such stars have a featureless spectrum in this wavelength region given by a black body curve at T = 25 500 K. The curve was normalized to the K magnitude of the star and subsequently used to divide the extracted and absorption corrected spectrum to obtain the flux calibration. The two dimensional spectra of all objects are shown in Fig. 1.
All candidates observed show an emission line, which means that our
selection was 100% efficient. We have fit Gaussian curves to the
emission lines applying a least squares method (see Fig. 2) in order to determine their central wavelength,
deconvolved FWHM (if resolved), flux and EW0 (Table 1).
Also presented in this table is the SFR derived from the H and UV luminosities (see discussion in Paper III), which is in the range 6-44
yr-1. We do not detect continuum emission in most
of the spectra, except in the co-added spectrum of candidate 131, 215
and the serendipitous object in slit one (
,
and
erg cm-2 s-1 Å-1, respectively), which
compare favourably with the line subtracted broad band fluxes measured
by imaging in Paper III (0.6, 1.0 and
erg cm-2 s-1 Å-1, respectively). The EW0 is therefore based on the
line flux measured in the spectra and the broad band magnitude
measured on the images.
The spectra of objects 131, 183 and 229 were fit by two Gaussians for
which the relative centers were fixed as for H and
[N II]
Å. Candidate 329 shows two emission peaks only
Å apart. Both lines could be H
emission from one
galaxy with two components, seperated by
km s-1 in velocity.
We have considered the possibility that the lines are due to
[O II] emission at
z = 4.5674 which would have a separation
of 15.0 Å. The emission line ratio of the supposed
[O II]
would be 0.7 which implies an
electron density
cm-3, normally only observed in the
central parts of nebulae (Osterbrock 1989). In addition, a faint (B =
27.0) counterpart in the B band, sampling a wavelength range below 912 Å for z = 4.6, makes the identification with [O II] improbable. Candidate 131 was included in two slits, with a position
angle difference of 40
.
The two fits to the spectra of
emitter 131 indicate a velocity difference of 180 km s-1 which can be
explained by the fact that different regions of the galaxy have been
sampled.
For how many of the emission lines can we be sure that the
identification with H at
is correct? For the three
objects with confirming [N II] lines, we can be certain. For the QSO object 215, [N II] is blended with the very broad H
line and
impossible to discern. Given a [N II]/H
line ratio of 1/3, we do not
expect to detect the [N II] line for objects 79, 207, 284, 329 above
the noise and the identification with H
is therefore consistent but
not 100% certain. An identification with H
for object 144 can only
be true if the [N II]/H
ratio is <1/6 (object 183 has an observed
ratio of 1/5.5). An alternative identification with
[O III]
5007 Å is improbable, as we do not detect its
counterpart [O III]
4958 Å at 2.045
.
We consider
this therefore a probable H
identification.
![]() |
Figure 2:
One dimensional spectra of the eight narrow line emitters (histograms)
with fits overlayed (solid lines). Object 329 has also an
alternative fit for [O II]. Units are in 10-18 erg cm-2 s-1 Å-1 and ![]() ![]() |
Open with DEXTER |
At least three of the emission lines are spatially extended. In two
cases, we detect unordered velocity structure, but the morphology of
the two dimensional spectrum of candidate 229 resembles a rotation
curve. A fit to this structure results in a rotation velocity of 50 km s-1 at 6 kpc radius, implying a dynamical mass of
.
Table 1: Properties of the observed emission lines.
In Fig. 3 we show the redshift distribution of the 9 emitters, assuming that all emission lines can be identified as H.
For this plot, we used the average redshift of the two lines of object 329. Also plotted is the sensitivity curve of the narrow band filter
used to select the candidates. A random distribution of emitter
redshifts would follow this curve. We have run a Monte Carlo
simulation with 10 000 realizations of a randomly sampled distribution
of emitters given this filter as selection criterium. Although the
mean of the measured distribution is only 0.18
away from the
mean of a random sample, its dispersion is much smaller, deviating by 1.75
from a random sample. The probability that the redshifts we measure are drawn from a random distribution is therefore 8%. Note
that the redshift of the radio galaxy is 2.156 (Röttgering et al. 1997), very
close to the mean of the selection filter (2.152), and to the mean of
the measured redshift distribution (2.154). The distribution is
consistent with a group of H
emitters associated with the radio
galaxy. The velocity dispersion of this group is 360 km s-1 (using the gapper sigma method, Beers et al. 1990), while the virial
radius of the nine emitters is 0.45 Mpc, implying a virial mass of
(assuming that all lines can be identified
with H
and taking the mean of the two redshifts for object 131).
This mass is merely illustrative as at this redshift it is improbable
that the structure is virialized. The velocity dispersion of the H
emitters is smaller than the velocity dispersion of the confirmed
Ly
emitters, both for the complete sample (1050 km s-1) and the nine
within the solid angle of the two ISAAC fields (760 km s-1, shaded part
of histogram in Fig. 3). There is no evidence for a
bimodal redshift distribution as observed for the Ly
emitters.
We can construct a complete sample out of the spectroscopic sample by
excluding the two objects with the lowest H flux and including the
radio galaxy. This collection represents all candidate H
emitters
with
erg cm-2 s-1 within 1
3
from the radio galaxy. The FWHM of the narrow band filter (
2.134 < z
< 2.174) and the solid angle given above define a comoving volume of 815 Mpc3, resulting in a volume density of 0.010 Mpc-3, which
is a factor four higher than the density of confirmed Ly
emitters
in this field. All star forming objects detected have line fluxes
lower than the high redshift H
surveys discussed in Paper III, but
we can compare the SFR density to the density at z = 2.23 derived
from H
emission in the HDF-N as measured by Iwamuro et al. (2000).
Following their cosmology (h0=0.5, q0=0.5) and procedure to
correct for the part of the H
luminosity function below the
detection limit, we obtain a SFR density of 0.48
yr-1 Mpc-3. This is 10 (5) times higher than the (reddening
corrected) value obtained for the HDF-N. Using the redshift range
defined by the H
emitters (
2.146 < z < 2.164) results in values
that are larger by a factor of two. Likewise, smaller SFR densities
would result if some of the H
lines have been misidentified.
The properties of the detected emission lines provide information
about the physical conditions in the galaxies. The FWHM of the narrow
nebular emission lines detected are in the range 40-360 km s-1 with
an average of 190 km s-1. These values are comparable to those found for
LBGs at
by Pettini et al. (2001). The [N II]/H
ratio can
be used to distinguish narrow-line active galaxies from H II region-like galaxies (Veilleux & Osterbrock 1987). The three emitters with detected [N II] have
[N II]/H
)
< -0.42, which
puts them among the star forming galaxies. In the absence of shock
excitation, the [N II]/H
ratio can also be used as
metallicity indicator. Using the empirical relation calibrated by
Denicoló et al. (2002), the average ratio of the three emitters implies
O/H
.
This value is comparable to the broad range of
values obtained for present-day spiral galaxies (van Zee et al. 1998). The EW0 of some detected narrow lines are surprisingly high, up to 1350 Å. This can be explained by very young stellar populations
where the continuum radiation around 6000 Å is still very weak.
EW0 values between 200 and 330 Å imply an age <100 Myr
(Leitherer et al. 1999). The moderately high metallicities found for the objects
with detected [N II] emission, however, require that the
galaxies are near the end of the star formation event. This
requirement seems to indicate that these emitters have undergone a
very similar evolution.
![]() |
Figure 3:
Redshift histogram for the H![]() ![]() ![]() |
Open with DEXTER |
Acknowledgements
We are grateful to the ESO VLT staff for excellent support during the observing run. We acknowledge fruitful discussions with B. Venemans and S. di Serego Alighieri. Comments of the anonymous referee have also helped to improve the manuscript. This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. We have also made use of NASA's Astrophysics Data System Bibliographic Services.