To distinguish between the different theories proposed to explain massive star formation, it is crucial to establish the distribution, the extinction, and the density of low-mass stars in massive star-forming regions. We analyzed deep X-ray observations of the Orion massive star-forming region using the Chandra Orion Ultradeep Project (COUP) catalog. We found that pre-main sequence (PMS) low-mass stars cluster toward the three massive star-forming regions: the Trapezium Cluster (TC), the Orion Hot Core (OHC), and OMC1-S. We derived low-mass stellar densities of 10^{5} stars pc^{-3} in the TC and OMC1-S, and of 10^{6} stars pc^{-3} in the OHC. The close association between the low-mass star clusters with massive star cradles supports the role of these clusters in the formation of massive stars. The X-ray observations show for the first time in the TC that low-mass stars with intermediate extinction are clustered toward the position of the most massive star, which is surrounded by a ring of non-extincted low-mass stars. Our analysis suggests that at least two basic ingredients are needed in massive star formation: the presence of dense gas and a cluster of low-mass stars. The scenario that better explains our findings assumes high fragmentation in the parental core, accretion at subcore scales that forms a low-mass stellar cluster, and subsequent competitive accretion.

Outflows arising from very young stars affect their surroundings and influence the star formation in the parental core. Multiple molecular outflows and Herbig-Haro (HH) objects have been observed in Orion, many of them originating from the embedded massive star-forming region known as OMC1-S. The detection of the outflow driving sources is commonly difficult, because they are still hidden behind large extinction, preventing their direct observation at optical and even near and mid-IR wavelengths. With the aim of improving the identification of the driving sources of the multiple outflows detected in OMC1-S, we used the catalog provided by deep X-ray observations, which have unveiled the very embedded population of pre-main sequence stars. We compared the position of stars observed by the Chandra Orion Ultra Deep project (COUP) in OMC1-S with the morphology of the molecular outflows and the directions of measured proper motions of HH optical objects. We find that 6 out of 7 molecular outflows reported in OMC1-S (detection rate of 86 %) have an extincted X-ray COUP star located at the expected position of the driving source. In several cases, X-rays detected the possible driving sources for the first time. This clustered embedded population revealed by Chandra is very young, with an estimated average age of few 10^{5} yr. It is also likely responsible for the multiple HH objects, which are the optical correspondence of flows arising from the cloud. We show that the molecular outflows driven by the members of the OMC1-S cluster can account for the observed turbulence at core-scales and regulate the star formation efficiency. We discuss the effects of outflow feedback in the formation of massive stars, concluding that the injected turbulence in OMC1-S is compatible with a competitive accretion scenario.

To distinguish between the different theories proposed to explain massive star formation, it is crucial to establish the distribution, the extinction, and the density of low-mass stars in massive star-forming regions. We analyze deep X-ray observations of the Orion massive star-forming region using the Chandra Orion Ultradeep Project (COUP) catalog. We studied the stellar distribution as a function of extinction, with cells of 0.03 pc x 0.03 pc, the typical size of protostellar cores. We derived stellar density maps and calculated cluster stellar densities. We found that low-mass stars cluster toward the three massive star-forming regions: the Trapezium Cluster (TC), the Orion Hot Core (OHC), and OMC1-S. We derived low-mass stellar densities of 10^{5} stars pc^{-3} in the TC and OMC1-S, and of 10^{6} stars pc^{-3} in the OHC. The close association between the low-mass star clusters with massive star cradles supports the role of these clusters in the formation of massive stars. The X-ray observations show for the first time in the TC that low-mass stars with intermediate extinction are clustered toward the position of the most massive star, which is surrounded by a ring of non-extincted low-mass stars. This 'envelope-core' structure is also supported by infrared and optical observations. Our analysis suggests that at least two basic ingredients are needed in massive star formation: the presence of dense gas and a cluster of low-mass stars. The scenario that better explains our findings assumes high fragmentation in the parental core, accretion at subcore scales that forms a low-mass stellar cluster, and subsequent competitive accretion. Finally, although coalescence does not seem a common mechanism for building up massive stars, we show that a single stellar merger may have occurred in the evolution of the OHC cluster, favored by the presence of disks, binaries, and gas accretion.

