PUBLICATIONS OF THE ASTRONOMICAL SOCIETY OF THE PACIFIC, 111 : 1392È1397, 1999 November ( 1999. The Astronomical Society of the PaciÐc. All rights reserved. Printed in U.S.A.
The Age of the Sculptor Dwarf Spheroidal Galaxy from Imaging with WFPC21 JACKIE MONKIEWICZ,2 JEREMY R. MOULD,2 JOHN S. GALLAGHER III,3 JOHN T. CLARKE,4 JOHN T. TRAUGER,5 CARL GRILLMAIR,6 GILDA E. BALLESTER,4 CHRISTOPHER J. BURROWS,7 DAVID CRISP,5 ROBIN EVANS,5 RICHARD GRIFFITHS,8 J. JEFF HESTER,9 JOHN G. HOESSEL,3 JON A. HOLTZMAN,10 JOHN E. KRIST,7 VICKI MEADOWS,5 PAUL A. SCOWEN,9 KARL R. STAPELFELDT,5 RAGVENDRA SAHAI,5 AND ALAN WATSON11 Received 1999 July 5 ; accepted 1999 July 22
ABSTRACT. From images taken with the Hubble Space T elescopeÏs WFPC2, we have obtained photometry of a Ðeld in the Sculptor dwarf spheroidal galaxy to 3 mag below the main-sequence turno†. We determine an age equal to that of the earliest globular clusters for the bulk of the stars in our Ðeld of view. We attempt to constrain the star formation history of the Sculptor dwarf by examining the main-sequence luminosity function. The presence of a half-dozen blue straggler candidates blueward of the turno† points to a possibly complex star formation history. However, the contribution of any intermediate-age population is difficult to measure conclusively, because of the uncertain origin of blue stragglers and the sparseness of the photometric sample.
L . Hierarchical models of formation depict the Galaxy as _ having formed by accreting smaller clumps of matter. Dwarf spheroidal galaxies may be the evolved remnants of these elementary building blocks and therefore tracers of the earliest epochs of star formation. Yet the star formation histories of the dwarf spheroidal satellites of the Milky Way are diverse (Da Costa 1998 ; Mateo 1998) : there are relatively simple systems, such as Draco, which is dominated by an old population of stars and is well-described by single burst of star formation, and there are systems with complex histories, such as Carina, which has a substantial intermediate-age population with a smaller fraction of older stars (Hurley-Keller, Mateo, & Nemec 1998 and references therein). Both Carina and Draco possess predominantly red horizontal branches, despite the fact that they are very metal poorÈthat is, they exhibit the second parameter e†ect. Although this is a clear characteristic, there are blue stars in Sculptor too (Kaluzny et al. 1995). The fact that the bulk of DracoÏs stars appears to be within 2 billion years of the oldest globular clusters (Grillmair et al. 1998) indicates that if age is the second parameter, horizontal-branch structure must depend very sensitively on it. The Sculptor dwarf spheroidal galaxy is a signiÐcantly Ñattened galaxy located toward the Milky WayÏs south pole. Stellar photometry in this galaxy reached the horizontal branch in the work of Kunkel & Demers (1977) and the main sequence in the Da Costa (1984) study. With a total luminosity exceeding 106 L , it is larger than either Carina _
1. INTRODUCTION Dwarf spheroidal galaxies are enigmatic objects, characterized by extremely low surface brightnesses and often high mass-to-light ratios. The Milky Way has some nine such companions, ranging in luminosity from 2 ] 105 to 2 ] 107 ÈÈÈÈÈÈÈÈÈÈÈÈÈÈÈ 1 Based on observations made with the NASA/ESA Hubble Space T elescope, obtained at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. 2 Research School of Astronomy and Astrophysics, Institute of Advanced Studies, Australian National University, Mount Stromlo Observatory, Private Bag, Weston Creek Post Office, ACT 2611, Australia ; jrm=mso.anu.edu.au. 3 Department of Astronomy, University of WisconsinÈMadison, 475 North Charter Street, Madison, WI 53706-1582 ; jsg=astro.wisc.edu. 4 Department of Atmospheric, Oceanic, and Space Sciences, University of Michigan, 2455 Hayward Street, Ann Arbor, MI 48109. 5 Jet Propulsion Laboratory, 4800 Oak Grove Drive, Mail Stop 183900, Pasadena, CA 91109. 6 IPAC, California Institute of Technology, Pasadena, CA 91125. 7 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218. 8 Department of Physics, Carnegie-Mellon University, Wean Hall, 5000 Forbes Avenue, Pittsburgh, PA 15213. 9 Department of Physics and Astronomy, Arizona State University, Tyler Mall, Tempe, AZ 85287. 10 Department of Astronomy, New Mexico State University, Box 30001, Department 4500, Las Cruces, NM 88003-8001. 11 Instituto de Astronomia-UNAM, 58090 Morelia, Michoacan, Mexico.
1392
AGE OF SCULPTOR DWARF SPHEROIDAL GALAXY
1393
two or more sizable clouds of H I inside the tidal radius and within the velocity range of the dwarf galaxy. Sculptor is apparently unique in this respect, although it should be noted that the H I velocity distribution of this region is quite complex (Putman et al. 1998 ; M. Putman 1999, private communication). The location of our Hubble Space T elescope Wide Field Planetary Camera 2 (WFPC2) Ðeld with respect to the Sculptor dwarf and CarignanÏs cloud is shown in Figure 1. It was chosen so that the photometry would go as deep as possible, undisturbed by confusion and crowding e†ects, rather than for proximity to this di†use gas. In this paper we apply the WFPC2 to Sculptor in order to measure the properties of the dwarf galaxyÏs mainsequence turno†. The photometry we obtain reaches to 3 or 4 mag below the turno† and allows us to compare the Sculptor dwarf spheroidal galaxy with its smaller northern counterpart in Draco.
FIG. 1.ÈThe WFPC2 Ðeld (dark gray) and the H I surface density contours of Carignan et al. (white) are overlaid on the STScI Digitized Sky Survey image of the Sculptor dwarf. This Ðeld is 40@ square.
or Draco, although Sculptor is metal poor (on average [Fe/H] \ [1.8 ^ 0.1 ; Mateo 1998). Unlike Carina, neither Sculptor nor Draco has stars signiÐcantly exceeding the luminosity of the red giant branch tip (Frogel et al. 1982). This lack of upper asymptotic giant branch (AGB) stars leads to the expectation that Sculptor also lacks a signiÐcant intermediate-age population. The carbon stars found by Aaronson & Mould (1985) and Azzopardi (1994) are at or below the red giant branch tip in luminosity. Sculptor is of particular interest because it may possess some neutral hydrogen. Carignan et al. (1998) discovered
TABLE 1
u4350201r u4350202r u4350203r u4350204r u4350205r u4350206r u4350207r u4350208r
........... ........... ........... ........... ........... ........... ........... ...........
Filter
HJD (start)
Exposure Length (s)
F555W F555W F555W F555W F814W F814W F814W F814W
2,450,795.7786 2,450,795.7946 2,450,795.8411 2,450,795.3578 2,450,795.9085 2,450,795.9252 2,450,795.9779 2,450,795.9946
1200 1200 1300 1300 1300 1300 1300 1300
NOTE.ÈUT date \ 1997 December 13. R.A. \ 01h00m07s. 42, decl. \ [33¡27@34A. 21. Roll angle(V3) \ 235¡.
1999 PASP, 111 : 1392È1397
A Ðeld well outside the core of Sculptor was selected in order to avoid image crowding problems. The images were taken on 1997 December 13 by the WFPC2 instrumentation aboard the Hubble Space T elescope. The observations comprise two 1200 s and two 1300 s images taken with the F555W Ðlter and four 1300 s images taken with the F814W Ðlter (see Table 1). The Ðeld of view is located
[email protected] from the center of the Sculptor dE, or a projected distance of 330 pc for a line-of-sight distance of 80 kpc. The raw observations were preprocessed via the WFPC2 standard pipeline calibration. Bad pixel and vignetting masks were then applied to the images, and intensity corrections were made for Ðeld distortions (Stetson et al. 1998).
3. PHOTOMETRY
LOG OF OBSERVATIONS FOR SCULPTOR dE
Archive Root Name
2. OBSERVATIONS
Apparent magnitudes were derived from the standard WFPC2 point-source functions developed for the H Key 0 Project using the DAOPHOT/ALLFRAME ensemble of programs (see Stetson 1994 for a complete description of this software). To reduce the e†ects of cosmic rays, the exposures in each Ðlter were combined together to produce the median of the four individual images in F555W and F814W. The initial star list was produced from by combining all eight frames. The WFPC2 mosaic of the combined observations is included in Figure 2. The detected stars were subtracted from the combined image, and a second detection pass was carried out to pick out point sources obscured
1394
MONKIEWICZ ET AL.
FIG. 2.ÈWFPC2 mosaic of our Sculptor Ðeld. This Ðeld is
[email protected] on a side.
by the skirts of brighter stars. After the convergence of the ALLFRAME program, obvious detections of background galaxies were removed from the Ðnal star list by hand. Stellar image quality constraints further eliminated spurious detections. Including contributions from Sculptor and the Galaxy, 400 stars were obtained from the three combined WF frames and 20 stars from the PC frame. Stars with V \ 27.5 are listed in Table 2. Fiducial stars for astrometry are listed in Table 3. In order to transform our photometry to the standard system, it is necessary to include the stellar light out to radial distances of 0A. 5 from the centroid of the point-source image. To this purpose, 20È30 bright, isolated stars were selected from each chip for the production of growth curves. The template growth curves and aperture corrections were constructed and applied to the rest of the photometry using DAOGROW and associated programs (Stetson 1990). Note that we have used standard transformations between F555W, F814W, and V , I (Hill et al. 1998). Holtzman et al.
(1995) found that synthetic photometry of models with metal abundances [Fe/H] \ [2.0 and 0.0 were indistinguishable for 0.0 \ V [I \ 1.0 in the two-color diagrams one can form from these magnitudes. Because the Ðeld is sufficiently uncrowded, aperture photometry was performed on a sample of bright stars from the Ðnal star list. Comparison with the ALLFRAME-derived magnitudes indicates that there is good chip-to-chip agreement between the two methods.
4. THE COLOR-MAGNITUDE DIAGRAM Figure 3 is a color-magnitude diagram like that of a globular cluster, and the natural comparison is with the metal-poor clusters of our Galaxy. The anticipated contamination of the Ðeld by foreground stars is negligible. The Galactic model of Bahcall and Soneira predicts some 20 foreground Galactic stars in our Ðeld for m \ 28 V 1999 PASP, 111 : 1392È1397
AGE OF SCULPTOR DWARF SPHEROIDAL GALAXY
1395
TABLE 2 STARS DETECTED IN WFPC2 SCULPTOR FIELD WITH m \ 27.5 V STAR
COORDINATES
MAGNITUDE
Chip
ID
x
y
V
V [I
PC . . . . . . PC . . . . . . PC . . . . . . PC . . . . . . PC . . . . . . PC . . . . . . PC . . . . . . PC . . . . . . PC . . . . . . PC . . . . . . PC . . . . . . PC . . . . . . PC . . . . . . PC . . . . . . PC . . . . . .
10325 800 76 502 708 564 419 584 594 703 917 157 267 340 420
671.02 230.69 364.99 714.24 671.93 230.22 330.33 308.73 666.03 484.24 469.75 528.85 635.96 280.15 353.34
568.76 705.78 120.33 446.68 630.23 496.94 373.15 520.38 525.79 625.43 792.42 192.03 264.90 317.87 375.35
19.33 26.33 26.02 26.66 26.05 24.62 24.70 24.52 24.65 25.41 25.32 22.90 24.25 24.51 23.55
1.09 0.86 1.14 1.07 0.88 0.70 0.72 0.70 0.61 0.84 0.82 0.75 0.72 0.71 0.61
NOTE.ÈTable 2 is published in its entirety in the electronic edition of the PASP. A portion is shown here for guidance regarding its form and content.
(Ratnatunga & Bahcall 1985). The majority of these will be halo M dwarfs and will lie redward of the Sculptor main sequence. In Figure 3 we overlay the Ðducial of the M68 main sequence, as derived from Walker (1994). We shift the Ðducial to the distance of Sculptor, taking distance moduli of 14.8 and 19.5 mag for M68 and Sculptor, respectively (Djorgovski 1992 ; Kunkel & Demers 1977). Using the TABLE 3 STELLAR POSITIONS IN WFPC2 SCULPTOR FIELD COORDINATES (J2000)
STAR Chip
ID
PC . . . . . . . . PC . . . . . . . . PC . . . . . . . . WF2 . . . . . . WF2 . . . . . . WF2 . . . . . . WF3 . . . . . . WF3 . . . . . . WF3 . . . . . . WF4 . . . . . . WF4 . . . . . . WF4 . . . . . .
10325 157 554 322 1014 775 406 865 592 652 787 1225
R.A. 1 1 1 1 1 1 1 1 1 1 1 1
00 00 00 00 00 00 00 00 00 00 00 00
09.97 08.54 09.34 09.89 11.65 10.76 05.14 03.52 04.75 05.50 08.11 05.85
Decl. [33 [33 [33 [33 [33 [33 [33 [33 [33 [33 [33 [33
27 27 27 27 28 28 28 28 27 26 26 26
04.9 08.3 28.2 53.2 39.1 22.6 23.9 00.3 59.9 47.2 45.4 14.4
NOTE.ÈUnits of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds.
1999 PASP, 111 : 1392È1397
FIG. 3.ÈColor-magnitude diagram of the WFPC2 Ðeld. We can Ðt the Sculptor main sequence by shifting the M68 main sequence (shown as solid symbols) by 4.66 mag in V and [0.04 in V [I.
method of Schlegel, Finkbeiner, & Davis (1998), we determined a reddening adjustment of E(V [I) \ 0.04 from the DIRBE sky maps. The metal abundance of M68 is [Fe/H] \ [2.09 (Djorgovski 1992). A change of 0.1 dex in [Fe/H] produces a change of 0.006 in V [I at M \ 4.0 in V the Yale isochrones (Green, Demarque, & King 1987). Although the overlaid Ðducial sequence for M68 lies toward the blue side of the Sculptor main sequence and giant branch base, the o†set is within the photometric calibration uncertainties of 0.03 mag (Stetson et al. 1998). These results are consistent with the proposition that there is no di†erence between the average ages of M68 and the stellar population of Sculptor. In order to put quantitative limits on what ““ no di†erence ÏÏ signiÐes, Yale isochrones for metallicity Z \ 2 ] 10~4 and helium abundance Y \ 0.23 were shifted to a common turno† color and magnitude and Ðtted to the Sculptor main sequence for ages of 13, 15, and 17 Gyr in Figure 4. According to Da Costa & Armandro† (1990), Yale isochrones are satisfactory Ðts to cluster giant branches, provided the Frogel, Persson, & Cohen (1983) horizontal-
1396
MONKIEWICZ ET AL.
FIG. 4.ÈColor-magnitude diagram of the WFPC2 Ðeld. We can Ðt the Sculptor main sequence by Yale isochrones with [Fe/H] \ [2.0, Y \ 0.23 for 13, 15, and 17 Gyr, registered to the bluest point on the isochrone and to the turno† of Sculptor, following VandenBerg, Bolte, & Stetson (1990). These isochrones represent the age limits of acceptable Ðts to the giant branch.
branch distance scale is adopted, which is the case here. Within the range of the small number statistics in our sample, the 13 and 17 Gyr isochrones lie outside the limits of probable Ðts to the red giant branch of Sculptor. Mainsequence Ðts, using VandenBerg et al. (1999) and Padova isochrones (Bertelli et al. 1994) of comparable metal and helium abundance, are generally in good agreement with this result. Including uncertainties in the metallicity determination, and requiring that our conclusion be independent of the choice of evolutionary model, we Ðnd an age of 15 ^ 2 Gyr for Sculptor. We emphasize that it is the ^2 that we are interested in determining in this way, rather than the absolute age. 5. DISCUSSION The age of the system is the most direct result of this new deep look at Sculptor, but we should keep in mind the
possibility that the stellar population varies radially (Da Costa et al. 1996), as is seen in the Small Magellanic Cloud (SMC ; Hatzidimitriou, Cannon, & Hawkins 1993), in NGC 147 (Han et al. 1997), and in Andromeda I (Armandro† et al. 1993). We have estimated an age of 15 ^ 2 Gyr for the bulk of the main-sequence stars. This age is comparable to the ages of the oldest Galactic globular clusters. However, the star formation histories of most dwarf galaxies are actually more complex (Carignan 1999) than the simple stellar population model of Renzini & Buzzoni (1986), and it is possible that this is true of Sculptor as well. The main-sequence luminosity function (MSLF) is a simple gauge of the lifetime star formation history, although in a composite population the MSLF can be degenerate over a very wide range of age/metallicity combinations. Kaluzny et al. (1995) show that there are no main-sequence stars with m \ 21.0 in the central area of the galaxy. CalcuV lations based on the Yale luminosity functions (Green et al. 1987) were carried out to compare the predictions of three possible star formation histories of Sculptor with the MSLF. The small number statistics of our present data permit a range of star formation histories to be taken as plausible Ðts to the data. The absence of stars brighter than V \ 21.0 rules out star formation more recent than 2 Gyr ago. A continuous star formation model is therefore inappropriate to the central regions of the galaxy at least. A spread in chemical composition in Sculptor (Da Costa 1984) would arise naturally from an extended history of star formation. While the placement of our Ðeld of view far from the center of Sculptor averts the problems associated with crowding, it severely restricts the number of stars brighter than the main-sequence turno†. The result is a sparsely populated giant branch and a statistically difficult number of stars in the blue straggler region of the color-magnitude diagram. This inhibits our ability to draw meaningful conclusions about the contribution of intermediate-age populations. We simply make the point that it is possible that star formation has been active in Sculptor steadily or sporadically over its lifetime, and not simply in its initial burst, presumably until the event which resulted in its gas dispersal. But we note that star formation activity inferred from the MSLF is an upper limit on the true level of activity ; the chance occurrence of a background quasar or a halo white dwarf in our Ðeld would lead to an overestimation of the amount of star formation during the last 5 Gyr. It is also possible that some or all of the stars in the 22d magnitude bins in Figure 3 are blue stragglers formed from stellar mergers rather than from a separate burst of star formation. If signiÐcantly more than 10% of the star formation in Sculptor occurred more recently than the SMC globular cluster NGC 121 was formed (11.5 Gyr ago according to Stryker, Da Costa, & Mould 1985), we would expect to see a few upper AGB carbon stars in the galaxy.
1999 PASP, 111 : 1392È1397
AGE OF SCULPTOR DWARF SPHEROIDAL GALAXY If the 22d magnitude main-sequence stars are in fact ordinary Sculptor main-sequence stars, the star formation history which is apparent in Sculptor is reminiscent in a very small way of what is played out fully in the Fornax dwarf (Grillmair et al. 1999 ; Stetson 1997). The range of star formation histories exhibited by dwarf spheroidal galaxies may run continuously from systems such as Draco, in which activity was strongly peaked in the past, to objects with a basically steady star formation rate which ceased recently, such as Fornax. Where exactly Sculptor is in this continuum remains to be seen. Clearly, more data are required to pursue these matters further. A detailed comparison of histories of Sculptor and
1397
Fornax (and other systems) is possible and well motivated. Such a study would require adequate spatial sampling of these galaxies, which are large compared with the WFPC2 Ðeld. Ground-based imaging is particularly important for capturing wide Ðelds where relatively rare younger stars are to be found. This research was carried out by the WFPC2 Investigation DeÐnition Team for JPL and was sponsored by NASA through contract NAS7-1260. We thank Don VandenBerg and Cesare Chiosi for communicating isochrones for use in this work. We are grateful to Gary Da Costa for reading and commenting on the manuscript.
REFERENCES Aaronson, M., & Mould, J. 1985, ApJ, 290, 191 Armandro†, T., Da Costa, G., Caldwell, N., & Seitzer, P. 1993, AJ, 106, 986 Azzopardi, M. 1994, in The Third ESO/CTIO Workshop on the Local Group, ed. A. Layden, R. Smith, & J. Storm (Garching : ESO), 129 Bertelli, G., Bressan, A., Chiosi, C., Fagotto, F., & Nasi, E. 1994, A&AS, 106, 275 Carignan, C. 1999, Publ. Astron. Soc. Australia, 16, 18 Carignan, C., Beaulieu, S., Cote, S., Demers, S., & Mateo, M. 1998, AJ, 116, 1690 Da Costa, G. S. 1984, ApJ, 285, 483 ÈÈÈ. 1998, in Stellar Astrophysics for the Local Group, ed A. Aparacio, A. Herrero, & F. Sanchez (Cambridge : Cambridge Univ. Press), 351 Da Costa, G. S., & Armandro†, T. 1990, AJ, 100, 162 Da Costa, G. S., Armandro†, T., Caldwell, N., & Seitzer, P. 1996, AJ, 112, 2576 Djorgovski, S. G. 1992, in ASP Conf. Ser. 50, Structure and Dynamics of Globular Clusters, ed. S. G. Djorgovski & G. Meylan (San Francisco : ASP), 373 Frogel, J., Blanco, V., McCarthy, M., & Cohen, J. 1982, ApJ, 252, 133 Frogel, J., Persson, S. E., & Cohen, J. 1983, ApJS, 53, 713 Green, E., Demarque, P., & King, C. 1987, The Revised Yale Isochrones and Luminosity Functions (New Haven : Yale Univ. Obs.)
1999 PASP, 111 : 1392È1397
Grillmair, C., et al. 1998, AJ, 115, 144 ÈÈÈ. 1999, AJ, 117, 167 Han, M., et al. 1997, AJ, 113, 1001 Hatzidmitriou, E., Cannon, R., & Hawkins, M. 1993, MNRAS, 261, 873 Hill, R. J., et al. 1998, ApJ, 496, 648 Holtzman, J., et al. 1995, PASP, 107, 1065 Hurley-Keller, D., Mateo, M., & Nemec, J. 1998, AJ, 115, 1840 Kaluzny, J., Kubiak, M., Szymanski, M., Udalski, A., Krzeminski, W., & Mateo, M. 1995, A&AS, 112, 407 Kunkel, W., & Demers, S. 1977, ApJ, 214, 21 Mateo, M. 1998, ARA&A, 36, 435 Putman, M. E., et al. 1998, Nature, 394, 752 Ratnatunga, K. U., & Bahcall, J. N. 1985, ApJS, 59, 63 Renzini, A., & Buzzoni, A. 1986, Spectral Evolution of Galaxies : Proceedings of the Fourth Workshop, Erice, Italy, ed. C. Choisi & A. Renzini (Dordrecht : Kluwer), 195 Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 Stetson, P. B. 1990, PASP, 102, 932 ÈÈÈ. 1994, PASP, 106, 250 ÈÈÈ. 1997, Baltic Astron., 6, 3 Stetson, P. B., et al. 1998, ApJ, 508, 491 Stryker, L., Da Costa, G., & Mould, J. 1985, ApJ, 298, 544 VandenBerg, D. A., Bolte, M., & Stetson, P. B. 1990, AJ, 100, 445 VandenBerg, D. A., Swenson, F., Rogers, F., Iglesias, C., & Alexander, D. 1999, in preparation Walker, A. 1994, AJ, 108, 555
