A&A 440, 61-66 (2005)
DOI: 10.1051/0004-6361:20034109
D. Rigopoulou1,2 - W. D. Vacca2,3 - S. Berta4 - A. Franceschini4 - H. Aussel5
1 - Astrophysics, University of Oxford, Denys Wilkinson Building, Keble Road, Oxford OX1 3RH, UK
2 - Max-Planck-Institut für extraterrestrische Physik (MPE), Postfach 1312, Garching 85741, Germany
3 - USRA, NASA Ames Research Center, MS 144-2 Moffett Field, CA 94035-1000, USA
4 - Dipartimento di Astronomia, Universita' di Padova, Vicolo dell'Osservatorio 2, 35122 Padova, Italy
5 - Institute for Astronomy, 2680 Woodlawn Drive, Honolulu, Hawaii 96822, USA
Received 24 July 2003 / Accepted 18 December 2004
Abstract
We present low resolution multi-object spectroscopy of an I-band magnitude limited (
23-23.5) sample of galaxies located in an area centered on the Hubble Deep Field-South (HDFS). The observations were obtained using the Focal Reducer/low dispersion Spectrograph (FORS) on the ESO Very Large Telescope. Thirty-two primary spectroscopic targets in the HST-WFPC2 HDFS were supplemented with galaxies detected in the Infrared Space
Observatory's survey of the HDFS and the ESO Imaging Deep Survey to comprise a sample of 100 galaxies for spectroscopic observations. Based on detections of several emission lines, such as
[OII]
3727, H
and [OIII]
5007, or of other spectroscopic features, we
measured accurate redshifts for 50 objects in the central HDFS and flanking fields. The redshift range of the current sample of galaxies is 0.6-1.2, with a median redshift of 1.13 (at I
23.5 not corrected for completeness). The sample is dominated by starburst galaxies with only a small fraction of ellipticals (
10%). For the emission line objects, the extinction corrected [OII]
3727 line strengths yield estimates of star formation rates in the
range 0.5-30
yr-1. We used the present data to derive the [OII]
3727 luminosity function up to redshift of 1.2. When combined with [OII]
3727 luminosity densities for the local and high redshift Universe, our results confirm the steep rise in the star formation rate (SFR) to
.
Key words: cosmology: observations - surveys - galaxies: luminosity function, mass function - galaxies: starburst
The Hubble Deep Field-South (HDFS, Williams et al. 2000; Casertano et al. 2000) has been the subject of extensive follow-up ground-based imaging observations at optical (ESO-Imaging Survey (EIS), DaCosta et al. 1998; Big Throughput Camera (BTC) survey, Teplitz et al. 2001), near-infrared (e.g DaCosta et al. 1998; Vanzella et al. 2001; Rudnick et al. 2001; Saracco et al. 2001) and mid-infrared (Infrared Space Observatory (ISO): Oliver et al. 2002; Mann et al. 2002; Aussel et al. 2005, in prep.) wavelengths. The field has also been observed by the Spitzer Telescope (Huang et al. 2005, in prep.). The unprecedented spatial resolution and depth of the Hubble Space Telescope (HST) images have allowed a reliable morphological classification of objects in the field (e.g. the SUNY collaboration, available at http://www.ess.sunysb.edu/astro/hdfs/wfpc2/wfpc2.html). In addition, many groups (e.g. Vanzella et al. 2001; Rudnick et al. 2001) have used various techniques to estimate photometric redshifts for large samples of galaxies detected in the HDFS region. With such a wealth of imaging observations already in place, a systematic spectroscopic follow-up of objects detected in the HDFS is imperative in order to fully exploit the data. To date, only a limited number of redshifts have been determined for selected objects in the HDFS: Lyman break galaxy candidates by Cristiani et al. (2000) and Vanzella et al. (2002); ISO-detected objects by Rigopoulou et al. (2000) and Franceschini et al. (2003); objects in the vicinity of the quasar J2233-606 by Tresse et al. (1999) and Bergeron et al. (1999). Recently, Vanzella et al. (2002) reported on low resolution spectroscopy of a sample of 65 galaxies from the HDFS targetting primarily high-redshift (z > 2) candidates while Sawicki & Mallen-Ornellas (2003) targeted lower redshift objects.
In this paper, we report on the optical spectroscopy of 100 objects located in
the HDFS and flanking fields. We have determined accurate redshifts for 50 objects. The sample selection, observations and data reduction are explained in
Sect. 2. The redshifts and redshift distribution are given in Sect. 3
while star formation rates and extinction estimates are discussed in Sect. 4.
Finally, in Sect. 5 the [OII]
3727 luminosity function and the
global star formation rate densities are presented. Throughout the paper we use
H0 = 70 km s-1 Mpc-3, q0 = 0.5 unless otherwise stated.
The WFPC2 observations of the HDFS resulted in the detection of a
number of objects down to a limiting magnitude
28 mag (Casertano et al. 2000). From this sample, galaxies were selected with
= 23.7 (corresponding to
= 23 mag), a limit imposed by the sensitivity of the FORS1 150I grism. This
criterion resulted in a primary sample of 32 objects all within the central WFPC2 field. The primary sample was supplemented with objects from the Flanking Fields, all of which have identifications in the ESO Deep Imaging Survey (EIS) and the BTC catalog. Additionally we included objects selected from the ISO HDFS detections (Aussel et al. 2005, in prep.). In all, our final sample included 100 objects.
The spectroscopic observations were carried out as part of ESO program ID 65.O-0418(A) from August 24-28, 2000, using the Focal Reducer Spectrograph (FORS1) (Appenzeller 2000), in the Multi-Object-Spectroscopy (MOS) mode, on the ANTU-ESO telescope (formerly UT1), on Paranal, Chile. The data presented here have been retrieved from ESO's VLT science archive.
In the FORS1/MOS mode, 19 slits are placed on a mask covering a
field of view (FOV) of 6.8
6.8 sq arcmin. To cover the 30 primary objects 5 MOS settings were needed. The supplemental objects were observed in the slits that were not occupied by the primary targets in each MOS setting. The slits were 1.2 wide and 22.5 arcsec long. In most cases the target galaxy was placed in the middle of the slit, so as to enable proper background substraction. In a few cases a second galaxy happened to fall in the same slit. Spectra of
these serendipitous galaxies were also analysed. The GRIS-300V, and GRIS-150I grisms were used, covering the 4450-8650 Å and 6000-11 000 Å ranges, respectively. The choice
of grism was based on previous photometric redshift
estimates for sample galaxies; for galaxies with
0.7 the
GRIS-300V grism was used, while the GRIS-150I grism was used for those
galaxies with
0.7. The pixel scale was 0.2 arcsec/pixel while the resolution was
= 440 and 260 for the GRIS-300V and GRIS-150I grisms,
respectively. At each MOS position, three exposures of 1900 s each were obtained, for a total integration time of 5700 s per position. The seeing during acquisition of the spectra varied between 0.8-1.5 arcsec. In Table 1 we present details of the observations.
Table 1: Log of the FORS1/MOS Observations.
For each field there were 3 separate exposures which were co-added before
spectral extraction.
We reduced the MOS spectra by treating each of the 19 slits separately
and using standard IRAF
routines. The
individual spectra were bias-subtracted and flat-fielded in the
standard manner. To remove the background sky emission we fit a
low-order polynomial in the spatial direction at each column and then
subtracted it for each individual exposure. Most of the cosmic ray
hits in the sky region of the spectral images were excluded from the
fits with a sigma-rejection algorithm. We then extracted a
variance-weighted one-dimensional spectrum. The final spectrum was
wavelength calibrated using arc lamp spectra and flux-calibrated using
the spectro-photometric standard stars included in the MOS frame. The
spectra of the standard stars were reduced in the same way as the
galaxy spectra.
In Table 2 we present an inventory of all objects targeted with FORS1/MOS. We list galaxy number, coordinates (in those cases where two galaxies are included in the slit we report only the position of the sample galaxy) and optical (UBVRI) and near-infrared (JHK) photometry. We give AB magnitudes, taken mostly from EIS but supplemented with data from the BTC survey available from http://hires.gsfc.nasa.gov/research/hdfs-btc. Additionally, we list HST identifications from the HDFS-WFPC2 catalogue (available from http://archive.stsci.edu/pub/hdf-south/version2) and the HDFS WFPC2 flanking fields catalogue (available from http://archive.stsci.edu/pub/hdf-south/version1). Finally, we also give ISO identifications (from Aussel et al. 2005, in prep.).
In Figs. 1a-1f and 2a-2b we show the final wavelength- and flux-calibrated spectra for each object in the M 2, M 4, M 5, M 10 (Figs. 1a-1f) and M 9 fields (Figs. 2a-2b). The spectra have
been grouped in the two figures according to the wavelength coverage:
spectra in Fig. 1 cover the 4450 to 8200 Å range (corresponding
to the GRIS-300V grism) while those in Fig. 2 cover the 6000-9200 Å range (corresponding to the GRIS-150I grism). Redshifts were determined based on the wavelengths of the detected emission lines
(usually the [OII]
3727 line originating in a starburst) or, in a few cases, absorption features such as the CaII H, K features (at 3968.5 and 3933.7 Å) prominent in the spectra of elliptical galaxies. When multiple emission lines were detected, we averaged the
redshifts from these lines. In the cases where we detect
absorption features we have not attempted to measure any properties of the
feature. The shape of the underlying continuum (especially in cases where emission
lines were present) was also taken into account and helped to
constrain the redshifts more efficiently.
In each spectrum we report the galaxy number (from Table 2) and
the corresponding redshift and mark the position of both detected and expected
emission/absorption lines. Note that emission features are not always seen at the expected locations of some lines (e.g., H
or [OIII]
5007).
In several spectra there are residual [O I] night-sky lines at
5577 and 6300-6364 that can sometimes lead to "spurious'' emission features in the spectra. These lines are highly variable and could not always be properly removed during the sky subtraction step.
In total we have obtained spectra for 100 objects. Some of the objects were determined to be stars. We have detected continuum in 78 objects and measured redshifts in 50 of them. We supplemented our redshift estimates with additional values taken from the literature (Rigopoulou et al. 2000; and Franceschini et al. 2003).
In Table 3 we give the redshifts and the emission line fluxes measured
from the present spectra. In the majority of the spectra we detect
[OII]
3727 and in some cases the [OIII]
4958
and
5007 lines as well. In half a dozen galaxies we measured
H
4861 line fluxes. Our measured redshifts are in
agreement with the redshifts of Sawicki & Mallen Ornellas (2003), although
the two lists are complementary as these authors target
lower redshift objects. Although we did not perform a detailed classification of
objects based on their spectral properties, it is evident that the majority
of the sample seem to be starburst galaxies in the redshift range
z = 0.6-1.2. We detected only a small fraction (less than 10%) of
elliptical galaxies
.
The highest redshift object in our sample (No. 48) is a QSO at
z = 1.561. Because of its redshift the Balmer lines are redshifted into
the near-infrared regime; in our optical spectra we detected
C III]
1909 and MgII
2798 Å. There is very good
agreement between the redshift determined from the present optical spectra and the near-infrared
spectra of Franceschini et al. (2003). Redshifts for objects Nos. 37 and 89 in our sample
were determined independently from near-infrared data as well, and the
agreement with our estimates is quite good. Based on the measured [OII]
3727 line fluxes we
have calculated the individual expected SFR ([OII]
3727) (see Sect. 4, for details). We note that the spectral profiles were examined carefully in order to exclude objects displaying broad lines (i.e. containing type 1 Seyferts and QSOs) from the further estimates of the
luminosity function (LF) and the total SFR density (see Sects. 5.1 and 5.2).
We estimated the median redshift as a function of the
magnitude using the values given in Table 2. A plot of the median redshift vs.
magnitude is presented in Fig. 3. For the plot we have used the entire FORS sample, i.e. objects from the main WFPC2 area, the Flanking Fields and the ISO selected targets (see Table 2). Although we have determined redshifts for the majority of the objects, the sample is rather incomplete at the fainter magnitudes.
Assuming that ni(z) and ni(tot) are, respectively, the number
of galaxies with spectroscopic redshift and the total number of galaxies (down
to
23.5), then the completeness function in a given magnitude interval
is defined as:
| (1) |
![]() |
Figure 3:
Median redshift as a function of
|
| Open with DEXTER | |
Finally, Fig. 4 shows the redshift distribution for our spectroscopic sample. The distribution shows an evident peak at z
0.58 which is most likely due to the large scale structure present
in the HDF-S region (e.g. Arnouts et al. 2002; Vanzella et al. 2002).
![]() |
Figure 4: Redshift distribution for all objects with confirmed spectroscopic redshifts. |
| Open with DEXTER | |
For ionization-bounded HII regions the Balmer emission line luminosities scale directly with the ionizing luminosities of the exciting stars and thus are proportional to the star formation rate (SFR). Therefore, it is possible to use the Balmer emission lines to derive
SFR in each galaxy. Traditionally, H
has been used to estimate SFR in galaxies in the redshift range up to z
0.4; beyond this redshift, H
moves
out of the optical window and can be accessed at z > 0.8 in the near-infrared regime.
At redshifts 0.5 < z < 1.4, the [OII]
3727 emission line
can be used to study the star formation rates.
The conversion factor between [OII] luminosity and SFR is computed
using an evolutionary synthesis model. To estimate SFR for the current
sample of HDFS galaxies we use the relationship given by Kennicutt (1998)
assuming a Salpeter IMF (0.1-100
), a constant star formation scenario and solar abundances:
![]() |
(2) |
[OII]
3727 luminosities were calculated from the measured
[OII]
3727 fluxes and the measured redshifts. We applied aperture
correction factors (with average value of 1.3) to account for the light missed by our 1.2'' slit. We calculated aperture corrections in two different ways: (1) we smoothed the WFPC2 F814 image to a
resolution of 1.2'' (to simulate the seeing during the present observations) and then measured the ratio between the total object counts (in the image) and the object counts through an
artificial 1.2'' slit; (2) for each galaxy we estimated the continuum flux at 3727 Å from the
magnitude (recall that for z
0.6 the V band filter samples the rest-frame
UV light) and calculated the [OII]
3727 flux from the Equivalent Width (EW, [OII]
3727) measured from our spectra. Both methods lead to similar correction factors of
1.3.
However, because the observed [OII]
3727 line flux is extremely
sensitive to extinction, an estimate of extinction should be made before we calculate SFRs.
We base our extinction estimates on the comparison between predicted and
observed V-K colour index (magnitudes taken from Table 2).
We use the Starburst99 code (Leitherer et al. 1999) for various star formation
histories (i.e., bursts of different durations and ages,
and continuous star formation) and calculate the range of intrinsic
colours. The intrinsic V-K colours predicted by the models
are in the range 1.1-1.5. We apply infrared (Poggianti 1997)
and optical (Coleman 1980) K-corrections for spiral (Sc) galaxies
(assuming the median redshift of our sample z
0.6).
Comparing the observed V-K colours of our galaxies with the predicted
ones we obtain a median extinction
of 0.9 assuming a screen model for the extinction (we assume no extinction in the K-band). This corresponds to an upwards correction factor of 3.5 in
the [OII]
3727 line fluxes. We list the final aperture-corrected and extinction-corrected
individual SFRs in Table 3.
A direct estimate of the present-day SFR activity in the Universe can be obtained by constructing the [OII]
3727 luminosity function (LF) from the present sample galaxies.
We used the the
method (e.g. Felten 1977) to estimate the number density of galaxies in the various luminosity bins. In computing
however, we have to take into account that the sample is incomplete especially towards the fainter magnitude bins. Given the completeness function (Eq. (1)) estimated in Sect. 3.2 the appropriate formula for each galaxy's
volume
is:
![]() |
(3) |
![]() |
(4) |
The LF constructed according to the
method is shown in Fig. 5.
We then fit the LF with a Schechter (1976) function, and the
best-fitting parameters (for a limiting [OII]
3727 flux of
1.02
10-17 erg cm-2 s-1) are:
= -1.24
0.12,
= 10
Mpc-3, L* = 10
erg s-1.
The solid line in Fig. 5 represents the Schechter function for the LF.
Since the luminosity function is well fitted by a Schechter function with
-2,
can be integrated over the whole range of luminosities:
![]() |
(5) |
In addition to the total [OII]
3727 luminosity density estimate
we divided our sample into three redshift bins at
= 0.4 and
= 0.8 and
= 1.2, while for each individual bin we estimated the
[OII]
3727 luminosity density. For each individual bin we used the sum of the individual densities for each galaxy (as described above) to find:
10
erg s-1 Mpc-3 at
= 0.4,
10
erg s-1 Mpc-3 at
= 0.8, and
10
erg s-1 Mpc-3 at
= 1.2.
The [OII] luminosity density as a function of median redshift is shown in Fig. 6.
Although these later estimates suffer from small number statistics (in each bin we have
10 < n < 20 galaxies), we note an evident trend towards higher
luminosity density with increasing redshift over the range 0.4 < z < 1.2.
![]() |
Figure 5:
[OII] luminosity function based on the
|
| Open with DEXTER | |
![]() |
Figure 6: [OII] luminosity density as a function of median redshift. |
| Open with DEXTER | |
Using Eq. (2) we can transform L[OII]
3727 into a global SFR density.
The total [OII]
3727 luminosity density we estimated in 5.1, translates into an SFR density of 0.015
0.03
yr-1 Mpc-3. If we assume a median extinction of
= 0.9 (see Sect. 4) then the total extinction-corrected [OII]
3727 luminosity density translates into SFR density of 0.11
0.02
yr-1 Mpc-3.
In Fig. 7 we plot SFR densities vs. redshifts from the present spectroscopic survey (for the redshift range 0.4-1.2) supplemented with data from other surveys in the
literature. We
restricted ourselves to SFR estimates based on [OII]
3727 surveys in order to avoid conversions from H
to [OII]
3727 luminosities.
For the local Universe we used the SFR estimated by
Gallego et al. (2002), and Hammer et al. (1997) for redshifts between 0.5 and 1.0. For redshifts beyond z
1.5 we used the results from Vanzella et al. (2002) and Steidel et al. (1997). Our results are in good agreement with the Canada-France-Redshift-Survey
(CFRS, Hammer et al. 1997) results for the same redshift range. Unfortunately, most surveys provide only the observed luminosities and not the extinction corrected values. For the present data we plot both observed and extinction-corrected values.
![]() |
Figure 7: SFR densities based on the present data, supplemented with literature values. Both the observed and extinction-corrected values are shown for the present sample and are represented as open and filled squares. Filled diamonds correspond to Gallego et al. (2002); filled triangles are for Hammer et al. (1997), filled pentagons for Vanzella et al. (2002) and filled stars are for Steidel et al. (1997). |
| Open with DEXTER | |
Our present data show evidence for strong evolution from the
local Universe out to z
1.2. Such a strong evolution of the SFR has
also been observed with the H
survey of Tresse et al. (2002)
and the CFRS (Hammer et al. 1997). Assuming that the evolution of the SFR density with redshift follows a relationship of the (1+z)n form, then the combination of all the SFR data points implies a value for
4. We caution, though,
that this result is to some degree dependent on the local SFR density
value used. As an example, the local Universe SFR density value
derived from the H
survey of Singleton et al. (2001) is lower
than the value of Gallego et al. (1995). Irrespective of the local value,
however, our present data show a clear rise in the SFR density by a factor
of 6 in the 0.4 < z < 1.2 regime.
We have presented new spectroscopic redshifts for 50 galaxies in the HDFS
and Flanking Fields. We determined redshifts based on the detections of
several emission lines such as [OII]
3727,
H
and [OIII]
5007 and the absorption features CaII H, K
3968.5 and 3933.7 A.
The redshift range of the present galaxy sample is
0.6 < z < 1.2 with a median redshift of 1.13 (at I
23.5 not
corrected for completeness). A large fraction of the galaxies detected are
starbursts while the remaining
10% are ellipticals.
The individual star formation rates have been estimated for the emission
line objects and the extinction correction rates range between
0.5-30
yr-1. Based on the present data we have estimated the global SFR density of the Universe out to z
1.3 and found it to be 0.11
0.02
yr-1 Mpc-3 using the total extinction corrected [OII]
3727 luminosity density.
Our data show evidence of a strong evolution from the local (z=0) Universe out to z
1.3, following (1+z)4. Finally, the availability of 8 m class telescopes allows us to sample
the evolution of the z
1 Universe precisely and improve on previous
determinations.
Acknowledgements
We acknowledge support of the EU TMR Network "Probing the Origin of the Extragalactic Background Radiation'' (HPRN-CT-2000-00138). We thank Michel Dennefeld and his team for acquiring the data and Niranjan Thatte for help with data analysis. S.B. thanks MPE for the hospitality during his visits.
![]() |
Figure 1:
Spectra of galaxies observed in the HDFS. The ordinate gives the
flux (in erg cm-2 s-1 Å-1). The wavelength coverage extends
between 4200 Å and 8200 Å. The numbers on the top left corner of
the spectra correspond to objects from Table 2 while the second number
denotes the redshift. Note that for objects 1-39 the resolution is
|
| Open with DEXTER | |
![]() |
Figure 1: continued. |
| Open with DEXTER | |
![]() |
Figure 1: continued. |
| Open with DEXTER | |
![]() |
Figure 1: continued. |
| Open with DEXTER | |
![]() |
Figure 1: continued. |
| Open with DEXTER | |
![]() |
Figure 1: continued. |
| Open with DEXTER | |
![]() |
Figure 1: continued. |
| Open with DEXTER | |
![]() |
Figure 2: Same as Fig. 1. The wavelength coverage extends between 6000-9200 Å. |
| Open with DEXTER | |
![]() |
Figure 2: continued. |
| Open with DEXTER | |