We present the results of a narrow-band near-infrared imaging survey for Molecular Hydrogen emission-line Objects (MHOs) toward 26 regions containing high-mass protostellar candidates and massive molecular outflows. We have detected a total of 236 MHOs, 156 of which are new detections, in 22 out of the 26 regions. We use H$_2$ 2.12-$\mu$m/H$_2$ 2.25-$\mu$m flux ratios, together with morphology, to separate the signatures of fluorescence associated with photo-dissociation regions (PDRs) from shocks associated with outflows in order to identify the MHOs. PDRs have typical low flux ratios of ~ 1.5 - 3, while the vast majority of MHOs display flux ratios typical of C-type shocks (~ 6-20). A few MHOs exhibit flux ratios consistent with expected values for J-type shocks (~ 3-4), but these are located in regions that may be contaminated with fluorescent emission. Some previously reported MHOs have low flux ratios, and are likely parts of PDRs rather than shocks indicative of outflows. We identify a total of 36 outflows across the 22 target regions where MHOs were detected. In over half these regions, MHO arrangements and fluorescent structures trace features present in CO outflow maps, suggesting the CO emission traces a combination of dynamical effects, which may include gas entrained in expanding PDRs as well as bipolar outflows. Where possible, we link MHO complexes to distinct outflows and identify candidate driving sources.

We have searched the L1340 A, B, and C clouds for shocks from protostellar outflows using the H$_2$ 2.122 $\mu$m near-IR line as a shock tracer. Substantial outflow activity has been found in each of the three regions of the cloud (L1340 A, L1340 B, & L1340 C). We find 42 distinct shock complexes (16 in L1340 A, 11 in L1340 B, and 15 in L1340 C). We were able to link 17 of those shock complexes in to 12 distinct outflows and identify candidate source stars for each. We examine the properties (A$_{V}$, T$_{bol}$, and L$_{bol}$) of the source protostars and compare that to the properties of the general population of Class 0/I and flat SED protostars and find that there is an indication, albeit at low statistical significance, that the outflow driving protostars are drawn from a population with lower A$_{V}$, higher L$_{bol}$, and lower T$_{bol}$ than the general population of protostars.

We present results of an infrared study of the molecular cloud Lynds 1340, forming three groups of low and intermediate-mass stars. Our goals are to identify and characterise the young stellar population of the cloud, study the relationships between the properties of the cloud and the emergent stellar groups, and integrate L1340 into the picture of the star-forming activity of our Galactic environment. We selected candidate young stellar objects from the Spitzer and WISE data bases using various published color criteria, and classified them based on the slope of the spectral energy distribution. We identified 170 Class II, 27 Flat SED, and Class 0/I sources. High angular resolution near-infrared observations of the RNO 7 cluster, embedded in L1340, revealed eight new young stars of near-infrared excess. The surface density distribution of young stellar objects shows three groups, associated with the three major molecular clumps of L1340, each consisting of less than 100 members, including both pre-main sequence stars and embedded protostars. New Herbig--Haro objects were identified in the Spitzer images. Our results demonstrate that L1340 is a prolific star-forming region of our Galactic environment in which several specific properties of the intermediate-mass mode of star formation can be studied in detail.

Yellowballs are a collection of approximately 900 compact, infrared sources identified and named by volunteers participating in the Milky Way Project (MWP), a citizen-science project that uses GLIMPSE/MIPSGAL images from Spitzer to explore topics related to Galactic star formation. In this paper, through a combination of catalog cross-matching and infrared color analysis, we show that yellowballs are a mix of compact star-forming regions, including ultra-compact and compact HII regions, as well as analogous regions for less massive B-type stars. The resulting MWP yellowball catalog provides a useful complement to the Red MSX Source (RMS) survey. It similarly highlights regions of massive star formation, but the selection of objects purely on the basis of their infrared morphology and color in Spitzer images identifies a signature of compact star-forming regions shared across a broad range of luminosities, and by inference, masses. We discuss the origin of their striking mid-infrared appearance, and suggest that future studies of the yellowball sample will improve our understanding of how massive and intermediate-mass star-forming regions transition from compact to more extended bubble-like structures.