The detection of narrow SiO thermal emission toward young outflows has been proposed to be a signature of the magnetic precursor of C-shocks. Recent modeling of the SiO emission across C-shocks predicts variations in the SiO line intensity and line shape at the precursor and intermediate-velocity regimes in only few years. We present high-angular resolution (3.8"x3.3") images of the thermal SiO J=2-1 emission toward the L1448-mm outflow in two epochs (November 2004-February 2005, March-April 2009). Several SiO condensations have appeared at intermediate velocities (20-50 km/s) toward the red-shifted lobe of the outflow since 2005. Toward one of the condensations (clump D), systematic differences of the dirty beams between 2005 and 2009 could be responsible for the SiO variability. At higher velocities (50-80 km/s), SiO could also have experienced changes in its intensity. We propose that the SiO variability toward L1448-mm is due to a real SiO enhancement by young C-shocks at the internal working surface between the jet and the ambient gas. For the precursor regime (5.2-9.2 km/s), several narrow and faint SiO components are detected. Narrow SiO tends to be compact, transient and shows elongated (bow-shock) morphologies perpendicular to the jet. We speculate that these features are associated with the precursor of C-shocks appearing at the interface of the new SiO components seen at intermediate velocities.

We present high angular resolution PdBI images (beam of ~0.33'') of the J=27-26 line from several vibrational levels (v7=1 and v6=1) of HC3N toward Cepheus A HW2. These images reveal the two main heating sources in the cluster: one centered in the disk collimating the HW2 radio jet (the HW2 disk), and the other associated with a hot core 0.3'' northeast HW2 (the HC). This is the first time that vibrationally excited emission of HC3N is spatially resolved in a disk. The kinematics of this emission shows that the HW2 disk rotates following a Keplerian law. We derive the temperature profiles in the two objects from the excitation of HC3N along the HW2 disk and the HC. These profiles reveal that both objects are centrally heated and show temperature gradients. The inner and hotter regions have temperatures of 350+-30K and 270+-20K for the HW2 disk and the HC, respectively. In the cooler and outer regions, the temperature drops to 250+-30K in the HW2 disk, and to 220+-15K in the HC. The estimated luminosity of the heating source of the HW2 disk is ~2.2E4Lo, and 3000Lo for the HC. The most massive protostar in the HW2 region is the powering source of the HW2 radio jet. We discuss the formation of multiple systems in this cluster. The proximity of the HC to HW2 suggest that these sources likely form a binary system of B stars, explaining the observed precession of the HW2 radio jet.

We present the first results of the expected variations of the molecular line emission arising from material recently affected by C-shocks (shock precursors). Our parametric model of the structure of C-shocks has been coupled with a radiative transfer code to calculate the molecular excitation and line profiles of shock tracers such as SiO, and of ion and neutral molecules such as H13CO+ and HN13C, as the shock propagates through the unperturbed medium. Our results show that the SiO emission arising from the early stage of the magnetic precursor typically has very narrow line profiles slightly shifted in velocity with respect to the ambient cloud. This narrow emission is generated in the region where the bulk of the ion fluid has already slipped to larger velocities in the precursor as observed toward the young L1448-mm outflow. This strongly suggests that the detection of narrow SiO emission and of an ion enhancement in young shocks, is produced by the magnetic precursor of C-shocks. In addition, our model shows that the different velocity components observed toward this outflow can be explained by the coexistence of different shocks at different evolutionary stages, within the same beam of the single-dish observations.

We present VLA and PdBI subarcsecond images (0.15"-0.6") of the radiocontinuum emission at 7 mm and of the SO2 J=19_{2,18}-18_{3,15} and J=27_{8,20}-28_{7,21} lines toward the Cepheus A HW2 region. The SO2 images reveal the presence of a hot core internally heated by an intermediate mass protostar, and a circumstellar rotating disk around the HW2 radio jet with size 600AUx100AU and mass of 1M_sun. Keplerian rotation for the disk velocity gradient of 5 kms-1 requires a 9 M_sun central star, which cannot explain the total luminosity observed in the region. This may indicate that the disk does not rotate with a Keplerian law due to the extreme youth of this object. Our high sensitivity radiocontinuum image at 7 mm shows in addition to the ionized jet, an extended emission to the west (and marginally to the south) of the HW2 jet, filling the south-west cavity of the HW2 disk. From the morphology and location of this free-free continuum emission at centimeter and millimeter wavelengths (spectral index of 0.4-1.5), we propose that the disk is photoevaporating due to the UV radiation from the central star. All this indicates that the Cepheus A HW2 region harbors a cluster of massive stars. Disk accretion seems to be the most plausible way to form massive stars in moderate density/luminosity clusters.

