A&A 404, 927-937 (2003)
DOI: 10.1051/0004-6361:20030524
L. Prisinzano 1 - G. Micela 2 - S. Sciortino 2 - F. Favata 3
1 - Dipartimento di Scienze Fisiche ed Astronomiche, Università
di Palermo, Piazza del Parlamento 1, 90134 Palermo, Italy
2 - INAF - Osservatorio Astronomico di Palermo, Piazza del Parlamento 1, 90134
Palermo, Italy
3 -
Astrophysics Division - Space Science Department of ESA,
ESTEC, Postbus 299, 2200 AG Noordwijk, The Netherlands
Received 6 May 2002 / Accepted 6 April 2003
Abstract
We present UBVRI photometry of the open cluster NGC 2422
(age
yr) down to a
limiting magnitude
.
These data are used
to derive the Luminosity
and Mass Functions and to study the cluster
spatial distribution. By considering the color-magnitude
diagram data and adopting a representative cluster main sequence,
we obtained
a list of candidate cluster members based on a photometric criterion.
Using a reference field region and an iterative procedure,
a correction for contaminating field stars has been derived in order
to obtain the Luminosity and the
Mass Functions in the
range.
By fitting the
spatial distribution, we infer that
a non-negligible number of cluster stars lies outside our investigated region.
We have
estimated a correction to the Mass Function of the cluster in order to take
into account the "missing'' cluster stars. The Present Day Mass
Function of NGC 2422 can be
represented by a power-law of index
(rms) - the Salpeter Mass Function in this
notation has index
- in the mass range
.
The index
and the total mass of the cluster
are very similar to those of the Pleiades.
Key words: open clusters and associations: individual: NGC 2422 - stars: luminosity function, mass function
Accurate determinations of the stellar Initial Mass Function together with the star formation rate are fundamental to understand star formation mechanisms and related astrophysical problems. Since Salpeter's estimation of the IMF for stars in the solar-neighborhood (Salpeter 1955), several investigations in Galactic and extra-galactic stellar systems seem to converge to a universal IMF described by a broken power-law (Scalo 1998). Many models of stellar population, chemical evolution and galactic evolution adopt a priori a single IMF assuming its universality, although their results are highly sensitive to uncertainties in the IMF (Kennicutt 1998).
Testing the universality of IMF is a challenge for astrophysicists because several indirect evidences suggest that the IMF ought to systematically vary with the time due to the different star forming conditions (Larson 1998). Nevertheless, no convincing proofs for a variable IMF still exist and evidences for uniformity of the IMF are deduced by estimates for different populations (Kroupa 2002; Scalo 1998,1986).
However, a large scatter in the logarithmic
power-law index, for stars more massive than ,
is evident. In order to
understand how large apparent IMF variations are due to uncertainties
inherent to any observational estimate of the IMF,
Kroupa et al. (2001) investigated the scatter, introduced by Poisson
noise and dynamical evolution of star clusters,
of the power-law indices inferred for
N-body model populations.
The resultant apparent variation of the IMF defines a "fundamental limit'' such that any true
variation in the IMF that is smaller than this fundamental limit is not
detectable. In addition,
determinations of the power-law indices are subject to systematic errors arising
from unresolved binaries.
Being systems of coeval and equidistant stars with the same chemical composition, open clusters are key samples in investigating the IMF and its possible spatial and temporal variations. Determination of the IMF of open clusters can however be challenging because of the contamination from background Galactic field stars. A further complication comes from the difficulty to transform the observed Present Day Mass Function (PDMF) into the IMF using proper assumptions on the stellar and dynamical evolution, mainly affecting high and low-mass stars, respectively.
In order to reduce these complications it is convenient to study young open clusters or star forming regions. In the low mass star range the IMF is more uncertain but stellar evolution effects are not important; in this mass range the Present Day Mass Function is representative of the IMF.
Additionally, for the question of the universality or variability of the IMF, it is convenient to compare clusters with same age in order to highlight other possible parametric dependence. For these reasons we choose to study the southern Galactic open cluster NGC 2422 which has an age comparable to the well studied Pleiades cluster.
The equatorial (J2000.0) and galactic coordinates of NGC 2422 are
RA
,
Dec = -14
30
(Lynga 1987),
l= 230
.97, b= 3
.13, respectively.
Estimates of the fundamental parameters of this cluster, such as age,
distance modulus
and reddening were given in the past
by various authors as summarized by Barbera et al. (2002).
The age of 108 yr was estimated using theoretical isochrones
(Rojo Arellano et al. 1997); the most recent value of the
distance
d=498+135-88 pc, corresponding to a distance modulus
(M-m)0=8.48+0.52-0.42,
was deduced by Hipparcos measurements on 4 stars
(Robichon et al. 1999) while the most recent
value of the reddening
was reported by
Dambis (1999).
Several papers have been devoted to NGC 2422 in the past:
the most recent of them report Strömgren photometry
(Shobbrook 1984; Rojo Arellano et al. 1997; Nissen 1988), while
the available
UBV photometric values, either photographic and photoelectric, extend only
down to
(Zug 1933; Hoag et al. 1961; Smyth & Nandy 1962; van Schewick 1966; Lynga 1959; Ishmukhamedov 1967).
Mean photometric data and spectral classification from the former papers
were compiled by Mermilliod (1986) in a catalog of 212 objects.
Using this catalog and 564 additional stars from Ishmukhamedov (1967),
Barbera et al. (2002) obtained a large literature-based compilation of measured data
of stars in the field of NGC 2422;
for some of these stars, X-ray
counterparts were found.
The layout of our paper is the following. We present, in Sect. 2, our photometry and astrometry and in Sect. 3 the method adopted to select the candidate cluster member sample. In Sect. 4 we describe how the Luminosity and Mass Functions of NGC 2422 and the spatial distribution of cluster members were obtained. Finally in Sect. 5 we summarize and discuss our results.
The data used in this paper are CCD images in the UBVRI pass-bands
collected at the 0.9-m telescope of the Cerro
Tololo Inter-American Observatory (CTIO) on January 30, 1997. The scale on the
sky of the instrument is 2.028 arcsec/pixel, for a total field of view
of
square degrees, making such an instrument suitable to cover most of the
apparent size of the cluster. The
observing log is summarized in Table 1.
While the quality of seeing was limited, all
images were collected
in photometric conditions. Possible effects
on the crowded field photometry, due to the limited
seeing, have been ruled out by the artificial
star test described in Sect. 2.4.
Table 1: Log-book of the NGC 2422 CCD observations.
All the images were pre-processed in a standard way with IRAF, using the sets of bias and sky flat-field images collected during the observing night. The instrumental magnitudes and the positions of the stars for each frame were derived by profile-fitting photometry with the package DAOPHOT II and ALLSTAR (Stetson 1987). Then we used ALLFRAME (Stetson 1994) to obtain the simultaneous PSF-fitting photometry of all the individual frames. In order to obtain the transformation equations relating the instrumental magnitudes to the standard UBV (Johnson), RI (Kron-Cousins) system, we also derived the instrumental profile-fitting photometry for the two Landolt (1992) fields of standard stars SA95 and SA98 observed during the same night.
In order to obtain the total integrated instrumental magnitudes we derived aperture photometry for the same stars that we used to define the PSF, after the digital subtraction of neighboring objects from the frames. We used DAOGROW (Stetson 1990) to obtain the aperture growth curves for each frame and to compute the aperture correction to the profile-fitting photometry.
The transformation coefficients to the standard system
were derived using
transformation equations of the form:
![]() |
|||
![]() |
|||
![]() |
(1) | ||
![]() |
|||
![]() |
Due to the limited number of standard star observations, we do not have a complete coverage of airmass and we have not been able to derive the extinction coefficients. In order to overcome this problem, we used typical values for CTIO available at http://www.noao.edu/scope/ccdtime/ctio.db; both the adopted extinction coefficients and the best fit values for the zero points and the color coefficients are summarized in Table 2. Second order color terms were tried and turned out to be negligible in comparison to their uncertainties.
Table 2: Coefficients of the calibration equations.
The calibration process we have adopted is based on two
sequential steps, following a procedure adopted by Stetson (private
communication). First, 45 program stars on the NGC 2422 images
were selected using the condition (following a criterion described in
Stetson 1993, Sect. 4.1) that each star has to be well
separated from its neighbors, observed in all frames, and with a
statistic index ,
relative to the goodness of the PSF-fitting
photometry, less than 1.5. The photometric calibration based on the
Landolt standards is applied to these selected stars only (which we
refer to as local standards), which were then used (together with
the Landolt standard stars) to calibrate the other program stars.
The astrometric solution has been computed using
as reference the recently released Guide Star catalog, version 2.2.01
(GSC 2.2).
At beginning, our pixel coordinate list was matched to the
celestial coordinate list of the GSC 2.2
by projecting the celestial coordinates onto
a plane and using as reference three stars for which
we had both the pixel and the celestial coordinates from the Hipparcos
catalog (Turon et al. 1993).
An initial estimate for the linear transformation has been computed using the reference coordinates of the 1350 matched stars. Then, a plate solution has been computed using the same matched pixel and celestial coordinates by fitting a power series polynomial (IRAF task CCMAP). The final accuracy is of 0.24 arcsec. Finally, the IRAF tasks CCSETWCS and SKYPIX were used to obtain the celestial coordinates of the total sample.
In order to obtain the color-magnitude diagram (CMD) of the cluster, a selection based on the sharp parameter was first done following Stetson (1987). The sharp parameter is related to the angular size of the astronomical object allowing to reject non-stellar objects as semi-resolved galaxies and unrecognized blended double stars, for which sharp is significantly greater than zero, or cosmic rays, for which sharp is significantly less than zero. We considered stellar objects those having the sharp parameter in the (-0.2-+0.2) range. In this way we obtained a list of 36 101 stellar objects for which we have Vand I magnitudes. For 35 732 of these we have VRI magnitudes, for 35 140 we have BVRI magnitudes and finally, for 33 494 we have UBVRI magnitudes.
The V vs. V-I color-magnitude diagram for all the stellar objects in the NGC 2422 field is shown in Fig. 1 (left upper panel). Horizontal bars indicate the median errors in color, while vertical bars indicate the median errors in magnitude for bins of one magnitude. Clearly, the diagram is heavily contaminated by background and foreground stars, as expected from the position of the cluster in the Galactic disk.
In order to determine the cluster main sequence minimizing
the contamination effects, we show
in Fig. 1 (right upper panel) the V vs. V-I color-magnitude diagram
approximately corresponding to the cluster core.
The centroid of this circular
region of 16 arcmin of radius
has been estimated using the stars with
.
We considered several theoretical isochrones available in literature
(Castellani et al. 1999; Girardi et al. 2000; Siess et al. 2000; Baraffe et al. 1998)
and we have found that the cluster main sequence is well fitted by
the 100 Myr and solar metallicity (Z = 0.02)
theoretical isochrone computed by
Siess et al. (2000, SDF00), for higher mass stars (
),
and the one computed
by Baraffe et al. (1998, BCAH98) with a general mixing length parameter
,
for lower mass stars (
).
Both isochrones were transformed into the V vs. V-I plane using
the apparent distance modulus
(m - M)V = 8.76, derived from the distance
modulus
(m-M)0 = 8.48 of Robichon et al. (1999) and the
interstellar absorption
AV= 0.273, from the standard relation
(Mathis 1990) where the reddening
is the value calculated by Dambis (1999).
The color excesses
,
,
were derived using the relations for
RV=3.1 given by Munari & Carraro (1996).
As shown in Fig. 1 (left lower panel), the resulting theoretical
isochrone is in good agreement with the apparent cluster main sequence,
confirming the literature cluster parameters, such as distance
and reddening.
Furthermore, using a set of theoretical isochrones by Siess et al. (2000) for
different ages, and the well defined main sequence of bright stars in
the B vs. U-B color-magnitude diagram, we verified that the NGC 2422 age is 108 yr, in agreement with the value given by Rojo Arellano et al. (1997).
Finally, we matched our photometric data with the X-ray sources of Barbera et al. (2002). The matched stars are plotted with large squares in Fig. 1 (right lower panel). There are two populations of X-ray sources: a blue faint population, dominated by background objects, and a group of sources belonging to the cluster main sequence. The photometric position of the X-ray members allow us to confirm our choice of the main sequence.
In Fig. 2, we show the
V vs. V-I , I vs. R-I , V vs. B-V and B vs. U-B color-magnitude diagrams for the stars in the NGC 2422 field.
The curve in the V vs. V-I and in the I vs. R-I color-magnitude diagrams
is the Siess et al. (2000) theoretical isochrone
extended to lower stars using the Baraffe et al. (1998)
theoretical isochrone as described above.
For the other two diagrams we only considered stars with
and
,
respectively, and
the Siess et al. (2000) theoretical isochrone because we use the V vs. B-V and the B vs. U-B color-magnitude diagrams to select only the bright
cluster members.
![]() |
Figure 1:
The V vs. V-I color-magnitude diagram for the full catalog ( left upper panel)
and for stars within ![]() |
Open with DEXTER |
![]() |
Figure 2:
The V vs. V-I , I vs. R-I , V vs. B-V and B vs.
U-B color-magnitude diagrams for the stars in the NGC 2422 field.
Horizontal bars indicate the median errors in color, while vertical bars
indicate the median errors in magnitude for bins of one magnitude. The solid
lines are the theoretical isochrones described in the text. In the V vs. B-V and B vs.
U-B color-magnitude diagrams only stars with
![]() ![]() |
Open with DEXTER |
Since we are interested in deriving the Luminosity Function from star counts and because of the limited seeing of our data, particular attention has been devoted to estimate the accuracy of the photometry and the completeness of the derived star list. For this, a list of artificial stars was created and added to the original frames in order to compare the photometric results of the recovered artificial stars and the input values.
In order to avoid overcrowding, the artificial stars were placed in a spatial grid (cf. Piotto & Zoccali 1999) such that the separation of the centers in each star pair was four PSF radii plus 1 pixel. Using a random-number generator for the V magnitudes, a list of 2470 stars (7% of the total sample) with a flat distribution of the instrumental V magnitudes between 11 and 21.5 was created. These limits were chosen on the basis of the instrumental V vs. V-I color-magnitude diagram. The V magnitudes for each artificial star were converted to the other filters using the fiducial lines representing the V vs. V-I , V vs. B-V , I vs. R-I and B vs. U-B color-magnitude diagrams.
DAOPHOT's ADDSTAR routine was used to add these artificial stars into copies of the original data frames, with the appropriate frame-to-frame shifts in their position and brightness. Calibrated magnitudes were derived using the same photometric parameters and the same procedure described in Sect. 2.1.
We estimated the completeness fraction as the ratio between the number of artificial stars recovered simultaneously in the V, I and R filters and the number of added stars per one magnitude bin. This condition was imposed because the photometric selection of low-mass candidate cluster members was based on the use of V vs. V-I and I vs. R-I color-magnitude diagrams.
We found that this ratio is equal to 1 down to V=18.26,
while it decreases to
to the limit
of our data. Therefore,
we can
conclude that our catalog is complete down to V=18.26 over the whole
V vs. V-I and I vs. R-I color-magnitude diagrams while it
is 94% complete to V=19, so that
a
6% correction will be required to the number of
faintest stars when computing the Luminosity Function.
Finally, photometric errors were estimated by the differences between
the "observed'' magnitudes and colors derived for the artificial stars and their
known input values. We defined the external error
for the V magnitudes and
for the V-I and R-I colors as in Stetson & Harris (1988); the results
are summarized
in Table 3.
Table 3: External photometric errors from the artificial star photometry.
In order to select possible NGC 2422 cluster members, we used only photometric informations as described in the following steps:
The selection procedure described above was chosen to take advantage of all the available photometric informations and to avoid excluding too many possible cluster members.
The photometric/astrometric catalog of the candidate cluster members,
containing 1277 stars, is given in Table 4
where we report RA and Dec (J2000.0) coordinates in decimal
degrees, an identification number for each star, the
U, B, V, R, I magnitudes and the associated uncertainties.
Figure 3 shows these cluster candidates in the color-magnitude diagrams so far considered. Clearly, the sample includes contaminating objects that do not belong to the cluster. The approach to estimate and to deal with such a contamination will be described in the following section.
![]() |
Figure 3: The color-magnitude diagrams for the possible photometric candidate members of NGC 2422. The sample has been obtained as described in Sect. 3. |
Open with DEXTER |
To construct the Luminosity and Mass Functions of the cluster we need to correct the luminosity distribution of our selected sample for field star contamination. To quantify the contamination we adopted an iterative procedure schematically summarized in the following steps:
As summarized above, a first approximation to the contamination
of foreground and background
sources has been obtained using the lower region of our field where we expect
to find only few cluster stars. In fact, as shown in Fig. 4,
the cluster is mainly concentrated in the upper region of our field where we
defined the "central cluster region'' corresponding to the circular area within the radius
of 27.07 arcmin.
A value for the cluster centroid position has been calculated as
the median value of the X and Y coordinates of the bright stars
(
). This position corresponds to the equatorial coordinates (J2000.0)
RA
Dec
.
We verified that the coordinates of
the centroid do not change significantly if the V magnitude limit
changes in the
range.
![]() |
Figure 4:
A map of the brightest selected stars in the total field of view
(
![]() ![]() ![]() |
Open with DEXTER |
We consider the rectangular region in Fig. 4, of
size (2000
),
as representative of the Galactic field (the "field region''), i.e.
we suppose that the observed luminosity distribution
in the rectangular region is given by
![]() |
(3) |
The NGC 2422 Luminosity Function,
,
has hence been
obtained subtracting the field
star luminosity distribution,
,
from the luminosity distribution
of all the stars in our
sample within the "central cluster region'',
.
Figure 5 shows the contaminated luminosity
distribution of the candidate members
(dotted line), the luminosity distribution of the field stars
(dashed line) and the
Luminosity Function of NGC 2422 obtained as described above and
corrected for incompleteness as derived in Sect. 2.4
(solid line).
We found that the contamination rate from field stars
is negligible for the brightest stars, but
starts to become significant at
where it is of the order of the 40%. At fainter magnitudes this value increases with a peak of the
order of
90% around
where the contribution of the field giant stars is dominant.
The total number of the cluster members estimated at this step is 347.
![]() |
Figure 5: Luminosity distribution of the candidate members without correction for contamination (dotted line), the luminosity distribution of the field stars (dashed line) and the Luminosity Function of NGC 2422 (solid line) obtained as difference between the last two and corrected for incompleteness, as described in Sect. 2.4. |
Open with DEXTER |
As already mentioned, due to dynamical evolution, the angular size of the cluster is almost certainly larger than the investigated "central cluster region'' hence this Mass Function is not representative of the whole cluster mass distribution. In order to estimate the correction to the Mass Function due to the dynamical evolution, we studied the spatial distribution of NGC 2422 as described in the following sections.
![]() |
Figure 6:
First approximation of the NGC 2422 Present Day Mass Function
obtained in our central cluster region and without any correction.
![]() |
Open with DEXTER |
The spatial distribution of the cluster members
was studied in order to investigate
the cluster dynamical evolution and to determine the cluster dynamical
parameters. As first approximation, we corrected the radial surface density
distribution of candidate cluster members for foreground
and background contaminating stars assuming that the field
stars are uniformly distributed, i.e.
![]() |
(4) |
Table 4:
Results of King's empirical profile fitting to the radial surface
density distribution. Columns 1 and 2
give the V magnitude ranges and the corresponding mass
ranges; Col. 3 the fitted normalization constant, Col. 4 the assumed
tidal radius, Col. 5 the fitted core radius, Col. 6 the
of the
fitting. Finally, Col. 7 lists the number of counted cluster members
in our investigated
region for each magnitude subset, while Cols. 8 and 9 are the total number of
cluster members predicted by the fitted model
within the investigated radius (
3.7 parsec) and
the tidal radius, respectively.
The values are given for each of the four step of the adopted iterative procedure.
In order to take into account the dynamical evolution and mass segregation effects due to the
energy equipartition (see Kroupa et al. 2001 and relative references),
we subdivided the stars into 4 bins of V magnitude as in Table 5.
In Fig. 7 we present the cluster surface density profiles in
as a function
of radius in parsec, assuming the distance value of 497.5 pc
by Robichon et al. (1999).
![]() |
Figure 7: The radial distribution of stellar density in NGC 2422, for 4 different bins of magnitude. |
Open with DEXTER |
The tidal radius was calculated as in Jeffries et al. (2001)
as
,
where
is the cluster
mass in solar masses.
In order to estimate the total mass of the cluster we integrated
the Mass Function obtained in Sect. 4.1
and found a value of
,
corresponding to the
tidal radius value
0.34 pc.
The results of King's empirical profile
fitting, with the
1-
uncertainty estimates for the parameters,
are reported in Table 5
and are shown in Fig. 7.
Comparing
the
values with the number of degrees of freedom (
),
we conclude
that the fits are all acceptable.
From the resulting distribution we can see that, while the distribution
of the brightest stars
vanishes to zero as the radius increases
(at the top on the left of Fig. 7),
the other profiles tend to non negligible values when fainter
V magnitudes are considered.
This result suggests an evidence of mass segregation
of high mass cluster stars toward the center of the cluster and
low mass cluster stars
out of our investigation region of 27 arcmin of radius.
This conclusion is confirmed by the increase of the core radius
as the mass decreases.
In order to estimate the fraction of low-mass stars beyond
the "central cluster region'',
we extrapolated the surface density profiles as far as the
tidal radius. Using the integral of Eq. (5) as given by King (1962),
we calculated, for each bin of magnitude,
the total number of cluster stars within the tidal
radius, indicated in Table 5 by
.
We compared the calculated total number
with the number of cluster members
,
detected in our survey, that is within
pc
and we found that while the stars with
are all in the investigated region,
the fraction of the cluster
stars lying outside our survey is
8% for
,
22% for
and 38% for
.
In order to verify the consistency of these
numbers we also integrated King's empirical profile within the radius
,
indicated by
,
finding consistent values with the number of the cluster stars found in the
central cluster region (
).
The results of the King's empirical profile integration also suggest that
the cluster
star contribution within the adopted "field region'',
is non-negligible and that the
estimated field star luminosity and spatial distributions have to be corrected for the
cluster star contamination. In order to estimate the correction factor, we split the
"field region'' in a grid of
subregions and we calculated
the King's profile
surface density,
,
where ri is
the radius of the centers of each subregion i. Therefore, the total
number of cluster stars within the "field region'', per magnitude bin, is
given by
![]() |
(6) |
Finally, using Eq. (2), we estimated a more
accurate field star distribution as
Using the value of
from Eq. (7) and the iterative approach
described in Sect. 4, we recomputed the Luminosity Function
of the cluster taking into account the presence of cluster stars in the "field region''.
This final Luminosity Function has been converted into a Mass
Function that has been integrated to compute a new value for the
total mass (
), corresponding to
pc. Dynamical parameters were also obtained in the
various iterations, and by integration of King's profiles we
interpolated the ratio
,
deriving a factor 1/R for the Mass Function to correct for the
cluster stars lying beyond the investigated cluster region. The
corrected Mass Function has been computed by the expression:
![]() |
(8) |
![]() |
Figure 8:
Comparison of the NGC 2422 Mass Function
corrected for dynamical evolution effects (filled points)
with the NGC 2422 Mass
Function without correction for dynamical evolution (empty triangles)
and with the Pleiades one (empty circles). The line is the power law fitting
the MF in the
![]() ![]() |
Open with DEXTER |
According to the recently proposed analytical IMF forms (Kroupa 2002),
we considered a
multi-part power-law IMF to fit our data.
We found that in the mass range
the value of
is
(rms) and
for NGC 2422 and
Pleiades, respectively.
The latest value is very close to the value 2.67 derived by
Barrado y Navascués et al. (2001) for the Pleiades,
using the Hambly et al. (1999) data in the
mass range.
The two values here derived
are in better agreement, within the errors, with the Scalo (1998) value
rather than with the Salpeter (1955) value
.
Using the final Mass Function, we also estimated the corrected cluster
total mass (
)
and we compared this value
with the NGC 2422 total mass estimated within the "central cluster region'',
.
In the same
mass range (
), the
total mass of the Pleiades is
,
a value very similar to the one of NGC 2422.
We note that our estimate of the total mass of NGC 2422 is, in any case,
a lower limit, since we do not have taken into account the presence of
companions in photometric binaries.
Using UBVRI images, covering a field of view of
square degrees, we extracted a
photometric and astrometric catalog of the stars in the field of NGC 2422 down to
(
).
Adopting a representative cluster main sequence to the
color-magnitude diagrams, a photometric criterion was defined to obtain
a list of candidate cluster members.
A test with artificial stars allowed us to verify that our data are
complete down to
V=18.26 while a correction of
6% was necessary for the lowest
V range
(
)
in deriving the Luminosity Function.
We have defined a "central cluster
region'' within a radius of 27.07 arcmin from the cluster centroid and
a "field region'' of
square degrees
to estimate the contamination of background and foreground stars.
An initial determination of the Luminosity and Mass Functions has been obtained
assuming that all the cluster stars lie in
the "central cluster region''. By applying an iterative procedure we have
estimated the number of cluster stars
within the "field region'' in order to obtain a more accurate correction
of the cluster
Luminosity Function for field star contamination.
Evidence for mass segregation and energy equipartition have been found from the
spatial distribution of the stars. Extrapolation of the King (1962)
empirical model
allows us to infer that, while all the stars with
are inside the "central cluster region'',
a non negligible fraction of lower mass stars lies outside.
In particular, by integrating the spatial distribution within the tidal radius
we are able to estimate the number of cluster stars lying beyond the "central
cluster region''
and within the tidal
radius. Therefore, we are able to estimate a correction to the Mass Function
in order to take into account the cluster dynamical evolution.
The corrected Present Day Mass Function was compared with that of the
Pleiades which is known down to
the brown dwarf limit (
).
We found that, in the mass range
,
the Mass Function of NGC 2422
can be represented by a power law of index
(rms)
comparable with the index
obtained
for the Pleiades in the
same mass range. The index is also consistent with
the data presented in the log M vs.
plot of Kroupa (2002).
By taking into account the correction due to the dynamical evolution, we computed a lower limit of
the total mass of NGC 2422
as
.
This value
is similar to that of the Pleiades.
A future spectroscopic study and a deep survey on a wider field of this cluster will allow us to find an independent membership criterion and to extend the Mass Function below the mass limit of the present survey.
Acknowledgements
We acknowledge P.B. Stetson for having made available to us his software and for useful discussions. This work has been partially supported by MIUR.