next previous
Up: Studies of the

3. Photometry

After the acquisition of the cluster sample our goal was to determine the apparent tex2html_wrap_inline1553-magnitude and other fundamental parameters for each of the program galaxies. Unfortunately, the du Pont film (and plate) was not suitable for this purpose due to the lack of sensitometer marks. Thus, we decided to do photometry with a first generation glass copy of the Schmidt plate SRC 323J from the archive of the ESO headquarters in Garching. The plate is equipped with 7 step sensitometer marks on its northern and eastern margin allowing to establish the density-to-intensity transformation.

3.1. Tracings

The part of the Schmidt plate 323J which coincides with the surveyed field was digitized with the Perkin Elmer PDS 1010A 2-D micro-densitometer system at the ESO headquarters. After the calibration of the instrument with an unexposed area (plate fog) at the edge of the plate, 16 subfields covering the survey region were scanned in the density mode with a 12 bit A/D converter. 25 tex2html_wrap_inline1555m were chosen for the size of the squared beam aperture as well as for the step size of the scan. With a plate scale of 67.5 arcsecs/mm and a seeing of 2.3arcsecs for the particular ESO Sky Survey plate, this led to a slight oversampling. The seeing value was derived from the faintest stars on the plate. The scanning velocity is inversely correlated to the maximum reliable density that can be registered by the system. We chose tex2html_wrap_inline1557 20mmstex2html_wrap_inline1559 enabling us to measure densities up to D=3.1. This density limit is somewhat higher than the expected saturation level due to photographic processes at tex2html_wrap_inline1563 or tex2html_wrap_inline156521.5Barcsecstex2html_wrap_inline1567. In this way, the dynamical range of the plate was covered best.

3.2. Characteristic curve

Additionally to the science field the two strips of sensitometer marks were fully traced. In Fig.2 (click here) we show plots of their smoothed profiles averaged over 50 scanning rows. While the northern strip yields very well defined density steps, with the exception of the highest, the average profile of the eastern strip is of only poor quality. Several density steps show gradients superposed by strong defects. No accurate measurements of these density values were possible. This circumstance together with the fact that no log(I)-values were available for the eastern projector, neither in the UKSTU-Handbook (Tritton 1983) nor from the Royal Observatory in Edinburgh, constrained us to work with the northern marks only.

To establish the characteristic curve of the plate we used the bi-logarithmic formula proposed by Llebaria & Figon (1981):
displaymath1569
where tex2html_wrap_inline1573 of the original formula is replaced by the more convenient quantity I=E/t. The different parameters have the following meaning: I and D are the relative intensity and the density, respectively. tex2html_wrap_inline1581 is the density at which the photographic plate becomes saturated, tex2html_wrap_inline1583 is the zero-point density of the photographic emulsion, and C is the zero point. The analytic curve was tex2html_wrap_inline1587-fitted at the points defined by the density values of the northern step marks and their corresponding log(I)-values (Tritton 1983). As already mentioned above, the highest step of the northern sensitometer mark was not very well defined. Consequently, a first fit with this original density value was unsatisfactory, i.e. the derived value of the fit parameter tex2html_wrap_inline1591 was lower than the highest scanned density found on the plate. To overcome this problem empirically we applied a 5% correction to the highest step value from 3.039 to 3.179 in order to reach the highest plate density. Figure3 (click here) shows the curve which fits best the finally used data points and in Table3 (click here) we list the best-fitting parameters.

  table325
Table 3: Parameters of the characteristic curve

  figure336
Figure 3: The scanned density values of the northern sensitometer marks versus the corresponding log(I)-values (plus sign). For the fit of the characteristic curve (line) we used the values indicated by circles. The highest measured density had to be slightly corrected for this purpose

For an external quality check of the density-to-intensity transformation Bender (from the CCD work of Bender 1994) kindly provided us with the major and minor-axis V-band surface-brightness profiles of the bright early-type galaxies NGC 4696, 4709, and 4729 in Centaurus. These profiles were compared with azimuthally averaged profiles in the B-band we derived for these galaxies (cf. next subsection). As no significant radial colour gradients are observed for early-type galaxies in general (e.g. Reid et al. 1994), our profiles should fit between Bender's major and minor-axis profiles after applying an individual zero-point offset correction. This is illustrated in Fig.4 (click here). Obviously a good transformation quality is achieved in the non-saturated regime from tex2html_wrap_inline1613Barcsectex2html_wrap_inline1615 (tex2html_wrap_inline1617) down to tex2html_wrap_inline1619Barcsectex2html_wrap_inline1621. The accuracy at a fainter level, between 25.0Barcsectex2html_wrap_inline1625 and tex2html_wrap_inline1627Barcsectex2html_wrap_inline1629 (which relies mostly on the accuracy of the sky background estimation) can be judged from the linear shapes of the surface-brightness profiles of faint dwarfs shown in Fig.5 (click here).

  figure349