Magnetic precursors of C-shocks accelerate, compress and heat molecular ions, modifying the kinematics and the physical conditions of the ion fluid with respect to the neutral one. Electron densities are also expected to be significantly enhanced in shock precursors. In this Letter, we present observations of strongly polar ion and neutral molecules such as SiO, H13CO+, HN13C and H13CN, which reveal the electron density enhancements associated with the precursor of the young L1448-mm outflow. While in the ambient gas the excitation of the ions and neutrals is explained by collisional excitation by H2 with a single density of 10E5cmE-3, H13CO+ shows an over excitation in the shock precursor component that requires H2 densities of a factor of >10 larger than those derived from the neutral species. This over excitation in H13CO+ can be explained if we consider an additional excitation by collisions with electrons and an electron density enhancement in the precursor stage by a factor of 500, i.e. a fractional ionization of 5x10E-5. These results show that multiline observations can be used to study the evolution of the ion and electron fluids at the first stages of the C-shock interaction.

We present the discovery of the first molecular hot core associated with an intermediate mass protostar in the CepA HW2 region. The hot condensation was detected from single dish and interferometric observations of several high excitation rotational lines (from 100 to 880K above the ground state) of SO2 in the ground vibrational state and of HC3N in the vibrationally excited states v7=1 and v7=2. The kinetic temperature derived from both molecules is 160K. The high-angular resolution observations (1.25'' x 0.99'') of the SO2 J=28(7,21)-29(6,24) line (488K above the ground state) show that the hot gas is concentrated in a compact condensation with a size of 0.6''(430AU), located 0.4'' (300AU) east from the radio-jet HW2. The total SO2 column density in the hot condensation is 10E18cm-2, with a H2 column density ranging from 10E23 to 6 x 10E24cm-2. The H2 density and the SO2 fractional abundance must be larger than 10E7cm-3 and 2 x 10E-7 respectively. The most likely alternatives for the nature of the hot and very dense condensation are discussed. From the large column densities of hot gas, the detection of the HC3N vibrationally excited lines and the large SO2 abundance, we favor the interpretation of a hot core heated by an intermediate mass protostar of 10E3 Lo. This indicates that the CepA HW2 region contains a cluster of very young stars.

The recent detection of shock-precursors toward the very young L1448-mm outflow offers us the possibility to study the grain chemistry during the first stages of the shock evolution, constraining the molecules ejected from grains and the species formed in gas phase. Observations of key molecules in the grain chemistry such as SiO, CH3OH, SO, CS, H2S, OCS, and SO2 toward this outflow are presented. The line profiles and the derived abundances show three distinct velocity regimes that trace the shock evolution: the preshock, the shock-precursor and the postshock gas. The SiO, CH3OH, SO, and CS abundances are enhanced with respect to the quiescent gas by 1 order of magnitude in the shock-precursor component, and by 3 orders of magnitude in the postshock gas. The derived SiO and CH3OH abundances are consistent with the recent ejection of these molecules from grains. Since H2S is only enhanced in the shock-precursor component, and OCS and SO2 are undetected, SO and CS are the most abundant sulfur-bearing species in the grain mantles of L1448-mm. The ejection of mainly SO and CS rather than H2S or OCS from grains, suggests that the sulfur chemistry will depend on the chemical "history" of the grain mantles in outflows and hot cores.

We present the detection of the SiO v=0 J=2-1 and J=3-2 lines, and of the HCO 101-000 J=3/2-1/2 F=2-1 line at ambient velocities towards the molecular outflows in L1448-mm and L1448-IRS3. This is the first detection of HCO in a dark cloud. We have also measured lines of H13CO+, H13CN, HN13C, CH3OH, and N2H+. While the HCO and the SiO lines have the narrowest profiles with linewidths of 0.5 km s-1, the other lines have widths of 1km s-1. Towards L1448-mm, all lines except those of SiO and HCO, show two distinct velocity components centered at 4.7 and 5.2 km s-1. HCO is only observed in the 4.7 km s-1 cloud, and SiO in the 5.2 km s-1 component. The SiO abundance is 10E-11 in the 5.2 km s-1 clouds, one order of magnitude larger than in the 4.7 km s-1 component and in other dark clouds. The HCO abundance is 10E-11, similar to that predicted by the ion-molecule reactions models for the quiescent gas in dark clouds. The large change in the SiO/HCO abundance ratio (>150) from the 4.7 to the 5.2 km s-1 component, and the distribution and kinematics of the SiO emission towards L1448-mm suggest that the ambient SiO is associated with the molecular outflows. We propose that the narrow linewidths and the abundances of SiO in the ambient gas are produced by the interaction of the magnetic and/or radiative precursors of the shocks with the clumpy pre-shocked ambient gas.

