A&A 442, 109-116 (2005)
S. A. Pustilnik1 - D. Engels2 - V. A. Lipovetsky1, - A. Y. Kniazev1,3 - A. G. Pramskij1 - A. V. Ugryumov1 - J. Masegosa4 - Y. I. Izotov5 - F. Chaffee6 - I. Márquez4 - A. L. Teplyakova7 - U. Hopp8 - N. Brosch9 - H.-J. Hagen2 - J.-M. Martin10
1 - Special Astrophysical Observatory, Nizhnij Arkhyz, Karachai-Circassia, 369167, Russia
2 - Hamburger Sternwarte, Gojenbergsweg 112, 21029 Hamburg, Germany
3 - European Southern Observatory, Karl-Schwarzschild-Strasse 2, 85748 Garching, Germany
4 - Instituto de Astrofisica de Andalucia, CSIC, Aptdo. 3004, 18080, Granada, Spain
5 - Main Astronomical Observatory, 27 Zabolotnoho str., Kyiv, 03680, Ukraine
6 - Keck Observatories, Hawaii, USA
7 - Sternberg Astronomical Institute, Moscow State University, Moscow, Russia
8 - Universitätssternwarte München, Scheiner Str. 1, 81679 München, Germany
9 - Wise Observatory, Tel-Aviv University, Tel-Aviv 69978, Israel
10 - GEPI, Observatoire de Paris, 92195 Meudon Cedex, France
Received 1 December 2004 / Accepted 15 June 2005
We present the sixth list with results of the Hamburg/SAO Survey for Emission-Line Galaxies. The final list resulted from follow-up spectroscopy conducted with the 4.5 m MMT telescope in 1996, and with 2.2 m CAHA and 6 m SAO telescopes in 2000 to 2003. The data of this snap-shot spectroscopy survey confirmed 134 emission-line objects out of 182 observed candidates and allowed their quantitative spectral classification and redshift determination. We classify 73 emission-line objects as definite or probable blue compact or H II galaxies (BCG), 8 as QSOs, 4 as Seyfert 1 and 2 galaxies. 30 low-excitation objects were classified as definite or probable starburst nuclei (SBN), 3 as dwarf amorphous nuclei starburst galaxies (DANS) and 2 as LINERs. Due to the low signal-to-noise ratio we could not classify 14 ELGs (NON). For another 9 galaxies we did not detect any significant emission lines. For 98 emission-line galaxies, the redshifts and/or line intensities are determined for the first time. For the remaining 28 previously-known ELGs we give either improved data the line intensities or some independent measurements. The detection rate of ELGs is 70%. This paper completes the classification of strong-lined ELGs found in the zone of the Hamburg/SAO survey. Together with previously known BCG/H II galaxies in this zone, this sample of 500 objects is the largest to date in a well bound region.
Key words: surveys - galaxies: fundamental parameters - galaxies: distances and redshifts - galaxies: starburst - galaxies: quasars: redshifts
The problem of creating large, homogeneous and deep samples of actively star-forming low-mass galaxies is very important for several applications in studies of galaxy evolution and spatial distribution. Several earlier projects, based on objective prism plates, like the Second Byurakan Survey (SBS) (Markarian et al. 1983; Stepanian 1994), the University of Michigan (UM) survey (e.g., Salzer et al. 1989) and the Case survey (Pesch et al. 1995, Salzer et al. 1995; Ugryumov et al. 1998), as well as some others (e.g., Kitt Peak International Spectral Survey - KISS, Salzer et al. 2000, based on CCD detector registration) identified several thousand emission-line galaxies. The Hamburg/SAO survey (HSS) creates a new very large homogeneous sample of such galaxies in the region of the Northern sky with an area of some 1700 square degrees. It was initiated, in particular, in order to close the gap between the sky regions of the SBS and the original Case survey, and as a result to get the combined sample of low-mass emission-line galaxies in a very large section of sky suitable for the study of their spatial distribution.
Table 1: Journal of observations.
Table 2: Summary of the samples observed and breakdown of the classifications after follow-up spectroscopy.
The basic outline of the HSS and its first results are described in Paper I (Ugryumov et al. 1999), while the additional results from the follow-up spectroscopy are given in Papers II, III, IV and V (Pustilnik et al. 1999; Hopp et al. 2000; Kniazev et al. 2001; Ugryumov et al. 2001). In this, the last paper, we present the results of the follow-up spectroscopy of another 182 objects selected on the Hamburg Quasar Survey (HQS) prism spectral plates as ELG candidates. In Table 2 we show the breakdown of these objects in the samples of the 1st and 2nd priority group, and the categories of detected objects as described below. Out of 134 emission-line objects (galaxies and QSOs) 69 were known as NED objects. For 28 of these galaxies either only redshift, or also some information on emission lines was known, mainly from the previous HSS papers. We included such objects in the presented list since we provide either significantly improved data or some independent measurements.
The article is organized as follows. In Sect. 2 we give the details of the spectroscopic observations and of the data reduction. In Sect. 3 the results of the observations are presented in several tables. In Sect. 4 we briefly discuss the new data and summarize the current state of the Hamburg/SAO survey. Throughout this paper a Hubble constant H0 = 75 km s-1 Mpc-1 is used.
The results presented here were obtained mostly in snap-shot observing mode during one run with the 4.5 m Multiple Mirror Telescope (MMT), two runs with the Calar Alto 2.2 m and seven runs with the SAO 6 m (BTA) telescopes (see Table 1).
The observations were carried out on May 20, 1996, with the Red Channel of the MMT Spectrograph through the long slit of . The 300 grooves mm-1 grating in first order provides a dispersion of 3.2 Å pixel-1, and a spectral resolution FWHM of about 10 Å. To avoid second-order contamination, a L-38 blocking filter was used. The total spectral range was 3700-7400 Å. The spectra were rebinned by a factor of 2 along the spatial axis. Hence, the spatial sampling was 0 6 pixel-1.
Short exposures (3-5 min) were taken in order to detect strong emission lines to allow redshift measurements and a crude classification. The slit was not oriented along the parallactic angle because of the snap-shot observing mode. Reference spectra of an Ar-Ne-He lamp were recorded to provide wavelength calibration. Spectrophotometric standard stars from Oke (1990) and Bohlin (1996) were observed at the beginning and at the end of the night for flux calibration. The dome flats, bias, dark and twilight sky frames were accumulated each night. The weather conditions were photometric, with seeing variations around 1 0 (FWHM).
Follow-up spectroscopy with the CAHA 2.2 m telescope was carried out during two runs (June-July 2000 and February 2002), using the Calar Alto Faint Object Spectrograph (CAFOS). During these runs a long slit of and a G-200 grism (187 Å mm-1, first order) were used. The spatial scale along the slit was pixel-1. A SITE 15 2K2K CCD was operated without binning. The wavelength coverage was 3700- 9500 Å with maximum sensitivity at 6000 Å. The spectral resolution was 12-16 Å (FWHM). The slit orientation was not aligned with the parallactic angle because of the snap-shot observing mode. The exposure times varied within 5-20 min depending on the object brightness and weather conditions. The observations were complemented by standard star flux measurements (Oke 1990; Bohlin 1996), reference spectra (Hg-Cd lamp) for wavelength calibration, dome flat, bias and dark frames. In the run of June-July 2000 the weather conditions were photometric most of the time with a seeing 1.5 (FWHM). During one night of this run, as well as during two nights in February 2002, the weather conditions were variable with a seeing of 3 -4 . The measurements in these nights are marked by "*'' in Table 4.
There was no order separation filter applied, therefore some second order contamination by the object UV light might be present at wavelengths longer than 7200 Å. However, as can be directly seen from the presented spectra, this effect is probably small, since it is undetectable in the continuum behavior around 7200 Å. In principle, one could expect an increase of the line fluxes at wavelengths longer than this due to the second order contamination in the spectra of the flux calibrating stars. There are seven objects whose emission line ratios could be potentially affected. These objects are listed and commented at the end of Sect. 3.1.
The observations with the 6 m telescope (BTA) of the Special Astrophysical Observatory of Russian Academy of Sciences (SAO RAS) were performed mainly as a back-up program. Therefore the weather conditions in most cases were rather poor. The seeing in the majority of the nights was in the range of 2 to 4 (FWHM) and/or the transparency was variable. Results obtained under non-photometric conditions are marked by "*'' in Table 4. In all cases we used the long slit spectrograph (LSS) in the BTA prime focus (Afanasiev et al. 1995) with a Photometrics 1K1K CCD detector with 24 m pixel size. The long slit of 120 was used with the slit width of either 1 5 or 2 0, depending on the seeing and grating. Three set-ups with the gratings of 325, 400 and 651 grooves mm-1 were used during various runs. The wavelength ranges of the spectra covered for different set-ups and their samplings in Å pixel-1 are given in Table 1. The respective effective resolutions were 14 Å, 11 Å and 7 Å.
Reference spectra of an Ar-Ne-He lamp were recorded before or after each observation to provide wavelength calibration. Spectrophotometric standard stars from Bohlin (1996) were observed for flux calibration. All observations were conducted mainly with the software package NICE in MIDAS, described by Kniazev & Shergin (1995).
The reduction of all data was performed at SAO using the standard reduction systems MIDAS and IRAF.
The MIDAS command FILTER/COSMIC was found to be a quite successful way to automatically remove all cosmic ray hits from the images. After that we applied the IRAF package CCDRED for bad pixel removal, trimming, bias-dark subtraction, slit profile and flat-field corrections.
To do accurate wavelength calibration, correction for distortion and tilt for each frame, sky subtraction and correction for atmospheric extinction, the IRAF package LONGSLIT was used.
To obtain an instrumental response function from observed spectrophotometric flux standards, we first extracted the apertures of standard stars. Then the determined sensitivity curve was applied to perform flux calibration for all object images. Finally we extracted one-dimensional spectra from the flux calibrated images. When more than one exposure was taken with the same setup for a given object, the extracted spectra were co-added and a mean vector was calculated. When several observations with different setups (telescopes or grisms) for the same object were obtained, the data were reduced and measured independently and the more accurate values were taken.
To speed-up and facilitate the line measurements we employed the dedicated command files created at SAO using the FIT context and MIDAS command language. The procedures for the measurements of line parameters and redshifts applied were also described in detail in Papers III, IV and in Kniazev et al. (2004).
In Table 2 we present the summary of the observation results. 182 candidates were selected from our first and second priority samples introduced in Paper IV.
Out of 76 first priority candidates (objects showing a clear density peak near 5000 Å and a blue continuum on the HQS prism spectral plates), 36 objects appeared in our list as new ones. 40 objects were listed in the NED as galaxies or objects from various catalogs and 4 of them already had information on emission lines and redshifts in earlier publications. Apart from these 4 objects, 24 more of the mentioned 40 NED galaxies have appeared in our previous HSS papers, but had data of rather low quality. All such objects were included in our observing program in order to improve spectral information. The comparison of our measured velocities with those of galaxies with already known redshift shows acceptable consistency for most objects in common within the uncertainties given. However, for five galaxies originally appearing in the HSS List II, the difference found is as high as 200-300 km s-1, which probably indicates the lower accuracy of some radial velocities from that list.
The remaining 106 observed objects were taken from the list of the second priority candidates, those with less prominent emission features on the high resolution spectra (HRS) obtained after scanning the original HQS objective prism plates. As described in Paper IV, we created from this list the "APM selected sample'', which uses additional information for the selection. The "APM selected'' sample comprises second priority candidates which are classified as non-stellar (at least in one of two filters) on Palomar Sky Survey plates (PSS) in the APM database, and have a blue colour according to the APM colour system ( ). Here we give the spectral data for 64 of them, that looked like ELGs or QSOs. 31 more 2nd priority candidates were classified as stars or galaxies without emission, and 12 objects with no emission lines were not classified at all due to poor S/N ratio spectra.
The parameters of the 126 observed emission line galaxies are listed in
Table 3, containing the following information:
Column 1: number in the table.
Column 2: the object's IAU-type name with the prefix HS.
Column 3: right ascension for equinox B1950.
Column 4: declination for equinox B1950. The coordinates were measured on direct plates of the HQS and are accurate to 2 (Hagen et al. 1995).
Column 5: heliocentric velocity and its rms uncertainty in km s-1.
Column 6: apparent B-magnitude obtained by calibration of the digitized photoplates with photometric standard stars (Engels et al. 1994), having an rms accuracy of for objects fainter than (Popescu et al. 1996). Since the algorithm to calibrate the objective prism spectra is optimized for point sources, the brightnesses of extended galaxies are underestimated. The resulting systematic uncertainties are expected to be as large as 2 mag (Popescu et al. 1996). For about 30% of our objects, B-magnitudes are unavailable at the moment. We present for them blue magnitudes obtained from the APM database. They are marked by a "*'' before the value in the corresponding column. According to our estimate they are systematically brighter by than the B-magnitudes obtained by calibration of the digitized photoplates (rms ).
Column 7: absolute B-magnitude, calculated from the apparent B-magnitude and the heliocentric velocity. No correction for galactic extinction is made because all objects are located at high galactic latitudes and the corrections are significantly smaller than the uncertainties in the magnitudes.
Column 8: preliminary spectral classification type according to the spectral data presented in this article. BCG means a galaxy possessing a characteristic H II-region spectrum with low enough luminosity ( ). SBN and DANS are galaxies of lower excitation with a corresponding position in the line ratio diagnostic diagrams, as discussed in Paper I. SBN are the brighter fraction of this type. Here we follow the notation of Salzer et al. (1989). The non-confident classification is followed by "?''. Three objects (HS 0807+4103, HS 1525+4344, HS 1627+3625) were recognized as Sy 1 galaxies due to the presence of broad Balmer lines and broad [Fe II] emission. HS 1644+3934 was recognized as a Seyfert 2 galaxy. The typical spectrum of low-ionization nuclear emission-line regions (LINERs) is identified for 2 galaxies. 14 ELGs are difficult to classify, mainly due to low S/N. They are coded as NON.
Column 9: one or more alternative names, according to the information from NED. References are given to the other sources of the redshift-spectral information indicating that a galaxy is an ELG.
The spectra of all emission-line galaxies are shown in Appendix A, which is available only in the electronic version of the journal.
The results of line flux measurements are given in Table 4
which contains the following information:
Column 1: number in the table.
Column 2: the object's IAU-type name with the prefix HS. Asterisks refer to the objects observed during non-photometric conditions.
Column 3: designation of the telescope with which the spectral data were obtained. "B'' means BTA, "C'' - Calar Alto 2.2 m telescope, and "M'' - MMT.
Column 4: observed flux (in 10-16 erg s-1 cm-2) of the H line. The accuracy of this and other parameters varies substantially over the whole table. We divided the relative errors into four intervals: 5%, (5-10)%, (10-20)% and (20-50)%. They are marked by the respective superscripts a, b, c and d right of each table entree. For about 40% of ELGs the Balmer absorptions from the underlying stellar population can somewhat affect the measured H emission flux and the related flux ratios. These objects are marked with "''. For several objects with non-detected H emission line, the fluxes are given for H and marked by a "''.
Columns 5-7: the observed flux ratios [O II]/H, [O III]/H and H/H.
Columns 8, 9: the observed flux ratios [N II] 6583 Å/H, and ([S II] 6716 Å + 6731 Å)/H.
Columns 10-12: equivalent widths of the lines [O II] 3727 Å, H and [O III] 5007 Å.
Below we give comments on some specific cases:
HS 1010+4907 and HS 1009+4906 comprise a compact group
(50 kpc in extent) with a fainter galaxy without evident emission
lines, namely HS 1010+4906 (see Table 6).
HS 1353+4706 was classified as an M-star in Paper I (Ugryumov et al. 1999). However, it was suspected that a wrong object had been observed about 0 3 away. This object will be referred to as HS 1353+4706A. The new observations indeed revealed that HS 1353+4706B is a very strong-lined BCG with very low metallicity (12+ (O/H ; see Pustilnik et al. 2004b). The M-dwarf HS 1353+4706A has B1950 coordinates 13 53 25.2 +47 06 46 and its brightness is B>18.7.
Seven galaxies observed with the Calar Alto 2.2 m telescope have H, [N II] or [S II] lines at Å. Their fluxes can be affected by the second order contamination as pointed out in Sect. 2.3. For these galaxies (1231+4349, 1235+4108, 1426+3658, 1437+3724, 1439+3704, 1525+4344 and 1614+4450) the affected parameters in Table 4 can be either the ratio F(H)/F(H), or the line flux F(H), if H was not detected. For the ratios of F([N II])/F(H) and F([S II])/F(H) the effect should be minor since these lines are close in wavelength.
|Figure 1: Distributions of the BCG sample in the zone of HSS of a) total B-magnitude, b) absolute B-magnitude, c) radial velocity and d) equivalent width of [O III 5007] line.|
|Open with DEXTER|
The main criteria applied to search for BCGs are blue continuum near 4000 Å and a strong emission line, the expected doublet [O III] 4959, 5007 Å, in the wavelength region between 5000 Å and the sensitivity break of the Kodak IIIa-J photoemulsion near 5400 Å (see Paper I). For this reason faint QSOs with Ly 1216 Å redshifted to , or with C IV 1549 Å redshifted to , or with Mg II 2798 Å redshifted to could be selected as BCG candidates. In Papers I-V we reported the discovery of a number of such faint QSOs. They were missed by the Hamburg Quasar Survey since it is restricted to bright QSOs ( ). Here we report the discovery of eight faint ( ) QSOs. For four of them we identified Ly 1216 Å redshifted to as the line responsible for its selection. Two objects (HS 1608+3546 and HS 1714+4202) show a broad emission line tentatively identified as Mg II 2798 Å at . Two more quasars with were selected due to the line C IV 1549 Å. Since for HS 1203+3811 only one broad line is seen in a rather poor S/N ratio spectrum, its identification as Ly should be considered as tentative. The data for all eight quasars are presented in Table 5. Finding charts and plots of their spectra can be found on the www-site of the Hamburg Quasar Survey (http://www.hs.uni-hamburg.de/hqs.html).
In total, for 49 candidates no (trustworthy) emission lines were detected. We divided them into three categories.
For nine non-emission line objects the signal-to-noise ratio of our spectra was sufficient to detect absorption lines, allowing the determination of redshifts. The data are presented in Table 6.
To separate the stars among the objects with no detectable emission lines, we cross-correlated a list of the most common stellar features with the observed spectra. In total, 26 objects with definite stellar spectra and redshifts close to zero were identified. All of them were crudely classified in categories from definite A-stars to G-stars, with most of them intermediate between A and F. The data for these stars are presented in Table 7.
It was not possible to classify 13 objects without emission lines. Their spectra have too low signal-to-noise ratio to detect trustworthy absorption features, or the EWs of their emission lines are too small.
As a result we have 182 observed candidates preselected on HQS objective prism plates, out of which 76 were first priority candidates and 106 were second priority. 134 objects (73% of the total) are found to be either ELGs (126) or quasars (8). 24 of these ELGs were presented in the previous HSS papers, and were reobserved in order to improve the data quality.
Seventy two out of 126 ELGs (57%) were classified based on the character of their spectra and their luminosity as H II/BCGs or probable BCGs.
14 ELGs are difficult to classify due to their poor signal-to-noise spectra. Six more ELGs were classified as Active Galactic Nuclei (AGN): 4 as Seyfert galaxies and 2 as LINERs. The remaining 33 ELGs are objects with low excitation: either starburst nuclei galaxies (SBN and probable SBN) or their lower mass analogs - dwarf amorphous nuclear starburst galaxies (DANS or probable DANS).
Summarizing the results of the Hamburg/SAO survey presented in Papers I through VI, we discovered altogether, from the 1st priority candidates, 463 new emission-line objects (26 of them are QSOs). For 100 more ELGs known from the literature (NED) we obtained quantitative data for their emission lines. The total number of confident or probable blue compact/H II-galaxies is 387. Relative to all observed 537 ELGs the fraction of BCGs is 72%.
42 more new BCGs and 56 other ELGs are found among the second priority candidates. Along with the BCGs selected from the HSS candidates, but not observed by us since they already were known from other surveys in this region, the total number of BCGs in the sky region covered by the HSS (1700 sq degrees of a single piece of sky) reaches 500. This constitutes the largest and deepest BCG sample in both hemispheres and will be presented elsewhere as a separate publication. The assembly and verifying of the whole HSS database is underway, and the most up-dated version of this will appear at: http://precise.sao.ru.
In Figures 1a-d we show the distributions of radial velocities, apparent and absolute magnitudes and EWs([O III] 5007) for all found BCGs or BCGs? in the zone of HSS (506 objects, including also the galaxies, selected by us as the HSS candidates, but not observed in this project due to lack of observing time; these objects were already known/classified and had redshift data). The latter group comprises about 70 galaxies which came mainly from the papers by Peimbert & Torres-Peimbert (1992), Vogel et al. (1993), Popescu et al. (1996,1997,1998), Ugryumov et al. (1998), Popescu & Hopp (2000) and some unpublished data on BCGs in the SBS zone, partly intersecting the HSS zone.
Table 8: Characteristics of the HSS BCG sample.
Since for most of the sample galaxies the total B-band magnitudes are not
available, we used their photographic blue magnitudes from the Palomar Sky
Survey, as provided in the APM database (Automatic Plate-measuring
at Cambridge, Irwin 1998), and calibrated them through the sample
of about a hundred BCGs with CCD-measured total B magnitudes.
We obtained the linear regression of:
The distributions shown in Fig. 1 have the characteristics presented in Table 8. A rough estimate of the completeness level from the shown distribution is . The radial velocity distribution peaks near V=7500 km s-1, while the distribution of EW(5007) peaks at the value of 60 Å.
Since the new BCG sample is the largest sample of low-mass galaxies and is well situated on the sky, it can be used for several purposes. First, as it was assumed in planning this survey, such a sample is suitable to address the problem of the spatial distribution of low-mass galaxies relative to the structures delineated by bright galaxies. A more detailed study of the already-known differences between luminous and faint galaxies (e.g., Salzer 1989; Pustilnik et al. 1995; Popescu et al. 1997) could help to gain a deeper understanding of the CDM structure N-body simulations (e.g., such as by Mathis & White 2002; and Gottlöber et al. 2003).
One of the aims of the HSS project was to close the gap between the sky regions of the SBS and Case survey. This goal has now been reached. Since the HSS has intersections with both, the possible differences in their selection functions and other sample characteristics can be quantified and accounted for. Thus, the useful ranges of galaxy parameters for which ELGs can be studied in the whole region of sky covered by the SBS, HSS and Case will be obtained.
Another important aspect of the BCG studies is related to the starburst triggering mechanisms. It can be addressed with the HSS sample to check the preliminary conclusion about the important role of galaxy interactions made, e.g., on the large sample of the SBS BCGs (Pustilnik et al. 2001). Similar studies have been conducted on the samples from the 2dF Survey project by Lambas et al. (2003) and Alonso et al. (2004).
One more interesting aspect of statistical studies of this BCG sample is related to the high S/N ratio spectra of the subsample of the most strong-lined ELGs. This allows us to determine in a large sample of galaxies the abundance of oxygen and other heavy elements by the classic method, and to use these data to compare the BCG properties with the models of galaxy chemical evolution. We already obtained such data for a significant fraction of the HSS BCG sample (15% of all BCG and 40% BCGs with the EW([O III]5007) > 150 Å, Pustilnik et al. 2004b). Some of the strong-lined HSS BCGs were used in the new primordial helium determination (Izotov & Thuan 2004).
One of the goals of the HSS project was to search for new extremely metal-deficient (XMD) BCGs (those with 12+ (O/H ), possible analogs of candidate young galaxies, like I Zw 18 and SBS 0335-052. Altogether, in addition to the two XMD galaxies in this zone known from previous studies (1415+437=CG 389 and 1224+3756=CG 1024), eight new such galaxies are found (see papers by Kniazev et al. 1998,2000a,b; Pustilnik et al. 2004a,b; Guseva et al. 2003). Thus, the fraction of XMD BCGs at the magnitude limit of the HSS is 2% (as already claimed by Pustilnik et al. 2003), about 1.5 times higher than the fraction found by Kniazev et al. (2003) for the Sloan Digital Sky Survey (SDSS). While we are dealing with small samples in either case, which makes this ratio uncertain, the ratio still indicates that we succeeded in creating a design for the HSS which is more sensitive in finding XMD BCGs than general galaxy surveys like the SDSS.
We performed the follow-up spectroscopy of the sixth and last list of candidates for ELGs (mainly of H II type) from the Hamburg/SAO Survey. Summarizing the results, the analysis of the spectral information and the discussion above we draw the following conclusions:
This work was supported by the grant of the Deutsche Forschungsgemeinschaft No. 436 RUS 17/77/94 and by the Russian Federal Program "Astronomy''. U.A.V. is very grateful to the staff of the Hamburg Observatory for their hospitality and kind assistance. Support by the INTAS grant No. 96-0500 is gratefully acknowledged. I.M. and J.M. acknowledge financial support by DGICyT grants AYA2001-2089 and AYA2003-00128 and the Junta de Andalucía. The authors thank the anonymous referee for useful comments and suggestions. The use of APM facility was very important for selection methods for additional candidates to BCGs from the 2nd priority list. This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. We have also used the Digitized Sky Survey, produced at the Space Telescope Science Institute under government grant NAG W-2166.