Hostname: page-component-7c8c6479df-p566r Total loading time: 0 Render date: 2024-03-29T12:36:22.030Z Has data issue: false hasContentIssue false

Complex Stellar System ESO65SC03: Open Cluster or Remnant?

Published online by Cambridge University Press:  09 June 2015

Gireesh C. Joshi*
Affiliation:
Centre of Advanced Study, Department of Physics, D.S.B. Campus, Kumaon University, Nainital 263002, Uttarakhand, India
Y. C. Joshi
Affiliation:
Aryabhatta Research Institute of Observational Sciences (ARIES), Manora peak, Nainital, Uttarakhand, India
S. Joshi
Affiliation:
Aryabhatta Research Institute of Observational Sciences (ARIES), Manora peak, Nainital, Uttarakhand, India
S. Chowdhury
Affiliation:
Aryabhatta Research Institute of Observational Sciences (ARIES), Manora peak, Nainital, Uttarakhand, India
R. K. Tyagi
Affiliation:
Department of Physics, Hemwati Nandan Bahuguna Govt. P.G. College, Khatima 262308, Uttarakhand, India
Rights & Permissions [Opens in a new window]

Abstract

We present a complete spatial and dynamical study of the poorly populated stellar system ESO65SC03. The radial distribution of the system gives a core and cluster radii of 1.10±0.63 and 5.36±0.24 arcmin, respectively. The SNDP does not show any clear enhancement of the surface stellar number density between the stars of the system and the field regions. We derive the optimum isochrone solution for a particular grid size in the CMD using the statistical cleaning procedure. Using the statistically cleaned CMDs, we find the distance modulus, (mM)0, and reddening, E(BV), of the system to be 11.8±0.2 and 0.45 mag, respectively. The mean proper motion of this system is − 5.37±0.81 mas yr−1 and 0.31±0.40 in RA and DEC directions, respectively. The mean proper motion of this system is found to be almost similar to the field region. The mass function for the brighter stars is found to be too high for the system to be an open cluster. These combined results place constraints on whether stellar system ESO65SC03 is a POCR or an Asterism. Our understanding is that the ESO65SC03 is in a stage of POCR by losing their main-sequence stars in the dynamic evolution processes.

Type
Research Article
Copyright
Copyright © Astronomical Society of Australia 2015 

1 INTRODUCTION

The comprehensive study of the Open Star Clusters (OCs) and their basic parameters provide a fundamental database to understand the Galactic evolution processes. Due to the dynamic evolution processes such as stellar encounters, interactions with the Galactic tidal field, etc., a significant fraction of the cluster stars are gradually lost to the field regions. The main effects of dynamical interactions are mass segregation and evaporation of low-mass members (Andersen & Nordström Reference Andersen, Nordström, Pallavicini, Micela and Scioortino2000; Patat & Carraro Reference Patat and Carraro1995; Carraro Reference Carraro2006). Studying the sparse region (the halo) along with the dense central region (the core) is therefore necessary to understand the evolution of any cluster. It is believed that the poorly populated OCs (PPOCs) do not survive longer than a few hundred Myr, whereas rich ones may survive for a longer period (Tanikawa & Fukushige Reference Tanikawa and Fukushige2005; Carraro et al. Reference Carraro, Geisler, Moitinho, Baume and Vázquez2005; Bonatto, Lima, & Bica Reference Bonatto, Lima and Bica2012). The photometric and structural properties of such PPOCs are similar to those of Asterisms which sometime creates conflicting classifications (Bica & Bonatto Reference Bica and Bonatto2011). The PPOCs may have a radial density profile (RDP) characterised essentially by a central overdensity together with significant noise outwards, typical features of an asterism. In these cases, only the colour–magnitude diagram (CMD) morphology may provide clues to the overdensity nature (Bonatto & Bica Reference Bonatto and Bica2008). In order to understand the difference between sparse clusters and asterisms as well as to find out any possible connection between them, we studied a poorly populated stellar system with a stellar enhanced region including PPOCs. For this, we searched PPOC systems in the WEBDA database and all available catalogues on VIZIER services. Kharchenko et al. (Reference Kharchenko, Piskunov, Schilbach, Roëser and Scholz2013) and Dias et al. (Reference Dias, Monteiro, Caetano, Lepine, Assafin and Oliveria2014) compiled several catalogues for optically visible and candidate OCs containing several basic spatial parameters of the cluster. We selected one such poorly populated stellar system ESO65SC03 for this study and extracted the data from various catalogues namely, Two Micron All Sky Survey (2MASS) (Skrutskie et al. Reference Skrutskie2006), Wide-field Infrared Survey Explorer (WISE) (Wright et al. Reference Wright, Eisenhardt and Mainzer2010), SPM4.0 (Girard et al. Reference Girard2011), and UCAC4 (Zacharias et al. Reference Zacharias, Finch, Girard, Henden, Bartlett, Monet and Zacharias2013).

