A&A 491, 767-780 (2008)
DOI: 10.1051/0004-6361:200810332
C. Bonatto - E. Bica
Universidade Federal do Rio Grande do Sul, Departamento de Astronomia CP 15051, RS, Porto Alegre 91501-970, Brazil
Received 5 June 2008 / Accepted 12 September 2008
Abstract
Context. Open clusters older than 4 Gyr are rare in the Galaxy. Affected by a series of mass-decreasing processes, the stellar content of most open clusters dissolves into the field on a time-scale shorter than
1 Gyr. In this sense, improving the statistics of old objects may provide constraints for a better understanding of the dynamical dissolution of open clusters.
Aims. Our main purpose is to investigate the nature of the globular cluster candidate FSR 1716, located at
and
.
We also derive parameters of the anti-centre open cluster Czernik 23 (FSR 834). Both objects have been detected as stellar overdensities in the Froebrich, Scholz & Raftery star cluster candidate catalogue.
Methods. The analyses are based on near-infrared colour-magnitude diagrams and stellar radial density profiles. The intrinsic colour-magnitude diagram morphology is enhanced by a field-star decontamination algorithm applied to the 2MASS J, H, and photometry.
Results. Isochrone fits indicate that FSR 1716 is more probably an old (7 Gyr) and absorbed (
0.2) open cluster, located
0.6 kpc inside the solar circle in a contaminated central field. However, we cannot rule out the possibility of a low-mass, loose globular cluster. Czernik 23 is shown to be an almost absorption-free open cluster,
5 Gyr old, located about 2.5 kpc towards the anti-centre. In both cases, solar and sub-solar (
)
metallicity isochrones represent equally well the stellar sequences. Both star clusters have a low mass content (
)
presently stored in stars. Their relatively small core and cluster radii are comparable to those of other open clusters of similar age. These structural parameters are probably a consequence of the several Gyr of mass loss due to stellar evolution, tidal interactions with the disk (and bulge in the case of FSR 1716), and possibly giant molecular clouds.
Conclusions. Czernik 23, and especially FSR 1716, are rare examples of extreme dynamical survivors. The identification of both as such represents an increase of
in the known population of open clusters older than
4 Gyr in the Galaxy.
Key words: Galaxy: open clusters and associations: general - Galaxy: open clusters and associations: individual: FSR 1716 - Galaxy: open clusters and associations: individual: Czernik 23
The Galaxy is an aggressive environment for star clusters in general, open clusters (OCs) in particular. These stellar systems continually undergo a series of dynamical processes such as mass loss associated with stellar evolution, mass segregation and evaporation, tidal interactions with the Galactic disk and bulge, and collisions with giant molecular clouds. Combined over time, such processes tend to accelerate the dynamical evolution, which produces significant changes in the cluster structure and eroded mass functions. Eventually, most OCs end up completely dissolved in the Galactic stellar field or as poorly-populated remnants (Pavani & Bica 2007, and references therein).
Theoretical (e.g. Spitzer 1958; Lamers & Gieles 2006), N-body (e.g. Baumgardt & Makino 2003; Goodwin & Bastian 2006; Khalisi et al. 2007), and observational (e.g. van den Bergh 1957; Oort 1958; von Hoerner 1958; Piskunov et al. 2007) evidence indicate that the disruption-time scale near the solar circle is shorter than 1 Gyr. Around this region, the disruption time-scale depends on mass as
(Lamers & Gieles 2006), which for clusters with mass in the range
corresponds to
.
In general, the effect of the relevant dynamical processes is stronger for the OCs more centrally located and the low-mass ones (see Bonatto & Bica 2007a, for a detailed discussion on disruption effects and time-scales). Indeed, OCs older than
1 Gyr are preferentially found near the solar circle and in the outer Galaxy (e.g. Friel 1995; Bonatto et al. 2006a), where the frequency of potentially damaging dynamical interactions with giant molecular clouds and the disk is lower (e.g. Salaris et al. 2004; Upgren et al. 1972). Disruption efficiency increases critically towards the Galactic
centre, to the point that the inner (
pc) tidal fields can dissolve a massive star cluster in
50 Myr (Portegies Zwart et al. 2002).
The above aspects considered, the natural expectation is that only a small fraction of the
OCs can reach old ages, and that the successful ones will be preferentially found at large
Galactocentric distances. In fact, of the 1000 OCs with known age listed in the
WEBDA
database, 180 are older than 1 Gyr, and only 18 (
)
are older than 4 Gyr (see also Ortolani et al. 2005a,b).
Not surprisingly, most of the OCs older than 1 Gyr identified so far are located outside the solar
circle (see, e.g., the spatial distribution of OCs of different ages in Fig. 1 of Bonatto & Bica 2007a).
In this context, it is naturally expected that the discovery and derivation of astrophysical parameters of old OCs will better define the OC parameter space. Thus, a better understanding of the dynamical survival rate of star clusters in the tidal field of the Galaxy can be reached. Such parameters, in turn, can be used in studies of star formation and evolution processes, dynamics of N-body systems, and the geometry of the Galaxy, among others.
In the present paper we study two stellar overdensities listed in the star cluster candidate catalogue of Froebrich et al. (2007a), which turn out to be very old star clusters. They are FSR 1716 and FSR 834. The present work employs near-IR J, H, and photometry obtained from the 2MASS
Point Source Catalogue (PSC). The spatial and photometric uniformity of 2MASS, which allow extraction of wide surrounding fields that provide high star-count statistics, are important to derive cluster parameters and probe the nature of stellar overdensities (e.g. Bica et al. 2008). For this purpose we have developed quantitative tools to statistically disentangle cluster evolutionary sequences from field stars in colour-magnitude diagrams (CMDs). We apply (i) field-star decontamination to quantify the statistical significance of the CMD morphology, which is important to derive reddening, age, and distance from the Sun; and (ii) colour-magnitude filters, which are essential for intrinsic stellar radial density profiles (RDPs), as well as luminosity and mass functions (MFs). In particular, field-star decontamination constrains the age and distance more, especially for low-latitude OCs (Bonatto et al. 2006a).
This paper is organised as follows. Section 2 contains basic properties and reviews literature data (when available) on both star cluster candidates. In Sect. 3 we present the 2MASS photometry and build the stellar surface-density distribution in the direction of both objects. In Sect. 4 we build CMDs, discuss the field-star decontamination algorithm, and provide tests of the old age of both clusters. In Sect. 5 we derive cluster fundamental parameters. Section 6 describes cluster structure by means of stellar RDPs. In Sect. 7 mass functions are built and cluster masses are estimated. In Sect. 8 aspects related to the structure and dynamical state of the present clusters are considered. Concluding remarks are given in Sect. 9.
The catalogue built by Froebrich et al. (2007a) includes 1021 star cluster candidates (hereafter FSR
objects) for Galactic latitudes
and all longitudes. The targets were selected by
an automated algorithm that identifies small-scale regions that present stellar overdensities,
applied to the 2MASS database. The overdensities are classified according to a quality flag,
``0'' and ``1'' for the most probable star clusters, and ``5'' and ``6'', which may be related to field fluctuations. Based on a combination of the number of cluster stars (corrected to a common
magnitude limit), the core radius and the central star density, Froebrich et al. (2007a) could discriminate
known globular clusters (GCs) from OCs. With this criterion applied to the overdensity catalogue,
they found 1012 OC and 9 GC candidates.
![]() |
Figure 1:
Left panel:
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
Table 1: General data on the clusters.
Several studies have explored the FSR catalogue with different approaches, with results that reflect the importance of such a catalogue. The recently discovered GCs FSR 1735 (Froebrich et al. 2007b) and FSR 1767 (Bonatto et al. 2007), and the probable GCs FSR 584 (Bica et al. 2007) and FSR 190 (Froebrich et al. 2008), are clear examples of the fundamental rôle played by the FSR catalogue in improving the statistics of very old star clusters. Indeed, FSR 1735 and FSR 1767 are the most recent additions to the Galactic GC population, currently of 160 members (e.g. Bonatto & Bica 2008a). On the other hand, the decontamination algorithm (Sect. 4) applied to the 2MASS photometry of some FSR GC candidates has shown that FSR 89 (Bonatto & Bica 2007a) and FSR 1603 (Bica & Bonatto 2008) are open clusters of age
1 Gyr, while FSR 1755 appears to be an embedded OC (Bica & Bonatto 2008).
During a close investigation of the overdensities, we noted that the CMDs of FSR 1716
and FSR 834 show features typical of old stellar systems. Indeed, FSR 1716 was classified as a GC candidate by Froebrich et al. (2007a). FSR 1716, at
and
is projected not far from the Galactic centre, which implies the presence of significant contamination by bulge stars. FSR 834, on the other hand, is a disk object projected towards the anti-centre (
and
). It has the poorly-studied OC Czernik 23 (hereafter Cz 23) as an optical counterpart. Czernik (1966) measured a diameter of 5
and estimated that Cz 23 contains about 40 member stars, and Ruprecht (1966) classified it as III 1p.
Figure 1 (left panel) shows a
2MASS
image of FSR 1716, where a significant concentration of stars is superimposed on a relatively crowded field, as expected from such a central direction. For Cz 23 we show in the right panel a
XDSS
R band image. In this case the cluster can be easily seen against a relatively clean field.
Table 1 provides fundamental data on both objects. Coordinates from Froebrich et al. (2007a) are given in Cols. 2 to 5; their quality flag in is Col. 6, and the core and tidal radii measured in the H band are given in Cols. 7 and 8. The age, central reddening, distance from the Sun, Galactocentric distance, and the components projected onto the Galactic plane derived in the present study (Sect. 5) are given in Cols. 9 to 16.
2MASS photometry in the J, H, and
bands was extracted in circular fields of extraction radius
centred on the coordinates of the objects (Table 1) by means of VizieR
. Previous studies with OCs in different environments (Sect. 1) revealed that in the absence of a neighbouring populous cluster and important differential absorption, wide extraction areas provide suitable statistics for a consistent characterisation of the field stars in terms of colour and magnitude. Based on this premise we adopted an extraction radius of
,
which is beyond the cluster radius (Sect. 6 and Col. 9 of Table 3) of the present objects. For decontamination purposes, comparison fields were extracted within wide rings located beyond the cluster radii. As a photometric quality constraint, the 2MASS extractions were restricted to stars (i) brighter than those of the 99.9% Point Source Catalogue completeness
limit
in the cluster
direction; and (ii) with errors in J, H, and
smaller than 0.3 mag. The 99.9% completeness limits refer to field stars, and depend on Galactic coordinates. In the present cases, the fraction of stars with J, H, and
uncertainties smaller than 0.06 mag is
.
A typical distribution of uncertainties as a function of magnitude, for objects projected towards the central parts of the Galaxy, can be found in Bonatto & Bica (2007b). Reddening transformations use the relations
,
,
,
and AJ=2.76
(Dutra et al. 2002), for a constant total-to-selective
absorption ratio
.
These ratios were derived from the extinction curve of Cardelli et al. (1989).
![]() |
Figure 2:
Stellar surface-density
![]() ![]() ![]() ![]() |
Open with DEXTER |
As a first step in understanding the nature of both overdensities, we show in Fig. 2 the
spatial distribution of the stellar surface-density as measured with 2MASS photometry. In both cases
we compute the surface density (,
in units of
)
in a rectangular
mesh with cells of dimensions
.
The mesh reaches total offsets of
with respect to the centre, in right ascension and
declination. Since the cluster radius of both objects is
(Table 3), most of the cluster structure should be contained in the central cell. Because of the significant contamination by field stars, the surface-density built with the observed (raw) photometry of the centrally projected OC FSR 1716 is very irregular (top-left panel), although an excess can be seen in the central cell, which corresponds to the overdensity detected by Froebrich et al. (2007a). Cz 23, on the other hand, clearly stands out in the central cell (top-right panel) against a less-irregular field. The bottom panels show the surface densities built with field-star decontaminated photometry (Sect. 4).
Especially in the case of FSR 1716, the stellar surface-density in the direction of both objects
clearly shows that field-star contamination should be taken into account. This assertion is further
supported by the
and
CMDs extracted from the
region of FSR 1716 (Fig. 3), and the
region of Cz 23 (Fig. 4). Features present in the central CMDs and in the respective comparison field (top and middle panels) show that field stars contribute in varying proportions to the CMDs, especially for the bulge-projected FSR 1716. Nevertheless, when contrasted with the equal-area comparison field extractions, sequences typical of old stellar systems are suggested in both cases - a conspicuous giant branch in FSR 1716 and a giant clump in Cz 23.
Field-star decontamination is a very important, yet difficult, step in the identification and characterisation of star clusters. Different approaches have been used for this, among them, those of Mercer et al. (2005) and Carraro et al. (2006). The first is based on spatial variations of the star-count density, but does not take into account CMD properties. In the latter, stars in a CMD extracted from an assumed cluster region are subtracted according to colour and magnitude similarity with the stars of an equal-area comparison field CMD.
In the present case, we apply the statistical algorithm described in Bonatto & Bica (2007b) to quantify the
field-star contamination in the CMDs. The algorithm makes use of both approaches above, in the sense that the relative star-count density together with colour/magnitude similarity between cluster and comparison field extractions are taken into account simultaneously. It measures the relative number densities of probable field and cluster stars in cubic CMD cells whose axes correspond to the J magnitude and the (J-H) and
colours. These are the 2MASS colours that provide the maximum variance among CMD sequences for OCs of different ages (e.g. Bonatto et al. 2004). The algorithm: (i) divides the full range of magnitude and colours covered by the CMD into a 3D grid, (ii) calculates the expected number density of field stars in each cell based on the number of comparison field stars with similar magnitude and colours as those in the cell; and (iii) subtracts the expected number of field stars from each cell. The algorithm is responsive to local variations of field-star contamination (Bonatto & Bica 2007b). Typical cell dimensions are
,
and
,
which are wide enough to
allow for sufficient star-count statistics in individual cells and narrow enough to preserve the
morphology of the CMD evolutionary sequences. The comparison fields are located within
(FSR 1716), and
(Cz 23). In
both cases, the inner boundary of the comparison field lies beyond the probable tidal radius
(Sect. 6), which minimises the probability of oversubtraction of member stars. We emphasise that the equal-area field extractions shown in the middle panels of Figs. 3 and 4 serve only for comparisons among the panels. Actually, the decontamination process is carried out with the wide surrounding area as described above. Further details on the algorithm, including discussions on subtraction efficiency and limitations, are given in Bonatto & Bica (2007b).
As extensively discussed in Bonatto & Bica (2007b), differential reddening between cluster and field stars
is critical for the decontamination algorithm. Large gradients would require wide cell sizes
or, in extreme cases, would preclude application of the algorithm altogether. It would
require, e.g. a
(0.25, in the present work) between cluster and comparison field for the differential reddening to affect the subtraction in a given cell. However, in both cases the CMDs extracted from the cluster region and comparison field (Figs. 3 and 4) indicate that the differential reddening, if present, is not important.
![]() |
Figure 3:
2MASS CMDs extracted from the
![]() ![]() ![]() ![]() |
Open with DEXTER |
Table 2: Statistics of the field-star decontamination discriminated by magnitude bins.
![]() |
Figure 4:
Same as Fig. 3 for the region
![]() ![]() |
Open with DEXTER |
The decontaminated stellar surface-density distributions (bottom panels of Fig. 2) are an indicator of the algorithm efficiency. In both cases the central cells, which correspond to the location of the overdensities, present conspicuous surface densities. Also, the residual surface-density around the centre was reduced to a minimum level. We note that the decontamination is essentially based on the colour-magnitude distributions of the stars in different spatial regions. The fact that the decontaminated surface density has a conspicuous excess only at the assumed cluster region (Fig. 2) implies significant differences among the spatial regions, both in the colour-magnitude distributions and in the number of stars within a given colour-magnitude bin. This is expected of star clusters, which can be characterised by a single stellar population, projected against a Galactic stellar field.
The decontaminated CMDs are shown in the bottom panels of Figs. 3 and 4. As expected, essentially all of the bulge and disk contamination in FSR 1716 is removed, leaving stellar sequences typical of a reddened old OC, with a rather populous main sequence (MS), well-defined main sequence turnoff (MSTO), and the giant branch. Alternatively, the decontaminated morphology might resemble that of a poorly-populated GC. We explore this possibility in Sect. 5.1. Of the 588 stars present in the
CMD, only
128 remain. A similar conclusion applies to Cz 23, in which the disk contamination has been subtracted, revealing a poorly-populated, old OC; 66 of the 228 stars remain in the CMD after decontamination.
The approximately central direction of FSR 1716 and the poorly-populated nature of Cz 23 require
additional statistical analysis. For this purpose, we present in Table 2 the full statistics of the decontaminated sequences and field stars, discriminated by magnitude bins. Statistically relevant parameters that characterise the nature of a star cluster are: (i)
which, for a given magnitude bin, corresponds to the ratio of the decontaminated number of stars to the
Poisson fluctuation of the number of observed stars; (ii)
,
which is related to the probability that the decontaminated stars result from the normal star count fluctuation in the comparison field and; (iii)
,
which measures the star-count uniformity of the comparison field. Properties of
,
,
and
,
measured in OCs and field fluctuations are discussed in Bica et al. (2008). Table 2 also provides integrated values of the above parameters, which correspond to the full magnitude range spanned by the CMD of each OC. The spatial regions considered here are those sampled by the CMDs shown in the top panels
of Figs. 3 and 4.
Star cluster CMDs should have integrated
values significantly higher than 1 (Bica et al. 2008),
a condition that is met by FSR 1716 (
)
and Cz 23 (
). Also, the number of decontaminated stars in each magnitude bin of FSR 1716 is higher (at the
level or higher) than what could be expected from field-star fluctuations. As a further test of the statistical significance of the above results we investigate star count properties of the field stars. First, the comparison field is divided into 8 sectors around the cluster centre. Next, we compute the parameter
,
which is the
Poisson fluctuation around the mean of the star counts measured in the 8 sectors of the comparison field (corrected for the different areas of the sectors and cluster extraction). In a spatially uniform comparison field,
is expected to be low. In this context, star clusters should have probable number of member stars (
)
higher than
,
to minimise the probability that
arises from fluctuations of a non-uniform comparison field. This condition is fully satisfied, in some cases reaching the level
.
The ratio decreases somewhat for Cz 23, probably because it has a low
number of member stars. We also provide in Table 2 the parameter
.
For a given magnitude bin we first compute the average number of stars over all sectors
and the corresponding
fluctuation
;
thus,
is defined as
.
Non-uniformities such as heavy differential reddening should result in high values of
.
Since we usually work with comparison fields wider than the possible-cluster extractions, the correction for the different spatial areas between field and cluster is expected to produce a fractional number of probable field stars (
)
in some cells. Before the cell-by-cell subtraction, the fractional numbers are rounded off to the nearest integer, but limited to the number of observed stars in each cell
,
where NI represents the nearest
integer. The global effect is quantified by means of the difference between the expected number
of field stars in each cell (
)
and the actual number of subtracted stars (
). Summed over all cells, this quantity provides an estimate of the total subtraction efficiency of the process,
The qualitative and quantitative expectations of the decontamination algorithm have been met by the output. On the one hand, the decontaminated photometry reveals a very high excess in the surface-density distributions with respect to the surroundings, in both cases (Fig. 2). In addition, CMDs extracted from the spatial regions where the excesses occur present statistically significant (Table 2) cluster-like features (Figs. 3 and 4).
Evidence drawn from the previous sections indicates that FSR 1716 is an old star cluster. However,
since it lies at a low latitude and is projected not far from the dense bulge stellar field, we provide an additional test to further probe its nature. A ``test-field'' of 20
of radius was extracted at
to the East of FSR 1716, with Galactic coordinates
and
,
thus projected somewhat closer to the centre than FSR 1716. The ``test-field'' photometry was analysed in the same way as that of FSR 1716. Similarly to FSR 1716, we considered a central region (
)
of the ``test-field'', extracted the equal-area comparison field near the border (
), and applied the decontamination algorithm. The results are shown in Fig. 5. Contrarily to the old-cluster decontaminated CMD morphology of FSR 1716 (Fig. 3), the decontaminated CMD of the ``test-field'' (bottom panels of Fig. 5) contains Poisson noise produced by the mutual subtraction of star samples that share similar distributions of magnitude and colours. The same conclusion
can be drawn from the featureless ``test-field'' RDP (Sect. 6; Fig. 10).
As discussed in Bica et al. (2008), Poisson fluctuations of the dense stellar field projected towards the bulge do not produce cluster-like CMDs and RDPs simultaneously. In particular, when field-fluctuation CMDs are field-star decontaminated, what results is represented by the above FSR 1716 ``test-field'' experiment (Fig. 5). In this sense, the striking differences exhibited by the CMD and RDP of FSR 1716, as compared to those of the ``test-field'', can be taken as robust evidence of the old-cluster nature of FSR 1716.
The CMD of Cz 23, especially the decontaminated one (Fig. 4), suggests a poorly-populated
old OC. Thus, we test whether the decontamination procedure, applied to random fields around Cz 23,
produces similar cluster-like sequences and stellar radial profiles. We consider
comparison fields, taken at 20
to each side of Cz 23. These fields have the same projected
area as that of Cz 23 shown in Fig. 4, and the centre to centre offsets correspond to
4 times the cluster radius of Cz 23 (Table 3). Field-star decontamination was applied to these CMDs using the same offset field as that employed for Cz 23 (Sect. 4), but restricted to
to avoid self subtraction.
The observed and decontaminated comparison field CMDs are shown in Fig. 6. In stark contrast to the decontaminated CMD of Cz 23 (Fig. 4), what remains in the
comparison field CMDs is a randomly scattered small number of stars, similar to the CMD of the FSR 1716 ``test-field'' (Sect. 4.1). To complete this analysis we also include in Fig. 6 the RDPs of the comparison fields. Both CMDs and RDPs of the comparison fields are typical of Poisson fluctuations of stellar fields projected towards the Galactic anti-centre (Bonatto & Bica 2008b).
As expected of field fluctuations, the CMDs and RDPs of the comparison fields (Fig. 6) are featureless, in contrast to the cluster CMD (Fig. 4) and RDP (Fig. 11) of Cz 23.
![]() |
Figure 5:
Same as Fig. 3 for the
![]() ![]() |
Open with DEXTER |
![]() |
Figure 6:
Left panels: observed and decontaminated CMDs extracted from
![]() ![]() |
Open with DEXTER |
![]() |
Figure 7:
Top panels: the field-decontaminated (
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 8:
Age and metallicity interpretation of FSR 1716. Top-left panel: isochrones of different ages and fixed solar metallicity. Top-right: same ages but with the sub-solar metallicity
![]() ![]() ![]() |
Open with DEXTER |
Having established that FSR 1716 and Cz 23 present old-age features, we proceed to derive
their fundamental parameters. We work with Padova isochrones (Girardi et al. 2002) of solar
and sub-solar metallicity, computed with the 2MASS J, H, and filters
. These isochrones are very similar to the Johnson-Kron-Cousins ones (e.g. Bessel & Brett 1988), with differences of at most 0.01 in (J-H) (Bonatto et al. 2004). We adopt
0.3 kpc (Bica et al. 2006a)
as the Sun's distance to the Galactic centre to compute the OC's Galactocentric distances. The value
0.3 kpc was derived by means of the GC spatial distribution. Other recent studies
gave similar results, e.g.
0.9 kpc (Eisenhauer et al. 2003),
0.32 kpc (Eisenhauer et al. 2005) and
0.10 kpc (Nishiyama et al. 2006), with different approaches. We make fits by eye, taking into account that, because of the presence of binaries, the isochrone should be shifted to the left of the MS fiducial line (e.g. Bonatto et al. 2005, and references therein).
In principle, the decontaminated CMD morphologies (bottom panels of Figs. 3 and 4) - as well as the additional tests discussed in Sect. 4.1 - should provide enough constraints to derive reliable cluster fundamental parameters. However, in view of the importance of identifying new old star clusters, we explore a range of possibilities in terms of age and metallicity. As discussed by, e.g. Friel (1995), OC metallicities in general range from solar (
,
or Z=0.019) to sub-solar (
,
Z=0.006) values.
First, we compare in Fig. 7 (top panels) the decontaminated CMD of FSR 1716 to that of
the 7 Gyr OC NGC 188 (Bonatto et al. 2005, and references therein). CMD morphologies in
both colours show a good agreement along the MS, MSTO, giant branch, and the MS width, which
suggests that both clusters are of similar age. In particular, both CMDs present a similar
smoothly-curved MSTO and a scarcely populated red clump.
Different ages and metallicities are tested in Fig. 8 solar-metallicity solutions
with isochrones of ages 5, 7, and 9 Gyr are shown in the top-left panel. The same ages but with
the low OC metallicity range (
,
which corresponds to
), are
considered in the top-right panels. As a fit constraint, the isochrones have in common a similar
description of the MS. In both panels, the youngest isochrone produces a relatively poorer fit,
especially to the giant branch. On the other hand, the older isochrones - of both metallicity
ranges - produce qualitatively similar fits. This not only confirms the old age of FSR 1716,
but, together with a few blue stars at
and
(Fig. 7),
which might suggest a blue horizontal branch (HB), raise the possibility of a poorly-populated,
metal-poor GC.
To test the GC hypothesis we compare the decontaminated CMD of FSR 1716 to those of the blue HB GCs NGC 6397 and M 4 (NGC 6121). Both GCs are among the nearest ones (
kpc),
which provides MS depth for 2MASS photometry, and are metal-poor, with
(NGC 6397) and
(M 4). We also use Padova isochrones of 12 Gyr to characterise the GC ages (e.g. Ortolani et al. 1995). The CMDs of NGC 6397 (extracted within
)
and M 4 (
)
were produced in the same way as those of FSR 1716 (Sect. 3). For a
quantitative derivation of fundamental parameters by means of a comparison with GCs, the isochrones
were first set to the GCs according to the fundamental data taken from Harris (1996, and the 2003 update
). Subsequently, we
searched for the best overall match between the
and FSR 1716 sequences.
The results are shown in the middle panels of Fig. 8. As a further test, we also
searched for solutions with the same isochrones and template GCs as above, but assuming a
brighter MSTO for FSR 1716, located at
.
These tentative fits are shown in the
bottom panels of Fig. 8
. Compared to NGC 6397 (bottom-left
panel), the general fit is acceptable, although the RGB of FSR 1716 is somewhat redder, which might suggest a higher metallicity. However, the fit with the more metal-rich GC M 4 (bottom-right) is
obviously poorer. Besides the red giants which are too red compared to M 4, the sub-giant
branch morphologies clearly are not the same. In the case of the M 4-like solution, FSR 1716
would be as close (
kpc) as an old OC. In the case of the NGC 6397 solution,
FSR 1716 would be the nearest (
kpc) GC so far discovered.
The above comparisons show that a very low metallicity, like that of NGC 6397, is not compatible
with the giant branch of FSR 1716, since it is steeper and bluer in NGC 6397. M 4, on the other
hand, provides a better match for the giant branch. What makes this possibility less likely is that
the few extreme blue stars in the CMD of FSR 1716 do not have a counterpart in either NGC 6397 or
M 4, whose blue HBs are redder and well distributed in colour. The CMD shape for
in FSR 1716 is more compatible with the MSTO and MS of an old OC than the matched subgiant branches of the two GCs. However, the evidence of a Palomar-like (i.e. low-mass and loose structure) GC
(Bonatto & Bica 2008a), although marginal, cannot be ruled out with the existing data.
We call attention to the fact that it is not an easy task to establish the nature of a low-mass GC or a low-mass old OC based on a few candidate HB or giant clump stars. For instance, NGC 188 (and other old OCs) presents only a few UV-bright stars, as shown by, e.g. Landsman et al. (1998). In fact, both low-mass GCs and low-mass old OCs hardly have any HBs or clump giants, including NGC 188 with only 5 clump giants (Dinescu et al. 1995; and Landsman et al. 1998). As examples of low-mass GCs with a small number of UV-bright stars we mention Palomar 13 (Côté et al. 2002) with only a few, and AM-4 (Inman & Carney 1987) with none.
![]() |
Figure 9:
Age and metallicity interpretation of Cz 23. Top panels: isochrones of different
ages and fixed solar metallicity. Bottom: same ages but with the sub-solar metallicity
![]() |
Open with DEXTER |
Given the above circumstances, deeper observations would be necessary for a more conclusive definition of the nature of FSR 1716 as a GC (as originally suggested by Froebrich et al. 2007a) or a very old OC, with a metallicity that can be sub solar.
The CMD features of Cz 23 (bottom panels of Fig. 7) indicate a younger age than FSR 1716. Indeed, its field-decontaminated CMD morphology is similar to that of the old OC M 67. The age of
M 67, as derived from near-infrared photometry, is 3.2 Gyr (Bonatto & Bica 2003, and references
therein). More recent determinations based on spectroscopy of individual stars put the age at
4 Gyr (e.g. Giampapa et al. 2006). In addition, the MSTO of Cz 23 seems to be at a slightly fainter magnitude range, which indicates an age older than that of M 67.
Thus, to probe the age and metallicity of Cz 23, we show in Fig. 9 fits with the 3, 5,
and 6 Gyr isochrones with a fixed solar metallicity (top panels), and the sub-solar metallicity
(bottom). The fits are required to provide a similar description of the MS. It is clear that, under this condition, the 3 Gyr isochrone fails to reproduce the giant clump, in both metallicity ranges. The 5 and 6 Gyr isochrones, on the other hand, provide acceptable MS and clump
fits, which confirms the old age of Cz 23. The 6 Gyr solutions are obtained with an essentially
zero reddening applied to the isochrones, which implies that older isochrones would require negative
correction. This argument constrains the age of Cz 23 to the range
4 to
6 Gyr.
Most of the arguments drawn above favour the old OC interpretation for the nature of FSR 1716. Since metal-poor OCs are preferentially distributed in the outer Galaxy (e.g. Friel 1995), and FSR 1716 is located inside the solar circle, its probable metallicity should be around the solar value.
Based on the above, we take the 7 Gyr, solar metallicity isochrone as representative of
the stellar sequences of FSR 1716, with an uncertainty of about 1 Gyr. The
corresponding best-fit was obtained with the reddening
0.02,
and observed distance modulus
0.2. This solution is shown in the top panels of Fig. 7. Taking into account fit uncertainties we derive the age 7
1 Gyr, reddening
0.1 or, equivalently (Sect. 3)
0.2, the absolute distance modulus
0.21, and the distance from the Sun
0.1 kpc. Thus, for
kpc, the Galactocentric distance
of FSR 1716 is
0.2 kpc, which puts it
0.6 kpc inside the solar circle.
Alternatively, we consider as well the less-likely possibility of a GC. In this case, the
best-fit of the 12 Gyr isochrone (based on the similarity with the M 4 CMD morphology -
Fig. 8) is obtained with
0.3,
0.04,
0.1 and
0.2,
0.32,
0.3 kpc, and
0.2 kpc, thus
1.8 kpc inside the solar circle.
With similar arguments as those used to estimate the age of FSR 1716, we take the 5 Gyr
(1 Gyr) isochrone to represent the age of Cz 23. Since both metallicity ranges
(
and 0.0) produce similar CMD fits (Fig. 9), for simplicity,
we also work with the solar metallicity isochrone.
Thus, fundamental parameters of Cz 23 are
0.01,
0.03 and
0.1,
0.1,
0.2 kpc, and
0.2 kpc; Cz 23 lies
2.5 kpc outside the solar circle. This solution is shown in the bottom panels of Fig. 7.
Structural parameters are derived by means of the RDPs, defined as the projected radial distribution of the number density of stars around the cluster centre.
Star clusters usually have RDPs that follow a well-defined analytical profile.
The most often used are the single mass, modified isothermal sphere (King 1966), the modified
isothermal sphere (Wilson 1975), and the power law with a core (Elson et al. 1987). Each function is characterised by different parameters that are related to cluster structure. However,
because the error bars in the present RDPs are significant (Fig. 10), we use the
analytical profile
,
where
is the residual background density,
is the central density of stars, and
is the core
radius. This function is similar to that used by King (1962) to describe the surface
brightness profiles in the central parts of globular clusters.
Table 3: Cluster structural parameters.
The RDPs are built with colour-magnitude filtered photometry to minimise noise. The most probable cluster sequences are isolated by means of colour-magnitude filters, which are used to exclude stars with colours different from those of the assumed cluster sequence. They are wide enough to include cluster MS and evolved star colour distributions, as well as the photometric uncertainties. Colour-magnitude filter widths should also account for formation or dynamical evolution-related effects, such as enhanced fractions of binaries (and other multiple systems) towards the central parts of clusters, since such systems
tend to widen the MS (e.g. Hurley & Tout 1998; Kerber et al. 2002; Bonatto & Bica 2007b; Bonatto et al. 2005). The
colour-magnitude filters for the present OCs are shown in the left panels of Fig. 7. Residual field stars with colours similar to those of the cluster are expected to remain inside the colour-magnitude filter region. They affect the intrinsic stellar RDP in a way that depends on the relative densities of field and cluster stars. The contribution of the residual contamination to the observed RDP is statistically subtracted by means of the field. In practical terms, the use of colour-magnitude filters in cluster sequences enhances the contrast of the RDP with respect to the background, especially for clusters in dense fields (e.g. Bonatto & Bica 2007b).
Oversampling near the centre and undersampling at large radii are avoided by using rings
of increasing width with distance from the centre. A typical set of ring widths is
,
and
,
respectively for
,
,
,
and
.
The number and
width of the rings can be set to produce RDPs with adequate spatial resolution and as low as
possible
Poisson errors. The residual background level of each RDP corresponds to the average number of colour-magnitude filtered stars measured in the field. The R coordinate (and uncertainty) of each ring corresponds to the average position and standard deviation of the
stars inside the ring.
![]() |
Figure 10:
Stellar RDPs (filled circles) of FSR 1716 ( left panels) and Cz 23 ( right). RDPs built with the observed (raw) photometry are shown in the top panels, while those built with colour-magnitude filtered photometry are in the bottom panels. Solid line: best-fit King-like profile. Horizontal shaded polygon: offset field stellar background level. Shaded regions: ![]() |
Open with DEXTER |
The RDPs of FSR 1716 and Cz 23 are given in Fig. 10. Because of the uncertainties
associated with the age (and metallicity) derivations (Sect. 5), which propagate to the absolute value of the structural parameters, RDPs in Fig. 10 are shown on an angular scale. Besides the RDPs built with the colour-magnitude filters, we show, for illustrative purposes, those produced with the observed (raw) photometry. Especially for FSR 1716, minimisation of the number of non-cluster stars by the colour-magnitude filter resulted in
a RDP with a higher contrast with respect to the background. Fits of the King-like profile
were performed with a non-linear least-squares fit routine that uses errors as weights. To minimise degrees of freedom,
and
were derived from the RDP fit, while
is measured in the field. These values are given in Table 3, and the best-fit solutions are
superimposed on the colour-magnitude filtered RDPs (Fig. 10). Table 3
also provides structural parameters in absolute units, computed with the cluster distances derived in Sect. 5.3. Because of the 2MASS photometric limit, which in most cases corresponds to a cutoff for stars brighter than
,
should be taken as a lower limit to the actual central number density. The adopted King-like function well describes the RDPs along the full radius range, within uncertainties.
To complete the structural description of the clusters we also estimate the cluster radius (
)
and
uncertainty by visually comparing the RDP level (and fluctuations) with the background. The cluster
radius corresponds to the distance from the cluster centre where the RDP and background are statistically indistinguishable (e.g. Bonatto & Bica 2005, and references therein). For practical purposes, most of the cluster stars are contained within
.
Note that
should not be mistaken for the tidal radius. Tidal radii are derived from King fits to RDPs, which depend on wide fields and adequate Poisson errors. For instance, in populous and relatively high Galactic latitude OCs such as M 67, NGC 188, and NGC 2477, cluster radii are a factor
0.5-0.7 of the respective tidal radii (Bonatto & Bica 2005). The cluster radii of the present objects are given in Cols. 6, 7 (angular scale) and 8, 9 (absolute scale) of Table 3. If cluster and tidal radii of FSR 1716 and Cz 23 are similarly related as in the bright OCs, the lower-limit of the radial range taken as a comparison field (Sect. 4) is located beyond the probable tidal radius. This, in turn, minimises the probability of cluster members at large radii, e.g. in the cluster halo, being considered as field stars by the decontamination algorithm.
Compared to the distribution of core radius derived for a sample of relatively nearby OCs by
Piskunov et al. (2007), FSR 1716 (especially the OC interpretation) and Cz 23 occupy the small-
tail. Assuming that the tidal radius is
2
,
Cz 23 and FSR 1716 (GC interpretation) would be around the median value of the distribution, while the OC interpretation of FSR 1716 again would be at the small tidal radius tail.
Table 3 (Col. 5) provides the density contrast parameter
.
Since
is measured in colour-magnitude-filtered RDPs, it may not correspond to the visual contrast produced by observed stellar distributions in images. FSR 1716 presents a low contrast in the 2MASS
image (Fig. 1) but, because most of the non-cluster stars have been excluded by the colour-magnitude filter, the corresponding RDP has a relatively high density contrast,
.
Obviously, the high-contrast OC Cz 23 (Figs. 1 and 10) is reflected in the high value of the density-contrast parameter
.
Interestingly, FSR 1716 and Cz 23 are projected against almost opposite regions (Table 1). Accordingly, the background (including foreground) stellar contribution turns out to be
in the direction of FSR 1716, and
towards Cz 23.
Finally, the RDP of the ``test-field'' of FSR 1716 is shown in the inset of panel (a) of Fig. 10. As expected from the CMD analysis (Sect. 4.1), the RDP is uniform from the centre to the borders of the field, characterised essentially by Poisson fluctuations. The featureless ``test-field'' RDP is low and flat, as compared the cluster RDP of FSR 1716 (panels a and c).
The cluster mass functions
are built following the methods presented in Bonatto & Bica (2005) (and references therein). We build them with colour-magnitude filtered photometry, the 3 2MASS bands separately, and the mass-luminosity relations obtained from the Padova isochrones and distances from the Sun adopted for FSR 1716 (Sect. 5.3.1) and Cz 23 (Sect. 5.3.2). Further details on MF construction are given in Bica et al. (2006b). The effective magnitude range over which MFs are computed is that where clusters present an excess of stars over the field. In both cases it begins below the MSTO and ends at a faint magnitude limit brighter than that stipulated by the 2MASS completeness limit (Sect. 3). The effective MS stellar mass ranges are
(FSR 1716), and
(Cz 23).
The MFs computed for the whole cluster region (
)
are shown in Fig. 11. Both cases suggest important depletion of the low-mass content, suggesting advanced dynamical evolution, especially in the case of FSR 1716. The drop in the MF of FSR 1716 begins at
,
which corresponds to
,
more than 1 mag brighter than the 2MASS 99.9% completeness limit (Sect. 3) at the cluster position. This suggests that the drop may be real, although we cannot rule out crowding to account for part of the important MF drop, since the 2MASS completeness limits refer to the field. Bearing in mind this caveat, we fit the MF of FSR 1716 with a two-segment function,
,
with
0.6 in the mass range
,
and
0.4 for
.
The high-mass range slope is, within the uncertainty, somewhat steeper than the
of the Salpeter (1955) Initial Mass Function (IMF). However, the low-mass range slope is much flatter than Salpeter's IMF. The drop in the MF of FSR 1716 (at
)
agrees with one of the breaks present
in the universal IMF of Kroupa (2001), which assumes increasing flattening towards low-mass stars. This IMF is described by the slopes
0.5 for the range
and
0.3 for
.
The low-mass range of the MF of FSR 1716 is flatter than Kroupa's IMF, which again points to advanced dynamical evolution, crowding or, more probably, a combination of both.
Although flat, the MF of Cz 23 is more monotonic than that of FSR 1716, and can be described by a
single power-law, with the slope
0.5, also flatter than Salpeter's or Kroupa's IMF.
![]() |
Figure 11:
Mass function computed for the whole region of FSR 1716 ( top panel) and Cz 23 ( bottom). The MF of FSR 1716 can be represented by the broken power-law
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
The number of observed MS and evolved stars in FSR 1716 is derived by counting the number
of stars (in the background-subtracted colour-magnitude filtered photometry) that lie in the
corresponding magnitude ranges,
for the MS and
for the evolved stars.
There are
5 and
7, MS and evolved stars, respectively; the corresponding mass values are
and
.
The evolved star mass corresponds to
multiplied by the stellar mass at the MSTO,
.
To estimate the total stellar
mass we extrapolate the observed MF down to the H-burning mass limit (
). Because of the very flat slope at the low-mass range, the extrapolated values are similar to the observed ones,
8 and
.
The number of MS and evolved stars in Cz 23 is
4 and
3, and the respective masses are
and
;
the extrapolated values are
8 and
.
Presently, both FSR 1716 and Cz 23 appear to be low-mass OCs, even when stars less massive than the observed range are taken into account. Most of the original mass must have been lost to the field because of dynamical effects during the several Gyr since their formation (Sect. 8).
The range of acceptable isochrone solutions (Sect. 5), together with the statistical tests
(Sect. 4) applied to FSR 1716 and Cz 23, i.e. (i) the decontamination algorithm; (ii) the integrated and per magnitude
parameter; and (iii) the ratio of
to
,
produce results consistent with both objects being old OCs. Examination of photometric and structural properties of the offset ``test-fields'' for FSR 1716 and Cz 23 also leads to the same conclusion.
With the analyses of the preceding sections, we provide fundamental and structural parameters, most of which are given for the first time for FSR 1716 and Cz 23. We now use these parameters to compare some of their properties with those of a set of well-studied OCs.
We compare the structural parameters, computed with the adopted isochrone solutions for FSR 1716
(Sect. 5.3.1) and Cz 23 (Sect. 5.3.2), with those of a reference sample of nearby OCs
with ages in the range 70-7000 Myr and masses within
400-5300
(Bonatto & Bica 2005).
To the original reference sample were added the young OCs NGC 6611 (Bonatto et al. 2006c) and NGC 4755
(Bonatto et al. 2006b). Clusters are distinguished according to total mass (lower or higher than 1000
). Bonatto & Bica (2005) discuss parameter correlations in the reference sample. For completeness, both sets of parameters for FSR 1716 (Table 3), corresponding to the old OC and GC possibilities, are considered separately.
Core and cluster radii of the reference sample relate as
(panel a), which suggests a similar scaling for both kinds of radii, at least for the sampled ranges of
and
.
Within the uncertainties, FSR 1716 (both interpretations) and Cz 23 fit tightly into that relation.
The reference OCs appear to follow the trend of increasing cluster radii with Galactocentric
distance (panel b), a dependence previously suggested by e.g. Lyngå (1982). To explain the
increase of GC radii with Galactocentric distance, van den Bergh et al. (1991) suggested that part of the
relation may be primordial, in the sense that the higher molecular gas density in central
Galactic regions may have produced clusters with small radii. After formation, mass loss
associated with stellar and dynamical evolution (such as mass segregation and evaporation),
together with tidal interactions with the Galactic potential and giant molecular clouds, also
contribute to the depletion of star clusters, especially the low-mass and centrally located
ones (Sect. 1). The cluster radius of FSR 1716 (GC interpretation) is consistent with
the assumed relation. On the other hand, Cz 23, and the OC interpretation of FSR 1716, appear
exceedingly small for their Galactocentric distance, which is consistent with the several Gyr
of depletion. A similar dependence on Galactocentric distance holds as well for ,
because of
the relation to
implied by panel (a).
In panels (c) and (d) of Fig. 12 we compare the cluster and core radii of FSR 1716 and
Cz 23 with those of the reference sample in terms of age. The locus occupied by Cz 23 in both
panels is consistent with the corresponding radii measured in low-mass, old OCs. A similar
conclusion applies to the GC interpretation of FSR 1716 (which again would be consistent with the
structural parameters of Palomar-like GCs - Bonatto & Bica 2008a). The significantly small radii
of the OC interpretation of FSR 1716, especially
,
are consistent with those of OCs that have
suffered severe depletion effects for long periods (e.g. Bonatto & Bica 2005; Bonatto & Bica 2007a,
and references therein).
![]() |
Figure 12:
Relations involving OC structural and fundamental parameters. Circles: nearby OCs.
Dotted circles: massive OCs (>
![]() |
Open with DEXTER |
![]() |
Figure 13: Galactic position of FSR 1716 and Cz 23 compared to the WEBDA OCs with ages younger (shaded circles) and older (filled circles) than 1 Gyr. The schematic projection of the Galaxy is seen from the North pole, with 7.2 kpc as the Suns distance to the Galactic centre. Main structures are identified. |
Open with DEXTER |
Further evidence in favour of dynamical effects affecting cluster size comes from the positions
of FSR 1716 and Cz 23 in the Galactic plane (Fig. 13). The Milky Way's spiral arm
structure is based on Momany et al. (2006) and Drimmel & Spergel (2001), derived from HII regions and molecular
clouds (e.g. Russeil 2003). The Galactic bar is shown with a 14 orientation and 3 kpc
in length (Freudenreich 1998; Vallée 2005). For comparison, we also include WEBDA clusters
younger and older than 1 Gyr. As discussed in Sect. 1, old OCs are distributed
preferentially outside the solar circle.
Interestingly, both clusters occur close to spiral arms, Cz 23 on Perseus, and FSR 1716 on
Carina-Sagittarius (OC interpretation) or Crux-Scutum (GC). Since they are located close to
the plane (Table 1), they may have interacted with the arms, especially by means of
encounters with giant molecular clouds. As discussed above, part of the
small sizes of FSR 1716 and Cz 23 can be accounted for by collisions with such structures
(e.g. Wielen 1971; Wielen 1991; Gieles et al. 2006; Gieles et al. 2007). A similar
effect was recently observed to occur with some small OCs located close to the Local and
Orion-Cygnus Arms (Bonatto & Bica 2008b).
Finally, both objects have MFs that appear to be very eroded and flat which, given the discussion on
structural parameters above, reflects the consequences of several Gyr of the relentless dynamical
effects and mass loss due to stellar evolution. Also, both objects are low-mass clusters, with
less than
presently stored in stars. They must have been formed as more massive open clusters to have survived for so long in the Galaxy.
In this paper we studied the nature of two stellar overdensities included in the catalogue of candidate star clusters of Froebrich et al. (2007a), FSR 1716 and FSR 834. The former was suggested to be a globular cluster candidate by Froebrich et al. (2007a), while the latter has the OC Cz 23 as optical counterpart. The analyses are based on field-star decontaminated 2MASS CMDs and stellar radial density profiles, with algorithms previously constructed by our group. Fundamental and structural parameters of the clusters are derived.
We present consistent evidence (e.g. CMD morphology, statistical tests, structural parameters, mass-function slope, and comparison with nearby OCs) that both objects are old star clusters.
FSR 1716 is significantly absorbed (
)
and projected not far from the bulge. Its
field-decontaminated CMD morphology is very similar to that of the
7 Gyr well-known OC NGC 188. Indeed, its CMD can be well represented by isochrones with ages older than
6 Gyr, both of solar and sub-solar (
)
metallicity. We adopted the former as the metallicity of FSR 1716, because of its relatively central location. Alternatively, we cannot rule out the possibility that FSR 1716 is a low-mass, loose (Palomar-like) globular cluster. FSR 1716 is located inside the solar circle,
0.6 kpc (in the case of an OC) or
1.8 kpc (GC).
The CMD morphology of Cz 23 is consistent with that of the 4 Gyr old OC M 67. Similarly to FSR 1716, it can be well represented by solar and sub-solar metallicity isochrones, but with ages
in the range 4-6 Gyr. With the
5 Gyr and solar metallicity solution, we find that Cz 23 is projected nearly towards the anti-centre, located
2.5 kpc outside the solar circle.
The core and cluster radii of FSR 1716 and Cz 23 are small when compared to a set of open clusters in the solar neighbourhood. However, such radii are comparable to those of other OCs of similar old
age. Also, the mass functions appear to be much flatter than Salpeter's IMF, especially FSR 1716,
which seems to present an increasing depletion in the number of low-mass stars. As a consequence, the total mass presently stored in stars in both clusters is lower than
.
Such a low value probably reflects the several Gyr-long period of mass loss due to stellar evolution, tidal interactions with the bulge (possibly in the case of FSR 1716), disk and giant molecular clouds. Because of its low mass content and flat mass function, Cz 23 may be evolving into an open cluster remnant (e.g. Pavani & Bica 2007).
Comprehensive catalogues of star cluster candidates, such as that of Froebrich et al. (2007a), should be
further explored with field-star decontamination algorithms and other tools, so that the nature
of the candidates can be probed and the age derived. It is remarkable how the decontamination
tool unveiled the intrinsic CMD sequences of FSR 1716, separating it from the crowded
field population. Consequently, the characterisation of FSR 1716 and Cz 23 as OCs older than
4 Gyr represents an important increase (
)
in the known population of such
objects in the Galaxy. In particular, FSR 1716 is the most recent addition to the 8 open clusters
older than 6 Gyr so far identified (WEBDA). In this sense, Cz 23 (FSR 834), and especially
FSR 1716, can be considered as rare examples of extreme dynamical survivors in disk-regions where
most open clusters are short-lived.
Acknowledgements
We thank the anonymous referee for helpful suggestions. This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Centre/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. This research has made use of the WEBDA database, operated at the Institute for Astronomy of the University of Vienna. We acknowledge support from the Brazilian Institution CNPq.