We report on the discovery of three candidate eruptive young stars, found during our comprehensive multi-wavelength study of the young stellar population of the dark cloud L1340. These stars are as follows. (1) IRAS 02224+7227 (2MASS 02270555+7241167, HH 487 S) exhibited FUor-like spectrum in our low-resolution optical spectra. The available photometric data restrict its luminosity to 23 L_sun < L_bol < 59 L_sun. (2) 2MASS 02263797+7304575, identified as a classical T Tauri star during our H alpha survey, exhibited an EXor type brightening in 2005 November, at the time of the SDSS observations of the region. (3) 2MASS 02325605+7246055, a low-mass embedded young star, associated with a fan-shaped infrared nebula, underwent an outburst between the DSS1 and DSS2 surveys, leading to the appearance of a faint optical nebula. Our [SII] and H alpha images, as well as the Spitzer IRAC 4.5 micron images revealed Herbig-Haro objects associated with this star. Our results suggest that amplitudes and time scales of outbursts do not necessarily correlate with the evolutionary stage of the stars.

We have discovered 12 new molecular hydrogen emission-line objects (MHOs) in the vicinity of the candidate massive young stellar object Mol 121, in addition to five that were previously known. H2 2.12-micron/H2 2.25-micron flux ratios indicate another region dominated by fluorescence from a photo-dissociation region (PDR), and one region that displays an anomalously low H2 2.12-micron/H2 2.25-micron flux ratio (<1) and coincides with a previously reported deeply embedded source (DES). Continuum observations at 3 mm reveal five dense cores; the brightest core is coincident with the DES. The next brightest cores are both associated with cm continuum emission. One of these is coincident with the IRAS source; the other lies at the centroid of a compact outflow defined by bipolar MHOs. The brighter of these bipolar MHOs exhibits [Fe II] emission and both MHOs are associated with CH3OH maser emission observed at 95 GHz and 44 GHz. Masses and column densities of all five cores are consistent with theoretical predictions for massive star formation. Although it is impossible to associate all MHOs with driving sources in this region, it is evident that there are several outflows along different position angles, and some unambiguous associations can be made. We discuss implications of observed H2 2.12-micron/H2 2.25-micron and [Fe II] 1.64-micron/H2 2.12-micron flux ratios and compare the estimated total H2 luminosity with the bolometric luminosity of the region. We conclude that the outflows are driven by massive young stellar objects embedded in cores that are likely to be in different evolutionary stages.

The Milky Way Project citizen science initiative recently increased the number of known infrared bubbles in the inner Galactic plane by an order of magnitude compared to previous studies. We present a detailed statistical analysis of this dataset with the Red MSX Source catalog of massive young stellar sources to investigate the association of these bubbles with massive star formation. We particularly address the question of massive triggered star formation near infrared bubbles. We find a strong positional correlation of massive young stellar objects (MYSOs) and H II regions with Milky Way Project bubbles at separations of < 2 bubble radii. As bubble sizes increase, a statistically significant overdensity of massive young sources emerges in the region of the bubble rims, possibly indicating the occurrence of triggered star formation. Based on numbers of bubble-associated RMS sources we find that 67+/-3% of MYSOs and (ultra)compact H II regions appear associated with a bubble. We estimate that approximately 22+/-2% of massive young stars may have formed as a result of feedback from expanding H II regions. Using MYSO-bubble correlations, we serendipitously recovered the location of the recently discovered massive cluster Mercer 81, suggesting the potential of such analyses for discovery of heavily extincted distant clusters.

We present a new catalogue of 5,106 infrared bubbles created through visual classification via the online citizen science website 'The Milky Way Project'. Bubbles in the new catalogue have been independently measured by at least 5 individuals, producing consensus parameters for their position, radius, thickness, eccentricity and position angle. Citizen scientists - volunteers recruited online and taking part in this research - have independently rediscovered the locations of at least 86% of three widely-used catalogues of bubbles and H ii regions whilst finding an order of magnitude more objects. 29% of the Milky Way Project catalogue bubbles lie on the rim of a larger bubble, or have smaller bubbles located within them, opening up the possibility of better statistical studies of triggered star formation. Also outlined is the creation of a 'heat map' of star-formation activity in the Galactic plane. This online resource provides a crowd-sourced map of bubbles and arcs in the Milky Way, and will enable better statistical analysis of Galactic star-formation sites.

