next previous
Up: Multi-spectral study of a


Subsections

3 Observations and data reduction

3.1 Observations

The observations were made at the Southern European Southern Observatory (La Silla, Chile) by Doublier and Caulet using the 1.54 m Danish Telescope, during two periods: April 29 - May 3 and December 20-23, 1995. In April-May, the direct camera was equipped with a thinned back-illuminated 1024$\times$1024 pixel Tektronix CCD (ESO #28) with a pixel scale of 0.36''. In December, DFOSC (Danish Faint Object Spectrograph and Camera) was equipped with the thinned Loral 2048$\times$2048 CCD#17 with a pixel scale of 0.38''.

Each galaxy field was exposed three times with a dithering step equal to at least the size of the galaxy in order to correct for blemishes on the CCD chip. Short exposures were taken to avoid saturation of bright stars in the fields that could create ghosts and uneven background, CCD saturation and bleeding. The data were bias subtracted, flat fielded and combined in the standard way, using the IRAF package.

During both runs, the weather was clear, and the seeing was between 0.9 and 1.1'' (FWHM measured on field stars close to the objects). Weather conditions did not vary much during the nights, the variations of airmass-dependent terms in Eqs. (1) and (2) were less than 0.005. In addition, we observed the galaxies close to meridian, to reduce the variation of atmospheric absorption. Therefore, seeing and/or absorption were not important effects when we combined the frames, even those taken during different nights.

We used Johnson B and R ESO filters in combination with the CCDs, producing instrumental systems that are easy to convert to the standard Cousins system. The color terms were found to be negligible for both filters. For the photometric calibration, some [Landolt(1992)] equatorial standard stars were followed throughout the nights at different airmasses. Table 2 summarizes the observation log, as well as the signal-to-noise ratio obtained at the surface brightness of 25 mag/arcsec2 in both bands.

  
Table 2: Observing log

\begin{tabular}
{l c c r@{~\&~}l}
\hline
\multicolumn{1}{l}{object Name}&
\multi...
 ...(25): Signal to noise ratio at $\mu$\space = 25 mag/arcsec$^2$.}\\ \end{tabular}

3.2 Standard stars

Aperture photometry was performed on standard Landolt stars. We obtained the following equations which are used to convert our instrumental magnitudes (br) and colors to the Johnson-Cousins system, the values in parenthesis correspond to those obtained during the December run:

\begin{eqnarray}
B = b_{\rm April, (December)} - 0.27 \times {\rm airmass} + 22....
 ...{(F_{ b,\rm cor} )} \\  
r = -2.5 \times \log {(F_{ r,\rm cor} )} \end{eqnarray} (1)
(2)
(3)
(4)
where $F_{ b,\rm cor}$ and $F_{ r,\rm cor}$ are the sky subtracted fluxes in counts per second, in B and R bands respectively.

The atmospheric absorption coefficients and zero points magnitudes were determined each night in a standard way. The values given in the above equations are average values for each observing period.

The Johnson R filter includes the H$_\alpha$ emission line in the rest frame; all objects have low redshift, and their H$_\alpha$emission is always included in the measured R flux. Over the galaxy, the H$_\alpha$ flux is less important than a substantial background from an evolved stellar population, but it may, in some cases, produce a local red excess inside the galaxies, especially in star forming complexes that may have a large equivalent width of H$_\alpha$emission.

3.3 Galaxy photometry

3.3.1 Data reduction

The surface photometry of the galaxies was done with our own MIDAS (see Paper I) macro-procedures, and the reader is referred to Paper I for a complete description of the methods used. We briefly summarize the method: an isophotal integration of the flux was done without setting any constrain on the geometry of the isophotes (equivalent profile method, [De Vaucouleurs1959]; [Fraser1977]). After a careful estimation of the local sky background, the following photometric parameters were derived (using the notation from [De Vaucouleurs1959]): the equivalent radius $r_{\rm eff}$ and r0.25, within which are contained half and a quarter of the total luminosity respectively; they were derived from the curve of growth of the instrumental magnitudes, before Galactic foreground extinction correction. The asymptotic magnitude was extrapolated from the last 4 points on the integrated magnitude curve of growth. The surface brightnesses at $r = r_{\rm eff}$ and r0.25, and the average surface brightness inside the effective radius, $<\mu_{\rm eff}\gt$, were also measured. We have checked that the variation of the derived parameters with different values of extinction is not large, considering the errors on the measurement of the radius, and the estimation of the asymptotic magnitude. For instance, adding an arbitrary value of AB = 0.3 mag to the galaxy leads to decrease the effective radius by 0.2 arcsec. This is not significant considering the average seeing of 1.1'' we had throughout our runs.

