About this Journal Submit a Manuscript Table of Contents
ISRN Astronomy and Astrophysics
Volume 2011 (2011), Article ID 127030, 9 pages
http://dx.doi.org/10.5402/2011/127030
Research Article

Luminosity Function of Some Open Clusters

1Departement of Astronomy, National Research Institute of Astronomy and Geophysics (NRIAG), Cairo 11421, Egypt
2Department of Physics, Faculty of Science, Northern Border University, Arar 1321, Saudi Arabia
3Department of Mathematics, Qassim University, Qassim 51452, Saudi Arabia

Received 11 January 2011; Accepted 21 February 2011

Academic Editors: G. Chernov and M. Marconi

Copyright © 2011 W. H. Elsanhoury et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

We investigated the luminosity function (LF) and initial mass function (IMF) of some open clusters having different ages. To calculate the LF, we followed the classical definition by van Rhijn (1936). Statistical investigation of the dispersion around a range of magnitudes concerning what is called Wielen dip revealed that the dip is unreal. To confirm the unreality of the dip, we computed the IMF for these open clusters, the statistical investigation of the IMF confirmed the results obtained using the LF, that is, there is no dip for these open clusters under study.

1. Introduction

Open clusters constitute an important part of a process transforming gas and dust into stars. They are observed as the most prominent parts in the regions of active star formation, or as tracers of the ceased star formation process in the general Galactic field. However, the role they are playing in this process has still not been fully understood. In spite of their prominence, there are indications that classical open clusters contribute only 10% or even less input [13] to the total stellar population of the Galactic disc [4].

The LF is among the most important tools to analyze the composition and evolution of stellar systems and is especially effective in studies of open star clusters, since the absolute magnitudes and ages of their constituent stars can be determined fairy accurately.

Wielen [5] derived the LF (𝜙) for solar neighborhood stars with aid of Gliese's catalogue [6, 7] by counting Gliese's stars in appropriate volumes of space. In this work, the unit of 𝜙 is stars per unit magnitude interval in a complete sphere of radius𝑟=20 pc.

Miller and Scalo [8] presented a comparison between McCuskey [9], Luyten [10], and Wielen [5] LFs along the LF adopted by him. The LF of Wielen and McCuskey shown a small dip around 𝑀𝑉=7. Miller and Scalo [8] attach there is no physical significance of this dip because (a) its size is within the errors of the determination of the LF and (b) this dip is not found by all investigators.

Lee and Chun [11] selected some of open the clusters (22 open clusters) and 3 associations on the following criteria: (i) correct photometric observations, (ii) easy selection of member stars, (iii) wide range of magnitude and mass, and (iv) large number of MS stars, taking into account the brightening effect of bright MS stars. Using these objects, they derived the LF and mass function MF. The results obtained were consistent with solar neighborhood stars by Wielen [5].

Lee and Chun [11] indicated the presence of a dip at 𝑀𝑉7 mag. or log0.2 in open clusters, which were studied in his paper, similar to what was called Wielen dip which appears in solar neighborhood stars in spite of the fact that Wielen, in his study [5] using solar neighborhood stars did not mention the existence of any dip. Moreover, Lee and Chun [11] mentioned that conclusive confirmation of the existence of this dip in the cluster initial mass function IMF must be made by using many LFs of open clusters older than Hyades because the possibility of the turnover of cluster LF at low mass (<1) has been suggested [12, 13].

Lee and Sung [14] studied the LF for the Pleiades cluster, showing what they called Wielen dip at 𝑀𝑉=7 mag. Also, Lee et al. [15] studied the LF for the Hyades and Praesepe clusters. According to this study, a Wielen dip was clearly seen at 𝑀𝑉=9 mag. (more precisely at𝑀𝑉=8.59 mag.). This dip is about 2 mag. fainter than the case of the Pleiades cluster whose Wielen dip position is consistent with that for the solar neighborhood field stars.

