R. Pelló1 - D. Schaerer2,1 - J. Richard1 - J.-F. Le Borgne1 - J.-P. Kneib3,1
1 - Laboratoire d'Astrophysique, UMR 5572,
Observatoire Midi-Pyrénées,
14 Avenue E. Belin, 31400 Toulouse, France
2 -
Observatoire de Genève,
51, Ch. des Maillettes, 1290 Sauverny, Switzerland
3 -
Caltech Astronomy, MC105-24, Pasadena, CA 91125, USA
Received 20 January 2004 / Accepted 13 February 2004
Abstract
We report the first likely spectroscopic confirmation
of a
galaxy from our ongoing search for distant
galaxies with ISAAC/VLT.
Galaxy candidates at
are selected from
ultra-deep JHKs images in the core of gravitational lensing
clusters for which deep optical imaging is
also available, including HST data.
The object reported here, found behind Abell 1835,
exhibits a faint emission line
detected in the J band, leading to z=10.0 when identified as Ly
,
in excellent agreement with the photometric redshift determination.
Redshifts z < 7 are very unlikely for various reasons we
discuss.
The object is located on the critical lines corresponding to z=9 to 11.
The magnification factor
ranges from 25 to 100.
For this object we estimate
yr-1 and
yr-1, both uncorrected for lensing. The steep UV slope
indicates a young object with negligible dust extinction.
SED fits with young low-metallicity stellar population models
yield (adopting
)
a lensing corrected stellar mass of
,
and luminosities of
,
corresponding to a dark
matter halo of a mass of typically
.
In general our observations show that under excellent conditions
and using strong gravitational lensing direct observations of
galaxies close to the "dark ages'' are feasible with ground-based
8-10 m class telescopes.
Key words: galaxies: high-redshift - galaxies: evolution - galaxies: starburst - galaxies: active - infrared: galaxies
In Sect. 2 we summarize the observational strategy adopted in the present study, and we present the photometric and spectroscopic characteristics of this object. The redshift identification is discussed in Sect. 3. The physical properties of this galaxy and implications of our finding are briefly discussed in Sect. 4.
To search for galaxy candidates at very high redshift (z > 7) we used the
following strategy. We apply the traditional Lyman drop-out
technique (Steidel et al. 1996)
to deep optical images.
At redshift
the main spectral discontinuity is
due to the Gunn-Peterson (GP) trough, i.e. the nearly complete
absorption of the flux shortward of Ly
(1216 Å), due to the large neutral H column density in the intergalactic medium (IGM).
In practice we search for objects which are non-detected in optical (VRI) bands. In
addition we require a detection in J, H, and Ks, a fairly red J-H due the appearance of
the GP trough in J, and a blue
colour corresponding to an intrinsically blue UV restframe spectrum such as expected for UV bright starbursts. These criteria, devised
from simulations based on synthesis models at different metallicities
(Schaerer & Pelló 2002; Schaerer 2003; Pelló et al. 2003) allow us to
select galaxies with redshifts between 7 and 10. Furthermore the detection in at least
two bands longward of Ly
and the combination with the above H-Ks colour
criterion allows us to avoid contamination by cool stars.
Table 1:
VRIJHKs photometry of #1916.
Near-IR magnitudes are measured within
1.5 arcsec apertures on seeing matched images.
Limiting magnitudes in the visible bands correspond to a
detection limit of
within the near-IR reference seeing
aperture, which is equivalent to a 3
detection on 4 ISAAC pixels.
Table entries denote the filter name (Col. 1), effective wavelength (2), total
exposure time (3), the image quality/seeing (4), and the observed
magnitudes in the Vega (5) and AB (6) systems.
Ultradeep JHKs imaging of the central arcmin2 of the lensing cluster Abell 1835 (z=0.253) was obtained at ESO/VLT with the near-IR imager and
spectrograph ISAAC in February 2003.
The total exposure times and the - excellent - seeing conditions are given in Table 1. Indeed, a large fraction of the total exposure time in H and Ks was obtained under excellent seeing conditions, with 74% of the data
below 0.6'' in H and 84% below 0.4'' in Ks.
Photometric data were complemented by deep VRI observations
taken at the Canada-France Hawaii Telescope (see Table 1)
and R band (F702W) WFPC2/HST images
(Programme ID: 8249, PI: Kneib, exposure time 7.5 ks).
Near-IR images were combined and stacked following a
procedure similar to the one used by Labbe et al. (2003),
which is well suited for point sources like the ones expected here.
Some modifications were adopted to improve the background subtraction within a
cluster of galaxies.
Photometry was performed on a 1.5 arcsec aperture
using the SExtractor package (Bertin &
Arnouts 1996). Instead of adopting the SExtractor error bars, which
are based on sky noise only, we have computed more
realistic background noise statistics by randomly sorting 200 blank regions of 1.5 arcsec aperture on the neighbouring sky area
10 arcsec
around the object. The resulting error estimates (adopted in Table 1)
are larger by 0.1 mag (J, H) and 0.2 mag (Ks) than the SExtractor
values.
The above procedure has also been used to derive the
random noise within 0.6'' apertures used to constrain the limiting
magnitudes in the filters where the object is not detected.
We have also checked that all the optical and near-IR colors for
cluster galaxies and, in particular the H-Ks versus J-H color-color diagram, are fully compatible with the expected values for
z=0.25 galaxies.
A more detailed report of these observations and
analysis of the photometric data will be published elsewhere (Richard
et al. in preparation).
![]() |
Figure 1: Composite JHKs ISAAC image showing: a) The core of the lensing cluster A1835, with the position of the long slit used during our spectroscopic survey with ISAAC/VLT, together with the location of objects #1916 and #2582, the reference acquisition star (circles), and a nearby field galaxy seen on the 2D spectra. Large and small dots show the position of the external and internal critical lines at z=10. b) Thumbnail image in the HST-WFPC2 R band, where the object is not detected. c), d), e) J, H, and Ks images respectively. f) Composite JHKs image. All thumbnail images are scaled to the original ISAAC pixel size, without smoothing, and are displayed with a linear scale. |
Open with DEXTER |
![]() |
Figure 2: True color composite images displaying the color difference between the lensed galaxy #1916 and other objects in this field: a) JHK and b) RJH. This object is blue in JHK as compared to cluster galaxies due to its intrinsically blue UV rest frame spectrum in HK, whereas it is red in RJH due to the presence of a break within the J band. |
Open with DEXTER |
Applying the above
search criteria has provided six z > 7 galaxy candidates in the observed
field. Among them, the best candidate (called #1916 hereafter) was
retained for subsequent follow-up spectroscopy. The coordinates of this
object are
:01:00.06,
:52:44.1. The location and
broad-brand spectral energy distribution (SED) of #1916 are shown in
Figs. 1 and 3; photometric data is given in Table 1.
The radius of the circles in the thumbnail images in Fig. 1 is similar to
the aperture used to compute magnitudes.
Limiting magnitudes adopted in Table 1 correspond to a
detection limit of 1
within the near-IR reference seeing
aperture (0.6''), which is roughly equivalent to a detection limit
of 3
within a region equivalent to 4 ISAAC pixels.
The source is not detected on the WPC2/HST F702W image with
R(F702W) > 27.3 mag
at 2
on 4 HST pixels.
Figure 2 shows the JHK and RJH true color
composite images displaying the color difference
between galaxy #1916 and other objects in this field, in particular
cluster galaxies.
![]() |
Figure 3:
Optical to near-IR spectral energy distribution of the
galaxy #1916. Shown are the broadband photometric measurements and
associated 1![]() ![]() ![]() ![]() |
Open with DEXTER |
The broad band SED of #1916 was used to constrain the source redshift
using our public photometric redshift code Hyperz (Bolzonella et al.
2000)
including numerous
empirical and theoretical spectral templates for both galaxies and AGNs,
and keeping internal extinction as a free parameter. The resulting
redshift probability distribution shows a clear maximum at redshift
.
This photometric redshift estimate is mainly corroborated
by the finding of a strong break of
3.1 to 3.7 AB mags between VRI
and H, which has a high significance independently of the definition
used for the limiting magnitudes. The low J-band flux is
attributed to the presence of the GP trough located within this filter
at these redshifts.
![]() |
Figure 4:
2D spectra showing the detected emission line
of #1916, as well the nearby field galaxy and the
[O III] ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
To search for
possible faint emission lines, we have systematically explored the J band, where
Ly
should be located for objects within the
8 < z < 10 redshift interval,
with
priorities set according to the photometric redshift probability distribution. For #1916, the following effective exposure times were obtained in the different adjacent
ISAAC bands:
m (3.6 ks), 1.255
m (6.4 ks),
1.315
m (10.8 ks) and 1.365
m (18.9 ks). Each band has
600 Å width.
Thus the observed spectral interval covers the range from 1.162 to 1.399
m;
in other words the "surveyed'' interval corresponds to redshifts
for Ly
.
With a spectral resolution for the sky lines of R=3100, corresponding to the instrumental 1 arcsec slit width, the fraction of spectral band lost
because of strong sky emission lines, mainly OH emission-lines, is of the order of 30%.
Spectroscopic data were reduced using IRAF procedures and conforming to the ISAAC
Data Reduction Guide 1.5
, using the
same procedure described in Richard et al. (2003).
The deep exposures resulted in the detection of one weak emission
line at the 4-5 level with an integrated flux
of
erg cm-2 s-1 at a wavelength of 1.33745
m (see
Figs. 4 and 5),
which appears on both the 1.315 and the 1.365
m bands.
If identified as Ly
the observed
wavelength translates to a redshift of
,
in excellent agreement with
the photometric redshift determination.
In the 1.315 m domain, we have also detected 3 emission lines on a secondary target
(called #2582)
included in the same slit as #1916 (Richard et al. 2003).
This galaxy, a faint R-band dropout at z=1.68,
shows [O III]
4959, 5007 Å and H
.
The weakest of them, H
,
has a flux
of
erg s-1 cm-2.
This provides an independent verification that emission
features as faint as
erg s-1 cm-2 could have been
detected in our spectra all the
way from 1.25
m to the end of the J band (1.399
m), with some gaps due to sky
features. At wavelengths below 1.25
m, the detection efficiency becomes smaller. For an
unresolved line of
erg s-1 cm-2 a S/N ratio between 2
and 3.5 (from blue to red)
would be expected at
m. If present, an unresolved line of 10-17 erg s-1 cm-2 would have been detected, with a S/N of 3.6 to 5.5, all the way from 1.162 to 1.399
m.
We used the good S/N of the [O III] 5007 line of the galaxy #2582 at z=1.68, which is seen in both the 1.315 and 1.365
m overlapping bands, to compute the offsets required to
stack both data sets. The reference star in the same slit, together with a
standard star, were used to carefully model the telluric absorptions
in these spectra. The correction applied on the relevant wavelength
domain is smooth. Figure 4 displays the 2D spectra obtained on
the two different overlapping and independent bands.
The combined 1D extracted spectrum around the emission-line is
shown in Fig. 5.
![]() |
Figure 5:
a) 1D spectrum of #1916, extracted from the composite 2D spectrum of the 1.315 and 1.365 ![]() ![]() |
Open with DEXTER |
How plausible is the observation of Ly
from such distant an object?
Despite the presence of neutral HI in the IGM at very high z, all galaxies with
have been found through their Ly
emission and several explanations exist for a
partial transmission of Ly
emission through the ambient ISM
and intervening IGM (Haiman 2002; Santos 2003) which are also applicable to higher redshifts.
In #1916 the comparison of UV and Ly
luminosities indicates
such a possible partial transmission (see Sect. 4).
The observed line is hardly resolved and barely
asymmetric (see Fig. 5). Best fits are obtained for rest
frame widths below 60 km s-1 using a simple
Gaussian to fit the red half of the Ly
emission convolved with the instrumental
profile. Line widths beyond 200 km s-1 (similar to those of the sample of
galaxies of Hu et al. 2003) are excluded with a high confidence level,
but intermediate widths of the order of 100 km s-1 cannot be excluded because
of the low S/N ratio in the wings of the line.
In any case, starburst galaxies show in general a wide range of Ly
profiles,
especially objects which are not selected by their line emission, like #1916.
In our case, a narrower Ly
profile could also be related to the rather small
inferred dark matter halo mass (cf. Sect. 4), corresponding to circular
velocities of the order of
km s-1.
Small or negligible line asymmetries as observed here can be explained
within simple models by a variety of combinations of parameters including
the intrinsic emission line width, the size of the surrounding H II region,
and galactic winds (cf. Haiman 2002; Santos 2003).
We therefore conclude that
the detection of Ly
emission from such a high redshift object and its
observed profile are not unexpected.
How secure is the identification of the single line as Ly,
i.e. the redshift z=10.0?
Although other line identifications corresponding to lower redshifts cannot be ruled out
in principle, they appear very unlikely for various reasons.
A rather "exotic'' solution at
,
albeit of lower probability,
can be found with an empirical template of a metal-poor blue compact galaxy
with strong emission lines such as SBS 0335-052.
In this case the observed strong H band flux is dominated by [O III]
4959, 5007, H
,
and H
emission lines, whereas the spectral range covered by Ks is
devoid a strong lines.
However, to decrease sufficiently the rest-frame UV flux between J and VRI an exceptionally large extinction of
is required.
In such a case the observed emission line should correspond to lines between [O II]
3727 and
4000 Å, including e.g. relatively strong lines
of [Ne III]
3868 or [Ne III]+H7
3968
(cf. Izotov & Thuan 1998). Whereas [O II]
3727 can be excluded (cf. above),
other identifications could be possible. However, more than one emission
line could then possibly have been detected in the observed wavelength
range (cf. above). Although quite unusual, such an explanation has
the advantage of being easily testable, as several very strong emission
lines would be expected in the H band.
For the reasons just discussed redshifts
seem very unplausible.
For
other possible identifications of lines sometimes observed in
ehigh excitation objects and Lyman break galaxies include lines at
Å such as C III]
1909 Å, He II
1640 Å, C IV
1550,
and N V
1240 Å.
However, also in this case, more than one emission line could possibly
have been detected in the observed wavelength range (cf. above).
According to the lensing model for Abell 1835 (Smith et al. 2002), within modelling errors the object #1916 is located on top of
the critical lines at
.
Therefore the magnification factor
is very large,
ranging with a large uncertainty from 25 to 100, i.e. 3.5 to 5 mag! As in principle the
likelihood to obtain the largest amplification is small we will adopt the minimum value
of
.
Once corrected for lensing, the
intrinsic AB magnitude of #1916 is
28.5 and 29 in the H and Ks bands respectively
(cf. Table 1). Clearly, without strong gravitational lensing
and excellent seeing conditions near-IR photometry and
spectroscopy of such a faint source would be impossible with current 8-10 m class
telescopes.
Because of its location close to the critical line the object
must be multiply imaged
by the cluster. However, the identification of the counter images is impossible with the
present data. Deeper observations with the Hubble Space Telescope in the near-IR
should allow to confirm our current interpretation.
Redshift 7-10 sources are expected to have virial radius of the
order of
1 kpc, corresponding to
0.1'',
and thus they are unresolved on our present images. HST images with
diffraction-limited resolution should thus allow to confirm the high-znature of this source (e.g. by the detection of strong tangential magnification or
multiple images), to constrain the lensing configuration through
morphological considerations, and to determine the physical scales
involved in star formation activities at such redshifts.
Estimates of the star formation rate (SFR) can be obtained from the Ly luminosity, the UV restframe flux, or SED fits. For a "standard'' Salpeter IMF from
1-100
and constant star formation at equilibrium we find a lower limit of
yr-1 and
yr-1,
both uncorrected for lensing, adopting published conversion factors (Schaerer 2003;
Kennicutt 1998). The lower SFR derived from Ly
is thought to reflect in large part the effects of photon
"destruction'' on the blue side of the Ly
emission line by scattering in the IGM
and losses due to other possible effects (dust, ISM geometry). SED fits with constant
star formation rate (e.g. shown in Fig. 3 for
and 3 Myr age) yield typically
yr-1 (0.8 to 4.8
yr-1 after lensing correction) depending on the age after the onset of star formation.
Lower SFR values would be obtained for "non standard'' IMFs favouring massive stars.
In general the shape of the overall UV SED does not allow one to distinguish age,
metallicity, and extinction as these parameters alter the UV spectrum in a degenerate
way. In the present case the situation is somewhat different, as the UV restframe
spectrum is very blue when compared to other known starbursts (Heckman et al. 1998) and to model
predictions (Leitherer et al. 1999; Schaerer 2003) -
the UV slope measured from the H and Ks-band flux is
.
Comparisons with observed UV spectra of galaxies and model predictions also show
that the observed UV is unlikely to be significantly influenced by line emission (if
present). From this we can quite safely conclude that the extinction must be small in this
z=10.0 galaxy.
SED modeling also allows us to estimate the total stellar mass and luminosity
involved in this star forming event. Scaling young burst models with metallicities
to the observed SED we obtain
(or lensing corrected values of
)
and luminosities of
(
lensing corrected)
for the above Salpeter IMF. Larger masses are obtained for older bursts. The
total current stellar mass estimated for constant star formation over
3 Myr is in good
agreement with the above estimate. Note that, as star formation is likely to continue for
longer time, the total mass of stars to be formed during this starburst event is likely
larger than the "current'' stellar mass.
The above estimates show that #1916 is comparable to or somewhat heavier than
the most massive old Globular Clusters or the most massive super star cluster observed
in nearby starbursts (Mandushev et al. 1991; Ho et al. 1996; Mengel et al. 2002).
Assuming typical values for the baryonic/dark matter content
and for the star formation efficiency (
)
this translates to a total mass of the dark matter halo of
.
Such "massive'' halos correspond to the collapse of
2
fluctuations
at redshift
,
where, if metals are absent or very
deficient, star formation is expected to occur thanks to cooling by atomic hydrogen
(Tegmark et al. 1997; Barkana & Loeb 2001).
"Massive'' starbursts like the one found here should contribute to the cosmic
reionization; they are possibly even the dominant contributors (Ciardi et al. 2003).
In short, the observed and derived properties of #1916 agree well with
expectations of a young protogalaxy experiencing a burst of star formation at z=10.
If we are witnessing the first star formation event, therefore potentially the formation of a
massive primordial (so-called Population III) star cluster, cannot be established or rejected
from the present data. Additional spectroscopic observations searching for
other possible emission line signatures (e.g. He II 1640, C IV
1548+51 Å or other metal lines) are necessary to answer this question.
In fact a measurement of such lines, although presumably intrinsically fainter than Ly
,
could well be feasible since the observed ("transmitted'') Ly
flux
represents only a relatively small fraction of the flux emitted from the source
(cf. above). The expected flux in other lines could therefore become comparable
to the observed emission line flux.
Finally, our observations show that under excellent conditions and using strong gravitational lensing direct observations of galaxies close to the "dark ages'' are possible with ground-based 8-10 m class telescopes. We are looking forward to the exploration of this yet unknown territory from the ground and with the forth-coming James Webb Space Telescope.
Acknowledgements
We are grateful to R. Behrend, B. Ciardi, T. Contini, A. Ferrara, Y. Izotov, M. Lemoine-Busserolle, D. Pfenniger, D. Valls-Gabaud, S. White and various other colleagues for useful comments and discussions. C. Charbonnel, G. Schaerer also contributed indirectly to the finalisation of this work. We thank the ESO Director General for a generous allocation of Director's Discretionary Time for ISAAC spectroscopy (DDT 271.A-5013). Also based on observations collected at the European Southern Observatory, Chile (70.A-0355), the NASA/ESA Hubble Space Telescope operated by the Association of Universities for Research in Astronomy, Inc., and the Canada-France-Hawaii Telescope operated by the National Research Council of Canada, the French Centre National de la Recherche Scientifique (CNRS) and the University of Hawaii. Part of this work was supported by the CNRS and the Swiss National Foundation. J.P.K. acknowledges support from Caltech and the CNRS.