The 2MASS (Skrutskie et al. Reference Skrutskie2006) has provided the photometric data of 410 million objects in J (1.25μm), H (1.65μm), and K (2.17μm) bands. The limiting sensitivity of J, H, K bands with 3σ are 17.1, 16.4, and 15.3 mag, respectively. The $\text{{\it JHK}--2MASS}$ system gives more accurate isochrone fitting for the lower portions of CMDs which improve the turnoff points, distance modulii, ages, and photometric membership of the clusters. Furthermore, this catalogue is a homogeneous database which allows us to observe young clusters in their dusty environments and reach out to their outer regions which is dominated by low-mass stars (Tadross Reference Tadross2011). The WISE (Wright et al. Reference Wright, Eisenhardt and Mainzer2010) is a NASA Medium Class Explorer mission that conducted a digital imaging survey of the entire sky in the mid-IR bands. The effective wavelength of mid-IR bands are 3.4 (W 1), 4.6 (W 2), 12 (W 3), and 22 μm (W 4). The SPM4.0 Catalogue provides position, absolute proper motions, and photographic BV photometry for over a hundred million stars and galaxies between the south celestial pole and − 20° declination and it is roughly complete down to V = 17.5. It is also based on the photographic and CCD observations taken with the Yale Southern Observatory's double-astrograph at Cesco Observatory in El Leoncito, Argentina (Girard et al. Reference Girard2011). UCAC4 (Zacharias et al. Reference Zacharias, Finch, Girard, Henden, Bartlett, Monet and Zacharias2013) is an all-sky star catalogue covering mainly 8 to 16 mag stars in a single passband between V and R. Furthermore, the positional errors are about 15 to 20 mas for stars in the 10 to 14 mag range. The proper motion in this catalogue have been derived for around 113 million stars utilising about 140 other star catalogues (including old SPM catalogue, Hipparcos and Tycho) with significant epoch difference to the UCAC4 CCD observations.

In the present paper, we derive the basic parameters and study the dynamic evolution of the stellar system ESO65SC03. This system is described as an optically visible sparse OC in the Southern hemisphere (Kharchenko et al. Reference Kharchenko, Piskunov, Schilbach, Roëser and Scholz2013). The basic parameters and the dynamical properties of this system has never been studied in detail. Along with studying the dynamic evolution of this system, we also aim to find out whether it is a Possible Open Cluster Remnant (POCR) or Asterism.

2 SPATIAL STRUCTURE OF SYSTEM

The classification of poorly populated stellar systems like ESO65SC03 is a challenging task due to the similarities in overdensity between sparse OCs and asterisms. The estimation of structural parameters and dynamical properties of poorly populated stellar systems can be used as an excellent probe for their classification and to find out any evolutionary connection between them. Using the estimated radius in Section 2.1, the extent of the system is drawn with a black circle in the W 1 band of WISE image as illustrated in Figure 1. The detailed estimation of the system parameters are discussed in the following sub-sections.

Figure 1. The finding chart $(20\,\text{arcmin}{\,\times\,} 20\,\text{arcmin})$ of stellar system ESO65SC03 in W 1 band of WISE taken from IRSA webportal. The field region starts after the radial distance of 5.36 arcmin (shown in black circle) and ends in 7.58 arcmin (shown in blue circle).

2.1 Radius of stellar system

A structural study of the system ESO65SC03 has been carried out using those stars which were detected in each of the 2MASS bands (J, H, and K) and their corresponding errors in magnitude < 0.1 mag. These stars were used to identify the system coordinates using maximum stellar counts per unit area. The radial stellar densities were estimated by taking the concentric rings from the centre of system with unit arcmin width, having maximum stellar counts and fitted with the King model (King Reference King1966). If the fit was not perfect, we selected the next unit cell with maximum stellar counts until we achieved the best fit. Finally, we determined the centre of the stellar system as $\alpha = 12^{\text{h}}:51^{\text{m}}:19.7^{\text{s}}, \delta = -69^{\circ} :43^{{\prime}}:21.6^{{\prime \prime}}$ and the equivalent Galactic coordinate as l = 302.91d, b = −6.84d which is consistent with the Kharchenko et al. (Reference Kharchenko, Piskunov, Schilbach, Roëser and Scholz2013).

The derived RDP of the system is shown in Figure 2. The modified king empirical formula [ρr − ρ0 = f 0(1 + (r/rc)2)−1, Kaluzny & Udalski (Reference Kaluzny and Udalski1992)] is drawn with the RDP of the cluster and it gives the peak stellar density (f 0), core-radius ($r_{\text{c}}$), and background stellar density (ρ0) as $2.3{\,\pm\,} 1.3\ \text{stars\,}\text{arcmin}^{-2}$, $1.1{\,\pm\,} 0.6\ \text{arcmin}$, and $3.2{\,\pm\,} 0.1\ \text{stars\,}\text{arcmin}^{-2}$, respectively. Here, the core radius of any stellar system is the radial distance from its centre where the enhanced stellar density above background level becomes half of that of the centre. The radius of the system is defined as the distance from the centre of the system to the point where the stellar enhanced density is not separated from the field region. After adding the 1σ correction in the background stellar density, we found the radius to be 5.36±0.24 arcmin.

Figure 2. The radial density profile (RDP) for the field of ES065SC03. The blue-dashed line represents the best fit empirical profile on the cluster RDP. The pink line represents ρ0 + σρ0, where σρ0 is the error in the estimation of field stellar density (ρ0).

Limiting radius of any stellar system is defined as a radial distance, after which no stars are gravitationally bound to the stellar system. It is independent of the background stellar density and only depends on peak stellar density, core radius, and estimation in the error of the background stellar density. The Bukowiecki et al. (Reference Bukowiecki, Maciejewski, Konorski and Strobel2011) relation ($r_{\text{lim}} = r_{\text{c}} \sqrt{f_{0}/3 \sigma _{\text{bg}} -1}$) is used to estimate the limiting radius ($r_{\text{limit}}$) of the cluster that comes out to be $2.96{\,\pm\,} 1.19\ \text{arcmin}$, which is significantly smaller than the radius of the stellar system. Since the system is gravitationally bound within ~ 3 arcmin and the model stellar radius is larger than expected gravitational attraction, therefore, the said stellar system may not be fall on the category of OCs. We found the spatial components of this system as X = 1.24 kpc, Y = − 1.92 kpc, and Z = − 0.27 kpc. The Galactocentric distance of the cluster is computed to be about 7.51 kpc. Our estimated core-radius and system-radius are in close agreement with those determined by the Kharchenko et al. (Reference Kharchenko, Piskunov, Schilbach, Roëser and Scholz2013).

