A publishing partnership

THE BLAST VIEW OF THE STAR-FORMING REGION IN AQUILA (ℓ = 45°, b = 0°)

, , , , , , , , , , , , , , , , , , , , , , , , , and

Published 2010 October 15 © 2010. The American Astronomical Society. All rights reserved.
, , Citation Alana Rivera-Ingraham et al 2010 ApJ 723 915 DOI 10.1088/0004-637X/723/1/915

0004-637X/723/1/915

ABSTRACT

We have carried out the first general submillimeter analysis of the field toward GRSMC 45.46+0.05, a massive star-forming region in Aquila. The deconvolved 6 deg2 (3° × 2°) maps provided by BLAST in 2005 at 250, 350, and 500 μm were used to perform a preliminary characterization of the clump population previously investigated in the infrared, radio, and molecular maps. Interferometric CORNISH data at 4.8 GHz have also been used to characterize the Ultracompact H ii regions (UCHiiRs) within the main clumps. By means of the BLAST maps, we have produced an initial census of the submillimeter structures that will be observed by Herschel, several of which are known Infrared Dark Clouds. Our spectral energy distributions of the main clumps in the field, located at ∼7 kpc, reveal an active population with temperatures of T∼ 35–40 K and masses of ∼103M for a dust emissivity index β = 1.5. The clump evolutionary stages range from evolved sources, with extended H ii regions and prominent IR stellar population, to massive young stellar objects, prior to the formation of an UCHiiR. The CORNISH data have revealed the details of the stellar content and structure of the UCHiiRs. In most cases, the ionizing stars corresponding to the brightest radio detections are capable of accounting for the clump bolometric luminosity, in most cases powered by embedded OB stellar clusters.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

The study of the earliest stages of star formation is intrinsically connected to the ability to probe the dusty regions within which stars are born. The short wavelength radiation from deeply embedded young stellar objects is absorbed by the dust and reaches the observer as re-processed far infrared (FIR) and submillimeter emission. The low optical depth at these longer wavelengths makes submillimeter studies the best tool to investigate the clumps and inner cores (∼1–0.1 pc; Williams et al. 2000) that produce new stars, especially the "coldest" structures prior to the onset of stellar activity.

The 2 m Balloon-borne Large-Aperture Submillimeter Telescope (BLAST; Pascale et al. 2008), with two consecutive and successful science missions in 2005 and 2006, mapped the sky simultaneously at the crucial wavelength bands of 250, 350, and 500 μm. The Galactic (e.g., Netterfield et al. 2009) and extragalactic (e.g., Devlin et al. 2009) data provided by the BLAST mission demonstrate the exciting science to be continued with the Herschel Space Observatory15 (hereafter Herschel; e.g., Elia et al. 2010) and its SPIRE instrument (Griffin et al. 2010; 18'' at 250 μm).

In this study, we present a preliminary analysis of the 6 deg2 field of a massive star-forming region in Aquila (ℓ =  45°, b = 0°), observed by BLAST during its 2005 flight. This work therefore complements the analysis of the other BLAST05 Galactic plane star-forming regions, including Vulpecula (Chapin et al. 2008) and Cygnus X (Roy et al. 2010b).

Aquila is well known for its prominent massive star formation (e.g., Rathborne et al. 2004), and is dominated by two main molecular complexes. These are host to numerous Ultracompact H ii regions (UCHiiRs; e.g., Wood & Churchwell 1989; [WC89]), maser emission, and outflows from young and highly embedded OB stellar clusters. The stellar activity of the Galactic Ring Survey Molecular Cloud (GRSMC) 45.60+0.30 has been extensively studied in the infrared and radio (e.g., Kraemer et al. 2003; Vig et al. 2006). However, although submillimeter data have also been used to complement the detailed study of this cloud (e.g., Vig et al. 2006), the overall population and properties of this entire region in Aquila are yet to be fully characterized at submillimeter wavelengths.

Previous BLAST studies of regions like Vulpecula (Chapin et al. 2008) and Vela (Netterfield et al. 2009) have provided robust samples of submillimeter sources in a range of evolutionary stages. By fitting modified blackbody spectral energy distributions (SEDs) to the BLAST measurements, and with the aid of ancillary data at infrared (IR) wavelengths for the Vulpecula sources (including IRAS, Spitzer MIPS, and the Midcourse Space Experiment, MSX), these studies constrained dust temperatures, masses, and bolometric luminosities for their respective samples. Their results aimed to characterize and investigate crucial properties of the submillimeter population, including their star formation stage, the (core) mass function, and core lifetimes (Vela).

In this paper, we provide the first general characterization of the submillimeter clump population of Aquila, which will serve as the initial census and preliminary analysis of those structures that will be observed in more detail by Herschel during its surveys of the massive star-forming regions in the Galactic plane (Hi-GAL: The Herschel infrared Galactic Plane Survey; Molinari et al. 2010). Like in previous BLAST studies, we provide estimates of physical parameters such as mass, dust temperature, and luminosity, although in this work this analysis is constrained to the most prominent clumps in the field, lying at the common distance of ∼7 kpc (Section 2.2). We have complemented our submillimeter analysis with a detailed study of the UCHiiRs and OB clusters within these main clumps with the aid of 4.8 GHz radio interferometry data provided by the Co-Ordinated Radio "N" Infrared Survey for High-mass star formation (CORNISH; Purcell et al. 2008). In Section 2, we describe the submillimeter and radio data used in this work, as well as the techniques required to produce our final results. We then analyze the major clouds and clumps in the region (Section 3), whose properties we discuss in Section 4. In Section 5, we give a brief summary of our results and the prospects for future studies of this massive star-forming region.

2. OBSERVATIONS AND DATA ANALYSIS

2.1. Submillimeter Data Products and Source Detection

BLAST observed Aquila during its first science flight in 2005 June, mapping a total of 6 deg2 (3°× 2°) at 250, 350, and 500 μm toward GRSMC 45.60+0.30. For the following analysis, we note that the three BLAST maps were made such that the average of these maps is zero (Patanchon et al. 2008). During the reduction and mapmaking process, the SANEPIC (Signal and Noise Estimation Procedure) algorithm removed low-frequency noise (mainly sky signal). Large-scale signals in the map were recovered, but the DC level was set to zero by application of a high-pass filter to the time-ordered data. The reduction process as well as the subsequent procedures to correct the unexpected optical degradation that occurred during the flight (Truch et al. 2008) have been described in detail in Chapin et al. (2008) and Roy et al. (2010a). The effective resolution is ∼40'', 50'', and 60'' at 250, 350, and 500 μm, respectively. The 250 and 350 μm maps were later convolved in order to produce BLAST maps with a common resolution of 60'' for photometry analysis.

As can be seen in Figure 1, the final products still suffered from artifacts (ripples) surrounding the brightest objects and the map edges arising from the deconvolution process. These introduce false detections when running the Interactive Data Language (idl16) routines designed to detect the main submillimeter peaks in the field, and which we therefore have to eliminate manually after visual inspection. For source identification we used the idl-implemented DAOPHOT "find" routine and a signal-to-noise ratio (S/N) detection limit of 3 at 250 μm, the band with higher S/N for a typical star-forming clump. We kept a total of 66 sources between ℓ = 44fdg6 and 46fdg8 (to avoid map edges) which also had a match at 350 μm within 20'' with an S/N >3. We also kept six additional sources (identified with * in Table 1) with a lower S/N at 250 and 350 μm (combined S/N between ∼2and3.5, visually prominent, and some with counterparts in the IRAS Point Source Catalog; PSC). Their lower S/N can be explained by their location in the local neighborhood of (clustered) bright sources, regions with difficult background estimation, and significantly above average noise. This source catalog (Table 1) is robust to ∼50 Jy at 350 μm (∼140 Jy at 250 μm). For comparison, the cirrus noise (3 σcirrus) at 250, 350, and 500 μm is estimated to be ∼20, 10, and 5 Jy for each map, respectively (Roy et al. 2010a). In the present work, we did not require completeness to faint flux densities, however, choosing to focus our analysis only on the most prominent submillimeter peaks of our sample.

Figure 1.

Figure 1. BLAST submillimeter map in Aquila combining three wavebands at 250, 350, and 500 μm, with color coding blue, green, and red, respectively.

Standard image High-resolution image

Table 1. Submillimeter Clumps in the BLAST Aquila Map

Source BLAST Namea b IRAS Nameb IRDCc Cloudd Clumpe D (kpc)f Assoc.g
A0 J191131+102635 44.502 0.347 19091+1021          
A1 J191124+102845 44.522 0.387 19090+1023          
A2 J191432+100817 44.577 −0.454   P3726 G044.34−00.21 c12 4.12 Y
A3 J191135+103138 44.585 0.370            
A4 J191148+103504 44.660 0.350 19094+1029          
A5 J191334+102308 44.686 −0.129 19112+1018   G044.49−00.16 c8 7.55  
            G044.69−00.11 c1 3.7 Y
A6 J191514+101158 44.711 −0.577     G045.09−00.51 c4 6.0  
            G044.74−00.56 c1,c7 ... Y
A7 J191537+101015 44.730 −0.675     G044.74−00.56 c2 ... Y
A8 J191514+101518 44.760 −0.551     G044.74−00.56 c1,c8 ... Y
A9 J191457+102119 44.816 −0.442 19125+1015   G045.09−00.51 c6 6.0  
A10 J191611+101914 44.927 −0.728            
A11 J191606+102731 45.039 −0.644            
A12 J191321+105049 45.070 0.132 19110+1045   G045.14+00.14 c1 4.53 Y
A13 J191538+103419 45.087 −0.492 19132+1029   G045.09−00.51 c1 6.0 Y
A14 J191528+103636 45.100 −0.435            
A15 J191421+104622 45.118 −0.119 19119+1041          
A16 J191327+105338 45.123 0.132 19111+1048   G045.14+00.14 c1 4.53 Y
A17 J191641+102952 45.140 −0.752 19143+1024          
A18 J191602+103743 45.183 −0.553     G045.39−00.76 c6 7.18 Y
A19 J191539+104116 45.190 −0.439 19132+1035   G045.09−00.51 c2 6.0 Y
A20 J191706+103032 45.198 −0.838            
A21 J191633+104658 45.378 −0.592 19142+1041   G045.74−00.26 c2 7.28 Y
            G045.39−00.76 c3,c10 7.18 ?
            G045.84−00.56 c4 3.62  
A22 J191706+104322 45.387 −0.738     G045.39−00.76 c1,c9,c11,c13 7.18 Y
A23 J191701+104447 45.399 −0.710     G045.39−00.76 c7,c11,c13 7.18 Y
A24 J191620+105035 45.407 −0.518 19139+1045   G045.84−00.56 c4 3.62  
A25* J191451+110244 45.416 −0.100 19125+1057          
A26* J191457+110324 45.437 −0.116            
A27 J191532+105933 45.447 −0.273            
A28 J191419+110902 45.449 0.063 19120+1103   G045.49+00.04 c1,c5,c7,c8 7.45 Y
            G045.14+00.14 c5 4.53 ?
            G045.24+00.19 c1 ...  
A28ah J191414+110809 45.425 0.077            
A29 J191425+110921 45.464 0.046 19120+1103i P3744 G045.49+00.04 c1,c4,c5,c7 7.45 Y
              c8    
            G045.14+00.14 c2 4.53 ?
