Optical and Near-infrared survey of the stellar contents associated with the star-forming Complex Sh2-252
We present the analyses of the stellar contents associated with the extended H ii region Sh2-252 using deep optical photometry, slit and slitless spectroscopy along with the near-infrared (NIR) data from 2MASS for an area 1 degree 1 degree. We have studied the sub-regions of Sh2-252 which includes four compact-H ii (CH ii ) regions, namely A, B, C and E and two clusters NGC 2175s and Teutsch 136 (Teu 136). Of the fifteen spectroscopically observed bright stars, eight have been identified as massive members of spectral class earlier than B3. From the spectro-photometric analyses, we derived the average distance of the region as 2.40.2 kpc and the reddening of the massive members is found to vary between 0.35 to 2.1 mag. We found that NGC 2175s and Teu 136, located towards the eastern edge of the complex are the sub-clusters of Sh2-252. The stellar surface density distribution in -band shows clustering associated with the regions A, C, E, NGC 2175s and Teu 136. We have also identified the candidate ionizing sources of the CH ii regions. 61 H emission sources are identified using slitless spectroscopy. The distribution of the H emission sources and candidate young stellar objects (YSOs) with IR excess on the colour magnitude diagram (CMD) shows that a majority of them have approximate ages between 0.1 - 5 Myr and masses in the range of 0.3 - 2.5 M. The optical CMDs of the candidate pre-main sequence (PMS) sources in the individual regions also show an age spread of 0.1 - 5 Myr for each of them. We calculated the -band luminosity functions (KLFs) for the sub-regions A, C, E, NGC 2175s and Teu 136. Within errors, the KLFs for all the sub-regions are found to be similar and comparable to that of young clusters of age 5 Myr. We also estimated the mass functions (MFs) of the PMS sample of the individual regions in the mass range of . In general, the slopes of the MFs of all the sub-regions are found comparable to the Salpeter value.
keywords:stars: formation stars: luminosity function, mass function stars: premainsequence H ii regions: individual: Sh2-252
It is known that most of the stars in the galaxy are formed in clusters or in OB associations where massive members of the region significantly affect their environment. This is possible in the form of strong stellar winds, intense UV radiation and eventually supernova explosions. These processes can ultimately cause their natal molecular cloud to be destroyed, which, in turn, will put an end to further star formation. However, it has also been observed that the above processes can trigger the birth of new generations of stars (e.g., Elmegreen & Lada 1977) and is likely to influence the key properties of the star formation process such as the initial mass function (IMF), star formation efficiency and the evolution of protostellar disks around the young stars (e.g., Clarke 2007). Keeping in mind the above points, it is important to study the large scale properties of OB associations so that we can have better understanding of the various processes governing the star formation.
As a continuation of our multi-wavelength analyses of star forming regions (SFRs; e.g., Pandey et al. 2008, Jose et al. 2008; 2011), an optically bright H ii region Sh2-252 (Sharpless 1959; = ; = ; l=190.04; b=+0.48) has been studied in this paper. It is located in the galactic anti-center direction and is a part of the Gemini OB1 association. Radio and CO surveys have already been done for this region by many authors (Felli et al. 1977, Lada & Wooden 1979, Kömpe et al. 1989; Szymczak et al. 2000 etc.). However, stellar contents of this region was poorly studied in optical and infrared wavebands. Also, authors have raised arguments regarding the presence of many ionizing sources in this complex (see Felli et al. 1977, Lada & Wooden 1979). However, the massive members of this region have not yet been characterized spectroscopically except three bright sources. In this paper, we use the deep optical data in bands, slit and slitless spectroscopic observations along with the data from 2MASS to identify and classify the massive members, to determine the reddening, distance of the region as well as the age and mass of candidate pre-main sequence (PMS) sources and finally to obtain the K-band luminosity function (KLF) and IMF. After a brief description of the previous studies on Sh2-252, section 2 describes various data sets used for the present study. Analyses and results including fundamental parameters of Sh2-252, identification and classification of massive stars in the associated clusters, its various properties such as age, distance, reddening, KLF and IMF are discussed in section 3 and section 4 summarizes the results.
1.1 Overview of previous studies on Sh2-252
The low density (n = 9 ) and the relatively large size (D 20 pc) of the extended region place Sh2-252 in the class of evolved H ii regions (Churchwell 1974). Main source of ionization of Sh2-252 is the central star HD 42088 with a spectral type O6.5 V (Conti & Alschular 1971; Walborn 1972), which is a member of the Gemini OB1 association. In Fig. 1, the DSS2-R band image of the region around Sh2-252 for an area of square degrees is given, which shows a bright nebula having a diameter 30 with a bright rim (Felli et al. 1977) at the west. See the colour composite image given in Fig. 2 for a close up view. Two stellar clusters were identified in this region optically, one is the dispersed cluster (namely NGC 2175) centered on the nebulosity and the other of smaller diameter (namely NGC 2175s) is located at the north-east of the nebula. Pismis (1970) estimated the distances to these clusters as 1900 pc and 3500 pc, respectively and hence NGC 2175s was thought to be a background object. Subsequently, using eight colour photometry, Chavarría et al. (1987) found that these two clusters are located at a distance of 2.0 kpc and hence claimed that both are associated with Sh2-252. A new near-infrared (NIR) cluster Teutsch 136 (here after Teu 136; see Fig. 1) has been identified towards the east of NGC 2175s by Koposov et al. (2008). However, they did not make any parameter determination for this cluster. Hence, the association of this object with Sh2-252 needs to be checked. In the literature, the distance to Sh2-252 from spectrophotometric and kinematic studies vary between 1.4 kpc to 2.6 kpc (e.g., Bonatto & Bica 2011, Pismis 1970, Reid et al. 2009 and Grasdalen & Carrasco 1975). Using the 5 GHz aperture synthesis observations, Felli et al. (1977) detected six extended radio sources towards Sh2-252 labeled as Sh2-252 A to F (see Fig. 1 in Lada & Wooden 1979). They have identified four of these sources (A, B, C and E) as compact H ii (CH ii ) regions, probably having a local source of ionization in each of them with a spectral type later than HD 42088 (i.e., O6.5V). The region F coincides with the optical bright rim (see Fig. 1, and see also Fig. 2 for a close up view) and no internal heating source is found within it (Lada & Wooden 1979). The molecular survey by Lada & Wooden (1979) shows that the most intense CO peak is located very close to Sh2-252A with water and methanol maser emissions within its proximity (Lada & Wooden 1979; Kömpe et al. 1989; Szymczak et al. 2000), which indicates the recent star formation activity towards this region. Most of the region around Sh2-252A seems to be optically obscured, however, there is an embedded cluster visible in NIR and recently, Tej et al. (2006) has studied this cluster in detail.
2 Observations and Data Reductions
2.1 Optical CCD imaging
The CCD observations of Sh2-252 were carried out using the 105-cm Schmidt Telescope of the Kiso observatory, Japan, on 2004 November 4. The 2K 2K CCD with each pixel corresponding to 1.5 arcsec covers a field of arcmin in the sky. The log of observations is tabulated in Table 1.
|Date of||Filter & Exp. time (sec) no. of frames|
|+20:28:10||2004.11.04||: 602; B: 202; V: 102; R: 102; I: 102|
|06:09:39||+20:29:12||2005.12.02||: 3602; B: 2402; V: 1202; R: 602; I: 602|
|06:10:09||+20:34:57||2006.10.19||: 6006; : 3006|
|06:09:14||+20:35:60||2006.10.23||: 6006; : 3006|
|06:10:09||+20:25:28||2006.10.23||: 6006; : 3006|
|06:09:18||+20:24:52||2006.12.06||: 6006; : 3006|
|06:08:06||+20:20:20||2007.01.15||: 6006; : 3006|
|06:08:06||+20:20:20||2007.01.15||: 6006; : 3006|
|06:10:56||+20:35:29||2007.01.15||: 6006; : 3006|
|06:09:11||+20:25:22||2007.11.06||: 6006; : 3006|
|06:10:54||+20:36:18||2007.11.06||: 6006; : 3006|
|06:08:51||+20:12:36||2008.01.08||: 6006; : 3006|
|06:11:56||+20:38:28||2008.01.08||: 6006; : 3006|
|06:10:25||+20:13:46||2008.04.15||: 6006; : 3006|
|06:11:31||+20:34:35||2008.04.15||: 6006; : 3006|
|06:08:31||+20:39:24||2008.12.18||: 6006; : 3006|
|06:09:29||+20:46:24||2008.12.18||: 6006; : 3006|
|06:10:24||+20:45:31||2008.12.18||: 6006; : 3006|
|06:09:42||+20:13:32||2008.12.18||: 6006; : 3006|
|06:11:09||+20:24:36||2008.12.19||: 6006; : 3006|
|06:08:23||+20:26:44||2008.12.19||: 6006; : 3006|
|06:11:16||+20:46:54||2008.12.19||: 6006; : 3006|
|06:09:21||+20:40:37||2005.01.07||Gr5/H-Br: 4503; H-Br: 902|
|06:09:21||+20:24:31||2005.01.08||Gr5/H-Br: 4503; H-Br: 902|
|06:09:60||+20:32:38||2005.01.09||Gr5/H-Br: 4503; H-Br: 902|
|06:09:60||+20:40:39||2005.01.09||Gr5/H-Br: 4503; H-Br: 902|
|06:10:37||+20:40:46||2005.01.09||Gr5/H-Br: 6003; H-Br: 1202|
|06:10:37||+20:32:55||2005.02.16||Gr5/H-Br: 6003; H-Br: 1202|
|06:10:37||+20:24:28||2005.02.16||Gr5/H-Br: 6003; H-Br: 1202|
|06:09:23||+20:48:58||2005.02.17||Gr5/H-Br: 6003; H-Br: 1202|
|06:09:23||+20:56:58||2005.02.17||Gr5/H-Br: 6003; H-Br: 1202|
|06:10:24||+20:15:44||2005.02.17||Gr5/H-Br: 6003; H-Br: 1202|
|06:09:23||+20:16:52||2005.02.17||Gr5/H-Br: 6003; H-Br: 1202|
|06:10:54||+20:36:23||2006.01.09||Gr5/H-Br: 6003; H-Br: 1502|
|06:08:43||+20:40:47||2007.01.25||Gr5/H-Br: 6003; H-Br: 1202|
|06:08:43||+20:32:52||2007.01.25||Gr5/H-Br: 6003; H-Br: 1202|
|06:08:43||+20:24:53||2007.01.25||Gr5/H-Br: 6003; H-Br: 1202|
|06:08:05||+20:32:57||2007.01.25||Gr5/H-Br: 4503; H-Br: 902|
|06:08:05||+20:40:56||2007.01.25||Gr5/H-Br: 4503; H-Br: 902|
|06:08:43||+20:16:58||2007.01.26||Gr5/H-Br: 4503; H-Br: 902|
|06:08:05||+20:25:01||2007.01.26||Gr5/H-Br: 4503; H-Br: 902|
|06:08:05||+20:33:01||2007.01.26||Gr5/H-Br: 4503; H-Br: 902|
|06:09:46||+20:38:44||2009.11.15||Gr5/H-Br: 6003; H-Br: 902|
|06:09:05||+20:38:44||2009.11.15||Gr5/H-Br: 6003; H-Br: 902|
|06:09:46||+20:28:53||2009.11.15||Gr5/H-Br: 6003; H-Br: 902|
|06:08:31||+20+53:16||2009.11.15||Gr5/H-Br: 6003; H-Br: 902|
|06:10:14||+20+52:54||2009.11.15||Gr5/H-Br: 6003; H-Br: 902|
105-cm Schmidt Telescope, Kiso, Japan
104-cm Sampurnanand Telescope, ARIES, Naini Tal
2-m Himalayan Chandra Telescope, IAO, Hanle
The CCD photometry of the central region ( = ; = ) of Sh2-252 was carried out using the 104-cm Sampurnand Telescope (ST) of ARIES, Nainital, India on 2005 December 2. The 2K 2K CCD with a plate scale of 0.37 arcsec pixel covers a field of arcmin in the sky. To improve the signal to noise ratio, the observations were carried out in a binning mode of pixels. The standard field SA 98 from Landolt (1992) was observed on the same night to apply the atmospheric and instrument corrections to the target field. We then performed PSF photometry with the DAOPHOT package in IRAF, and derived the extinction coefficients and colour coefficients using the observations of standard field SA 98. The central region of Sh2-252 observed with ST was calibrated by applying these coefficients. The calibration uncertainties between the standard and transformed magnitudes and , , , colours were 0.03 mag. The selected isolated stars of this central calibrated region were used to calibrate the wide field observations taken with Kiso Schmidt Telescope.
We further carried out deep observations in and bands using ST for 20 overlapping sub-regions of Sh2-252, to cover an area 1 degree 1 degree around it. The observations were conducted from 2006 October to 2008 December in good photometric conditions and the log of observations is tabulated in Table 1. Each region was exposed for 60 min in and 30 min in -band. The combined area of these 20 overlapping regions is shown in Fig. 1 using solid lines. The secondary standards from the Kiso observations were used to calibrate these individual regions. The final catalog is made from these three sets of observations, where, the magnitudes of the bright ( 15 mag) and faint ( 15 mag) objects were taken from Kiso observations and ST observations, respectively. To ensure good photometric quality, we selected only those sources having uncertainty 0.2 mag in -band to make the final catalog. Thus we obtained photometry of 8791 sources detected at least in and -bands with a limiting magnitude of 23 mag.
In order to check the photometry accuracy, we compared our photometry with the photometries available in the literature by Grasdalen & Carrasco (1975) and Haikala (1994). The mean difference between our photometry and theirs were 0.03 mag in and , which shows that our photometry is in agreement with the previous studies.
2.2 slitless grism spectroscopy
An emission line survey of Sh2-252 was conducted between 2005 January and 2009 November. The observations were done in the slitless mode with a grism as the dispersing element using the HFOSC of Himalayan Chandra Telescope (HCT) of IAO, Hanle, India. The central 2K 2K pixels of the 2K 4K CCD were used for data acquisition. The CCD with an image scale of 0.296 arcsec pixel covers an area of 10 10 arcmin in the sky. A combination of broad-band filter (-Br; 6100 - 6740 Å) and Grism 5 were used. The resolution of Grism 5 is 870. We observed 25 overlapping sub-regions around Sh2-252 and the area covered for this observations is shown in Fig. 1 using a dot-dashed box. Multiple frames were taken to ensure the presence of emitting sources. The log of the observations are given in Table 1. emission line stars with an enhancement over the continuum at the line were visually identified. We obtained 61 emitting sources within our surveyed area. The detection limit of the H-alpha survey is about 3 Å in terms of equivalent width, or V 22 in terms of magnitude. In Table 2, we have provided the J2000 coordinates, -band magnitudes and colours for these emission line sources, if available.
2.3 Slit spectroscopy
Low resolution optical spectroscopic observations of fifteen optically bright sources in Sh2-252 were made using the HFOSC of HCT. These sources are the probable massive members of Sh2-252 and some of them are the candidate ionizing sources of the sub-H ii regions (see section 3.4.1). The locations of these sources are marked in Fig. 2. The observations were taken on 16 and 17 November 2009. The spectra in the wavelength range 3800-6840 Å with a dispersion of 1.45 Å pixel were obtained using the low resolution Grism 7 with a slit of 2 arcsec width and exposure time of 900s. One-dimensional spectra were extracted from the bias-subtracted and flat-field corrected images by using the optimal extraction method in IRAF. Wavelength calibration of the spectra was done by using an FeAr lamp source. Spectro photometric standard star (G191B2B) was observed on these two nights to give the flux calibration to the target spectra.
2.4 Near-infrared data from 2MASS
NIR data for point sources within a radius of 30 arcmin around Sh2-252 have been obtained from 2MASS point source catalog (PSC) (Cutri et al. 2003). To improve photometric accuracy, we used photometric quality flag (phqual = AAA) which gives a S/N 10 and a photometric uncertainty 0.10 mag. This selection criterion ensures best quality detection in terms of photometry and astrometry (cf. Lee et al. 2005). The data were transformed from the 2MASS system to the California Institute of Technology (CIT) system by using the relations given by Carpenter (2001).
2.5 Completeness of the data
To study the luminosity function (LF) / mass function (MF), it is necessary to take into account the incompleteness of the data that could occur due to various factors (e.g., crowding of the stars). We used ADDSTAR routine of DAOPHOT to determine the completeness factor (CF). The procedure has been outlined in detail in our earlier works (see e.g., Pandey et al. 2001). Briefly, we randomly added artificial stars to both and images taken with ST in such a way that they have similar geometrical locations but differ in brightness according to the mean colour ( mag) of the data sample. The luminosity distribution of artificial stars was chosen in such a way that more number of stars were inserted towards the fainter magnitude bins. The frames were reduced using the same procedure used for the original frames. The ratio of the number of stars recovered to those added in each magnitude interval gives the CF as a function of magnitude. The minimum value of the CF of the pair (i.e., - and -bands) was used to correct the data incompleteness. It is to be noted that the spatially varying nebular background (see. Fig. 2) and stellar crowding characteristics (see Fig. 9) could affect the local completeness limit. In order to account for this, we estimated the CFs of two representative sub-regions Sh2-252C and NGC 2175s and the values obtained are given in Table 3. The CFs of both the regions seem to be agreeing well with each other except for the faint magnitude bin, where the variation seems to be significant.
We also estimated the CFs of the 2MASS data using the -band images from archive around the regions Sh2-252C and NGC 2175s. Using the DAOPHOT ALLSTAR package, we detected all the stars in the 2MASS PSC, with photometry accuracy better that 0.05 mag. We then performed the completeness analysis using ADDSTAR routine as mentioned above. The CFs thus obtained for the two sub-regions are given in Table 3. The CFs of both the regions seem to be more or less similar in all the magnitude bins.
|mag. range||Sh2-252C||NGC 2175s||Sh2-252C||NGC 2175s|
|10 - 11||1.00||1.00||1.00||1.00|
|11 - 12||1.00||1.00||0.99||0.99|
|12 - 13||1.00||1.00||0.97||0.97|
|13 - 14||1.00||1.00||0.92||0.95|
|14 - 15||1.00||1.00||0.85||0.88|
|15 - 16||1.00||1.00||0.67||0.70|
|16 - 17||1.00||1.00||-||-|
|17 - 18||0.99||0.99||-||-|
|18 - 19||0.96||0.98||-||-|
|19 - 20||0.89||0.92||-||-|
|20 - 21||0.80||0.85||-||-|
|21 - 22||0.59||0.71||-||-|
3.1 Spectral classification of the optically bright sources in Sh2-252
The targets for the low resolution spectroscopy were selected on the basis of their brightness, their location within 2 arcmin radius of the sub-regions of Sh2-252 (see Fig. 2) as well as their location on the optical colour-magnitude diagram (CMD; Fig. 7). The coordinates of these sources are given in Table 4. To determine the spectral type, we extracted the low-resolution, one dimensional spectrum of each source. The flux calibrated, normalized spectrum of these sources are shown in Fig. 3.
In order to classify the stars earlier than B3, we used the criteria given by Walborn & Fitzpatrick (1990). For later type stars we used the criteria by Jacoby et al. (1984) and Torres-Dodgen & Weaver (1993). Spectra of O and B stars have the features of hydrogen, helium and other atomic lines (e.g., O ii, C iii, Si iii, Si iv and Mg ii etc.). Hydrogen and helium lines are usually seen in absorption for dwarfs whereas they may be in emission in super giants. In the case of early type stars, the ratio of He i 4471/He ii 4542 is a primary indicator of the spectral type and the ratio is greater than 1 for spectral type later than O7. The line strength of He ii gets weaker for late O type stars and He ii 4686 is last seen in B0.5 type stars (Walborn & Fitzpatrick 1990). If the spectrum displays He ii line at 4200 Å, along with the O ii/C iii blend at 4650 Å, it indicates that the spectral type of the star must be earlier than B1. The absence of He ii 4200, He ii 4686 and Mg ii 4481 and the weak features of silicon lines along with the weak O ii/C iii blends at 4070 and 4650 are the supporting criteria in the B1-B2 range. For type B2, He i is in its maximum and for later types Si ii 4128-4130 and Mg ii 4481 appears stronger (Walborn & Fitzpatrick 1990). The presence of the He i lines in absorption constrains the spectral types to be earlier than B5-B7. Presence of the spectral lines Mg ii 4481 and Si iii 4552 is an indication of evolved early B type stars. The spectral lines such as Na i 5893, He i, H, Ca i 6122, 6162, Fe ii 6456 are used to classify the late type sources (Torres-Dodgen & Weaver 1993). Finally, the spectral type was assigned to each source by visual comparison to the standard library spectra given in the literature (Jacoby et al. 1984) and our classification is given in Table 4. However, because of low resolution observations an uncertainty of 0.5 to 1 in the sub-class identification is expected.
Eight out of the fifteen stars observed are found to have spectral types earlier than B3V. The sub-H ii regions of Sh2-252 (i.e., A, B, C and E) are found to have at least one star of spectral type early than B3V (i.e., star 1, 3, 5 and 11). It is quite interesting to note that the remaining four early type sources (i.e., star 6, 7, 8 and 9) belong to the small cluster NGC 2175s. The star 2 shows H 6563 line in strong emission (not shown in Fig. 3). This star also shows strong NIR excess emission in the colour-colour (C-C) diagram (see Fig. 5) most likely originating from the extensive, hot circumstellar material around the star. The emission line spectra as well as the NIR excess emission are the primary indicators of an young stellar object (YSO) and hence we tentatively put this object in the Herbig Be category.
The spectral types of some of these sources have already been mentioned in the literature. On the basis of the 1280 MHz radio continuum flux, Tej et al. (2006) predicted the spectral type of the ionizing source of Sh2-252A as B0-B0.5V. The star 3, located within 2 arcmin of the radio peak of Sh2-252A, is the brightest source in the optical and NIR wavelengths. This could be the ionizing source of Sh2-252A. We classify this source as an O9.5V type star, which is in agreement with the classification by Tej et al. (2006). The star 5, the probable exciting source of Sh2-252E, had been assigned a spectral type between O8 to B0.5V by Garnier & Lortet-Zuckermann (1971) and B1V by Grasdalen & Carrasco (1975). Our classification (B1V) is in good agrrement with that given by Grasdalen & Carrasco (1975). The star 11, the brightest source of Sh2-252B was classified as B2V by Grasdalen & Carrasco (1975). However, our classification (B0V) leads to a class slightly earlier than their study. The star 2 had been identified as a probable Herbig Be source of spectral type B9 by Chavarría et al. (1989), which seems to match with our classification. Thus within errors the present spectral classifications are in accordance with those reported in the literature.
For 16 data is taken from Pismis (1977)
3.2 Reddening, membership and distance of the bright sources in Sh2-252
The extinction in a star forming region arises due to two distinct sources: (1) the general ISM in the foreground of the region , and (2) the localized cloud associated with the region [ = ], where is the observed reddening of the star embedded in the local cloud. The presence of the molecular cloud (Lada & Wooden 1979) and the nebulosity associated with Sh2-252 as evident in Fig. 2 would give rise to variable extinction within this region. To estimate the extinction towards the region, we used the C-C diagram. Fig. 4 shows the C-C distribution of all the stars detected in our Kiso observations, where the spectroscopically observed sources are numbered in red. It is to be noted that some of the spectroscopically observed stars are not covered/detected in the and bands, hence they are not included in Fig. 4. The star 16, which is the ionizing source of Sh2-252 (i.e., HD 42088), was saturated in our observations. Hence the data for this star is taken from Pismis (1977). The zero-age-main-sequence (ZAMS) by Girardi et al. (2002) that is reddened by = 0.23 mag along the reddening vector having a normal slope of = 0.72 seems to match with the first sequence of distribution in Fig. 4. However, majority of the spectroscopically identified stars earlier than spectral type B3 show in the range 0.35 - 0.75 mag. The stars 10, 13 and 15 are found to be late type stars and hence may not be associated with Sh2-252. Assuming that uncertainties in the and colours mainly originate from their calibration (ref. section 2.1), the typical error in the estimation of would be 0.04 mag. A careful inspection of the C-C diagram indicates the presence of further reddened population at 0.65 - 0.90 mag. Some of them may be the reddened members of the Sh2-252 region. The population having 1.00 - 1.25 mag could be part of the background Norma-Cygnus arm as discussed by Carraro et al. (2005) and Pandey et al. (2006).
We also estimated the individual reddening of spectroscopically classified stars. We used the relation = , where is the observed colour and is the intrinsic colour taken from Schmidt-Kaler (1982). Some of the stars do not have the data (i.e., stars 3, 6, 7 and 14) and for these stars was calculated by using the relations = and = /0.11 (Cohen et al. 1981), where, colours were taken from 2MASS and intrinsic colours were taken from Koornneef (1983). The thus calculated are given in Table 4. The range of reddening of the spectroscopically identified stars shown in Fig. 4 agree well with that given in Table 4.
In Fig. 5 we have shown the C-C distribution of all the spectroscopically identified stars in Sh2-252. The sources falling towards the right of the reddening vector are considered to have NIR excess (see sect. 3.4.3 for details). The star 2 falls in the location where Herbig Ae/Be stars are generally found (Hernandez et al. 2005) and this justifies our spectral classification of this source as a Herbig Be type (see sect. 3.1). The figure also shows that star is highly reddened when compared to other massive stars in the field.
The estimated spectral classification, apparent magnitudes and extinction of the massive members of Sh2-252 (Table 4) enable us to derive their distances. Using the present spectral classification along with and values from Table 4 and values from the table of Schmidt-Kaler (1982), we estimated the intrinsic distance modulus () of each source and the same is given in Table 4. The average value of the intrinsic distance modulus is found to be 11.91 0.18 mag, which corresponds to a distance of pc. We also calculated the intrinsic distance modulus () of each source based on the observed colours from 2MASS data and intrinsic colours and absolute magnitudes from Koornneef (1983) and Schmidt-Kaler (1982), respectively. The average value of the intrinsic distance modulus using 2MASS data was found to be 11.80 0.10 mag, which is comparable with that of the intrinsic distance modulus estimated from optical data within errors. Fig. 6 shows dereddened CMD for the probable main sequence (MS) members of Sh2-252 identified spectroscopically. In this figure, we have also plotted the theoretical ZAMS locus for solar metalicity () by Girardi et al. (2002), adjusted for the distance of 2.4 kpc, which seems to be matching well with the distribution of the probable MS members. The present distance estimate is in good agreement with that of Georgelin & Georgelin (1970; 2.4 kpc), Neckel et al. (1980; 2.3 kpc) and Chavarría et al. (1987; 2.3 kpc). At the same time, it is also in agreement with the distance estimates by Grasdalen & Carrasco (1975; 2.6 kpc) and Reid et al. (2009; 2.1 0.03 kpc) within errors.
An important issue to be discussed here is the association of NGC 2175s and Teu 136 with Sh2-252. NGC 2175s is a small, optically visible cluster located at the east of Sh2-252 (at = ; = ; see Figs. 1 and 2). The photoelectric study by Pismis (1970) put this cluster at a distance of 3500 pc. However, in the present spectroscopic observations, we identified four early B type members of this cluster (i.e., stars 6, 7, 8 and 9, see Table 4) and the distance estimates of these four sources are in good agreement with that of the average distance of Sh2-252. Hence, our spectroscopic analysis supports the association of NGC 2175s with the Sh2-252 complex. Similarly, the cluster Teu 136 was identified as a NIR embedded cluster by Koposov et al. (2008) in a survey of the 2MASS catalog designed to identify new open clusters. Located at coordinates = ; = , 15 NE of NGC 2175s, this cluster lies close to the eastern edge of the Sh2-252 complex. We obtained spectra of two stars located close to this cluster (i.e., stars 14 and 15) and the distance of star 14 (see Table 4) seems to match well with that of Sh2-252. Also, the spectral types of two bright stars at the center of this cluster, which are saturated in our photometric observations are given as B1.5 and B2.5, respectively in Reed (2003). The average distance modulus of these two stars comes out to be 11.74 mag which seems to match with the average distance modulus of the Sh2-252 complex. Therefore we presume that Teu 136 could be associated with the Sh2-252 complex, too. Hence, in this paper we report that the NIR cluster Teu 136 as a sub-cluster of the Sh2-252 complex.
In Fig. 7, we show the CMD for all the stars detected in our optical photometry. Spectroscopically identified massive MS members are shown with large filled circles in red, where the brightest star is HD 42088. The dashed curves are the ZAMS by Girardi et al. (2002) shifted for the distance of 2.4 kpc and reddening =0.35 and =0.75, respectively. The average photometric uncertainty in magnitudes and colours estimated from the difference between the input and recovered magnitudes of the artificial stars (see Sect. 2.5) are shown in the left side of the figure. The CMD clearly shows a broad MS band down to 14 mag with a significant number of sources located within the reddening strips (i.e., = 0.35 to 0.75 mag). In the absence of spectroscopic information, it is difficult to assign membership for these bright sources. Since Sh2-252 is located close to the Galactic plane (l=190.04; b=+0.48) at a distance of 2.4 kpc, it is quite natural that the region could be significantly contaminated by Galactic field MS and giant populations. The extinction variation inside the Sh2-252 region may cause a spread in the MS both for the member sources as well as for the background sources. A significant number of sources can also be noticed towards the right side of the ZAMS, part of which could be PMS stars.
3.3 Age and mass distribution of candidate YSOs in Sh2-252
In the absence of proper motion data, the membership of Sh2-252 can be established by identifying stars with indicators of youth. YSOs emit excess radiation in the infrared in comparison to MS stars due to the thermal emission from their circumstellar material. So YSOs can be identified on the basis of their infrared excess emission. We have also carried out a YSO survey of the Sh2-252 region using the 2MASS and the Spitzer-IRAC data sets (3.6, 4.5, 5.8 and 8.0 m) and the detailed analyses will be presented in a forthcoming paper. A total of 96 Class I and 400 Class II YSOs are identified within this complex on the basis of IR colour excess criteria given by Gutermuth et al. (2009). The possible contamination such as PAH-emitting galaxies, broad-line AGNs etc. have been excluded using the criteria given by Gutermuth et al. (2009). Similarly, the presence of H emission line is considered as a significant characteristic of a YSO with ongoing disk accretion process (e.g., Dahm 2005). We identified 61 H emission line stars in our slitless spectroscopy survey of the Sh2-252 region (see section 2.2). Our deep photometry provided and counterparts for 15 Class I, 179 Class II and 54 H emission line sources. Of these 54 H sources 4 are Class I and 33 are Class II YSOs. The distribution of these candidate YSOs i.e., Class I, Class II and H emission line sources on the CMD is shown in Fig. 8 by using red squares, green triangles and blue triangles, respectively. The distance to Sh2-252 is adopted as 2.4 kpc as estimated on the basis of the spectroscopic analyses mentioned in section 3.2. In Fig. 8 the ZAMS (thick solid curve) by Girardi et al. (2002) and the PMS isochrones (dashed curves) by Siess et al. (2000) for age 0.1 and 5 Myr are also shown. The ZAMS and isochrones are shifted for the distance of 2.4 kpc and average reddening i.e., = 0.5 mag. The average photometric uncertainty in magnitudes and colours estimated from the completeness simulation (see Sect. 2.5) close to region C are shown in the right side of the figure. The distribution of YSOs on the CMD is an ideal tool to estimate the approximate ages of YSOs. It is evident from this figure that a majority of the YSOs in Sh2-252 are distributed between the PMS isochrones of age 0.1 and 5 Myr, which shows that the region Sh2-252 is significantly comprised of young sources. However, the photometric uncertainty, differential reddening, binarity, different evolutionary stages of YSOs etc. can cause a spread in the CMD. The spectroscopic analysis of the member stars shows that the reddening in the region varies between 0.9 to 6.5 mag (see Table 4). Also, the NIR colour-colour analysis of the candidate PMS sources shows that they are reddened up to 8 mag (see Sect. 3.4.3), which suggest that the region is significantly affected by the differential reddening. Similarly, the photometric uncertainties in the case of YSOs (see Fig. 8) may be relatively large due to crowding and variable background which may cause the scatter in the distribution of the Class I and Class II sources in the CMD. Hence the differential reddening and photometric uncertainty would have an impact on the observed age spread in the CMD. The detailed spatial distribution and evolutionary status of the YSOs in the Sh2-252 region will be analyzed in the forthcoming paper. In Fig. 8 we have also shown the PMS evolutionary tracks by Siess et al. (2000) (thin solid curves) for various mass bins (the value of mass for each track is given towards its right), which indicate that a majority of the YSOs have masses in the range between 0.3 - 2.5 .
The ages of the young clusters/H ii regions can also be derived from the post-main-sequence age of their most massive members. The most massive member of Sh2-252 is HD 42088, an O6.5 star which is still in the MS (Conti & Alschular 1971; Walborn 1972) and is considered to be the main ionizing source of Sh2-252. Also, none of the spectroscopically identified OB stars of Sh2-252 are found to be evolved (see section 3.1). Hence, the age of the H ii region Sh2-252 should be younger or of the order of the MS life time of the O6.5V star. i.e., 4.4 Myr (Meynet et al. 1994). Thus we can put an upper limit to the age of Sh2-252 as 4 Myr. However, this has to be taken as an approximate estimate as the low mass members of the region could have been formed prior to HD 42088.
3.4 Stellar contents of the sub-regions in Sh2-252
3.4.1 Ionizing sources of the CH ii regions
As discussed in section 1.1, the radio observation of this region by Felli et al. (1977) identified six compact thermal radio sources, of which the regions A, B, C and E have been claimed as CH ii regions having its own ionizing sources in each of them. In order to search for the ionizing sources of these regions, we obtained spectra of optically bright sources within 2 radius around each region. The spectral types of these bright sources have been discussed in section 3.1 and each region is found to have at least one source with spectral type earlier than B3V. Their spectral type can also be obtained by measuring the number of Lyman continuum photons (N ) emitted per second by the star that would be responsible for the ionization of each region. From the radio continuum flux at 1415 MHz by Felli et al. (1977) together with our estimated distance of 2.4 kpc and assuming an electron temperature of 10 K, we calculated log N using the relation given by Martín-Hernández et al. (2003). The log N have been obtained as 46.73, 47.52, 46.57 and 46.55, respectively for regions A, B, C and E. From these log N values we derived the spectral type of the ionizing sources (Panagia 1973) as B0-B0.5, B0, B0.5 and B0.5, respectively, for regions A, B, C and E. These spectral types are consistent with the earliest spectral types derived by our optical observations for regions A, B and E (i.e., star 3, 11 and 5, see Table 4) within the errors. The 1280 MHz observations by Tej et al. (2006) also estimated the spectral type of the ionizing source of region A as B0-B0.5V which is consistent with our estimate. Slight difference in the spectral type estimate is expected due to uncertainties in the effective temperature and ionizing flux of massive stars, for a given spectral type. Thus we conclude that the stars 3, 11 and 5 are the ionizing sources of the CH ii regions A, B and E, respectively. However in the present observation, the earliest spectral type derived in region C is B2.5V (see Table 4), which does not seem to be early enough to explain the radio continuum flux. The expected spectral type from the radio flux is much earlier than B2.5V. The study by Kömpe et al. (1989) showed that the probable spectral type of the IRAS point source (IRAS 06063+2040) located at the center of the nebulosity in region C is B0V, which could be the candidate ionizing source of this CH ii region. Since this source is deeply embedded in the cloud (=9.16, 8 mag), we could not carry out optical spectroscopy of this source and future spectroscopy would help us to confirm this result. In Fig. 9 we have shown the -band images of the sub regions A, B, C, E, NGC 2175s and Teu 136 for an area 8 8. The probable ionizing sources of the regions A, B, C and E are marked with a white star symbol.
3.4.2 Stellar surface density distribution in Sh2-252
Though we have the optical and NIR observations for an area 1 degree 1 degree around the Sh2-252 region, however a majority of the members of Sh2-252 might be confined within its sub-regions (see section 1.1). In order to define the extent of each sub-region and to study the morphology of the region, we generated isodensity contours for stellar population detected in 2MASS -band. The two dimensional stellar surface density distribution of the Sh2-252 region manifests prominent stellar density enhancements at the locations of A, C, E, NGC 2175s and Teu 136. This indicates the presence of clustering associated with these sub-regions. On the other hand, no clustering is apparent towards the sub-regions B and F. In Fig. 9 the stellar surface density contours above 3 of the background level are over plotted on the -band images of regions A, B (no contours although), C, E, NGC 2175s and Teu 136 . We will discuss the properties of these five sub-regions i.e., A, C, E, NGC 2175s and Teu 136 individually in ensuing sections. The extent of these sub-regions is estimated on the basis of the surface density distribution. The contour above the 3 background level around each region is considered as its extent. Thus we estimated the radius as 3.5, 3.5, 3.0, 3.0, 3.0, respectively, for the region A, C, E, NGC 2175s and Teu 136. We used these extents for further analysis of these sub-regions.
3.4.3 NIR colour-colour diagrams
NIR C-C diagrams are ideal tools to identify the candidate PMS sources having NIR excess (Hunter et al. 1995; Haisch et al. 2000; 2001; Sugitani et al. 2002; Devine et al. 2008; Chavarria et al. 2010). The / C-C diagrams for the sub-regions A, C, E, NGC 2175s and Teu 136 and a nearby control field are shown in Fig. 10. We have considered an area of 3.0 radius centered at = , = as the control field. This region is devoid of nebulosity in the -band (see Fig. 1), hence can be considered as a field region. The thin and thick solid curves are the locations of unreddened MS and giant stars (Bessell Brett 1988), respectively. The dotted lines represent the locus of intrinsic and reddened ( = 4.0 mag; 8.0 mag) classical T Tauri stars (CTTSs; Meyer et al. 1997). The two long parallel dashed lines are the reddening vectors for the early MS and giant type stars (drawn from the base and tip of the two branches). One more reddening vector is plotted from the tip of the unreddened CTTS locus. The crosses on the reddening vectors are separated by an value of 5 mag. The extinction ratios, and , are adopted from Cohen et al. (1981). The magnitudes, colours of the stars and the curves are in the CIT system.
Following Ojha et al. (2004a), we classified sources according to their locations in C-C diagrams. The ‘F’ sources are those located between the reddening vectors projected from the intrinsic colours of MS and giant stars. These sources are reddened field stars (MS and giants) or Class III/Class II sources with little or no NIR excess (viz., weak-lined T Tauri sources (WTTSs) but some CTTSs may also be included). The sources located redward of region ‘F’ are considered to have NIR excess. Among these, the ‘T’ sources are located redward of ‘F’ but blueward of the reddening line projected from the red end of the CTTS locus. These sources are considered to be mostly CTTSs (Class II objects) with large NIR excesses (Lada & Adams 1992). There may be an overlap in NIR colours of Herbig Ae/Be stars and T Tauri stars in the ‘T’ region (Hillenbrand et al. 1992). The ‘P’ sources are those located in the region redward of region ‘T’ and are most likely Class I objects (protostellar-like) showing large amount of NIR excess. Here it is worthwhile to mention that Robitaille et al. (2006) have shown that there is a significant overlap between protostellar-like objects and CTTSs in the C-C diagram.
A comparison of the distribution of sources in the sub-regions and the control field in Fig. 10 suggests that there is an appreciable difference between them. A significant fraction of sources within the sub-regions are concentrated between the intrinsic and reddened CTTS locus (i.e., within the ‘F’ and ‘T’ regions), whereas a majority of the sources in the control field are mainly concentrated towards the left of the ‘F’ region (i.e., 0.35 mag). Average value of the uncertainty in colour is 0.05 mag. Taking this uncertainty into consideration and on the basis of the comparison of the C-C diagrams, we can safely assume that a majority of the sources in the sub-regions located between the intrinsic and reddened CTTS locus with 0.40 mag and lying in the ‘F’, ‘T’ and ‘P’ regions are most likely PMS members. Some of the above sources lying in the ‘F’ region could be the reddened field stars, however a majority of them are likely candidate WTTSs or CTTSs with little or no NIR excess. Fig. 10 shows that the clusters NGC 2175s and Teu 136 do not have many PMS members with 4.0 mag, whereas a significant fraction of the PMS members of the sub-regions A, C and E are reddened up to 8.0 mag. This indicates that the clusters NGC 2175s and Teu 136, which are located towards the east of Sh2-252, are less reddened in comparison to other sub-regions of Sh2-252 located towards the west.
3.4.4 Optical colour-magnitude diagrams
The CMDs for the sources within the sub-regions of Sh2-252 and the nearby control field are shown in Fig. 11. The ZAMS (thick solid curve) by Girardi et al. (2002) and PMS isochrones (dashed curves) by Siess et al. (2000) for age 0.1 and 5 Myr are also shown. The ZAMS and isochrones are shifted for the distance of 2.4 kpc and reddening. The encircled sources are those located above the intrinsic CTTS locus in the C-C diagram (Fig. 10) and with 0.4 mag. A majority of these sources are found to have ages 5 Myr in all the sub-regions, suggesting that the sources with 0.4 mag could be probable PMS members as discussed in section 3.4.3. A comparison of these CMDs with that of the control field also supports the notion that the sources lying between 0.1 and 5 Myr isochrones could be the PMS sources associated with the sub-regions. Similarly, the CMDs of the sub-regions resemble the CMD for the candidate YSOs shown in Fig. 8. It further supports that the majority of the sources identified on the basis C-C diagrams (Fig. 10) in the sub-regions could be the PMS members and that they have an age spread of 5 Myr. A detailed analysis on the evolutionary status of the YSOs within the individual regions will be done in the forthcoming paper.
3.5 -band luminosity functions of sub-regions in Sh2-252
The luminosity function in the - band is frequently used in studies of young clusters and star forming regions as a diagnostic tool of the mass function and the star formation history of their stellar populations. Pioneering work on the interpretation of KLF was presented by Zinnecker et al. (1993). During the last decade several studies have been carried out with the aim of determining the KLF of young clusters (e.g., Muench et al. 2000; Lada & Lada 2003; Ojha et al. 2004b; Sanchawala et al. 2007; Pandey et al. 2008; Jose et al. 2008; 2011). We used the 2MASS - band data to study the KLF of the Sh2-252 region. Because this region shows clusterings around the regions A, C, E, NGC 2175s and Teu 136 (see section 3.4.2), we estimated the KLF for each region independently. The KLF is calculated within the area of each region as estimated on the basis of stellar surface density distribution discussed in section 3.4.2.
KLF is the number of stars as a function of - band magnitude. In order to convert the observed KLF to the true KLF, it is necessary to account for the data incompleteness as well as the background and foreground source contamination. The regions C and NGC 2175s were corrected for their respective CFs given in Table 3 and the average value of the CFs for regions C and NGC 2175s was applied to regions A, E and Teu 136 to correct for their data incompleteness. The data from the control field region (see section 3.4.3) has been used to remove the field star contribution. Since this control field is located off the Sh2-252 region, its background population will be affected by a smaller interstellar reddening. Whereas, the background field stars of the sub-regions of Sh2-252 are seen through a larger reddening due to the matter associated with it. Hence the population in the background of sub-regions will be more reddened as compared to the control field. Therefore a direct subtraction of the observed control field from the sub-regions of Sh2-252 will yield incorrect LF. To account for the higher extinction towards the embedded sub-regions, the field star population towards the direction of the control field is predicted on the basis of the Besançon Galactic model of stellar population synthesis (Robin et al. 2003) by using a similar procedure as described by Ojha et al. (2004b). An advantage of using this model is that we can simulate foreground (d 2.4 kpc) and the background (d 2.4 kpc) field star populations separately. The foreground population was simulated by using the model with the extinction = 1.1 mag ( = 0.35 mag; cf. section 3.2) and d 2.4 kpc. As discussed in section 3.4.3, a majority of the PMS members of regions A, C and E are reddened up to =8.0 mag, whereas in NGC 2175s and Teu 136 they are reddened up to =4.0 mag. Hence the background population in these regions are seen through clouds with extinction of 8.0 and 4.0 magnitudes, respectively. Hence, we simulated two sets of background population (d 2.4 kpc) with values of 4.0 and 8.0 mag, respectively. Then we determined the fractions of the contaminating stars (foreground background) over the total model counts. The fraction was calculated separately for both sets of the simulated background population. This fractions were used to scale the nearby observed control field and subsequently the star counts of the modified control field were subtracted from the KLFs of the sub-regions of Sh2-252 to obtain their final corrected KLFs.
The KLFs of young embedded clusters are known to follow power-law shapes (Lada et al. 1991; 1993), which is expressed as:
where is the number of stars per 1 mag bin and is the slope of the power law. We counted the number stars per unit magnitude interval after correcting for the field star contribution. We then fit a power law to the data. The resulting field star subtracted KLFs for the five sub-regions of Sh2-252 are shown in Fig. 12. The slopes of the KLFs of regions A, C, E, NGC 2175s and Teu 136 have been obtained as 0.310.03, 0.340.03, 0.250.05, 0.200.05 and 0.360.02, respectively. Within errors, the KLFs for all the five regions seem to match with each other.
KLFs of different ages are known to have different peak magnitudes and slopes (Muench et al. 2000) and hence the KLF slope could be an age indicator of young clusters. For clusters up to 10 Myr old, the KLF slope gets steeper as the cluster gets older (Ali & Depoy 1995; Lada & Lada 1995). However, there exists no precise age - KLF relationship in the literature due to huge uncertainties in their correlation (Devine et al. 2008). There are many studies on KLF of young clusters. The studies by Blum et al. (2000); Figuerêdo et al. (2002); Leistra et al. (2005; 2006) and Devine at al. (2008) indicate that the KLF slope varies from 0.2 - 0.4 for clusters younger than 5 Myr, which is in agreement with the KLFs obtained for the sub-regions of Sh2-252. The KLFs of the sub-regions of Sh2-252 are worth comparing with the recent studies of star forming regions viz; NGC 1893 (Sharma et al. 2007), Be 59 (Pandey et al. 2008), Stock 8 and NGC 1624 (Jose et al. 2008; 2011), since all the KLFs are obtained by using a similar technique. The slopes of the KLFs obtained for the sub-regions in Sh2-252 are comparable with those obtained for NGC 1893 (), Stock 8 () Be 59 () and NGC 1624 ().
3.6 Initial mass functions of sub-regions in Sh2-252
One of the most fundamental disciplines of astrophysical research is the origin of stars and stellar masses. The distribution of stellar masses that form in a star formation event in a given volume of space is called IMF and is one of the most important measurable quantities in star formation studies. Together with star formation rate, the IMF dictates the evolution and fate of star clusters and galaxies (Kroupa 2002). Young clusters and star forming regions are important objects to study the IMF since their mass functions (MFs) can be considered as IMFs as they are too young to loose a significant number of the members either by dynamical or stellar evolution.
The MF is often expressed by the power law, and the slope of the MF is given as
where is the number of stars per unit logarithmic mass interval. For the mass range , the classical value derived by Salpeter (1955) is .
One of the main objectives of this work is to measure the MFs of the sub-regions of Sh2-252. As discussed in section 3.4.4, the majority of the sources lying above the 5 Myr isochrone of the CMDs could be probable PMS members of the region and most of them have a mass range of 0.3 - 2.5M. We calculated the MFs of these PMS sources falling above 5 Myr and mass range of 0.3 - 2.5M for the five sub-regions. To remove the contamination due to field stars from the PMS sample, we statistically subtracted the contribution of field stars from the observed CMDs of sub-regions (see Fig. 11) using the procedure described in our earlier studies (e.g., Pandey et al. 2008, Jose et al. 2008; 2011). After statistically subtracting the field star contribution, we used all the candidate PMS sources lying above the 5 Myr isochrone and having mass in the range of 0.3 - 2.5M to calculate the MFs. The MFs for the candidate PMS sources were obtained by counting the number of stars in various mass bins, shown as evolutionary tracks by Siess et al. (2000) in Fig. 8. Necessary corrections for data incompleteness as a function of magnitude were taken into account to calculate the MFs. The CFs given in Table 3 for regions C and NGC 2175s were used to correct for their data incompleteness whereas the average value of the CFs for regions C and NGC 2175s was applied to regions A, E and Teu 136 to correct for their data incompleteness. The resulting field star decontaminated and completeness corrected MFs of the PMS sources of the five sub-regions of Sh2-252 region are shown in Fig. 13. The slopes, , of the MFs in the mass range are found to be -1.330.51, -2.390.15, -1.400.30, -1.650.28, -1.820.26, respectively, for regions A, C, E, NGC 2175s and Teu 136. The slopes of the MFs of the regions A and E are in good agreement with that of the Salpeter (1955) value, and in view of large errors, those in the case of NGC 2175s and Teu 136 can also be considered comparable to the Salpeter value. However, in the case of the region C, the MF seems to be steeper than the Salpeter value. In order to check this we also calculated the MF using -band luminosity function (JLF). The procedure for this calculation has been mentioned in our earlier paper (Jose et al. 2011). The MF thus obtained for region C from JLF is shown in Fig. 13 as dashed line and the slope obtained is -2.100.54. Considering the large uncertainty in the IMF calculation, we conclude that all the regions have IMF slopes, which do not significantly deviate from the salpeter value.
The shape of the stellar IMF and whether it is universal or not are key issues in astrophysics. For clusters within 2 kpc, there is no compelling evidence for variations in the stellar IMF with respect to their environment (e.g. Meyer et al. 2000; Kroupa 2002; Chabrier 2005). We have been pursuing studies of various star forming regions in the galaxy, hence a comparative study of their IMFs obtained using similar techniques will give useful information on IMFs. Our recent analysis of young clusters (age 2-4 Myr), viz., NGC 1893 (Sharma et al. 2007), Stock 8, NGC 1624 (Jose et al. 2008, 2011) and W5 E (Chauhan et al. 2011) have yielded values similar to the Salpeter value. In the present study, the MFs of all the sub-regions of Sh2-252 are found to be in general comparable to the Salpeter value.
In this paper we studied the stellar contents of the H ii region Sh2-252 which itself contains four CH ii regions such as A, B, C and E and two clusters NGC 2175s and Teu 136. We used the deep optical photometry in bands, slit and slitless spectroscopic observations along with the data from 2MASS for our analyses. An attempt has been made to identify and classify the massive members, to determine the fundamental parameters of Sh2-252 as well as to obtain the age and mass distribution of candidate PMS sources and finally to constrain the KLFs and IMFs of the sub-regions of Sh2-252.
We have carried out optical spectroscopy of fifteen bright sources of the region, out of which, eight have been identified as massive members of spectral type earlier than B3. The CH ii regions A, B, C and E are found to have at least one star of spectral type earlier than B3, whereas, the small optically bright cluster NGC 2175s is found to have four stars of spectral type earlier than B3 around its center. We have also identified the probable candidate ionizing sources of the CH ii regions A, B, C and E. From the spectro-photometric analysis, we derived the average distance of the region as 2.4 0.2 kpc. The distance estimated from the four bright stars of NGC 2175s is found to be in good agreement with the average distance of Sh2-252, which supports the association of NGC 2175s with the Sh2-252 complex. In this paper, we also report that the NIR embedded cluster Teu 136, located at the eastern edge of the complex as a sub-cluster of Sh2-252. Based on the optical colour-colour diagram and spectroscopic properties, the reddening of the massive members of the region is found to vary between 0.35 mag to 2.1 mag.
Using slitless spectroscopy survey, we have identified 61 emission sources in the region. We obtained optical photometry for 211 candidate YSOs of the region. The distribution of these candidate YSOs on the CMD shows that a majority of them have an age distribution between 0.1 - 5 Myr and masses in the range of 0.3 - 2.5 M. The stellar surface density distribution in -band shows that there is at least five sub-clusters, each associated with the regions A, C, E, NGC 2175s and Teu 136. The C-C diagrams of the individual regions show that the regions A, C and E have more reddened population when compared to the clusters NGC 2175s and Teu 136, which are located towards the east of Sh2-252. The optical CMDs of the individual regions also show that the candidate PMS sources are distributed between 0.1 - 5 Myr, suggesting an age spread for them.
We also calculated the KLFs of the sub-regions A, C, E, NGC 2175s and Teu 136 and have estimated as 0.310.03, 0.340.03, 0.250.05, 0.200.05 and 0.360.02, respectively. Within errors, the KLFs of all the sub-regions are similar and comparable to that of young clusters of age 5 Myr. We also derived the MFs of PMS sample of the individual regions in the mass range of . The slopes of the MFs of all the sub-regions are found to match with the Salpeter value.
Authors are thankful to the anonymous referee for the useful comments which has improved contents and presentation of the paper significantly. We thank the staff of the Kiso Observatory, Japan, ARIES, Naini Tal, India and IAO, Hanle and its remote control station at CREST, Hosakote for their assistance during the observations. This publication also makes use of data from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. JJ and NC are thankful for the financial support for this study through a stipend from CSIR and DST, India. AKP is thankful to DST (India) and JSPS (Japan) for providing funds to visit KISO observatory to carry out the observations.
Ali B., Depoy D. L., 1995, AJ, 109, 709
Bessell M., Brett J. M., 1988, PASP, 100, 1134
Blum R. D., Conti P. S., Damineli A., 2000, AJ, 119, 1860
Bonatto C., Bica E., 2011, MNRAS, 414, 3769
Carpenter J. M., 2001, AJ, 121, 2851
Carraro G., Vazquez R. A., Moitinho A., Baume G., 2005, ApJ, 630, L153
Chabrier G., 2005, The Initial Mass Function 50 Years Later, 327, 41
Chauhan N., Pandey A. K., Ogura K., Jose J., Ojha D. K., Samal M. R., Mito H., 2011, MNRAS, 415, 1202
Conti P.S., Alschuler W. R., 1971, ApJ, 170, 325
Cohen J. G., Frogel J. A., Persson S. E., Ellias J. H., 1981, ApJ, 249, 481
Cutri R. M., Skrutskie M. F., van Dyk S. et al., 2003, The IRSA 2MASS All Sky Point Source Catalog, NASA/IPAC Infrared Science Archive, http://irsa.ipac.caltech.edu/applications/Gator/
Chavarría-K. C., de Lara E., Hasse I., 1987, A&A, 171, 216
Chavarría-K. C. et al., 1989, A&A, 215, 51
Chavarria L., Mardones D., Garay G., Escala A., Bronfman L., Lizano S., 2010, ApJ, 710, 583
Churchwell E., 1974, in proceedings of the Second European Regional Meeting in Astron. Mem. Soc. Astron. Ital. 45, 259
Clarke C. J., 2007, MNRAS, 376, 1350
Dahm S., 2005, AJ, 130, 1805
Devine K. E., Churchwell E. B., Indebetouw R., Watson C., Crawford S. M., 2008, AJ, 135, 2095
Elmegreen B.G., Lada C.J., 1977, ApJ, 214, 725
Felli M., Habing H. J., Israël F.P., 1977, A&A, 59, 43
Figuerêdo E., Blum R. D., Damineli A., Conti P. S., 2002, AJ, 124, 2739
Garnier R., Lortet-Zuckermann M. C., 1971, A&A, 31, 41
Georgelin Y. P., Georgelin Y. M., 1970, A&A, 6, 349
Girardi L., Bertelli G., Bressan A., Chiosi C., Groenewegen M. A. T. et al., 2002, A&A, 391, 195
Grasdalen G. L., Carrasco L., 1975, A&A, 43, 259
Gutermuth R. A., Megeath S. T., Myers P. C., Allen L. E., Pipher J. L., Fazio G. G., 2009, ApJS, 184, 18
Haikala L. K., 1994, A&A, 108, 643
Hernandez J., Calvet N., Hartmann L., Briceno C., Sicilia-Aguilar A., Berlind P., 2005, AJ, 129, 856
Haisch K. E., Lada E. A., Lada C. J., 2000, AJ, 120, 1396
Haisch K. E., Lada E. A., Lada C. J., 2001, AJ, 121, 2065
Hillenbrand L. A., Strom S. E., Vrba F. J., Keene J., 1992, ApJ, 397, 613
Hunter T. R., Testi L., Taylor G. B., Tofani G., Felli M., Phillips T. G., 1995, A&A, 302, 249
Jacoby G. H., Hunter D. A., Christian C. A., 1984, ApJS, 56, 257
Jose, J. et al., 2008, MNRAS, 384, 1675
Jose, J. et al., 2011, MNRAS, 411, 2530
Kömpe C., Joncas G., Baudry A., Wouterloot J. G. A., 1989, A&A, 221, 295
Koposov S.E., Glushkova E.V., Zolotukhin I.Yu., 2008, A&A, 486, 771
Koornneef J., 1983, A&A, 128, 84
Kroupa P., 2002, SCIENCE, 295, 82
Lada C. J., Wooden D., 1979, ApJ, 232, 158
Lada C. J., Lada E. A., 1991, in ASP Conf. Ser. 13, The Formation and Evolution of Star Clusters, ed. K. Janes (San Francisco: ASP), 3
Lada C. J., Adams F. C., 1992, ApJ, 393, 278
Lada C. J., Young T., Greene T., 1993, ApJ, 408, 471
Lada E . A ., Lada C . J., 1995, AJ , 109, 1682
Lada C. J., Lada E. A., 2003, ARA&A, 41, 57
Landolt A.U., 1992, AJ, 104, 340
Lee H. T., Chen W. P., Zhang Z. W., Hu J.Y., 2005, ApJ, 624,808
Leistra A., Cotera A. S., Leibert J., Burton M., 2005, AJ, 130, 1719
Leistra A., Cotera A. S., Liebert J., 2006, AJ, 131, 2571
Martín-Hernández N. L., van der Hulst J. M., Tielens A. G. G. M., 2003, A&A, 407, 957
Meyer M., Calvet N., Hillenbrand L. A., 1997, AJ, 114, 288
Meyer M. R., Adams F. C., Hillenbrand L. A., Carpenter J. M., Larson R. B., 2000, Protostars and Planets IV, 121
Meynet G., Maeder A., Schaller G., Schaerer D., Charbonnel C., 1994, A&AS, 103, 97
Muench A. A., Lada E.A., Lada C.J., 2000, ApJ, 553, 338
Neckel T., Klare G., Sarcander M., 1980, Bull. Inform. Centre, des Donees Stellaires, 19
Ojha D. K., Tamura M., Nakijama Y. et al., 2004a, ApJ, 608, 797
Ojha D. K., Tamura M., Nakajima Y. et al., 2004b, ApJ, 616, 1042
Panagia, N., 1973, AJ, 78, 929
Pandey A.K., Nilakshi, Ogura K., Sagar R., Tarusawa K., 2001, A&A, 374, 504
Pandey A. K., Sharma S., Ogura K., 2006, MNRAS, 373, 255
Pandey A. K., Sharma S., Ogura K., Ojha D. K., Chen W. P. et al., 2008, MNRAS, 383, 1241
Pismis P., 1970, Boletin Observatorio Tonantzintla Tacubaya, 5, 219
Pismis P., 1977, Rev. Mex. Astron. Astrofis., 2, 59
Reed B. C., 2003, AJ, 125, 2531
Reid M. J., Menten K. M., Brunthaler A., Zheng X. W., Moscadelli L., Xu Y., 2009, ApJ, 693, 397
Robin A. C., Reyle C., Derriere S., Picaud S., 2003, A&A, 409, 523
Robitaille T. P., Whitney B. A., Indebetouw R., Wood K., Denzmore P., 2006, ApJS, 167, 256
Salpeter E. E., 1955, ApJ, 121, 161
Sanchawala K., Chen W. P., Ojha D. et al., 2007, ApJ, 667, 963
Schmidt-Kaler Th. 1982, in Landolt-Bornstein, Vol. 2b, ed. K. Schaifers, H. H. Voigt, & H. Landolt (Berlin: Springer), 19
Sharpless S., 1959, ApJS, 4, 257
Sharma S., Pandey A. K., Ojha D. K., Chen W. P., Ghosh S. K., Bhatt B. C., Maheswar G., Sagar R., 2007, MNRAS, 380, 1141
Siess L., Dufour E., Forestini M., 2000, A&A, 358, 593
Sugitani K. et al., 2002, ApJ, 565, L25
Szymczak M., Hrynek G., Kus A. J., 2000, A&AS, 143, 269
Tej A., ojha D. K., Ghosh S. K., Kulkarni V. K., Verma R. P., Vig S., Prabhu T. P., 2006, A&A, 452, 203
Torres-Dodgen Ana V., Weaver W. B., 1993, PASP, 105, 693
Walborn N.R., 1972, AJ, 77, 312
Walborn N. R., Fitzpatrick E. L., 1990, PASP, 102, 379
Zinnecker H., McCaughrean M. J., Wilking B. A., 1993, in Protostars and Planets III, ed. E. Levy & J. Lunine (Tucson: Univ. Arizona Press), 429