2.2 Surface number density profile

The number of stars per magnitude interval in unit area is defined as the surface number density (SND), σ(r, m), and is a function of stellar magnitude and radius from the centre (Sung et al. Reference Sung, Lee, Lee and Ann1996). We computed the σ(r, m) for both the system and field regions as drawn in Figure 3. We see that the σ(r, m) of each region increases towards the higher magnitude end of each region. The SND profile (SNDP) of the system indicates that the system region have enhanced brighter stars compared to the fainter stars. The system stars are further divided into two groups through their SNDPs namely (a) brighter stars group (BSG), (b) fainter stars group (FSG). A separation at about 13 mag has been taken as cut-off magnitude between BSG and FSG group of stars. A vertical line separating the two groups for both the regions is shown in Figure 3. The SNDPs of system and field regions give different slopes ($m=\frac{\Delta \sigma (r,m)}{\Delta M_{H}}$) for BSG and FSG, which is estimated as 0.40±0.05 and − 0.28±0.13, respectively. For the field region, the slope of BSG and FSG are estimated as 0.60±0.11 and − 0.07±0.09, respectively. The slope of the BSG for the field region is found to be steeper than the system region while opposite is found for FSG.

Figure 3. The distribution of logarthim value of stellar number density versus J magnitude. The black points represent the logarithm SND values for the system region while blue points represent same for the field region. The solid vertical line has been used to separate the BSG and FSG stars of the stellar system.

2.3 Colour–magnitude diagram (CMD)

There is a possibility of some stars getting rejected from the system membership, if the chosen field region has more area than the system region. Whereas, if the area of the chosen field is less, the probability of the field stars in the stellar system is increased. In either case, the observed CMD must get deformed, thereby increasing the inaccuracy in estimation of the system parameters. Therefore, we have considered the area of the field region as 90.29 $\text{arcmin}^2$ which is equal to the area of the stellar system. Our field region begins after the radial distance of 5.36 arcmin from the centre as shown in Figure 1. The field star decontamination is necessary for reducing the colour-scattering on the CMDs. The field stars are removed by comparing the magnitude and colour similarities of stars between the stellar system to the field region. Based on the following criteria, the stars are excluded from the stellar enhanced region:

  • The obtained J-band magnitude difference between star in the system is not more than of 0.20 mag from the star of the field region at the same place in the CMD.

  • The colour difference (JH) of stars between system and field regions is not more than 0.05 mag.

  • Similarly, the colour difference (HK) of stars between system and field regions is not more than 0.05 mag.

We define a term $D_{\text{mc}}= \sqrt{\Delta _{J}^2+ \Delta _{JH}^2 + \Delta _{HK}^2}$ which is the summation of colour and magnitude differences of the field star and the star in the system, which fall within the grid of above specifications. We removed the star from the grid whose $D_{\text{mc}}$ is minimum. After iterating the procedure for each field region star, we found 158 stars on the statistically cleaned CMDs.

The statistical approach depends on the colour and magnitude grid size. If the boundaries of the grid size is increased, the system stars may get excluded from the membership of the stellar system and if the grid size is decreased, there is a possibility for the field stars to remain in the stellar system. Thus, the optimum grid size is chosen in such a way that the main-sequence (MS) pattern appears in the observed CMD. Before the field star decontamination, no such pattern was visible while after the decontamination of the field stars, we see a clear MS pattern (Figure 4(b)). We noticed an almost vertical sequence of stars with (JK)~0.8 and J fainter than 13 mag in the $2\text{MASS}$ decontamination CMD. This vertical sequence represents residual field star contamination. These residual field stars are removed from the decontaminated CMD by taking grid size as $\Delta {J}=0.5,\ \Delta {(J-H)}=0.2,\ {\text{and}}\ \Delta {(H-K)}=0.2,$ but exact cluster evolutionary sequence is not occurred in the decontaminated CMD (Figure 4(c)). In this scenario, the estimation of system parameters may not be possible, therefore, we have derived parameters through decontaminated CMD as obtained by applying the earlier grid size. The limitation of this approach, however, is that the cluster characteristics depend on the size of the grid.

Figure 4. J vs. (J–K) CMD for the system. (a) The cluster and the field stars are represented by red and green dots, respectively. (b) CMD for the remaining stars after the field star subtraction through statistical method (grid size, ΔJ = 0.20, Δ(JH) = 0.05, Δ(HK) = 0.05). The solid line represents the cluster evolutionary sequence while dotted lines represent the sequences beyond which CMDs may be contaminated by field stars. The black points represent possible residual field stars. (c) CMD for the remaining stars after the field star subtraction through statistical method with increased grid size (ΔJ = 0.50, Δ(JH) = 0.20, Δ(HK) = 0.20) (see text for detail).

The observed CMDs of the stellar system are used to determine the model distance, age, and reddening of the system by comparing them with the theoretical stellar isochrones. To test a suitable number of possible isochrone solutions, we first determined initial values for each physical parameter by visual inspection [log(τ)i, (mM)i, (JK)i] and subsequently, we investigated a large parameter space in an iterative manner. Our isochrone solutions typically span the physical ranges as log(τ) = log(τ)i ± 0.05,(mM) = (mM)i ± 0.50, and (JK) = (JK)i ± 0.10. We used the solar metallicity isochrones of log(age)=8.75 of UBVRIJHK photometric system (Marigo et al. Reference Marigo, Girardi, Bressan, Groenewegen, Silva and Granato2008) in J vs. (JK) CMD (Figure 4(b)). Through these CMDs, we found the apparent distance modulus of J-band i.e., (mM)J and log(age) of the system as 12.06 mag and 7.75 Myr, respectively. Furthermore, the colour-excess, E(JK), is estimated as 0.22 mag. The obtained log(age) value is in close agreement with that given by Kharchenko et al. (Reference Kharchenko, Piskunov, Schilbach, Roëser and Scholz2013). The obtained age of stellar system is however higher than that of sparse OCs (Bonatto et al. Reference Bonatto, Lima and Bica2012; Tanikawa & Fukushige Reference Tanikawa and Fukushige2005).