Belikov et al. [16] studied the fine structure of Pleiades LF and premain sequence evolution, they found that there are three features (dips) in the observed cluster LF in a magnitude range 𝑀𝑉=5 mag.–12 mag. The first dip together with the adjacent maximum has a pre-MS nature. In their study of the theoretical LFs construction and comparing with the observations, they argued that since the LF depends on both IMF and mass-luminosity relation MLR derivative and since MLR of pre-MS stars differs from that of MS stars, the pre-MS branch (in the cluster where it is present) should affect the LF. Especially, a sharp contrast is expected at (or just below) the turn on the point at the transition between MS and pre-MS branches. According to Piskunov et al. [17], this transition produces a local “bump” in the LF followed by a dip (herein they call this “bump” as “H-peak” and the effect of peak and dip existence as “H-feature”: they chose this abbreviation to stress that the feature is related to the beginning of hydrogen burning as a star approaches the ZAMS) while the two dips (at 𝑀𝑉=7.5 mag. and 𝑀𝑉=9.5 mag.) are assumed to be field LF features: Wielen and Kroupa et al. (WK for short) dips, into which theoretical models fail to reproduce them.

In the study of the low-luminosity stellar MF, Kroupa et al. [18] discussed possible mechanisms causing the WK dips. They proposed that both details of the LF are the result of mass-luminosity relation fine structure generated by the equation of state and the opacity law. They proposed this idea by causing low-mass ZAMS models and showed that the Kroupa et al. dip is caused by the influence of H2 molecules on the equation of state. They attributed the Wielen dip to the increasing importance of H ions as an opacity source. Belikov et al. in [16] using the more sophisticated models of D’Antona and Mazzitelli [19], concluded that these dips did not appear in their results.

In the present paper, we tried to study the LF of the some open clusters of different ages and magnitudes, with concentration on magnitude range 6𝑀𝑉9 (the range of the so called Wielen dip) in order to(i)test the reality of the presence of this dip (i.e., Wielen dip) in open cluster,(ii)to find the correct position of that dip, if present, and its dependence on cluster age.

2. Observational Materials

If the number of stars is increased in the study of luminosity function of open clusters, the statistical weight can be increased, reducing the possible selection effort. This general approach can be achieved only if the each cluster has statistically enough member stars. Accordingly, the criteria of the selection process of the open clusters to be included in this study to meet the goals that must be achieved were as follows: (a) the open cluster must include enough numbers of most probable members that cover the range of mag. (6𝑀𝑉9), and (b) the most probable member stars in the cluster are determined from proper motion and position.

In our study of the local population of Galactic open clusters, we obtained data that now allow a reliable construction of their luminosity and mass functions. Our cluster sample identified in the all-sky compiled catalogue ASCC-2.5 [20]. For each star projected on a cluster area, a membership probability was determined in an iterative that takes spatial, photometric, and proper motion distributions of stars into account within the corresponding area on the sky [21]. At the end of the iterations for each cluster, we obtained new coordinates of the cluster center, the cluster size, the mean proper motion, the distance from the Sun, reddening, and age [22, 23]. These parameters were determined with data on the most probable members, for example, stars having kinematic membership probability higher than 50%. The results are included in the catalogue of open cluster data (COCD) and its extension [22, 23]. This completeness limit corresponds to a range of magnitude6𝑀𝑉9. Table 1 shows the basic data of 9 open clusters selected according to the above criteria.

tab1
Table 1: Basic data.

3. Method of Analysis

The LF in our analysis was determined using the classical method of Van Rhijn [24]. It is defined as the number of stars in certain magnitude interval per unit volume. The volume and magnitude elements were determined by dividing the cluster into shells with certain width depending on the number of the stars in each cluster.

Double frequency tables of the number of stars in each radius interval (shell) corresponding to each magnitude interval were established. In order to determine the LF from the double frequency table, we counted the number of stars which fall in the magnitude interval corresponding to each distance interval, then we divided the counts in each shell by the volume of this shell (i.e. density) and constructed a table representing the mean absolute magnitude and density𝜙(𝑀𝑉). Table 2 represents a sample of double frequency table for Stock 2 as an example, into which the number of members is 85, distributed into equal magnitude intervals (i.e., absolute magnitude) and distance intervals (measuring in parsec unit).

tab2
Table 2: Stock 2 double frequency table.