A30 J191407+111222 45.475 0.133 19117+1107   G045.49+00.04 c2,c3,c5-c8 7.45 Y
            G045.24+00.19 c1 ...  
A30a J191409+111245 45.485 0.128            
A31 J191350+111548 45.493 0.222     G045.24+00.19 c1 ...  
A32 J191536+110259 45.506 −0.262            
A33 J191558+110113 45.521 −0.354            
A34 J191608+110059 45.536 −0.391 19138+1055          
A35 J191450+111127 45.543 −0.029     G045.49+00.04 c4 7.45 Y
A36 J191446+111159 45.544 −0.012 19124+1106   G045.49+00.04 c4 7.45 Y
A37 J191512+111007 45.565 −0.119     G045.49+00.04 c4 7.45  
            G045.64-00.01 c2 11.12 Y
A38 J191459+111541 45.623 −0.030     G045.49+00.04 c4 7.45  
            G045.64-00.01 c1,c3,c4 11.12 Y
A39 J191458+111639 45.635 −0.018     G045.49+00.04 c4 7.45  
            G045.64-00.01 c1,c4 11.12 Y
A40 J191358+112532 45.651 0.270     G045.49+00.04 c9 7.45  
            G045.89-00.36 c4 6.72  
            G045.64+00.29 c1,c5,c6 1.88  
A41 J191709+110129 45.660 −0.608     G045.39-00.76 c2 7.18 Y
            G045.94-00.56 c5 4.78 ?
A42* J191553+111347 45.696 −0.237     G046.19-00.56 c4 3.75  
            G045.74-00.26 c1,c10,c11 7.28 Y
              c13-c15,c23,c24    
            G045.89-00.36 c2 6.72 ?
            G045.84-00.56 c3 3.62  
A43 J191705+110509 45.707 −0.566            
A44 J191609+111423 45.736 −0.291 19137+1108   G046.19-00.56 c4 3.75  
            G045.74-00.26 c1,c8,c21 7.28 Y
            G045.89-00.36 c2 6.72  
            G046.14-00.21 c2 1.8  
            G045.84-00.56 c3 3.62  
A45 J191604+111725 45.772 −0.251     G046.19-00.56 c4 3.75  
            G045.74-00.26 c1,c12,c21 7.28 Y
            G045.89-00.36 c2 6.72 ?
            G045.84-00.56 c3 3.62  
            G046.09+00.24 c2 ...  
A46 J191630+111606 45.802 −0.355 19141+1110   G045.74-00.26 c16,c19 7.28 Y
            G045.89-00.36 c2 6.72 ?
            G046.09+00.24 c2 ...  
A47 J191618+111909 45.823 −0.286 19139+1113   G046.19-00.56 c4 3.75  
            G045.74-00.26 c1 7.28 ?
            G045.89-00.36 c2 6.72 Y
            G045.84-00.56 c10 3.62  
            G046.09+00.24 c2 ...  
A48 J191713+111605 45.884 −0.511   P3758 G045.89-00.36 c1 6.72 Y
A49 J191447+113705 45.915 0.181 19124+1131   G045.89-00.36 c6 6.72  
            G045.64+00.29 c1 1.88  
A50 J191655+112153 45.935 −0.401 19145+1116   G045.89-00.36 c1,c2 6.72 Y
            G045.94-00.56 c2 4.78 ?
A51 J191514+113632 45.959 0.078     G045.89-00.36 c6 6.72  
            G046.14-00.21 c6 1.8  
A52 J191456+114614 46.068 0.218     G045.89-00.36 c6 6.72  
            G045.64+00.29 c1 1.88  
            G046.09+00.24 c1 ... Y
A53 J191450+114805 46.084 0.254 19125+1142   G045.89-00.36 c6 6.72  
            G045.64+00.29 c1 1.88  
            G046.09+00.24 c1 ... Y
A53a J191448+114910 46.096 0.270            
A54 J191413+115330 46.094 0.429     G046.14+00.39 c1 4.25 Y
            G045.64+00.29 c1 1.88  
A55 J191423+115342 46.115 0.396 19120+1148   G046.14+00.39 c1 4.25 Y
            G045.64+00.29 c1 1.88  
A56 J191617+114240 46.169 −0.102 19139+1137   G045.89-00.36 c6 6.72  
            G046.14-00.21 c1 1.8  
A57 J191750+113102 46.174 −0.527   P3762 G046.19-00.56 c1 3.75 Y
A58 J191400+120052 46.177 0.535 19116+1155          
A59 J191402+120157 46.197 0.536            
A59a   46.212 0.551            
A60 J191819+113119 46.233 −0.628     G046.24-00.66 c1 4.07 Y
A61 J191851+112750 46.244 −0.773            
A62 J191702+114640 46.314 −0.234 19146+1141   G046.19-00.56 c2 3.75 ?
            G046.34-00.21 c2,c3,c5,c7 4.12 Y
              c9-c11    
A63* J191703+115012 46.367 −0.208     G046.19-00.56 c2 3.75 ?
            G046.34-00.21 c2,c5,c9,c12 4.12 Y
              c13    
A64* J191709+114956 46.374 −0.231     G046.19-00.56 c2 3.75 ?
            G046.34-00.21 c2,c5,c12 4.12 Y
A65 J191655+115407 46.410 −0.150 19146+1148 P3773 G046.19-00.56 c2 3.75 Y
            G046.34-00.21 c2,c5,c12 4.12  
A66 J191715+115225 46.423 −0.235 19149+1147 P3776 G046.19-00.56 c2 3.75 ?
            G046.34-00.21 c1,c2,c5,c6 4.12 Y
              c12    
A67* J191711+115548 46.465 −0.194   P3778 G046.19-00.56 c2 3.75 ?
            G046.34-00.21 c1,c5,c6,c12 4.12 Y
            G046.79+00.04 c3 3.6  
A68 J191631+120108 46.468 −0.009            
A69 J191630+120314 46.496 0.013            
A70 J191849+115057 46.581 −0.586 19164+1145   G046.19-00.56 c3 3.75 ?
            G046.89-00.36 c1 8.05 Y
A71 J191537+121816 46.618 0.319 19132+1212   G046.59+00.34 c1 0.65 Y
            G046.59+00.19 c2 0.8 ?

Notes. *Detection S/N <3. aSource names include the prefix "BLAST" as part of the standard BLAST nomenclature. bClosest IRAS match within 30'' in the Point Source Catalog (PSC). cClosest IRDC match within 30'' in the catalog from Peretto & Fuller (2009). dGRS name of the extended clouds from Rathborne et al. (2009) which contain the molecular clump counterpart(s) to the BLAST source. A blank implies the absence of a clump counterpart in the catalog. eClump(s) from Rathborne et al. (2009) containing the BLAST peak within its/their ℓ and b FWHM. fDistance from Roman-Duval et al. (2009). gCloud association: (Y) - BLAST clump is associated with the cloud. Cloud has a molecular clump whose velocity FWHM contains the strongest molecular peak of spectrum, measured at the peak of the submillimeter emission. Molecular clump also has peak velocity closest to that measured for the BLAST source. (?) - Cloud membership unclear. The main molecular emission in the BLAST spectrum is within the velocity FWHM of this molecular cloud/clump, but there is still another cloud with a clump with peak velocity closer to that of the strongest line in the BLAST spectrum. Contributions from different clouds are possible. A blank implies that the velocity of main BLAST molecular emission is not compatible with that assigned to the cloud/clump, despite the spatial agreement and apparent location of the submillimeter peak within the spatial FWHM of clump(s) associated with that cloud. hSource designations ending in "a" refer to additional Gaussians required during the fit of the main clump. iPosition beyond 30'', but visually within the IRAS emission.

Download table as:  ASCIITypeset images: 1 2 3

2.2. Distance Estimation

In order to investigate the velocities and distances toward the main clumps in the field we used the 13CO molecular spectra from the BU-FCRAO Galactic Ring Survey (GRS; Jackson et al. 2006). The data sets cover the region between 18°<ℓ <55fdg7 with a velocity resolution of 0.2 km s−1. The spatial resolution is 46'' with an angular sampling of 22''. Due to the distribution of the clumps throughout the field, this spectral analysis is essential in order to estimate and justify the use of a suitable common distance to the most active regions in Aquila. A comparison of BLAST with the molecular emission is shown in Figure 2. To aid qualitative intercomparison, BLAST contours are presented in this and a number of later figures, from images in which an approximate DC level has been restored (300, 160, and 50 MJy sr−1 at 250, 350, and 500 μm, respectively).

Figure 2.

Figure 2. Gray scale: 13CO molecular emission at ∼60 ± 10 km s−1 with labels marking position and names (GRSMC) of clouds in this range from Rathborne et al. (2009). Contours are BLAST 500 μm surface densities. White: contours from 10% to 90% of map peak value of 2153 MJy sr−1 in 5% steps. Black: 5% contour level (see the text).

Standard image High-resolution image

For the spectral analysis we first extracted and fitted the 13CO spectra at the position of each submillimeter peak. Molecular counterparts from the catalogs of Roman-Duval et al. (2009) and Rathborne et al. (2009) were assigned to each BLAST peak if the peak was located within the spatial FWHM of the molecular clump. The velocity we obtained from the spectral fitting was then compared with the clump velocities in the catalog. This allowed us to distinguish the appropriate molecular counterpart to the BLAST clump in those cases where more than one molecular counterpart (but with different local standard of rest (LSR) radial velocity) was found for a particular BLAST clump (e.g., Table 1). Our analysis and line fitting of the spectra is in good agreement with the significant dispersion of velocities (e.g., Figure 3) existing in the recent cloud/clump catalogs of Rathborne et al. (2009) and Roman-Duval et al. (2009). The main activity of the field is located in the main broad high-velocity peak in Figure 3. The smaller peak at lower velocities could arise from clumps located either much closer or much farther away (∼1 or ∼11 kpc; Roman-Duval et al. 2009) than the complexes forming the high-velocity peak (which cluster around the tangent point of ∼6 kpc; e.g., Simon et al. 2001).

Figure 3.

Figure 3. Histogram of 13CO radial velocities of 68 BLAST detections with molecular clump counterpart (BLAST peak lying within spatial FWHM of clump) in catalogs of Roman-Duval et al. (2009) and Rathborne et al. (2009). The main stellar activity has velocities within the broad (high velocity) peak.

Standard image High-resolution image

Using the Clemens (1985) rotation curve of the Milky Way, Roman-Duval et al. (2009) estimated distances to GRS molecular clouds (e.g., Rathborne et al. 2009) with the aid of H i self-absorption (HISA) and continuum absorption, both in data from the Very Large Array Galactic Plane Survey (VGPS; Stil et al. 2006). There are four clouds in their sample within the defined BLAST map limits (Section 2.1) with an assigned distance between ∼7 ± 0.5 kpc (LSR velocity within the second peak in Figure 3: v ∼60 km s−1): G045.49+00.04, G045.74−00.26, G045.89−00.36, and G045.39−00.76. However, one of the major clouds in the BLAST field, G045.14+00.14 (vLSR ≈ 60 km s−1), has instead been assigned its near distance of ∼4.5 kpc.