We have carried out a high-sensitivity search for circumstellar disks around Herbig Be stars in the continuum at 1.4mm and 2.7mm using the IRAM interferometer at the Plateau de Bure (PdBI) . In this letter, we report data on three well studied B0 stars, MWC 1080, MWC 137 and R Mon. The two latter have also been observed in the continuum at 0.7 cm and 1.3 cm using the NRAO Very Large Array (VLA) . We report the detection of circumstellar disks around MWC 1080 and R Mon with masses of Md ~ 0.003 and 0.01 Msun, respectively, while for MWC 137 we estimate a disk mass upper limit of 0.007 Msun. Our results show that the ratio Md/M* is at least an order of magnitude lower in Herbig Be stars than in Herbig Ae and T Tauri stars.

We report the first detection of the reactive ions CO+ and HOC+ towards ultracompact (UC) HII regions, particularly in Mon R2 and G29.96-0.02. We have observed two positions in Mon R2, namely the peak of the UC HII region and a position in the high density molecular cloud which bounds it. CO+ and HOC+ were detected at the UC HII region but not at the molecular cloud, as expected if the CO$^+$ and HOC$^+$ emissions arise in the PDR surrounding the \uch. The measured CO$^+$ and HOC$^+$ column densities are of the order of 10$^{11}$ cm$^{-2}$ in both sources, which yields a strikingly low [HCO$^+$]/[HOC$^+$] abundance ratio of 460 in Mon R2. These values are similar to those found in some other well-known PDRs, like NGC 7023 and the Orion Bar. We briefly discuss the chemical implications of these results.

We present images of the J=10-9 rotational lines of HC3N in the vibrationally excited levels 1v7, 1v6 and 1v5 of the hot core (HC) in Orion KL. The images show that the spatial distribution and the size emission from the 1v7 and 1v5 levels are different. While the J=10-9 1v7 line has a size of 4''x 6'' and peaks 1.1'' NE of the 3 mm continuum peak, the J=10--9 1v5 line emission is unresolved (<3'') and peaks 1.3'' south of the 3 mm peak. This is a clear indication that the HC is composed of condensations with very different temperatures (170 K for the 1v7 peak and $>230$ K for the 1v5 peak). The temperature derived from the 1v7 and 1v5 lines increases with the projected distance to the suspected main heating source I. Projection effects along the line of sight could explain the temperature gradient as produced by source I. However, the large luminosity required for source I, >5 10^5 Lsolar, to explain the 1v5 line suggests that external heating by this source may not dominate the heating of the HC. Simple model calculations of the vibrationally excited emission indicate that the HC can be internally heated by a source with a luminosity of 10^5 Lsolar, located 1.2'' SW of the 1v5 line peak (1.8'' south of source I). We also report the first detection of high-velocity gas from vibrationally excited HC3N emission. Based on excitation arguments we conclude that the main heating source is also driving the molecular outflow. We speculate that all the data presented in this letter and the IR images are consistent with a young massive protostar embedded in an edge-on disk.

We present HHT observations in the N=3-2 rotational transition of the CN radical toward selected positions of the Trapezium region and of the molecular Ridge in the Orion molecular cloud. Two of the positions in the Ridge were also observed in the N=2-1 line of CN and 13CN. The N=3-2 CN lines have been combined with observations of the N=2-1 and N=1-0 transitions of CN, and of the N=2-1 of 13CN to estimate the physical conditions and CN abundances in the molecular gas. We analyze in detail the excitation of the CN lines and find that the hyperfine ratios of the N=3-2 line are always close to the Local Thermodynamic Equilibrium (LTE) optically thin values even in the case of optically thick emission. This is due to different excitation temperatures for the different hyperfine lines. From the line intensity ratios between the different CN transitions we derive H_2 densities of ~10^5cm^-3 for the molecular Ridge and of ~3x10^6cm^-3 for the Trapezium region. The CN column densities are one order of magnitude larger in the Ridge than in the Trapezium region, but the CN to H_2 ratio is similar both in the Trapezium and in the Ridge. The combination of the low CN column densities, high H_2 densities and relatively high CN abundances toward the Trapezium region requires that the CN emission arises from a thin layer with a depth along the line of sight of only \~5x10^15cm. This high density thin layer of molecular gas seems to be related with material that confines the rear side of the HII region Orion A. However the molecular layer is not moving as expected from the expansion of the HII region, but it is ``static'' with respect to the gas in the molecular cloud. We discuss the implication of a high density ``static'' layer in the evolution of an HII region.

