H2 emission from shocks in molecular outflows: the significance of departures from a stationary state

We have computed the time dependence of the H₂ rovibrational emission spectrum from molecular outflows. This emission arises in shock waves generated by the impact of jets, associated with low-mass star formation, on molecular gas. The shocks are unlikely to have attained a state of equilibrium, and so their structure will exhibit both C- and J-type characteristics. The rotational excitation diagram is found to provide a measure of the age of the shock; in the case of the outflow observed in Cepheus A West by the ISO satellite, the shock age is found to be approximately 1.5×10³ yr. Emission by other species, such as NH₃ and SiO, is also considered, as are the intensities of the fine-structure transitions of atoms and ions.     

Unlocking the Keyhole: H2 and PAH emission from molecular clumps in the Keyhole Nebula

To better understand the environment surrounding CO emission clumps in the Keyhole Nebula, we have made images of the region in H₂ 1–0 S(1) (2.122-μm) emission and polycyclic aromatic hydrocarbon (PAH) emission at 3.29 μm. Our results show that the H₂ and PAH emission regions are morphologically similar, existing as several clumps, all of which correspond to CO emission clumps and dark optical features. The emission confirms the existence of photodissociation regions (PDRs) on the surface of the clumps. By comparing the velocity range of the CO emission with the optical appearance of the H₂ and PAH emission, we present a model of the Keyhole Nebula whereby the most negative velocity clumps are in front of the ionization region, the clumps at intermediate velocities are in it and those which have the least negative velocities are at the far side. It may be that these clumps, which appear to have been swept up from molecular gas by the stellar winds from η  Car, are now being overrun by the ionization region and forming PDRs on their surfaces. These clumps comprise the last remnants of the ambient molecular cloud around η Car.     

The contributions of J-type shocks to the H2 emission from molecular outflow sources

We present the results of modelling of the H₂ emission from molecular outflow sources, induced by shock waves propagating in the gas. We emphasize the importance of proper allowance for departures from equilibrium owing to the finite flow velocity of the hot, compressed gas, with special reference to the excitation, dissociation and reformation of H₂. The salient features of our computer code are described. The code is applied to interpreting the spectra of the outflow sources Cepheus A West and HH43. Particular attention is paid to determining the cooling times in shocks whose speeds are sufficient for collisional dissociation of H₂ to take place; the possible observational consequences of the subsequent reformation of H₂ are also examined. Because molecular outflow sources are intrinsically young objects, J-type shocks may be present in conjunction with magnetic precursors, which have a C-type structure. We note that very different physical and dynamical conditions are implied by models of C- and J-type shocks which may appear to fit the same H₂ excitation diagram.     

On the photodissociation of H2 by the first stars

The first star formation in the Universe is expected to take place within small protogalaxies, in which the gas is cooled by molecular hydrogen. However, if massive stars form within these protogalaxies, they may suppress further star formation by photodissociating the H₂. We examine the importance of this effect by estimating the time-scale on which significant H₂ is destroyed. We show that photodissociation is significant in the least massive protogalaxies, but becomes less so as the protogalactic mass increases. We also examine the effects of photodissociation on dense clumps of gas within the protogalaxy. We find that while collapse will be inhibited in low-density clumps, denser ones may survive to form stars.     

Fabry–Perot imaging of NGC 6090 in H2 emission

We have obtained an H₂ v =1–0 S(1) image of a merging galaxy system, NGC 6090, by using a Fabry–Perot imager. The H₂ emission originates between the double nuclei of pre-merger galaxies, and exhibits an arc-like or ring-like structure almost connecting the double nuclei. This structure is similar to that suggested for Arp 220 from the velocity field measured by CO radio emission. The separation of the double nuclei in NGC 6090 is 5–6 arcsec, corresponding to a projected distance of 3–4 kpc. This is much larger than that of Arp 220 and suggests that the molecular gas distribution can form an organized shape between the nuclei, such as a ring, in a rather early phase of merging.     

The ortho- to para- ratio of H2 in the starburst of NGC 253

We present measurements of several near-infrared emission lines from the nearby galaxy NGC 253. We have been able to measure four H₂ lines across the circumnuclear starburst, from which we estimate the ortho- to para- ratio of excited H₂ to be ∼2. This indicates that the bulk of the H₂ emission arises from photodissociation regions (PDRs), rather than from shocks. This is the case across the entire region of active star formation.
As the H₂ emission arises from PDRs, it is likely that the ratio of H₂ to Brγ (the bright hydrogen recombination line) is a measure of the relative geometry of O and B stars and PDRs. Towards the nucleus of NGC 253 the geometry is deduced to be tightly clustered O and B stars in a few giant H  II regions that are encompassed by PDRs. Away from the nuclear region, the geometry becomes that of PDRs bathed in a relatively diffuse ultraviolet radiation field.
The rotation curves of 1–0 S(1) and Brγ suggest that the ionized gas is tracing a kinetic system different from that of the molecular gas in NGC 253, particularly away from the nucleus.     

A simple model for H2 line profiles in bow shocks

We present a model for empirically reproducing line profiles of molecular hydrogen emission in bow shocks. The model takes into account bow velocity, dissociation limit, a cooling function, viewing angle, bow shape and a limited form of extinction. Our results show that both geometrical factors and shock physics can significantly affect the profile morphology. In a companion paper we will apply this model to Fabry–Perot observations of bow shocks in the Orion BN–KL outflow.     