Taking into account a maximum standard error of 10% for the distances (Roman-Duval et al. 2009), we estimate that a common distance of ∼7 kpc would be a suitable approximation for the analysis of all the clouds with velocities of ∼60 km s−1 and located in this distance range. Furthermore, this estimate is in good agreement with the distances derived by Pandian et al. (2009), who also resolved the kinematic distance ambiguities using H i self-absorption with data from the VGPS. These authors calculated, by means of methanol masers and a flat rotation curve, a range of distances for the clumps within G045.49+00.04 between 6.9 and 7.4 kpc, which again agrees well with the estimate from Roman-Duval et al. (2009). However, they assigned the far distance of 7.4 kpc to the clump in G045.14+00.14, which disagrees with Roman-Duval et al. (2009). This coincides with previous studies (e.g., Simon et al. 2001) that placed this cloud at the same kinematic distance as its neighbor G045.49+00.04. In this paper, we adopt the far distance estimate of Pandian et al. (2009) for G045.14+00.14, and we restrict our main analysis to four clouds (G045.49+00.04, G045.74−00.26, G045.89−00.36, and G045.14+00.14), which comprise the most active complexes of the field. Our main BLAST detections corresponding to G045.39−00.76 are associated with a bubble (Churchwell et al. 2006), separate from the main clouds of the field, and we therefore omit these sources from the main analysis. A common distance of ∼7 kpc (1° ∼0.1 kpc scale) is assumed in the present study (in contrast with those used by previous authors of ∼6 and ∼8 kpc; e.g., Simon et al. 2001; Kraemer et al. 2003; Paron et al. 2009).

2.3. Ancillary Data

To create our SEDs, we used the BLAST fluxes estimated using our own Gaussian-fitting routine (Section 2.4). Whenever possible, we also included data from the IRAS Galaxy Atlas (IGA; Cao et al. 1997: HIRES; Aumann et al. 1990) at 60 and 100 μm and BOLOCAM (Galactic Plane Survey) data at 1.1 mm (Rosolowsky et al. 2010). IRAS flux densities were estimated individually using our Gaussian-fitting routine, while for the convolved (60'') BOLOCAM images we used aperture photometry (single sources: an aperture centered on the BLAST source with 1farcm8 radius; an estimate for the background was obtained from an annulus formed by the source aperture and another aperture with radius 1.3 times larger than that of the source). The measured BOLOCAM fluxes were found to be within the error range of those provided by the GATOR service,17 and the latter values were therefore used to constrain the SEDs at longer wavelengths if a counterpart to the BLAST clump was found within 30'' (we found 39 sources in Table 1 with BOLOCAM counterparts, including all those with measured SEDs, with the exception of A22). We note that IRAS was found to be saturated for the brightest clumps, for which we adopted the upper limits estimated by Kraemer et al. (2003). Most of our measured IRAS fluxes were also found to be higher than those from the IRAS PSC. This effect, combined with the lower resolution of IRAS and the photometry of sometimes blended sources, implies that these values should also be generally considered upper limits in our fits.

To investigate the stellar content and the morphology of the emission at shorter wavelengths we used the Infrared Array Camera (IRAC) on Spitzer. Of the four available IRAC bands (3.6, 4.5, 5.8, and 8 μm), we focused on the longest wavelength in our morphology analysis. For our study of the stellar population we used data from the Galactic Legacy Infrared Midplane Survey Extraordinaire (GLIMPSE; Benjamin et al. 2003) through the online gator query system. The GLIMPSE I online data sets, covering about 220 deg2 of the Galactic plane (ℓ = 10°–65°) in all four IRAC bands, include a highly reliable (GR) catalog and a more complete, but less reliable one (GC). Both were used in the present work.

2.4. Photometry: Spectral Energy Distributions and Parameter Estimation

Similar to other Galactic BLAST studies (e.g., Chapin et al. 2008), our SED analysis was carried out by estimating flux densities with our own idl spatial Gaussian-fitting routine, and subsequently fitting them with a modified blackbody of the form

Equation (1)

where β is the dust emissivity index, (ν/ν0)β is the emissivity factor normalized at a fixed frequency ν0, and Bν(T) is the Planck function. The amplitude A can be expressed as

Equation (2)

where D is the distance to the source, Mc is the total clump mass, κ is the dust mass absorption coefficient, evaluated at ν0 = c/250 μm, and R is the gas-to-dust ratio. A two-dimensional Gaussian is fitted to the three bands simultaneously with the amplitudes, centroid, orientation angle, and major and minor FWHM as free parameters. The fit is carried out using a nonlinear least-squares minimization idl routine (based on mpfit; Markwardt 2009), which also allowed for simultaneous multiple Gaussian fitting in crowded regions of the field. The background was estimated by fitting a fourth-order polynomial to the regions around the clump.

Following Chapin et al. (2008), estimates and uncertainties for the dust temperatures, masses, and luminosities were obtained through a Monte Carlo analysis performed on our SED models, keeping the dust emissivities and R as fixed parameters (β = 1.5, R = 100, κ = 10 cm2 g−1). Examples of SEDs are shown in Figures 4 and 5. A flux density consistency check between all convolved and deconvolved maps shows flux conservation to better than 5%.

Figure 4.

Figure 4. SED of A16 (Table 1). Symbols are BLAST and ancillary data (IRAS; BOLOCAM). Upper and lower curves are the 1σ envelope of the dust emission model (Chapin et al. 2008).

Standard image High-resolution image
Figure 5.

Figure 5. Same as Figure 4, but for A50 (Table 1).

Standard image High-resolution image

The bolometric luminosities derived from the fits for each evolved clump were subsequently used to estimate equivalent single zero-age-main-sequence star (ZAMS) spectral types using the tables of Panagia (1973). This luminosity value was used as a general measure of the energy output of each submillimeter detection. We note that the BLAST sources at a distance of ∼7 kpc are parsec-sized, and therefore "clump-like" objects (e.g., Bergin & Tafalla 2007). We expect them to be powered by several embedded stars, with the luminosity dominated by that of the most massive source(s) within the clump. In addition, we also include the equivalent luminosity class V star spectral types from the more recent stellar calibration provided by Martins et al. (2005), which takes into account the effects of line-blanketing and winds in the non-LTE atmosphere models of O stars. For this case, we used the calibration based on the observational Teff scale; these authors note that the theoretical and observational scales are in good agreement for early-type dwarfs and supergiants.

2.5. 4.8 GHz Radio Interferometry Data and Analysis

For our radio analysis we used the CORNISH interferometry data sets, a 4.8 GHz Very Large Array (VLA) survey covering the northern GLIMPSE region with a resolution of ∼1''. Our analysis of these data was carried out using the same procedure and Gaussian-fitting codes used for the submillimeter. The 4.8 GHz fluxes from our fitting were then used to estimate the total number of ionizing photons per unit time for each radio source (e.g., Panagia 1973; Martins et al. 2005):

Equation (3)

where ν0 is the frequency at the Lyman edge.