This paper studies the global ISM content in a sample of 104 accreting galaxies, including counterrotators and polar rings, which spans the entire Hubble sequence. The molecular, atomic and hot gas content of accretors is compared to a newly compiled sample of normal galaxies. We present results of a small survey of the J=1--0 line of {12}CO with the 15m SEST telescope on a sample of 11 accretors (10 counterrotators and 1 polar ring). The SEST sample is enlarged with published data from 48 galaxies, for which observational evidence of counterrotation in the gas and/or the stars has been found. Furthermore, the available data on a sample of 46 polar ring galaxies has been compiled. It was found that the normalized content of cold gas (M_gas/L_B) in polar rings is ~1 order of magnitude higher than the reference value derived for normal galaxies. The inferred gas masses are sufficient to stabilize polar rings through self-gravity. In contrast, it was found that the cold gas content of counterrotators is close to normal for all galaxy types. Although counterrotators and polar rings probably share a common origin, the gas masses estimated here confirm that light gas rings accreted by future counterrotators may have evolved faster than the self-gravitating structures of polar rings. In this scenario, the transformation of atomic into molecular gas could be enhanced near the transition region between the prograde and the retrograde disks, especially in late-type accretors characterized by a high content of primordial gas. This is tentatively confirmed in this work: the measured H_2/HI ratio seems larger in counterrotators than in normal or polar ring galaxies for types later than S0s.

We present high sensitivity and high angular resolution images of the high velocity (vLSR>30kms^-1) CO emission in the J=1-0 and J=2-1 lines of the Orion KL region. These results reveal the morphology of the high-velocity CO emission at the most extreme velocities. High velocity emission have been only detected in two regions: BN/KL (IRc2/I) and Orion-S. The Orion-S region contains a very young (dynamical age of 10^3years), very fast (~110kms^-1) and very compact (<0.16pc) bipolar outflow. From the morphology of the high-velocity gas we estimate that the position of the powering source must be ~20'' north of FIR4. For the IRc2/I molecular outflow the morphology of the moderate velocity (<60kms^-1) gas shows a weak bipolarity around IRc2/I. The gas at the most extreme velocities does not show any bipolarity around IRc2/I, if any, it is found ~30'' north from these sources. The blue and redshifted gas at moderate velocities shows similar spatial distribution with a systematic trend for the size of the high-velocity gas to decrease as the terminal radial velocity increases. The size-velocity relationship is fitted with a simple velocity law which considers a highly collimated jet and entrained material outside the jet moving in the radial direction. We also find that most of the CO outflowing at moderate velocities is located at the head of the jet. Our results and the spatial distribution and kinematics of the shock tracers in this outflow can be explained if the IRc2/I outflow is driven by a precessing jet oriented along the line of sight.

We have observed the S(0), S(1), S(2), S(3), S(4) and S(5) rotational lines of molecular hydrogen (H2) towards the peak of the photodissociation region (PDR) associated with the reflection nebula NGC 7023. The observed H2 line ratios show that they arise in warm gas with kinetic temperatures ~300 - 700 K. However, the data cannot be fitted by an ortho- to para- (OTP) ratio of 3. An OTP ratio in the range ~1.5 - 2 is necessary to explain our observations. This is the first detection of a non-equilibrium OTP ratio measured from the H2 pure-rotational lines in a PDR. The existence of a dynamical PDR is discussed as the most likely explanation for this low OTP ratio.

We present high sensitivity maps of the High Velocity (HV) CO emission toward the molecular outflows around IRc2 and Orion-S in the Orion A molecular cloud. The maps reveal the presence of HV bullets in both outflows with velocities between 40-80kms-1 from the ambient gas velocity. The blue and redshifted CO HV bullets associated with the IRc2 outflow are distributed in thin (12"-20", 0.02-0.04pc) elliptical ring-like structures with a size of ~10"x50" (0.02x0.1pc). The CO emission at the most extreme blue and redshifted velocities (EHV) peaks 20" north of source I, just inside the rings of the HV bullets. The low velocity H2O masers and the H2* bullets around IRc2 are located at the inner edges of the ring of CO HV bullets and surrounding the EHV CO emission. Furthermore, the high velocity H2O masers are very well correlated with the EHV CO emission. This morphology is consistent with a model of a jet driven molecular outflow oriented close to the line of sight. In the Orion-S outflow, the morphology of the CO HV bullets shows a bipolar structure in the southeast-northwest direction, and the H2O masers are found only at low velocities in the region between the exciting source and the CO HV bullets. The morphology of the CO HV bullets, the radial velocities and the spatial distribution of the H2O masers in both outflows, as well as the H2* features around IRc2, are consistent with a model in which these outflows are driven by a jet variable in direction.In this scenario, the large traverse velocity measured for the H2O masers in the IRc2 outflow, ~18kms-1, supports the evolutionary connection between the jet and the shell-like outflows.