Australia Telescope Compact Array Survey of Candidate Ultra-Compact and Buried HII Regions in the Magellanic Clouds R´emy Indebetouw, Kelsey E. Johnson University of Wisconsin Astronomy Department, 475 N. Charter St., Madison, WI 53705 ([email protected], [email protected]) and Peter Conti JILA, University of Colorado, 440-UCB, Boulder, CO, 80309 ([email protected]) ABSTRACT We present a systematic survey for ultracompact H II (UCH II ) regions in the Magellanic Clouds. Understanding the physics of massive star formation (MSF) is a critical astrophysical problem. The study of MSF began in our galaxy with surveys of UCH II regions, but before now this has not been done for other galax- ies. We selected candidates based on their Infrared Astronomical Satellite (IRAS) colors and imaged them at 3 and 6 cm with the Australia Telescope Compact Array (ATCA). Nearly all of the observed regions contain compact radio sources consistent with thermal emission. Many of the sources are related to optically visible H II regions, and often the radio emission traces the youngest and densest part of the H II region. The luminosity function and number distribution of Ly- man continuum fluxes of the compact radio sources are consistent with standard arXiv:astro-ph/0407516v1 23 Jul 2004 stellar and cluster initial mass functions. This type of systematic assessment of IRAS diagnostics is important for interpreting Spitzer Space Telescope data, which will probe similar physical scales in nearby galaxies as IRAS did in the Magellanic Clouds. Subject headings: stars: formation — HII regions — Magellanic Clouds 1. Introduction Massive stars play a major role in the dynamical evolution of galaxies: they are respon- sible for the ionization of the interstellar medium, their stellar winds and supernovae are –2– dominant sources of mechanical energy, their ultraviolet radiation powers far-infrared lumi- nosities through the heating of dust, and they are a main driver of chemical evolution in the universe via supernova explosions at the end of their lives. However, despite the significant role of massive stars throughout the universe, their birth is not well understood, and we are only beginning to piece together a scenario for the youngest stages of massive star evolution. We have made some progress understanding the early stages of massive star formation in the Galaxy, but the current knowledge about the early stages of massive star evolution in other environments is mediocre at best. The reasons for this dearth of information about extremely young massive stars in other galaxies are predominantly twofold: (1) the earliest stages of massive star evolution are deeply enshrouded and inaccessible in the optical and near-infrared regimes; and (2) high spatial resolutions are necessary to disentangle the individual massive stars from their surrounding environment and background contamination. Radio observations using the Aus- tralia Telescope Compact Array (ATCA) are uniquely capable of overcoming these obstacles for the galaxies nearest to our own, the Large and Small Magellanic Clouds. Synthesis imag- ing observations at centimeter wavelengths are sensitive to the free-free emission emitted by the densest compact ionized regions, and the centimeter spectral index can be used to differentiate those extremely dense regions suffering from self-absorption. We have conducted a survey of candidate UCH II regions in the Magellanic Clouds with high (1–2′′) spatial resolution at 3 and 6 cm using the ATCA. Section 2 contains a description of the FIR-based candidate selection, the radio observing strategy, data reduction, and lists the detected sources. Individual compact H II regions with interesting radio morphologies are described in the Appendix, and the statistical properties of the population as a whole are analyzed in section 3. 2. Observations 2.1. Far Infrared Target Selection One goal of this project was to test the effectiveness of using far-infrared (FIR) colors at fairly low spatial resolution to discover massive protostellar objects. This exercise using IRAS data in the Magellanic Clouds is particularly interesting because the Spitzer Space Telescope probes similar spatial scales in nearby galaxies as IRAS did in the Magellanic Clouds. The primary target selection criterion for this survey was thus FIR color. Wood & Churchwell (1989b, hereafter WC89) determined that the majority of Galactic ultra-compact –3– H II (UCH II ) regions have particular mid- to far-infrared colors which distinguish them from other astronomical objects. They indicated (their Figure 1b) that most sources with IRAS colors Fν(60µm)/Fν(12µm) ≥ 20. and Fν(25µm)/Fν(12µm) ≥ 3.7 are UCH II s (these colors correspond to black-body temperatures .105 and .160 K, respectively). More importantly in the context of this work, they found that there was little color contamination – compared to a random selection of Galactic IRAS sources, most of the objects in that region of color-color space were UCH II regions. To exploit this fact, we first combined the IRAS point source catalog with a multi- wavelength Parkes radio continuum survey (Filipovi´c, Haynes, White, & Jones 1998, and references therein). The Parkes surveys have a resolution of 2.8 arcmin at 8.55 Ghz (which corresponds to a linear distance of ∼ 49 pc in the SMC and ∼ 41 pc in the LMC), so we took as a match to each radio source the brightest IRAS source within 2.5 arcminutes. We then examined the FIR to radio spectral energy distribution of all the objects with slightly broader colors than the WC89 criteria: F24µm/F12µm > 1 and F60µm/F12µm > 10. This sample contained 98 objects in the LMC (47 with the stricter WC89 colors) and 22 in the SMC (8 with WC89 colors). Finally, we chose those targets which were consistent with having flat or inverted radio emission at λ ∼1–5 cm. Thus our sample is FIR-color selected, but objects which are highly nonthermal at cm wavelengths at the resolution of the Parkes survey are excluded. This last criterion was imposed to bias against supernova remnants and other sources of nonthermal emission. Sources consistent with optically thin thermal α emission (slightly negative spectral index α ∼ −0.1 where Fν ∝ ν ) were kept as candidates under the presumption that H II regions observed at low resolution will be dominated by the optically thin envelopes, but may still contain compact, more optically thick sources. Figure 1 shows our target objects in an IRAS color-color space identical to that of WC89’s Figure 1b. The two targets which fall farthest outside of the color selection box were cases for which the IRAS flux densities or identification were ambiguous. 2.2. Radio Observations and Data Reduction Candidate fields were observed with the ATCA during two runs in May 2001 and May 2002. Simultaneous observations were obtained at 3 cm (8.6 GHz) and 6 cm (4.8 GHz) in the highest resolution configurations (6F in 2001, 6A in 2002). Each field was observed for roughly 3 hours, with some time added for a few of the faintest fields. We observed 19 of the 29 selected fields in the LMC, and 7 of the 8 selected fields in the SMC. For some of the fields, archival data was available, which allowed compact sources to be tabulated for 3 additional candidates. Our observations are summarized in Table 1, and the archival –4– Fig. 1.— Candidate UCH II regions in the Magellanic Clouds. All IRAS sources with Parkes radio detections in the Clouds are marked (crosses), our candidates in bold symbols (filled diamonds observed 5/2001, squares observed 5/2002, filled circles existing archival observations). The color selection box of WC89 is marked; we used a slightly relaxed color selection criterion in order to test the validity of that box at the lower spatial resolution available in the Magellanic Clouds. –5– observations which we completely re-reduced for consistency, are summarized in Table 2. The receivers simultaneously detect 4800 MHz (∼6 cm) and 8640 MHz (∼3 cm), with primary beams (field of view) of 10 and 5 arcminutes, respectively. Each frequency had a bandwidth of 128 MHz, and the resulting theoretical RMS noise for these observations is ≈ 0.07 mJy/beam. The absolute flux density scale was determined by assuming a flux density for PKSB1934-638 of 5.83 Jy at 4.8 GHz and 2.84 at 8.6 GHz, and we estimate that the resulting flux density scale is uncertain by less than ∼ 5%. Calibration was carried out using the Miriad data reduction package, including gain, phase, and polarization calibrations of the flux density calibrator and phase calibrators. Following the calibration, data were loaded into the Astronomical Image Processing System (AIPS) where they were inverted and cleaned using the task imagr. In order to optimize the synthesized beam at 4.8 GHz, a fairly uniform weighting was used (robust=0), while a more natural weighting was using at 8.6 GHz (robust=3) in order to increase the sensitivity. In addition, the uv range of the 4.8 GHz data was limited in order to better match the shortest uv coverage of the 8.6 GHz data. It should be noted that this does not create matched beams at 3cm and 6cm, but simply mitigates the influence of extended structure at 6cm; the 6cm observations still have denser uv coverage at short spacings, and the 3cm observations still have uv coverage at long spacings not available in the 6cm observations. Based on the resulting uv coverage, the largest angular scales to which these data are sensitive are . 30′′ at 6 cm and . 20′′ at 3 cm. The images at both frequencies were then corrected for the primary beam sensitivity. Radio flux densities were measured in AIPS++ using two methods. The first method used the viewer program to place identical apertures around the 3 cm and 6 cm images, and background levels were estimated using surrounding annuli.