L. Verdes-Montenegro 1 - J. Sulentic 2 - U. Lisenfeld 1,3 - S. Leon 1 - D. Espada 1 - E. Garcia 1 - J. Sabater 1 - S. Verley 1,4
1 - Instituto de Astrofísica de Andalucía, CSIC, Apdo. 3004, 18080 Granada, Spain
2 -
Department of Astronomy, Univ. of Alabama, Tuscaloosa, USA
3 -
Dept. Física Teórica y del Cosmos, Universidad de Granada,
Spain
4 -
LERMA - Observatoire de Paris, 61 Av. de l'Observatoire, 75014 Paris, France
Received 29 October 2004 / Accepted 16 February 2005
Abstract
The AMIGA project (Analysis of the Interstellar Medium of
Isolated Galaxies) is compiling a multiwavelength database of isolated
galaxies that includes optical (B and H
), infrared (FIR and
NIR) and radio (continuum plus HI and CO lines) properties. It involves a
refinement of the pioneering Catalog of Isolated Galaxies.
This paper is the first in a series and begins
with analysis of the global properties of the nearly redshift-complete
CIG with emphasis on the Optical Luminosity Function (OLF) which we
compare with other recent estimates of the OLF for a variety of
environments. The CIG redshift distribution for n= 956 galaxies
re-enforces the evidence for a bimodal structure seen earlier in smaller
samples. The peaks at redshift near 1500 and 6000 km s-1, corresponding
respectively to galaxies in the local supercluster and those in more
distant large-scale components (particularly Perseus-Pisces). The two
peaks in the redshift distribution are superimposed on 50% or more of
the sample that is distributed in a much more homogeneous way. The CIG
probably represents the most homogeneous local field example that has
ever been compiled. Our derivation of the CIG OLF is consistent with
other studies of the OLF for lower density environments. This
comparison via the Schechter parameter formalization shows that: 1) M* increases with galaxy surface density on the sky and 2)
shows a weaker tendency to do the same. The CIG represents the largest and most complete foundation for studies of isolated galaxies and is likely
as close as we can come to a field sample.
Key words: galaxies: evolution - galaxies: interactions - galaxies: luminosity function, mass function - surveys
The evolutionary history of galaxies is thought to be strongly conditioned by the environment. Evidence has emerged for interaction-induced emission enhancements (e.g. Sulentic 1976; Larson & Tinsley 1978; Joseph & Wright 1985; Bushouse 1987; Xu & Sulentic 1991, hereafter XS91) and interaction-driven secular evolutionary effects (e.g. Moore et al. 1996; Verdes-Montenegro et al. 2001) in galaxies that are members of binaries and dense groups. The observational evidence is sometimes weak or unclear. Part of the difficulty lies in the confusion between the roles of one-on-one interactions vs. more general correlations with average galaxy environmental density. Many of the uncertainties, both of the amplitude of enhancements and the connection between environment and parameters, reflect a lack of suitable control samples to which interacting sample properties can be compared. Ideally this would involve samples of isolated galaxies. Samples of isolated pairs and compact groups provide the parameters to quantify effects of 2 and n body interactions. Isolated galaxy samples should provide the baseline for interpreting the strength and properties of interaction-induced effects. Awareness of these two effects (one-on-one vs. local overdensity) is important in compiling isolated galaxy samples.
The most common reference or control samples found in the literature can be described as either "field'' or "normal''. The former refer to the most isolated galaxies while the latter refer to galaxies which show none of the generally accepted signs of interaction-induced activity. A field sample (e.g. Kennicutt & Kent 1983) might include any galaxy not belonging to a cluster, so galaxies in pairs, triplets and loose/compact groups would not necessarily be excluded. Normal galaxy samples would be defined in terms of specific parameters such as HI content (Boselli et al. 2001) or a specified level of nuclear activity. Study of a selected quantity as a function of the environment is then one way to quantify the level of environmentally induced activity.
The alternative approach involves sample selection using an isolation criterion. In the case of isolated binaries this likely reduces the interaction equation to the effects of one-on-one encounters. Studies of isolated galaxies usually involve from 10s to 100-200 objects (e.g. Huchra & Thuan 1977; Vettolani et al. 1986; Marquez & Moles 1999; Marquez et al. 2000; Colbert et al. 2001; Pisano et al. 2002; Varela et al. 2004). The largest samples of isolated galaxies in the literature involve, in most cases, monochromatic observations of subsamples from the Catalog of Isolated Galaxies (CIG: Karachentseva 1973, also referred as K73 in SIMBAD and KIG in NED databases; see Sect. 2) (Adams et al. 1980; Haynes & Giovanelli 1980; Sulentic 1989; Young et al. 1986; XS91; Hernandez-Toledo et al. 1999; Perea et al. 1997; Sauty et al. 2003).
Previous work suggests that small samples of isolated galaxies have limited statistical value. Ideally we seek a sample large enough to isolate a significant population of the most isolated galaxies. This motivated us to use the CIG as the basis for a large, well-defined and statistically significant multiwavelength database that can serve as a comparison template for the study of galaxies in denser environments. CIG galaxies were selected to be free of equal mass perturbers but hierarchical pairs and groups could not be removed without reducing the sample to negligible size. A large sample like CIG can be refined and quantified in terms of degree of isolation. It can then be correlated with multiwavelength interstellar medium (ISM) properties. The result can be a sample large enough to characterize the low density tail of the two-point correlation function. The result will also tell us: a) if truly isolated galaxies exist; b) in what numbers; and c) the environmental level where the onset of interaction-induced activity can first be detected. This study constitutes the AMIGA project (Analysis of the Interstellar Medium of Isolated Galaxies). AMIGA is compiling data that will characterise the all phases of the ISM and it is being released and periodically updated at http://www.iaa.csic.es/AMIGA.html.
This paper studies the optical properties of the entire CIG as the
foundation sample for AMIGA. These properties can be compared with
future refinements to monitor changes and possible biases in the
statistical properties. The CIG is now almost redshift-complete
allowing study of both apparent and distance-dependent properties. We
first consider the distribution of CIG galaxies in 2 and 3 dimensions.
This allows us to decide which local large scale-components contribute
most heavily to the sample and how close the CIG comes to showing
homogeneity. We then analyze sample completeness via the
test
followed by derivation of the optical luminosity function (OLF). This is a
much more reliable derivation than previous estimates because of the
near completeness of redshift measures. Finally, we compare the OLF
derived for the CIG sample with those of other samples representing different
environments.
Statistical studies of isolated galaxies require a large, preselected
and morphologically diverse sample. The tendency for galaxies to
aggregate in multiple systems and clusters at all scales might suggest
that such a sample is difficult to find. However Karachentseva (1973)
compiled the Catalog of Isolated Galaxies (CIG) which includes 1051
objects. All of the CIG objects are found in the Catalog of Galaxies
and Clusters of Galaxies (Zwicky et al. 1961-1968; CGCG) with
< 15.7 and
,
3% of the CGCG). Only one of the compiled objects is not a galaxy, but a globular cluster (CIG 781 = Palomar 15, Nilson 1973), so the size of the sample considered in the rest of this paper is n = 1050. We used the CIG as our starting point
because it has a number of strengths:
The CIG can be improved in several ways that take advantage of the digitized sky surveys (POSS1 and POSS2). Our two largest refinements include uniform reevaluation of morphology and isolation degree. This is being done for the entire sample except for the nearest dwarf galaxy subpopulation where numerous sources of distances and morphologies now exist.
The first papers discussing CIG morphology, and isolated galaxy
morphology in general, are roughly contemporary with the appearance of
the catalog in 1973. Galaxy classification data for CIG galaxies is
non-uniform and often contradictory. According to the NED and LEDA
databases the CIG is composed of
20% early types (E+S0), however the
distribution of the individual mophologies shows a large discrepancy,
as shown in Fig. 1.
We re-evaluated CIG morphologies using the POSS2 images (Sulentic et al.
in preparation) and find it possible to obtain reliable galaxy types
for 80% or more of the sample. The population of luminous
isolated spirals are the easiest to classify: near face-on spirals
could be easily recognized beyond 10 000 km s-1. The remaining 20% of the
sample are being supplemented with archival data (e.g. SDSS; James et al. 2004) or new CCD images on 1-2 m class telescopes. POSS2 provides the higher spatial resolution necessary to distinguish between basic subtypes.
![]() |
Figure 1: Histogram of the morphological types of the full CIG sample obtained from NED ( bottom) and LEDA ( top). |
| Open with DEXTER | |
Comparing CIG positions in the SIMBAD database and the Updated Zwicky Catalogue (UZC;
Falco et al. 2000) we found differences of up to several tens of arcsec
for some galaxies, large enough to make accurate telescope pointings or cross
correlations with on-line databases impossible. This motivated us to
systematically revise all of the CIG positions using SExtractor on the
images of the digitized sky surveys (Leon & Verdes-Montenegro 2003).
We found differences between old and new positions of up to 38.9
with a mean value of 2.4
for both SIMBAD and UZC.
The fraction of CIG with measured redshift has almost doubled in the
past 15 years. Our archival and bibliographic search reveals data for
almost the entire sample (956/1050 galaxies
). The redshift measures are given in
Table 1. About one half of the redshifts were found in NED
and we compiled the rest from 37 different sources. This includes 10 new HI observations from Nancay and Green Bank (see footnote to Table 1). Our search increased the CIG redshift sample by n = 489 relative to the recent studies by Hernandez-Toledo et al. (1999,
2001). Redshift distances were derived for all galaxies with V > 1000 km s-1 and are expressed as D =
where
is the velocity after the 3K correction (as given in Table 1) and assuming H0 = 75 km s-1 Mpc-1. 3K corrected velocities are
computed in the reference frame defined by the 3K cosmological background
radiation. They are corrected for local velocity inhomogeneities due
to the Local Group and Virgo Cluster. The velocity conversion is made
with the standard correction as defined in Courteau & Van den Bergh
(1999). The velocity and apex directions of the Sun relative to the
comoving frame have been derived from an analysis of the FIRAS data
(Fixsen et al. 1996) with
and
.
Redshift-independent distance
estimates and references are provided for galaxies with V < 1000 km s-1 (Table 2), except for CIG 506 (V = 998 km s-1), 657 (V = 626 km s-1), 711 (V = 976 km s-1),
748 (V = 961 km s-1) and 753 (V = 851 km s-1) listed in Table 1 since only redshift-derived distances could be obtained.
Table 1:
Recession velocities for the CIG sample (
> 1000 km s-1)1.
Table 2:
Distances for the CIG sample (
< 1000 km s-1).
Figure 2 shows the distribution of the CIG sample on the sky in
3000 km s-1 velocity intervals. This velocity segmentation makes it
easier to recognize concentrations associated with major components of
large-scale structure in the local Universe. The core of the Virgo cluster is indicated in the first segment with a circle of D = 12
.
Other Abell clusters in the same redshift range and with richness classes 1 or 2 are indicated with circles corresponding to
their core radius. As expected we see little correspondence between
the positions of the nearby cluster cores and CIG galaxies. Of course
some correspondence with more complex local large-scale structure
components has been found (Haynes & Giovanelli 1983). The 2-point
correlation function for the CIG (Vettolani et al. 1986) also shows
evidence for weak clustering.
![]() |
Figure 2:
Aitoff projection in right ascension and declination coordinates
showing the distribution on the sky of CIG galaxies in the following velocity ranges.
Galaxies in the 0-3000 km s-1 velocity interval are shown in a),
3000-6000 km s-1 in b), 6000-9000 km s-1 in c)
and 9000-12 000 km s-1 in d). The core of the Virgo cluster is
indicated by a circle with
|
| Open with DEXTER | |
Figure 3 shows the CIG redshift distribution which can be compared
with earlier studies involving smaller parts of the sample (Haynes &
Giovanelli 1983; XS91) when far fewer redshifts were available. The
earlier studies commented on evidence for large scale structure
components including the local (
1500 km s-1) and Pisces-Perseus
(
5-6000 km s-1) superclusters as well as the local void (
3000 km s-1) surrounding
the local supercluster. It is not clear that the latter void is actually seen. It is more appropriate to say that the level of the curve in the 3000 km s-1 region can be used to place an upper
limit on any quasi-homogeneous component of the CIG. It is clear that
the question of the existence of a galaxy "field'' component is out of
date. Large-scale structure dominates the distribution of galaxies and
one can only try to isolate the lowest density regions of this
structure. Figure 3 shows a comparison between the CIG redshift distribution and a corresponding homogeneous distribution of the same sample size with the same redshift distribution and Schechter luminosity function (see Sect. 4.3). It is clear that the fit is
not satisfactory for V < 6000 km s-1 due to the above mentioned
structure components. Removing these structures would provide an
estimate of the fraction of CIG galaxies that is homogeneously
distributed, at least in 2D. In order to estimate this number we have
assumed that the CIG is composed of both homogeneous and inhomogeneous
populations. The latter is dominated by peaks at 1500 km s-1 and 5000 km s-1. We can fit a homogeneous distribution to the population
underlying these peaks by scaling the solid curve downwards by a factor
of 0.6 (dashed curve). Thus about half of the CIG sample can be argued
to be reasonably homogeneous (see Figs. 2a,b). The solid curve fit to
the complete sample approaches homogeneity at
6500 km s-1 corresponding to a volume of a radius of about 90 Mpc). The higher velocity
part of the CIG samples a large enough volume to make sure that details
of individual large-scale structure components have little effect on
the velocity distribution. The residuals after subtraction of the
underlying homogeneous (dashed curve) component show two peaks
corresponding to the local and Pisces-Perseus Superclusters. Since the
total CIG comprises about 3% of the CGCG this means that 1-2% of
the CGCG can be argued to show homogeneity. This is about the same
population fraction as the dense isolated compact groups (Mendes de Oliveira & Hickson 1991, hereafter MH91; SR94) that lie at the other end of the "field'' clustering spectrum. This fractional similarity is probably reasonable because both CIG (densest regions and least dense
i.e. voids) and (e.g.) HCG share a similar avoidance of the most
clustered regions via an isolation selection criterion (Sulentic 1987).
The CIG is likely as close as we can hope to come to a local
homogeneous component of the galaxy distribution.
![]() |
Figure 3: Histogram of the optical heliocentric velocities of the 956 CIG galaxies with redshift data. Only CIG 402 is out of the plot, with V = 40 658 km s-1. The solid line corresponds to a homogeneous redshift distribution of the same sample size, velocity distribution and Schechter function. The dashed line has been obtained by scaling down the previous distribution by a factor of 0.6. |
| Open with DEXTER | |
We compiled blue magnitudes (mB) from the CGCG for all CIG galaxies and applied the following corrections.
![]() |
(1) |
| (2) |
Table 3: Optical magnitudes and luminosities of the CIG sample.
![]() |
Figure 4:
a) Zwicky-corrected magnitudes versus
|
| Open with DEXTER | |
Figure 5 shows the distribution of apparent
corrected magnitudes for the complete sample. We find only a few
galaxies (n = 19) brighter than
= 11.0. This weak tail
extending to
= 8.5 represents galaxies that are in a sense
interlopers to the CIG. They are the few large and bright
galaxies in the nearby Universe that escaped deletion by the isolation
criterion. They are almost certainly less isolated than the bulk of
the CIG.
![]() |
Figure 5: Distribution of the corrected apparent magnitudes (n = 1050). |
| Open with DEXTER | |
Table 4:
for different limiting magnitudes and different
absolute magnitude bins.
We used the
test (Schmidt 1968) to evaluate CIG sample
completeness. We calculate for each object the volume V contained
in a sphere whose radius is the distance to the object and the maximum
volume
contained in a sphere whose radius is the maximum
distance at which the galaxy would still be visible given the magnitude
limit of the CIG. We then calculate the average of the objects
brighter than the magnitude limit. This calculation is sensitive to
the choice of the brightest magnitudes included. We choose to neglect
galaxies brighter than 11 mag because they are obviously incomplete due
to their small number per magnitude bin (see Fig. 4). This was also
noted in the
value when including them. It was not necessary
to exclude galaxies fainter than 11 mag as they proved to be reasonably
complete. Figure 6 shows the cumulative
distribution as a function of limiting apparent magnitude. Results are
presented for the complete sample (n = 1031) and the subsample for
which recession velocities are available (n = 937), where galaxies
brighter than 11 mag have been excluded. The difference between the two samples is small and only visible at fainter magnitudes where the sample becomes slightly more complete. Our test suggests that the CIG is surprisingly complete (between 80-95%) brighter than 15.0. The
sample becomes rapidly less complete at fainter magnitudes. We
therefore adopted m = 15.0 (
= 0.41 for the sample with
recession velocities) as the cutoff for inclusion in the sample used to
derive the OLF. Hence the OLF has been calculated using n = 725 galaxies which corresponds to 734 galaxies with known distance and magnitudes in the range 11-15 mag minus 9 galaxies with very high or low luminosity excluded since they were scattered in bins containing a
low number of galaxies. In Table 4 we tabulate
as a function of both apparent and absolute magnitudes for
the redshift-complete sample. The high and low luminosity bins contain
few sources because the former are rare and the latter are restricted
to the extreme low redshift part of our sample. The depression centered
at about limiting magnitude
13.3 reflects the gap between our
local supercluster and large-scale structure components beyond.
![]() |
Figure 6:
The
|
| Open with DEXTER | |
The surprising level of completeness found for the CIG highlights its
strength as mentioned earlier. The completeness correction is done by
calculating
in intervals of 0.1 mag and adding the
necessary number of galaxies in order to obtain for
a value
of
0.5, characteristic of a complete homogeneous sample (see
e.g. Huchra & Sargent 1973). The galaxies added in each bin are taken
into account when calculating
for the next fainter bin.
The final result depends somewhat on the bin size because a change in
this parameter affects the assumed magnitude distribution of the added
galaxies. Here we choose to make the smallest bin size commensurate
with the precision of the data which corresponds roughly to the
expected error in the adopted apparent magnitudes.
Table 5 gives the results. We see that we have
added 716 galaxies in order to make the sample complete to
= 15.0. This yields a correction factor
= of 2.0. We can
estimate the uncertainty of this value due to the adopted bin size. If
we had chosen a bin size of 0.2 or 0.05 mag the number of
added galaxies would have changed to 532 and 851 respectively. These
correction factors would have been 1.7 and 2.2 respectively. We
estimate an uncertainty of about 15%.
Table 5: Correction for incompleteness.
Table 6: Optical luminosity function.
The differential OLF
estimating the number
of galaxies per unit volume and per unit absolute magnitude level,
is estimated from (Felten 1976)
![]() |
(4) |
![]() |
(5) |
![]() |
Figure 7:
The
|
| Open with DEXTER | |
![]() |
Figure 8:
Optical luminosity function for CIG galaxies for which
velocity information exists and with apparent magnitudes
between 11 and 15 (n = 725). The Schechter fit to this sample is
plotted as a solid line. The dashed line corresponds to a fit to the same sample when
galaxies with |
| Open with DEXTER | |
Table 7: Optical luminosity function for the CIG sample.
We fit the OLF with a Schechter function:
![]() |
(6) |
![]() |
Figure 9: All three OLF estimations for the CIG involving n=725 galaxies (our sample, filled triangles), 450 galaxies (CIG sample of XS91, open circles) and 294 galaxies (AIG sample of XS91, crosses). The CIG data are given in corrected magnitudes while the data of XS91 refer to uncorrected magnitudes. |
| Open with DEXTER | |
![]() |
Figure 10: The same as in Fig. 9 but calculated using uncorrected magnitudes for the CIG galaxies, as was done for the other two samples by XS91. |
| Open with DEXTER | |
Table 8: Optical luminosity function for the samples from the bibliography.
The
parameter shows less difference to previous estimates.
The underfit at the faint end results from the "interloper''
population of dwarf galaxies mentioned earlier. All of the galaxies in
the faintest two bins (n = 21) lie within a recession velocity
= 2000 km s-1 and 13 of them within
= 1000 km s-1. Only a very
small local volume is sensitive to such intrinsically faint galaxies.
However many were found in this volume because of the ineffectiveness
of the isolation criterion for galaxies within a few Mpc. We
recalculated the OLF without local galaxies using several velocity cuts
in addition to the restriction in magnitude range to 11-15. When we
remove galaxies with
<
km s-1, where
ranges
from 500 to 1500 km s-1,
changes systematically from -1.3 to -0.8. Since all except 3 of the n = 59 removed galaxies with
< 1500 km s-1 are fainter than
= -20 there is no change in the OLF for magnitudes brighter than -20.
However the fit to the
bright part of the OLF changes because we are truncating our reasonably complete sample, resulting in a failure of the Schechter function to model the OLF properly. The only interesting result to
emerge from this truncation involves the decrease in the
parameter. The fit obtained for
= 1500 is shown with a dashed line in Fig. 8 and the fit parameters are given in Table 7.
In this section we compare our CIG OLF with results of other samples
involving a range of environments (Table 8). We have
also included the results from the compilation of Binggeli et al. (1988)
for the OLF of field galaxies. We selected several
samples for comparison and concentrate on the shape of the OLF rather than
the actual space density (
)
which is an extremely difficult thing to
compare. The absolute magnitude range covered by each sample is given in
Table 8. All published values of M* have been
reduced to H0 = 75 km s-1 Mpc-1 and transformed to Zwicky
corrected magnitudes (
)
using the appropriate relation:
=
,
=
,
=
,
=
,
obtained from the relation given in Eq. (3) combined with the
transformations from Liske et al. (2003).
The comparison samples include:
![]() |
Figure 11: The Schechter function M* parameter as a function of the environment, ordered in an approximate way. The label "all'' indicates that the M* value has been obtained for all galaxies in the sample independently of the environment. |
| Open with DEXTER | |
Figure 11 plots M* for each sample ordered roughly by environmental density. The sequence indicates reasonably clearly the change in OLF parameters as one proceeds from higher to lower density samples. The former show an excess of high luminosity galaxies as inferred from M*. Our CIG value for M* betterfits somewhat denser environment than voids but with a lower density than some field estimates. This is consistent with the fact that the isolation criterion mitigates any possible bias towards inclusion of void members in the CIG. Void galaxies are often not isolated in projection but only in 3D. Even void samples contain interacting pairs (e.g. Grogin & Geller 2000) and M* will be affected depending on their fractional representation in the sample.
We also find a possible environmental trend (albeit with larger
scatter) for the
parameter (Fig. 12) in the sense
that it becomes more negative for denser environments. The location of
the CIG in this plot obviously depends on the inclusion or exclusion of
the local part of the sample (
< 1500 km s-1) which is dominated
by low luminosity dwarf galaxies (M> -19). Our results are consistent
with Marinoni et al. (1999) when local dwarfs are included in the CIG
which is reasonable since they sample galaxies with absolute
magnitudes down to MB = -15.2. Results for the HCG are
controversial since a significant dwarf population is only found if one
increases the diameter of the groups as defined for the high luminosity
members. While low luminosity CPG pairs are found, very few compact
groups composed entirely of low luminosity members are found, for example,
in the HCG (see SR94 Fig. 1).
![]() |
Figure 12:
The same as in Fig. 11 for the |
| Open with DEXTER | |
The CIG sample is the basis of the AMIGA project. It has many
advantages as a source of galaxies in low density environments, not the
least of which is its relatively large size. This means that it can be
refined without reducing the final sample population below a size that
would be statistically useful. We find that its 2D distribution is reasonably
homogeneous as we would expect for a distribution sampling, predominantly, the
peripheries of large-scale structure features. It is affected by the
local and Pisces-Perseus superclusters in 3D. The former because we are
inside it and the latter because it is rather large and diffuse.
Underlying these two bumps in the redshift distribution we again find
evidence that 50% or more of the sample shows a quasi-homogeneous
redshift distribution, motivating us to suggest that CIG is as close as
we can hope to come towards achieving a local "field'' population. A
test confirms the completeness of the CIG and a comparison of the OLF of the
CIG with that of other samples re-enforces the credibility of the idea that CIG OLF is
representative of the lower density parts of the galaxy environment.
Care must be taken with the local supercluster contribution to the CIG
because it samples the OLF to much lower luminosities than the rest
of the sample.
Acknowledgements
L.V.-M., U.L., D.E., S.L., S.V. and E.G. are partially supported by DGI (Spain) AYA 2002-03338 and Junta de Andalucía TIC-114 (Spain).