The Fiorucci & Munari (Reference Fiorucci and Munari2003) relations ($\frac{E({J-H})}{E({B-V})}= 0.309 {\,\pm\,} 0.130, \ \frac{E({J-K})}{E({B-V})} = 0.485 {\,\pm\,} 0.150$) have been used to estimate the reddening i.e. E(BV) and E(JH), which are estimated as 0.45 and 0.14 mag, respectively. The corresponding absolute distance modulus [(mM)0 = (mM)J − 1.125E(JK)] and distance obtained as 11.81±0.25 mag and 2.30±0.35 kpc, respectively. The error in these values is estimated by considering the errors of Fiorucci & Munari (Reference Fiorucci and Munari2003) relations. Kharchenko et al. (Reference Kharchenko, Piskunov, Schilbach, Roëser and Scholz2013) estimated the apparent distance modulus, reddening, and distance as 12.6 mag, 0.56 mag, and 3.05 kpc, respectively which is slightly different than our obtained values.

2.4 Two-colour diagrams (TCDs)

TCD are useful to derive the relationships among various colours. The colour ratio, $\frac{J-K}{V_{\text{ph}}-K}$ and $\frac{J-K}{(B-V)_{\text{ph}}}$, are obtained as 0.19±0.01 mag and 0.30±0.03 mag, respectively and the resultant plots are shown in Figures 6(a) and (b). The $(V_{\text{ph}}-J)$ vs. (JK) TCD are utilised to estimate the interstellar extinction in the near-IR, as shown in Figure 5. The best fit in this TCD gives a colour-excess of $E(V_{\text{ph}}-K)=0.90{\,\pm\,} 0.05$ mag and E(JK) = 0.17±0.19 mag. We used Whittet & van Breda (Reference Whittet and van Breda1980) relation to estimate the reddening vector, $R_{\text{cluster}}$=1.1E(VK)/E(BV) in the direction of system. We have found this vector as 2.2 which is far-off from the normal value of 3.1. Dutra, Santiago, & Bica (Reference Dutra, Santiago and Bica2002) relation $\frac{A_K}{A_V}=0.112$ yields the AK as 0.11. In fact, normal reddening law (RV = 3.1) is expected in the line of sight of the cluster in the absence of dust and intermediate stellar gases (Snedden et al. Reference Snedden, Gehrz, Hackwell, York and Snow1978). So, we have concluded that system region is affected by interstellar gases and dust.

Figure 5. The ($V_{\text{ph}}-K$) vs. (JK) TCD for statistically cleaned members of system. The green-dotted line shows the solar metallicity isochrone for log(age) = 8.75, while the two solid arrows indicate the direction of the normal reddening vector. The blue-dotted line is obtained by using reddening as E(VK) = 0.90 mag and E(JK) = 0.17 mag.

Figure 6. The ($V_{\text{ph}}-K$) vs. (JK) and $(B-V)_{\text{ph}}$ vs. (JK) (a) and (b) TCDs for the system's members which are found SPM4.0 catalogue. The dotted straight line represents the best linear fit. The lower panel represents the $(V_{\text{ph}}-\lambda)/(B-V)$ diagrams where λ have indicated any one photometric band from J, H, K, W 1, and W 2 bands.

Further, the TCD are used to examine the nature of variation of colour-excess with increasing wavelength and normal colour ratio values in the direction of system. For this purpose, the $(V_{\text{ph}}-\lambda)$ vs. $(B-V)_{\text{ph}}$ (λ is one of J, H, K, W 1, and W 2 bands.) plots are constructed in Figure 6(c) and the resultant colour ratios ($m_{\text{cluster}}$) are listed in Table 1. This clearly indicates that the value of $m_{\text{cluster}}$ increases from J to W 1. To derive the normal values for various colour-ratios, we have used the following approximate Neckel & Chini (Reference Neckel and Chini1981) relation:

\begin{equation*} m_{\text{normal}}=m_{\text{cluster}}{\times} \frac{R_{\text{normal}}}{R_{\text{cluster}}}. \end{equation*}

The present estimated normal values are close to known values. The variation of these values may be occurred because photographic plates are less sensitive to measurement of stellar magnitudes and have more errors.

Table 1. The various colour ratios are given in the table. The third column represents the estimated normal colour ratios through $R_{\text{cluster}}$ and Neckel and Chini relation (Neckel & Chini Reference Neckel and Chini1981) while the value in fourth column represents the expected normal colour ratios.

2.5 Proper motion of the system