Shocked gas around Cepheus A: evidence for multiple outflows from H2S and SO2 observations

The Cepheus A star-forming region has been investigated through a multiline H₂S and SO₂ survey at millimetre wavelengths. Large-scale maps and high-resolution line profiles reveal the occurrence of several outflows. Cep A East is associated with multiple mass-loss processes: in particular, we detect a 0.6-pc jet-like structure which shows for the first time that the Cep A East young stellar objects are driving a collimated outflow moving towards the south.
The observed outflows show different clumps associated with definitely different H₂S/SO₂ integrated emission ratios, indicating that the gas chemistry in Cepheus A has been altered by the passage of shocks. H₂S appears to be more abundant than SO₂ in high-velocity clumps, in agreement with chemical models. However, we also find quite small H₂S linewidths, suggestive of regions where the evaporated H₂S molecules had enough time to slow down but not to freeze out on to dust grains. Finally, comparison between the line profiles indicates that the excitation conditions increase with the velocity, as expected for a propagation of collimated bow shocks.     

Testing the circumstellar disc hypothesis: a search for H2 outflow signatures from massive young stellar objects with linearly distributed methanol masers

The results of a survey searching for outflows using near-infrared imaging are presented. Targets were chosen from a compiled list of massive young stellar objects associated with methanol masers in linear distributions. Presently, it is a widely held belief that these methanol masers are found in (and delineate) circumstellar accretion discs around massive stars. If this scenario is correct, one way to test the disc hypothesis is to search for outflows perpendicular to the methanol maser distributions. The main objective of the survey was to obtain wide-field near-infrared images of the sites of linearly distributed methanol masers using a narrow-band 2.12-μm filter. This filter is centred on the  H₂ v = 1–0 S(1)  line; a shock diagnostic that has been shown to successfully trace CO outflows from young stellar objects. 28 sources in total were imaged of which 18 sources display H₂ emission. Of these, only two sources showed emission found to be dominantly perpendicular to the methanol maser distribution. Surprisingly, the H₂ emission in these fields is not distributed randomly, but instead the majority of sources are found to have H₂ emission dominantly parallel to their distribution of methanol masers. These results seriously question the hypothesis that methanol masers exist in circumstellar discs. The possibility that linearly distributed methanol masers are instead directly associated with outflows is discussed.     

H3+ in diffuse interstellar gas

A model is constructed of the material in front of the star Cygnus OB2 no. 12 in which dense cores are embedded in diffuse clumps of gas. The model reproduces the measured abundances of C₂ and CO, and predicts a column density of 910¹⁰ cm² for HCO⁺.     

The effect of grain size distribution on H2 formation rate in the interstellar medium

The formation of molecular hydrogen  (H₂)  in the interstellar medium takes place on the surfaces of dust grains. Hydrogen molecules play a role in gas-phase reactions that produce other molecules, some of which serve as coolants during gravitational collapse and star formation. Thus, the evaluation of the production rate of hydrogen molecules and its dependence on the physical conditions in the cloud are of great importance. Interstellar dust grains exhibit a broad size distribution in which the small grains capture most of the surface area. Recent studies have shown that the production efficiency strongly depends on the grain composition and temperature as well as on its size. In this paper, we present a formula that provides the total production rate of  H₂  per unit volume in the cloud, taking into account the grain composition and temperature as well as the grain size distribution. The formula agrees very well with the master equation results. It shows that for a physically relevant range of grain temperatures, the production rate of  H₂  is significantly enhanced due to their broad size distribution.     

Formation of H2 on an olivine surface: a computational study

The formation of H₂ on a pristine olivine surface [forsterite (010)] is investigated computationally. Calculations show that the forsterite surface catalyzes H₂ formation by providing chemisorption sites for H atoms. The chemisorption route allows for stepwise release of the reaction exothermicity and stronger coupling to the surface, which increases the efficiency of energy dissipation. This suggests that H₂ formed on a pristine olivine surface should be much less rovibrationally excited than H₂ formed on a graphite surface. Gas-phase H atoms impinging on the surface will first physisorb relatively strongly  ( E _phys= 1240 K)  . The H atom can then migrate via desorption and re-adsorption, with a barrier equal to the adsorption energy. The barrier for a physisorbed H atom to become chemisorbed is equal to the physisorption energy, therefore there is almost no gas-phase barrier to chemisorption. An impinging gas-phase H atom can easily chemisorb  ( E _chem= 12 200 K)  , creating a defect where a silicate O atom is protonated and a single electron resides on the surface above the adjacent magnesium ion. This defect directs any subsequent impinging H atoms to chemisorb strongly (39 800 K) on the surface electron site. The two adjacent chemisorbed atoms can subsequently recombine to form H₂ via a barrier (5610 K) that is lower than the chemisorption energy of the second H atom. Alternatively, the adsorbed surface species can react with another incoming H atom to yield H₂ and regenerate the surface electron site. This double chemisorption 'relay mechanism' catalyzes H₂ formation on the olivine surface and is expected to attenuate the rovibrational excitation of H₂ thus formed.