Table 2 contains the best fit model for each image. Together with
the galaxy name (1), the filter (2),
and the referring image (3), with integration time,
and run ID, we list the inclination (4), the best fitting
function for the z-distribution (5), the calibration index
(6) (Ref. Sect. 2.4), and the central surface
brightness of the model (7), without correcting for
inclination. According to the distance tabulated in Table 1,
the cut-off radius
(8) is given in kpc and arcsec
as well as the scalelength h (9) and the vertical
scaleheight z0 (10) which is normalised to the isothermal
case, being two times an exponential scaleheight hz.
For the seven galaxies with available images in more than one filter, we
do not see any correlation of fitted parameters with different wavelength,
although we find the same inclination angle for the best fitted disk within
the range of the errors.
Appendix A shows the best fitting model as an overlay to selected
radial profiles for each image. The subsequent analysis of the distribution
for the different parameters concerning the formation and evolution of
galaxies will be given in forthcoming papers.
galaxy | filter | image | i | f(z) | cali. | h | z0 | |||||
[ ] | mag | [ ] | [kpc] | [ ] | [kpc] | [ ] | [kpc] | |||||
(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) | |||
ESO 112-004 | r | 40E3 | 87.5 | sech | l | 21.14 | 45.0 | -- | 22.7 | -- | 2.8. | -- |
ESO 150-014 | r | 20E3 | 90.0 | sech | l | 22.00 | 64.1 | 33.27 | 23.4 | 12.14 | 5.5 | 2.85 |
NGC 585 | R | 20L1 | 88.0 | sech | e | 21.76 | 70.0 | 24.49 | 32.7 | 11.44 | 10.6 | 3.71 |
ESO 244-048 | r | 15E3 | 87.0 | l | 20.54 | 45.4 | 19.17 | 13.7 | 5.78 | 7.2 | 3.04 | |
NGC 973 | R | 10L1 | 89.5 | e | 20.83 | 105.0 | 33.74 | 51.5 | 16.55 | 12.4 | 3.99 | |
UGC 3326 | R | 30L1 | 88.0 | sech2 | e | 21.50 | 101.5 | 28.45 | 70.4 | 19.74 | 4.7 | 1.32 |
UGC 3425 | R | 30L1 | 87.0 | sech2 | e | 21.01 | 80.5 | 22.26 | 29.2 | 8.08 | 7.4 | 2.05 |
NGC 2424 | R | 15L1 | 86.5 | e | 20.52 | 112.0 | 23.40 | 31.7 | 6.62 | 11.9 | 2.49 | |
IC 2207 | R | 10L2 | 86.5 | e | 21.04 | 56.0 | 17.67 | 38.2 | 12.06 | 6.7 | 2.12 | |
ESO 564-027 | r | 30E2 | 88.0 | sech | l | 21.11 | 140.4 | 18.33 | 50.9 | 6.65 | 6.6 | 0.86 |
ESO 436-034 | g | 60E3 | 88.0 | sech | l | 20.96 | 82.1 | 18.32 | 22.6 | 5.04 | 6.8 | 1.52 |
ESO 319-026 | g | 30E3 | 86.5 | sech | e | 21.51 | 60.1 | 13.21 | 14.2 | 3.12 | 3.1 | 0.68 |
ESO 319-026 | r | 30E2 | 86.0 | i | 21.14 | 63.0 | 13.84 | 14.7 | 3.23 | 3.2 | 0.70 | |
ESO 319-026 | r | 30E3 | 88.0 | sech2 | i | 21.94 | 61.9 | 13.60 | 13.2 | 2.90 | 2.8 | 0.62 |
ESO 319-026 | i | 30E3 | 88.0 | sech | l | 21.01 | 64.8 | 14.24 | 14.3 | 3.14 | 3.4 | 0.75 |
ESO 321-010 | g | 30E3 | 88.0 | sech | l | 20.11 | 64.4 | 12.34 | 19.9 | 3.81 | 4.7 | 0.90 |
ESO 321-010 | r | 30E3 | 88.0 | sech | l | 19.54 | 64.8 | 12.42 | 21.6 | 4.14 | 5.0 | 0.96 |
NGC 4835A | r | 40E3 | 85.5 | l | 20.92 | 90.0 | 18.56 | 40.4 | 8.33 | 7.9 | 1.63 | |
ESO 575-059 | r | 15E2 | 87.0 | sech2 | l | 20.78 | 60.5 | 17.53 | 22.6 | 6.55 | 6.4 | 1.86 |
ESO 578-025 | g | 30E2 | 86.5 | l | 21.58 | 50.7 | 20.63 | 17.0 | 6.92 | 7.5 | 3.05 | |
ESO 578-025 | g | 30E3 | 86.5 | l | 21.69 | 50.7 | 20.63 | 17.0 | 6.92 | 7.5 | 3.05 | |
ESO 578-025 | r | 30E2 | 86.0 | l | 21.02 | 50.4 | 20.51 | 14.9 | 6.06 | 7.4 | 3.01 | |
ESO 578-025 | i | 30E2 | 86.0 | l | 20.27 | 47.5 | 19.33 | 20.3 | 8.26 | 7.6 | 3.09 | |
ESO 446-018 | r | 30E2 | 86.5 | sech | l | 20.45 | 75.6 | 22.78 | 23.7 | 7.14 | 3.9 | 1.18 |
IC 4393 | r | 30E2 | 87.0 | l | 20.30 | 75.6 | 12.88 | 29.8 | 5.08 | 5.8 | 0.99 | |
ESO 581-006 | r | 30E3 | 86.5 | l | 21.33 | 55.8 | 11.07 | 18.2 | 3.61 | 5.6 | 1.11 | |
ESO 583-008 | r | 30E3 | 87.0 | i | 20.84 | 56.2 | 26.70 | 14.0 | 6.65 | 3.8 | 1.81 | |
UGC 10535 | r | 25E2 | 88.0 | sech2 | i | 21.50 | 41.4 | 20.57 | 11.4 | 5.66 | 4.3 | 2.14 |
NGC 6722 | r | 10E3 | 86.5 | sech | l | 19.47 | 86.0 | 30.77 | 21.6 | 7.73 | 6.2 | 2.22 |
ESO 461-006 | r | 60E3 | 87.5 | l | 20.96 | 61.6 | 23.37 | 21.1 | 8.01 | 3.8 | 1.44 | |
IC 4937 | g | 20E1 | 88.5 | sech2 | l | 22.73 | 74.9 | 10.41 | 93.8 | 13.04 | 5.6 | 0.78 |
IC 4937 | r | 30E3 | 88.0 | sech2 | l | 21.29 | 83.5 | 11.61 | 32.1 | 4.46 | 5.9 | 0.82 |
IC 4937 | i | 20E1 | 88.5 | sech | e | 19.92 | 78.1 | 10.86 | 38.6 | 5.37 | 7.2 | 1.00 |
ESO 528-017 | g | 30E3 | 86.5 | l | 21.98 | 57.6 | 22.54 | 20.5 | 8.02 | 3.5 | 1.37 | |
ESO 528-017 | r | 60E3 | 86.5 | l | 21.28 | 55.1 | 21.56 | 19.9 | 7.79 | 2.7 | 1.06 | |
ESO 528-017 | i | 30E3 | 86.5 | sech | l | 20.89 | 50.0 | 19.57 | 22.6 | 8.85 | 3.2 | 1.25 |
ESO 187-008 | r | 30E3 | 85.5 | l | 20.77 | 50.0 | 13.65 | 15.1 | 4.12 | 4.4 | 1.20 | |
ESO 466-001 | i | 40E3 | 87.0 | e | 19.52 | 52.6 | 23.76 | 13.0 | 5.87 | 8.2 | 3.70 | |
ESO 189-012 | g | 60E3 | 87.0 | l | 21.50 | 56.2 | 29.81 | 26.8 | 14.21 | 3.8 | 2.02 | |
ESO 189-012 | r | 30E3 | 86.5 | l | 20.71 | 56.5 | 29.97 | 20.6 | 10.93 | 3.4 | 1.80 | |
ESO 189-012 | i | 20E3 | 87.0 | sech | l | 20.44 | 54.7 | 29.01 | 22.7 | 12.04 | 3.2 | 1.70 |
ESO 533-004 | r | 20E1 | 88.0 | l | 20.32 | 68.4 | 11.12 | 33.1 | 5.38 | 7.7 | 1.25 | |
IC 5199 | g | 30E3 | 86.5 | l | 21.55 | 64.1 | 20.44 | 19.9 | 6.35 | 5.6 | 1.79 | |
IC 5199 | i | 30E3 | 86.5 | l | 19.97 | 57.6 | 18.37 | 19.9 | 6.35 | 4.9 | 1.56 | |
ESO 604-006 | r | 30E3 | 90.0 | sech | l | 21.22 | 70.6 | 34.54 | 27.9 | 13.65 | 3.8 | 1.86 |
The different fitting methods were independently developed within two diploma
theses (Lütticke 1996; Schwarzkopf 1996).
The quality of the data basis for each project was the same. From the sample
presented here there were five objects in common. These are used to compare
the two methods and determine the quantitative difference of the derived
parameters.
galaxy | f(z) | i | z | h | |
[ ] | [ ] | [ ] | [ ] | ||
ESO 321-010 r { | sech | 88.0 | 5.0 | 21.6 | 64.8 |
88.0 | 5.2 | 25.9 | 64.1 | ||
ESO 461-006 r { | 87.5 | 3.8 | 21.1 | 61.6 | |
87.5 | 3.6 | 16.2 | 61.2 | ||
IC 4937 r { | 86.5 | 6.1 | 26.8 | 83.5 | |
sech | 89.5 | 7.0 | 27.4 | 75.6 | |
ESO 189-012 r { | 86.5 | 3.2 | 19.3 | 56.2 | |
88.0 | 3.6 | 21.3 | 56.9 | ||
ESO 604-006 r { | sech | 88.5 | 4.1 | 39.8 | 67.7 |
sech2 | 89.5 | 3.2 | 24.5 | 73.4 |
A subsequent analysis shows that it is not possible to ascribe the sometimes quite large discrepancies to the quality of the individual method. It turns out, that the main problem is the non-uniform determination of the fitting area. The intrinsic asymmetric variations of a real galactic disk compared to the model enforce a more subjective restriction of the galaxy image to the fitting region, whereby for example one has to exclude the bulge area and the dust lane.
This finding is in agreement with the study of Knapen & van der Kruit ([1991]) who compared published values of the scalelength and find an average value of 23% for the discrepancy between different sources. As already mentioned by Schombert & Bothun ([1987]) the limiting factor for accuracy of the decomposition is not the typical S/N from the CCD-telescope combination nor the errors in the determination of the sky background, but the deviation of real galaxies from the standard model.
In our former study (Paper I) with an earlier method to adapt Eq. (4), 20 of our 45 galaxy images have already been used. We decided to re-use them in this study to get models for as many galaxies as possible in a homogeneous way. Additionally, Paper I only presents the best fit values for the isothermal model, and uses a different definition of the cut-off radius.
Only three galaxies are in common with the sample of de Grijs ([1998]): ESO 564-027, ESO 321-010, and ESO 446-018. The mean difference for the scalelength is 10.4% (ranging from ) and for the scale height (normalised to the isothermal case) 4.0% . For the remaining galaxies there are no models in the literature.
Our model only represents a rather simple axisymmetric three-dimensional model for a galactic disk, consisting of an one component radial exponential disk with three different laws for the density distribution in the z-direction and a sharp outer truncation. Therefore it does not include additional components, such as bulges, bars, thick disks, or rings, and cannot deal with any asymmetries. Features like spiral structure or warps are not included, whereas Reshetnikov & Combes ([1998]) multiply their exponential disk by a spiral function introducing an expression to characterise an intrinsic warp depending on the position angle outside of a critical radius.
The choice of our fitting area tries to avoid the dust lane, possible only for almost edge-on galaxies, as a first step to account for the dust influence (cf. Sect. 4.4.2). Examples of models including a radiative transfer with an extinction coefficient can be found in Xilouris et al. ([1999]). However, introducing more and more new components and features automatically increases the amount of free parameters. Therefore we restricted our model to the described six parameters, to obtain statistically meaningful characteristics for galactic disks.
In the following we demonstrate that a simple disk model omitting the bulge component and the dust lane give indeed reasonable parameters.
We have studied the influence of the bulge for some of our objects including
the earliest type galaxy in our sample (ESO 575-059) presented here.
We have subtracted our derived disk model from the galaxy and then tried to
find the best representation for the remaining bulge by a de Vaucouleurs
r1/4 or an exponential model.
Taking the slope of the vertical profile at
and a fixed axis ratio,
we have constructed the 2-dimensional model of the bulge.
In agreement with Andredakis et al. ([1995]) we find, that bulges of
early type galaxies are better fitted by an exponential profile than by a
r1/4.
Figure 3 shows the resulting vertical and radial cuts for
ESO 575-059 together with the models.
Figure 3: Minor and major axes profile for ESO575059r15E2 (solid lines) together with the best disk model (dotted), the best bulge model (dashed), and the resulting combined profile (dotted-dashed) |
Therefore we conclude, that it is possible to nearly avoid any influence of the central component by fitting outside the clearly visible bulge region.
We checked the influence of the dust on our determined parameter
set, by studying simulated galaxy images with three different dust
distributions.
These dusty galaxies were kindly provided by Simone Bianchi
who calculated images with known input parameters for the star and dust
distribution with his Monte Carlo radiative transfer method (Bianchi et al.
[1996]).
We have defined a worst, best, and transparent case, according to
the dust distributions presented by Xilouris et al. ([1999]), of our
mean stellar disk (
,
,
and f1(z)).
The worst case is calculated with
,
,
and
,
the best case with
,
,
and
,
and a transparent case
without dust.
To be comparable we used the same method for selecting the fitting region by
masking the "visible'' dust lane and reserving a typical area for a possible
bulge component using mean values for the transparent case.
In contrast to our standard procedure we do not restrict the inclination
range from the appearance of the dust lane, but specify the best SQ model
in the range
in Table 4. For the models
marked with a "'' we pretend the correct input inclination.
Table 4 demonstrates, that even for the worst case we are able to
reproduce the input parameters within the range of the typical 20% error
discussed in Sect. 4.2. It should be mentioned that for each case
we overestimate the input scalelength and -height, whereas the determination
of
does not depend on the dust distribution. The implication on the
distribution of the ratio
will be discussed in a forthcoming
paper.
case | fn | ||||||
[ ] | [ ] | [ ] | [ ] | [%] | [%] | ||
trans. | 87.5 | 1 | 87.5 | 8.5 | 34.0 | 1.2 | 0.9 |
best | 85.0 | 1 | 84.5 | 8.5 | 35.0 | 1.2 | 3.9 |
best | 87.5 | 1 | 86.5 | 8.6 | 35.5 | 2.4 | 5.3 |
best | 90.0 | 1 | 87.0 | 8.7 | 36.0 | 3.6 | 6.8 |
best | 90.0 | 1 | 90.0 | 8.8 | 35.0 | 4.8 | 3.9 |
worst | 85.0 | 2 | 85.0 | 8.9 | 39.9 | 6.0 | 18.4 |
worst | 87.5 | 2 | 86.0 | 8.8 | 39.1 | 4.8 | 16.0 |
worst | 90.0 | 2 | 88.0 | 9.7 | 39.4 | 15.5 | 16.9 |
worst | 90.0 | 2 | 90.0 | 9.8 | 39.1 | 16.7 | 16.0 |
Trying to adapt a simple, perfect, and exact symmetric model to real galaxies always implies a compromise between the degree of any deviation and the final model (Sect. 4.6). The following list will provide some typical caveats found during the fit procedure which will characterise the quality of the specified model for individual galaxies.
ESO 112-004: warped, asymmetric, central part slightly tilted compared to disk, after fitting still remaining residuals.
ESO 150-014: slightly warped, minor flatfield problems.
NGC 585: remaining residuals.
ESO 244-048: possible two component system, slope of inner radial profile significantly higher than of an outer one, final model fits the inner parts.
NGC 973: one side disturbed by stray light of nearby star, seems to be radially asymmetric, remaining residuals.
UGC 3425: superimposed star on one edge.
NGC 2424: model does not fit very well without obvious reason.
ESO 436-034: strong bulge component, possibly barred, hard to pinpoint final model, remaining residuals.
ESO 319-026: outer parts show u-shaped behaviour, remaining residuals, therefore large ( ) difference in inclination angle.
ESO 321-010: u-shaped, no clear major axis visible, therefore uncertain rotation angle, bar visible, bulge rotated against disk.
NGC 4835A: strong residuals.
ESO 446-018: the different sides of the disk are asymmetric visible in radial profiles and on the contour plot.
IC 4393: similar to NGC 4835A.
ESO 581-006: galaxy shows typical late type profile, questionable, but nevertheless final model seems to fit well.
ESO 583-008: disturbed by superimposed star, shows warp feature and a bar structure, questionable, remaining residuals.
UGC 10535: one side slightly extended.
NGC 6722: only one side observed, bulge rotated against disk, barred, strongly disturbed by dust absorption, radial extension visible, therefore should be treated with caution.
ESO 461-006: minor flatfield problem seems to cause asymmetry, although model looks fine.
IC 4937: similar to NGC 6722, dominating bulge, small disk, model significantly different compared to the i and r image, model possibly hampered by strong dust lane.
ESO 578-025: bar visible.
ESO 466-001: maybe two components, final model represents only inner part, outer part clearly different from normal disk component.
ESO 189-012: slightly warped.
ESO 533-004: similar to NGC 4835A, model fits the whole galaxy, leaving more or less no bulge component.
IC 5199: slightly radial asymmetric.
ESO 604-006: only one side observed, bar structure visible.
The model limitations described above constrain the application of our fitting process. Therefore we had to exclude about 20 galaxies from our original sample. They all show significant deviations from the simple geometry and an inclusion of their parameters obtained by forcing the model to fit the data will spoil the resulting parameter distribution. One larger group classified mainly as S0 galaxies (e.g. NGC 2549, ESO 376-009, NGC 7332, ESO 383-085, ESO 506-033) shows a completely different behaviour of the luminosity distribution in the outer parts compared to the other galaxies. They all show an additional component, mainly characterised as an elliptical envelope. This is already visible in the contour plot, but becomes even more evident in a radial cut parallel to the major axis. In these cases the usual common curved decline of the profile (e.g. ESO 578-025) is missing, and is replaced by a more or less straight decline into the noise level, sometimes even by an upwards curved profile. Fitting these luminosity distribution by our one component exponential disk with cut-off, will therefore naturally provide parameters qualitatively different compared to late type disks. This will be discussed in detail in a forthcoming paper.
Another group consists of galaxies dominated mainly by their bulges, whereas the disk is only an underlying component, partly characterized as having thick boxy bulges (Dettmar & Lütticke [1999]), e.g. IC 4745, ESO 383-005, although there are also pure elliptical bulges (e.g. ESO 445-049, NGC 6948).
In the case of ESO 383-048 and ESO 510-074 the radial profiles clearly indicate that a more complex model will be needed to fit these kind of multicomponent galaxies. Galaxies like UGC 7170 or ESO 113-006 were excluded due to their strong warps, which made it impossible to fit the model in a consistent way. Mainly late type galaxies such as ESO 385-008, IC 4871, UGC 1281, or ESO 376-023 show a too patchy and asymmetric light distribution, that any attempt to fit the profiles will give only very crude, low quality parameters. UGC 11859 and UGC 12423 were rejected due to their thin faint disks, which will maybe overcome by taking new images with longer integration times to get a higher signal-to-noise ratio, whereas NGC 5193A is completely embedded into the surface brightness distribution of its near companion.
Acknowledgements
This work was supported by the Deutsche Forschungsgemeinschaft, DFG. 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 made use of the LEDA database (www-obs.univ-lyon1.fr). The authors wish to thank Simone Bianchi, who kindly provided us dusty-galaxy images produced with his radiative transfer code.
Copyright The European Southern Observatory (ESO)