It was found that 80 members are included in the distance interval (4.370; 546.223 pc), the remaining 5 members were distributed in other intervals as shown in the table, according to the statistical distribution of the data. Using the double frequency table, we calculated the LF as it is represented in Table 3 according to Van Rhijn’s [24] definition of LF.

tab3
Table 3: Stock 2 luminosity function calculations.

4. Results

4.1. The Luminosity Function

Depending on the previously mentioned statistical method, the LF of the studied clusters of our interest were calculated. The LF of the clusters studied are tabulated in Tables 4, 5, 6, 7, 8, 9, 10, 11, and 12.

tab4
Table 4: Stock 2 LF.
tab5
Table 5: NGC 869 (h-Per) LF.
tab6
Table 6: NGC 884 (χ-Per) LF.
tab7
Table 7: NGC 2168 (M35) LF.
tab8
Table 8: NGC 2264 (S Mon) LF.
tab9
Table 9: NGC 2516 LF.
tab10
Table 10: NGC 2632 (Praesepe, M44) LF.
tab11
Table 11: NGC 3532 LF.
tab12
Table 12: Ruprecht 147 LF.

NGC 2264, 2168, and Praesepe had been studied before by Lee et al. in 1984, the data used were selected on basis of photometric observations and proper motion studies. The LF and MF of these clusters are different with cluster age. On the other hand, other clusters (except Ruprecht 147) had been also studied by Lee and Chun in 1988. According to these studies, the combined LFs of open clusters and associations, the initial luminosity function ILF, and initial mass function IMF of these clusters were derived. They show slight differences in the bright and faint parts as compared with the field star ILF and IMF. However, the general shapes of the cluster ILF and IMF agree well with those for field stars in the solar neighborhood. We note that, all of these above studies [11, 25] depend on counting the number of stars within certain magnitude interval, regardless of the classical definition adopted by Van Rhijn [24].

Recently, Kalirai et al. [26] studied luminosity functions for NGC 2168: the LF rises until 𝑀𝑉=5 and then dips down at 𝑀𝑉=7 and rises again. This dip was most likely the result of poor statistics or cluster stars being subtracted off in the blank field and not physical. In that sense, Kalirai et al. [26]: stated that we do note that the luminosity function given in Barrado Y Navascués et al. [27] does not show a dip at this magnitude. This was constructed by using the I band LF in Table 3 of Barrado et al. [27] and correcting to V by using the empirical fiducial in Table 1 of that paper.

Sanner and Geffert [28] studied the fields of nine open clusters (e.g., Stock 2, Pleiades, and Praesepe), based on Tycho-2 catalogue. They determined membership probabilities for the stars in the cluster fields from the stellar proper motions and used the Tycho-2 photometry to compute the initial mass function (IMF) for the clusters from the main sequence turn-off point down to approximately. 1𝑀.

In this work, the computed LFs of the studied open clusters are represented in Figure 1. From the general appearance of the curves in this figure, it is so difficult to conclude whether there is a dip or not.

127030.fig.001
Figure 1: LFs of the present studied open clusters.

In order to have an accurate decision for the presence of the dip, we calculated the dispersion (i.e., standard deviation) of LF of the whole curve and the dispersion around the region of the dip which was suggested by Lee et al., 1997.

The results of these calculations are shown in Table 13. From this table, it is clear that the dispersion around the region where Wielen dip position was suggested is much smaller than the dispersion of the whole cluster, which assures the absence of any dip in the region of 7 mag. or 9 mag., except in Praesepe where the dispersion around the dip was slightly larger than the dispersion of the whole curve. The results obtained in the present work are in agreement with those obtained by Miller and Scalo [8], where they concluded that there is no physical significance of this dip because its size is within the errors of the determination of the LF, and the dip is not found by all investigators.

tab13
Table 13: The dispersions of the IMF of the whole cluster compared with that of the dispersions around the dip region.

The absence of any dip in the LF of the studied clusters in the present work makes us more interested to use Lee’s original data [15] for calculating the general dispersion of the figure and the dispersion around the suggested magnitude intervals of the dip.