We found a total of 158 stars in statistical cleaned J vs (J–K) CMD. We have also drawn a colour–magnitude filter around the best fitted isochrone as shown in Figure 4(b) by dotted lines. A total of 73 stars are found inside the colour–magnitude filter, which are defined as the MS stars and considered for the proper motion study (depicted by red dots in Figure 4(b)). The mean proper motion of this stellar system is estimated using two different catalogues. Using a criterion of 1-arcsec uncertainity, a total of 62 and 69 MS stars matched with the UCAC4 and SPM4.0 catalogues, respectively. The stars which do not fall within 3σ value of the mean are rejected. Using proper motions values from UCAC4 catalogue for 59 remaining MS stars after 3σ clipping, the mean proper motion values of the stellar system is found to be ${\bar{\mu}}_{x}$ = − 2.42±0.99 mas yr−1 and ${\bar{\mu}}_{y}$ = − 2.18±0.78 mas yr−1 in RA and DEC directions, respectively. Similarly, the mean proper motion values of remaining 61 MS stars, as found through SPM4.0 catalogue, are found to be − 5.37±0.81 mas yr−1 and 0.31±0.40 mas yr−1. One can see that both the catalogues have different number of stars and give different values of mean proper motion (shown in Figure 7). Since, the uncertainties in the stellar proper motion values of the SPM4.0 catalogue is less than UCAC4 catalogue, therefore, we considered the mean proper motion of the system ${\bar{\mu}}_{x}$ = − 2.42±0.99 mas yr−1 and ${\bar{\mu}}_{y}$ = − 2.18±0.78 mas yr−1 in RA and DEC directions, respectively as obtained through the former catalogue. Our obtained proper motion do not match with that estimated by Kharchenko et al. (Reference Kharchenko, Piskunov, Schilbach, Roëser and Scholz2013), using the PPMXL catalogue (${\bar{\mu}}_{x}$= − 14.54 mas yr−1, ${\bar{\mu}}_{y}$=4.63 mas yr−1).

Figure 7. The open circles are depicted the location of stars at $\mu _{x}\text{--}\mu _{y}$ plane based on the proper motion of the system in SPM4.0 catalogue for the member of the stellar system. The filled circles represent the same using the UCAC4 catalogue.

When we consider all the stars in Figure 4(b) without any colour–magnitude filter and adopting the same approach as above using the SPM4.0 catalogue, we found a mean proper motion of ${\bar{\mu}}_{x} = -4.96 {\,\pm\,} 0.56$ and ${\bar{\mu}}_{y} = 0.43{\,\pm\,} 0.32$ for the system. This value is similar to the mean proper motion ${\bar{\mu}}_{x} = -4.45{\,\pm\,} 0.32$ mas yr−1 and ${\bar{\mu}}_{y} = 0.48{\,\pm\,} 0.23$ mas yr−1 estimated for the field region. For gravitationally bound cluster system, this result is unexpected. Hence, we could not separate out the stars between the system and field regions on the basis of proper motion membership probabilities.

3 DYNAMICAL STUDIES OF SYSTEM

The dynamical properties provide the information about the member evaporation rate, moving direction of cluster members, and the strength of gravitational interaction within systems. In order to understand the dynamical behaviour of OCs, the mass function (MF) and mass segregation within the cluster are looked into. We studied the dynamical properties of the system by assuming two sample of stars. The first SAMPLE-A contains those stars which are found on statistically cleaned near-IR CMDs. The second SAMPLE-B consists of those members which are common among SAMPLE-A, 2MASS catalogue (near-IR survey), and WISE catalogue (mid-IR survey) within 1-arcsec uncertainty. Generally, the brighter stars are found in almost all filters but detection of fainter stars varies from filter to filter. This is because the completeness factor of stars gradually decreases towards the fainter end of stellar magnitude, as indicated by several completeness studies for different filters in the past (Joshi et al. Reference Joshi, Balona, Joshi and Kumar2014). By studying two different samples, we can understand the effect of stellar incompleteness on dynamical behaviour of the system. Table 2 lists all the stars of SAMPLE-A and SAMPLE-B in various H-band magnitude bins. The analysis of dynamical study of these samples is described below.

Table 2. The MF of the cluster ESO65SC03 derived from cleaned (H–K)/H CMD through (a) stars from 2MASS catalogue and (b) common stars between WISE and 2MASS (represented by suffix A and B, respectively).

3.1 Mass function

During the process of star formation, the distribution of stellar mass per unit volume is defined as the initial mass function (IMF) which determines the subsequent evolution of the cluster (Kroupa Reference Kroupa2002). The direct estimation of the IMF is not possible due to the dynamical evolution of the cluster. However, we can derive the MF, which is the relative number of stars per unit mass and can be expressed by a power law N(logM)∝M Γ. The slope, Γ, of the MF can be determined from

\begin{equation*} \Gamma = \frac{\text{d}\log N(\log {\it m})}{\text{d}\log {\it m}}, \end{equation*}

where Nlog(m) is the number of stars per unit logarithmic mass. The masses of probable cluster members can be determined by comparing observed magnitudes with those predicted by a stellar evolutionary model if the age, reddening, distance, and metallicity are known. The MF values estimated for the system are listed in Table 2, while the variation of MF with average mass in different magnitude bins are shown in Figure 8.

Figure 8. The variation of logarithm MF with the logarithm of average mass of members. These quantities have been computed using members as leftover after field star decontamination through the statistical procedure. The Figure A and B are represented the said plot for members of SAMPLE-A and SAMPLE-B.

