Received 5 June 2015; accepted 25 September 2015; published 28 September 2015

ABSTRACT

In this paper we present astrophysical parameters ofthe very old open star cluster Berkeley 39 using 2MASS JHK data and observation ofBV CCD photometry (V ~ 21.5 mag) in Johnson-Cousin filters. We havederived reddening, distance, metallicity and age etc. using two color diagrams in near-infrared and optical bands.We have found E_B-V = 0.099±0.02, E_J-H = 0.018 ± 0.001 and E_J-Ks = 0.016 ± 0.001 mag, and the distance modulus of (m-M)₀ = 13.24±0.03 mag, which give a distance of 4.2±0.04 Kpc, and the age of the cluster is 5.6 Gyr with metallicity; Z= 0.008 (obtained from fitting of the observed CMDs with Padovaisochrones). The radius of the cluster obtained from the analysis of 2MASS data is equal 9.0±0.13 arc minutes; it reflects that the cluster have large radius with new cluster center. Luminosity and mass functions have been investigated as well as the dynamical relaxation time of the cluster Berkeley 39.

Keywords:

Galaxy, Open Clusters, Associations, Individual, Berkeley 39

1. Introduction

Open star clusters are the excellent tool to trace the properties of star formation, the formation and the evolution of the disc of our galaxies and structural points ofview. Studying the old open star cluster trace back the first moments of the disc lifetime and putting crucial constraints on its formation mechanism and early evolution of the open cluster system, all this allow us to explore the history of star and cluster formation over wide range of the galaxy.

