near - infrared surveys prompted the discovery of a population of objects with very red optical - infrared colors ( extremely red objects , eros hereafter ; e.g. elston et al .
1988 , mccarthy et al .
1992 , hu & ridgway 1994 ; thompson et al . 1999 ; yan et al . 2000 ; scodeggio & silva 2000 ) . in general ,
objects have been classified as eros when they had redder colors than late type galaxies ( with negligible dust extinction ) at any redshift .
however , depending on the depth of the photometry and on the available filters , different selection criteria have been used to select eros . in this paper
, eros are defined as objects with @xmath55 ( see sect . 5 for more details on this choice ) .
the red colors of eros are consistent with two classes of galaxies : they could be old , passively evolving elliptical galaxies at @xmath10 which are so red because of the large @xmath7correction .
eros may also be strongly dust  reddened star  forming galaxies or agn .
the observational results of the last few years showed that both classes of galaxies are indeed present in the ero population : on one hand , a few objects were spectroscopically confirmed to be @xmath11 ellipticals ( dunlop et al .
1996 , spinrad et al .
1997 , liu et al .
2000 , and marginally , soifer et al .
1999 ) , or to have surface brightness profiles consistent with being dynamically relaxed early type galaxies ( e.g. stiavelli et al .
1999 , benitez et al .
1999 ) . on the other hand
, other eros have been detected in the sub
mm ( cimatti et al . 1998 , dey et al . 1999 , smail et al . 1999 ,
andreani et .
al . 2000 ) , thus providing examples of high  redshift starburst galaxies reddened by strong dust extinction and characterized by high star formation rates .
the relative contribution of the two classes of objects to the whole ero population is still unknown , but there are preliminary indications , based on near - ir and optical spectroscopy and on surface brightness analysis , that ellipticals may represent the largest fraction of this population ( e.g. cimatti et al . 1999 , 2000 ; liu et al . 2000 ; moriondo et al . 2000 ) . a small fraction of low - mass - stars and brown dwarfs among eros is also expected in case of unresolved objects ( e.g. thompson et al . 2000 , cuby et al . 1999 ) .
the importance of studying eros is clear especially for the clues that they could provide on the formation and evolution of elliptical galaxies .
for instance , existing realizations of hierarchical models of galaxy formation predict a significant decline in the comoving density of the ellipticals with @xmath9 , as they should form through merging at @xmath12 ( kauffmann 1996 , baugh et al .
1996a ) , so that a measure of a decline of their comoving density would provide a stringent proof of these models .
conflicting results have been found so far about such issue : some works claim the detection of a deficit of @xmath11 ellipticals ( e.g. kauffmann et al .
1996 , zepf 1997 , franceschini et al .
1998 , barger et al . 1999 ,
menanteau et al . 1999 ) , whereas others find that a constant comoving density of ellipticals up to @xmath13 is consistent with the data ( totani & yoshii 1997 , benitez et al . 1999 , broadhurst & bowens 2000 , schade et al . 1999 ) .
a potentially serious problem in these studies is suspected to be the influence of the field - to - field variations in the density of eros due to the small fields of view usually covered in the near - infrared .
for instance , barger et al . ( 1999 ) found a very low surface density of eros in a 60 arcmin@xmath2 survey to @xmath14 , while on a similar area mccracken et al .
( 1999 ) observed a density three times larger at the same @xmath7 level .
similar uncertainties have been found in the attempts of deriving the fraction of high redshift galaxies among ir selected samples , which is believed to be a stringent test for the formation of the massive galaxies ( broadhurst et al .
1992 , kauffmann & charlot 1998 ) .
fontana et al .
( 1999 ) , using photometric redshifts , found that the fraction of high@xmath9 galaxies in a collection of small deep fields complete to @xmath15 was low and comparable with the predictions of the cold dark matter ( cdm ) models , but not with passive evolution ( ple ) models .
the preliminary results of eisenhardt et al .
( 2000 ) , suggesting a much higher fraction of high@xmath9 galaxies , are instead consistent with both cdm and ple models .
the main aim of our survey was to encompass the difficulties induced by the cosmic variance , obtaining a sample of eros on a large area , at moderately deep @xmath7 levels , in order to minimize and possibly to detect the effects of their clustering , and to compare the observed density with that expected in the case of passive evolution of ellipticals .
our survey is larger by more than a factor of four than the thompson et al .
( 1999 ) survey , and by more than an order of magnitude than all the other previous surveys for eros , at the same limiting magnitudes . with
the large area covered we aimed also to detect a sample , or place limits to the surface density of the very rare class of extreme eros with @xmath16 .    in this paper , the observational results of this survey and their main implications are presented .
a more detailed interpretation of our findings will be presented in a forthcoming paper ( daddi et al .
, in preparation ) .
the paper is organized as follows : we first describe the data reduction and analysis , then we present the counts of field galaxies . in sect .
[ sec : eros ] the sample of eros is described .
[ sec : clustering ] contains the analysis of the clustering of field galaxies and eros .
the main implications of our findings are discussed in sect .
[ sec : implic ] .
km  s@xmath18  mpc@xmath18 throughout the paper .
the @xmath0 observations were made with the eso ntt 3.5 m telescope in la silla , during the nights of 2730 march 1999 , using the sofi camera ( moorwood et al .
1998 ) with a field of view of about 5@xmath195@xmath20 .
sofi is equipped with a hawaii hgcdte 1024x1024 array , with a scale of 0.29@xmath21/pixel .
the @xmath0 filter has @xmath22 and @xmath23 and it is slightly bluer than the standard @xmath7 filter in order to reduce the thermal background .
the center of the observed field is at @xmath24 = 14@xmath2549@xmath2629@xmath27 and @xmath28 = 09@xmath29 ( j2000 ) .
the observed field is one of the fields described in yee et al .
( 2000 ) to which we refer for details about its selection .
the main criteria were not to have any apparent nearby clusters and to be at high galactic latitude .
the images were taken with a pattern of fixed offsets of 144@xmath21 ( about half of the sofi field of view ) over a grid of 9@xmath3013 pointings .
the total area covered by the observations was about 24@xmath3034 arcmin , with a local integration time of 12 minutes in the central deepest region of the field . in the shallower region , the effective integration time is reduced to not less than 6 minutes .
the total amount of time required to cover the whole field was about 5.5 hours .
the data reduction was carried out using the iraf software .
the images were flat - fielded with twilight flats .
the sky background was estimated and subtracted for each frame using a clipped average of 68 adjacent frames ( excluding the central frame itself ) .
the photometric calibration was achieved each night with the observation of 5 - 7 standard stars taken from persson et al .
the zero - points have a scatter of @xmath31 0.015 magnitudes in each night and a night - to - night variation within 0.02 magnitudes .
each frame was scaled to the same photometric level correcting for the different zero - points and airmasses .
accurate spatial offsets were measured for each frame using the area in common with the adjacent frames .
the images were then combined , masking the known bad pixels , in order to obtain the final mosaic .
the cosmic rays were detected and replaced by the local median using the task _ cosmicrays _ of the iraf package _
ccdred_. the effective seeing of the final coadded mosaic ranges from 0.9@xmath21 to 1.1@xmath21 .
the @xmath1-band data were taken in may 1921 1998 with the 4.2 m william herschel telescope on la palma .
the observations were done using the prime focus camera , equipped with a thinned 2048x4096 pixels eev10 chip , with a scale of 0.237@xmath21/pixel .
this gives a field of view of about 8.1@xmath32 .
a standard johnson @xmath1filter was used .
the whole field has been covered by a mosaic of 6 pointings . each pointing consisted of at least 3 exposures of 1200s taken with small offsets .
the total integration times per pointing was therefore 3600s , with the exceptions of two pointings with 4800s and 6000s .
the photometric calibration was achieved with standard stars taken from landolt ( 1992 ) with a scatter in the zeropoints , from the different stars used , below 0.01 magnitudes .
the images were de - biased and then flatfielded , using a master flatfield constructed from the science exposures , scaled to the same zeropoint and then combined .
the seeing of the final @xmath1-band mosaic was between 0.7@xmath21 and 0.8@xmath21 .
the software sextractor ( bertin & arnouts 1996 ) was run on the @xmath0 mosaic with a background weighted threshold in order to take into account the depth variations across the area , as defined by sextractor . among the detected objects ,
all those with s / n@xmath33 in a 2@xmath34 circular aperture ( twice the average seeing fwhm ) were selected and added to the catalog .
a few spurious detections ( e.g. close to image defects ) have been excluded after a visual inspection of the image .
the final catalog includes 4585 objects . in the central deepest region ,
the 5@xmath35 limiting aperture magnitude is @xmath0(2@xmath34)=19.6 , whereas in the remaining area such limit is @xmath0(2@xmath34)19.2 because of the reduced integration time .
isophotal magnitudes were measured with a limiting threshold of about 0.7@xmath36 corresponding in the central area to a surface brightness limit of about @xmath37 mag arcsec@xmath38 . the aperture correction from 2@xmath34 to total magnitudes
was estimated throughout the area by measuring the difference between the isophotal and the 2@xmath21 aperture magnitudes for the stars with @xmath3916 .
a differential correction , in the range of 0.160.30 magnitudes , was measured for different regions of the mosaic with a typical scatter less than 0.03 magnitudes . for the bright objects the isophotal magnitudes were on average consistent with the kron automatic aperture magnitudes .
however , we adopted the isophotal magnitudes because the kron magnitudes are rather unstable at faint flux levels , where the low signal often does not allow to define the correct automatic aperture .
the total @xmath0 magnitudes were then defined as the brightest between the isophotal and the corrected aperture magnitude .
this allowed to safely assign a total magnitude for both the faint and the bright objects .
the typical @xmath0 magnitude where the corrected aperture magnitude begins to be adopted as the total one is @xmath0@xmath3118 in the central deepest region .    the completeness of our catalog has been estimated by adding artificial objects to the @xmath0 mosaic in empty positions , using the iraf package _
artdata_. point - like sources as well as objects with de vaucouleurs and exponential profiles ( convolved with the seeing psf ) were simulated , and sextractor was run with the same detection parameters as for the real data .
the 85% completeness magnitude for the deepest area is @xmath40 for point - like sources .
the completeness decreases to @xmath31  70% for the worst case that we have tested , i.e. for exponential galaxies with 0.7@xmath21 half  light radius . in the shallower area
the corresponding limiting magnitude is @xmath3 .
most of the @xmath41 galaxies are anyway expected to be only barely resolved with the @xmath0 1@xmath21 seeing ( saracco et al .
1997 ) , and this certainly occurs for the distant @xmath10 ellipticals , and thus their completeness limits can be assumed to be similar to those for stars .      in order to recover the @xmath1-band counterparts of the @xmath0-selected objects , a coordinate mapping between the @xmath0 and the @xmath1 images
was derived .
sextractor was then run in assoc mode with a search box of 2@xmath30 fwhm@xmath42 .
the regions around bright stars or defects in the @xmath1 and @xmath0 band images were excluded from this analysis . the final effective area is 701 arcmin@xmath2 at @xmath3 and 447.5 arcmin@xmath2 at @xmath43 .    whenever an object had s / n@xmath443 in the @xmath1-band image 3@xmath35 limits were assigned .
the 3@xmath35 limiting magnitude in a 2@xmath21 diameter aperture is @xmath45 for most of the area , reaching @xmath46 in the deepest pointing . when s / n@xmath473 , 2@xmath21 diameter corrected aperture magnitudes were assigned to each object .
the aperture correction was derived in the same way as for the @xmath0-band , with slightly smaller corrections because of the better @xmath1-band seeing .
the magnitudes were dereddened for galactic extinction . at the galactic coordinates of the center of our field ( @xmath48 , @xmath49 ) , the extinction coefficients from burstein & heiles ( 1982 ) and from schlegel et al .
( 1998 ) are @xmath50 and @xmath51 respectively .
since the two values are derived in different ways , neither of the two can be discarded .
the average was therefore adopted , obtaining a correction of 0.052 magnitudes in @xmath1 and negligible in @xmath0 .
this introduces an uncertainty of @xmath31 0.03 magnitudes in the dereddened @xmath1 magnitudes .
finally , the @xmath8 colors are defined for all the objects as the difference between the @xmath1 and @xmath0 corrected aperture magnitudes .
thanks to the depth of the @xmath1-band data , colors as red as @xmath52 could be measured down to the @xmath0 magnitude limits of our survey . because of the long integrations used in the @xmath1-band the objects with @xmath53 were saturated .
this effect is obviously more important for stars than for galaxies ; moreover , since we are interested in the extremely red galaxy population , the effects due to saturation have no impact on our results .      the star
 galaxy ( s / g ) separation was done by means of the sextractor class_star parameter in both the @xmath1 and @xmath0 band .