To determine the MF slope, we have only used those stars which lie in the MS of the stellar system. The mostly brighter stars (H < 12 mag) may be giants, therefore, they have excluded in the estimation of MF slope. The MF slope value Γ is found to be 0.07±0.42 for the mass range from 2.44 $\text{M}_{\odot} $ to 0.85 $\text{M}_{\odot} $ which is low numerical value compared to the Salpeter value [Salpeter (Reference Salpeter1955) gave the classical slope of − 1.35 for the stars of mass range 0.4 $\text{M}_{\odot} $ to 10 $\text{M}_{\odot} $] and have positive slope. Such flatness/positive value of MF slope may be due to the dynamical evolution, rather than the intrinsic difference in the IMF (Nilakshi & Sagar Reference Nilakshi and Sagar2002). The stellar incompleteness increases toward the fainter end of photometric magnitude and consequently reduce the MF value thereby increasing the flatness of the MF slope. The MF values have also been computed for SAMPLE-B members of stellar system. The purpose of computing the MF values for SAMPLE-B members was to find out any correlation between MF values and the common stellar fraction of stars in various photometric regions. The MF value of the stars from SAMPLE-B increases as the average mass of the stars increases from 0.85 $\text{M}_{\odot} $ to 2.57 $\text{M}_{\odot} $. Brighter stars are vital members to identify any cluster. However, even after the cleaning, the field stars may still be present in the fainter end. The MF value depends on the number of stellar members, hence the MF value should increase in presence of the residual field stars. However, the stellar incompleteness gives the opposite result, which is much effective and prominent compared to the effect of residual field stars. The corresponding MF slope value for SAMPLE-B is 1.68±0.67 which is numerically comparable with the Salpeter value but opposite in slope. Since the stars of SAMPLE-B are found in both 2MASS and WISE catalogues, the common stellar fraction of stars in both catalogues gradually decreases towards the end of fainter stars. Due to the incompleteness of the data, in terms of the number of common stars, we expect the MF slope increases towards the positive slope value.

3.2 Mass segregation

Through the process of mass segregation, the stellar encounters of the cluster members are gradually increased. Furthermore, the higher-mass cluster members gradually sink towards the cluster centre by transferring their kinetic energy to the more numerous lower-mass stellar members of the cluster (Yadav, Sariya, & Sagar Reference Yadav, Sariya and Sagar2013). For the investigation of the dynamical evolution and mass segregation phenomenon due to the energy equipartition, the cluster stars are divided into two mass range of $2.567 \ {\le} \text{M}_{\odot} < 2.438$ and $2.438 \ {\le} \text{M}_{\odot} < 0.850$. The cumulative radial stellar distribution of stars for various mass ranges are shown in Figure 9. This clearly indicates that the mass segregation effect is present in the stellar system, and the radial commutative distributions of various mass groups are seen separately from system centre to a distance of about 4.2 arcmin from the centre. The Kolmogorov–Smirnov (K–S) test have been performed for verifying the same kind of mass distribution. Using the relation $\chi ^2=\frac{4 D^2 n_1 n_2}{n_1 + n_2}$, we have also estimated the χ2 value of mass segregation for the system, where D is the maximum distance between two cumulative radial distribution of stars and n 1 and n 2 are the number of stars found for various mass range. For the SAMPLE-A & SAMPLE-B, the χ2 value comes out to be 4.14 (confidence level is 87%) and 3.38 (confidence level is 81%), respectively. On the basis of this K–S test, we conclude that mass segregation phenomenon is more for the stars of SAMPLE-A than the other sample.

Figure 9. The radial density cumulative distribution of stars for various mass groups. The mass segregation phenomenon for the star groups are shown. The red and green lines represent the mass range: $2.567 \ {\le} \text{M}_{\odot} < 2.438$ and $2.438 \ {\le} \text{M}_{\odot} < 0.850$, respectively.

4 DISCUSSION

The PPOCs and Asterisms have similar structural properties (Bica & Bonatto Reference Bica and Bonatto2011) but only differ in their comprehensive structural and dynamical studies. The overdensity profile does not provide any clue for their type. However, SNDP may provide us with an in-depth knowledge of their characteristics. The Figures 4(b) and (c) show the dependency of observed MS on grid size of statistically cleaned procedure. The smaller grid size of ΔJ = 0.20, ΔJH = 0.05, ΔHK = 0.05 provides the MS pattern for the stellar system as well as the isochrone solution. We do not find any MS pattern on the observed CMD when the grid size is increased to ΔJ = 0.50, ΔJH = 0.20, ΔHK = 0.20 which is less than the typical grid sizes of sparse clusters. Thus, the grid size did not provide any useful information about the MS characteristics of the stellar system.

Using two different catalogues [PPMXL: Roëser, Demleitner, & Schilbach (Reference Roëser, Demleitner and Schilbach2010) and UCAC4: Zacharias et al. (Reference Zacharias, Finch, Girard, Henden, Bartlett, Monet and Zacharias2013)], Joshi et al. (Reference Joshi, Balona, Joshi and Kumar2014) found that the proper motion of the OC NGC 559 is similar to each other within their given uncertainties. Due to the low subtraction efficiency of the field decontamination procedure, the mean proper motion of the system may get shifted towards the proper motion centre of the field regions, as expected. The mean proper motion value of the field and the system regions, as estimated by us, are almost similar. As a result, the low subtraction efficiency of the field decontamination have negligible effect on the mean proper motion of the system. The stellar density of the system will increase due to the low subtracting efficiency, leading to the increased MF but the slope value of MF decreases due to the stellar incompleteness. In our study, the effect of latter is more dominant than the low subtracting efficiency. However, if the corrected completeness factor is applied to the stellar system, the slope of the MF may increase. Due to fainter members of ESO65SC03, we also found the flat MF slope indicating more evolved members. In addition, the MF values of brighter members of the system is too low for it to be an OC. The typical MF value for brighter members of any OC reported is much lower than our value (Lata et al. Reference Lata, Pandey, Sharma, Bonatto and Yadav2014; Joshi et al. Reference Joshi, Balona, Joshi and Kumar2014; Glushkova et al. Reference Glushkova, Zabolotskikh, Koposov, Spiridonova, Leonova, Vlasyuk and Rastorguev2013; Yadav et al. Reference Yadav, Sariya and Sagar2013; Tadross et al. Reference Tadross, El-Bendary, Osman, Ismail and Bakry2012; Tadross Reference Tadross2011; Bhatt, Sagar, & Pandey Reference Bhatt, Sagar and Pandey2012; Sharma et al. Reference Sharma, Pandey, Ogura, Aoki, Pandey, Sandhu and Sagar2008). There is also a possibility that the lighter members of this system has already evaporated from its central region as expected by the phenomenon of mass segregation.