Old open star cluster allow us to investigate the chemical evaluation of the galaxies and stage of formation the galactic disk. Many detailed studies of the old open clusters[1] -[6] . The oldest open star cluster Berkeley 39 is located in the outer part of the disc, about 11 kpc from the galactic center, [5] and reference there in, at (α2000.0 = 07h 47m 18s,δ2000.0 = −04˚44', l =223˚, b =10˚). [1] find that the total absorption in the direction of the cluster Berkeley 39E_(B-V) = 0.15 ± 0.02anddated that the cluster age about 8 Gyr with E_B-V= 0.12 and (m-M)_v = 13.4 mag,based on BV CCD observation. The mean velocity of Berkeley 39 is +55± 10 as determined by [7] . [3] , determined the datasets for Berkeley 39 are: age is 6±2 Gyr and (m-M)_o = 12.9±0.2 mag with adapted the reddening E_B-V =0.11±0.09.

The database of cluster Berkeley 39 published in the WEBDA is Reddening = 0.12 mag., log Age = 9.9 Gyr, with Distance modulus = 13.8 mag, Metallicity=−0.17 and Distance = 4.7 kPc. [8] find the reddening of E_(B-V)= 0.11 ± 0.09 using Hβ line for blue stars in the cluster field. Also [9] determines the chemical composition of Berkeley 39 from high-resolution spectra of giants and finds that the cluster is slightly metal-poor ([Fe/H] =−0.21±0.01). 2MASS JHK bands beside CCD photometry (V ~ 21.5 mag) is taken in BV Johnson-Cousin filters which enable us to properly investigate photometric and structure parameters of theclusters.

Old open star cluster is very important objects, however, Berkeley 39 is important to re-investigate itsposition and age, which are estimated in a precision way from their Color-Magnitude Diagrams using the most recent Near?Infrared photometric data from the 2-Micron All Sky Survey (2MASS). Other motivation of this project is to study the useful parameters to improve our knowledge about the distance scale and the kinematic properties and structures of the Milky Way Galaxy.

The paper is organized as follows: the observations and data reduction are described in Section 2. The analysis of the color-magnitude diagrams and the derived photometric parameters are given in Section 3. The Cluster Membership is presented in Section 4. Radial density distribution is presented in Section 5.The luminosity function and mass function are described in Section 6, and mass segregation and dynamical relaxation time are described in Section 7. The conclusion and summary are presented in last section.

2. Observations and Data Reduction

The cluster was observed with the direct camera of the 188 cm reflector telescope at Newtonian focus (f/4.9) during one night on 2013, March 19, in BV CCD images.The telescope is operating by Kottamia observatory, astronomy department of NRIAG, Egypt.It was equipped with a 2k × 2k EEV 42−40 CCD camera, its pixel size of 13.5 µm. The telescope and this CCD camera give us a square field of view of about 10 × 10 arcmin² on the sky with a scale of 0.30 arc sec per pixel, more information about the capabilities of Kottamia telescope are presented in [10] . Log of observation gives in Table 1, the numbers of observations for each filter, exposer times and seeing are listed. Bias and twilight flat field frames were taken in the same period and same filters of observation. The observations of the standard stars in the selected area SA 107, [11] have been performed to calibrate the target observations also in the same period of observations [12] .

All CCD reductions for the cluster and standard stars were done using IRAF. These processes include overscan, bias subtraction, flat field corrections, and removal of cosmic rays in the standard way. 20 sky flat field were used to flat field the object and standard fields. The instrumental magnitude processes of the cluster were obtained using point spread function using IRAF photometry tasks.The magnitude transformations to the standard system (Londolt 1992) also have been done using IRAF tasks [12] . The transformation equations are in the following form:

whereB, V, and b, v, are the standard and the instrumental magnitudes respectively. X is the airmass, while z_b, z_v, and k_b, k_v, are the photometric zero points and the extinction coefficients in B, V filters respectively. The a_b,a_v, are the color coefficients. The values of extinction coefficients, the color coefficients and the photometric zero points in their corresponding filters are given in Table 2. The near-infrared JHKs filters (2MASS) photometric data were extracted from Vizier web page for the star cluster Berkeley 39.

3. The Color Magnitude Diagrams

The most important parameter in the cluster is the age and distance. Color-Magnitude-Diagram (CMD) has been used to determine the age, reddening and distance.So, the optical (V, B-V), and infrared (J, J-H), (H, J-H), and

(K_s, J-H) CMDs for the total number of the observed stars in the decontaminated cluster region (rB-V = 0.099±0.02, and the color excess values for 2MASS E_J-H= 0.018 ± 0.001 and E_J-Ks = 0.016 ± 0.001 mag, the distance modulus of (m-M)₀ = 13.24±0.03 mag, which gives a distance of 4.2±0.04 Kpc.

Some of these results agree well with which given by [1] , dated that the total absorption in the direction of the cluster Berkeley 39E(B-V) = 0.15 ± 0.02anddated that the cluster age about 8 Gyr with E_B-V= 0.12 and (m-M)v = 13.4 mag.,While [2] , determined the datasets for Berkeley 39 are: age is 6±2 Gyr and (m-M)o = 12.9±0.2 mag, The reddening E_B-V =0.11± 0.09 determined by[2] [7] found the reddening of E(B-V) = 0.11 ± 0.09 using Hβ line for blue stars in the cluster field.

4. Photometric Selection of Candidate Cluster Members

Determination membership of stars cluster from field stars is an essential task to obtain an accurate photometric study for star clusters. To decontaminate the membership stars in a cluster, we taking in account the area of the optical field of view is small comparing to 2MASS area, we used a technique given by [14] as described in the following steps:

a) In the V vs. B-V color magnitude diagramin Figure 1, the stars have been selected as probable members all those stars which according to their σV and σB-V errors, belong to a well-defined strip in the CMD. The envelope of this strip follows the representative main sequence for the cluster.

b) We considered at first the stars chosen above and have 2MASS photometry available along an analogous strip, in the J vs. J-H color―magnitude diagram, and rejected from our sample of initial possible candidate members stars which, have distance from the center less than cluster region (r

c) Finally, we considered as possible cluster candidate members those belonging to the strip defined in the V vs. B-V and J vs.J-H colors magnitude diagram.The star member in the 2MASS data more than the optical one due to the difference in field of view.

5. Cluster Center & Radial Density Distribution

Due to the small area of the optical field of view, we investigate the cluster structure based only on the 2MASS data, starting with determination of the cluster center. The cluster center determination by constructing a rectangular strip around the cluster center (visually determined) in α andδ direction with a side of one degree, we count the stars in each strip to build the frequency distribution in both directions. The histograms were fitted by a Gaussian function, whereas the location of maximum number of stars (peak) indicates the new cluster center,

