Shocked molecular gas towards the supernova remnant G359.1–0.5 and the Snake

We have found a bar of shocked molecular hydrogen (H₂) towards the OH(1720 MHz) maser located at the projected intersection of supernova remnant (SNR)  G359.1–0.5  and the non-thermal radio filament known as the Snake. The H₂ bar is well aligned with the SNR shell and almost perpendicular to the Snake. The OH(1720 MHz) maser is located inside the sharp western edge of the H₂ emission, which is consistent with the scenario in which the SNR drives a shock into a molecular cloud at that location. The spectral line profiles of ¹²CO, HCO⁺ and CS towards the maser show broad-line absorption, which is absent in the ¹³CO spectra and most probably originates from the pre-shock gas. A density gradient is present across the region and is consistent with the passage of the SNR shock, while the H₂ filament is located at the boundary between the pre-shock and post-shock regions.     

Chemical evolution ahead of Herbig–Haro objects

Compact regions of enhanced HCO⁺ and NH₃ emission have been detected close to a number of Herbig–Haro objects. An interpretation of these detections is the following: a transient clump within the molecular cloud has been irradiated by the shock that generates the Herbig–Haro object. The irradiation releases icy mantles from the grains within the transient clump and initiates a photochemistry. On the basis of this picture, we have developed an extensive chemical model which predicts that a wide range of species, other than NH₃ and HCO⁺, should also be detectable. These include CH₃OH, H₂S, C₃H₄, H₂CO, SO, SO₂, H₂CS and NS. The chemical effects should last ∼  10⁴ yr  .     

A source of extended HCO+ emission in young stellar objects

Anomalous molecular line profile shapes are the strongest indicators of the presence of the infall of gas that is associated with star formation. Such profiles are seen for well-known tracers, such as HCO⁺, CS and H₂CO. In certain cases, optically thick emission lines with appropriate excitation criteria may possess the asymmetric double-peaked profiles that are characteristic of infall. However, recent interpretations of the HCO⁺ infall profile observed towards the protostellar infall candidate B335 have revealed a significant discrepancy between the inferred overall column density of the molecule and that which is predicted by standard dark cloud chemical modelling.
This paper presents a model for the source of the HCO⁺ emission excess. Observations have shown that, in low-mass star-forming regions, the collapse process is invariably accompanied by the presence of collimated outflows; we therefore propose the presence of an interface region around the outflow in which the chemistry is enriched by the action of jets. This hypothesis suggests that the line profiles of HCO⁺, as well as other molecular species, may require a more complex interpretation than can be provided by simple, chemically quiescent, spherically symmetric infall models.
The enhancement of HCO⁺ depends primarily on the presence of a shock-generated radiation field in the interface. Plausible estimates of the radiation intensity imply molecular abundances that are consistent with those observed. Further, high-resolution observations of an infall-outflow source show HCO⁺ emission morphology that is consistent with that predicted by this model.     

An ultra-high-resolution study of the interstellar medium towards Orion

We report ultra-high-resolution observations     of Na  i , Ca  ii , K  i , CH and CH⁺ for interstellar sightlines towards 12 bright stars in Orion. These data enable the detection of many more absorption components than previously recognized, providing a more accurate perspective on the absorbing medium. This is especially so for the line of sight to the Orion nebula, a region not previously studied at very high resolution. Model fits have been constructed for the absorption-line profiles, providing estimates for the column density, velocity dispersion and central velocity for each constituent velocity component. A comparison between the absorption occurring in sightlines with small angular separations has been used, along with comparisons with other studies, to estimate the line-of-sight velocity structures. Comparisons with earlier studies have also revealed temporal variability in the absorption-line profile of ζ Ori, highlighting the presence of small-scale spatial structure in the interstellar medium on scales of ≈10 au. Where absorption from both Na⁰ and K⁰ is observed for a particular cloud, a comparison of the velocity dispersions measured for each of these species provides rigorous limits on both the kinetic temperature and turbulent velocity prevailing in each cloud. Our results indicate the turbulent motions to be subsonic in each case.     abundance ratios are derived for individual clouds, providing an indication of their physical state.     

Molecular hydrogen in damped Lyα systems: clues to interstellar physics at high redshift