The sky background was checked to be flat in the vicinity of the objects. Since our procedures to derive the surface photometry take into account the local sky background, we wanted to avoid sky background subtraction on the image. In the cases of sky gradients under the objects, we subtracted a two-dimensional fit of the sky from the images.

As in Paper I, to increase the signal-to-noise ratio of the faint outermost isophotes, we have performed the isophotal integration in two steps. First, we used the direct image to produce a surface brightness profile, valid at bright intensity levels, then we smoothed the images with a Gaussian filter of a few pixels of width that matched the seeing disk of the images, and we built a second profile. Finally, the two surface brightness profiles were merged.

3.3.2 Errors

The errors in magnitudes and surface brightness were estimated assuming Poisson statistics, and taking into account the errors on the sky background estimation (based on local 1-$\sigma$ variation). We used the relation given by [Saglia et al.(1997)].

The errors on the asymptotic magnitudes are more difficult to estimate accurately because pure Poisson statistics do not apply, and due to the extrapolation of a magnitude based on the last four points of the photometric profile. We estimate that the errors are of the order of those made on the photometric calibration, e.g. 1%.

 
\begin{figure}
\includegraphics [clip]{figure1a.ps}
\end{figure} Figure 1:  a) Atlas of isophotal maps and surface brightness profile distributions. This atlas presents the isophotal B maps and surface brightness distribution in the B and R bands. On the top right corner, the B-R color distribution is shown. The scale bar at the bottom of the map represents 1 kpc. The threshold brightness is indicated, the isophotal step on the map is 0.5 magnitude. Orientation is North up East left for the galaxies taken in December 1995, and North down and East right for the images taken in April-May 1995

 
\begin{figure}
\includegraphics [clip]{figure1b.ps}
\end{figure} Figure 1: b) continued

 
\begin{figure}
\includegraphics [clip]{figure1c.ps}
\end{figure} Figure 1: c) continued

 
\begin{figure}
\includegraphics [clip]{figure1d.ps}
\end{figure} Figure 1: d) continued

 
\begin{figure}
\includegraphics [clip]{figure1e.ps}
\end{figure} Figure 1: e) continued

 
\begin{figure}
\includegraphics [clip]{figure1f.ps}
\end{figure} Figure 1: f) continued

 
\begin{figure}
\includegraphics [clip]{figure1g.ps}
\end{figure} Figure 1: g) continued

 
\begin{figure}
\includegraphics [clip]{figure1h.ps}
\end{figure} Figure 1: h) continued

 
\begin{figure}
\includegraphics [clip]{figure1i.ps}
\end{figure} Figure 1: i) continued

 
\begin{figure}
\includegraphics [clip]{figure1j.ps}
\end{figure} Figure 1: j) continued

 
\begin{figure}
\includegraphics [clip]{figure1k.ps}
\end{figure} Figure 1: k) continued

 
\begin{figure}
\includegraphics [clip]{figure1l.ps}
\end{figure} Figure 1: l) continued


  
Table 3: Results of the photometry

\begin{tabular}
{l ccccccc}
\hline
\multicolumn{8}{l}{}\\ [-3mm]
\multicolumn{1}...
 ...ce colors, defined as the $r_{\rm eff} / r_{0.25}$\space ratio.}\\ \end{tabular}

3.3.3 Photometric results

The results of the photometry are presented in Fig. 1 and Table 3. For each galaxy, Fig. 1 shows the isophotal map in B, the surface brightness distribution in the B and R bands, and the B-R distribution. The orientations are from Haro 14 to Tol 3 North to the top and East to the left, and from Fairall 6 to Tol 1937-423 North down and East right. Table 3 collects the observable parameters derived from the photometry.

The data in Table 3 have not been corrected for Galactic foreground extinction. Extinction has been applied when deriving the integrated parameters such as the absolute magnitudes in B and R (MB and MR resp.). We used the E(B-V) maps of [Burstein & Heiles1982] and derived the absorption coefficients using Savage & Mathis (1979) relations: AB = 4.0 E(B-V) and AR = 2.52 E(B-V).