Figure 4: The major and minor-axis surface-brightness profiles in the V-band for NGC 4969, 4709, and 4729 are shown as lines. Our azimuthally averaged surface-brightness profiles are added as dotted circles. The best-fitting offset from the B-band to the V-band magnitude system was determined individually for each galaxy

Each individual zero-point offset is the sum of the calibration constant of our instrumental magnitude to the B-magnitude system and the B-V galaxy colour which could be used to calibrate our magnitude system in principal. But neither tex2html_wrap_inline1641 nor B-V are very well known for the three galaxies (Sadler 1984; Poulain 1988; Dressler et al. 1991; Longo & de Vaucouleurs 1983; de Vaucouleurs et al. 1991; Prugniel et al. 1993).

3.3. Image reduction

The neighbourhood of each cluster galaxy image was disentangled interactively from surrounding stars. To do this we used procedures developed in the image processing program MIDAS. If possible, stars projected on a galaxy surface were eliminated by taking advantage of the symmetric property of the galaxy. The pixels of such a star were substituted by the pixels lying on the point symmetrical opposite of the galactic centre. In cases where this technique was not applicable, e.g. the opposite pixels were also contaminated by stars or the galaxy showed asymmetrical appearance, the underlying galactic ground was modelled by fitting locally a second order 2D-polynomial.

After the galaxy image was star subtracted we determined for the galaxy a growth curve based on a visually selected centre. The pixel intensities were integrated in concentric round apertures of increasing radius in steps of 1 pixel (=1.68'') outwards. Simultaneously, the sky intensity was subtracted so the growth curve became asymptotically flat for sufficiently large radii and a correct sky background determination. This first growth curve was used to define the maximum radius where all light of the galaxy is included and the noise of the background starts to dominate. A new galaxy centre was defined by the luminosity-weighted first moment of all pixels within the maximum aperture and the second (final) growth curve tex2html_wrap_inline1655 was established for the galaxy.

Apart from few cases this growth curve was smooth enough to derive the corresponding instrumental surface-brightness profile:


displaymath1649

The profile was approximated analytically in three possible ways depending on the Hubble type of the galaxy. Classically, the light profiles of ellipticals or the spheroidal components of disk galaxies were described by the tex2html_wrap_inline1657-law (de Vaucouleurs 1948). A more general approach represents the generalized exponential law (Sérsic 1968): tex2html_wrap_inline1659 or tex2html_wrap_inline1661, where the exponent remains free instead of being fixed at n=1/4. It has been shown recently (Caon et al. 1993; Graham et al. 1996) that this formula offers a much better approach to the observed profiles of elliptical galaxies. We were thus motivated to fit most of these profiles by this analytic form. The shape parameters n are listed in the CCC.

S0's and spirals consist of two different profile components. The inner bulge follows the de Vaucouleurs law as the ellipticals. The outer disk part shows a linear decay of the surface brightness corresponding to an exponent n=1 in the generalized exponential law. Consequently we fitted a two-component model at the profiles of these galaxies.

Profiles of some bright nucleated dwarf galaxies (dE, N and dS0, N) exhibit a clear bipartition. In these cases we approximated the profile by a two-component exponential model (Binggeli & Cameron 1991). Generally, a classical exponential law was used for the dE's, dS0's, Im's, and BCD's as well as for some dE, N's and dS0, N's. Note that the nuclei of dE, N's or dS0, N's in Centaurus have a typical size of tex2html_wrap_inline1669. Therefore, the seeing effect and the discreteness of the scanning process can easily destroy the central light excess in the empirical surface-brightness profile.

The innermost radius limit for all fits was 3arcsecs and a lower limit in instrumental surface-brightness was applied equivalent to 27Barcsectex2html_wrap_inline1673. For some bright galaxies the centre of the image was saturated due to photographical limitation. Here, the fit was restricted to the undisturbed outer part of the profile starting at an equivalent surface brightness of 21.5Barcsectex2html_wrap_inline1677.

The model-dependent parameters tex2html_wrap_inline1679 and tex2html_wrap_inline1681 were derived from the pure exponential outer part of a profile. As the surface brightness scale is logarithmic, an exponential profile appears as a straight line. The linear extrapolation of the best-fitting line into the centre yields the central exponential surface brightness tex2html_wrap_inline1683. The exponential scale length tex2html_wrap_inline1685 of the galaxy corresponds to the slope of this line.