this classification was found to be reliable for objects with @xmath5417.5 and @xmath5524 . because of the seeing variations through the area in both bands , a variable class_star threshold was adopted in different subareas .
given the better seeing of the @xmath1-band , the classification was based mostly on that band , switching to the @xmath0 class_star for objects close to the saturation level in @xmath1 .
objects which are fainter than @xmath0 = 17.5 and @xmath1 = 24 have not been classified and have been considered to be galaxies .
this means that a s / g separation can be provided only for objects with colors @xmath565 at @xmath0=19 , @xmath566 at @xmath0=18 , and so on ( see fig .
[ fig : zoom ] and the upper panel of fig .
[ fig : ccdiag ] ) . from the small number of objects classified as stars which are close to the diagonal straight line which indicates our s / g classification limit in the color - magnitude plane
, we can safely conclude that our inability to properly classify very red , faint objects has almost no effect on the total star number counts and , therefore , on the galaxy number counts .
in addition , it is possible that near the faint limit of our survey a very small fraction of very compact , blue galaxies , such as for instance agn or compact narrow emission line galaxies ( e.g. koo & kron 1988 ) , could have been incorrectly classified as stars .
galaxy number counts in the @xmath0 band can provide more advantages in studying galaxy evolution and cosmological geometry than optical counts because they are much less sensitive to the evolution of stellar population and to the dust extinction .
our survey , which covers the magnitude range @xmath57 , represents the widest among the previous deep surveys at levels fainter than @xmath58 .
[ tab : xx ]    cccc & & & + 11.7  12.2 & 701 & - & 4 + 12.2  12.7 & `` & - & 7 + 12.7  13.2 & '' & - & 9 + 13.2  13.7 & `` & - & 17 + 13.7  14.2 & '' & 4 & 21 + 14.2  14.7 & `` & 6 & 38 + 14.7  15.2 & '' & 16 & 37 + 15.2  15.7 & `` & 30 & 62 + 15.7  16.2 & '' & 74 & 73 + 16.2  16.7 & `` & 100 & 88 + 16.7  17.2 & '' & 178 & 128 + 17.2  17.7 & `` & 372 & 127 + 17.7  18.2 & '' & 633 & 144 + 18.2  18.7 & `` & 892 & 185 + & & & + 18.7  18.8 & '' & 200 & 32 + 18.8  19.2@xmath59}$ ] & 447.5 & 628 & 84 +    table [ tab : xx ] summarizes the number of galaxies and stars detected in each @xmath0 bin and fig .
[ fig : stars1 ] shows the corresponding differential counts in 0.5 magnitude bins .
no correction for incompleteness was applied .
galaxies start to dominate over stars at @xmath60 and their surface density is about a factor of 8 higher than the stellar surface density at @xmath61    the slopes of the galaxy number counts were derived over the magnitude range covered by our survey . at bright magnitudes a slope of @xmath62
is found in the range @xmath63 .
we confirm that the @xmath7-band galaxy counts show a flattening at @xmath64 , where the best fit slope changes from @xmath65 to @xmath66 ( see fig .
[ fig : counts_gala ] ) .
the leveling off of the counts below a slope of 0.4 indicates that the differential contribution to the extragalactic background light ( ebl ) in the @xmath7 band peaks at @xmath67 and then starts to decrease at fainter fluxes .
the contribution to the ebl over the magnitude range @xmath68 sampled by our survey is about @xmath69 nw / m@xmath2/sr , which constitutes about @xmath70 of the estimated ebl from discrete sources in the @xmath7-band ( cf .
pozzetti et al .
1998 , madau & pozzetti 2000 ) .
[ fig : counts_gala ] shows the differential galaxy number counts in our survey compared with a compilation of @xmath7-band published surveys .
no attempt was made to correct for different filters .
as shown in the figure , our counts are in very good agreement with the average counts of previous surveys ( hall & green 1998 ) .
[ tab : median ]    ccc & & + 16.5  17.0 & 143 & 3.54 + 17.0  17.5 & 277 & 3.80 + 17.5  18.0 & 522 & 3.92 + 18.0  18.5 & 759 & 4.08 + 18.5  18.8 & 613 & 4.11 + 18.8  19.2 & 628 & 4.04 +
[ tab : eros ]    c|c|ccc|ccc|ccc|ccc & & & & & + & & & & & & & & & & & & & +    @xmath7117.0 & 701 & 2 & 0.006 & 0.003 &  &  &  &  &  &  &  &  &  + @xmath7117.5 & 701 & 15 & 0.025 & 0.02 & 5 & 0.008 & 0.007 &  &  &  &  &  &
+ @xmath7118.0 & 701 & 58 & 0.051 & 0.08 & 19 & 0.017 & 0.027 &  &  &  &  &  &  + @xmath7118.5 & 701 & 158 & 0.084 & 0.23 & 75 & 0.040 & 0.11 & 7 & 0.004 & 0.01 &  &  &  + @xmath7118.8 & 701 & 279 & 0.111 & 0.40 & 150 & 0.060 & 0.21 & 30 & 0.012 & 0.04 & 2 & 0.0008 & 0.003 + @xmath7119.0 & 447.5 & 220 & 0.116 & 0.49 & 133 & 0.070 & 0.30 & 33 & 0.017 & 0.07 & 4 & 0.0021 & 0.009 + @xmath7119.2 & 447.5 & 281 & 0.126 & 0.63 & 173 & 0.079 & 0.39 & 44 & 0.020 & 0.10 & 5 & 0.0023 & 0.011 + @xmath7119.2 & 701 + 447.5 & 279 + 119 & 0.127 & 0.67 & 150 + 87 & 0.076 & 0.40 & 30 + 27 & 0.018 & 0.10 & 2 + 4 & 0.0019 & 0.012 +    in fig .
[ fig : ccdiag ] the @xmath8 vs. @xmath0 color - magnitude diagram is plotted for both stars and galaxies in our sample .
the diagonal straight line in the upper panel indicates our s / g classification limit in this plane .
because objects above this line are not classified , no star can appear in the upper right corner .
however , this figure shows that there are very few stars with a color redder than @xmath72 , even in the region of this color - magnitude diagram where our morphological classification is still reliable .
this suggests that very few stars should be present in the sample of objects for which no morphological classification was possible .
[ fig : histoc ] shows the distribution of colors for galaxies fainter than @xmath73 , from which the median color at different @xmath0 magnitude levels have been calculated ( see table [ tab : median ] ) . in spite of the presence of objects with color upper and lower limits , the use of the median colors ( instead of the mean colors ) allows unbiased estimates in the @xmath0 range we are considering .
the faintest bin with 18.8@xmath7419.2 includes only the objects detected in the deeper region . the median @xmath8 color of the galaxies increases by 0.5 magnitudes from @xmath73 to @xmath75 and
then it remains almost constant up to the limits of our survey ( @xmath0=19.2 ) .
this trend is similar to what is found by saracco et al .
( 1999 ) for the median @xmath76 galaxies color , which also reaches a maximum at @xmath77 1819 and then it becomes bluer , while the median @xmath78 color gets significantly bluer at brighter ( @xmath79 ) magnitudes ( gardner et al .
1993 ) .
our wide - field survey allows us to select a statistically significant and complete sample of eros which can provide stringent constraints on the density of high-@xmath9 ellipticals ( see sect .
1 ) . for this reason
, eros are defined as objects with @xmath80 because this corresponds to select passively evolving ellipticals at @xmath81 ( see fig .
[ fig : color_fig ] ) . table [ tab : eros ] shows the results of the selection for different magnitude and color thresholds . among the others ,
the threshold @xmath82 is used because it corresponds to the selection of @xmath83 elliptical galaxies ( see fig .
[ fig : color_fig ] ) .
we detected 279 eros with @xmath3 from the whole area and 119 eros with @xmath43 in the deeper area , yielding a total sample of 398 objects ( see fig .
[ fig : zoom ] and table [ tab : eros ] ) .
this is by far the largest sample of eros obtained to date
. a small but complete sample of eros with @xmath16 has also been selected , and we estimate for the first time their surface density to be @xmath84 arcmin@xmath38 at @xmath85 .    a comparison of the surface densities of the eros in our sample with those in the thompson et al .
( 1999 ) survey can be directly done after taking into account the different filters used in the two surveys .
for the @xmath86 filter used by thompson et al .
( 1999 ) , @xmath0 @xmath31 @xmath86 - 0.2 ( adopting @xmath87 ) , and therefore their limit at @xmath88 corresponds to @xmath89 which is the shallower limit of our survey , while their @xmath90 color is bluer than our @xmath8 by about 0.1 magnitudes for the redder objects ( thompson , private communication ) . at the level of @xmath3
we find a density of @xmath91 arcmin@xmath38 for eros with @xmath92 , to be compared with the value of 0.039@xmath930.016 that they find ( these errors are poissonian ) .
thus , the two surface densities are in excellent agreement with each other .
we also verified that the average @xmath8 color of all our objects with @xmath94 ( @xmath95 , determined with a kaplan - meier estimator ) , is in good agreement with the average @xmath96 in @xmath97 by thompson et al .
statistical measurements of the clustering of faint galaxies are important for studying the evolution of galaxies and the formation of structures in the universe .
in fact , the amplitude of clustering in 2d space is a useful probe of the underlying 3d structure ( e.g. connolly et al .
@xcite , efstathiou et al .
1991 , magliocchetti & maddox 1999 ) .
the clustering of galaxies on the sky has been studied extensively especially in the optical , but also in the near - infrared ( e.g. roche et al . 1998 and 1999 , postman et al .
1998 , baugh et al .
our survey , as noted before , is the widest at the limits of @xmath85 .
it is therefore interesting to estimate the clustering of our sample of galaxies .
the angular two - point correlation function @xmath98 is defined as the excess probability ( over a poissonian distribution ) of finding galaxies separated by the apparent distance @xmath99 : @xmath100d\omega_{\rm 1}d\omega_{\rm 2 } \label{eq : dp}\ ] ] where n is the mean density per steradian ( groth & peebles , 1977 ) .
several methods for estimating @xmath98 from a set of object positions have been proposed and used , but the most bias  free and suitable for faint galaxies samples resulted to be the landy & szalay technique ( landy & szalay 1993 , see also kerscher et al
this technique ( adopted for the calculations in this paper ) consists in deriving the counts of objects binned in logarithmic distance intervals , for the data - data sample @xmath101 $ ] , the data - random sample @xmath102 $ ] and the random  random sample @xmath103 $ ] .
these counts have to be normalized , i.e. divided for the total number of couples in each of the 3 samples . from them we can estimate @xmath104 as : @xmath105 - 2[dr]+[rr]}{[rr ] } \label{eq : w_b}\ ] ] which is biased to lower values with respect to the real correlation function @xmath98 : @xmath106 where @xmath107 is the `` integral constraint '' ( groth & peebles , 1977 ) : @xmath108 assuming that the angular correlation function @xmath98 can be described by a power law of the form @xmath109 , then , following roche et al .
( 1999 ) , we can extimate the ratio between @xmath107 and the amplitude @xmath110 using the random
random sample : @xmath111    the amplitude of the real two - point correlation function @xmath98 can then be estimated by fitting to the measured @xmath104 the function : @xmath112    the errors can be estimated , following baugh et al .
( 1996b ) , as : @xmath113{(1+w_{\rm b}(\theta))/dd ) } \label{eq : error}\ ] ] where @xmath114 is the non normalized histogram of @xmath101 $ ] .
( [ eq : error ] ) is equivalent to assuming 2@xmath35 poissonian errors for the correlations , and it gives estimates that are comparable to the errors obtained with the bootstrap technique ( ling , frenk & barrow 1986 ) .
this is necessary because it is known that , as the counts in the different bins are not completely independent , assuming the 1@xmath35 poissonian errors would result in an underestimate of the true variance of the global parameters of the angular correlation ( see mo et al .
1992 ) .    in case of the presence of a randomly distributed spurious component among the analyzed sample of objects (
an example of this case would be a residual stellar component among the galaxy sample ) , the resulting amplitudes are apparently reduced by a factor @xmath115 , where @xmath116 is the fraction of the randomly distributed component ( see e.g. roche et al . 1999 ) , and the corresponding correction should be applied .
[ tab : clust_all ]    cccccc & & & & & + 18.0 & 701 & 1131 & 1.3@xmath930.5 & 0.8 & 4.55 + 18.5 & 701 & 1890 & 1.6@xmath930.3 & 0.8 & 4.55 + 18.8 & 701 & 2503 & 1.5@xmath930.2 & 0.8 & 4.55 + 19.2 & 447.5 & 2222 & 1.6@xmath930.2 & 0.8 & 5.16 +    the random samples used in our analysis were obtained using the pseudo  random number generator routine of the c library function _
drand48_. random samples with up to 200000 objects were used .
typically the number of objects in the random samples were a factor of 100200 larger than the number of observed objects .
the random sample was generated with the same geometrical constraints as the data sample , avoiding for instance to place objects in the regions excluded around the brightest stars .      in our analysis
a fixed slope of @xmath117  0.8 was assumed , as this is consistent with the typical slopes measured in both faint and bright surveys ( e.g. baugh et al .
1996b , roche et al .
1996 , maddox et al . 1990 ) , and because it gives us the possibility to directly compare our results with the published ones that are typically obtained adopting such a slope .
the factor @xmath118 was estimated ( with eq .
( [ eq : c ] ) ) for both the whole and the deeper areas , turning out to be 4.55 and 5.16 respectively ( the angles are expressed in degrees , if not differently stated ) . in fig .
[ fig : all_clust ] the observed , bias corrected , two - point correlation functions @xmath98 are shown ; the bins have a constant logarithmic width ( @xmath119 ) , with the bin centers ranging from 3.6@xmath21 to 15@xmath20 .
we clearly detect a positive correlation signal for our sample with an angular dependence broadly consistent with the adopted slope @xmath117  0.8 , even if the measurements show some deviations , in particular for the brightest samples .
a few cluster candidates are present in our survey .
these possible clusters include galaxies with @xmath120 , and are therefore expected to be at @xmath121 .
a detailed analysis of the cluster candidates will be given in a forthcoming paper .
for the purpose of the present work , we tested that the measured clustering amplitudes are stable in case of removal of the galaxies of the most evident cluster from the sample .
however , the presence of such clusters , most of which happen to be in the shallower area , could partly explain the observed deviations from the fitted @xmath122 power laws for the three brightest samples .
the derived clustering amplitudes are presented in table [ tab : clust_all ] .
the amplitude errors are obtained from the fit assuming eq .
( [ eq : error ] ) .
no correction for the stellar contamination was applied . in fig .
[ fig : roche ] the clustering amplitudes of our samples are compared with other published measurements .
the data are shown together with a number of ple models with different clustering evolutions , which are described in detail in roche et al .
( 1998 , 1999 ) .
our measurements are in good agreement with both the models and the previous estimates of roche et al .
( 1999 ) , except for the point with limiting magnitude @xmath123 , which is however the most uncertain of our data points .    as a check
it was verified that the correlations of the stellar sample are consistent with zero , within the measured errors , at all the scales .
this is a confirmation that the stars are homogeneously distributed on the field ( as they should be ) and the seeing variations across the area did not cause a detectable bias in our classification .
[ tab : clust_eros ]    cccccccc & & & & & + & & & & & & & + 18.0 & 701 & 58 & 24@xmath9310 &  &  & 0.8 & 4.55 + 18.25 & 701 & 106 & 25@xmath935 &  &  & 0.8 & 4.55 + 18.5 & 701 & 158 & 22@xmath933 &  &  & 0.8 & 4.55 + 18.8 & 701 & 279 & 14@xmath932 & 150 & 14@xmath933.4 & 0.8 & 4.55 + 19.2 & 447.5 & 281 & 13@xmath931.5 & 173 & 12@xmath932.3 & 0.8 & 5.16 +    the large sample of eros derived from our survey allowed us for the first time to estimate their clustering properties .
even a simple visual inspection of the sky distribution of the objects with @xmath80 ( see fig . [
fig : distrib_sky ] ) shows that the eros have a very inhomogeneous distribution .
the results of the quantitative analysis of the clustering are shown in fig .
[ fig : eros_clust ] , where the observed , bias  corrected angular correlations @xmath98 of the objects with @xmath8@xmath1245 are displayed .
a strong clustering is indeed present at all the scales that could be measured , and its amplitudes ( table [ tab : clust_eros ] ) are about an order of magnitude higher than the ones of the field population at the same @xmath0 limits .
the correlations are very well fitted by a @xmath125 power law .
no attempt was made to correct the amplitudes for the stellar contamination ( see sect . [ calculus ] ) , and we stress that such corrections would increase them . adopting the errors derived from the fits , our detections are significant at more than 7@xmath35 level for the samples with @xmath80 and @xmath0 limits equal to or fainter than @xmath126 .
the amplitudes shown in fig .
5 suggest a possible trend of decreasing strength of the clustering for fainter eros : the @xmath7119.2 eros are less clustered than the ones with @xmath7118.5 and the difference is significant at @xmath127 , based on the derived errors .
the significance of this effect would , however , decrease if the contamination by randomly distributed stars increases towards the limit of our survey .
although difficult to be quantified , this is likely to happen because , differently from the brightest eros , only a small fraction of the eros fainter than @xmath0 = 18.5 could be morphologically classified ( see fig .
[ fig : zoom ] ) .    defining redder thresholds drastically reduces the number of eros and it is not possible to estimate with sufficient accuracy how the amplitudes change for objects with even redder @xmath8 colors .
we could only verify that the sample of eros with @xmath82 has clustering amplitudes consistent with those of the @xmath80 samples ( see table [ tab : clust_eros ] ) . to measure the clustering of the @xmath56 eros , an area at least 10 times larger than ours ( i.e. @xmath312 square degrees ) at @xmath128
would be needed , assuming that their clustering amplitudes are similar to those of the eros with @xmath80 .
finally , it was studied if and how the clustering amplitude changes as a function of @xmath8 for the @xmath7118.8 sample ( see fig .
[ fig : col_var ] ) . a clear increase of @xmath110 with @xmath8 is present for colors @xmath8@xmath1243.5 , while the @xmath53 sample has an amplitude that is consistent with that of the whole sample of field galaxies .
the variation of @xmath110 can be described with a power law in the range of @xmath129 .
previous efforts to disentangle the clustering properties of the red and blue populations in faint @xmath7selected samples probably failed because the ero population was not sufficiently sampled .
for instance , kummel & wagner ( 2000 ) did not find significant differences in the clustering of objects with color bluer or redder than @xmath8=3.49 for their @xmath13017 sample .
this is not surprising since at k@xmath4417 the ero population is almost absent ( see table [ tab : eros ] and fig .
[ fig : ccdiag ] ) .    to check for the stability of these results , possible systematics that could produce a bias in our work were analyzed .
first af all , as the clustering of our @xmath0selected galaxies is in good agreement with the literature data ( fig .
[ fig : roche ] ) , we can exclude the presence of measurable biases coming from the selection of the sample .    regarding the color measuraments ,
since eros are the tail of objects in the @xmath8 color distribution , systematic variations of the photometric zeropoints across the area could have the effect of creating artificial ero overdensities and voids . to exclude this possibility we verified that the blue tail of the @xmath8 distribution is homogeneously spread across our survey , with a very low clustering amplitude . in case of zeropoints variations these should produce the same effect in both the tails of the color distribution .
moreover , to test the reality of the large void of eros clearly seen in the bottom region of our survey ( see fig .
[ fig : distrib_sky ] ) , the @xmath8 color distribution of the galaxies inside and outside this large void were compared by means of a kolmogorov - smirnov test , selecting only the galaxies with @xmath1314 in both regions .
the probability that the two distributions are extracted from the same population is 43% .
thus , the two regions are fully consistent with each other with respect to the color distribution of the blue population .
this , together with the fact that no underdensity is present in this region when only the bluer galaxies are considered , shows that the void of eros should be considered a real feature .
all these tests strongly suggest that the inhomogeneous ero distribution is a real effect .
the strong clustering signal that we find to increase with the @xmath8 threshold and to reach very high values for the eros is potentially capable of giving insight about the nature of these objects .
the main possible source of objects that may contribute to the ero population , as discussed in the introduction , are old passively evolving @xmath10 ellipticals , dust - reddened starburst galaxies and , in case of unresolved objects , low - mass stars or brown dwarfs . our field being at high galactic latitude ( @xmath49 )
, stars are expected to have no clustering , and to be homogeneously distributed , and certainly not to give the strong clustering signal detected . as for the starburst galaxies
, it must be noted that in such galaxies the red colors are mainly driven by the amount of dust extinction and not by the redshift , as in the case of ellipticals ( see fig .
[ fig : color_fig ] ) , and therefore a wide redshift distribution is expected which should dilute their intrinsic clustering .
moreover , it is known that the iras  selected galaxies ( which are typically star  forming galaxies ) have very low intrinsic clustering ( e.g. fisher et al .
we can therefore reasonably conclude that the observed signal is due to the clustering of high redshift ellipticals .
this is also suggested by studies of the local universe which have shown that early - type galaxies are much more clustered than late - type galaxies ( e.g. guzzo et al .
1997 , willmer et al .
1998 ) . in this regard the results plotted in fig .
[ fig : col_var ] could be qualitatively explained by noting that in selecting redder samples the fraction of early type galaxies increases ( the color of local ellipticals is just around @xmath132 ) and by assuming that these galaxies are intrinsically more clustered , while such plot would be difficult to understand if mainly driven by the strongly reddened starburst galaxies .
these considerations strongly suggest that eros are mainly composed by @xmath10 ellipticals , confirming the previous indications that had been found on this issue .    as the elliptical galaxies are the dominant population of galaxy clusters , we investigated the possibility that the detected clustering of eros could be the result of a few massive clusters at @xmath133 present in our field .
for example , in the region inside the circle in fig .
[ fig : distrib_sky ] , a large ero overdensity is found , that one could suspect to be due to a high-@xmath9 cluster of galaxies . however , fig .
[ fig : eros_cm ] shows that there is no clear color - magnitude sequence among the @xmath134 objects inside that region , suggesting that they do not all belong to a single cluster . in case of a cluster , even at high@xmath9 ,
a well defined color - magnitude sequence is in fact generally observed ( e.g. stanford et al .
1998 ) . in the last years
a few examples of massive @xmath1351 clusters of galaxies have been discovered ( e.g. stanford et al .
1997 , rosati et al .
1999 ) , with x ray luminosity of @xmath136 erg s@xmath18 .
a crude estimate of the number of structures of this sort that could be observed in our survey can be derived by calculating the number of high@xmath9 clusters with l@xmath137 erg s@xmath18 expected in the volume we are sampling . from the x ray luminosity function of such structures at @xmath138 ( rosati et al .
2000 ) we estimate that the expected number of massive clusters in our field in the redshift interval @xmath139 is only @xmath140 ( for @xmath141 ) .
moreover , the detection of the ero positive correlation , following a @xmath125 power law on all the scales from 10@xmath21 to 15@xmath20 ( corresponding to @xmath318 mpc at @xmath138 ) suggests that the clustering signal does not come from a few possible clusters detected in our field , but rather from the whole large scale structure traced by the elliptical galaxies .
our results on the clustering of eros have important consequences on the problem of estimating the density of high-@xmath9 ellipticals ( see sect .
1 ) .    the existence of an ero angular correlation with @xmath125 and with a high amplitude implies significant surface density variations around the mean value even for relatively large areas . in the presence of a correlation with amplitude
a , the rms fluctuations of the counts around the mean value @xmath142 is ( see for example roche et al .
1999 ) : @xmath143 the factor @xmath118 is the same as in eq .
( [ eq : c ] ) and , by applying eq .
( [ eq : c ] ) for several areas , it was found that it can be approximated as : @xmath144 if the area is expressed in arcmin@xmath2 and @xmath125 .
the validity of such an approximation has been tested for square regions and for areas not larger than the ones of our survey . with eq .
( [ eq : sigma ] ) and ( [ eq : c2 ] ) the expected variations of the ero number counts can be calculated , once their clustering amplitude is known .    to verify the consistency of this picture , we derived the distribution of the number of eros ( with @xmath8@xmath1245 and @xmath7118.8 , i.e. those in fig .
[ fig : distrib_sky ] ) that can be recovered in our area by sampling it with a field of view of 5@xmath1455@xmath20 , which is the typical field of view of a near - infrared imager such as sofi . in fig .
[ fig : sofi_f ] the observed frequencies of the number of eros recovered in this counts - in - cell analysis is plotted . as the mean expected
number of eros is about 10 , the poissonian fluctuations would be @xmath1463.2 , while fluctuations with @xmath1475.4 are actually observed . applying eq .
( [ eq : sigma ] ) , the measured clustering amplitude @xmath148 implies @xmath149 , in excellent agreement with the measured @xmath35 value .
we also note that the distribution of the numbers of eros in fig .
[ fig : sofi_f ] is not only asymmetric , but also very broad , ranging from @xmath150=0 to @xmath150=30 . in 29% of the cases
the number of eros recovered is @xmath1515 , corresponding to a surface density half of the real one , while only in 19% of the cases the observed number is @xmath15215 .
this shows that , on average , it is more probable to underestimate the real surface density of these objects .
this is a clear property of the sky distribution that we observe , as the voids extend on a large fraction of the surveyed area .
these results show how strong the effects of the field - to - field variations are in the estimate of the sky surface density of eros . in this respect , it should be noted here that all previous estimates of the number density of high@xmath9 ellipticals were based on surveys made with small fields of view , typically ranging from 1 arcmin@xmath2 in the case of the nicmos hdf - s ( benitez et al . 1999 ) to 60 arcmin@xmath2 in the case of barger et al .
( 2000 ) .
the selection of galaxies with colors @xmath134 can be used to search for elliptical galaxies at @xmath153 ( see fig . [
fig : color_fig ] ) , and to study their evolution by comparing their observed surface densities with those expected from ple or hierarchical models of massive galaxy evolution . in this respect ,
very discrepant results have been obtained so far , making the formation of spheroids one of the most controversial problems of galaxy evolution ( see the introduction ) .
our results on the ero clustering clearly show that for such a comparison to be reliable , both a wide field survey ( resulting in a large number of eros ) and a consistent estimate of their surface density fluctuations are necessary before reaching solid conclusions on the evolution of elliptical galaxies .    in this section , with the main aim to show the effect of the increased uncertanties due to the clustering , a preliminary comparison is presented between the sky density of eros observed in our survey ( table [ tab : eros ] ) and the predictions of an extreme ple model similar to that used by zepf ( 1997 ) . in this model , ellipticals formed at @xmath154=5 and their star formation rate ( @xmath155 ) is characterized by an exponentially decaying burst with @xmath156 , with @xmath157 gyr . adopting the markze et al .
( 1994 ) local luminosity function of ellipticals , and the bruzual & charlot ( 1997 ) models with solar metallicity and salpeter imf , the expected surface densities of passively evolving ellipticals with @xmath92 was calculated for different limiting @xmath0 magnitudes .
[ fig : eros_dens ] shows the comparison between the expected and the observed densities of eros with @xmath92 ( such a color threshold should select passively evolving galaxies at @xmath158 ) . for each data point
we show three different error bars , which are actually the region of confidence in the poissonian case ( at 1@xmath35 ) and in the true ( i.e. clustering corrected ) case ( at 1@xmath35 and 2@xmath35 ) .
such confidence regions have been estimated , following the prescriptions of eq .
( [ eq : sigma ] ) , by finding the range of values for the true average counts @xmath142 for which the observed @xmath150 would represent a deviation of the required number of @xmath35 from the real density . in other words such ranges
are defined from the two solutions of the equation : @xmath159 where @xmath24 is the number of @xmath35 considered .
the amplitudes of the angular correlation function used for the @xmath92 eros are those derived for the @xmath80 eros , which is likely to be a conservative assumption as the amplitudes of redder samples should be higher , as suggested by fig .
[ fig : col_var ] .    fig .
[ fig : eros_dens ] shows that the observed eros densities are indeed lower than the predictions of this particular ple model .
however , even for the most deviant point , this ple model can be rejected at only the 2.5@xmath35 and 2.3@xmath35 level for @xmath141 and @xmath160 respectively , if we use the `` correct '' error bars .
note that the different points plotted in that figure are not statistically independent because they are partially obtained with the same objects ( they are cumulative values ) .
it should be recalled here that the observed ero densities plotted in fig .
14 are an overestimate of the true density of the ellipticals because of the contamination by dusty starbursts ( see cimatti et al .
1998 , 1999 ; dey et al . 1999 ; smail et al . 1999 ) and by field low - mass stars .
the fraction of dusty starbursts in complete ero samples is not known yet , as discussed in the introduction , but our results show that they should not be the dominant population . for instance , assuming that the fraction of dusty starbursts and low - mass stars is 20% and 10% of the ero population respectively , this would decrease the observed densities plotted in fig .
[ fig : eros_dens ] accordingly , but it would also increase by a factor of 2 the clustering amplitudes of the high-@xmath9 ellipticals ( see sect .
[ calculus ] ) , and hence the error bars related to those points . as a consequence , the statistical significance of the difference between data and model in this case would be only at the 2@xmath35 level .
it is relevant to mention that the predictions of the ple models depend very strongly on many parameters that have to be adopted _ a priori _ such as @xmath161 , @xmath162 , the local lf of ellipticals ( uncertain by up to a factor of 2 ) , the redshift of formation @xmath154 , the history of star formation , the metallicity , the imf , the spectral synthesis models .
for instance , even just a decrease of @xmath154 , or a small residual star formation at z @xmath31 1.5 ( menanteau et al .
98 , jimenez et al .
99 ) , would decrease the predicted numbers of eros making them more consistent with our data .
we therefore conclude that it seems premature to reject even extreme ple models at a high level of statistical significance on the basis of these data .
a preliminary comparison of our results can be made with some aspects of the hierarchical models of galaxy formation .
first of all , our findings could qualitatively fit into the predictions of such models , where high-@xmath9 ellipticals should be very clustered ( kauffmann et al .
1999 ) because they are expected to be linked to the most massive dark matter haloes which are strongly clustered at high@xmath9 .
the indication ( marginally significant at @xmath163 level ) of a decrease of the clustering amplitude of the eros with the @xmath0 magnitude ( see sect . [
sec : clust_eros ] ) , if mainly due to the mass of the galaxies , could also fit well in this framework because smaller objects should be connected to smaller dark matter haloes which are expected to be less correlated . on the other hand ,
our results seem to conflict with the predictions made by kauffmann & charlot ( 1998 ) on the fraction of @xmath7-selected galaxies with @xmath164 .
in fact , the fraction of galaxies observed to have color @xmath82 ( which corresponds to the selection of @xmath83 ellipticals ) is about 7% of the total in our survey ( see table [ tab : eros ] ) , to be compared with the 23% of @xmath11 galaxies with k@xmath16519 expected in the kauffmann & charlot ( 1998 ) hierarchical model .
this result on the fraction of @xmath10 galaxies in our @xmath7-selected sample broadly agree with the finding of eisenhardt et al .
the main results of this work are :    * we have presented a survey which covers 701 arcmin@xmath2 and is 85% complete to @xmath3 over the whole area and to @xmath4 over 447.5 arcmin@xmath2 ; the r - band limit is @xmath166 at the 3@xmath35 level . *
the observed galaxy counts are derived over the largest area so far published in the range of @xmath167 .
such counts are in excellent agreement with other published data . *
the median @xmath8 color of field galaxies increases by 0.5 mags from @xmath73 to @xmath168 , and it remains constant to @xmath40 . * a sample of 398 eros has been selected .
this sample is the largest published to date and is characterized by an area larger by about four times than previous surveys .
the ero counts and surface densities have been derived for several color thresholds and @xmath0 limiting magnitudes . in particular
, we find @xmath169 ( poissonian ) eros arcmin@xmath38 with @xmath80 and @xmath170 eros arcmin@xmath38 with @xmath92 at @xmath4 . *
the surface density of eros with @xmath16 has been estimated for the first time to be of the order of @xmath84 arcmin@xmath38 at @xmath85 .
* the angular correlation function of field galaxies , fitted with a fixed slope @xmath171 , has an amplitude @xmath172 at @xmath173 , in agreement with previous measurements .
* for the first time , we detected the clustering of eros , with an amplitude @xmath174 for the objects with @xmath80 , in the range @xmath173 ) which is about a factor of ten higher than that of field galaxies .
the ero two point correlations are very well fitted by a @xmath171 power law . *
the clustering amplitude of the galaxies increases with the @xmath8 color threshold following the relation  @xmath175 , for @xmath176 at @xmath89 . * the strong clustering of eros
is shown to be a direct evidence that a large fraction of these objects are indeed high@xmath9 ellipticals .
our result is therefore the first detection of the large scale structure traced by the elliptical galaxies at @xmath138 . *
the ero clustering explains the conflicting results obtained so far on the density of high-@xmath9 ellipticals in terms of strong field - to - field variations affecting the surveys based on small fields of view ( e.g. @xmath177 arcmin ) . *
taking into account the clustering of eros , even the predictions of extreme ple models for the comoving density of high@xmath9 ellipticals can not be rejected at much more than 2@xmath35 significance level .
we would like to thank nathan roche for providing his models in digital form , gustavo bruzual and stephane charlot for their synthetic stellar population models .
we also thank leonardo vanzi for his assistance during the ntt observations and the anonymous referee for useful comments .
lp acknowledges the support of cnaa during the realization of this project .