at (α2000.0 = 07h 46m 42s, δ2000.0 = −04˚39', l =223.29˚, b =10.02˚), asshown in Figure 3.

The internal and external dynamical process is affected of stellar density distribution of a star in the cluster. A good way to investigate the evidence of extra tidal extension in the cluster is the radial surface density profile of the cluster (the spatial distribution of stars around the cluster center) [15] . We examined the radial surface density distribution of stars as a function of distance from the center of the cluster. The number of stars per unit area at different radius from cluster’s center outward is the Radial Density Profile (RDP). The radial density profile is determined through counting the stars in concentric rings around the cluster’s center. We calculated the density of each ring by dividing the number of stars in the ring by its area. The ring radius was chosen to be ranged from 0.25 to 15 arc-minutes. We applied the empirical King model [16] to fit the observed RDP. The equation of King model for the density function is in the following form:

wherer_core, f₀,andf_bgare the core radius, the central surface density and background surface density and their values of r_core = 3 ±0.09 arcmin, f₀=0.32.7± 0.20stars/arcmin² andf_bg = 4.1 ± 0. 21stars/arcmin²respectively.

The radius at which the central surface density drops to its half value is called the core radius. The density profile must represent the well exponential decreasing of the cluster surface density outward from its center, as well as the constancy of the dense background after the cluster limit as it has been shown on Figure 4. We suppose that the background star density lies within the lines represent the standard deviation of the background density level, therefore we define the cluster limits (cluster limited radius; r_lim = 9.2 ± 0.17 arcmin) at the intersection of the King profile with the upper limit of the density standard deviation line, see Figure 4, the background star density at this point (), where is uncertainty of.

6. The Luminosity and Mass Function

The 2MASS uniform sky coverage has been shown to be an ideal tool for studies the cluster area. In this section we analyze the observed stars counts as a function of magnitude to obtain the Luminosity Function (LF) and Mass Function (MF), as well as their spatial dependence. Near-infrared surveys are very useful for the investigation of clusters. It is relatively less affected by high reddening from the Galactic plane. However, due to the small area of the optical field of view, we investigate the Luminosity Function (LF) based on the 2MASS data, it

is given in terms of the absolute magnitude (MJ), obtained after applying the distance modulus derived in Section 3, the frequency distribution of the J absolute magnitude has been obtained, as shown in Figure 5, and the bin used is ∆MJ = 1 mag. LFs inside Berkeley 39 have been built following the structures present in the radial density profile (Figure 4).

The mass function (MF) is obtunding from (LF) using the theoretical evolutionary tracks and their isochrones with different ages selected [13] . Masses of possible cluster members were derived from the polynomial expression developed by [17] with solar metallicity.The absolute magnitude (MJ) and the masses for the adopted isochrones with metallicity Z = 0.008 are used to construct the relation between M_ʘ and absolute magnitude. The relation is a polynomial function, used to determine masses of the cluster from the observed absolute magnitude (MJ), then the number of stars as a function of mass interval was preformed and presented in Figure 6. It is clear that the MFs shown in Figure 6, being flat in the core of the cluster region and very steep in the outer cluster region log (M/M_ʘ) ≥ 1. The mass frequency distribution is linearly fitted and it is noted that the investigated MF slope value of −4.808, the properties and evolution of a star are closely related to its mass.

7.Mass Segregation and Dynamical Relaxation Time

Any stellar system may be characterized by the dynamical time scales (the two body relaxation time) and mass segregation of star cluster, which is an important tool to investigate the star forming process and dynamical evolution of star clusters. The mass segregation resulting from star formation, fast dynamical evolution of star from their parent molecular cloud, i.e. the distribution of masses of the stars along the radius of the cluster showed a clear segregation of mass in the cluster, with increased concentration of massive stars towards the centre. The relaxation time is the most important vector in the dynamics of any stellar system. After a long enough time, the star’s perpendicular speed will (on average) grow to equal its original speed, which is called the relaxation time. We have to estimate the relaxation time to decide whether the mass segregation is primordial or due to dynamical relaxation. To compute the dynamical relaxation time for the clusters, T_E, we have used the relation given by [18] .

Where R_h is the radius containing half the cluster mass, N is the number of cluster members, and is the average star mass in the cluster. The calculated dynamical relaxation time (T_E = 66.22 Myr) with R_h equal to half of the cluster radius in linear units. We have been comparing the values of T_E with the cluster ages, we find Age/T_E = 78.4given the fact that therelaxation times are shorter than the cluster age due to the evaluation of the system and the cluster is more relaxed and stable.

8. Conclusion

In the present work, the observation of open star clusters Berkeley 39 in optical BV and JHKs 2MASS (infrared) data is used to determine the cluster structure and photometric parameters. The optical and infrared dataset reveals that the cluster is very old open star cluster, and rich star cluster and the cluster in stability and relaxed star cluster. The astrophysical parameters have been estimated some of these parameters spatial infrared parameters for the first time for these clusters. The main parameters are summarized and listed in Table 3.

Acknowledgements

This paper is a part of the project No. STDF-1335; funded by Science & Technology Development Fund (STDF) under the Egyptian Ministry for Scientific Research. The project team expresses their deep appreciation to the

administrators of STDF and its organization.This publication also makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation”.

Cite this paper

A.Essam,I. M.Selim, (2015) Optical and Infrared Photometric Study of Old Open Star Cluster Berkeley 39. International Journal of Astronomy and Astrophysics,05,173-181. doi: 10.4236/ijaa.2015.53022

References