<<

arXiv:astro-ph/9810199v2 26 Apr 1999 aais vlto omlg:osrain gravita- – lensing observations tional Cosmology: – evolution : words: Key and counts number dependence. wavelength faint the their compute discuss the we for Finally, curves analysis. our depletion evo- by con- density favored number strong is includes a lution which as model scenarios. used hierarchical evolution A be galaxy different can disentangle partic- tail to and straint arclets, high of its the distribution ularly estimate redshift We the arclets. properly of that to distribution show redshift order and in number account, statistical sources into background taken of number par- is distribution In arclets ellipticity image. the HST the the ticular, predict in distribution to redshift prescrip- and evolution counts spectrophotometric galaxy a for new with tion images. the model multiple combine mass we identified cluster (1998a), newly al. of presented, B´ezecourt et number is Following halos) a galaxy on clus- cluster the based includes of (that model core mass ter refined A predictions. earlier with inlpi Kebe l 93,addtriearedshift gravita- a determine B2–B3 and the 1993), of al. of et nature (Kneib con- image pair tional spectroscopically multiple faint we of the properties particular, firm statistical In the galaxies. on observa- lensed light These and new galaxies. image shed faint F675W tions of HST/WFPC2 spectroscopy deep 3.6m a ESO 370: Abell lens Abstract. Accepted , Received eecp fteErpa otenOsraoy aSilla, 3.6m La Kea, USA. Observatory, the Mauna Hawaii, at Telescope Southern with Canada-France-Hawaii the European and and Chile the Facility of Coordinating Telescope European scope Telescope to ⋆ requests offprint Send b of B´ezecourtJ. distribution redshift and counts number sources the Modeling I. 370 Abell in galaxies Lensed 1 2 3 2(10. bl 7;1.52 20.;12.07.1) 12.03.3; 11.05.2; 370; Abell (11.03.4 12 are: codes thesaurus later) Your hand by inserted be (will no. manuscript A&A bevtieMidi-Pyr´en´ees,Observatoire d’Astrophysi Laboratoire aty nttt,Psbs80 70A rnne,TeNeth The Groningen, AV 9700 800, Postbus Institute, Kapteyn nttt fAtooy aige od abig B 0HA, CB3 Cambridge Road, Madingley Astronomy, of Institute z ae nosrain ihteNASA/ESA the with observations on Based 0 = . bandfo h aaaciea h pc Tele- Space the at archive data the from obtained 806 epeetnwosrain ftecluster- the of observations new present We ⋆ ± 1 aais lse:idvda:Ael30– 370 Abell individual: cluster: Galaxies: , 2 0 ..Kneib J.P. , . 0 hc si eygo agreement good very in is which 002 .B´ezecourt, [email protected] J. : 1 .Soucail G. , ubeSpace Hubble 1 n ...Ebbels T.M.D. and , u,UR57,1 vneE ei,F340Tuos,France Toulouse, F-31400 Belin, E. Avenue 14 5572, UMR que, te eydsatlne ore ( sources lensed distant very Other h otdsatglx eetdi rvttoa arc gravitational a is of one detected Indeed, galaxy at galaxies. distant distant most probe the to tool are clusters powerful galaxy a amplification, gravitational to Thanks Introduction 1. eidAel21.Ms fte ofimteln redshift lens the objects confirm them lensed of of Most 2218. 19 Abell measured behind have (1998) Ebbels particular, al. and In et method undertaken. arclet been this gravitational have of validate spectroscopy an programs extensive to of biases, be shape order its study the can In accurately redshift arclet. measure a individual can that one demonstrated if distribution, (1996) estimated mass al given 1998, et clus- a al. For Kneib et scale of 1998). Natarajan Schneider large 1996, contribution & al. and Geiger the et smooth (Kneib and halos dis- a galaxy are ter cluster) of shear distributions sum (the weak mass the component by cluster and described that images best shown multiple has of tortions analysis the h ako togseta etrsi h pia domain optical the 1 in with features and galaxies spectral for 1998), strong al. bright- of et surface lack Ebbels low the 1997, their Soucail the (B´ezecourt to & ness due despite dif- a task But is observational redshifts 1991). ficult arclets al. measuring et magnification, Mellier ( cluster radius 1988, arc an al. the an- et giving within the Soucail hence estimate fixes configuration, mass it optical cluster as absolute the importance of of great redshift scales of the gular is Determining arclets redshift. and evolution large arcs and at to history galaxies formation leading of the 1998) on Broadhurst results & important Frye A2390 and 1998a, 1997) al. al. (Pell´o et et (Trager Cl0939+4713 behind tified erlands UK z cuaemdln fcutrptnil ae on based potentials cluster of modeling Accurate 49 ntecutrC15+2(rn ta.1997). al. et (Franx Cl1358+62 cluster the in =4.92 3 < z < 2 . 5. > z ASTROPHYSICS )hv eniden- been have 3) ackground ASTRONOMY 20.3.2018 AND e.g. 2 J. B´ezecourt et al.: Lensed galaxies in prediction, and allow the accuracy of the mass model to age combinations after centering and cosmic ray cleaning. be improved. Similarly, B´ezecourt and Soucail (1997) have The absolute photometry was obtained using magnitude started the spectroscopy of arclets in Abell 2390, which zero-points given in Holtzman et al. (1995). The photom- has been used to constrain the mass distribution in this etry of the field was obtained with the Sextractor package cluster with a great accuracy (Kneib et al. 1998). (Bertin & Arnouts 1996). A criterion of 12 contiguous pix- Using these accurate cluster mass models and a spec- els above the given detection threshold was used to iden- trophotometric description of galaxy evolution (Bruzual tify an object. The 1 σ detection limit is R675W = 24.9 2 & Charlot 1993, Pozzetti et al. 1996), B´ezecourt et mag arcsec− . From the magnitude histogram obtained al. (1998a) have predicted the expected arclet num- from the catalog, we estimate the completeness limit to ber counts and their redshift distribution. This model be R675W = 25.5. presents many improvements with respect to previous We also built a sample of arclets for the purpose of work (e.g. Nemiroff & Dekel 1989, Wu & Hammer 1993, our study of their photometric and statistical properties. Grossman & Saha 1994, Hattori et al. 1997) as it includes To define an arclet we imposed the following selection cri- many observational limits such as magnitude ranges, sur- teria: at least 12 contiguous pixels above 2σ of the local face brightness cut-off or a choice of the optical waveband, sky level, an axis ratio greater than 2, a central surface −2 and this for any mass distribution, regardless of its com- brightness lower than µ675W = 25.5 mag. arcsec and plexity. a magnitude range 20 < R675W < 26. The final sample Abell 2218 is the first cluster where the number contains 81 arclets and their magnitude histogram is dis- counts and redshift distribution of the background arclets cussed in Section 4.2. have been examined in detail (B´ezecourt et al. 1998a, Ebbels et al. 1998). We propose in this paper to further 2.2. Identification of multiple images and arclets extend this study to another well known cluster lens, namely Abell 370. Its mass distribution was first accu- Abell 370 (z=0.37) is a rich optical cluster dominated by rately derived by Kneib et al. (1993) [hereafter K93] who two giant elliptical galaxies identified as #20 and #35, showed that the mass model has to be bimodal in shape in following the numbering of Mellier et al. (1988). A set order to accommodate the gravitational pair B2–B3. This of multiple image candidates and gravitational arcs are was later confirmed by the HST/WFPC1 observation de- identified on the WFPC2/F675W image and are displayed scribed in Smail et al. (1996) and the X-ray map of the in Figure 1. Their photometric and geometrical properties cluster, displaying a bimodal shape compatible with the are also summarized in Table 1. We now discuss them in lens model (Fort & Mellier 1994). detail: We present new observations of the cluster Abell 370 in A0 : Near galaxy #35 is located the spectacular gi- Section 2: a deep HST/WFPC2 image and spectroscopic ant arc (z=0.724) initially detected by Soucail et al. data. Section 3 discusses the new lensing mass model. In (1987). From the high resolution WFPC2 image more section 4, we use an improved version of the code devel- details on its morphology are clearly visible, suggesting oped by B´ezecourt et al. (1998a) to determine the ex- that it is a gravitationally lensed image of a pected counts and redshift distribution of arclets in Abell (Soucail et al. 1998a). 370. Our analysis explores two different scenarios of galaxy A1 to A6 : A first set of arclets (labelled A1 to A6) was evolution to study their differences, and compute the de- detected from ground based images by Fort et al. (1988) pletion curves of background number counts at different and further discussed by Kneib et al. (1994). Most of them wavebands. Discussion, conclusions and further prospects are blue objects, but none have a spectroscopic redshift to constrain galaxy evolution through lenses are presented yet. A5 is the most extended one and presents very blue in Section 5. Throughout the paper, we use a Hubble con- colors and a strong dimming in the reddest bands, sug- −1 −1 stant of H0 = 50kms Mpc , with Ω0= 0 or 1 and gesting a young forming galaxy. Despite deep spectro- ΩΛ = 0. scopic data for A5 (Mellier et al. 1991) no significant emis- sion line has been identified suggesting that 1

Fig. 1. Detailed view of the multiple image candidates detected in the WFPC2/F675W image. B2–B3–B4 is a triple image configuration, as well as C1–C2–C4. R1/R2 is a radial arc. E1–E2–E3 is also a triple configuration with a clear inversion of parity between E2 and E1/E3 (see text for more details). multiple image systems. A redshift estimate based on the mass model was proposed for each of the pairs: zB =0.865, zC =0.81 and zD =0.95. The reality of B4 was confirmed Table 1. Main geometrical and photometric properties of the a posteriori by the HST/WFPC1 data (Smail et al. 1996). multiple images candidates identified on the F675W HST im- C3 is likely to be a wrong identification of the counter age. The origin of the coordinates is the center of galaxy #35 image of C1–C2, and we denote C4 the correct counter and the XY orientation is the CCD one, i.e. North is to the image used in the present model. top with an angle of 85◦ 41′ clockwise from the X axis and The radial arc R : In the HST/WFPC1 image, Smail East is to the left. Coordinates are in arcseconds. e is the el- et al. (1996) discovered a radial arc candidate R (R1-R2) lipticity of the objects (i.e. 1 – b/a where b/a is the axis ratio located close to galaxy #35. They modeled it as a 5-image of the isophotes) and θ is the orientation of the isophote with configuration, and predicted a redshift of z ∼ 1.3 using the respect to the X axis. Surface photometry was computed on ellipse K93 model. This radial arc is well identified in the F675W the two brightest using the package in STSDAS. The two galaxies are well fitted by a de Vaucouleurs image and is clearly the merging of 2 images. −1 −1 law with Re = 32.6h50 kpc for #20 and Re = 44h50 kpc for New multiple images : The detailed insight of the mass #35. model and the exquisite HST resolution allow the identifi- cation of other multiple image candidates. We only discuss ′′ ′′ here the E1–E2–E3 multiple configuration as it presents a Object X ( ) Y( ) e θ R675W characteristic inversion of parity as expected from lensing A1 16.35 57.30 0.74 -23.8 24.83 theory (see also Smail et al. 1995 and Colley et al. 1996 A2 5.56 60.43 0.75 -4.6 24.38 for other spectacular examples). Other multiple image B2 8.46 21.14 0.76 -0.2 23.08 candidates will be discussed in a forthcoming paper. B3 13.80 20.58 0.74 -15.7 23.23 B4 -19.45 20.98 0.42 -41.6 23.82 C1 6.46 6.50 0.17 4.6 24.05 2.3. Spectroscopic Observations C2 10.19 4.96 0.10 63.2 24.00 Spectroscopic data were acquired at the 3.6m telescope of C3 -21.82 2.78 0.14 -60.9 23.99 La Silla (ESO) with EFOSC on October 1996. A long slit R1/R2 5.27 7.03 0.76 76.4 23.98 of 1.5′′ width was positioned along the two objects B2 and E1 32.86 18.62 0.51 74.6 24.83 E2 -31.22 19.15 0.68 64.4 24.53 B3 for a total integration time of 3 hours with an average E3 0.54 22.21 0.66 -69.8 24.44 seeing of 1.9′′. Following the predictions of K93, the [Oii] emission line was expected to lie at λ ∼ 6950A.˚ Thus the # 20 3.09 37.85 0.36 -70±3 16.9 R300 grism was used, providing a useful spectral range # 35 0.00 0.00 0.45 -77±4 17.1 6000A–9000˚ A˚ and a dispersion of 7.5A/pixel.˚ The data were reduced with standard procedures for flat fielding, wavelength and flux calibration in the IRAF environment. 4 J. B´ezecourt et al.: Lensed galaxies in Abell 370

Sky subtraction was performed on the 2D image, and the 2.4. A closer look at the B2–B3–B4 multiple image final spectra were extracted with an optimal extraction algorithm.

Fig. 2. Two dimensional spectrum of objects B2 and B3. The [Oii] emission line is visible for the spectrum of B2 (bottom) and B3 (top) in the center of the image.

Fig. 4. Multicolor photometry of object B3 in U,B,R,I,J and K′. The best fit of the data points corresponds to a star forming galaxy with a constant star formation rate seen at t = 1.7 Gyr, with Z = Z⊙ and a null absorption. Magni- tudes are: U = 22.43 ± 0.1, B = 24.58 ± 0.1, R = 23.84 ± 0.2, J = 22.08 ± 0.2, K′ = 21.26 ± 0.5.

This multiple image is one of the bluest arclets ob- served in this field. Accurate photometry is available in many filters extending from U to K′ with data coming from a U HST image with the F336W filter (ID: 5709, P.I. J.M Deharveng, B´ezecourt et al. 1998b, hereafter Pa- per II), B, R and I images from CFHT (Kneib et al. 1994) and unpublished CFHT infrared images in J and K′ bands (Le Borgne and Picat, private communication). B2 lies too Fig. 3. Spectra of objects B2 and B3 co-added and flux cali- close to galaxy #20 and is contaminated by its envelope brated. The ordinate is Fλ in arbitrary units with zero level at at long wavelengths, hence we only concentrate on B3. the bottom of the graph. The wavelength scale is the observed We compute the spectral energy distribution (SED) one at the bottom, and the rest frame one at the top. The two of B3, corrected from its redshift, in order to study the emission lines detected in the spectra are quoted with their stellar content of the galaxy. We compare and fit the B3 identification, as well as the location of the absorption H and SED with synthetic ones from the Bruzual and Charlot K lines, not clearly detected. spectrophotometric evolutionary code (GISSEL96). This code can provide SEDs for different ranging from Z⊙/50 to 5 Z⊙, different spectroscopic galaxy type Both objects B2 and B3 show a prominent emission determined by the star formation history (burst, ellipti- line at λ = 6728A˚ and λ = 6727A˚ respectively on the two cal, spiral and irregular), and different internal absorp- dimensional spectra, which we identify with [Oii] 3727A˚ tion by dust can be studied. The extinction law used is (Figure 2 and 3). Another feature appears at λ = 9045A˚ the one proposed by Calzetti (1997) which is similar to and λ = 9042A˚ respectively, identified with [Oiii] 5007A˚ what is observed in our galaxy without the 2175A˚ bump. thus confirming the redshifts. We find zB2 = 0.8058 and We explore the 3-dimensional parameter space (spectro- zB3 = 0.8054, clearly demonstrating the similarity of the scopic type, and absorption, E(B − V ) from 0 two objects within the error bars, and in good agreement to 0.3 mag.) to fit the B3 photometry by the model SEDs. with the z =0.865 prediction. This definitely confirms the The main constraint comes from the high UV flux nature of these objects as a multiply imaged galaxy. which requires a weak internal absorption (E(B − V ) < Although, arclets A1 and A2 were observed during 4.5 0.2mag) and a low metallicity (Z

1 L 2 r0 = r0∗ . (3) L∗  The scaling relations adopted are motivated by the prop- erties of the fundamental plane (FP) and are similar to the one used by Brainerd et al. (1996). In particular, the exponent 0.8 in eq. (2) leads to a total mass-to-light ra- tio that scales with L0.3 in agreement with the observed correlation of the FP (e.g. Jorgensen, Franx & Kjaergaard 1996). Brighter galaxies have more extended and therefore more massive dark halos. The orientation and ellipticity of the galaxy halos are taken from the observed values of the light distri- bution while σ0∗ and rt∗ are optimized. r0∗ is fixed at 0.15 kpc. To model the cluster components we considered Fig. 5. Multicolor photometry of object A5 in U,B,R,I, and J two ‘large scale’ mass distributions (represented by 2 trun- (undetected in K′). The best fit of the data points corresponds cated PIEMD mass distribution) centered on galaxies #20 to a burst of star formation redshifted at z = 1.3 and seen at and #35. Their orientations, ellipticities, velocity disper- − t = 0.3 Gyr, with Z = Z⊙/50 and E(B V ) = 0.1 mag. Magni- sions and core radius are left as free parameters. ± ± ± tudes are: U = 21.70 0.15, B = 22.90 0.04, R = 22.26 0.06, The constraints imposed in the optimization procedure ± ± I = 21.67 0.09, J = 20.07 0.11. are the multiple images described in section 2.2, along with the redshifts of A0 and of the B2–B3–B4 triplet. The best fit parameters are summarized in Table 2 (χ2=4.5). These two objects B3 and A5 are typical examples of We also derive a “lensing” redshift estimate for the mul- the population revealed by UV imaging. The UV selec- tiple images. The C system is at zC =0.75 ± 0.1, D is at tion exhibits young objects with low metallicity and low zD =0.85 ± 0.1 and E is at zE =1.3 ± 0.1. R is moved to absorption in the redshift range [0.5,2.0]. larger redshift with zR =1.7 ± 0.2 and A1/A2 could form Such detailed stellar population study, could in prin- a gravitational pair at zA1/A2 = 1.4 ± 0.2. We computed ciple be applied on any other arclets in the field, for the total projected mass in different apertures centered which redshifts have been spectroscopically confirmed or at the barycentre of the #20 and #35 galaxies. We find are strongly constrained by lensing and photometric red- that within 75, 150 and 300 kpc the total mass is, respec- shift techniques (Pell´oet al. 1998a). tively, for this modeling and the previous one (K93): 0.5 6 J. B´ezecourt et al.: Lensed galaxies in Abell 370

Fig. 6. Central part of Abell 370 seen with WFPC2 in F675W. Overlaid is the mass distribution (black thin lines) and the critical lines at z = 0.806 (thick white lines) and z = 4 (thick black lines). The most important arclets and multiple images are shown.

14 14 (0.45) ±0.05 10 M⊙, 1.8 (1.6) ±0.1 10 M⊙ and 4.8 4. Statistical properties of faint galaxies 14 (4.3) ±0.15 10 M⊙. Furthermore, 5% of the total mass is retained in cluster galaxy halos. The total mass-to-light ratio is ∼180 (M/LV )⊙ (out to 300kpc, 160(M/LV )⊙ for 4.1. Modeling the number counts of lensed galaxies the K93 modeling), and ∼9 (M/LV )⊙ for a L∗ galaxy ha- los. These results are similar to the one found in other cluster lenses like Abell 2218, 2390 and AC114 (Kneib et The number counts of gravitational arclets in the field of a al 1996, 1998; Natarajan et al 1998). massive cluster of galaxies are a competition between the magnification of the luminosity by the cluster potential that makes more objects visible and the surface dilution that decreases the surface density of arclets by the same factor as for the magnification. J. B´ezecourt et al.: Lensed galaxies in Abell 370 7

the complex relation given by Kneib et al. (1996) that re- Table 2. Model parameters for A370 potential lates the image deformation τ I as a function of the source deformation τ S and the strength of the potential τ pot: id ε θ σ0 r0 rt a2−b2 −1 −1 ℜ ∗ a2+b2 kms kpc kpc sgn(det A )τ I = τ S + τ pot δS + (τ Sgpot) (6) 35 0.23 -80 1050 75 800  2 20 0.12 -57 1100 91 800 where δ = (1 + τ )=1+ gτ ∗ (see Kneib et al 1996 for ℜ ∗ galaxies – – σ0∗ = 125 r0∗ = 0.15 rt∗ = 15 a full descriptionp of this formalism). The term (τ Sgpot) is a correction for strong lensing only. A lower threshold in the observed axis ratio (a/b > 2 here) corresponds to a lower threshold in the deformation (τI > 0.75) what- The number of arcs brighter than magnitude m with ever the position angle of the image. Within this limit, we an axis ratio greater than q and a surface brigthness min then scan all the allowed solutions for τ S and the source brighter than µ0 within the field of a cluster of galaxies is: position angle θS. The fraction F(xI ,yI ,z) of galaxies at the location (xI ,yI ) in the image plane and at redshift z N(m, qmin,µ0)= which fill the condition is computed using the probability (4) distribution of equation (4) for τ S and an average over zmax ∞ Lmax dV i zl qmin S(q,z) Lmin Φi(L,z) dL dq dz dz the position angle θS. R R R P 2π The sum is over the different morphological types i. z 1 l F(xI ,yI ,z)= p(τS)dτS dθS (7) is the lens redshift and zmax(µ0,i) is the redshift cutoff 2π Z0 Zτ S | τI >0.75 corresponding to the limit in central surface brightness µ0. S(q,z) is the angular area in the source plane (at redshift 4.2. Comparison between the observed and predicted ar- z) that gives arcs with an axis ratio between q and q + dq. clets number counts Φi(L,z) is the luminosity function at redshift z for each morphological type. It was required as a preliminary step that counts in wavebands from U to K and redshift distributions in B and I in empty field are correctly reproduced. We used the model for galaxy evolution of Bruzual and Charlot (1993) with the prescriptions of Pozzetti et al. (1996) for the op- timisation of the parameters representing the field galaxy distributions (see B´ezecourt et al. 1998a for details). The weak point in the B´ezecourt et al. 1998a model, is the approximation of circular sources used for the cal- culation of arclets axis ratio. Indeed, this may underesti- mate the number of galaxies with axis ratio larger than a certain threshold especially at low redshifts and then bias the predicted redshift distribution. To provide a more accurate and reliable redshift distribution, we now take into account the galaxy ellipticity distribution as given by Ebbels (1998) (this distribution is based on the analysis of HST/MDS-like images, and is therefore appropriate to Fig. 7. Number counts of arclets in A370 in the F675W our study): HST image with the following selection criteria: a/b > 2 and µ675W < 25.5. Observed counts (◦) are compared to two num- τ 0.54 ber counts models with respectively q0 = 0.5 (solid line) and p(τ) ∝ τ exp − (5)  0.036  q0 = 0 (dashed line). Error bars are statistical poissonian un-   certainties.

2 2 a −b where τ = 2 a b for an elliptical object with semi ma- jor axis a and semi minor axis b. The complex deformation The magnitude histogram of the arclets observed τ is defined by: τ = τ cos(2θ)+ iτ sin(2θ) where θ is the in A370 is shown in Figure 7. As already stressed in orientation of the major axis. At each point in the image B´ezecourt et al. (1998a), we find that the addition of plane (and at each redshift step for the sources), we es- galaxy scale components and the effect of the source ellip- timate the fraction of the background galaxies which can ticity distribution significantly increase arclets counts by give a lensed image more distorted than a fixed limit. We a factor of 1.4 (q0 = 0.5) and 1.7 (q0 = 0) with respect first consider the effect of the lens deformation through to a K93 model with the assumption of circular sources. 8 J. B´ezecourt et al.: Lensed galaxies in Abell 370

Such increase in the arclets counts was already observed in In order to see how stable the redshift distribution is B´ezecourt et al. (1998a) when using a detailed mass distri- with respect to the ellipticity distribution, another p(τ) is bution but only with circular sources. Considering an ellip- considered following Ebbels (1998): ticity distribution for the galaxies also increases the arclet τ 0.85 counts and strongly modifies the redshift distribution as p(τ) ∝ τ exp − (8) discussed below. With this refinement, the arclet counts  0.20  prediction is now more consistent with the data over the whole magnitude range. Clearly, detailed lens models and representative of objects with 18 2.5) have already been found (Ebbels et al. 1996, Trageretal.1997, Franxetal.1997, Frye & Broadhurst 1998, Pell´oet al. 1998a). – The two galaxy evolution models presented in this pa- per (q0 = 0 and q0 =0.5) have been validated in empty fields where the number counts and redshift distribu- tion of galaxies were compatible with the observed ones. However they present significant differences at high redshift for arclets. The q0 = 0 model with “Pure Luminosity Evolution” (PLE) predicts more high red- shift arclets than the q0 = 0.5 model where number density evolution is included. This different behavior Fig. 8. Redshift distribution of arclets in A370 for the follow- of the redshift distribution at z > 2 is encouraging ing selection criteria: R675W < 23.5, µR < 25.5 and a/b > 2. as it could be a way to distinguish the two scenarios The solid line corresponds to the model q0 = 0.5 and the by analyzing the redshift distribution of arclets with dashed line is for q0 = 0. The dotted line corresponds to the z > 2. Although it is presently difficult to speculate on ellipticity distribution of equation 8 with q0 = 0.5. the true fraction of very high redshift objects in well defined samples of arclets, this population of very high redshift arclets does not seem to dominate the actual The next step is to consider the redshift distribution samples. Hence we favor the q0 = 0.5 model corre- of the arclets. Figure 8 shows the predicted one with sponding to a merging scenario of galaxy formation the selection criteria close to the observational limits for and evolution. faint object spectroscopy (R675 < 23.5, µR < 25.5 and – Because the high redshift domain is sensitive to many a/b > 2). The comparison with real data is not easy at uncertainties included in the evolutionary code (dust present because no well defined sample of arclets, complete obscuration in the UV, uncertainties in the slope of in magnitude, has yet been spectroscopically explored. the IMF for massive or in their UV tracks, influ- The redshift distributions displayed in Figure 8 present ence of the short time scale phenomena), we must be a prominent peak at z ≃ 0.6 and a secondary one at z > 2 cautious on the above conclusions. Furthermore, the due essentially to the contribution of bright and young galaxy evolution models assume that all the galaxies Elliptical-type galaxies in their process of strong star for- form instantaneously at a given redshift, which is a mation. This second peak is strongly attenuated com- clear limitation of these models (we have zform = 4.5 pared to what was discussed in B´ezecourt et al. (1998a) for q0 = 0, which gives for the present day galaxies an where there was clearly a strong excess in the high redshift age of 16 Gyr, and zform = 5 for q0 =0.5, correspond- tail. Taking into account the ellipticity distribution of the ing to an age of 12.2 Gyr). A more phenomenological sources has fixed this problem. model of galaxy formation that follows in detail the J. B´ezecourt et al.: Lensed galaxies in Abell 370 9

mass evolution of galaxy sub units (Baugh et al. 1998) is also very sensitive to the magnitude threshold (see the I may improve this simple description, and would likely curve in Figure 9). UV counts have a slope larger than 0.4 be the best to be compared with the current data. at bright magnitudes (U < 23) but it quickly decreases at fainter levels. This fast flattening is due to the Lyman 4.4. Depletion curves in Abell 370 break that goes through the red limit of the U filter at redshift z ≃ 2.5. Hence, the lack of objects at faint magni- tudes in U produces the depletion curve shown in Figure 9 (see Paper II for more details of the UV modeling). The detection of this magnification bias is highly de- pendant on the poissonian noise of the background sources and the contamination by cluster members and foreground objects. This is critical in cluster cores where the den- sity of objects is very low, of the order of a few units per arcmin2. Hence, these very poor statistics cannot bring valuable information on the sources redshift distribution as most of the information is expected in the dip of the de- pletion curve. Only massive clusters are able to show such a depletion curve but the recovery of the sources redshift distribution seems out of reach with present datasets.

5. Conclusions In this paper, we have presented a new analysis of the cluster lens Abell 370. Thanks to the measurement of Fig. 9. Depletion of lensed object in different filters and lim- iting magnitudes: U < 23 (– · · · –), B < 24.5 (—), R < 24 the redshift for the B2–B3 gravitational pair (z = 0.806) (– – –), I < 23 (– · –) and K < 21 (· · ·) for q0 = 0. The and the identification of several new multiple arclets in small excess near r = 40′′ is due to the local magnification the WFPC2 image, a more accurate and well constrained enhancement of an individual cluster galaxy. model of the mass distribution in Abell 370 is proposed, including galaxy scale components. We studied the spec- tral energy distribution (SED) of arclets B3 and A5 which For a simple description of the background popula- are found to be both young, low metallicity star forming tion, it is easy to show that the number density of objects objects without strong interstellar extinction. brighter than a given magnitude m behind a lens is: The lens model has been used to study the background 2.5α−1 population of galaxies. Taking the galaxy ellipticity dis- N(< m,A)= N0( 2) galaxies, compared to nounced at longer wavelengths (Fort et al. 1996, Taylor et what is currently observed. If this effect is observationally al. 1998). confirmed it may constrain the number density evolution With our model, it is possible to compute radial de- of galaxies in order to interpret the redshift distribution pletion curves instead of global number counts, for any of arclets. filter or magnitude range. A clear illustration of the wave- To understand further the properties of the back- length dependence of the magnification bias is given by ground population through cluster lenses, deep dedicated the ratio of the number of lensed objects expected in a cluster surveys are needed in order to enlarge the num- given area in the field of Abell 370 over the number of ber of arclets and have significant statistical properties. A field galaxies in the same area (Figure 9). One can note multi-color approach is favored as it can help in constrain- that the predicted intensity of the depletion is higher at ing the redshift distribution. This is particularly true for longer wavelengths, as expected, because of the shallower those unresolved sources for which no shape information is slope of the field galaxies counts. Moreover, because of the available and only the spectral energy distribution can dis- flattening of the counts at faint magnitudes, the depletion criminate between background objects and cluster mem- 10 J. B´ezecourt et al.: Lensed galaxies in Abell 370 bers. In addition, detailed study of their spectral energy Kneib J.–P., Pell´oR., Fort B., Mellier Y., Soucail G., Arag´on- distribution will be useful to characterize properties like Salamanca A., Ellis R.S., Smail I.R., Miralda-Escud´eJ., the luminosity function, dust extinction and star forma- 1998, in preparation tion histories of distant galaxies (e.g. Pell´oet al. 1998b). Mellier Y., Soucail G., Fort B., Mathez G., 1988, A&A, 199, 13 Mellier Y., Fort B., Soucail G., Mathez G., Cailloux M., 1991, Acknowledgements. We thank R. Pell´o, R.S. Ellis and Y. ApJ, 380, 334 Mellier for many fruitful discussions and encouragements. Natarajan P., Kneib J.–P., Smail I., Ellis R.S., 1998, ApJ, 499, This research has been conducted under the auspices of 600 a European TMR network programme made possible via Nemiroff R.J. & Dekel A., 1989, ApJ, 344, 51 generous financial support from the European Commission Pell´oR., Kneib J.P., Le Borgne J.F., B´ezecourt J., Ebbels (http://www.ast.cam.ac.uk/IoA/lensnet/). This work was T.M., Tijera I., Bruzual G.A., Miralles J.M., Smail I., Sou- also supported by the Programme National de Cosmologie and cail G., Bridges T.J., 1999, A&A in press, astro-ph/9810390 CNRS. Pell´oR., Kneib J.P., Bruzual G.A., Miralles J.M., 1998b, in “The Birth of Galaxies”, 10th Rencontres de Blois, June 1998, Editions Fronti`eres, astro-ph/9809151 References Pozzetti L., Bruzual G.A., Zamorani G., 1996, MNRAS, 281, AbdelSalam H.M., Saha P., Williams L.L.R, 1998, MNRAS, 953 294, 734 Smail I., Couch W.J., Ellis R.S., Sharples R.M., 1995, ApJ 440, Baugh C.M., Cole S., Frenk C.S., Lacey C.G., 1998, ApJ 498, 501 504 Smail, I., Dressler, A., Kneib, J.-P., Ellis, Richard S., Couch, Bertin E. & Arnouts S., 1996, A&AS, 117, 393 Warrick J., Sharples, R.M., Oemler, A., 1996, ApJ, 469, B´ezecourt J. & Soucail G., 1997, A&A, 317, 661 508. B´ezecourt J., Pell´oR., Soucail G., 1998a, A&A, 330, 399 Soucail G., Fort B., Mellier Y., Picat J.P., 1987, A&A 187, L1 B´ezecourt J., Soucail G., Ellis R.S., Kneib J.P., 1998b, in Soucail G., Mellier Y., Fort B., Mathez G., Cailloux M., 1988, preparation A&A, 191, L19 Brainerd, T., Blandford, R., Smail, I., 1996, ApJ, 466, 623 Soucail G., Kneib J.P., Picat J.P., Ellis R.S., 1998a, in “The Bruzual G.A. & Charlot S., 1993, ApJ, 405, 538 Young Universe – Galaxy Formation and Evolution at Inter- Calzetti D., 1997, Proceedings of the Conference ”The Ultravi- mediate and High Redshift”, Monteporzio, September 1997. olet Universe at Low and High Redshift”, eds. W. H. Waller, A.S.P. Conference Series, Eds S. D’Odorico, A. Fontana, E. M. N. Fanelli, and A. C. Danks (New York:AIP), astro- Giallongo. ph/9706121 Taylor A.N., Dye S., Braodhurst T.J., Ben´ıtez N., Van Kampen Colley W.N., Tyson J.A., Turner E.L., 1996, ApJ 461, L83 E., 1998, ApJ, 501, 539 Ebbels T.M.D., Le Borgne J.F., Pell´oR., Kneib J.–P., Smail Trager S.C., Faber S.M., Dressler A., Oemler A., 1997, ApJ, I.R., Sanahuja B., 1996, MNRAS 281, L75 485, 92 Ebbels T.M.D., 1998, PhD Thesis, U. of Cambridge, UK Wu X.P. & Hammer F., 1993, MNRAS, 262, 187 Ebbels T.M.D., Ellis, R.S., Kneib J.–P., Le Borgne J.F., Pell´o R., Smail I.R., Sanahuja B., 1998, MNRAS 295, 75 Fort B., Prieur J.L., Mathez G., Mellier Y., Soucail G., 1988, A&A, 200, L17 Fort B. & Mellier Y., 1994, A&AR, 5, 239 Fort B., Mellier Y., Dantel–Fort M., 1996, A&A, 321, 353 Franx M., Illingworth G.D., Kelson D.D., van Dokkum P.G., Tran K-V., 1997, ApJ, 486, L75 Frye B., Broadhurst T., 1998, ApJ, 499, L115 Geiger, B. & Schneider, P., 1998, MNRAS 295, 497 Grossman S.A. & Saha P., 1994, ApJ 431, 74 Hattori M., Watanabe K., Yamashita K., 1997, A&A, 319, 764 Hjorth J., Kneib J.–P., 1998, ApJ submitted Holtzman J.A., Burrows C.J., Casertano S., Hester J.J., Trauger J.T., Watson A.M., Worthey G., 1995, PASP 107, 1065 Jorgensen, I., Franx, M. & Kjaergaard, P., 1996, MNRAS, 280, 167. Kassiola A. & Kovner I., 1993, ApJ 417, 450 Kneib J.–P., Mellier Y., Fort B., Mathez G., 1993, A&A, 273, 367 [K93] Kneib J.–P., Mathez G., Fort B., Mellier Y., Soucail G., Lon- garetti P.–Y., 1994, A&A, 286, 701 Kneib J.–P., Ellis R.S., Smail I., Couch W.J., Sharples R.M., This article was processed by the author using Springer-Verlag 1996, ApJ, 471, 643 LATEX A&A style file L-AA version 3.