THE SKY AS LABORATORY – 2008/2009 1 Photometric analysis of the open cluster NGC 2420 and of the NGC 6229

Boldrin Elena(1), Bombieri Davide(1), Ercole Luca(2)

(1)Liceo “Giuseppe Berto”, sez. Scientifico, Mogliano Veneto (2) Liceo “Giordano Bruno”, sez. Scientifico, Mestre

ABSTRACT under 500 Myr, and reducing to the higher spectral types. We studied the photometry of the globular cluster NGC Obviously, when the age of the cluster grows, all the 6229, and the open cluster NGC 2420, by using g, r features typical of the diagram become evident and the and i bands. By transforming g ,r, i into BVR Johnson- in each region have a certain mass. Cousins system and constructing CMD diagrams, we The point in the diagram where the stars go out the MS determined distance modulus (DM) by ZAMS fitting for is called turn-off (TO), and can be used to calculate the NGC 2420 and Mv(RR Lyrae) vs V(RR Lyrae) in HB age of the cluster, whether globular or open. for NGC 6229. Finally we determined the age of both Since every of the cluster is close to the others, NGC 2420 and NGC 6229 by fitting different sets of relating to the distance of the cluster from Earth, the isochrones, for different Z, onto CMD difference between apparent and absolute magnitudes diagrams. The age of NGC 6229 has also been is called distance modulus, DM, and is constant. estimated by determining difference MV(TO)-MV(HB). Therefore a color-magnitude diagram is also a HR diagram and this allows us to estimate the distance if we have the DM. I. INTRODUCTION A is a group of stars born in the same age from the same molecular cloud with the same chemical II. OBSERVATIONAL DATA composition, kept together by gravitational force. We can distinguish between globular and open clusters. We produced the photometry of the globular cluster The first ones are generally bigger and composed of a NGC 6229 (fig.1) high number of old population II stars, and are distributed in the whole galactic halo. They have a visible spherical shape due to the strong gravity, that keeps a high stellar density in their centre. They do not have any gas or dust. Open clusters have a lower number of stars, so their gravitational bound is weaker. They are generally very bright, because they contain young stars at a high temperature. Analysing the light coming from an open cluster we can estimate its age. An abundance of blue stars means that the cluster is young. The stellar density can be extremely variable, sometimes similar to that of the field stars, making the cluster hard to recognize; anyway it is always lower than the density of the central areas of a globular cluster. Open clusters can be found only in the galactic disk of spiral and in irregular ones. Generally, since all the stars in a cluster have about the same age, chemical composition and distance, every difference in the stellar parameters is caused just by the initial mass of each star, that determines a different Fig. 1 NGC 6229 evolution. When we create a CMD or a H-R diagram of the stars and of the open cluster NGC 2420 (fig. 2) in a cluster, its shape depends on the age of the cluster, showing an almost complete MS for the young ones, THE SKY AS LABORATORY – 2008/2009 2 subtracting off an average sky count, determined in a circle ring around the object, to obtain the effective flux. In this case this technique was not suitable because of the high stellar density, that made inaccurate, or even impossible for the globular cluster, the average sky count determination. Therefore we applied the PSF photometry (Point Spread Function). This is a mathematical model of the stellar profile which describes the average distribution on the CCD surface of the photons coming from a single star. Because of the Earth atmosphere effect (seeing), this distribution has a Gaussian shape, whose typical width (FWHM). To calculate the PSF we chose some stars which showed a regular distribution of light and presented the characteristic trend of the cluster stars. These stars ought to be quite isolated, in order Fig. 2 NGC 2420 not to have the tail end of the profiles overlapping. Besides, they had to be regularly distributed throughout The globular cluster NGC 6229, situated in the the whole image. IRAF determines, by applying the psf at RA = 16h 46m 58s.9 DEC = instruction, an average profile model which, once +47° 31’ 40” (J2000), has a size of about 4,5’ and is of applied to every star of the image, allows to define the class 4, while the open cluster NGC 2420, in the exact position of the star and the quantity of received Gemini constellation at RA = 07h 38m 23.9s DEC= photons, subtracting at the same time the light intensity +21° 34′ 27″ (J2000), has a size of about 10’ and is of of the background sky and dividing the stellar flux by class I 1 r. the exposition time. As a result, we obtained the stellar The photometry has been realized with the data instrumental magnitudes. We applied this procedure for extracted from the SDSS (Sloan Digital Sky Survey) the images in g, r and i band of both clusters. Data Release 6, whose photometric survey has been made with a 120 Mpx CCD (Charge Coupled Device) Since the instrumental magnitude is only an indication of the 2.5-meters telescope in Apache Point (New of the photon flux coming to the stars, we had to bring Mexico, USA). back to a calibrated magnitude, which can be We have done the photometry in the three photometric determined with the following equation: bands g, r and i, and then we have converted them into the more common BVR bands of Johnson-Cousins m=m0+(ms-25)-kx system.