5 CONCLUSION

The stellar enhancement over the field region is clearly visible in the RDP of ESO65SC03 and the enhancement in the SNDP is only due to the brighter members of the system. Since the SNDPs of the field and system regions are intersecting each other at several points, the system does not seem to be an OC. On the other hand, the limiting radius is found to be significantly smaller than the cluster radius which is unexpected for an OC system. We estimated the age, heliocentric distance, galactocentric distance, and reddening of the system as 562 Myr, 2.30±0.35 kpc, 7.51 kpc, and 0.45 mag, respectively, as listed in Table 3.

Table 3. The structural and dynamical parameters of stellar system have been listed here. In this table, G.D. and H.D. represent Galactocentric and Heliocentric distances, respectively. Other notations are described in the text.

The TCD for this system indicate that the colour excess values [E(Vph − λ)] increases from J-band to W 1-band. Our photographic colour-ratio values for the stellar system is less than the typical [$E({V_{\text{CCD}}-\lambda})$] values found through CCD observations. The MF value for the brighter stars is too low for the system to be an OC. However for stars with mass 2.57 < M < 0.85, the MF slope is found to be 0.07±0.42 indicating that most of the stars are evolved from the MS. The MF slope for the common stars between WISE and 2MASS catalogue is 1.68±0.67 which indicates that the MF slope increases towards the positive value due to the incompleteness effect. The system also shows the mass segregation phenomenon. The number of stars in the system are found to be about 6.5% more than that of the field region. Based on all the above said factors and considering the criteria suggested by Pavani et al. (Reference Pavani, Kerber, Bica and Maciel2011), we conclude that ESO65SC03 may not be a normal star cluster but a POCR. Since we get an isochrone solution for system by applying limiting grid size in statistical field star decontamination procedure, therefore, the possibility of system as a asterism is ruled out. We conclude that ESO65SC03 is reached at a stage of POCR by losing their MS stars in the dynamic evolution processes through the tidal effects and stellar encounters. A further multi-bands analysis is required to conclude the exact nature of this system.

ACKNOWLEDGEMENTS

GCJ is thankful to University Grants Commission (New Delhi) for the financial assistance in terms of an RFSMS fellowship. SC acknowledges the Indo–Russian (RFBR) project INT/RFBR/P-118 for the financial support. Authors are also thankful to the referee for his constructive comments and suggestions which helped us to improve the scientific contents. This publication made use of data products from the 2MASS, 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 publication used the data products from the WISE, which is a joint project of the University of California, Los Angeles and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration.

References

REFERENCES

Andersen, J., & Nordström, B. 2000, in ASP Conf. Ser. 198, Stellar Clusters and Associations, ed. Pallavicini, R., Micela, G., & Scioortino, S. (San Francisco: Astron. Soc. Pac.), 171Google Scholar
Bhatt, H., Sagar, R., & Pandey, J. C. 2012, NewA, 17, 160CrossRefGoogle Scholar
Bica, E., & Bonatto, C. 2011, A&A, 530, A32Google Scholar
Bonatto, C., & Bica, E. 2008, A&A, 485, 81Google Scholar
Bonatto, C., Lima, E. F., & Bica, E. 2012, A&A, 540, 137Google Scholar
Bukowiecki, L., Maciejewski, G., Konorski, P., & Strobel, A. 2011, AcA, 61, 231Google Scholar
Carraro, G. 2006, BASI, 34, 153Google Scholar
Carraro, G., Geisler, D., Moitinho, A., Baume, G., & Vázquez, R. A. 2005, A&A, 442, 917Google Scholar
Dias, W. S., Monteiro, H., Caetano, T. C., Lepine, J. R. D., Assafin, M., & Oliveria, A. F. 2014, A&A, 564A, 79Google Scholar
Dutra, C. M., Santiago, B. X., & Bica, E. 2002, A&A, 381, 219Google Scholar
Fiorucci, M., & Munari, U. 2003, A&A, 401, 781Google Scholar
Girard, T. M., et al. 2011, AJ, 142, 15CrossRefGoogle Scholar
Glushkova, E. V., Zabolotskikh, M. V., Koposov, S. E., Spiridonova, O. I., Leonova, S. I., Vlasyuk, V. V., & Rastorguev, A. S. 2013, MNRAS, 429, 1102CrossRefGoogle Scholar
Joshi, Y. C., Balona, L. A., Joshi, S., & Kumar, B. 2014, MNRAS, 437, 804CrossRefGoogle Scholar
Kaluzny, J., & Udalski, A. 1992, AcA, 42, 49Google Scholar
King, I. 1966, AJ, 71, 64CrossRefGoogle Scholar
Kharchenko, N. V., Piskunov, A. E., Schilbach, E., Roëser, S., & Scholz, R.-D. 2013, A&A, 558A, 53Google Scholar
Kroupa, P. 2002, Sci, 295, 82CrossRefGoogle Scholar
Lata, S., Pandey, A. K., Sharma, S., Bonatto, C., & Yadav, R. K. 2014, NewA, 26, 77CrossRefGoogle Scholar
Marigo, P., Girardi, L., Bressan, A., Groenewegen, M. A. T., Silva, L., & Granato, G. L. 2008, A&A, 482, 883Google Scholar
Neckel, T., & Chini, R. 1981, A&A, 45, 451Google Scholar
Nilakshi, & Sagar, R. 2002, A&A, 381, 65Google Scholar
Patat, F., & Carraro, G. 1995, A&ASS, 114, 281Google Scholar
Pavani, D. B., Kerber, L. O., Bica, E., & Maciel, W. J. 2011, MNRAS, 412, 1611CrossRefGoogle Scholar
Roëser, S., Demleitner, M., & Schilbach, E. 2010, AJ, 139, 2440CrossRefGoogle Scholar
Salpeter, E. E. 1955, ApJ, 121, 161CrossRefGoogle Scholar
Sharma, S., Pandey, A. K., Ogura, K., Aoki, T., Pandey, K., Sandhu, T. S., & Sagar, R. 2008, AJ, 135, 1934CrossRefGoogle Scholar
Skrutskie, M., et al. 2006, AJ, 131, 1163CrossRefGoogle Scholar
Snedden, C., Gehrz, R. D., Hackwell, J. A., York, D. G., & Snow, T. P. 1978, ApJ, 223, 168CrossRefGoogle Scholar
Sung, H., Lee, S., Lee, M. G., & Ann, H. B. 1996, JKAS, 29, 269Google Scholar
Tadross, A. L. 2011, JKAS, 44, 1Google Scholar
Tadross, A. L., El-Bendary, R., Osman, A., Ismail, N., & Bakry, A. 2012, RAA, 12, 75Google Scholar
Tanikawa, A., & Fukushige, T. 2005, PASJ, 57, 155CrossRefGoogle Scholar
Whittet, D. C. B., & van Breda, I. G. 1980, MNRAS, 192, 467CrossRefGoogle Scholar
Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868CrossRefGoogle Scholar
Yadav, R. K. S., Sariya, D. P., & Sagar, R. 2013, MNRAS, 430, 3350CrossRefGoogle Scholar
Zacharias, N., Finch, C. T., Girard, T. M., Henden, A., Bartlett, J. L., Monet, D. G., & Zacharias, M. I. 2013, AJ, 145, 44Google Scholar
Figure 0