We redraw Lee's data of 1997 presented in Table 14 as shown in Figure 2 for the three open clusters, Hyades (closed triangle), Praesepe (opened square), and Pleiades (opened circle) with solar neighborhood (closed circle), keeping Lee’s scale in our drawings. Comparing Figure 2 with Lee’s Figure 3 which is copied from the paper published by Lee et al. (1997). It was found that there is a mistake in his Figure 3 regarding to the Pleiades cluster. Into which the LF of the Pleiades cluster lies above the LF of the solar neighborhood stars and below the LF of the Hyades cluster while in Figure 2 the LF of the Pleiades cluster is above all the LF curves of the other open clusters in most parts, and it lies between Hyades and Praesepe clusters starting from 𝑀𝑉15 mag.

tab14
Table 14: Luminosity functions for three open clusters and solar neighborhood stars due to Lee. et al., 1997.
127030.fig.002
Figure 2: LFs for three open clusters and solar neighborhood stars after redrawing using Lee’s data.
127030.fig.003
Figure 3: LFs for three open clusters and solar neighborhood stars, the original graph, copied from Lee et al. 1997.

We calculate the dispersions of the LF curve of the whole cluster and the region around the suggested dip for these three open clusters given by Lee et al. in 1997. The results are shown in Table 15, from which it is clear that the dispersion around the region of the dip suggested by them is much smaller than the dispersion in the LF curve of the whole cluster, also its amplitude was approximately measured based on shape of the curve around the dip, it is clear that these amplitude values are within the dispersion around the dip region as shown in Table 15. This makes us on the side of the nonreality of the suggested dip. Especially, if we know that Lee et al. have included white dwarfs in their LF calculations, and LF as defined in the literature was not followed exactly by them since the volume was not considered in their calculations of LF.

tab15
Table 15: The dispersions in comparison with absolute value of the dip in each cluster.
4.2. The Present Day Mass Function (PDMF) and IMF

The PDMF 𝜙ms(log) is defined as the number of main sequence stars per unit logarithmic mass interval per square parsec in the solar neighborhood. The PDMF of main sequence field star is related to the LF of field stars 𝜙(𝑀𝑉) [8] by𝜙ms𝑀(log)=𝜙𝑉||||𝑑𝑀𝑉||||𝑀𝑑log2𝐻𝑉𝑓ms𝑀𝑉,(1) where 𝑑𝑀𝑉/𝑑log is the slope of the [(absolute magnitude 𝑀𝑉, mass -) relation] and converts the LF to an MF. The term 2𝐻(𝑀𝑉) is the result of integrating the LF perpendicular to the plane of the Galaxy, assuming an exponential distribution with scale height 𝐻(𝑀𝑉). The factor 𝑓ms(𝑀𝑉) gives the fraction of stars at a given magnitude which is on the MS. The values of the absolute magnitudes, masses, 𝑑𝑀𝑉/𝑑log, scale heights, and the fraction 𝑓ms adopted in the present investigations were given by Miller and Scalo [8].

The definition of the field star IMF in terms of all stars ever formed may seem confusing at first. However, in the case of a time-constant IMF (as assumed here), the IMF at any given time has the same shape as the IMF at any other time and, therefore, the same shape as the IMF of all stars ever formed. Note that the field star IMF as defined here and the IMF of stars in open clusters are similar quantities in that both give the mass distribution of stars at birth.

Stars with MS lifetimes greater than the age of the Galaxy will be found on the MS today regardless of when they were formed. For these stars (e.g., open clusters), the PDMF (𝜙ms(log)) and IMF (𝜉(log)) are identical [8, 15]:𝜙ms(log)=𝜉(log),𝑇ms𝑇0.(2)

Now, we aimed to study the IMF of the above list of clusters (where the data were collected from ASCC-2.5 catalogue) to show if there is a Wielen dip at region log0.25(0.56) for Praesepe cluster and log0.12(0.76) for other open clusters [15]. Now, referring to (1) to calculate IMF for the above list of clusters, we get the data in Tables 16, 17, 18, 19, 20, 21, 22, 23, and 24.