Moreover, we determined three model-independent parameters for each galaxy: the instrumental total magnitude tex2html_wrap_inline1687, the effective radius tex2html_wrap_inline1689 (radius containing half of the total light), and the mean effective surface brightness tex2html_wrap_inline1691. For those galaxies with no saturation problems in the surface-brightness profile we analysed the observed growth curve for this purpose. In particular, the asymptotic limit of the growth curve at maximum aperture defined the total instrumental intensity of the galaxy. For all other cases, the parameters were derived from the growth curve corresponding to the best-fitting surface brightness model.

  figure383
Figure 5: Selected non-smoothed mean surface-brightness profiles in order of decreasing total apparent magnitudes are shown. The CCC numbers are listed, along with the galaxy classifications. The points are plotted at 1'' intervals. Notice the saturation effect in the inner part of the CCC130 profile and its correction by the best-fitting model

As already mentioned above, there were few galaxies, mainly faint early-type dwarf galaxies, where the growth curve was so badly affected by sky background noise that a derivation of the surface-brightness profile would require to run a smoothing algorithm first. All structure parameters have been derived from the best-fitting model of an integrated pure exponential law into the noisy growth curve starting at r>3'':


displaymath1650
with tex2html_wrap_inline1697, I the intensity, and tex2html_wrap_inline1701 the exponential scale length.

For illustration we present in Fig.5 (click here) an arbitrarily chosen collection of eight calibrated surface-brightness profiles. In the upper left diagram the saturation effect due to the limits of the photographic plate is clearly visible. A fit of a generalized exponential law to the undisturbed outer part of the profile was used to extrapolate into the central region.

3.4. Vignetting correction

All our magnitudes were corrected for the geometrical vignetting effect of the UKSTU Schmidt telescope. We used the experimentally determined vignetting function (cf. Dawe & Metcalfe 1982; Tritton 1983). In general, the correction factor varies on the photographic plate from 1.0 at the centre to 1.25 at 4 degrees off axis. At the position of the Centaurus cluster on plate 323J a typical correction of 3% had to be applied.

