UBC Theses and Dissertations

UBC Theses Logo

UBC Theses and Dissertations

CCD photometry of NGC 3201 Brewer, James Philip 1991

Your browser doesn't seem to have a PDF viewer, please download the PDF to view this item.

Item Metadata

Download

Media
831-UBC_1991_A6_7 B74.pdf [ 9.29MB ]
Metadata
JSON: 831-1.0053473.json
JSON-LD: 831-1.0053473-ld.json
RDF/XML (Pretty): 831-1.0053473-rdf.xml
RDF/JSON: 831-1.0053473-rdf.json
Turtle: 831-1.0053473-turtle.txt
N-Triples: 831-1.0053473-rdf-ntriples.txt
Original Record: 831-1.0053473-source.json
Full Text
831-1.0053473-fulltext.txt
Citation
831-1.0053473.ris

Full Text

CCD PHOTOMETRY OF NGC 3201 by JAMES PHILIP BREWER B.Sc, The University of Leeds, 1988 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF T H E REQUIREMENTS FOR THE DEGREE OF MASTER OF SCIENCE in T H E FACULTY OF GRADUATE STUDIES (Department of Geophysics and Astronomy)  We accept this thesis as conforming to the required standard  T H E UNIVERSITY OF BRITISH COLUMBIA September 1991 © James Philip Brewer, 1991  In  presenting this  degree at the  thesis  in  University of  partial  fulfilment  of  of  department  this thesis for or  by  his  or  requirements  British Columbia, I agree that the  freely available for reference and study. I further copying  the  representatives.  an advanced  Library shall make it  agree that permission for extensive  scholarly purposes may be her  for  It  is  granted  by the  understood  that  head of copying  my or  publication of this thesis for financial gain shall not be allowed without my written permission.  Department The University of British Columbia Vancouver, Canada  Date  DE-6 (2/88)  ABSTRACT CCD observations of the halo globular cluster NGC 3201, acquired at the Dupont 2.5m telescope at Las Campanas, are reduced and analyzed. The data are from two separate observing runs. During the first run deep data were secured for a field at seven core radii from the centre, whilst during the second run, shallow data were secured for the central field of the cluster. Colour-magnitude diagrams are presented, and a two-colour diagram is constructed using data from the latter run. From the colour-magnitude diagram and colour-colour diagrams, estimates for reddening, metallicity, distance modulus and age are derived, and consistency with values in the literature is shown. Blue stragglers were observed in the central field, and are shown to be more centrally concentrated than subgiants in an equal magnitude interval. A U V bright star was observed in the central field, and magnitudes and colours are reported for it. The possible nature of this object is discussed. Mass functions were derived using deep V and / band data, by binning both on mass and on luminosity. Discrepancies between the two mass functions are discussed and resolved. The counts in a background field are compared to those expected from the Galactic model of Bahcall and Soneira (1980, 1984), and a discrepancy is found at fainter magnitudes. Mass-luminosity relationships (MLRs) in V and / were used to generate an isochrone which showed good agreement with the data, lending credence to the MLRs. Finally, the slope of the mass function is measured for stars with M < 0.4.M©, and values are obtained for the half-mass relaxation time and 'destruction time' of NGC 3201. These values show consistency with the results of Richer et al. (1991).  ii  CONTENTS ABSTRACT  ii  LIST OF FIGURES  vi  LIST OF TABLES  viii  ACKNOWLEDGMENTS  ix  Chapter 1  1  PRELUDE  1.1  Globular Clusters, an Overview  1  1.2  NGC 3201, Previous Studies  3  1.3  Objectives of Present Study  4  Chapter 2  OBSERVATIONS AND REDUCTIONS  6  2.1  Charge-Coupled Devices  2.2  IRAF  .7  2.3  Data  7  2.4  Preprocessing of Data  8  2.4.1  8  2.5  TI Data  2.4.2 TEK Data  9  Production of Program Frames  9  2.5.1  9  Classification of Observations  2.5.2 TI Data 2.5.3 2.6 Chapter 3  6  9  TEK Data  10  DAOPHOT  10  CALIBRATION  20  3.1  Modus Operandi  20  3.2  Aperture Photometry Transformations for the TI Data  21  3.3  ALLSTAR Photometry Transformations for the TI Background Field . . .  3.4  22  ALLSTAR Photometry Transformations for the TI Program Field . 25 3.4.1  Secondary Standards  25  3.4.2  Calibration of the TI Program Field Frames  25  iii  3.4.3  Calibration of the Binned TI Program Field Frames  26  3.5  Error Propagation in the Transformation Equations  27  3.6  Aperture Photometry Transformations for the TEK Data  28  3.7  ALLSTAR Photometry Transformations for the TEK Data  29  3.7.1  Secondary Standards  29  3.7.2  Calibration of the TEK Program Field Frames  Chapter 4  :  29  COLOUR-MAGNITUDE AND COLOUR-COLOUR DIAGRAMS 43  4.1  Data Cleaning  43  4.2  Production of Diagrams  43  4.2.1  TIData  43  4.2.2  TEK Data ; . .  44  4.3  .Two-Colour Diagrams  44  4.3.1  Reddening  45  4.3.2  Metallicity  46  4.4  Distance Modulus  48  4.5  Isochrone Fitting  48  4.5.1  Synthetic Isochrones  48  4.5.2  Deep CMDs  49  4.5.3  Age Determination from Deep CMDs  49  4.5.4  Shallow CMDs  50  4.5.5  Combined CMDs  50  4.6  Binaries  51  4.7  Blue Stragglers  52  Chapter 5  LEE'S UNUSUAL STAR  71  5.1  Observations  71  5.2  Nature of the Ultraviolet Excess Object  72  iv  Chapter 6  STAR COUNTS IN NGC 3201  76  6.1  Motivation  76  6.2  Production of Frames for ADDSTAR Tests  77  6.3  Synthetic Star Tests  77  6.4  Robustness of the Technique  78  6.5  Comparison with the Bahcall and Soneira Galaxy Model  79  6.6  Mass-Luminosity Relationships  81  6.7  Error Propagation in ADDSTAR Tests  81  6.8  The Mass Function  83  6.9  Implications of the Mass Function  86  Chapter 7  CONCLUSIONS  109  BIBLIOGRAPHY  110  v  LIST O F FIGURES 2.1  Location of the TI and TEK fields in NGC 3201  12  3.1  Aperture photometry transformations for the night of May 22/23 1988 . . . 31  3.2  Aperture photometry transformations for the night of May 23/24 1988 . . . 32  3.3  Aperture photometry transformations for the nights of May 22/23 and 23/24 1988  33  3.4  TI backgroundfieldALLSTAR zero points (short exposure frames)  . . . . 34  3.5  TI backgroundfieldALLSTAR zero points (long exposure frames)  3.6  Comparison of aperture photometry between frames with long and short  35  exposure times  36  3.7  TI program field ALLSTAR zero points (short exposure frames)  37  3.8  TI program field ALLSTAR zero points (long exposure frames)  38  3.9  Colourless TI programfieldALLSTAR zero points (long exposure frames)  39  3.10  Colourless TI programfieldALLSTAR zero points (short exposure frames)  40  3.11  Aperture photometry transformations for the TEK data set  41  3.12  TEK program field ALLSTAR zero points  42  4.1  Shallow CMD, with slight cleaning, showing divisions made to distinguish classes of star  54  4.2  Shallow CMD with cleaning  55  4.3  Deep (B - V) CMD  56  4.4  Deep (B - I) CMD  57  4.5  Deep (V - J) CMD  58  4.6  Two-colour diagram for the TEK program field, with luminosity class V fiducials  4.7  59  Two-colour diagram for the TEK program field, with minimum and maximum abundance  fiducials  60  4.8  Deep (B - V) CMD of NGC 3201 with a subdwarffiducialimposed . . . 61  4.9  Deep (B - V) CMD with VB85 isochrones imposed  62  4.10  Deep (V - I) CMD withVB85 isochrones imposed  63  4.11  Deep (V - I) CMD for a backgroundfieldnear NGC 3201  64  4.12  Deep (B - I) CMD with VB85 isochrones imposed  65  4.13  The (B - V) CMD of ALA89  66  4.14  Shallow (B - V) CMD with VB85 isochrones imposed  67  vi  4.15  Combined (B - V) CMD for TI and TEK data sets  68  4.16  Location of blue stragglers on the TEK program  4.17  Cumulative distributions of blue stragglers and subgiants in the TEK program  field  69  field  70  5.1  Lee's unusual star  73  5.2  Finding chart for the TEK program  5.3  Comparison with the photometry of Lee (1977)  6.1  Error diagram for synthetic stars added to the binned TI program field  field  74 75  frames  88  6.2  Varying crowding conditions on the I programfield:  89  6.3  Uncorrected luminosity functions for various crowding conditions  90  6.4  Recovery efficiency for various crowding conditions  91  6.5  Corrected luminosity functions for various crowding conditions  92  6.6  Comparison of frames with 5% and 100% add levels  93  6.7  Stars found in two colours, with bins and incompleteness shown  94  6.8  V mass-luminosity relationship  95  6.9  I mass-luminosity relationship  96  6.10  The effect of cluster age on the V MLR  97  6.11  Mass function for the V TI program  field  98  6.12  Mass function for the I TI program  field  99  6.13  Mass functions for V and I TI program  6.14  An isochrone derived from the V and I MLRs  6.15  Uncorrected mass function for the V TL program  6.16  Uncorrected mass function for the I TI program  6.17  Corrected mass function for the V TI program  field  104  6.18  Corrected mass function for the I TI program  field  105  6.19  A 'magnitude-magnitude diagram' for the TI program  6.20  A 'magnitude-magnitude' diagram for the TI background  6.21  Plots of half-mass relaxation time and 'time until destruction' against mass function slope below OAM0 .  vii  fields  100 101  field field  field field  102 103  106 107 108  LIST O F T A B L E S 1.1  Some parameters of NGC 3201  5  2.1  Characteristics of the TEK and TI CCDs  8  2.2  Characteristics of the TI program field frames  10  2.3  Characteristics of the TEK programfieldframes  10  2.4  Observing log for the TI data set  13  2.5  Observing log for the TEK data set  17  4.1  Corrections of 6 to 6(0.6)  46  4.2  Metallicity values for NGC 3201 from previous studies  47  6.1  Characteristics of frames used for ADDSTAR tests  77  6.2  Comparison of counts with the Bahcall and Soneira model  80  6.3  Properties of clusters with derived mass functions extending well below O.4M  87  0  viii  ACKNOWLEDGMENTS Thanks to life, for being such a wonderful and truly joyous experience.  "We are all in the gutter, but some of us are looking at the stars." Oscar Wilde (1854-1900)  ix  Chapter  1  P R E L U D E  1.1 Globular Clusters, an Overview Globular clusters are systems of stars which are gravitationally bound and long-lived, as opposed to Galactic clusters which tend to be loosely bound and short-lived. Most globular clusters appear spherical, their axial ratios lying between 0.9:1 and 1:1, though a few clusters have axial ratios of up to 0.6:1 in their outer regions. N-body modelling of globular cluster dynamics is not realizable with current computing technology, but their kinematics can be described by King Models (King 1966) which are spatially limited, steady state solutions of the Fokker-Planck equation. The Fokker-Planck equation is an approximation to the collisionless Boltzmann-equation. Globular clusters are Population II objects which, in the Milky Way, are found to have a bimodal metallicity distribution. There are both halo and disk globular clusters, the halo clusters having a low systemic rotation and a high random velocity dispersion, and the disk clusters having a larger rotation and a smaller random velocity dispersion (Zinn 1985). The Milky Way has associated with it a system of 150 to 200 globular clusters, each of which contains approximately 10 stars. Globular clusters, per se, make a negligible contribution to 5  the mass of the Galaxy. Over half of the globular clusters are at a distance of less than 10 kpc from the Galactic centre, the radius of the disk being 25 - 30 kpc (Mihalas and Binney 1981). The distribution extends far out into the halo, for example the globular cluster NGC 2419 is at a Galactocentric distance of 99.4 kpc (Webbink 1985). The distribution above the Galactic plane does not differ significantly from that in the two orthogonal directions; Shapley recognized this, and was able to determine the dimensions of the Milky Way, concluding that the sun is not near the centre of the Galaxy. This is contrary to what had been previously thought by Kapetyn, who had failed to allow for the absorption of starlight by interstellar material (Mihalas and Binney 1981). Synthetic isochrone fitting to globular cluster colour-magnitude diagrams (CMDs) shows them to be approximately 15 gigayears (Gyr) old, with a definite age spread. These ages, on the scale defined by VandenBerg and Bell (1985, henceforth VB85), are hard to reconcile with the age of the universe derived using Hubble's constant, which is believed to lie in the range HQ=50 - 100 k m s ^ p c , implying an age of between 7 and 13 Gyr, assuming a £=0, A=0 —  - 1  l  2  Chapter 1. PRELUDE  Friedrnann-Lemaitre model (Borner 1988). The values may be reconciled if it is found that the absolute ages of the isochrones, which are very uncertain, are somewhat overestimated. A globular cluster will show a large range in its surface brightness. Two parameters are used to characterize the surface brightness distribution, they are •  r =Core radius; the radius at which the surface brightness is half of its central value. c  r =Tidal radius; the radius at which the surface brightness reaches zero, often obtained by t  extrapolation of the surface brightness profile. Of the Galactic globular clusters, 5/6 of them have (0.3 < r < 10) pc and (10 < rt/r < 100). c  c  Inspection of Table 1.1 shows NGC 3201 to be typical in its surface brightness parameters. Most globular clusters have a core where the surface brightness changes slowly, i.e. they have an isothermal velocity dispersion, while for a few, the surface brightness keeps increasing down to the smallest radius which can be resolved in the cluster. These clusters are said to possess 'cusps' and are believed to be in the post core-collapse stage of evolution (Spitzer 1987). The distribution of the integrated luminosity of globular clusters has a mean value of 4.8 with a dispersion of 0.5 in log ^j,"*'"). The central velocity dispersion, which can be obtained using line widths of integrated spectra or measurement of individual giants, combined with the virial theorem, allows for the total mass of a cluster to be estimated (this is model dependant). Using this method, it is found that the mass to light ratio, in solar units, is in the range 1 to 4, though recent studies by Richer and Fahlman (1991) and Richer et al. (1991), where data were fitted to King models, indicate values of around 6. Knowing this, and applying it to fainter clusters from which integrated spectra cannot be obtained, it would appear that the masses of globular clusters lie in the range 10 to 10 4  6  Mq.  The integrated colour of globular clusters, (B — V ) , is found to lie in the range 0  0 ™4 < (B - V ) < 0 ™8, with the distribution peaking at (B - V ) = 0 ™57. This shows that 0  0  the light from globular clusters is dominated by fight from stars cooler than the sun, i.e. of a later spectral type than G5. A n integrated colour is not a good parameter for classification; the colour of a star depends on two factors, viz.: temperature and metallicity. Clusters can be classified into concentration classes, the classes being designated I to XII. In this system, devised by Shapley and Sawyer, class I clusters show the highest central concentration of stars whilst class XII clusters the lowest. The divisions were made such that the clusters divide almost equally between the twelve classes.  3  Chapter 1. PRELUDE  1.2 N G C 3201, Previous Studies  NGC 3201 is in Vela, a constellation in the southern hemisphere near Crux. The constellation of Vela lies partly in the plane of the Milky Way. NGC 3201 is a halo globular cluster, with fairly typical properties, as shown in Table 1.1. One striking thing about NGC 3201 is its high heliocentric radial velocity, which is +493 k m s  (Kinman 1959). This is the largest radial velocity, in the absolute sense, of  - 1  any cluster in the Milky Way. If the effect of the sun's motion is removed from the radial velocity of NGC 3201, the radial velocity is approximately 240 k m s . With this velocity, the - 1  cluster probably cannot escape from the Milky Way, but will be in a retrograde orbit (White and Kraft 1972). Menzies (1967) was the first person to produce a colour-magnitude diagram (CMD) on the BV system for NGC 3201. He measured 25 stars photoelectrically and 300 stars photographically, producing a CMD with a limiting magnitude of V=15 ™6. The results of this work are published in Alcaino's (1974) Atlas of Globular Clusters. From his work, Menzies determined a distance of 6.9 kpc (Alcaino and Liller 1981) and a reddening of E(B - F ) = 0 ™ 0 7 (Alcaino et al. 1989) for NGC 3201. Alcaino (1976) carried out an investigation of NGC 3201, using B and V band photometry, measuring 18 stars photoelectrically and 313 stars photographically. The limiting magnitude of his survey was V=16 ™30. Alcaino derived (m - M) =14 V  ™10 and E(B - V)=0 ™14, though  he advises the use of caution with the reddening value. Alcaino noted that the horizontal branch was richly populated with blue, red and RR-Lyrae stars. Lee (1977) performed UBV photometry on NGC 3201, measuring 78 stars photoelectrically and 649 photographically. In this study, U, B, and V magnitudes were obtained for all stars having V < 15 ™8, within 5.8 arcmin of the cluster centre. Lee noted that the CMD of NGC 3201 resembles that of M 3 , the main difference between the clusters being that NGC 3201 has a concentration class of X , while M 3 has a concentration class of V I (Lang 1980). Using a two-colour diagram Lee deduced E(B — V ) = 0 ™21 + 0 ™02, and the true distance modulus to be (m — M ) = 1 3 ™57, assuming that My(horizontal branch)=0 ™6. Using metal abundance 0  indicators, Lee deduced the metallicity of NGC 3201 as being —1.2 + 0.1 dex relative to the Hyades. Lee (1977) suggests that he may have blue stragglers in his photometry, though he stresses that his photometry is not reliable at the magnitude of the candidate blue stragglers.  4  Chapter 1. PRELUDE  Alcaino and Liller (1981) used a Pickering-Racine wedge, a device which produces two images with a constant magnitude difference between them for each star, to perform photographic photometry on 1452 stars in NGC 3201. The magnitude limit of their survey was V=19 ™6, and was the first study which reached the main sequence turn off.  Alcaino and  Liller determined the turn off to lie at a magnitude of V=18 ™0 ± 0 ™18, and a colour of (B - V)=0™61 ± 0™03. They deduced a distance modulus of (m - M ) = 1 4 ™ 0 7 ± 0™15 F  and a reddening of E(B — F ) = 0 ™ 2 2 + 0™03. Alcaino and Liller also estimated the age of the cluster as lying between 7.9 and 12.3 Gyr by several means, including using the formula of Iben and Rood (1970), fitting Demarque and McClure's (1977) isochrones and Carney's (1980) relationship, based on Demarque and McClure's (1977) isochrones. Penny (1984), working with Griffiths, was the first to produce CMDs for NGC 3201 by the use of a CCD (see §2.1) detector. Using B, V and / filters, they produced CMDs which extend below the main sequence turn off, and are well defined down to V= 21 ™0. Metallicity and reddening were not determined in this study, but VandenBerg (1983) isochrones were fitted to the data, and agreement was found between the morphologies in the colour-magnitude plane. Alcaino, Liller and Alvarado 1989 (henceforth ALA89) used 67 CCD frames obtained at the 1.54 metre Danish telescope at L a Silla, European Southern Observatory, to perform BVRI photometry on NGC 3201.  They found the main sequence turn off at a magnitude  of V=18 ™20 ± 0 VM, and a colour of (B - V)=0 ™57 ± 0 ™02. They report the magnitude difference between the turnoff and the horizontal branch as being 3 ™44 + 0™12. The CMDs were fitted to VB85 isochrones, a best fit being found for Y=0.2, [Fe/H]=-1.27, a=1.65, (m - M ) ^ = 1 4 ™ 0 0 , and an age of 16 ± 2 Gyrs. 1.3 Objectives of Present Study The main objective of the present study is to determine the mass function of NGC 3201, which has not been published. The data are also used to investigate the reddening, metallicity, distance modulus and age of the cluster, as well as looking for a binary sequence and the concentration of the blue stragglers.  Chapter!.  TABLE  5  PRELUDE  1.1: Some parameters  of NGC 3201  Value 10 15 34 -46° 09' 42" Declination, 61950 277 °229 Galactic Longitude, / 8°641 Galactic Latitude, b 0™21 Foreground Reddening, E(B — V) 5.0 Kpc Heliocentric Distance, Rs -7.47 Integrated Absolute Magnitude, My 14™75 Magnitude of Horizontal Branch, V H B Colour of Subgiant Branch at Horizontal OTO Branch Level, (B - V) X Concentration Class 65:60:51 Ratio of Stars in Horizontal Branch,  Quantity Right Ascension, 01950  h  m  s  Reference Webbink (1985) Webbink (1985) Webbink (1985) Webbink (1985) Lee (1977) Webbink (1985) Webbink (1985) Lee (1977) Lee (1977)  Q  N B H : N A : N R H , for r < V  5'.8  Galactocentric Distance, R G C Heliocentric Radial Velocity, v  r  Metallicity, [m/H] Core Radius, r Limiting Radius, r Core Concentration Parameter c  t  Lang (1980) Lee (1977)  9.5 Kpc 493 km s — 1.6 dex 1.58 pc 52.4 pc 1.55  -1  Webbink (1985) Kinman (1959) Webbink Webbink Webbink Webbink  (1985) (1985) (1985) (1985)  Chapter OBSERVATIONS  A N D  2 REDUCTIONS  2.1 Charge-Coupled Devices  The mid-1970's were the dawn of a new age for observational astronomy with the invention of the charged-coupled device (CCD). A CCD is near to being the perfect detector; it has a high quantum efficiency (QE), is linear in response to near its saturation limit, and has a large dynamic range. CCDs were originally designed for television applications, having an imaging and read out area, this type of CCD is known as a 'frame transfer device'. The Texas Instruments 800 x 800 CCD, made for the Hubble Space Telescope wide-field/planetary camera, and the Tektronix CCD, both used for the observations in this thesis, are for non-television applications and have only an imaging area. The principle of operation of a CCD is quite simple: a photon incident upon the silicon of a CCD will cause the production of an electron/hole pair. Electrodes are used to create potential wells in the silicon in which the photoelectrons are trapped. By use of driving pulses the charge can be transferred along the CCD. An analogue to digital converter is used to read out the accumulated charge, and the actual values obtained are in terms of analogue to digital units (ADUs). The read out process has read out noise (RON) associated with it. This noise is reduced by cooling, but at the expense of sensitivity loss at the red end of the spectrum.  A bias frame is a frame containing only the pedestal which was imposed on the CCD before an exposure. This helps to alleviate the low level charge transfer inefficiency. The CCDs used for the observations in this thesis were 'injected' with charge, to create a pedestal, before the actua exposure. MacKay (1986) states that "The only uniform CCD is a dead CCD." The response of any contemporary CCD is not uniform; individual pixels vary in their QEs, and there is often a large scale structure over the chip. To overcome the non-linearities, 'flatfields',images of an area with uniform illumination, are taken using the appropriatefilterand normalized. These are then divided into an image to normalize the response of the pixels. This is not an ideal solution, non uniformities are colour dependent and so it is impossible to correct thefieldin a perfect manner as the objects in thefieldwill have varying colours. Either domeflats,using artificial illumination on the inside of the dome, skyflats,which are images of the dusk, dawn or dark sky are used. MacKay (1986) says offlatfieldingthat it " .. . is at best a compromise and in no way a substitute for intrinsic device uniformity."  6  Chapter 2. OBSERVATIONS AND REDUCTIONS  7  2.2 IRAF  IRAF (Image Reduction and Analysis Facility) is a general purpose software system, containing a suite of programs for astronomical applications, which include routines capable of handling CCD data frames. The IRAF project was initiated in the autumn of 1981 at the Kitt Peak National Observatory, the project later being taken over and further developed by the National Optical Astronomy Observatories. The version of IRAF that was mainly used was 2.9.1, though earlier versions were utilized.  2.3 Data  In this thesis two data sets are analyzed. The deeper of the two data sets was secured using the 2.5 metre Dupont telescope at Las Campanas Observatory in Chile, on the nights of May 22/23 and May 23/24 1988. The observers were G.G. Fahlman and LB. Thompson. This data set will be referred to as the TI data, and observation numbers will have the prefix nee, an acronym for nineteen eighty eight, for the rest of this thesis. The detector used for the observations was a Texas Instruments CCD, which had the characteristics shown in Table 2.1. The telescope configuration gave a plate scale of 0.162 arcsec per pixel, which gives a total field area of approximately 2' 9" x 2' 9". The clusterfieldis located at «i988=10 16 26. and <$i988=—46° 22' 15", and lies at a distance of approximately 7 core radii from the cluster centre (Figure 2.1). The observing log for the TI data is given in Table 2.4. Observations wer secured usingfilterswhich emulate the Johnson B and Vfilters,and the Kron-Cousins / filter. h  m  The second data set analyzed, also secured at the 2.5 metre Dupont telescope at Las Campana Observatory, was obtained on the nights of May 22/23, May 27/28 and May 29/30 1990. The observers were again G.G. Fahlman and I.B. Thompson. This data set will be referred to as t TEK data, and the observations will have the prefix nn, standing for nineteen ninety, for the rest of this thesis. A C100 Tektronix CCD was used as a detector, which had the characteristics shown in Table 2.1. The telescope configuration gave a plate scale of 0.29 arcsec per pixel, giving a totalfieldarea of approximately 2' 29" x 2' 29". The TEK clusterfieldis located at the centr of NGC 3201 (Figure 2.1), and has coordinates a o=10 15 34 and <$io=-46 09' 42". The observing log for the TEK data is given in Table 2.5. Observations were secured using filters that emulate the Johnson U, B, and V filters. h  195  m  s  0  9 5  Chapter 2. OBSERVATIONS AND REDUCTIONS  TABLE 2.1: Characteristics  8  of the TEK and TI CCDs  Property Manufacturer Pixel Size Dimensions Read Out Noise 2  1  Gain  1  Limit of Linearity  1  TI Texas Instruments 15 microns 800 x 800 pixels  TEK Tektronix 27 microns 512 x 512 pixels  12 e pixel  -1  8 e pixel  2.1 e ADU 25000 ADU  -1  2 e ADU 32000 ADU  -1  -1  Notes to Table 2.1: 1) G.G. Fahlman, Private Communication 2) CTIO Observer's manual  Throughout this thesis I adopt the convention that U, B and V represent the Johnson system while / represents the Kron-Cousins system. By 'standard system' I am referring to either the Johnson UBV and/or the Kron-Cousins / system. 2.4 Preprocessing of Data 2.4.1 TI Data  The CCD frames had been previously debiased, but not flatfielded.Three dome flat-fie were produced by combining several individualflats,see observations neeYll —» nee\9\ and nee216 —> nee29?>. Theflatswere combined, using a medianfilterto aid the removal of ray strikes, and normalized to a value of unity. The data frames were divided through by the appropriateflatfield,and trimmed to a size of 798 x 794 pixels to remove overscan regions. The B and V frames showed localized regions where theflatfieldinghad been poor; the mean sk level being 918 in thefieldand 754 in the worst region of the B frame. This may be attribut to problems with the CCD temperature at the start andfinishof the run on the night of May 22/23; the CCD may not have fully cooled, thus causing the effect. (The quantum efficiency a CCD is very sensitive to temperature, and a non uniform temperature distribution could have caused this effect.) It is also possible that the colour-temperature of the illumination, as discuss in §2.1, may have had an effect, the spectral sensitivity changing rapidly in the poor regions.  Fringing is an effect which manifests itself as a ripple pattern on a frame due to interferen in the silicon layer of the chip. No fringing was apparent on any of the frames.  Chapter 2. OBSERVATIONS AND REDUCTIONS  9  2.4.2 T E K Data  The CCD frames had not been debiased and no bias frames had been taken. The frames were debiased by evaluating the mean ADU value in the overscan region of the frame and subtracting this valuefromevery pixel in the frame. It is not expected that the bias will show appreciabl structure over the frame, so this method of debiasing should be valid. A U domeflatwas constructed, as above, using seven frames. Dome flats for B and V had already been constructed at the telescope by using the median offiveframes. The data frames were divided through by the appropriateflatfieldand trimmed to a size of 503 x 511 pixels to remove overscan regions. No fringing was apparent on any of the frames. 2.5 Production of Program Frames 2.5.1 Classification of Observations  The observations consist of essentially three different types offields.What I shall refer to as a programfieldis afieldin the cluster, whereas a backgroundfieldis afieldlocated near the cluster, but devoid of cluster stars. The third type offieldI refer to is a standardfield,which contains stars of known magnitudes which are used in the primary calibration. 2.5.2 TI Data  Where possible, frames were averaged together using the AVSIGCLIP routine in IRAF. This will average frames together, rejecting pixels which deviate from the mean by more than a specified number of standard deviations. Obviously, when only two frames were averaged, it was not possible to apply the AVSIGCLIP routine. Combining frames changes the read out noise and gain of the newframe.The DAOPHOT user's manual (Stetson 1987a) gives the following formulae for when frames are averaged together (2.1)  Gains and read out noise were calculated for each of the new program frames, by the use of equation 2.1. The IRAF task IMEXAMJNE was used to calculate the FWHM of the stellar images on the frames, and the standard deviation of a sample of pixels constituting the sky. These results, for the TI programfieldframes, are summarized in Table 2.2.  Chapter 2. OBSERVATIONS AND REDUCTIONS  TABLE  2.2: Characteristics  No. of Constituent Frames 2 3 5  B V I  10  of the TI program  field  Exp. Time per RON Frame (e pixel ) (seconds) 8.49 900 600 6.93 5.37 300 -1  frames  Gain (e A D U  Sky a - 1  )  4.2 6.3 10.5  (ADU)  26 19 16  FWHM (arcsecs) 1.2 1.3 0.9  2.5.3 T E K Data  All frames were first debiased andflatfielded.Frames were combined together to increase the signal to noise ratio, as for the TI data, where possible. Characteristics of the three TEK programfieldframes built are shown in Table 2.3. TABLE  2.3: Characteristics  of the TEK program  No. of Constituent Frames  Exp. Time per Frame  1 7 9  300 15 15  u B V  RON (c pixel" )  field  frames  Sky  Gain (cADU  - 1  )  rr  (ADU)  FWHM (arcsecs)  (seconds)  8.0 3.02 2.67  2.0 14 18  9 5.3 9.2  1.4 1.2 1.2  2.6 D A O P H O T  DAOPHOT (Stetson 1987b) is a suite of routines written specifically for the reduction of CCD frames containing crowded images of starfields.Each star will produce a point-spread function (PSF) on the frame, the PSF of a star being much larger than its theoretical diffraction pattern (King 1971). The PSF is determined by the phenomena of atmospheric refraction, diffraction and scattering, as well as instrumental diffraction and scattering, it may also depend on telescope guiding errors. Due to the many variable factors, the PSF from one frame to the next will be different, making invalid the use of a universal PSF within a data set. A PSF has to be built for each frame on which profile-fitting photometry is required. DAOPHOT models the PSF by firstfittingthe central regions of a stellar profile with an analytic bivariate Gaussian function in 1  1  The point-spread function is the two-dimensional brightness distribution produced on the focal plane of the detector by an  unresolved source.  11  Chapter 2. OBSERVATIONS AND REDUCTIONS  which the Gaussian acts as a first order approximation to the profile, and then a look up table of residuals from the Gaussian is calculated. Other programs along similar lines (e.g. Penny and Dickens 1986, Peterson et al. 1978, and Buonanno et al. 1983) model a PSF by using a flexible analytical function. There are advantages and disadvantages to both methods, and I refer the interested reader to Stetson (1987b) for a fuller discussion of these matters. A typical reduction of aframeusing DAOPHOT involvesfirst'building' a PSF, using the least crowded stars on the frame as prototypes. The stars on the frame are then 'found' and split into groups, this being necessary because the reduction time for a group of N stars is 0(N ), where as if the group is split up into v subgroups, the reduction time is 0(N /^ ) (Stetson 1987b). DAOPHOT will simultaneouslyfitthe stars within a group, using the PSF as thefittingprofile. The stars thusfittedare subtracted from the frame, leaving stars which were 'missed' on thefirstfind due, probably, to crowding. Thefindingroutine is again employed tofindthe missed stars, and they are grouped andfitted.The twofilescontaining the stars are concatenated and a final simultaneousfitis performed. This general procedure is what I will refer to as a standard reduction for the rest of this thesis. After a standard reduction all of the recovered stars have profile-fitting magnitudes, which I shall refer to as ALLSTAR magnitudes for the rest of this thesis. The name ALLSTAR comes from a routine in DAOPHOT which groups andfitsprofiles to stars. The ALLSTAR magnitudes need to be calibrated to a standard system for much of the analysis. 3  3  2  Chapter 2. OBSERVATIONS AND  12  REDUCTIONS  FIGURE 2.1. Location of the TI and TEK fields in NGC 3201 The location and relative sizes of the TI and TEK fields are shown. Concentric annuli indicate core radii. The core radius is 1.58 pc which, at the distance of NGC 3201, is equivalent to an angular width of 64 arcsec. The limiting radius, r , of NGC 3201 is 52.4 pc, which is equivalent to 33 core radii. t  13  Chapter 2. OBSERVATIONS AND REDUCTIONS  TABLE 2.4: Observing log for the TI data set  Object  M  '  Filter  Exposure (seconds)  Universal Time  Niuht of Ma\ 22 2 * B dome flats  nee 111 —• nee 191  B  90  < 23:17  Program field  nee195  I  20  23:17  Program field  nee196  I  600  23:20  Program field  nee197  I  600  23:27  Program field  nee198  I  600  23:34  Program field  nee199  I  600  23:45  Program field  nee200  I  600  23:52  Program field  nee201  I  600  23:58  A & L Stnds.  nee202  I  20  00:05  A & L Stnds.  nee203  I  20  00:09  A & L Stnds.  nee204  B  20  00:30  A & L Smds.  nee205  B  60  00:32  Program field  nee206  B  60  00:36  Program field  nee207  B  900  00:43  Program field  nee208  B  900  01:00  Program field  nee209  V  60  01:29  Program field  nee210  V  600  01:31  Program field  nee211  V  600  01:43  Program field  nee212  V  600  01:54  PG1323 - A & B  nee213  V  5  02:16  PG1323 - C  nee214  V  20  02:26  PG1323  nee215  V  20  02:33  PG1323  nee216  B  20  02:37  PG1323  nee217  I  20  02:40  PG1323 - C  nee218  I  20  02:44  PG1323 - C  nee219  B  20  02:48  PG1323 - A & B  nee220  B  20  02:56  PG1323 - A & B  nee221  I  20  02:59  SA110 AM/AN/AO  nee250  B  20  08:19  SA110 AM/AN/AO (Focus)  nee251  B  20  08:23  SA110 AM/AN/AO  nee252  V  20  08:27  14  Chapter 2. OBSERVATIONS AND REDUCTIONS  TABLE 2.4: (Continued) Object  Observing log for the TI data set Observation Number  Filler  Exposure (seconds)  Universal Time  SA 110 AM/AN/AO  nee253  B  20  08:31  SA110 AM/AN/AO  nee25A  I  20  08:35  SA110 A L  nee255  I  20  08:38  SA110 A L  nee256  I  10  08:43  SA110 A L  nee251  B  20  08:45  SA110 A L  nee258  V  20  08:48  SA110 CD & C C  nee259  V  15  09:21  SA110 CD & CC  nee26Q  I  10  09:25  SA110 CD & CC  nee26l  B  20  09:30  SA110 CD & C C  nee262  V  20  09:33  SA110 CD & CC  nee263  /  10  09:37  SA110 CE & CF  nee264  /  10  09:43  SA110 CE & CF  nee265  /  5  09:45  SA110 CE & CF  nee266  B  20  09:47  SA110 CE & CF  neeUbl  V  20  09:50  SA113 492/493/495  nee268  V  10  09:54  SA113 492/493/495  nee269  V  10  09:59  SA113 492/493/495  nee21Q  I  5  10:01  SA113 492/493/495  nee21\  B  10  10:04  V dome flats  nee216 -> nee284  V  30  >10:04  / dome flats  nee2S5 ->• nee293  I  10  >10:04  Nifihi r,f May 23 24 PG0918 - D  nee294  I  4  22:51  PG0918 - B  nee295  I  5  22:53  PG0918 - C  nee296  I  5  22:55  PG0918 - A  nee291  I  5  22:56  PG0918  nee29S  I  5  22:58  PG0918  nee299  B  20  23:01  PG0918 - D  nee300  B  20  23:04  PG0918 - B  nee30l  B  20  23:06  PG0918 - C  nee302  B  20  23:08  PG0918 - A  nee303  B  40  23:10  PG0918 - A  nee304  V  30  23:13  PG0918  nee305  V  20  23:16  Chapter 2. OBSERVATIONS AND REDUCTIONS  TABLE 2.4: (Continued) Object  15  Observing log for the TI data set Observation Number  Filler  Exposure (seconds)  Universal Time  PG0918 - D  nee306  V  20  23:18  PG0918 - B  nee307  V  30  23:20  PG0918 - C  nee308  V  30  23:22  Program field  nee309  I  30  23:50  Program field  nee3l0  I  300  23:52  Program field  nee3\\  I  300  00:08  Program field  nee3\2  I  300  00:15  Program field  nee3\3  I  300  00:21  Program field  «ce314  I  300  00:28  Program field  «ee315  I  300  00:34  Program field  nee316  V  300  00:42  Program field  nee3\l  B  300  00:50  Background field  nee31S  B  120  01:07  Background field  nee3\9  B  900  01:11  Background field  nee320  V  120  01:29  Background field  nee32\  V  600  01:36  Background field  nee322  I  60  01:56  Background field  nee323  I  450  02:00  Background field  nee32A  I  450  02:08  Background field  nee325  I  450  02:17  Background field  nee326  I  450  02:26  PG1323 - A & B  nee336  B  20  04:00  PG1323  nee331  B  20  04:04  PG1323  nee338  V  15  04:07  PG1323 - A & B  nee339  V  15  04:09  PG1323 - A & B  nee3A0  I  5  04:12  PG1323  nee34\  I  5  04:14  SA110 AM/AN/AO  nee366  I  5  08:31  SA110 A L  nee367  I  5  08:34  SA110 A L  nee368  V  20  08:37  SA11G AM/AN/AO  nee369  V  20  08:39  SA110 AM/AN/AO  nee370  B  20  08:42  SA110AL  nee31\  B  20  08:44  SA110 CE, CF  nee312  B  20  08:48  Chapter 2. OBSERVATIONS AND REDUCTIONS  TABLE 2.4: (Continued)  16  Observing log for the TI data set Observation Number  Filter  Exposure (seconds)  Universal Time  SA110 CC, CD  nee373  B  20  08:51  SA110 CC, CD  nee374  V  20  08:53  SA110 CE, CF  nee31S  V  20  08:55  SA110 CE, CF  nee376  I  10  08:58  SA110 CC, CD  neehll  I  10  09:00  SA112 805  nee378  I  10  09:06  SA112 805  nee379  V  20  09:09  SA112 805  nee3S0  V  10  09:12  SA112 805  «e<»381  B  10  09:15  Object  Notes to Table 2.4: A & L Stnds: Alcaino and Liller (1984) Standard stars u, w, t and x PG0918: Palomar Green Survey, Observatory standard (I. Thompson, Private Communication) PG1323: Palomar Green Survey, Observatory standard (I. Thompson, Private Communication) PG1323 - A: Noted to be variable; not used in the calibration. SA110: Standard area, Landolt (1973) SA112: Standard area, Landolt (1973) SA113: Standard area, Landolt (1973)  Chapter 2. OBSERVATIONS AND REDUCTIONS  17  TABLE 2.5: Observing log for the TEK data set  Object  Observation Number  l-'ilter  Exposure (seconds)  Universal 1  May 22:23 Median of 5 B dome flats  None  B  60  < 22:2!  Median of 5 V dome flats  None  V  10  < 22:2!  nnl —• nnl  U  10, 15, 20, 25, 35, 50, 65  U dome flats  22:25 ->2;  \i:>hl of May 27 2H A & L Stnds.  ««299  V  60  23:38  A & L Stnds.  nn300  V  300  23:44  A & L Stnds.  nn301  V  300  23:52  A & L Stnds.  n«302  V  30  23:59  A & L Stnds.  nn303  B  30  00:02  A & L Stnds.  nn304  B  300  00:04  A & L Stnds.  nn305  B  300  00:12  A & L Stnds.  nn306  U  600  00:29  Program field  nn307  U  300  00:46  Program field  nn308  B  15  00:54  Program field  «n309  B  15  00:56  Program field  nn310  B  15  00:57  Program field  nn311  B  15  00:58  Program field  nnhYl  B  15  00:59  Program field  ««313  B  15  01:00  Program field  nn314  B  15  01:00  Program field  n«315  B  15  01:01  Program field  n/i316  B  15  01:01  Program field  nn317  B  15  01:02  Program field  nn318  B  15  01:02  Program field  n«319  B  15  01:03  Program field  n«320  V  15  01:05  Program field  «n321  V  15  01:08  Program field  nn322  V  15  01:08  Program field  n/i323  V  15  01:09  Program field  n«324  V  15  01:10  Program field  nn325  V  15  01:11  Chapter 2. OBSERVATIONS AND REDUCTIONS  TABLE 2.5: (Continued) Object  18  Observing log for the TEK data set '  .  Filter  Exposure (seconds)  Universal Time  Program field  nn326  V  15  01:20  Program field  ««327  V  15  01:21  Program field  «/i328  V  15  01:21  Program field  nn329  V  15  01:22  Program field  «n330  V,  15  01:23  Nif-hi of May 2V Jo G162-66  ««355  U  10  23:13  G162-66  nn356  U  10  23:14  G162-66  «n357  B  5  23:16  G162-66  n«358  B  5  23:17  G162-66  ««359  V  5  23:19  G162-66  n«360  V  5  23:21  Program field  ««361  V  20  23:43  Program field  nn362  V  20  23:51  Program field  ««363  V  20  23:53  Program field  nn364  B  20  23:55  Program field  /m365  B  20  23:56  Program field  n«366  B  20  23:58  Program field  «n367  B  20  23:59  Program field  ««368  B  20  00:00  Program field  nn369  U  300  00:03  Program field  nn370  u  300  00:10  Cen A 3/6/9  /i«378  V  5  01:24  Cen A 3/6/9  ww379  V  5  01:26  Cen A 3/6/9  nn380  B  5  01:27  Cen A 3/6/9  nn381  B  5  01:28  Cen A 3/6/9  n«382  U  100  01:28  Cen A 3/6/9  nn383  u  100  01:33  Cen A 2/5/8  nn384  u  100  01:39  Cen A 2/5/8  /in385  u  100  01:45  Cen A 2/5/8  n«386  B  5  01:49  Cen A 2/5/8  nn387  B  5  01:51  Cen A 2/5/8  nn388  V  5  01:52  Cen A 2/5/8  nn389  V  5  01:54  19  Chapter!. OBSERVATIONS AND REDUCTIONS  TABLE 2.5: (Continued)  Observing log for the TEK data set Observation Number  Filter  Exposure (seconds)  Universal Time  PG1323  nn390  V  20  02:02  PG1323  nn391  V  20  02:04  PG1323  nn392  B  20  02:06  PG1323  «n393  B  20  02:08  PG1323  n«394  U  200  02:10  PG1323  nn395  u  150  02:15  SA113 492/493/495  nn435  u  100  09:41  SA113 492/493/495  nn436  u  100  09:43  SA113 492/493/495  nn437  B  5  09:47  SA113 492/493/495  /w438  B  5  09:48  SA113 492/493/495  «n439  V  5  09:49  SA113 492/493/495  «n440  V  5  09:50  Feige 108  n«441  V  5  09:53  Feige 108  nn442  V  5  09:54  Feige 108  ««443  B  5  09:55  Feige 108  «n444  B  5  09:58  Feige 108  nn445  u  20  10:01  Feige 108  nn446  u  20  10:02  SA112 805 & 822  nn447  u  20  10:09  SA112 805 & 822  ««448  u  40  10:10  SA112 805 & 822  «/i449  B  5  10:11  SA112 805 & 822  nn450  B  5  10:12  SA112 805 & 822  n/j451  V  5  10:13  SA112 805 & 822  n«452  V  4  10:14  SA112 805 & 822  n«453  V  4  10:15  Object  Notes to Table 2.5: A & L Stnds: Alcaino and Liller (1984), Standard stars u, w, t and x G162-66: Giclas star, see Landolt (1973) Cen A (NGC 5128), See Graham (1981) PG1323: Palomar Green Survey, Observatory standard (I. Thompson, Private Communication) Feige 108: See Landolt (1973) SA113: Standard area, Landolt (1973) SA112: Standard area, Landolt (1973)  Chapter  3  CALIBRATION 3.1 Modus  Operandi  To place the photometry system of this thesis onto a standard system, transformation equations need to be derived. If afiltersystem had been used for the observations which closely resembled the standard system, then one would expect the magnitudes and colours of the systems to disagree by only zero point offsets. As thefiltersdiffer, it is necessary to calculate colour terms. Below is a synopsis of the photometry calibration.  The DAOPHOT task PHOT is used to evaluate aperture magnitude values of primary standards, a primary standard being a star with a hitherto known magnitude. PHOT locates th centre of a stellar image, and integrates thefluxwithin a radius, returning afistof magnitudes evaluated at different radii. Errors aside, the magnitude should converge as the radius increase and all thefluxin the PSF is counted. The convergence value is the aperture magnitude of the s Having aperture magnitudes for standard stars allows for transformation equations to be derived to convert aperture magnitudes to a standard system. In deriving the transformation equations i is important to have standards with a wide range in colour to help tie down the colour term the transformation equations. The next step in the calibration is to perform aperture photometry on secondary standards; uncrowded stars in thefieldto be calibrated. For this a single fram is used so as to avoid any ambiguity over corrections for air mass and exposure time. Having aperture magnitudes for the secondary standards will then allow the standard magnitudes of these stars to be derived.  There are both advantages and disadvantages to using long or short exposure frames to obtai secondary standards. Short exposure frames will be less crowded, but at the expense of having lower signal to noise ratio. Long exposure frames have a higher signal to noise ratio, but are m crowded. It is not possible to increase the signal to noise in the data frames, but it is possibl reduce the effects of crowding. This can be achieved by subtracting out non-secondary standards from the frame. The secondary standards are then used to determine the new zero points in t transformation equations for profile-fitting magnitudes to standard system magnitudes. N.b. The colour terms in the transformation equations for aperture and profile-fitting magnitudes will not change, only the offsets change as there is a different normalization. These equations can the be applied to all the profile-fitting magnitudes to determine standard magnitudes. 20  Chapter 3. CALIBRATION  21  3.2 Aperture Photometry Transformations for the T I Data  Primary standards were observed on both the nights of May 22/23 and May 23/24 1988, ref to Table 2.4. These stars were used to derive the parameters in the transformation equations -V) = a -v ) = ap  (V  + /3(b - v )  ap  ap  7 a p  -0 — €a  + S(B-V)  ~\~  P  ap  C(^ap  (3.1)  ^ap)-  Many of the frames used had 'multiple exposures', where the star(s) had been observed, the telescope shifted, and then reobserved on the same frame. Where this had been done, aperture magnitude values measured were in excellent agreement with each other, and the value adopted was the mean of the measured values. The aperture magnitude values for the primar standards were corrected to an exposure time of one second and an air mass of zero atmosphere This is the normalization for exposure and air mass that I adopt for the remainder of this thes The extinction coefficients used to correct for air mass were those as used by Landolt (1983) the Cerro-Tololo Inter-American Observatory (CTIO), these values should be appropriate to Las Campanas due to the similar climatic conditions and altitudes of the observatories. The values used were Q  y  h  = 0.172 =  0.111  (3.2)  kc, = 0.087,  where Q , k\ and £9 are the extinctions in V, (B — V) and (V — I) respectively. y  The parameters and errors of thefitsin this chapter were determined using OPDATA (Bennett 1989). No significant second order colour terms were found in thefits,and such terms were droppedfromfurther consideration. After an initialfitto the data, the deviation fro thefitwas found for each point. A cleaning procedure was employed where any point having a deviation greater than twice the standard deviation of the distribution of all the points was rejected, on the premise that the point was an outlier and did not belong to the distribution. the distribution is normal, it is expected that approximately one point in twenty will be unfairly rejected by this procedure. The standard deviation of the distribution is calculated with outliers included, and so is larger than the intrinsic standard deviation of the distribution, thus decreasin even further the chance of a point being unfairly rejected. The data werefittedagain, excludin  22  Chapter 3. CALIBRATION  the outliers, which are shown as open circles (Figures 3.1 and 3.2). Thefitsobtained are (73 _ V)  =  0.583 ± 0.004 + 1.211 ± 0.006(b  (V - vap) = -(8.6 ± 1.5) x 10  -2  (V — J)  =  -  ap  v ) ap  + 9 x 10~ ± 1.25 x 10~ (B - V) 2  4  0.876 ± 0.008 + 0.97 ± 0.01(v  -  ap  (3.3)  i ) ap  for the night of May 22/23, and (B -V) = 0.56 ± 0.01 + 1.205 ± 0.015(6 - v ) ap  op  (V - v ) = -(9.36 ± 1) x 10 - 7.1 x 10 ± 1 x 10 (5 - V) -2  -3  (3.4)  _2  ap  (V -I) = 0.885 ± 0.008 + 0.98 ± 0.01(u - i ) ap  ap  for the night of May 23/24.  The parameters in thefitsare consistent between the two nights at a one sigma level, with th exception of the zero point in the (B — V) transformation, which is only marginally inconsistent. As the programfieldframes being used for the secondary standards came from both nights of observation (§2.5.2), I decided to combine all of the primary standards together and obtain transformation equations applicable to aperture photometry from both nights. An initialfitwas performed on the combined data set and the aforementioned cleaning procedure employed. A fit was performed on the cleaned data (Figure 3.3), and thefitsobtained were (B-V)  = 0.58 ± 0.01 + 1.19 ± 0M(b  ap  (V - v ) = ap  -(9 ± 1) x lO  -  v ) ap  (3.5)  -2  (V — I) = 0.880 ± 0.006 + 0.975 ± 0.008(u - i ), ap  ap  with the colour term being dropped from the (V — v ) transformation; it being consistent with zero when the errors are considered. That the colour term is small is an indication that the V filter used for these observations is congruous with the Johnson system V filter. ap  3.3 A L L S T A R Photometry Transformations for the T I Background Field V Calibration  To calibrate the TI V backgroundfield,a colourless transformation equation was employed, viz.  (V-v ) als  =  7 a l s  + 6(B-V),  (3.6)  Chapter 3. CALIBRATION  23  with 8 being set equal to zero, in accordance with equations 3.5. I determined 7^ in two ways. Thefirstwas to select nineteen secondary standards in the backgroundfieldand perform aperture photometry on them using the shorter exposure B and V frames, «ee318 and nee320. The backgroundfieldwas uncrowded, and it was deemed unnecessary to subtract out neighbours. Aperture magnitude values obtained were normalized, the extinction coefficients used being those given in equation 3.2. Applying equations 3.5 to the aperture magnitudes, the Johnson B and V magnitudes for the nineteen secondary standard stars were obtained. The nineteen secondary standards were identified in the ALLSTARfileresulting from the standard reduction on the TI V backgroundfieldwhich had been binned 2 x 2 (see §6.2) and a plot was made of (V — v i ) against (B — V), see Figure 3.4. Two points were arbitrarily rejected on the premise that they were considered to be outliers, and not to belong to the intrinsic distribution values. The mean of the remaining seventeen points gave a 3  7ais = (V-v ) ah  = 4.81 ± 0 . 0 1 ,  (3.7)  with the error being the standard error, given by standard error =—.  ,  (3.8)  where a is the standard deviation of the distribution and n is the number of points in the distribution. The second way in which j i was determined was to repeat the above procedure, but using the longer exposure B and V frames, nee3\9 and nee321, for the aperture photometry. Again, it was deemed unnecessary to subtract out neighbours. A plot of (V — v i ) against (B — V) is shown in Figure 3.5. Two points were dropped from consideration in evaluating the mean value of 7 /j. The mean of the remaining seventeen values gave a s  a s  a  7 a / s  = (V-  the error being the standard error.  v ) = 4.787 ± 0.006, als  (3.9)  Chapter 3. CALIBRATION  24  The two values for 7^ are in reasonable agreement, as anticipated. The latter value of 7^ has considerably less scatter in it, as expected from signal to noise considerations, and I adopte this value for the calibration of the binned background V fields. / Calibration  The / backgroundfieldgoes deeper than the V backgroundfield.To calibrate the / background field a colourless transformation equation, 1  = iais + Vats,  (3.10)  is required, otherwise only stars found on both the V and / frames can be calibrated by use of equations 3.5. Aperture photometry magnitudes were evaluated on both the short and long exposure / frames, nee322 and «ee323 respectively, for the same nineteen secondary standard as used in the V backgroundfieldcalibration. The frames were deemed sufficiently uncrowded as not to require neighbour subtraction. Aperture magnitude values obtained were normalized, using the extinction coefficients given in equation 3.2. The aperture photometry transformation equations for the TI data, equations 3.5, were applied to obtain / magnitudes for the nineteen seconda standards, with the short and long exposure aperture magnitudes being kept distinct.  The nineteen secondary standards were identified in the ALLSTARfilewhich had resulted from a standard reduction on the TI / backgroundfield,which had been binned 2x2 (see §6.2), and a plot was made of (/ — i i ) against (V — I), see Figures 3.4 and 3.5. a  s  A least squaresfitto the data points showed that there was an insignificant colour term. Tha the colour term was negligible indicated that thefilterused for the observations was congruous with the Kron-Cousins / filter.  Using the short exposure (60 second) backgroundfieldI obtained, after the rejection of thre outliers, V a l s  = (I - i ) = als  3.74 ±0.01,  (3.11)  the error being the standard error. Similarly, the long exposure (450 seconds) TI / backgroun field led, after the rejection of three outliers, to  25  Chapter 3. CALIBRATION  rials = (I-iais)  =  (3.12)  3-756 ± 0 . 0 0 5 .  The two values were consistent with each other at the one sigma level, but the latter, obtained using stars with a higher signal to noise, has less scatter. I adopted the latter value for the calibration of the binned TI / background fields. 3.4 A L L S T A R Photometry Transformations for the T I Program Field 3.4.1 Secondary Standards  Seventeen uncrowded stars were selected on the TI programfieldframes to be secondary standards. The calibration of these stars, to a standard system, was achieved in two ways; by using both long and short exposurefieldframes.  To reduce the effects of crowding, a standard reduction was performed on each of the above frames to produce frames that contained only the seventeen secondary standards and anything which had been missed by the standard reduction. One must realize that profilefittingis never perfect: shot noise will be left where a star has been subtracted out, due to the increased number of photons in a stellar image. The method of neighbour subtraction does have its limitations.  Aperture photometry was performed on the secondary standard stars, and the aperture magnitudes thus obtained were normalized, the extinction coefficients used being those given in equation 3.2. Two stars were lost in the short exposure calibration, due to the B frame being displaced with respect to the V and / frames. A comparison of the aperture magnitude values obtained for thefifteenstars in common is shown in Figure 3.6. One sees good agreement for the brighter secondary standards, with the scatter increasing at the fainter end. Standard magnitudes were calculated by use of equations 3.5. 3.4.2 Calibration of the TI Program Field Frames  The seventeen secondary standards from §3.4.1 were identified in the ALLSTARfilesresulting from standard reductions on the TI B, V and / programfields.The transformation equations for ALLSTAR photometry are V) = a  als  (V  v  (V  I) = tais  als)  +  lA9±0.01(b -v ) als  als  (3.13)  — Tals  + 0.975 ± 0.008(y  o/s  -  i  a l s  ).  26  Chapter 3. CALIBRATION  Values for a i , j i and were determined for each secondary standard star, and plots made of the value against the appropriate colour index, (Figures 3.7 and 3.8). Outliers deemed not to belong to the distribution, and excluded from the fit, are indicated with open circles. One expects no colour gradient if the colour coefficient is correct; this is the case. a  s  a s  Values obtained using the short exposure secondary standards are a^7  = 2.13 ±0.01  7 ^ 7 = 4.83 ± 0.02  (3.14)  e^7 = 0.17 ±0.005, whilst using the long exposure secondary standards gives OaTs^ 2.13 ±0.01 77J7  = 4.82 ± 0.005  (3.15)  tai; = 0.18 ±0.01.  With errors taken into account, these two sets of values are consistent with each other at the o sigma level. I adopted the latter set of values due to the generally smaller errors, which result in the transformation equations (B-V)  = 2.13 ± 0.01 + 1.188 ± 0.01(6 - v ) a/a  a}s  (y-v )  = 4.820 ±0.005  (V  0.18 ± 0.01 + 0.975 ± 0.008(u, -  als  — I) =  (3.16) a s  i ) als  for ALLSTAR photometry to standard system photometry on the TI program frames. 3.4.3 Calibration of the Binned TI Program Field Frames  To calibrate the binned V and / frames, which were used to construct the luminosity function colourless transformation equations, see §3.3, were used. The seventeen secondary standards from §3.4.1 were identified in the ALLSTAR list which had resulted from a standard reduction on the V and / frames that had been binned 2x2 (see §6.2). Figures 3.9 and 3.10 show plots (I - ials) against (V - I) and (V - v i ) against (B - V). Open circles indicate points which were arbitrarily considered to be outliers, and were dropped from consideration. a  s  27  Chapter 3. CALIBRATION  Using the shorter exposure programfields,yields the transformation equations 7afc = U = 3-30 ± 0.01 rials = (V - v i ) = 4.444 ± 0.006,  \D.LI)  a s  while using the longer exposure programfieldsgives ~f ls = (I-^als)  = 3.274 ± 0 . 0 0 4  a  (3.1s)  rials = (V - v ) = 4.425 ± 0.005. als  The values of j i  a s  and riai are reasonably consistent with each other at around the two s  sigma level. The standard deviations of the distributions obtained using the longer exposure frames are smaller than those for the shorter ones, and therefore have better defined means. I adopted the values obtained using the longer exposure frames for the calibration of the binned program field frames.  3.5 Error Propagation in the Transformation Equations Using the error returned along with the ALLSTAR magnitude in DAOPHOT, it is possible to propagate errors and determine the errors in the calibrated colours and magnitudes. Writing (B-V)  = a + 3(b -v ), als  (3.19)  aU  we have, by the addition of errors in quadrature,  = l + l(b -v ) a  '°IB-V)  a  als  aU  •  ( - °) 3  2  Now, Bevington (1969) gives that if x = ±auv,  (3.21)  then  a = a v al + a VCTJ + 2a^uva\ 2  2  2  x  with a  uv  v  (3.22)  being the covariance, defined as o  lim  1 -T[( -u)( -v)]. Ul  n —> oo n  Vt  (3.23)  Writing x = B(b u - v ) a  aU  (3.24)  28  Chapter 3. CALIBRATION  and a=1 (3.25)  u= B V  = (b  als  -  V l ), a s  the covariance in the error term can be ignored as no correlation between  /3 and ( b i a  s  —  vi) a  s  is expected. Hence  Using the addition of errors in quadrature I obtain °(B-V)  =  +P  (°L + O + (bais ~ vah) ai,  2  (3-27)  2  which is readily evaluated. The errors for the other magnitudes and colours can be propagated similarly. The errors shown in the colour-magnitude diagrams and colour-colour diagrams of chapter 4 are calculated as above. They essentially reflect two components, one component being a calibration error, which can be regarded as an error in coordinate labels, the other error being an intrinsic error which, as expected, increases for fainter stars.  3.6 Aperture Photometry Transformations for the T E K Data Observations of primary standard stars were made on the night of May 29/30 1990 (see Table 2.5). Using these standards, transformation equations from aperture photometry to Johnson photometry were derived. This was achieved by performing aperture photometry on the standard stars and normalizing for air mass and exposure time. The extinction coefficients used were  Q = 0.172 y  k\ = 0.111  (3.28)  fa = 0.290, where Q , k\ and £3 are the extinctions in V, (B — V) and (U — B) respectively. Thefirsttwo of y  these values are from Landolt (1983), while the third is from Landolt (private communication). The (U — B) extinction coefficient was adopted after finding a large discrepancy between the value reported by Landolt (1983) and that in the CTIO user manual. An initialfitwas performed upon the data; points which had a deviation from the initialfitgreater than two standard deviations  Chapter 3. CALIBRATION  29  were deemed to be outliers and dropped from consideration (Figure 3.11) . The transformation equations thus derived are ( B - V ) = 0.379 ± 0.007 + 1.26 ± 0.02(b  ap  - v ) a p  (V - v ) = -0.80 ± 0.01 + -0.01 ± 0.01(5 - V ) a p  (U - B ) = —1.55 ± 0.02 + 0.90 ± 0.01(u  ap  (3.29)  -b ). a p  The colour term in the (V — v ) transformation equation was considered small enough to be a p  neglected. This is a result of having a Vfilterthat is congruous with the Johnson Vfilter.The colour term was dropped and the mean of the points, excluding the outliers, was found, giving the transformation equation (V - v ) = -0.81 ± 0.01. a p  (3.30)  3.7 A L L S T A R Photometry Transformations for the T E K Data 3.7.1 Secondary Standards Frames /m362, nn366 and n«370 are short exposure TEK program field frames, acquired on the same night as the standardfields,so it is legitimate to apply equations 3.29 to normalized aperture photometry on these frames. Thirteen secondary standard stars were chosen which were common to all three short exposure program field frames, as well as to the program field frames from the night of May 27/28 1990, which were used to build the program frames. The frames were deemed to be too crowded for aperture photometry and so a standard reduction was performed on each frame to subtract neighbours. Johnson magnitudes were derived for these stars using the transformation equations derived in §3.6.  3.7.2 Calibration of the T E K Program Field Frames To derive the ALLSTAR to Johnson transformation, it is necessary to redetermine the zero points in the aperture transformation equations, as in §3.4.2. The thirteen secondary standards were identified in ALLSTAR file resulting from a standard reduction on the frame and values were derived for ot i , ^ah a s  and e i and were plotted against a colour index (Figure 3.12). Open a s  circles indicate points deemed as not belonging to the distribution, which were dropped from  Chapter 3.  CALIBRATION  30  r  consideration. Taking the mean and standard error of the three distributions gives ( B - V ) = -0.28 ± 0.02 + 1.26 ± 0.02(b  - v )  al3  (V - v  a h  a h  ) = 1.32 ±0.01  (3.31)  ( U - B ) = 0.63 ± 0.02 + 0.90 ± 0.01(««j - b ) a  als  Errors for Johnson magnitudes and colours were propagated and evaluated as outlined in  §3.5.  Chapter 3. CALIBRATION  31  I 03  J  I  I  l_  1.5  1  r  n — n  1  1  1—-i  I  I  I  r  i  i  i  i  i  i  r  -  -0.05 I  -0.1 -0.15 -1  I L _  1  0.5  1.5  (B-V) l  1  I  \  !  I  L_  I  _J  (v -i ) a p  a p  FIGURE 3.1. Aperture photometry transformations for the night of May 22/23 1988 Open circles indicate points which had deviations greater than 2a from the fit to all the points, and were dropped from consideration, while the solid lines indicate the fits to the remaining points. The subscript ap refers to normalized aperture photometry magnitudes, while upper case denotes the standard system.  32  Chapter 3. CALIBRATION  l  2  l  I  1.5 I  1 0.5 0 _J  |  I  |_  0.5  H  1  1  T  -i  1  1  1  i  i  i  1.5  1  1  1  r  i  ,  r  n  -0.05  -0.1 -0.15 _i  i  i  i  i  i  J_  i  i  i  i  0.5  i  _i  i  i  i  1.5  i  i  i  (B-V) ~ i  1  1  r  I  -1  J | | L_  _|  ( ^ap  |  |  i  i ap )  FIGURE 3.2. Aperture photometry transformations for the night of May 23/24 1988 Open circles indicate points which had deviations greater than 2a from the fit to all the points, and were dropped from consideration while the solid lines indicate the fits to the remaining points. The subscript ap refers to normalized aperture photometry magnitudes, while upper case denotes the standard system.  Chapter 3. CALIBRATION  33  (v -i ) a p  a p  FIGURE 3.3. Aperture photometry transformations for the nights of May 22/23 and 23/24 1988 Open circles indicate points which had deviations greater than 2a from the initial fit to all the points, and were dropped from consideration, while the solid lines indicate the fits to the remaining points. The subscript ap refers to normalized aperture photometry magnitudes, while upper case denotes the standard system.  34  Chapter 3. CALIBRATION ~i  1  _j  i  r  4.85  ri ^  4.8  4.75 J  |  |_  I  I  L_  0.8  0.6  1.2  i_  1.4  (B-V) 3.85  -i—n  r  ~i  1  r  _i  i  i_  3.8  3.75  0.8  1.2  1.4  1.6  (V-I)  FIGURE 3.4. TI background field ALLSTAR zero points (short exposure frames) Open circles indicate points dropped from consideration, and the dashed line indicates the mean value of the remaining points. The aperture magnitude values, used to derive the standard magnitudes, were determined using the frames nee318, nee320 and nee322. The subscript als refers to profile-fitting magnitudes, while upper case denotes the standard system.  Chapter 3.  35  CALIBRATION  4.85 4.8 4.75 -1  I  -J  l_  I  L_  —J  I  L.  0.8  0.6  —I  I  1.2  1_  1.4  (B-V) 3.85  ~i  1  r  3.8  3.75  0.8  1.2  1.4  1.6  (V-I)  FIGURE 3.5. TI background field ALLSTAR zero points (long exposure frames) Open circles indicate points dropped from consideration, and the dashed line indicates the mean value of the remaining points. The aperture magnitude values, used to derive the standard magnitudes, were determined using the frames nee319, nee321 and nee323. The subscript als refers to profile-fitting magnitudes, while upper case denotes the standard system.  36  Chapter 3. CALIBRATION 0.1  —j  1  1  r  -J  I  I  l_  ~i  1  1  -I  I  i  r  cr(Ab )=0.037 ap  0.05 0 -0.05 -0.1  -J  I  I  L  17  17.5  _l  '  18  I  I  18.5  •ap n  0.1  1  1  r  n  i  i  I  i  i  i  \  i  i  p  ~i  1  r  i  i  i_  cr(Av- ) = 0.041 ap  0.05 0 -0.05 -0.1  _L  16.5  17  — i — i — i — i — i  17.5  18  i  i  18.5  Vap ~i  0.1  1—  i  cr(Ai ) = 0.050 ap  0.05 0 -0.05 -o.i  F.j  16  i  17  i_  18 l  ap  FIGURE 3.6. Comparison of aperture photometry between frames with long and short exposure times Aperture magnitudes were normalized to an exposure time of 1 second and an air mass of zero. The difference between the magnitudes, A , is in the sense (long - short). The dashed lines indicate the mean values of the distributions, while their standard deviations are given in each panel.  37  Chapter 3. CALIBRATION i  2.25  n  r  i  r  r  2.2  3  2.15 2.1 2.05 _i  0.8  0.6  i  1.2  (B-V) 4.95  i  ~\  r  r  ~i  r  4.9 4.85 4.8 4.75  _|  L.  _J  0.6  I  |_  0.8  1.2  (B-V) -i  r  0.6 0.4 (JJ  0.2  0.8  1.2  (V-I) FIGURE 3.7. TI program field ALLSTAR zero points (short exposure frames)  Open circles indicate points deemed as outliers and dropped from consideration, and the dashed lines indicate mean values of the remaining points in each of the three panels. The zero points were derived using secondary standards on the TI program field frames nee206, nee209 and nee309. Upper case letters denote the standard system.  38  Chapter 3. CALIBRATION  2.25  ~i  r  _i  1  -1  1  r  ~i  r  _J  I  L_  _j  i_  2.2 "3  2.15 2.1 2.05 L_  1.2  0.8  0.6  (B-V) 4.95 4.9  J  -1  r  _J  |_  1  1  -  4.85 4.8 4.75 0.6  _|  L_  0.8  1.2  (B-V) r  1  1  r  0.6 0.4 0.2  0.8  1.2  (V-I) FIGURE 3.8. TI program field ALLSTAR zero points (long exposure frames) Open circles indicate points deemed as outliers and dropped from consideration, and the dashed lines indicate mean values of the remaining points in each of the three panels. The zero points were derived using secondary standards on the program field frames nee207, nee211 and nee311. Upper case letters denote the standard system.  39  Chapter 3. CALIBRATION 3.45  "i  r  ~i  i  r  3.4 »  3.35  1,  3.3 3.25 3.2  ° I  1.2  0.8 (v-i)  4.7  n  -i  r  r  4.6 «  |  4.5 4.4  _I  0.6  L_  0.8  (B-V)  FIGURE 3.9. Colourless TI program field ALLSTAR zero points (long exposure frames)  Open circles were deemed as outliers and were dropped from consideration, and the dashed lines indicate mean values of the remaining points in each of the three panels. Aperture magnitude values were determined using the frames nee2Ql, nee2\ 1 and nee311. The subscript als refers to profile-fitting magnitudes, while upper case is used to denote the standard system.  40  Chapter 3. CALIBRATION 3.45  -|  ~i  ~r  r  ~i  r  3.4 « i,  3.35 3.3 3.25 3.2  _J  L _  1.2  0.8  (V-I) 4.7  ~i  r  ~i  r  ~ i — —  r  _i  i  4.6 4.5 4.4 0.6  i_  0.8  (B-V)  FIGURE 3.10. Colourless TI program field ALLSTAR zero points (short exposure frames)  Open circles were deemed as outliers and were dropped from consideration, and the dashed lines indicate mean values of the remaining points. Aperture magnitude values were obtained using the frames nee206, nee209 and nee309. The subscript als refers to profile-fitting magnitudes, while upper case is used to denote the standard system.  41  Chapter 3. CALIBRATION \  r  I 05  _|  1  |_  0.5  -0.5  (b -v ) ep  i 0  — .,_ 1  1  i  i  -0.8  — -  1  1  1  1  1  1  0  1  1  1  0.5  |  1  1  1  -  • "" r  1  1  0  1  •  O t  1  • ' _ • • ••  *  -0.82  i  1  i  -0.78  n.  i?  i  ap  1  *  —  •—-  ,  ,  i  (B-V)  cq  I  _J  L_  _j  |  |  L  J  |_  (u -b ) ap  ap  FIGURE 3.11. Aperture photometry transformations  for the TEK data set  Open circles had deviations greater than 2a from the fit to all the points, and were dropped from consideration, while the solid lines indicate the fits to the remaining points. The dashed line in the middle panel indicates the mean value of (V - v ) after the colour term is dropped in the fit. The subscript ap refers to normalized aperture photometry values, while upper case denotes the Johnson system. a p  42  Chapter 3. CALIBRATION -0.2  r  1  ~i  n  1  1  r  _|  |  |  L _  -0.4 -0.6 -0.8 _i  i  i  i_  0.5  (B-V)  c5  (B-V) -i  1  r  i  i  n  1  |  i  1  1  TV  r  0.9 0.8 0.7 0.6 j  I  |  1  0.2  i  0.4 T  i  1  0.6  |  |  L_  0.8  ( U—B)  FIGURE 3.12. TEK program field ALLSTAR zero points Open circles were deemed as outliers and dropped from consideration. The dashed line indicates the mean value of the remaining points in each of the three panels. Upper case denotes the standard system.  Chapter 4 COLOUR-MAGNITUDE AND COLOUR-COLOUR DIAGRAMS 4.1 Data Cleaning Profile-fitting photometry can produce spurious results, e.g. a stellar image may have an unidentified neighbour in its profile which could cause the fitting to be poor. ALLSTAR returns, along with a magnitude, an estimate of the error in the magnitude. As one would generally expect, this error increases for fainter magnitudes, as the signal to noise ratio drops. If the data are cleaned by imposing an upper cut off for the error, it is found that fainter stars, with errors in an acceptable range, are rejected, while brighter stars, with relatively large errors, are accepted. To circumvent this, I employed PHOTSELECT, a program written by G. Drukier at the University of British Columbia. PHOTSELECT works by fitting a curve to the lower envelope of an error diagram, binning in magnitude, and evaluating the mean error in each bin. A 'mean multiplier' is chosen, and stars are rejected if their errors are larger than the mean multiplier multiplied by the mean error in the bin. Figures 4.1 and 4.2 show the effect of cleaning data using PHOTSELECT; note how much better delineated the blue stragglers are in the cleaner diagram.  4.2 Production of Diagrams 4.2.1 TI Data The three program field frames which had been built for the TI data (see §2.5.2) had standard reductions performed upon them to yield ALLSTAR magnitudes for the stars. The finding thresholds used were determined by plotting a graph of detected objects against finding threshold, a steady rise in the number of detected objects is seen as the finding threshold is decreased, to a point where the shot noise in the sky starts being detected and the number of detections takes an abrupt upturn. The threshold values were chosen at the point where the sky shot-noise starts to be detected. Using a standard reduction, 544, 660 and 1229 stars were fitted on the TI B, V and / frames respectively. The frames had linear displacements applied to translate them to a common origin. (N.b. As the telescope and its configuration were the same for all observations, one would expect linear 43  44  Chapter 4. COLOUR-MAGNITUDE AND COLOUR-COLOUR DIAGRAMS  shifts in the X and Y coordinates only, i.e. scale changes or rotations are not expected, and were not allowed for.) The photometry files were cleaned to varying degrees (see §4.1) and local programs were utilized to match stars between frames and producefileswith calibrated magnitudes and errors. Only stars which were found on all three frames were considered as real helping to reduce spurious detections. Using mean multipliers of 2, 4 and 8 in the PHOTS EL routine, the number of stars common to the three frames were 301, 335 and 364 respective One always expects fewer stars to be matched between frames; at faint magnitudes, where th incompleteness is high, only the intersection of the samples will be counted as real. CMDs for (B — V), (B — I) and (V — I) versus V, were plotted, along with the errors in photometry (Figures 4.3, 4.4 and 4.5). It is apparent that the data just reaches the main sequ turn off, with the turnoff being poorly defined due to a lack of turn off stars and subgiants. With this caveat noted, the main sequence turn off is measured as lying at Vxo=18.2 ± 0.2, (B - V) =0.65 ± 0.03, (B - J) =1.49 ± 0.02 and (V - I) =0.82 ± 0.03, as indicated by arrows in Figures 4.3, 4.4 and 4.5. TO  4.2.2  TO  TO  T E K Data  The programfieldU, B and V frames were reduced in a similar manner to the TI program field frames. Initially a lot of faint stars were left in the frames. To help clean up the subtr frame, thefindingthresholds were dropped and the frames reduced again. This provided a cleane final frame and should improve the photometry by ensuring that neighbours do not contamina the stars beingfitted.The resulting ALLSTARfileswere processed as above and, using mean multipliers of 1, 4 and 8, the number of stars matched between the frames were 1022, 222 2553 respectively. A low value for the mean multiplier was found to give a well defined CMD Using a low mean multiplier is valid as I expect a lot of contamination from poorfits,due to t crowding on the frame. A (B — V) versus V colour-magnitude diagram is shown in Figure 4.2. It is apparent that the limit of reliable photometry is at approximately 17 5, the width of the main sequence increasing rapidly after this. Using the CMD to measure the main sequence tur off yields VTO—IS.3 ± 0.3 and (B — V) =0.6 ± 0.03, as indicated by the arrows in Figure 4.2. The magnitude and colour of the turn off are poorly defined due to the large scatter. ni  TO  4.3 Two-Colour Diagrams  The two-colour diagram considered in this work is a plot of the colour index (U — B) agains the colour index (B — V). If the stars in the diagram were perfect black bodies, at different  Chapter 4. COLOUR-MAGNITUDE AND COLOUR-COLOUR DIAGRAMS,  45  temperatures, then they would define a diagonal sequence, similar to the reddening locus, in t two-colour plane. The 'S' shape in the two-colour diagram is caused mainly by absorption in the Balmer continuum, which causes the stars to be non-blackbody radiators. 4.3.1 Reddening  For spectral types O through B9, the reddening vector in the two-colour plane is given by E(U E(B  -B) - V)  0.72 + 0.005E(J5 -  V)  (4.1)  (Mihalas and Binney (1981), page 186). The colour term on the right hand side of the equati makes only a small contribution and is usually neglected.for modest values of reddening. The cluster interstellar reddening can be determined by shifting the clusterfiducialalong the reddening vector until the early spectral types coincide with an unreddened luminosity class V fiducial. Early spectral types are used as these stars are sufficiendy hot that the metallicity does not affe their positions in the two-colour plane, unlike cooler stars. The stars which meet this requiremen in the current data are the horizontal branch (HB) stars. Although not luminosity class V, th HB stars suffice for the purpose. Figure 4.6 shows data from the TEK data andfiducialswhich have been reddened by E(B - V)=0™16, 0™21 and 0™25. There is a lot of scatter in the diagram and the HB stars do not appear to lie along a particular line. It has been suggested by Cacciari (1984) and Da Costa et al. (1981) that the reddening varies across the face of the cluster, but the magnitu of the effect proposed, AE(B — V) ~0 ™03, is too small to explain the scatter seen in the two-colour diagram.  Looking at Figure 4.2 one sees that the blue end of the HB starts to droop, and breaks of into what appear to be two branches. The HB of NGC 3201 does droop, as is clearly seen in Lee's Figure 7. With the paucity of stars on the HB, little can be said on its morphology. T bluest horizontal branch stars are critical in the determination of the reddening, and the spread in the HB causes the determination to be poor. No explanation is offered for the cause of thi effect; some of the horizontal branch stars were blended, though they all subtracted from the frame cleanly. To elucidate the reddening via this method, afieldis needed containing more horizontal branch stars; the smallfieldsof CCDs are not congenial to providing this.  The reddening may also be determined by the use offieldstars. Field stars that lie outside the galactic plane suffer the same absorption as the cluster. If thefieldstars are plotted on a t  Chapter 4. COLOUR-MAGNITUDE  AND COLOUR-COLOUR  46  DIAGRAMS  colour diagram, the reddening can be determined by working out the trajectories to translate the stars back onto the Population I luminosity class Vfiducial,this of course assumes that the stars are Population I luminosity class V stars. There were no obviousfieldstar candidates appearing on the CMD (Figure 4.2), and it was concluded that this method would not be of worth. The value adopted for the interstellar reddening for the rest of this thesis was 0 ^21 in accordance with the value given by Lee (1977), this value not being inconsistent with the determination of the reddening in this thesis. 4.3.2 Metallicity  Population JJ stars have a lower metallicity than Population I stars. The difference in metallicity means that the blue and ultraviolet atmospheric opacity in later type stars is less for Population II than for Population I. Metal-poor stars will consequently be found to he above the metal rich stars in the two-colour diagram. The position of the clusterfiducialin the twocolour diagram will thus allow for an estimation of the metallicity. Figure 4.7 shows the main sequence and sub giant stars from the TEK frame, as defined in Figure 4.1, plotted in the twocolour plane. The data were corrected for the effects of interstellar reddening, and fiducials for metal abundances of zero and the Hyades abundance plotted. The quantity 8(U — B) is defined as being the height at which the data lies above the Hyadesfiducial,at the intrinsic colour of (B — I0o=0 ™6. The data shown in Figure 4.7 are of poor quality. An attempt was made to determine the metallicity from the available data by determining the excess at other colours, and then interpolating in the table given by Sandage (1969) to determine the excess at (73-y)o=O 6. A mean line through the data was drawn by eye, this was judged as being superior to simply finding the mean in a bin, due to the large number of outliers scattered to higher values of (U — B) . Values for 6 were calculated and corrected as shown in Table 4.1 below. 6 0  ni  0  TABLE 4.1: Corrections of 8 to 8(0.6) 8  0.35 0.40 0.45 0.50  0.12 0.12 0.14 0.10  M 0.54 0.54 0.61 0.42  Correction 1.19 1.19 1.18 1.09  8(0.6)  0.14 0.14 0.17 0.11  A mean value of <5(0.6)=0.14 is obtained. To convert 8(0.6) to a value for metallicity, a calibration curve is required. I use the calibration curve as given by Richer and Fahlman (1986),  Chapter 4. COLOUR-MAGNITUDE AND COLOUR-COLOUR DIAGRAMS  47  which results in the value [Fe/H]iR=—0.83. The IR subscript indicates that the metallicity is on the scale as defined by Frogel, Cohen and Personn (1983). The derived value must be taken cum gran salis: this method of metallicity determination depends vitally on a good calibration and a reliable measure of reddening. A small error in 8(0.6) can lead to a large error in the metallicity determination due to the large gradient of the calibration curve. To obtain a metallicity of —1.4, a value of <S(0.6)=0.2 would be required, which is not inconsistent with the data. Other values determined for the metallicity of NGC 3201 are summarized in Table 4.2. TABLE 4.2: Metallicity values for NGC 3201 from previous studies  I.Fe/H| -1.2 + 0.1 -1.0 -1.40 ±0.06 -1.51 ± 0.02  Reference Lee (1977) Pilachowski et al. (1980) Zinn (1980a) Zinn (1980b)  Method Two-Colour Diagram High Dispersion Spectra Q39: Integrated Light Index Q39 (Calibration modified c.f. Zinn 1980a) S(W) and CMDs High Dispersion Spectra High Dispersion Spectra  -1.4 ±0.2 -1.19 + 0.05 -0.95 ± 0.2 -1.41 -1.34 -0.83  Da Costa et al. (1981) Gratton (1982) Pilachowski et al. (1983) Zinn and West (1984) Q39 Grattton and Ortolani (1989) High Dispersion Spectra This Work Two-Colour Diagram  The metal abundance of globular clusters was recently a topic of debate in the literature. Resultsfromearly photographic high dispersion spectra (Pilachowski et al. 1983) and photometric indices (Zinn 1980a and 1980b, Zinn and West 1984) were in disagreement. The source of the discrepancy, in the case of NGC 3201, was the definition of 'metallicity', which includes numerous metals, as compared to specific quantities like [Fe/H] (Pilachowski, private communication). Recent high dispersion spectra, obtained using CCDs (e.g. Gratton and Ortolani (1989)), give more accurate results, which are in good agreement with results from photometry and low dispersion spectra. For subsequent work I adopt the canonical value of [Fe/H]=—1.27 + 0.1 (Z=0.001), in accordance with ALA89. This value is in reasonable agreement with the results from contemporary high-dispersion spectra studies and photometry indices.  Chapter 4. COLOUR-MAGNITUDE  AND COLOUR-COLOUR  48  DIAGRAMS  4.4 Distance Modulus  To determine the distance modulus of a globular cluster, it is possible to make a comparison withfieldsubdwarfs which have measured parallaxes. Lutz, Hanson and Van Altena (1988), using 50fieldsubdwarfs with accurate trigonometric parallaxes and spectroscopic [Fe/H] determinations, derive the formula M (TT) = V  1.41 + 5.17(B - V ) - 0.94[Fe/H]  (4.2)  for a Population IIfiducial.The distance modulus of NGC 3201 was determined by translating the subdwarffiducialin V magnitude until it coincides with the cluster main sequence (Figure 4.8). Thefiducialwasfittedto the lower part of the NGC 3201 main sequence by eye, so as to avoid problems with outliers,fieldstars and the main-sequence turn off. A goodfitwas found for a distance modulus of (m — M) = 14.2 ± 0.15, in agreement with the value derived by Lee (1977). This value is consistent with the assumption that My(HB)=06; as shown in Figure 4.2. V  m  4.5 Isochrone Fitting 4.5.1 Synthetic Isochrones  The isochrones used in this study are those of VB85 and are not enhanced for oxygen. By mass, oxygen constitutes half of all metals, even at solar abundance, and is an important source opacity. Pilachowski, Sneden and Canterna (1980) report oxygen as being enhanced by a factor of two in NGC 3201 over normal solar abundances. Oxygen enhancement in globular clusters is believed to be a result of enrichment from Type II supernovae, possibly from Population III stars. At the time of writing there is some controversy over whether the use of oxygen enhanced isochrones will make any difference, as the other a-elements, i.e. Ne, Mg, Si, S and Ca may also have enhanced abundances. The use of oxygen enhanced isochrones leads to lower derived ages, whereas the use of isochrones, enhanced for all the a-elements, does not (Chieffi, Straniero and Salaris 1991). The isochrones used in the subsequentfittingsare those of VB85 and use canonical values of Y=0.2 for the helium abundance and a=1.6 for the mixing length parameter. The isochrones are applicable to a metallicity of [Fe/H]=—1.27, see §4.3.2, and are scaled to solar abundance.  Chapter 4. COLOUR-MAGNITUDE AND COLOUR-COLOUR DIAGRAMS  49  4.5.2 Deep CMDs  VB85 isochrones were superimposed on the CMDs, Figures 4.9, 4.10 and 4.12, the isochr being corrected for a cluster having a reddening of E(B - V)=0 ™21, and a V distance modulus of 14.2, in accordance with §4.3.1 and §4.4. The value of total to selective absorption used w 3.1 in accordance with Savage and Mathis (1979). The equivalent reddenings for other colour were taken from Fahlman et al. (1989) and are given by E(V — I) = IME{B  V)  E(B - I) = 2.36E(5  V).  (4.3)  To achieve coincidence of the synthetic isochrone with the zero age main sequence (ZAMS) the (B — V), V plane, (Figure 4.9), an offset in colour of 0 ™07 is required. There is a bo evidence (Hesser etal. 1987) to suggest that the VB85 isochrones are too blue by approximate 0 ™03 to 0 ™04, which may be attributable to the difficulties encountered in predicting (B — from the temperature. The value adopted for the reddening is also open to question. Cacciari (1984) and Da Costa, Frogel and Cohen (1981) report a variable reddening across the face o NGC 3201 of the order AE(B — V) ~ 0.03. This may, in part, account for some of the l colour offset required to obtain afit.Also, the possibility of undetected systematic calibration errors, such as caused by non-photometric conditions, can not be precluded. Figure 4.10 shows the synthetic isochrones in the (V — I), V plane. To obtain afitto th ZAMS, the synthetic isochrones were corrected for distance modulus and reddening as above, except the equivalent reddening in (V — I) was derived from equation 4.3 above. In this case a colour offset of — 0 ™01 in (V — I) was applied to the isochrones.  Figure 4.12 shows the synthetic isochrones on the (B — I), V plane. The isochrones we corrected for the distance modulus and had an equivalent reddening in (B — I) applied to the and a colour offset in (B — I) of 0 ™05 was applied to make the isochrones coincide with t ZAMS. 4.5.3 Age Determination from Deep CMDs  Cluster dating via isochronefittingis an art with many caveats. To be able to date a clus one needs to know accurately the distance modulus, the cluster reddening, the cluster metallici as well as having data with a good calibration. The input to the synthetic models is still very uncertain; no one is sure what abundance ratios to use, or whether the mixing length parame is correct for stars unlike the sun, and there appears to be problems in translating the isochro  Chapter 4. COLOUR-MAGNITUDE AND COLOUR-COLOUR DIAGRAMS  50  from the theoretical to observational planes as mentioned above. Given all of this, a quixotic attitude is also required. With the above pitfalls noted, I state that NGC 3201 appears to be around 15 ± 2 Gyr ol though I place little confidence in this value given the above and noting the lack of subgiant a main sequence turn off stars to allow for discrimination between the isochrones. The subgiant branch of the synthetic isochrones do not appear tofitthe data at all, though the brighter st probably are not all subgiants; if they were they would form a tight sequence instead of havin the large spread in colour and magnitude which is seen. The stars seen are probablyfieldstars. back up this claim by presenting the CMD for a backgroundfieldin Figure 4.11. The backgroun field is the one as discussed in chapter 6. The studies of ALA89 indicate an age of 16 + 3 G for NGC 3201. Figure 4.13 shows points defining the V, (B — V) main sequence of NGC 3201 from the studies of ALA89. It is immediately obvious that the morphologies of the isochrone and data are not similar, and the diagram clearly demonstrates the subjective nature of isochron fitting. For example, at V < 17 ™5 thefiducialof ALA89 appears tofitthe 16 Gyr isochro well, while for 17 ™5 < V < 18 ?5 thefiducialfits the 14 Gyr isochrone. For V > 19 ™0 the fiducial does notfitthe isochrones. In conclusion, I state that my age estimate based on synthetic isochronefittingis 15 ± 2 Gy This is in agreement with the value reached by ALA89, though I feel that the ALA89 data ma be questionable given the morphology of the derived CMD. I can state that the morphology of the CMDs from the present study agrees with the synthetic isochrones. 4.5.4 Shallow C M D s  VB85 isochrones were superimposed on the (B — V), V CMD resulting from the TEK data, (Figure 4.14). The isochrones were corrected for a cluster having a reddening of E(B — F)=0 21 and a V distance modulus of 14.2, in accordance with §4.3.1 and §4.4. To obtain coincidence of the giant branch of the synthetic isochrone with the data, it was foun necessary to apply a colour offset to the isochrones of 0 ™04. The possible reasons for this discussed in §4.5.2. It is hard to draw any solid conclusions from the data, due to the ill-defin main sequence. It appears that the age is around 16 ± 2 Gyr; this being consistent with the value derived in §4.5.3. ni  4.5.5 Combined C M D s  For completeness, a montage of the data in Figures 4.2 and 4.3 was plotted, (Figure 4.15) No allowance was made for possible zero point errors in the calibration or possible differential  Chapter 4. COLOUR-MAGNITUDE AND COLOUR-COLOUR DIAGRAMS  51  reddening effects; there is a possibility that thefiducialsequences are misaligned. No stars wer common to the two frames (Figure 2.1), and so it was not possible tofinda common zero po VandenBerg, Bolte and Stetson (1990) discuss and implement a technique for determining relative ages of clusters with similar metallicities. The technique does not require the cluster reddening or distance modulus to be known, it only requires that a well definedfiducialexist for the cluster. The technique works by utilizing the fact that the difference in colour betwee the main sequence turn off and the giant branch is an age sensitive parameter. One chooses fiducial cluster, and then aligns the clusterfiducialin colour by matching the main sequence offs, and in magnitude by matching a point slightly redder than the turn off. Differences in th position of the giant branch can then be interpreted in terms of age differences.  Figure 4.3 shows the cleaned CMD derived from the TI data. As can be seen, the main sequence turn off is barely defined, and there is no giant branch. Figure 4.2 shows the cle CMD derived from the TEK data. Although the giant branch is reasonably well defined, the main sequence turn off is very noisy and it would be foolhardy to extract afiducialsequence from it. An attempt was made to implement the technique, but was abandoned as being futi due to problems in defining the main sequence turn off. VandenBerg et al. (1990) state that uncertainty in (B — V) of + 0 ™02 leads to an uncertainty in age of ± 2 Gyrs, and so, cle a better defined sequence is required to implement this technique with any degree of credulity. For the rest of this thesis, an age of 16 Gyr is adopted for NGC 3201. 4.6 Binaries  Evidence, from photometry techniques, for binary stars in globular clusters has been found, e.g. McClure et al. (1985) found photometric evidence for binaries in E 3, as did Bolte (198 and Sato, Richer and Fahlman (1989) for NGC 288 and M12 respectively. On the other hand some studies show absolutely no evidence for the existence of binaries in other clusters. Binar stars are important in the dynamical consideration of globular clusters as they can provide a heating mechanism to reexpand the cluster core after gravothermal collapse.  From Figures 4.2 and 4.3 there is no obvious evidence for an equal-mass binary sequence Binary stars would manifest themselves as a sequence displaced upwards by approximately 0 ™7 magnitudes (for stars of equal mass). Whether one would expect tofindbinaries in the outer field is doubtful; mass segregation would tend to cause the heavier binaries to migrate to th cluster core. It has been shown by Spitzer and Mathieu (1980) that if initially 10 per cent o the stars in a globular cluster are binaries, then the density of binaries in the cluster core wou  Chapter 4. COLOUR-MAGNITUDE  AND COLOUR-COLOUR  DIAGRAMS  52  rapidly rise to 50 —> 100 per cent. Hence, to investigate whether NGC 3201 has a large nu of binaries, a study of the cluster centre, extending down the main sequence, would be benefici Photometric techniques for detecting binaries will underestimate the binary abundance as systems with large mass ratios, or with sub-luminous companions will be undiscovered. Griffi et al. (1988) found numerous binary systems, recognized by radial velocity variations, lying along the main sequence in the Hyades. Griffin et al. (1988)findthat approximately one third of the binaries are well separated from the main sequence, indicating that the detection efficiencies for binaries via photometry will be similar. It can be demonstrated (Richer 1991, preprint) that if the cluster were to be entirely composed of binaries, with a spectrum of mas ratios, then it may be hard to spot the binary candidates from the errors in the photometry. Gi this, it may be hasty to reject the presence of binaries in NGC 3201 based solely on the curr data. A more thorough study looking carefully at the main-sequence width and errors is needed 4.7 Blue Stragglers  Figure 4.1 clearly shows a sequence of blue stragglers. The blue stragglers are seen as extension of the zero age main sequence beyond the main sequence turn off. Many explanatio have been proposed for the formation of blue stragglers, and I point the interested reader to Leonard (1989), and references therein, for a contemporary discussion of proposed formation mechanisms. Figure 4.16 shows the location of the blue stragglers on the field.  If blue stragglers are coalesced binaries then they would be expected to be more massive th other cluster members and therefore maybe more centrally concentrated, due to mass segregatio effects, provided that their ages are longer than the cluster relaxation time. To test this hypothes a comparison of the central concentration of blue stragglers and subgiants is made on the TE frame. The subgiants are chosen to have the same magnitude range as the blue stragglers s that the two sets of stars are similarly incomplete, and any bias due to magnitude dependent crowding (Nemec and Harris 1987) is avoided. Using the selection criterion for blue stragglers that (0 ™2 < (B - V) < 0 6) and (16 < V < 18), there were 18 candidates, while 141 subgia were selected with the criteria that (0 ™6 < (B - V) < 1 ™0) and (16 < V < 18). A graph of cumulative fraction against the logarithm of the radial distance is shown in Figure 4.17. It appears, from Figure 4.17, that the blue stragglers are more centrally concentrated than the subgiants. To quantify this hypothesis the Kolomogorov-Smirnov (K-S) test is used. The K-S statistic will test the null hypothesis that the distributions are the same, against the alternativ that they are different. The K-S statistic is given by ni  Chapter 4. COLOUR-MAGNITUDE  D  AND COLOUR-COLOUR  max — OO < X < OO  DIAGRAMS  |s 0) - s Nl  for two cumulative distributions, where S N ^ Z ) and S N functions (Press et al. 1987). The function  2  N2  (z)|,  (Z)  53  (4.4)  are the two cumulative distribution  OO  (4.5) j=i  allows for the calculation of the significance of the result. The significance level of an observed value of D (as a disproof that the distributions are the same) is (4.6) where N is the number of points in the distributions (Press et al. 1987). For the distributions under consideration, D was calculated to be 0.407. The significance level of this result was calculated to be ~0.01, i.e. there is a one per cent chance that the two samples were drawn from the same distribution. From this I conclude that the blue stragglers in the TEK field are more centrally concentrated, at a confidence level exceeding 99 per cent, than the subgiants.  Chapter 4. COLOUR-MAGNITUDE  AND COLOUR-COLOUR  •  DIAGRAMS  54  (B-V)  FIGURE 4.1. Shallow CMD, with slight cleaning, showing divisions made to distinguish classes of star The upper left region defines the horizontal branch stars, the upper right region the giants and the lower division is for the main sequence and subgiant branch stars. Dashed lines indicate the RR-Lyrae gap as defined by Lee (1977). The data were cleaned using PHOTSELECT (§4.1) with a mean multiplier of 8. The square indicates Lee's unusual star (see chapter 5).  55  Chapter 4. COLOUR-MAGNITUDE AND COLOUR-COLOUR DIAGRAMS  NGC 3201 -  14  -  16  -  18  -  20  i  i  i  i  i  J  0 _  1  1  I  I  I  0.5 1  I  I  L  1.5  i  i  i  1  i  i  i  i  1  i  i  i  i  i  i  1  i  i  i  i  1  i  i  i  1  1  .  1  1  .  1  .  .  (B-V) FIGURE  4.2. Shallow CMD with cleaning  The data were cleaned using a mean multiplier of 1 in the PHOTSELECT routine (see §4.1). Arrows indicate the values chosen for the main sequence turn off. The dashed line indicates a magnitude of 14 8 , the expected magnitude of the horizontal branch, assuming My(HB)=0 ™6 and (m - M ) = 1 4 ™2. m  v  Chapter 4. COLOUR-MAGNITUDE AND COLOUR-COLOUR DIAGRAMS  56  (B-V) FIGURE  4.3. Deep (B - V ) CMD  The data were cleaned using a mean multiplier of 2 in the PHOTSELECT routine (see §4.1). Arrows indicate the values chosen for the main sequence turn off.  Chapter 4. COLOUR-MAGNITUDE AND COLOUR-COLOUR DIAGRAMS  —  -i-  |  1  1  p—i  -  - NGC 3201 16. -  -  18 -  1  1  1  1  1111  57  1  1  1 1r  i•  (B-I) FIGURE  4.4. Deep (B - I) CMD  The data were cleaned using a mean multiplier of 2 in the PHOTSELECT routine (see §4.1). Arrows indicate the values chosen for the main sequence turn off.  Chapter 4.  COLOUR-MAGNITUDE  AND COLOUR-COLOUR  DIAGRAMS  58  (V-I) FIGURE  4.5. Deep (V - I) CMD  The data were cleaned using a mean multiplier of 2 in the PHOTSELECT routine (see §4.1). Arrows indicate the values chosen for the main sequence turn off.  Chapter 4. COLOUR-MAGNITUDE  '• 1  1  AND COLOUR-COLOUR  ~~\—i—i—i—i—i—i—i—r  1 _  i—r  NGC 3201  -i -  —  59  DIAGRAMS  ° are HB stars o are giants  — -0.5 -  0I £3  0.5 -  •  _B  _ 1 -  , i , r 0 0.1 0.2 _:  1.5  y~  ^0.5  0  0.5  °"(C7-B)  ~i  i  i  i  1.5 i  i  i  i  I  i  i  i  r  ~i—r  J  L  (B-V) FIGURE 4.6. Two-colour diagram for the TEK program field, with luminosity class V fiducials Horizontal branch stars, giants and main sequence/subgiant branch stars, as defined in Figure 4.1, are indicated. The solid curve is the unreddened Population I luminosity class V two-colour fiducial (Mihalas and Binney 1981), whilst dashed curves represent the same fiducial reddened by the amounts 0 16, 0™21 and 0P25. The data were cleaned using a mean multiplier of 1 in the PHOTSELECT routine (see §4.1). m  60  Chapter 4. COLOUR-MAGNITUDE AND COLOUR-COLOUR DIAGRAMS  ~i—I—r  —i—i—|—i—i—i—i—|—r  NGC 3201 E(B-V) = .21  -0.5  0  o  i  la  0.5  1  i  0  0.1 0.2  0  i  J  i  L  0.5  (u-B)  1.5  a  a  i  1—i  r  ~i  1  ~\ 1 1  r  J  (B-V)  "i—r  I  I  I  r—r  l_  J  L  0  FIGURE 4.7. Two-colour diagram for the TEK program field, with minimum and maximum abundance fiducials The diagram contains main sequence and subgiant stars only, as defined in Figure 4.1. The solid line is the Hyades fiducial sequence, while the dashed line represents the minimum abundance fiducial (Sandage 1969). The colours of the stars have been corrected for reddening by amounts E ( B - F)=0 ™21 and E( U - B)=0 ™153. The data were cleaned using a mean multiplier of 1 in the PHOTSELECT routine (see §4.1).  Chapter 4. COLOUR-MAGNITUDE AND COLOUR-COLOUR DIAGRAMS  61  (B-V) FIGURE 4.8. Deep (B - V ) CMD of NGC 3201 with a subdwarf fiducial imposed  The subdwarf fiducial is from Lutz, Hanson and Van Altena (1988). The fiducial, corrected for a V distance modulus of 14.2, is shown as a solid line, while fiducials with distance moduli of 14.2 + 0 1 5 are shown as dashed lines. The subdwarf fiducial was appropriate for a metallicity of [Fe/H]=-1.27 and had been reddened by E(B - V)=0 ™21. m  62  Chapter 4. COLOUR-MAGNITUDE AND COLOUR-COLOUR DIAGRAMS  16  ~i  1  1  1  1  r  NGC 3201  E(B-V)=.21 Ay=0.651 14 16 18 Gyrs (B-V)  offset=.07  (m-¥) =13.55 0  18  20  -  22 -  24,  J  I  L  J  !  0.5 n  i  I  r  L  1.5 i  1  n— —r  r  :  0.1 0  J  I  L  J  L  (B-V) FIGURE 4.9. Deep (B - V ) CMD with VB85 isochrones imposed The isochrones are for stellar models with a=1.6, [Fe/H]=-1.27, and Y=0.20 and have been reddened by E ( B - V)=0 ™21 and corrected for the distance modulus of NGC 3201. The data were cleaned using a mean multiplier of 2 in the PHOTSELECT routine (see §4.1).  Chapter 4. COLOUR-MAGNITUDE  AND COLOUR-COLOUR  DIAGRAMS  63  (V-I) FIGURE 4.10. Deep (V-I)  CMD withVB85 isochrones imposed  The isochrones are for stellar models with a=1.6, [Fe/H]=-1.27, and Y=0.20 and have been reddened by E( V -1)=0 ™286 and corrected for the distance modulus of NGC 3201. The data were cleaned using a mean multiplier of 2 in the PHOTSELECT routine (see §4.1).  64  Chapter 4. COLOUR-MAGNITUDE AND COLOUR-COLOUR DIAGRAMS  16  i  11111  r  i  -  18  -  20  -  22  J  I  I  L_  J  I  [_  r  J  I  L  1.5 ~  i  I  i  r  ~i  i  r  j  i  L  (V-I) FIGURE 4.17. Deep (V - I) CMD for a background field near NGC 3201  Background fields in Vand /were reduced, and stars found were calibrated and matched. The errors shown are those returned by ALLSTAR, being added in quadrature for (V - I).  Chapter 4. COLOUR-MAGNITUDE  AND COLOUR-COLOUR  1  — i  -  1  11  n  1  1  r E(B-J)=0.4956  r  NGC 3201  16  65  DIAGRAMS  Ay=0.651 14 16 18 Gyrs (B-I) offset=.05 (m-M) =13.55 0  18 -  20 -  22 -  j  i_  1.5  2.5 i  J  L  j  i  i  L  j  F  i  r  i  i  i_  (B-I)  FIGURE 4.12. Deep (B - I) CMD with VB85 isochrones imposed The isochrones are for stellar models with «=1.6, [Fe/H]=-1.27, and Y=0.20 and have been reddened by E ( B - J)=0 ™496 and corrected for the distance modulus of NGC 3201. The data were cleaned using a mean multiplier of 2 in the PHOTSELECT routine (see §4.1).  Chapter 4. COLOUR-MAGNITUDE AND COLOUR-COLOUR DIAGRAMS  FIGURE  66  4.13. The (B - V ) CMD of ALA89  The open points are the mean values of colour taken at every 0 ™20. The solid lines are VB85 isochrones which have had the reddening and distance modulus adopted by ALA89 applied to them.  67  Chapter 4. COLOUR-MAGNITUDE AND COLOUR-COLOUR DIAGRAMS  1  12  i  1  i  i  i  i  i  i  i  i  i  i  i  i  1  i  i  i  NGC 3201 14  E(B-V)=.21 ~ A =0.651 14 16 18 Gyrs (B-V) offset=.04 (m-M) =13.55  \  v  J• •  -?<:  0  16 •t  i . - JL \  -  - |  18  —  —  -  '% 20  0  -  • '>i .''-'}^\:/-''-'- • • • x,  —  22^  ••  i  i  i  i  i  1  i  i  i  i  i  i  i  i  i  0.5  0 i  1  i  i  i  1.5 i  i  i  r  "l  r  J  L  1 I  r  I  I  I  I  (B-V) FIGURE 4.14. Shallow (B - V ) CMD with VB85 isochrones imposed The isochrones are for stellar models with a=1.6, [Fe/H]=-1.27, and Y=0.20 and have been reddened by E ( B - V)-0 ™21 and corrected for the distance modulus of NGC 3201. The data were cleaned using a mean multiplier of 1 in the PHOTSELECT routine (see §4.1).  68  Chapter 4. COLOUR-MAGNITUDE AND COLOUR-COLOUR DIAGRAMS  12  ~ i — r  n — i — i — | — i — i — i — r  ~i—1—r  NGC 3201 -  14  16  18 -  20 .1.  22 -  24.)  i  i  i  0.1 /  I  i  i  i  0  i  I  J  I  I  I  I  I  1  ~ i — r  j  I  0.5  I  I  I  I  I  1.5  2  ~i—i—i—i—i—r  i  I  i  i  L  _J  I  1  I  I  I  I  J  1  I  I  I  I  L  (B-V) FIGURE 4.15. Combined (B - V ) CMD for TI and TEK data sets The diagram is a montage of the data from Figures 4.2 and 4.3. The TEK data were only included for V magnitudes brighter than 18 5. ni  Chapter 4. COLOUR-MAGNITUDE AND COLOUR-COLOUR DIAGRAMS  69  FIGURE 4.16. Location of blue stragglers on the TEK program field  The circles in the diagram are centred on blue stragglers. The image is a V image of the TEK program field.  Chapter 4. COLOUR-MAGNITUDE  1  -  AND COLOUR-COLOUR  70  DIAGRAMS  i  1  1  r  J  I  I  L  NGC 3201  0.8  o • rH  o  u  0.6  «rH  > • I—I  H-J  0.4 3  0.2  0 _i  I  I  L  0.5 log R  0 (arcmin)  FIGURE 4.17. Cumulative distributions of blue stragglers and subgiants in the TEK program field The ordinate is the cumulative fraction, while the abscissa is the logarithm of the radial distance in arcmins, the cluster centre being defined according to Lee (1977), i.e. the centre of the concentric circles on Figure 5.2. The solid line represents the blue stragglers and the dashed line is for the selected subgiants.  Chapter 5 LEE'S UNUSUAL STAR 5.1 Observations  Lee (1977) reports an unusual star, with a large UV excess, located towards the centre of NGC 3201. He measures the star as having V=15 "89, (B - V)=0 ™78 and (U - B)=-Q ™90. Lee's Unusual Star (LUS) was located on the TEKfield,and is shown in U, B and V band images in Figure 5.1. Close examination of the U band image showed LUS to be blended with a neighbour, which can be visually resolved on the B and V frames. Analysis of the U image showed that LUS corresponded to the lower of the two stars; the centre of the combined images corresponding to the centre of the lower star. My photometry on LUS yields the values V=17 ™08, (B - V)=-0 ?07 and (U - B)=-0 ™94. Figure 4.1 shows the location of LUS in the colour-magnitude plane, it is not seen in Figure 4.2 due to rejection by the PHOTSELECT routine. The errors returned by ALLSTAR were around 0 ™01 for each of the three magnitudes, and LUS was likely rejected due to its larger errors compared with other objects in the same magnitude range.  My photometry is in stark contrast to the values as reported by Lee, except for the (U — B) colour. To see whether systematic differences between the photometry could explain the anomaly, stars common to both data sets were identified (Figure 5.2) and the difference in colour/magnitude were plotted against colour/magnitude (Figure 5.3). Figure 5.3, shows clearly that the discrepancies are not attributable to systematic differences between the photometry. It is interesting to note the relationships between Lee's photometry and that of the present study. The differences may be attributable to Lee having difficulty in obtaining a good sky estimate in the crowded frame of the cluster centre, the effect of the sky becoming more pronounced at fainter magnitudes, as seen in the lower panel of Figure 5.3. At V=17™0 with my system, I expect, by extrapolation, Lee to obtain V ~16 7 with his system. Lee reports having a plate limit of V ~16 3 and therefore I would not have expected him to have detected LUS. From a close examination of Lee's (1977) Figure 9 showing LUS, it is not at all clear which star is bein selected, and I believe that the discrepancy in the photometry may be due to a mismatch by Lee. m  ni  71  72  Chapter 5. LEE'S UNUSUAL STAR  5.2 Nature of the Ultraviolet Excess Object  Assuming the object to lie out of the galactic plane, and applying reddening corrections applicable to NGC 3201, intrinsic colours of (B - V) =-0™28, (U - B) =-l ™09 and (U — V)Q=— 1 ™37 are obtained. Assuming LUS to be a cluster member, and correcting for the adopted distance modulus of NGC 3201, an absolute magnitude of My=2 ™88 is derived. The (75 — V) and (U — B) colours correspond to a spectral class of approximately BO, and hence an effective temperature of around 30,500 °k (Mihalas and Binney 1981). It is interesting to speculate on the nature of LUS. The colours of LUS are consistent with it being a white dwarf. If the object is assumed to be a cluster member, then it is too luminous to be a white dwarf, though the possibility of LUS being afieldwhite dwarf cannot be excluded. Hot subdwarfs have been seen in the globular clusters M4 and M71 by Drukier et al. (1989). Hot subdwarfs, along with the central stars of planetary nebula, appear to be the progenitors of white dwarfs. The origin of hot subdwarfs is still speculative, and I point the interested reader to Heber (1986) for further discussion. The possibility of LUS being a hot subdwarf seems a reasonable one: Greenstein and Sargent (1974) determined a mean absolute magnitude for sdB's of My=+3 ™1 with a dispersion of less than a magnitude, while Heber (1986) determined an average absolute magnitude of My=+4 ™1 ± 0 8 for sdB/sdOB stars. Downes (1986) determined absolute magnitudes of My=(5 ™0/5 ™3) + 1 4 for sdB's, depending on the intrinsic colour adopted. Some other less likely possibilities are listed below. 0  0  Q  0  nl  m  1. ) Quasar: LUS is probably too blue to be a quasar. Richer and Fahlman (1988) state that mean QSO colours are (U - V)=-0 .6 with very few bluer than (U - V)=-0™8. 2. ) Blue drooper: The horizontal branch of NGC 3201 does appear to drop at an uncorrected colour of (B — V)=0 ™2 (Figure 4.1). Whether LUS is a blue drooper is open to debate: the paucity of stars between it and the horizontal branch provides a strong case against the notion. 3. ) Blue straggler: LUS is too blue, and appears quite distinct from the observed blue stragglers r n  To summarize, LUS is probably either afieldwhite dwarf or a hot subdwarf. To make the distinction, a spectrum of LUS is needed.  73  Chapter 5. LEE'S UNUSUAL STAR  I FIGURE  B  I  5.1. Lee's unusual star  The star in question is shown between the inequality signs. The letters below the images indicate the filter, and are U, V and B, clockwise from top left, respectively.  Chapter 5. LEE'S UNUSUAL STAR  74  FIGURE 5.2. Finding chart for the TEK program field The notation in the diagram is similar to that of Lee (1977), with the exception of the pixel coordinates. The dashed line indicates the extent of the TEK program field, while the dotted line indicates the location of the region in Figure 5.1. The numbered stars are used in the comparison of the current photometry with that of Lee's (1977).  75  Chapter 5. LEE'S UNUSUAL STAR  0.3 F  ~i  1  1  r  n  r  J  I  1  L  J  I  ~]  1  r  I  L  0.2 I  0.1 0 0  I  L  J  I  0.5  1.5  (B-V) "i  0.2  r  This work  i  i  r  0.1 I  0  <f  - 0.1  •  •  0.2 J  ,  I  I  L  0  J  I  \  I  L  0.5  1.5  (U—B i  J  I  1  1  1  )  T h i s  work  n  r  1  1  r  0 •••*  •0.2  • x  •0.4 i •  12  13  i  J  i  I  14  L  J  15 V  FIGURE 5.3. Comparison  I  I  I  U  16  This work  with the photometry of Lee (1977)  A comparison between the magnitudes and colours derived from the TEK data for stars in common with the photometry of Lee (1977). The difference, A , is in the sense (Lee value - This work).  Chapter 6 STAR COUNTS IN NGC  3201  6.1 Motivation  So it would seem with galaxies anyway. There is now a growing body of evidence, mainly from galactic rotation curves, that a major part of the mass in the Galaxy is locked up in non-luminous matter, generically known as dark matter. If the mass function of a globular cluster can be traced to the hydrogen burning limit and is found to be still steeply rising, then it can be inferred that a large amount of the cluster mass is locked up in brown dwarfs. It has been shown (Fich and Tremaine 1990) that the dark matter needed to obtain the Galactic rotation curve can be modeled as a spherical halo of radius 35 kpc. It is believed (Larson 1976, Searle and Zinn 1978) that most globular clusters were destroyed at an early epoch of the Galaxy, their stars now constituting part of the halo. It is tempting to hypothesize that the halo mass function will be similar to cluster mass functions, both the halo and globular clusters belonging to Population II, and also that clusters contain brown dwarfs. Brown dwarfs could then constitute the dark matter needed to account for the galactic rotation curve, and hence aid in demystifying the dark matter enigma. 'The  whole is more than the sum  of its (luminous) parts'.  It has been shown (Richer and Fahlman 1991) that at large core radii in a globular cluster, the mass function of the less massive stars, M < OAMQ, will have a slope resembling that of the initial mass function (IMF), provided the cluster has not been extensively tidally stripped and that King models are reasonable approximations to globular cluster kinematics, as appears to be the case. The IMF, in the models of Larson (1986), depends critically on the temperature-density relationship of the gas out of which the stars formed. Knowledge of cluster IMFs may thus help place constraints on the physical properties of the gas in the early universe. If a relationship between cluster metallicity and mass function is proved to exist (McClure et al. 1986), then this would demonstrate that it is easier to form low mass stars in metal-poor environments, though contemporary thinking errs away from such a relationship.  The luminosity function, which is afirststep in deriving the mass function, helps provide a check on the theory of low mass stars. The interior physics of low mass stars is very sensitive to the atmosphere, but the input physics of cool atmospheres is very uncertain. It is believed that a steepening of the mass-luminosity relationship will occur for stars with M < 0.5^0. This 76  Chapter 6. STAR COUNTS IN NGC 3201  77  will be reflected in the luminosity function of the cluster, and hence the luminosity function will provide a check on the theory of cool stellar atmospheres. 6.2 Production of Frames for ADDSTAR Tests  The V and I programfieldsfrom the TI data were binned 2x2, increasing the signal to noise ratio in the images. When binning one must ensure that the stellar profiles are at least still critically sampled, otherwise problems will be encountered with profile-fitting. Standard reductions were performed upon the two binned frames, and IMEXAMINE was used to determine the standard deviation of the sky. Thefindingthreshold was initially set to approximately 3.5 times the sky sigma, but was dropped to a lower value for two of thefieldswhen it was found that fainter stars were not being recovered. Frames of the backgroundfield,located at ai 88=10 38 08 and <S 8=-46 22' 23", h  m  s  0  9  198  were taken in V and / (see observations nee320 —> nee326). Observation nee32\ was binned 2 x 2 to produce the V background frame, whilst observations nee323 —> nee326 were registere combined and binned 2 x 2 to produce the / background frame. The parameters for these frames are summarized below. TABLE 6.1: Characteristics  of frames used for ADDSTAR tests  RON (c  V Program I Program V Background I Background  pixel  3.46 2.69 6 3  Gain - 1 1  (c  A D U  25.2 42 8.4 33.6  Sky o -  1  )  (ADU)  5 4.3 10.5 7  Threshold (ADU)  12 14 35 20  FWHM (arcbccs)  1.6 0.9 1.3 1.1  6.3 Synthetic Star Tests  The efficiency with which stars are recovered from crowded data frames is not known a priori, it will depend, inter alia, on the crowding in the frame, the PSF profile (and hence seeing and guiding) and routines employed in the reduction. These factors affect the efficiency in ways that are hard to determine analytically. To overcome this lack of knowledge, a Monte Carlo technique is employed, where synthetic stars are added into a data frame, and then the data frame is reduced to determine the efficiency with which stars are found. The DAOPHOT routine, ADDSTAR, will generate synthetic stars, with the appropriate Poisson shot noise, and  Chapter  6. STAR COUNTS  IN NGC 3201  78  add them to a data frame. The luminosity function (LF) of the stars used for the ADDSTAR tests should be as similar as possible to the true LF, so that the correction factors calculated will be applicable to the real stars in the frame. Consider, for example, what will happen if real stars are biased to one side of a bin, and a correction is calculated with an even distribution of synthetic stars across the bin, with thefindingefficiency rapidly changing, viz.: the correction applied will be erroneous. Ideally, one would guess the form of the LF, add n per cent of this LF into theframe,evaluate the corrections and deduce a corrected LF. One would then iterate on this method, adding in n per cent of the latest corrected LF until convergence of the LF was found. This method was deemed inappropriate, due to the exorbitant amount of computing time required. In lieu of iteration, a "Guessed Luminosity Function" (GLF) was produced. The GLF was created by performing a standard reduction on theframeunder question, and producing its uncorrected LF. For the brighter stars on theframe,it is a good assumption that they are all recovered. The turnovers seen in uncorrected LFs (see, for example, Figure 6.3) were assumed, and later found, to be due to the inefficiency in the recovery of the fainter stars. To generate the GLF I extrapolate the LF of the bright end to the faint end assuming, a priori, that the turnover is due to incompleteness in the recovery of the fainter stars. Having the GLF then allows a file to be generated containing stars at art n per cent level of the GLF, with a random spatial distribution on the frame. ADDSTAR is then used to add the stars into a frame, and the new frame has a standard reduction performed upon it to determine the efficiency with which the synthetic stars are recovered.  When a synthetic star is recovered, it is possible to make a comparison of its recovered magnitude with its added magnitude. This was done for both the V and I binned program fields, and is shown in Figure 6.1. One notices the scatter to be fairly symmetrical at brighter magnitudes, but at fainter magnitudes there is a tendency for the recovered magnitude to be larger, i.e. fainter than the added magnitude. This is contrary to what is expected, one would expect the stars to be recovered at brighter magnitudes, due to them being on noise spikes. No explanation is offered, though it does seem to be the case for stars with I > 23 ™0. A star was only counted as a recovery if the difference between its added and recovered magnitudes was smaller than 0 ™7, in accordance with Fahlman et al. (1989). 6.4 Robustness of the Technique  To determine the robustness of the Monte Carlo technique, it was decided to see what effect adding stars at varying levels would have on the corrected LF. A GLF was generated for the  Chapter 6. STAR COUNTS IN NGC 3201  79  binned / frame, and ADDSTAR tests were performed at the 5, 25, 75 and 100 per cent leve adding. Figure 6.2 shows the I frame with these various levels of stars added, while Figure shows the derived uncorrected luminosity functions for these frames. Fewer of the fainter star are found on the more crowded frame, as can be seen from the middle panel of Figure 6.6. that in the middle panel there is a point for which more stars are found on the crowded fram attribute this to counting noise: at the faint magnitude associated with the point, only about six stars are recovered. Figure 6.4 shows the recovery efficiency for stars on the frames. Compari the 5 per cent and the 100 per cent add levels, see Figure 6.6, it is seen that in the less c frame stars are recovered with a higher efficiency, as expected. Figure 6.5 shows the correc luminosity functions derived using the frames. The 75 per cent and 100 per cent diagrams sho a sharper upturn at faint, I > 18 ^0, magnitudes where the incompleteness is high, viz: The corrected luminosity functions for the crowded frames contain more faint stars than the corrected luminosity functions for the less crowded frames, as shown in the lower panel of Figure 6.6  I propose that this is due to noise spikes being detected; the probability of detecting a noise spike is not that offindinga star, the stellar image being a distributed function. As the highe add levels are more incomplete (Figure 6.4) a larger incompleteness correction is applied to t same noise spikes, hence causing the sharper upturn. Comparing the 5 per cent and the 100 cent corrected luminosity function graphs, it is seen that the less crowded frame goes deeper than the more crowded frame. The derived LFs for these two frames are in good agreement to where the completeness drops below 50 per cent.  From these tests it was concluded that adding stars at a 10 per cent level will not adverse affect the derived LF, the differences between 5 per cent and 100 per cent only becoming manif at faint magnitudes where the incompleteness is high. All subsequent results were derived by adding synthetic stars to the frames at the 10 per cent level. 6.5 Comparison with the Bahcall and Soneira Galaxy Model  The Bahcall and Soneira (1980, 1984) Galaxy model consists of an exponential disk, representing Population I stars, and a de Vaucouleurs spheroid for modeling the Population II stars. Symmetry is assumed about a north-south direction and also the direction of the galactic centre. Ratnatunga and Bahcall (1985) calculated the expected number offieldstars in the direction of 76 globular clusters with \b\ > 10° using the above model. NGC 3201, with galactic latitude of 8 ? 64 is not included, and the authors give the below warning.  Chapter 6. STAR COUNTS IN NGC 3201  80  "The estimates given for the 32 globular clusters with 20° > b > 10° and for all fields with V > 22 must be used with special caution because these regions are outside the galactic latitud or apparent magnitude range over which the galaxy model has been tested."  With these facts noted, a comparison is made. The corrected counts are shown in Table 6.2 The counts for the clusters NGC 2808 and NGC 6101 were multiplied by 4.6, so as to have the same sky coverage as the TI frame. Figure 6.7 shows stars which were matched between the TI backgroundfieldI and V frames, along with the magnitude bins as used by Ratnatunga and Bahcall, and the incompleteness in V each bin. The completeness corrections were calculated using PSFs, and so should be applied to stars only, not noise. To try and ensure this was happening, the stars were counted in two ways. The column headed unmatched is the corrected number of stars counted in each V bin, and is susceptible to noise contamination at the faint end. The matched column is the number of stars in each V bin which were also found on the I frame. By setting this criterion, a lot of noise is rejected and the counts should be more reliable. The I data are fairly complete to around 22 ™0; the last V bin may have too few stars due to incompleteness in I. The stars in the last V bin span a range of 20 ™5 < I < 25 ™0, and it is not obvious how to correct for this incompleteness. The matched and unmatched counts show good agreement in all but the last two bins. The discrepancy is vindicative of noise contamination in the last two bins. TABLE 6.2: Comparison of counts with the Bahcall and Soneira model  Range in V 13 < V< 15 15 < V< 17 17 < V < 19 19 < V< 21 21 < V < 23 23 < V < 25  NGC 3201, Background Field  NGC 2808  7=280 ?34, 6=10 °54  unmatched  matched  0 4 19 26 79 1806  0 4 19 26 53 230  7=282°, b=-U  0.92 3.588 11.5 22.08 33.58 59.8  ".3  NGC 6101 7=317°, £ = - 1 5 °  1.012 3.864 11.96 23.26 37.26 64.4  Comparing the counts for NGC 2808 and NGC 6101, one sees that they are similar. When symmetry considerations are taken into account, the Galactic coordinates of the TI background field differ from those of NGC 2808 less than NGC 2808 differs from NGC 6101. With this observation the assumption is made that the NGC 3201 counts should be similar to those of NGC 2808. This is approximately seen for 13 0 < V < 21 0, but divergence is seen thereafter. One m  m  Chapter 6. STAR COUNTS IN NGC 3201  81  has to be careful about drawing conclusions as the model was never meant for the low galactic latitudes of NGC 3201, though I tentatively suggest that the Bahcall and Soneira model predicts too few faint stars at low galactic latitudes. Fahlman, Richer and VandenBerg (1985) found a similar effect for M15, with /=25° and b—21°, where too few faint stars were predicted by the Bahcall-Soneira standard Galaxy model. The divergence from the models is a clear warning not to use the model to correct for faint background stars.  6.6 Mass-Luminosity Relationships A mass-luminosity relationship (MLR) is a relationship between magnitude, in this case V or I, and stellar mass. The MLR, like an isochrone, is dependent on the stellar model, i.e. the luminosity is not only a function of mass. The MLR will change as the cluster ages, and will also be a function of metallicity, mixing length parameter, helium content and other less significant factors. The MLR I adopt for the Johnson V band in this work is that of Drukier et al. (1988), which was obtained using a 16 Gyr VB85 isochrone, and points which extend the relationship down to M = 0.15.Mo (Figure 6.8). The MLR is applicable to stars with a metallicity of [Fe/H]=-1.27.  The I band MLR was derived by using the MLR of Fahlman et al. (1989) and the models of VB85. Points were read off the Fahlman et al. MLR, and splined to points from the VB85 16 Gyr isochrone. It was found that to ensure a smooth transition between the two sets of points, the brighter end of the Fahlman et al. relationship had to be truncated, see Figure 6.9. The Fahlman et al. relationship was derived from LFT (Luyten Five Tenths Catalogue, Luyten 1955) subdwarfs, observed by Rodgers and Eggen (1974), and the models of D'Antona (1987). I refer the interested reader to Fahlman et al. for more details. The I band MLR adopted was applicable to Population II stars of intermediate metallicity, and an age of 16 Gyr. Figure 6.10 shows the effect of having alternate ages for the MLR, and it is seen that changing the age of the cluster will not have a large effect for stars less massive than 0.7MQ. The MLRs are applicable only to main sequence stars. 6.7 Error Propagation in ADDSTAR Tests  To derive the mass function in the programfieldit is necessary to count the stars in magnitude bins. The stars actually counted on a frame have to be corrected for incompleteness. Denoting No as the number of cluster members in a particular bin, we have No = Np — N R with Np and N B being thefieldand background counts in the same magnitude bin, after being corrected for  82  Chapter 6. STAR COUNTS IN NGC 3201  incompleteness. Denoting IC as the incompleteness correction and RC as the raw counts, N  R  = IC x R C F  F  (6.1) N  B  = IC x RC . B  B  Propagation of errors leads to the expression  for either thefieldor background frames. The error in the raw counts is given by the Poissoni counting error,  The statistics of ADDSTAR tests are described well by the binomial distribution (Bolte 1987). The standard deviation of the binomial distribution is given by VV(I-P),  (6.3)  where n is the number of trials, equivalent to number of stars added, and p is the probabilit of success, equivalent to the number of stars recovered divided by the number of stars added. Defining N A as the number added, and NR, as the number recovered then yields -(N .) = R  N R  A  / ( ^ - i - )  .  (6.4)  Defining IC as  ( ^-Y=(^)\(^) a  2  .  (6.5)  The number of stars added to the frame is a known quantity and so a A will be zero. Havin the error in the incompleteness corrections will then allow the errors in the background and fi counts to be determined. These errors are added in quadrature to determine the error in NoThe error in is evaluated by dividing the error in the true counts by the width of the mass bin. The MLR is assumed to have no error, and so dM is also assumed to be error fr  Chapter  6. STAR COUNTS  IN NGC 3201  83  6.8 The Mass Function The derived LFs from the program fields consist of two components; a cluster component and afieldcomponent. Thefieldcomponent consists of background and foreground objects not associated with the cluster. To determine the true cluster LF, it is necessary to produce a LF for the backgroundfieldsand then, assuming the background distribution to be homogeneous and hence the background contribution to the program frame to be similar to the background frame, subtract the background LF from the programfieldLF. This procedure should automatically make corrections for galaxy contamination; the samefindingroutines are employed on both frames. Taking E(B - V)=0 ™21 and (m - M ) V = U ™2 for NGC 3201, the absorption in the V band will be  A  V  =R  V  xE(B - V) = 0.65,  (6.6)  where R is taken to be 3.1, in accordance with Savage and Mathis (1979). This yields a true distance modulus of 13 ™55 for NGC 3201. Cohen et al. (1981) give that A / = 0.6 x A y ,  the application of which then gives (m — M)j—13  nl  (6.7)  94.  Applying the colourless transformation equations (§3.4.3) to ALLSTAR magnitudes from the V and I program fields and correcting the values for the appropriate distance modulus yields absolute magnitudes for the cluster members. The V and J background frames were reduced using standard reductions, and GLFs were built. The frames had ADDSTAR tests, at a 10 per cent level of the GLF, performed upon them. The ALLSTAR magnitudes of objects detected on the background frame were calibrated to standard magnitudes (§3.3), and corrected for the appropriate distance modulus, V or I, of the cluster. This ensured that the background stars appeared in appropriate magnitude bins to allow the calibrated program field counts to be corrected. Cluster counts were derived by applying incompleteness corrections to the magnitude bins, and then subtracting the background counts from the field counts. At bright magnitudes, where there is a paucity of stars, it is possible for the cluster count to be negative due to the effects of small number statistics.  Chapter  6. STAR COUNTS  IN NGC 3201  84  The quantity <f> is defined as  with N representing number of stars and M representing mass. The quantity <f> is a conti function, though I only derive the mean value in a bin. To obtain <f> for a bin, the mas bin ends is evaluated, and the change in mass is derived. The number of stars in the bin is t divided by the change in mass to yield a value for <f>. Mass functions were derived from the V and I frames in two ways. Thefirstway was to count the stars in magnitude bins, the bins having a width of 0 5 and an interval of 0™25 so that the bins overlap. A bin width of 0 5 was chosen as a compromise between resolution and Poisson-counting noise. This produced the mass functions shown in Figures 6.11 and 6.12. Figure 6.13 shows the two mass functions superimposed. These mass functions have been corrected for background counts, and points were only plotted if the uncorrected number of cluster stars was larger than ten. Comparing the mass functions derived using the V and / frames, one notices considerable differences between them, for example the two minima seen in the 'I mass function', located at approximately —0.2 and —0.4 in logio^^^^, which do not appear in the 'V mass function'. One would intuitively not expect there to be a difference since the same stars are being counted. A possible explanation of this anomaly may be that the star colours are inconsistent with the MLRs. To check this, a synthetic isochrone was generated by evaluating V and I magnitudes at equal masses on the two MLRs. The isochrone was plotted in the (V — I), I plane, along with the calibrated data from a standard reduction on a binned V frame and a binned / frame, see Figure 6.14. One sees that the synthetic isochrone is a goodfitto the data. The fact that t isochronefitsthe data to its limiting magnitude lends credence to the MLRs of Fahlman et al. (1989) and Drukier et al. (1988), especially as the extensions to the VB85 isochrones, for eac of the MLRs, were derived in alternate ways. m  nl  To make a more direct comparison between the two mass functions, it was decided to bin on the mass interval, as opposed to the magnitude interval so that, bin jumping aside, exactly the same stars would be counted in the same bins. The magnitude intervals corresponding to mass intervals of 0.05.Mo were calculated using the V and / MLRs. Stars on the program and background frames were counted in these bins, and the applicable corrections were derived by considering the stars added and recovered from  Chapter 6. STAR COUNTS IN NGC  3201  85  each of these bins. Figures 6.15 and 6.16 show the two mass functions derived in this way before any background corrections, and Figures 6.17 and 6.18 show the mass functions after the background corrections. One immediately notices that the two minima seen in the 7 mass function' are still present and still do not appear in the 'V mass function' for both the corrected and uncorrected mass functions. To explain these anomalies between the mass functions, plots were made of I magnitude against V magnitude of stars in both thefieldand background frames (Figures 6.19 and 6.20). The results from a single reduction are expected to be very similar to the average of the results of twenty reductions, the number of ADDSTAR tests used to derive the mass function. These graphs show the magnitude equivalent of equally spaced mass bins, and the completeness corrections applied in each bin. The stars found in boxes running from bottom left to top right are stars which are located in the correct bins in both colours, whereas stars located in the off-diagonal boxes have 'bin jumped'. Bin jumping is inevitable using bins, and will not have an adverse effect on the results provided that the scatter is symmetrical about bin intersections. A calibration error may cause the points to be asymmetric. It is hard to decide whether such an error exists in the present data, but it appears not. Considering the third bin, which shows a discrepancy between the mass functions, one sees that, assuming all the completeness corrections to be unity, N  V :  N  = C + 20  0  d  = B + 3  B  d  (6.9) J  No = C + 14  :  d  N  =  B  B  d  + 14.  In the above notation, N represents the total number of stars in the cluster or background bin, while (C/B)a is the number of stars in the 'on-diagonal' box, indicated by the intersection of the arrows in Figures 6.19 and 6.20. Before correcting for background counts, the third I bin has a deficit of six stars as compared to the third V bin, this is sufficient to explain the discrepancies seen between the mass functions infigures6.15 and 6.16. After correcting for background counts, the number of stars in the third mass bin, derived using the V data will be (C - B ) + d  d  20 - 3,  while for the I data it is (C d  B) d  + 14 -  14.  Chapter  6. STAR COUNTS  IN NGC 3201  86  The anomaly, far from being corrected, is now worse. The difference between the two values is sufficient to explain the lower value seen in the third bin of the background corrected mass function that was derived using the I data. The anomaly between the two mass functions is caused by asymmetrical scattering and stars with 'wrong' colours, for which the MLRs are not applicable. If the stars on the 'off-diagonal' bins were not to be counted, on the premise that their colours are inconsistent with the isochrone, and therefore the MLR is not applicable to them, then stars with genuine scatter would be lost, especially at bin intersections and at fainter magnitudes. If one is to count the off diagonal stars then the validity of the MLR to these stars must be questioned. For this thesis, it was decided that off-diagonal stars would be counted. 6.9 Implications of the Mass Function  The mass function of NGC 3201 is found to have a similar morphology to previously derived mass functions (e.g. see Richer et al. 1991, henceforth RFBFST). In common with RFBFST, an upturn is found in the mass function for stars with M < 0.4.M©. Fitting the three points corresponding to masses less than 0AM®, for which the completeness is greater than 50 per cent, yields a value forjc of approximately 1.9. If the last point shown in Figure 6.18 were to be taken into consideration, the mass function would be found to be more steeply rising, but the validity of this last point must be questioned due to its large incompleteness. No evidence whatsoever is seen for a flattening of the mass function at lower masses. Assuming the last point to be valid, I can state that no downturn is seen in the mass function down to a mass of approximately 0.2MQ. An extension of the mass function into the brown dwarf region is then suggested by continuity arguments; there being no known reason for a cut off in the star formation at the hydrogen burning limit. 1  Table 6.3 is from RFBFST, with the addition of NGC 3201. The value of x is for M < 0.4.M©, while the third column is the half-mass relaxation time, T1/2, from Binney and Tremaine (1987), and thefifthand sixth columns are the concentration parameter, c, and metal abundance, [M/H], taken from Webbink (1985). The entry in the fourth column, To, is what RFBFST term as the "time until destruction" and is the inverse of the "destruction rate" for the cluster, as calculated by Aguilar, Hut and Ostriker (1988). In the Population II regime, the hydrogen burning limit is believed to lie in the range (D' Antona 1987).  Q.WMQ  <  Muiu  < O.II-MQ  Chapter 6. STAR COUNTS IN NGC 3201  87  TABLE 6.3: Properties of clusters with derived mass functions extending well below OAMQ  Cluster  X  2.7 2.5 2.1  M13 NGC 6752 (jj Cen  NGC  3201  M5 NGC 6397 M 71  I  B  1.3 0.9 0.5  T 9.40 8.94 9.96 9.14 9.44 8.47 8.35  log  1/2  T 7.1 7.1 6.7 6.25 4.5 2.5 2.9 D  c 1.4 1.6 1.2 1.55 1.8 1.6 1.1  [M/H] -1.6 -1.6 -1.6 -1.6 -1.6 -2.0 -0.5  RFBFST conclude from the six clusters in their study that there is no correlation between th mass function slope and the concentration parameter or the metallicity. This conclusion weighs heavily against the McClure (1986) effect where a correlation between x and [Fe/H] is proposed RFBFST suggest that there is a weak correlation between x and the half-mass relaxation time, and a strong correlation between x and the destruction time. It is seen that the results for NG 3201 provide no contradiction between these correlations. Figure 6.21, adapted from RFBFST, shows plots of the destruction time and the half-mass relaxation times against x. A possible interpretation of the results is that all globular clusters start with mass function slopes of approximately three, which are modified due to dynamical processes (RFBFST). Mass segregation is faster in clusters with short relaxation-times, resulting in an enhancement of lowmass stars in the outer parts of the cluster. The lower-mass stars are preferentially stripped wh the cluster transits the Galactic plane. This can explain the positive correlation seen between half-mass relaxation time and mass-function slope; clusters with larger half-mass relaxation times have lost fewer of their low-mass stars. A stronger correlation between 'time until destruction' and x would then be anticipated, as T D takes account of all the processes thought to be importan in modifying the MF of the cluster.  The derived MF slopes of globular clusters can be taken as lower limits on the slope of t IMF. If the halo mass-function is similar to that of the globular clusters, then it will also have steep slope at low masses. Extrapolation of the halo mass-function into the brown dwarf regio suggests a large population of brown dwarfs in the halo, and hence a large mass-to-light ratio A large mass-to-light ratio in the halo would have a significant bearing on Galactic structure.  Chapter 6? STAR COUNTS IN NGC 3201  "i  88  r  n  r  L  j  i  0.5 0-  T3  d  -0.5 J  I  i_  20  18  24  22 add  ~i  1  i  r  1  r  0.5 I  0 -0.5 J  18  20  22  l_  24  add  FIGURE frames  6.1. Error diagram for synthetic stars added to the binned TI program field  The subscripts 'add' and 'rec' refer to the magnitudes at which the synthetic star was added and recovered. Stars below the zero line on the ordinate were recovered at a fainter magnitude at which they were added, and vice versa.  Chapter 6. STAR COUNTS IN NGC  3201  89  FIGURE 6.2. Varying crowding conditions on the I program field  From top left, going clockwise, the percentage add levels are 0, 25, 75 and 100 per cent.  Chapter  6. STAR COUNTS  n  IN NGC 3201  90  | i i i | i i i | i i i | i  i  r  25% Add level  5% Add level ii  i  iii  i  ****  1  1 t „l  ii  1.5  i—I—i—i—i—i—i—i—r  a  Ml  o  i  -  1  cd  w 0.5  M i l l  I I I  CD  100% Add level  r£3  i  5 ?J*J  i i  I  I  I  75% Add level  TT  i  3 iz;  «  tiO  o  ii 0.5  J  i  S Ti  1.5 -  I i i i I i i i L _ i i i I i i il I i I i i i I i i i  14  16  18  20  14  16  I  :  i i  18  20  KC  FIGURE 6.3. Uncorrected luminosity functions for various crowding conditions The above plots were generated by adding stars to the TI binned / program field at varying levels. Differences in the morphologies only become manifest after the turn over, the fainter stars being lost on the more crowded frames. (N.b. only real stars are counted.)  91  Chapter 6. STAR COUNTS IN NGC 3201  1  -i  1  1  1  1  1  1  1  1  1  1  1  1  r  5% Add level  I  I  I  T"  ~|  i  i  \  |  i  i  i  25% Add level  Si* ! 1  0.8  0.6  0.4 id  0> rH  0.2  CD >  o o  CD  0  rH  1  H—I—I—I—I—I—I—I—I—I—I—h  o o cd  100% Add l e v e l -  ********  • rH H->  H—I——I—I—I——I—h 75% Add level  0.8  0.6  0.4  0.2  '  i  I  16  I  I  I  I  18  I  I  L.  20  1_  i  22  i  i  i  i  16  i  18  i  •••••(  i J  20  i  i  i _  22  -^KC  FIGURE 6.4. Recovery efficiency for various crowding conditions The above plots were generated by adding stars to the TI binned / program frame at varying levels. At 100 per cent, more stars are 'lost' and hence the luminosity function cannot be determined to such faint magnitudes, this effect is attributable to crowding.  Chapter 6. STAR COUNTS IN NGC 3201  ~i  r  r • i 1  i  I  i  i  i  I  92  ~p  i  "T  5% Add level  i  I  i  i  i  i  i  r  i  I  i  r  25% Add level 5$$  o  .II*  = 11'  w rH  CO  w «fH  o rH  CD  a  r Q  I I I  Tj  M i l l  I I I  I  100% Add level  CD H->  I I I  I I I  I I I !  75% Add level  O CD rH rH  o 00  o  ***  ii  ii  r J  ii  1  .1* i  14  i  i  i  i  16  1_  I  I  I  18  I  I  I  I  20  I  I  L  j  i  14  i  i  i  i  16  i  i  i  i  i  i  18  i  I  i  i  20  FIGURE 6.5. Corrected luminosity functions for various crowding conditions The above plots were generated by adding stars to the TI binned / program frame at varying levels. The open circles indicate bins in which the completeness was less than 50 per cent. The difference in the morphologies is discussed in the text.  i  Chapter 6. STAR COUNTS IN NGC 3201  ~\  1  O CD rH  r1  _•_ _ — =-  1  • • •  cd rH  1  93  0.8  1  —•  ~t  1  —  1  r  H  h  r  0.6 0.4  H->  cd  PH  0.2 0  P>H  0  X*  -0.2  -I  h  H  1  h  =1  h  H  1  h  H  h  -I  h  hJ  CD  H->  O CD  -0.4  SH SH  o o  -0.6  ro  -0.8  PI  -1 En rj X) CD H-> O CD  H  1  h  _1  I  L.  H  h  J  I  -F  0.6 0.4 0.2  rH rH  o o  0 -0.2 14  16  18  1_  20  ^KC  FIGURE  6.6. Comparison  of frames with 5% and 100% add  levels  The ratio of the recovery fractions (in the sense 100% 5%) are shown in the upper panel. The lower panels show differences in the uncorrected and corrected luminosity functions (see Figures 6.3 and 6.5). The difference, A, is in the sense (100% -5%) -=- 5%. Open points in the bottom panel indicate that the recovery fraction was less than 50 per cent, and dashed lines are to guide the eye.  Chapter 6. STAR COUNTS IN NGC  ~\  3201  1  94  1  19  r  1  ~i  26  1  46  r  30  25  • »•  • •• •  * »> •  20  15 0.87  0.98 _J  15  I  0.13 I  L  20  25  V  FIGURE 6.7. Stars found in two colours, with bins and incompleteness  shown  The stars shown were found on both the TI Vand / background fields. The incompleteness of each Vbin is indicated in the lower part of the panel and the uncorrected number of stars in each bin in the upper part of the panel.  FIGURE 6.8. V mass-luminosity relationship The relationship is from Drukier et al. (1988). Circles denote values from VB85 isochrones, while squares indicate the Drukier et al. extension of the relationship to fainter magnitudes.  Chapter 6. STAR COUNTS IN NGC 3201  FIGURE  96  6.9. I mass-luminosity relationship  Circles denote points taken from VB85 isochrones, while squares denote points from the Fahlman et al. (1989) MLR, which was adopted for / > 6 P5.  Chapter  6. STAR COUNTS  ~i  i  IN NGC  i  I  i  i  3201  i  97  I  i  i  i  I  I  I  1  I  1 1 1 1 1 1 1 1 1 1  I  I  r  A  6  ^  8 -  M = 0.7Ji . r  10  12  V 14  0  -0.2  J  1  •0.4  -0.6  I  I  I  -0.8  I  I  I  I  I  I  I  1.2  log (JL/Jl ) 10  o  FIGURE 6.10. The effect of cluster age on the V MLR The three curves are VB85 model isochrones with «=1.6, [Fe/H]=-1.27, and Y=0.20 and ages of 14, 16 and 18 Gyr from left to right respectively. The point on the abscissa corresponding to a mass of 0.7M is indicated. Q  Chapter 6. STAR COUNTS IN NGC 3201  ~i  NGC  •%  98  1  n  r  I  L  r  3201  3.5  O  J  0  i  i  i  i  0.2  0.4 log (Jl/Jl ) 10  o  J  J  -0.6  L  -0.8  (  FIGURE 6.11. Mass function for the V TI program field  The stars were binned into bins of width 0 ™5 and interval 0 ^25. The mass function has been corrected for background objects. Open circles indicate bins in which the completeness was less than 0.5.  Chapter 6. STAR COUNTS IN NGC 3201  1  1  1  |  NGC  99  ,  l  1  i  i  i  i  i  i  i  i  3201  >  I  —  if  nj  nn  fi  1  u II  i I  -  <  i  0  i  1  ,  ,  i , -0.4  ,  -0.2  ,  ,  i -0.6  i  i  i  -0.8  log (Ji/Jl ) 1 0  e  FIGURE 6.12. Mass function for the I TI program field The stars were binned into bins of width 0 ™5 and interval 0 25. The mass function has been corrected for background objects. Open points indicate bins in which the completeness was less than 0.5. ni  Chapter 6. STAR COUNTS IN NGC  I  0  i  ik_i  i  100  3201  I  i  i  i  I i  i  i  -0.4  -0.2  i  i  -0.6  i  i  i  i  -0.8  log (Jt/it ) 1 0  FIGURE 6.13. Mass functions forV  o  and I TI program fields  Squares are for the 'V mass function' and Circles are for the 7 mass function'. The stars were binned into bins of width 0 5 and interval 0 ™25. The mass functions were corrected for background objects. Open points indicate bins in which the completeness was less than 0.5. ni  101  Chapter 6. STAR COUNTS IN NGC 3201  (V-I) FIGURE 6.14. An isochrone derived from the V and I MLRs The points shown were derived by matching stars found on a standard reduction of the binned Vand /fields. The stars were calibrated to absolute magnitudes by use of colourless transformations, equations 3.18, and the errors are those returned by ALLSTAR, added in quadrature for (V-I). The isochrone was reddened by E(V — I)=0 286 and corrected for the / distance modulus of NGC 3201. ni  102  Chapter 6. STAR COUNTS IN NGC 3201  i  i  i  i  i  i  i  1  1  i  i  i  i  3201  NGC  4 >  <>  -  3.5  —  I I  ;  ti  [11  ]  3  i  0  i  i  i  -0.2  i  i  i  -0.4  log (JL/Ji ) 10  i  i  -0.6  i  i  -0.8  o  FIGURE 6.15. Uncorrected mass function for the V TI program field The mass function was not corrected for background objects. The bins used are of equal width in mass. Open circles indicate bins in which the completeness was less than 0.5.  103  Chapter 6. STAR COUNTS IN NGC 3201  ~i  NGC  3.5  ~i  r  1  1  1  1  1  r  I  L  3201  -  o  J  0  i  i_  J  J  L.  -0.2  -0.4  J  -0.6  l_  -0.8  log (JA/Jlo) 10  FIGURE 6.16. Uncorrected mass function for the I TI program field  The mass function was not corrected for background objects. The bins used are of equal width in mass. Open circles indicate bins in which the completeness was less than 0.5.  104  Chapter 6. STAR COUNTS IN NGC 3201  T  i  i  NGC  •%  |  i  i  —|  i  i  i  i  J  I  L  |  i  i  r  3201  3.5  O  0  -0.2  -0.4  J  -0.6  L  -0.8  \og (M/M ) 10  Q  FIGURE 6.17. Corrected mass function for the V TI program field  The mass function was corrected for background objects. The bins used are of equal width in mass. Open circles indicate bins in which the completeness was less than 0.5.  Chapter 6. STAR COUNTS IN NGC 3201  105  FIGURE 6.18. Corrected mass function for the I TI program field. The mass function was corrected for background objects. Bins used are of equal width in mass. Open circles indicate bins in which the completeness was less than 0.5. The solid lines indicate relationships of the form log\o<f> = login  with x having the values indicated.  M  -(l+x)-  Chapter 6. STAR COUNTS  1  10  Bin  8  6  1  -  i  IN NGC 3201  |  i  >  i  106  |  i  1  1  1  1  1  1  1  1  1  1  1  1  12  0.21  -  11 10 9 8 7 6  0.52 0.66 0.79 0.86 0.87 0.9  -  5  0.89  -  4  0.94  —  3  0.96  2  0.94  1  0.96  -  3«-  -  .V " 1  it  '.  ./'  1\  i  0  1  2  ' i  i  -i  1 i  4  i  0.94  o cn c6  cn d  0.79  0.69  0.48 0.28 0.13  —  -  CV2  CO  •*  in>  co  co cn  1  1  6  1  1  1  8  t  TOO  0.96  -  0  1  Completeness  4  2  1  o  o  —  o 1  10  1  1  1  12  V  FIGURE 6.19. A 'magnitude-magnitude  diagram' for the TI program field  The diagram was generated by matching stars on the binned V and / program fields. The magnitude equivalents of equal mass divisions (in 0.05M ) are shown. The bins are numbered, and the recovery efficiencies in the bins are indicated. Q  Chapter 6. STAR COUNTS IN NGC 3201  1  10  1  Bin  8  -  1  |  i  1  1  |  107  i  1  1  1  12 11 10 9 8 7 6 5 4 3 2  0.53 0.93 0.9 0.98 0.93 0.93 0.98 0.98 0.97 0.93 1  1  1  |  1  1  !  1  1  — 3»-  A  —  CO 05 o  ^H  C\2 i  0  1  j  _  -  0  1  Completeness  -  2  1  1  2  i  i  i  !  4  i  \  CO  1  1  6  1  CO 05 O  CM OS O  in  CD  1  8  CO CO in CO CO . - H q o o o d  in  CD 05 I  I  o  o  o  -rH  1  —  o  1  10  C\2  1  1  1  12  V  FIGURE 6.20. A 'magnitude-magnitude' diagram for the TI background field The diagram was generated by matching stars on the binned Vand / background fields. The magnitude equivalents of equal mass divisions (in 0.05.M©) are shown. The bins are numbered, and the recovery efficiencies indicated.  Chapter 6. STAR COUNTS IN NGC 3201  108  i i i i i i i i i i i u i i i i i i i ii  i i i i I i i i i Ii I i i Ii i i i I i  10  7  6  9.5 CD  w  6  rH  CD \  o&o 9  r—I  X a  8.5 i i i i i i i i i i i i i i i i i i i i i i  0.5 1  1.5 x  2 2.5  i i i i i i i i i i i i i i i i i i i i i  0.5  1 1.5 2 2.5 x  FIGURE 6.21. Plots of half-mass relaxation time and 'time until destruction' against mass function slope below 0.4M Q  Circles represent data from RFBFST, while the square represents data for NGC 3201, the value of x being determined in this study.  Chapter 7 CONCLUSIONS Both deep and shallow CMDs were produced for NGC 3201. The shallow data allowed for estimates of the cluster reddening and metallicity. The values obtained were not as tightly constrained as previous studies, but did show general agreement. Matching the deep CMD to a subdwarf fiducial allowed for a value of the distance modulus to be derived. A value of (m — M)^=14 *P-2 was found, in agreement with other values in the literature. Synthetic isochronefittingto the deep CMDs implied an age of 15 + 2 Gyrs on the scale of VB85. Isochronefittingto the shallow data set produced consistent results, though in both cases thefittingwas poor, and the derived ages must be taken as being inconclusive. No obvious equal-mass binary sequence was seen in either the TI or TEK field, though it was stressed that this does not preclude the possibility of binaries in NGC 3201. The CMD for the shallower data set had a group of blue stragglers. The blue stragglers were shown to be more centrally concentrated than subgiant stars in the same magnitude range. This provides strong evidence for mass segregation in NGC 3201. An unusual star,firstobserved by Lee (1977), was noted. New values were derived for the U, B and V magnitudes of the star, the new values being contrary to those of Lee, due to a possible mismatch by him. The possible nature of the object was speculated on. Using Monte Carlo techniques, the luminosity function of the cluster in both the V and the I band was derived. The I band MLR from Fahlman et al.. (1989) was adapted slightly, and then used to generate the mass function of NGC 3201. A V band MLR from Drukier et al. (1988) was similarly used to generate a mass function. Differences between the two mass functions are discussed, and a reason proposed for the differences seen. The two MLRs were used to generate a synthetic isochrone, and this isochrone was shown to follow the morphology of the data. This shows that the two MLRs are of the right form, unless errors had conspired to give the apparent match. Finally, the gradient of the mass function was determined for the lower mass stars. The mass function for the less massive stars was found to have a slope of x=\.9. The 'time until destruction' and the 'destruction rate' for the cluster are found, and the data are found to fit into the picture of RFBFST well, implying that many dynamical processes are at work in the destruction of globular clusters. 109  BIBLIOGRAPHY  110  BIBLIOGRAPHY  Aguilar, L., Hut, P., and Ostriker, J.P. 1988, Ap. J. 335, 720 Alcaino, G. 1974, in Atlas of catolica de Chile), p. 26 Alcaino, G. 1976, Astron.  Globular Clusters  (Talleres de la imprenta de la universid  Astrophys. Suppl. 28, 409  Alcaino, G., and Liller, W. 1981, A. J. 86, 1480 Alcaino, G., and Liller, W. 1984, A. J. 89, 1712 Alcaino, G, Liller, W., and Alvarado, F. 1989, Astron.  Astrophys. 216,  Bahcall, J.N., and Soneira, R.M. 1980, Ap.  J. Suppl., 44, 73  Bahcall, J.N., and Soneira, R.M. 1984, Ap.  J. Suppl., 55, 67  Bennett, P.D. 1989, Document)  OPDATA Reference Manual  Bevington, P.R. 1969, Data Reduction York: McGraw-Hill), chapter 4  68 (ALA89)  (Vancouver: TRIUMF Computing  and Error Analysis for the Physical Sciences  (N  Binney, J., and Tremaine, S. 1987, Galactic Dynamics (Princeton: Princeton University Press), p. 514 Bolte, M. 1988, PA.S.P.  100, 869  Bolte, M. 1989, Ap. J. 341, 168 Borner, G. 1988,  The Early Universe  (Berlin Springer-Verlag)  Buonanno, R., Buscerna, G., Corsi, C.E., Ferraro, I., and Iannicola, G. 1983, Astr. 278 Cacciari, C. 1984, A. J. 89, 231 Carney, B.W. 1980, Ap.  J. Suppl. 42, 481  A  BIBLIOGRAPHY  111  Chieffl, A., Straniero, O., and Salaris, M. 1991, in ASP Conference Series Volume Formation and Evolution of Star Clusters, ed. Kenneth James (BookCrafters inc.), p. 219  13,  Cohen, J.G., Frogel, J.A., Personn, S.E., and Elias, J.H. 1981, Ap. J. 249, 481 Da Costa, G.S., Frogel, J.A., and Cohen, J.G. 1981, Ap. J. 248, 612 D'Antona, F. 1987, Ap. J. 320, 653  Demarque, P., and McClure, R. 1977, In The Evolution of Galaxies and Stellar Populatio eds. B. Tinsley and R.B. Larson (Yale University Observatory, New Haven) Downes, R.A. 1986,  Ap. J. Suppl. 61,  569  Drukier, G.A., Fahlman, G.G., Richer, H.B., and VandenBerg, D.A. 1988, A. J. 95, 14 Drukier, G.A., Fahlman, G.G., and Richer, H.B. 1989, ApJ. 342, L27 Fahlman, G.G., Richer, H.B., and VandenBerg, D.A. 1985, Ap. J.  Suppl. 58,  225  Fahlman, G.G., Richer, H.B., Searle, L., and Thompson, LB. 1989, Ap. J. 343, L49 Fich, M., and Tremaine, S. 1990, Preprint Frogel, J.A., Cohen, J.G., and Personn, S.E. 1983, Ap. J. 275, 773 Graham, J.A. 1981,  PA.S.P. 93,  291  Gratton, R.G. 1982, Astron. Astrophys.  115,  171  Gratton, R.G., and Ortolani, S. 1989, Astron. Astrophys. Ill, 41 Greenstein, J.L., and Sargent, A.I. 1974, Ap. J. Suppl. 28, 157 Heber, U. 1986, Astron. Astrophys. 155, 33 Iben, I., and Rood, R.T. 1970, Ap. J. 159, 605 King, I.R. 1966, A. J. 71, 64  BIBLIOGRAPHY  King, I.R. 1971,  112  PA.S.P. 83, 199  Kinman, T.D. 1959,  MNJIA.S.  Landolt, A.U. 1973,  Astron. J. 78,  959  Landolt, A.U. 1983,  Astron. J. 88,  439  119,  157  Lang, K.R. 1980, in Astrophysical Formula p. 518  (second edition)  (New York: Springer-Verlag),  Larson, R.B. 1976, in Galaxies, eds. L. Martinet and M. Mayer (Geneva: Geneva Observatory), p. 67 Larson, R.B. 1986, in  IAU Symposium 126, The Harlow-Shapley Symposium on Globular Cluster Systems in Galaxies, eds. J.E. Grindlay and A.G. Davis Philip (Dordrecht:  Kluwer), p. 311 Lee, S.W.  1977, Astron. Astrophys. Suppl. 28, 409  Leonard, P.J.T. 1989, A. J. 98, 217 Lutz, T.E., Hanson, R.B., and Van Altena, W.F. 1988, in The Extragalactic Distance Scale, eds. S. van den Bergh and C.J. Pritchett (Utah: Brigham Young University Print Services) Luyten, W.J. 1955, A Catalogue of 1849 Stars with Proper Motions Exceeding 0."5 Annually (Minneapolis: Lund Press) MacKay, CD. 1986,  Ann. Rev. Astron. Astrophys. 24, 255  McClure, R.D., Hesser, J.E., Stetson, P.B. and Stryker, L.L. 1985,  PA.S.P. 97, 665  McClure, R.D., VandenBerg, D.A., Smith, G.H., Fahlman, G.G., Richer, H.B., Hesser, J.E., Harris, W.E., Stetson, P.B., and Bell, R.A. 1986, Ap. J. 307, L49 Menzies, J. 1967,  Ph. D. thesis  (Australian National University)  Mihalas, D., and Binney, J. 1981, Freeman)  Galactic Astronomy (second edition)  (San Francisco:  BIBLIOGRAPHY  113.  Nemec, J.M., and Harris, H.C. 1987, Ap. J. 316, 172  Penny, A. J. 1984, in IAU symposium 105, Observational Tests of the Stellar Evolution Theo eds. A. Maeder and A. Renzini (Dordrecht: Reidel), p. 157 Penny, A.J., and Dickens, RJ. 1986,  M.N.RA.S. 220,  845  Peterson, B.A., Murdin, P., Wallace, P.T., Manchester, R.N., Penny, A.J., Jordan, Hartley, K.F., and King, D. 1978, Nature 276, 475  Pilachowski, C.A., Sneden, C, and Canterna, R. 1980, in IAU symposium 85, Star Clus ed. J. Hesser (Dordrecht: Reidel), p. 467 Pilachowski, C.A., Sneden, C, and Wallerstein, G. 1983, Ap. J. Suppl. 52, 241  Press, W.H., Flannery, B.P., Teukolsky, S.A., and Vetterling, W.T. 1987, in Numeric Recipes. The Art of Scientific Computing (Cambridge: Cambridge University Press) Richer, H.B., and Fahlman, G.G. 1986, Ap. J. 304, 273 Richer, H.B., and Fahlman, G.G. 1988, Ap. J. 325, 218 Richer, H.B. 1991, Preprint  Richer, H.B., and Fahlman, G.G., 1991, in ASP Conference Series Volume 13, The Forma and Evolution of Star Clusters, ed. Kenneth James (BookCrafters inc.), p. 219  Richer, H.B., Fahlman, G.G., Buonanno, R., Fusi Pecci, F., Searle, L., and Thompson 1991, Ap. J. 381 in press (RFBFST) Rodgers, A.W., and Eggen, J. 1974,  PA.S.P. 86,  742  Sandage, A. 1969, Ap. J. 158, 1115 Sato, T., Richer, H.B., and Fahlman, G.G. 1989, A. J. 98, 1335 Savage, B.D., and Mathis, J.S. 1979, Ann. Rev. Searle, L., and Zinn, R. 1978, Ap. J. 225, 357  Astron. Astrophys. 17,  73  BIBLIOGRAPHY  114  Spitzer, L. (Jr.), and Mathieu, R.D. 1980, Ap. J. 241, 618 Spitzer, L. (Jr.) 1987, in Dynamical Evolution of University Press)  Globular Clusters (Princeton:  Stetson, P.B. 1987a, DAOPHOT Observatory)  (Victoria: Dominion Astrophysical  Users Manual  Stetson, P.B. 1987b, P.A.S.P.  99,  VandenBerg, D.A. 1983, Ap.  J. Suppl. 51,  Princeton  191 29  VandenBerg, D. A., and Bell, R.A. 1985, Ap.  J. Suppl. 58,  561 (VB85)  VandenBerg, D.A., Bolte, M., and Stetson, P.B. 1990, A. J. 100, 445 Webbink, R.F. 1985, in IAU symposium 113, Dynamics of Star and P. Hut (Dordrecht: Reidel), p. 541 White, R.E., and Kraft, T.T. 1972, PA.S.P. Zinn, R.J. 1980a, Ap.  J. Suppl. 42,  84,  298  19  Zinn, R.J. 1980b, Ap. J. 241, 602 Zinn, R.J., and West, M.J. 1984, Ap. Zinn, R.J. 1985, Ap. J. 293, 424  J. Suppl. 55,  45  Clusters,  eds. J. Goodman  

Cite

Citation Scheme:

        

Citations by CSL (citeproc-js)

Usage Statistics

Share

Embed

Customize your widget with the following options, then copy and paste the code below into the HTML of your page to embed this item in your website.
                        
                            <div id="ubcOpenCollectionsWidgetDisplay">
                            <script id="ubcOpenCollectionsWidget"
                            src="{[{embed.src}]}"
                            data-item="{[{embed.item}]}"
                            data-collection="{[{embed.collection}]}"
                            data-metadata="{[{embed.showMetadata}]}"
                            data-width="{[{embed.width}]}"
                            async >
                            </script>
                            </div>
                        
                    
IIIF logo Our image viewer uses the IIIF 2.0 standard. To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.831.1-0053473/manifest

Comment

Related Items