Figure 1. The finding chart $(20\,\text{arcmin}{\,\times\,} 20\,\text{arcmin})$ of stellar system ESO65SC03 in W1 band of WISE taken from IRSA webportal. The field region starts after the radial distance of 5.36 arcmin (shown in black circle) and ends in 7.58 arcmin (shown in blue circle).

Figure 1

Figure 2. The radial density profile (RDP) for the field of ES065SC03. The blue-dashed line represents the best fit empirical profile on the cluster RDP. The pink line represents ρ0 + σρ0, where σρ0 is the error in the estimation of field stellar density (ρ0).

Figure 2

Figure 3. The distribution of logarthim value of stellar number density versus J magnitude. The black points represent the logarithm SND values for the system region while blue points represent same for the field region. The solid vertical line has been used to separate the BSG and FSG stars of the stellar system.

Figure 3

Figure 4. J vs. (J–K) CMD for the system. (a) The cluster and the field stars are represented by red and green dots, respectively. (b) CMD for the remaining stars after the field star subtraction through statistical method (grid size, ΔJ = 0.20, Δ(JH) = 0.05, Δ(HK) = 0.05). The solid line represents the cluster evolutionary sequence while dotted lines represent the sequences beyond which CMDs may be contaminated by field stars. The black points represent possible residual field stars. (c) CMD for the remaining stars after the field star subtraction through statistical method with increased grid size (ΔJ = 0.50, Δ(JH) = 0.20, Δ(HK) = 0.20) (see text for detail).

Figure 4

Figure 5. The ($V_{\text{ph}}-K$) vs. (JK) TCD for statistically cleaned members of system. The green-dotted line shows the solar metallicity isochrone for log(age) = 8.75, while the two solid arrows indicate the direction of the normal reddening vector. The blue-dotted line is obtained by using reddening as E(VK) = 0.90 mag and E(JK) = 0.17 mag.

Figure 5

Figure 6. The ($V_{\text{ph}}-K$) vs. (JK) and $(B-V)_{\text{ph}}$ vs. (JK) (a) and (b) TCDs for the system's members which are found SPM4.0 catalogue. The dotted straight line represents the best linear fit. The lower panel represents the $(V_{\text{ph}}-\lambda)/(B-V)$ diagrams where λ have indicated any one photometric band from J, H, K, W1, and W2 bands.

Figure 6

Table 1. The various colour ratios are given in the table. The third column represents the estimated normal colour ratios through $R_{\text{cluster}}$ and Neckel and Chini relation (Neckel & Chini 1981) while the value in fourth column represents the expected normal colour ratios.

Figure 7

Figure 7. The open circles are depicted the location of stars at $\mu _{x}\text{--}\mu _{y}$ plane based on the proper motion of the system in SPM4.0 catalogue for the member of the stellar system. The filled circles represent the same using the UCAC4 catalogue.

Figure 8

Table 2. The MF of the cluster ESO65SC03 derived from cleaned (H–K)/H CMD through (a) stars from 2MASS catalogue and (b) common stars between WISE and 2MASS (represented by suffix A and B, respectively).

Figure 9

Figure 8. The variation of logarithm MF with the logarithm of average mass of members. These quantities have been computed using members as leftover after field star decontamination through the statistical procedure. The Figure A and B are represented the said plot for members of SAMPLE-A and SAMPLE-B.

Figure 10

Figure 9. The radial density cumulative distribution of stars for various mass groups. The mass segregation phenomenon for the star groups are shown. The red and green lines represent the mass range: $2.567 \ {\le} \text{M}_{\odot} < 2.438$ and $2.438 \ {\le} \text{M}_{\odot} < 0.850$, respectively.

Figure 11

Table 3. The structural and dynamical parameters of stellar system have been listed here. In this table, G.D. and H.D. represent Galactocentric and Heliocentric distances, respectively. Other notations are described in the text.