Characteristic wavelengths of the photometric ugriz m0 = constant depending on the considered filter bands are: ms = instrumental magnitude u 3551 Ǻ k = atmospheric extinction coefficient x = airmass g 4686 Ǻ 25 = constant applied by the program to the r 6165 Ǻ instrumental magnitude values, to make them positive. i 7481 Ǻ z 8931 Ǻ The parameters values, acquired by the SDSS catalogue for our fields, are:

NGC 2420 G r i III. DESCRIPTION OF THE WORK m0 24,45 24,07 23,74 k 0,14 0,09 0,03 To perform the photometric analysis of the FITS files we used IRAF software (Image Reduction and x 1,06 1,06 1,06 Analysis Facility), in particular DAOPHOT and IMAGES packages and their daofind, phot, psf, allstar NGC 6229 g r i and tv mark, imexamine tasks. The data are divided into m0 24,58 24,12 23,75 rectangular fields, partially overlapping to allow a k 0,20 0,11 0,06 precise combination, based on the common elements of x 1,18 1,18 1,18 two or more adjacent fields. One of the two clusters we studied was entirely contained in just one field, while We obtained the g, r and i magnitudes of 915 stars in for the other one we matched the images of two the open cluster NGC 2420 and 1970 in the globular adjacent fields within the three studied bands, by using cluster NGC 6229, and then we converted them in B, V the imcombine function. and R by using the transformation equations elaborated The techniques we used to make the photometric by Jordi et al. (2006) [1]. analysis were the aperture photometry and the PSF We used these values to elaborate the color-magnitude fitting. The first one consists of adding up the pixel diagrams (B-V vs V), from which we determined the counts within circles centred on each object and THE SKY AS LABORATORY – 2008/2009 3 distance, the color- (B-V vs Mv) diagrams and the age of the two clusters. To determine the age of NGC 2420 we fitted In case of NGC 2420 we corrected the B-V color isochrones with different metallicity Z obtained from indices on the basis of a mean color excess =0.05, published in scientific literature for this bin/cmd_2.1 getting the best result with ages between cluster [2][3][4] , according to the relation: 1,26 Gyr and 1,59 Gyr and metallicity Z=0,01, as fig. 4 shows.

(B-V) = (B-V)0 + The correction process of the B-V and V values has not and the V magnitudes for the interstellar extinction been done for the globular cluster NGC 6229 since its according to the relation: E(B-V) proves to be equal to 0 or at most 0,01 from literature. V=V0+R We built up the (B-V)0 vs V diagram (fig. 5),

where we put R=3.1. (B-V)0 - V

15,0

(B-V)0 - V0 E(B-V)=0,05 R=3,1 16,0

13,60

17,0 AGB

HB 18,0 15,60 RR Lyrae 19,0 V RGB

Mv+13,3 20,0 17,60 0 V 21,0 TO

22,0 19,60

Mv+13 Mv+13,6 23,0

-0,3 -0,2 -0,1 0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 1,1 1,2 1,3 1,4 1,5 21,60 (B-V)0

Fig.5 (B-V)0 vs V diagram 0,10 0,30 0,50 0,70 0,90 1,10 1,30 1,50