tab16
Table 16: The Stock 2 IMF.
tab17
Table 17: The NGC 869 (h-Per) IMF.
tab18
Table 18: The NGC 884 (χ-Per) IMF.
tab19
Table 19: The NGC 2168 (M35) IMF.
tab20
Table 20: The NGC 2264 (S Mon).
tab21
Table 21: The NGC 2516 IMF.
tab22
Table 22: The NGC 2632 (Praesepe, M44) IMF.
tab23
Table 23: The NGC 3532 IMF.
tab24
Table 24: The Ruprecht 147 IMF.

Again, we computed the dispersion of IMF and also the dispersion around the region of the dip (as a function of mass) which was suggested by Lee et al. [15].

The results of these calculations are shown in Table 25. Here, it is clear that the dispersion around the region of the Wielen dip which was suggested by Lee et al. [15] is much smaller than the dispersion of the whole cluster, which assure the absence of any dip in the region of log0.25(0.56) for Praesepe cluster and log0.12(0.76) for other open clusters. It is worthy to note that the dispersion of IMF of whole Praesepe cluster is much larger than the dispersion around the dip.

tab25
Table 25: The dispersions of the LF of the whole cluster compared with that of the dispersions around the dip region.

5. Conclusion

Finally, we concluded that the examination of IMF revealed the same results obtained with the LF, for example, there is no dip in the distribution of both LFs and IMFs. Also, Wielen [29] stated that most people explain the so-called Wielen dip in the LF of field stars as a consequence of a special property of the relation between stellar mass and stellar luminosity in the region of the dip. If the IMF would be universal and valid for clusters too, then we should see the dip also in the LF of clusters. However, there are strong indications that the IMF in clusters differ from that of the field stars. In this case, the dip may not be strongly visible in the cluster LF.

Acknowledgment

The authors would like to thank Professor Dr. R. Wielen, for his communication and comments on this paper.

