EXTINCTION MAP OF BAADE'S WINDOW 1 K. Z. Stanek Princeton University Observatory, Princeton, NJ 08544{1001 e-mail: [email protected] ABSTRACT Recently Wo(cid:19)zniak & Stanek (1996) proposed a new method to investigate 6 9 interstellar extinction, based on two band photometry, which uses red clump 9 stars as a means to construct the reddening curve. I apply this method to 1 n the color-magnitude diagrams obtained by the Optical Gravitational Lensing a 0 2 Experiment (OGLE) to construct an extinction map of (40) region of Baade's J 4 Window, with resolution of (cid:24) 30 arcsec. Such a map should be useful for studies of this frequently observed region of the Galactic bulge. The map and 7 software useful for its applications are available via anonymous ftp. The total 3 0 2 1 extinction AV varies from 1:26 mag to 2:79 mag within the (40) (cid:12)eld of view 2 1 centered on ((cid:11)2000;(cid:14)2000) = (18:03:20.9,{30:02:06), i.e. (l;b) = (1:001;(cid:0)3:885). 5 The ratio AV=E(V (cid:0)I) = 2:49(cid:6)0:02 is determined with this new method. 9 / h p - Subject headings: Galaxy: general { stars: Hertzsprung-Russell diagram { stars: o r statistics t s a 1. INTRODUCTION The Optical Gravitational Lensing Experiment (OGLE, Udalski et al. 1993; 1994) is an extensive photometric search for the gravitational microlensing events towards the Galactic bulge. However, the wealth of data it provides (Szyman(cid:19)ski & Udalski 1993) may be applied to a host of separate problems. Recently Wo(cid:19)zniak & Stanek (1996) proposed a new method for investigating the value of the coe(cid:14)cient of the selective extinction. This method uses red clump stars as a means to construct the reddening curve. Red clump stars are the equivalent of the horizontal-branch stars in metal-rich populations, i.e. they are relatively low-mass stars burning helium in their cores. From observations and also from 1 On leave fromN. Copernicus AstronomicalCenter, Bartycka 18, Warszawa 00{716,Poland { 2 { stellar evolution theory (Castellani, Chie(cid:14), & Straniero 1992) we expect the bulge red clump stars to be relatively bright and have a narrow luminosity distribution with weak dependence on the metallicity. Therefore red clump stars form a suitable population with which to investigate the e(cid:11)ects of the interstellar extinction. I apply the method of Wo(cid:19)zniak & Stanek to the OGLE CMDs to map the interstellar extinction for Baade's Window region of the Galactic bulge. To my knowledge, this is the (cid:12)rst attempt to map in detail the extinction in this frequently observed region of the Galactic bulge. Blanco et al. (1984) divided their photographic plate of Baade's Window into four sub(cid:12)elds according to the surface density of stars. However, their division was done in rather arbitrary fashion, by visual inspection of the plate. The method I use in this paper, based on Wo(cid:19)zniak and Stanek (1996), is entirely quantitative. Terndrup et al. (1995) determined the reddening values for regions de(cid:12)ned by Blanco et al. (1984) using medium-resolution spectra of K giants (which, from their color-magnitude diagram, seem overlap largely with what we call \red clump" stars). In Section 2, I discuss the data used in this paper. In Section 3, I use the method of Wo(cid:19)zniak & Stanek to construct the extinction map for nine (cid:12)elds in Baade's Window region of the Galactic bulge. In Section 4, I apply di(cid:11)erent checks to the results and discuss possible applications of the extinction map. 2. THE DATA I use two somewhat separate data products created by the OGLE collaboration. Szyman(cid:19)ski & Udalski (1993) constructed a database of photometric measurements in V and I bands for the (cid:12)elds observed by the OGLE project. Udalski et al. (1993) present color-magnitude diagrams of 13 (cid:12)elds in the direction of the Galactic bulge. In this paper 0 2 I discuss nine (cid:12)elds in Baade's window (BW), which together cover (cid:24) (40) and contain 5 (cid:24) 5 (cid:2)10 stars. All observations were made using the 1-meter Swope telescope at the Las Campanas Observatory (operated by the Carnegie Institution of Washington) and a 0 0 2048(cid:2)2048 pixel Ford/Loral CCD detector with the pixel size 0:44 arcsec covering 15 (cid:2)15 (cid:12)eld of view. CMDs for the (cid:12)elds analyzed in this paper can be seen in Udalski et al. (1993). Most of each diagram is dominated by bulge stars, with distinct red clump, red giant and turn-o(cid:11) point stars. The part of the diagram dominated by the disk stars has been analyzed by Paczyn(cid:19)ski et al. (1994). Stanek et al. (1994) has used red clump dominated parts of the CMDs to (cid:12)nd a strong signature of the Galactic bar in the OGLE data. { 3 { Fig. 1.|Schematicmapof Baade'sWindow1deg(cid:2)1deg regionof theGalacticbulge. Large overlapping squares represent nine OGLE BW (cid:12)elds. Small dots show the grids on which the extinctionwas measured. The big circlecentered on the large black dot (globular cluster NGC 6522) represents the regions A, B, C observed by Blanco et al. (1984). The shaded annulus centered on NGC 6522 represents the region analyzed by Terndrup et al. (1995). The smaller black dot represents the globular cluster NGC 6528. The coordinates are given in epoch 2000:0. 3. THE ANALYSIS Wo(cid:19)zniak & Stanek (1996) describe in detail the method I use in this paper. To summarize, the database (Szyman(cid:19)ski & Udalski 1993) of V measurements is used to construct star density maps of each of nine Baade's Window (cid:12)elds, BW1-8 and BWC (see Fig.1). Assuming that the extinction a(cid:11)ects signi(cid:12)cantly the number of stars detected in V band in a given sub(cid:12)eld, arranging the sub(cid:12)elds in order of increasing number of database stars per sub(cid:12)eld should re(cid:13)ect a decreasing extinction in a sub(cid:12)eld. This indeed has been shown the case by Wo(cid:19)zniak & Stanek. The density maps are then used for separating CMDs of sub(cid:12)elds with di(cid:11)erent extinction for a given (cid:12)eld. The resulting CMDs are used to obtain the quantitative values of the o(cid:11)set on the CMD between the di(cid:11)erent sub(cid:12)elds, caused by di(cid:11)erential extinction. Wo(cid:19)zniak & Stanek use the red clump dominated parts of CMDs for determining the o(cid:11)sets { the clump is seen at fainter magnitudes and redder { 4 { Fig. 2.| Map of the interstellar extinction A in Baade's Window. The two black dots V correspond to the globular clusters NGC 6522 (center) and NGC 6528 (left). For detailed discussion see text. Fig. 3.| Image of the same 1deg(cid:2)1deg region as in Fig.1 and Fig.2 adapted from the Digitized Sky Survey available on the CD ROMs. Dark regions corresponds to regions of low density of stars, caused by higher extinction. The two globular clusters marked as black dots in Fig.2 can be seen near the center (NGC 6522) and to the left (NGC 6528) of the image. { 5 { colors in sub(cid:12)elds with higher extinction. This method does not depend on knowing the true luminosity function of the clump giants, the only assumption is that it is the same 0 0 within an 15 (cid:2) 15 OGLE (cid:12)eld. As a result, for each of the nine (cid:12)elds I obtain two independent 30(cid:2)30 maps of the relative extinction and reddening within this (cid:12)eld, with pixel size corresponding to 64 (cid:2)0:44 arcsec = 28:16 arcsec. The boundary of 64 CCD pixels was removedin each (cid:12)eld due to the shifts between V and I frames, resulting in some empty boundary areas. From the nine 30 (cid:2)30 maps of the di(cid:11)erential extinction and reddening I then create two common maps. First, the individual maps are re-mapped onto one common, uniform grid. As seen in Fig.1, there are overlap regions between the individual (cid:12)elds. These regions are used for matching the values of di(cid:11)erential extinction between di(cid:11)erent (cid:12)elds. They also provide an independent measure of error in our extinction values { one can check how well the extinction and reddening agree in the matching points. Except for the expected zero-point o(cid:11)set the matches are very good, with the rms scatter corresponding to 0:1 mag for the extinction and 0:04 mag for the reddening. This is about twice the formal errors obtained using Wo(cid:19)zniak and Stanek method. So far I have described a method for constructing maps of the di(cid:11)erential extinction and reddening within Baade's Window. To construct the (cid:12)nal maps, I need zero-points o(cid:11)sets of these variables for some part of the maps. Terndrup et al. (1995) have recently determinedE(V (cid:0)I) for sample of stars in BW and found that the region \A" seen in Fig.1 has an average of E(V (cid:0)I) = 0:59(cid:6)0:08 mag. Using matching points between my map and their region \A" I apply this value to the reddening map. In order to apply the o(cid:11)set in A , for every point in my maps I obtain A =E(V (cid:0)I) ratio (this is discussed in the next V V section). The averageratio of 2:49 is then used to multiplythevalue of E(V (cid:0)I)= 0:59 mag of Terndrup et al. (1995) and the resulting zero-point o(cid:11)set of A = 1:47 mag is then V;0 applied to the extinction map. The resulting map of extinction is shown as a greyscale plot in Fig.2. This map can be compared visually with the Palomar Survey plate of the same region seen in Fig.3. The correspondence of the high extinction regions in my map to low density of stars regions on the Palomar plate is striking. In Fig.4 I show the histogram of A values for the points on the common grid. Most of BW has the extinction below V A = 2:0 mag, but there is a tail of high extinction values, mostly in the BW3 (cid:12)eld, as sen V in Fig.2. The numerical values of both A extinction and E(V (cid:0)I) reddening maps are V available through the anonymous ftp serviceor by request. For the anonymous ftp service, use astro.princeton.edu server, go to stanek/Extinction subdirectory and retrieve the README (cid:12)le with the instructions. { 6 { 4. DISCUSSION In the previous section I have obtained maps of the interstellar extinction A and V 0 0 reddening E(V (cid:0)I) in roughly 40 (cid:2)40 region of Baade's Window. I will now apply several checks to see if these maps are indeed representative of the interstellar extinction in this region. As a (cid:12)rst check, for every point on the common grid (before applying the zero-point o(cid:11)sets) I plot the value of A as a function of E(V (cid:0)I). This is shown in Fig.5. By (cid:12)tting V a straight line to these points, one can obtain the value of A =E(V (cid:0)I), which in our case V is A =E(V (cid:0)I) = 2:49(cid:6)0:02. This is close to the value 2.6 of Dean, Warren, & Cousins V (1978) and Walker (1985). An additional check comes from the overlap regions between the nine BW (cid:12)elds. The (cid:12)elds are matched so that BWC is used to match other four adjacent (cid:12)elds, which in turn are used to match \corner" (cid:12)elds (see Fig.1). So, for each of the corner (cid:12)elds the two overlap regions can be used to see if, for example, matching the BW3 through the BW6 or through the BW7 is producing similar results. I have found no signi(cid:12)cant di(cid:11)erence in the resulting maps regardless of the way the matching procedure went. I have also used my extinction map to check what is the di(cid:11)erence in A extinction V between regions A, B and C of Blanco et al. (1984). I have found that A is (cid:24) 0:2 mag V higher in their region C than in their regions A and B, but there is a (cid:24) 0:5 mag scatter of A in all these regions. One should be cautious when using average extinction { V examination of the plates of BW shows that there may be structure in the reddening on scales even smaller than our \pixel" size, which is (cid:24) 30 arcsec. Another warning for users of the reddening and extinction maps { the regions close to the twoglobular clusters, NGC 6522 and NGC 6528, inBaade's Window are \contaminated" by the cluster's stars, which may result in extinction values too low for these regions. Indeed, examination of Fig.2 shows somewhat lower extinction around both clusters. Also, there is an artifact of CCD bleeding column visible in the BW8 (cid:12)eld as a vertical strip of higher extinction. One should employ caution when using the extinction map in these regions, and perhaps use average extinction from the neighboring regions. The FORTRAN program map.f available on anonymous ftp warns the user when the coordinates fall into these \bad" regions. I will now turn to some of the possible applications of the extinction and reddening maps produced in this paper. One such application was already explored by Wo(cid:19)zniak and Stanek { correcting the CMDs for the di(cid:11)erences in the extinction within the (cid:12)eld (see their Fig.5). Their dereddened CMD is less di(cid:11)use and shows both a much more compact red { 7 { Fig. 4.| Histogram of A values for the points on the common grid. V Fig. 5.| Correlation betweendi(cid:11)erentialreddening E(V (cid:0)I) and di(cid:11)erentialextinctionA . V { 8 { clump and narrower band of red giants, demonstrating that their method can be used to remove the e(cid:11)ects of the di(cid:11)erential extinction within the observed (cid:12)eld. Another possible application involves removing the e(cid:11)ects of interstellar extinction from samples of variable stars. It should be stressed here that my maps show only the total extinction and reddening towards Baade's Window, without any information on how they change with the distance. This problem was addressed to some extent by Arp (1965) and Paczyn(cid:19)ski et al. (1994). Another possible application of this extinction map would involve measuring the statistical properties of interstellar extinction distribution. In this paper I have constructed extinction and reddening maps of Baade's Window region of the galactic bulge, using OGLE CMD data. These maps show highly irregular extinction, with a di(cid:11)erence between the lowest and highest extinction regions of more than 1:5 mag in A . The errors of the maps are dominated by the error of the zero-point V o(cid:11)set, taken from Terndrup et al. (1995). When using the map to correct for the di(cid:11)erential extinction within BW region, the errors are (cid:24) 0:1 mag in A and (cid:24) 0:04 mag in E(V (cid:0)I). V The numerical values of both A extinction and E(V (cid:0)I) reddening maps are available V through the anonymous ftp service or by request. I would like to thank the OGLE collaboration for making their data available to me and Michal Szyman(cid:19)ski for his invaluable help with the OGLE software. I thank Przemek Wo(cid:19)zniak for some of the software created during our earlier collaboration and used for this project. I also thank Bohdan Paczyn(cid:19)ski and Slavek Rucin(cid:19)ski for indicating the usefulness of this project and comments on the paper. Wes Colley and the referee of the paper Donald Terndrup has also provided me with very useful comments on this paper. This project was supported with the NSF grant AST 9216494 and also with the NAS Grant-in-Aid of Research through Sigma Xi, The Scienti(cid:12)c Research Society. REFERENCES Arp, H., 1965, ApJ, 141, 43 Blanco, V. M., McCarthy, M. F., & Blanco, B. M., 1984, AJ, 89, 636 Castellani, V., Chie(cid:14), A., & Straniero, O., 1992, ApJS, 78, 517 Dean, J. F., Warren, P. R., & Cousins, A. W. J., 1978, MNRAS, 183, 569 Paczyn(cid:19)ski, B., Stanek, K. Z., Udalski, A., Szyman(cid:19)ski, M., Ka luz_ny, J., Kubiak, M., & Mateo, M., 1994, AJ, 107, 2060 { 9 { Stanek, K. Z., Mateo, M., Udalski, A., Szyman(cid:19)ski, M., Ka luz_ny, J., & Kubiak, M., 1994, ApJ, 429, L73 Szyman(cid:19)ski, M., & Udalski, A., 1993, Acta Astron., 43, 91 Terndrup, D. M., Sadler, E. M., & Rich, R. M., 1995, AJ, 110, 1774 Udalski, A., Szyman(cid:19)ski, M., Ka luz_ny, J., Kubiak, M., & Mateo, M., 1993, Acta Astron., 43, 69 Udalski, A., Szyman(cid:19)ski,M., Stanek, K. Z., Ka luz_ny, J., Kubiak, M., Mateo, M., Krzemin(cid:19)ski, W., Paczyn(cid:19)ski, B., & Venkat, R., 1994, Acta Astron., 44, 165 Walker, A. R., 1985, MNRAS, 213, 889 Wo(cid:19)zniak, P. R., & Stanek, K. Z., 1996, ApJ, in press This preprint was prepared with the AAS LATEX macros v4.0.