In order to interpret H₂ quasar absorption-line observations of damped Lyα systems (DLAs) and subDLAs, we model their H₂ abundance as a function of dust-to-gas ratio, including H₂ self-shielding and dust extinction against dissociating photons. Then, we constrain the physical state of the gas by using H₂ data. Using H₂ excitation data for DLAs with H₂ detections, we derive a gas density  1.5 ≲ log n (cm⁻³) ≲ 2.5  , temperature  1.5 ≲ log T (K) ≲ 3  , and an internal ultraviolet (UV) radiation field (in units of the Galactic value)  0.5 ≲ log χ≲ 1.5  . We then find that the observed relation between the molecular fraction and the dust-to-gas ratio of the sample is naturally explained by the above conditions. However, it is still possible that H₂ deficient DLAs and subDLAs with H₂ fractions less than  ∼10⁻⁶  are in a more diffuse and warmer state. The efficient photodissociation by the internal UV radiation field explains the extremely small H₂ fraction  (≲10⁻⁶)  observed for  κ≲ 1/30  (κ is the dust-to-gas ratio in units of the Galactic value); H₂ self-shielding causes a rapid increase in, and large variations of, H₂ abundance for  κ≳ 1/30  . We finally propose an independent method to estimate the star formation rates of DLAs from H₂ abundances; such rates are then critically compared with those derived from other proposed methods. The implications for the contribution of DLAs to the cosmic star formation history are briefly discussed.     

New, high-resolution, near-infrared observations of HH 1

We present new, high-resolution, near-infrared images of the HH 1 jet and bow shock. H₂ and [Fe  ii ] images are combined to trace excitation changes along the jet and across the many shock features in this flow. Echelle spectra of H₂ profiles towards a few locations in HH 1 are also discussed. Gas excitation in oblique, planar C-type shocks best explains the observations, although J-type shocks must be responsible for the observed [Fe  ii ] emission features. Clearly, no single shock model can account for all of the observations. This will probably be true of most, if not all, Herbig–Haro flows.     

The thermal balance of the first structures in the primordial gas

The contraction of matter in the primordial medium, to form the first gravitationally bound structures, was mediated by radiative cooling of the gas by H₂ and HD. We have computed the initial phases of free-fall collapse, incorporating the results of quantum mechanical calculations of rate coefficients for collisional excitation of H₂ and HD by the principal perturbers, H, He, H₂ and H⁺. The structure of shock waves produced when the collapse speed exceeds the local sound speed is determined. In the post-shock gas, radiative cooling by H₂ exceeds that by HD, but by a factor of only 4. The intensities of the strongest emission lines of H₂– rotational transitions within the vibrational ground state – are calculated. Even with coarse spectral and angular resolution, these transitions might be observable as inhomogeneities in the cosmic background radiation.     

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.     

Excitation and kinematics in H2 bow shocks: near-infrared observations of HH 99 and VLA 1623A (HH 313)

We present a comprehensive near-infrared study of two molecular bow shocks in two protostellar outflows, HH 99 in R Coronae Australis and VLA 1623A (HH 313) in Rho Ophiuchi. New, high-resolution, narrow-band images reveal the well-defined bow shock morphologies of both sources. These are compared with two-dimensional MHD modelling of molecular bows from which we infer flow inclination angles, shock speeds and the magnetic field in the pre-shock gas in each system. With combined echelle spectroscopy and low-resolution K -band spectra we further examine the kinematics and excitation of each source. Bow shock models are used to interpret excitation (CDR) diagrams and estimate the extinction and, in the case of VLA 1623, the ortho–para ratio associated with the observed H₂ population. For the first time, morphology, excitation and kinematics are fitted with a single bow shock model.
Specifically, we find that HH 99 is best fitted by a C-type bow shock model (although a J-type cap is probably responsible for the [Fe  ii ] emission). The bow is flowing away from the observer (at an angle to the line of sight of ∼45°) at a speed of roughly 100 km s⁻¹. VLA 1623A is interpreted in terms of a C-type bow moving towards the observer (at an angle to the line of sight of ∼75°) at a speed of ∼80 km s⁻¹. The magnetic field associated with HH 99 is thought to be orientated parallel to the flow axis; in VLA 1623A the field is probably oblique to the flow axis, since this source is clearly asymmetric in our H₂ images.     

A textbook case of bow shock entrainment in a YSO outflow