References

  1. R. Wielen, “The age distribution and total lifetimes of galactic clusters,” Astron & Astrophys, vol. 13, p. 309, 1971.
  2. R. E. Miller and J. M. Scalo, “On the birthplaces of stars,” Physical Audio Signal Processing, vol. 90, p. 506, 1978.
  3. A. E. Piskunov, N. V. Kharchenko, S. Röser, E. Schilbach, and R. D. Scholz, “Revisiting the population of Galactic open clusters,” Astronomy and Astrophysics, vol. 445, no. 2, pp. 545–565, 2006. View at Publisher · View at Google Scholar
  4. A. E. Piskunov, N. V. Kharchenko, E. Schilbach, S. Röser, R. D. Scholz, and H. Zinnecker, “The initial luminosity and mass functions of the Galactic open clusters,” Astronomy and Astrophysics, vol. 487, no. 2, pp. 557–566, 2008. View at Publisher · View at Google Scholar
  5. R. Wielen, “The kinematics and ages of stars in Gliese's catalogue,” in IAU Highlight of Astronomy, G. Contopoulos and G. Contopoulos, Eds., vol. 3, p. 395, Reidel, Dordrecht, The Netherlands, 1974.
  6. W. Gliese, “Katalog der Sterne näher als 20 Parsek für 1950.0,” Mitt. Astron. Rechen-Institut, Heidelberg Serie A, no. 8, pp. 1–89, 1957.
  7. W. Gliese, “Catalogue of nearby stars. Edition,” Veroffentl. Astron. Rechen-Institut, Heidelberg, no. 22, 1969.
  8. G. E. Miller and J. M. Scalo, “The initial mass function and stellar birthrate in the solar neighborhood,” Astrophysical Journal Supplement Series, vol. 41, pp. 513–547, 1979.
  9. S. W. McCuskey, “The stellar luminosity function,” Vistas in Astronomy, vol. 7, pp. 141–171, 1966.
  10. W. J. Luyten, “A new determination of the luminosity function,” Monthly Notices of the Royal Astronomical Society, vol. 139, pp. 221–224, 1968.
  11. S.-W. Lee and M. Y. Chun, “Initial luminosity function and initial mass function for open clusters and associations,” Journal of Korean Astronomical Society, vol. 21, p. 37, 1988.
  12. S. van den Bergh, “The luminosity function of star formation,” The Astrophysical Journal, vol. 134, p. 553, 1961.
  13. R. B. Larson, “Turbulence and star formation in molecular clouds,” Monthly Notices of the Royal Astronomical Society, vol. 194, p. 809, 1981.
  14. S.-W. Lee and H. Sung, “The luminosity function and initial mass function for the pleiades cluster,” Journal of Korean Astronomical Society, vol. 28, p. 45, 1995.
  15. S.-W. Lee, H. Sung, and D.-H. Cho, “The luminosity function and mass function for the praesepe and hyades clusters,” Journal of Korean Astronomical Society, vol. 30, p. 181, 1997.
  16. A. N. Belikov, S. Hirte, H. Meusinger, A. E. Piskunov, and E. Schilbach, “The fine structure of the Pleiades luminosity function and pre-main sequence evolution,” Astronomy and Astrophysics, vol. 332, no. 2, pp. 575–585, 1998.
  17. A. E. Piskunov and A. N. Belikov, “Fine structure of the luminosity functions of young galactic clusters,” Astronomy Letters, vol. 22, no. 4, p. 522, 1996.
  18. P. Kroupa, C. A. Tout, and G. Gilmore, “The low-luminosity stellar mass function,” Monthly Notices of the Royal Astronomical Society, vol. 244, p. 76, 1990.
  19. F. D'Antona and I. Mazzitelli, “New pre-main-sequence tracks for M ≤ 2.5 M as tests of opacities and convection model,” Astrophysical Journal, Supplement Series, vol. 90, no. 1, pp. 467–500, 1994.
  20. N. V. Kharchenko, “All-sky compiled catalogue of 2.5 million stars,” Kinematics and Physics of Celestial Bodies, vol. 17, p. 409, 2001.
  21. N. V. Kharchenko, A. E. Piskunov, S. Röser, E. Schilbach, and R. D. Scholz, “Astrophysical supplements to the ASCC-2.5. II. Membership probabilities in 520 Galactic open cluster sky areas,” Astronomische Nachrichten, vol. 325, no. 9, pp. 740–748, 2004. View at Publisher · View at Google Scholar
  22. N. V. Kharchenko, A. E. Piskunov, S. Röser, E. Schilbach, and R. D. Scholz, “Astrophysical parameters of Galactic open clusters,” Astronomy and Astrophysics, vol. 438, no. 3, pp. 1163–1173, 2005. View at Publisher · View at Google Scholar
  23. N. V. Kharchenko, A. E. Piskunov, S. Röser, E. Schilbach, and R. D. Scholz, “109 new Galactic open clusters,” Astronomy and Astrophysics, vol. 440, no. 1, pp. 403–408, 2005. View at Publisher · View at Google Scholar
  24. P. J. van Rhijn, “The absorption of light in interstellar galactic space and the galactic density distribution,” Publications of the Kapteyn Astronomical Laboratory Groningen, vol. 47, pp. 1–34, 1936.
  25. S.-W. Lee and H. S. Lee, “Luminosity functions and mass functions of eight open clusters,” Journal of Korean Astronomical Society, vol. 17, p. 51, 1984.
  26. J. S. Kalirai, G. G. Fahlman, H. B. Richer, and P. Ventura, “The CFHT open star cluster survey. IV. Two rich, young open star clusters: NGC 2168 (M35) and NGC 2323 (M50),” Astronomical Journal, vol. 126, no. 3, pp. 1402–1414, 2003. View at Publisher · View at Google Scholar
  27. D. Barrado Y Navascués, J. R. Stauffer, J. Bouvier, and E. L. Martín, “From the top to the bottom of the main sequence: a complete mass function of the young open cluster M35,” Astrophysical Journal, vol. 546, no. 2, pp. 1006–1018, 2001. View at Publisher · View at Google Scholar
  28. J. Sanner and M. Geffert, “The IMF of open star clusters with Tycho-2,” Astronomy and Astrophysics, vol. 370, no. 1, pp. 87–99, 2001.
  29. R. Wielen, “Private communication,” 2009.