We have discovered two compact sources of shocked H2 2.12-micron emission coincident with Mol 160 (IRAS 23385+6053), a massive star-forming core thought to be a precursor to an ultracompact HII region. The 2.12-micron sources lie within 2" (0.05 pc) of a millimeter-wavelength continuum peak where the column density is >= 10e24 cm$^{-2}$. We estimate that the ratio of molecular hydrogen luminosity to bolometric luminosity is > 0.2%, indicating a high ratio of mechanical to radiant luminosity. CS J=2-1 and HCO$^+$ J=1-0 observations with CARMA indicate that the protostellar molecular core has a peculiar velocity of ~ 2 km s$^{-1}$ with respect to its parent molecular cloud. We also observed 95 GHz CH3OH J=8$-7 Class I maser emission from several locations within the core. Comparison with previous observations of 44-GHz CH3OH maser emission shows the maser sources have a high mean ratio of 95-GHz to 44-GHz intensity. Our observations strengthen the case that Mol 160 (IRAS 23385+6053) is a rapidly accreting massive protostellar system in a very early phase of its evolution.

We present ground-based mid-infrared observations of Class 0 protostars in LDN 1448. Of the five known protostars in this cloud, we detected two, L1448N:A and L1448C, at 12.5, 17.9, 20.8, and 24.5 microns, and a third, L1448 IRS 2, at 24.5 microns. We present high-resolution images of the detected sources, and photometry or upper limits for all five Class 0 sources in this cloud. With these data, we are able to augment existing spectral energy distributions (SEDs) for all five objects and place them on an evolutionary status diagram.

We present far-infrared, submillimetre, and millimetre observations of bright IRAS sources and outflows that are associated with massive CS clumps in the Monoceros OB1 Dark Cloud. Individual star-forming cores are identified within each clump. We show that combining submillimetre maps, obtained with SCUBA on the JCMT, with HIRES-processed and modelled IRAS data is a powerful technique that can be used to place better limits on individual source contributions to the far-infrared flux in clustered regions. Three previously categorized "Class I objects" are shown to consist of multiple sources in different evolutionary stages. In each case, the IRAS point source dominates the flux at 12 and 25 microns. In two cases, the IRAS point source is not evident at submillimetre wavelengths. The submillimetre sources contribute significantly to the 60 and 100 micron fluxes, dominating the flux in the 100 micron waveband. Using fluxes derived from our technique, we present the spectral energy distribution and physical parameters for an intermediate-mass Class 0 object in one of the regions. Our new CO J=2-1 outflow maps of the three regions studied indicate complex morphology suggestive of multiple driving sources. We discuss the possible implications of our results for published correlations between outflow momentum deposition rates and "source" luminosities, and for using these derived properties to estimate the ratio of mass ejection rates to mass accretion rates onto protostars.

We present sensitive, large-scale maps of the CO J=1-0 emission of the L1448 dark cloud. These maps were acquired using the On-The-Fly capability of the NRAO 12-meter telescope. CO outflow activity is seen in L1448 on parsec-scales for the first time. Careful comparison of the spatial and velocity distribution of our high-velocity CO maps with previously published optical and near-infrared images and spectra has led to the identification of six distinct CO outflows. We show the direct link between the heretofore unknown, giant, highly-collimated, protostellar molecular outflows and their previously discovered, distant optical manifestations. The outflows traced by our CO mapping generally reach the projected cloud boundaries. Integrated intensity maps over narrow velocity intervals indicate there is significant overlap of blue- and red-shifted gas, suggesting the outflows are highly inclined with respect to the line-of-sight, although the individual outflow position angles are significantly different. The velocity channel maps also show that the outflows dominate the CO line cores as well as the high-velocity wings. The magnitude of the combined flow momenta, as well as the combined kinetic energy of the flows, are sufficient to disperse the 50 solar mass NH3 cores in which the protostars are currently forming, although some question remains as to the exact processes involved in redirecting the directionality of the outflow momenta to effect the complete dispersal of the parent cloud.