Near-infrared images in H₂ line emission and submillimetre maps in CO J  = 3–2 emission illustrate the remarkable association between a molecular bow shock and the redshifted molecular outflow lobe in W75N. The flow lobe fits perfectly into the wake of the bow, as one would expect if the lobe represented swept-up gas. Indeed, these observations strongly support the 'bow shock' entrainment scenario for molecular outflows driven by young stars.   The characteristics of the bow shock and CO outflow lobe are compared with those of numerical simulations of jet-driven flows. These models successfully reproduce the bulge and limb-brightening in the CO outflow, although the model H₂ bow exhibits more structure extending back along the flow axis. We also find that the size of the flow, the high mass fraction in the flow at low outflow velocities (low γ values) and the high CO/H₂ luminosity ratio indicate that the system is evolved. We also predict a correlation, in evolved systems, between outflow age and the CO/H₂ luminosity ratio.     

Molecular hydrogen formation on porous dust grains

Recent laboratory experiments on interstellar dust analogues have shown that H₂ formation on dust-grain surfaces is efficient in a range of grain temperatures below 20 K. These results indicate that surface processes may account for the observed H₂ abundance in cold diffuse and dense clouds. However, high abundances of H₂ have also been observed in warmer clouds, including photon-dominated regions (PDRs), where grain temperatures may reach 50 K, making the surface processes extremely inefficient. It was suggested that this apparent discrepancy can be resolved by chemisorption sites. However, recent experiments indicate that chemisorption processes may not be efficient at PDR temperatures. Here we consider the effect of grain porosity on H₂ formation, and analyse it using a rate-equation model. It is found that porosity extends the efficiency of the recombination process to higher temperatures. This is because H atoms that desorb from the internal surfaces of the pores may re-adsorb many times and thus stay longer on the surface. However, this porosity-driven extension may enable efficient H₂ formation in PDRs only if porosity also contributes to significant cooling of the grains, compared to non-porous grains.     

Probing the surfaces of interstellar dust grains: the adsorption of CO at bare grain surfaces

A solid-state feature was detected at around 2175 cm⁻¹ towards 30 embedded young stellar objects in spectra obtained using the Infrared Spectrometer and Array Camera at the European Southern Observatory Very Large Telescope. We present results from laboratory studies of CO adsorbed at the surface of zeolite wafers, where absorption bands were detected at 2177 and 2168 cm⁻¹ (corresponding to CO chemisorbed at the zeolite surface) and 2130 cm⁻¹ (corresponding to CO physisorbed at the zeolite surface), providing an excellent match to the observational data. We propose that the main carrier of the 2175-band is CO chemisorbed at bare surfaces of dust grains in the interstellar medium. This result provides the first direct evidence that gas–surface interactions do not have to result in the formation of ice mantles on interstellar dust. The strength of the 2175-band is estimated to be  ∼4 × 10⁻¹⁹ cm  molecule⁻¹. The abundance of CO adsorbed at bare grain surfaces ranges from 0.06 to 0.16 relative to H₂O ice, which is, at most, half of the abundance (relative to H₂O ice) of CO residing in H₂O-dominated ice environments. These findings imply that interstellar grains have a large (catalytically active) surface area, providing a refuge for interstellar species. Consequently, the potential exists for heterogeneous chemistry to occur involving CO molecules in unique surface chemistry pathways not currently considered in gas grain models of the interstellar medium.     

The excitation and kinematical properties of H2 and [Fe ii] in the HH 46/47 bipolar outflow

Long-slit spectra of the molecular outflow Herbig–Haro (HH) 46/47 have been taken in the J and K near-infrared bands. The observed H₂ line emission confirms the existence of a bright and extended redshifted counter-jet outflow south-west of HH 46. In contrast with the optical appearance of this object, we show that this outflow seems to be composed of two different emission regions characterized by distinct heliocentric velocities. This implies an acceleration of the counter-jet.
The observed [Fe  ii ] emission suggests an average extinction of 7–9 visual magnitudes for the region associated with the counter-jet.
Through position–velocity diagrams, we show the existence of different morphologies for the H₂ and [Fe  ii ] emission regions in the northern part of the HH 46/47 outflow. We have detected for the first time high-velocity (−250 km s⁻¹) [Fe  ii ] emission in the region bridging HH 46 to HH 47A. The two strong peaks detected can be identified with the optical positions B8 and HH 47B.
The H₂ excitation diagrams for the counter-jet shock suggest an excitation temperature for the gas of T _ex≈2600 K . The lack of emission from the higher energy H₂ lines, such as the 4–3 S(3) transition, suggests a thermal excitation scenario for the origin of the observed emission. Comparison of the H₂ line ratios with various shock models yielded useful constraints about the geometry and type of these shocks. Planar shocks can be ruled out whereas curved or bow shocks (both J- and C-type) can be parametrized to fit our data.     

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 possibility of observing H2 emission from primordial molecular cloud kernels