3.4 B-R color profiles

The method used to trace the B-R color profiles is the same as in Paper I. We interpolated the surface brightness (SB) profiles in B and R because the sampling of the profiles are different in each band, and subtracted the resulting profiles from one another. The results are in excellent agreement with those obtained directly by performing the photometry on B-R two-dimensional maps of several galaxies (Mk 996, Haro 14 and Fairall 301), only of better signal-to-noise toward the faint regions. Since we did not correct neither the profiles, nor the images from the seeing effects, the profiles are not to be seriously considered within the central first 1'' radius.

The errors are calculated assuming that the errors in B and R are independent, using:
\begin{eqnarray}
\sigma_{B-R} = \sqrt{\sigma_{B}^2 + \sigma_{R}^2}.\end{eqnarray} (5)
Here, we should emphasize that, in the central part of the BCDGs, the B-R colors could be contaminated by nebular emission such as H$_\alpha$ present in the R band filter. The equivalent width of this line reaches generaly a few hundreds of Angström, therefore it can contribute significantly within the R filter to the total R emission. The outer parts of the BCDGs, however, do not suffer from such contamination, since the emission can be supposedly attributed to the stellar component only.

3.5 Comments on individual objects

Most of our objects show several star forming regions or knots with B-R colors different from those of their host galaxies. We note that the knots get bluer away from the photometric center of the galaxy. This behaviour could reflect a variation of internal reddening due to dust, or more likely an age spread. However, those objects do not show obvious presence of dust as the ratio between the recombination lines of hydrogen H$\alpha$/H$\beta$ remains close to a value of 3. Some galaxies have a single compact star forming knot ($\le$ 5-10''), whereas , in other galaxies, "multi-knots'' regions are distributed over more than 2/3 of the galaxy's surface.

As in Paper I, the morphology of several BCDGs in the present sample departs conspicuously from axisymmetric shapes, for instance Tol 1937-423, UM 461, in the present sample. Because irregular morphology in star forming galaxies can be explained by mergers and dynamical interaction, and because our galaxies are not associated with known companions (except UM 465), we have searched for faint galaxy companions close to our BCDGs on our deep optical images. The results are presented in Sect. 4.3. We discuss now in detail the properties of individual galaxies. Other distortions (boxy isophotes, isophote apparent major axis rotation) are noted in the inner regions of several galaxies.

A summary of the anomalies noted in the BCDGs studied in our sample is given in Table 4 at the end of this section.

3.6 The surface brightness distributions

In Paper I, we concluded that the 23 BCDGs from the Byurakan Surveys could be subdivided in three groups according to their surface brightness distribution: 10 objects clearly dominated by a "spheroidal'' component fitted by a dominant de Vaucouleurs (r1/4) law, 7 objects clearly dominated by an exponential brightness distribution, and 6 composite or unclassifiable profiles. With this limited northern sample, we found that about 3/4 of our objects belong to the two first categories.

In the present sample, we find a similar repartition of our BCDGs in the three groups just defined: 12 objects are dominated by an exponential component (7 are fitted by a pure exponential law), 12 are dominated by a r1/4 law (5 are fitted by a pure r1/4 law). By "pure'' we mean that no other significant component has to be added to account for all the light of the galaxy. For the total sample of 44 Northern and Southern galaxies (1 companion included), we find that: 11 galaxies (25%) have SB profiles following closely a pure exponential distribution, and 8 (18%) galaxies with composite profiles in which the exponential component dominates; 9 galaxies (20%) have SB profiles following closely a pure r1/4 distribution, and 12 galaxies (27%) with composite profiles in which the r1/4 component dominates; finally, 3 galaxies are unclassifiable (7%), those peculiar galaxies are Mk 1499, Mk 1131 and SBS 1331+493 (Paper I).

Therefore, we support that the global structure of BCDGs, traced by their projected luminosity density, is not different, in terms of dynamical components, from the global structure in normal galaxies where two main stellar populations dominate. Those components are thick disks (whose true axial ratio distribution remains to be determined in the disk dominated BCDG galaxian population, see Sung et al. 1998), and spheroidal components obeying the r1/4 law (implying a significant degree of relaxation). The present 20 disk-dominated BCDGs will be used to derive tight constraints on the true axial ratios of the disks, this will be treated in a following paper. In any case, kinematical data are also necessary, in addition to photometry of a larger sample of objects. Regarding the dynamics of the BCDGs, both exponential and r1/4, we insist that a deeper analysis requires kinematical data from the gas and stellar populations.