(B-V)0 and we deduced the DM by using the mean HB visual Fig. 3 ZAMS fitting magnitude and in particular its zone where there are the pulsating variable stars RR Lyrae. These stars have a To find the cluster distance we translated the ZAMS, typical B-V color index between 0,2 and 0,4, and their obtained with Johnson’s data [5] until we fitted the absolute magnitude is known, since it depends on correct diagram, determining the distance modulus DM [Fe/H]. It has been estimated in various ways for and obtaining (fig. 3) globular clusters that these stars have Mv=0,60±0,12 at [Fe/H]= –1,5, in particular after the studies carried out DM(NGC 2420)= 13,3± 0,2 with HIPPARCOS satellite [6]. The HB, and therefore the RR Lyrae, have a visual magnitude of +18,1, therefore by comparing the (B-V)0 - Mv E(B-V)=0,05 R=3,1 isochrones Z=0,01 observed visual magnitude value with the absolute one 0,00 we can obtain a DM= V - Mv = 18,10-0,60= 17,50.

1,00 Afterwards we plotted also for this cluster the diagram

experimental data in absolute magnitude (B-V vs Mv) and we fitted 2,00 1,259 Gyr different-metallicity isochrones to establish its age. We 1,413 Gyr 3,00 obtained the best fitting with isochrones having 1,585 Gyr Z=0,001, and age between 10,00 Gyr e 14,13 Gyr. 4,00 V M

5,00 Another method we applied to calculate the age is based on the difference between the stars which are on 6,00 the HB and the ones on the TO according to the

7,00 following equation:

8,00 Log t9=0,37 ΔV -0,03

9,00 0,10 0,30 0,50 0,70 0,90 1,10 1,30 1,50 where ΔV =MV(TO)-MV(HB)=3,2. (B-V)0

Fig. 4 CMD of NGC 2420 and Isochrones THE SKY AS LABORATORY – 2008/2009 4 [2] THE GAP IN THE COLOR-MAGNITUDINE (B-V)0 - MV isochrones Z=0,001 DIAGRAM OF NGC 2420: A TEST OF -3,0 CONVECTIVE OVERSHOOT AND CLUSTER AGE Pierre Damarque, Ata Sarajedini and X.-J. Guo -2,0 experimental data The Astrophysical Journal, 426:165-169, 1994 May 1 10,00 Gyr 12,59 Gyr 14,13 Gyr -1,0 [3] UBVI CCD PHOTOMETRY OF THE OPEN CLUSTER NGC 2420

0,0 Sang Hyun Lee, Yong-Woo Kang, and Hong Bae Ann Korean Astronomical Society 14:61-67, 1999 1,0 V M ΔV=3,8-0,6 [4] GALACTIC GLOBULAR AND OPEN 2,0 CLUSTERS IN THE SLOAN DIGITAL SKY SURVEY. I. CROWDED FIELD PHOTOMETRY 3,0 AND CLUSTER FIDUCIAL SEQUENCES IN ugriz AAVV 4,0 APJS

5,0 [5] published in: “ Basic Astronomical Data “ Chicago U.P. page 216 chapter 11 “Photometric

6,0 System.” -0,3 -0,2 -0,1 0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 1,1 1,2 1,3 1,4 1,5 Johnson - Iriarte (B-V)0 Fig. 6 CMD of NGC 6229, Isochrones and ΔV [6] AN ASTROMETRIC CALIBRATION OF RR LYRAE PERIOD-LUMINOSITY-METALLICITY RELATIONS George Fritz Benedict IV. RESULTS The distance of the two clusters is calculated starting from the DM through the equation:

M–m=5–5Log d which gives, for the open cluster, a distance of 4,6 kpc, a measure which remarkably differs from the officially accepted values. Likewise the age determination to 1,4 Gyr does not agree with the one reported by literature, although the sources report meaningfully different data, and this is due to the peculiar characteristics of the cluster. In this study probably there is a systematic error in the acquisition of the photometrical data, which causes the translation of B-V values of about -0,2, so that ZAMS are fitted to a higher DM, consequently increasing the estimated distance and the isochrones fitting, therefore the valued distance. As regards the globular cluster our results totally agree with literature; actually, with the isochrones-fitting method, the distance proves to be 31,6 kpc and the age 12,6 Gyr. This is substantially confirmed by the age determination of 14,2 Gyr, obtained with the HB-TO distance method.

[1] EMPIRICAL COLOR TRANSFORMATIONS BETWEEN SDSS PHOTOMETRY AND OTHER PHOTOMETRIC SYSTEMS K. Jordi⋆, E.K. Grebel, and K. Ammon Astronomy & Astrophysics manuscript no. 6082 February 5, 2008