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 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.

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.