Because of the difficulties to quantify the effect of the desensitization of hypersensitized plates ("Malin effect''), no correction was applied for this additional effect. It has been studied in detail by Campbell (1982) and a maximum error of less than tex2html_wrap_inline1705 is quoted by Dawe & Metcalfe (1982).

3.5. Magnitude zero point determination

The goal was to transform the instrumental magnitudes tex2html_wrap_inline1709 into CCD B-magnitudes. For this purpose two sets of independent data were used from the literature. A selected sample of bright galaxies was taken from the Surface Photometry Catalogue of the ESO-Uppsala Galaxies (Lauberts & Valentijn 1989 hereafter LV). These magnitudes proved to be in good agreement with the standard RC3 system (Paturel et al. 1994). The selection criteria we applied to find those galaxies which are most reliable for the calibration were twofold: First, the central part of the surface-brightness profile must not be affected by saturation effects. Second, our derived value for tex2html_wrap_inline1713 should not differ too much from the effective diameter tex2html_wrap_inline1715 given by LV. The motivation for the second constraint is based on an existing discrepancy between LV and other authors about the tex2html_wrap_inline1717 values of an early-type galaxy sample (cf. Fig. 9 in LV). Based on this uncertainty we arbitrarily fixed the allowed maximum deviation between tex2html_wrap_inline1719 and 2tex2html_wrap_inline1721 at 15%. By this selection we expect to take into consideration only the uncontroversial galaxy data of both studies. The two constraints were met by 10 galaxies.

LV offers only few galaxies fainter than tex2html_wrap_inline1723. Thus, we supplemented the calibrator sample with a CCD-based data set of 22 dwarf galaxies (Bothun et al. 1989 hereafter BCS89). This second sample covers well the magnitude interval tex2html_wrap_inline1725 and will improve the calibration accuracy at the faint end. The two data sets of LV and BCS89 galaxies are collected in Table4 (click here).

 
 table413

Table 4: Calibration data

  figure431
Figure 6: The upper graph shows the calibration diagram based on 32 selected galaxies in common with LV (circle) and BCS89 (triangle). The indicated line is not a fit, but a line with unity slope representing a perfect relation between two magnitude systems. Obviously, this line is very well approximated by the data points. The residual diagram for our tex2html_wrap_inline1745-magnitudes is given in the bottom panel

  figure436
Figure 7: tex2html_wrap_inline1747-magnitudes of all galaxies in common with LV. Filled circles are galaxies used for the magnitude calibration. Open circles are galaxies required individual corrections for central saturation effects

We converted the CCD B-magnitudes from the BCS89 sample to the bandpass IIIa-J (tex2html_wrap_inline1751) using the colour equation tex2html_wrap_inline1753 as determined by Metcalfe et al. (1995). The approximately equivalent relation tex2html_wrap_inline1755 was applied on the tex2html_wrap_inline1757-magnitudes of LV. The calibration diagram is shown in Fig.6 (click here). A linear fit quantifies the transformation rule to be tex2html_wrap_inline1759. If restricted to the magnitude range of particular interest between tex2html_wrap_inline1761 and 21, the fitted line differs at most tex2html_wrap_inline1765 from a transformation rule with a unity slope. Based on this good agreement we assumed an underlying relation tex2html_wrap_inline1767 and determined the zero point with ZP=-23.23. This value was applied to transform our instrumental magnitudes into tex2html_wrap_inline1771-magnitudes. The rms scatter in Fig.6 (click here) is tex2html_wrap_inline1773 which translates to an uncertainty of the calibration zero point of tex2html_wrap_inline1775. Subsequently all tex2html_wrap_inline1777-magnitudes were converted into B-magnitudes using the colour equations given above where the measured colours were available. The average colour of the BCS sample tex2html_wrap_inline1781 was used where colours had not been measured. The last value is about equivalent to B-R=1.49 and comparable to the mean B-R colour of tex2html_wrap_inline1787 for all LV galaxies in common (not only the calibrators).

  figure456
Figure 8: Comparison of the model- depending parameters tex2html_wrap_inline1789 and tex2html_wrap_inline1791 for a sample of early-type dwarfs in common with BCS89. The solid lines represent perfect agreement between the two data sets

A weakness of photographic plates is the relatively small dynamic range they cover. Long-time exposures are leading to saturation effects on central parts of bright galaxies (cf. Fig.4 (click here)). The surface-brightness profile of such a galaxy had to be extrapolated to the centre in a way as described before. To get an idea how well this procedure worked, we compare in Fig.7 (click here) the tex2html_wrap_inline1793-magnitudes of the 10 uncorrected LV calibrators with those of other 28 galaxies in common with LV which required a saturation correction. We found good statistical agreement between the two samples, tex2html_wrap_inline1795 and tex2html_wrap_inline1797 with only a slight systematic error for the corrected sample of tex2html_wrap_inline1799. This result shows the reliability of our correction procedure within the expected magnitude error as described below. Furthermore, we note that the complete LV sample covers a wider colour range 0.8<(B-R)<1.7 than the calibrators. Apparently the transformation rule between the instrumental and B-magnitude systems seems to work well for all these galaxies.

3.6. Accuracy of galaxy parameters

The photometric accuracy of the instrumental magnitudes was estimated from the residuals illustrated in the bottom diagram of Fig.7 (click here). First of all it is satisfactory to see that there are no systematic magnitude differences between our magnitudes and those of the bright and faint calibrator samples, respectively: tex2html_wrap_inline1817 and tex2html_wrap_inline1819. We conclude that the global errors of our tex2html_wrap_inline1821-magnitudes for brighter galaxies are comparable to the typical uncertainty tex2html_wrap_inline1823 for LV magnitudes (Paturel et al. 1994). We further propose an error of tex2html_wrap_inline1825 for galaxies fainter than tex2html_wrap_inline1827.

BCS89 published the model parameters tex2html_wrap_inline1829 and tex2html_wrap_inline1831 for their 22 faint galaxies which can be used to estimate the accuracy of other structure parameters. In Fig.8 (click here) we plot the scale length and central exponential surface brightness data of our studies. Obviously, good agreement is found for tex2html_wrap_inline1833. The typical error is tex2html_wrap_inline1835 and no systematic offset is evident in the data. In the case of tex2html_wrap_inline1837 we find a mean difference of tex2html_wrap_inline1839 with a standard deviation of tex2html_wrap_inline1841.

Most of the dwarf elliptical galaxies exhibit in first order a pure exponential surface-brightness profile (e.g. Binggeli & Cameron 1991) which relate the effective radius and the scale length by the equality tex2html_wrap_inline1843. Taking advantage of this formula we estimate the error for tex2html_wrap_inline1845 to be roughly the same as that of tex2html_wrap_inline1847. The error of tex2html_wrap_inline1849 may be tex2html_wrap_inline1851 combining the uncertainties of tex2html_wrap_inline1853 and tex2html_wrap_inline1855.

Acknowledgements

This paper is part of the PhD thesis of one of the authors (HJ). He would like to thank his supervisors G.A. Tammann and Bruno Binggeli for their interesting ideas and invaluable advice. We are most grateful to Ralf Bender who made available to us galaxy surface-brightness profiles as well as to the ESO head quarter in Garching for giving access to the PDS plate-scanning machine. We thank Ken Freeman and the referee Dr. E. Bertin for having read the manuscript and having improved the origin version by many suggestions. HJ thanks the Swiss National Science Foundations for financial support.


next previous
Up: Studies of the

Copyright by the European Southern Observatory (ESO)
web@ed-phys.fr