We study the prospects for observing H₂ emission during the assembly of primordial molecular cloud kernels. The primordial molecular cloud cores, which resemble those at the present epoch, can emerge around  1+ z ∼20  according to recent numerical simulations. The kernels form inside the cores, and the first stars will appear inside the kernels. A kernel typically contracts to form one of the first generation stars with an accretion rate that is as large as ∼0.01 M_⊙ yr⁻¹. This occurs owing to the primordial abundances, which result in a kernel temperature of order 1000 K, and the collapsing kernel emits H₂ line radiation at a rate ∼10³⁵ erg s⁻¹. Predominantly   J =5-3   ( v =0)  rotational emission of H₂ is expected. At redshift  1+ z ∼20  , the expected flux is ∼0.01 μJy for a single kernel. While an individual object is not observable by any facilities available in the near future, the expected assembly of primordial star clusters on subgalactic scales can result in fluxes at the sub-mJy level. This is marginally observable with ASTRO-F and ALMA. We also examine the rotational   J =2-0   ( v =0)  and vibrational   δv =1  emission lines. The former may possibly be detectable with ALMA.     

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

On the detection of interstellar homonuclear diatomic molecules

The pure rotational spectrum of homonuclear diatomic molecules in the interstellar medium is strongly forbidden, and no such spectrum has been detected. In regions of high excitation, vibrational emission may occur, as is widely detected in the case of H₂ in interstellar shocks and photon-dominated regions. However, it is of considerable interest to know the abundance of homonuclear diatomics in quiescent regions. We propose that vibrational emission from homonuclear diatomic molecules in cold clouds may be detectable, where the excitation is mainly through collisions with non-thermal electrons arising from the cosmic-ray ionization of H₂. As an example, we estimate the intensity of emission from N₂ in cold, dark interstellar clouds. We show that such emission is at the limit of detectability with current technology. Other excitation mechanisms may also contribute and enhance this emission.     

Thermal desorption of water ice in the interstellar medium

Water (H₂O) ice is an important solid constituent of many astrophysical environments. To comprehend the role of such ices in the chemistry and evolution of dense molecular clouds and comets, it is necessary to understand the freeze-out, potential surface reactivity and desorption mechanisms of such molecular systems. Consequently, there is a real need from within the astronomical modelling community for accurate empirical molecular data pertaining to these processes. Here we give the first results of a laboratory programme to provide such data. Measurements of the thermal desorption of H₂O ice, under interstellar conditions, are presented. For ice deposited under conditions that realistically mimic those in a dense molecular cloud, the thermal desorption of thin films (≪50 molecular layers) is found to occur with zeroth-order kinetics characterized by a surface binding energy, E _des, of 5773 ± 60 K, and a pre-exponential factor, A , of 10^30 ± 2 molecules cm⁻² s⁻¹. These results imply that, in the dense interstellar medium, thermal desorption of H₂O ice will occur at significantly higher temperatures than has previously been assumed.     

Interaction between the north-eastern boundary of Sgr A East and giant molecular clouds

We have detected the   v = 1 → 0 S(1) (λ= 2.1218 μm)  and   v = 2 → 1 S(1) (λ= 2.2477 μm)  lines of H₂ in the Galactic Centre, in a  90 × 27 arcsec²  region between the north-eastern boundary of the non-thermal source Sgr A East, and the giant molecular cloud (GMC)  M−0.02 − 0.07  . The detected  H₂ v = 1 → 0  S(1) emission has an intensity of  1.6–21 × 10⁻¹⁸ W m⁻² arcsec⁻²  and is present over most of the region. Along with the high intensity, the large linewidths  (FWHM = 40–70 km s⁻¹)  and the  H₂ v = 2 → 1 S(1)  to   v = 1 → 0 S(1)  line ratios (0.3–0.5) can be best explained by a combination of C-type shocks and fluorescence. The detection of shocked H₂ is clear evidence that Sgr A East is driving material into the surrounding adjacent cool molecular gas. The H₂ emission lines have two velocity components at ∼+50 and  ∼0 km s⁻¹  , which are also present in the NH₃(3, 3) emission mapped by McGary, Coil & Ho. This two-velocity structure can be explained if Sgr A East is driving C-type shocks into both the  GMC M−0.02 − 0.07  and the northern ridge of McGary et al.