The use of the light distribution to trace the mass distribution can be done to some extent, keeping in mind that the mass-to-light ratio usually varies within the galaxies, as indicated by the color gradient in BCDGs ([Papaderos et al.1996b]; [Doublier et al.1997] and this work). The central starburst dominated regions have mass-to-light ratio close to 0.1 ([Charlot et al.1996]) in the visual, while the outer old star population dominated regions ([Thuan1983]; [Hunter & Gallagher1985]; [Doublier et al.1999]) would have a mass-to-light ratio larger than 1. As a result, the mass distribution would flatten towards the center compared to the light distribution. Nevertheless, studies of the light profile in the near-infrared ([Doublier et al.1999]) shows that BCDGs with "optical'' r1/4 law profiles display r1/4 law profiles in the K band where the old stellar population dominates without a light excess in the central parts. This leads us to believe that the r1/4 law is intrinsic, rather than being due to the light excess caused by the presence of the star formation regions.


  
Table 6: r1/4 galaxies; north and south samples


\begin{tabular}
{l r@{~$\pm$~}l r@{~$\pm$~}l c c c }
\hline
\multicolumn{1}{l}{O...
 ... 10.49&C \\  
Mk 1426 & --& 8.92&0.09 & -- & 11.82& \\ \hline\hline\end{tabular}


In Tables 5 and 6, we summarize the parameters derived from our photometric analysis. The surface brightness profiles were fitted uniquely either by a r1/4 law or an exponential law. If neither case applies, a note in the last column of Tables 5 and 6 mentions that the surface brightness is "composite''.

The parameters of the exponential law and of the de Vaucouleurs (r1/4) law were estimated using a crossed-linear regression applied to the regions that were not affected by the central "bulge-like'' component (in case of the exponential law) and seeing (in case of the de Vaucouleurs law), and excluding the outermost isophotes where the signal-to-noise ratio is too low. In a following section, we discuss the "extra'' components such as the central excess of light, or any possible excess of light present in the low S/N regions of the galaxies.

The excess of light in the central regions is probably due to the star forming regions: the method described by De Vaucouleurs (1959) we used to derive our surface photometry defines the center as the photometric center, while for ellipses fitting of the isophotes the center of the galaxy is defined as the center of last inner ellipse. We have tried the ellipses fitting on some of the most clumpy BCDGs in our sample, and the results are sensibly in agreement (within a few percent) in the outer regions, while the fit diverge considerably in the central regions.

The excess of light compared to a pure exponential law, or to the de Vaucouleurs law, seen in the outer parts of some BCDGs is not likely due to large errors (see surface brightness distributions in Fig. 1). In some cases, it could be attributed to an old, very low surface brightness component, whether this component could be exponential or r1/4 is not clear. This component is detected in few cases (UM 465A, Fairall 301) in the near IR ([Doublier et al.1999]).

Figure 2 displays the surface brightness distributions of the r1/4 BCDGs in the [$\mu$, $(r / r_{\rm eff})^{1/4}$] plane, where $\mu$ is the surface brightness in mag/arcsec2, and $r_{\rm eff}$ the effective radius defined in Sect. 3.3.1 (a similar plot can be found in Paper I, for the northern sample). The plotted line represents the relations between the reduced radius and the surface brightness for a pure de Vaucouleurs law. Most of the r1/4 profiles show a much steeper slope than the de Vaucouleurs relation. We discuss this difference in terms of the surface brightness gradient in the next section.

 
\begin{figure}
\includegraphics [width=8.5cm,clip]{figure3.ps}

\includegraphics [width=8.5cm,clip]{figure3bc.ps}
\end{figure} Figure 3:  Relationship between the scale-length in B and R band. a) Relationship between de Vaucouleurs law slopes. The solid line represents the line of equal slopes in B and R. The relationship is obviously steeper than 1, the coefficients are systematically larger in B than in R, and systematically larger than the usual value 8.327. The large star represents the mean values of the subsample (10.0 in B and 9.7 in R) while the large triangle shows the empirical values of 8.327 in both bands. b) Relationship between the scale-lengths of the disk galaxies (or composite when the disk was easy to separate) in B and R bands. The solid line represents the line of equal scale lengths in B and R. The scale-lengths in both colors are clearly correlated. b, c) Variation of the $\chi^2$ as a function of the Equivalent. Radius. b) Variation of the $\chi^2$ for the elliptical-like galaxy: II SZ 34. The $\chi^2$ was obtained by dividing the surface brightness distribution modeled by the best fitting r1/4 law, by the observed surface brightness distribution. Note the local maxima, in B and R, corresponding to inner features of the galaxy. c) $\chi^2$ variation for the disk galaxy: UM 465, obtained in the same way as above. Note the local maxima of the $\chi^2$ that are, unlikely to the elliptical case, behaving differently in B and R