These values were also used in conjunction with the tables from Panagia (1973) and Martins et al. (2005) to estimate the equivalent single-star radio spectral types for each detection. The positions and parameters derived using these data have been included in Tables 24. The target S/N of our sample is >3. As with the BLAST sources, we have also kept a few radio detections with lower S/N. These were identified and required by our Gaussian-fitting routine, and are almost exclusively secondary peaks neighboring the brightest sources (clusters: A16–C24, A28–C15), where the S/N is greatly deteriorated (>1.5). This effect has already been discussed in Purcell et al. (2008), and these sources will therefore require further examination as soon as additional/better data become available. Extra detections with low S/N have been identified (*) in Table 2. Nomenclature of radio sources is formed by the name of the BLAST clump, followed by that of the CORNISH detection (A#−C#).

Table 2. Coordinates of CORNISH Sources in the BLAST Aquila Field

BLAST CORNISH R.A. Decl. BLAST CORNISH R.A. Decl.
    (J2000) (J2000)     (J2000) (J2000)
    (h m s) (° ' '')     (h m s) (° ' '')
A12 A12−C0 19 13 19.063 +10 51 26.521 A28 A28−C0 19 14 20.710 +11 09 13.540
  A12−C1* 19 13 21.617 +10 50 56.756   A28−C1 19 14 20.714 +11 08 57.156
  A12−C2 19 13 21.886 +10 50 49.704   A28−C2 19 14 20.784 +11 09 17.737
  A12−C3 19 13 22.130 +10 50 52.163   A28−C3 19 14 20.918 +11 09 15.926
  A12−C4 19 13 23.834 +10 51 42.624   A28−C4 19 14 20.957 +11 09 01.728
  A12−C5 19 13 23.954 +10 51 41.497   A28−C5 19 14 21.002 +11 09 16.578
  A12−C6 19 13 25.548 +10 51 07.567   A28−C6 19 14 21.002 +11 09 20.246
  A12−C7 19 13 25.584 +10 51 07.535   A28−C7 19 14 21.120 +11 09 08.550
  A12−C8 19 13 25.601 +10 51 09.320   A28−C8 19 14 21.190 +11 09 07.268
  A12−C9 19 13 25.697 +10 50 57.977   A28−C9 19 14 21.226 +11 09 02.110
A16 A16−C0 19 13 24.283 +10 53 28.525   A28−C10 19 14 21.250 +11 09 14.915
  A16−C1* 19 13 24.374 +10 53 27.247   A28−C11 19 14 21.276 +11 09 11.534
  A16−C2 19 13 24.446 +10 53 25.919   A28−C12 19 14 21.338 +11 09 22.162
  A16−C3 19 13 24.994 +10 53 13.895   A28−C13 19 14 21.362 +11 09 10.084
  A16−C4 19 13 25.116 +10 54 16.153   A28−C14 19 14 21.386 +11 09 11.242
  A16−C5 19 13 25.205 +10 53 14.932   A28−C15 19 14 21.408 +11 09 13.309
  A16−C6 19 13 27.055 +10 53 19.946   A28−C16 19 14 21.578 +11 09 04.824
  A16−C7 19 13 27.218 +10 53 22.794   A28−C17 19 14 21.617 +11 09 15.264
  A16−C8 19 13 27.331 +10 53 46.489   A28−C18 19 14 21.720 +11 08 55.424
  A16−C9* 19 13 27.516 +10 53 46.558   A28−C19 19 14 21.854 +11 08 54.413
  A16−C10 19 13 27.545 +10 53 30.196   A28−C20 19 14 21.862 +11 09 30.956
  A16−C11* 19 13 27.583 +10 53 46.237   A28−C21 19 14 21.893 +11 09 17.791
  A16−C12* 19 13 27.583 +10 53 49.636   A28−C22 19 14 22.022 +11 09 45.385
  A16−C13 19 13 27.586 +10 53 27.920   A28−C23 19 14 22.099 +11 09 45.972
  A16−C14 19 13 27.602 +10 53 25.483   A28−C24 19 14 22.272 +11 09 37.760
  A16−C15 19 13 27.703 +10 53 29.249   A28−C25 19 14 22.481 +11 09 31.090
  A16−C16 19 13 27.710 +10 53 29.807   A28−C26 19 14 22.615 +11 09 43.690
  A16−C17 19 13 27.763 +10 53 38.990   A28−C27 19 14 22.776 +11 08 57.970
  A16−C18 19 13 27.773 +10 53 34.897 A29 A29−C0 19 14 25.690 +11 09 25.330
  A16−C19 19 13 27.797 +10 53 32.194   A29−C1 19 14 26.222 +11 08 35.488
  A16−C20 19 13 27.866 +10 53 30.181   A29−C2 19 14 26.942 +11 09 14.897
  A16−C21 19 13 27.888 +10 53 40.402   A29−C3 19 14 27.079 +11 09 12.456
  A16−C22 19 13 27.912 +10 53 33.659   A29−C4 19 14 27.559 +11 09 56.617
  A16−C23 19 13 27.929 +10 53 27.013 A30 A30−C0 19 14 07.330 +11 12 41.742
  A16−C24 19 13 27.965 +10 53 35.675   A30−C1* 19 14 08.227 +11 12 22.684
  A16−C25 19 13 28.008 +10 53 45.150   A30−C2 19 14 08.232 +11 12 36.929
  A16−C26 19 13 28.022 +10 53 28.806   A30−C3 19 14 08.602 +11 12 26.327
  A16−C27 19 13 28.030 +10 53 24.598   A30−C4* 19 14 08.830 +11 12 26.586
  A16−C28 19 13 28.034 +10 53 42.302   A30−C5 19 14 09.046 +11 12 25.661
  A16−C29 19 13 28.128 +10 53 31.708   A30−C6 19 14 09.094 +11 12 22.075
  A16−C30 19 13 28.212 +10 53 37.547 A30a A30a−C0 19 14 08.604 +11 13 19.175
  A16−C31 19 13 28.222 +10 53 34.260 A47 A47−C0 19 16 17.645 +11 18 59.670
  A16−C32 19 13 28.226 +10 53 35.905   A47−C1 19 16 17.746 +11 19 14.509
  A16−C33 19 13 28.380 +10 53 32.298   A47−C2 19 16 19.330 +11 19 14.992
  A16−C34 19 13 29.842 +10 53 34.318 A50 A50−C0 19 16 55.730 +11 21 48.355
  A16−C35 19 13 31.090 +10 53 45.031        

Note. *Detection S/N <3.

Download table as:  ASCIITypeset image

Table 3. Spectral Types and IR Counterparts for the BLAST Peaks A16 and A12 and the Radio Sources Detected in the CORNISH Data Within ∼1' of the Submillimeter Peaks

Source IR+submillimetera IR+submillimetera IRb TFT1c TFT2c
  (ZAMS) (Class V)   (6 cm) (3.6 cm)
  log Q0d 4.8 GHz:ZAMS/Class V IRe 3.6 cme GLIMPSEf
  (s−1)        
A16 O6−O5.5 O6−O4 O5 O6 O6
A16−C0 46.86+0.01−0.01 B0.5−B0/−      
A16−C1 46.96+0.16−0.26 B0.5−B0/−      
A16−C2 46.97+0.01−0.01 B0/−      
A16−C3 47.36+0.00−0.00 B0/−      
A16−C4 46.70+0.01−0.01 B0.5/−      
A16−C5 46.83+0.01−0.01 B0.5−B0/−      
A16−C6 47.06+0.17−0.27 B0.5−B0/−     G045.1163+00.1327
A16−C7 47.03+0.01−0.01 B0/−     G045.1163+00.1327
A16−C8 46.68+0.04−0.05 B0.5/−      
A16−C9 46.71+0.05−0.06 B0.5/−      
A16−C10 47.58+0.01−0.01 B0/−      
A16−C11 46.51+0.01−0.01 B0.5/−      
A16−C12 46.73+0.01−0.01 B0.5/−      
A16−C13 47.33+0.12−0.16 B0/−      
A16−C14 47.73+0.08−0.09 O9.5/−      
A16−C15 47.54+0.05−0.05 B0/−      
A16−C16 47.63+0.16−0.25 B0−O9.5/−      
A16−C17 48.30+0.03−0.03 O8.5−O8/O8.5     G045.1221+00.1322
A16−C18g 47.74+0.44−47.74 <O8.5/O8.5     G045.1221+00.1322
A16−C19 47.89+0.01−0.01 O9.5/O9.5     G045.1221+00.1322
A16−C20 47.79+0.01−0.01 O9.5/−      
A16−C21 48.25+0.20−0.38 O9.5−O7.5/O8     G045.1221+00.1322
A16−C22g 48.16+0.01−0.01 O8.5/O9−O8.5     G045.1221+00.1322
A16−C23 46.88+0.02−0.02 B0.5−B0/−      
A16−C24 49.14+0.00−0.00 O6/O5.5 O5 O6 G045.1221+00.1322
A16−C25 47.53+0.07−0.08 B0/−      
A16−C26 47.90+0.05−0.06 O9.5/O9.5 B0.5 B0  
A16−C27 46.59+0.04−0.05 B0.5/−      
A16−C28 47.98+0.02−0.02 O9.5−O9/O9.5−O9      
A16−C29 47.87+0.01−0.01 O9.5/−      
A16−C30 48.09+0.16−0.26 O9.5−O8.5/O8.5      
A16−C31 47.93+0.23−0.52 B0−O8.5/O9      
A16−C32 48.29+0.01−0.01 O8/O8.5      
A16−C33 47.58+0.29−1.18 B0.5 − 09.5/−      
A16−C34 47.06+0.01−0.01 B0/−      
A16−C35 47.65+0.00−0.00 O9.5/−     G045.1302+00.1223
A12 O6 O5 O6 O8.5h O8h
A12−C0 47.32+0.04−0.05 B0/−     G045.0729+00.1484
A12−C1 46.58+0.00−0.00 B0.5/−      
A12−C2 46.49+0.00−0.00 B0.5/− B0.5 B0−O9.5 G045.0712+00.1321
A12−C3 47.88+0.00−0.00 O9.5/− O6.5 O9 G045.0712+00.1321
A12−C4 47.53+0.07−0.08 B0−O9.5/−      
A12−C5 47.99+0.00−0.00 O9/O9      
A12−C6 46.47+0.02−0.02 B0.5/−      
A12−C7g 48.15+0.2−0.39 O9.5−O9/<O8.5      
A12−C8 47.03+0.09−0.11 B0/−      
A12−C9 47.32+0.03−0.03 B0/−     G045.0788+00.1191

Notes. aZAMS and class V spectral types from Panagia (1973) and Martins et al. (2005) for D = 7 kpc. If the luminosity is not within the range of the calibration scale, no estimate (−) is provided. bZAMS spectral type derived by Kraemer et al. (2003) from the total source flux in the IR (IRAS) for D = 6 kpc. cZAMS spectral type of the counterpart of the brightest CORNISH source within BLAST clump in Testi et al. (1999) for D = 7 kpc. dTotal number of ionizing photons per unit time for a distance of 7 kpc and Te ∼ 8000 K, assuming optically thin emission. eZAMS spectral type derived by Kraemer et al. (2003) for D = 6 kpc. Values are provided only for sources lying closer than 5'' to the CORNISH source. fClosest IRAC GLIMPSE I Complete Catalog match within ∼5'' of the CORNISH peak. gExtra source required during Gaussian fitting. hEstimate includes total flux of C0 and C1, as both are unresolved at 6 cm (Testi et al. 1999).

Download table as:  ASCIITypeset image

Table 4. Spectral Types and IR Counterparts for the BLAST Peaks A28, A29, A30, and A30a and the Radio Sources Detected in the CORNISH Data Within ∼1' of the Submillimeter Peaks

Source IR+submillimeter IR+submillimeter IR TFT1 TFT2
  (ZAMS) (Class V)   (6 cm) (3.6 cm)
  log Q0 4.8 GHz:ZAMS/Class V IR 3.6 cm GLIMPSE
  (s−1)        
AA29 O6.5−O5.5 O7−O4.5      
A28     O5.5   O6.5a−O5.5b
Complex 1 49.25+0.00−0.00 O5.5/O5     G045.4509+00.0570c
A28−C1 47.35+0.04−0.04 B0/−      
A28−C18 46.79+0.15−0.23 B0.5−B0/−      
A28−C19 47.19+0.01−0.01 B0/−      
A28−C20 47.11+0.11−0.16 B0/−      
A28−C22 46.38+0.21−0.42 B0.5/−      
A28−C23 47.13+0.15−0.23 B0.5−B0/−      
A28−C24 46.69+0.03−0.04 B0.5/−      
A28−C25 46.62+0.02−0.02 B0.5/−      
A28−C26 46.91+0.02−0.02 B0.5−B0/−     G045.4640+00.0574
A28−C27 46.69+0.02−0.02 B0.5/−      
A29         O9.5
A29−C0 47.84+0.00−0.00 O9.5/−     G045.4661+00.0457(GR)
          G045.4658+00.0452 (GC)
A29−C1 46.44+0.02−0.02 B0.5−B0/−     G045.4551+00.0359
A29−C2 47.17+0.01−0.01 B0/−      
A29−C3 47.28+0.00−0.01 B0/−      
A29−C4 47.29+0.02−0.02 B0/−      
AA30 B0−O6.5 <O6.5      
A30     O6.5−O6   O6.5
Complex 2 48.47+0.00−0.00 O7.5/O8 O6.5 O7  
A30−C0 46.78+0.04−0.04 B0.5−B0/−     G045.4790+00.1365
A30−C1 47.09+0.06−0.06 B0/−      
A30−C2 47.27+0.11−0.16 B0/−      
A30a          
A30a−C0 46.49+0.03−0.03 B0.5/−     G045.4909+00.1376(GR)
          G045.4904+00.1367(GC)
A46 B0.5−B0        
A47 O8−O7.5 O9−O8      
A47−C0 46.93+0.26−0.74 B0.5−B0/−     G045.8202−00.2866
A47−C1 46.61+0.06−0.08 B0.5/−     G045.8232−00.2842 (GR)
A47−C2 46.73+0.08−0.09 B0.5−B0/−     G045.8266−00.2892 (GR)
A50 O9.5−O8.5 <O9.5      
A50−C0 46.97+0.09−0.11 B0.5−B0/−     G045.9344−00.4016

Notes. aCompact source only, ∼6'' in size, and therefore smaller than our aperture used for photometry. bCompact source including the extended emission. cThe closest CORNISH counterpart to this source is C54 (Table 2).

Download table as:  ASCIITypeset image

Considering the uncertainties in our measurements and the above procedure for spectral type estimation, assuming optically thin emission, we estimate an uncertainty for these measurements of ∼0.5 in spectral type (for a given calibration scale of the stellar parameters). This uncertainty is expected to be larger than any spectral type error that could arise when comparing our estimates (for a distance of ∼7 kpc) and those obtained by previous authors with a distance of ∼6 kpc (Tables 3 and 4).

3. DISCUSSION OF INDIVIDUAL SOURCES

In the following sections, we apply the results of our analysis to describe the main submillimeter regions (G045.49+00.04, G045.74−00.26, G045.89−00.36, and G045.14+00.14), adopting the designations from Rathborne et al (2009; Figure 2).

3.1. GRSMC G045.14+00.14

3.1.1. G45.12

GRSMC 45.122+0.132 (G45.12; IRAS 19111+1048; A16 in the present work) contains a cometary UCHiiR (WC89) classified as a massive young stellar object (MYSO) by Chan et al. (1996). It shows extended submillimeter emission toward the north (upper left of Figure 6), observed in the MIR as a bright rimmed H ii region. This structure is likely produced by the activity of the young massive population embedded in the innermost regions of the clump, as detected by Vig et al. (2006) using the Giant Metrewave Radio Telescope (GMRT) at 1200 and 610 MHz. The 19 embedded OB stars observed by these authors are located within and surrounding the peak of our BLAST clump, mainly within the H ii region. Their brightest source agrees well with our submillimeter peak (≲5''), and although it remained unresolved in their data sets, they suspected it to be a cluster responsible for the extended morphology in the radio. A possible cluster powering this region has also been suggested through [Ne ii] observations by Takahashi et al. (2000), and later by Zhu et al. (2008), who detected three additional peaks (G45.12N, G45.12SE, and G45.12SW) nearby the main UCHiiR from WC89.

Figure 6.

Figure 6. Gray-scale IRAC 8 μm image of IRAS 19111+1048 (A16) and IRAS 19110+1045 (A12) with BLAST 500 μm contours superimposed. Contours are from 10% to 90% of map peak value 2153 MJy sr−1 in 10% steps.

Standard image High-resolution image

Our analysis of the 4.8 GHz CORNISH data confirms such suspicions, with the detection of 28 radio sources in a region ∼0.8 pc in radius centered on (and including) the central brightest member (A16−C24; Table 3 and Figure 7). These objects are believed to be forming the unresolved source observed by Vig et al. (2006; S14 in their sample). None have a radio ZAMS spectral type later than B0 (with the exception of A16−C18, which was not initially detected by our extraction routine).

Figure 7.

Figure 7. Gray-scale CORNISH 4.8 GHz image of central cluster (S14 in Vig et al. 2006) in IRAS 19111+1048 (A16) with numbered emission peaks (Table 3). White contours are from 2.5% to 12.5% of map peak value 0.7 Jy beam−1 in 1% steps. Gray contours are from 15% to 85% in 10% steps.

Standard image High-resolution image

We find that 22 of these radio sources lie within ∼0.5 pc of A16−C24. This object has a bolometric luminosity equivalent to an O6 ZAMS star, which is in good agreement with previous results for a distance of ∼6 kpc (e.g., Kraemer et al. 2003). The position of A16−C24 also agrees with the UCHiiR of WC89 and Testi et al. (1999). The two Two Micron All Sky Survey (2MASS) IR sources detected by Vig et al. (2006) within their S14 source (IR4 and IR5) agree well with A16−C24 and A16−C26 and the two peaks (KJK1 and KJK2) detected by Kraemer et al. (2003). The G45.12SE peak from Zhu et al. (2008) traces the region of A16−C26; G45.12SW is formed by the conglomeration of radio sources toward A16−C13, and G45.12N is centered on the arc-shaped radio emission around A16−C21. The estimated spectral types of all radio sources within ∼1' of the BLAST coordinates have been included in Table 3 for an electron temperature similar to that estimated by Vig et al. (2006; ∼8000 K, for their source containing the cluster S14).

Taking into account the 3σ positional error of our CORNISH source A16−C24 (∼0farcs5) and that expected for the OH maser detection of Baart & Cohen (1985), ∼10 mas, the separation of ∼3'' between the two may be physically significant. Indeed, our projected separation of ∼0.1 pc positions the maser source at the edge of the UCHiiR, and agrees well with the displacement measured by Baart & Cohen (1985), of the order of ∼0.03 pc (with respect to the coordinates from Matthews et al. 1977). This distance, as pointed out by Baart & Cohen (1985), is well within the upper limit for the radius at which OH masers are produced in a shell-like front surrounding a massive young star. This suggests that these may be dense fragments of the expanding shell. The blueshift for the maser clusters detected by these same authors (∼6 km s−1) is also similar to that observed for the dense gas traced with ammonia by Hofner et al. (1999), whose absorption peak coordinates are also offset from our radio peak in the same direction as the maser source. This supports the real physical origin of the observed displacement.

Although further investigation is required, our analysis of the CORNISH data and the evidence in the mid-infrared (MIR) and submillimeter suggest that A16−C24 may be the initial trigger and power source of G45.12. The massive outflow detected by Hunter et al. (1997) appears to follow the overall direction of the distribution of CORNISH sources well, which could indicate that A16−C24 has induced the formation of its lower mass companions. In addition, the distribution of CORNISH sources appears to have a slightly curved, bow-shock-like morphology pointing toward the submillimeter comma-like structure extending downward in Figure 8. This structure also coincides with the strongest emission in 13CO and CS. Considering its central position, A16−C24 appears as a suitable candidate for the main triggering/compressing source.

Figure 8.

Figure 8. Gray-scale IRAC 8 μm image of IRAS 19111+1048 (A16) with 250 μm BLAST contours overlaid. Cross marks the position of submillimeter peak and circles the positions of main CORNISH detections (central cluster). Contours are from 5% to 85% of the map peak value of 31,500 MJy sr−1 in 10% steps.

Standard image High-resolution image

3.1.2. G45.07

Just like its neighbor, IRAS 19110+1045 (A12; GRSMC 45.073+0.129 [G45.07]) also shows submillimeter, IR, and radio structure signaling the presence of massive star formation, albeit significantly less pronounced than in G45.12. Only nine possible CORNISH detections are found within ∼1' of the BLAST position, with radio ZAMS/class V spectral types ranging between B0.5 and O9 for an electron temperature of 8000 K (Table 3). The lower number of stars and their later apparent radio spectral type is consistent with the more compact size of the resulting ionized region compared to G45.12.

Two CORNISH sources, A12−C2 and A12−C3 (Figure 9), match (≲3farcs5) those detected by Vig et al. (2006) and the two embedded IR sources (KJK1, KJK3) from Kraemer et al. (2003). Despite the positional accuracy of our CORNISH measurements (<1farcs5 difference between our two coordinates and those given by Testi et al. 1999 and WC89), we find no indication of the third IR source (KJK2) detected by Kraemer et al. (2003). This source is located ∼2'' from A12−C3 and should have been identifiable at the resolution limit of the CORNISH data.

Figure 9.

Figure 9. Gray-scale CORNISH image of three central sources likely powering IRAS 19110+1045 (A12). Contours are from 2% to 9% of map peak value 0.7 Jy beam−1 in 1% steps.

Standard image High-resolution image

As with G45.12, G45.07 also contains a UCHiiR (spherical/unresolved UCHiiR; WC89) with a CO outflow (Hunter et al. 1997) and masers. Our brightest radio source lies at <1farcs5 from the coordinates given by WC89. The masers show a more significant displacement (∼2'') and appear to extend along the overall outflow direction. In addition to hydroxyl masers, G45.07 also shows several water maser sources (e.g., Hofner & Churchwell 1996) and 6.7 GHz methanol emission (Pandian et al. 2007), with distances ranging between ∼1''–3farcs5 and ∼10farcs5 from our brightest CORNISH radio peak, respectively.

Most water masers align well with our radio source A12−C3 and the overall outflow direction. This, together with the clear separation from our radio peak, again supports the scenario where they are energized by the bipolar outflow (e.g., Hofner & Churchwell 1996; De Buizer et al. 2005). The presence of these water masers also implies a less advanced evolutionary stage than that of its neighbor G45.12 (in agreement with previous studies: e.g., Vig et al. 2006). This is supported by the simpler overall morphology, the lower stellar content, and the much earlier spectral type inferred from the IR/submillimeter for A12−C3 and for the total clump with respect to those derived from the radio, suggesting a much more deeply embedded stage with significant radio optical depth.

While there appears to be extended MIR emission in the GLIMPSE images, consistent with the outflow orientation, there is also considerable emission in an almost orthogonal direction to its main axis (western upper right corner in Figure 10). The peak of the 21 cm radio continuum is displaced toward this region. The overall morphology appears to be similar to the extended H ii in G45.12, albeit in a smaller scale, which could be due to its younger age and/or the lower stellar content in the innermost central region.

Figure 10.

Figure 10. Gray-scale IRAC 8 μm image of IRAS 19110+1045 (A12). Symbols and contours like Figure 8.

Standard image High-resolution image

3.1.3. Submillimeter Analysis

Our analysis of the SEDs yields total system masses for G45.12 and G45.07 of ∼3550 M and ∼1400 M, respectively, for a distance of 7 kpc (Table 6). This differs from the estimates given by Hunter et al. (1997), who obtained (with temperature, dust emissivity index, and optical depth as free parameters) masses ∼7 and ∼4 times larger than our estimated masses, respectively (when scaled to their distance of 8.3 kpc). Vig et al. (2006) obtained, through SED fitting and radiative transfer modeling (assuming a spherically symmetric homogeneous cloud of gas and dust and a distance of 6 kpc), masses of 5000 M (∼20 M of dust and a larger gas-to-dust ratio of 250) and 13,000 M (for a dust mass of 30 M and gas-to-dust ratio of 450) for G45.12 and G45.07, respectively. Scaling to the same distance, we obtain Md ∼ 26 and Md ∼ 10.5 M. While our dust mass estimate for G45.12 is in good agreement with their estimate from the SCUBA maps, the mass for G45.07 is in better agreement with the mass these authors estimated from their 130 μm map (obtained with the TIFR 1 m balloon-borne telescope; Ghosh et al. 1988) of ∼13 M.

We find temperatures for G45.12 and G45.07 of ∼40 K, and clump luminosities equivalent to an ∼O6 ZAMS spectral type star (∼O5 class V star). This is in good agreement with those estimated from the IR (e.g., Kraemer et al. 2003) for the brightest sources within the clumps. This suggests that A16−C24 and A12−C3 are the main sources responsible for the dust heating in G45.12 and G45.07, respectively. We note, however, that the significantly later spectral type derived from the radio for A12−C3 appears to suggest a much younger and highly embedded stage than A16−C24.

3.2. GRSMC G045.49+00.04

This cloud complex is dominated by the other two main submillimeter peaks in the BLAST region, IRAS 19120+1103 (GRSMC 45.453+0.060; [G45.45]; a MYSO candidate in the Red MSX Source (RMS) Survey; e.g., Mottram et al. 2010) and IRAS 19117+1107 (GRSMC 45.478+0.131; [G45.48]). Both sources show significant extended submillimeter emission (Figure 11), and required additional Gaussians to model the emission local to the main submillimeter peaks. At IRAS resolution, the extended emission of both sources was not resolved, and so our BLAST flux densities were obtained with aperture photometry combining A28 with A29 (AA29; Figure 12, comprising G45.47+0.05 and G45.45) and A30 with A30a (AA30).

Figure 11.

Figure 11. Gray-scale IRAC 8 μm image of IRAS 19120+1103 and IRAS 19117+1107. Contours like Figure 6.

Standard image High-resolution image
Figure 12.

Figure 12. Gray-scale IRAC 8 μm image of IRAS 19120+1103 (A28−A29) with 250 μm BLAST contours overlaid. Crosses mark the position of main GR sources. Contours are from 5% to 55% of the map peak value of 31,500 MJy sr−1 in 10% steps. Other symbols like Figure 8.

Standard image High-resolution image

The presence of UCHiiRs was detected by WC89 and Testi et al. (1999) at 3.6 cm. Chan et al. (1996) classified AA29 as a MYSO. The activity and complexity of this region is particularly prominent in the IRAC images, as can be observed in Figures 1113.

Figure 13.

Figure 13. Gray-scale IRAC 8 μm image of most central regions of A28 and A29 in Figure 12 with CORNISH contours. Contours are from 10% to 90% of the map peak value of 0.07 Jy beam−1 in 10% steps.

Standard image High-resolution image

3.2.1. G45.47+0.05

A29 (G45.47+0.05) is the most prominent of all submillimeter peaks in this complex, and yet lacks significant emission in the MIR (e.g., Figures 12 and 13). The submillimeter peak is located in the neighborhood of 14 GR sources (five of which are associated with bright IR sources; Figure 12). All these sources are located within ∼30'' of the BLAST peak. The appearance in the MIR may suggest that the submillimeter emission is tracing a particularly dense part of the molecular clump. This material could be confining the expansion of the photodissociation region (PDR) as seen in IRAC, which also appears to be in a very young and embedded stage of star formation. Pillar-like structures seen in extinction and extending from A29 into G45.45 (Figure 12), and the main five IR detections embedded within the BLAST source (as suggested by their flux increasing with increasing wavelength) support this scenario.

The brightest CORNISH source in the neighborhood of the submillimeter peak, A16−C0 (Figure 14), is a known UCHiiR (WC89; Testi et al. 1999 center agreement <1''). The flux from our radio analysis is equivalent to an O9.5 radio ZAMS star at 7 kpc, which agrees with the spectral type we derived using the 6 cm flux from Urquhart et al. (2009). This source is also one of the most complex of the clump, and contains H2O (e.g., Forster & Caswell 1989), OH (e.g., Forster & Caswell 1989; Argon et al. 2000), and methanol (e.g., Menten 1991) masers. The most prominent OH sources are located closer (<0farcs5) to the radio peak than the methanol and H2O masers (∼1''), and they all lie preferentially toward the northern (upper left direction; Figure 14) rim of the radio emission and toward the extended MIR region.

Figure 14.

Figure 14. Gray-scale IRAC 8 μm image of region around the main UCHiiR (A29−C0; C0 in the image) within A29 (Figure 12) with CORNISH contours superimposed. Contours are from 10% to 60% of map peak value of 0.07 Jy beam−1 in 10% steps. The figure shows closest IR counterparts retrieved from the GLIMPSE I complete catalog (GC), labeled as G1, G2, and G3.

Standard image High-resolution image

In the IRAC maps, the GR source (G045.4661+00.0457; G3 in Figure 14) is at the peak of the IR emission, at the opposite corner of the MIR elongation to the CORNISH source (at ∼3''). Other GR sources are also present nearby the radio peak. Such astrometry differences between the radio and IR detections have been previously observed (De Buizer et al. 2005), which confirms this to be a true physical separation. The MIR elongation may be tracing the outflow morphology detected by Hunter et al. (1997) and the NH3 emission (e.g., Hofner et al. 1999). Most maser sources are indeed located between the MIR peak (G3) and the radio peak, along the elongation, with a few OH masers located between the radio source and G2 (Figure 14). At the northern extreme edge of the elongation and closer to the main MIR peak, there are more detections of methanol maser signatures (Kurtz et al. 2004), as well as SiO and NH3 emission peaks (∼10''–15'' from the radio peak; Hofner et al. 1999). The presence of shocked gas along the northern MIR emission, the outflow, and the location of the UCHiiR opposite the main MIR peak may suggest the scenario where the outflow/jet originating at the UCHiiR interacts with its local environment, producing the observed structure in the MIR. An alternative interpretation claims this elongation may arise from accretion (De Buizer et al. 2005) or collapse (e.g., Cesaroni et al. 1992; Hofner et al. 1999) onto the UCHiiR. This could explain the emission and redshifted absorption observed in NH3. More recent studies have not confirmed the presence of infall signatures (e.g., Wilner et al. 1996; Klaassen & Wilson 2007). Further analysis of the innermost regions is needed to solve this controversy.

The maser emission (especially methanol and H2O) suggests that this source is one of the youngest in this complex. It has also been classified as an Extended Green Object (EGO) and a likely MYSO outflow candidate (Cyganowski et al. 2008). This supports the HCO analysis of Wilner et al. (1996), who suggested that this source is in the early stages of forming an OB star cluster.

3.2.2. G45.45

A28 is close to the most active region in the cloud as observed in the MIR (G45.45+0.06 [G45.45]), although from its position and extended morphology the submillimeter emission may also be tracing the compressed and extended molecular gas. G45.45 (also a MYSO candidate in the RMS Survey; e.g., Mottram et al. 2010) has been observed to lie in the border of a larger and fainter H ii region ∼3' in size (G45L; Paron et al. 2009). Although this structure is also traced by the BLAST images, it has not been included in our photometry/radio analysis.

The IRAC images show a complex morphology, with the strongest emission in the MIR and 21 cm coincident with the main nebula-like radio structure observed in the CORNISH data (e.g., Figures 12 and 13). In Figure 13, the MIR emission appears to follow a horseshoe shape around a cavity-like structure (∼1 pc in size at 7 kpc), within which no stars are detected in the IRAC images. Although this structure can be observed at 1.3 cm in the analysis of Mooney et al. (1995) at ∼12'' resolution, the details and morphology of the dust in this region remain unresolved in our BLAST maps.

In the radio, the CORNISH data reveal a similarly complex substructure for G45.45 (Figure 15). There are 18 radio peaks detected within and surrounding the radio "nebula," coincident with the brightest area in the GLIMPSE images. Within this extended radio emission, called "The Orion nebula's younger brother" by Feldt et al. (1998), numerous emission peaks are noticeable. While we cannot easily distinguish or fit these peaks, we detect at least eight within an estimated size of 0.35 × 0.2 pc at 7 kpc (Figure 16). This supports previous studies that classify this source as a cluster of H ii regions (e.g., Garay et al. 1993; Testi et al. 1999). Aperture photometry on the whole structure (including extended emission) yields a radio flux density comparable to that of an ∼O5.5 ZAMS spectral type (∼O5 class V) star, which is in agreement with previous estimates (e.g., Garay et al. 1993). Using the 6 cm flux from Urquhart et al. (2009) for a source major axis of ∼8'', we obtain an estimate of ∼O6.5. This is comparable to the results we obtain using the flux estimates from Testi et al. (1999; Table 4).

Figure 15.

Figure 15. Detailed gray-scale CORNISH image of central UCHiiR near A28 (part of IRAS 19120+1103; Figures 12 and 13), with numbered emission peaks (Table 2). Contours like Figure 13.

Standard image High-resolution image
Figure 16.

Figure 16. More detailed gray-scale CORNISH image of innermost regions of the central UCHiiR in Figure 15 with numbered emission peaks (Table 2). Letters are NIR detections from Feldt et al. (1998) falling within VLA map (Table 5). Roman numbers are other NIR detections (Feldt et al. 1998). Contours like Figure 13.

Standard image High-resolution image

Regarding the stellar content, Feldt et al. (1998) detected up to 14 near-infrared (NIR) sources (from a to o) possibly associated with this complex. They also detected additional MIR sources outside the main radio nebula (Figure 16). The roman numbers indicate other point sources in their NIR images outside their VLA map. We confirm their observations that the majority of these objects lie on the sharp and distorted upper edge of the radio emission (Figure 16). This may support the scenario proposed by these authors where a central cluster of OB stars, deeply embedded within the radio complex, produces a front that has triggered further star formation. In their analysis, they claim that IR sources l, m, n, o (Figure 16) may be the source of the UCHiiR. Although these objects do lie close to two bright radio peaks within the cluster (A28−C13, A28−C11; Table 5), their source b is the closest to the brightest peak in the radio (A28−C15 in Figure 16). This CORNISH source lies at less than ∼1'' from both the WC89 coordinates and source b.

Table 5. 4.8 GHz Radio (CORNISH) Sources with IR Counterparts (Feldt et al. 1998) Within 1''

IR Radio
b A28−C15
d A28−C10
f A28−C11
l A28−C11
m A28−C11
n A28−C11
o A28−C13
MIR1 A28−C10

Download table as:  ASCIITypeset image

At a resolution of 0farcs4, WC89 defined this source as a cometary UCHiiR. Indeed, two parallel extensions appear to emerge from the central source A28−C15. These structures form a tail-like morphology extending toward the right in Figure 16, within which A28−C13, A28−C14, and A28−C11 are detected. In total, we find that eight of these IR sources have CORNISH counterparts within ∼1''. Sometimes one radio peak is found between several in the NIR (Table 5 and Figure 16).

Blum & McGregor (2008) classified the main radio structure as a group of MYSOs with late and early B stars surrounding the central source, probably similar to the central cluster powering G45.12 within A16. The offset between the NIR and radio peaks could be explained if the radio detections are dense clumps ionized by the surrounding massive stars (Feldt et al. 1998). Alternatively, we do not discard the possibility that the NIR stars may have been produced by the central and embedded star/cluster, whose expanding H ii region could have triggered the formation of these sources by compression of the dense material traced in the submillimeter. This last scenario is supported by the overall morphology in the CORNISH images; the observed emission supports A28−C15, the brightest radio source, as being responsible for the origin of l, m, n, o, rather than the latter sources being the actual "ground zero" for the triggering mechanism (as suggested by Feldt et al. 1998). Further investigation at high resolution is required to corroborate the triggering hypothesis and to fully distinguish between the different cases.

The candidate counterpart to MIR1 (A28−C10) is clearly visible in the CORNISH data, albeit relatively weak compared to the most central sources. This weak emission is in contrast with the strong IR luminosity of this source. In combination with the nearby OH maser (Argon et al. 2000), this may again suggest that at least some structures are actual stars still embedded within their natal cocoons (e.g., Kraemer et al. 2003).

The remaining possible radio detections within ∼1' and not lying within the main complex have been included in Table 4. Due to the complexity of the radio structures around this complex, we do not discard the presence of additional sources in this region (e.g., Figure 15).

3.2.3. G45.48

AA30 (G45.48; north of AA29 in Figure 11) also shows complex structure in the submillimeter and IR. Our two submillimeter peaks lie at opposite sides of the main emission visible in the IRAC images. As in previous cases, they may be tracing relatively cold, dense gas that is confining the expansion of the H ii region.

WC89 classified this UCHiiR as irregular, but suggested that the elongated form could be due to an imaging artifact. Testi et al. (1999) identified no compact sources, only extended emission. The elongation of this source, as well as several embedded peaks, are evident at 4.8 GHz (Figure 17), and the structure is identified in the MIR images as the equatorial ridge of the bipolar-like MIR structure (Figure 18).

Figure 17.

Figure 17. Gray-scale CORNISH image of elongated (and brightest) radio emission in A30−A30a (IRAS 19117+1107; Complex 2 in Table 4); with numbered emission peaks (Table 2). Contours are from 10% to 25% of map peak value of 0.07 Jy beam−1 in 3% steps.

Standard image High-resolution image
Figure 18.

Figure 18. Gray-scale IRAC 8 μm image of the H ii region within IRAS 19117+1107. Numbers as in Figure 17. Crosses mark the position of submillimeter peaks A30 and A30a.

Standard image High-resolution image

Just like in G45.45, the MIR and main radio emission appear to wrap around a cavity-shaped structure, within which no IRAC sources are visible. The BLAST emission includes this "cavity," but any such structure would be unresolved in our images. The adjacent radio/IR ridge could be compatible with illumination from the south, although the main source powering the region could in fact be located more toward the north, closer to source A30−C0 (Figure 18), the northern radio source detected by Garay et al. (1993), and the diffuse emission in the IR (KJK1) of Kraemer et al. (2003). Only one GC was found within ∼5'' of A30−C0 (G045.4790+00.1365). KJK1 does appear to be coincident with an IRAC emission peak, even though there is no counterpart in the GLIMPSE I catalogs. The irregular radio structures, and the lack of a significant radio counterpart for KJK1, could imply the presence of an embedded star/cluster around this position, which also coincides with an opacity peak and is estimated to be much colder than the southern counterpart (Kraemer et al. 2003).

The multi-peaked radio morphology of the elongated radio structure (coincident with the southernmost source detected by Garay et al. 1993 and the IR source KJK2 from Kraemer et al. 2003) clearly indicates the presence of embedded sources. Our Gaussian routine detects at least four main peaks in a filamentary structure about ∼0.2 pc in length (Figure 17). Aperture photometry on the entire filament yields a total radio flux equivalent to an O7.5 radio ZAMS (O8 class V) star. This is significantly later than the IR estimate (Table 4), perhaps an indication of self-absorption and a young embedded stage of evolution, despite our earlier spectral type estimate obtained using the 3.6 cm flux from Testi et al. (1999). This structure is also coincident with methanol masers (Menten 1991), recognized as a signature of highly embedded stages of massive star formation. In addition, we observe (weaker) non-Gaussian, filamentary structures to the right of the main structure in Figure 17, also likely associated with methanol emission. A prominent cluster of methanol, OH, and water masers has been detected about ∼15''–20'' southwest from A30−C0 and A30−C1. The presence of a GC/GR counterpart (G045.4725+00.1335) for this source and the lack of an equivalent CORNISH detection again argues in favor of the young stage of the stellar population in this region.

3.2.4. BLAST A35 and A36

Away from the main complex, A35 and A36 (Figure 19) also show some signatures of stellar activity, as suggested by their morphology.

Figure 19.

Figure 19. Gray-scale IRAC 8 μm image showing cometary shape of A35+A36 with BLAST 500 μm contours superimposed. Contours are from 8% to 12% of map peak value of 2153 MJy sr−1 in 2% steps.

Standard image High-resolution image

Although our BLAST peaks generally overlap well with the emission observed in IRAC, IRAS 100 μm shows poor correlation. Only at 60 μm do we observe a reasonable match with A36, the only one of these two sources with a PSC IRAS match within ∼30''. This IR source, IRAS 19124+1106, was estimated to be at a distance of 6.5 kpc, and has been suggested to be one of the impact sites of the jets emanating from the microquasar and superluminal source GRS 1915+105 (Kaiser et al. 2004). The other impact site was identified as our BLAST source A19, a candidate MYSO from the RMS Survey (Mottram et al. 2010). GRS 1915+105 is believed to be a binary system consisting of an early-type K giant, accreting via a Roche lobe overflow onto a black hole (e.g., Greiner 2001). The MIR structure of A36, if associated with this IRAS source, does not resemble what one would expect to see if caused by a jet originating at the position of GRS 1915+105. Furthermore, the new distance estimate for GRS 1915+105 of 11 kpc (Zdziarski et al. 2005) is now much larger than the kinematic distance to this IR/submillimeter emission.

Our analysis of the CORNISH images reveals that, with the exception of some weak structures near the emission peak of A36 (likely tracing the ionized region as observed in the VGPS 21 cm maps), there are no major radio sources near these objects. This may suggest very early or very late/isolated low-mass star-forming activity.

3.2.5. Submillimeter Analysis

Estimated BLAST parameters for the main objects in this complex have been included in Table 6. The temperatures are warm, characteristic of active star-forming clumps with ongoing stellar activity. The total mass of AA29 is ∼4800 M for a temperature of ∼35 K, which contrasts with the mass estimate of 26,000 M obtained by Mooney et al. (1995) for a distance of 8.1 kpc and a much lower fixed dust temperature of 20 K. Scaling to this distance we obtain a total mass of ∼6300 M, and ∼34,000 M for a fixed temperature of 20 K and free β (best-fit β = 2.2). This shows the importance of actually determining the temperature using the BLAST observations.

Table 6. Masses and Temperatures of the Main BLAST Sources at ∼7 kpc

Source T Mc Lbol Single-star Spectral Type
  (K) (102 M) (103 L)  
A12 42.5 ± 0.5 14.2 ± 1.0 282.8 ± 0.6 O6a/O5b
A16 39.9 ± 4.6 35.4 ± 8.0 340.9 ± 131.7 O6−O5.5/O6−O4
A22c 16.2 ± 4.6 65.7 ± 47.6 7.5 ± 6.6 <B0.5/−
A28−A29 36.7 ± 4.8 47.8 ± 10.3 262.3 ± 121.6 Combined due to IRAS unresolved (and saturated) emission
        O6.5−O5.5/O7−O4.5
A30−A30a 27.4 ± 5.8 33.4 ± 9.1 89.1 ± 67.2 Combined due to IRAS unresolved (and saturated) emission
        B0−O6.5/<O6.5
A35 33.9 ± 10.9 4.3 ± 1.5 62.0 ± 144.9 ∼O8/∼O9
A36 41.6 ± 1.4 1.7 ± 0.3 30.1 ± 0.4 B0/−
A46 32.5 ± 1.6 4.5 ± 0.7 19.1 ± 2.7 B0.5−B0/−
A47 37.9 ± 2.5 7.9 ± 1.3 76.0 ± 14.7 O8−O7.5/O9−O8
A50 35.0 ± 2.0 6.9 ± 1.2 43.6 ± 6.6 O9.5−O8.5/<O9.5

Notes. aSpectral type from Lbol for a ZAMS star (Panagia 1973). bSpectral type from Lbol for a luminosity class V star (Martins et al. 2005). cNo PSC counterpart within 30''. No BOLOCAM data available. Possibly associated (visually) with an IRAS source, but this source covers nearby BLAST sources as well and is therefore ignored in the fits. The rising SED in the submillimeter and the lack of data at shorter wavelengths result in large error ranges.

Download table as:  ASCIITypeset image

Our estimated IR+submillimeter single-star spectral type for AA29 of ∼O5.5 is in excellent agreement with the radio spectral type of the radio nebula within G45.45 (Complex 1 in Table 4), the most luminous system within the BLAST emission formed by A28 and A29 (AA29).

The IR+submillimeter spectral type of A30+A30a is in agreement with the IR spectral type estimate of Complex 2 (the elongated radio emission; Table 4), but the CORNISH radio spectral type of this substructure is considerably later.

Although the radio flux estimates from Testi et al. (1999) at 3.6 cm and the uncertainties in the SED bolometric luminosities could be compatible with a relatively optically thin state for this source, the CORNISH measurements suggest a more embedded (and younger) stage for AA30 than for AA29. Just like the clumps within GRSMC G045.14+00.14, the BLAST sources in GRSMC G045.49+00.04 appear to be powered by the central OB young stellar clusters, which could also be responsible for the additional star formation detected within their parent clumps.

3.3. GRSMC G045.74−00.26 and GRSMC G045.89−00.36

This complex is dominated by three main (sub)millimeter sources (Figure 20) within an extended and largely filamentary molecular region.

Figure 20.

Figure 20. Gray-scale IRAC 8 μm image of IRAS 19141+1110 (A46), IRAS 19139+1113 (A47), and IRAS 19145+1116 (A50), with BLAST 500 μm contours superimposed. Contours are from 10% to 20% of map peak value of 2153 MJy sr−1 in 5% steps. White crosses are CORNISH detections.

Standard image High-resolution image

3.3.1. BLAST A50

In the submillimeter, A50 (IRAS 19145+1116) is largely structureless and symmetric, possibly with some weak emission surrounding the central peak. In the MIR, this H ii region reveals a clumpy structure and numerous MIR peaks clustered within the main submillimeter emission. There are seven GC sources found at less than ∼20'' from our BLAST peak. Among the closest of these we find a string of three sources, G045.9352−00.4017 (the central one), G045.9362−00.4017 (the brightest), and G045.9344−00.4016, which constitute the most prominent structure in the GLIMPSE images.

Overall, the extended structure in the MIR and at 21 cm is rather circular, although emission arcs are visible within the region in the MIR. We detect only one weak 4.8 GHz source, found at less than ∼3'' from the GLIMPSE source G045.9344−00.4016, with a flux corresponding to a B0.5−B0 ZAMS star. The IRAS clump was reported to satisfy the color criteria of WC89 for UCHiiRs (e.g., Bronfman et al. 1996), and was reported by Codella et al. (1995) to be possibly associated with water masers. This would argue against the more evolved state of the system indicated by its "disturbed" morphology, although from the coordinates provided by these authors we suspect that this maser emission could also be associated with a later generation of star formation within an older clump. The lack of a CORNISH or IR counterpart for this maser may also be indicating that the structure is due to a low-mass star.

3.3.2. BLAST A46

While relatively structureless in the submillimeter, the GLIMPSE images of A46 (IRAS 19141+1110) show a bow-shock cometary-like structure, with the front oriented toward the south (lower right of Figure 20). Some possibly associated sources toward the south and southeast of the main structure are also detected. Despite the presence of numerous GC/GR sources, the bow shock itself may be produced by a highly embedded source, closest to G045.8062−00.3518 (Figure 20).

Even though A46 also has IRAS, MSX, and IRAC counterparts, there appears to be negligible 21 cm emission and no NVSS, 4.8 GHz emission, H2O or OH masers, or any other tracers that may be indicative of ongoing star formation. An exception is the methanol maser emission ∼25'' eastward (lower left) from the BLAST peak detected by Pandian et al. (2007). This is the 51 − 60A+ line of methanol at 6.7 GHz, the strongest of class II methanol masers, seen exclusively near massive stars. The authors suggest a very early stage of evolution, in which methanol activity has just turned on, before the onset of additional signs of star formation. Class II maser sources are radiatively pumped but are believed to be quenched by the formation of an UCHiiR, which is consistent with the absence of 4.8 GHz CORNISH sources and the "young" age of the system (e.g., Minier et al. 2005; Ellingsen et al. 2007). This scenario would also be in agreement with Cyganowski et al. (2008), who classified this source as a "possible" MYSO. The observed IRAC structures and methanol emission are all within the prominent BLAST submillimeter emission, which also suggests the presence of local high density molecular material.

3.3.3. BLAST A47

A47 (IRAS 19139+1113) is a prominent MYSO (e.g., Chan et al. 1996) showing clear signatures of active star formation. A prominent carved structure (∼2 pc in size at 7 kpc) is observed toward the north (Figure 20). This H ii region is the most extended feature of this cloud, and is detectable in the IR, submillimeter, millimeter, and 21 cm images. The submillimeter and radio peaks appear to trace a region of diffuse emission within this structure, close to three bright IRAC sources (G045.8263−00.2820, G045.8212−00.2846, and G045.8291−00.2858). The location of these IR objects, close to the center of curvature of some of the ionization fronts forming the edge of the ionized region in the molecular cloud, suggests that the combined contribution from these objects results in the sinusoidal boundary of the overall region.

We find three likely sources at 4.8 GHz within ∼1' of our BLAST source, all with GC counterparts within ∼5'', and indicative of ongoing massive star formation. The extended morphology and the advanced state of stellar activity within the main H ii would support the "older" age suggested by the lack of maser emission (e.g., van der Walt et al. 1995; Szymczak et al. 2000).

3.3.4. Radio and Submillimeter Analysis

From their measured fluxes, we estimate ZAMS spectral types for the radio peaks in these clumps of later than ∼B0 stars. Our SED fits show that, overall, none of the host clumps have an equivalent single-star ZAMS spectral type earlier than an ∼O8 star. In the case of A47, it is possible that additional sources are contributing to the total luminosity, in addition to those traced in the radio.

The detection of 4.8 GHz emission confirms previous studies that classified A47 and A50 as UCHiiRs from the WC89 color criteria (Bronfman et al. 1996). Although the radio peaks are not coincident with the brightest IRAC sources, they clearly indicate that massive stars are indeed being formed in this cloud, albeit at a lower level than that observed in the two main clouds described above.

4. AQUILA IN PERSPECTIVE: THE STAR-FORMING ACTIVITY IN THE BLAST 6 deg2 MAP

4.1. The Submillimeter Population and Structural Properties of the Aquila Field

This region of Aquila consists principally of three main active complexes at ∼ 7 kpc (G045.49+00.04, G045.14+00.14, and the neighboring region of G045.74−00.26, and G045.89−00.36), as well as a quiescent cloud (GRSMC 45.60+0.30) at 1.8 kpc (e.g., Simon et al. 2001). There also exists (at various distances) additional dispersed (clump) star-forming activity throughout the field, two cataloged (and two candidate) supernova remnants (Green 2009), and a microquasar (GRS 1915+105).

In the present work, we have used our BLAST submillimeter data to identify and characterize the global parameters of the clumps at ∼7 ± 0.5 kpc, in which the molecular emission (and star formation activity) is the strongest. These clumps have been further characterized in the MIR and the radio with IRAC GLIMPSE and 21 cm continuum maps. The stellar content as well as the molecular structure have been investigated with CORNISH interferometry data and GRS 13CO/CS emission.

At ∼60 ± 10 km s−1 the 13CO maps reveal a dispersed and filamentary molecular structure covering the size of the BLAST map. Two of the main complexes are located at the central "knots" of an overall heart-shaped-like morphology (Figure 2). From the catalog of Rathborne et al. (2009) and Roman-Duval et al. (2009) only five clouds are at a distance of ∼7 ± 0.7 kpc (six including G045.14+00.14), but these contain the brightest objects in the field.

The submillimeter emission follows the molecular emission well (Figure 2). Clumps with typical star-forming characteristics (T ∼ 30–35 K for β = 1.5) are associated mainly with the three main complexes. Outside of the central regions the star-forming activity as traced in the submillimeter is less organized, although still following the overall contours of the molecular emission. More "isolated" clumps with IRAS counterparts (not all part of our final 7 kpc sample) show typical star-forming clump temperatures of ∼30 K and disturbed environments, with their stellar content generally resolved, and in an apparently less embedded stage of evolution.

Relatively cold structures in Aquila (T <20 K for β = 1.5 as measured in this analysis) have counterparts in major Infrared Dark Cloud (IRDC) surveys. Simon et al. (2006) detected 20 IRDC candidates in absorption with MSX, and Peretto & Fuller (2009) found 79 candidates using Spitzer data. We have checked for IRDC matches within 30'' of our submillimeter peaks in the entire field and found seven examples, all from the sample of Peretto & Fuller (2009). We note that several IRDCs, although not located within 30'' of our submillimeter peaks, are frequently clustered in their immediate neighborhood; these may have not be resolved in our data, or they may be confused by deconvolution structures. One prominent example of this is G045.49+00.04. Several IRDCs are observed in the local regions surrounding the brightest IRAS source formed by A28 and A29, especially toward A29. The intricate structures observed in the IRAC images, such as the pillar-like structures coincident with this BLAST source and the "cavity," will require submillimeter data with higher resolution and sensitivity in order to be properly characterized. We find a similar situation for the submillimeter complex labeled in Figure 2 as G046.34−00.21, where the majority of the submillimeter sources are identified as dark features in the IRAC maps. Many of these objects appear to be part of globular-like structures "blown" by the stellar activity located to the northeast. The characteristics of these and similar sources outside the distance range chosen for the present work (which focuses on the most prominent and active clumps of the field) will be investigated in detail with the release of new Herschel data from Hi-Gal.

Although a more in-depth study is required to confirm our results, the distribution and characteristics of the population within the main clumps suggest that triggering mechanisms (new star formation induced/initiated by external agents; e.g., Elmegreen 1992) are likely taking place in the region. Our multiwavelength (MIR–radio) analysis is compatible with induced activity from the innermost region to the farthest (parsec) extent of the parent clump.

4.2. Evolution and Massive Star Formation in Aquila

Within G045.14+00.14, our analysis in the IR, submillimeter, and radio, the molecular data, and the maser and stellar activity all support the scenario where A16 is in a more evolved stage of evolution than its neighbor A12. Figure 21 shows the measured bolometric luminosity for our BLAST sources relative to the combined total number of ionizing photons per unit time for all CORNISH sources possibly associated with each submillimeter source. The apparent underionizing state of A12 suggests optically thick radio emission, and therefore a highly embedded stage of the central sources.

Figure 21.

Figure 21. Comparison of Lbol (Table 6) and total Q0 for BLAST sources with analyzed CORNISH sources (Tables 3 and 4). Results are compared with theoretical curves for ZAMS (Panagia 1973) and class V (Martins et al. 2005) stars. An example of expected effect of increasing (decreasing) the distance by 1 kpc is shown by the diagonal line of source A12.

Standard image High-resolution image

A28 likely contains the "older" stellar activity of the observed clumps, while A29 and A30 appear to have characteristics of a relatively "younger" population. The later spectral type obtained from the CORNISH data for the main complex in A30 (the filamentary structure) relative to its estimate from the IR would support the embedded stage of the main star-forming activity in this region. Although our total luminosity measurement (and the radio estimates from Testi et al. 1999) is consistent with the theoretical curve in Figure 21, the large error range in luminosity for AA30 relative to the uncertainty in Q0 could still support this conclusion.

Despite the apparent agreement of AA29 (A28+A29) with the theoretical models in Figure 21, the more compact radio emission of A29 and the weak MIR emission (compared to the prominent submillimeter peak) would argue that A29 is the youngest active part of the entire cloud. Paron et al. (2009) suggested that G45.45 itself could have been triggered by three possible O-type stars ionizing the extended H ii region G45L. Although not analyzed in the present work, this is also a possibility for the original site ("first generation") of the activity within the BLAST clump.

A35 and A36 also have SED parameters consistent with ongoing stellar activity, but lack obvious 4.8 GHz detections. Neither of these two sources show maser emission in the catalogs used in this work, and both have weaker 13CO and CS antenna temperature. These sources may thus be in a relatively young and still embedded stage, or the star-forming activity may be low and/or without significant massive star formation.

Of the three main submillimeter peaks in the complex formed by G045.74−00.26 and G045.89−00.36, A50 has the strongest molecular emission. It is also possibly associated water masers, and has a relatively symmetric IR, submillimeter, and radio configuration. This source appears to have at least one possibly associated 4.8 GHz source, which would agree with previous studies that classified this source as an UCHiiR from the color criteria of WC89.

A46 shows the weakest molecular and radio emission, and no CORNISH counterparts or prominent signatures of star-forming activity. However, the presence of methanol emission and the possible association with an EGO suggest that this clump is in a very early and deeply embedded stage of massive star formation.

A47 shows not only more extended emission at all wavelengths, but it also appears to contain a few associated CORNISH sources within the main H ii region. The lack of any obvious maser emission, its morphology, the high dust temperatures, and weak molecular emission indicate that this clump has the most evolved stellar activity of the entire complex. The weak 13CO and CS antenna temperatures for this source could also be explained by the more evolved stage of the system.

From the results in Figure 21, both A47 and A50 appear to be underionizing for their luminosity. This could be an indication that the total luminosity has significant contribution from several "cooler" stars. However, the presence of water masers and CORNISH detections may still suggest that unaccounted very optically thick radio emission, from a new embedded generation of star formation in these clumps, could also be a possibility.

5. CONCLUSION

By means of the submillimeter maps provided by BLAST at 250, 350, and 500 μm, together with IRAC imaging, 21 cm continuum emission, 4.8 GHz CORNISH radio interferometry data, and GRS FCRAO 13CO/CS molecular data cubes, we have characterized the main clump population in the region of Aquila.

Our SED fitting of the most prominent star-forming sources at a distance of ∼7 kpc has revealed clumps with temperatures ranging between T ∼ 35 K and 40 K (for β = 1.5). Their total bolometric luminosities have equivalent single-star ZAMS spectral types earlier than ∼B0 stars.

The 4.8 GHz interferometry maps have been used to investigate the UCHiiRs within the submillimeter clumps. The ZAMS spectral types estimated from the CORNISH data are overall in good agreement with previous results. The "later" types obtained for some sources relative to their IR estimates support our BLAST and maser/outflow analysis, which suggest that highly embedded young OB stellar clusters are powering the most active clumps in the region.

Ongoing massive star formation is most prominent in G045.49+00.04 and G045.14+00.14 (containing our BLAST sources A12, A16, and A28, A29, A30, respectively). We confirm the presence of an OB stellar cluster deep within IRAS 19111+1048 (A16), containing a bright ∼O6 star surrounded by numerous late O and early B stars. This is in contrast with the scarcer maser and 4.8 GHz emission toward the complex formed by G045.74−00.26 and G045.89−00.36, in the field containing A46, A47, and A50. Despite the "disturbed" environments observed in the MIR, the evolutionary stage of the three main clumps in this complex appear to range from quite evolved (with extended and bright rimmed H ii regions and a prominent IR stellar population), to the earliest stages of massive star formation detected so far in this field (A46, prior to the onset of significant maser emission and UCHiiRs). The parameters derived from these SEDs are typical of active clumps and are indeed suggestive of massive star formation, although at a lower level to that observed in the previous two clouds.

The enhanced data sets that will be provided by Herschel (Hi-GAL) and SCUBA-2 at JCMT (JCMT Galactic Plane Survey; Moore et al. 2005) will be crucial to fully characterize and complete the census of the clump population in this field, especially the fainter and the coldest structures. In addition, the higher resolution of these instruments will help to probe down to core scales and to investigate the dusty structures within the clumps, which appear to be channeling the stellar activity and shaping the ionized regions in the neighborhood of the embedded OB stellar clusters.

The BLAST collaboration acknowledges the support of NASA through grants NAG5-12785, NAG5-13301, and NNGO-6GI11G, the Canadian Space Agency (CSA), the UK Particle Physics and Astronomy Research Council (PPARC), the Canada Foundation for Innovation (CFI), the Ontario Innovation Trust (OIT), and Canada's Natural Sciences and Engineering Research Council (NSERC). We also thank the Columbia Scientific Balloon Facility (CSBF) staff for their outstanding work. We also thank the referee for very useful suggestions and improvements to our paper, and Mubdi Rahman for useful discussions.

Footnotes

Please wait… references are loading.
10.1088/0004-637X/723/1/915