3.6.1 "De Vaucouleurs'' profiles (23%)

Figure 3 (top panel) shows the relation between the observed slopes of the r1/4 laws in reduced coordinates in R and B. The slope values are systematically larger than the canonical value of 8.327 for the law. Furthermore, the slopes are clearly larger in B than in R. This was also observed in our northern sample ([Doublier et al.1997]) for which we suggested either having systematically overestimated the effective radii or the presence of an homogeneous component in addition to the dominant spheroid. The systematic excess over the 8.327 value and the fact that this excess is larger in B than in R leads us to conclude that the presence of the starburst "component'' is probably responsible for this effect.

The most natural explanation is that the effective radius is overestimated. The overestimation of $r_{\rm eff}$ could result from errors in extrapolating the asymptotic magnitudes; however, we performed the surface photometry several times showing that the internal consistency of the asymptotic magnitude values is good. Then the excess $r_{\rm eff}$ values could result from the presence of a small compact component of a small scale length in the central part of the BCDG.

To investigate this hypothesis, we have simulated 3 types of galaxies: a pure elliptical, a dominant elliptical with a small additional r1/4 central component, and a dominant elliptical with a small exponential central component. Seeing effects have been added by means of a convolution of the model galaxies with a Gaussian. The model images have subsequently been submitted to the surface photometry procedure, yielding a clear tendency to overestimate the resulting $r_{\rm eff}$ of the "Host galaxy + Starburst component'' whatever the nature of the starburst component.

To check the consistency of the fit to the r1/4 law, we have performed $\chi^2$ tests by subtracting the model r1/4 SB distribution that best fits from the observed SB. Figures 2b and 3c displays, as an example, the results on Mk 996 and II SZ 34 in B and R. The value of $\chi^2$ for II SZ 34 shows local maxima at 5, 15 and 30'' from the center, implying local significant departures from the pure r1/4 law in these regions. These local departures originate from the fact that not all the star formation regions are concentrated in the central parts of the galaxies. This is commonly encountered in almost all "spheroidal'' BCDGs of our sample.

3.6.2 Exponential profiles (25%)

The mean value of the scale length parameter is 1.63 kpc-1 without taking into account the scale length value of the first component of UM 461 (7.57 $\pm$ 0.1 (kpc-1)) that may be due to the two very bright nuclei.

Figure 3 (bottom panel) shows the relation between the scale length in B and R. The relation shows that the scale lengths in B and R are well correlated. The disks should have similar stellar distribution in R and in B. In some cases, the B scale length is significantly larger than the R scale length.

Figures 3b and 3c show the $\chi^2$ variation for UM 465A, as a function of the equivalent radius, in B and R. As for the de Vaucouleurs profiles, the $\chi^2$ show small variations larger than the mean scale length that account for the existence of small structures within the disk of the galaxy: off-centered star formation regions, foreground faint stars, and globular clusters.

The disks and de Vaucouleurs parameters describing the galaxies are listed in Tables 5 and 6 respectively.

3.6.3 Composite profiles

We call "composite'' SB profiles that cannot be fitted by a pure exponential or a pure r1/4 law. The shapes of these profiles (in B and R) can be very different from one galaxy to another. In most cases, two components are easily seen, the observed profiles being the sums of an exponential and a r1/4, as in giant spirals, with one of these components being more or less dominant. However, for some other galaxies, the profiles displayed in the [$r_{\rm e}$, $\mu$] plane show a convex shape but when displayed in the [r1/4, $\mu$]plane, the profiles exhibit a concave shape. Also, these galaxies show very disturbed isophotes. These systems may be relaxing merger remnants, or they could have a thick disk not yet settled inside a host elliptical galaxy.


next previous
Up: Multi-spectral study of a

Copyright The European Southern Observatory (ESO)